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

Transcription

1 Mathematical Problems in Engineering Volume 06, Article ID , 8 pages 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 604, China Correspondence should be addressed to Kai Yang; kyang@dlut.edu.cn Received 7 December 05; Revised February 06; Accepted 7 February 06 Academic Editor: Ivano Benedetti Copyright 06 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. 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 twodimensional 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.. Introduction Boundary element method (BEM) [] is a very promising approach for solving various engineering problems, such as heat transfer problems [ 5] in energy science and engineering, in which accurate evaluation of boundary integrals [6 0] is required. Weak, strong, or hypersingularities are involved in these boundary integrals, if the source point is located on the element under integration []. These issues have been focused on by researchers, and many techniques, such as Gauss logarithmic quadrature rule for weak singularities [], rigid body displacement algorithm for strong singularities [], and regularization for hypersingularities [3], have been proposed to remove the singularities. Each method has its advantage and disadvantage [ 9]. A comprehensive review has been given in [] on the treatment of singularities. In general, the efficient methods for singular curved line and surfaceintegralsarethosebasedontheoperationoverthe intrinsic coordinate system [4 6]. Recently, Gao et al. [] proposed another efficient direct method for numerically evaluating all kinds of singular curved boundary integrals for two-dimensional (D) 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 [] 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 []. 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 [0] is introduced to determine the third-order derivatives. For D boundary integral, the third-order derivatives are analytically derived, which are used to verify the CVDM. The paper is organized as follows. In Section, review and modification of the direct method in [] are provided. In Section 3, numerical examples are given to show the

2 Mathematical Problems in Engineering effectiveness and the accuracy of the present method. Finally, conclusions are drawn from computational results and discussions. x x p Q ρ P. Development of the Direct Method for Evaluating Arbitrary High Order Singular Curved Boundary Integrals ξ Q r P.. Review of Singular Integrals []. Singular integrals over an element can be classified into the following form: I e (x p )= Γ e f(x p, x) r β (x p, x) dγ (x) () in which x p and x represent the source and field points, respectively; Γ e is the boundary element under integration, which is a curved line or a curved surface for D or 3D problems, respectively; and r is the distance between the source and field points. The term f(x p, x) in () is a regular function and β is the order of the singularity. In D problems, the following type of singular line integrals is also frequently encountered: J e (x p )= Γ e f(x p, x) ln [r (x p, x)] dγ (x). () r β (x p, x) x x Figure : The curved line boundary and the projection line. In D problems, the boundary element is a curved line as shown in Figure. The local orthogonal coordinate system is denoted by (x,x ) and the projection line is along the axis x. There is only one independent variable for a line, and all other geometry quantities can be expressed by x.similarto the treatment in [], the coordinate x at the field point over thecurvedlinecanbeexpressedbyx using Taylor series. Different from the treatment in [], the third-order derivative is considered in the present work: It is assumed that the integrals in () and () always exist: that is, each integration has a finite value... Review of Expansions of Geometry Quantities on the Projection Line/Plane []. In [], a projection line for D 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: x i =L ij (x j x o j ), (3) x i =x o i +L ji x j, (4) L ij is the direction cosine of the local coordinate axes with respect to the global one. The repeated subscripts represent summation; x o i 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. x = x x x + x x x x x 3 x 3 =a (x )x +a (x )x +a 3 (x )x 3, a (x )= x x = x x, a (x )= x x = x a 3 (x )= 3 x 6 x 3 = 6 3 x 3 + x 3 ξ. x 3 = 6 x x x x x x ( x ) x +3 x + x x ξ x ξ, x (5) (6).3. High Order Expansion of Geometry Quantities in terms of Variables on the Projection Line for D Problems. For D and 3D problems, derivations for the modification are different, which are separately described in detail. It can be seen that the expressions are more complex. For D problem, the corresponding derivatives are analytically derived, and the detailed derivation process is omitted here.

