Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 9523405, 8 pages
http://dx.doi.org/10.1155/2016/9523405

Research Article

High Order Projection Plane Method for Evaluation of Supersingular Curved Boundary Integrals in BEM

Miao Cui, Wei-zhe Feng, Xiao-wei Gao, and Kai Yang

School of Aeronautics and Astronautics, Faculty of Vehicle Engineering and Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China

Received 7 December 2015; Revised 12 February 2016; Accepted 17 February 2016

Academic Editor: Ivano Benedetti

Copyright © 2016 Miao Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Boundary element method (BEM) is a very promising approach for solving various engineering problems, in which accurate evaluation of boundary integrals is required. In the present work, the direct method for evaluating singular curved boundary integrals is developed by considering the third-order derivatives in the projection plane method when expanding the geometry quantities at the field point as Taylor series. New analytical formulas are derived for geometry quantities defined on the curved line/plane, and unified expressions are obtained for both two-dimensional and three-dimensional problems. For the two-dimensional boundary integrals, analytical expressions for the third-order derivatives are derived and are employed to verify the complex-variable-differentiation method (CVDM) which is used to evaluate the high order derivatives for three-dimensional problems. A few numerical examples are given to show the effectiveness and the accuracy of the present method.

1. Introduction

Boundary element method (BEM) [1] is a very promising approach for solving various engineering problems, such as heat transfer problems [2–5] in energy science and engineering, in which accurate evaluation of boundary integrals [6–10] is required. Weak, strong, or hypersingularities are involved in these boundary integrals, if the source point is located on the element under integration [11]. These issues have been focused on by researchers, and many techniques, such as Gauss logarithmic quadrature rule for weak singularities [12], rigid body displacement algorithm for strong singularities [12], and regularization for hypersingularities [13], have been proposed to remove the singularities. Each method has its advantage and disadvantage [11–19]. A comprehensive review has been given in [11] on the treatment of singularities. In general, the efficient methods for singular curved line and surface integrals are those based on the operation over the intrinsic coordinate system [14–16].

Recently, Gao et al. [11] proposed another efficient direct method for numerically evaluating all kinds of singular curved boundary integrals for two-dimensional (2D) and three-dimensional (3D) BEM analyses, and numerical examples were given to demonstrate the stability of the proposed method. However, only quadratic terms were truncated when expanding the coordinate of the field point as Taylor series, which means the method will be invalid if the curved boundary is beyond the quadratic, such as the cubed curve.

In the present work, the direct method in [11] for evaluating arbitrary high order singular curved boundary integrals is developed, by considering the third-order derivatives when expanding the geometry quantities at the field point as Taylor series. Corresponding analytical expressions are derived. The first innovation is that the present method will extend the application range of the method in [11]. That is, the singular integrals over higher order elements, such as cubed curve, can be directly evaluated. Another innovation is that the complex-variable-differentiation method [20] is introduced to determine the third-order derivatives. For 2D boundary integral, the third-order derivatives are analytically derived, which are used to verify the CVDM.

The paper is organized as follows. In Section 2, review and modification of the direct method in [11] are provided. In Section 3, numerical examples are given to show the effectiveness and the accuracy of the present method. Finally, conclusions are drawn from computational results and discussions.

2. Development of the Direct Method for Evaluating Arbitrary High Order Singular Curved Boundary Integrals

2.1. Review of Singular Integrals [11]

Singular integrals over an element can be classified into the following form:

in which and represent the source and field points, respectively; is the boundary element under integration, which is a curved line or a curved surface for 2D or 3D problems, respectively; and is the distance between the source and field points.

The term in (1) is a regular function and is the order of the singularity.

In 2D problems, the following type of singular line integrals is also frequently encountered:

It is assumed that the integrals in (1) and (2) always exist: that is, each integration has a finite value.

2.2. Review of Expansions of Geometry Quantities on the Projection Line/Plane [11]

