Weakly Imposed Symmetry and Robust Preconditioners for Biot’s Consolidation Model

Trygve Bærland; Jeonghun J. Lee; Kent-Andre Mardal; Ragnar Winther

doi:10.1515/cmam-2017-0016

Artikel

Weakly Imposed Symmetry and Robust Preconditioners for Biot’s Consolidation Model

Trygve Bærland , Jeonghun J. Lee , Kent-Andre Mardal und Ragnar Winther

Veröffentlicht/Copyright: 17. Juni 2017

Veröffentlicht von

Veröffentlichen auch Sie bei De Gruyter Brill

Manuskript einreichen Informationen für Autor*innen Erkunden Sie dieses Fachgebiet

Aus der Zeitschrift Computational Methods in Applied Mathematics Band 17 Heft 3

Abstract

We discuss the construction of robust preconditioners for finite element approximations of Biot’s consolidation model in poroelasticity. More precisely, we study finite element methods based on generalizations of the Hellinger–Reissner principle of linear elasticity, where the stress tensor is one of the unknowns. The Biot model has a number of applications in science, medicine, and engineering. A challenge in many of these applications is that the model parameters range over several orders of magnitude. Therefore, discretization procedures which are well behaved with respect to such variations are needed. The focus of the present paper will be on the construction of preconditioners, such that the preconditioned discrete systems are well-conditioned with respect to variations of the model parameters as well as refinements of the discretization. As a byproduct, we also obtain preconditioners for linear elasticity that are robust in the incompressible limit.

Keywords: Poroelasticity; Finite Elements; Preconditioning

MSC 2010: 65N30; 65N55; 74S05

Funding statement: The research leading to these results has received funding the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 339643. The work of Kent-Andre Mardal has also been supported by the Research Council of Norway through grant no. 179578/F50.

A Appendix: Proofs of (3.3) and (3.5)

Here we provide proofs of inequalities (3.3) and (3.5).

Proof of (3.3).

Fix τ∈Σ and recall that |Γt|>0 and τ⋅ν^=0 on Γt. By the pointwise decomposition

τ=PD⁢τ+(I-PD)⁢τ,

and the fact that (I-PD)⁢τ=1n⁢tr⁡τ⁢𝕀, it suffices to show that

∥tr⁡τ∥02≤C⁢((12⁢μ⁢PD⁢τ,PD⁢τ)+∥𝐝𝐢𝐯⁡τ∥02)

for some constant C independent of τ. To prove this, we use a well-known result for the right inverse of the divergence operator: There exists ϕ∈HΓd1⁢(Ω;𝕍):={φ∈H1⁢(Ω;𝕍):φ|Γd=0} such that

(A.1)

div⁡ϕ=tr⁡τ,∥ϕ∥1≤C⁢∥tr⁡τ∥0

with C>0 independent of τ, cf. Appendix B. We then have that

∥tr⁡τ∥02=(tr⁡τ,div⁡ϕ)=(tr⁡τ⁢𝕀,𝐠𝐫𝐚𝐝⁡ϕ).

Since tr⁡τ⁢𝕀=n⁢(τ-PD⁢τ), we get

∥tr⁡τ∥02=n⁢(τ,𝐠𝐫𝐚𝐝⁡ϕ)-n⁢(PD⁢τ,𝐠𝐫𝐚𝐝⁡ϕ)

=-n⁢(𝐝𝐢𝐯⁡τ,ϕ)-n⁢(PD⁢τ,𝐠𝐫𝐚𝐝⁡ϕ),

where the first term of the final form is a result of integration by parts. Next, we may use the Cauchy–Schwarz inequality, which results in

∥tr⁡τ∥02≤n⁢(∥𝐝𝐢𝐯⁡τ∥0⁢∥ϕ∥0+∥PD⁢τ∥0⁢∥grad⁡ϕ∥0)

≤n⁢(∥𝐝𝐢𝐯⁡τ∥02+∥PD⁢τ∥02)12⁢∥ϕ∥1

≤C⁢(∥𝐝𝐢𝐯⁡τ∥02+∥PD⁢τ∥02)12⁢∥tr⁡τ∥0,

and so the result follows after dividing by ∥tr⁡τ∥0. ∎

When |Γt|=0, i.e., Γd=∂⁡Ω, (A.1) can only hold if tr⁡τ has mean value zero. However, with this constraint, we can prove (3.5) with almost the same argument as above.

Proof of (3.5).

Fix any τ∈Σ. From the decomposition P0⁢τ=PD⁢τ+(P0-PD)⁢τ, it suffices to prove the estimate for (P0-PD)⁢τ component. Denoting the mean value of the trace by

