Splitting methods in the design of coupled flow and mechanics simulators
|
|
- Elijah Page
- 5 years ago
- Views:
Transcription
1 Splitting methods in the design of coupled flow and mechanics simulators Sílvia Barbeiro CMUC, Department of Mathematics, University of Coimbra PhD Program in Mathematics Coimbra, November 17, 2010 Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 1 / 36
2 Introduction Modeling porous media systems aims to describe the changes in a permeable media (made of somewhat rigid material containing open spaces) due to changes of the fluid it held. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 2 / 36
3 Introduction The changes in pressure caused by motion, removal or addition of liquid or gas may cause deformations in the structure holding the fluid. In the past, modelers tended to ignore the geomechanics in their calculations. Side effects of drilling: consolidation - reduction in volume due to fluid extraction compaction - reduction in volume due to air removal subsidence - collapse Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 3 / 36
4 Motivation Subsidence of Ekofisk Oil Field was a side effect of drilling. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 4 / 36
5 Motivation They subsided so much they had to go in and raise the platforms, costing them several billion dollars. If they d known ahead of time, they could have built their platforms taller Rick Dean, in Modeling complex, multiphase porous media systems, Siam News, April Photo: Norwegian Petroleum Museum/ConocoPhillips Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 5 / 36
6 Motivation Subsidence of Venice cased by the extraction of the aquifer. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 6 / 36
7 Modeling Poroelasticity refers to fluid flow within a deformable porous medium under the assumption of relatively small deformations. Examples of poroelastic structures soil, rock, cartilage, brain, heart, bone Various environmental, energy industry and biomechanics applications Subsidence, reservoir compaction Well stability, sand production Waste disposal Sequestration of carbon in saline aquifers Estimate tumor-induced stress levels in the brain Development of prosthetic devices for cartilage, bone, heart valves Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 7 / 36
8 Biot s consolidation model Karl von Terzaghi (October 2, October 25, 1963) Austrian civil engineer and geologist. Frequently called the father of soil mechanics. Maurice Anthony Biot (May 25, September 12, 1985) Belgian-American physicist Founder of the theory of poroelasticity. Karl von Terzaghi Maurice Anthony Biot Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 8 / 36
9 Biot s consolidation model primary variables: displacement u = (u 1, u 2, u 3 ) and fluid pressure p Balance of linear momentum (σ(u) αpi ) = f σ(u) - effective stress tensor, linear elastic, f - body force, α - Biot-Willis constant where σ(u) = 2µɛ(u) + λtr(ɛ(u))i, ɛ(u) = 1 2 ( grad u + (grad u) t), λ, µ - Lamé constants The momentum conservation equation is very similar to the equation governing linear elasticity, the exception is the addition of the term involving pressure. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 9 / 36
10 Biot s consolidation model Mass conservation η - fluid content, v f - flux of fluid, s f η t = v f + s f - volumetric fluid source term Equation for fluid content η in terms of fluid pressure p and material volume u η = c 0 p + α u c 0 - constrained specific storage coefficient Darcy s law for the flux of fluid v f = 1 µ f K( p ρ f g) K - permeability tensor, µ f - fluid viscosity, g - gravity Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 10 / 36
11 Summary of equations Coupling equations In the domain, at any time t Initial conditions at time t = 0 (σ(u) αpi ) = f t (c 0p + α u) 1 µ f K( p ρ f g) = s f Plus adequate boundary conditions p(0) = p 0, u(0) = u 0 For the mixed formulation for the flow, we introduce the variable z = 1 µ f K( p ρ f g). Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 11 / 36
12 Variational problem Integrating by parts over Ω, we obtain the variational problem Find u, p and z such that a u (u, v) α(p, v) = ( p ) ( ) c 0 t, w + α t u, w Ω f v + r N v, Γ N s f w, Ω p D s η Γ p + ( z, w) = (µ f K 1 z, s) (p, s) = ρ f g s holds for all (v, w, s) and t [0, T ], where a u (u, v) = (2µ(ɛ(u) : ɛ(v)) + λ( u)( v)) dx. Ω Ω Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 12 / 36
13 Time discretization Using the backward Euler method. Let n denote the time step, and t the time increment a u (u n+1, v) α(p n+1, v) = ( ) ( c 0 (p n+1 p n ), w + α (µ f K 1 z n+1, s) (p n+1, s) = holds for all (v, w, s). Ω f v + r N v, Γ N ) + t( z n+1, w) = t ρ f g s p D s η Ω Γ p (u n+1 u n ), w Ω s f w, Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 13 / 36
14 Time discretization Operator splitting α u = c r p + c r α σ FLOW ( ) (c 0 + c r )(p n+1,k+1 p n ), w + t( z n+1,k+1, w) = t s f w Ω ) α ( σn+1,k σ n ), w, ( cr (µ f K 1 z n+1,k+1, s) (p n+1,k+1, s) = p D s η + Γ p MECHANICS a u (u n+1,k+1, v) = Ω f v + r N v + α(p n+1,k+1, v) Γ N Ω ρ f g s Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 14 / 36
15 Splitting: iterative coupling Time step loop Time iteration n = n + 1 k = k + 1 FLOW σ n,k 1 MECHANICS No Converges? Yes Convergence criterion σ n+1,k σ n+1,k 1 L < Tol Iterative coupling is stable and accurate. [Wheeler and Gai, Numer. Meth. PDEs, 2007] Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 15 / 36
16 Fully, iteratively, explicit and loosely coupled Fully coupled The coupled governing equations of flow and geomechanics are solved simultaneously at every time step. Iteratively coupled Either the flow, or mechanical, problem is solved first, then the other problem is solved using the intermediate solution information. This sequential procedure is iterated at each time step until the solution converges to within an acceptable tolerance. The converged solution is identical to that obtained using the fully coupled approach. Explicitly coupled This is a special case of the iteratively coupled method, where only one iteration is taken. Loosely coupled The coupling between the two problems is resolved only after a certain number of flow time steps. This method can save computational cost compared to the other strategies, but it is less accurate and requires reliable estimates of when to update the mechanical response. [Kim, Tchelepi, Juanes, SPE, 2009] Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 16 / 36
17 Space discretization A very simple example u (x) + u(x) = (1 + π 2 ) sin(πx), x (0, 1), u(0) = 0, u(1) = 0 Variational formulation Multiplying the equation by any arbitrary weight function v H0 1 (0, 1) and integrating over the interval (0, 1) 1 0 u (x)v(x) dx + Integrating by parts 1 0 u(x)v(x) dx = 1 0 (1 + π 2 ) sin(πx)v(x) dx. 1 u (x)v (x) dx u(x)v(x) dx = 1 0 (1 + π 2 ) sin(πx)v(x) dx +u (1)v(1) u (0)v(0). Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 17 / 36
18 Finite element method The finite element method supplies an approximation to the analytical solution in the form of a piecewise polynomial function, defined over the entire computational domain. Example: The simplest case of linear splines. 1 ϕ i (x x i 1 )/h i, x i 1 x x i, ϕ i (x) = (x i+1 x)/h i+1, x i x x i+1, 0, otherwise. x i 1 x i x i+1 x A piecewise linear finite element basis function ϕ i (hat functions). Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 18 / 36
19 Finite element method x 0 = a x 1 x 2 x 3 x 4 x 5 x 6 x 7... A piecewise linear function. x n = b x Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 19 / 36
20 Finite element method Matrix form Since the aim of Finite Element Method method is the production of a linear system of equations, we build its matrix form, which can be used to compute the solution by a computer program. n We expand u n in respect to this basis, u n = c j ϕ j to obtain j=1 n c j a(ϕ j, ϕ i ) = f (ϕ i ) i = 1,..., n, j=1 which is a linear system of equations AU = F, where a ij = a(ϕ j, ϕ i ), F i = f (ϕ i ). Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 20 / 36
21 Finite element method For the finite element method the important property of the basis functions ϕ i, 1 i n is that they have local support, being nonzero only in one pair of adjacent intervals (x i 1, x i ] and [x i+1, x i ). This means that, A ij = 0 if i j > 0. = The matrix A is symmetric, positive definite and tridiagonal, and the associated system of linear equations can be solved very efficiently. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 21 / 36
22 Finite element method Back to the very simple example (Exact solution: u(x) = sin(πx)) Uniform subdivision of [0, 1] of spacing h = 1/n. n = 2 n = 4 Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 22 / 36
23 Finite element method Uniform subdivision of [0, 1] of spacing h = 1/n. n = 6 n = 100 In the last figure, the approximation error is so small that u and u n are indistinguishable. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 23 / 36
24 Finite element method Basis function in 2D Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 24 / 36
25 Finite element method Find u n V n such that v n V n, a(u n, v n ) = f (v n ). Galerkin orthogonality The key property of the Galerkin approach is that the error is orthogonal to the chosen subspaces. Since V n V, we can use v n as a test vector in the original equation. Subtracting the two, we get the Galerkin orthogonality relation for the error, e n = u u n and a(e n, v n ) = a(u, v n ) a(u n, v n ) = f (v n ) f (v n ) = 0 a(e n, e n ) = a(u u n, u u n ) = min v n V n a(u v n, u v n ). Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 25 / 36
26 Finite element method Galerkin orthogonality u ϕ H0 1 (a, b) u u n = (u ϕ) (u n ϕ) V n 0 u n ϕ Energy norm: v E = a(v, v) 1/2 u n is the best approximation from V n to the weak solution u H0 1 (a, b), when we measure the error of the approximation in the energy norm. Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 26 / 36
27 Spatial discretization E h and E H be two nondegenerate partitions of the polyhedral domain Ω with maximal element diameter h and H, respectively. Mixed spaces for flow variables Examples of mixed spaces with the needed properties are the Raviart-Thomas-Nedelec spaces. Example: Lowest order Raviart-Thomas 2D triangles: in each element p =const, z = (a + bx, c + by) t quadrilaterals: in each element p =const, z = (a + bx, c + dy) t Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 27 / 36
28 Mandel s problem Mandel s solution has been used as a benchmark problem for testing the validity of numerical codes of poroelasticity. [Mandel, 1953] - analytical solution for pressure [Abousleiman et al., 1996] - analytical solution for displacement and stress 2F y x 2F Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 28 / 36
29 Mandel s problem (λ + µ) ( u) µ 2 u + α p = 0 in Ω (0, T ] t (c 0p + α u) 1 µ f K p = 0 in Ω (0, T ] boundary conditions p = 0, x = a, 1 µ f K p η = 0, x = 0, y = 0, y = b, u 1 = 0, x = 0, u 2 = 0, y = 0, u y = 0, y = b, x ση = ( F /a)η, y = b, ση = 0, x = 0, x = a, y = 0, u x = 0 z x = 0 u y =0 z y =0 p = 0 computational domain Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 29 / 36
30 Mandel s problem surface: p, arrows: u 1 and u 2 T = Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 30 / 36
31 Mandel s problem surface: p, arrows: u 1 and u 2 T = 0.1 Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 31 / 36
32 Mandel s problem surface: p, arrows: u 1 and u 2 T = 0.2 Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 32 / 36
33 Mandel s problem surface: p, arrows: u 1 and u 2 T = 0.5 Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 33 / 36
34 Mandel s problem Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 34 / 36
35 Mandel s problem Numerical results H u ū H 1 rate p p L 2 rate z z L 2 rate 5.000e e e e e e e e e e e e e e e e e e e e e e e e-2 - Convergence rates: bilinear elements for u, lowest order Raviart-Thomas space for p and z Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 35 / 36
36 Important Questions Is the method stable? Is the method convergente? Next seminar... I really enjoy developing efficient and accurate solutions to real-world problems, while maintaining a solid theoretical base. MFW Sílvia Barbeiro (CMUC) Splitting methods for flow and mechanics simulators 36 / 36
Multiphysics modelling of the Mandel-Cryer effect
Int. Jnl. of Multiphysics Volume 10 Number 1 2016 11 Multiphysics modelling of the Mandel-Cryer effect E Holzbecher* German University of Technology in Oman (GUtech) Muscat, Oman ABSTRACT In porous medium
More informationA SIMULATOR WITH NUMERICAL UPSCALING FOR THE ANALYSIS OF COUPLED MULTIPHASE FLOW AND GEOMECHANICS IN HETEROGE-
A SIMULATOR WITH NUMERICAL UPSCALING FOR THE ANALYSIS OF COUPLED MULTIPHASE FLOW AND GEOMECHANICS IN HETEROGE- NEOUS AND DEFORMABLE POROUS AND FRACTURED MEDIA A Dissertation by DAEGIL YANG Submitted to
More informationIntro to Soil Mechanics: the what, why & how. José E. Andrade, Caltech
Intro to Soil Mechanics: the what, why & how José E. Andrade, Caltech The What? What is Soil Mechanics? erdbaumechanik The application of the laws of mechanics (physics) to soils as engineering materials
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationICES REPORT October Sanghyun Lee, Mary F. Wheeler, Thomas Wick
ICES REPORT 16-23 October 2016 Iterative coupling of flow, geomechanics and adaptive phase-field fracture including level-set crack width approaches by Sanghyun Lee, Mary F. Wheeler, Thomas Wick The Institute
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationGuaranteed and computable bounds of errors for semi-discrete approximations of the Biot problem
www.oeaw.ac.at Guaranteed and computable bounds of errors for semi-discrete approximations of the Biot problem K. Kumar, S. Matculevich, J.M. Nordbotten, S. Repin RICAM-Report 018-8 www.ricam.oeaw.ac.at
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationLab 5 - Aquifer Elasticity and Specific Storage. Due October 12 14, 2010
Lab 5 - Aquifer Elasticity and Specific Storage Due October 12 14, 2010 The goal of this experiment is to measure the specific storage S s of a balloon, which simulates aquifer elasticity. The experiment
More informationSEQUENTIAL METHODS FOR COUPLED GEOMECHANICS AND MULTIPHASE FLOW
SEQUENTIAL METHODS FOR COUPLED GEOMECHANICS AND MULTIPHASE FLOW A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN
More informationThe Effect of Stress Arching on the Permeability Sensitive Experiment in the Su Lige Gas Field
The Effect of Stress Arching on the Permeability Sensitive Experiment in the Su Lige Gas Field Fanliao Wang, Xiangfang Li, Gary Couples, Mingchuan Wang, Yiqun Zhang and Jingjing Zhao THE EFFECT OF STRESS
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationJ. Kim, SPE; H. A. Tchelepi, SPE, Stanford U.; and R. Juanes, SPE, Massachusetts Institute of Technology
SPE 9084 Stability, Accuracy and Efficiency of Sequential ethods for Coupled Flow and Geomechanics J. Kim, SPE; H. A. Tchelepi, SPE, Stanford U.; and R. Juanes, SPE, assachusetts Institute of Technology
More informationPre- and Post-Fracturing Analysis of a Hydraulically Stimulated Reservoir
ARMA 18 503 Pre- and Post-Fracturing Analysis of a Hydraulically Stimulated Reservoir Erfan Sarvaramini Department of Civil and environmental Engineering, University of Waterloo, Waterloo, ON, Canada Robert
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationINL Capabilities and Approach to CO 2 Sequestration. 4 th U.S.-China CO2 Emissions Control Science & Technology Symposium
www.inl.gov INL Capabilities and Approach to CO 2 Sequestration 4 th U.S.-China CO2 Emissions Control Science & Technology Symposium Travis McLing, Ph.D. Energy Resources Recovery and Sustainability September
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationHow are calculus, alchemy and forging coins related?
BMOLE 452-689 Transport Chapter 8. Transport in Porous Media Text Book: Transport Phenomena in Biological Systems Authors: Truskey, Yuan, Katz Focus on what is presented in class and problems Dr. Corey
More informationChapter 1 INTRODUCTION
CONSTITUTIVE MODELS AND CHALK Chapter 1 INTRODUCTION Computing power has advanced significantly in the years since soil mechanics and rock mechanics first became mature fields. Probably the single most
More informationReservoir Simulator Compaction Modelling: A Predictor for Accelerated Coupled Rock Mechanics -- Reservoir Simulation
Reservoir Simulator Compaction Modelling: A Predictor for Accelerated Coupled Rock Mechanics -- Reservoir Simulation by Øystein Pettersen Centre for Integrated Petroleum Research, Bergen, Norway ECMOR
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationDefmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationA NOVEL FULLY COUPLED GEOMECHANICAL MODEL FOR CO 2 SEQUESTRATION IN FRACTURED AND POROUS BRINE AQUIFERS
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champagne June 17-22, 2012 A NOVEL FULLY COUPLED GEOMECHANICAL MODEL FOR CO 2 SEQUESTRATION IN FRACTURED AND POROUS
More informationSolving an elasto-plastic model using DOLFIN
Solving an elasto-plastic model using DOLFIN TTI 2005 Johan Jansson johanjan@math.chalmers.se Chalmers University of Technology Solving an elasto-plastic model using DOLFIN p. Overview Motivation Previous
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationPractical Geomechanics
www.bakerhughes.com Practical Geomechanics Baker Hughes - RDS Geomechanics Services 2015 Baker Hughes Incorporated. All rights reserved Copyright By accepting these materials you agree that all materials
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationModeling two-phase flow in strongly heterogeneous porous media
Presented at the COMSOL Conference 2010 China COMSOL 2010 I«^rc Modeling two-phase flow in strongly heterogeneous porous media Zhaoqin Huang Research Center for Oil & Gas Flow in Reservoir, School of Petroleum
More informationSHAPED CHARGE PENETRATION INTO STRESSED ROCK
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 SHAPED CHARGE PENETRATION INTO STRESSED ROCK B. Grove 1, J. Heiland 1, and I. Walton 1 1 Schlumberger Reservoir Completions
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationPhase-field modeling of a fluid-driven fracture in a poroelastic medium
Phase-field modeling of a fluid-driven fracture in a poroelastic medium Université Lyon 1, Institut Camille Jordan, FRANCE joint work with M.F. Wheeler (UT Austin), T. Wick (Ecole Polytechnique) Multiphysics,
More informationFinite Element Methods in Soil Mechanics
Aalto University Department of Mathematics and Systems Analysis Graduate School in Engineering Mechanics Annual Seminar 28.-29.1.2010 Background Master's Thesis late 2007 Topic: composite Reissner-Mindlin
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationA phase-field method for propagating fluid-filled fractures coupled to a surrounding porous medium
Accepted for publication in SIAM Multiscale Modeling and Simulation, 215 A phase-field method for propagating fluid-filled fractures coupled to a surrounding porous medium Andro Mikelić Mary F. Wheeler
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationBiological Transportation Networks PDE Modeling and Numerics
MÜNSTER Biological Transportation Networks PDE Modeling and Numerics joint work with Martin Burger () Jan Haskovec (KAUST) Peter Markowich (KAUST/Cambridge) Benoît Perthame (LLJLUPMC) Data-Rich Phenomena
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationA New Method for Calculating Oil-Water Relative Permeabilities with Consideration of Capillary Pressure
A Ne Method for Calculating Oil-Water Relative Permeabilities ith Consideration of Capillary Pressure K. Li, P. Shen, & T. Qing Research Institute of Petroleum Exploration and Development (RIPED), P.O.B.
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationICES REPORT A Phase-Field Method For Propagating Fluid-Filled Fractures Coupled To A Surrounding Porous Medium
ICES REPORT 14-8 April 214 A Phase-Field Method For Propagating Fluid-Filled Fractures Coupled To A Surrounding Porous Medium by Andro Mikelic, Mary F. Wheeler, Thomas Wick The Institute for Computational
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationVISCOUS FLUID FLOW IN AN INCLINED CHANNEL WITH DEFORMABLE POROUS MEDIUM
International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue, January 08, pp. 970 979 Article ID: IJMET_09_0_04 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationA Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach
SPE-172990-MS A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach Benjamin Ganis, Kundan Kumar, and Gergina Pencheva, SPE; Mary F. Wheeler, The
More informationractical Geomechanics for Unconventional Resources
P ractical Geomechanics for Unconventional Resources 24-26 October 2012, Calgary, Canada Practical Geomechanics for Unconventional Resources Nowadays, unconventional resources have been brought into the
More informationWELL-POSEDNESS OF A FULLY COUPLED THERMO-CHEMO-POROELASTIC SYSTEM WITH APPLICATIONS TO PETROLEUM ROCK MECHANICS
Electronic Journal of Differential Equations, Vol. 017 017), No. 137, pp. 1. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu WELL-POSEDNESS OF A FULLY COUPLED THERMO-CHEMO-POROELASTIC
More informationOn a variational inequality of Bingham and Navier-Stokes type in three dimension
PDEs for multiphase ADvanced MATerials Palazzone, Cortona (Arezzo), Italy, September 17-21, 2012 On a variational inequality of Bingham and Navier-Stokes type in three dimension Takeshi FUKAO Kyoto University
More informationA High Order Method for Three Phase Flow in Homogeneous Porous Media
A High Order Method for Three Phase Flow in Homogeneous Porous Media Justin Dong Advisor: B eatrice Rivi`ere February 5, 2014 Abstract The modeling of three-phase fluid flow has many important applications
More informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationDIMENSIONAL MODEL REDUCTION FOR FLOW THROUGH FRACTURES IN POROELASTIC MEDIA
MOX-Report No. 45/6 DIMENSIONAL MODEL REDUCTION FOR FLOW THROUGH FRACTURES IN POROELASTIC MEDIA Bukac, M.; Yotov, I.; Zunino, P. MOX, Dipartimento di Matematica Politecnico di Milano, Via Bonardi 9-33
More informationICES REPORT January Saumik Dana and Mary. F. Wheeler
ICES REPORT 18-01 January 2018 Convergence analysis o ixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter by Saumik Dana and Mary. F. Wheeler The Institute or Computational
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationINTEGRATION OF RESERVOIR SIMULATION AND GEOMECHANICS
INTEGRATION OF RESERVOIR SIMULATION AND GEOMECHANICS by Nan Zhao A dissertation submitted to the faculty of The University of Utah in partial fulfillment of the requirements for the degree of Doctor of
More informationClick to add title. Continuum mechanics of two-phase porous media
Click to add title Continuum mechanics of two-phase porous media Ragnar Larsson Division of Material and Computational Mechanics Department of Applied Mechanics Chalmers University of Technology S-412
More informationUniversitat Politècnica de Catalunya BARCELONATECH Escola Tècnica Superior d Enginyers de Camins, Canals i Ports. Soil Mechanics.
Universitat Politècnica de Catalunya BARCELONATECH Escola Tècnica Superior d Enginyers de Camins, Canals i Ports Soil Mechanics Chapter 8 Consolidation Chapter 6 1. One-dimensional consolidation theory.
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationThis paper was prepared for presentation at the Unconventional Resources Technology Conference held in Denver, Colorado, USA, August 2013.
1 SPE 168873 / URTeC 1606914 Numerical Upcaling of Coupled Flow and Geomechanics in Highly Heterogeneous Porous Media Daegil Yang, Texas A&M University; George J. Moridis, Lawrence Berkeley National Laboratory;
More informationUpscaling of coupled geomechanics, flow, and heat in a poro-elastic medium in the quasi-static situation
Upscaling of coupled geomechanics, flow, and heat in a poro-elastic medium in the quasi-static situation Mats K. Brun, Florin A. Radu, Inga Berre, J. M. Nordbotten November 14, 2017 Abstract Motivated
More informationModeling Multiphase Flow in Porous Media with Complementary Constraints
Background Modeling Multiphase Flow in Porous Media with Complementary Constraints Applied Math, Statistics, and Scientific Computation, University of Maryland - College Park October 07, 2014 Advisors
More informationDerivation of the fractional flow equation for a one-dimensional oil-water system. Consider displacement of oil by water in a system of dip angle α
TPG45 Reservoir Recovery Techniques 27 /9 BUCKLEY-LEVERETT ANALYSIS Derivation of the fractional flow equation for a one-dimensional oil-water system Consider displacement of oil by water in a system of
More informationComputational Techniques for Efficient Solution of Discretized Biot s Theory for Fluid Flow in Deformable Porous Media. Im Soo Lee
Computational Techniques for Efficient Solution of Discretized Biot s Theory for Fluid Flow in Deformable Porous Media by Im Soo Lee Dissertation submitted to the Faculty of the Virginia Polytechnic Institute
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationICES REPORT February Saumik Dana and Mary F. Wheeler
ICE REPORT 18-2 February 218 A framework for arriving at bounds on effective moduli in heterogeneous poroelastic solids by aumik Dana and Mary F. Wheeler The Institute for Computational Engineering and
More informationractical Geomechanics for Oil & Gas Industry
P ractical Geomechanics for Oil & Gas Industry Practical Geomechanics for Oil and Gas Industry The integrity of the wellbore plays an important role in petroleum operations including drilling, completion
More informationPORE PRESSURE EVOLUTION AND CORE DAMAGE: A COMPUTATIONAL FLUID DYNAMICS APPROACH
SCA211-41 1/6 PORE PRESSURE EVOLUTION AND CORE DAMAGE: A COMPUTATIONAL FLUID DYNAMICS APPROACH I. Zubizarreta, M. Byrne, M.A. Jimenez, E. Roas, Y. Sorrentino and M.A. Velazco. Senergy. Aberdeen, United
More information