A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
|
|
- Adele Bond
- 6 years ago
- Views:
Transcription
1 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 (LJLL), Sorbonne U. Scientific Computing and Cardiovascular Simulation (SCCS) Lab, NTU FreeFem++ days 2017
2 2 / 41 Outline
3 Formulation of coupled system A general fluid-structure coupled system General Laws of Continuum Mechanics ( ) u ρ t + u ( u) = (σ) + f (1) For Newtonian incompressible fluid : σ f = p f I + µd For hyperelastic incompressible material : σ s = p s I + F ΨF T For hyperelastic compressible material : σ s = J 1 F Ψ F T D = ( u + u ) T X: Ω 0 (0, T ) Ω t : X ( x 0, t ) : Lagrangian position at t of x 0 d = X ( x 0, t ) x 0 : displacement F i,j = x 0 j X i ; x 0 j : original position J = det F, the Jacobian of the deformation 3 / 41
4 4 / 41 Hyperelastic material model Mooney Rivlin model Ψ(F) = c 1 tr F T F + c 2 ( ) tr (F T F) tr 2 2 F T F F ΨF T = (2c 1 4c 2 tr B ) B + 4c 2 B 2, where B = FF T. (2) S t Venant Kirchhoff model Ψ (F) = λ s 2 tr2 E + µ s tr E 2, E = 1 ( F T F I ) 2 [ ( F ΨF T trb = λ s 2 3 ) ] µ s B + µ s B 2 (3) 2
5 5 / 41 Lagrangian points ( x 0, t ) to Eulerian points (x, t) F T ( = x 0X = x 0 d + x 0 ) = x 0d + I = F T d + I F = (I d) T C = I B = Dd d T d Cayley Hamilton theorem For a general n n invertible matrix A, i.e., one with nonzero determinant, A 1 can thus be written as a n th order polynomial expression in A n = 3, A 3 tr A A (tr2 A tr A 2)A det A I = 0 n=3 B = tr B I 1 ( ) tr 2 2 B tr B 2 B 1 + J 2 B 2, B 2 = 1 [ ( ) tr 2 2 B + tr B 2 I + J 2 1 ( ) ] (4) tr 3 2 B tr B tr B 2 B 1 + tr B J 2 B 2
6 Stress tensor for hyperelastic models (n=3) Mooney Rivlin model F ΨF T = 2c 1 ( Dd d T d ) 2 + 2c3 ( Dd d T d ) + αi (5) c 3 := c 1 2 S t Venant Kirchhoff model ( tr 2 B tr B 2 4 ) 2c 2 F ΨF T = αi + β ( Dd d T d ) + γ ( Dd d T d ) 2 (6) { α = λ s 4 (tr B 3) ( 2tr B tr 2 B + tr B 2 + 2J 2) + µs 2 tr ( B 2trB + tr B 2 tr 2 B 2 + 2J 2) β = λs 4 (tr B 3) ( tr 2 B tr B 2 4J 2) ( + µs 2 tr 3 B tr B tr B 2 tr 2 B + tr B 2 + 2J 2 4tr B J 2) γ = λs 2 (tr B 3) J 2 + µ s (tr B 1) J 2 6 / 41
7 7 / 41 Variational formulation (n=3) (ρ f D tu û + µ ) 2 Du: Dû p û ˆp u Ω t f + + β 2 For simplicity, Ω t s [ ρ sd tu û + γ ( 2 Dd d d) T : Dû 2 ( ) Dd d T d ] : Dû + α û = (7) f û Ω t D td = u, (8) homogeneous boundary conditions on Γ Ω homogeneous Neumann conditions on Ω t \ Γ. So, given Ω 0 f, Ω0 s, d, and u at t = 0, we must find ( ) u, p, d, Ω t f, Ωt s with u Γ = 0 and equation (7) and (8)
8 8 / 41 A first order in time consistent approximation Find Ω n+1, u n+1 H 1 0(Ω n+1 ), p n+1 L 2 0(Ω n+1 ) such that for all û n+1 H 1 0(Ω n+1 ), ˆp n+1 L 2 0(Ω n+1 ), n+1 u u n Y n+1 ρn+1 Ωn+1 δt + δt Ω n+1 s û + Ω n+1 f ( p n+1 û ˆp u µ ) 2 Dun+1 : Dû [γ n+1 ( Du n+1 u n+1 T dn dn T u n+1 ) ( D d n d n T dn ) + βn+1 ( Du n+1 u n+1 T T n+1 )] dn dn u : Dû 2 [ γ n+1 + (D d n Ω n+1 d n T ) 2 β n+1 ( dn : Dû + D d n d n T ) dn : Dû s 2 2 +α n+1 ] û = f û Ω n+1 (9) where d n stands for d n (Y n+1 ) and where d n is updated by d n+1 = d n Y n+1 + δtu n+1 (10)
9 Solution algorithm Fixed point iteration 1 Set ρ = ρ n, α = α n, β = β n, γ = γ n, Ω = Ω n, u = u n h, Y (x) = x uδt, 2 Solve equation (9), 3 Set u = u n+1 h, Y (x) = x uδt, Ω r = Y 1 (Ω n r ), r = s, f, update α, β, γ and ρ. 4 If not converged return to step 2. 9 / 41
10 10 / 41 Outline A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin 1 2 A thin elastic plate clamped into a small rigid square body immersed in a flowing fluid A thin elastic plate clamped into a rigid cylinder immersed in a flowing fluid 3
11 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin A thin elastic plate clamped into a small rigid square body immersed in a flowing fluid Table: Material parameters [ ρ f g cm 3] [ ν f g cm 1 s 1] ρ s [g cm 3] E [g cm 1 s 2] U in [cm s 1] u=uin v=0 v=0 u y= u x= v x=0 12 v=0 u y= / 41
12 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin A thin elastic plate clamped into a small rigid square body immersed in a flowing fluid Author present, hmin = 0.02 present, hmin = Hubner (a) magnitude frequency (b) Figure: Contours of velocity magnitude at (a) t = 12.5 and (b) t = / 41
13 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin A thin elastic plate clamped into a rigid cylinder immersed in a flowing fluid Poiseuille profile inflow U (y) = 1.5U in y(h y) (0.5H) 2 and H = 0.41 For fluid, ρ f = 1000 and ν f = Two steady cases, g s = 2, 4 and U in = 0 (structure ρ s = 1000 and E = ) FSI1, U in = 1, ρ s = 1000 and E = FSI2 U in = 2, ρ s = and E = u=u(y) v=0 0.2 u=v= u=v=0 2.3 σ n= / 41
14 14 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin g s=2 g s=4 d x d y d x d y Present h min = Dunne [?] Hron and Turek [?] Richter [?] Wick [?] (a) (b) Figure: (a) displacement on x and y with respect to time, (b) centerlines position.
15 15 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin
16 16 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin FSI2 FSI3 amplitude frequency amplitude frequency Present h min = s s 1 Present h min = s s 1 Turek [?] Dune [?] s s 1 Thomas [?] s 1 - -
17 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin U in = 1.07 Whole solid system flaps due to instability of wake region. Front cylinder rotates at frequency of 6.38 Hz and 7.35 Hz, experimentally and numerically determined value. u=uin v= u=v= σ n= Origin of coordinate is placed at center of cylinder. Measure points A and B locate at (0.082, 0) and (0.082, 0.04) u=v= / 41
18 18 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin
19 19 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin
20 20 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin Relatively good agreement in displacement, frequency and velocity at measure points elements (P1b,P1), and first order in time.
21 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin [ ] [ ] An elastic plate h 2, h b 2 2, b [0, L] 2 is clamped to a rectangular tube [ 5b, 16b] [ 8b, 8b] [0, 17b] h/b = 0.2 and L/b = 5 A constant inflow U in = 1 along x direction at Γ in Re = U inb ν = 1600 and ρ f = 1 ρ s = 0.678, E = , and ν s = 0.4 z x Buoyancy force of solid f b = b Γ in Γ out h L 5b 16b 21 / 41
22 22 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin (a) (b) Figure:. (a) Experimental results [?], (b) computational results.
23 23 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin C D D x /b D z /b Present Luhar and Nepf [?] 1.15 (with 10% error) Tian and Dai etc [?] Table: Comparison of drag coefficient C D and deflection in x and z direction with respect to referenced data.
24 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin A cylindrical chamber with r 2 = 76.2 and length of 200 Two inlet pipes with r 1 = 21.0 and length of 20 A silicon filament with h = 2 thick, b = 11 wide, and l = 65 long ρ f = and ν f = ρ s = , Young s modulus E s = , and Poisson ratioν s = Gravity g = 9810 along y direction r 1 r 1 r 2 b (0,0,0) h l h y r 2 r 1 y r 1 z x 24 / 41
25 Initial configuration due to buoyancy force A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin Figure: Initial deflection of silicone filament because of buoyancy force. 25 / 41
26 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin Case I : steady deformation Inlet velocity are Poiseuille profile with peak velocityof 630 and 615 for upper and lower inlet respectively. (a) (b) Figure: Computational results for phase I (a) position of center line along z direction, (b) velocity profiles at symmetric plane x = / 41
27 27 / 41 Case II : starting flow A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin n 1 n 2 n 3 d 0 b 1 b 2 b 3 b 4 Î k ˆv x [0, 4.07] ˆv y [0, 5.51] ˆv z [0, 5.27] Table: For case II, curve fitting coefficients of inlet peak velocity for ˆv k (t) = Σ 3 i=1 n it i /Σ 4 j=0 b j t j with ˆv k = 0 for t I\I k. Note, that flow in y direction is applied only in the upper inlet.
28 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin (a) (b) (c) (d) (e) (f) Figure: Deflection of the silicone filament in phase II (a) t = 0.073, (b) t = 0.721, (c) t = 1.153, (d) t = 1.585, (e) t = 2.017, (f) t = / 41
29 29 / 41 A thin elastic plate clamped into a small rigid square body immersed A thin elastic plate clamped into a rigid cylinder immersed in a flowin (a) (b) (c) (d) Figure: At midplane x = 0 mm, flow velocity component u z at (a) z 0.20l, (b) z 0.51l, (c) z 0.82l, (d) z 1.12l ( l = 65 mm, length of the silicon filament ).
30 30 / 41 Variational inequality/constrain for contact Valved veins Outline Variational inequality/constrain for contact Valved veins 4
31 31 / 41 Variational inequality/constrain for contact Valved veins Variational inequality/constrain for contact A (u, v) + Λ (u, v) = b (v) where Λ is a variational constrain for contact { λ (x, t) 0, x, Λ (u, v) = λn v, Ω t d S t (x) λ (x, t) = 0, S t i i To each simply onnected disjointed part Si t, i = 1 ns is associated a signed distance function x d S t measuring distance of x to S t i i. d S t i d S t i d S t i < 0, structure region = 0, contact point/line/surface > 0, fluid region for some x Ω/S t i
32 32 / 41 Variational inequality/constrain for contact Valved veins Semi Smooth Newton Method In k th iteration for each time step n + 1 ( λ n+1,k+1 min {0, λ n+1,k + c 0 d S x + δtu n+1,k)}. x S ti n+1 i Original equations can be recast into G ( u n+1,k+1, v ) = ns A ( u n+1,k+1, v ) + τ ck G ( u n+1,k+1, v ) = b (v) i=1 S i n+1 {x : λ n+1,k +c 0 d S i n+1 (x+δtu n+1,k ) } g ( u n+1,k+1, v ) where function f is associated to rebound velocity and variant with different model.
33 Variational inequality/constrain for contact Valved veins Rebound velocity For example g ( u n+1,k+1) = u n+1,k+1 f ( u n+1,k n ) n f(u n )=-βu n, 0<β f(u n )=ε closer away 33 / 41
34 34 / 41 Variational inequality/constrain for contact Valved veins Free falling ball hit a slope
35 35 / 41 Variational inequality/constrain for contact Valved veins d h w=d/20 h T=0.92d w=0.625d h B=0.58d L v 0.8d h v=0.01d (a) (b) Figure: Shape and size of one valved vein
36 36 / 41 Variational inequality/constrain for contact Valved veins
37 37 / 41 Variational inequality/constrain for contact Valved veins
38 38 / 41 Outline
39 Modeling In/out-let boundary conditions of blood vessels Anisotropic hyperelastic material model for vascular walls and valves Model for agregation and deformation of red blood capsules Simulation Parallize whole computational process with domain decomposition method Vascular model and simulation in 3D Comparison with experimental data 39 / 41
40 40 / 41 Thanks for your listening.
41 41 / 41 Flow rate and volume (a) (b)
A 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 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 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 informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
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 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 informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
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 informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
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 information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
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 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 informationSebastian Gjertsen Master s Thesis, Spring 2017
Development of a Verified and Validated Computational Framework for Fluid-Structure Interaction: Investigating Lifting Operators and Numerical Stability Sebastian Gjertsen Master s Thesis, Spring 2017
More informationNumerical Benchmarking of Fluid-Structure Interaction: A comparison of different discretization and solution approaches
Numerical Benchmarking of Fluid-Structure Interaction: A comparison of different discretization and solution approaches Stefan Turek, Jaroslav Hron, Mudassar Razzaq, Hilmar Wobker, and Michael Schäfer
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 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 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 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 informationPart II Fundamentals of Fluid Mechanics By Munson, Young, and Okiishi
Part II Fundamentals of Fluid Mechanics By Munson, Young, and Okiishi WHAT we will learn I. Characterization of Fluids - What is the fluid? (Physical properties of Fluid) II. Behavior of fluids - Fluid
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
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 informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
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 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 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 informationR09. d water surface. Prove that the depth of pressure is equal to p +.
Code No:A109210105 R09 SET-1 B.Tech II Year - I Semester Examinations, December 2011 FLUID MECHANICS (CIVIL ENGINEERING) Time: 3 hours Max. Marks: 75 Answer any five questions All questions carry equal
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationPREDICTION OF PULSATILE 3D FLOW IN ELASTIC TUBES USING STAR CCM+ CODE
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More information1. The Properties of Fluids
1. The Properties of Fluids [This material relates predominantly to modules ELP034, ELP035] 1.1 Fluids 1.1 Fluids 1.2 Newton s Law of Viscosity 1.3 Fluids Vs Solids 1.4 Liquids Vs Gases 1.5 Causes of viscosity
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 informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
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: v1 [math.na] 19 Dec 2013
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH HYPERELASTIC MODELS I: A PARTITIONED APPROACH arxiv:1312.5561v1 [math.na] 19 Dec 2013 ULRICH LANGER AND HUIDONG YANG Abstract. In this
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 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 informationAn Added-Mass Partition Algorithm for Fluid-Structure Interactions of Compressible Fluids and Nonlinear Solids
An Added-Mass Partition Algorithm for Fluid-Structure Interactions of Compressible Fluids and Nonlinear Solids J. W. Banks a,,2, W. D. Henshaw a,,3,, A. K. Kapila a,4, D. W. Schwendeman a,,3,4 a Department
More information2007 Problem Topic Comment 1 Kinematics Position-time equation Kinematics 7 2 Kinematics Velocity-time graph Dynamics 6 3 Kinematics Average velocity
2007 Problem Topic Comment 1 Kinematics Position-time equation Kinematics 7 2 Kinematics Velocity-time graph Dynamics 6 3 Kinematics Average velocity Energy 7 4 Kinematics Free fall Collisions 3 5 Dynamics
More informationLaw of behavior very-rubber band: almost incompressible material
Titre : Loi de comportement hyperélastique : matériau pres[...] Date : 25/09/2013 Page : 1/9 Law of behavior very-rubber band: almost incompressible material Summary: One describes here the formulation
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 informationMECHANICAL PROPERTIES OF FLUIDS:
Important Definitions: MECHANICAL PROPERTIES OF FLUIDS: Fluid: A substance that can flow is called Fluid Both liquids and gases are fluids Pressure: The normal force acting per unit area of a surface is
More informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 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 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 informationFluid Mechanics. Jim Emery 9/3/ Bernoulli s Law 2. 2 Torricelli s Law 3. 3 Time to Empty a Tank 3. 4 Viscosity 4. 5 The Acceleration 5
Fluid Mechanics Jim Emery 9/3/2010 Contents 1 Bernoulli s Law 2 2 Torricelli s Law 3 3 Time to Empty a Tank 3 4 Viscosity 4 5 The Acceleration 5 6 The Equation of Motion Due to Cauchy 7 7 The Deformation
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 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 informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationReview of Fluid Mechanics
Chapter 3 Review of Fluid Mechanics 3.1 Units and Basic Definitions Newton s Second law forms the basis of all units of measurement. For a particle of mass m subjected to a resultant force F the law may
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 informationUnified Lagrangian formulation for solid and fluid mechanics and FSI problems
Unified Lagrangian formulation for solid and fluid mechanics and FSI problems Alessandro Franci 1,2, Eugenio Oñate 1,2, Josep Maria Carbonell 1,2 1 International Center for Numerical Methods in Engineering
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 information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
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 informationThe Non-Linear Field Theories of Mechanics
С. Truesdell-W.Noll The Non-Linear Field Theories of Mechanics Second Edition with 28 Figures Springer-Verlag Berlin Heidelberg NewYork London Paris Tokyo Hong Kong Barcelona Budapest Contents. The Non-Linear
More informationNonlinear Modeling of Fiber-Reinforced Elastomers and the Response of a Rubber Muscle Actuator
Nonlinear Modeling of Fiber-Reinforced Elastomers and the Response of a Rubber Muscle Actuator Larry D. Peel, Ph.D.* Department of Mechanical & Industrial Engineering Texas A&M Univ. - Kingsville David
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 informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationarxiv: v1 [physics.flu-dyn] 3 Nov 2017
This draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics arxiv:7.65v [physics.flu-dyn] 3 Nov 27 On the Mechanism of Large Amplitude Flapping of Inverted Foil in a
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 informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationThis introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference.
Chapter 1 Introduction to Elasticity This introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference. 1.1 Kinematics of finite deformations We
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationOptimized Surface Acoustic Waves Devices With FreeFem++ Using an Original FEM/BEM Numerical Model
Optimized Surface Acoustic Waves Devices With FreeFem++ Using an Original FEM/BEM Numerical Model P. Ventura*, F. Hecht**, Pierre Dufilié*** *PV R&D Consulting, Nice, France Laboratoire LEM3, Université
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 informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
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 informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
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 informationMultiplicative Jacobian Energy Decomposition Method for Fast Porous Visco-Hyperelastic Soft Tissue Model
Multiplicative Jacobian Energy Decomposition Method for Fast Porous Visco-Hyperelastic Soft Tissue Model Stéphanie Marchesseau 1, Tobias Heimann 1, Simon Chatelin 2, Rémy Willinger 2, and Hervé Delingette
More informationFluid Mechanics. Spring 2009
Instructor: Dr. Yang-Cheng Shih Department of Energy and Refrigerating Air-Conditioning Engineering National Taipei University of Technology Spring 2009 Chapter 1 Introduction 1-1 General Remarks 1-2 Scope
More informationA MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 A MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD James E. HILTON
More information1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts)
1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts) Introduction, general information Course overview Fluids as a continuum Density Compressibility Viscosity Exercises: A1 Fluid mechanics Fluid
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationStress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation
Energy Research Journal 1 (2): 165-170, 2010 ISSN 1949-0151 2010 Science Publications Stress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation P. Khamdaengyodtai, P. Sakulchangsatjatai
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 informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
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 informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationA numerical study of vortex-induced vibrations (viv) in an elastic cantilever
A numerical study of vortex-induced vibrations (viv) in an elastic cantilever Lorentzon, Johan; Austrell, Per Erik; Szász, Robert-Zoltán; Revstedt, Johan 2009 Link to publication Citation for published
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 informationSTRUCTURAL ANALYSIS OF A WESTFALL 2800 MIXER, BETA = 0.8 GFS R1. By Kimbal A. Hall, PE. Submitted to: WESTFALL MANUFACTURING COMPANY
STRUCTURAL ANALYSIS OF A WESTFALL 2800 MIXER, BETA = 0.8 GFS-411519-1R1 By Kimbal A. Hall, PE Submitted to: WESTFALL MANUFACTURING COMPANY OCTOBER 2011 ALDEN RESEARCH LABORATORY, INC. 30 Shrewsbury Street
More informationExam paper: Biomechanics
Exam paper: Biomechanics Tuesday August 10th 2010; 9.00-12.00 AM Code: 8W020 Biomedical Engineering Department, Eindhoven University of Technology The exam comprises 10 problems. Every problem has a maximum
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 informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics
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 informationSamantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2
Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force
More informationTraction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment
Traction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment Colangeli E., Repetto R., Tatone A. and Testa A. Grenoble, 24 th October 2007 Table of contents
More information