A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

Solution of Poisson’s equation by analytical boundary element integration

The solution of Poisson’s equation is essential for many branches of science and engineering such as fluid-mechanics, geosciences, and electrostatics. Solution of two-dimensional Poisson’s equations is carried out by BEM based on Galerkin Vector Method in which the integrals that appear in the boundary element method are expressed by analytical integration. In this paper, the Galerkin vector method is developed for more general case than presented in the previous works. The integrals are computed for constant and linear elements in BEM. By employing analytical integration in BEM computation, the numerical schemes and coordinate transformations can be avoided. The presented method can also be used for the multiple domain case. The results of the analytical integration are employed in BEM code and the obtained analytical expression will be applied to several examples where the exact solution exists. The produced results are in good agreement with the exact solution.     

An analytical contour integration method for handling body forces in elasticity

In the boundary element and superposition methods for two-dimensional elasticity problems, the presence of body forces requires the integration of the basic point load solution against the body force field over a specified area. For polynomial body forces, x^myⁿ, these area integrals may be transformed into contour integrals which have closed form analytic representations when the area is a polygon. This paper shows that commonly used numerical integration techniques developed for area integrals are unsuitable for handling body forces using the above named methods. Not only are these techniques inaccurate but expensive in terms of execution time in comparison to the analytical algorithm developed here. A computer program and four example problems are included.     

Numerical integration of logarithmic and nearly logarithmic singularity in BEMs

A new semi-analytical integration scheme is proposed for evaluation of logarithmically singular and/or nearly singular integrals occurring in 2D BEM formulations. Extensive numerical experiments are performed to study the accuracy by the proposed and other schemes of numerical integration. The accuracy of the numerical integration is almost exact even for curvilinear elements and the accuracy of numerical computation is determined by the accuracy of the approximation of the boundary density and geometry.     

A general algorithm for the numerical evaluation of nearly singular integrals on 3D boundary element

A general numerical method is proposed to compute nearly singular integrals arising in the boundary integral equations (BIEs). The method provides a new implementation of the conventional distance transformation technique to make the result stable and accurate no matter where the projection point is located. The distance functions are redefined in two local coordinate systems. A new system denoted as (α,β) is introduced here firstly. Its implementation is simpler than that of the polar system and it also performs efficiently. Then a new distance transformation is developed to remove or weaken the near singularities. To perform integration on irregular elements, an adaptive integration scheme is applied. Numerical examples are presented for both planar and curved surface elements. The results demonstrate that our method can provide accurate results even when the source point is very close to the integration element, and can keep reasonable accuracy on very irregular elements. Furthermore, the accuracy of our method is much less sensitive to the position of the projection point than the conventional method.     

Some algorithms of the projective-Newton method for the finite element space

The paper presents the projective-Newton method in the finite element space with curved finite elements along the boundary. The properties of this method are studied when the integrals are calculated by numerical quadrature. Some algorithms are derived when ω = {x ? R²: 0 < x₁, x₂ < 1}.     

Evaluation of regular and singular domain integrals with boundary-only discretization—theory and Fortran code

In this paper, a set of boundary integrals are derived based on a radial integration technique to accurately evaluate two dimensional (2D) and three dimensional (3D), regular and singular domain integrals. A self-contained Fortran code is listed and described for numerical implementation of these boundary integrals. The main feature of the theory is that only the boundary of the integration domain needs to be discretized into elements. This feature cannot only save considerable efforts in discretizing the integration domain into internal cells (as in the conventional method), but also make computational results for singular domain integrals more accurate since the integrals have been regularized. Some examples are provided to verify the correctness of the presented formulations and the included code.     

On numerical cubatures of singular surface integrals in boundary element methods

Summary We present and analyze methods for the accurate and efficient evaluation of weakly, Cauchy and hypersingular integrals over piecewise analytic curved surfaces in ³.The class of admissible integrands includes all kernels arising in the numerical solution of elliptic boundary value problems in three-dimensional domains by the boundary integral equation method. The possibly not absolutely integrable kernels of boundary integral operators in local coordinates are pseudohomogeneous with analytic characteristics depending on the local geometry of the surface at the source point. This rules out weighted quadrature approaches with a fixed singular weight.For weakly singular integrals it is shown that Duffy's triangular coordinates leadalways to a removal of the kernel singularity. Also asymptotic estimates of the integration error are provided as the size of the boundary element patch tends to zero. These are based on the Rabinowitz-Richter estimates in connection with an asymptotic estimate of domains of analyticity in ².It is further shown that the modified extrapolation approach due to Lyness is in the weakly singular case always applicable. Corresponding error and asymptotic work estimates are presented.For the weakly singular as well as for Cauchy and hypersingular integrals which e.g. arise in the study of crack problems we analyze a family of product integration rules in local polar coordinates. These rules are hierarchically constructed from finite part integration formulas in radial and Gaussian formulas in angular direction. Again, we show how the Rabinowitz-Richter estimates can be applied providing asymptotic error estimates in terms of orders of the boundary element size.Partially supported by the Priority Research Programme Boundary Element Methods of the German Research Foundation DFG under Grant No. We 659/16-1 (guest programme) and under AFOSR-grant 89-0252.     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

A product quadrature algorithm by Hermite interpolation

This paper is concerned with numerical integration of ∫¹₋₁f(x)k(x)dx by product integration rules based on Hermite interpolation. Special attention is given to the kernel k(x) = e^iτx, with a view to providing high precision rules for oscillatory integrals. Convergence results and error estimates are obtained in the case where the points of integration are zeros of p_n(W; x) or of (1 − x²)p_n−2(W; x), where p_n(W; x), n = 0, 1, 2…, are the orthonormal polynomials associated with a generalized Jacobi weight W. Further, examples are given that test the performance of the algorithm for oscillatory weight functions.     

