Home About us Contact
|
Home About us Contact | ||
|
|
||
Singular Kernel (singular + kernel)
Selected AbstractsSingularity extraction technique for integral equation methods with higher order basis functions on plane triangles and tetrahedraINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2003Seppo Järvenpää Abstract A numerical solution of integral equations typically requires calculation of integrals with singular kernels. The integration of singular terms can be considered either by purely numerical techniques, e.g. Duffy's method, polar co-ordinate transformation, or by singularity extraction. In the latter method the extracted singular integral is calculated in closed form and the remaining integral is calculated numerically. This method has been well established for linear and constant shape functions. In this paper we extend the method for polynomial shape functions of arbitrary order. We present recursive formulas by which we can extract any number of terms from the singular kernel defined by the fundamental solution of the Helmholtz equation, or its gradient, and integrate the extracted terms times a polynomial shape function in closed form over plane triangles or tetrahedra. The presented formulas generalize the singularity extraction technique for surface and volume integral equation methods with high-order basis functions. Numerical experiments show that the developed method leads to a more accurate and robust integration scheme, and in many cases also a faster method than, for example, Duffy's transformation. Copyright © 2003 John Wiley & Sons, Ltd. [source] On the spectrum of the electric field integral equation and the convergence of the moment methodINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 1 2001Karl F. Warnick Abstract Existing convergence estimates for numerical scattering methods based on boundary integral equations are asymptotic in the limit of vanishing discretization length, and break down as the electrical size of the problem grows. In order to analyse the efficiency and accuracy of numerical methods for the large scattering problems of interest in computational electromagnetics, we study the spectrum of the electric field integral equation (EFIE) for an infinite, conducting strip for both the TM (weakly singular kernel) and TE polarizations (hypersingular kernel). Due to the self-coupling of surface wave modes, the condition number of the discretized integral equation increases as the square root of the electrical size of the strip for both polarizations. From the spectrum of the EFIE, the solution error introduced by discretization of the integral equation can also be estimated. Away from the edge singularities of the solution, the error is second order in the discretization length for low-order bases with exact integration of matrix elements, and is first order if an approximate quadrature rule is employed. Comparison with numerical results demonstrates the validity of these condition number and solution error estimates. The spectral theory offers insights into the behaviour of numerical methods commonly observed in computational electromagnetics. Copyright © 2001 John Wiley & Sons, Ltd. [source] Blow up for a Cauchy viscoelastic problem with a nonlinear dissipation of cubic convolution typeMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 15 2009Shengqi Yu Abstract In this paper, we consider a Cauchy viscoelastic problem with a nonlinear source of polynomial type and a nonlinear dissipation of cubic convolution type involving a singular kernel. Under suitable conditions on the initial data and the relaxation functions, it is proved that the solution of this particular problem blows up in finite time. Copyright © 2009 John Wiley & Sons, Ltd. [source] A new integral equation approach to the Neumann problem in acoustic scatteringMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 16 2001P. A. Krutitskii We suggest a new approach of reduction of the Neumann problem in acoustic scattering to a uniquely solvable Fredholm integral equation of the second kind with weakly singular kernel. To derive this equation we placed an additional boundary with an appropriate boundary condition inside the scatterer. The solution of the problem is obtained in the form of a single layer potential on the whole boundary. The density in the potential satisfies a uniquely solvable Fredholm integral equation of the second kind and can be computed by standard codes. Copyright © 2001 John Wiley & Sons, Ltd. [source] Electrostatic BEM for MEMS with thin beamsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 6 2005Zhongping Bao Abstract Micro-electro-mechanical (MEM) and nano-electro-mechanical (NEM) systems sometimes use beam- or plate-shaped conductors that can be very thin,with h/L,,,(10,2,10,3) (in terms of the thickness h and length L of a beam or the side of a square pate). Conventional boundary element method (BEM) analysis of the electric field in a region exterior to such thin conductors can become difficult to carry out accurately and efficiently,especially since MEMS analysis requires computation of charge densities (and then surface tractions) separately on the top and bottom surfaces of such objects. A new boundary integral equation (BIE) is derived in this work that, when used together with the standard BIE with logarithmically singular kernels, results in a powerful technique for the BEM analysis of such problems with thin beams. This new approach, in fact, works best for very thin beams. This thin beam BEM is derived and discussed in this work. Copyright © 2005 John Wiley & Sons, Ltd. [source] Boundary element formulation for 3D transversely isotropic cracked bodiesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2004M. P. Ariza Abstract The boundary traction integral representation is obtained in elasticity when the classical displacement representation is differentiated and combined according to Hooke's law. The use of both traction and displacement integral representations leads to a mixed (or dual) formulation of the BEM where the discretization effort for crack problems is much smaller than in the classical formulation. A boundary element analysis of three-dimensional fracture mechanics problems of transversely isotropic solids based on the mixed formulation is presented in this paper. The hypersingular and strongly singular kernels appearing in the formulation are regularized by using two terms of the displacement series expansion and one term of the traction expansion, at the collocation point. All the remaining integrals are analytically evaluated or transformed by means of Stokes' theorem into regular or weakly singular integrals, which are numerically computed. The method is general and can be used for elements of any shape including quarter-point crack front elements. No change of co-ordinates is required for the integration. The formulation as presented in this paper is something as clear, general and easy to handle as the classical BE formulation. It is used in combination with three-dimensional quadratic and quarter-point elements to obtain accurate results for several different crack problems. Cracks in boundless and finite transversely isotropic domains are studied. The meshes are simple and include only discretization of the crack and the external boundary. The obtained results are in good agreement with those existing in the literature. Copyright © 2004 John Wiley & Sons, Ltd. [source] Singularity extraction technique for integral equation methods with higher order basis functions on plane triangles and tetrahedraINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2003Seppo Järvenpää Abstract A numerical solution of integral equations typically requires calculation of integrals with singular kernels. The integration of singular terms can be considered either by purely numerical techniques, e.g. Duffy's method, polar co-ordinate transformation, or by singularity extraction. In the latter method the extracted singular integral is calculated in closed form and the remaining integral is calculated numerically. This method has been well established for linear and constant shape functions. In this paper we extend the method for polynomial shape functions of arbitrary order. We present recursive formulas by which we can extract any number of terms from the singular kernel defined by the fundamental solution of the Helmholtz equation, or its gradient, and integrate the extracted terms times a polynomial shape function in closed form over plane triangles or tetrahedra. The presented formulas generalize the singularity extraction technique for surface and volume integral equation methods with high-order basis functions. Numerical experiments show that the developed method leads to a more accurate and robust integration scheme, and in many cases also a faster method than, for example, Duffy's transformation. Copyright © 2003 John Wiley & Sons, Ltd. [source] | ||||||||||