In [11], a projection line for 2D or a projection plane for 3D problems, which is the tangential line/plane of the element to the origin of the intrinsic coordinate system, was introduced. A local orthogonal coordinate system was established on the projection line/plane. The coordinate transformation between the local and global systems can be performed using the following relationships:

where is the direction cosine of the local coordinate axes with respect to the global one. The repeated subscripts represent summation; is the global coordinates of the origin of the local coordinate system.

The original curved element can be projected onto the projection line or plane to form a straight or a flat projection element by (3), and then all geometry quantities can be expressed in terms of variables defined on the projection line/plane.

2.3. High Order Expansion of Geometry Quantities in terms of Variables on the Projection Line for 2D Problems

For 2D and 3D problems, derivations for the modification are different, which are separately described in detail.

In 2D problems, the boundary element is a curved line as shown in Figure 1. The local orthogonal coordinate system is denoted by () and the projection line is along the axis .

Figure 1: The curved line boundary and the projection line.

There is only one independent variable for a line, and all other geometry quantities can be expressed by . Similar to the treatment in [11], the coordinate at the field point over the curved line can be expressed by using Taylor series. Different from the treatment in [11], the third-order derivative is considered in the present work:

where

It can be seen that the expressions are more complex. For 2D problem, the corresponding derivatives are analytically derived, and the detailed derivation process is omitted here.

In (5) and (6),

In (7a)–(8c),

Define a local distance:

where is the local coordinate of point which is the projection point of the source point and is the local coordinate of point which is the projection point of the field point .

The derivative of with respect to local coordinate is

In (12a), or .

Then, we can obtain

By substituting (12a) into (5), we can obtain

where

In (15),

Similarly, we can express the intrinsic coordinate in terms of the local distance as

where and can be determined by replacing with in (13)–(16c).

Substituting (12b) and (13) into (4), we can obtain

Then, substituting (15) into (18), we can obtain

where

Referring to Figure 1, the global distance can be easily expressed as

where

In (22), , , , , and .

It can be easily derived that

2.4. The Complex-Variable-Differentiation Method (CVDM)

From the above derivations, it can be seen that the corresponding expressions have been changed, by considering the third-order derivatives when expanding the geometry quantities at the field point as Taylor series. The third-order derivatives should be first determined.

In the present work, the complex-variable-differentiation method (CVDM) [20] is introduced to determine the high order derivatives. For 2D analysis, the third-order derivatives are determined explicitly by (7c), (8c), (9), (10a), and (10b), which are also calculated by the CVDM.

In the CVDM, the variable of a real function is replaced by a complex one , with the imaginary part being very small (usually = 10⁻²⁰). The function can be expanded into Taylor’s series as

Since is small enough, (25) can be obtained:

where Im denotes the imaginary part.

The CVDM is a very promising method, since the derivatives only require direct evaluations [21]. The CVDM has been validated to be an efficient and accurate approach for determination of derivatives [22–24].

For example, considering , we can add the very small imaginary part to ; then the real value of yields . In addition, the value of can be easily obtained using the analytical expressions in the previous work [11]. The CVDM and the derived expressions in (7c)–(8c) can be validated by each other.

2.5. Expansion of Geometry Quantities in terms of Variables on the Projection Plane for 3D Problems

In 3D problems, the boundary element is a curved surface marked as in Figure 2. A local orthogonal coordinate system () is established on the projection plane with its origin being at the point (, ), in which axes and are located within the plane. The axis is along direction and axis is along the outward normal direction to the element. By (4), we can project the original curved surface element onto the projection plane, to form a flat projection element marked as in Figure 2.

Figure 2: The curved surface boundary and the projection plane.

In a plane, there are only two independent variables. We choose and as the independent variables, and over the curved surface can be expressed in terms of and . Compared with the treatment in [11], is expanded as Taylor series by and , by considering the third-order derivatives:

where , , and are integers varying from 1 to 2, and

The treatment is the same as that in [11], to determine in (27a) and in (27b). For in (27c), we use the CVDM, to avoid rather complicated analytical derivations.

