NONCONFORMING SPECTRAL ELEMENT METHOD FOR ELASTICITY INTERFACE PROBLEMS
NONCONFORMING SPECTRAL ELEMENT METHOD FOR ELASTICITY INTERFACE PROBLEMS
Journal of Applied Mathematics & Informatics. 2014. Sep, 32(5_6): 761-781
• Received : September 25, 2013
• Accepted : March 18, 2014
• Published : September 30, 2014
PDF
e-PUB
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
N. KISHORE KUMAR

Abstract
An exponentially accurate nonconforming spectral element method for elasticity systems with discontinuities in the coefficients and the flux across the interface is proposed in this paper. The method is leastsquares spectral element method. The jump in the flux across the interface is incorporated (in appropriate Sobolev norm) in the functional to be minimized. The interface is resolved exactly using blending elements. The solution is obtained by the preconditioned conjugate gradient method. The numerical solution for different examples with discontinuous coefficients and non-homogeneous jump in the flux across the interface are presented to show the efficiency of the proposed method. AMS Mathematics Subject Classification : 65N35, 65Y05, 74B05.
Keywords
1. Introduction
The elliptic interface problems arises in many engineering problems, for example, in heat conduction or elasticity problems where the domain of definition is composed of different materials. In this paper we study the nonconforming spectral element method for elasticity interface problems. These problems have wide applications in continuum mechanics, multi-phase elasticity problems, etc.
There exists several methods in the literature to solve elliptic interface problems (see [26] ). There are two types of finite element methods for elliptic interface problems: fitted and unfitted finite element methods. Fitted finite element method is a common approach where the mesh is fitted to the interface, so the interface conditions are satisfied in the weak formulation. The interface is approximated by the sides of isoparametric elements in the discretization. The accuracy of the method depends on the approximation quality of the interface. In this case, the method converges with optimal rates in h . If the mesh is not fitted to the interface, suboptimal convergence behavior will occur or may be the method does not converge at all. To avoid this difficulty, in [2] Babuska have formulated an equivalent minimization problem with all boundary and jump conditions incorporated in the cost functions.
Unfitted finite element methods are based on a mesh which is independent of interface. In [3] unfitted finite element method based on penalized problem as in [2] has been proposed. With an appropriate choice of penalty term the approximation converges to the solution at optimal rate (in h ) in H 1 norm. Only lower order finite elements have been studied in the literature (for more details see [27] ). In [27] , a conforming higher order finite element method has been analyzed for elliptic interface problems. One can also look for different formats of finite element methods for interface problems in [24 , 33 , 35 , 41] . Many nonconforming approaches are available in the literature like discontinuous Galerkin methods, Mortar finite element methods, etc (see [26] ). In [8] a priori and a posteriori error estimates have been derived for discontinuous Galerkin method. Quasi optimal a priori estimates for interface problems even with lower smoothness conditions on the solution were derived.
This problem is also studied in the framework of least-squares finite element method [4 , 5 , 6 , 11] . In these formulations, the given differential equation is converted into first order partial differential system and a suitable least-squares formulation is applied. Optimal convergence rates in h have been shown. The first order system least-squares method (FOSLS) for linear elasticity problems has been proposed in [9 , 10] . Least-squares spectral element method has been proposed in [19 , 20] . In [31 , 32] iterative substructuring methods for spectral element discretizations of elliptic systems have been proposed. The method provides an efficient preconditioner with an optimal condition number. The extended finite element method (XFEM) or generalized finite element method (GFEM) is a useful method for approximating the solutions with singularities and solutions of interface problems. This method extends the FEM approach by enriching the solution space. The approximation consists of standard finite element approximation and the enrichment through the partition of unity concepts [17 , 18] .
Immersed interface method has been widely studied for elliptic interface problems [28] . Finite difference based explicit jump immersed interface method for elasticity systems was described in [25] . An immersed finite element method for elasticity equations with interfaces has been studied in [29 , 39 , 40] . In this method the mesh is independent of the interface and basis functions are chosen such that they satisfies the interface conditions. Optimal convergence rates in h have been derived. In [30] , linear and bilinear immersed finite elements for planar elasticity interface problems have been discussed. 2 D linear, bilinear immersed finite elements which satisfy the interface jump conditions were used. Optimal convergence rates in h were shown in L 2 and semi H 1 norms. Details and complete citation list on immersed finite element methods can be found in [38] . In [12] an adaptive immersed interface finite element method for elasticity interface problem was presented.
In [22] Nitsche’s method has been described. In [23] a finite element solution of elliptic interface problem using an approach due to Nitsche has been proposed. The method allows for discontinuities, internal to the elements, in the approximation across the interface and it was shown to be second order accurate (in h ) in L 2 norm. In this method the interface conditions are satisfied weakly by means of variant of Nitsche’s method. A hp Nitsche’s method for interface problems with nonconforming unstructured finite element meshes have been proposed and error estimates with optimal bound in h and suboptimal bound in p by degree p 1/2 were obtained in [13] .
In this paper, we propose a least-squares spectral element method for elasticity interface problems based on the method proposed in [26] . In the least-squares formulation of the method, a solution is sought which minimizes the sum of the squares of a squared norms of the residuals in the partial differential equation and the sum of the residuals in the boundary conditions in fractional Sobolev norms and the sum of the jumps in the displacement and the flux across the interface in appropriate fractional Sobolev norms and enforce the continuity along the inter element boundaries by adding a term which measures the sum of the squares of the jump in the function and its derivatives in fractional Sobolev norms. The proposed numerical formulation is based on the regularity estimate for the interface problems stated in [7] and the stability estimate proved in [26] .
This method is nonconforming (in terms of approximation). This formulation is different from the standard techniques in LSFEM used to convert the second order elliptic equations into first order system. The interface is resolved completely using blending elements [21] . Higher order spectral elements are used to approximate the solution. The spectral elements are the sum of tensor products of the polynomials of degree W in each variable. The solution is obtained using preconditioned conjugate gradient method (PCGM) without storing the stiffness matrix and load vector. Even though we do not store the matrix, the added advantage of the proposed method is the resulting stiffness matrix is symmetric and positive definite. The integrals involved in the residual computations are obtained efficiently and inexpensively [36] (a brief description is given in the Appendix).
The rest of the paper is organized as follows: In Section 2 the elasticity interface problem is defined. The discretization of the domain is given in Section 3 and the numerical scheme is described. Finally in Section 4 numerical results are presented for various examples.
2. Elasticity interface problem
In this section we state the elasticity interface problem on a domain Ω ⊆ ℝ 2 . First we define the function spaces which we need in the latter sections.
Denoting Hk (Ω), the usual Sobolev space of integer order k with the norm ∥.∥ k,Ω as given below,
PPT Slide
Lager Image
Further, let
PPT Slide
Lager Image
denote the fractional Sobolev norm of order s , where 0 < s < 1. Here J denotes an interval contained in ℝ.
We denote vectors and vector spaces by bold characters. For example, u = ( u 1 , u 2 ) T , H k (Ω) = Hk (Ω) × Hk (Ω), etc. The norms are given by
PPT Slide
Lager Image
for
PPT Slide
Lager Image
etc.
2.1. Linear elasticity system. Let x = ( x , y ) be a point in space, u = ( u 1 ( x , y ), u 2 ( x , y )) T be the displacement vector and ϵ = ( ϵij ) be the strain tensor. If u 1 , u 2 are the two planar displacement components, then the straindisplacement relation is given by
PPT Slide
Lager Image
The relation between stresses and strains (from the Hooke’s law) is given by,
PPT Slide
Lager Image
where λ and μ are the Lame coefficients, and
PPT Slide
Lager Image
Let σ = ( σij ) be the stress tensor, f ( x ) = ( f 1 , f 2 ) T be the applied body forces, then the stress tensor satisfies the following partial differential equations,
PPT Slide
Lager Image
From the above equations, we can re-write the above system as the system of plane elasticity equations of the following,
PPT Slide
Lager Image
The Lame coefficients λ and μ are given by
PPT Slide
Lager Image
where E is the Young’s modulus and ν is the Poisson’s ratio. Let us assume that the constants λ and μ have finite jumps across the interface; so does the flux σ n . Now the elasticity interface problem is defined below.
2.2. The interface problem. Let Ω and
PPT Slide
Lager Image
be open bounded domains with boundaries
PPT Slide
Lager Image
respectively. Assume that the boundary Γ 0 is sufficiently smooth. Further, let
PPT Slide
Lager Image
Let
PPT Slide
Lager Image
Now the elasticity interface problem can be written as follows:
PPT Slide
Lager Image
where L u = ( L 1 u , L 2 u ) T , f = ( f 1 , f 2 ) T , q = ( q 1 , q 2 ) T , g = ( g 1 , g 2 ) T are known vector functions. n is the unit outward normal to the interface Γ 0 . The jump [.] is defined as the difference of the limiting values from the outside of the interface to the inside. The coefficients λ and μ are piecewise constant, i . e .,
PPT Slide
Lager Image
3. Discretization and Numerical Scheme
Considered the circular domain Ω 1 such that
PPT Slide
Lager Image
where Ω is square whose boundary is
PPT Slide
Lager Image
as shown in Figure 1 , for brevity. Let
PPT Slide
Lager Image
and the interface is Γ 0 which is smooth as shown in Figure 1 . The results presented are applicable to arbitrary smooth interfaces also.
PPT Slide
Lager Image
The domain Ω and discretization.
Now the domain Ω 1 and Ω 2 are partitioned into finite number of quadrilateral subdomains (elements)
PPT Slide
Lager Image
such that the subdomain divisions match on the interface. The interface is completely resolved using blending elements [21] .
Each element is mapped to the master square S = (−1, 1) 2 . Define an analytic map
PPT Slide
Lager Image
from the master square
PPT Slide
Lager Image
by (see [1 , 21] )
PPT Slide
Lager Image
A brief description of the map is given in the Appendix A1. Here and in the rest of this section l = 1, ..., p for i = 1 and l = 1, ..., q for i = 2.
Define the spectral element functions 1
PPT Slide
Lager Image
as the tensor product of polynomials of degree W in each variable ξ and η as
PPT Slide
Lager Image
Then 2
PPT Slide
Lager Image
are given by
PPT Slide
Lager Image
Now
PPT Slide
Lager Image
Here
PPT Slide
Lager Image
is the Jacobian of the mapping
PPT Slide
Lager Image
from
PPT Slide
Lager Image
Define 3
PPT Slide
Lager Image
(the differential operator in the transformed coordinates in the domains Ω 1 and Ω 2 respectively). Then
PPT Slide
Lager Image
Define
PPT Slide
Lager Image
Let
PPT Slide
Lager Image
for i = 1, 2. Define 4
PPT Slide
Lager Image
Let γs be a side common to the two adjacent elements
PPT Slide
Lager Image
(as shown in Fig. 2 (a) for i = 2). Assume that γs is the image of η = −1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
and also the image of η = 1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
By chain rule
PPT Slide
Lager Image
PPT Slide
Lager Image
Elements with common edges
Then the jumps along the inter-element boundaries are defined as
PPT Slide
Lager Image
Here and in what follows, I is an interval (−1, 1).
As the division of the domain into subdomains match along the interface, we define the jump across the interface by taking it (a part of interface) as the common edge. Consider the elements
PPT Slide
Lager Image
(as shown in Fig. 2 (b)) which have the common edge γs ⊆ Γ 0 . Let γs be the image of ξ = 1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
and also the image of ξ = −1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
Define
PPT Slide
Lager Image
where
PPT Slide
Lager Image
are the tangential derivatives of
PPT Slide
Lager Image
respectively.
Now along the boundary
PPT Slide
Lager Image
(for some j ) be the image of ξ = 1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
Then
PPT Slide
Lager Image
As defined earlier u 1 = u | Ω1 and u 2 = u | Ω2 , so the boundary condition u = g on Γ in the discrete form will be u 2 = g on
PPT Slide
Lager Image
be the image of the mapping
PPT Slide
Lager Image
corresponding to the side ξ = 1 and
PPT Slide
Lager Image
where −1 ≤ η ≤ 1.
On the interface Γ 0 we have [ u ] = 0 and [ σ n ] = q . Let γs ⊆ Γ 0 be the image of ξ = 1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
and also the image of ξ = −1 under the mapping
PPT Slide
Lager Image
which maps
PPT Slide
Lager Image
Let
PPT Slide
Lager Image
Let
PPT Slide
Lager Image
the space of spectral element functions. Define the functional
PPT Slide
Lager Image
The approximate solution is chosen as the unique
PPT Slide
Lager Image
which minimizes the functional
PPT Slide
Lager Image
over all
PPT Slide
Lager Image
The minimization problem leads to a symmetric and positive definite lineear system A Z = b . Where Z be a vector assembled from the values of
PPT Slide
Lager Image
at the Gauss-Lobatto-Legendre points arranged in lexicographic order for 1 ≤ k p , 1 ≤ l q . The solution is obtained by preconditioned conjugate gradient method. The action of a matrix on a vector in each iteration is obtained efficiently and inexpensively without storing the matrix A (since PCGM requires the action of a matrix on a vector). The details are shown in Appendix A2.
We used a preconditioner which was proposed in [15] . The preconditioner 5 is block diagonal matrix, where each diagonal block corresponds to the H 2 norm of the spectral element function representation of each component of the vector on a particular element which is mapped onto the master square S . The obtained solution of the preconditioned system is nonconforming. A set of corrections are made to the solution so that the corrected solution is conforming (see A3).
Let u be the exact solution and z be the approximate solution which is conforming. Let e = u z . Then for W large enough we have the following error estimate H 1 norm (since the jump in displacement across the interface is zero, u H 1 (Ω))
PPT Slide
Lager Image
holds, where C and b are constants. Proof is very similar to the one proven in [26] .
The displacement components in the transformed coordinates, where i gives the domain status (Ω1 or Ω2) and l gives the number of the element in that domain.
is vector u in element
4. Numerical Results
To prove the effectiveness of the method we present the numerical results for the problem defined in Section 3. The relative error ∥ e ER is defined as
PPT Slide
Lager Image
In all the examples, degree of the approximation polynomial is denoted by 6 W, ’DOF’ means the number of degrees of freedom and ’Iters’ means the total number of iterations required to compute the solution using PCGM. We have used the relative residual norm as a stopping criteria in PCGM. That is, the iteration process is stopped when the relative residual norm
PPT Slide
Lager Image
( r i is the residual in ith iteration, ∥.∥ 2 is vector norm) is less than the tolerance ϵ .
Example 1. Interface problem with homogeneous jump conditions: Consider the linear elasticity interface problem (5) (plane strain, see (4)) stated in Section 2 on a square domain [−0.75,0.75] 2 with a circle centered at the origin of radius s as the interface as shown in Fig. 1 . The coefficients λ, μ are given by
PPT Slide
Lager Image
Chosen the data such that the given interface problem has the exact solution u = ( u 1 , u 2 )
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Here we choose the radius of the circle
PPT Slide
Lager Image
Note that the solution u satisfies homogeneous jump conditions across the interface.
This problem have been studied in [39] . We discretized the domain into 9 quadrilateral elements as shown in Fig. 1 . The conforming numerical solution has been obtained for various values of b (see (8)) for different degree of the approximating polynomial W . The relative error ∥ e ER in percent, the iterations are tabulated for different values of b (see (8)) in Table 1 and Table 2 .
The relative error in percent and iterations for differentW
PPT Slide
Lager Image
The relative error in percent and iterations for different W
The relative error in percent and iterations for differentW
PPT Slide
Lager Image
The relative error in percent and iterations for different W
For smaller values of b (see (8)) the iteration count is less but the iteration count is large when b b is large. More efficient preconditioner is under investigation. The log of relative error against the degree of the approximating polynomial W is drawn in Fig. 3 for b = 50, 100. The relation is almost linear. This shows the exponential accuracy of the method. The error decays exponentially for all values of b .
PPT Slide
Lager Image
Log of the relative error against W for b = 50, 100
Example 2. Interface problem with nonhomogeneous jump in the flux: Consider the linear elasticity problem as defined in (5) (plane strain problem) on the same domain as considered in Example 1 with a circle centered at the origin of radius s as the interface. The coefficients λ, μ are given by
PPT Slide
Lager Image
Chosen the data such that the given interface problem has the exact solution u = ( u 1 , u 2 )
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Here we choose the radius of the circle
PPT Slide
Lager Image
For any b , the solution u is continuous across the interface and [ σ n ] = q = ( q 1 , q 2 ),
PPT Slide
Lager Image
The domain is discretized as in the above example and the conforming solution is obtained for b =10, 0.1. The relative error ∥ e ER in percent, the iterations are tabulated in Table 3 . The log of relative error against W is drawn in Fig. 4 for b = 10, 0.1. Here one can see that the iteration count is large compared to the count in previous example (look at Table 1 for b = 10).
The relative error in percent and iterations for differentW
PPT Slide
Lager Image
The relative error in percent and iterations for different W
PPT Slide
Lager Image
Log of the relative error against W for b = 0.1, 10
For the simplicity of the programming W is assumed to be uniform in all elements.
5. Conclusions
The proposed method is nonconforming and exponentially accurate. The interface is resolved exactly using blending elements. A small data has to be interchanged in between the elements for each iteration of the PCGM and the residuals in the normal equations can be obtained efficiently and inexpensively. The proposed method is efficient even when the jump in the coefficient is large. The numerical results shows that large differences in the coefficients leads to increase in the number of iterations of the PCGM. A more efficient preconditioner is under investigation. This method is also efficient on parallel computers. The method is applicable to arbitrary smooth interfaces too and the method can be extended to the singular case which is ongoing work.
BIO
N. Kishore Kumar received PhD from Indian Institute of Technology, Kanpur, India. Since December 2013 he has been at BITS-Pilani Hyderabad Campus, India. His research interests include numerical solutions to PDE and Tensor decompositions.
e-mail:thimmaki@gmail.com
References
Babuska I. (1998) The hp version of the finite element method for domains with curved boundaries SIAM J. Numer. Anal. 25 (4) 837 - 861    DOI : 10.1137/0725048
Babuska I. (1970) The finite element method for elliptic equations with discontinuous coefficients Computing 5 207 - 213    DOI : 10.1007/BF02248021
Barret J.W. , Elliott C.M. (1987) Fitted and unfitted finite element methods for elliptic equations with smooth interfaces IMA Journal of Num. Anal. 7 283 - 300    DOI : 10.1093/imanum/7.3.283
Berndt M. , Manteuffel T.A. , McCormick S.F. , Starke G. (1995) Analysis of first-order system least-squares (FOSLS) for elliptic problems with discontinuous coefficients Part I, SIAM J. Numer. Anal. 43 (1) 386 - 408    DOI : 10.1137/S0036142903427688
Berndt M. , Manteuffel T.A. , McCormick S.F. (1995) Analysis of first-order system leastsquares (FOSLS) for elliptic problems with discontinuous coefficients Part II, SIAM J. Numer. Anal. 43 (1) 409 - 436    DOI : 10.1137/S003614290342769X
Bochev P.V. , Gunzburger M.D. 2009 Least-squares finite element methods Springer
Bramble J.H. , King J.T. (1996) finite element method for interface problems in domains with smooth boundaries and interfaces Adv. Com. Math. 6 109 - 138    DOI : 10.1007/BF02127700
Cai Z. , Ye X. , Shun Z. (2011) Discontinuous Galerkin finite element methods for interface problems; a priori and a posteriori error estimations SIAM Jour. of Num. Anal. 49 (5) 1761 - 1787    DOI : 10.1137/100805133
Cai Z. , Manteuffel T.A. , McCormick S.F. , Parter S. (1998) First-order system least squares for the pure traction problem in planar linear elasticity SIAM J. Numer. Anal. 35 320 - 335    DOI : 10.1137/S0036142995294930
Cai Z. , Lee C.O. , Manteuffel T.A. , McCormick S.F. (2000) First-order system least squares for planar linear elasticity: numerical results SIAM J. Sci. Comp. 21 1706 - 1727    DOI : 10.1137/S1064827598338640
Cao Y. , Gunzburger M.D. (1998) Least-square finite element approximations to solutions of interface problems SIAM. J. Numer. Anal. 35 (1) 393 - 405    DOI : 10.1137/S0036142996303249
Chang Y. (2012) The adaptive immersed interface finite element method for elasticity interface problem Journal of Computational Mathematics 30 (6) 629 - 642    DOI : 10.4208/jcm.1203-m3869
Chernov A. , Hansbo P. (2011) A hp Nitsche’s method for interface problems with nonconforming unstructured finite element meshes Lec. Notes in Comp. Sci. and Eng. 76 153 - 161
Dutt P.K. , Bedekar S. (2001) Spectral methods for hyperbolic initial value problems on parallel computers J. of Comp. and Appl. Math. 134 165 - 190    DOI : 10.1016/S0377-0427(00)00535-5
Dutt P. , Biswas P. , Naga Raju G. (2008) Preconditioners for spectral element methods for elliptic and parabolic problems J. Comput. Appl. Math. 215 (1) 152 - 166    DOI : 10.1016/j.cam.2007.03.030
Dutt P.K. , Singh A.K. (1994) The Galerkin-Collocation method for hyperbolic initial boundary value problems Jour. Comp. Physics. 112 211 - 225    DOI : 10.1006/jcph.1994.1093
Fries T.P. , Belytschko T. (2006) The intrinsic XFEM, a method for arbitrary discontinuous without additional unknowns Int. J. Num. Meth. Engg. www.xfem.rwth-aachen.de. 68 1358 - 1385    DOI : 10.1002/nme.1761
Fries T.P. , Belytschko T. (2010) The extended/generalized FEM; An overview of the method and its applications Int. J. Num. Meth. Engg. 84 253 - 304
Galvao A. , Gerritsma M.I. , Maercshalck De (2008) hp Adaptive least-squares spectral element method for hyperbolic partial differential equations J. Comp. Appl. Math. 215 (2) 409 - 418    DOI : 10.1016/j.cam.2006.03.063
Gerritsma M.I. , Proot M.J. (2002) Analysis of a discontinuous least-squares spectral element method J. Sci. Com. 17 (1-4) 297 - 306    DOI : 10.1023/A:1015173203136
Gordan W.J. , Hall C.A. (1973) Transfinite element methods: Blending-function interpolation over arbitrary curved element domains Numer. Math. 21 (2) 109 - 129    DOI : 10.1007/BF01436298
Hansbo P. (2006) Nitche’s method for interface problems in computational mechanics GAMMMitt 28 (2) 183 - 206    DOI : 10.1002/gamm.201490018
Hansbo A. , Hansbo P. An unfitted finite element method based on Nitsche’s method for elliptic interface problems Com. Meth. Appl. Mech. and Eng. 191 (47-48) 5537 - 5552    DOI : 10.1016/S0045-7825(02)00524-8
Hou S. , Li Z. , Wang L. , Wang W. (2012) A numerical method for solving elasticity equations with interfaces Comm. in Comp. Physics 12 (2) 595 - 612
Rutka V. , Wiegmann A. (2006) Explicit jump immersed interface method for virtual material design of the effective elastic moduli of composite materials Numer. Algo. 43 (4) 309 - 330    DOI : 10.1007/s11075-007-9063-9
Kishore Kumar N. , Naga Raju G. (2012) Nonconforming least-squares method for elliptic partial differential equations with smooth interfaces Jour. of Scientific Computing 53 (2) 295 - 319    DOI : 10.1007/s10915-011-9572-5
Li J. , Melenk J.M. , Wohlmuth B. , Zou J. (2010) Optimal a priori estimates for higher order finite element methods for elliptic interface problems Applied Num. Math. 60 19 - 37    DOI : 10.1016/j.apnum.2009.08.005
Li Z. , Ito K. 2006 The immersed interface method: Numerical solutions of PDEs involving interfaces and irregular domains Frontiers Appl. Math. SIAM Philadelphia 33
Li Z. , Yang X. 2005 An immersed finite element method for elasticity equations with interfaces Contemporary Mathematics 383 285 - 298
Lin T. , Zhang X. Linear and bilinear immersed finite elements for planar elasticity interface problems Jour. Com. App. Math. 236 (18) 4681 - 4699    DOI : 10.1016/j.cam.2012.03.012
Pavarino L. , Widlund O. Iterative substructuring methods for spectral element discretizations of elliptic systems. I: Compressible linear elasticity SIAM J. Num. Anal. 37 (2) 353 - 374    DOI : 10.1137/S003614299732824X
Pavarino L. , Widlund O. Iterative substructuring methods for spectral element discretizations of elliptic systems. II: Mixed methods for elasticity and Stokes flow SIAM J. Num. Anal. 37 (2) 375 - 402    DOI : 10.1137/S0036142998333092
Peichl G.H. , Touzani R. (2007) An optimal order finite element method for elliptic interface problems Proc. Appl. Math. Mech. 7 1025403 - 1025404    DOI : 10.1002/pamm.200701144
Schwab Ch. 1998 p and h - p Finite Element Method Clarendon Press Oxford
Sinha R.K. , Deka B. (2007) An unfitted finite element method for elliptic and parabolic interface problems IMA Jour. of Num. Ana. 27 529 - 549    DOI : 10.1093/imanum/drl029
Tomar S.K. (2006) h - p spectral element method for elliptic problems on non-smooth domains using parallel computers Computing 78 117 - 143    DOI : 10.1007/s00607-006-0176-0
Tomar S.K. , PhD thesis 2001 h - p spectral element for elliptic problems on non-smooth domains using parallel computers IIT Kanpur India PhD thesis http://doc.utwente.nl/65818/ Reprint available as Tec. Rep. No. 1631, Department of applied Mathematics, University of Twente, The Netherlands
He Xiaoming , PhD thesis Bilinear immersed finite elements for interface problems Virginia polytechnic institute and state university PhD thesis
Xie H. , Li Z. , Qiao Z. A finite element method for elasticity interface problems with locally modified triangles Int. Jour. of Nume. Anal. and Mod. 8 (2) 189 - 200
Yang X. , Li B. , Li Z. 2003 The immersed interface method for elasticity problems with interfaces, Dynamics of continuous, Discrete and impulsive systems 10 (5) 783 - 808
Wu H. , Xiao Y. An unfitted hp-interface penalty finite element method for elliptic interface problems http://arxiv.org/pdf/1007.2893v1.pdf