tr⁡τ¯:=1|Ω|⁢∫Ωtr⁡τ⁢d⁢x,

we have that (I-PD)⁢P0⁢τ=(P0-PD)⁢τ=1n⁢(tr⁡τ-tr⁡τ¯)⁢𝕀, and so it is sufficient to show that

∥tr⁡τ-tr⁡τ¯∥02≤C⁢((12⁢μ⁢PD⁢τ,PD⁢τ)+∥𝐝𝐢𝐯⁡τ∥02).

Since tr⁡τ-tr⁡τ¯ is mean-value zero, there exists ϕ∈{φ∈H1⁢(Ω;𝕍):φ|∂⁡Ω=0} such that

div⁡ϕ=tr⁡τ-tr⁡τ¯,∥ϕ∥1≤C⁢∥tr⁡τ-tr⁡τ¯∥0

with C>0 independent of τ (cf. [14, Theorem 5.1]). The rest of the proof is completely analogous to the proof of (3.3) above. ∎

B Appendix: Right Inverse of Divergence Operator

A result for the right inverse of the divergence operator, as expressed by (A.1), is closely related to the inf-sup condition for the Stokes problem, and therefore well-known. However, we are not aware of a proper reference for the case when |∂⁡Ω|>|Γt|>0, i.e., for the case when |Γd|>0, but Γd is not all of ∂⁡Ω. Therefore, for completeness, we include a proof here.

Lemma B.1.

Assume that |Γt|>0 and set HΓd1⁢(Ω;V)={ϕ∈H1⁢(Ω;V):ϕ|Γd=0}. Then there is a constant C>0 so that for every f∈L2⁢(Ω) there is a ϕ∈HΓd1⁢(Ω;V) so that

div⁡ϕ=f,∥ϕ∥1≤C⁢∥f∥0.

Proof.

Take any f∈L2⁢(Ω). We first decompose f into its mean value zero- and mean value part as f=f0+fc, where f0∈L02⁢(Ω) and fc=af⁢1Ω for af∈ℝ. Further, we can decompose HΓd1⁢(Ω;𝕍)=H01⁢(Ω;𝕍)⊕V1, where

V1:={ϕ∈HΓd1⁢(Ω;𝕍):(𝐠𝐫𝐚𝐝⁡ϕ,𝐠𝐫𝐚𝐝⁡ψ)=0⁢ for all ⁢ψ∈H01⁢(Ω;𝕍)}.

Consider then the problem of finding ζ∈V1 so that

(B.1)(𝐠𝐫𝐚𝐝⁡ζ,𝐠𝐫𝐚𝐝⁡ψ)=(𝕀,𝐠𝐫𝐚𝐝⁡ψ) for all ⁢ψ∈V1.

By the Lax–Milgram lemma (cf. e.g. [8, Theorem 4.1.6]), problem (B.1) has a unique solution ζ and ∥ζ∥1≤C1 for some constant C1>0 depending on Ω. Taking ψ=ζ in (B.1), we obtain

∫Ωdiv⁡ζ⁢d⁢x=∥𝐠𝐫𝐚𝐝⁡ζ∥02.

Therefore, if we set

ω=af∥𝐠𝐫𝐚𝐝⁡ζ∥02⁢ζ,

we have

∫Ωdiv⁡ω⁢d⁢x=af

and ∥ω∥1≤C⁢∥fc∥0 for some constant C depending on ζ. It follows that f-div⁡ω∈L02⁢(Ω), i.e., f-div⁡ω has mean value zero. From the theory of Stokes equation, we can thus find a ω0∈H01⁢(Ω;𝕍) so that

(B.2)div⁡ω0=f-div⁡ω,∥ω0∥1≤C2⁢∥f-div⁡ω∥0,

where the constant C2 is independent of f-div⁡ω (cf. [14, Theorem 5.1]). We set ϕ=ω0+ω, and it follows from (B.2) that div⁡ϕ=f. Using the triangle inequality, (B.2) and the properties of ω, we estimate ∥ϕ∥1 as

∥ϕ∥1≤∥ω0∥1+∥ω∥1≤C⁢(∥f-div⁡ω∥0+∥fc∥0)≤C⁢(∥f∥0+∥ω∥1)≤C⁢∥f∥0,

which completes the proof. ∎

References

[1] D. N. Arnold, F. Brezzi, J. Douglas, Jr., PEERS: A new mixed finite element for plane elasticity, Japan J. Appl. Math. 1 (1984), no. 2, 347–367. 10.1007/BF03167064Suche in Google Scholar

[2] D. N. Arnold, J. Douglas, Jr. and C. P. Gupta, A family of higher order mixed finite element methods for plane elasticity, Numer. Math. 45 (1984), no. 1, 1–22. 10.1007/BF01379659Suche in Google Scholar