The local distance projected from the global distance onto the projection plane and its derivative are introduced as follows:

The local coordinates at the field point can be expressed as

Similar to the 2D analysis, substituting (29) into (26) yields

where

In (32),

Similar to the 2D analysis, intrinsic coordinates can be expressed as

where and can be determined by replacing with in (31)–(33c).

Substituting (29) and (30) into (4) yields

Substituting (32) into (35), we can obtain

where

The expression of is as follows:

where

In (39), , , , , and .

It can be easily derived that has the same expression as that in 2D analysis.

Considering that and or for 2D problems, it can be seen that unified expressions are obtained for both two-dimensional and three-dimensional problems.

Then, we can handle the singular boundary integrals in (1) and (2), with the above modifications. The 2D problem is taken as an example to describe the procedure to determine the integral in (1). The differential relationship between the local distance and the real curved line can be written as

where and are the outward normals to the tangential lines of passing the origin of the local coordinate system and the field point, respectively, and is the angle between the two outward normals.

Substituting (40) and (38) into (1) leads to

where is the distance from the projected source point to one of the ends of the projected line element and is the regular part: that is,

In order to integrate (41), the nonsingular part is expanded as a power series in , such that

in which is the order of the power series, usually taking a value between 2 and 7 depending on the size of , and are constants which are determined by collocating points over the integration region . In the paper, equally spaced points are used: that is, . The coefficient for the first point is and other coefficients can be solved using the following equation set:

where is a square matrix with the order of :

and are vectors as follows:

Solving (44) for coefficient vector and then substituting (43) into (41) yield

where

Then, singularities can be analytically removed.

3. Numerical Examples

Based on the Fortran code for singular integrals in [11], a new Fortran code has been developed with the above modification presented in this paper. It should be emphasized that the functions of the old Fortran code are fully included in the new Fortran code. If all the third-order derivatives are set to be zero, the functions are the same with those of the method in [11]. Besides, the new code has its unique characteristics.

For 2D analysis, the third-order derivatives are computed by both CVDM and analytical expression. For 3D analysis, the high order derivatives are computed by the CVDM, because analytical derivations are very difficult and complicated to be obtained.

In this section, the first example is given to show the advantage of CVDM over a difference method. Then, two singular integrals are calculated, to show the effectiveness and the accuracy of the present method.

3.1. The Advantage of the CVDM

A singular function is selected:

Then, the exact derivative with respect to the variable can be obtained:

For comparison, the central difference formulation is adopted: that is,

When , , the exact first-order derivative is 2 × 10⁹. The dependence of the derivatives on for the CVDM and the difference method is examined, and the results are listed in Table 1.

Table 1: The first-order derivatives with different h for .

CVDM
Difference method

It can be seen that the difference method can give satisfactory results only for between 10⁻⁵ and 10⁻¹⁰. For smaller , its result becomes unacceptable. However, the result of the CVDM tends to be independent of as the step size becomes sufficiently small. This clearly verifies the advantage of the CVDM.

3.2. A Weakly Singular Line Integral over a Cubed Curve

In this example, a singular line integral for a 2D BEM problem is chosen:

The integration line is a cubed curve as shown in Figure 3.

Figure 3: The cubed curve for integration.

The four typical points are 1 , 2 , 3 , and 4 . The computational point (source point) is at 1. The curve is taken as one boundary element. It should be emphasized here that the direct method in [11] has been invalid, which cannot be used for the singular line integration in (52) over the cubed curve in Figure 4, with only one element, because only quadratic terms were truncated when expanding the geometry quantities at the field point as Taylor series.

Figure 4: The cylindrical surface.

Firstly, the accuracy of the CVDM for the proposed technique is validated. We use the derived analytical expressions in (7c)–(8c) to compute the third-order derivatives. Then, the CVDM is employed to determine the same derivatives.

Using the derived analytical expressions in this paper, the third-order derivatives can be obtained:

The results of the CVDM are exactly the same with that in (53). It can be seen that the CVDM and the analytical expressions in (7c)–(8c) are validated by each other.

