Skip to main content
×
Home
    • Aa
    • Aa

Polyharmonic homogenization, rough polyharmonic splines and sparse super-localization

  • Houman Owhadi (a1), Lei Zhang (a2) and Leonid Berlyand (a3)
Abstract

We introduce a new variational method for the numerical homogenization of divergence form elliptic, parabolic and hyperbolic equations with arbitrary rough (L) coefficients. Our method does not rely on concepts of ergodicity or scale-separation but on compactness properties of the solution space and a new variational approach to homogenization. The approximation space is generated by an interpolation basis (over scattered points forming a mesh of resolution H) minimizing the L2 norm of the source terms; its (pre-)computation involves minimizing 𝒪(Hd) quadratic (cell) problems on (super-)localized sub-domains of size 𝒪(H ln(1/H)). The resulting localized linear systems remain sparse and banded. The resulting interpolation basis functions are biharmonic for d ≤ 3, and polyharmonic for d ≥ 4, for the operator −div(a∇·) and can be seen as a generalization of polyharmonic splines to differential operators with arbitrary rough coefficients. The accuracy of the method (𝒪(H) in energy norm and independent from aspect ratios of the mesh formed by the scattered points) is established via the introduction of a new class of higher-order Poincaré inequalities. The method bypasses (pre-)computations on the full domain and naturally generalizes to time dependent problems, it also provides a natural solution to the inverse problem of recovering the solution of a divergence form elliptic equation from a finite number of point measurements.

We introduce a new variational method for the numerical homogenization of divergence form elliptic, parabolic and hyperbolic equations with arbitrary rough (L) coefficients. Our method does not rely on concepts of ergodicity or scale-separation but on compactness properties of the solution space and a new variational approach to homogenization. The approximation space is generated by an interpolation basis (over scattered points forming a mesh of resolution H) minimizing the L2 norm of the source terms; its (pre-)computation involves minimizing 𝒪(Hd) quadratic (cell) problems on (super-)localized sub-domains of size 𝒪(H ln(1/H)). The resulting localized linear systems remain sparse and banded. The resulting interpolation basis functions are biharmonic for d ≤ 3, and polyharmonic for d ≥ 4, for the operator −div(a∇·) and can be seen as a generalization of polyharmonic splines to differential operators with arbitrary rough coefficients. The accuracy of the method (𝒪(H) in energy norm and independent from aspect ratios of the mesh formed by the scattered points) is established via the introduction of a new class of higher-order Poincaré inequalities. The method bypasses (pre-)computations on the full domain and naturally generalizes to time dependent problems, it also provides a natural solution to the inverse problem of recovering the solution of a divergence form elliptic equation from a finite number of point measurements.

Copyright
Linked references
Hide All

This list contains references from the content that can be linked to their source. For a full set of references and notes please see the PDF or HTML where available.

A. Abdulle and M.J. Grote , Finite element heterogeneous multiscale method for the wave equation. Multiscale Model. Simul. 9 (2011) 766792.

A. Abdulle and Ch. Schwab , Heterogeneous multiscale FEM for diffusion problems on rough surfaces. Multiscale Model. Simul. 3 (2004) 195220.

G. Allaire and R. Brizzi , A multiscale finite element method for numerical homogenization. Multiscale Model. Simul. 4 (2005) 790812.

T. Arbogast and K.J. Boyd . Subgrid upscaling and mixed multiscale finite elements. SIAM J. Numer. Anal. 44 (2006) 11501171.

T. Arbogast , C.-S. Huang and S.-M. Yang , Improved accuracy for alternating-direction methods for parabolic equations based on regular and mixed finite elements. Math. Models Methods Appl. Sci. 17 (2007) 12791305.

S.N. Armstrong and P.E. Souganidis , Stochastic homogenization of Hamilton-Jacobi and degenerate Bellman equations in unbounded environments. J. Math. Pures Appl. 97 (2012) 460504.

I. Babuška , G. Caloz and J.E. Osborn , Special finite element methods for a class of second order elliptic problems with rough coefficients. SIAM J. Numer. Anal. 31 (1994) 945981.

I. Babuška and R. Lipton , Optimal local approximation spaces for generalized finite element methods with application to multiscale problems. Multiscale Model. Simul. 9 (2011) 373406.

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) 510536.

I. Babuška and J.E. Osborn , Can a finite element method perform arbitrarily badly? Math. Comput. 69 (2000) 443462.

G. Bal and W. Jing , Corrector theory for MSFEM and HMM in random media. Multiscale Model. Simul. 9 (2011) 15491587.

G. Ben Arous and H. Owhadi , Multiscale homogenization with bounded ratios and anomalous slow diffusion. Comm. Pure Appl. Math. 56 (2003) 80113.

L. Berlyand and H. Owhadi , Flux norm approach to finite dimensional homogenization approximations with non-separated scales and high contrast. Arch. Rational Mech. Anal. 198 (2010) 677721.

X. Blanc , C. Le Bris and P.-L. Lions , Une variante de la théorie de l’homogénéisation stochastique des opérateurs elliptiques. C. R. Math. Acad. Sci. Paris 343 (2006) 717724.

X. Blanc , C. Le Bris and P.-L. Lions , Stochastic homogenization and random lattices. J. Math. Pures Appl. 88 (2007) 3463.

L.A. Caffarelli and P.E. Souganidis , A rate of convergence for monotone finite difference approximations to fully nonlinear, uniformly elliptic PDEs. Comm. Pure Appl. Math. 61 (2008) 117.