[3] D. N. Arnold, R. S. Falk and R. Winther, Preconditioning in H⁢(div) and applications, Math. Comp. 66 (1997), no. 219, 957–984. 10.1090/S0025-5718-97-00826-0Suche in Google Scholar

[4] D. N. Arnold, R. S. Falk and R. Winther, Mixed finite element methods for linear elasticity with weakly imposed symmetry, Math. Comp. 76 (2007), no. 260, 1699–1723. 10.1090/S0025-5718-07-01998-9Suche in Google Scholar

[5] O. Axelsson, R. Blaheta and P. Byczanski, Stable discretization of poroelasticity problems and efficient preconditioners for arising saddle point type matrices, Comput. Vis. Sci. 15 (2012), no. 4, 191–207. 10.1007/s00791-013-0209-0Suche in Google Scholar

[6] L. Berger, R. Bordas, D. Kay and S. Tavener, Stabilized lowest-order finite element approximation for linear three-field poroelasticity, SIAM J. Sci. Comput. 37 (2015), no. 5, A2222–A2245. 10.1137/15M1009822Suche in Google Scholar

[7] D. Boffi, F. Brezzi and M. Fortin, Reduced symmetry elements in linear elasticity, Commun. Pure Appl. Anal. 8 (2009), no. 1, 95–121. 10.3934/cpaa.2009.8.95Suche in Google Scholar

[8] D. Boffi, F. Brezzi and M. Fortin, Mixed Finite Element Methods and Applications, Springer Ser. Comput. Math. 44, Springer, Heidelberg, 2013. 10.1007/978-3-642-36519-5Suche in Google Scholar

[9] J. H. Bramble, Multigrid Methods, Pitman Res. Notes in Math. 294, Longman Scientific & Technical, Harlow, 1993. Suche in Google Scholar

[10] Y. Chen, Y. Luo and M. Feng, Analysis of a discontinuous Galerkin method for the Biot’s consolidation problem, Appl. Math. Comput. 219 (2013), no. 17, 9043–9056. 10.1016/j.amc.2013.03.104Suche in Google Scholar

[11] B. Cockburn, J. Gopalakrishnan and J. Guzmán, A new elasticity element made for enforcing weak stress symmetry, Math. Comp. 79 (2010), no. 271, 1331–1349. 10.1090/S0025-5718-10-02343-4Suche in Google Scholar

[12] O. Coussy, Poromechanics, John Wiley & Sons, Hoboken, 2004. 10.1002/0470092718Suche in Google Scholar

[13] M. Farhloul and M. Fortin, Dual hybrid methods for the elasticity and the Stokes problems: A unified approach, Numer. Math. 76 (1997), no. 4, 419–440. 10.1007/s002110050270Suche in Google Scholar

[14] V. Girault and P.-A. Raviart, Finite Element Methods for Navier–Stokes Equations, Springer Ser. Comput. Math. 5, Springer, Berlin, 1986. 10.1007/978-3-642-61623-5Suche in Google Scholar

[15] J. Gopalakrishnan and J. Guzmán, A second elasticity element using the matrix bubble, IMA J. Numer. Anal. 32 (2012), no. 1, 352–372. 10.1093/imanum/drq047Suche in Google Scholar

[16] J. B. Haga, H. Osnes and H. P. Langtangen, A parallel block preconditioner for large-scale poroelasticity with highly heterogeneous material parameters, Comput. Geosci. 16 (2012), no. 3, 723–734. 10.1007/s10596-012-9284-4Suche in Google Scholar

[17] R. Hiptmair and J. Xu, Nodal auxiliary space preconditioning in H⁢(curl) and H⁢(div) spaces, SIAM J. Numer. Anal. 45 (2007), no. 6, 2483–2509. 10.1137/060660588Suche in Google Scholar

[18] J. J. Lee, Robust error analysis of coupled mixed methods for Biot’s consolidation model, J. Sci. Comput. 69 (2016), no. 2, 610–632. 10.1007/s10915-016-0210-0Suche in Google Scholar

[19] J. J. Lee, Towards a unified analysis of mixed methods for elasticity with weakly symmetric stress, Adv. Comput. Math. 42 (2016), no. 2, 361–376. 10.1007/s10444-015-9427-ySuche in Google Scholar

[20] J. J. Lee, K.-A. Mardal and R. Winther, Parameter-robust discretization and preconditioning of Biot’s consolidation model, SIAM J. Sci. Comput. 39 (2017), no. 1, A1–A24. 10.1137/15M1029473Suche in Google Scholar