Secondly, the singular integral is evaluated using the present method. The result is 0.229935, with the Gauss order being 6.

This integral is also evaluated using the Gauss quadrature rule for logarithmic function [12]. The result is . It can be seen that the present method is effective and with high accuracy, for higher order curved boundary.

3.3. A Supersingular Integral over a Curved Surface

In this example, the accuracy of the present method in handling singular integrals over a highly curved surface is examined. The boundary element is a cylindrical surface with radius 1 and the rotational angle is . The geometrical dimensions are shown in Figure 4.

The boundary integral is as follows:

Two source points, and , are located at positions with intrinsic coordinates of and .

Guiggiani et al. [14] computed this integral for the case of . Table 2 lists the results using both methods. For the present method, the Gauss order is 4.

Table 2: Results for two points a and b ().

Point a
Point b

Method in [14]
−0.3438
−0.4971

Present method
−0.3458
−0.4981

From Table 2, it can be seen that he results using the present method agree well with the Guiggiani et al.’s results, which indicate the high accuracy of the present method for singular integrals over a highly curved surface.

4. Conclusion and Future Work

In the present work, high order projection plane method for evaluation of supersingular curved boundary integrals in BEM has been developed, by considering the third-order derivatives when expanding the geometry quantities at the field point as Taylor series. Corresponding analytical expressions are derived, and unified expressions are obtained for both two-dimensional and three-dimensional problems.

The present method extends the application range of the method in [11]. That is, the singular integrals over higher order elements, such as cubed curve, can be directly evaluated. The complex-variable-differentiation method is introduced to determine the third-order derivatives, which effectively avoids rather complicated analytical derivations for the three-dimensional analysis. For the two-dimensional boundary integral, the third-order derivatives are analytically derived and calculated, which are used to verify the complex-variable-differentiation method.

The results show that the present method is effective and with high accuracy for high order singular curved boundary integrals.

Future work is to apply the mathematical derivations and equations to evaluate the singularities involved in radial integration boundary element method in solving nonhomogeneous thermoelastic problems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the National Natural Science Foundation of China (51576026 and 11572075) and the Fundamental Research Funds for the Central Universities (DUT14LK03).

References

C. F. Rowlatt and T. N. Phillips, “A spectral element formulation of the immersed boundary method for Newtonian fluids,” Computer Methods in Applied Mechanics and Engineering, vol. 298, pp. 29–57, 2016.

H. Qi, L. M. Ruan, H. C. Zhang, Y. M. Wang, and H. P. Tan, “Inverse radiation analysis of a one-dimensional participating slab by stochastic particle swarm optimizer algorithm,” International Journal of Thermal Sciences, vol. 46, no. 7, pp. 649–661, 2007.

S. Zhang, Y. Zhu, B. Qiao, and Z. Zhang, “System model of heat and mass transfer process for mobile solvent vapor phase drying equipment,” Mathematical Problems in Engineering, vol. 2014, Article ID 267276, 11 pages, 2014.

T. Lu, W. W. Han, P. X. Jiang, Y. H. Zhu, J. Wu, and C. L. Liu, “A two-dimensional inverse heat conduction problem for simultaneous estimation of heat convection coefficient, fluid temperature and wall temperature on the inner wall of a pipeline,” Progress in Nuclear Energy, vol. 81, pp. 161–168, 2015.

H. Qi, Y. B. Qiao, S. C. Sun, Y. Yao, and L. Ruan, “Image reconstruction of two-dimensional highly scattering inhomogeneous medium using MAP-based estimation,” Mathematical Problems in Engineering, vol. 2015, Article ID 412315, 9 pages, 2015.

J. Sladek, V. Sladek, and S. N. Atluri, “Local boundary integral equation (LBIE) method for solving problems of elasticity with nonhomogeneous material properties,” Computational Mechanics, vol. 24, no. 6, pp. 456–462, 2000.