3 Mathematical Problems in Engineering 3 In (5) and (6), x Nnode = N α α= x α, (7a) p(x )=a x p +a (x )x p x p (4) +a 3 (x )x p x p x p =x p, x N node = α= 3 x N node 3 = α= N α x α, (7b) 3 N n 3 x α, (7c) x = m, (8a) ξ x = m 4 ( x x + x x ), (8b) 3 ξ x 3 = x ( ξ x )= m B 4C. (8c) m 7 In (7a) (8c), m = ( x ) +( x ), (9) B=( x ) +( x ) + x 3 x 3 + x 3 x 3, C=( x ) ( x ) + x x x x +( x ) ( x ). Define a local distance: (0a) (0b) ρ= x x p, () x p is the local coordinate of point P which is the projection point of the source point P and x is the local coordinate of point Q which is the projection point of the field point Q. The derivative of ρ with respect to local coordinate x is In (a), ρ, =or. Then, we can obtain ρ, = ρ x = x x p. (a) ρ x =x p +ρ,ρ. (b) By substituting (a) into (5), we can obtain x =p(x )+q(x,ρ)ρ, (3) In (5), q (x ) q(x,ρ)=q (x )+q (x )ρ+q 3 (x )ρ. (5) =[a (x )+a (x )x p +3a 3 (x )x p ]ρ,, (6a) q (x )=a (x )ρ, +3a 3 (x )x p (6b) =a (x )+3a 3 (x )x p, q 3 (x,ρ)=a 3 (x )ρ,. (6c) Similarly, we can express the intrinsic coordinate in terms of the local distance as ξ=p(ξ) +q(ξ, ρ) ρ, (7) p(ξ) and q(ξ, ρ) can be determined by replacing x with ξ in (3) (6c). Substituting (b) and (3) into (4), we can obtain x i =x o i +L ji x j =xo i +L i x +L ix =x o i +L i (x p +ρ,ρ) + L i [x p +q(x,ρ)ρ] =x p i +L i ρ, ρ+l i q(x,ρ)ρ. Then, substituting (5) into (8), we can obtain (8) x i =x p i +(b i +c i ρ+d i ρ )ρ, (9) b i =L i ρ, +L i q (y ), c i =L i q (y ), d i =L i q 3 (y ). (0a) (0b) (0c) Referring to Figure, the global distance r canbeeasily expressed as r= (x i x p i )(x i x p i ) = g (ρ) ρ, () g(ρ)= (b i +c i ρ+d i ρ )(b i +c i ρ+d i ρ ) = G 0 +G ρ+g ρ +G 3 ρ 3 +G 4 ρ 4. () In (), G 0 =b i b i, G =b i c i, G =c i c i +b i d i, G 3 =c i d i,and G 4 =d i d i. Itcanbeeasilyderivedthat r,i = r = x i x p i x i r = b i +c i ρ+d i ρ. (3) g(ρ)

4 4 Mathematical Problems in Engineering.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) [0] is introduced to determine the high order derivatives. For D analysis, the third-order derivatives are determined explicitly by (7c), (8c), (9), (0a), and (0b), which are also calculated by the CVDM. In the CVDM, the variable x of a real function f(x) is replaced by a complex one x+ih, with the imaginary part h being very small (usually h =0 0 ). The function f(x + ih) can be expanded into Taylor s series as x C C B B x x 3 o n 0 A A D D f (x+ih) =f(x) +ihf (x) h f (x) +o(h 3 ). (4) Since h issmallenough,(5)canbeobtained: f (x) = Im (f (x+ih)), (5) h Im denotes the imaginary part. The CVDM is a very promising method, since the derivatives only require direct evaluations []. The CVDM has been validated to be an efficient and accurate approach for determination of derivatives [ 4]. For example, considering 3 x / x 3 = ( / x )( x / x ), we can add the very small imaginary part h to x ; then the real value of 3 x / x 3 yields Im( x / x )/h. In addition, the value of x / x can be easily obtained using the analytical expressions in the previous work []. The CVDM and the derived expressions in (7c) (8c) can be validatedbyeachother..5. Expansion of Geometry Quantities in terms of Variables ontheprojectionplanefor3dproblems. In 3D problems, the boundary element is a curved surface marked as ABCD in Figure. A local orthogonal coordinate system (x,x,x 3 )is established on the projection plane with its origin being at the point (ξ =0, ξ =0), in which axes x and x are located within the plane. The axis x is along ξ direction and axis x 3 is along the outward normal direction to the element. By (4), wecanprojecttheoriginalcurvedsurfaceelementontothe projection plane, to form a flat projection element marked as A B C D in Figure. In a plane, there are only two independent variables. We choose x and x as the independent variables, and x 3 over the curved surface can be expressed in terms of x and x. Compared with the treatment in [], x 3 is expanded as Taylor series by x and x, by considering the third-order derivatives: x 3 =a I (x 3 )x I +a IJ (x 3 )x I x J +a IJK (x 3 )x I x J x K, (6) Figure : The curved surface boundary and the projection plane. I, J,andK are integers varying from to, and a I (x 3 )= x 3 x I a IJ (x 3 )= x 3 xi x J = x 3 K, (7a) K x I = x 3 K L K L xi xj + x 3 K a IJK (x 3 )= 3 x 3 6 xi = a IJ (x 3 ) x J x K 3 xk = a IJ (x 3 ) L. 3 L x K ξ K xi, x J (7b) (7c) The treatment is the same as that in [], to determine a I (x 3 ) in (7a) and a IJ(x 3 ) in (7b). For a IJK(x 3 ) in (7c), we use the CVDM, to avoid rather complicated analytical derivations. The local distance ρ projected from the global distance r ontotheprojectionplaneanditsderivativeρ,i are introduced as follows: as ρ= (x x p ) +(x x p ), ρ,i = ρ x I = x I x p I. ρ (8) Thelocalcoordinatesatthefieldpointcanbeexpressed x I =x p I +ρ,i ρ. (9) Similar to the D analysis, substituting (9) into (6) yields x 3 =p(x 3 )+q(x 3,ρ)ρ, (30)