Extrapolation methods to compute hypersingular integral in boundary element methods

The composite trapezoidal rule for the computation of Hadamard finite-part integrals in boundary element methods with the hypersingular kernel 1/sin 2(x-s) is discussed,and the main part of the asymptotic expansion of error function is obtained.Based on the main part of the asymptotic expansion,a series is constructed to approach the singular point.An extrapolation algorithm is presented and the convergence rate is proved.Some numerical results are also presented to confirm the theoretical results and show the efficiency of the algorithms.     

Semi-analytic integration of hypersingular Galerkin BIEs for three-dimensional potential problems

An accurate and efficient semi-analytic integration technique is developed for three-dimensional hypersingular boundary integral equations of potential theory. Investigated in the context of a Galerkin approach, surface integrals are defined as limits to the boundary and linear surface elements are employed to approximate the geometry and field variables on the boundary. In the inner integration procedure, all singular and non-singular integrals over a triangular boundary element are expressed exactly as analytic formulae over the edges of the integration triangle. In the outer integration scheme, closed-form expressions are obtained for the coincident case, wherein the divergent terms are identified explicitly and are shown to cancel with corresponding terms from the edge-adjacent case. The remaining surface integrals, containing only weak singularities, are carried out successfully by use of standard numerical cubatures. Sample problems are included to illustrate the performance and validity of the proposed algorithm.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method

A systematic treatment of the three-dimensional Poisson equation via singular and hypersingular boundary integral equation techniques is investigated in the context of a Galerkin approximation. Developed to conveniently deal with domain integrals without a volume-fitted mesh, the proposed method initially converts domain integrals featuring the Newton potential and its gradient into equivalent surface integrals. Then, the resulting boundary integrals are evaluated by means of well-established cubature methods. In this transformation, weakly-singular domain integrals, defined over simply- or multiply-connected domains with Lipschitz boundaries, are rigorously converted into weakly-singular surface integrals. Combined with the semi-analytic integration approach developed for potential problems to accurately calculate singular and hypersingular Galerkin surface integrals, this technique can be employed to effectively deal with mixed boundary-value problems without the need to partition the underlying domain into volume cells. Sample problems are included to validate the proposed approach.     

Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method

A systematic treatment of the three-dimensional Poisson equation via singular and hypersingular boundary integral equation techniques is investigated in the context of a Galerkin approximation. Developed to conveniently deal with domain integrals without a volume-fitted mesh, the proposed method initially converts domain integrals featuring the Newton potential and its gradient into equivalent surface integrals. Then, the resulting boundary integrals are evaluated by means of well-established cubature methods. In this transformation, weakly-singular domain integrals, defined over simply- or multiply-connected domains with Lipschitz boundaries, are rigorously converted into weakly-singular surface integrals. Combined with the semi-analytic integration approach developed for potential problems to accurately calculate singular and hypersingular Galerkin surface integrals, this technique can be employed to effectively deal with mixed boundary-value problems without the need to partition the underlying domain into volume cells. Sample problems are included to validate the proposed approach.     

Interface problem in holonomic elastoplasticity

The three-dimensional interface problem with the homogeneous Lamé system in an unbounded exterior domain and holonomic material behaviour in a bounded interior Lipschitz domain is considered. Existence and uniqueness of solutions of the interface problem are obtained rewriting the exterior problem in terms of boundary integral operators following the symmetric coupling procedure. The numerical approximation of the solutions consists in coupling of the boundary element method (BEM) and the finite element method (FEM). A Céa-like error estimate is presented for the discrete solutions of the numerical procedure proving its convergence.     

Singularity cancellation method for time-domain boundary element formulation of elastodynamics: A direct approach

In the implementation of time-domain boundary element method for elasto-dynamic problems, there are two types of singularities: the wave front singularity arising when the product of wave velocity and time is equal to the distance between the source point and the field point, and the spatial singularity arising when the source point coincides with the field point. In this paper, the singularity of the first type in the integrand is eliminated by an analytical integration over time, Cauchy principal value and Hadamard finite part integral. Four types of singularities with different orders appear in the integrand after analytical time integration. In order to accurately calculate the integral, in which the integrand is piecewise continuous, the integral domain is subdivided into several patches based on the relation between the product of wave velocity and time and the distance. In singular patches, the integrands are separated into a regular part and a singular part. The regular part can be computed by traditional numerical integration method such as Gaussian integration, while the singular part can be analytically integrated. Using the proposed method, the spatial singular integrals can be calculated directly. Numerical examples using various kinds of elements are presented to verify the proposed method.     

Clifford algebra valued boundary integral equations for three-dimensional elasticity

Applications of Clifford analysis to three-dimensional elasticity are addressed in the present paper. The governing equation for the displacement field is formulated in terms of the Dirac operator and Clifford algebra valued functions so that a general solution is obtained analytically in terms of one monogenic function and one multiple-component spatial harmonic function together with its derivative. In order to solve numerically the three-dimensional problems of elasticity for an arbitrary domain with complicated boundary conditions, Clifford algebra valued boundary integral equations (BIEs) for multiple-component spatial harmonic functions at an observation point, either inside the domain, on the boundary, or outside the domain, are constructed. Both smooth and non-smooth boundaries are considered in the construction. Moreover, the singularities of the integrals are evaluated exactly so that in the end singularity-free BIEs for the observation point on the boundary taking values on Clifford numbers can be obtained. A Clifford algebra valued boundary element method (BEM) based on the singularity-free BIEs is then developed for solving three-dimensional problems of elasticity. The accuracy of the Clifford algebra valued BEM is demonstrated numerically.