J. T. Chen and H.-K. Hong, “Review of dual boundary element methods with emphasis on hypersingular integrals and divergent series,” Applied Mechanics Reviews, vol. 52, no. 1, pp. 17–33, 1999.

Ch. Zhang, M. Cui, J. Wang, X. W. Gao, J. Sladek, and V. Sladek, “3D crack analysis in functionally graded materials,” Engineering Fracture Mechanics, vol. 78, no. 3, pp. 585–604, 2011.

Y. Ochiai, “Three-dimensional thermo-elastoplastic analysis by tripe-reciprocity boundary element method,” Engineering Analysis with Boundary Elements, vol. 35, pp. 478–488, 2011.

J. Bremer and Z. Gimbutas, “On the numerical evaluation of the singular integrals of scattering theory,” Journal of Computational Physics, vol. 251, pp. 327–343, 2013.

X.-W. Gao, W.-Z. Feng, K. Yang, and M. Cui, “Projection plane method for evaluation of arbitrary high order singular boundary integrals,” Engineering Analysis with Boundary Elements, vol. 50, pp. 265–274, 2015.

J. L. Tsalamengas, “Quadrature rules for weakly singular, strongly singular, and hypersingular integrals in boundary integral equation methods,” Journal of Computational Physics, vol. 303, pp. 498–513, 2015.

V. Sladek and J. Sladek, Singular Integrals in Boundary Element Methods, WIT Press, Southampton, UK, 1998.

M. Guiggiani, G. Krishnasamy, T. J. Rudolphi, and F. J. Rizzo, “A general algorithm for the numerical solution of hypersingular boundary integral equations,” Journal of Applied Mechanics, vol. 59, no. 3, pp. 604–614, 1992.

P. Bettini, M. F. Palumbo, and R. Specogna, “A boundary element method for eddy-current problems in fusion devices,” Fusion Engineering and Design, vol. 96-97, pp. 620–623, 2015.

A. Frangi and M. Guiggiani, “Direct approach for boundary integral equations with high-order singularities,” International Journal for Numerical Methods in Engineering, vol. 49, no. 7, pp. 871–898, 2000.

S. Mukherjee and Y. X. Mukherjee, “The hypersingular boundary contour method for three-dimensional linear elasticity,” Journal of Applied Mechanics, vol. 65, no. 2, pp. 300–309, 1998.

W. Chen, J.-Y. Zhang, and Z.-J. Fu, “Singular boundary method for modified Helmholtz equations,” Engineering Analysis with Boundary Elements, vol. 44, pp. 112–119, 2014.

M. H. Bazyar and A. Talebi, “Scaled boundary finite-element method for solving non-homogeneous anisotropic heat conduction problems,” Applied Mathematical Modelling, vol. 39, no. 23-24, pp. 7583–7599, 2015.

J. N. Lyness and C. B. Moler, “Numerical differentiation of analytic functions,” SIAM Journal on Numerical Analysis, vol. 4, pp. 202–210, 1967.

M. Cui, X. Gao, and H. Chen, “Inverse radiation analysis in an absorbing, emitting and non-gray participating medium,” International Journal of Thermal Sciences, vol. 50, no. 6, pp. 898–905, 2011.

M. Cui, X. W. Gao, and H. G. Chen, “A new approach for the equivalent gray radiative property of a non-gray medium using a modified zonal method and the complex-variable-differentiation method,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 112, no. 8, pp. 1336–1342, 2011.

M. Cui, Q. H. Zhu, and X. W. Gao, “A modified conjugate gradient method for transient nonlinear inverse heat conduction problems: a case study for identifying temperature-dependent thermal conductivities,” ASME Journal of Heat Transfer, vol. 136, no. 9, Article ID 091301, 2014.

M. Cui, W.-W. Duan, and X.-W. Gao, “A new inverse analysis method based on a relaxation factor optimization technique for solving transient nonlinear inverse heat conduction problems,” International Journal of Heat and Mass Transfer, vol. 90, pp. 491–498, 2015.