ICES REPORT Multirate Undrained Splitting for Coupled Flow and Geomechanics in Porous Media
|
|
- Edith Floyd
- 6 years ago
- Views:
Transcription
1 ICES REPORT 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 Te University o Texas at Austin Austin, Texas Reerence: Tameem Almani, Kundan Kumar, Mary F. Weeler, "Multirate Undrained Splitting or Coupled Flow and Geomecanics in Porous Media," ICES REPORT 16-13, Te Institute or Computational Engineering and Sciences, Te University o Texas at Austin, May 2016.
2 Multirate Undrained Splitting or Coupled Flow and Geomecanics in Porous Media T. Almani 1, K. Kumar 2, 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 a multirate iterative sceme or te uasi-static Biot euations modelling te coupled low and geomecanics in a porous medium. Te iterative sceme is based on undrained splitting were te low and mecanics euations are decoupled wit te mecanics solve ollowed by te pressure solve. Te multirate sceme proposed ere uses dierent time steps or te two euations, tat is, uses low steps or eac coarse mecanics step and may be interpreted as using a regularization parameter or te mecanics euation. We prove te convergence o te sceme and te proo reveals te appropriate regularization parameter and also te eect o te number o low steps witin coarse mecanics step on te convergence rate. Keywords. poroelasticity; undrained split iterative coupling; multirate sceme; contraction mapping 1 Introduction Coupling o geomecanics and low in poroelastic media as many important applications suc as subsidence events, carbon seuestration, ground water remediation, ydrocarbon production, enanced geotermal systems, solid waste disposal, and biomedical modeling. Starting rom te pioneering work o Terzagi and Biot [1], tere as been active investigation into te coupled geomecanics and low problems [12]. Te Biot model consists o a geomecanics euation coupled to a low model wit te displacement, pressure and low velocity as unknowns. Tere is a uge literature on Biot euations and tey ave been analyzed by a number o autors wo establised existence, uniueness, and regularity, see Sowalter [14], Pillips and Weeler [9] and reerences terein. In contrast to a ully implicit sceme or solving te coupled model o Biot euations, iterative metods are oten employed in practice [2 4]. Te iterative scemes allow te decoupling o te low and mecanics euations and tus oer several attractive eatures suc as use o existing low and mecanics codes, use o appropriate pre-conditioners and solvers or te two models, and ease o implementation. Te design o iterative scemes owever is an important consideration or an eicient, convergent, and robust algoritm. 1
3 a Single Rate b Multirate Figure 1.1: Flowcart o te undrained split single rate and multirate iterative coupling algoritms. In te single rate sceme, bot low and mecanics sare te exact same time step t. In te multirate sceme, te low iner time step is t, and mecanics coarser time step is t. In addition, oten we can take a coarser time step or te mecanics euation tan or te low. Here, we consider one o te iterative scemes oten used in practice: undrained splitting and propose a multirate iterative sceme, as illustrated by Figure 1.1b. Tis sceme considers a iner time step or te low model and a coarser time step or mecanics low steps or eac mecanics step and ten perorms an iteration between te mecanics and iner low steps. In contrast, in te single rate sceme, as illustrated by Figure 1.1a, te low and mecanics problems sare te exact same time step, and te coupling iteration continues until convergence. Bot scemes are iterative in te sense tat or eac coarse mecanics time step, we solve or low iner time steps ollowed by a mecanics step and we urter repeat te process. Details about convergence criteria can be ound in [7]. Te converged solutions solve te coupled time-discrete system consisting o low solves and one mecanics solve. Te low iner solve uses te mecanics at te coarse step and ence, te coupled system is ully implicit. Since te cost o mecanics is oten muc more tan te low, a less mecanics solves leads to computational savings. Our work is motivated by te recent work o Mikelić and Weeler [7,8] were tey ave considered dierent iterative scemes or low and mecanics couplings and establised contractive results in suitable norms, see also [6] or studying te von Neumann stability o iterative algoritms, [13] or multirate scemes or Darcy-Stokes, and [10,11] or relationsip o tese iterative metods to te linearization procedures. 2
4 2 Model euations, discretization and splitting algoritm 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. In general, we assume tat te boundary is decomposed into Diriclet boundary Γ D, and Neumann boundary Γ N, associated wit Diriclet and Neumann boundary conditions respectively, suc tat Γ D Γ N = Ω. Our model assume a linear, elastic, omogeneous, and isotropic poro-elastic medium Ω R d, in wic te reservoir is saturated wit a sligtly compressible viscous luid. 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 uasi-static Biot model [1] will be employed in tis work. Te model reads: Find u and p satisying te euations below or all time t ]0, T [: Flow Euation: t 1 M + c ϕ 0 p + α u Mecanics Euations: div σ por u, p = in Ω, Boundary Cond.: 1 K p ρ,r g η = in Ω σ por u, p = σu α p I in Ω, σu = λ ui + 2 Gεu in Ω u = 0 on Ω, K p ρ,r g η n = 0 on Γ N, p = 0 on Γ D, Initial Cond. t=0: 1 M + c ϕ 0 p + α u 0 = 1 M + c ϕ 0 p0 + α u 0. were: g is te gravitational constant, η is te distance in te direction o gravity assumed to be constant in time, ρ,r > 0 is a constant reerence density relative to te reerence pressure p r, ϕ 0 is te initial porosity, M is te Biot constant, = ρ,r were 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 involving α. 2.1 Mixed variational ormulation We use a mixed inite element ormulation or low and a conormal Galerkin ormulation or mecanics or te spatial discretization and a backward-euler or te 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 3
5 velocities luxes: V = {v H 1 Ω d ; T T, v T P 1 d, v Ω = 0}, Q = {p L 2 Ω ; T T, p T P 0 }, Z = { Hdiv; Ω d ; T T, T P 1 d, n = 0 on Ω}. We also assume tat te iner time step is given by: t = t k t k 1. I we denote te total number o timesteps by N, ten te total simulation time is given by T = t N, and t i = i t, 0 i N denote te discrete time points. Te proo presented ere can be easily extended to oter mixed metod approaces e.g., [15] or Conormal Galerkin discretizations. Notation: 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, t stands or te unit iner time step, and is te ixed number o local low time steps per coarse mecanics time step. We start by analyzing te ully discrete single-rate undrained-split coupling algoritm, in wic bot low and mecanics sare te exact same time step. Tis can be viewed as an extension o te work o Mikelić & Weeler [8], in wic conormal Galerkin is used or low discretization versus mixed orm or te low in tis work, and te contraction coeicient obtained ere is optimized. In addition, understanding te strategy o te single-rate proo serves as a good introduction beore tackling te more complicated multirate case. 2.2 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. Deinition 2.1 mecanics euation Find u k+1 V suc tat, v V, 2G εu k+1, εv + λ u k+1, v α p k+1, v =, v. 2.1 and low euation ind p k+1 Q, and z k+1 Z suc tat,, θ θ Q, Z, 1 1 t α t M + c ϕ 0 p k+1 u k+1 K 1 z k+1, = p k, θ + u k p k+1, + wit te initial conditions or te irst discrete time step, + 1 z k+1, θ =, θ, 2.2 ρ,r g η,, 2.3 p 0 = p 0, u 0 = u
6 2.3 Fully discrete sceme or multirate Te weak ormulation o a multirate sceme reads as ollows. Deinition 2.2 mecanics euation Find u k+ V suc tat, v V, 2G εu k+ and low euation For 1 m, ind p m+k θ Q, Z, 1 1 t, εv + λ u k+, v α p k+, v =, v. 2.5 α t M + c ϕ 0 p m+k u k+ K 1 z m+k, = Q, and z m+k, θ p m 1+k, θ + u k p m+k, + wit te initial conditions or te irst discrete time step, Z suc tat, + 1 z m+k, θ =, θ, 2.6 ρ,r g η,, 2.7 p 0 = p 0, u 0 = u Single Rate Formulation - Undrained Split Iterative Metod Te sceme starts by solving te mecanics problem ollowed by te low problem, and iterates between te two until convergence is acieved. Te iteration assumes a constant luid mass during te deormation o te structure and can be interpreted as a regularization o mecanics euation. We start by presenting te sceme. In te above, we ave used β = 1 M + c ϕ 0 or te notational convenience. L is a regularization parameter and te corresponding term vanises in te case o convergence. 3.1 Analysis o te single rate undrained sceme Te convergence proo is based on studying te dierence o two successive iterates and deriving te contraction o appropriate uantities in suitable norms. Accordingly, we deine: δξ n,k = ξ n+1,k ξ n,k, were ξ = p, u, or z. Te uantity to be contracted on is a composite one consisting o bot pressure p n,k and volumetric strain terms u n,k. For a particular coupling iteration, n 1, and or time step t k, we deine te uantity to be contracted on as: m n,k = L γ un,k + α γ pn,k, were γ is an adjustable coeicient tat will be selected careully suc tat te sceme acieves contraction on m. Te presence o γ does not alter te contractivity, owever, it simpliies te algebra and provides a systematic tecniue or obtaining similar results or oter problems. 5
7 Algoritm 1: Undrained Split Single Rate Algoritm 1 or k = 0, 1, 2, 3,.. do /* mecanics time step iteration index */ 2 or n = 1, 2,.. do /* coupling iteration index */ 3 First Step: Mecanics euations 4 Given p n,k+1, solve or u n+1,k+1 satisying assuming an initial value is given or te irst iteration: p 0,k+1 : 2G εu n+1,k+1 λ + L u n+1,k+1 I = α Second Step: Flow euations 5 Given u n+1,k+1, solve or p n+1,k+1 and z n+1,k+1 p n+1,k+1 β t p k pn,k+1 I L u n,k+1 I z n+1,k+1 = α z n+1,k+1 satisying: u n+1,k+1 t u k = K p n+1,k+1 ρ,r g η Weak ormulation o dierence o two successive iterates Considering te dierence between two consecutive iterative coupling iterations, te weak ormulation o euations corresponding to 3.1, 3.2, and 3.3 can be written as ollows. Mecanics step: V suc tat, Given δp n,k+1 rom te previous coupling iteration, ind δu n+1,k+1 v V, 2G εδu n+1,k+1, εv + λ + L δu n+1,k+1, v = α δp n,k+1, v + L δu n,k+1, v, 3.4 Flow step: θ Q, Given δu n+1,k+1, ind δp n+1,k+1 Q, δz n+1,k+1 Z suc tat: β δp n+1,k+1 t, θ + 1 δz n+1,k+1 α, θ = t δu n+1,k+1, θ 3.5 Z, K 1 δz n+1,k+1, = δp n+1,k+1, 3.6 Step 1: Mecanics euation First, we analyze te mecanics euation. Testing 3.4 wit v = δu n+1,k+1, we get: 2G εδu n+1,k λ + L δu n+1,k+1 2 = αδp n,k+1 + L δu n,k+1, δu n+1,k+1 = γδ m n,k+1, δu n+1,k+1 ε δu n+1,k δ + 2 2ε γ2 m n,k+1 2 6
8 by Young s ineuality. For ε = λ + L, we obtain, 4G λ + L εδu n+1,k+1 Step 2: Flow euation 2 + δu n+1,k+1 2 γ 2 δ m n,k+1 λ + L Testing 3.5, wit θ = δp n+1,k+1, and multiplying by t, we get: recall β = 1 M + c ϕ 0 βδp n+1,k t δz n+1,k+1, δp n+1,k+1 = α δu n+1,k+1, δp n+1,k Testing 3.6 wit = δz n+1,k+1, we get and taking te dierence between tem, we get K 1/2 δz n+1,k+1 Substituting 3.9 into 3.8, we ave β δp n+1,k+1 2 = δp n+1,k t K 1/2 δz n+1,k α Step 3: Combining Mecanics and Flow, δz n+1,k δu n+1,k+1, δp n+1,k+1 = Multiplying 3.10 by anoter ree parameter c 2 and adding 3.7, we obtain 4G εδu n+1,k+1 λ + L { 2 + c 2 βδp n+1,k c 2 α δu n+1,k+1, δp n+1,k+1 + } δu n+1,k c2 t K 1/2 δz n+1,k+1 2 γ 2 δ m n,k+1 λ + L Step 4: Identiying te parameters Below we provide te procedure or determining te tree adjustable parameters c 2, γ, and L yielding a contraction. Tese parameters sould be cosen suc tat te terms on te let and side o 3.11 remain positive, and te sceme acieves contraction on m. Clearly, δ m n+1,k+1 2 L 2 = δu n+1,k+1 γ α2 δp n+1,k+1 γ αL γ 2 δp n+1,k+1, δu n+1,k+1. Matcing coeicients by comparing wit te terms in te curly brackets in 3.11 provides us wit te ollowing conditions: L 2 γ 2 = 1, α 2 γ 2 = c2 β, 2αL γ 2 = c 2 α. Tis gives, L = γ = α2 2β, c2 = 2 L. Substituting in 3.11 leads to a contraction actor o te orm L 2, λ+l wic is strictly less tan one. 7
9 3.2 Main result: Optimized Contraction or Single rate Our main result summarises te above contraction result. Teorem 3.1 Wit L = α2 2β and c2 = 4β, te undrained single rate iterative sceme deined α 2 by is a contraction given by 4G εδu n+1,k+1 λ + L 2 + c2 t K 1/2 δz n+1,k δ m n+1,k+1 2 L 2 δ m n,k+1 2. λ + L 3.12 Furtermore, te seuences deined by tis sceme converge to te uniue solution o te weak ormulation Remark 3.2 Te above contraction result implies tat te composite uantity m n+1,k+1, symmetric strain εu n+1,k+1, and lux z n+1,k+1 converge at a geometric rate. Relatively straigtorward arguments tat include standard mixed metod or controlling pressure by lux, Korn s ineuality to control te H 1 norm by te L 2 norm o te symmetric strain tensor, imply te uniue convergence o p n+1,k+1, u n+1,k+1 in L 2 and H 1 norms respectively. Remark 3.3 Te above remark 3.2 can also be used to strengten te contraction coeicient in teorem 3.1. Using triangle s ineuality, δ m n,k L γ δun,k + α γ δpn,k and using standard mixed metod to estimate te pressure by te lux and te Korn ineuality, we obtain δ m n,k C n,k δz εδu + n,k Now consider 3.12 in te teorem 3.1 and denote te irst two terms on te let by I n+1,k+1, tat is, I n+1,k+1 = 4G εδu n+1,k+1 λ + L 2 + c2 t K 1/2 δz n+1,k+1 2. Te above ineuality 3.13 can be rewritten as using a generic C, δ m n+1,k+1 2 C 2 I n+1,k Te ineuality 3.12 takes te orm C δ m n+1,k+1 2 L 2 δ m n,k+1 2 λ + L yielding an improved contraction constant, δ m n+1,k+1 2 C L 2 δ m n,k C + 1 λ + L In practice, it is diicult to estimate C, owever, te above computations sow te relative contributions o te extra positive terms in 3.12 aect te contraction result observed in practice. 8
10 Remark 3.4 Te single rate undrained split iterative coupling sceme as been rigorously analyzed by Mikelić and Weeler [8]. In teir analysis, a contraction on te luid mass per bulk volume, deined by: m = m 0 + ρ 0 α u + ρ 0 M p p 0, as been obtained or a continuous in time and space ormulation. Te value o te introduced ree parameter L as been cosen a priori to be L = Mα 2. Following teir approac, and adapting to our ully discrete ormulation mixed orm or low, and conormal Galerkin or mecanics, we ave te ollowing result: Teorem 3.5 [Mikelić & Weeler [8]] For L = Mα 2, ϕ n,k := ϕ 0 +α u n,k + 1 M pn,k p 0, and δ denoting te dierence o two successive iterates, te undrained split sceme as given in algoritm 1 is a contraction given by 4Gα 2 εδu n+1,k λ+mα t K 1/2 M δz n+1,k 2 + 2c ϕ 0 M δp n+1,k 2 + δϕ n+1,k+1 2 Mα 2 2 δϕ n,k+1 λ+mα 2 2. Furtermore, te seuences deined by tis sceme converge to te uniue solution o te weak ormulation We note tat our proo is optimized in te sense tat it reveals te optimal value o te parameter L = α2 2 2β, wic leads to a sarper contraction estimates: L λ+l = α 2 M 2. 2λ+2λMc ϕ 0 +α 2 M 9
11 4 Multirate Formulation - Undrained Split Iterative Metod Te sceme starts by solving te mecanics problem ollowed by a seuence o low problems, and iterates between te two until convergence is acieved. Te iteration assumes a constant luid mass during te deormation o te structure and can be interpreted as a regularization o mecanics euation. We start by presenting te sceme. Algoritm 2: Undrained Split Multirate Algoritm 1 or k = 0,, 2, 3,.. do /* mecanics time step iteration index */ 2 or n = 1, 2,.. do /* coupling iteration index */ 3 First Step: Mecanics euations 4 Given p n,k+, solve or u n+1,k+ satisying assuming an initial value is given or te irst iteration: p 0,k+ : 2G εu n+1,k+ λ + L u n+1,k+ I = α pn,k+ I L u n,k+ I Second Step: Flow euations 5 Given u n+1,k+ 6 or m = 1, 2,.., do /* low iner time steps iteration index */ 7 Solve or p n+1,m+k and z n+1,m+k satisying: p n+1,m+k β p n+1,m 1+k + 1 t z n+1,m+k z n+1,m+k α = u n+1,k+ u k t = K p n+1,m+k ρ,r g η 4.3 In te above, we ave used β = 1 M + c ϕ 0 or te notational convenience. L is a regularization parameter and te corresponding term vanises in te case o convergence. 4.1 Analysis o te multirate undrained sceme Following a similar approac to tat o te single rate case, te convergence proo is based on studying te dierence o two successive iterates and deriving te contraction o appropriate uantities in suitable norms. We recall tat δξ n,k = ξ n+1,k ξ n,k, were ξ = p, u, or z. It is interesting tat te contracting uantity is a composite one consisting o bot pressure p n+1,k+m and volumetric strain terms u n+1,k+. For a particular coupling iteration, n 1, and between two coarse mecanics time steps t k and t k+, we deine te uantity to be contracted on as: m n+1,k+m = L γ un+1,k+ + α γ pn+1,k+m p n+1,k+m 1, or 1 m, 10
12 were γ is an adjustable coeicient tat will be selected careully suc tat Banac contraction olds on m. As in te single-rate case, te contractivity o te sceme is not altered by te presence o γ, as it only simpliies te algebra, and scales te te uantity o contraction is a way suc tat te contraction coeicient is most optimized Weak ormulation o dierence o two successive iterates Considering te dierence between one local low iteration and its corresponding local low iteration in te previous coupling iteration, and te dierence between two consecutive mecanics coupling iterations, te weak ormulation o euations corresponding to 4.1, 4.2, and 4.3 can be written as ollows. Mecanics step: V suc tat, Flow step: tat: Given δp n,k+ rom te previous coupling iteration, ind δu n+1,k+ v V, 2G εδu n+1,k+, εv + λ + L δu n+1,k+, v = α δp n,k+, v + L δu n,k+, v, 4.4 Given δu n+1,k+, or 1 m, ind δp n+1,m+k Q, δz n+1,m+k Z suc δp n+1,m+k θ Q, β, θ t + 1 δz n+1,m+k, θ = α t δu n+1,k+, θ 4.5 Z, K 1 δz n+1,m+k, = δp n+1,m+k, 4.6 Step 1: Mecanics euation First, we analyze te mecanics euation. Testing 4.4 wit v = δu n+1,k+ by noting tat 2G εδu n+1,k+ = αδp n,k+ = ε λ + L δu n+1,k+ 2 + L δu n,k+, δu n+1,k+ α δp n,m+k δp n,m 1+k δu n+1,k ε γ2 δp n,m+k δp n,m 1+k L + δun,k+ δ m n,k+m ε = λ + L, we obtain ater some simpliications, 4G λ + L εδu n+1,k+ = δp n,k+ 2 + δu n+1,k+ 2 2, δu n+1,k+, we get: and using Young s ineuality. For γ 2 λ + L 2 δ m n,k+m
13 Step 2: Flow euation Testing 4.5, wit θ = δp n+1,m+k β = 1 M + c ϕ 0 β δp n+1,m+k α δu n+1,k+, δp n+1,m+k, and multiplying by t, we get: recall 2 + t δz n+1,m+k, δp n+1,m+k = 4.8 Now, consider 4.6 or two consecutive local low iner time steps, t = t m+k, and t = t m 1+k, and test wit = δz n+1,m+k and taking te dierence between tem, we get K 1 δz n+1,m+k δz n+1,m 1+k,δz n+1,m+k = δp n+1,m+k, δz n+1,m+k. 4.9 Substituting 4.9 into 4.8, we ave β δp n+1,m+k 2 + t K 1 δz n+1,m+k α δu n+1,k+ By Young s ineuality, wit urter simpliications, δz n+1,m 1+k, δp n+1,m+k βδp n+1,m+k 2 + α δu n+1,k+, δp n+1,m+k t K 1/2 δz n+1,m+k 2 2 K 1/2 δz n+1,m 1+k 2 + K 1/2 δz n+1,m+k, δz n+1,m+k =. δz n+1,m 1+k 2 = 0. Summing or local low time steps and ater some simpliications telescopic cancellations togeter wit te act tat δz n+1,k = 0, we get β δp n+1,m+k t K 1/2 δz n+1,k t α δu n+1,k+ K 1/2 δz n+1,m+k, δp n+1,m+k δz n+1,m 1+k 2 =
14 Step 3: Combining Mecanics and Flow Multiplying 4.10 by anoter ree parameter c 2 and adding 4.10, we obtain 4G λ + L + c2 α εδu n+1,k+ 2 + δu n+1,k+ { c 2 β, δp n+1,m+k δp n+1,m+k 2 + δu n+1,k+ 2} + c2 t K 1/2 δz n+1,k c2 t K 1/2 δz n+1,m+k 2 δz n+1,m 1+k 2 γ 2 λ + L 2 δ m n,k+m Step 4: Identiying te parameters Note tat we ave tree ree parameters: c 2, γ, and L. Below we provide te procedure or determining tese parameters yielding a contraction. Tese parameters sould be cosen suc tat te terms on te let and side o 4.11 remain positive, and te sceme acieves contraction on m. Clearly, δ m n+1,k+m 2 = L2 δu n+1,k+ 2 γ α2 p n+1,k+m γ 2 p n+1,k+m 1 + 2αL γ 2 p n+1,k+m p n+1,k+m 1, δu n+1,k+. Matcing coeicients by comparing wit te terms in te curly brackets in 4.11 provides L us 2 α = 1, 2 c 2 β, and 2αL = c2 α α2 2 γ 2 γ 2 γ 2. Tis gives, L = γ, L 2β and since te contraction actor is monotone wit respect to L, its minimum is acieved wen we coose, L = α2 2β implying γ = α2 2β and c2 = 42 β α 2. 2 L Using above in 4.11 we note tat te contraction actor is 2 and is smaller wen 2 λ+l 2 is larger. Also, wen = 1, te above contraction rate reduces to tat o te single rate case [8] wen te time steps or te mecanics and low are te same. 4.2 Main result: Contraction or Multirate Our main result summarises te above contraction result. Teorem 4.1 Wit L = α2 2β and c2 = 42 β, te undrained multirate iterative sceme α 2 13
15 deined by is a contraction given by c 2 t K 1/2 δz n+1,k c2 t 2 + δ m n+1,k+m 2 + 4G λ + L εδu n+1,k+ K 1/2 δz n+1,m+k 2 L 2 2 λ + L 2 δz n+1,m 1+k 2 δ m n,k+m Furtermore, te seuences deined by tis sceme converge to te uniue solution o te weak ormulation Convergence to discrete multirate ormulation From te result obtained above, we establis convergence o te seuences generated by te multirate undrained split algoritm and sow tat te converged uantities satisy te weak ormulation Te proo uses te matematical induction or te iner low euations combined wit te contraction estimates obtained above. Lemma 4.2 For every coarser mecanics time step, t = t k, tere exist a limit unction u k suc tat u n,k u k strongly in H 1 Ω d. Proo. Te contraction result in teorem 4.1 implies tat or a coarser time step t = t k, εδu n+1,k converges geometrically to zero. Using Korn s ineuality, tis implies tat xi u n+1,k, i = 1, 2, 3 is a Caucy seuence converging geometrically to a uniue limit in L 2 Ω d. It ollows immediately tat u n+1,k is a Caucy seuence converging geometrically to a uniue limit in H 1 Ω d, being a Hilbert space. Lemma 4.3 For every two consecutive coarser mecanics time steps, t = t k, and t = t k+, and or every 1 m, tere exist limit unctions p m+k, z m+k suc tat p n,m+k p m+k in L 2 Ω, z n,m+k z m+k in Hdiv, Ω d, wit strong convergence in te norms o te above spaces. Proo. Te contraction result in teorem 4.1 implies tat te uantities K 1/2 δz n+1,m+k δz n+1,m 1+k 2, and δ m n+1,m+k 2 converge geometrically to zero. It ollows tat or 1 m, K 1/2 δz n+1,m+k δz n+1,m 1+k 2, and δ m n+1,m+k 2 converge ge- ometrically to zero. Moreover, by 4.3, and Poincaré ineuality, K 1/2 δp n+1,m+k δp n+1,m 1+k 2 and δp n+1,m+k 2 converge geometrically to zero, respectively. Tis implies tat or every 1 m, te iner time step dierences p n,m+k 14
16 p n,m 1+k, z n,m+k z n,m 1+k, and te uantity deined by m n,m+k are Caucy seuences in L 2 Ω. We will sow strong convergence o te pressure seuence by induction. Te proo o strong convergence o te lux seuence ollows in te same way. Given an initial pressure value or t = t 0 : p n,0 = p 0, rom te above discussion, p n,1 p 0 is a Caucy seuence in L 2 Ω, and, in turn, p n,1 is a Caucy seuence in te complete space L 2 Ω, and as a uniue limit. Tis completes te base case or induction. For te inductive ypotesis, we assume tat or any coarser mecanics time step t = t k, and or any 1 m, p n,m+k is a Caucy seuence converging to a uniue limit in L 2 Ω: p n,m+k p m+k in L 2 Ω. We will sow tat p n,m+k+2 is also a Caucy seuence converging to a uniue limit in L 2 Ω. Consider te two Caucy seuences in L 2 Ω: p n,m+k+2 p n,m+k+1 and p n,m+k+1 p n,m+k. Let or a, b L 2 Ω. It ollows tat Tus, p n,m+k+2 p n,m+k+2 p n,m+k+1 p n,m+k+1 p n,m+k+2 p n,m+k+1 a in L 2 Ω, p n,m+k b in L 2 Ω, + p n,m+k+1 p n,m+k a + b + p m+k a + b in L 2 Ω. in L 2 Ω, by te inductive ypotesis. Tis completes te inductive step. Tereore, we obtain tat or all coarser mecanics time steps t = t k, and or 1 m, p n,m+k, z n,m+k are Caucy seuences converging geometrically to uniue limits in L 2 Ω. For te divergence o te lux, we note tat 4.5 amounts to te ollowing euality a.e. in L 2 Ω: δz n+1,m+k Te convergence o z n+1,m+k p n+1,m+k p n+1,m+k p n+1,m+k p n+1,m+k = β t δpn+1,m+k α t δun+1,k+. and m n+1,m+k and m n+1,m+k in L 2 Ω ollows rom te convergence o te dierence in L 2 Ω, establised above te convergence o implies te convergence o u n+1,k+ by deini- and z n,m+k converging geometrically to uniue limits tion. Tus, we ave bot z n,m+k in L 2 Ω, and ence z n+1,k converges to a uniue limit in Hdiv, Ω d. It remains to pass to te limit in Tis is straigtorward since te euations are linear and all operators involved are continuous in te spaces invoked in te statements o Lemmas 4.2 and 4.3. Moreover te convergences are strong. Tereore, we easily retrieve te ully discrete multirate ormulation. Remark 4.4 As in te single rate case discussed in remark 3.3, in te multirate case too te contraction in teorem 4.1 can be improved. Te above proo already provides te arguments 15
17 reuired. Note tat using triangle s ineuality, δ m n+1,k+m L γ δun+1,k+ + α γ δpn+1,k+m δp n+1,k+m 1, and using standard mixed metod to estimate te pressure by te lux and te Korn ineuality, we obtain δ m n+1,k+m 2 C εδu n+1,k+ 2 + C εδu n+1,k+ 2 + K 1/2 δz n+1,m+k K 1/2 δz n+1,m+k Now consider 3.12 in te teorem 4.1 and denote I n+1,k+ I n+1,k+ = 4G λ + L εδu n+1,k+ 2 + c2 t 2 δz n+1,m 1+k 2 δz n+1,m 1+k 2 + K 1/2 δz n+1,k K 1/2 δz n+1,m+k + c2 t K 1/2 δz n+1,k Te above ineuality 4.13 can be rewritten as using a generic C, Te ineuality 4.12 takes te orm C δz n+1,m 1+k 2 δ m n+1,k+m 2 C 2 I n+1,k δ m n+1,k+m 2 L 2 δ m n,k+m λ + L yielding an improved contraction constant, δ m n+1,k+1 2 C L 2 δ m n,k C + 1 λ + L In practice, it is diicult to estimate C, owever, te above computations sow te relative contributions o te extra positive terms in 3.12 aect te contraction result observed in practice. 5 Conclusions and outlook We ave considered single rate and multirate iterative coupling scemes or te seuential coupling o low wit mecanics based on te undrained split iterative coupling algoritm. 2 16
18 For bot scemes, we ave proved Banac ixed-point contraction, and convergence to te weak solution o te corresponding ully discrete sceme ollows immediately. Te multirate is a natural extension o te single rate sceme, and contracts on a composite uantity consisting o pressure and volumetric strain terms. Contraction proos are optimal in te sense tat contraction uantities are scaled suc tat more terms on te let and side are absorbed. Compared to previously obtained results [8], our derived contraction coeicients are saper. To te best o our knowledge, tis is te irst time a contraction result as been rigorously obtained or te multirate undrained split iterative coupling sceme. However, It sould be noted tat our analysis limits to one coarser time step and we ave not considered and investigated te propagation o error due to spatial and temporal discretizations. Tese error estimates providing te convergence rate can be studied and analyzed, or example, in te spirit o [5]. Furter, te nonlinear extensions o tese scemes, teir matematical analyses and computational perormance are interesting uestions tat will be addressed in uture work. Moreover, te perormance o tese algoritms sould be investigated numerically, and appropriate convergence stopping criteria or sould be devised accordingly. 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] M. A. Biot. General teory o tree-dimensional consolidation. J. Appl. Pys., 122: , [2] 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, [3] X. Gai, S. Sun, M. F. Weeler, and H. Klie. A timestepping sceme or coupled reservoir low and geomecanics on nonmatcing grids. In SPE Annual Tecnical Conerence and Exibition, SPE [4] V. Girault, K. Kumar, and M. F. Weeler. Convergence o iterative coupling o geomecanics wit low in a ractured poroelastic medium. Ices report 15-05, Institute or Computational Engineering and Sciences, Te University o Texas at Austin, Austin, Texas, [5] 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,
19 [6] J. Kim, H. A. Tcelepi, and R. Juanes. Stability, accuracy, and eiciency o seuential metods or coupled low and geomecanics. In Te SPE Reservoir Simulation Symposium, Houston, Texas, February 2-4, SPE [7] A. Mikelić, B. Wang, and M. F. Weeler. Numerical convergence study o iterative coupling or coupled low and geomecanics. Computational Geosciences, 18: , [8] A. Mikelić and M. F. Weeler. Convergence o iterative coupling or coupled low and geomecanics. Computational Geosciences, 17: , [9] P. J. Pillips and M. F. Weeler. A coupling o mixed and continuous Galerkin inite element metods or poroelasticity. I. Te continuous in time case. Comput. Geosci., 112: , [10] I. S. Pop, F. Radu, and P. Knabner. Mixed inite elements or te ricards euation: Linearization procedure. Journal o Computational and Applied Matematics, : , July [11] F. A. Radu, J. M. Nordbotten, I. S. Pop, and K. Kumar. A robust linearization sceme or inite volume based discretizations or simulation o two-pase low in porous media. Journal o Computational and Applied Matematics, 289: , [12] 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, [13] Li San, Haibiao Zeng, and William J. Layton. A decoupling metod wit dierent subdomain time steps or te nonstationary stokes?darcy model. Numerical Metods or Partial Dierential Euations, 292: , [14] R. E. Sowalter. Diusion in poro-elastic media. J. Mat. Anal. Appl., 2511: , [15] M. F. Weeler and I. Yotov. A multipoint lux mixed inite element metod. SIAM Journal o Numerical Analysis, 44: ,
ICES 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationE E B B. over U is a full subcategory of the fiber of C,D over U. Given [B, B ], and h=θ over V, the Cartesian arrow M=f
ppendix : ibrational teory o L euivalence E E onsider ibrations P, P, and te category Fib[,] o all maps E @E 2 =(, ):@ :: P { P 2 () @ o ibrations; Fib[,] is a ull subcategory o [,] ; see [3]. Fib[,] is
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 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 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 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 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 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 informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
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 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 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 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 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 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 informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
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 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 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 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] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
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 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 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 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 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 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 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 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 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 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 informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
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 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 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 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 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 information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
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 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 informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
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 informationAnalysis of time-dependent Navier-Stokes flow coupled with Darcy
Analysis of time-dependent Navier-Stokes flow coupled wit Darcy flow Ayçıl Çeşmelioğlu and Béatrice Rivière Abstract Tis paper formulates and analyzes a weak solution to te coupling of time-dependent Navier-Stokes
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 information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
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 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 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 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 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 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 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 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 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 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 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 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 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 informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
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 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 informationFinite Difference Formulae for Unequal Sub- Intervals Using Lagrange s Interpolation Formula
Int. Journal o Mat. Analyi, Vol., 9, no. 7, 85-87 Finite Dierence Formulae or Unequal Sub- Interval Uing Lagrange Interpolation Formula Aok K. Sing a and B. S. Badauria b Department o Matematic, Faculty
More informationDifferentiation Rules c 2002 Donald Kreider and Dwight Lahr
Dierentiation Rules c 00 Donal Kreier an Dwigt Lar Te Power Rule is an example o a ierentiation rule. For unctions o te orm x r, were r is a constant real number, we can simply write own te erivative rater
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 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 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 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 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 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 informationREGULARITY OF SOLUTIONS OF PARTIAL NEUTRAL FUNCTIONAL DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY
Proyecciones Vol. 21, N o 1, pp. 65-95, May 22. Universidad Católica del Norte Antofagasta - Cile REGULARITY OF SOLUTIONS OF PARTIAL NEUTRAL FUNCTIONAL DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY EDUARDO
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 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 informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationAn Implicit Finite Element Method for the Landau-Lifshitz-Gilbert Equation with Exchange and Magnetostriction
SMA An Implicit Finite Element Metod for te Landau-Lifsitz-Gilbert Equation wit Excange and Magnetostriction Jonatan ROCHAT under te supervision of L. Banas and A. Abdulle 2 ACKNOWLEDGEMENTS I would like
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
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 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 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 information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
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 informationDifferentiation. introduction to limits
9 9A Introduction to limits 9B Limits o discontinuous, rational and brid unctions 9C Dierentiation using i rst principles 9D Finding derivatives b rule 9E Antidierentiation 9F Deriving te original unction
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationKrazy Katt, the mechanical cat
Krazy Katt, te mecanical cat Te cat rigting relex is a cat's innate ability to orient itsel as it alls in order to land on its eet. Te rigting relex begins to appear at 3 4 weeks o age, and is perected
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More information