Coupling Stokes and Darcy equations: modeling and numerical methods
|
|
- Randolf Barton
- 5 years ago
- Views:
Transcription
1 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), Marie Curie grant
2 Motivation We want to model the filtration of a free fluid through a porous medium: Environmental application Biomedical simulations [simulation by E. Miglio] [simulation by P. Zunino] Industrial applications: filters, porous foams, fuel cells... We must introduce a suitable modeling framework and set up effective numerical methods.
3 Modeling filtration (i) To model the filtration process three different approaches are generally considered: I - Microscale II - Mesoscale III - Macroscale [Chandesris and Jamet, Transp. Porous Med., 29]
4 Modeling filtration (ii) I - Microscale approach The geometry of the porous medium must be taken into account and the Navier-Stokes equations must be solved in the whole domain. [R. Hilfer, Lecture Notes in Physics 554, p (2)]
5 Modeling filtration (iii) II - Mesoscale approach A transition region is considered between the porous medium and the free fluid [Whitaker and Ochoa-Tapia (995); Jackson et al (22)]. Penalization approach [Iliev et al. (24,27); Angot (999)]: divσ(u,p)+(u )u+ nν K u+ C F K u u = f divu = in Ω f Ω d The penalization terms are set to zero in Ω f.
6 Modeling filtration (iv) III - Macroscale approach Sharp interface Γ with continuity conditions across it. Γ Ω f Ω d Free fluid Porous medium Free flow: Navier-Stokes equations divσ(u f,p f )+(u f )u f = f f divu f = in Ω f where σ(u f,p f ) = 2µ s u f p f I is the Cauchy stress tensor. Porous media flow: Darcy equation u d + K µ ( p ( ) d f d ) = K div divu d = µ ( p d f d ) = in Ω d
7 Modeling filtration (v) Continuity conditions (I) Continuity of normal velocities (due to incompressibility) u f n = u d n on Γ 2 Continuity of normal stresses n σ(u f,p f ) n = p d on Γ (pressures may be discontinuous across Γ: p f p d on Γ); 3 Beavers Joseph Saffman condition τ σ(u f,p f ) n = α d µ τ K τ u f τ [Beavers and Joseph (967); Saffman (97); Jäger and Mikelić (996,2,2) Miglio et al (22); Layton et al (23)]
8 Modeling filtration (vi) Continuity conditions (II) Continuity of normal velocities (due to incompressibility) u f n = u d n on Γ 2 Continuity of pressures p f = p d on Γ [Ene and Sanchez-Palencia (975); Levy and Sanchez-Palencia (975)]
9 Modeling filtration (vii) The microscale approach is very demanding: careful reconstruction of the geometry and solution of the Navier-Stokes equations. The macroscale approach allows using ad-hoc models in each subdomain. Once the coupling conditions are set, a global well-posed problem may be set up and non-overlapping domain decomposition methods can be used to solve the corresponding interface problem (substructuring techniques). The mesoscale approach requires solving a more expensive problem in Ω f Ω d than the macroscale approach.
10 Coupling through a sharp interface: domain decomposition methods without overlap (macroscale setting)
11 Weak form of the Stokes-Darcy problem Find u f H (Ω f ), p f L 2 (Ω f ), p d H (Ω d ): 2µ s u f : s α d µ v+ (u f ) τ v τ p f divv Ω f Γ τ K τ Ω f + p d (v n) = f f v Γ Ω f q f divu f = Ω f q d K Ω d µ p K d (u f n)q d = Γ Ω d µ f d q d
12 Algebraic form A suitable conforming finite element approximation of the coupled problem leads to the linear system: (A f +ᾱm Γ ) II (A f ) IΓ D T I (A f ) ΓI (A f ) ΓΓ D T Γ C T Γ D I D Γ C Γ A d (u f ) I (u f ) n p f = b p d with (u f ) n vector of the nodal values of u f n on Γ
13 An equivalent interface equation The Stokes-Darcy problem can be formulated in terms of the solution λ (normal velocity across Γ) of the interface problem S f fluid operator: S f : λ (normal velocities on Γ) S f λ+s d λ = χ on Γ solve n σ(u f,p f ) n (normal stresses on Γ). Stokes S d porous medium operator: S d : λ (normal velocities on Γ) solve Darcy p d Γ (pressure p d on Γ). Discrete interface equation for the normal velocity: (Σ f +Σ d )(u f ) n = χ s +χ d [MD, Quarteroni (23)]
14 Dimensional analysis of the Stokes-Darcy problem The Stokes-Darcy problem can be characterized in terms of three dimensionless numbers: the Reynolds number Re = U fx f ρ µ the Euler number Eu = P f ρu 2 f the Darcy number [Bear (99)] Da = K X 2 f [MD (in preparation)]
15 Dimensionless Schur complement system The Schur complement system with respect to (u f ) n becomes: (Σ s +Σ d )(u f ) n = χ s +χ d with Σ f = (ReEu) ˆΣf Σ d = (ReEuDa) ˆΣd and ˆΣ f = (A f ) ΓΓ (A f ) ΓI (A f +αda /2 M Γ ) II (A f IΓ ( ) + (A f ) ΓI (A f +αda /2 M Γ ) II D T I D T Γ ( ) D I (A f +αda /2 M Γ ) II D T I ( ) D I (A f +αda /2 M Γ ) II (A f ) IΓ D Γ ˆΣ d = C T Γ A d C Γ
16 If we multiply by (ReEuDa), we obtain: (Da ˆΣ f + ˆΣ d )(u f ) n = ˆχ f + ˆχ d Remark that we have independence of Re. We can prove that cond((da ˆΣ f + ˆΣ d )) c k + k Dak +h Dak + + k k+ h h Re = 2 Da: h = / Re = h = / h = / e+3.289e+3 h = / e e+3 h = / h = / e+3.289e+3 h = / h = / e e+3
17 Dirichlet-Neumann type preconditioners (i) It can be proved that cond((da ˆΣ f ) (ˆΣ f + ˆΣ d )) C + k Da + k h 2 k + k + Da from which we can see that if Da h 2, then cond((da ˆΣ f ) (ˆΣ f + ˆΣ d )) if Da h 2, then cond((da ˆΣ f ) (ˆΣ f + ˆΣ d )) h 2 PCG iterations, tol =, Re = Da: h = / h = / h = / h = / no prec. Da: h = / h = / h = / h = / P = (Da ˆΣ f )
18 Dirichlet-Neumann type preconditioners (ii) It can be proved that cond(ˆσ d (ˆΣ f + ˆΣ + k + k + Da k d )) C h 2 +k Da from which we can see that if Da h 2, then cond(ˆσ d (ˆΣ f + ˆΣ d )) h 2 if Da h 2, then cond(ˆσ d (ˆΣ f + ˆΣ d )) PCG iterations, tol =, Re = Da: h = / h = / h = / h = / no prec. Da: h = / h = / h = / h = / P = (ˆΣ d )
19 Comments Using (ˆΣ f ) requires solving a Stokes problem: y = ˆΣ f w y = u f n on Γ where Stokes(u f,p f ) = in Ω f n σ(u f,p f ) n = w on Γ Using (ˆΣ d ) requires solving a Darcy problem: y = ˆΣ d w y = K p d µ n div( K µ p d on Γ where ) = in Ω d p d = w on Γ
20 Numerical results (i) We consider the following setting proposed in [Hanspal et al, 29]: Ω f =.5.5 m, Ω d =.5.25 m - Inflow with max horizontal velocity. m/s, Re = 5 K: m 2 2 m 2 4 m 2 m 2 2 m 2 4 m 2 Da: h = / h = / h = / >25 >25 no prec. P = (Da ˆΣ f ) m 2 2 m 2 4 m P = (ˆΣ d )
21 Numerical results (ii) Following [Zunino (22)] we consider an artery with - radius =.3 cm - wall thickness =.3 cm - blood density ρ b =.4 g/cm 3 - blood (kinematic) viscosity ν b =.33 cm 2 /s - wall permeability K = 2 4 cm 2 This corresponds to - Re = Da = Dofs (u f + p f + p d ) h No prec. PCG with ˆΣ d / / / /
22 Numerical results (iii) Ω f = (,) (,2), Ω d = (,) (,) Re = Variable permeability 8 K 6 m 2 (Da = 3.3 ) Iterations Relative residual (PCG with ˆΣ d ) (tol = ) Velocity field Vertical velocity Hydrodynamic pressure x
23 Coupling through a transition region: interface control domain decomposition methods (mesoscale setting) Joint work with: Paola Gervasio (Università di Brescia) Alfio Quarteroni (EPF Lausanne)
24 The ICDD method for the Stokes-Darcy problem (i) Similarly to the mesoscale approach, we introduce a transition region overlapping Ω f and Ω d : Γ 2 Ω f Ωfd Γ Stokes divσ(u f,p f ) = f f in Ω f div u f = in Ω f Φ f (u f,p f ) = λ on Γ Ω d Darcy ( ) K div µ ( p d f d ) = in Ω d Φ d (p d ) = λ 2 on Γ 2 λ, λ 2 control interface functions
25 The ICDD method for the Stokes-Darcy problem (ii) Consider controls corresponding to the traces of the Stokes velocity and the Darcy pressure on the interfaces: u f = λ on Γ and p d = λ 2 on Γ 2. We set the minimization problem: inf λ,λ 2 J(λ,λ 2 ) with J(λ,λ 2 ) = 2 u f u d 2 L 2 (Γ ) + 2 p f p d 2 L 2 (Γ 2 ) [MD, Gervasio, Quarteroni (submitted 24)]
26 The ICDD method for the Stokes-Darcy problem (iii) The corresponding optimality system becomes: state problems: divσ(u f,p f ) = f f in Ω ( ) f K div u f = in Ω div f µ ( p d f d ) = in Ω d u f = λ on Γ p d = λ 2 on Γ 2 pseudo-adjoint problems: divσ(w f,q f ) = in Ω ( ) f K div w f = in Ω div f µ ψ d = in Ω d w f = u f u d on Γ ψ d = p d p f on Γ 2 pseudo-euler-lagrange equations: (u f u d )+ψ d = on Γ (p d p f )+q f = on Γ 2.
27 Numerical results (i) BiCGStab iterations: Dependence on K K iter inf J 9.6e-2 6.8e e-25.3e-25 7.e e-3.5e-24 u f u d L 2 (Ω fd ) 3.9e-6.7e-6.9e-6.9e-6.9e-6.9e-6.9e-6 p f p d L 2 (Ω fd ) 4.5e- 2.7e-9 4.7e-9 5.2e-9 5.3e-9 5.3e-9 5.3e-9 Dependence on δ δ iter (K = 8 ) 4 4 iter (K = 2 ) 2 2 Physical considerations lead to choose δ K. Dependence on h h iter (K = 8 ) iter (K = 2 ) 2 2 2
28 Comparison with the sharp-interface solution (K = 7 m 2 ) ICDD: Sharp interface (Beavers-Joseph-Saffman): Horizontal velocity Vertical velocity Hydrodynamic pressure
29 Comparison with the sharp-interface solution (K = 2 m 2 ) ICDD: Sharp interface (Beavers-Joseph-Saffman): Horizontal velocity Vertical velocity Hydrodynamic pressure
30 Comparison with the sharp-interface solution Stokes domain K = 7 m y 2 2 x y 2 x y K = 2 m 2 x y 2 2 x y 2 x y u ICDD u BJS u BJS v ICDD v BJS v BJS p ICDD p BJS p BJS x 2 3
31 Numerical results (ii) K =, K 2 = 4, K 3 = 6 m 2 δ =.5%, convergence in 4 iterations (tol = 9 ) ICDD: Velocity field Hydrodynamic pressure Sharp interface (Beavers-Joseph-Saffman):
32 Numerical results (iii) K = 7, K 2 = m 2, Re =. δ =.5%, convergence in 2 iterations (tol = 9 ) ICDD: Velocity field Hydrodynamic pressure Stokes Hydrodynamic pressure Darcy This case is not meaningful using Beavers-Joseph-Saffman: ICDD can treat more general fluid regimes.
33 Conclusions Sharp-interface approach: An interface problem and ad-hoc preconditioners for the Schur complement system (Steklov-Poincaré equation). Transition-region approach: A minimization problem which takes into account the physical properties of the problem without relying on any specific interface conditions.
34
35 Optimized Schwarz (Robin-Robin type) preconditioner (i) How to deal with intermediate values of Da? Considering the Schur complement of the linear system associated to the Stokes-Darcy problem with respect to both (u f ) n and p d Γ, we obtain the following augmented interface system: where Σ p Σ d. ( Σs C T )( ) Γ (uf ) n C Γ Σ p p d Γ ( ) χs = χ p [MD (23)]
36 Optimized Schwarz (Robin-Robin type) preconditioner (ii) Consider the Robin-Robin method: Solve the Stokes problem in Ω f with interface conditions on Γ: τ σ(u (k) f,p (k) f ) n = α dµ τ K τ u (k) f τ n σ(u (k) f,p (k) f ) n α f u (k) f n = p (k ) d α f u (k ) d n Solve the Darcy problem in Ω d with interface condition on Γ: p (k) d +α p u (k) d n = n σ(u (k) f,p (k) f ) n+α p u (k) f n where the parameters α f and α p must be chosen in a suitable way to ensure (optimal) convergence.
37 Optimized Schwarz (Robin-Robin type) preconditioner (iii) The Robin-Robin method can be proved to be equivalent to a Gauss-Seidel scheme for the augmented interface Schur complement system multiplied by the preconditioner: P α p I OS = ( I ) αf I I which can be used in the framework of Krylov methods like GMRES or BiCGStab. The theory of Optimized Schwarz methods allows to characterize α f and α p to maximize the convergence rate of the method. [MD, Gerardo-Giorda (in preparation)]
38 Numerical results (i) Domain Ω f = (,) (,2) and Ω d = (,) (,). Solution u f = ((y ) 2 +(y )+ Da,x(x )) p f = 2 (x +y ) ReEu ( p d = (x( x)(y )+(y ReEu Da )3 /3)+2x ) We set Re = 3 and Da = 6, 8,, 2 m 2. Using optimized Schwarz methods we can estimate α f = , α p =.528. Da: h = / h = / h = / BiCGStab iterations with tol =
39 .5.5 Stokes velocity field Numerical results (iii) Setting as in [Hanspal et al, 29] but with discontinuous permeability (reduced by one order of magnitude) and inflow from the left (no-slip conditions on the remaining fluid boundary): Darcy pressure Here Re = 5 and Da between.34 5 and Optimal coefficients: α f = 3.24 and α p =.24. h: 3/ 3/2 3/4 elements iterations BiCGStab iterations with tol =
40 The ICDD method ICDD is inspired to both the Virtual Control Method by J.L. Lions, O. Pironneau (998) and the Least Squares Method by R. Glowinski, Q. Dinh, J. Periaux (983). VCM/LSM ICDD J(λ) = 2 u u 2 L 2 (Ω 2 ) J(λ) = 2 u u 2 L 2 (Γ Γ 2 ) distributed observation interface observation classical Optimality System non-standard OS (integration by parts) Interface observation makes ICDD more suitable for heterogeneous problems and it leads to more efficient algorithms than distributed observation.
Iterative 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 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 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 informationKey words. domain decomposition, optimal control, elliptic boundary value problems, nonconforming discretizations, heterogeneous coupling, ICDD
THE INTERFACE CONTROL DOMAIN DECOMPOSITION (ICDD) METHOD FOR ELLIPTIC PROBLEMS MARCO DISCACCIATI, PAOLA GERVASIO, AND ALFIO QUARTERONI Abstract. Interface controls are unknown functions used as Dirichlet
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 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 domain decomposition framework for fluid-structure interaction problems
A domain decomposition framework for fluid-structure interaction problems Simone Deparis 1, Marco Discacciati 1, and Alfio Quarteroni 1,2 1 Institute of Analysis and Scientific Computing, CMCS-IACS, EPF
More informationAn Introduction to Fluid-Porous Interface Coupling
Bergische Universität Wuppertal Fachbereich Mathematik und Naturwissenschaften Lehrstuhl für Angewandte Mathematik und Numerische Mathematik Lehrstuhl für Optimierung und Approximation Preprint BUW-AMNA-OPAP
More informationCoupled free-flow and porous media flow: a numerical and experimental investigation
Coupled free-flow and porous media flow: a numerical and experimental investigation Master s Thesis Pavan Cornelissen 3863514 Supervisors: Kilian Weishaupt, MSc prof. dr. ir. Rainer Helmig prof. dr. ir.
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 informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationSome examples in domain decomposition
Some examples in domain decomposition Frédéric Nataf nataf@ann.jussieu.fr, www.ann.jussieu.fr/ nataf Laboratoire J.L. Lions, CNRS UMR7598. Université Pierre et Marie Curie, France Joint work with V. Dolean
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 informationON THE QUASISTATIC APPROXIMATION IN THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS
ON THE QUASSTATC APPROXMATON N THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS Marina Moraiti Department of Mathematics University of Pittsburgh Pittsburgh, PA 156, U.S.A. Tel: +1 (41) 889-9945
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 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 informationOn Limiting Behavior of Contaminant Transport Models in Coupled Surface and Groundwater Flows
Axioms 25, xx, -x; doi:.339/ Article OPEN ACCESS axioms ISSN 275-68 www.mdpi.com/journal/axioms On Limiting Behavior of Contaminant Transport Models in Coupled Surface and Groundwater Flows Vincent J.
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 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 informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
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 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 informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationPreprint Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction
Fakultät für Mathematik und Informatik Preprint 2015-15 Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction ISSN 1433-9307
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 45, No. 5, pp. 1801 1813 c 2007 Society for Industrial and Applied Mathematics A TWO-GRID METHOD OF A MIXED STOKES DARCY MODEL FOR COUPLING FLUID FLOW WITH POROUS MEDIA FLOW MO
More informationExtended variational formulation for heterogeneous partial differential equations
To appear in Computational Methods in Applied Mathematics (CMAM) Extended variational formulation for heterogeneous partial differential equations Pablo Blanco, Paola Gervasio, Alfio Quarteroni Abstract
More informationStokes-Darcy coupling for periodically curved interfaces
Stokes-Darcy coupling for periodically curved interfaces Sören Dobberschütz sdobber@nano.ku.dk February 21, 2014 Abstract We investigate the boundary condition between a free fluid and a porous medium,
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 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 informationMagnetic Drug Targeting in Cancer Therapy
Magnetic Drug Targeting in Cancer Therapy SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior
More informationMORTAR 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 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 informationDouble penalisation method for porous-fluid problems and numerical simulation of passive and active control strategies for bluff-body flows
Double penalisation method for porous-fluid problems and numerical simulation of passive and active control strategies for bluff-body flows Iraj Mortazavi Institut Polytechnique de Bordeaux INRIA Bordeaux
More informationNumerical Simulation of Newtonian and Non-Newtonian Flows in Bypass
Numerical Simulation of Newtonian and Non-Newtonian Flows in Bypass Vladimír Prokop, Karel Kozel Czech Technical University Faculty of Mechanical Engineering Department of Technical Mathematics Vladimír
More informationNumerical Analysis of Laminar flow of Viscous Fluid Between Two Porous Bounding walls
Numerical Analysis of Laminar flow of Viscous Fluid Between Two Porous Bounding walls Ramesh Yadav Department of Mathematics Babu Banarasi Das National Institute of Technology & Management Lucknow Uttar
More informationOne-Shot Domain Decomposition Methods for Shape Optimization Problems 1 Introduction Page 565
1 One-Shot Domain Decomposition Methods for Shape 2 Optimization Problems 3 Rongliang Chen 1 and Xiao-Chuan Cai 2 4 1 College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082,
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationPressure jump interface law for the Stokes-Darcy coupling: Confirmation by direct numerical simulations
Under consideration for publication in J. Fluid Mech. Pressure jump interface law for the Stokes-Darcy coupling: Confirmation by direct numerical simulations C a r r a r o T., G o l l C., M a r c i n i
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 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 informationReduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling Toni Lassila, Andrea Manzoni, Gianluigi Rozza CMCS - MATHICSE - École Polytechnique Fédérale
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationarxiv: v3 [physics.flu-dyn] 8 Nov 2016
Under consideration for publication in J. Fluid Mech. 1 arxiv:1604.02880v3 [physics.flu-dyn] 8 Nov 2016 A framework for computing effective boundary conditions at the interface between free fluid and a
More informationPrince Chidyagwai, Scott Ladenheim, and Daniel B. Szyld. Report July 2015
CONSTRAINT PRECONDITIONING FOR THE COUPLED STOKES-DARCY SYSTEM Prince Chidyagwai, Scott Ladenheim, and Daniel B. Szyld Report 5-07-9 July 05 This is a report of the Department of Mathematics at Temple
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationPrince Chidyagwai, Scott Ladenheim, and Daniel B. Szyld. Report July 2015, Revised October 2015
CONSTRAINT PRECONDITIONING FOR THE COUPLED STOKES-DARCY SYSTEM Prince Chidyagwai, Scott Ladenheim, and Daniel B. Szyld Report 5-07-9 July 05, Revised October 05 This is a report of the Department of Mathematics
More informationTectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop
.. Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Politecnico di Milano - MOX, Dipartimento di
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
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 informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
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 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 informationUNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS
UNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS MAXIM A. OLSHANSKII, JÖRG PETERS, AND ARNOLD REUSKEN Abstract. We consider
More information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationAn Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media
An Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media D.J. Lopez Penha a, B.J. Geurts a, S. Stolz a,b, M. Nordlund b a Dept. of Applied Mathematics, University of
More informationForchheimer law derived by homogenization of gas flow in turbomachines
Journal of Computational and Applied Mathematics 215 (2008) 467 476 www.elsevier.com/locate/cam Forchheimer law derived by homogenization of gas flow in turbomachines Gottfried Laschet ACCESS e.v., IntzestraYe
More informationINVERSE VISCOSITY BOUNDARY VALUE PROBLEM FOR THE STOKES EVOLUTIONARY EQUATION
INVERSE VISCOSITY BOUNDARY VALUE PROBLEM FOR THE STOKES EVOLUTIONARY EQUATION Sebastián Zamorano Aliaga Departamento de Ingeniería Matemática Universidad de Chile Workshop Chile-Euskadi 9-10 December 2014
More informationFluid flow over and through a regular bundle of rigid fibres
Under consideration for publication in J. Fluid Mech. 1 Fluid flow over and through a regular bundle of rigid fibres Giuseppe A. Zampogna and Alessandro Bottaro DICCA, Scuola Politecnica, Università di
More informationDomain Decomposition, Operator Trigonometry, Robin Condition
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03039-9 Domain Decomposition, Operator Trigonometry, Robin Condition Karl Gustafson 1. Introduction The purpose of this paper is to bring to the
More informationConstraint interface preconditioning for the incompressible Stokes equations Loghin, Daniel
Constraint interface preconditioning for the incompressible Stokes equations Loghin, Daniel DOI: 10.1137/16M1085437 License: Other (please specify with Rights Statement) Document Version Publisher's PDF,
More informationPhase-field modeling of a fluid-driven fracture in a poroelastic medium
Phase-field modeling of a fluid-driven fracture in a poroelastic medium Université Lyon 1, Institut Camille Jordan, FRANCE joint work with M.F. Wheeler (UT Austin), T. Wick (Ecole Polytechnique) Multiphysics,
More informationVerification of the boundary condition at the porous medium fluid interface
EPJ Web of Conferences 114, 02125 (2016) DOI: 10.1051/ epjconf/ 2016114 02125 C Owned by the authors, published by EDP Sciences, 2016 Verification of the boundary condition at the porous medium fluid interface
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationMéthode de quasi-newton en interaction fluide-structure. J.-F. Gerbeau & M. Vidrascu
Méthode de quasi-newton en interaction fluide-structure J.-F. Gerbeau & M. Vidrascu Motivations médicales Athérosclérose : Sistole Diastole La formation des plaques est liée, entre autre, la dynamique
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
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 informationFinite Elements for Magnetohydrodynamics and its Optimal Control
Finite Elements for Magnetohydrodynamics and its Karl Kunisch Marco Discacciati (RICAM) FEM Symposium Chemnitz September 25 27, 2006 Overview 1 2 3 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD)
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationIntroduction to Domain Decomposition Methods in the numerical approximation of PDEs. Luca Gerardo-Giorda
Introduction to Domain Decomposition Methods in the numerical approximation of PDEs Luca Gerardo-Giorda Plan 0 Motivation Non-overlapping domain decomposition methods Multi-domain formulation Variational
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationDIMENSIONAL MODEL REDUCTION FOR FLOW THROUGH FRACTURES IN POROELASTIC MEDIA
MOX-Report No. 45/6 DIMENSIONAL MODEL REDUCTION FOR FLOW THROUGH FRACTURES IN POROELASTIC MEDIA Bukac, M.; Yotov, I.; Zunino, P. MOX, Dipartimento di Matematica Politecnico di Milano, Via Bonardi 9-33
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationExtending the theory for domain decomposition algorithms to less regular subdomains
Extending the theory for domain decomposition algorithms to less regular subdomains Olof Widlund Courant Institute of Mathematical Sciences New York University http://www.cs.nyu.edu/cs/faculty/widlund/
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationMultilevel and Adaptive Iterative Substructuring Methods. Jan Mandel University of Colorado Denver
Multilevel and Adaptive Iterative Substructuring Methods Jan Mandel University of Colorado Denver The multilevel BDDC method is joint work with Bedřich Sousedík, Czech Technical University, and Clark Dohrmann,
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
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 informationKey words. Convection-Diffusion, Domain Decomposition, Preconditioning, Spectral Element Method
University of Maryland Department of Computer Science TR-C4931 University of Maryland Institute for Advanced Computer Studies TR-2009-03 FAST ITERATIVE SOLVER FOR CONVECTION-DIFFUSION SYSTEMS WITH SPECTRAL
More informationNumerical modelling of shear-thinning non-newtonian flows in compliant vessels
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2007; 00:1 [Version: 2002/09/18 v1.01] Numerical modelling of shear-thinning non-newtonian flows in compliant vessels M.
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners for
More information