ICES REPORT January Saumik Dana and Mary. F. Wheeler
|
|
- Kristopher Byrd
- 5 years ago
- Views:
Transcription
1 ICES REPORT 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 Engineering and Sciences The University o Texas at Austin Austin, Texas Reerence: Saumik Dana and Mary. F. Wheeler, "Convergence analysis o ixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter," ICES REPORT 18-01, The Institute or Computational Engineering and Sciences, The University o Texas at Austin, January 2018.
2 Convergence analysis o ixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter Saumik Dana a, Mary. F. Wheeler a a Center or Subsurace Modeling, Institute or Computational Engineering and Sciences, UT Austin, Austin, TX Abstract We perorm a convergence analysis o the ixed stress split iterative algorithm or the Biot system modeling coupled low and deormation in anisotropic poroelastic media with ull Biot tensor. The ixed stress split iterative scheme solves the low subproblem with all components o the stress tensor rozen using a multipoint lux mixed inite element method, ollowed by the poromechanics subproblem using a conorming Galerkin method in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed or time marching. The convergence analysis is based on studying the equations satisied by the dierence o iterates to show that the ixed stress split iterative scheme or anisotropic poroelasticity with Biot tensor is contractive. Keywords: Anisotropic poroelastic medium, Biot tensor, Fixed stress split iterative coupling, Contractivity 1. Introduction Staggered solution algorithms are those in which operator splitting strategies are used to decompose coupled problems into subproblems which are then solved sequentially in successive iterations until a convergence criterion is met at each time step (Felippa et al. [13], Armero and Simo [3], Turska and Schreler [25], Schreler et al. [21]). The ixed stress split is one such operator splitting strategy or isotropic poroelasticity with scalar Biot parameter in which the low problem is solved irst while reezing the volumetric mean stress o the porous solid. Such an approach has been shown to enjoy numerical and theoretical convergence addresses: saumik@utexas.edu (Saumik Dana), mw@ices.utexas.edu (Mary. F. Wheeler)
3 New time step New coupling iteration Solve low with stress tensor ixed Solve poromechanics yes converged? no Figure 1: Fixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter (Kim et al. [18], Mikelić and Wheeler [19], Castelletto et al. [10], Almani et al. [2]). The subject o this paper is the extension o the decoupling assumption to the case o anisotropic poroelasticity with tensor Biot parameter. The micromechanical analyses or the case o anisotropic poroelasticity (see Carroll [8], Carroll and Katsube [9], Katsube [17]) revealed that the modiication to the stress applied to the porous solid due to the presence o pore luid pressure is not hydrostatic, as it is in the case o isotropic poroelasticity (see Nur and Byerlee [20]). Further, unlike in case o isotropic poroelasticity where the solid-luid coupling parameter is a scalar (see Biot [4], Geertsma [14], Skempton [23], Nur and Byerlee [20]), the coupling parameter in case o anisotropic poroelasticity is a tensor (see Biot [5], Coussy [12]). Intuitively, this implies that the decoupling assumption o reezing the hydrostatic part o the stress tensor during the low solve in every coupling iteration or the case o isotropic poroelasticity requires a generalization or the case o anisotropic poroelasticity. As we shall show in the ensuing convergence analysis, the decoupling assumption o reezing all components o the stress tensor during the low solve in every coupling iteration does enjoy theoretical convergence or the case o anisotropic poroelasticity with tensor Biot parameter. The convergence analysis is motivated by the previous work o Mikelić and Wheeler [19] and Almani et al. [2]. Mikelić and Wheeler [19] proved that the standard ixed stress split scheme or isotropic poroelasticity with scalar Biot parameter is a contraction 2
4 map with respect to appropriately chosen metrics. Almani et al. [2] extended those results to establish convergence o the multirate ixed stress split scheme in which the subproblems are resolved on dierent time scales with the poromechanics subproblem being resolved on a larger time scale. To the best o our knowledge, this is the irst time a contraction map has been rigorously derived or the ixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter. This paper is structured as ollows: Section 2 presents the model equations or low and poromechanics. Section 3 presents the statement o contraction o the ixed stress split iterative scheme or anisotropic poroelasticity with tensor Biot parameter. Section 4 presents the convergence criterion or the algorithm Preliminaries Given a bounded convex domain Ω R 3, we use P k (Ω) to represent the restriction o the space o polynomials o degree less that or equal to k to Ω and Q 1 (Ω) to denote the space o trilinears on Ω. For the sake o convenience, we discard the dierential in the integration o any scalar ield χ over Ω as ollows χ(x) χ(x) dv ( x Ω) Ω Ω Sobolev spaces are based on the space o square integrable unctions on Ω given by L 2 (Ω) { θ : θ 2 Ω := θ 2 < + }, The inner product o two second order tensors S and T is given by (see Gurtin et al. [15]) S : T = S ij T ij (i, j = 1, 2, 3) A ourth order tensor is a linear transormation o a second order tensor to a second order tensor in the ollowing manner (see Gurtin et al. [15]) PS = T P ijkl S kl = T ij (i, j, k, l = 1, 2, 3) Ω 2. Model equations 2.1. Flow model Let the boundary Ω = Γ D Γ N where Γ D is the Dirichlet boundary and Γ N is the Neumann boundary. The luid mass conservation equation (2.1) in the presence o deormable 3
5 and anisotropic porous medium with the Darcy law (2.2) and linear pressure dependence o density (2.3) with boundary conditions (2.4) and initial conditions (2.5) is ζ t + z = q (2.1) z = K µ ( p ρ 0g) = κ( p ρ 0 g) (2.2) ρ = ρ 0 (1 + c (p p 0 )) (2.3) p = g on Γ D (0, T ], z n = 0 on Γ N (0, T ] (2.4) p(x, 0) = p 0 (x), ρ(x, 0) = ρ 0 (x), φ(x, 0) = φ 0 (x) ( x Ω) (2.5) where p : Ω (0, T ] R is the luid pressure, z : Ω (0, T ] R 3 is the luid lux, ζ is the luid content, n is the unit outward normal on Γ N, q is the source or sink term, K is the uniormly symmetric positive deinite absolute permeability tensor, µ is the luid viscosity, ρ 0 is a reerence density, φ is the porosity, κ = K µ is a measure o the hydraulic conductivity o the pore luid, c is the luid compressibility and T > 0 is the time interval Poromechanics model Let the boundary Ω = Γ p D Γp N where Γp D is the Dirichlet boundary and Γp N is the Neumann boundary. Linear momentum balance or the anisotropic porous solid in the quasistatic limit o interest (2.6) with small strain assumption (2.8) with boundary conditions (2.9) and initial condition (2.10) is σ + = 0 (2.6) = ρφg + ρ r (1 φ)g (2.7) ɛ(u) = 1 2 ( u + ( u)t ) (2.8) u n 1 = 0 on Γ p D [0, T ], σt n 2 = t on Γ p N [0, T ] (2.9) u(x, 0) = 0 x Ω (2.10) where u : Ω [0, T ] R 3 is the solid displacement, ρ r is the rock density, is the body orce per unit volume, n 1 is the unit outward normal to Γ p D, n 2 is the unit outward normal to Γ p N, t is the traction speciied on Γp N, ɛ is the strain tensor, σ is the Cauchy stress tensor given by the generalized Hooke s law (see Coussy [12]) σ = σ 0 + Mɛ α(p p 0 ) (2.11) 4
6 where σ 0 is the in situ stress, M is the ourth order anisotropic elasticity tensor and α is the Biot tensor. The inverse o the generalized Hooke s law (2.11) is given by (see Cheng [11]) ɛ = Cσ + 1 CBp (2.12) 3 where C is the ourth order anisotropic compliance tensor, C(> 0) is a generalized Hooke s law constant and B is a generalization o the Skempton pore pressure coeicient B (see Skempton [22]) or anisotropic poroelasticity Fluid content The luid content ζ is given by (see Cheng [11]) ζ = Cp CB : σ 1 M p + α : ɛ (2.13) where M(> 0) is a generalization o the Biot modulus (see Biot and Willis [6]) or anisotropic poroelasticity. 3. Statement o contraction o the ixed stress split scheme or anisotropic poroelasticity with Biot tensor We use the notations ( ) n+1 or any quantity ( ) evaluated at time level n + 1, ( ) m,n+1 or any quantity ( ) evaluated at the m th coupling iteration at time level n + 1, δ (m) ( ) or the change in the quantity ( ) during the low solve over the (m + 1) th coupling iteration at any time level and δ (m) ( ) or the change in the quantity ( ) over the (m + 1) th coupling iteration at any time level. Let T h be inite element partition o Ω consisting o distorted hexahedral elements E where h = max E T h diam(e). The details o the inite element mapping are given in Appendix A. The discrete variational statements in terms o coupling iteration dierences is : ind δ (m) p h W h, δ (m) z h V h and δ (m) u h U h such that C(δ (m) p h, θ h ) Ω + t( δ (m) z h, θ h ) Ω = C 3 (B : δ(m 1) σ, θ h ) Ω (3.1) (κ 1 δ (m) z h, v h ) Ω = (δ (m) p h, v h ) Ω (3.2) (δ (m) σ : ɛ(q h )) Ω = 0 (3.3) 5
7 where the inite dimensional spaces W h, V h and U h are W h = { } θ h : θ h P 0 (E) E T h V h = { v h : v h E ˆv Ê ˆV(Ê) E T h, v h n = 0 on Γ } N U h = { q h = (u, v, w) E Q 1 (E) E T h, q h = 0 on Γ p } D and the details o ˆV(Ê) are given in Appendix B. The equations (3.1), (3.2) and (3.3) are the discrete variational statements (in terms o coupling iteration dierences) o (2.1), (2.2) and (2.6) respectively. The details o (3.1) and (3.2) are given in Appendix C whereas the details o (3.3) are given in Appendix D. Theorem 3.1. The ixed stress split iterative coupling scheme or anisotropic poroelasticity with Biot tensor in which the low problem is solved irst by reezing all components o the stress tensor is a contraction given by B : δ (m) σ 2 Ω + 2 δ (m) p h 2 Ω + 6 t C κ 1/2 δ (m) z h 2 Ω + 6 C 2 δ(m) ζ δ (m) ζ 2 Ω B : δ (m 1) σ 2 Ω 0 6 C (δ(m) σ : Cδ (m) σ) Ω + 3 C 2 δ(m) ζ 2 Ω where the quantity δ (m) ζ δ (m) ζ 2 Ω is driven to a small value via the convergence criterion or the algorithm. Proo. Step 1: Flow equations Testing (3.1) with θ h δ (m) p h, we get C δ (m) p h 2 Ω + t( δ (m) z h, δ (m) p h ) Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.4) Testing (3.2) with v h δ (m) z h, we get From (3.4) and (3.5), we get κ 1/2 δ (m) z h 2 Ω = (δ (m) p h, δ (m) z h ) Ω (3.5) C δ (m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.6) Step 2: Poromechanics equations Testing (3.3) with q h δ (m) u h, we get (δ (m) σ : δ (m) ɛ) Ω = 0 (3.7) 6
8 We now invoke (2.12) to arrive at the expression or change in strain tensor over the (m+1) th coupling iteration as ollows δ (m) ɛ = Cδ (m) σ + C 3 Bδ(m) p h (3.8) Substituting (3.8) in (3.7), we get (δ (m) σ : Cδ (m) σ) Ω + C 3 (B : δ(m) σ, δ (m) p h ) Ω = 0 (3.9) Step 3: Combining low and poromechanics equations Adding (3.6) and (3.9), we get C δ (m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + C 3 (B : δ(m) σ, δ (m) p h ) Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.10) Step 4: Variation in luid content In lieu o (2.13), the variation in luid content in the (m + 1) th coupling iteration can be written as As a result, we can write δ (m) ζ = Cδ (m) p h + C 3 B : δ(m) σ (3.11) 1 2C δ(m) ζ 2 Ω C 2 δ(m) p h 2 Ω C 18 B : δ(m) σ 2 Ω = C 3 (B : δ(m) σ, δ (m) p h ) Ω (3.12) From (3.10) and (3.12), we get C δ (m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + 1 2C δ(m) ζ 2 Ω C 2 δ(m) p h 2 Ω C 18 B : δ(m) σ 2 Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.13) Adding and subtracting C 6 B : δ(m) σ 2 Ω to the LHS o (3.13) results in C 6 B : δ(m) σ 2 Ω + C 2 δ(m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + 1 2C δ(m) ζ 2 Ω C 9 B : δ(m) σ 2 Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.14) In lieu o (2.13) and the ixed stress constraint during the low solve, the variation in luid content during the low solve in the (m + 1) th coupling iteration is given by δ (m) ζ = Cδ (m) p h + C 3 B : δ (m) σ 0 7
9 Further, since the pore pressure is rozen during the poromechanical solve, we have δ (m) p h = δ (m) p h. As a result, we can write Subtracting (3.15) rom (3.11), we can write which implies that In lieu o (3.16), we can write (3.14) as δ (m) ζ = Cδ (m) p h (3.15) δ (m) ζ δ (m) ζ = C 3 B : δ(m) σ 1 C δ(m) ζ δ (m) ζ 2 Ω = C 9 B : δ(m) σ 2 Ω (3.16) C 6 B : δ(m) σ 2 Ω + C 2 δ(m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + 1 2C δ(m) ζ 2 Ω 1 C δ(m) ζ δ (m) ζ 2 Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.17) Step 5: Invoking positive-semideiniteness o the compliance tensor The ourth order anisotropic compliance tensor C is positive-semideinite (see Gurtin et al. [15]) in the sense that it obeys S : CS 0 or all symmetric tensors S. In lieu o the above, since the Cauchy stress tensor is symmetric, the ollowing holds true In lieu o (3.18), we can write (3.17) as δ (m) σ : Cδ (m) σ 0 (3.18) 0 C 6 B : δ(m) σ 2 C Ω + 2 δ(m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + 1 2C δ(m) ζ 2 Ω 1 C δ(m) ζ δ (m) ζ 2 Ω = C 3 (B : δ(m 1) σ, δ (m) p h ) Ω (3.19) The quantity δ (m) ζ δ (m) ζ 2 Ω on the LHS o (3.19) is driven to a small value via the convergence criterion or the algorithm as explained in Section 4. Step 6: Invoking the Young s inequality Since the sum o the terms on the LHS o (3.19) is nonnegative, the RHS is also nonnegative. We invoke the Young s inequality (see Steele [24]) ab a2 2 + b2 2 ( a, b R) 8
10 or the RHS o (3.19) as ollows C 3 (B : δ(m 1) σ, δ (m) p h ) Ω C 3 ( 1 2 B : δ(m 1) σ 2 Ω + 1 ) 2 δ(m) p h 2 Ω (3.20) In lieu o (3.20), we write (3.19) as 0 C 6 B : δ(m) σ 2 C Ω + 2 δ(m) p h 2 Ω + t κ 1/2 δ (m) z h 2 Ω + (δ (m) σ : Cδ (m) σ) Ω + 1 C δ(m) ζ δ (m) ζ 2 Ω C 6 B : δ(m 1) σ 2 Ω + C 6 δ(m) p h 2 Ω which can also be written as B : δ (m) σ 2 Ω + 2 δ (m) p h 2 Ω C δ(m) ζ 2 Ω 6 t C κ 1/2 δ (m) z h 2 6 Ω + C (δ(m) σ : Cδ (m) 3 σ) Ω + C 2 δ(m) ζ 2 Ω 6 C 2 δ(m) ζ δ (m) ζ 2 Ω B : δ (m 1) σ 2 Ω (3.21) 4. Convergence criterion or the ixed stress split algorithm or anisotropic poroelasticity with Biot tensor In lieu o (2.13), the variation in luid content δ (m) ζ measured during the low solve in the (m + 1) th coupling iteration at any time step is δ (m) ζ = 1 M δ(m) p h + α : δ (m) ɛ = ζp m+1 ζ m (4.1) where ζ m is the luid content at the end o the previous or m th coupling iteration and ζ m+1 p serves as the predictor to the luid content at the end o the current or (m + 1) th coupling iteration. Similarly, the variation in luid content δ (m) ζ over the (m+1) th coupling iteration (including the low solve and poromechanics solve) at any time step is δ (m) ζ = 1 M δ(m) p h + α : δ (m) ɛ = ζ m+1 c ζ m (4.2) where ζc m+1 serves as the corrector to ζp m+1. Subtracting (4.1) rom (4.2), we get δ (m) ζ δ (m) ζ ζc m+1 ζp m+1 = α : (δ (m) ɛ δ (m) ɛ) = α : (ɛ m+1 c ɛ m+1 p ) 9
11 which implies that the dierence between the predicted value and the corrected value o the luid content at the end o the (m + 1) th coupling iteration is equal to the dierence between the predicted value (ɛ m+1 p ) and the corrected value (ɛ m+1 c ) o the strain tensor at the end o the (m + 1) th coupling iteration scaled by the Biot tensor α. In lieu o the above, the stopping criterion or coupling iterations at any time step is δ (m) ζ δ (m) ζ m+1 c where T OL > 0 is a pre-speciied tolerance. ζ ζ c m+1 ζ m+1 p T OL L L ζ m+1 c Appendix A. Finite element mapping ˆr 8 ˆr 7 r 8 r 7 ˆr 5 ˆr 6 ˆr 4 ˆr 3 F E r r 6 5 r 4 r 3 ˆr 1 ˆr 2 r 1 r 2 Figure A.2: Trilinear mapping F E : Ê E or 8 noded distorted hexahedral elements. The aces o E can be non-planar. Let r i, i = 1,.., 8 be the vertices o E. Now consider a reerence cube Ê with vertices ˆr 1 = [0 0 0] T, ˆr 2 = [1 0 0] T, ˆr 3 = [1 1 0] T, ˆr 4 = [0 1 0] T, ˆr 5 = [0 0 1] T, ˆr 6 = [1 0 1] T, ˆr 7 = [1 1 1] T and ˆr 8 = [0 1 1] T as shown in Figure A.2. Let ˆx = (ˆx, ŷ, ẑ) Ê and x = (x, y, z) E. The unction F E (ˆx) : Ê E is F E (ˆx) = r 1 (1 ˆx)(1 ŷ)(1 ẑ) + r 2ˆx(1 ŷ)(1 ẑ) + r 3ˆxŷ(1 ẑ) + r 4 (1 ˆx)ŷ(1 ẑ) + r 5 (1 ˆx)(1 ŷ)ẑ + r 6ˆx(1 ŷ)ẑ + r 7ˆxŷẑ + r 8 (1 ˆx)ŷẑ Denote Jacobian matrix by DF E and let J E = det(df E ). Deining r ij r i r j, we have DF E (ˆx)= r 21 + (r 34 r 21 )ŷ + (r 65 r 21 )ẑ + ((r 21 r 34 ) (r 65 r 78 ))ŷẑ; r 41 + (r 34 r 21 )ˆx + (r 85 r 41 )ẑ + ((r 21 r 34 ) (r 65 r 78 ))ˆxẑ; r 51 + (r 65 r 21 )ˆx + (r 85 r 41 )ŷ + ((r 21 r 34 ) (r 65 r 78 ))ˆxŷ 10
12 Denote inverse mapping by F 1 1, its Jacobian matrix by DF such that E E and let J F 1 E = det(df 1 E ) DF 1 E (x) = (DF E) 1 (ˆx); J F 1(x) = (J E ) 1 (ˆx) Let φ(x) be any unction deined on E and ˆφ(ˆx) be its corresponding deinition on Ê. Then we have E φ = (DF 1 E )T (x) ˆ ˆφ = (DF E ) T (ˆx) ˆ ˆφ (A.1) Appendix B. Enhanced BDDF 1 spaces ˆv 51 ˆv 52 ˆv 11 ˆv 53 ˆv 83 ˆv 81 ˆv 82 ˆv 41 ˆv 42 ˆv 43 ˆv 62 ˆv 63 ˆv 73 ˆv61 ˆv 21 ˆv ˆv ˆv 22 ˆv 23 ˆv 72 ˆv 32 ˆv 31 ˆv 33 ˆv 71 F E v 51 v 52 v 41 v 43 v 83 v 81 v 82 v 53 v 61 v 62 v 42 v 73 v 63 v 71 v 11 v 12 v 13 v22 v 23 v 21 v 72 v 32 v 33 v 31 Figure B.3: Degrees o reedom and basis unctions or the enhanced BDDF 1 velocity space on hexahedra. For the sake o clarity, we provide a brie description o the mixed inite element spaces used in the low model. Let V h W h be the lowest order BDDF 1 MFE spaces on hexahedra (see Brezzi et al. [7]). On the reerence unit cube these spaces are deined as ˆV (Ê) = P 1(Ê) + r 0 curl(0, 0, ˆxŷẑ) T + r 1 curl(0, 0, ˆxŷ 2 ) T + s 0 curl(ˆxŷẑ, 0, 0) T + s 1 curl(ŷẑ 2, 0, 0) T + t 0 curl(0, ˆxŷẑ, 0) T + t 1 curl(0, ˆx 2 ẑ, 0) T Ŵ (Ê) = P0(Ê) with the ollowing properties ˆ ˆV (Ê) = Ŵ (Ê), and ˆv ˆV (Ê), ê Ê, ˆv ˆn ê P 1 (ê) The multipoint lux approximation procedure requires on each ace one velocity degree o reedom to be associated with each vertex. Since the BDDF 1 space V h has only three 11
13 degrees o reedom per ace, it is augmented with six more degrees o reedom (one extra degree o reedom per ace). Since the constant divergence, the linear independence o the shape unctions and the continuity o the normal component across the element aces are to be preserved, six curl terms are added (Ingram et al. [16]). Let V h W h be the enhanced BDDF 1 spaces on hexahedra. On the reerence unit cube these spaces are ˆV(Ê) = ˆV (Ê) + r 2 curl(0, 0, ˆx 2 ẑ) T + r 3 curl(0, 0, ˆx 2 ŷẑ) T + s 2 curl(ˆxŷ 2, 0, 0) T + s 3 curl(ˆxŷ 2 ẑ 2, 0, 0) T + t 2 curl(0, ŷẑ 2, 0) T + t 3 curl(0, ˆxŷẑ 2, 0) T Ŵ (Ê) = P0(Ê) with the ollowing properties ˆ ˆV(Ê) = Ŵ (Ê), and ˆv ˆV(Ê), ê Ê, ˆv ˆn ê Q 1 (ê) where Q 1 is the space o bilinear unctions. Since dim Q 1 (ê) = 4, the dimension o ˆV(Ê) is 24 as shown in Figure B.3. Appendix C. Discrete variational statements or the low subproblem in terms o coupling iteration dierences Beore arriving at the discrete variational statement o the low model, we impose the ixed stress constraint on the strong orm o the mass conservation equation (2.1). In lieu o (2.13), we write (2.1) as t (Cp + C 3 B : σ) + z = q C p t + z = q C 3 B : σ t (C.1) Using backward Euler in time, the discrete in time orm o (C.1) or the m th coupling iteration in the (n + 1) th time step is written as C 1 t (pm,n+1 p n ) + z m,n+1 = q n+1 1 C t 3 B : (σm,n+1 σ n ) where t is the time step and the source term as well as the terms evaluated at the previous time level n do not depend on the coupling iteration count as they are known quantities. The ixed stress constraint implies that σ m,n+1 gets replaced by σ m 1,n+1 i.e. the computation 12
14 o p m,n+1 and z m,n+1 is based on the value o stress updated ater the poromechanics solve rom the previous coupling iteration m 1 at the current time level n + 1. The modiied equation is written as C(p m,n+1 p n ) + t z m,n+1 = tq n+1 C 3 B : (σm,n+1 σ n ) As a result, the discrete weak orm o (2.1) is given by C(p m,n+1 h p n h, θ h) Ω + t( z m,n+1 h, θ h ) Ω = t(q n+1, θ h ) Ω C 3 (B : (σm 1,n+1 σ n ), θ h ) Ω Replacing m by m + 1 and subtracting the two equations, we get C(δ (m) p h, θ h ) Ω + t( δ (m) z h, θ h ) Ω = C 3 (B : δ(m 1) σ, θ h ) Ω The weak orm o the Darcy law (2.2) or the m th coupling iteration in the (n + 1) th time step is given by (κ 1 z m,n+1, v) Ω = ( p m,n+1, v) Ω + (ρ 0 g, v) Ω v V(Ω) (C.2) where V(Ω) is given by V(Ω) H(div, Ω) { v : v n = 0 on Γ } N and H(div, Ω) is given by H(div, Ω) { v : v (L 2 (Ω)) 3, v L 2 (Ω) } We use the divergence theorem to evaluate the irst term on RHS o (C.2) as ollows ( p m,n+1, v) Ω = (, p m,n+1 v) Ω (p m,n+1, v) Ω = (p m,n+1, v n) Ω (p m,n+1, v) Ω = (g, v n) Γ (p m,n+1, v) Ω D (C.3) where we invoke v n = 0 on Γ N. In lieu o (C.2) and (C.3), we get (κ 1 z m,n+1, v) Ω = (g, v n) Γ + (p m,n+1, v) Ω + (ρ 0 g, v) Ω D Replacing m by m + 1 and subtracting the two equations, we get (κ 1 δ (m) z h, v h ) Ω = (δ (m) p h, v h ) Ω 13
15 Appendix D. Discrete variational statement or the poromechanics subproblem in terms o coupling iteration dierences The weak orm o the linear momentum balance (2.6) is given by where U(Ω) is given by ( σ, q) Ω + ( q) Ω = 0 ( q U(Ω)) (D.1) U(Ω) { q = (u, v, w) : u, v, w H 1 (Ω), q = 0 on Γ p } D where H m (Ω) is deined, in general, or any integer m 0 as H m (Ω) { w : D α w L 2 (Ω) α m }, where the derivatives are taken in the sense o distributions and given by We know rom tensor calculus that D α α w w = x α, α = α α n, 1.. xαn n ( σ, q) Ω (, σq) Ω (σ : q) Ω (D.2) Further, using the divergence theorem and the symmetry o σ, we arrive at (, σq) Ω (q, σn) Ω (D.3) We decompose q into a symmetric part ( q) s 1 2( q + ( q) T ) ɛ(q) and skewsymmetric part ( q) ss and note that the contraction between a symmetric and skewsymmetric tensor is zero to obtain From (D.1), (D.2), (D.3) and (D.4), we get σ : q σ : ( q) s + 0 σ : ( q) ss = σ : ɛ(q) (σn, q) Ω (σ : ɛ(q)) Ω + (, q) Ω = 0 (D.4) which, ater invoking the traction boundary condition, results in the discrete weak orm or the m th coupling iteration as (t n+1, q h ) Γ p N (σm,n+1 : ɛ(q h )) Ω + ( n+1, q h ) Ω = 0 Replacing m by m + 1 and subtracting the two equations, we get (δ (m) σ : ɛ(q h )) Ω = 0 14
16 Reerences [1] Y. Abousleiman, A. H. D. Cheng, L. Cui, E. Detournay, and J. C. Roegiers. Mandel s problem revisited. Géotechnique, 46(2): , [2] T. Almani, K. Kumar, A. Dogru, G. Singh, and M. F. Wheeler. Convergence analysis o multirate ixed-stress split iterative schemes or coupling low with geomechanics. Computer Methods in Applied Mechanics and Engineering, 311: , [3] F. Armero and J. C. Simo. A new unconditionally stable ractional step method or non-linear coupled thermomechanical problems. International Journal or Numerical Methods in Engineering, 35(4): , [4] M. A. Biot. General theory o three dimensional consolidation. Journal o Applied Physics, 12: , [5] M. A. Biot. Theory o elasticity and consolidation or a porous anisotropic solid. Journal o Applied Physics, 26(2): , [6] M. A. Biot and D. G. Willis. The elastic coeicients o the theory o consolidation. Journal o Applied Mechanics, 24: , [7] F. Brezzi, J. Douglas, R. Durán, and M. Fortin. Mixed inite elements or second order elliptic problems in three variables. Numerische Mathematik, 51(2): , [8] M. M. Carroll. An eective stress law or anisotropic elastic deormation. Journal o Geophysical Research, 84(B13), [9] M. M. Carroll and N. Katsube. The role o terzaghi eective stress in linearly elastic deormation. Journal o Energy Resources Technology, 105(4): , [10] N. Castelletto, J. A. White, and H. A. Tchelepi. Accuracy and convergence properties o the ixed-stress iterative solution o two-way coupled poromechanics. International Journal or Numerical and Analytical Methods in Geomechanics, 39(14): ,
17 [11] A. H. D. Cheng. Material coeicients o anisotropic poroelasticity. International Journal o Rock Mechanics and Mining Sciences, 34: , [12] O. Coussy. Poromechanics. Wiley, 2nd edition, [13] C. A. Felippa, K. C. Park, and C. Farhat. Partitioned analysis o coupled mechanical systems. Computer Methods in Applied Mechanics and Engineering, 190(24): , [14] J. Geertsma. The eect o luid pressure decline on volumetric changes o porous rocks. SPE, 210: , [15] M. E. Gurtin, E. Fried, and L. Anand. The Mechanics and Thermodynamics o Continua. Cambridge University Press, 1 edition, [16] R. Ingram, M. F. Wheeler, and I. Yotov. A multipoint lux mixed inite element method on hexahedra. SIAM Journal o Numerical Analysis, 48(4): , [17] N. Katsube. The constitutive theory or luid-illed porous materials. Journal o Applied Mechanics, 52(1): , [18] J. Kim, H. A. Tchelepi, and R. Juanes. Stability, accuracy and eiciency o sequential methods or coupled low and geomechanics. SPE Journal, 16(2): , [19] A. Mikelić and M. F. Wheeler. Convergence o iterative coupling or coupled low and geomechanics. Computational Geosciences, 17(3): , [20] A. Nur and J. D. Byerlee. An exact eective stress law or elastic deormation o rock with luids. Journal o Geophysical Research, 76(26): , [21] B. A. Schreler, L. Simoni, and E. Turska. Standard staggered and staggered newton schemes in thermo-hydro-mechanical problems. Computer Methods in Applied Mechanics and Engineering, 144(1-2):93 109, [22] A. W. Skempton. The pore-pressure coeicients a and b. Géotechnique, 4(4): ,
18 [23] A. W. Skempton. Signiicance o terzaghi s concept o eective stress (terzaghi s discovery o eective stress). In From theory to practice in soil mechanics (L. Bjerrum and A. Casagrande and R. B. Peek and A. W. Skempton, eds.). New York-London: John Wiley and Sons 1960, [24] J. M. Steele. The Cauchy-Schwarz Master Class: An Introduction to the Art o Mathematical Inequalities. Maa Problem Books Series. Cambridge University Press, [25] E. Turska and B. A. Schreler. On convergence conditions o partitioned solution procedures or consolidation problems. Computer Methods in Applied Mechanics and Engineering, 106(1-2):51 63,
ICES REPORT April Saumik Dana and Mary F. Wheeler
ICES REPORT 18-05 April 2018 Convergence analysis o ixed stress split iterative scheme or small strain anisotropic poroelastoplasticity: A primer by Saumik Dana and Mary F. Wheeler The Institute or Computational
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 informationICES REPORT May Saumik Dana and Mary F. Wheeler
ICES REPORT 8-08 May 208 Te correspondence between Voigt and Reuss bounds and te decoupling constraint in a two-grid staggered solution algoritm or coupled low and deormation in eterogeneous poroelastic
More informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationSplitting methods in the design of coupled flow and mechanics simulators
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
More informationA Multiscale Mortar Multipoint Flux Mixed Finite Element Method
A Multiscale Mortar Multipoint Flux Mixed Finite Element Method Mary F. Wheeler Guangri Xue Ivan Yotov Abstract In this paper, we develop a multiscale mortar multipoint flux mixed finite element method
More informationA solid-fluid mixture theory of porous media
A solid-luid mixture theory o porous media I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro Abstract The theories o mixtures in the ramework o continuum mechanics have been developed
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationMultiphysics 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 informationICES REPORT May Tameem Almani, Kundan Kumar, Gurpreet Singh, Mary F. Wheeler
ICES REPORT 16-12 May 2016 Stability o Multirate Explicit Coupling o Geomecanics wit Flow in a Poroelastic Medium by Tameem Almani, Kundan Kumar, Gurpreet Sing, Mary F. Weeler Te Institute or Computational
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationHEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS
HEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS RICHARD S. FAL, PAOLO GATTO, AND PETER MON Abstract. We study the approximation properties of some finite element subspaces of H(div; Ω) and H(curl ; Ω) defined
More informationarxiv: v1 [math.na] 2 Nov 2018
A multipoint stress mixed finite element method for elasticity II: Quadrilateral grids Ilona Ambartsumyan Eldar Khattatov Jan Nordbotten Ivan Yotov October 30, 2018 arxiv:1811.01928v1 [math.na] 2 Nov 2018
More informationBIOT COEFFICIENTS OF UNIDIRECTIONAL FIBRE PACKED ASSEMBLIES: EXPERIMENTAL AND THEORETICAL
BIOT COEFFICIENTS OF UNIDIRECTIONAL FIBRE PACKED ASSEMBLIES: EXPERIMENTAL AND THEORETICAL Thang Tran, Sébastien Comas-Cardona,*, Nor-Edine Abria 2, Christophe Binétruy Polymers and Composites Technology
More informationMary Fanett Wheeler 1, Guangri Xue 2 and Ivan Yotov 3
ESAIM: M2AN 46 (2012) 759 796 DOI: 10.1051/m2an/2011064 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org A MULTISCALE MORTAR MULTIPOINT FLUX MIXED FINITE ELEMENT METHOD Mary Fanett
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationDouble porosity modeling in elastic wave propagation for reservoir characterization
Stanord Exploration Project, Report 98, September 2, 1998, pages 1 243 Double porosity modeling in elastic wave propagation or reservoir characterization James G. Berryman 1 and Herbert F. Wang 2 keywords:
More informationON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS
ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS VINCENT FONTAINE 1,2 AND ANIS YOUNES 2 1 Laboratoire de Physique des Bâtiments et des Systèmes, Université de La Réunion, 15 avenue
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationOn the relative importance of global and squirt flow in cracked porous media
On the relative importance o global and squirt low in cracked porous media Aamir Ali 1 and Morten Jakobsen 1, 1 Department o Earth Sciences, University o Bergen, Allegaten 41, 5007 Bergen, Norway. Centre
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
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 information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut ür Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 0946 8633 On iterative subdomain methods or the Stokes Darcy problem Alonso Caiazzo
More informationTHE NEAR-WELLBORE PRESSURE CALCULATION MODEL INCORPORATING THERMOCHEMICAL EFFECT
HERMAL CIENCE: Year 2018, Vol. 22, No. 1B, pp. 623-630 623 HE NEAR-WELLBORE PREURE CALCULAION MODEL INCORPORAING HERMOCHEMICAL EFFEC by Zhiqiang ANG, Qian LI *, and Hu YIN Petroleum and Natural Gas Engineering
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More informationСИБИРСКИЕ ЭЛЕКТРОННЫЕ МАТЕМАТИЧЕСКИЕ ИЗВЕСТИЯ
... S e MR ISSN 1813-3304 СИБИРСКИЕ ЭЛЕКТРОННЫЕ МАТЕМАТИЧЕСКИЕ ИЗВЕСТИЯ Siberian Electronic Mathematical Reports http://semr.math.nsc.ru Том 13, стр. 75 88 (2016) УДК 517.95;534.21 DOI 10.17377/semi.2016.13.006
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationDERIVATION OF THE DOUBLE POROSITY MODEL OF SINGLE PHASE FLOW VIA HOMOGENIZATION THEORY*
SIAM J. MATH. ANAL. c 1990 Society or Industrial and Applied Mathematics Vol. 21, No. 4, pp. 823 836, July 1990 001 DERIVATION OF THE DOUBLE POROSITY MODEL OF SINGLE PHASE FLOW VIA HOMOGENIZATION THEORY*
More informationSCIENCE CHINA Physics, Mechanics & Astronomy
SCIENCE CHINA Physics, Mechanics & Astronomy Article June 014 Vol.57 No.6: 1068 1077 doi: 10.1007/s11433-014-5435-z Wave dispersion and attenuation in viscoelastic isotropic media containing multiphase
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139
MASSACHUSTTS INSTITUT OF TCHNOLOGY DPARTMNT OF MATRIALS SCINC AND NGINRING CAMBRIDG, MASSACHUSTTS 0239 322 MCHANICAL PROPRTIS OF MATRIALS PROBLM ST 4 SOLUTIONS Consider a 500 nm thick aluminum ilm on a
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 informationICES REPORT Multirate Undrained Splitting for Coupled Flow and Geomechanics in Porous Media
ICES REPORT 16-13 May 2016 Multirate Undrained Splitting or Coupled Flow and Geomecanics in Porous Media by Tameem Almani, Kundan Kumar, Mary F. Weeler Te Institute or Computational Engineering and Sciences
More informationEFFECTS OF LATERAL MOVEMENTS ON THE DEFORMATION OF SEDIMENTARY BASINS: A NUMERICAL APPROACH
EFFECTS OF LATERAL MOVEMENTS ON THE DEFORMATION OF SEDIMENTARY BASINS: A NUMERICAL APPROACH André Brüch Department o Civil Engineering Federal University o Rio Grande do Sul, Porto Alegre RS, Brazil andrebruch@hotmail.com
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 informationMcMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007
McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007 RANDOM POLYCRYSTALS OF GRAINS WITH CRACKS: MODEL OF ELASTIC BEHAVIOR FOR FRACTURED SYSTEMS James G. Berryman Earth Sciences Division
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 informationSIMULATION OF A FLOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT
SIMULATION O A LOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT R. M. Saldanha da Gama a, J. J. Pedrosa ilho a, and M. L. Martins-Costa b a Universidade do Estado do Rio de Janeiro Departamento
More informationA Total Lagrangian Based Method for Recovering the Undeformed Configuration in Finite Elasticity
A Total Lagrangian Based Method or Recovering the Undeormed Coniguration in Finite Elasticity Grand Roman Joldes*, Adam Wittek and Karol Miller Intelligent Systems or Medicine Laboratory, School o Mechanical
More informationPARCYCLIC: FINITE ELEMENT MODELING OF EARTHQUAKE LIQUEFACTION RESPONSE ON PARALLEL COMPUTERS
13 th World Conerence on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 361 PARCYCLIC: FINITE ELEMENT MODELING OF EARTHQUAKE LIQUEFACTION RESPONSE ON PARALLEL COMPUTERS Jun PENG
More informationStrain and Stress Measurements with a Two-Dimensional Detector
Copyright ISSN (C) 97-, JCPDS-International Advances in X-ray Centre Analysis, or Volume Diraction 4 Data 999 5 Strain and Stress Measurements with a Two-Dimensional Detector Baoping Bob He and Kingsley
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
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 informationAlgorithms for fluid-structure interaction problems
9 Algorithms or luid-structure interaction problems Miguel A. Fernández, Jean-Frédéric Gerbeau Fluid-structure coupled problems and their analysis have been addressed in the previous chapters particularly
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationEXISTENCE OF SOLUTIONS TO SYSTEMS OF EQUATIONS MODELLING COMPRESSIBLE FLUID FLOW
Electronic Journal o Dierential Equations, Vol. 15 (15, No. 16, pp. 1 8. ISSN: 17-6691. URL: http://ejde.math.txstate.e or http://ejde.math.unt.e tp ejde.math.txstate.e EXISTENCE OF SOLUTIONS TO SYSTEMS
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More information(One Dimension) Problem: for a function f(x), find x 0 such that f(x 0 ) = 0. f(x)
Solving Nonlinear Equations & Optimization One Dimension Problem: or a unction, ind 0 such that 0 = 0. 0 One Root: The Bisection Method This one s guaranteed to converge at least to a singularity, i not
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam Abstract. A reormulation o the planetary
More informationChapter 3: Newtonian Fluid Mechanics QUICK START W
Chapter 3: Newtonian Fluid Mechanics TWO GOAL Derive governing equations (mass and momentum balances olve governing equations or velocity and stress ields QUICK TART W First, beore we get deep into derivation,
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationS. Srinivasan, Technip Offshore, Inc., Houston, TX
9 th ASCE Specialty Conerence on Probabilistic Mechanics and Structural Reliability PROBABILISTIC FAILURE PREDICTION OF FILAMENT-WOUND GLASS-FIBER Abstract REINFORCED COMPOSITE TUBES UNDER BIAXIAL LOADING
More informationLocalization in Undrained Deformation
Localization in Undrained Deformation J. W. Rudnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 621-319 John.Rudnicki@gmail.com Fourth Biot Conference on
More informationUnjacketed bulk compressibility of sandstone in laboratory experiments. R. M. Makhnenko 1 and J. F. Labuz 1
481 Unjacketed bulk compressibility of sandstone in laboratory experiments R. M. Makhnenko 1 and J. F. Labuz 1 1 Department of Civil Engineering, University of Minnesota, Minneapolis, MN 55455; PH (612)
More information3.5 Analysis of Members under Flexure (Part IV)
3.5 Analysis o Members under Flexure (Part IV) This section covers the ollowing topics. Analysis o a Flanged Section 3.5.1 Analysis o a Flanged Section Introduction A beam can have langes or lexural eiciency.
More informationRobin-Robin preconditioned Krylov methods for fluid-structure interaction problems
MOX Report No. 19/2008 Robin-Robin preconditioned Krylov methods or luid-structure interaction problems Santiago Badia, Fabio Nobile, Christian Vergara MOX, Dipartimento di Matematica F. Brioschi Politecnico
More informationOpen Access Establishment of Mathematical Model and Sensitivity Analysis of Plugging Capacity of Multi-Component Foam Flooding
Send Orders or Reprints to reprints@benthamscience.ae 84 The Open Chemical Engineering Journal, 2015, 9, 84-89 Open Access Establishment o Mathematical Model and Sensitivity Analysis o Plugging Capacity
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 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 informationTwo-phase flow in a fissurized-porous media
Two-phase low in a issurized-porous media P. Ø. Andersen and S. Evje Citation: AIP Con. Proc. 479, 234 (22); doi:.63/.4756663 View online: http://dx.doi.org/.63/.4756663 View Table o Contents: http://proceedings.aip.org/dbt/dbt.jsp?key=apcpcs&volume=479&issue=
More informationAn Integrated Discrete Fracture Model for Description of Dynamic Behavior in Fractured Reservoirs
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanord University, Stanord, Caliornia, February 24-26, 2014 SGP-TR-202 An Integrated Discrete Fracture Model or Description o Dynamic
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
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 informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationBANDELET IMAGE APPROXIMATION AND COMPRESSION
BANDELET IMAGE APPOXIMATION AND COMPESSION E. LE PENNEC AND S. MALLAT Abstract. Finding eicient geometric representations o images is a central issue to improve image compression and noise removal algorithms.
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationINFLUENCE OF POROSITY AND RADIATION IN A VISCO ELASTIC FLUID OF SECOND ORDER FLUID WITHIN A CHANNEL WITH PERMEABLE WALLS
International Journal o Physics and Mathematical Sciences ISSN: 77- (Online) Vol. () January-March, pp.8-9/murthy and Rajani INFLUENCE OF POROSITY AND RADIATION IN A VISCO ELASTIC FLUID OF SECOND ORDER
More informationAn Iterative Substructuring Method for Coupled Fluid-Solid Acoustic Problems 1
An Iterative Substructuring Method or Coupled Fluid-Solid Acoustic Problems 1 Jan Mandel Department o Mathematics and Center or Computational Mathematics, University o Colorado at Denver, Denver, CO 8217-3364,
More informationA Numerical Study on the Effect of Anisotropy on Hydraulic Fractures
Rock Mechanics and Rock Engineering (219) 52:591 69 https://doi.org/1.17/s63-17-1362-4 ORIGINAL PAPER V. Valliappan 1 J. J. C. Remmers 1 A. Barnhoorn 2 D. M. J. Smeulders 1 Received: 3 April 217 / Accepted:
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More information3D Numerical Modelling of Convective Heat Transfer through Two-sided Vertical Channel Symmetrically Filled with Metal Foams
P Periodica Polytechnica Mechanical Engineering P 60(4), pp. 193-202, 2016 DOI: 10.3311/PPme.8511 Creative Commons Attribution b 3D Numerical Modelling o Convective Heat Transer through Two-sided Vertical
More informationPAijpam.eu NEW H 1 (Ω) CONFORMING FINITE ELEMENTS ON HEXAHEDRA
International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 609-617 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 10.12732/ijpam.v109i3.10
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationShort Note. A simple derivation of the effective stress coefficient for seismic velocities in porous rocks. Boris Gurevich. n = 1 K 0 /K s, (4)
GEOPHYSICS, VOL. 69, NO. 2 (MARCH-APRIL 2004); P. 393 397, 1 FIG. 10.1190/1.1707058 Short Note A simple derivation of the effective stress coefficient for seismic velocities in porous rocks Boris Gurevich
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationGeneral method for simulating laboratory tests with constitutive models for geomechanics
General method for simulating laboratory tests with constitutive models for geomechanics Tomáš Janda 1 and David Mašín 2 1 Czech Technical University in Prague, Faculty of Civil Engineering, Czech Republic
More informationCategories and Natural Transformations
Categories and Natural Transormations Ethan Jerzak 17 August 2007 1 Introduction The motivation or studying Category Theory is to ormalise the underlying similarities between a broad range o mathematical
More informationFILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS
The 8 th International Conerence on Flow Processes in Composite Materials (FPCM8) FILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS S.T.Jespersen
More informationPulling by Pushing, Slip with Infinite Friction, and Perfectly Rough Surfaces
1993 IEEE International Conerence on Robotics and Automation Pulling by Pushing, Slip with Ininite Friction, and Perectly Rough Suraces Kevin M. Lynch Matthew T. Mason The Robotics Institute and School
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationNatural convection in a vertical strip immersed in a porous medium
European Journal o Mechanics B/Fluids 22 (2003) 545 553 Natural convection in a vertical strip immersed in a porous medium L. Martínez-Suástegui a,c.treviño b,,f.méndez a a Facultad de Ingeniería, UNAM,
More informationICES REPORT February Sanghyun Lee, Mary F. Wheeler, Thomas Wick
ICES REPORT 16-04 February 2016 Pressure and fluid-driven fracture propagation in porous media using an adaptive finite element phase field model by Sanghyun Lee, Mary F. Wheeler, Thomas Wick The Institute
More informationSTAT 801: Mathematical Statistics. Hypothesis Testing
STAT 801: Mathematical Statistics Hypothesis Testing Hypothesis testing: a statistical problem where you must choose, on the basis o data X, between two alternatives. We ormalize this as the problem o
More informationFiltration. Praktikum Mechanical Engineering. Spring semester ML F16 Tel.:
Praktikum Mechanical Engineering Spring semester 2018 Filtration Supervisor: Davide Stucchi ML F16 stucchid@ptl.mavt.ethz.ch Tel.: 044 632 25 05 1 1 Table o Contents 1 TABLE OF CONTENTS... 2 2 INTRODUCTION...
More information