[21] K.-A. Mardal and R. Winther, Preconditioning discretizations of systems of partial differential equations, Numer. Linear Algebra Appl. 18 (2011), no. 1, 1–40. 10.1002/nla.716Suche in Google Scholar

[22] M. A. Murad and A. F. D. Loula, Improved accuracy in finite element analysis of Biot’s consolidation problem, Comput. Methods Appl. Mech. Engrg. 95 (1992), no. 3, 359–382. 10.1016/0045-7825(92)90193-NSuche in Google Scholar

[23] M. A. Murad and A. F. D. Loula, On stability and convergence of finite element approximations of Biot’s consolidation problem, Internat. J. Numer. Methods Engrg. 37 (1994), no. 4, 645–667. 10.1002/nme.1620370407Suche in Google Scholar

[24] P. J. Phillips and M. F. Wheeler, A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: The continuous in time case, Comput. Geosci. 11 (2007), no. 2, 131–144. 10.1007/s10596-007-9045-ySuche in Google Scholar

[25] P. J. Phillips and M. F. Wheeler, A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: The discrete-in-time case, Comput. Geosci. 11 (2007), no. 2, 145–158. 10.1007/s10596-007-9044-zSuche in Google Scholar

[26] P. J. Phillips and M. F. Wheeler, A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity, Comput. Geosci. 12 (2008), no. 4, 417–435. 10.1007/s10596-008-9082-1Suche in Google Scholar

[27] M. B. Reed, An investigation of numerical errors in the analysis of consolidation by finite elements, Internat. J. Numer. Analyt. Methods Geomech. 8 (1984), no. 3, 243–257. 10.1002/nag.1610080304Suche in Google Scholar

[28] S. Rhebergen, G. N. Wells, A. J. Wathen and R. F. Katz, Three-field block preconditioners for models of coupled magma/mantle dynamics, SIAM J. Sci. Comput. 37 (2015), no. 5, A2270–A2294. 10.1137/14099718XSuche in Google Scholar

[29] R. E. Showalter, Diffusion in poro-elastic media, J. Math. Anal. Appl. 251 (2000), no. 1, 310–340. 10.1006/jmaa.2000.7048Suche in Google Scholar

[30] J. H. Smith and J. A. Humphrey, Interstitial transport and transvascular fluid exchange during infusion into brain and tumor tissue, Microvasc. Res. 73 (2007), no. 1, 58–73. 10.1016/j.mvr.2006.07.001Suche in Google Scholar PubMed

[31] R. Stenberg, A family of mixed finite elements for the elasticity problem, Numer. Math. 53 (1988), no. 5, 513–538. 10.1007/BF01397550Suche in Google Scholar

[32] K. H. Støverud, M. Alnæs, H. P. Langtangen, V. Haughton and K.-A. Mardal, Poro-elastic modeling of syringomyelia – A systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord, Comput. Methods Biomech. Biomed. Engin. 19 (2016), no. 6, 686–698. 10.1080/10255842.2015.1058927Suche in Google Scholar PubMed

[33] P. A. Vermeer and A. Verruijt, An accuracy condition for consolidation by finite elements, Internat. J. Numer. Analyt. Methods Geomech. 5 (1981), no. 1, 1–14. 10.1002/nag.1610050103Suche in Google Scholar

[34] H. F. Wang, Theory of Linear Poroelasticity, Princeton Ser. Geophys., Princeton University Press, Princeton, 2000. Suche in Google Scholar

[35] S.-Y. Yi, A coupling of nonconforming and mixed finite element methods for Biot’s consolidation model, Numer. Methods Partial Differ. Equ. 29 (2013), no. 5, 1749–1777. 10.1002/num.21775Suche in Google Scholar

[36] S.-Y. Yi, Convergence analysis of a new mixed finite element method for Biot’s consolidation model, Numer. Methods Partial Differ. Equ. 30 (2014), no. 4, 1189–1210. 10.1002/num.21865Suche in Google Scholar

[37] O. C. Zienkiewicz and T. Shiomi, Dynamic behaviour of saturated porous media; the generalized Biot formulation and its numerical solution, Internat. J. Numer. Analyt. Methods Geomech. 8 (1984), no. 1, 71–96. 10.1002/nag.1610080106Suche in Google Scholar

Received: 2017-3-15

Revised: 2017-5-15

Accepted: 2017-5-16

Published Online: 2017-6-17

Published in Print: 2017-7-1

Sie haben derzeit keinen Zugang zu diesem Inhalt.

Artikel in diesem Heft

https://doi.org/10.1515/cmam-2017-0016

Schlagwörter für diesen Artikel

Poroelasticity; Finite Elements; Preconditioning