5 Mathematical Problems in Engineering 5 In (3), p(x 3 )=a I (x 3 )x p I +a IJ (x 3 )x p I x p J +a IJK (x 3 )x p I x p J x p K =x p 3, (3) q(x 3,ρ)=q (x 3 )+q (x 3 )ρ+q 3 (x 3 )ρ. (3) q (x 3 )=a I (x 3 )ρ,i +a IJ (x 3 )x p I ρ,j +a IJ (x 3 )x p J ρ,i +a IJK (x 3 )x p I x p K ρ,j +a IJK (x 3 )x p J x p K ρ,i (33a) In (39), G 0 =b i b i, G =b i c i, G =c i c i +b i d i, G 3 =c i d i,and G 4 =d i d i. Itcanbeeasilyderivedthatr,i has the same expression as that in D analysis. Considering that I=J=K=and ρ, =or for D 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 () and (), with the above modifications. The D problem is taken as an example to describe the procedure to determine the integral in (). The differential relationship between the localdistanceandtherealcurvedlinecanbewrittenas dρ = dγ cos φ=dγn 0 i n i, (40) +a IJK (x 3 )x p I x p J ρ,k, q (x 3 )=a IJ (x 3 )ρ,iρ,j +a IJK (x 3 )x p K ρ,iρ,j +a IJK (x 3 )x p I ρ,jρ,k +a IJK (x 3 )x p J ρ,iρ,k, q 3 (x 3 )=a IJK (x 3 )ρ,iρ,j ρ,k. (33b) (33c) n 0 i and n i are the outward normals to the tangential linesofpassingtheoriginofthelocalcoordinatesystemand the field point, respectively, and φ is the angle between the two outward normals. Substituting (40) and (38) into () leads to I e (x p )=lim ρ E ε 0 ε F(ρ) dρ, (4) ρβ Similar to the D analysis, intrinsic coordinates can be expressed as ξ K =p(ξ K )+q(ξ K,ρ)ρ, (34) p(ξ K ) and q(ξ K,ρ)can be determined by replacing x 3 with ξ K in (3) (33c). Substituting (9) and (30) into (4) yields x i =x o i +L Ji x J +L 3ix 3 =x o i +L Ji (x p J =x p i +L Ji ρ,j ρ+l 3i q(x 3,ρ)ρ. +ρ,jρ) + L 3i [x p 3 +q(x 3,ρ)ρ] Substituting (3) into (35), we can obtain (35) x i =x p i +(b i +c i ρ+d i ρ )ρ, (36) b i =L Ji ρ,j +L 3i q (x 3 ), c i =L 3i q (x 3 ), d i =L 3i q 3 (x 3 ). The expression of r is as follows: (37a) (37b) (37c) r= (x i x p i )(x i x p i ) = g (ρ) ρ, (38) g(ρ)= (b i +c i ρ+d i ρ )(b i +c i ρ+d i ρ ) = G 0 +G ρ+g ρ +G 3 ρ 3 +G 4 ρ 4. (39) ρ E isthedistancefromtheprojectedsourcepointp to one of the ends of the projected line element and F(ρ) is the regular part: that is, F(ρ)= f(xp, x (ρ)) g β (ρ) n 0 i n i (ρ). (4) In order to integrate (4), the nonsingular part F is expanded as a power series in ρ,suchthat F(ρ)= N n=0 B (n) ρ n (43) in which N is the order of the power series, usually taking a value between and 7 depending on the size of ρ E,andB (n) are constants which are determined by collocating N+points over the integration region (0, ρ E ).Inthepaper,N+equally spaced points are used: that is, (0, ρ,...,ρ N ).Thecoefficient for the first point (n = 0) is B (0) = F(0) and other coefficients can be solved using the following equation set: [R] B = Y, (44) [R] is a square matrix with the order of N: ρ ρ N ρ ρ N [R] = [... (45). ] [ ρ N ρ N N ]

