A new multiscale finite element method for high-contrast elliptic interface problems
HTML articles powered by AMS MathViewer
- by C.-C. Chu, I. G. Graham and T.-Y. Hou PDF
- Math. Comp. 79 (2010), 1915-1955 Request permission
Abstract:
We introduce a new multiscale finite element method which is able to accurately capture solutions of elliptic interface problems with high contrast coefficients by using only coarse quasiuniform meshes, and without resolving the interfaces. A typical application would be the modelling of flow in a porous medium containing a number of inclusions of low (or high) permeability embedded in a matrix of high (respectively low) permeability. Our method is $H^1$- conforming, with degrees of freedom at the nodes of a triangular mesh and requiring the solution of subgrid problems for the basis functions on elements which straddle the coefficient interface but which use standard linear approximation otherwise. A key point is the introduction of novel coefficient-dependent boundary conditions for the subgrid problems. Under moderate assumptions, we prove that our methods have (optimal) convergence rate of $O(h)$ in the energy norm and $O(h^2)$ in the $L_2$ norm where $h$ is the (coarse) mesh diameter and the hidden constants in these estimates are independent of the “contrast” (i.e. ratio of largest to smallest value) of the PDE coefficient. For standard elements the best estimate in the energy norm would be $\mathcal {O}(h^{1/2-\varepsilon })$ with a hidden constant which in general depends on the contrast. The new interior boundary conditions depend not only on the contrast of the coefficients, but also on the angles of intersection of the interface with the element edges.References
- Loyce Adams and Zhilin Li, The immersed interface/multigrid methods for interface problems, SIAM J. Sci. Comput. 24 (2002), no. 2, 463–479. MR 1951051, DOI 10.1137/S1064827501389849
- Mark Ainsworth, Robust a posteriori error estimation for nonconforming finite element approximation, SIAM J. Numer. Anal. 42 (2005), no. 6, 2320–2341. MR 2139395, DOI 10.1137/S0036142903425112
- Ivo Babuška, The finite element method for elliptic equations with discontinuous coefficients, Computing (Arch. Elektron. Rechnen) 5 (1970), 207–213. MR 277119, DOI 10.1007/bf02248021
- Ivo Babuška, Gabriel Caloz, and John E. Osborn, Special finite element methods for a class of second order elliptic problems with rough coefficients, SIAM J. Numer. Anal. 31 (1994), no. 4, 945–981. MR 1286212, DOI 10.1137/0731051
- I. Babuška and J. E. Osborn, Generalized finite element methods: their performance and their relation to mixed methods, SIAM J. Numer. Anal. 20 (1983), no. 3, 510–536. MR 701094, DOI 10.1137/0720034
- John W. Barrett and Charles M. Elliott, Fitted and unfitted finite-element methods for elliptic equations with smooth interfaces, IMA J. Numer. Anal. 7 (1987), no. 3, 283–300. MR 968524, DOI 10.1093/imanum/7.3.283
- C. Bernardi and R. Verfürth, Adaptive finite element methods for elliptic equations with non-smooth coefficients, Numer. Math. 85 (2000), no. 4, 579–608 (English, with English and French summaries). MR 1771781, DOI 10.1007/PL00005393
- I-Liang Chern and Yu-Chen Shu, A coupling interface method for elliptic interface problems, J. Comput. Phys. 225 (2007), no. 2, 2138–2174. MR 2349699, DOI 10.1016/j.jcp.2007.03.012
- Zhiming Chen and Jun Zou, Finite element methods and their convergence for elliptic and parabolic interface problems, Numer. Math. 79 (1998), no. 2, 175–202. MR 1622502, DOI 10.1007/s002110050336
- Maksymilian Dryja, On discontinuous Galerkin methods for elliptic problems with discontinuous coefficients, Comput. Methods Appl. Math. 3 (2003), no. 1, 76–85. Dedicated to Raytcho Lazarov. MR 2002258, DOI 10.2478/cmam-2003-0007
- Miguel A. Dumett and James P. Keener, An immersed interface method for solving anisotropic elliptic boundary value problems in three dimensions, SIAM J. Sci. Comput. 25 (2003), no. 1, 348–367. MR 2047209, DOI 10.1137/S106482750240697X
- Yalchin Efendiev and Thomas Y. Hou, Multiscale finite element methods, Surveys and Tutorials in the Applied Mathematical Sciences, vol. 4, Springer, New York, 2009. Theory and applications. MR 2477579
- Ronald P. Fedkiw, Tariq Aslam, Barry Merriman, and Stanley Osher, A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method), J. Comput. Phys. 152 (1999), no. 2, 457–492. MR 1699710, DOI 10.1006/jcph.1999.6236
- David Gilbarg and Neil S. Trudinger, Elliptic partial differential equations of second order, Grundlehren der Mathematischen Wissenschaften, Vol. 224, Springer-Verlag, Berlin-New York, 1977. MR 0473443, DOI 10.1007/978-3-642-96379-7
- I. G. Graham, P. O. Lechner, and R. Scheichl, Domain decomposition for multiscale PDEs, Numer. Math. 106 (2007), no. 4, 589–626. MR 2317926, DOI 10.1007/s00211-007-0074-1
- I. G. Graham and R. Scheichl, Robust domain decomposition algorithms for multiscale PDEs, Numer. Methods Partial Differential Equations 23 (2007), no. 4, 859–878. MR 2326197, DOI 10.1002/num.20254
- P. Grisvard, Elliptic problems in nonsmooth domains, Monographs and Studies in Mathematics, vol. 24, Pitman (Advanced Publishing Program), Boston, MA, 1985. MR 775683
- Thomas Y. Hou, Zhilin Li, Stanley Osher, and Hongkai Zhao, A hybrid method for moving interface problems with application to the Hele-Shaw flow, J. Comput. Phys. 134 (1997), no. 2, 236–252. MR 1458828, DOI 10.1006/jcph.1997.5689
- Thomas Y. Hou and Xiao-Hui Wu, A multiscale finite element method for elliptic problems in composite materials and porous media, J. Comput. Phys. 134 (1997), no. 1, 169–189. MR 1455261, DOI 10.1006/jcph.1997.5682
- Thomas Y. Hou, Xiao-Hui Wu, and Zhiqiang Cai, Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients, Math. Comp. 68 (1999), no. 227, 913–943. MR 1642758, DOI 10.1090/S0025-5718-99-01077-7
- Jianguo Huang and Jun Zou, Some new a priori estimates for second-order elliptic and parabolic interface problems, J. Differential Equations 184 (2002), no. 2, 570–586. MR 1929889, DOI 10.1006/jdeq.2001.4154
- T. G. Liu, B. C. Khoo, and C. W. Wang, The ghost fluid method for compressible gas-water simulation, J. Comput. Phys. 204 (2005), no. 1, 193–221. MR 2121909, DOI 10.1016/j.jcp.2004.10.012
- B.C. Khoo, T.G. Liu and K.S. Yeo, Ghost fluid method for strong shock impacting on material interface, J. Comput. Phys. 190 (2003), 651-681.
- Randall J. LeVeque and Zhi Lin Li, The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM J. Numer. Anal. 31 (1994), no. 4, 1019–1044. MR 1286215, DOI 10.1137/0731054
- J. Li, J.M. Melenk, B. Wohlmuth and J. Zou, Optimal a priori estimates for higher order finite element methods for elliptic interface problems, Applied Numerical Mathematics (in press) (2009).
- Zhilin Li, A fast iterative algorithm for elliptic interface problems, SIAM J. Numer. Anal. 35 (1998), no. 1, 230–254. MR 1618460, DOI 10.1137/S0036142995291329
- Zhilin Li and Kazufumi Ito, Maximum principle preserving schemes for interface problems with discontinuous coefficients, SIAM J. Sci. Comput. 23 (2001), no. 1, 339–361. MR 1860918, DOI 10.1137/S1064827500370160
- Zhilin Li, Tao Lin, and Xiaohui Wu, New Cartesian grid methods for interface problems using the finite element formulation, Numer. Math. 96 (2003), no. 1, 61–98. MR 2018791, DOI 10.1007/s00211-003-0473-x
- J.L. Lions and E. Magenes, Non-Homogeneous Boundary Value Problems and Applications, Springer-Verlag, Berlin, 1972.
- Xu-Dong Liu, Ronald P. Fedkiw, and Myungjoo Kang, A boundary condition capturing method for Poisson’s equation on irregular domains, J. Comput. Phys. 160 (2000), no. 1, 151–178. MR 1756763, DOI 10.1006/jcph.2000.6444
- Xu-Dong Liu and Thomas C. Sideris, Convergence of the ghost fluid method for elliptic equations with interfaces, Math. Comp. 72 (2003), no. 244, 1731–1746. MR 1986802, DOI 10.1090/S0025-5718-03-01525-4
- Hyeonbae Kang and Jin Keun Seo, The layer potential technique for the inverse conductivity problem, Inverse Problems 12 (1996), no. 3, 267–278. MR 1391539, DOI 10.1088/0266-5611/12/3/007
- Hyeonbae Kang, Jin Keun Seo, and Dongwoo Sheen, Numerical identification of discontinuous conductivity coefficients, Inverse Problems 13 (1997), no. 1, 113–123. MR 1435871, DOI 10.1088/0266-5611/13/1/009
- William McLean, Strongly elliptic systems and boundary integral equations, Cambridge University Press, Cambridge, 2000. MR 1742312
- Houman Owhadi and Lei Zhang, Metric-based upscaling, Comm. Pure Appl. Math. 60 (2007), no. 5, 675–723. MR 2292954, DOI 10.1002/cpa.20163
- Clemens Pechstein and Robert Scheichl, Analysis of FETI methods for multiscale PDEs, Numer. Math. 111 (2008), no. 2, 293–333. MR 2456834, DOI 10.1007/s00211-008-0186-2
- Charles S. Peskin, Numerical analysis of blood flow in the heart, J. Comput. Phys. 25 (1977), no. 3, 220–252. MR 490027, DOI 10.1016/0021-9991(77)90100-0
- Charles S. Peskin, The immersed boundary method, Acta Numer. 11 (2002), 479–517. MR 2009378, DOI 10.1017/S0962492902000077
- Michael Plum and Christian Wieners, Optimal a priori estimates for interface problems, Numer. Math. 95 (2003), no. 4, 735–759. MR 2013126, DOI 10.1007/s002110200395
- Michael E. Taylor, Partial differential equations, Texts in Applied Mathematics, vol. 23, Springer-Verlag, New York, 1996. Basic theory. MR 1395147, DOI 10.1007/978-1-4684-9320-7
- Andrea Toselli and Olof Widlund, Domain decomposition methods—algorithms and theory, Springer Series in Computational Mathematics, vol. 34, Springer-Verlag, Berlin, 2005. MR 2104179, DOI 10.1007/b137868
- S.O. Unverdi and G. Tryggvason, A front-tracking method for viscous, incompressible, multi-fluid flows, J. Comput. Phys. 100 (1992), 25-37.
- T. von Petersdorff, Randwertprobleme der Elastizitätstheorie für Polyeder Singularitätaten und Approximation mit Randelementmethoden, Ph.D. Thesis, T.U. Darmstadt, German, 1989.
- Martin Vohralík, Two types of guaranteed (and robust) a posteriori estimates for finite volume methods, Finite volumes for complex applications V, ISTE, London, 2008, pp. 649–656. MR 2451464
- Y. C. Zhou, Shan Zhao, Michael Feig, and G. W. Wei, High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources, J. Comput. Phys. 213 (2006), no. 1, 1–30. MR 2186592, DOI 10.1016/j.jcp.2005.07.022
Additional Information
- C.-C. Chu
- Affiliation: Department of Applied and Computational Mathematics, California Institute of Technology, Pasadena, California 91125
- Address at time of publication: Department of Mathematics, University of Texas at Austin, 1 University Station C1200, Austin, Texas 78712
- Email: ccchu@acm.caltech.edu
- I. G. Graham
- Affiliation: Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom
- MR Author ID: 76020
- Email: I.G.Graham@bath.ac.uk
- T.-Y. Hou
- Affiliation: Department of Applied and Computational Mathematics, California Insitute of Technology, Pasadena, California 91125
- Email: hou@acm.caltech.edu
- Received by editor(s): February 24, 2009
- Published electronically: May 25, 2010
- Additional Notes: The authors thank Rob Scheichl and Jens Markus Melenk for useful discussions. The second author acknowledges financial support from the Applied and Computational Mathematics Group at California Institute of Technology. The research of the third author was supported in part by an NSF Grant DMS-0713670 and a DOE Grant DE-FG02-06ER25727.
- © Copyright 2010 American Mathematical Society
- Journal: Math. Comp. 79 (2010), 1915-1955
- MSC (2010): Primary 65N12, 65N30
- DOI: https://doi.org/10.1090/S0025-5718-2010-02372-5
- MathSciNet review: 2684351