Space-Time Domain Decomposition Methods for Transport Problems in Porous Media
|
|
- Andrew Wade
- 5 years ago
- Views:
Transcription
1 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 Maison de la Simulation October 16, 2015 Enit Lamsin, Tunis
2 2 / 49 OUTLINE 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
3 3 / 49 Outline Introduction 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
4 4 / 49 Introduction Objective: to formulate numerical methods for flow and transport in heterogeneous porous media Examples of heterogeneous media:
5 4 / 49 Introduction Objective: to formulate numerical methods for flow and transport in heterogeneous porous media Examples of heterogeneous media: porous media around underground nuclear waste deposit sites
6 5 / 49 Introduction Heterogeneities mean difficulties for simulation Deep underground repository (High-level waste) A repository 2km 2km
7 5 / 49 Introduction Heterogeneities mean difficulties for simulation Deep underground repository (High-level waste) Different materials strong heterogeneity, different time scales. Large differences in spatial scales. Long-term computations. A repository 2km 2km
8 6 / 49 Introduction Objective: to formulate numerical methods for flow and transport in heterogeneous porous media Examples of heterogeneous media: porous media around underground nuclear waste deposit sites porous media with fractures
9 7 / 49 Introduction Difficulty for modeling flow in media with fractures A problem requiring multi-scale modelling
10 Introduction Difficulty for modeling flow in media with fractures A problem requiring multi-scale modelling Fractures represent heterogeneities in porous media Usually of much higher permeability than surrounding medium May be of much lower permeability so that they act as a barrier 7 / 49
11 Introduction Difficulty for modeling flow in media with fractures A problem requiring multi-scale modelling Fractures represent heterogeneities in porous media Usually of much higher permeability than surrounding medium May be of much lower permeability so that they act as a barrier Fracture width much smaller than any reasonable parameter of spatial discretization. 7 / 49
12 Introduction Difficulty for modeling flow in media with fractures A problem requiring multi-scale modelling Fractures represent heterogeneities in porous media Usually of much higher permeability than surrounding medium May be of much lower permeability so that they act as a barrier Fracture width much smaller than any reasonable parameter of spatial discretization. Different types of models for flow in fractures double continuum models. discrete fracture networks (DFN s) (no exchange with surrounding matrix rock) reduced fracture models (with exchange with matrix rock) 7 / 49
13 8 / 49 Introduction Objective here: to formulate methods for subdomain time-stepping More specifically: develop and compare two different space-time (global in time) domain decomposition methods for the linear transport problem in mixed formulation. extend these methods to the case of a domain with a discrete fracture extend these method to two phase flow models, with discontinuous capillary pressure (in progress)
14 9 / 49 Introduction Domain decomposition (DD) methods Domain decomposition in space t y x
15 9 / 49 Introduction Domain decomposition (DD) methods Domain decomposition in space t y x
16 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space t y x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain.
17 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space Space-time domain decomposition t t y y x x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain.
18 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space Space-time domain decomposition t t y y x x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain. Solve time-dependent problems in the subdomains Exchange information through the space-time interface Enable local discretizations both in space and in time local time stepping
19 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space t Space-time DD with Time windows t y y T1 x 0 x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain.
20 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space t Space-time DD with Time windows t y y T2 x T1 0 x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain.
21 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space t Space-time DD with Time windows t y T3 y T2 x T1 0 x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain.
22 Introduction 9 / 49 Domain decomposition (DD) methods Domain decomposition in space t Space-time DD with Time windows t y T y T2 T1 x 0 x Discretize in time and apply DD algorithm at each time step: Solve stationary problems in the subdomains Exchange information through the interface Use the same time step on the whole domain. Perform few iterations per window Use different space-time grids in each window Use the solution in the previous window to calculate a good initial guess on the interface.
23 10 / 49 Outline Pure diffusion problems 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
24 11 / 49 Model problem Pure diffusion problems Multi-domain mixed formulations Transport of a contaminant in a porous medium under the effect of diffusion, written in mixed form: L(c, r) := φ t c + div r = f in Ω (0, T ), M(c, r) := D 1 r + c = 0 in Ω (0, T ), c = 0 on Ω (0, T ), c(, 0) = c 0 in Ω, c concentration of a contaminant dissolved in a fluid, r diffusive flux. φ porosity; D symmetric, positive definite, time-independent diffusion tensor.
25 11 / 49 Model problem Pure diffusion problems Multi-domain mixed formulations Transport of a contaminant in a porous medium under the effect of diffusion, written in mixed form: L(c, r) := φ t c + div r = f in Ω (0, T ), M(c, r) := D 1 r + c = 0 in Ω (0, T ), c = 0 on Ω (0, T ), c(, 0) = c 0 in Ω, c concentration of a contaminant dissolved in a fluid, r diffusive flux. φ porosity; D symmetric, positive definite, time-independent diffusion tensor. Existence and uniqueness If D L (Ω), f L 2 (0, T ; L 2 (Ω)) and c 0 H0 1 (Ω) then problem above has a unique weak solution ( ) (c, r) H 1 (0, T ; L 2 (Ω)) L 2 (0, T ; H(div, Ω)) L (0, T ; L 2 (Ω)). Moreover, if D W 1, (Ω), f H 1 (0, T ; L 2 (Ω)) and c 0 H 2 (Ω) H0 1 (Ω), then ( ) (c, r) W 1, (0, T ; L 2 (Ω)) L (0, T ; H(div, Ω)) H 1 (0, T ; L 2 (Ω)). Proof. Galerkin s method and a priori estimates.
26 Pure diffusion problems Multi-domain mixed formulations Multi-domain problem t T y Γ Ω 1 Ω 2 x 0 Equivalent multi-domain problem: L(c 1, r 1 ) = f, on Ω 1 (0, T ), M(c 1, r 1 ) = 0, on Ω 1 (0, T ), c 1 (, 0) = c 0, in Ω 1, L(c 2, r 2 ) = f, on Ω 2 (0, T ), M(c 2, r 2 ) = 0, on Ω 2 (0, T ), c 2 (, 0) = c 0, in Ω 2, 12 / 49
27 Pure diffusion problems Multi-domain mixed formulations Multi-domain problem t T y Γ Ω 1 Ω 2 x 0 Equivalent multi-domain problem: L(c 1, r 1 ) = f, on Ω 1 (0, T ), M(c 1, r 1 ) = 0, on Ω 1 (0, T ), c 1 (, 0) = c 0, in Ω 1, L(c 2, r 2 ) = f, on Ω 2 (0, T ), M(c 2, r 2 ) = 0, on Ω 2 (0, T ), c 2 (, 0) = c 0, in Ω 2, together with the transmission conditions on the space-time interface c 1 = c 2 r 1 n 1 + r 2 n 2 = 0 on Γ (0, T ). 12 / 49
28 13 / 49 An overview Pure diffusion problems Multi-domain mixed formulations Different (equivalent) transmission conditions (TCs)
29 13 / 49 An overview Pure diffusion problems Multi-domain mixed formulations Different (equivalent) transmission conditions (TCs) GTP Schur Physical TCs + N-N preconditioner GTO Schwarz More general TCs with optimized parameters accelerate the convergence rate. Robin TCs Ventcell TCs
30 13 / 49 An overview Pure diffusion problems Multi-domain mixed formulations Different (equivalent) transmission conditions (TCs) GTP Schur Physical TCs + N-N preconditioner GTO Schwarz More general TCs with optimized parameters accelerate the convergence rate. Robin TCs Ventcell TCs Substructuring technique: Space-time interface problem
31 13 / 49 An overview Pure diffusion problems Multi-domain mixed formulations Different (equivalent) transmission conditions (TCs) GTP Schur Physical TCs + N-N preconditioner GTO Schwarz More general TCs with optimized parameters accelerate the convergence rate. Robin TCs Ventcell TCs Substructuring technique: Space-time interface problem Iterative solvers (GMRES, Richardson iteration)
32 T.T.P. Hoang, J. Jaffré, C. Japhet, M.K, and J. E. Roberts. Space-time domain decomposition methods for diffusion problems in mixed formulations. SIAM J. Numer. Anal., 51(6): , / 49 An overview Pure diffusion problems Multi-domain mixed formulations Different (equivalent) transmission conditions (TCs) GTP Schur Physical TCs + N-N preconditioner GTO Schwarz More general TCs with optimized parameters accelerate the convergence rate. Robin TCs Ventcell TCs Substructuring technique: Space-time interface problem Iterative solvers (GMRES, Richardson iteration)
33 14 / 49 Pure diffusion problems Multi-domain mixed formulations Time-dependent Steklov-Poincaré operator Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, on Ω i (0, T ), M(c i, r i ) = 0, on Ω i (0, T ), c i = λ, on Γ (0, T ), c i (, 0) = c 0, in Ω i.
34 14 / 49 Pure diffusion problems Multi-domain mixed formulations Time-dependent Steklov-Poincaré operator Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, on Ω i (0, T ), M(c i, r i ) = 0, on Ω i (0, T ), c i = λ, on Γ (0, T ), c i (, 0) = c 0, in Ω i. Space-time interface problem: S1 DtN (λ, f, c 0 ) + S2 DtN (λ, f, c 0 ) = 0, 2 2 Si DtN (λ, 0, 0) = Si DtN (0, f, c 0 ), i=1 i=1 Sλ = χ, on Γ (0, T ).
35 14 / 49 Pure diffusion problems Multi-domain mixed formulations Time-dependent Steklov-Poincaré operator Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, on Ω i (0, T ), M(c i, r i ) = 0, on Ω i (0, T ), c i = λ, on Γ (0, T ), c i (, 0) = c 0, in Ω i. Space-time interface problem: S1 DtN (λ, f, c 0 ) + S2 DtN (λ, f, c 0 ) = 0, 2 2 Si DtN (λ, 0, 0) = Si DtN (0, f, c 0 ), i=1 i=1 Sλ = χ, on Γ (0, T ). Neumann-Neumann preconditioner with weights: ( ) σ 1 S1 NtD + σ 2 S2 NtD Sλ = χ, on Γ (0, T ), where σ i : Γ (0, T ) [0, 1] such that σ 1 + σ 2 = 1.
36 15 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Robin transmission conditions Equivalent Robin TCs on Γ (0, T ): for α 1, α 2 > 0 r 1 n 1 + α 1 c 1 = r 2 n 1 + α 1 c 2, r 2 n 2 + α 2 c 2 = r 1 n 2 + α 2 c 1,
37 15 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Robin transmission conditions Equivalent Robin TCs on Γ (0, T ): for α 1, α 2 > 0 r 1 n 1 + α 1 c 1 = r 2 n 1 + α 1 c 2, r 2 n 2 + α 2 c 2 = r 1 n 2 + α 2 c 1, Robin-to-Robin operators, for i = 1, 2 and j = 3 i: where (c i, r i ), i = 1, 2, is the solution of S RtR i : (ξ i, f, c 0 ) ( r i n j + α j c i ) Γ, L(c i, r i ) = f, on Ω i (0, T ), M(c i, r i ) = 0, on Ω i (0, T ), r i n i + α i c i = ξ i, on Γ (0, T ), c i (, 0) = c 0, in Ω i.
38 15 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Robin transmission conditions Equivalent Robin TCs on Γ (0, T ): for α 1, α 2 > 0 r 1 n 1 + α 1 c 1 = r 2 n 1 + α 1 c 2, r 2 n 2 + α 2 c 2 = r 1 n 2 + α 2 c 1, Robin-to-Robin operators, for i = 1, 2 and j = 3 i: where (c i, r i ), i = 1, 2, is the solution of S RtR i : (ξ i, f, c 0 ) ( r i n j + α j c i ) Γ, L(c i, r i ) = f, on Ω i (0, T ), M(c i, r i ) = 0, on Ω i (0, T ), r i n i + α i c i = ξ i, on Γ (0, T ), c i (, 0) = c 0, in Ω i. Space-time interface problem with two Lagrange multipliers: or equivalently, ξ 1 ξ 2 = S2 RtR (ξ 2, f, c 0 ), = S1 RtR (ξ 1, f, c 0 ), on Γ (0, T ), S R ( ξ1 ξ 2 ) = χ R, on Γ (0, T ).
39 16 / 49 Pure diffusion problems OSWR algorithm with Robin TCs Multi-domain mixed formulations OSWR iterative algorithm: at the k th iteration, for i = 1, 2 and j = (3 i) φ i t ci k + div r k i = f in Ω i (0, T ) D 1 i r k i + ci k = 0 in Ω i (0, T ) r k i n i + α i ci k ( = r k 1 j n i + α i c k 1 j on Γ (0, T ), ) for given initial guess g i = r 0 i n i + α i ci 0, i = 1, 2.
40 16 / 49 Pure diffusion problems OSWR algorithm with Robin TCs Multi-domain mixed formulations OSWR iterative algorithm: at the k th iteration, for i = 1, 2 and j = (3 i) φ i t ci k + div r k i = f in Ω i (0, T ) D 1 i r k i + ci k = 0 in Ω i (0, T ) r k i n i + α i ci k ( = r k 1 j n i + α i c k 1 j on Γ (0, T ), ) for given initial guess g i = r 0 i n i + α i ci 0, i = 1, 2. Theorem (Convergence of OSWR algorithm in mixed formulation) ( ) If the algorithm above is initialized by (g i ) H 1 0, T ; L 2 (Γ) for i = 1, 2, then ( ) ( ) a sequence of iterates ci k, r k i H 1 0, T ; L 2 (Ω i ) L 2 (0, T ; H (div, Ω i )) is well-defined 2 i=1 ( ) ci k c Ωi H1 (0,T ;L 2 (Ω i )) + rk i r Ωi 2 k L 2 0. (0,T ;H(div,Ω i )) where H (div, Ω i ) := { } v H (div, Ω i ) : (v n i ) Γ L 2 (Γ). Remark. The proof is carried out for the multiple subdomain case.
41 17 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Ventcell transmission conditions With sufficient regularity equivalent Ventcell transmission conditions
42 17 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Ventcell transmission conditions With sufficient regularity equivalent Ventcell transmission conditions In primal form: on Γ (0, T ): r 1 n 1 + α 1 c 1 + = r 2 n 1 + α 1 c 2 + r 2 n 2 + α 2 c 2 + = r 1 n 2 + α 2 c 1 +
43 Pure diffusion problems Multi-domain mixed formulations 17 / 49 GTO Schwarz: Ventcell transmission conditions With sufficient regularity equivalent Ventcell transmission conditions In primal form: on Γ (0, T ): ( r 1 n 1 + α 1 c 1 + β 1 φ2 t c 1 + div τ ( D 2,Γ τ c 1 ) ) = r 2 n 1 + α 1 c 2 + ( β 1 φ2 t c 2 + div τ ( D 2,Γ τ c 2 ) ), ( r 2 n 2 + α 2 c 2 + β 2 φ1 t c 2 + div τ ( D 1,Γ τ c 2 ) ) = r 1 n 2 + α 2 c 1 + ( β 2 φ1 t c 1 + div τ ( D 1,Γ τ c 1 ) ). α i, β i : positive constants to be optimized to accelerate convergence rate.
44 17 / 49 Pure diffusion problems Multi-domain mixed formulations GTO Schwarz: Ventcell transmission conditions With sufficient regularity equivalent Ventcell transmission conditions In primal form: on Γ (0, T ): r 1 n 1 + α 1 c 1 + β 1 ( φ2 t c 1 + div τ ( D 2,Γ τ c 1 ) ) = r 2 n 1 + α 1 c 2 + r 2 n 2 + α 2 c 2 + β 2 ( φ1 t c 2 + div τ ( D 1,Γ τ c 2 ) ) = r 1 n 2 + α 2 c 1 + β 1 ( φ2 t c 2 + div τ ( D 2,Γ τ c 2 ) ), β 2 ( φ1 t c 1 + div τ ( D 1,Γ τ c 1 ) ). α i, β i : positive constants to be optimized to accelerate convergence rate. In mixed form: introduce Lagrange multipliers on the interface, c i,γ and r Γ,i, for i = 1, 2, ( ) r i n i + α i c i,γ + β i φj t c i,γ + div τ r Γ,i = r j n i + α i c j,γ + ( ( )) β i φ j t c j,γ + div τ D j,γ D 1 i,γ r Γ,j, D 1 j,γ r Γ,i + τ c i,γ = 0. c i,γ : pressure trace on Γ. r Γ,i := D j,γ τ c i,γ : NOT the tangential trace of r i on Γ (0, T ).
45 18 / 49 Pure diffusion problems Nonconforming discretizations in time Nonconforming discretizations in time t T t 1 n t 2 m 0 Ω 1 Ω 2 x
46 18 / 49 Pure diffusion problems Nonconforming discretizations in time Nonconforming discretizations in time t T t 1 n t 2 m 0 Ω 1 Ω 2 x Time discretization: non-conforming time grids T 1, T 2 ; discontinuous Galerkin with piecewise polynomials of degree 0. Projections: Π ji is an L 2 projection from piecewise constant functions on T i onto piecewise constant functions on T j. Ex: (Π 21 (λ 1 )) m = 1 M 1 Jm 2 λ 1, for m = 1,, M 2. n=1 J 1 n J 2 m
47 19 / 49 Pure diffusion problems Nonconforming discretizations in time Semi-discrete transmission conditions with nonconforming time grids ( ) For GTP Schur method: take λ = λ 1,, λ M 1 piecewise constant on ( ) Jn 1 = t1 n, t n+1 1, for n = 0,, M 1 1. Continuity of concentration: c 1 = Π 11 (λ) and c 2 = Π 21 (λ). Conservation of the flux over each time subinterval (Π 11 (r 1 n 1 ) + Π 12 (r 2 n 2 )) dt = 0, for n = 0,, M 1 1. J 1 n Γ
48 19 / 49 Pure diffusion problems Nonconforming discretizations in time Semi-discrete transmission conditions with nonconforming time grids ( ) For GTP Schur method: take λ = λ 1,, λ M 1 piecewise constant on ( ) Jn 1 = t1 n, t n+1 1, for n = 0,, M 1 1. Continuity of concentration: c 1 = Π 11 (λ) and c 2 = Π 21 (λ). Conservation of the flux over each time subinterval (Π 11 (r 1 n 1 ) + Π 12 (r 2 n 2 )) dt = 0, for n = 0,, M 1 1. J 1 n Γ For GTO Schwarz method: conservation of the two Robin (Ventcell) conditions across the interface over each time subinterval [( r 1 n 1 + α 1 c 1 ) Π 12 ( r 2 n 1 + α 1 c 2 )] dt = 0, n = 0,, M 1 1, J 1 n J 2 m Γ Γ [( r 2 n 2 + α 2 c 2 ) Π 21 ( r 1 n 2 + α 2 c 1 )] dt = 0, m = 0,, M 2 1. Convergence of semi-discrete, nonconforming in time, OSWR algorithm
49 20 / 49 Outline Advection-diffusion problems 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
50 21 / 49 Advection-diffusion problems Operator splitting Extension to advection-diffusion problems Linear advection-diffusion equation: Operator splitting φ tc + div (uc) + div r = f in Ω (0, T ), c + D 1 r = 0 in Ω (0, T ), c = 0 on Ω (0, T ), c(, 0) = c 0 in Ω. Advection eq.: explicit Euler + upwind, cell-centered finite volumes. Diffusion eq.: implicit Euler + mixed finite elements. CFL condition: sub-time steps for the advection. T = N 1 t 1 = N 2 t 2 t T L 1 t 1,a = t 1 t 2 = L 2 t 2,a 0 Ω 1 Ω 2 x
51 Advection-diffusion problems Discrete interface problems Operator splitting GTP Schur method: S h ( λa λ ) = χ h, on Γ (0, T ). = Generalized Neumann-Neumann preconditioner GTO Schwarz method with Robin TCs: λ a S R,h ξ 1 = χ R,h, on Γ (0, T ). ξ 2 = Optimized Robin parameters for the diffusion eq. only fully implicit scheme. Remark. λ a Λ N L h while λ, ξ 1, ξ 2 Λ N h. T.T.P. Hoang, J. Jaffré, C. Japhet, M.K., and J. E. Roberts. Proc. Mamern 2015, in preparation. 22 / 49
52 Advection-diffusion problems Operator splitting Test case 2: A near-field simulation (project PAMINA ) Performance Assessment Methodologies IN Application to Guide the Development of the Safety Case Parameters of the simulation Material Permeability (m.s 1 ) Porosity Diffusion (m 2. s 1 ) Host rock EDZ Vitrified waste m 0.6m 100m m / 49
53 Advection-diffusion problems Operator splitting Advection field: Darcy flow div u = 0 in Ω, u = K p in Ω. BCs: Homogeneous Neumann at x = 0 and x = 10, Dirichlet conditions with p = 100 Pa at y = 0 and p = 0 at y = / 49
54 25 / 49 Advection-diffusion problems Operator splitting Transport problem: time windows and decomposition Ω 7 Ω 8 Ω 9 Ω 4 Ω 5 Ω 6 Ω 1 Ω 2 Ω 3 Final time: T f = s ( years) 200 time windows with size T = 10 9 s. Decomposition into 9 subdomains. Nonconforming time grids: Diffusion step: t i = T /500, i = 5, t i = T /100, i 5. calculation and its decomposition. Diffusion-dominated: Pe L t a,i = t i. Non-uniform mesh in space: uniform mesh in the repository (10 by 10), then progressively coarser with a factor of 1.05.
55 26 / 49 Advection-diffusion problems Evolution of concentration field Operator splitting 100 years 5000 years years years In the repository In the host rock
56 27 / 49 Advection-diffusion problems Operator splitting Performance of one time window Log of error in c No Precond. Schur Precond. Schur Opt. Schwarz Error in c Error in r Number of subdomain solves Convergence with GMRES t Error in c with nonconforming time grids t Error in r with nonconforming time grids. Time grid 1: t i = T /500, i Time grid 2: t 5 = T /500, t i = T /100, i 5 Time grid 3: t 5 = T /100, t i = T /500, i 5 Time grid 4: t i = T /100, i
57 A subsurface waste storage simulation Zone Hydraulic conductivity Porosity Molecular diffusion K (m/year) φ d m (m 2 /year) terrain dalleprotec/dalleobtur voile remplissage conteneur1/conteneur dechet1/dechet radier drainant dalleobtur dalleprotec terrain voile remplissage conteneur2 conteneur1 dechet1 9.6 m dechet2 drainant forme radier drain 12 m
58 Advection-diffusion problems Operator splitting 29 / 49 Darcy flow div u = 0 in Ω, u = K h in Ω. BCs: Homogeneous Neumann at x = 0 and x = 12m, Dirichlet conditions with h = 9.998m at y = 0 and h = 10m at y = 9.6m Hydraulic head field
59 30 / 49 Advection-diffusion problems Operator splitting Concentration field after 500 years
60 31 / 49 Outline Extension to two-phase flow 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
61 Extension to two-phase flow 32 / 49 Model problem: Two phase immiscible flow Mathematical model t (ωρ α S α ) + div (ρ α u α ) = q α mass conservation u α = k rα µ α K ( p α ρ α g) Darcy s law S n + S w = 1 p n p w = π(s w ) capillary pressure Phase α = w water, n gas or oil. π(s w ) increasing function on [0, 1] (extend coninuously to R). ω porosity S α phase saturation u α phase velocity k rα relative permeability K permeability p α : phase pressure ρ α phase density µ α viscosity
62 33 / 49 Simplified model Extension to two-phase flow Enchery et al. (06), Cances (08), Brenner et al. (13), no gravity 1 Global pressure (Chavent) P g (S) = p w + 2 Kirchhoff transformation : φ(s) = S 0 K S 0 k rn (u)/µ n π (u) du, k rn(u) k rw (u) µ n + µ w k rn (u)k rw (u) µ n k rw (u) + µ w k rn (u) π (u) du.
63 33 / 49 Simplified model Extension to two-phase flow Enchery et al. (06), Cances (08), Brenner et al. (13), no gravity 1 Global pressure (Chavent) P g (S) = p w + 2 Kirchhoff transformation : φ(s) = S 0 K S 0 k rn (u)/µ n π (u) du, k rn(u) k rw (u) µ n + µ w k rn (u)k rw (u) µ n k rw (u) + µ w k rn (u) π (u) du. µ w k rn (S) Transformed system : f (S) = µ w k rn (s) + µ n k rw (S), λ(s) = k rn(s) + k rw(s). µ n µ w { ω t S + div (f (S)q T ) φ(s) = 0 in Ω [0, T ] div q T = 0, q T = K λ(s) grad P g
64 33 / 49 Simplified model Extension to two-phase flow Enchery et al. (06), Cances (08), Brenner et al. (13), no gravity 1 Global pressure (Chavent) P g (S) = p w + 2 Kirchhoff transformation : φ(s) = S 0 K S 0 k rn (u)/µ n π (u) du, k rn(u) k rw (u) µ n + µ w k rn (u)k rw (u) µ n k rw (u) + µ w k rn (u) π (u) du. µ w k rn (S) Transformed system : f (S) = µ w k rn (s) + µ n k rw (S), λ(s) = k rn(s) + k rw(s). µ n µ w { ω t S + div (f (S)q T ) φ(s) = 0 in Ω [0, T ] div q T = 0, q T = K λ(s) grad P g Simplified system: neglect advection ω t S φ(s) = 0 in Ω [0, T ] Nonlinear (degenerate) diffusion equation
65 Extension to two-phase flow Discontinuous capillary pressure: transmission conditions Two subdomains Ω = Ω 1 Ω 2, Ω 1 Ω 2 =. Γ = Ω 1 Ω 2 P_c2(1) P_c1(1) Capillary pressure P_c P_c2(0) P_{c1}(0) 0 s_2 s_1 1 Saturation Transmission conditions on the interface Continuity of capillary pressure π 1 (S 1 ) = π 2 (S 2 ) on Γ Continuity of the flux φ 1 (S 1 ).n 1 = φ 2 (S 2 ).n 2 on Γ Chavent Jaffré (86), Enchéry et al. (06), Cances (08), Ern et al (10), Brenner et al. (13). 34 / 49
66 35 / 49 Extension to two-phase flow Non-linear Schwarz algorithm Robin transmission conditions φ 1 (S 1 ).n 1 + β 1 π 1 (S 1 ) = φ 2 (S 2 ).n 2 + β 1 π 2 (S 2 ) φ 2 (S 2 ).n 2 + β 2 π 2 (S 2 ) = φ 1 (S 1 ).n 1 + β 2 π 1 (S 1 ) Schwarz algorithm Given S 0 i, iterate for k = 0,... Solve for S k+1 i, i = 1, 2, j = 3 i ω t S k+1 i φ i (S k+1 i ) = 0 in Ω i [0, T ] φ i (S k+1 i ).n i + β i π i (S k+1 i ) = φ j (S k j ).n j + β i π j (S k j ) on Γ [0, T ], (β 1, β 2 ) are free parameters chosen to accelerate convergence Basic ingredient: subdomain solver with Robin bc.
67 36 / 49 Extension to two-phase flow Finite volume scheme (1) Extension to Robin bc of cell centered FV scheme by Enchéry et al. (06). Triangulation T, cells K T, boundary faces σ Γ. Unknowns : cell values (S K ) K T, boundary face values (S σ ) σ EΓ Notations: K L = edge between K and L, τ K L = average). m(k L) K K L (eg harmonic
68 37 / 49 Extension to two-phase flow Finite volume scheme (2) Interior equation m(k ) Sn+1 K S n K δt + L N (K ) + σ E Γ E K τ K,σ ( ) τ K L φ(s n+1 K ) φ(s n+1 L ) ( ) φ(s n+1 K ) φ(s n+1 σ ) = 0, K T. Robin BC for boundary faces ( ) τ K,σ φ(s n+1 K ) φ(s n+1 σ ) + βm(σ)π(s n+1 σ ) = g σ, σ E Γ Implemented with Matlab Reservoir Simulation Toolbox (K. A. Lie et al. (14)) Solver with automatics differentiation : no explicit computation of Jacobian
69 Extension to two-phase flow 10 2 alpha=0.009,beta=0.008 Numerical example Homogeneous medium, Ω 1 = (0, 100) 3, Ω 2 = (100, 200) (0, 100) 2. Mobilities λ 0 (S) = S, S [0, 1], Capillary pressure π(s) = 5S 2, S [0, 1] 10 0 alpha=8, beta=0.008 alpha=4,beta=4 alpha=0.008, beta= Convergence history for various parameters 38 / 49
70 39 / 49 Extension to two-phase flow Evolution of the concentration
71 40 / 49 Outline Extension to reduced fracture models 1 Introduction 2 Pure diffusion problems Multi-domain mixed formulations Nonconforming discretizations in time 3 Advection-diffusion problems Operator splitting 4 Extension to two-phase flow 5 Extension to reduced fracture models
72 t p i + div u i = q i in Ω i (0, T ), 41 / 49 and we denote by n i the outward unit normal vector field on Ω i. We use the convention Extension that to reduced for anyfracture scalar, models vector or tensor valued function φ defined on Ω, φ i denotes the A reduced model: δ interface-fracture Alboin-Jaffré-Roberts-Serres (2002) n 1 Martin-Jaffré-Roberts (2005) Knabner-Roberts (2014) γ 1 (Forchheimer γ 2 flow) Ω 1 Ω f Ω 2 In this work: assumenthat 2 D/δ large concentration continuity across the fracture In the subdomains n 1 Ω 1 Ω 2 Figure φ 1. Left: The domain Ω with the fracture Ω f. Right: The domain Ω with the i t c i + div r i = f i in Ω interface-fracture i (0, T ), γ. r i = D i c i in Ω i (0, T ), restriction of φ to Ω i, i = 1, c2, i f. We = 0rewrite problem on Ω (1) i as Ω the following (0, T ), transmission for i = problem: 1, 2, c i = c γ on γ (0, T ), c t p i + div u i = q i in Ω i (0, T ), i = 1, 2, f, i (, 0) = c 0,i in Ω i, u i = K i p i in Ω i (0, T ), i = 1, 2, f, p i = 0 on Ω i Ω (0, T ), i = 1, 2, f, p i = p f on γ i (0, T ), i = 1, 2, u i n i = u f n i on γ i (0, T ), i = 1, 2, p i (, 0) = p 0,i in Ω i, i = 1, 2, f. 2. Reduced fracture model In this section, the model where the fracture is treated as an interface γ between subdomains Ω 1 and Ω 2 [1,4] is considered. We use the notation τ and div τ for the tangential gradient and divergence, respectively. A reduced model consists of the equations in the subdomains n 2 γ (2)
73 t p i + div u i = q i in Ω i (0, T ), 41 / 49 and we denote by n i the outward unit normal vector field on Ω i. We use the convention Extension that to reduced for anyfracture scalar, models vector or tensor valued function φ defined on Ω, φ i denotes the A reduced model: δ interface-fracture Alboin-Jaffré-Roberts-Serres (2002) n 1 Martin-Jaffré-Roberts (2005) Knabner-Roberts (2014) γ 1 (Forchheimer γ 2 flow) Ω 1 Ω f Ω 2 In this work: assumenthat 2 D/δ large concentration continuity across the fracture In the subdomains n 1 Ω 1 Ω 2 Figure φ 1. Left: The domain Ω with the fracture Ω f. Right: The domain Ω with the i t c i + div r i = f i in Ω interface-fracture i (0, T ), γ. r i = D i c i in Ω i (0, T ), restriction of φ to Ω i, i = 1, c2, i f. We = 0rewrite problem on Ω (1) i as Ω the following (0, T ), transmission for i = problem: 1, 2, c i = c γ on γ (0, T ), and in the fracture c t p i + div u i = q in Ω (0, T ), i = 1, 2, f, i (, = c 0) u i 0,i i in = K i p i Ω i, i in Ω i (0, T ), i = 1, 2, f, p i = 0 on Ω i Ω (0, T ), i = 1, 2, f, p i = p f on γ i (0, T ), i = 1, 2, (2) u i n i i = 1, 2, φ γ t c γ + div p i (, τ 0) r γ ini γ = 1, (0, 2, f. T ), = u f n i on γ i (0, T ), = f γ p 0,i + ( ) r 1 nin 1 γ Ω i, + r 2 n 2 γ r γ = D γδ τ c γ in γ (0, T ), c γ 2. = Reduced 0 fracture modelon γ (0, T ), In this section, the model c γ(, where 0) the= fracture c 0,γ is treated as an interface γ between in γ. subdomains Ω 1 and Ω 2 [1,4] is considered. We use the notation τ and div τ for the tangential gradient and divergence, respectively. A reduced Communication model consists between of the equations the fracture in the and subdomains the rock matrix. n 2 γ
74 42 / 49 Extension to reduced fracture models Formulation as an interface problem (GTP Schur) The same (as with simple DD) Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), c i = λ, on γ (0, T ), c i (, 0) = c 0, in Ω i.
75 42 / 49 Extension to reduced fracture models Formulation as an interface problem (GTP Schur) The same (as with simple DD) Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of Different space-time interface problem: 2 i=1 S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), c i = λ, on γ (0, T ), c i (, 0) = c 0, in Ω i. instead of 2 Si DtN (λ, 0, 0) = Si DtN (0, f, c 0 ), i=1 Sλ = χ, on γ (0, T ).
76 42 / 49 Extension to reduced fracture models Formulation as an interface problem (GTP Schur) The same (as with simple DD) Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of Different space-time interface problem: S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), c i = λ, on γ (0, T ), c i (, 0) = c 0, in Ω i. L γ(λ, r γ) + Sλ = χ + f γ, in γ (0, T ), M γ(λ, r γ) = 0 in γ (0, T ), λ(, 0) = c 0,γ, in γ.
77 42 / 49 Extension to reduced fracture models Formulation as an interface problem (GTP Schur) The same (as with simple DD) Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of Different space-time interface problem: Two possible preconditionners: S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), c i = λ, on γ (0, T ), c i (, 0) = c 0, in Ω i. L γ(λ, r γ) + Sλ = χ + f γ, in γ (0, T ), M γ(λ, r γ) = 0 in γ (0, T ), λ(, 0) = c 0,γ, in γ. a Neumann-Neumann preconditionner with weights
78 Extension to reduced fracture models Formulation as an interface problem (GTP Schur) The same (as with simple DD) Dirichlet-to-Neumann operators, for i = 1, 2: where (c i, r i ), i = 1, 2, is the solution of Different space-time interface problem: Two possible preconditionners: S DtN i : (λ, f, c 0 ) (r i n i ) Γ, L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), c i = λ, on γ (0, T ), c i (, 0) = c 0, in Ω i. L γ(λ, r γ) + Sλ = χ + f γ, in γ (0, T ), M γ(λ, r γ) = 0 in γ (0, T ), λ(, 0) = c 0,γ, in γ. a Neumann-Neumann preconditionner with weights a local preconditioner (coming from the observation that the interface problem is dominated by the 2nd order operator, Amir, MK, Martin, Robert, Arima 06) 42 / 49
79 43 / 49 Extension to reduced fracture models Transmission conditions for a GTO Schwarz method Taking a linear combination of the transmission conditions for the GTP Schur method we obtain: r 1 n 1 + α 1 c 1,γ + φ γ t c i,γ + div τ r γ,1 = r 2 n 1 + α 1 c 2,γ + f γ r γ,1 = D γδ τ c 1,γ r 2 n 2 + α 2 c 2,γ + φ γ t c 2,γ + div τ r γ,2 = r 1 n 2 + α 2 c 1,γ + f γ r γ,2 = D γδ τ c 2,γ
80 Extension to reduced fracture models 44 / 49 Formulation as an interface problem (GTO Schwarz) We use Ventcell to Robin operators, for i = 1, 2: S VtR i : (θ i, f, c 0, f γ, c 0,γ ) ( r i n j + αc i ) Γ, where ( c i, r i, c i,γ, r γ,i ), i = 1, 2, is the solution of L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), r i n i + αc i,γ + φ γ t c i,γ + div τ r γ,i = θ i on γ (0, T ), r γ,i + D γδ τ c i,γ = 0, on γ (0, T ), c i (, 0) = c 0, in Ω i c i,γ (, 0) = c 0,γ, in γ.
81 44 / 49 Extension to reduced fracture models Formulation as an interface problem (GTO Schwarz) We use Ventcell to Robin operators, for i = 1, 2: S VtR i : (θ i, f, c 0, f γ, c 0,γ ) ( r i n j + αc i ) Γ, where ( c i, r i, c i,γ, r γ,i ), i = 1, 2, is the solution of L(c i, r i ) = f, in Ω i (0, T ), M(c i, r i ) = 0, in Ω i (0, T ), r i n i + αc i,γ + φ γ t c i,γ + div τ r γ,i = θ i on γ (0, T ), r γ,i + D γδ τ c i,γ = 0, on γ (0, T ), c i (, 0) = c 0, in Ω i c i,γ (, 0) = c 0,γ, in γ. Space-time interface problem: θ 1 θ 2 = S2 VtR (θ 2, f, c 0, f γ, c 0,γ ) + f γ, on γ (0, T ), = S1 VtR (θ 1, f, c 0, f γ, c 0,γ ) + f γ, on γ (0, T ).
82 45 / 49 Numerical results Extension to reduced fracture models 1 c=0 c=0 0 c= c=1 Geometry and boundary conditions. Isotropic coefficients: D i = 1, i = 1, 2, and D γ = 1/δ = 1000.
83 45 / 49 Numerical results Extension to reduced fracture models 1 c=0 c=0 0 c= c=1 Geometry and boundary conditions. Isotropic coefficients: D i = 1, i = 1, 2, and D γ = 1/δ = Zero source terms and initial condition. Spatial discretization: uniform rectangular mesh h = 1/100 mixed FE with the lowest-order Raviart-Thomas spaces. Time discretization (case 1): conforming grids t m = t γ = T /300 with T = 0.5.
84 46 / 49 Extension to reduced fracture models Snapshots of solution - concentration field c t= t t=t/4 t=t/2 t=t
85 47 / 49 Extension to reduced fracture models Snapshots of solution - diffusive flux r 1 1 t= t t=t/ t=t/ t=t
86 Extension to reduced fracture models Convergence - GMRES Number of Iterations Number of Iterations Figure 4. Convergence curves for the compressible flow: errors in p (on the lef TITLE TITLE WILL WILL BE SET BE BY SETTHE BY PUBLISHER THE PUBLISHER (on the right) - Method 1 with 13no preconditioner 13 (blue), Method 1 with local pr (green), Method 1 with Neumann-Neumann preconditioner (cyan) and Metho Error in p Error in p Error in u Error in u Error in u Number Number of Iterations of Iterations Number Number of Iterations of Iterations Figure 5. The errors (in logarithmic scale) after 10 Jacobi iterations for variou Figure Figure 4. Convergence 4. curves curves for the forcompressible the flow: flow: errors errors in p in (onp (onthe left) the Robin left) andparameter. in andu in uthe red star shows the optimized parameters computed by (on the (onright) the right) - Method - Method 1 with 1 with no preconditioner no (blue), (blue), Method Method 1 with 1 with local local preconditioner minimizing the continuous convergence factor. (green), (green), Method Method 1 with 1 with Neumann-Neumann preconditioner (cyan) (cyan) and Method and Method 2 (red). 2 (red). T.T.P. Hoang, J. Jaffré, 8 8 C. Japhet, M.K., and J. E. Roberts. Space-time Domain Error in u Error in u Anal., to appear, α L 2 concentration errors (c) L 2 flux errors (r) L 2 error versus α GT Schur with no preconditioner GTP Schur with local preconditioner GTP Schur with NN preconditioner GTO Schwarz method the number of iterations for different schemes are shown: errors in p (on the left) and in see that Method 2 is superior to the other schemes. For Method 1, the Neumann-Neumann a better acceleration than the local one, and both are better than without a precondition For Method 2, we vary Robin parameters and plot the logarithmic scale of the error iterations. In Figure 5, we see that the optimized Robin parameter (the red star), w numerically minimizing the convergence factor, are located close to those giving the sm same number of iterations. Decomposition and Mixed Formulation for reduced fracture models. SIAM J. Numer / 49
87 49 / 49 Extension to reduced fracture models Conclusions perspectives Space time DD method with Robi TC for diffusion and advection diffusion Extension to fractured media Convergence for GTP Schur (Gander et al. for homogeneous media) Convergence for fractured media Influence of Robin parameter β, find optimal parameter Study interface problem for non-linear case, Jacobi (SWR) vs Newton Extension to full two-phase model Convergence of Schwarz alg. for nonlinear case Large scale parallel solver (MdS)
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure Elyes Ahmed, Caroline Japhet, Michel Kern INRIA Paris Maison de la Simulation ENPC University Paris
More informationA posteriori error estimates for space-time domain decomposition method for two-phase flow problem
A posteriori error estimates for space-time domain decomposition method for two-phase flow problem Sarah Ali Hassan, Elyes Ahmed, Caroline Japhet, Michel Kern, Martin Vohralík INRIA Paris & ENPC (project-team
More informationANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow
ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern
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 informationA reduced fracture model for two-phase flow with different rock types
A reduced fracture model for two-phase flow with different rock types Elyes Ahmed, Jérôme Jaffré, Jean E. Roberts To cite this version: Elyes Ahmed, Jérôme Jaffré, Jean E. Roberts. different rock types.
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 informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationPreconditioning the Newton-Krylov method for reactive transport
Preconditioning the Newton-Krylov method for reactive transport Michel Kern, L. Amir, A. Taakili INRIA Paris-Rocquencourt Maison de la Simulation MoMaS workshop Reactive Transport Modeling in the Geological
More informationAN UNFITTED METHOD FOR TWO-PHASE FLOW IN FRACTURED POROUS MEDIA
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 AN UNFIED MEHOD FOR WO-PHASE FLOW IN FRACURED POROUS MEDIA A. Fumagalli and A. Scotti
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
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
More informationReactive transport in porous media: formulations, non-linear solvers and preconditioners
Reactive transport in porous media: formulations, non-linear solvers and preconditioners Michel Kern with L. Amir, B. Gueslin, A. Taakili Institut National de Recherche en Informatique et Automatique High
More informationSchur Complement Technique for Advection-Diffusion Equation using Matching Structured Finite Volumes
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 8, Number 1, pp. 51 62 (213) http://campus.mst.edu/adsa Schur Complement Technique for Advection-Diffusion Equation using Matching Structured
More informationCoupling transport with geochemistry Applications to CO2 sequestration
Coupling transport with geochemistry Applications to CO2 sequestration Michel Kern with L. Amir, B. Gueslin, A. Taakili INRIA Paris Rocquencourt 3ème Conférence Internationale de la Société Marocaine de
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationMixed-hybrid finite element method for modelling two-phase flow in porous media
Journal of Math-for-Industry, Vol. 3 (211C-2), pp. 9 19 Mixed-hybrid finite element method for modelling two-phase flow in porous media Radek Fučík and Jiří Mikyška Revised on September 23, 211 Abstract.
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 informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
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 informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationA posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media
A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media D. A. Di Pietro, M. Vohraĺık, and S. Yousef Université Montpellier 2 Marseille,
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationParallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media
Parallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media Jakub Šístek1 joint work with Jan Březina 2 and Bedřich Sousedík 3 1 Institute of Mathematics of the AS CR Nečas Center
More informationA High Order Method for Three Phase Flow in Homogeneous Porous Media
A High Order Method for Three Phase Flow in Homogeneous Porous Media Justin Dong Advisor: B eatrice Rivi`ere February 5, 2014 Abstract The modeling of three-phase fluid flow has many important applications
More informationNONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS MEDIA APPLICATIONS
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira Eds Lisbon, Portugal,14-17 June 2010 NONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationUncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media
Uncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media M.Köppel, C.Rohde Institute for Applied Analysis and Numerical Simulation Inria, Nov 15th, 2016 Porous Media Examples: sponge,
More informationTIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA
XVIII International Conference on Water Resources CMWR 2 J. Carrera (Ed) c CIMNE, Barcelona, 2 TIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA Florian Haeberlein and Anthony Michel
More informationDownloaded 08/04/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 4, pp. B593 B618 c 2016 Society for Industrial and Applied Mathematics ACTIVE-SET REDUCED-SPACE METHODS WITH NONLINEAR ELIMINATION FOR TWO-PHASE FLOW PROBLEMS IN POROUS
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 informationSafety assessment for disposal of hazardous waste in abandoned underground mines
Safety assessment for disposal of hazardous waste in abandoned underground mines A. Peratta & V. Popov Wessex Institute of Technology, Southampton, UK Abstract Disposal of hazardous chemical waste in abandoned
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 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 informationFrom Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D
From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct
More informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationPRECONDITIONING A COUPLED MODEL FOR REACTIVE TRANSPORT IN POROUS MEDIA
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 16, Number 1, Pages 18 48 c 2019 Institute for Scientific Computing and Information PRECONDITIONING A COUPLED MODEL FOR REACTIVE TRANSPORT
More informationA Newton Krylov method for coupling transport with chemistry in porous media
A Newton Krylov method for coupling transport with chemistry in porous media Laila Amir 1,2 Michel Kern 2 1 Itasca Consultants, Lyon 2 INRIA Rocquencourt Joint work with J. Erhel (INRIA Rennes /IRISA)
More informationA Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach
SPE-172990-MS A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach Benjamin Ganis, Kundan Kumar, and Gergina Pencheva, SPE; Mary F. Wheeler, The
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
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 informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
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 informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationAgain we will consider the following one dimensional slab of porous material:
page 1 of 7 REVIEW OF BASIC STEPS IN DERIVATION OF FLOW EQUATIONS Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations,
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
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 finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
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 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 informationCoupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationSeminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011
Seminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011 Modeling Two-Phase-Two-Component Processes in Porous Media Using an Upscaling-Multi-Scale Method Elin Solberg 27 January 2011 Contents
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 informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationA MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS
A MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS JØRG E. AARNES AND YALCHIN EFENDIEV Abstract. methods are currently under active investigation for the simulation
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationThermo hydro mechanical coupling for underground waste storage simulations
Thermo hydro mechanical coupling for underground waste storage simulations Clément Chavant - Sylvie Granet Roméo Frenandes EDF R&D 1 Outline Underground waste storage concepts Main phenomena and modelisation
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 informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
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 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 informationEfficient and Accurate simulation of nonlinearly coupled multiphase flow in porous media. 25 August 2015 Meeneesh Kardale
Efficient and Accurate simulation of nonlinearly coupled multiphase flow in porous media 25 August 2015 Meeneesh Kardale Title : Efficient and accurate simulation of nonlinearly coupled multiphase flow
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 informationMortar upscaling for multiphase flow in porous media
Computational Geosciences 6: 73 100, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. Mortar upscaling for multiphase flow in porous media Małgorzata Peszyńska a, Mary F. Wheeler a andivanyotov
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationMMsFEM and Streamlines
MMsFEM and Streamlines Combining a mixed multiscale FEM with streamline simulation for enhanced reservoir performance prediction on large grid models. Jørg E. Aarnes, Vegard Kippe and Knut-Andreas Lie
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationOn efficient implicit upwind schemes
On efficient implicit upwind schemes Jostein R. Natvig Knut Andreas Lie June 30, 2008 Abstract Although many advanced methods have been devised for the hyperbolic transport problems that arise in reservoir
More informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationMulti-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media Jan Tecklenburg a,, Insa Neuweiler a, Jesus Carrera b, Marco Dentz b a Institute of Fluid Mechanics
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationElectronic Transactions on Numerical Analysis Volume 49, 2018
Electronic Transactions on Numerical Analysis Volume 49, 2018 Contents 1 Adaptive FETI-DP and BDDC methods with a generalized transformation of basis for heterogeneous problems. Axel Klawonn, Martin Kühn,
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationOptimized Schwarz methods with overlap for the Helmholtz equation
Optimized Schwarz methods with overlap for the Helmholtz equation Martin J. Gander and Hui Zhang Introduction For the Helmholtz equation, simple absorbing conditions of the form n iω were proposed as transmission
More informationA mixed finite element method for deformed cubic meshes
A mixed finite element method for deformed cubic meshes Inria Project Lab 2S@Exa, technical meeting at ANDRA, hâtenay-malabry Nabil Birgle POMDAPI project-team Inria Paris-Rocquencourt, UPM With : Martin
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationICES REPORT October Sanghyun Lee, Mary F. Wheeler, Thomas Wick
ICES REPORT 16-23 October 2016 Iterative coupling of flow, geomechanics and adaptive phase-field fracture including level-set crack width approaches by Sanghyun Lee, Mary F. Wheeler, Thomas Wick The Institute
More informationOn Application of the Weak Galerkin Finite Element Method to a Two-Phase Model for Subsurface Flow
J Sci Comput (2016) 66:225 239 DOI 10.1007/s10915-015-0021-8 On Application of the Weak Galerkin Finite lement Method to a Two-Phase Model for Subsurface Flow Victor Ginting 1 Guang Lin 2 Jiangguo Liu
More informationMixed and Mixed Hybrid Finite Element Methods: Theory, Implementation and Applications
Mixed and Mixed Hybrid Finite Element Methods: Theory, Implementation and Applications Max-Planck Institute for Mathematics in the Sciences Leipzig, Germany mailto:fradu@mis.mpg.de Introduction Objective
More informationThree-dimensional Modelling of Reactive Solutes Transport in Porous Media
151 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 41, 214 Guest Editors: Simonetta Palmas, Michele Mascia, Annalisa Vacca Copyright 214, AIDIC Servizi S.r.l., ISBN 978-88-9568-32-7; ISSN 2283-9216
More informationRecovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements
Recovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements Zhiqiang Cai Shun Zhang Department of Mathematics Purdue University Finite Element Circus, Fall 2008,
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationA study of the modelling error in two operator splitting algorithms for porous media flow
Computational Geosciences 2 (1998) 23 36 23 A study of the modelling error in two operator splitting algorithms for porous media flow K. Brusdal a,h.k.dahle a, K. Hvistendahl Karlsen a and T. Mannseth
More informationHigh performance computing for neutron diffusion and transport equations
High performance computing for neutron diffusion and transport equations Horizon Maths 2012 Fondation Science Mathématiques de Paris A.-M. Baudron, C. Calvin, J. Dubois, E. Jamelot, J.-J. Lautard, O. Mula-Hernandez
More information