6 6 Mathematical Problems in Engineering B and Yare vectors as follows: B = B B. B N, [F(ρ ) B (0) ] ρ [F(ρ ) B (0) ] Y = ρ... [F(ρ N ) B (0) ] ρ N (46) Solving (44) for coefficient vector B and then substituting (43) into (4) yield E n I e (x p )= N n=0 ρ E ε 0 ε B (n) lim ρ n β dρ = N n=0 B (n) E n, (47) n β+ ( ρ lim β n ρ = E ε 0ρ ) (n =β ) β n ε ln ρ E lim ln ρ ε (n = β ). ρε 0 Then, singularities can be analytically removed. 3. Numerical Examples (48) BasedontheFortrancodeforsingularintegralsin[],anew Fortrancodehasbeendevelopedwiththeabovemodification 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 bezero,thefunctionsarethesamewiththoseofthemethod in []. Besides, the new code has its unique characteristics. For D 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 complicatedtobeobtained. In this section, the first example is given to show the advantageofcvdmoveradifferencemethod.then,two singular integrals are calculated, to show the effectiveness and the accuracy of the present method. 3.. The Advantage of the CVDM. A singular function is selected: f (x, x p ) = (x x p ). (49) Table : The first-order derivatives with different h for f(x, x p )= /(x x p ). h CVDM Difference method.00e E E E E E E E E E E E E E E E E E E E E E E E E E E + 00 Then, the exact derivative with respect to the variable x p can be obtained: f x p = (x x p ) 3. (50) For comparison, the central difference formulation is adopted: that is, f x p = f(x,xp +h) f(x,x p h). (5) h When x =.00, x p =, the exact first-order derivative is 0 9. The dependence of the derivatives on h for the CVDM and the difference method is examined, and the results are listed in Table. It can be seen that the difference method can give satisfactory results only for h between 0 5 and 0 0.For smaller h, its result becomes unacceptable. However, the result of the CVDM tends to be independent of h as the step size becomes sufficiently small. This clearly verifies the advantage of the CVDM. 3.. A Weakly Singular Line Integral over a Cubed Curve. In this example, a singular line integral for a D BEM problem is chosen: J= ln (r) dγ. (5) Γ The integration line is a cubed curve as shown in Figure 3. The four typical points are (, ), (, ), 3 ( /3, /7), and 4(/3, /7). The computational point (sourcepoint)isat.thecurveistakenasoneboundary element. It should be emphasized here that the direct method in [] has been invalid, which cannot be used for the singular line integration in (5) 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. 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.

7 Mathematical Problems in Engineering 7 Table : Results for two points a and b (β =3). Point a Point b Method in [4] Present method z Figure 3: The cubed curve for integration. a b The boundary integral is as follows: I(x p )= Γ 4πr [3r r β,3 n n 3]dΓ. (54) Two source points, a and b, are located at positions with intrinsic coordinates of (0, 0) and (0.66, 0). Guiggiani et al. [4] computed this integral for the case of β=3. Table lists the results using both methods. For the present method, the Gauss order is 4. From Table, 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. Figure 4: The cylindrical surface. Using the derived analytical expressions in this paper, the third-order derivatives can be obtained: y 3 x x 3 = , 3 ξ x 3 = x (53) TheresultsoftheCVDMareexactlythesamewiththat 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 , with the Gauss order being 6. This integral is also evaluated using the Gauss quadrature rule for logarithmic function []. The result is J = It can be seen that the present method is effective and with high accuracy, for higher order curved boundary 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 and the rotational angle is 90. The geometrical dimensions are shown in Figure Conclusion and Future Work Inthepresentwork,highorderprojectionplanemethod forevaluationofsupersingularcurvedboundaryintegrals in BEM has been developed, by considering the thirdorder 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 []. 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 complexvariable-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 ( and 57075) and

