ICES REPORT May Tameem Almani, Kundan Kumar, Gurpreet Singh, Mary F. Wheeler
|
|
- Edith Georgina Anthony
- 6 years ago
- Views:
Transcription
1 ICES REPORT 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 Engineering and Sciences Te University o Texas at Austin Austin, Texas Reerence: Tameem Almani, Kundan Kumar, Gurpreet Sing, Mary F. Weeler, "Stability o Multirate Explicit Coupling o Geomecanics wit Flow in a Poroelastic Medium," ICES REPORT 16-12, Te Institute or Computational Engineering and Sciences, Te University o Texas at Austin, May 2016.
2 Stability o multirate explicit coupling o geomecanics wit low in a poroelastic medium T. Almani 1, K. Kumar 2, G. Sing 1, M. F. Weeler 1 1 Center or Subsurace Modeling, ICES, UT Austin, USA 2 Matematics Institute, University o Bergen, Norway {tameem,gurpreet,mw}@ices.utexas.edu, kundan.kumar@uib.no May 11, 2016 Abstract We consider single rate and multirate explicit scemes or te Biot system modeling coupled low and geomecanics in a poro-elastic medium. Tese scemes are te most widely used in practice tat ollows a sequential procedure in wic te low and mecanics problems are ully decoupled. In suc a sceme, te low problem is solved irst wit time-lagging te displacement term ollowed by te mecanics solve. Te multirate explicit coupling sceme exploits te dierent time scales or te mecanics and low problems by taking multiple iner time steps or low witin one coarse mecanics time step. We provide ully discrete scemes or bot te single and multirate approaces tat use Backward Euler time discretization and mixed spaces or low and conormal Galerkin or mecanics. We perorm a rigorous stability analysis and derive te conditions on reservoir parameters and te number o iner low solves to ensure stability or bot scemes. Furtermore, we investigate te computational time savings or explicit coupling scemes against iterative coupling scemes. Keywords. poroelasticity; Biot; iterative and explicit coupling; multirate sceme; mixed ormulation 1 Introduction Te coupling between subsurace low and reservoir geomecanics plays a critical role in obtaining accurate results or models involving reservoir deormation, surace subsidence, well stability, sand production, waste deposition, ydraulic racturing, CO 2 sequestration, and ydrocarbon recovery [2, 8, 14]. Te quasi-static Biot equations are used to model te subsurace coupled low and mecanics and consists o a system o two coupled linear partial dierential equations, eac o wic is typically associated wit te low and mecanics, respectively. Quite oten in practice, te geomecanics problem as a muc slower evolution tan tat o te low problem. In suc cases, te mecanics problem can cope wit a muc coarser time step compared to te low problem. Te multirate sceme exploits tis dierence in te two equations and allows te low to take several iner time steps beore updating te mecanics and is a natural candidate in tis setting. Figures 1.1a and 1
3 1.1b illustrate te dierences between single rate versus multirate explicit coupling scemes. Figure 1.1a represents a typical single rate sceme, in wic te low and mecanics problems sare te exact same time step. In contrast, Figure 1.1b demonstrates a typical multirate sceme, in wic te low problem takes multiple iner local time steps witin one coarser mecanics time step. Te explicit coupling approac tat we consider ere is te most widely utilized sceme in practice. Te decoupling o te two equations makes it easy to implement and te time marcing witout any iterations leads to a lower computational cost. Te drawback is tat tis sceme is only conditionally stable. For te single rate sceme, te rigorous stability properties ave been investigated in te work o Mikelić and Weeler [15]. However, in te case wen te multiple low time steps are taken or one mecanics time step, it is unclear ow tese stability properties cange. In tis work, we ocus our attention on te explicit coupling approac, establis its stability teoretically or bot ully discrete single rate and multirate scemes, and investigate its computational time savings numerically. Moreover, in contrast to te explicit coupling approac, te iterative coupling approac as been investigated in te past. In tis approac, te two coupled subsystems are solved iteratively by excanging te values o te sared state variables in an iterative manner. Te procedure is iterated at eac time step until te solution is obtained wit an acceptable tolerance [5, 13, 14]. Multirate iterative coupling scemes, extending te ixed-stress split coupling algoritm, ave been rigorously designed and analyzed in [1, 2]. Unconditional stability o suc scemes ollows immediately by establising teir Banac contraction properties. In addition, multirate iterative coupling scemes, based on te undrained-split coupling algoritm, are sown to be contractive, and tus unconditionally stable [12]. In tis work, we consider explicit coupling scemes and rigorously establis teir stability properties. In addition, we perorm numerical computations on ield scale problems to compare te eiciency and computational perormance o tese two approaces. Te coupled low and geomecanics problem as been intensively investigated in te past pioneered by Terzagi [21] and Biot [3, 4]. Terzagi was te irst to propose an explanation o te soil consolidation process, in wic e assumed tat grains orming te soil are bound togeter by some molecular orces resulting in te ormation o te porous material wit elastic properties. It is te success o Terzag s teory in predicting te settlement o dierent types o soils tat led to te creation o te science o soil mecanics [4]. Biot ten extended Terzag s one dimensional work to te tree-dimensional case, and presented a more rigorous generalized teory o consolidation [4]. A compreensive treatment o poromecanics and te teory o mecanics o porous continua can be ound in [6] by Coussy. Oter nonlinear extensions on te teory o poroelasticity can be ound in [7, 17, 20]. Recently, te work o Mikelić and Weeler [15] establises stability and geometric convergence contraction wit respect to appropriately cosen metrics or dierent low and geomecanics coupling scemes. In addition, Kim et al. [10, 11] ave used von Neumann stability analysis to study te stability and convergence o similar scemes. Te multirate scemes or te non-stationary Stokes-Darcy model ave been investigated in [18, 24]. In tis work o multirate explicit coupling o low wit geomecanics, we establis stability results or bot te single rate and multirate scemes, and investigate teir accuracies and computational time savings numerically. To te best o our knowledge, tis is te irst 2
4 a Single Rate b Multirate Figure 1.1: Flowcart or explicit single rate and multirate time stepping or coupled geomecanics and low problem analysis o te multirate explicit coupling sceme or Biot equations. Te paper is structured as ollows. We present te model and discretization in Section 2. Te single rate and multirate explicit coupling scemes are introduced and analyzed in Section 3. Numerical results are sown in Section 4. Conclusions and outlook are discussed in Section Preliminaries Let Ω be a bounded domain open and connected o IR d, were te dimension d = 2 or 3, wit a Lipscitz continuous boundary Ω, and let Γ be a part o Ω wit positive measure. Wen d = 3, we assume tat te boundary o Γ is also Lipscitz continuous. As usual, we denote by H 1 Ω te classical Sobolev space equipped wit te semi-norm and norm: We also deine: H 1 Ω = {v L 2 Ω ; v L 2 Ω d }, v H 1 Ω= v L 2 Ω d, v H 1 Ω= v 2 L 2 Ω + v 2 H 1 Ω 1/2. H 1 0 Ω = {v H 1 Ω ; v Ω = 0}, and or te divergence operator, we sall use te space equipped wit te norm Hdiv; Ω = {v L 2 Ω d ; v L 2 Ω}, v Hdiv;Ω = v 2 L 2 Ω + v 2 L 2 Ω 1/2. 3
5 For a vector v in IR d, recall te strain or symmetric gradient tensor εv: εv = 1 2 v + vt. In te sequel we sall use Poincaré s and Korn s inequalities. Poincaré s inequality in H 1 0 Ω reads: Tere exists a constant P Ω depending only on Ω suc tat v H 1 0 Ω, v L 2 Ω P Ω v H 1 Ω. 1.1 Next, recall Korn s irst inequality in H 1 0 Ωd : Tere exists a constant C κ depending only on Ω suc tat v H 1 0 Ω d, v H 1 Ω C κ εv L 2 Ω Model equations, discretization and splitting algoritm We assume a linear, elastic, omogeneous, and isotropic poro-elastic medium Ω R d, d = 2 or 3, in wic te reservoir is saturated wit a sligtly compressible viscous luid. 2.1 Assumptions Te luid is assumed to be sligtly compressible and its density is a linear unction o pressure, wit a constant viscosity µ > 0. Te reerence density o te luid ρ > 0, te Lamé coeicients > 0 and G > 0, te dimensionless Biot coeicient α, and te pore volume ϕ are all positive. Te absolute permeability tensor, K, is assumed to be symmetric, bounded, uniormly positive deinite in space and constant in time. A quasi-static Biot model [4, 8] will be employed in tis work. Te model reads: Find u and p satisying te equations below or all time t ]0, T [: Flow Equation: t 1 M + c ϕ 0 p + α u Mecanics Equations: Boundary Conditions: 1 µ K p ρ,r g η = q in Ω div σ por u, p = in Ω, σ por u, p = σu α p I in Ω, σu = ui + 2 Gεu in Ω u = 0, K p ρ,r g η n = 0 on Ω Initial Condition t = 0 : 1 M + c ϕ 0 p + α u 0 = 1 M + c ϕ 0 p 0 + α u 0. were: g is te gravitational constant, η is te distance in te vertical direction assumed to be constant in time, ρ,r > 0 is a constant reerence density relative to te reerence 4
6 pressure p r, ϕ 0 is te initial porosity, M is te Biot constant, q = were q is a mass source or sink term taking into account injection into or out o te reservoir. We remark tat te above system is linear and coupled troug te Biot coeicient terms. 2.2 Mixed variational ormulation A mixed inite element ormulation or low and a conormal Galerkin ormulation or mecanics will be used. Te mixed ormulation is a locally mass conservative sceme, and allows or explicit lux computation. Te lux is deined as a separate unknown and te low equation is rewritten as a system o irst order equations. Accordingly, or te ully discrete ormulation discrete in time and space, we use a mixed inite element metod or space discretization and a backward-euler time discretization. Let T denote a regular amily o conorming triangular elements o te domain o interest, Ω. Using te lowest order Raviart-Tomas RT spaces, we ave te ollowing discrete spaces V or discrete displacements, Q or discrete pressures, and Z or discrete velocities luxes: q ρ,r V = {v H 1 Ω d ; T T, v T P 1 d, v Ω = 0} 2.1 Q = {p L 2 Ω ; T T, p T P 0 } 2.2 Z = {q Hdiv; Ω d ; T T, q T P 1 d, q n = 0 on Ω} 2.3 Te space o displacements, V, is equipped wit te norm: v V = d v i 2 Ω 1/2. i=1 We also assume tat te iner time step is given by: = t k t k 1. I we denote te total number o timesteps by N, ten te total simulation time is given by T = N, and t i = i, 0 i N denote te discrete time points. For te ully discrete sceme, we ave cosen te Raviart-Tomas spaces or te mixed inite element discretization. However, te proo extends to oter coices or te mixed spaces, and we will state te results or Multipoint Flux Mixed Finite Element MFMFE spaces [23] in Remark 3.6. Remark 2.1 Notation: We adopt te ollowing notations,k denotes te coarser time step iteration index or indexing mecanics coarse time steps, m is te iner local time step iteration index or indexing low ine time steps, stands or te unit iner time step, and q is te ixed number o local low time steps per coarse mecanics time step. 2.3 Fully discrete sceme or single rate As discussed above, using te mixed inite element metod in space and te backward Euler inite dierence metod in time, te weak ormulation o te single rate sceme reads as ollows. 5
7 Deinition 2.2 low equation Find p k+1 Q, and z k+1 Z suc tat, uk uk 1 θ Q, 1 p k+1 M + c ϕ 0 p k, θ + 1 µ z k+1, θ + α q Z, K 1 z k+1, q = p k+1, q + ρ,r gη, q, θ = q, θ Deinition 2.3 mecanics equation Find u k+q V suc tat, v V, 2Gεu k+1, εv + u k+1, v αp k+1, v = k+1, v Fully discrete sceme or multirate Te weak ormulation o te multirate sceme reads as ollows. Deinition 2.4 low equation For 1 m q, ind p m+k Q, and z m+k Z suc tat, 1 θ Q, 1 M + c ϕ 0 p m+k, θ + 1 z m+k µ, θ = α u k q uk q, θ + q, θ, 2.7 q Z, K 1 z m+k, q = p m+k, q + Deinition 2.5 mecanics equation Find u k+q v V, 2Gεu k+q, εv + u k+q ρ,r g η, q, 2.8 V suc tat,, v αp k+q, v = k+q, v. 2.9 wit te initial condition or te irst discrete time step or bot single rate and multirate scemes, p 0 = p Note tat or te multirate sceme, te pressure unknowns p and lux unknowns z are being solved at iner time steps t k+m, m = 0,..., q wereas te mecanics variables u are being solved at t iq, i N. Tereore, or eac mecanics time step o size q, tere are q low solves justiying te nomenclature o multirate. Moreover, te above system o PDEs is linear, decoupled and te inormation excange taking place at te coarse time steps. 3 Analysis o Explicit Coupling Scemes In tis section, we perorm a matematically rigorous analysis o te stability o te single rate and multirate explicit coupling scemes. Recall tat in te single rate case Figure 1.1a, te low and mecanics problems sare te exact same time step. In contrast, in te multirate case Figure 1.1b, te low problem takes q iner local time steps witin one coarser mecanics time step. As as been stated above, te explicit coupling approac is a sequential procedure in wic te low or te mecanics problem is solved irst ollowed by te oter. Tere is no coupling iteration between te two problems. 6
8 3.1 Single Rate Formulation: We start by analyzing te single-rate explicit coupling algoritm, in wic bot low and mecanics sare te same time step. To te best o our knowledge, tis is te irst rigorous matematical analysis o te ully discrete single-rate explicitly coupled Biot system. In addition, te analysis reveals a more general stability condition compared to te one obtained in [16] by elementary means. Te algoritm is given as ollows: Algoritm 1: Single Rate Explicit Coupling Algoritm 1 Given initial conditions u 0 and p0, solve ully implicitly or p1, u1 satisying Biot model 2 or k = 1, 2,... do /* time step index */ 3 First Step: Flow equations 4 Given u k and uk 1 : 5 Solve or p k+1 and z k+1 satisying deinition Second Step: Mecanics equations 7 Given p k+1 and, z k+1 : 8 Solve or u k+1 satisying deinition 2.3 Note tat we begin wit k = 1 and we require bot u 1 and u0 or obtaining p2. In te irst step, we use a ully implicit metod to solve or p 1, u1. Alternatively, to keep te problem decoupled, we can use iterative tecniques suc as ixed stress splitting or undrained splitting [15] Assumptions For notational convenience, we deine β = 1 M + c ϕ 0. For stability to old, we assume te ollowing: Result A 1 β > α2. Our results make explicit te dependence o te stability on te dierence o te above quantities. we ave te ollowing stability result. Teorem 3.1 [Single rate] Under te Assumption A 1 above, te ollowing stability result olds or te single rate explicit coupling sceme or time steps t 0 t k t J : K 1/2 z J+1 µ 2 + K 1/2 z k+1 z k 2 + 2G εu k+1 u k 2 k=1 k=1 + u J+1 u J P 2 C + 2β α2 q + Ω C2 κ 2G k=1 k+1 k=1 k 2 7
9 3.1.3 Stability Analysis: Te proo o te above teorem is carried out in tree steps by considering te low equation, te mecanics equation and ten combining te two togeter. Recall tat β = 1 M + c ϕ 0. Proo. Step 1: Flow equations Testing 2.4 wit θ = p k+1 p k, we obtain β 1 p k+1 p k z k+1 µ, p k+1 p k + α u k uk 1, p k+1 = p k q, p k+1 p k 3.1 Next, we consider te lux equation 2.5. Taking te dierence o two consecutive time steps t = t k+1 and t = t k and testing wit q = z k+1, we obtain: K 1 z k+1 z k, zk+1 = p k+1 p k, zk Substituting 3.2 into 3.1, ater some algebraic manipulations o te resulting term using: aa b = 1 2 a2 b 2 + a b 2, we derive β p k+1 p k K 1/2 z k+1 2µ 2 K 1/2 z k 2 + K 1/2 z k+1 z k 2 + α u k uk 1, p k+1 p k = q, p k+1 p k 3.3 Step 2: Elasticity equation Considering 2.6 or te dierence o two consecutive time steps, t = t k+1 and t = t k, and testing wit v = uk+1 u k, we obtain 2G εu k+1 u k 2 + u k+1 u k 2 α p k+1 = 1 p k, uk+1 u k k+1 k, uk+1 u k 3.4 8
10 Step 3: Combining low and elasticity equations Combining 3.3 wit 3.4 yields β p k+1 p k K 1/2 z k+1 2µ 2 K 1/2 z k 2 + K 1/2 z k+1 z k 2 + 2G εu k+1 u k 2 + u k+1 u k 2 α = u k uk 1, p k+1 p k }{{} + α p k+1 p k, uk+1 u k + q, p k+1 p k }{{}}{{} R 2 R 3 + R 1 1 k+1 k, uk+1 u k }{{} R Denoting by R 1, R 2, R 3, and R 4 te terms on te rigt and side, togeter wit Poincaré s, Korn s, and Young s inequalities, we estimate R 1 α 1 u k 2ɛ u k α ɛ 1 p k p k 2 R 2 α u k+1 2ɛ u k 2 αɛ 2 + p k p k 2 R ɛ 3 q + p k+1 2ɛ 3 2 p k 2 R 4 1 k+1 2ɛ k 2 + ɛ 4 u k u k 2 1 k+1 2ɛ k 2 + ɛ 4PΩ 2 C2 κ εu k u k 2. or ɛ 1, ɛ 2, ɛ 3, and ɛ 4 > 0. Coosing ɛ 1 = ɛ 2 = α, ɛ 3 = 2 α2 β, ɛ 4 = multiplying 3.5 by 2, we derive K 1/2 z k+1 µ 2 K 1/2 z k 2 + K 1/2 z k+1 z k 2 + 2G + u k+1 u k 2 u k u k P + q + Ω 2 C2 κ 2β 2α 2 2G 2G P 2 Ω C2 κ and εu k+1 u k 2 k+1 k Summing up 3.6 or 1 k J, or J time steps, wit telescopic cancellations, we get: K 1/2 z J+1 µ 2 + K 1/2 z k+1 z k 2 + 2G k=1 εu k+1 u k 2 k=1 + u J+1 u J 2 u 1 u 0 2 K + 1/2 z µ + P2 Ω C2 κ 2G 9 k+1 k=1 2 2β α2 2 q k=1 k 2, 3.7
11 Recall tat u 1, z1 ave been computed using te ully implicit time discretization. Using standard a priori estimates or te coupled Biot model, we conclude tat u 1 u0 C 2 and K 1/2 z 1 2 C. Tis completes te derivation. Remark 3.2 Te above proo also provides a way to devise an explicitly coupled algoritm tat is unconditionally stable. For te single rate algoritm, we replace 2.4 by te ollowing equation: low equation Find p k+1 Q and z k+1 Z suc tat, θ Q, 1 M + c ϕ 0 + α2 p k+1 p k, θ + 1 z k+1 µ, θ + α uk uk 1, θ = 3.8 Note tat te stabilisation term α2 pk+1 p k as been added above in contrast to 2.4. Te stability result is ten obtained wit te assumption A 1 relaxed. Te consistence error is expected to be o order O wic is also expected or te sceme. To see te unconditional stability o te new sceme, consider te analog o 3.5 and proceed as in te previous case, β + α2 p k+1 p k K 1/2 z k+1 2µ 2 K 1/2 z k 2 + K 1/2 z k+1 z k 2 + 2G εu k+1 u k 2 + u k+1 u k 2 α = u k uk 1, p k+1 p k }{{} + α p k+1 p k, uk+1 u k + q, p k+1 p k }{{}}{{} R 2 R 3 + R 1 1 k+1 k, uk+1 u k. }{{} R Denoting by R 1, R 2, R 3, and R 4 te terms on te rigt and side, togeter wit Poincaré s, Korn s, and Young s inequalities, we estimate R 1 α 1 u k 2ɛ u k α ɛ 1 p k p k 2 R 2 α u k+1 2ɛ u k 2 αɛ 2 + p k p k 2 R ɛ 3 q + p k+1 2ɛ 3 2 p k 2 R 4 1 k+1 2ɛ k 2 + ɛ 4 u k u k 2 1 k+1 2ɛ k 2 + ɛ 4PΩ 2 C2 κ εu k u k 2. q, θ. 10
12 or ɛ 1, ɛ 2, and ɛ 4 > 0. Coosing ɛ 1 = α, ɛ 2 = α, ɛ 3 = 2β, and ɛ 4 = 2G and multiplying PΩ 2 C2 κ 3.9 by 2, we derive K 1/2 z k+1 µ 2 K 1/2 z k 2 + K 1/2 z k+1 z k 2 + 2G εu k+1 u k 2 + u k+1 u k 2 u k u k P 2 q + Ω C2 κ k+1 2β 2G k and rest o te steps proceeds as ollows. 3.2 Multirate Formulation: Recall tat in te multirate explicit coupling approac, te low problem is solved q times wit a iner time step witin a coarser mecanics time step. Algoritm 2: Multirate Explicit Coupling Algoritm 1 Given initial conditions u 0 and p0, solve implicitly or um, pm, zm, m = 1, 2,..., q satisying ully coupled multirate Biot model 2 or k = q, 2q, 3q,.. do /* mecanics time step iteration index */ 3 First Step: Flow equations 4 Given u k 5 or m = 1, 2,.., q do /* low iner time steps iteration index */ 6 Solve or p m+k and z m+k satisying deinition Second Step: Mecanics equations 8 Given p k+q and, z k+q 9 Solve or u k+q satisying deinition Assumptions Te stability assumption in te multirate case takes te orm: A q β > q + q α2 or q 1, were q is te number o low iner time steps witin one coarse mecanics time step. As in te single rate case, we need to prepare te initial data or starting te time stepping. Accordingly, in te irst step o te multirate algoritm Algoritm 2, or k = 0, and m = 1, 2,..., q, te initial conditions are computed by solving te coupled Biot system wit ully implicit time discretization wit a time step o size or te q coupled solves. Alternatively, decoupled iterative scemes [2, 12] suc as ixed stress iterative single rate sceme can be used to compute u m, pm, zm, m = 1, 2,..., q. Note tat i q = 1, te multirate condition A q is identical to te single rate condition A 1. Our main result is te ollowing stability estimate. 11
13 Teorem 3.3 [Multirate] Under te assumption A q, te ollowing stability result olds or te multirate explicit coupling sceme or mecanics time steps t 0 t k t J, k = q, 2q,..: 2G εu k+q u k 2 K + 1/2 z J+q µ 2 + k=q + u J+q u J 2 q 2 C + 2 β α2 q + q 3.3 Stability Analysis k=q m=1 K 1/2 z m+k 2 P 2 q + Ω C2 κ 2G k=q z m 1+k k+q k=q 2 k Te proo or te stability analysis ollows te same ideas as in te single rate proo, owever te use o multiple time steps requires additional estimates. We ollow te same principle o estimating te low equation ollowed by mecanics equation and ten combining te two togeter to obtain te stability estimates. Proo. Step 1: Flow equations Testing 2.7 wit θ = p m+k β p m+k z m+k µ + α u k q uk q, we get, p m+k, p m+k = q, p m+k 3.12 In te lux equation 2.8, considering te dierence or two consecutive iner time steps t = t m+k and t = t m 1+k, and testing wit q = z m+k, we obtain K 1 z m+k z m 1+k, z m+k = p m+k, z m+k Substituting 3.13 into 3.12, we derive βp m+k 2 + K 1 z m+k µ = α u k q uk q z m 1+k, p m+k, z m+k + q, p m+k Summing across low iner time steps 1 m q, we get use aa b = 1 2 a2 b 2 + a b 2 and te telescopic cancellations β + m=1 m=1 p m+k K 1/2 z m+k 2 + K 1/2 z k+q 2µ 2 K 1/2 z k z m 1+k 2 = α q + q, u k uk q, p m+k m=1 2 p m+k m=
14 Step 2: Elasticity equation Considering 2.9 or te dierence o two consecutive mecanics time steps, t = t k and t = t k+q, and testing wit v = u k+q u k, we obtain 2Gεu k+q u k 2 + u k+q u k 2 αp k+q p k, uk+q u k = Step 3: Combining low and elasticity equations Combining 3.14 wit 3.15 gives β m =1 p m+k 2 + 2Gεu k+q + K 1/2 z k+q 2µ 2 K 1/2 z k 2 + m=1 = α u k q uk q, p m+k m=1 }{{} R 1 + αp k+q p k, uk+q u k }{{ } R 3 + u k 2 + u k+q u k 2 K 1/2 z m+k + q, k+q k, uk+q u k }{{} k+q k, uk+q u k 3.15 z m 1+k p m+k m=1 2 }{{} R 4. R Denoting by R 1 and R 2 te irst two terms on te rigt and side, Youngs s and triangle s inequalities give R 1 α q ɛ1 2 ɛ2 R 2 2 m=1 m=1 p m+k p m+k Using te act tat p k+q p k = m=1 p m+k q 2ɛ 1 u k u k q q 2 q 2ɛ 2 2, togeter wit Young s and triangle s inequalities, te tird term on te rigt and side o 3.16, denoted by R 3, can be written as R 3 αɛ 3 2 m=1 p m+k 2 + qα u k+q 2ɛ u k
15 By Poincaré s, Korn s, and Young s inequalities, te last term on te rigt and side o 3.16, denoted by R 4, can be written as R 4 1 k+q 2ɛ 4 Coosing ɛ 1 = α, ɛ 2 = 2 2 we derive 1 k+q 2ɛ 4 β q K 1/2 z k+q µ 2 K 1/2 z k G k 2 + ɛ 4 u k+q 2 u k 2 k 2 + ɛ 4PΩ 2 C2 κ εu k+q 2 u k 2. α2 + q, ɛ 3 = qα, ɛ 4 = 2G, and multiplying by PΩ 2 C2 κ m=1 K 1/2 z m+k z m 1+k 2 εu k+q u k 2 + u k+q u k 2 u k u k q + q 2 2 β q + q α2 2 P 2 q + Ω C2 κ 2G k+q 2 k We need to impose te ollowing condition: β α2 q + q > 0, wic is noting but te Assumption A q. Summing up equation 3.18 or q k J k is a multiple o q, tat is, k = q, 2q,.., we write 2G εu k+q u k 2 K + 1/2 z J+q µ 2 + K 1/2 z m+k z m 1+k 2 k=q m=1 + u J+q u J 2 K 1/2 z q µ 2 + u q u0 2 k=q q β α2 q + q 2 P 2 q + Ω C2 κ 2G k=q k+q k=q k To estimate te irst two terms on te rigt and side, we need to obtain a priori estimates or te ully implicit sceme or te multirate Biot. Tis a priori estimate is obtained by a sligt variation o te tecnique rom te single rate sceme and yields u q u0 2 Cq 2 2 and K 1/2 z q C. We spare te details o obtaining tese a priori estimates. Putting togeter, we conclude te result. Remark 3.4 As in te single rate case in remark 3.2, te multirate case can also be made unconditionally stable by adding a stabilisation term. In te deinition 2.4, we modiy te low equation 2.7 by adding a stabilisation term γ α2 pm+k p m 1+k, were γ = q +q. Te modiied equation reads: 14
16 low equation For 1 m q, ind p m+k Q, and z m+k Z suc tat, θ Q, 1 1 M + c ϕ 0 + γα2 α q u k uk q p m+k, θ +, θ + 1 µ z m+k, θ = q, θ Te proo or te unconditional stability ollows te same ideas as in te single rate case and is skipped ere. Remark 3.5 For te numerical simulations we will be using te multipoint lux mixed inite element metod MFMFE [22,23] or te low discretization. All our obtained results remain valid or tis discretization. Indeed, or suc a sceme, te stability results 3.19 translates to, µ + 2G K 1 z J+q εu k+q k=q, z J+q + Q µ k=q m=1 u k 2 + u J+q + u q u0 2 q β α2 q + q K 1 z m+k z m 1+k, z m+k z m 1+k Q u J 2 K 1 z q µ, zq Q 2 P 2 q + Ω C2 κ 2G k=q k+q k=q k 2, 3.21 were K 1.,. Q is te quadrature rule deined in [23] or te MFMFE corresponding spaces. It was sown by Weeler and Yotov in [23], and ten extended to distorted quadrilaterals and exaedra in [22], tat or any z Z, C 1 z 2 K 1 z, z Q C 2 z 2, or a constant C 1, C 2 > 0. Tis immediately leads to a similar stability result. Te same argument olds or single rate case. Remark 3.6 Te well source/sink term q can be assumed to be varying wit discrete ine/coarse time steps, and all obtained results remain valid. 4 Numerical Results 4.1 Iterative vs. Explicit Coupling Scemes In tis section, we compare single rate and multirate explicit coupling scemes versus iterative coupling scemes. Bot scemes are implemented in te Integrated Parallel Accurate Reservoir Simulator IPARS on top o a single-pase low model coupled wit a linear poroelasticity model. Te Multipoint Flux Mixed Finite Element Metod MFMFE is used or low discretization and Conormal Galerkin is used or elasticity discretization. Mikelić and Weeler [15] ave analyzed dierent iterative coupling scemes, and ave sown tat te two oten used tecniques known as te ixed-stress split and te undrained-split coupling algoritms are unconditionally stable. Te numerical computations in [14] sow te relative perormances o te two metods wit ixed stress splitting perorming better. 15
17 Total Simulation time: 64.0 days k xx Range: , md Finer Unit time step: 1.0 days k yy Range: , md Number o grids: grids k zz Range: , md Possion Ratio, ν 0.35 Biot Modulus, M 1.E15 psi Biot s constant, α 0.9 = Eν 1+ν1 2ν 4.32E7 psi Initial porosity, ϕ Flow Boundary Conditions: zero low B.C. Fluid viscosity, µ w 1.0 cp Mecanics B.C.: Initial luid pressure, p psi X+ boundary σ xx = σ n x = 10, 000psi luid compressibility c w : 1.E-6 1/psi X- - boundary u = 0, zero displacement Rock compressibility: 1.E-6 1/psi Y+ - boundary u = 0, zero displacement Rock density: lb m/t 3 Y- - boundary σ yy = σ n y = 2000psi Initial luid density, ρ o lb m/t 3 Z+ - boundary u = 0, zero displacement Young s Modulus E 5.0E7 psi Z- - boundary σ zz = σ n z = 1000psi Table 1: Input Parameters or Brugge Field Model % o Reduction in: q = 1 q = 2 q = 4 q = 8 CPU run time 47.34% 47.62% 44.97% 51.69% Number o low linear iterations 48.50% 52.92% 52.16% 53.55% Number o mecanics linear iterations 52.21% 52.20% 51.50% 52.89% Table 2: Computational savings o explicit coupling scemes versus iterative coupling scemes or dierent values o q te number o low ine time steps witin one coarse mecanics time step. In te multirate case te unconditional stability o tese two scemes ave been studied in [2, 12]. For our numerical tests, we consider te iterative ixed-stress coupling algoritm wen comparing te eiciency o te iterative coupling scemes versus explicit coupling scemes Brugge Fileld Model We consider te Brugge ield model [19] or comparing te accuracy and eiciency o iterative versus explicit coupling scemes. Te model consists o a general exaedral elements capturing te ield geometry, wit 30 bottom-ole pressure speciied wells, 10 o wic are injectors at a pressure o 4600 psi, and 20 are producers at a pressure o 1200 psi. Producers are located at a lower elevation compared to injectors. No low boundary condition is enorced across all external boundaries. For te mecanics model, we apply a mixture o zero displacement and traction boundary conditions. we also include te eects o gravity. Detailed speciications o te input parameters can be ound in Table 1. We note ere tat assumptions A 1 and A q are bot satisied or te single rare and multirate explicit coupling cases q = 1, 2, 4, and 8, respectively. 16
18 Accumulated CPU Time minutes Accumulated CPU Run Time vs Simulation Period Iterative Coupling - Single Rate Iterative Coupling - Multirate q = 2 Iterative Coupling - Multirate q = 4 Iterative Coupling - Multirate q = 8 Explicit Coupling - Single Rate Explicit Coupling - Multirate q = 2 Explicit Coupling - Multirate q = 4 Explicit Coupling - Multirate q = 8 Accumulated Mecanics Linear Iterations Accumulated # o Mecanics Linear Itrns vs Simulation Period Iterative Coupling - Single Rate Iterative Coupling - Multirate q = 2 Iterative Coupling - Multirate q = 4 Iterative Coupling - Multirate q = 8 Explicit Coupling - Single Rate Explicit Coupling - Multirate q = 2 Explicit Coupling - Multirate q = 4 Explicit Coupling - Multirate q = Simulation Period days a CPU Run Time vs Simulation Days Simulation Time days b Total Number o Mecanics Linear Iterations vs Simulation Days Accumulated Flow Linear Iterations Accumulated # o Flow Linear Itrns vs Simulation Period Iterative Coupling - Single Rate Iterative Coupling - Multirate q = 2 Iterative Coupling - Multirate q = 4 Iterative Coupling - Multirate q = 8 Explicit Coupling - Single Rate Explicit Coupling - Multirate q = 2 Explicit Coupling - Multirate q = 4 Explicit Coupling - Multirate q = Simulation Time days c Total Number o Flow Linear Iterations vs Simulation Days Figure 4.1: Brugge Field Model Numerical Results 17
19 a Pressure ield at 64.0 days Iterative Coupling b Displacement ield at 64.0 days Iterative Coupling c Pressure ield at 64.0 days Explicit Coupling d Displacement ield at 64.0 days Explicit Coupling Figure 4.2: Iterative vs Explicit Coupling Results: Background: Pressure Proile, Arrows: Mecanical Displacements 18
20 4.1.2 Results & Discussion Figure 4.1a sows te accumulated CPU run time or te single rate case q = 1, and or multirate cases: q = 2, 4, and 8, or bot iterative and explicit coupling scemes. In general, or a ixed q, explicit coupling scemes are more eicient, compared to teir counterpart iterative coupling scemes. Tis is expected as explicit scemes eliminate any coupling iteration between te two problems. Tis results in a uge reduction in te total number o low and mecanics linear iterations or explicit coupling scemes, as sown in Figures 4.1c, and 4.1b respectively. Te results obtained sow tat explicit coupling scemes can reduce te accumulative number o low linear iterations or te wole simulation run by almost 50.0% compared to iterative coupling scemes. In addition, te accumulative number o mecanics linear iterations is reduced as well wen comparing an explicit coupling sceme to an iterative sceme or a ixed value o q. As sown in igure 4.1b, te single rate iterative coupling sceme results in te igest number o total mecanics linear iterations. In contrast, te multirate explicit coupling sceme q=8 results in te lowest number o mecanics linear iterations or te wole simulation run. Computational savings o explicit coupling scemes versus iterative coupling scemes are sown in Table 2. Figures 4.2a and 4.2b sow te pressure and displacement ields or te iterative coupling sceme ater 64.0 days o simulation o te Brugge ield case. Figures 4.2c and 4.2d sow te corresponding ields or te explicit coupling sceme. Te solutions or bot te approaces are airly close wit a sligt dierence between te iterative and explicit coupling being more apparent or pressure ields. Te dierences in displacement ields or bot scemes are negligible. 4.2 Validating Teoretical Assumptions In tis section, we try to validate our teoretically induced assumptions or te single rate and multirate explicit coupling scemes against te Frio ield model. Located on te Gul Coast, near Dayton, Texas, at Sout Liberty oil ield, te Frio ield model is a ield-scale problem wit a geometrically callenging geological ormation [9]. Te ield is curved in te dept direction, wit several tin curved aults [9]. In tis work, we only consider te callenging geometry o te ield, and its real permeability distribution. Gravity eects are included in tis model. Oter input speciications are sown in Table Results & Discussion We recall tat or te single rate case, te stability assumption is 1 M + c ϕ 0 > α2 and in te multirate case it reads 1 M + c ϕ 0 > α2 q + q. We consider a particular coice or q = 2 and or te parameters sown in Table 3, our assumption requires 1 M + c ϕ 0 > For te numerical test cases, we consider two dierent compressibility values corresponding to 1 satisying te stability condition and 2 te stability assumption is violated. In te irst case, we coose c = satisying te stability assumption. Te pressure proile ater 4010 simulation days is sown in igure 4.3a. Resulting pressures lie 19
21 Wells: 3 production wells, 6 injection well Injection well 1: Pressure speciied, psi Injection well 2: Pressure speciied, psi Injection well 3: Pressure speciied, psi Injection well 4: Pressure speciied, psi Injection well 5: Pressure speciied, psi Injection well 6: Pressure speciied, psi Production well 1: Pressure speciied, psi Production well 2: Pressure speciied, psi Production well 3: Pressure speciied, psi Total Simulation time: days Finer low Unit time step: 1.0 days Coarse mecanics time step: 2.0 days q = 2 Number o grids: 891 grids Permeabilities: k xx, k yy, k zz igly varying, range: 5.27E-10, 3.10E+3 md Initial porosity, ϕ 0: 0.2 Fluid viscosity, µ : 1.0 cp Initial pressure, p 0: psi Fluid compressibilities: Case 1, condition is satisied, c : 1.E-4 1/psi Case 2, condition is not satisied, c : 1.E-13 1/psi Case 3, condition is not satisied, c : 1.E-8 1/psi Rock compressibility: 1.E-6 1/psi Rock density: lb m/t 3 Initial luid density, ρ : lb m/t 3 Young s Modulus E: 1.E5 psi Possion Ratio, ν: 0.3 Biot s constant, α: 0.7 Biot Modulus, M: 1.0E16 psi Eν = 1+ν1 2ν psi Flow Boundary Conditions: no low boundary condition on all 6 boundaries Mecanics B.C.: X+ boundary EBCXX1: σ xx = σ n x = 10, 000psi, overburden pressure X- - boundary EBCXXN1: u = 0, zero displacement Y+ - boundary EBCYY1: u = 0, zero displacement Y- - boundary EBCYYN1: σ yy = σ n y = 2000psi Z+ - boundary EBCZZ1: u = 0, zero displacement Z- - boundary EBCZZN1: σ zz = σ n z = 1000psi Table 3: Input Parameters or Frio Field Model 20
22 a Pressure proile wen te compressibility o te luid satisies te derived stability condition c = Results are pysically correct, and lie between te expected range o values. b Pressure proile wen te compressibility o te luid does not satisy te derived stability condition c = Results are not pysically correct, as pressure values drop below zero. Figure 4.3: Pressure proiles o te multirate explicit coupling sceme q=2 or te Frio ield model. in te expected range o values, based on wells injection and production rates speciied in table 3. Next, we consider te case wen we coose c = , tat strongly violates te stability condition. In tis case, te coupling iteration did not converge, as a result o producing extremely ig pressure values in magnitudes, and tat, in turn, triggered te pore-volume values o grid blocks to exceed teir corresponding bulk-volume values, wic is pysically meaningless. To urter test te eect o compressibility, we increase te compressibility and coose c = still violating te stability condition. In tis case, te pressure proiles ater 4010 simulation days are sown or two compressibilities in igure 4.3b. It is clear rom te igure tat te pressure proiles are unpysical since te pressure values drop below zero. Given te values o te initial pressure, and wells injection and production rates speciied in table 3, tis is a non-pysical solution. 5 Conclusions and outlook In tis work, we considered single rate and multirate explicit coupling scemes or coupling low wit geomecanics in poro-elastic media. We derived stability criteria or bot multirate and single rate scemes and derived te assumptions on reservoir parameters or te stability to old. In addition, we perorm te numerical experiments were we compare te time savings in te explicit coupling scemes compared to te iterative ixed stress scemes. Te multirate iterative scemes ave been proven to be geometrically convergent. Our computational results sow tat, i te parameters satisy te stability condition, explicit coupling scemes reduce CPU run time eiciently as compared to iterative scemes. 21
23 Acknowledgements TA is unded by Saudi Aramco. We tank Paulo Zunino and Ivan Yotov or elpul discussions. KK would like to acknowledge te support o Statoil Akademia Grant Bergen. Te autors would like to acknowledge te CSM Industrial Ailiates program, DOE grant ER25617, and ConocoPillips grant UTA Reerences [1] T. Almani, A. H. Dogru, K. Kumar, G. Sing, and M. F. Weeler. Convergence o multirate iterative coupling o geomecanics wit low in a poroelastic medium. Saudi Aramco Journal o Tecnology, Spring [2] T. Almani, K. Kumar, A. H. Dogru, G. Sing, and M. F. Weeler. Convergence analysis o multirate ixed-stress split iterative scemes or coupling low wit geomecanics. Ices report 16-07, Institute or Computational Engineering and Sciences, Te University o Texas at Austin, Austin, Texas, [3] M. A. Biot. Consolidation settlement under a rectangular load distribution. J. Appl. Pys., 125: , [4] M. A. Biot. General teory o tree-dimensional consolidation. J. Appl. Pys., 122: , [5] L. Y. Cin, L. K. Tomas, J. E. Sylte, and R. G. Pierson. Iterative coupled analysis o geomecanics and luid low or rock compaction in reservoir simulation. Oil and Gas Science and Tecnology, 575: , [6] O. Coussy. Mecanics o Porous Continua. Wiley, West Sussex PO19 1UD, England, [7] X. Gai, R. H. Dean, M. F. Weeler, and R. Liu. Coupled geomecanical and reservoir modeling on parallel computers. In Te SPE Reservoir Simulation Symposium, Houston, Texas, Feb. 3-5, [8] V. Girault, M. F. Weeler, B. Ganis, and M. Mear. A lubrication racture model in a poro-elastic medium. Tecnical report, Te Institute or Computational Engineering and Sciences, Te University o Texas at Austin, [9] M. Juntunen and M. F. Weeler. Two-pase low in complicated geometries - modeling te rio data using improved computational meses. Computational Geosciences, 17: , [10] J. Kim, H. A. Tcelepi, and R. Juanes. Stability and convergence o sequential metods or coupled low and geomecanics: ixed-stress and ixed-strain splits. Comput. Metods Appl. Mec. Engrg., : ,
24 [11] J. Kim, H. A. Tcelepi, and R. Juanes. Stability, accuracy, and eiciency o sequential metods or coupled low and geomecanics. In Te SPE Reservoir Simulation Symposium, Houston, Texas, February 2-4, SPE [12] K. Kumar, T. Almani, G. Sing, and M. F. Weeler. Multirate undrained splitting or coupled low and geomecanics in porous media. In ENUMATH 2015 Proceedings. European Conerence on Numerical Matematics and Advanced Applications, September 14-18, submitted. [13] M. Mainguy and P. Longuemare. Coupling luid low and rock mecanics: ormulations o te partial coupling between reservoir and geomecanics simulators. Oil and Gas Science and Tecnology - Rev. IFP, 574: , [14] A. Mikelić, B. Wang, and M. F. Weeler. Numerical convergence study o iterative coupling or coupled low and geomecanics. Computational Geosciences, 18: , [15] A. Mikelić and M. F. Weeler. Convergence o iterative coupling or coupled low and geomecanics. Computational Geosciences, 17: , [16] P. Samier, A. Onaisi, and S. d. Gennaro. A practical iterative sceme or coupling geomecanics wit reservoir simulation. SPE Reservoir Evaluation and Engineering, 11: [17] A. Settari and F. M. Mourits. Coupling o geomecanics and reservoir simulation models. In Siriwardane and Zema, editors, Comp. Metods and Advances in Geomec., pages , Balkema, Rotterdam, [18] L. San, H. Zeng, and W. J. Layton. A decoupling metod wit dierent subdomain time steps or te nonstationary stokes-darcy model. Numerical Metods or Partial Dierential Equations, 292: , [19] G. Sing. Coupled Flow and Geomecanics Modeling or Fractured Poroelastic Reservoirs. PD tesis, Te University o Texas at Austin, Austin, Texas, [20] J. C. Small, J. R. Booker, and E. H. Davis. Elasto-plastic consolidation o soil. Int. J. Solids Struct., 126: , [21] K. V. Terzagi. Teoretical Soil Mecanics. Wiley, New York, [22] M. F. Weeler, G. Xue, and I. Yotov. A multipoint lux mixed inite element metod on distorted quadrilaterals and exaedra. Numerisce Matematik, 1211: , [23] M. F. Weeler and I. Yotov. A multipoint lux mixed inite element metod. SIAM Journal o Numerical Analysis, 44: , [24] X. Xiong. Partitioned metods or coupled luid low and turbulence low problems. PD tesis, University o Pittsburg, Pittsburg, Pennsylvania,
ICES 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 informationICES REPORT Convergence and Error Analysis of Fully Discrete Iterative Coupling Schemes for Coupling Flow with Geomechanics
ICES REPORT 16-4 October 016 Convergence and Error Analysis of Fully Discrete Iterative Coupling Scemes for Coupling Flow wit Geomecanics by Tameem Almani, Kundan Kumar, Mary F. Weeler Te Institute for
More informationICES REPORT February Tameem Almani, Kundan Kumar, Ali H. Dogru, Gurpreet Singh, Mary F. Wheeler
ICES REPORT 16-07 February 2016 Convergence Analysis of Multirate Fixed-Stress Split Iterative Scemes for Coupling Flow wit Geomecanics by Tameem Almani, Kundan Kumar, Ali H. Dogru, Gurpreet Sing, Mary
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 informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationHydraulic validation of the LHC cold mass heat exchanger tube.
Hydraulic validation o te LHC cold mass eat excanger tube. LHC Project Note 155 1998-07-22 (pilippe.provenaz@cern.c) Pilippe PROVENAZ / LHC-ACR Division Summary Te knowledge o te elium mass low vs. te
More informationA Two-field Finite Element Solver for Poroelasticity on Quadrilateral Meshes
A Two-field Finite Element Solver for Poroelasticity on Quadrilateral Meses Graam Harper 1, Jiangguo Liu 1, Simon Tavener 1, and Zuoran Wang 1 Colorado State University, Fort Collins, CO 80523, USA {arper,liu,tavener,wangz}@mat.colostate.edu
More informationContinuity and Differentiability
Continuity and Dierentiability Tis capter requires a good understanding o its. Te concepts o continuity and dierentiability are more or less obvious etensions o te concept o its. Section - INTRODUCTION
More informationOn One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations
Pure and Applied Matematics Journal 7; 6(5: 74 ttp://wwwsciencepublisinggroupcom/j/pamj doi: 648/jpamj765 ISSN: 6979 (Print; ISSN: 698 (Online On One Justiication on te Use o Hybrids or te Solution o First
More informationICES 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 informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
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 informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
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 informationJournal of Applied Science and Agriculture. The Effects Of Corrugated Geometry On Flow And Heat Transfer In Corrugated Channel Using Nanofluid
Journal o Applied Science and Agriculture, 9() February 04, Pages: 408-47 AENSI Journals Journal o Applied Science and Agriculture ISSN 86-9 Journal ome page: www.aensiweb.com/jasa/index.tml Te Eects O
More informationNumerical Analysis of the Double Porosity Consolidation Model
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real 21-25 septiembre 2009 (pp. 1 8) Numerical Analysis of te Double Porosity Consolidation Model N. Boal
More informationVariational Methods for Time Dependent Wave Propagation Problems
Variational Metods or Time Dependent Wave Propagation Problems Patrick Joly INRIA Rocquencourt, BP105 Le Cesnay, France (patrick.joly@inria.r) 1 Introduction Tere is an important need or numerical metods
More informationTHE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION
THE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION I.V.Vodolazskaya, Yu.Yu.Tarasevic Astrakan State University, a Tatiscev St., Astrakan, 41456, Russia e-mail: tarasevic@aspu.ru e ave proposed and
More informationA Two-Field Finite Element Solver for Poroelasticity on Quadrilateral Meshes
A Two-Field Finite Element Solver for Poroelasticity on Quadrilateral Meses Graam Harper, Jiangguo Liu, Simon Tavener, and Zuoran Wang (B) Colorado State University, Fort Collins, CO 80523, USA {arper,liu,tavener,wangz}@mat.colostate.edu
More informationPARTITIONED METHODS FOR COUPLED FLUID FLOW PROBLEMS
PARTITIONED METHODS FOR COUPLED FLUID FLOW PROBLEMS by Hoang A. Tran B.S. in Matematics and Computer Sciences, University o Natural Sciences, Ho Ci Min City, Vietnam, 006 France-Vietnam M.S. in Applied
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationOnline Appendix for Lerner Symmetry: A Modern Treatment
Online Appendix or Lerner Symmetry: A Modern Treatment Arnaud Costinot MIT Iván Werning MIT May 2018 Abstract Tis Appendix provides te proos o Teorem 1, Teorem 2, and Proposition 1. 1 Perect Competition
More informationAn Efficient Multigrid Solver for a Reformulated Version of the Poroelasticity System
An Efficient Multigrid Solver for a Reformulated Version of te Poroelasticity System F.J. Gaspar a F.J. Lisbona a, C.W. Oosterlee b P.N. Vabiscevic c a Departamento de Matemática Aplicada, University of
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationUNIT #6 EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS
Answer Key Name: Date: UNIT # EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS Part I Questions. Te epression 0 can be simpliied to () () 0 0. Wic o te ollowing is equivalent to () () 8 8? 8.
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
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 informationChapter 2 Limits and Continuity. Section 2.1 Rates of Change and Limits (pp ) Section Quick Review 2.1
Section. 6. (a) N(t) t (b) days: 6 guppies week: 7 guppies (c) Nt () t t t ln ln t ln ln ln t 8. 968 Tere will be guppies ater ln 8.968 days, or ater nearly 9 days. (d) Because it suggests te number o
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationENTROPY GENERATION IN RECTANGULAR DUCTS WITH NONUNIFORM TEMPERATURE ON THE CONTOUR
ENROPY GENERAION IN RECANGULAR UCS WIH NONUNIFORM EMPERAURE ON HE CONOUR Cesar A. Marcelino Matoso Instituto ecnológico de Aeronáutica cesar_matoso@yaoo.ca Ezio Castejon Garcia Instituto ecnológico de
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationHOMOGENIZATION OF IMMISCIBLE COMPRESSIBLE TWO-PHASE FLOW IN DOUBLE POROSITY MEDIA
Electronic Journal o Dierential Equations, Vol. 2016 (2016), No. 52, pp. 1 28. ISSN: 1072-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu tp ejde.mat.txstate.edu HOMOGENIZATION OF IMMISCIBLE
More informationEssential Microeconomics : EQUILIBRIUM AND EFFICIENCY WITH PRODUCTION Key ideas: Walrasian equilibrium, first and second welfare theorems
Essential Microeconomics -- 5.2: EQUILIBRIUM AND EFFICIENCY WITH PRODUCTION Key ideas: Walrasian equilibrium, irst and second welare teorems A general model 2 First welare Teorem 7 Second welare teorem
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationA New Projected Gradient Algorithm for Total Variation Based Image Denoising
Noname manuscript No. (will be inserted by te editor) A New Projected Gradient Algoritm or Total Variation Based Image Denoising Ming-Jun Lai Leopold Matamba Messi te date o receipt and acceptance sould
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationSolving Second Order Linear Dirichlet and Neumann Boundary Value Problems by Block Method
IAENG International Journal o Applied Matematics 43: IJAM_43 04 Solving Second Order Linear Diriclet and Neumann Boundar Value Problems b Block Metod Zanaria Abdul Majid Mod Mugti Hasni and Norazak Senu
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationDIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS
Journal o Applied Analysis Vol. 14, No. 2 2008, pp. 259 271 DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS B. BELAÏDI and A. EL FARISSI Received December 5, 2007 and,
More informationMATH1901 Differential Calculus (Advanced)
MATH1901 Dierential Calculus (Advanced) Capter 3: Functions Deinitions : A B A and B are sets assigns to eac element in A eactl one element in B A is te domain o te unction B is te codomain o te unction
More informationConductance from Transmission Probability
Conductance rom Transmission Probability Kelly Ceung Department o Pysics & Astronomy University o Britis Columbia Vancouver, BC. Canada, V6T1Z1 (Dated: November 5, 005). ntroduction For large conductors,
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationNonlinear Darcy fluid flow with deposition
Nonlinear Darcy fluid flow wit deposition V.J. Ervin Hyesuk Lee J. Ruiz-Ramírez Marc 23, 2016 Abstract In tis article we consider a model of a filtration process. Te process is modeled using te nonlinear
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
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 informationCALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY
CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationUNIFORMLY-STABLE FINITE ELEMENT METHODS FOR DARCY-STOKES-BRINKMAN MODELS *
Journal o Computational Matematics, Vol.26, No.3, 2008, 437 455. UNIFORMLY-SABLE FINIE ELEMEN MEHODS FOR DARCY-SOKES-BRINKMAN MODELS * Xiaoping Xie Yangtze Center o Matematics, and Scool o Matematics,
More informationRecent Advances in Time-Domain Maxwell s Equations in Metamaterials
Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University,
More informationUNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics
UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):
More informationAn operator splitting approach for the interaction between a fluid and a multilayered poroelastic structure
An operator splitting approac or te interaction between a luid and a multilayered poroelastic structure M. Bukač I. Yotov P. Zunino Abstract We develop a loosely coupled luid-structure interaction inite
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationTheoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump
TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing
More informationContinuous Hybrid Multistep Methods with Legendre Basis Function for Direct Treatment of Second Order Stiff ODEs
American Journal o Computational and Applied Matematics 06, 6(): 8-9 DOI: 0.59/j.ajcam.06060.0 Continuous Hybrid Multistep Metods wit Legendre Basis Function or Direct Treatment o Second Order Sti ODEs
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationEvaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir
American Journal of Environmental Engineering and Science 2016; 3(2): 68-74 ttp://www.aascit.org/journal/ajees ISSN: 2381-1153 (Print); ISSN: 2381-1161 (Online) Evaluation and Accurate Estimation from
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationObservability Analysis of Nonlinear Systems Using Pseudo-Linear Transformation
9t IFAC Symposium on Nonlinear Control Systems Toulouse, France, September 4-6, 2013 TC3.4 Observability Analysis o Nonlinear Systems Using Pseudo-Linear Transormation Yu Kawano Tosiyui Otsua Osaa University,
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
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 informationA symmetric mixed finite element method for nearly incompressible elasticity based on. biorthogonal system
A symmetric mixed finite element metod for nearly incompressible elasticity based on biortogonal systems Bisnu P. Lamicane and Ernst P. Stepan February 18, 2011 Abstract We present a symmetric version
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA
ROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA Y. EFENDIEV, J. GALVIS, S. KI KANG, AND R.D. LAZAROV Abstract. In tis paper, we study robust iterative solvers for
More information