C.-C. Chu , I.G. Graham and T.Y. Hou , A new multiscale finite element method for high-contrast elliptic interface problems. Math. Comput. 79 (2010) 19151955.

Y. Efendiev , J. Galvis and X. Wu , Multiscale finite element and domain decomposition methods for high-contrast problems using local spectral basis functions. J. Comput. Phys. 230 (2011) 937955.

Y. Efendiev , V. Ginting , T. Hou and R. Ewing , Accurate multiscale finite element methods for two-phase flow simulations. J. Comput. Phys. 220 (2006) 155174.

Y. Efendiev and T. Hou , Multiscale finite element methods for porous media flows and their applications. Appl. Numer. Math. 57 (2007) 577596.

B. Engquist and P.E. Souganidis , Asymptotic and numerical homogenization. Acta Numerica 17 (2008) 147190.

B. Engquist , H. Holst and O. Runborg , Multi-scale methods for wave propagation in heterogeneous media. Commun. Math. Sci. 9 (2011) 3356.

A. Gloria , Analytical framework for the numerical homogenization of elliptic monotone operators and quasiconvex energies. SIAM MMS 5 (2006) 9961043.

A. Gloria , Reduction of the resonance error-Part 1: Approximation of homogenized coefficients. Math. Models Methods Appl. Sci. 21 (2011) 16011630.

A. Gloria and F. Otto , An optimal error estimate in stochastic homogenization of discrete elliptic equations. Ann. Appl. Probab. 22 (2012) 128.

L. Grasedyck , I. Greff and S. Sauter , The al basis for the solution of elliptic problems in heterogeneous media. Multiscale Modeling and Simulation 10 (2012) 245258.

M. Grüter and K. Widman , The green function for uniformly elliptic equations. Manuscripta Math. 37 (1982) 303342.

R.L. Harder and R.N. Desmarais , Interpolation using surface splines. J. Aircr. 9 (1972) 189191.

T.Y. Hou , X.-H. Wu and Z. Cai , Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients. Math. Comput. 68 (1999) 913943.

T.Y. Hou and X.H. Wu , A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys. 134 (1997) 169189.

E. Kosygina , F. Rezakhanlou and S.R.S. Varadhan , Stochastic homogenization of Hamilton-Jacobi-Bellman equations. Comm. Pure Appl. Math. 59 (2006) 14891521.

O. Kounchev and H. Render , Polyharmonic splines on grids Z× aZn and their limits. Math. Comput. 74 (2005) 18311841.

P.-L. Lions and P.E. Souganidis , Correctors for the homogenization of Hamilton-Jacobi equations in the stationary ergodic setting. Comm. Pure Appl. Math. 56 (2003) 15011524.

W.R. Madych and S.A. Nelson , Multivariate interpolation and conditionally positive definite functions. II. Math. Comput. 54 (1990) 211230.

W.R. Madych and S.A. Nelson , Polyharmonic cardinal splines. J. Approx. Theory 60 (1990) 141156.

W.R. Madych and S.A. Nelson , Polyharmonic cardinal splines: a minimization property. J. Approx. Theory 63 (1990) 303320.

J.M. Melenk , On n-widths for elliptic problems. J. Math. Anal. Appl. 247 (2000) 272289.

P. Ming and X. Yue , Numerical methods for multiscale elliptic problems. J. Comput. Phys. 214 (2006) 421445.

F.J. Narcowich , J.D. Ward and H. Wendland , Sobolev bounds on functions with scattered zeros, with applications to radial basis function surface fitting. Math. Comput. 74 (2005) 743763.

G. Nguetseng , A general convergence result for a functional related to the theory of homogenization. SIAM J. Math. Anal. 20 (1989) 608623.

J. Nolen , G. Papanicolaou and O. Pironneau , A framework for adaptive multiscale methods for elliptic problems. Multiscale Model. Simul. 7 (2008) 171196.

H. Owhadi and L. Zhang , Metric-based upscaling. Comm. Pure Appl. Math. 60 (2007) 675723.

H. Owhadi and L. Zhang . Localized bases for finite dimensional homogenization approximations with non-separated scales and high-contrast. SIAM Multiscale Model. Simul. 9 (2011) 13731398. arXiv:1011.0986.

H. Owhadi , Anomalous slow diffusion from perpetual homogenization. Ann. Probab. 31 (2003) 19351969.

H. Owhadi , Averaging versus chaos in turbulent transport? Comm. Math. Phys. 247 (2004) 553599.

Ch. Rabut , High level m-harmonic cardinal B-splines. Numer. Algorithms 2 (1992) 6384.

G. Stampacchia , Le problème de dirichlet pour les équations elliptiques du second ordre à coefficients discontinus. Ann. Inst. Fourier (Grenoble) 15 (1965) 189258.

J. Vybiral , Widths of embeddings in function spaces. J. Complexity 24 (2008) 545570.

X.H. Wu , Y. Efendiev and T.Y. Hou , Analysis of upscaling absolute permeability. Discrete Contin. Dyn. Syst. Ser. B 2 (2002) 185204.

Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

ESAIM: Mathematical Modelling and Numerical Analysis
  • ISSN: 0764-583X
  • EISSN: 1290-3841
  • URL: /core/journals/esaim-mathematical-modelling-and-numerical-analysis
Please enter your name
Please enter a valid email address
Who would you like to send this to? *
×

Keywords:

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 1 *
Loading metrics...

Abstract views

Total abstract views: 12 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 23rd September 2017. This data will be updated every 24 hours.