8 8 Mathematical Problems in Engineering the Fundamental Research Funds for the Central Universities (DUT4LK03). 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. 98, pp. 9 57, 06. []H.Qi,L.M.Ruan,H.C.Zhang,Y.M.Wang,andH.P.Tan, Inverse radiation analysis of a one-dimensional participating slab by stochastic particle swarm optimizer algorithm, International Thermal Sciences,vol.46,no.7,pp , 007. [3] S.Zhang,Y.Zhu,B.Qiao,andZ.Zhang, Systemmodelofheat andmasstransferprocessformobilesolventvaporphasedrying equipment, Mathematical Problems in Engineering, vol.04, Article ID 6776, pages, 04. [4] 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.8,pp.6 68,05. [5] H.Qi,Y.B.Qiao,S.C.Sun,Y.Yao,andL.Ruan, Imagereconstruction of two-dimensional highly scattering inhomogeneous medium using MAP-based estimation, Mathematical Problems in Engineering, vol. 05, Article ID 435, 9 pages, 05. [6] 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.4,no.6,pp ,000. [7] 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. 5, no., pp. 7 33, 999. [8]Ch.Zhang,M.Cui,J.Wang,X.W.Gao,J.Sladek,andV. Sladek, 3D crack analysis in functionally graded materials, Engineering Fracture Mechanics,vol.78,no.3,pp ,0. [9] Y. Ochiai, Three-dimensional thermo-elastoplastic analysis by tripe-reciprocity boundary element method, Engineering Analysis with Boundary Elements,vol.35,pp ,0. [0] J. Bremer and Z. Gimbutas, On the numerical evaluation of the singular integrals of scattering theory, Computational Physics,vol.5,pp ,03. [] 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 , 05. [] J. L. Tsalamengas, Quadrature rules for weakly singular, strongly singular, and hypersingular integrals in boundary integral equation methods, Computational Physics, vol. 303, pp , 05. [3] V.SladekandJ.Sladek,Singular Integrals in Boundary Element Methods, WIT Press, Southampton, UK, 998. [4] M. Guiggiani, G. Krishnasamy, T. J. Rudolphi, and F. J. Rizzo, A general algorithm for the numerical solution of hypersingular boundary integral equations, Applied Mechanics,vol. 59,no.3,pp ,99. [5] P. Bettini, M. F. Palumbo, and R. Specogna, A boundary element method for eddy-current problems in fusion devices, Fusion Engineering and Design, vol , pp , 05. [6] 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 , 000. [7] S. Mukherjee and Y. X. Mukherjee, The hypersingular boundary contour method for three-dimensional linear elasticity, Applied Mechanics, vol. 65, no., pp , 998. [8] W. Chen, J.-Y. Zhang, and Z.-J. Fu, Singular boundary method for modified Helmholtz equations, Engineering Analysis with Boundary Elements,vol.44,pp. 9,04. [9] 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. 3-4, pp , 05. [0] J. N. Lyness and C. B. Moler, Numerical differentiation of analytic functions, SIAMJournalonNumericalAnalysis,vol. 4, pp. 0 0, 967. [] M. Cui, X. Gao, and H. Chen, Inverse radiation analysis in an absorbing, emitting and non-gray participating medium, International Thermal Sciences,vol.50,no.6,pp , 0. [] 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-variabledifferentiation method, Quantitative Spectroscopy and Radiative Transfer,vol.,no.8,pp ,0. [3] 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 Heat Transfer,vol.36, no.9,articleid0930,04. [4] 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 Heat and Mass Transfer,vol.90,pp , 05.

9 Advances in Operations Research Advances in Decision Sciences Applied Mathematics Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Discrete Mathematics Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis International Stochastic Analysis Optimization