Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels Mini-Workshop I
|
|
- Paul Harrington
- 5 years ago
- Views:
Transcription
1 Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels INRIA Rocquencourt, projet REO CRM Spring School, Montréal May 16, 2005
2 Motivations It is generally admitted that local blood flow behavior (shear stress) plays a key role in the understanding of several cardiovascular diseases Artery walls are compliant We aim at simulating the mechanical interaction between the blood flow and the arterial wall [press] [defor] [caro] [anev] Difficulties: Large displacements: geometrical non-linearities, moving computational domain Numerical stability: blood and wall densities are close Computational cost: strong research effort to improve the coupling algorithms FSI cost fluid cost + solid cost
3 Outline 1 Motivations 2 Part I: Model Equations The Structure Equations The Fluid Equations: ALE formulation The Coupled Problem Global Energy Balance 3 Part II: Explicit and Implicit Coupling Explicit Coupling Numerical Stability: Added Mass Effect Implicit Coupling Fixed-Point and Newton Methods 4 Part III: Semi-Implicit Coupling Remarks on Explicit and Implicit Coupling A New Coupling Scheme Numerical Stability
4 Part I Model Equations: Kinematics
5 Problem Setting Ω s (t) Γ in Γ w (t) Ω f (t) Γ out Ω f (t): fluid domain (filled with a viscous fluid, blood), Ω s (t): solid current configuration (arterial wall) Γ w (t): fluid-solid interface, where we enforce continuity of velocity equilibrium of stress
6 Problem Setting Ω s (t) Γ in Γ w (t) Ω f (t) Γ out Ω f (t): fluid domain (filled with a viscous fluid, blood), Ω s (t): solid current configuration (arterial wall) Γ w (t): fluid-solid interface, where we enforce continuity of velocity equilibrium of stress The problem Determine Ω f (t), velocity and stress within the fluid and the solid.
7 Solid Kinematics: Lagrange Formulation (I) Let Ω s 0 be reference material configuration, for instance Ωs 0 = Ωs (0). x s (, t) Ω s 0 Γ w 0 Γ w (t) Ω s (t) We introduce the solid deformation x s as an application that maps Ω s 0 into Ωs (t), i.e. x s (Ω s 0 ) = Ωs (t). x s : Ω s 0 R+ R 3. x def = x s (x 0, t) gives the position, at time t, of the material point x 0 Ω s 0. In other words, x s tracks the material points motion.
8 Solid Kinematics: Lagrange Formulation (II) We introduce the following notations: displacement: d s (x 0, t) def = x s (x 0, t) x 0 velocity: acceleration: d s def = ds t d s def = 2 d s t 2 deformation gradient: F (d s ) def = 0 x s = I + 0 d s jacobian: J(d s ) = det F (d s ) solid density in the reference configuration: ρ s (x 0 ) second Piola-Kirchoff stress tensor S(d s ), related to d s through a appropriate consitutive law. For instance, for an hyper-elastic material we have S(d s ) = W (E(ds )), E(d s ) def = 1 F (d s ) T F (d s ) I, E 2 with W a given elastic energy density and E(d s ) the Green-Lagrange stress tensor. unit normal vector on Ω s 0 : ns 0
9 Solid Equations: Lagrange Formulation From the Momentum conservation, we have that d s solves: ρ s ds div 0 (F (d s )S(d s )) = 0, in Ω s 0, Boundary conditions: d s = 0, on Γ D 0, F (d s )S(d s )n s 0 = 0, on ΓN 0, F (d s )S(d s )n s 0 = g fluid, on Γ w 0 (equilibrium of stress). Γ N 0 Γ D 0 Ω s 0 Γ w 0 Γ D 0 Γ N 0
10 Solid Equations: Variational Formulation Test functions space: V s def = H 1 Γ D (Ω s 0 ) def = 0 n v s H 1 (Ω s 0 ) : vs Γ D 0 o = 0 Solid variational formulation Find d s (t) V s such that ρ s d s v s + F (d s )S(d s ) : 0 d s = Ω s 0 Ω s 0 Γ w 0 g fluid v s, v s V s. Idea of the proof. Multiply by v s V s, integrate by parts and take into account the boundary conditions.
11 Fluid Equations: Some Notation We assume the fluid to be homogeneous, Newtonian and incompressible. We introduce the notations: density: ρ f dynamic viscosity: µ velocity: u(x, t), x Ω f (t) pressure: p(x, t) Cauchy stress tensor: σ(u, p), given by strain rate tensor: ε(u), defined as unit normal vector on Ω f (t): n unit normal vector on Ω f 0 : n0 σ(u, p) = pi + 2µε(u) ε(u) def = 1 h u + ( u) Ti 2
12 Fluid Equations: Eulerian Formulation (I) Γ in Ω f (t) Γ w (t) Γ out Incompressible Navier-Stokes equations in Eulerian formulation: «u ρ f t + u u 2µ div ε(u) + p = 0, in Ω f (t), Boundary conditions: σ(u, p)n = g in out, on Γ in out, div u = 0, in Ω f (t). u = d s, on Γ w (t) (continuity of velocities),
13 Fluid Equations: Eulerian Formulation (II) Γ w (t n+1 ) Γ w (t n ) j i i j Ω f (t n ) Remarks Ω f (t n+1 ) In the sequel we will deal with fluid meshes following the interface motion in order (accurate interface fluid load). In such case, the approximation of u (x, t), t becomes troublesome. This can be overcome by the introduction of a change of variables (the ALE map)
14 Fluid Domain Description (I) We set Ω f 0 def = Ω f (0) as a reference fluid domain. x f (, t) Ω f 0 Γ w 0 Ω f (t) Γ w (t) Assume we know the interface motion through the solid deformation x s Γ w. 0 We introduce the fluid domain deformation x f as an injective application x f : Ω f 0 R + R 3, such that x f Γ w 0 = x s Γ w 0, xf (Ω f 0 ) = Ωf (t). In other words x f = Ext(x s Γ w ), with Ext a given extension operator. 0 Remark Inside Ω f 0 the map xf is arbitrary
15 Fluid Domain Description (II) Thus, x f tracks the motion of the fluid domain. But, in general, it doesn t track the motion of the fluid particles. It is useful to introduce: the fluid domain displacement: d f (x 0, t) def = x f (x 0, t) x 0 Example of extension operator: the fluid domain velocity: w, defined by 8 d >< = 0, in Ω f 0, Harmonic lift: d f = d s, on Γ w 0 >: d f = 0, on Γ in out, w def = df t x 0 = xf t x 0 the deformation gradient: F (d f ) def = 0x f = I + 0d f and the jacobian: J(d f ) = det F (d f )
16 Fluid Domain Description (III) Remarks In general, w(x 0, t) u(x f (x 0, t), t) for x 0 Ω f 0. However, since by definition x f = x s on Γ w 0, we have and, therefore, w = u on Γ w 0 d f t x 0 = ds t x 0, on Γ w 0,
17 Time Derivatives: Eulerian vs. ALE (I) Let u = u(x, t) be a vector field defined on Ω f (t), t > 0. For x Ω f (t), the Arbitrary Lagrange-Euler (ALE) derivative of u is defined by with x 0 = ˆx f (, t) 1 (x), i.e. x = x f (x 0, t). u t x 0 (x, t) def = t u(xf (x 0, t), t), Proposition (ALE vs. Eulerian) Let x Ω f (t) with x = x f (x 0, t) and x 0 Ω f 0. Then with w def = d f. u u (x, t) = w(x 0, t) u(x, t) + (x, t), t x 0 t {z } {z } ALE derivative Eulerian derivative Idea of the proof. Use chain rule.
18 Why the terminology ALE? If w = 0, i.e. the fluid domain does not changes in time, we have u u (x, t) = (x, t), t x 0 t {z } {z } ALE derivative Eulerian derivative hence we recover the classical Eulerian derivative if w = u, i.e the map x f tracks the fluid particles motion, we have u (x, t) = u(x, t) u(x, t) + u (x, t), t x 0 t {z } {z } ALE derivative material or Lagrangian derivative hence we recover the classical Lagrangain or material derivative
19 Fluid Equations: ALE Formulation Γ in Ω f (t) Γ w (t) Γ out Incompressible Navier-Stokes equations in ALE formulation: «u ρ f + (w u) u 2µ div ε(u) + p = 0, t x 0 in Ω f (t), Boundary conditions: σ(u, p)n = g in, on Γ in σ(u, p)n = g out, on Γ in u = w, on Γ w (t) div u = 0, in Ω f (t).
20 Fluid Equations: ALE Variational Formulation (I) Test functions space o V f (t) = nv f : Ω f (t) R 3 : v f = ˆv [x f (, t)] 1, ˆv HΓ 1 w0 (Ωf 0 ) In particular, for each v f V f we have v f t x 0 = 0. Fluid variational formulation Find u(t) H 1 (Ω f (t)) with u Γ w (t) = w Γ w (t) and p(t) L 2 (Ω f (t)) such that d d t + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f u v f Ω f (t) ε(u) : ε(v f ) q div u = 0, ρ f (div w)u v f + Ω f (t) q L 2 (Ω f (t)), Ω f (t) p div v f = g v f, Γ in out ρ f (w u) u v f v f V f (t),
21 Fluid Equations: ALE Variational Formulation (II) Idea of the proof. By multiplying by v f V f, integrating by parts and taking into account the boundary conditions. Here we only consider the mass term: ρ f u v f = ρ f J(d f ) u v f t x 0 t x 0 Ω f (t) = Ω f 0 Ω f 0 = d d t = d d t ρ f J(d f ) v f J(d f ) ρ f u v f t x 0 Ω f t x 0 0 J(d f )ρ f u v f J(d f )(div w)ρ f u v f Ω f 0 Ω f (t) Ω f 0 ρ f u v f (div w)ρ f u v f. Ω f (t) where we used the fact that J(d f ) t x 0 = J(d f ) div w.
22 Fluid Load at the Interface: Variational Residual Proposition Let (u, p) solution of the previous fluid sub-problem, v s V s and L : H 1 2 (Γ w 0 ) H1 Ω f (Ω f 0 Γw 0 ), 0 a given lift operator. Let v f = L(v s ) [x f (, t)] 1 (obviously v f / V f (t)), then Γ w 0 g fluid v s = σ(u, p)n v s = Γ w (t) Idea of the proof. Integration by parts. Remarks: d ρ f u v f d t Ω f (t) + ρ f (w u) u v f + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f (div w)u v f ε(u) : ε(v f ) We solve the fluid equations with test functions vanishing on the interface. Ω f (t) We compute the forcing term with test functions not vanishing on the interface. This explains the terminology variational residual This approach has a straightforward discrete counterpart. p div v f.
23 The Coupled Problem: Strong Formulation Fluid: Solid: «u ρ f + (u w) u 2µ div ε(u) + p = 0, in Ω f (t), t x 0 div u = 0, in Ω f (t), σ(u, p)n = g in out, on Γ in out, ρ s d s div 0 (F (d s )S(d s )) = 0, in Ω s 0, Coupling conditions: geometry, velocity and stress d s = 0, on Γ D 0, F (d s )S(d s )n s 0 = 0, on ΓN 0, d f = Ext(d s Γ w 0 ), in Ωf 0, Ωf (t) = (I + d f )(Ω f 0 ), u = w(d f ), on Γ w (t), F (d s )S(d s )n 0 = J(d f )σ(u, p)f (d f ) T n 0, on Γ w 0.
24 The Coupled Problem: Variational Formulation Fluid domain motion: d f = Ext(d s Γ w 0 ), w = d f, Ω f (t) = (I + d f )(Ω f 0 ). Fluid: Find u(t) H 1 (Ω f (t)) with u Γ w (t) = w Γ w (t) and p(t) L 2 (Ω f (t)) such that d d t + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f u v f Ω f (t) ε(u) : ε(v f ) q div u = 0, ρ f (div w)u v f + Ω f (t) q L 2 (Ω f (t)), Ω f (t) p div v f = g v f, Γ in out ρ f (w u) u v f v f V f (t), Solid: Find d s (t) V s such that ρ s ds v s + F (d s )S(d s ) : 0 d s = Ω s 0 Ω s 0 Γ w 0 g fluid v s, v s V s.
25 The Coupled Problem: Variational Formulation Fluid domain motion: d f = Ext(d s Γ w 0 ), w = d f, Ω f (t) = (I + d f )(Ω f 0 ). Fluid: Find u(t) H 1 (Ω f (t)) with u Γ w (t) = w Γ w (t) and p(t) L 2 (Ω f (t)) such that d d t + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f u v f Ω f (t) ε(u) : ε(v f ) q div u = 0, ρ f (div w)u v f + Ω f (t) q L 2 (Ω f (t)), Ω f (t) p div v f = g v f, Γ in out ρ f (w u) u v f v f V f (t), Solid: Find d s (t) V s such that ρ s ds v s + F (d s )S(d s ) : 0 d s = R f`u, p, d f ; L(v s ), v s V s. Ω s 0 Ω s {z } 0 variational fluid residual
26 Global Energy Equality Proposition Assume that u = 0 on Γ in out and that the solid is an hyper-elastic material. Then, we have the following energy equality d ρ f ρ s! d t Ω f (t) 2 u 2 + Ω s 2 d s 2 + W (E(d s )) + 2µ ε(u) 2 = 0 0 Ω s 0 Ω f (t)
27 Global Energy Equality Proposition Assume that u = 0 on Γ in out and that the solid is an hyper-elastic material. Then, we have the following energy equality d d t ρ f ρ s Ω f (t) 2 u 2 + Ω s 2 d s 2 0 {z } Kinetic energy + Ω s 0 W (E(d s )) {z } Elastic energy «+ Ω f (t) 2µ ε(u) 2 = 0 {z } Dissipated power
28 Global Energy Equality Proposition Assume that u = 0 on Γ in out and that the solid is an hyper-elastic material. Then, we have the following energy equality d d t ρ f ρ s Ω f (t) 2 u 2 + Ω s 2 d s 2 0 {z } Kinetic energy + Ω s 0 W (E(d s )) {z } Elastic energy «+ Ω f (t) 2µ ε(u) 2 = 0 {z } Dissipated power Idea of the proof. Multiply the fluid and solid equations by u(t) and d s, respectively, integrate by parts and take into account the boundary conditions.
29 Part II Time Discretization: Explicit and Implicit Coupling
30 Time Discretization Fluid domain motion: d f = Ext(d s Γ w 0 ), w = d f, Ω f (t) = (I + d f )(Ω f 0 ). Fluid : Find u(t) H 1 (Ω f (t)) with u Γ w (t) = w Γ w (t) and p(t) L 2 (Ω f (t)) such that d d t + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f u v f Ω f (t) ε(u) : ε(v f ) q div u = 0, ρ f (div w)u v f + Ω f (t) q L 2 (Ω f (t)), Ω f (t) p div v f = g v f, Γ in out ρ f (w u) u v f v f V f (t), Solid : Find d s (t) V s such that ρ s ds v s + F (d s )S(d s ) : 0 d s = R f `u, p, d f ; L(v s ), v s V s. Ω s 0 Ω s 0
31 Time Discretization Fluid domain motion: d f,n+1 = Ext( d sn+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ). Fluid : Find u(t) H 1 (Ω f (t)) with u Γ w (t) = w Γ w (t) and p(t) L 2 (Ω f (t)) such that d d t + 2µ Ω f (t) Ω f (t) Ω f (t) ρ f u v f Ω f (t) ε(u) : ε(v f ) q div u = 0, ρ f (div w)u v f + Ω f (t) q L 2 (Ω f (t)), Ω f (t) p div v f = g v f, Γ in out ρ f (w u) u v f v f V f (t), Solid : Find d s (t) V s such that ρ s ds v s + F (d s )S(d s ) : 0 d s = R f `u, p, d f ; L(v s ), v s V s. Ω s 0 Ω s 0
32 Time Discretization Fluid domain motion: d f,n+1 = Ext( d sn+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ). Fluid (Implicit Euler): Find u n+1 H 1 (Ω f,n+1 ) with u n+1 Γw,n+1 = wn+1 Γw,n+1 and p n+1 L 2 (Ω f,n+1 ) such that 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out Solid : Find d s (t) V s such that ρ s ds v s + F (d s )S(d s ) : 0 d s = R f `u, p, d f ; L(v s ), v s V s. Ω s 0 Ω s 0
33 Time Discretization Fluid domain motion: d f,n+1 = Ext( d sn+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ). Fluid (Implicit Euler): Find u n+1 H 1 (Ω f,n+1 ) with u n+1 Γw,n+1 = wn+1 Γw,n+1 and p n+1 L 2 (Ω f,n+1 ) such that 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out Solid (Leap-Frog): Ω s 0 Find d s,n+1 V s such that ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
34 Coupling Strategies: Explicit/Implicit Key point: Give a prediction d sn+1 Γ w 0 of the interface displacement.
35 Coupling Strategies: Explicit/Implicit Key point: Give a prediction d sn+1 Γ w 0 of the interface displacement. Two possibilities:
36 Coupling Strategies: Explicit/Implicit Key point: Give a prediction d sn+1 Γ w 0 of the interface displacement. Two possibilities: Explicit coupling: d sn+1 Γ w 0 time-steps: d s,n, d s,n 1,... 0th-order interface prediction: 1st-order interface prediction:... given in terms of the interface displacement at the previous d sn+1 Γ w 0 = d s,n Γ w 0 d sn+1 Γ w 0 = d s,n Γ w. 0 d s,n Γ w 0 + d s,n 1 Γ w 0.
37 Coupling Strategies: Explicit/Implicit Key point: Give a prediction d sn+1 Γ w 0 of the interface displacement. Two possibilities: Explicit coupling: d sn+1 Γ w 0 time-steps: d s,n, d s,n 1,... 0th-order interface prediction: 1st-order interface prediction:... given in terms of the interface displacement at the previous d sn+1 Γ w 0 = d s,n Γ w 0 d sn+1 Γ w 0 Implicit coupling: No prediction, we simply set d sn+1 Γ w 0 = d s,n Γ w. 0 d s,n Γ w 0 + = d s,n+1 Γ w. 0 d s,n 1 Γ w 0.
38 Explicit Coupling: Detailed Algorithm 1 Compute fluid domain displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ). 2 Solve fluid in the new domain Ω f,n+1 : 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out u n+1 Γw,n+1 = wn+1 Γ w,n+1. 3 Solve solid with the new fluid forcing term R f Ω s 0 4 Next time step ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
39 Explicit Coupling (II) Some observations are in order:
40 Explicit Coupling (II) Some observations are in order: This algorithm is very attractive: FSI cost fluid cost + solid cost
41 Explicit Coupling (II) Some observations are in order: This algorithm is very attractive: FSI cost fluid cost + solid cost We have a fully uncoupled problem.
42 Explicit Coupling (II) Some observations are in order: This algorithm is very attractive: We have a fully uncoupled problem. FSI cost fluid cost + solid cost Velocity continuity is not exactly enforced (due to extrapolation) u n+1 Γ w 0 = d s,n Γ w 0 d s,n 1 Γ w 0 {z } extrapolated interface velocity d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity.
43 Explicit Coupling (II) Some observations are in order: This algorithm is very attractive: We have a fully uncoupled problem. FSI cost fluid cost + solid cost Velocity continuity is not exactly enforced (due to extrapolation) u n+1 Γ w 0 = d s,n Γ w 0 d s,n 1 Γ w 0 {z } extrapolated interface velocity d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity. Two major questions: Is energy well balanced at the discrete level?. Is this scheme stable?. [movie1] [movie2]
44 Explicit Coupling: A Priori Energy Estimate (I) Proposition Assume that u n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 + " + 1 W (E(d s,n+1 )) Ω s 0 Ω s 0 Ω s 0 ρ s 2 d s,n+1 d s,n 2 # W (E(d s,n )) q with q t x 0 = 0. ρ s Ω s µ ε(u n+1 ) 2 Ω f,n+1 d s,n d s,n 1 2 # R f `u n+1, p n+1, d f,n+1 ; L(u n+1 ds,n+1 d s,n ) 0
45 Explicit Coupling: A Priori Energy Estimate (I) Proposition Assume that u n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 + " + 1 W (E(d s,n+1 )) Ω s 0 Ω s 0 Ω s 0 ρ s 2 d s,n+1 d s,n 2 # W (E(d s,n )) q with q t x 0 = 0. ρ s Ω s µ ε(u n+1 ) 2 Ω f,n+1 d s,n d s,n 1 2 # R f `u n+1, p n+1, d f,n+1 ; L(u n+1 ds,n+1 d s,n ) 0 {z } Γw,n+1 σ(u n+1, p n+1 )n (u n+1 ds,n+1 d s,n )
46 Explicit Coupling: A Priori Energy Estimate (I) Proposition Assume that u n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds q with q t x 0 = 0. " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 ρ s + d s,n+1 d s,n 2 ρ s Ω s 2 d s,n d s,n 1 2 # 0 Ω s 2 0 " + 1 # W (E(d s,n+1 )) W (E(d s,n )) + 2µ ε(u n+1 ) 2 Ω s 0 Ω s 0 Ω f,n+1 j artificial R f `u n+1, p n+1, d f,n+1 ; L(u n+1 ds,n+1 d s,n ) 0 power {z } Γw,n+1 σ(u n+1, p n+1 )n (u n+1 ds,n+1 d s,n )
47 Explicit Coupling: A Priori Energy Estimate (II) Fluid domain displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), w n+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ). Fluid subproblem: Find u n+1 H 1 (Ω f,n+1 ) with u n+1 Γw,n+1 = wn+1 Γ and w,n+1 pn+1 L 2 (Ω f,n+1 ) such that 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n+1 Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out Solid subproblem: Find d s,n+1 V s such that Ω s 0 ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
48 Explicit Coupling: A Priori Energy Estimate (III) Idea of the proof. We take appropriate test functions and then we add the fluid and the solid variational formulations. In particular, for the solid discrete variational formulation: v s = ds,n+1 d s,n, Remark v s Γ w 0 = 0 and hence v s V s. for the fluid discrete variational formulation: v f = u n+1 L(v s Γ w 0 ) L(un+1 Γ w 0 v s Γ w 0 ) V f (t n+1 ). Remark 1 v f Γ w 0 = 0 and hence v f V f. 2 Let us recall that u n+1 Γ w 0 v Γ s w. 0
49 Explicit Coupling: Let s Summarize
50 Explicit Coupling: Let s Summarize Observations: Explicit coupling is cheap FSI cost fluid cost + solid cost
51 Explicit Coupling: Let s Summarize Observations: Explicit coupling is cheap FSI cost fluid cost + solid cost Numerical experiments (haemodynamics) show hat explicit coupling can be unstable Energy estimate with an artificial power term Γw,n+1 σ(u n+1, p n+1 )n (u n+1 ds,n+1 d s,n ).
52 Explicit Coupling: Let s Summarize Observations: Explicit coupling is cheap FSI cost fluid cost + solid cost Numerical experiments (haemodynamics) show hat explicit coupling can be unstable Energy estimate with an artificial power term Γw,n+1 σ(u n+1, p n+1 )n (u n+1 ds,n+1 d s,n ). Explicit coupling is stable, and widely used, in aeroelasticity (see Farhat, Piperno)
53 Explicit Coupling: Let s Summarize Observations: Explicit coupling is cheap FSI cost fluid cost + solid cost Numerical experiments (haemodynamics) show hat explicit coupling can be unstable Energy estimate with an artificial power term Γw,n+1 σ(u n+1, p n+1 )n (u n+1 ds,n+1 d s,n ). Explicit coupling is stable, and widely used, in aeroelasticity (see Farhat, Piperno) Major questions: What is the source of instabilities? Is there a hidden stability condition?
54 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out
55 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out Solid: string model (small displacements) ρ s h s η + Lη = p Γ w, in Γ w,
56 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out Solid: string model (small displacements) ρ s h s η + Lη = p Γ w, in Γ w, with η: vertical displacement h s : thickness L: linear operator (for instance Lη = aη b 2 η x )
57 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out Solid: string model (small displacements) ρ s h s η + Lη = p Γ w, in Γ w, Fluid: we keep a fixed fluid domain and we neglect viscous and convective effects 8 ρ f u >< t + p = 0, in Ωf, div u = 0, in Ω f, u n = η, on Γ >: w, p = p in out, on Γ in out.
58 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out Solid: string model (small displacements) ρ s h s η + Lη = p Γ w, in Γ w, Fluid: we keep a fixed fluid domain and we neglect viscous and convective effects 8 ρ f u >< t + p = 0, in 8 Ωf, p = 0, in Ω f, >< div u = 0, in Ω f, p u = = ρf u n = η, on Γ >: w, n t n = ρf η, on Γ w, >: p = p in out, on Γ in out. p = p in out, on Γ in out.
59 A 2D Simplified Model Γ w η Γ in Ω f symmetry axis Γ out Solid: string model (small displacements) ρ s h s η + Lη = p Γ w, in Γ w, Fluid: we keep a fixed fluid domain and we neglect viscous and convective effects 8 ρ f u >< t + p = 0, in 8 Ωf, p = 0, in Ω f, >< div u = 0, in Ω f, p u = = ρf u n = η, on Γ >: w, n t n = ρf η, on Γ w, >: p = p in out, on Γ in out. p = p in out, on Γ in out. Why a new model? Physics: reproduces propagation phenomena [movie] Numerics: explicit coupling unstable [movie]
60 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w,
61 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Definition (Inverse Steklov-Poincaré operator) The operator M A : H 1 2 (Γ w ) H 1 2 (Γ w ) is defined as follows: for each g H 1 2 (Γ w ) we define M A (g) = q Γ w, with q H 1 (Ω f ) the solution of 8 q = 0, in Ω f, >< q n = g, on Γw, >: q = 0, on Γ in out.
62 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Definition (Inverse Steklov-Poincaré operator) The operator M A : H 1 2 (Γ w ) H 1 2 (Γ w ) is defined as follows: for each g H 1 2 (Γ w ) we define M A (g) = q Γ w, with q H 1 (Ω f ) the solution of Using this definition we have: 8 q = 0, in Ω f, >< q n = g, on Γw, >: q = 0, on Γ in out. p Γ w = M A ( ρ f η).
63 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Condensed solid problem: we have p Γ w = M A ( ρ f η), thus ρ s h s η + Lη = M A ( ρ f η).
64 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Condensed solid problem: we have p Γ w = M A ( ρ f η), thus `ρs h s + ρ f M A η + Lη = 0, in Γ w. (2)
65 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Condensed solid problem: we have p Γ w = M A ( ρ f η), thus `ρs h s + ρ f M A {z } added mass η + Lη = 0, in Γ w. (2)
66 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Condensed solid problem: we have p Γ w = M A ( ρ f η), thus `ρs h s + ρ f M A {z } added mass η + Lη = 0, in Γ w. (2) Remarks: This new equation looks like a structure equation, except for the extra term ρ f M A The fluid-structure coupling can be condensed into an extra mass acting on the structure (hence the terminology added mass effect )
67 The Added Mass Effect The simplified coupled problem: Fluid: 8 p = 0, in Ω f, >< p n = ρf η, on Γ w, >: p = 0, on Γ in out. (1) Solid: ρ s h s η + Lη = p Γ w, in Γ w, Condensed solid problem: we have p Γ w = M A ( ρ f η), thus `ρs h s + ρ f M A {z } added mass η + Lη = 0, in Γ w. (2) Remarks: This new equation looks like a structure equation, except for the extra term ρ f M A The fluid-structure coupling can be condensed into an extra mass acting on the structure (hence the terminology added mass effect ) Major question: What kind of time integration scheme of (2) arises from the a explicit coupling of (1)?
68 Time Discretization: Explicit Coupling (I) Fluid: 8>< >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2.
69 Time Discretization: Explicit Coupling (I) Fluid: 8>< Solid: >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2. ρ s h s ηn+1 2η n + η n 1 () 2 + Lη n+1 = p n+1 Γ w.
70 Time Discretization: Explicit Coupling (I) Fluid: 8>< Solid: >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n Condensed solid problem: ρ s h s ηn+1 2η n + η n 1 () 2 p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2. ρ s h s ηn+1 2η n + η n 1 () 2 + Lη n+1 = p n+1 Γ w. + ρ f M A η n 2η n 1 + η n 2 () 2 + Lη n+1 = 0.
71 Time Discretization: Explicit Coupling (I) Fluid: 8>< Solid: >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n Condensed solid problem: p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2. ρ s h s ηn+1 2η n + η n 1 () 2 + Lη n+1 = p n+1 Γ w. ρ s h s ηn+1 2η n + η n 1 () 2 +ρ f η n 2η n 1 + η n 2 M A () 2 + Lη n+1 = 0. {z } implicit
72 Time Discretization: Explicit Coupling (I) Fluid: 8>< Solid: >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n Condensed solid problem: p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2. ρ s h s ηn+1 2η n + η n 1 () 2 + Lη n+1 = p n+1 Γ w. ρ s h s ηn+1 2η n + η n 1 () 2 +ρ f η n 2η n 1 + η n 2 M A () 2 +Lη n+1 = 0. {z } {z } implicit explicit
73 Time Discretization: Explicit Coupling (I) Fluid: 8>< Solid: >: + p n+1 = 0, 8 >< div u n+1 = 0, = >: u n+1 n = ηn η n 1. ρ f un+1 u n Condensed solid problem: p n+1 = 0, p n = ρf ηn 2η n 1 + η n 2 () 2. ρ s h s ηn+1 2η n + η n 1 () 2 + Lη n+1 = p n+1 Γ w. ρ s h s ηn+1 2η n + η n 1 () 2 +ρ f η n 2η n 1 + η n 2 M A () 2 +Lη n+1 = 0. {z } {z } implicit explicit Observation: Explicit coupling leads to a explicit discretization of the added mass term.
74 Times discretization: Explicit Coupling (II) Theorem (Causin-Gerbeau-Nobile 04) Let λ max be the for largest eigenvalue of M A and assume that Lη = aη. Then, the previous explicit coupling scheme is unconditionally unstable whenever ρ f λ max ρ s h s 1. (3)
75 Times discretization: Explicit Coupling (II) Theorem (Causin-Gerbeau-Nobile 04) Let λ max be the for largest eigenvalue of M A and assume that Lη = aη. Then, the previous explicit coupling scheme is unconditionally unstable whenever ρ f λ max ρ s h s 1. (3) Remarks: The stability of the coupling scheme does not depend on the time step. The instability condition (3) confirms empirical observations: instabilities may occur when the structure is light, as in haemodynamic applications ρ f ρ s in aeroelasticity ρ f ρ s, hence explicit coupling is stable
76 Times discretization: Explicit Coupling (II) Theorem (Causin-Gerbeau-Nobile 04) Let λ max be the for largest eigenvalue of M A and assume that Lη = aη. Then, the previous explicit coupling scheme is unconditionally unstable whenever ρ f λ max ρ s h s 1. (3) Remarks: The stability of the coupling scheme does not depend on the time step. The instability condition (3) confirms empirical observations: Idea of the proof. instabilities may occur when the structure is light, as in haemodynamic applications ρ f ρ s in aeroelasticity ρ f ρ s, hence explicit coupling is stable Note that M A is a compact, postive and self-adjoint operator on L 2 (Γ w ). Expand η n+1, η n, η n 1 and η n 2 on a orthonormal basis made of eigen-vectors of M A. Conclude by proving that characteristic polynomial of the resulting difference equation, has a root with magnitude larger than 1 under condition (3).
77 Time Discretization: Implicit Coupling Fluid: Solid: Ω s 0 d f,n+1 = Ext(d s,n+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ), 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out u n+1 Γw,n+1 = wn+1 Γ w,n+1 ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
78 Time Discretization: Implicit Coupling Fluid: d f,n+1 = Ext(d s,n+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ), 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out u n+1 Γw,n+1 = wn+1 Γ w,n+1 We formally note this problem as (d f,n+1, u n+1, p n+1 ) = F(d s,n+1 ). Solid: Ω s 0 ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
79 Time Discretization: Implicit Coupling Fluid: d f,n+1 = Ext(d s,n+1 Γ w ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 0 ), 1 ρ f u n+1 v f 1 ρ f u n v f + ρ f (div w n+1 )u n+1 v f Ω f,n+1 Ω f,n Ω f,n+1 + ρ f (w n+1 u n+1 ) u n+1 v f + 2µ ε(u n+1 ) : ε(v f ) Ω f,n+1 Ω f,n+1 p n+1 div v f = g n+1 v f, v f V f (t n+1 ), Ω f,n+1 Γ in out u n+1 Γw,n+1 = wn+1 Γ w,n+1 We formally note this problem as (d f,n+1, u n+1, p n+1 ) = F(d s,n+1 ). Solid: Ω s 0 ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 We formally note this problem as d s,n+1 = S(d f,n+1, u n+1, p n+1 ). = R f`u n+1, p n+1, d f,n+1 ; L(v s ), v s V s.
80 Explicit Coupling (II) Some observations are in order:
81 Explicit Coupling (II) Some observations are in order: Implicit coupling leads to stable schemes. Velocity continuity is exactly enforced u n+1 Γ w 0 = d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity We then recover an energy estimate without artificial power at the interface.
82 Explicit Coupling (II) Some observations are in order: Implicit coupling leads to stable schemes. Velocity continuity is exactly enforced u n+1 Γ w 0 = d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity We then recover an energy estimate without artificial power at the interface. Highly non-linear coupled problem (geometrical non-linearities): (d f,n+1, u n+1, p n+1 ) = F(d s,n+1 ) d s,n+1 = S(d f,n+1, u n+1, p n+1 )
83 Explicit Coupling (II) Some observations are in order: Implicit coupling leads to stable schemes. Velocity continuity is exactly enforced u n+1 Γ w 0 = d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity We then recover an energy estimate without artificial power at the interface. Highly non-linear coupled problem (geometrical non-linearities): Thus, by composition: (d f,n+1, u n+1, p n+1 ) = F(d s,n+1 ) d s,n+1 = S(d f,n+1, u n+1, p n+1 ) Fixed-point problem: d s,n+1 = S F(d s,n+1 )
84 Explicit Coupling (II) Some observations are in order: Implicit coupling leads to stable schemes. Velocity continuity is exactly enforced u n+1 Γ w 0 = d s,n+1 Γ w 0 d s,n Γ w 0 {z } interface velocity We then recover an energy estimate without artificial power at the interface. Highly non-linear coupled problem (geometrical non-linearities): Thus, by composition: (d f,n+1, u n+1, p n+1 ) = F(d s,n+1 ) d s,n+1 = S(d f,n+1, u n+1, p n+1 ) Fixed-point problem: d s,n+1 = S F(d s,n+1 ) Root-finding problem: R(d s,n+1 ) def = d s,n+1 S F(d s,n+1 ) = 0
85 Solution Methods: Fixed-Point Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Relaxation step: d s k+1 ω kd s k+1 + (1 ω k)d s k
86 Solution Methods: Fixed-Point Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Relaxation step: d s k+1 ω kd s k+1 + (1 ω k)d s k Acceleration via dynamic relaxation: Aitken s method (Mok-Wall-Ramm 01, Gerbeau-Vidrascu 03) ω k = (ds k ds k 1 ) (ds k S F(ds k ) ds k 1 + S F(ds k 1 )) d s k S F(ds k ) ds k 1 + S F(ds k 1 ) 2
87 Solution Methods: Fixed-Point Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Relaxation step: d s k+1 ω kd s k+1 + (1 ω k)d s k Acceleration via dynamic relaxation: Aitken s method (Mok-Wall-Ramm 01, Gerbeau-Vidrascu 03) ω k = (ds k ds k 1 ) (ds k S F(ds k ) ds k 1 + S F(ds k 1 )) d s k S F(ds k ) ds k 1 + S F(ds k 1 ) 2 Reducing the fluid computational cost via transpiration conditions (Deparis-MF-Formaggia 03)
88 Solution Methods: Fixed-Point Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Relaxation step: d s k+1 ω kd s k+1 + (1 ω k)d s k Acceleration via dynamic relaxation: Aitken s method (Mok-Wall-Ramm 01, Gerbeau-Vidrascu 03) ω k = (ds k ds k 1 ) (ds k S F(ds k ) ds k 1 + S F(ds k 1 )) d s k S F(ds k ) ds k 1 + S F(ds k 1 ) 2 Reducing the fluid computational cost via transpiration conditions (Deparis-MF-Formaggia 03) Observation Too expensive and may fail to converge!
89 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds
90 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds The bottleneck of Newton s method For a given solid variation z, how to compute ˆ D d s R(d s k ) z?
91 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds The bottleneck of Newton s method For a given solid variation z, how to compute ˆ D d s R(d s k ) z? Finite difference Jacobian approximations (Matthies-Steindorf 01, Heil 03), ˆ Dd s R(d s k ) z R(ds k + ɛz) R(ds k ), 0 < ɛ 1 ɛ
92 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds The bottleneck of Newton s method For a given solid variation z, how to compute ˆ D d s R(d s k ) z? Finite difference Jacobian approximations (Matthies-Steindorf 01, Heil 03), too expensive and may fail to converge!
93 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds The bottleneck of Newton s method For a given solid variation z, how to compute ˆ D d s R(d s k ) z? Finite difference Jacobian approximations (Matthies-Steindorf 01, Heil 03), too expensive and may fail to converge! Physically based Jacobian approximations (Tezduyar 01, Gerbeau-Vidrascu 03, Deparis-Gerbeau-Vaseur 04),
94 Solution Methods: Newton Based Methods 1 Initial guess d s 0, ds def 0 = d s,n 2 For k = 0, 1,... until convergence (a) Solve fluid and solid: d s k+1 = S F(ds k ) (b) Evaluate residual: R(d s k ) def = d s k S F(ds k ) = ds k ds k+1 (c) Solve tangent problem ˆ Dd s R(d s k ) δd s = R(d s k ) (solved via GMRES) (d) Update rule: d s k+1 ds k + δds The bottleneck of Newton s method For a given solid variation z, how to compute ˆ D d s R(d s k ) z? Finite difference Jacobian approximations (Matthies-Steindorf 01, Heil 03), too expensive and may fail to converge! Physically based Jacobian approximations (Tezduyar 01, Gerbeau-Vidrascu 03, Deparis-Gerbeau-Vaseur 04), cheaper but still may fail to converge!
95 Exact Jacobian Evaluation Given a solid displacement z we need to evaluate ˆ D d s R( d s ) z with R(d s ) def = d s S F(d s ).
96 Exact Jacobian Evaluation Given a solid displacement z we need to evaluate ˆ D d s R( d s ) z with R(d s ) def = d s S F(d s ). One application of the chain rule yields ˆ Dd s R( d s ) z = z ˆ D (d s,u,p) S(F( d s )) ˆ D d s F( d s ) z.
97 Exact Jacobian Evaluation Given a solid displacement z we need to evaluate ˆ D d s R( d s ) z with R(d s ) def = d s S F(d s ). One application of the chain rule yields ˆ Dd s R( d s ) z = z ˆ D (d f,u,p) S(F( d s )) ˆ D d s F( d s ) z {z } δu f = (δu, δp, δd f ).
98 Exact Jacobian Evaluation Given a solid displacement z we need to evaluate ˆ D d s R( d s ) z with R(d s ) def = d s S F(d s ). One application of the chain rule yields ˆ Dd s R( d s ) z = z ˆ D (d f,u,p) S(F( d s )) ˆ D d s F( d s ) z {z } δu f = (δu, δp, δd f ) {z } δz.
99 Fluid Tangent Problem Using shape-derivative calculus (Sokolowski-olesio 91, Delfour-olesio 01), (δu, δp) solves the following linear problem: 8 ρ h δu + ρ div δu (u n w g( d i f )) σ(δu, δp) = ρ (div δdf )ū nh div ρū (u n w g( d i f )) σ(ū, p) hi div δd f ( δd f ) Tio >< + ρ n h div(ū δdf ) div µ ū δd f + ( δd f ) T ( ū) Tio, in Ω f ( d), div δu = div n ū h I div δd f ( δd f ) Tio, in Ω f ( d f ), δu = δdf, on Γw ( d f ), h >: σ(δu, δp)n = µ ū δd f + ( δd f ) T ( ū) Ti n, on Γ in out, with δd f = Ext ( d s )z (MF-Moubachir 03,04).
100 Fluid Tangent Problem Approximations An approximate fluid tangent problem can be derived by neglecting the shape derivative terms, yielding 8 ρ h δu + ρ div δu (u n w g( d i f )) σ(δu, δp) = 0, >< div δu = 0, in Ω f ( d f ), in Ω f ( d f ), δu = δdf, on Γw ( d f ), >: σ(δu, δp)n = 0, on Γ in out.
101 Fluid Tangent Problem Approximations An approximate fluid tangent problem can be derived by neglecting the shape derivative terms, yielding 8 ρ h δu + ρ div δu (u n w g( d i f )) σ(δu, δp) = 0, in Ω f ( d f ), >< div δu = 0, in Ω f ( d f ), δu = δdf, on Γw ( d f ), >: σ(δu, δp)n = 0, on Γ in out. On the other hand, neglecting here the convective and diffusive terms (Gerbeau-Vidrascu 2003), we get the following problem 8 ρ δu + δp = 0, in Ωf ( d f ), >< div δu = 0, in Ω f ( d f ), δu = δdf, on Γw ( d f ), >: δpn = 0, on Γ in out ( d f ).
102 Iterations History ( = 10 4 s) Newton: exact Newton: not exact FP-Aitken number of iterations time
103 Computational Cost ( = 10 4 s) Algorithm CPU time (dimensionless) Picard-Aitken 1.00 Quasi-Newton 0.55 Newton 0.60
104 Computational Cost ( = 10 4 s) Algorithm CPU time (dimensionless) Picard-Aitken 1.00 Quasi-Newton 0.55 Newton 0.60 Question What about increasing the time step?
105 Residuals History ( = 10 3 s) 0.01 Newton: exact Newton: not exact residual 1e-05 1e-06 1e-07 1e number of iterations
106 Some Remarks on Implicit Schemes Newton based methods perform better than fixed-point based methods Approximated Jacobians may lead to divergent quasi-newton methods Exact Jacobians may improve the convergence of the Newton s loop (moderate time steps) but are too expensive Realistic simulations with this kind of methods are very expensive
107 Part III Semi-Implicit Coupling
108 Let s Summarize Explicit coupling is cheap but unstable
109 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive
110 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive Added-mass effect seems to be the source of numerical instabilities in explicit coupling
111 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive Added-mass effect seems to be the source of numerical instabilities in explicit coupling Geometrical non-linearities (moving domains), convective and viscous effects do not seem to affect the stability of a coupling algorithm. However, they are implicitly treated in fully implicit schemes (very expensive!)
112 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive Added-mass effect seems to be the source of numerical instabilities in explicit coupling Geometrical non-linearities (moving domains), convective and viscous effects do not seem to affect the stability of a coupling algorithm. However, they are implicitly treated in fully implicit schemes (very expensive!) Major question: How to balance stability and computational cost?
113 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive Added-mass effect seems to be the source of numerical instabilities in explicit coupling Geometrical non-linearities (moving domains), convective and viscous effects do not seem to affect the stability of a coupling algorithm. However, they are implicitly treated in fully implicit schemes (very expensive!) Major question: How to balance stability and computational cost? First idea: semi-implicit coupling Treat implicitly the added-mass effect Treat explicitly the geometrical non-linearities and the convective and viscous effects
114 Let s Summarize Explicit coupling is cheap but unstable Implicit coupling is stable but too expensive Added-mass effect seems to be the source of numerical instabilities in explicit coupling Geometrical non-linearities (moving domains), convective and viscous effects do not seem to affect the stability of a coupling algorithm. However, they are implicitly treated in fully implicit schemes (very expensive!) Major question: How to balance stability and computational cost? First idea: semi-implicit coupling Treat implicitly the added-mass effect Treat explicitly the geometrical non-linearities and the convective and viscous effects Second idea: semi-implicit coupling Perform this using a projection scheme (Chorin-Teman) within the fluid
115 The Chorin-Teman Projection Scheme Main feature Incompressibility and viscous/convective effects are decoupled Set u 0 = u 0, for n 0 compute (ũ n+1, u n+1, p n+1 ) Elliptic step: Projection step: 8 >< ρ f ũn+1 u n >: + ũ n+1 ũ n+1 «2µ div ε(ũ n+1 ) = f n+1, 8 >< ρ f un+1 ũ n+1 >: + p n+1 = 0, u n+1 = 0, u n+1 n Γ = 0. (see Guermond-Minev-Shen 04, Rannacher 91, Shen 92) ũ n+1 Γ = 0.
116 The Chorin-Teman Projection Scheme Main feature Incompressibility and viscous/convective effects are decoupled Set u 0 = u 0, for n 0 compute (ũ n+1, u n+1, p n+1 ) Elliptic step: 8 >< ρ f ũn+1 u n >: + ũ n+1 ũ n+1 «2µ div ε(ũ n+1 ) = f n+1, ũ n+1 Γ = 0. Projection step: 8 >< ρ f un+1 ũ n+1 >: + p n+1 = 0, u n+1 = 0, u n+1 n Γ = 0. 8 >< p n+1 = ρf ũn+1, = p >: n+1 = 0. n Γ (see Guermond-Minev-Shen 04, Rannacher 91, Shen 92)
117 Semi-Implicit Coupling Algorithm: Explicit Part Update fluid domain (mesh) displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ),
118 Semi-Implicit Coupling Algorithm: Explicit Part Update fluid domain (mesh) displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ), Fluid: Elliptic step in the new domain Ω f,n+1 : ρ f ũn+1 u n + (ũ n+1 w n+1 ) ũ n+1 «2µ div ε(ũ n+1 ) = 0, ũ n+1 Γ w,n+1 = wn+1 Γ w,n+1.
119 Semi-Implicit Coupling Algorithm: Explicit Part Update fluid domain (mesh) displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ), Fluid: Elliptic step in the new domain Ω f,n+1 : 1 ρ f ũ n+1 v f Ω f,n ρ f u n v f Ω f,n 1 + ρ f (div w n+1 )ũ n+1 v f Ω f,n ρ f (w n+1 ũ n+1 ) ũ n+1 v f + 2µ ε(ũ n+1 ) : ε(v f Ω f,n+1 Ω f,n+1 1 ) = 0, vf 1 V f (t n+1 ), ũ n+1 Γ w,n+1 = wn+1 Γ w,n+1.
120 Semi-Implicit Coupling Algorithm: Explicit Part Update fluid domain (mesh) displacement and velocity: d f,n+1 = Ext(d s,n Γ w 0 ), wn+1 = df,n+1 d f,n, Ω f,n+1 = (I + d f,n+1 )(Ω f 0 ), Fluid: Elliptic step in the new domain Ω f,n+1 : 1 ρ f ũ n+1 v f Ω f,n ρ f u n v f Ω f,n 1 + ρ f (div w n+1 )ũ n+1 v f Ω f,n ρ f (w n+1 ũ n+1 ) ũ n+1 v f + 2µ ε(ũ n+1 ) : ε(v f Ω f,n+1 Ω f,n+1 1 ) = 0, vf 1 V f (t n+1 ), ũ n+1 Γ w,n+1 = wn+1 Γ w,n+1. Observations Geometrical non-linearities, viscous and convective effects are out of inner iterations
121 Semi-Implicit Coupling Algorithm: Implicit Part Fluid: Projection step in Ω f,n+1 (known domain) 8 >< ρ f un+1 ũ n+1 + p n+1 = 0, in Ω f,n+1 >: u n+1 = 0, in Ω f,n+1, Solid: u n+1 n = ds,n+1 d s,n n, on Γ w,n+1. 8 >< ρ s ds,n+1 2d s,n + d s,n 1 () 2 div 0 `F (d s,n+1 )S(d s,n+1 ) = 0, in Ω s 0, >: F (d s,n+1 )S(d s,n+1 )n 0 = J(d f,n+1 )σ(ũ n+1, p n+1 )F (d f,n+1 ) T n 0, on Γ w 0.
122 Semi-Implicit Coupling Algorithm: Implicit Part Fluid: Projection step in Ω f,n+1 (known domain) 1 ρ f u n+1 v f Ω f,n ρ f ũ n+1 v f Ω f,n+1 2 p n+1 div v f Ω f,n q div u n+1 = 0, Ω f,n+1 v f 2 H Γ w,n+1(div; Ωf,n+1 ), q L 2 (Ω f,n+1 ), Solid: u n+1 n = ds,n+1 d s,n n, on Γ w,n+1. 8 >< ρ s ds,n+1 2d s,n + d s,n 1 () 2 div 0 `F (d s,n+1 )S(d s,n+1 ) = 0, in Ω s 0, >: F (d s,n+1 )S(d s,n+1 )n 0 = J(d f,n+1 )σ(ũ n+1, p n+1 )F (d f,n+1 ) T n 0, on Γ w 0.
123 Semi-Implicit Coupling Algorithm: Implicit Part Fluid: Projection step in Ω f,n+1 (known domain) 1 ρ f u n+1 v f Ω f,n ρ f ũ n+1 v f Ω f,n+1 2 p n+1 div v f Ω f,n q div u n+1 = 0, Ω f,n+1 Solid: Ω s 0 v f 2 H Γ w,n+1(div; Ωf,n+1 ), q L 2 (Ω f,n+1 ), u n+1 n = ds,n+1 d s,n n, on Γ w,n+1. ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = p n+1 v s n 2µ ε(ũ n+1 )n v s, v s V s. Γ w,n+1 Γ w,n+1
124 Semi-Implicit Coupling Algorithm: Implicit Part Fluid: Projection step in Ω f,n+1 (known domain) 1 ρ f u n+1 v f Ω f,n ρ f ũ n+1 v f Ω f,n+1 2 p n+1 div v f Ω f,n q div u n+1 = 0, Ω f,n+1 Solid: Ω s 0 v f 2 H Γ w,n+1(div; Ωf,n+1 ), q L 2 (Ω f,n+1 ), u n+1 n = ds,n+1 d s,n n, on Γ w,n+1. ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f p`pn+1 ; L 2 (v s n ) Γ w R f u`ũn+1, d f,n+1 ; L 0 1 (v s ) Γ {z } w, v s V s. 0 {z } Projection residual Elliptic residual
125 Semi-Implicit Coupling Algorithm: Implicit Part Fluid: Projection step in Ω f,n+1 (known domain) 1 ρ f u n+1 v f Ω f,n ρ f ũ n+1 v f Ω f,n+1 2 p n+1 div v f Ω f,n q div u n+1 = 0, Ω f,n+1 Solid: Ω s 0 v f 2 H Γ w,n+1(div; Ωf,n+1 ), q L 2 (Ω f,n+1 ), u n+1 n = ds,n+1 d s,n n, on Γ w,n+1. ρ s ds,n+1 2d s,n + d s,n 1 () 2 v s + F (d s,n+1 )S(d s,n+1 ) : 0 d s Ω s 0 = R f p`pn+1 ; L 2 (v s n ) Γ w R f u`ũn+1, d f,n+1 ; L 0 1 (v s ) Γ {z } w, v s V s. 0 {z } Projection residual Elliptic residual Observations Projection sub-step in a fixed fluid domain (fixed matrix) Linear coupled problem if linear solid This coupled problem requires cheap inner iterations
126 Semi-Implicit Coupling: A Priori Energy Estimate (I) Proposition Assume that ũ n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 + " + 1 W (E(d s,n+1 )) Ω s 0 Ω s 0 Ω s 0 ρ s 2 d s,n+1 d s,n 2 # W (E(d s,n )) q with q t x 0 = 0. ρ s Ω s µ ε(ũ n+1 ) 2 Ω f,n+1 d s,n d s,n 1 2 # R f `ũn+1 u, d f,n+1 ; L 2 (ũ n+1 ds,n+1 d s,n ) 0
127 Semi-Implicit Coupling: A Priori Energy Estimate (I) Proposition Assume that ũ n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 + " + 1 W (E(d s,n+1 )) Ω s 0 Ω s 0 Ω s 0 ρ s 2 d s,n+1 d s,n 2 # W (E(d s,n )) q with q t x 0 = 0. ρ s Ω s µ ε(ũ n+1 ) 2 Ω f,n+1 d s,n d s,n 1 2 # R f `ũn+1 u, d f,n+1 ; L 1 (ũ n+1 ds,n+1 d s,n ) 0 {z } Γw,n+1 2µ ε(ũ n+1 )n (ũ n+1 ds,n+1 d s,n )
128 Semi-Implicit Coupling: A Priori Energy Estimate (I) Proposition Assume that ũ n+1 = 0 on Γ in out, the structure is an hyper-elastic material and» 1 q q = div w n+1 q, Ω f,n+1 Ω f,n Ω f,n+1 Then, the following energy inequality holds q with q t x 0 = 0. " 1 ρ f Ωf,n Ω f,n+1 2 un+1 2 ρ f 2 un 2 ρ s + d s,n+1 d s,n 2 ρ s Ω s 2 d s,n d s,n 1 2 # 0 Ω s 2 0 " + 1 # W (E(d s,n+1 )) W (E(d s,n )) + 2µ ε(ũ n+1 ) 2 Ω s 0 Ω s 0 Ω f,n+1 j artificial R f `ũn+1 viscous power u, d f,n+1 ; L 1 (ũ n+1 ds,n+1 d s,n ) 0 {z } Γw,n+1 2µ ε(ũ n+1 )n (ũ n+1 ds,n+1 d s,n )
129 Semi-Implicit Coupling: A Priori Energy Estimate (II) Idea of the proof. We take appropriate test functions and then we add the fluid and the solid variational formulations. In particular, for the solid: for the fluid (explicit part): for the fluid (implicit part): v f 1 = ũn+1 L 1 (v s Γ w 0 v s = ds,n+1 d s,n ) L 1 (ũ n+1 Γ w 0 V s. v Γ s w ) V f (t n+1 ). 0 v f 2 = un+1 L 2 (v s n Γ w 0 ) H Γ w,n+1(div; Ω f,n+1 ). Remarks: 1 v f 1 Γ w 0 = 0 and hence v f V f 2 v f 2 n Γ w 0 = 0 and hence v f 2 H Γ w,n+1(div; Ωf,n+1 ) 3 Let us recall that ũ n+1 Γ w 0 v s Γ w 0
Mé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 informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationFluid-structure algorithms based on Steklov-Poincaré operators
Fluid-structure algorithms based on Steklov-Poincaré operators S. Deparis, M. Discacciati, G. Fourestey, A. Quarteroni April 25, 2005 Abstract In this paper we review some classical algorithms for fluid-structure
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 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 informationA semi-implicit approach for fluid-structure interaction based on an algebraic fractional step method
A semi-implicit approach for fluid-structure interaction based on an algebraic fractional step method Annalisa Quaini, Alfio Quarteroni 3rd October 6 Chair of Modelling and Scientific Computing (CMCS),
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationA modular, operator splitting scheme for fluid-structure interaction problems with thick structures
Numerical Analysis and Scientific Computing Preprint Seria A modular, operator splitting scheme for fluid-structure interaction problems with thick structures M. Bukač S. Čanić R. Glowinski B. Muha A.
More informationStability of the kinematically coupled β-scheme for fluid-structure interaction problems in hemodynamics
Numerical Analysis and Scientific Computing Preprint Seria Stability of the kinematically coupled β-scheme for fluid-structure interaction problems in hemodynamics S. Čanić B. Muha M. Bukač Preprint #13
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationAn introduction to the Immersed boundary method and its finite element approximation
FE An introduction to te and its finite element approximation Collaborators: Daniele offi, Luca Heltai and Nicola Cavallini, Pavia Pavia, 27 maggio 2010 FE History Introduced by Peskin Flow patterns around
More informationFourier Series Development for Solving a Steady Fluid Structure Interaction Problem
Int. J. Contemp. Math. Sciences, Vol. 8, 203, no. 2, 75-84 HIKARI Ltd, www.m-hikari.com Fourier Series Development for Solving a Steady Fluid Structure Interaction Problem Ibrahima Mbaye University of
More informationarxiv: v2 [math.ap] 4 Mar 2016
Existence of a weak solution to a fluid-elastic structure interaction problem with the Navier slip boundary condition Boris Muha Sunčica Čanić arxiv:155.446v [math.ap] 4 Mar 16 Abstract We study a nonlinear,
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 informationA QUASI-NEWTON ALGORITHM BASED ON A REDUCED MODEL FOR FLUID-STRUCTURE INTERACTION PROBLEMS IN BLOOD FLOWS
ESAIM: M2AN Vol. 37, N o 4, 2003, pp. 631 647 DOI: 10.1051/m2an:2003049 ESAIM: Mathematical Modelling and Numerical Analysis A QUASI-NEWTON ALGORITHM BASED ON A REDUCED MODEL OR LUID-STRUCTURE INTERACTION
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationStability and convergence analysis of the kinematically coupled scheme for the fluid-structure interaction
Stability and convergence analysis of the kinematically coupled scheme for the fluid-structure interaction Boris Muha Department of Mathematics, Faculty of Science, University of Zagreb 2018 Modeling,
More informationAlgorithms for fluid-structure interaction problems
9 Algorithms or luid-structure interaction problems Miguel A. Fernández, Jean-Frédéric Gerbeau Fluid-structure coupled problems and their analysis have been addressed in the previous chapters particularly
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
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 informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationAdvances in the mathematical theory of the finite element immersed boundary method
Advances in the mathematical theory of the finite element immersed boundary method Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 12, 2014
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 informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationExtrapolation methods for accelerating unsteady partitioned fluid-structure interaction simulations
Int. Jnl. of Multiphysics Volume 5 Number 4 211 287 Extrapolation methods for accelerating unsteady partitioned fluid-structure interaction simulations Stephen Sachs*, Marcel Streitenberger, Dörte C. Sternel**,
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationSTABILITY OF THE KINEMATICALLY COUPLED β-scheme FOR FLUID-STRUCTURE INTERACTION PROBLEMS IN HEMODYNAMICS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Number, Pages 54 8 c 25 Institute for Scientific Computing and Information STABILITY OF THE KINEMATICALLY COUPLED β-scheme FOR FLUID-STRUCTURE
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 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 informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationFSI with Application in Hemodynamics Analysis and Simulation
FSI with Application in Hemodynamics Analysis and Simulation Mária Lukáčová Institute of Mathematics, University of Mainz A. Hundertmark (Uni-Mainz) Š. Nečasová (Academy of Sciences, Prague) G. Rusnáková
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationA monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
1 / 41 A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion Chen-Yu CHIANG O. Pironneau, T.W.H. Sheu, M. Thiriet Laboratoire Jacques-Louis Lions
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
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 multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction
A multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction J Degroote, S Annerel and J Vierendeels Department of Flow, Heat and Combustion Mechanics, Ghent University,
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationarxiv: v1 [math.na] 3 May 2018
; :1 25 Fluid-structure interaction simulations of venous valves: a monolithic ALE method for large structural displacements arxiv:1805.01321v1 [math.na] 3 May 2018 S. Calandrini 1 *, E. Aulisa 1 1 Department
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
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 informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
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 informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationStability and Geometric Conservation Laws
and and Geometric Conservation Laws Dipartimento di Matematica, Università di Pavia http://www-dimat.unipv.it/boffi Advanced Computational Methods for FSI May 3-7, 2006. Ibiza, Spain and Introduction to
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationCoupling techniques for partitioned fluid structure interaction simulations with black box solvers
Coupling techniques for partitioned fluid structure interaction simulations with black box solvers Joris Degroote 1, Robby Haelterman 2, Sebastiaan Annerel 1, Jan Vierendeels 1 1 Ghent University, Department
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
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 informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationMultivariable calculus: What is it good for?
: What is it good for? Department of Mathematics and Statistics, UMBC http://www.math.umbc.edu/ rouben/ January 2011 Purpose The purpose of these slides is to show that the material that one learns in
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
More informationFluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers
Fluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers R. Sanchez 1, D. Thomas 2, R. Palacios 1, V. Terrapon 2 1 Department of Aeronautics, Imperial College London 2 Department
More informationNumerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell
J.S.Howell Research Statement 1 Numerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell Fig. 1. Plot of the magnitude of the velocity (upper left) and
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationBarotropic geophysical flows and two-dimensional fluid flows: Conserved Quantities
Barotropic geophysical flows and two-dimensional fluid flows: Conserved Quantities Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 2016 Advanced Topics in Applied Math Di Qi,
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
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 informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 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 informationA Fully Implicit Domain Decomposition Based ALE Framework for Three-dimensional Fluid-structure Interaction with Application in Blood Flow Computation
A Fully Implicit Domain Decomposition Based ALE Framework for Three-dimensional Fluid-structure Interaction with Application in Blood Flow Computation Yuqi Wu a, Xiao-Chuan Cai b a Department of Applied
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More information