Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
|
|
- Horace Leonard
- 5 years ago
- Views:
Transcription
1 Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh KAUST WEP Workshop January 3 February, 2 Joint work with Vivette Girault, Paris VI, and Danail Vassilev, University of Pittsburgh Department of Mathematics, University of Pittsburgh
2 Outline Mathematical model for the coupled flow problem Interface conditions Existence and uniqueness for a global weak formulation Equivalence to a domain decomposition weak formulation Finite element discretizations Conforming Stokes elements and mixed finite elements for Darcy Mortar finite elements on all subdomain interfaces Discrete inf-sup condition Existence and uniqueness of a discrete solution Convergence analysis Non-overlapping domain decomposition - reduction to a mortar interface problem Coupling of Stokes-Darcy flow with transport Department of Mathematics, University of Pittsburgh 2
3 Applications of coupled Stokes and Darcy flows Interactions between surface water and groundwater flows Flows in fractured porous media Flows through vuggy rocks Flows through industrial filters Fuel cells Blood flows Glaucoma Department of Mathematics, University of Pittsburgh 3
4 Model problem Γ Γ Ω = fluid region (Stokes flow) n 2 n Γ Ω = saturated porous Γ medium (Darcy s law) Γ I Γ Γ 2 u j : Ω j R d, fluid velocity in Ω j, p j : Ω j R, fluid pressure in Ω j. Deformation tensor D and stress tensor T in Ω s : D(u s ) := 2 ( u s + ( u s ) T ), T(u s, p s ) := p s I + 2µD(u s ). Department of Mathematics, University of Pittsburgh 4
5 Flow equations Stokes flow on Ω s : div T(u s, p s ) 2µ div D(u s ) + p s = f s in Ω (conservation of momentum), div u s = in Ω s (conservation of mass), u s = on Γ s (no slip). Darcy flow on Ω d : µk u d + p d = f d in Ω d (Darcy s law), div u d = q d in Ω d (conservation of mass), u d n d = on Γ d (no flow), Solvability condition: Ω d q d dx = Department of Mathematics, University of Pittsburgh 5
6 Mass conservation across Γ sd : Continuity of normal stress on Γ d : Interface conditions u s n s + u d n d = on Γ sd. n s T n s p s 2µn s D(u s ) n s = p d on Γ sd. Slip with friction interface condition: (Beavers-Joseph (967),Saffman (97), Jones (973), Jäger and Mikelić (2)) n s T τ j 2µn s D(u s ) τ j = µα Kj u τ j, j =, d, on Γ sd, where K j = τ j K τ j. Department of Mathematics, University of Pittsburgh 6
7 Global variational formulation For u d, v d in L 2 (Ω d ) d and u s, v s in H (Ω s ) d, a(u, v) = µ K u d v d +2 µ Ω d D(u s ) : D(v s )+ Ω s d j= For v d H(div; Ω d ), v s H(div; Ω s ) and q L 2 (Ω), b(v, q) = q div v d Ω d q div v s. Ω s µα (u s τ j )(v s τ j ). Γ sd Kj = {v H(div; Ω) ; v s H (Ω d ) d, v Γs =, (v n) Γd = }, v, v = ( v 2 H(div;Ω) + v s 2 H (Ω s )) /2. Find (u, p) L 2 (Ω) such that v, a(u, v) + b(v, p) = r L 2 (Ω), b(u, r) = Ω Ω f v, r q d. Department of Mathematics, University of Pittsburgh 7
8 Existence and uniqueness of a weak solution Lemma: q L 2 (Ω), sup v b(v, q) v β q L 2 (Ω) Lemma: a(v, v) γ v 2 v V = {v : div = } Proof: By Korn s inequality, a(v, v) C( v s 2 H (Ω s ) + v d 2 L 2 (Ω d ) ). Coercivity now follows from v V, v s L 2 (Ω s ) C ( v s 2 H (Ω s ) + v d 2 /2. L 2 (Ω d )) Lemma: The variational problem has a unique solution. Department of Mathematics, University of Pittsburgh 8
9 Domain decomposition variational formulation Let Ω s = Ω s,i, Ω d = Ω d,i. Interface conditions Stokes-Stokes interfaces: Darcy-Darcy interfaces: [v s ] =, [T n] = on Γ ss [u d n] =, [p d ] = on Γ dd = {v L 2 (Ω) d : v Ωs,i H (Ω s,i ) d, v Ωd,i H(div, Ω d,i ), v n L r ( Ω d,i )(r > ), v Γs =, (v n) Γd = }, W = L 2 (Ω), Λ = {ξ : ξ Γij H /2 (Γ ij ) Γ ij Γ ss, ξ Γij H /2 (Γ ij ) Γ ij Γ dd Γ sd }. Department of Mathematics, University of Pittsburgh 9
10 Domain decomposition variational formulation, cont. a s,i (u s,i, v s,i ) = 2 µ D(u s,i ) : D(v s,i )+ Ω s,i a d,i (u d,i, v d,i ) = µ K u d,i v d,i, Ω d,i d j= Ω s,i Γ sd µα Kj (u s,i τ j )(v s,i τ j ) b i (v i, w i ) = w i div v i Ω i a(, ) = a s,i (, ) + a d,i (, ), b(, ) = b i (, ) b Γ (v, λ) = [v] λ + [v n] λ + [v n] λ Γ ss Γ dd Γ sd Find u, p W, λ Λ: v, a(u, v) + b(v, p) + b Γ (v, λ) = w W, b(u, w) = w q d, Ω ξ Λ, b Γ (u, ξ) =. Ω f v, Department of Mathematics, University of Pittsburgh
11 Equivalence Lemma: The two variational formulations are equivalent. Department of Mathematics, University of Pittsburgh
12 Partition T h i on Ω i ; T h i Finite element discretization and T h j need not match at Γ ij. Stokes elements h s,i W h in Ω s,i : MINI (Arnold-Brezzi-Fortin), Taylor- Hood, Bernardi-Raugel; contain at least polynomials of degree r s and r s resp. Pressure Velocity Mixed finite elements d,i h W d,i h in Ω d,i : RT, BDM, BDFM, BDDF; contain at least polynomials of degree r d velocity pressure Lagrange multplier T H ij h := h i, W h := W h i L 2 (Ω) - partition of Γ ij, possibly different from the traces of T h i and T h j Λ H ij : continuous or discontinuous piecewise polynomials of degree at least m s on Γ ss or m d on Γ dd and Γ sd Λ H := Λ H ij Nonconforming approximation: Λ h Λ Department of Mathematics, University of Pittsburgh 2
13 Mortar finite element method Find u h h, p h W h, λ H Λ H : v h h, a(u h, v h ) + b(v h, p h ) + b Γ (v h, λ H ) = w h W h, b(u h, w h ) = w h q d, Equivalently, letting ξ H Λ H, b Γ (u h, ξ H ) =. Ω Ω f v h, V h = {v h h : b Γ (v h, ξ H ) = ξ H Λ H }, Find u h V h, p h W h : v h V h, a(u h, v h ) + b(v h, p h ) = f v h, w h W h, b(u h, w h ) = Ω Ω w h q d. Department of Mathematics, University of Pittsburgh 3
14 Mortar compatibility conditions On Γ ij Γ dd Γ sd, i < j (assume that Ω i is a Darcy domain), ξ H Λ H, sup ϕ h i V h i n ϕ h i, ξh Γij ϕ h i L 2 (Γ ij ) β d ξ H L 2 (Γ ij ). On Γ ij Γ ss, i < j, ξ H Λ H, sup ϕ h i V h i Γ ij H /2 (Γ ij) d ϕ h i, ξh Γij ϕ h i H /2 (Γ ij ) β s ξ H H /2 (Γ ij ). These are much more general than the mortar conditions used in [Bernardi-Maday- Patera], [Ben Belgacem] for Laplace and Stokes, and in [Layton-Schieweck-.], [Riviere-.], [Galvis-Sarkis] for Stokes-Darcy. The above condition on Γ dd is similar to the one used in [Arbogast-Cowsar-Wheeler-.] and [Arbogast-Pencheva-Wheeler-.]. Department of Mathematics, University of Pittsburgh 4
15 Π h s v H (Ω s,i) C v H (Ω s,i) Π h d v H(div;Ωs,i) C v H (Ω s,i) Weakly continuous interpolants V h s = {v h h s : [v h ], ξ H Γss = ξ H Λ H } V h d = {v h h d : [v h n], ξ H Γdd = ξ H Λ H } There exisit Π h s : H (Ω s ) V h d such that (div(v Π h sv), w h ) Ωs,i = w h W h, Ω s,i Ω s, There exist Π h d : H (Ω d ) V h s such that (div(v Π h dv), w h ) Ωd,i = w h W h, Ω d,i Ω d, Department of Mathematics, University of Pittsburgh 5
16 Discrete inf-sup condition Lemma: w h W h, sup b(vh, wh ) v h V h v Proof: Given w h W h, construct v H (Ω): β w W div v = w h, v H (Ω) C w h L 2 (Ω). Then construct an interpolant Π h : H (Ω) V h such that b(π h v v, q h ) = q h W h, Π h v C v H (Ω). Department of Mathematics, University of Pittsburgh 6
17 Existence and uniqueness of a discrete solution Lemma: There exists a unique solution to: Find u h V h, p h W h : v h V h, a(u h, v h ) + b(v h, p h ) = (f, v h ), w h W h, b(u h, w h ) = (q d, w h ). Proof: Coercivity: The discrete Korn s inequality gives which, combined with a s,i (v h, v h ) α s,i v h 2 H (Ω s,i ), v h V h, v h 2 L 2 (Ω s,i ) v h 2 H (Ω s,i ), implies as,i(v h, v h ) α s v h 2 H (Ω s,i ). On Ω d, ad,i (v h, v h ) α d v h 2 H(div;Ω d,i) vh : (div v h, w h ) = w h W h. Coercivity and inf-sup condition imply existence and uniqueness of a solution. Department of Mathematics, University of Pittsburgh 7
18 Approximation properties of V h Construct an interpolant I h : H (Ω) V h, I h v = (I h s v, I h d v), v I h s v H (Ω s,i) C h r s v H rs+ (Ω s,i), v I h d v H(div;Ωs,i) C (h r d+ v H rd + (Ωd,i) + hr s v H rs+ (Ω s,i)). Department of Mathematics, University of Pittsburgh 8
19 Convergence analysis Theorem: Proof: u u h + p p h W C(h r s + h r d+ + H m s+/2 + H m d+/2 ) u u h + p p h W C( inf v h V h u vh + where R h a(u, = sup vh ) + b(vh, p) (f, vh ) v h V h v h Bound on the non-conforming error: Θ(v h ) := a(u, v h ) + b(v h, p) (f, v h ) = b Γ (v h, λ) inf p w h W ) + R h, w h W h [v h n], λ Γsd CH md+ v h [v h n], λ Γdd CH md+/2 v h [v h ], λ Γss CH ms+/2 v h Department of Mathematics, University of Pittsburgh 9
20 Example : Smooth solution Ω = Ω Ω 2, where Ω = [, ] [ 2, ] and Ω 2 = [, ] [, 2 ] Taylor-Hood on triangles in Ω ; lowest order Raviart-Thomas on rectangules in Ω 2. u = u 2 = [ sin( x G + ω)ey/g cos( x G + ω)ey/g ], p = ( G K µ G )cos( x G + ω)e/(2g) + y.5, p 2 = G K cos( x G + ω)ey/g, where µ =., K =, α =.5, G = µk, and ω =.5. α Department of Mathematics, University of Pittsburgh 2
21 Convergence rates for Example mesh u u,h,ω rate p p,h,ω rate 4x4 9.74e e-3 8x8 2.4e e x6 6.5e e x32.59e e x64 4.5e e-5.9 Table : Numerical errors and convergence rates in Ω. mesh u 2 u 2,h,Ω2 rate p p 2,h,Ω2 rate 4x4 3.8e-2 8.e-3 8x8 9.2e e-3.7 6x6 2.48e e x e e x64.72e e-5.99 Table 2: Numerical errors and convergence rates in Ω 2. Department of Mathematics, University of Pittsburgh 2
22 Example 2: discontinuous tangential velocity u = u 2 = [ (2 x)(.5 y)(y ξ) y3 3 + y2 2 (ξ +.5).5ξy.5 + sin(ωx) ], [ ω cos(ωx)y χ(y +.5) + sin(ωx) p = sin(ωx) + χ + µ(.5 ξ) + cos(πy), 2K p 2 = χ (y +.5) 2 sin(ωx)y K 2 K, ], where ω = 6 and the other parameters are defined as in Example. Department of Mathematics, University of Pittsburgh 22
23 Computed pressure and velocity for Example 2 U V Left: horizontal velocity; Right: vertical velocity Department of Mathematics, University of Pittsburgh 23
24 Convergence rates for Example 2 mesh u u,h,ω rate p p,h,ω rate 4x4 3.54e- 3.e-2 8x8 8.6e e x6 2.5e e x e e x64.4e e-4.99 Table 3: Numerical errors and convergence rates in Ω. mesh u 2 u 2,h,Ω2 rate p p 2,h,Ω2 rate 4x4 2.6e-.8e- 8x8 5.79e e x6.47e e x32 3.7e e x e e-4 2. Table 4: Numerical errors and convergence rates in Ω 2. Department of Mathematics, University of Pittsburgh 24
25 Solution of the algebraic flow system A B L B T L T u p λ = f f 2 ( R L L T ) ( x λ ) = ( f ) x - subdomain unknowns; λ - interface unknowns Matrix is symmetric but indefinite. Form the Shur complement system: L T R Lλ = L T R f () Iterative method for () (e.g. CG). Each iteration requires evaluating R = R... R n, i.e., solving subdomain problems. Advantages: subdomain solves in parallel; reuse existing Stokes and Darcy codes Department of Mathematics, University of Pittsburgh 25
26 Domain decomposition Subdomain problems with specified interface normal stress (R Lλ): find (u h,i (λ), p h,i (λ)) Vi h W h i such that a i (u h,i(λ), v i ) + b i (v i, p h,i(λ)) = λ v i n i, v i V i h, Γ I b i (u h,i(λ), w i ) =, w i Wh. i λ Complementary subdomain problems (R f): find (ū h, p h ) Vh i W h i such that BC a i (ū h, v i ) + b i (v i, p h ) = f i v i, v i V i BC h, BC q Ω i q 2 BC b i (ū h, w i ) = q i w i, w i Wh. i BC BC Ω i Interface problem: find λ h Λ h such that s h (λ h, µ) b I (u h(λ h ), µ) = b I (ū h, µ), µ Λ h, Recover global velocity and pressure: u h = u h (λ h) + ū h, p h = p h (λ h) + p h. Department of Mathematics, University of Pittsburgh 26
27 Subdomain problems Both types of local problems are well posed due to the local discrete inf-sup conditions. Neumann-Robin BC on Γ I for the Stokes region: n T n = λ, n T τ j µα Kj u τ j =. Dirichlet BC on Γ I for the Darcy region: p 2 = λ. Interface (Steklov Poincaré) operator: S h : Λ h Λ h, (S h λ, µ) = s h (λ, µ) λ, µ Λ h ; S h : λ [u h(λ) n] Interface problem: Find λ h Λ h such that S h λ h = g h. Department of Mathematics, University of Pittsburgh 27
28 Condition number estimate Lemma: For all λ M h = {µ M h : Γ I µ = }, C, h λ 2 Γ I s (λ, λ) C,2 λ 2 Γ I (S : H /2 (Γ I ) H /2 (Γ I )) C 2, λ 2 Γ I s 2 (λ, λ) C 2,2 h λ 2 Γ I (S 2 : H /2 (Γ I ) H /2 (Γ I )) Theorem: cond(s h ) = O(h ) Proof: (C, h + C 2, ) λ 2 Γ I s(λ, λ) (C,2 + C 2,2 h ) λ 2 Γ, Department of Mathematics, University of Pittsburgh 28
29 Algorithm on non-matching grids Apply the Conjugate Gradient method for S h λ h = g h. Computing the action of the operator (needed at each CG iteration): Given mortar data λ h Λ h, project onto subdomain grids: λ h Q hi λ h Solve local Stokes and Darcy problems in parallel with boundary data Q hi λ h Project fluxes onto the mortar space and compute the jump: u h,i n i Q T h i u h,i n i, S h λ h = Q T h u h, n + Q T h 2 u h,2 n 2 Department of Mathematics, University of Pittsburgh 29
30 Multiple subdomains Three types of interfaces: Stokes Stokes, Stokes Darcy, and Darcy Darcy. The algorithm above works for Stokes Darcy, and Darcy Darcy interfaces: only u n is continuous. On Stokes Stokes interfaces u τ also has to be continuous. Augment S h with the tangential stress as primary variable: S h : ( ) n T n n T τ j where [ ] is the jump operator. ( [u n ] [u τ j ] S h is a Neumann-to-Dirichlet operator. Different number of primary variables on different interfaces. ) or, equivalently, S h : n T [u ], Department of Mathematics, University of Pittsburgh 3
31 Condition number for multiple subdomains Γ S : Stokes-Stokes interfaces Γ D : Stokes-Darcy and Darcy-Darcy interfaces C (h λ 2 Γ S + λ 2 Γ D ) s(λ, λ) C 2 ( λ 2 Γ S + h λ 2 Γ D ) C h λ 2 Γ s(λ, λ) C 2 h λ 2 Γ Theorem: cond(s h ) = O(h 2 ) Department of Mathematics, University of Pittsburgh 3
32 Condition number studies: two subdomains /h eigmin eigmax iter Department of Mathematics, University of Pittsburgh 32
33 Condition number studies: multiple subdomains 4 subdomains /h eigmin eigmax iter Varying number of subdomains /H eigmin eigmax iter Department of Mathematics, University of Pittsburgh 33
34 Coupling Stokes-Darcy flow with transport concentration of a chemical c(x, t) porosity of the medium φ(x) diffusion - dispersion tensor D t (x, t) (symmetric and positive definite) source Q φc t + (uc D t c) = φq, (x, t) Ω (, T ) IC: c(x, ) = c (x), x Ω BC: { (uc Dt c) n = uc I n on Γ in (D t c) n = on Γ out where Γ in := {x Ω : u n < } and Γ out := {x Ω : u n }. Department of Mathematics, University of Pittsburgh 34
35 Stability and accuracy of a discontinuous Galerkin method z = D t c; (c h, z h ) : approximation of (c, z); Norm for (c h, z h ): (c h, z h ) 2 := φ /2 c h (T ) 2,Ω + + T { u h n, (c h )2 Γ + l T where [c h ] is the jump across the interior boundary edge γ l. Theorem: (Vassilev and.) D /2 t z h 2,Ω dt u h n l, [c h ] 2 γl } dt (c h, z h ) { φ /2 c 2,Ω + T c I u n /2 2,Γ in dt} /2 + T φ /2 Q,Ω dt Theorem: (Vassilev and.) (c c h, z z h ) C(h k + u u h ) C(h k + h r + h r 2+ ) Department of Mathematics, University of Pittsburgh 35
36 Convergence test U V D = 3 I D = mesh c c h L (L 2 ) rate z z h L 2 (L 2 ) rate c c h L (L 2 ) rate 4x4.99e+ 8.95e-3 2.7e+ 8x8 3.27e e e x6 8.48e e e x e e e x64 5.6e e e-3.74 Department of Mathematics, University of Pittsburgh 36
37 Contaminant transport example Permeability perm V Computed velocity Department of Mathematics, University of Pittsburgh 37
38 Transport: initial plume.2.2 Time = C Time = Department of Mathematics, University of Pittsburgh 38
39 Transport: initial plume.2.2 Time = 3 C Time = 3 Department of Mathematics, University of Pittsburgh 39
40 Transport: initial plume.2.2 Time = 5 C Time = 5 Department of Mathematics, University of Pittsburgh 4
41 Transport: initial plume.2.2 Time = 9 C Time = 9 Department of Mathematics, University of Pittsburgh 4
42 Transport: initial plume.2.2 Time = 6 C Time = 6 Department of Mathematics, University of Pittsburgh 42
43 Transport: inflow.2.2 Time = 2 C Time = 2 Department of Mathematics, University of Pittsburgh 43
44 Transport: inflow.2.2 Time = C Time = Department of Mathematics, University of Pittsburgh 44
45 Transport: inflow.2.2 Time = 7 C Time = 7 Department of Mathematics, University of Pittsburgh 45
46 Summary Well-posed mathematical flow model based on continuity of flux and normal stress, and the Beavers-Joseph-Saffman condition Error analysis for mortar multiscale finite element approximations with very general mortar conditions Non-overlapping domain decomposition: s.p.d. interface problem; local Stokes and Darcy solves required at every CG iteration Condition number of interface problem: O(h ) for two subdomains; O(h 2 ) for multiple subdomains Conforming and DG Stokes elements and mixed finite elements for flow are coupled with LDG for transport Department of Mathematics, University of Pittsburgh 46
47 Current and future work Balancing preconditioner for the interface problem A posteriori error estimation and adaptive mesh refinement Extensions to Navier-Stokes and Brinkman Polygonal elements using MFD and DG P.2 K Permeability Pressure and velocity Department of Mathematics, University of Pittsburgh 47
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationCoupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements V.J. Ervin E.W. Jenkins S. Sun Department of Mathematical Sciences Clemson University, Clemson, SC, 29634-0975, USA. Abstract We study
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationInterior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids
Interior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids Gergina Pencheva a,, Ivan Yotov b a Center for Subsurface Modeling, Institute for Computational Engineering
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 40, No. 6, pp. 195 18 c 003 Society for Industrial and Applied Mathematics COUPLING FLUID FLOW WITH POROUS MEDIA FLOW WILLIAM J. LAYTON, FRIEDHELM SCHIEWECK, AND IVAN YOTOV Abstract.
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationMixed Multiscale Methods for Heterogeneous Elliptic Problems
Mixed Multiscale Methods for Heterogeneous Elliptic Problems Todd Arbogast Abstract We consider a second order elliptic problem written in mixed form, i.e., as a system of two first order equations. Such
More informationA Multiscale Mortar Multipoint Flux Mixed Finite Element Method
A Multiscale Mortar Multipoint Flux Mixed Finite Element Method Mary F. Wheeler Guangri Xue Ivan Yotov Abstract In this paper, we develop a multiscale mortar multipoint flux mixed finite element method
More informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationNumerical Simulations of Fluid Pressure in the Human Eye
Numerical Simulations of Fluid Pressure in the Human Eye T.R. Crowder V.J. Ervin Department of Mathematical Sciences Clemson University Clemson, SC, 2963-0975 USA. January 22, 2013 Abstract In this article
More informationBalancing domain decomposition for mortar mixed nite element methods
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2003; 10:159 180 (DOI: 10.1002/nla.316) Balancing domain decomposition for mortar mixed nite element methods Gergina Pencheva and
More informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationSpace-Time Domain Decomposition Methods for Transport Problems in Porous Media
1 / 49 Space-Time Domain Decomposition Methods for Transport Problems in Porous Media Thi-Thao-Phuong Hoang, Elyes Ahmed, Jérôme Jaffré, Caroline Japhet, Michel Kern, Jean Roberts INRIA Paris-Rocquencourt
More informationTexas at Austin 2 Nazarbayev University. January 15, 2019
Recovery of the Interface Velocity for the Incompressible Flow in Enhanced Velocity Mixed Finite Element Method arxiv:1901.04401v1 [math.na] 14 Jan 2019 Yerlan Amanbek 1,2, Gurpreet Singh 1, and Mary F.
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationApproximation of the Stokes-Darcy system by optimization
Approximation of the Stokes-Darcy system by optimization V.J. Ervin E.W. Jenkins Hyesuk Lee June 2, 202 Abstract A solution algorithm for the linear/nonlinear Stokes-Darcy coupled problem is proposed and
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationA Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 UNCORRECTED PROOF
1 A Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 Herbert Egger 1 and Christian Waluga 2 4 1 Center of Mathematics, Technische Universität München, Boltzmannstraße 3, 85748 5 Garching
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationOutline: 1 Motivation: Domain Decomposition Method 2 3 4
Multiscale Basis Functions for Iterative Domain Decomposition Procedures A. Francisco 1, V. Ginting 2, F. Pereira 3 and J. Rigelo 2 1 Department Mechanical Engineering Federal Fluminense University, Volta
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS
ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS VINCENT FONTAINE 1,2 AND ANIS YOUNES 2 1 Laboratoire de Physique des Bâtiments et des Systèmes, Université de La Réunion, 15 avenue
More informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
More informationOn some numerical convergence studies of mixed finite element methods for flow in porous media
On some numerical convergence studies of mixed finite element methods for flow in porous media Gergina Pencheva Abstract We consider an expanded mixed finite element method for solving second-order elliptic
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationA multiscale preconditioner for stochastic mortar mixed finite elements
A multiscale preconditioner for stochastic mortar mixed finite elements Mary F. Wheeler b, Tim Wildey b, and Ivan Yotov a a Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260 b The
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationIFE for Stokes interface problem
IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationITERATIVE SOLUTION OF A COMBINED MIXED AND STANDARD GALERKIN DISCRETIZATION METHOD FOR ELLIPTIC PROBLEMS
ITERATIVE SOLUTION OF A COMBINED MIXED AND STANDARD GALERKIN DISCRETIZATION METHOD FOR ELLIPTIC PROBLEMS RAYTCHO D. LAZAROV, JOSEPH E. PASCIAK, AND PANAYOT S. VASSILEVSKI Abstract. In this paper, we consider
More informationarxiv: v1 [math.na] 25 Aug 2016
A conservative local multiscale model reduction technique for Stoes flows in heterogeneous perforated domains Eric T. Chung Maria Vasilyeva Yating Wang arxiv:1608.07268v1 [math.na] 25 Aug 2016 October
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT DISCRETIZATION FOR STEADY STOKES FLOWS WITH THRESHOLD SLIP BOUNDARY CONDITION
DISCONTINUOUS GALRKIN FINIT LMNT DISCRTIZATION FOR STADY STOKS FLOWS WITH THRSHOLD SLIP BOUNDARY CONDITION J.K. Djoko Department of Mathematics and Applied Mathematics, University of Pretoria, Private
More informationBasics and some applications of the mortar element method
GAMM-Mitt. 28, No. 2, 97 123 (2005) Basics and some applications of the mortar element method Christine Bernardi 1, Yvon Maday 1, and Francesca Rapetti 2 1 Laboratoire Jacques-Louis Lions, C.N.R.S. & université
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More information