An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
|
|
- Lorena Sanders
- 5 years ago
- Views:
Transcription
1 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 Joint work with: S. Canic and A. Quaini Acknowledgements: NSF DMS , DMS
2 Motivation Outline Motivation for Applications in Hemodynamics
3 Motivation Main motivation We are interested in study the blood flow inside the Abdominal Aortic Aneurysm with stent-graft treated. Provide new information and give insight in new treatment Serve as a potential predictor of AAA arterial wall rupture Why Discontinuous Galerkin? High order and less dissipative solution Ideal for problems contain discontinuous solution
4 Fluid model Thin structure model The coupling conditions Fluid model We treat blood as an incompressible, viscous Newtonian fluid, and its motion is described by the Navier-Stokes equations. tu + (u )u = 1 ρ ( σ) = 1 p + ν u ρ u = 0 in Ω(t), in Ω(t), where u is the fluid velocity, σ is the Cauchy stress tensor, ρ is the fluid density, and ν is the fluid kinematic viscosity. The Cauchy stress tensor is given by σ = pi + 2µD(u), where p is the pressure and D(u) = ( u + ( u) T )/2 is the strain rate tensor.
5 Fluid model Thin structure model The coupling conditions Thin structure model The elastodynamics of the thin structure will be described in general by ρ sh s ttη + Lη + γv tη = f s in Γ 0, Operator L is a linear operator acting on vector function η describing the elastic properties of the structure, and V is a linear operator acting on vector function tη describing the viscoelastic properties of the structure. Γ 0 denote the fixed reference structure domain, ρ s is the structure density, h s is the structure thickness, γ is a viscoelastic parameter, and f s is the force acting on the structure. We consider two different forms of operators L and V.
6 Fluid model Thin structure model The coupling conditions Thin structure model - I We consider the generalized string model for an elastic tube with the reference configuration as a cylinder of radius R 0 and length L (Quarteroni 2000). And assuming that the longitudinal displacement η x and the azimuthal displacement η z are negligible. ( 2 η y Lη = 0, kgh s x 2 + Ehs (1 νs 2 ) η y R 2 0 ) T, 0, V η ( ) T t = 0, 3 η y x 2 t, 0, where k is the Timoshenko shear correction factor, G is the shear modulus, E is the Young s modulus and ν s is the Poisson ratio of the structure. Structure model is endowed with absorbing boundary conditions of first order: η y t η y t + kg η y x ρ s kg η y ρ s x = 0 at x = 0, = 0 at x = L.
7 Fluid model Thin structure model The coupling conditions Thin structure model - II The second model we consider involves no spatial derivatives and no viscoelasticity. Operators L and V are given by: Lη = (C x η x, C y η y, 0) T, V = (0, 0, 0) T, where C x and C y are spring constants. Notice that the displacement components are independent. The structure model is supplemented with homogeneous Dirichlet boundary conditions: η = 0 on Γ 0.
8 Fluid model Thin structure model The coupling conditions The coupling conditions The fluid and the structure are coupled via two boundary conditions: - Kinematic coupling condition: it describes the continuity of velocity at the fluid-structure interface (no-slip condition) + u A t = tη on Γ 0 (0, T ); - Dynamic coupling condition: it describes the continuity of the stress at the fluid-structure interface + J σn Γ(t) = f s on Γ 0 (0, T ), where J denotes the Jacobian of the transformation from Eulerian to Lagrangian coordinates, and σn Γ(t) denotes the normal fluid stress at the deformed fluid-structure interface, evaluated with respect to the reference configuration. Vector n is the outward normal to the deformed fluid domain.
9 Fluid model Thin structure model The coupling conditions The coupling conditions Coupling conditions above can be written in the equivalent form: α f u A t J σn Γ(t) = α f tη f s on Γ 0 (0, T ), α su A t + J σn Γ(t) = α s tη + f s on Γ 0 (0, T ), where α f > 0 and α s > 0 (α f α s) are constants.
10 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Partitioned methods for the FSI problem The FSI problem will be solved using two different partitioned strategies based on Domain Decomposition methods: the Dirichlet-Neumann (DN) and the Robin-Neumann (RN) algorithms. (A. Quarteroni and A. Valli 1999) The Dirichlet-Neumann: Dirichlet boundary conditions are imposed at the interface Γ(t) for the fluid problem, whereas Neumann boundary conditions are supplemented for the structure problem. The DN algorithm iterates over these two sub-problems until meeting the convergence criteria for the interface displacement. When the density of structure has the same magnitude as fluid density, relaxation step on the structure displacement is required (Nobile 2001). The Robin-Neumann: When the FSI problem contains a closure structure, the Dirichlet-Neumann algorithms fails due to the fact that the coupling conditions are not implemented synchronously. Robin-Neumann algorithm always converges without any relaxation and its convergence is almost insensitive to the added-mass effect.
11 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise We consider a smooth mapping: A t : Ω 0 Ω(t), A t : x 0 x, where x and x 0 are the coordinates in the physical domain Ω(t) and the reference domain Ω 0, respectively. The domain velocity w is given by w(t, ) = dat dt (t, At(t, ) 1 ). The ALE time derivative of the fluid velocity is defined as: tu x0 = D tu(t, A t(x 0 )) = tu(t, x) + w(t, x) u(t, x), for x = A t(x 0 ), x 0 Ω 0.
12 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise The incompressible Navier-Stokes equations in the ALE formulation: tu x0 + (u w) u = 1 p + ν u ρ u = 0 in Ω(t), in Ω(t),
13 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise We introduce a notation of F stands for the nonlinear flux vector. Given the velocity components u = (u, v) T for d = 2 and u = (u, v, w) T for d = 3, F is defined as: F = [ ] u 2 uv vu v 2 if d = 2, or u 2 uv uw vu v 2 vw if d = 3. wu wv w 2 Let J At = det( x x 0 ) be the Jacobian of the deformation gradient: tj At x0 = J At w. We can derive the fully conservative form of the NS equation in ALE: t(j At u) + J At (F wu) = J A t x0 ρ p + J A t ν u.
14 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Harmonic extension Describing the ALE map in terms of the fluid domain displacement d(x 0, t) = A t(x 0 ) x 0. Let η be the displacement of the deformable boundary Γ(t). Inside Ω 0 the displacement d is arbitrary: it can be any reasonable extension of η over Ω 0. A classical choice is to consider a harmonic extension in the reference domain, that is: with boundary conditions: d = 0, in Ω(t), d = η d = 0 on Γ(t), on Ω(t)/Γ(t). The fluid domain velocity w is obtained by differentiation: w = td.
15 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Temporal discretization of the fluid problem The second order semi-implicit splitting scheme is used to discretize NS equation (Karniadakis 1991). The NS equation in the reference domain Ω 0 is splitted to three sub-equations. Step 1: Using the Backward Difference Formula (BDF2) to obtain the first intermediate velocity field J n u k+ 1 by the explicit approximation of the convective 3 component (F wu): (3J n u k+ 1 4J n u n + J n 1 u n 1 ) 3 = 2J n (F wu) n + J n 1 (F wu) n 1. 2 t Step 2: Known the J n u k+ 1 3, we solve for the second intermediate velocity field J n u k+ 2 3 by including the pressure gradient term: J n u k+ 2 J n u 3 3 k t = 1 ρ Jn p k+1. We impose the divergence free requirement to J n u k+ 2 3, (J n International Workshop on Fluid-Structure Interaction Problemsu k+ 2 IMS ) = National 0. University of Singapore 30 May - 3 June 2016
16 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Temporal discretization of the fluid problem Taking the divergence of both sides of the equation: 1 ρ Jn p k+1 = 3 2 t (Jn u k+ 1 ), 3 and solve for the pressure J n p k+1. After obtaining the pressure field J n p k+1, the second intermediate velocity field J n u k+ 2 is updated from: 3 J n u k+ 2 3 = J n u k t 3ρ Jn p k+1. Step 3: Known the second intermediate velocity field J n u k+ 2 3, we implicitly integrate the viscous component ( u) to solve a Helmholtz problem for J n u k+1 : J n u k+1 J n u k t = νj n u k+1.
17 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Discontinuous Galerkin method Let us consider the local polynomial approximation of fluid velocity component of u(r, s): N N u(r, s) = û i Φ i (r, s) = u j l j (r, s), i=1 j=1 Here, we are assuming that the number of nodes and number of modes are the same, and are equal to N. û i is the expansion coefficient associating with the orthonormal modal base Φ i (r, s), while u j is the value of u(r, s) at grid point of (r, s) j and l j (r, s) is the corresponding nodal base. In the 2D case, we have N = (n + 1)(n + 2)/2 which represents the least number of terms required to have the completeness of expansion using the n th order polynomial.
18 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Discontinuous Galerkin method Modal bases: A set of orthonormal polynomial modes defined in standard isosceles right triangle Ω st (Gottlieb 1977, Dubiner 1991,Hesthaven 2008): Φ i (r, s) = 2P 0,0 k (2 1 + r 1 s 1)P(2k+1,0) l (s)(1 s) k, in which r and s are the coordinates of two short sides of an isosceles right triangle. Subscript i stands for the i th mode of the basis functions, and i = (N + 1)k + l + 1 k (k 1), (k, l) 0; k + l N. P(α,β) 2 k (x) is Jacobi polynomials. Nodal bases: Lagrangian polynomials built on Fekete points (GKarniadakis 2005)
19 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Discontinuous Galerkin method Nodal bases are expressed in terms of modal bases to simplify the evaluation of the numerical integration. And the generalized Vandermonde matrix V works as a role to connect the nodal bases with modal bases, and it is defined as: V ij = Φ j (r, s) i. Between modal Φ j (r, s) and nodal l j (r, s) bases, one can be transformed from the other via Vandermonde matrix: l(r, s) = (V T ) 1 Φ(r, s). In other words, nodal basis l i (r, s), it can be expressed as: N l i (r, s) = (V T ) 1 ik Φ k(r, s). k=1
20 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Discontinuous Galerkin method In an standard isosceles triangle, the local mass matrix is: M ij = (VV T ) 1 ij. And the local matrices for convective term are: where D r and D s can be computed from: S rij = (MD r ) ij, S sij = (MD s) ij. D r = V r V 1, D s = V sv 1. And V r and V s denote two Vandermonde matrices as below: V rij = dφ j (r, s) dr (r,s)i, V sij = dφ j (r, s) ds (r,s)i.
21 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Spatial discretization of the fluid problem Step 1: For the convective equation, we have: 3J n u k+ 1 4J n u n + J n 1 u n 1 M 3 e 2 t = 2J n S : (F wu) n + J n 1 S : (F wu) n 1 2J n n ((F wu) (F wu) ) n l(x, y)dxdy E + J n 1 n ((F wu) (F wu) ) n 1 l(x, y)dxdy. E Local Lax-Friedrichs flux is used and denoted by (F wu) : (F wu) = (F wu)+ + (F wu) + τ(u + u ), 2 in which the coefficient τ is the local maximum velocity magnitude.
22 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Spatial discretization of the fluid problem Step 2: For the pressure possion equation, we convert it into a system of first order equations (Cockburn and Shu 1998): J n ρ p k+1 = q q = 3 2 t (Jn u k+ 1 3 ) In each triangle E, we discretize (J n p k+1, q) as ((J n p k+1 ) e, q x, q y ) and approximate it by the n th order basis function l(x, y) constructed on that triangle. Internal penalty fluxes as below are used: q = Jn ( (pk+1 ) + e + (p k+1) e ρ 2 p k+1 = p+ k+1 + p k+1 2 where the τ is penalty parameter. τ [ n (p k+1 ) e. + n+ (p k+1 ) + ] ) e,
23 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Spatial discretization of the fluid problem Then we have: J n ( (S x ) ep k+1 ρ E J n ( (S y ) ep k+1 ρ E (S x ) eq x n (q x q )l(x, y)dxdy E + (S y ) eq y E ( n x pk+1 p ) ) k+1 l(x, y)dxdy = M eq x, ( n y pk+1 p ) ) k+1 l(x, y)dxdy = M eq y. n (q y q )l(x, y)dxdy = 3 2 t S (Jn u k+ 1 3 ).
24 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Temporal discretization of the structure problem Consider the general string model for the structure, we introduce a new variable of U: U = ρ sh s η y t γ 2 η y x 2, which combines the temporal derivative terms. Such that the problem can be rewritten as: ( ρ sh s η y t t γ 2 η y x 2 ) = U t = 2 η y kghs x 2 Ehs (1 ν 2 ) η y R J σn Γ(t). We also introduce new variables W denote the second order spatial derivative term 2 η y x 2.
25 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Temporal discretization of the structure problem As a result, the original problem is transformed into a system of ODEs: η y ρ sh s = U + γw, t U t = kghsw Ehs η y (1 ν 2 ) R0 2 + J σn Γ(t). The initial Conditions of η and U can be reasonably obtained as following: η y t=0 = 0, η y t t=0 = 0, U t=0 = ρ sh s η y t t=0 γw t=0 = 0.
26 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Temporal discretization of the structure problem The above system can be discretized by using forward difference scheme. Given η n y, U n, W n and σ k+1, find the value of (η y ) k+1, U k+1 and W k+1, so that: (η y ) k+1 ηy n ρ sh s = U k+1 + γw n, t U k+1 U n = kgh sw n Ehs ηy n t 1 ν 2 R0 2 + J n σ k+1 n Γ(t n ), W k+1 = 2 (η y ) k+1 x 2.
27 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Spatial discretization of the structure problem We introduce an intermediate variable q to split W into a system of two first order differential equations which can later be easily solved by Local Discontinuous Galerkin (LDG) (Cockburn and Shu 1998): q = ηy x W = q x By choosing the fluxes as q = q and η y = η + y, we can write out the matrix form of eq. (1) as below: M eq = (S x ) eη y n x (η y η y )l(x, y)dxdy, E M ew = (S x ) eq n x (q q )l(x, y)dxdy; E.
28 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Element-wise The internal penalty flux also is used in solving the displacement equation in the mapping, our approach does not force two neighbor triangles to have the same displacement at the common edge but limit the differences between them and disallow large gaps in the displacement. We reconstruct harmonic extension equation as a system of two first order equations for each component of displacement vector: { d = q q = 0 where d corresponds to each component of the displacement vector. We then replace q with the following internal penalty fluxes: q = d+ e + d e 2 τ[n d + n + d + ], d = d+ e + d e 2 Here, d correspond to the vector of each component of displacement. Parameter τ is the penalty parameter..
29 Partitioned methods for the FSI problem Harmonic extension Discretization of the fluid problem Discretization of the structure problem Element-wise Element-wise Then the martix form are: M eq x = (S x ) ed n x (d d )l(x, y)dxdy, E M eq y = (S y ) ed n y (d d )l(x, y)dxdy; E (S x ) eq x + (S y ) eq y n (q x q )l(x, y)dxdy E n (q y q )l(x, y)dxdy = 0. E
30 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 1: Compliant vessel Inlet Γ Γ Outlet (a) Computational domain (b) Inlet data over time (c) Mesh Figure: Test 1: (a) computational domain at the initial time, (b) the inlet pressure pulse, and (c) unstructured triangle mesh (fine mesh).
31 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 1: Compliant vessel (a) t = 2 ms (b) t = 4 ms (c) t = 6 ms (d) t = 8 ms (e) t = 10 ms (f) t = 12 ms Figure: Test 1: pressure at different times: (a) t = 2 ms, (b) t = 4 ms, (c) t = 6 ms, (d) t = 8 ms, (e) t = 10 ms, (f) t = 12 ms. The unit for the pressure is dynes/cm 2.
32 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 1: Compliant vessel (a) t = 2 ms (b) t = 4 ms (c) t = 6 ms (d) t = 8 ms (e) t = 10 ms (f) t = 12 ms Figure: Test 1: pressure along the centerline at different times: (a) t = 2 ms, (b) t = 4 ms, (c) t = 6 ms, (d) t = 8 ms, (e) t = 10 ms, (f) t = 12 ms. In the legend, A. Quaini denotes the results in (Quaini 2009), DG-IP HI and DG-IP LO denote the results obtained with our DG solver with IP fluxes with the fine and coarse mesh, respectively.
33 Introduction Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 2: Oscillating immersed linearly elastic 1D membrane (a) coarse mesh (b) fine mesh Figure: Test 3: (a) coarse mesh and (b) fine mesh. International Workshop on Fluid-Structure Interaction Problems IMS National University of Singapore 30 May - 3 June 2016
34 Introduction Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 2: Oscillating immersed linearly elastic 1D membrane (a) t = 0.1 s (b) t = 0.5 s (c) t = 1 s (e) t = 2 s (f) t = 8 s (d) t = 1.5 s Figure: Test 3: pressure contour and velocity vector field at time (a) t = 0.1 s, (b) t = 0.5 s, (c) t = 1 s, (d) t = 1.5 s, (e) t = 2 s, (f) t = 8 s. International Workshop on Fluid-Structure Interaction Problems IMS National University of Singapore 30 May - 3 June 2016
35 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Benchmark test 2: Oscillating immersed linearly elastic 1D membrane (a) maximum x-coordinate over time (b) difference Figure: Test 3: (a) maximum x-coordinate over time, (b) difference in absolute value between the curves in (a).
36 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Section of an Abdominal Aortic Aneurysm We consider a section of a patient-specific descending aorta affected by AAA with compliant wall. Inlet Triple point 4.03 cm Γ stent 11.2 cm Γ Outlet Outlet (a) Geometry (b) Mesh Figure: (a) Geometry of the section of aorta affected by AAA under consideration and (b) mesh used for the computations.
37 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation (a) Inflow Velocity (b) Outlet Pressure Figure: Rescaled (a) experimental average velocity (cm 2 /s) at the inlet section and (b) experimental pressure (Pa) at the outlet section.
38 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation (a) Pressure (b) Velocity magnitude (c) Displacement magnitude
39 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation (d) Pressure (e) Velocity magnitude (f) Displacement magnitude Figure: (a) and (d) Pressure, (b) and (e) velocity magnitude and vector field, (c) and (f) displacement magnitude on a section of aorta affected by AAA without (top) and with (bottom) stent-graft at time t = 2.3 s.
40 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation (a) Pressure (b) Velocity (c) Displacement
41 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation (d) Pressure (e) Velocity (f) Displacement Figure: (a) and (d) Pressure, (b) and (e) velocity magnitude and vector field, (c) and (f) displacement magnitude on a section of aorta affected by AAA without (top) and with (bottom) stent-graft at time t = 2.85 s.
42 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Figure: Maximum displacement magnitude over time (3 cycles): blue line corresponds the case without stent-graft, red line corresponds the case with the stent-graft.
43 Benchmark test 1 Benchmark test 2 Section of an Abdominal Aortic Aneurysm Preliminary 3D simulation Preliminary 3D simulation (a) Inlet data over time Figure: The inlet pressure pulse. Pressure Propagation Vortex
44 We proposed an DG ALE approach for FSI problem. Benefit from allowing an element-wise mapping, the implementation of DG in ALE frame becomes straight forward. For certain FSI which contains discontinuity solution our DG-ALE approach will be a feasible approach. Our solver is a serial version only at the moment, parallelization is desired.
45 THANK YOU FOR YOUR ATTENTION!
Stability 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 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 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 informationFluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement
Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement Martina Bukač a, Sunčica Čanića, Roland Glowinski a,b, Josip Tambača c, Annalisa Quaini a a Department of
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 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 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 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 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 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 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 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 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 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 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 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 informationA NONLINEAR FLUID-STRUCTURE INTERACTION PROBLEM IN COMPLIANT ARTERIES TREATED WITH VASCULAR STENTS
A NONLINEAR FLUID-STRUCTURE INTERACTION PROBLEM IN COMPLIANT ARTERIES TREATED WITH VASCULAR STENTS MARTINA BUKAČ, SUNČICA ČANIĆ, AND BORIS MUHA Abstract. We study a nonlinear fluid-structure interaction
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 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 informationPULSE WAVE PROPAGATION IN LARGE BLOOD VESSELS BASED ON FLUID-SOLID INTERACTIONS METHODS
PULSE WAVE PROPAGATION IN LARGE BLOOD VESSELS BASED ON FLUID-SOLID INTERACTIONS METHODS Tomohiro Fukui 1,, Kim H. Parker 2 and Takami Yamaguchi 3 1. Department of Mechanical and System Engineering, Kyoto
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationBME 419/519 Hernandez 2002
Vascular Biology 2 - Hemodynamics A. Flow relationships : some basic definitions Q v = A v = velocity, Q = flow rate A = cross sectional area Ohm s Law for fluids: Flow is driven by a pressure gradient
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 informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
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 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 multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More informationMathematical Models and Numerical Simulations for the Blood Flow in Large Vessels
Mathematical Models and Numerical Simulations for the Blood Flow in Large Vessels Balazs ALBERT 1 Titus PETRILA 2a Corresponding author 1 Babes-Bolyai University M. Kogalniceanu nr. 1 400084 Cluj-Napoca
More informationLarge Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes
Large Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes Thomas Gallinger * Alexander Kupzok Roland Wüchner Kai-Uwe Bletzinger Lehrstuhl für Statik Technische Universität München
More informationA monolithic approach to fluid-composite structure interaction
Noname manuscript No. (will be inserted by the editor) A monolithic approach to fluid-composite structure interaction Davide Forti Martina Buac Annalisa Quaini Suncica Canic Simone Deparis Received: date
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 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 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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 informationFlow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems
Flow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems Martin Engel and Michael Griebel Institute of Numerical Simulation, University of Bonn, Wegelerstr.
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 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 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 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 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 informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
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 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 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 information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
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 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 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 informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationAn Operator Splitting Approach to the Solution of Fluid-Structure Interaction Problems in Hemodynamics
An Operator Splitting Approach to the Solution of Fluid-Structure Interaction Problems in Hemodynamics Martina Bukač, Sunčica Čanić, Boris Muha, and Roland Glowinski Abstract We present a loosely coupled
More informationNavier-Stokes Flow in Cylindrical Elastic Tubes
Navier-Stokes Flow in Cylindrical Elastic Tubes Taha Sochi University College London, Department of Physics & stronomy, Gower Street, London, WC1E 6BT Email: t.sochi@ucl.ac.uk. bstract nalytical expressions
More informationDetermination of pressure data in aortic valves
Determination of pressure data in aortic valves Helena Švihlová a, Jaroslav Hron a, Josef Málek a, K.R.Rajagopal b and Keshava Rajagopal c a, Mathematical Institute, Charles University, Czech Republic
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationNumerical Study of Blood Flow through Symmetry and Non- Symmetric Stenosis Artery under Various Flow Rates
IOSR Journal of Dental and Medical Sciences (IOSR-JDMS) e-issn: 2279-0853, p-issn: 2279-0861.Volume 16, Issue 6 Ver. I (June. 2017), PP 106-115 www.iosrjournals.org Numerical Study of Blood Flow through
More informationNumerical Simulations of the Fluid-Structure Coupling in Physiological Vessels Mini-Workshop I
Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels INRIA Rocquencourt, projet REO miguel.fernandez@inria.fr CRM Spring School, Montréal May 16, 2005 Motivations It is generally
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
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 informationChalmers University of Technology
2016-11-16 1 What is an Abdominal Aortic Aneurysm(AAA)and how is it rectified? AAA's are caused due to local dilation of the abdominal aortic vessel common among males 60 years or over How do we treat
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationAccepted Manuscript. A study of self-propelled elastic cylindrical micro-swimmers using modeling and computation
Accepted Manuscript A study of self-propelled elastic cylindrical micro-swimmers using modeling and computation Lingling Shi, Sunčica Čanić, Annalisa Quaini, Tsorng-Whay Pan PII: S1-9991(16)143-1 DOI:
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 informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 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 informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
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 informationNumerical study of blood fluid rheology in the abdominal aorta
Design and Nature IV 169 Numerical study of blood fluid rheology in the abdominal aorta F. Carneiro 1, V. Gama Ribeiro 2, J. C. F. Teixeira 1 & S. F. C. F. Teixeira 3 1 Universidade do Minho, Departamento
More informationModeling of non-newtonian Blood Flow through a Stenosed Artery Incorporating Fluid-Structure Interaction
Modeling of non-newtonian Blood Flow through a Stenosed Artery Incorporating Fluid-Structure Interaction W. Y. Chan Y.Ding J. Y. Tu December 8, 2006 Abstract This study investigated fluid and structural
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationAnalysis and Simulation of Blood Flow in the Portal Vein with Uncertainty Quantification
Analysis and Simulation of Blood Flow in the Portal Vein with Uncertainty Quantification João Pedro Carvalho Rêgo de Serra e Moura Instituto Superior Técnico Abstract Blood flow simulations in CFD are
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 informationProposal for numerical benchmarking of fluid-structure interaction between an elastic object and laminar incompressible flow
Proposal for numerical benchmarking of fluid-structure interaction between an elastic object and laminar incompressible flow Stefan Turek and Jaroslav Hron Institute for Applied Mathematics and Numerics,
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 informationA nonlinear, 3D fluid-structure interaction problem driven by the time-dependent dynamic pressure data: a constructive existence proof
umerical Analysis and Scientific Computing Preprint Seria A nonlinear, 3D fluid-structure interaction problem driven by the time-dependent dynamic pressure data: a constructive existence proof B. Muha
More informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION REP- RESENTED BY HUMAN VOCAL FOLD
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION REP- RESENTED BY HUMAN VOCAL FOLD J. Valášek 1, P. Sváček 1 1 Department of Technical Mathematics, Faculty of Mechanical Engineering, Czech Technical
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 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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationTECHNISCHE UNIVERSITEIT EINDHOVEN Department of Biomedical Engineering, section Cardiovascular Biomechanics
TECHNISCHE UNIVERSITEIT EINDHOVEN Department of Biomedical Engineering, section Cardiovascular Biomechanics Exam Cardiovascular Fluid Mechanics (8W9) page 1/4 Monday March 1, 8, 14-17 hour Maximum score
More informationTWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM
Proceedings of ALGORITMY 212 pp. 85 94 TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM JAROSLAVA HASNEDLOVÁ-PROKOPOVÁ,, MILOSLAV FEISTAUER, JAROMÍR HORÁČEK, ADAM KOSÍK, AND VÁCLAV
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More information