Applications of Poromechanics to Energy Engineeering

Total Page:16

File Type:pdf, Size:1020Kb

Applications of Poromechanics to Energy Engineeering Applications of Poromechanics to Energy Engineeering Chloé Arson Assistant Professor Texas A&M University, Department of Civil Engineering C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 1 / 66 Introduction Modeling Thoughts... « Toute loi physique, étant une loi approchée, est à la merci d’un progrès qui, en augmentant la précision des expériences, rendra insuffisant le degré d’approximation que comporte cette loi ; elle est essentiellement provisoire. L’appréciation de sa valeur varie d’un physicien à l’autre, au gré des moyens d’observation dont ils disposent et de l’exactitude que réclament leurs recherches ; elle est essentiellement relative. » “Every law of Physics, being approximated, is at the mercy of progress, which, by increasing the accuracy of experiments, will make the degree of approximation of this law insufficient ; [every law of Physics] is essentially temporary. The appreciation of its value varies from one physicist to the other, depending on the means of observation that they have, and on the exactness required by their research work ; [every law of Physics] is essentially relative.” Pierre Duhem, La Théorie Physique. Son Objet - Sa Structure. Chap. 5, x 3, Revue de Philosophie, 1906-1914. C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 2 / 66 Introduction 1 Thermodynamic Framework of Poromechanics 2 Performance Assessment of Heat Exchanger Piles 3 Modeling Damage in Porous Media 4 Study of the EDZ around Nuclear Waste Disposals 5 Prediction of Permeability in Fractured Rock C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 3 / 66 Poromechanics 1 Thermodynamic Framework of Poromechanics General Thermodynamic Framework Application to Constitutive Modeling of Porous Media 2 Performance Assessment of Heat Exchanger Piles A Review of Geothermal Systems Preliminary Results and Research Plans 3 Modeling Damage in Porous Media State of the Art Limitations and Challenges An Alternative for Saturated Rock : Double Effective Stress 4 Study of the EDZ around Nuclear Waste Disposals Outline of the THHMD Model Study of Nuclear Waste Disposals with the THHMD Model 5 Prediction of Permeability in Fractured Rock A New Model of Permeability for Cracked Porous Rock Results : Simulation of Triaxial Compression Tests C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 4 / 66 Poromechanics General Thermodynamic Framework Thermodynamics of Open Systems [Coussy, 2004] First Law of Thermodynamics (porous solid filled with a non reactive fluid mixture) _ _ _ _ _ K + Eint = Etot = Pmec + Echem + Q _ _ X _ _ Pmec = Pdefo + K ) Eint = Pdefo + µj Nj + Q j _ _ _ X _ _ Ψ = Eint − TS ) Ψ + T S + ST = Pdefo + µj Nj + Q j Second Law of Thermodynamics _ S_ Q > T Inequality of Clausius-Duhem X _ _ _ Φ = Pdefo + µj Nj − ST − Ψ > 0 j Thermodynamic potentials depend on state variables and internal variables. C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 5 / 66 Poromechanics General Thermodynamic Framework State Variables for Non-Isothermal Unsaturated Porous Media 2 miscible pure fluids (liquid, gas), small deformation, absence of body forces. Inequality of Clausius-Duhem [Coussy, 2004] Φ = Φs + Φl + Φg + ΦT > 0 s s s d d φl d φg dT d Ψs Φs = σ : dt + pl dt + pg dt − Ss dt − dt r Φl = [−rX pl ] · φl Vl r Φg = [−rX pg ] · φg Vg Q ΦT = − T · rX T “deformation” “stress” e σ e φl pl e φg pg TSs C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 6 / 66 Poromechanics General Thermodynamic Framework Thermodynamic Conjugation Relationships The free energy of the solid skeleton is sought in the form : e e e Ψs = Ψs( ; φl ; φg; T ; χ) χ : vector containing all internal variables of interest, e.g. damage (Ω), hardening variables such as the equivalent plastic strain (γp) Mechanical dissipation inequality (in the absence of plastic porosity changes): p e ! e ! dφe d @Ψs d @Ψs dφl @Ψs g @Ψs dT @Ψs dχ σ : + σ − + pl − + pg − − Ss + − > 0 @e @φe @φe @ @χ dt dt l dt g dt T dt dt In the absence of irreversible microstructure change... e ! e ! e @Ψ d @Ψ dφ @Ψ dφg @Ψ dT σ − s + p − s l + p − s − S + s = 0 e l e g e s @ dt @φl dt @φg dt @T dt ... from which thermodynamic conjugation relationships are deduced... @Ψ @Ψ @Ψ @Ψ σ = s ; p = s ; p = s ; S = − s e l e g e s @ @φl @φg @T ... and the reduced dissipation inequality is obtained : dp dχ @Ψ σ : + ξ · 0 ; ξ = − s dt dt > @χ C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 7 / 66 Poromechanics Application to Constitutive Modeling of Porous Media Tangent Properties [Coussy, 2004] Expression of the skeleton free energy ) REV’s tangent properties. Example : thermo-poro-elasticity (saturated porous medium) e with a free energy sought in the form Ψs = Ψs( ; p; T ) : 2 2 2 σ_ = @ Ψs : _e + @ Ψs p_ + @ Ψs T_ @e2 @e@p @e@T σ_ = D : _e − Bp_ − α (D : δ) T_ D : stiffness tensor ; B : Biot’s tensor ; α : thermal expansion coefficient p_ φ_ = b : _e + − 3α T_ N φ 1=N : inverse of Biot’s modulus ; αφ : coefficient for volumetric thermal dilation due to porosity changes only T_ S_ = (D : δ) α : _e − 3α p_ + C s φ T C : skeleton volumetric heat capacity C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 8 / 66 Poromechanics Application to Constitutive Modeling of Porous Media Flow Properties [Coussy, 2004] Conduction laws ensure the positivity of the fluid and thermal dissipations in the absence of TH couplings : “Darcy’s” law : v = φ Vr = k l :[−r p ] l l l ηl X l ) Φ = [−r p ] · v = k l :[r p ]:[r p ] ≥ 0 l X l l ηl X l X l v = φ Vr = k g :[−r p ] g g g ηg X g ) Φ = [−r p ] · v = k g :[r p ]:[r p ] ≥ 0 g X g g ηg X g X g Fourier’s law : Q = −λrX T λ ) ΦT = T rX T · rX T ≥ 0 with positive permeabilities and thermal conductivity C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 9 / 66 Poromechanics Application to Constitutive Modeling of Porous Media General case of an isothermal solid saturated with water and air Typical Constitutive Relationships (hyp. incompressible solid grains, σnet = σ − pg δ, s = pg − pl ) _ Stress state law : σ_ net = D : _ − α (D : δ) T − βs_ Heat flow equation : ∗ QT = − λT rT + hfg ρw V vap + ρvap V a ∗ + ρw CPw V w + ρw CPvap V vap + ρa CPa V a (T − T0) Moisture transfer equation : Q wtot = Q w + Q vap = ρw (−DT · rT + DP · rs − Kw · rz) Air flow equation : 1 @pa(x;T (x)) pa V a = − Ka · rT (x) − Ka · r − Ka · rz γa @T (x) γa Balance equations Solid skeleton balance equation : r · σ + f = 0 @ET Energy conservation : @t + r · QT = 0; ET = CT (T − T0) + n (1 − Sw ) ρvap hfg @ρm ∗ Moisture mass conservation : @t + r · ρw V w + V vap = 0 Air mass conservation : @ ∗ @t [n ρa (1 − Sw + HSw )] + r · [ρa V a] + r · [ρa H V w ] − r · ρw V vap = 0 C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 10 / 66 Poromechanics Application to Constitutive Modeling of Porous Media Choice of the stress state variables for unsaturated porous media number of stress variables ? A long-standing debate [Fredlund and Morgenstein 1977, Houlsby 1997] Example for incompressible solid grains : 0 Bishop’s effective stress : σ = (σ − pg δ) + χ (pg − pl ) δ independent state variables : (σ − pg δ; pg − pl ) or (pg − pl ; σ − pl δ) or (σ − pl δ; σ − pg δ) nature of the state variables ? No unique formulation [Coussy, 2004 ; Dangla, 2010]. Example for elastic isothermal unsaturated porous media, based on the separation of energies : e e Ψs ; φ, Sl = s ; φ + φ U (φ, Sl ) e 8 @ s( ,φ) σ = > @e > > < @ e,φ @(φ U(φ,Sl )) s( ) π = Sl pl + Sg pg − @φ = @φ > > > e > @U( ,φl ,φg ) : pc = − @Sl C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 11 / 66 Geothermal Piles 1 Thermodynamic Framework of Poromechanics General Thermodynamic Framework Application to Constitutive Modeling of Porous Media 2 Performance Assessment of Heat Exchanger Piles A Review of Geothermal Systems Preliminary Results and Research Plans 3 Modeling Damage in Porous Media State of the Art Limitations and Challenges An Alternative for Saturated Rock : Double Effective Stress 4 Study of the EDZ around Nuclear Waste Disposals Outline of the THHMD Model Study of Nuclear Waste Disposals with the THHMD Model 5 Prediction of Permeability in Fractured Rock A New Model of Permeability for Cracked Porous Rock Results : Simulation of Triaxial Compression Tests C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 12 / 66 Geothermal Piles A Review of Geothermal Systems Forms of Geothermal Energy [Lund, 2007] comes “from the decay of the naturally occurring isotopes of uranium, thorium and potassium in the earth” natural thermal reservoirs above 10km depth : 1:3 × 1027J ' 3 × 1017 oil barrels geothermal resources > 150oC : electrical power plants geothermal resources < 150oC : direct use for heating and cooling Geothermal Resources Direct Use of Geothermal Energy (2005) C. Arson (Texas A&M University) Poromechanics & Energy Engineering January 6th, 2012 13 / 66 Geothermal Piles A Review of Geothermal Systems GCHP with Heat Exchanger Piles ground = heat reservoir high thermal conductivity of ground and grout ) foundation piles used as heat exchangers ) thermo-mechanical behavior of the foundation ? heat stored in summer, pumped in winter ) annual energy balance ? ) design and performance ? C.
Recommended publications
  • Fully Coupled Two-Phase Flow and Poromechanics Modeling Of
    1Fully Coupled Two-Phase Flow and Poromechanics Modeling of Coalbed Methane 2Recovery: Impact of Geomechanics on Production Rate 3 4Tianran Ma a, b, c, Jonny Rutqvist c, Curtis M. Oldenburg c, Weiqun Liu a, b, Junguo Chen d, a 5 Key Laboratory of Coal-based CO2 Capture and Geological Storage, China University of 6Mining and Technology, Xuzhou, Jiangsu, China 7b State Key Laboratory for Geomechanics and Deep Underground Engineering, China 8University of Mining and Technology, Xuzhou, Jiangsu, China 9c Lawrence Berkeley National Laboratory, Earth Sciences Division, Berkeley, CA, USA 10d College of Mining and Safety Engineering, Shandong University of Science and 11Technology, Qingdao, Shandong, China 12 13 14 15 Submitted manuscript accepted for publication in 16 Journal of Natural Gas Science and Engineering 17 https://doi.org/10.1016/j.jngse.2017.05.024 18 19 Final Version Published as: 20 Ma T., Rutqvist J., Oldenburg C.M., Liu W. and Chen J., Fully coupled two-phase flow and 21 poromechanics modeling of coalbed methane recovery: Impact of geomechanics on 22 production rate. Journal of Natural Gas Science and Engineering 45, 474- 486 (2017). 23 24 25 26 27 28 29 1 1 30Abstract 31This study presents development and application of a fully coupled two-phase (methane and 32water) and poromechanics numerical model for the analysis of geomechanical impact on 33coalbed methane (CBM) production. The model considers changes in two-phase fluid flow 34properties, i.e., coal porosity, permeability, water retention, and relative permeability curves 35through changes in cleat fractures induced by effective stress variations and desorption- 36induced shrinkage. The coupled simulator is first verified for poromechamics coupling and 37simulation parameters of a CBM reservoir model are calibrated by history matching against 38one year of CBM production field data from Shanxi Province, China.
    [Show full text]
  • Coupled Multiphase Flow and Poromechanics: 10.1002/2013WR015175 a Computational Model of Pore Pressure Effects On
    PUBLICATIONS Water Resources Research RESEARCH ARTICLE Coupled multiphase flow and poromechanics: 10.1002/2013WR015175 A computational model of pore pressure effects on Key Points: fault slip and earthquake triggering New computational approach to coupled multiphase flow and Birendra Jha1 and Ruben Juanes1 geomechanics 1 Faults are represented as surfaces, Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, capable of simulating runaway slip USA Unconditionally stable sequential solution of the fully coupled equations Abstract The coupling between subsurface flow and geomechanical deformation is critical in the assess- ment of the environmental impacts of groundwater use, underground liquid waste disposal, geologic stor- Supporting Information: Readme age of carbon dioxide, and exploitation of shale gas reserves. In particular, seismicity induced by fluid Videos S1 and S2 injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present a new computational approach to Correspondence to: model coupled multiphase flow and geomechanics of faulted reservoirs. We represent faults as surfaces R. Juanes, embedded in a three-dimensional medium by using zero-thickness interface elements to accurately model [email protected] fault slip under dynamically evolving fluid pressure and fault strength. We incorporate the effect of fluid pressures from multiphase flow in the mechanical
    [Show full text]
  • Theoretical and Computational Poromechanics from 14:30 to 16:30 Via Ferrata 3, Pavia Carlo Callari, Università Del Molise
    DIPARTIMENTO DI INGEGNERIA CIVILE E ARCHITETTURA DOTTORATO IN INGEGNERIA CIVILE E ARCHITETTURA UNIVERSITÀ DEGLI STUDI DI PAVIA SHORT COURSE ON JUNE, 22, 23, 24, 25, 26 THEORETICAL AND COMPUTATIONAL POROMECHANICS FROM 14:30 TO 16:30 VIA FERRATA 3, PAVIA CARLO CALLARI, UNIVERSITÀ DEL MOLISE Detailed Outline MOTIVATIONS: The role of poromechanics in civil, environmental and medical engineering: problem analysis and material modeling PART 1: Mechanics of fully saturated porous media with compressible phases Biot's thermodynamics and constitutive equations: • Thermodynamics of barotropic and inviscid fluids • Fluid mass balance in the porous medium • Local forms of first and second principles • Clausius-Duhem inequality • Transport laws • General form of constitutive laws for the porous medium • Linear poroelasticity and extension to poro-elastoplasticity • Fluid-to-solid and solid-to-fluid limit uncoupled influences • Strain-dependent permeability models Finite element formulations for porous media: • Boundary conditions for mechanical and fluid-flow fields • Formulation of mechanical and fluid-flow problems • Extension to non-linear response Research applications to saturated porous media: • 2D response of dam and rock mass to reservoir operations • Tunnel face stability • Poroelastic damage and applications to hydrocarbon wells PART 2: Extension to multiphase fluids Biot's thermodynamics and poroelastic laws: • General hyperelastic laws for three-phase porous media • Generalized Darcy law • Constitutive equations from the theory of mixtures
    [Show full text]
  • POROMECHANICS of POROUS and FRACTURED RESERVOIRS Jishan Liu, Derek Elsworth KIGAM, Daejeon, Korea May 15-18, 2017
    POROMECHANICS OF POROUS AND FRACTURED RESERVOIRS Jishan Liu, Derek Elsworth KIGAM, Daejeon, Korea May 15-18, 2017 1. Poromechanics – Flow Properties (Jishan Liu) 1:1 Reservoir Pressure System – How to calculate overburden stress and reservoir pressure Day 11 1:2 Darcy’s Law – Permeability and its changes, reservoir classification 1:3 Mass Conservation Law – flow equations 1:4 Steady-State Behaviors – Solutions of simple flow problems 1:5 Hydraulic Diffusivity – Definition, physical meaning, and its application in reservoirs 1:6 Rock Properties – Their Dependence on Stress Conditions Day 1 ……………………………………………………………………………………………………………... 2. Poromechanics – Fluid Storage Properties (Jishan Liu) 2:1 Fluid Properties – How they change and affect flow Day 2 2:2 Mechanisms of Liquid Production or Injection 2:3 Estimation of Original Hydrocarbons in Place 2:4 Estimation of Ultimate Recovery or Injection of Hydrocarbons 2:5 Flow – Deformation Coupling in Coal 2:6 Flow – Deformation Coupling in Shale Day 2 ……………………………………………………………………………………………………………... 3. Poromechanics –Modeling Porous Medium Flows (Derek Elsworth) 3:1 Single porosity flows - Finite Element Methods [2:1] Lecture Day 3 3:2 2D Triangular Constant Gradient Elements [2:3] Lecture 3:3 Transient Behavior - Mass Matrices [2:6] Lecture 3:4 Transient Behavior - Integration in Time [2:7] Lecture 3:5 Dual-Porosity-Dual-Permeability Models [6:1] Lecture Day 3 ……………………………………………………………………………………………………………... 4. Poromechanics –Modeling Coupled Porous Medium Flow and Deformation (Derek Elsworth) 4:1 Mechanical properties – http://www.ems.psu.edu/~elsworth/courses/geoee500/GeoEE500_1.PDF Day 4 4:2 Biot consolidation – http://www.ems.psu.edu/~elsworth/courses/geoee500/GeoEE500_1.PDF 4:3 Dual-porosity poroelasticity 4:4 Mechanical deformation - 1D and 2D Elements [5:1][5:2] Lecture 4:5 Coupled Hydro-Mechanical Models [6:2] Lecture Day 4 ……………………………………………………………………………………………………………..
    [Show full text]
  • Poromechanics V Proceedings of the Fifth Biot Conference on Poromechanics
    POROMECHANICS V PROCEEDINGS OF THE FIFTH BIOT CONFERENCE ON POROMECHANICS July 10–12, 2013 Vienna, Austria SPONSORED BY Vienna University of Technology/Technische Universität Wien Engineering Mechanics Institute of the American Society of Civil Engineers EDITED BY Christian Hellmich Bernhard Pichler Dietmar Adam Published by the American Society of Civil Engineers Preface We are extremely honored, and very happy, to cordially welcome you in the cosmopolitan and culturally rich City of Vienna and at the Vienna University of Technology (TU Wien), the oldest of its type in Central Europe, for the Fifth Biot Conference on Poromechanics (BIOT-5). In 2013, this year marking the 50th anniversary of the death of Karl von Terzaghi (1883-1963), we hold BIOT-5 in memory of this “Grandfather of Poromechanics,” the founder of consolidation experiments and theoretical soil mechanics in general (see the Introduction for further historical remarks). We are very grateful that a significant number of brilliant scientists were ready to organize mini-symposia on topics that they rated as most interesting and stimulating. These mini- symposia are the natural backbone of BIOT-5’s scientific program, their titles and organizer affiliations representing the scientific and geographic breadth of BIOT-5 (see the Introduction for all names and further details). The design and set-up of BIOT-5 was strongly supported by the Local Organizing Committee, the International Scientific Committee, and the Advisory Committee, comprising all the earlier Biot conference hosts, in particular by Franz Ulm (M.I.T.), as well as by Alex Cheng (University of Mississippi) and Younane Abousleiman (University of Oklahoma), the latter two having been so instrumental in launching, with the help of Madame Nadine Biot, the Biot conference series, and staying the ever-motivated “good spirits” behind the scenes.
    [Show full text]
  • From Linear to Nonlinear Poroelasticity and Poroviscoelasticity E
    Poromechanics: from Linear to Nonlinear Poroelasticity and Poroviscoelasticity E. Bemer, M. Boutéca, O. Vincké, N. Hoteit, O. Ozanam To cite this version: E. Bemer, M. Boutéca, O. Vincké, N. Hoteit, O. Ozanam. Poromechanics: from Linear to Nonlinear Poroelasticity and Poroviscoelasticity. Oil & Gas Science and Technology - Revue d’IFP Energies nou- velles, Institut Français du Pétrole, 2001, 56 (6), pp.531-544. 10.2516/ogst:2001043. hal-02053960 HAL Id: hal-02053960 https://hal-ifp.archives-ouvertes.fr/hal-02053960 Submitted on 1 Mar 2019 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Oil & Gas Science and Technology – Rev. IFP, Vol. 56 (2001), No. 6, pp. 531-544 Copyright © 2001, Éditions Technip Poromechanics: From Linear to Nonlinear Poroelasticity and Poroviscoelasticity E. Bemer1, M. Boutéca1, O. Vincké1, N. Hoteit2 and O. Ozanam2 1 Institut français du pétrole, 1 et 4, avenue de Bois-Préau, 92852 Rueil-Malmaison Cedex - France 2 ANDRA, Parc de la Croix-Blanche, 92298 Châtenay-Malabry Cedex - France e-mail: [email protected] - [email protected] - [email protected] - [email protected] - [email protected] Résumé — Poromécanique : de la poroélasticité linéaire à la poroélasticité non linéaire et la poroviscoélasticité — Compte tenu des répercussions sur les calculs de productivité et de réserves, une modélisation fiable du comportement des roches est essentielle en ingénierie de réservoir.
    [Show full text]
  • Energy Approach to Poromechanics Energy Approach To
    Biot Medal Lecture 2003 Energy Approach to Poromechanics Olivier COUSSY NAVIER INSTITUTE Energy Approach to Poromechanics • drying and shrinkage • rocks, soils • wetting and swelling • concretes • temperature rise and build-up of pressure • woods • hydraulic diffusion and subsidence •foams • freezing and spalling • bones • capillarity and cracking • living tissues • electro-osmosis and consolidation •etc. •etc. Outlines : 1. Revisiting the Energy Approach of Saturated Poroelasticity 2. Energy Approach to Partially Saturated Poroelasticity 3. From Poromechanics to the Mechanics of Colloids 1 Revisiting Biot’s Energy Approach to Saturated Poroelasticity Open system: Element of sk fl m fl dΩ0 Solid Saturating Porous = + fluid mass material Skeleton fluid content Assuming no p dissipation σ dε + f + fl dm −dF = 0 ij ij fl fl ρfl Strain work Helmoltz Free energy infinitesimal change of Free energy Work to make the the open system supplied mass enter the by convection porous element g fldm fl Gibbs potential Separate equality for the skeleton Mesoscopic σij dεij + g fl dm fl −dF = 0 information dpfl m fl = φ ρfl dg fl = ρfl Fsk = F − ffl m fl Fluid state equations Skeleton free energy σij dεij + pfl dφ −dFsk = 0 Strain work Related to the skeleton only Lagrangean porosity 2 Linear Poroelasticity and Poroplasticity (saturated) σij dεij + pfl dφ −dFsk = 0 ∂Fsk ∂Fsk State equations: σij = pfl = ∂εij ∂φ LINEAR POROELASTICITY σ = Kε −bpfl φ − φ0 = pfl /N +b ε sij = 2G eij POROPLASTICITY p εij → εij − εij p ∂f p ∂f f (σij ,pfl ) ≤ 0 dε = dλ
    [Show full text]
  • Coupling Biot and Navier-Stokes Problems for Fluid-Poroelastic
    Coupling Biot and Navier-Stokes problems for fluid-poroelastic structure interaction Santiago Badiaa,1, Annalisa Quainib, Alfio Quarteronic,b aInternational Center for Numerical Methods in Engineering (CIMNE), Universitat Polit`ecnica de Catalunya, Jordi Girona 1-3, Edifici C1, 08034 Barcelona, Spain. bChair of Modelling and Scientific Computing (CMCS), Institute of Analysis and Scientific Computing (IACS), Ecole´ Polytechnique F´ed´erale de Lausanne, CH-1015, Lausanne, Switzerland. cMOX– Modellistica e Calcolo Scientifico Dipartimento di Matematica “F. Brioschi” Politecnico di Milano via Bonardi 9, 20133 Milano, Italy. Abstract The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier-Stokes equations for the fluid with the Biot system for the structure. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we design residual-based stabilization tech- niques for the Biot system, that have been motivated using the variational multiscale approach. Then, we state the monolithic Navier-Stokes/Biot system with the appro- priate transmission conditions on the interface. We consider both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are used and their convergence has been analyzed for a simplified problem. The efficiency of these methods has been compared for a test problem motivated from hemodynamics applications. Key words: Stabilized finite elements, Darcy’s problem, Biot system, Poromechanics, Hemodynamics, Fluid-structure interaction 2000 MSC: 65M12, 65M60 1. Introduction The arterial wall is made of three layers of tissue called intima, media, and adven- titia. From the arterial lumen (where blood flows), the blood enters the intima and, after crossing the media, serves the adventitia and muscle cells.
    [Show full text]
  • 3RD BIOT POROMECHANICS CONFERENCE University of Oklahoma, Norman, OK, May 24–27, 2005
    3RD BIOT POROMECHANICS CONFERENCE University of Oklahoma, Norman, OK, May 24–27, 2005 BIOT LECTURE: TWENTY-FIVE YEARS OF THE SLOW WAVE James G. Berryman University of California Lawrence Livermore National Laboratory Livermore, CA BRIEF HISTORY OF POROELASTICITY 1941 – Biot (quasi-statics, theory and analysis) 1944 – Frenkel (electroseismics and waves) 1951 – Gassmann (undrained behavior, theory) 1954 – Skempton (undrained behavior, experiments in soil) 1956 – Biot (waves, Lagrangian app., prediction of the slow wave) 1957 – Biot and Willis (coefficients from experiment) 1962 – Biot (reformulation of wave theory, Hamiltonian app.) 1976 – Rice and Cleary (quasi-statics, numerical methods) 1980 (Feb. 15) – Plona (slow wave first observed! but no theory!!) 1980 (Aug. 15) – Berryman (effective mass and “structure factor”) 1980 (August) – Brown (elect. tortuosity and “structure factor”) 1980 (Dec. 15) – Johnson (liquid He, 4th sound, and the slow wave) OUTLINE What is the slow wave? Why is it important? • Methods of measuring the tortuosity factor • Brief review of experimental vs. theoretical status of slow wave • Summary of the main themes • Conclusions • Biot’s (1962) Strain Energy Functional 2 2 2E = He 2Ceζ + Mζ 4µI2 − − where H, C, M, and µ are poroelastic constants, ϕ = porosity, ~u = solid frame displacement, ~uf = pore fluid displacement, ~w = ϕ(~u ~u) = relative displacement, and f − e = ~u = frame dilatation, ∇· ζ = ~w = increment of fluid content, −∇ · 1 2 I2 = e e + e e + e e (γ + ...) = a strain invariant. x y y z z x − 4 x Biot’s Equations of Dynamic Poroelasticity After first Fourier transforming from time to frequency: ω2ρ~u +(H µ) e + µ 2~u = ω2ρ ~w + C ζ, − ∇ ∇ − f ∇ ω2q(ω)~w M ζ = ω2ρ ~u C e, − ∇ − f − ∇ where ω = 2πf = angular frequency, ρ = ϕρ + (1 ϕ)ρ = the average density, fluid − solid q(ω)= ρf [α/ϕ + iF (ξ)η/κω], and α = tortuosity, κ = permeability, F (ξ)η = dynamic viscosity, p = M ~w C ~u = Mζ Ce = fluid pressure.
    [Show full text]
  • Large Deformation Poromechanics with Local Mass Conservation: an Enriched Galerkin Finite Element Framework
    Received: 24 March 2018 Revised: 17 June 2018 Accepted: 19 June 2018 DOI: 10.1002/nme.5915 RESEARCH ARTICLE Large deformation poromechanics with local mass conservation: An enriched Galerkin finite element framework Jinhyun Choo Department of Civil Engineering, The University of Hong Kong, Pokfulam, Summary Hong Kong Numerical modeling of large deformations in fluid-infiltrated porous media Correspondence must accurately describe not only geometrically nonlinear kinematics but also Jinhyun Choo, Department of Civil fluid flow in heterogeneously deforming pore structure. Accurate simulation of Engineering, The University of Hong fluid flow in heterogeneous porous media often requires a numerical method Kong, Pokfulam, Hong Kong. Email: [email protected] that features the local (elementwise) conservation property. Here, we introduce a new finite element framework for locally mass conservative solution of cou- Funding information Seed Fund for Basic Research of The pled poromechanical problems at large strains. At the core of our approach is the University of Hong Kong, Grant/Award enriched Galerkin discretization of the fluid mass balance equation, whereby Number: 201801159010 elementwise constant functions are augmented to the standard continuous Galerkin discretization. The resulting numerical method provides local mass conservation by construction with a usually affordable cost added to the con- tinuous Galerkin counterpart. Two equivalent formulations are developed using total Lagrangian and updated Lagrangian approaches. The local mass conserva- tion property of the proposed method is verified through numerical examples involving saturated and unsaturated flow in porous media at finite strains. The numerical examples also demonstrate that local mass conservation can be a crit- ical element of accurate simulation of both fluid flow and large deformation in porous media.
    [Show full text]
  • Prediction of Permeability and Formation Factor of Sandstone With
    International Journal of Rock Mechanics and Mining Sciences 106 (2018) 269–277 Contents lists available at ScienceDirect International Journal of Rock Mechanics and Mining Sciences journal homepage: www.elsevier.com/locate/ijrmms Prediction of permeability and formation factor of sandstone with hybrid T lattice Boltzmann/finite element simulation on microtomographic images ⁎ WaiChing Suna, , Teng-fong Wongb a Department of Civil Engineering and Engineering Mechanics, Columbia University, New York, USA b Earth System Science Programme, Faculty of Science, The Chinese University of Hong Kong, Hong Kong ARTICLE INFO ABSTRACT Keywords: In Fontainebleau sandstone, the evolution of transport properties with porosity is related to changes in both the Digital rock physics size and connectivity of the pore space. Microcomputed tomography can be used to characterize the relevant Computed tomography geometric attributes, with the resolution that is sufficiently refined for realistic simulation of transport properties Hydraulic permeability based on the 3D image. In this study, we adopted a hybrid computation scheme that is based on a hierarchical Formation factor multi-scale approach. The specimen was partitioned into cubic sub-volumes for pore-scale simulation of hy- draulic permeability and formation factor using the lattice Boltzmann method. The pore-scale results were then linked with finite element simulation in a homogenized scheme to compute and upscale the transport properties to specimen scale. The simulated permeability and formation factor have magnitude and anisotropy that are in good agreement with experimental rock physics data. Together with simulated and measured values of con- nected porosity and specific surface area, they provide useful insights into how pore geometry controls the evolution of the transport properties.
    [Show full text]
  • Poromechanics
    Poromechanics Poromechanics Olivier Coussy Laboratoire Central des Ponts et Chauss´ees Institut Navier, Marne-la-Vall´ee, France Copyright c 2004 John Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex PO19 8SQ, England Telephone (+44) 1243 779777 Email (for orders and customer service enquiries): [email protected] Visit our Home Page on www.wileyeurope.com or www.wiley.com All Rights Reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic, mechanical, photocopying, recording, scanning or otherwise, except under the terms of the Copyright, Designs and Patents Act 1988 or under the terms of a licence issued by the Copyright Licensing Agency Ltd, 90 Tottenham Court Road, London W1T 4LP, UK, without the permission in writing of the Publisher. Requests to the Publisher should be addressed to the Permissions Department, John Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex PO19 8SQ, England, or emailed to [email protected], or faxed to (+44) 1243 770620. This publication is designed to provide accurate and authoritative information in regard to the subject matter covered. It is sold on the understanding that the Publisher is not engaged in rendering professional services. If professional advice or other expert assistance is required, the services of a competent professional should be sought. Other Wiley Editorial Offices John Wiley & Sons Inc., 111 River Street, Hoboken, NJ 07030, USA Jossey-Bass, 989 Market Street, San Francisco, CA 94103-1741, USA Wiley-VCH Verlag GmbH, Boschstr. 12, D-69469 Weinheim, Germany John Wiley & Sons Australia Ltd, 33 Park Road, Milton, Queensland 4064, Australia John Wiley & Sons (Asia) Pte Ltd, 2 Clementi Loop #02-01, Jin Xing Distripark, Singapore 129809 John Wiley & Sons Canada Ltd, 22 Worcester Road, Etobicoke, Ontario, Canada M9W 1L1 Wiley also publishes its books in a variety of electronic formats.
    [Show full text]