Solving Fluid Structure Interaction Problems on Overlapping Grids
|
|
- Ella Allison
- 5 years ago
- Views:
Transcription
1 Solving Fluid Structure Interaction Problems on Overlapping Grids Bill Henshaw Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA, USA. Virginia Tech, Blacksburg Virginia, October 2, Henshaw (LLNL) FSI on overlapping grids VaTech 1 / 46
2 Acknowledgments. Primary collaborators for this effort: Jeff Banks (LLNL), Kyle Chand (LLNL), Don Schwendeman (RPI), Björn Sjögreen (LLNL). Supported by: ASCR Department of Energy, Office of Science, ASCR Applied Math Program. LDRD LLNL: Laboratory Directed Research and Development (LDRD) program. NSF National Science Foundation. Henshaw (LLNL) FSI on overlapping grids VaTech 2 / 46
3 Outline 1 Background: overlapping grids, Overture and CG. 2 Incompressible flows and rigid bodies (split-step schemes). 1 high-order accurate factored scheme and boundary conditions. 2 matrix-free multigrid. 3 moving grid generation. 3 Compressible flow and deforming bodies (partitioned schemes). 1 Deforming Composite Grids (DCG). 2 a one-dimensional model problem. 3 an interface projection scheme (fluid-solid Riemann problem). 4 Compressible flow and light rigid bodies (partitioned schemes). 1 a one-dimensional model problem. 2 the added-mass Newton Euler equations. 3 overcoming the added-mass instability for light solids. Henshaw (LLNL) FSI on overlapping grids VaTech 3 / 46
4 The Overture project is developing PDE solvers for a wide class of continuum mechanics applications. Overture is a toolkit for solving PDE s on overlapping grids and includes CAD, grid generation, numerical approximations, AMR and graphics. The CG (Composite Grid) suite of PDE solvers (cgcns, cgins, cgmx, cgsm, cgad, cgmp) provide algorithms for modeling gases, fluids, solids and E&M. Overture and CG are available from We are looking at a variety of applications: wind turbines, building flows (cgins), explosives modeling (cgcns), fluid-structure interactions (e.g. blast effects) (cgmp+cgcns+cgsm), conjugate heat transfer (e.g. NIF holhraum) (cgmp+cgins+cgad), damage mitigation in NIF laser optics (cgmx). Henshaw (LLNL) FSI on overlapping grids VaTech 4 / 46
5 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
6 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
7 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
8 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
9 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
10 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
11 What are overlapping grids and why are they useful? Overlapping grid: a set of structured grids that overlap. Overlapping grids can be rapidly generated as bodies move. High quality grids under large displacements. Cartesian grids for efficiency. Smooth grids for accuracy at boundaries. Efficient for high-order methods. Asymptotic Performance Principle for overlapping grids As grids are refined, total CPU/memory usage can approach that of a Cartesian grid. Henshaw (LLNL) FSI on overlapping grids VaTech 5 / 46
12 Components of an overlapping grid G 1 G 1 G 2 G 2 interpolation ghost unused Left: an overlapping grid consisting of two structured curvilinear component grids, x = G 1 (r) and x = G 2 (r). Middle and right: component grids for the square and annular grids in the unit square parameter space r. Grid points are classified as discretization points, interpolation points or unused points. Ghost points are used to apply boundary conditions. Henshaw (LLNL) FSI on overlapping grids VaTech 6 / 46
13 Incompressible flows and moving bodies. Challenges: efficient solvers (esp. since grids change at every time step). high-order accuracy including the boundary. efficient grid generation. Henshaw (LLNL) FSI on overlapping grids VaTech 7 / 46
14 Incompressible Navier-Stokes. u t +(u )u + p ν u f = 0, t > 0, x Ω u = 0 t > 0, x Ω Divergence damping term: α u is important. Wall boundary conditions: u = 0, u = 0, (pressure BC) x Ω, with numerical boundary condition: p n = n ( ν u )+n f. Use u instead of u for implicit time-stepping. Henshaw (LLNL) FSI on overlapping grids VaTech 8 / 46
15 Incompressible Navier-Stokes. Split-step, velocity-pressure formulation: u t +(u )u + p ν u f = 0, t > 0, x Ω p u : u α u f = 0, t > 0, x Ω Divergence damping term: α u is important. Wall boundary conditions: u = 0, u = 0, (pressure BC) x Ω, with numerical boundary condition: p n = n ( ν u )+n f. Use u instead of u for implicit time-stepping. Henshaw (LLNL) FSI on overlapping grids VaTech 8 / 46
16 Incompressible Navier-Stokes. Split-step, velocity-pressure formulation: u t +(u )u + p ν u f = 0, t > 0, x Ω p u : u α u f = 0, t > 0, x Ω Divergence damping term: α u is important. Wall boundary conditions: u = 0, u = 0, (pressure BC) x Ω, with numerical boundary condition: p n = n ( ν u )+n f. Use u instead of u for implicit time-stepping. Henshaw (LLNL) FSI on overlapping grids VaTech 8 / 46
17 Incompressible Navier-Stokes. Split-step, velocity-pressure formulation: u t +(u )u + p ν u f = 0, t > 0, x Ω p u : u α u f = 0, t > 0, x Ω Divergence damping term: α u is important. Wall boundary conditions: u = 0, u = 0, (pressure BC) x Ω, with numerical boundary condition: p n = n ( ν u )+n f. Use u instead of u for implicit time-stepping. Henshaw (LLNL) FSI on overlapping grids VaTech 8 / 46
18 Incompressible Navier-Stokes. Split-step, velocity-pressure formulation: u t +(u )u + p ν u f = 0, t > 0, x Ω p u : u α u f = 0, t > 0, x Ω Divergence damping term: α u is important. Wall boundary conditions: u = 0, u = 0, (pressure BC) x Ω, with numerical boundary condition: p n = n ( ν u )+n f. Use u instead of u for implicit time-stepping. Henshaw (LLNL) FSI on overlapping grids VaTech 8 / 46
19 PDE based numerical boundary conditions (NBCs) High-order FD schemes with wide stencils require additional NBCs. One-sided difference approximations can be less stable and accurate. NBC Principle Derive NBCs from the PDE and physical boundary conditions. Example (heat equation): u t = ν(u xx + u yy )+f (x, y, t), 0 < x < 1, < y <, u(0, y, t) =g(y, t), u x (1, y, t) =k(y, t). Compatibility conditions are used as NBCs: (Dirichlet boundary condition), (Neumann boundary condition), g t (y, t) =ν(u xx (0, y, t)+g yy (0, y, t)) + f (0, y, t), k t (y, t) =ν(u xxx (1, y, t)+k yy (1, y, t)) + f x (1, y, t). Normal mode theory shows why compatibility conditions can be applied to lower order accuracy. Henshaw (LLNL) FSI on overlapping grids VaTech 9 / 46
20 PDE based numerical boundary conditions (NBCs) High-order FD schemes with wide stencils require additional NBCs. One-sided difference approximations can be less stable and accurate. NBC Principle Derive NBCs from the PDE and physical boundary conditions. Example (heat equation): u t = ν(u xx + u yy )+f (x, y, t), 0 < x < 1, < y <, u(0, y, t) =g(y, t), u x (1, y, t) =k(y, t). Compatibility conditions are used as NBCs: (Dirichlet boundary condition), (Neumann boundary condition), g t (y, t) =ν(u xx (0, y, t)+g yy (0, y, t)) + f (0, y, t), k t (y, t) =ν(u xxx (1, y, t)+k yy (1, y, t)) + f x (1, y, t). Normal mode theory shows why compatibility conditions can be applied to lower order accuracy. Henshaw (LLNL) FSI on overlapping grids VaTech 9 / 46
21 Cgins: an efficient solver for the incompressible Navier-Stokes equations A parallel split-step solver is being developed based on: 1 Fourth-order accurate approximate-factored/compact time-stepping scheme for the momentum equations. 2 Fourth-order accurate multigrid solver for the pressure equation. 3 Fast overlapping grid generation for moving geometry. WDH., A Fourth-Order Accurate Method for the Incompressible Navier-Stokes Equations on Overlapping Grids, J. Comput. Phys, (1994). WDH, N.A. Petersson, A Split-Step Scheme for the Incompressible Navier-Stokes Equations, WDH., K. K. Chand, A Composite Grid Solver for Conjugate Heat Transfer in Fluid-Structure Systems, J. Comput. Phys, Henshaw (LLNL) FSI on overlapping grids VaTech 10 / 46
22 Approximate factorization & compact discretizations - a key point is maintaining accuracy at boundaries. Approximate factorization (AF) schemes offer larger timesteps with second order accuracy in time: becomes (I + t 2 (A + B))Un+1 =(I t (A + B))Un 2 (I + t t A)(I B)Un+1 = (I t t A)(I 2 2 B)Un Compact schemes can be integrated into the AF solves Special combined compact schemes have been developed: reduce the number of factors (a r + b 2 rr )u P 1 (D r a + D rr b)u + corrections preserve accuracy at boundaries 4 th and 6 th order accuracy with a 5 point stencil special penta-diagonal solvers that handle wider boundary stencils Henshaw (LLNL) FSI on overlapping grids VaTech 11 / 46
23 Smallest Scale LES (SSLES) model Our high-order accurate schemes all use central differences. Nonlinear dissipation is used to stablize the schemes. It can be proved [HKR] that the minimum scale of the INS locally satisfies ν λ min u loc + c. Setting h = λ min leads to the non-linear dissipation operators: D 2 (u i )=(a 21 + a 22 h u i )h 2 h u i, D 4 (u i )=(a 41 + a 42 h u i )h 4 ( 2 h)u i. 1 These serve as simple LES models (Smagorinsky type). 2 These terms stabilize the scheme even if ν = 0. [HKR1] On the smallest scale for the incompressible Navier-Stokes equations, WDH, H.-O. Kreiss, L.G.M. Reyna, Theoret. Comput. Fluid Dynamics, [HKR2] Smallest scale estimates for the incompressible Navier-Stokes equations, WDH, H.-O. Kreiss, L.G.M. Reyna, Arch. Rational Mech. Anal., Henshaw (LLNL) FSI on overlapping grids VaTech 12 / 46
24 Code verification: an integral part of our development Twilight zone solutions (a.k.a. manufactured solutions) are extensively supported by the Overture framework. h max e p e u e v e w u 1.34e e e e e e e e-02 (15.9) 1.84e-02 (18.6) 1.01e-02 (21.2) 6.43e-03 (19.1) 2.71e-01 (10.3) 3.34e e-03 (11.8) 1.49e-03 (12.3) 7.50e-04 (13.5) 3.01e-04 (21.4) 3.10e-02 (8.7) 1.67e e-05 (12.9) 1.25e-04 (11.9) 6.55e-05 (11.5) 1.74e-05 (17.3) 4.00e-03 (7.8) rate Henshaw (LLNL) FSI on overlapping grids VaTech 13 / 46
25 AFS performance: flow past two cylinders. Grid tcilce64, N g = 22M grid-points Method P-solver N p TPSM TTS RPG PC24 BICGS(5) PC24 MG AFS24 MG Speed-up from baseline PC24-BICGS(5) to AFS24-MG is 16000/ Memory is reduced by a factor of 205/32 6. N p : number of processors. N g : number of grid-points. TPS : CPU (s) time-per-step. TPSM : CPU time (s) per-step, per-million-grid-points, per-processor = N p TPS/(N g/10 6 ). TTS : normalized time-to-solution = TPSM / CFL-number. RPG : memory usage in real-per-grid-point. PC24: predictor-corrector, 2nd-order time, 4th-order space. AFS24: approximate-factored scheme, 2nd-order time, 4th-order space. BICGS(5) : bi-conjugate-gradient-stabilized with ILU(5) preconditionner. Henshaw (LLNL) FSI on overlapping grids VaTech 14 / 46
26 Fourth-order Accurate Parallel Multigrid Solver for the Pressure Equation. Henshaw (LLNL) FSI on overlapping grids VaTech 15 / 46
27 The Ogmg overlapping grid multigrid solver has been extended to 4th-order accuracy and parallel. Ogmg is many times faster than Krylov methods. matrix-free; optimized for Cartesian grids. automatic coarse grid generation. adaptive smoothing variable sub-smooths per component grid. interpolation-boundary smoothing (IBS). Galerkin coarse grid operators (operator averaging). PDE-based numerical boundary conditions for Dirichlet and Neumann problems. WDH., On Multigrid For Overlapping Grids, SIAM J. Sci. Comput. 26, no. 5, (2005) Henshaw (LLNL) FSI on overlapping grids VaTech 16 / 46
28 Automatic coarse grid generation is a key feature. overlap increases interpolation accuracy reduced Henshaw (LLNL) FSI on overlapping grids VaTech 17 / 46
29 Local Fourier analysis significantly improves convergence rates. Over-relaxed Red-Black smoothers and Galerkin coarse grid operators D Poisson, Order 4, RB GS ρ V[1,1] 3G ρ V[1,1] Galerkin 3G ρ 3G V[2,1] ρ 3G V[2,1] Galerkin default rates (ω = 1) optimal rates 3D Poisson, Order 4, RB GS 0.35 ρ 3G V[1,1] ρ 0.3 3G V[1,1] Galerkin ρ 3G V[2,1] 0.25 ρ V[2,1] Galerkin 3G ω ω Three-grid multigrid convergence rates as a function ω. ω : relaxation parameter in Red-Black Gauss-Seidel smoother. ρ 3G : convergence rate per cycle for a 3 grid (i.e. 3 level) MG. V [m, n] : MG V-cycle, m pre-smooths and n-post smooths. Galerkin : Galerkin coarse grid operators. Henshaw (LLNL) FSI on overlapping grids VaTech 18 / 46
30 Accuracy and convergence of the new fourth-order accurate parallel version of Ogmg. Maximum error Dirichlet Neumann Sphere in a Box, Errors versus h Slope 4 maximum residual 10 4 Sphere in a Box, Order 4, V[2,1] CR=0.246, ECR=0.74 CR=0.058, ECR=0.61 no IBS IBS h multigrid cycle Henshaw (LLNL) FSI on overlapping grids VaTech 19 / 46
31 Multigrid is much faster than Krylov based methods. And uses much less memory. Poisson equation, SECOND-order accuracy Poisson equation, FOURTH-order accuracy Ogmg vs Krylov sq2048 cic64 tcilc64 box256 sphere4 sphere8 Ogmg vs Krylov sq2048 cic64 tcilc64 box256 sphere4 sphere8 speed-up memory-factor speed-up memory-factor Performance of Ogmg for solving Poisson s equation to on various grids. Ogmg is compared to a Krylov solver (bi-cg stab, ILU(1)) from PETSc. Henshaw (LLNL) FSI on overlapping grids VaTech 20 / 46
32 Rapid grid generation for moving geometry is critical. Ogen is the overlapping grid generator in Overture. Overlapping grid generation consists of two major steps: 1 construct component grids (Mappings). 2 grid connectivity: cut holes and determine interpolation information using Ogen (this is the step that requires most of the CPU time). In recent work changes have been made to support the generation of large (billion point +) grids. parallel moving-grid flow simulations (Ogen is called at each time step). Grid order of grid points processors cpu (s) accuracy (nodes p/n) Sphere in a box billion 16 (8 2) 136 Re-entry vehicle million 128 (16 8) 1990 Henshaw (LLNL) FSI on overlapping grids VaTech 21 / 46
33 Rapid grid generation for moving geometry is critical. Ogen is the overlapping grid generator in Overture. Overlapping grid generation consists of two major steps: 1 construct component grids (Mappings). 2 grid connectivity: cut holes and determine interpolation information using Ogen (this is the step that requires most of the CPU time). In recent work changes have been made to support the generation of large (billion point +) grids. parallel moving-grid flow simulations (Ogen is called at each time step). Grid order of grid points processors cpu (s) accuracy (nodes p/n) Sphere in a box billion 16 (8 2) 136 Re-entry vehicle million 128 (16 8) 1990 Henshaw (LLNL) FSI on overlapping grids VaTech 21 / 46
34 FSI for compressible fluids and deforming elastic solids (partitioned algorithms). Overcoming the added-mass instablity. Henshaw (LLNL) FSI on overlapping grids VaTech 22 / 46
35 Deforming composite grids (DCG) for Fluid-Structure Interactions (FSI) Goal: To perform coupled simulations of compressible fluids and deforming solids. A mixed Eulerian-Lagrangian approach: Fluids: general moving coordinate system with overlapping grids. Solids : fixed reference frame with overlapping-grids (later: unstructured-grids, or beam/plate models). Boundary fitted deforming grids for fluid-solid interfaces. Strengths of the approach: maintains high quality grids for large deformations/displacements. efficient structured grid methods (AMR) optimized for Cartesian grids. Henshaw (LLNL) FSI on overlapping grids VaTech 23 / 46
36 Deforming Composite Grids for FSI t=1.0 t=1.5 t=2.0 solid fluid interfaces solid Solve Euler equations in the fluid domains on moving grids. Solve equations of linear elasticity in the solid domains. Fluid grids at the interface deform over time (hyperbolic grid generator). Henshaw (LLNL) FSI on overlapping grids VaTech 24 / 46
37 The traditional partitioned algorithm for FSI 1 Advance solid using the stress from the fluid. 2 Advance fluid using the solid position/velocity. Note: traditional scheme has trouble for light solids (e.g. becomes unstable or requires a smaller time-step). this is known as an added mass instability. Added mass: The force required to accelerate a body in a fluid is greater than that in a vacuum (since both body and fluid must be accelerated). Henshaw (LLNL) FSI on overlapping grids VaTech 25 / 46
38 The elastic piston - a 1D FSI model problem solid fluid ρ : density ū : displacement v : velocity σ : stress c p : speed of sound z = ρc p : impedance ρ : density v : velocity σ = p : stress, pressure a : speed of sound z = ρa : impedance Henshaw (LLNL) FSI on overlapping grids VaTech 26 / 46
39 The elastic piston x = 0 x = F (t) t C + x = c pt ū 0 ( x) v 0 ( x) σ 0 ( x) solid (F (τ),τ) (x, t) fluid x = a 0 t ρ 0 v 0 p 0 x The governing equations for the solid and fluid are ū t v = 0 v t σ x / ρ = 0, for x < 0, σ t ρc p 2 v x = 0 where ρe = p/(γ 1)+ρv 2 /2. The interface conditions are v(0, t) =v(f (t), t), ρ t +(ρv) x = 0 (ρv) t +(ρv 2 + p) x = 0, for x > F (t), (ρe) t +(ρev + pv) x = 0 σ(0, t) =σ(f (t), t) p(f (t), t)+p e. Henshaw (LLNL) FSI on overlapping grids VaTech 27 / 46
40 The elastic piston x = 0 x = F (t) t C + x = c pt ū 0 ( x) v 0 ( x) σ 0 ( x) solid (F (τ),τ) (x, t) fluid x = a 0 t ρ 0 v 0 p 0 x The governing equations for the solid and fluid are ū t v = 0 v t σ x / ρ = 0, for x < 0, σ t ρc p 2 v x = 0 where ρe = p/(γ 1)+ρv 2 /2. The interface conditions are v(0, t) =v(f (t), t), ρ t +(ρv) x = 0 (ρv) t +(ρv 2 + p) x = 0, for x > F (t), (ρe) t +(ρev + pv) x = 0 σ(0, t) =σ(f (t), t) p(f (t), t)+p e. Henshaw (LLNL) FSI on overlapping grids VaTech 27 / 46
41 The solution of the fluid-solid Riemann (FSR) problem defines our interface projection. x = c pt v0 σ 0 x = 0 x = F (t) v1 σ 1 solid t ρ 1 v 1 p 1 fluid x =(v 1 + a 1 )t x = St ρ 0 v 0 p 0 x =(v 0 + a 0 )t special case of elastic piston problem - constant initial conditions in fluid and solid. the fluid may have a shock or expansion fan on the C + characteristic. x Henshaw (LLNL) FSI on overlapping grids VaTech 28 / 46
42 Solution of the linearized fluid-solid Riemann problem Characteristic relations: Solid: z v σ = z v 0 σ 0, on d x/dt = ±c p, Fluid: zv σ = zv 0 σ 0, on dx/dt = v 0 ± a 0, where z = ρc p and z = ρa 0 are the acoustic impedances. The state next to the interface is an impedance weighted average of the fluid and solid states: v 1 = v 1 = z v 0 + zv 0 z + z + σ 0 σ 0 z + z, σ 1 = σ 1 = z 1 σ 0 + z 1 σ 0 z 1 + z 1 + v 0 v 0 z 1 + z 1. The solution to the full nonlinear problem can also be determined. Henshaw (LLNL) FSI on overlapping grids VaTech 29 / 46
43 The FSI-DCG time-stepping algorithm 1 Advance solid to give provisional interface values ( v 0, σ 0 ). 2 Advance fluid to give provisional interface values (v 0,σ 0 ). 3 project the interface values: Notes: v I = z v 0 + zv 0 z + z + σ 0 σ 0 z + z, σ I = z 1 σ 0 + z 1 σ 0 z 1 + z 1 + v 0 v 0 z 1 + z 1 One can prove that the new scheme remains stable under the usual CFL time step restriction, even for very light solids [Banks/Sjögreen 2011]. The standard FSI scheme uses the heavy solid limit, z z, v I = v 0 σ I = σ 0 and is unstable for light solids (i.e. requires a small time step). Henshaw (LLNL) FSI on overlapping grids VaTech 30 / 46
44 The FSI-DCG time-stepping algorithm 1 Advance solid to give provisional interface values ( v 0, σ 0 ). 2 Advance fluid to give provisional interface values (v 0,σ 0 ). 3 project the interface values: Notes: v I = z v 0 + zv 0 z + z + σ 0 σ 0 z + z, σ I = z 1 σ 0 + z 1 σ 0 z 1 + z 1 + v 0 v 0 z 1 + z 1 One can prove that the new scheme remains stable under the usual CFL time step restriction, even for very light solids [Banks/Sjögreen 2011]. The standard FSI scheme uses the heavy solid limit, z z, v I = v 0 σ I = σ 0 and is unstable for light solids (i.e. requires a small time step). Henshaw (LLNL) FSI on overlapping grids VaTech 30 / 46
45 The FSI-DCG time-stepping algorithm 1 Advance solid to give provisional interface values ( v 0, σ 0 ). 2 Advance fluid to give provisional interface values (v 0,σ 0 ). 3 project the interface values: Notes: v I = z v 0 + zv 0 z + z + σ 0 σ 0 z + z, σ I = z 1 σ 0 + z 1 σ 0 z 1 + z 1 + v 0 v 0 z 1 + z 1 One can prove that the new scheme remains stable under the usual CFL time step restriction, even for very light solids [Banks/Sjögreen 2011]. The standard FSI scheme uses the heavy solid limit, z z, v I = v 0 σ I = σ 0 and is unstable for light solids (i.e. requires a small time step). Henshaw (LLNL) FSI on overlapping grids VaTech 30 / 46
46 Verification: elastic piston numerical results Computed solution for a smoothly receding piston. Fluid Solid grid N ρ r v r p r ū r v r σ r G e-3 2.7e-3 1.6e-3 4.9e-4 9.3e-5 1.2e-4 G e e e e e e G e e e e e e G e e e e e e rate Max-norm errors for very light solid: ρ = 10 5, z/ z = Fluid Solid grid N ρ r v r p r ū r v r σ r G e-4 9.6e-4 8.5e-4 3.0e-4 1.3e-5 1.1e-5 G e e e e e e G e e e e e e G e e e e e e rate Max-norm errors for very heavy solid: ρ = 10 5, z/ z = Max-norm errors are converging to second-order accuracy. Scheme is stable for large and small impedance ratios. Henshaw (LLNL) FSI on overlapping grids VaTech 31 / 46
47 The superseismic shock FSI problem y shock w 1 ξ w 0 fluid θ S solid x interface w 2 s-wave η s w 1 η p p-wave w 0 a fluid shock moving from left to right deflects a solid interface. an FSI problem with an analytic traveling wave solution. Henshaw (LLNL) FSI on overlapping grids VaTech 32 / 46
48 Superseismic shock: grids and computed solution ρ v 2 The red grid for the solid domain is shown adjusted for the displacement. Henshaw (LLNL) FSI on overlapping grids VaTech 33 / 46
49 L 1 -norm errors and convergence rates for the superseismic shock Solid Fluid Grid ū r v r σ r ρ r v r T r G (4) ss 8.6e-5 2.8e-3 7.9e-3 4.6e-3 2.2e-2 1.1e-2 G (8) ss 3.6e e e e e e G (16) ss 1.5e e e e e e G (32) ss 5.6e e e e e e rate Due to the discontinuities, the L 1 -norm errors do not converge at second-order. The solid variables v and σ converge at the expected rates of 2/3 (due to spreading of the linear discontinuities). In isolation the fluid domain should converge at first order. Henshaw (LLNL) FSI on overlapping grids VaTech 34 / 46
50 FSI for compressible fluids and light rigid bodies (partitioned algorithms). Overcoming the added-mass instability. Henshaw (LLNL) FSI on overlapping grids VaTech 35 / 46
51 Light rigid body in a fluid - model problem w b /2 0 w b /2 x fluid (acoustics) solid (rigid body) fluid (acoustics)... 3 vl vr σ L Fluid domains (acoustics): t i vk σ k v b σ R 0 1 ρ vk ρc 2 = 0, k=l,r, 0 x σ k Rigid body: (Question: what happens when m b 0?) i Interface conditions: m b v b = F, F = σ R x=wb /2 σ L x= wb /2 v L x= wb /2 = v b, v R x=wb /2 = v b. Henshaw (LLNL) FSI on overlapping grids VaTech 36 / 46
52 The fluid-solid Riemann problem for a rigid body. x = r b (t) t v b,σ b C : σ + zv = σ0 + zv0 f b body fluid v 0,σ 0 x From characteristics we obtain the fluid force on the body: σ(r b, t) =σ(r b +, t )+z v(r b +, t ) v b (t). (1) where z = ρc is the fluid impedence. Equation (1) defines our stress projection. Note the dependence of the stress on the velocity of the body v b. Henshaw (LLNL) FSI on overlapping grids VaTech 37 / 46
53 Light rigid body in a fluid - algorithm First order upwind scheme n+1 vk σ k i = n vk σ k i n + tr k Λ vk k R 1 k D σ k i n + tr k Λ + vk k R 1 k D +, σ k i Partitioned time-stepping algorithm (new scheme: α = z, traditional: α = 0), n+1 vk 1 Advance fluid domans: σ k i 2 Set F n+1 = σ n+1 R,1 n+1 + α(vr,1 v n+1 b ) σ n+1 n+1 L, 1 α(vb v n+1 L, 1 ), 3 Solve for rigid body: m b v n+1 b = m b v n b + tf n+1, 1 v n+1 b = m b + 2 tα m b vb n + t σ n+1 R,1 σn+1 n+1 L, 1 + α(vr,1 + v n+1 L, 1. ) 4 Set ghost values from projection. Note the added mass term 2 tα. Henshaw (LLNL) FSI on overlapping grids VaTech 38 / 46
54 Normal mode stability analysis Theorem The first-order upwind scheme with projection parameter α = z is stable for λ = c t x < 1 and any m b 0. Theorem The first-order upwind scheme with projection parameter α = 0 has the time-step restriction λ<1 and t < m b (4 λ)/(zλ). Note 1: The projection scheme (α = z) remains stable with the usual time-step restriction for any m b 0. Note 2: The traditional scheme (α = 0) is formally stable for m b > 0 but the time-step t goes to zero at m b 0. Henshaw (LLNL) FSI on overlapping grids VaTech 39 / 46
55 The Newton-Euler equations for rigid body dynamics The equations of motion for the rigid body are the Newton-Euler equations: ẋ b = v b, m b v b = F, A ω = WAω + T, Ė = WE. m b : mass of the body, x b (t) R 3 : position of the center of mass, v b (t) R 3 : is the velocity of the center of mass, ω(t) R 3 : is the angular velocity, A(t) R 3 3 : moment of inertia matrix, E(t) R 3 3 : matrix with columns being the principle axes of inertia, W (t) R 3 3 : angular velocity matrix, W a = ω a, F(t) R 3 : resultant force on the body, T (t) R 3 : resultant torque the body. Henshaw (LLNL) FSI on overlapping grids VaTech 40 / 46
56 Extension to 3D: added mass matrices For a point r on the rigid body in 3D, the interface projection is p(r, t)n = p f n + z f n T v f v(r, t) n. Force and torque: (y = r x b (t), v(r, t) =v b (t) y(r, t) ω(t)) F = p(r, t)n ds = A vv v b A vω ω + F, r B T = y ( p(r, t)n) ds = A ωv v b A ωω ω + T r B where A ij are the added-mass matrices A vv = z f nn T ds, A vω = z f n(y n) T ds, (2) B B A ωv = z f (Y n)n T ds A ωω = z f Y n(y n) T ds, (3) B B and where Y is the matrix corresponding to y. Henshaw (LLNL) FSI on overlapping grids VaTech 41 / 46
57 Extension to 3D: added mass matrices For a point r on the rigid body in 3D, the interface projection is p(r, t)n = p f n + z f n T v f v(r, t) n. Force and torque: (y = r x b (t), v(r, t) =v b (t) y(r, t) ω(t)) F = p(r, t)n ds = A vv v b A vω ω + F, r B T = y ( p(r, t)n) ds = A ωv v b A ωω ω + T r B where A ij are the added-mass matrices A vv = z f nn T ds, A vω = z f n(y n) T ds, (2) B B A ωv = z f (Y n)n T ds A ωω = z f Y n(y n) T ds, (3) B B and where Y is the matrix corresponding to y. Henshaw (LLNL) FSI on overlapping grids VaTech 41 / 46
58 Multi-dimensional algorithm Scheme: 1 second-order accurate Godunov based scheme for the (reactive) Euler equations. 2 Implicit Runge-Kutta (DIRK) time-stepping for added-mass rigid-body equations. Notes: 1 Second-order accurate in the max norm for smooth solutions (with limiter turned off). 2 First-order accurate in the discrete L 1 -norm for problems with shocks. 3 Adaptive mesh refinement. 4 Equations are solved in the moving grid coordinate system. Henshaw (LLNL) FSI on overlapping grids VaTech 42 / 46
59 Shock driven zero mass ellipse t = t = t = p p p Henshaw (LLNL) FSI on overlapping grids VaTech 43 / 46
60 Shock driven zero mass ellipse 2.5 x 1 Shock driven ellipse, M=0 2 x 2 v 1 v w t Time histories of rigid-body properties. Colour: grid G re (8), black: fine grid G re (32) Grid G (j) h j e ρ (j) r e (j) u r e (j) v r e (j) p r G re (8) 1/40 2.1e-3 9.3e-4 9.6e-4 2.1e-3 G re (16) 1/80 9.9e e e e G re (32) 1/ e e e e rate A posteriori estimated errors (L 1 -norm) and convergence rates for the accelerated ellipse at t = 1.0. Henshaw (LLNL) FSI on overlapping grids VaTech 44 / 46
61 Shock driven zero mass ellipse Coarse grid G re (32) (left), medium grid G re (16 4) (middle), fine (AMR) grid G re (8 4 4) (right). Henshaw (LLNL) FSI on overlapping grids VaTech 45 / 46
62 Summary We have been developing efficient algorithms for modeling 1 incompressible flows and rigid body motion, 2 compressible flows and rigid/deforming bodies. The approach is based on 1 overlapping grids for flexible representation of geometry, 2 high-order accurate algorithms, 3 accurate treatment of interfaces and boundary conditions. 4 solutions to fluid-solid Riemann problems to overcome added-mass instabilities. Henshaw (LLNL) FSI on overlapping grids VaTech 46 / 46
Solving Fluid Structure Interaction Problems on Overlapping Grids
Solving Fluid Structure Interaction Problems on Overlapping Grids Bill Henshaw Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA, USA. Massachusetts Institute
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 informationRecent Developments in Overture
Recent Developments in Overture Bill Henshaw Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA, USA. 11th Symposium on Overset Grids and Solution Technology,
More informationDeforming Composite Grids for Solving Fluid Structure Problems
Deforming Composite Grids for Solving Fluid Structure Problems Jeffrey W. Banks a,1, William D. Henshaw a,1,, Donald W. Schwendeman b,2 a Centre for Applied Scientific Computing, Lawrence Livermore National
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 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 informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
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 informationNumerical Methods for Solid Mechanics on Overlapping Grids: Linear Elasticity
Manuscript Click here to view linked References Numerical Methods for Solid Mechanics on Overlapping Grids: Linear Elasticity Daniel Appelö a,, Jeffrey W. Banks b,, William D. Henshaw b,,, Donald W. Schwendeman
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 informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
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 informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
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 informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
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 informationNavier Stokes Solvers. Length (Log meters) Particle Methods. Density (Log n/n 0. Knudsen number =1/10 2. Dense Gas
Mesh and Algorithm Reænement Alejandro L. Garcia Center for Comp. Sciences and Engineering Lawrence Berkeley National Laboratory and Dept. Physics, San Jose State University Collaborators: Berni Alder
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
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 informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationIterative Solvers in the Finite Element Solution of Transient Heat Conduction
Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Mile R. Vuji~i} PhD student Steve G.R. Brown Senior Lecturer Materials Research Centre School of Engineering University of
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
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 information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
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 informationA multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction
A multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction J Degroote, S Annerel and J Vierendeels Department of Flow, Heat and Combustion Mechanics, Ghent University,
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
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 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 informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
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 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 Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
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 informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
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 informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationEfficient wave propagation on complex domains
Center for Turbulence Research Annual Research Briefs 2006 223 Efficient wave propagation on complex domains By K. Mattsson, F. Ham AND G. Iaccarino 1. Motivation and objectives In many applications, such
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationTAU Solver Improvement [Implicit methods]
TAU Solver Improvement [Implicit methods] Richard Dwight Megadesign 23-24 May 2007 Folie 1 > Vortrag > Autor Outline Motivation (convergence acceleration to steady state, fast unsteady) Implicit methods
More informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More informationSplitting based Implicit solvers for compressible fluid models
Splitting based Implicit solvers for compressible fluid models D. Coulette 3, E. Franck 1, M. Gaja 2, P. Helluy 3, J. Lakhlili 2, M. Mazza 2, M. Mehrenberger 3, A. Ratnani 2, S. Serra-Capizzano 4, E. Sonnendrücker
More informationModeling and numerical simulations of swimmers
Plafrim, may 31, 2011 p. 1 Modeling and numerical simulations of swimmers Michel Bergmann INRIA Bordeaux Sud-Ouest, project-team MC2 Institut de Mathématiques Appliquées de Bordeaux 33405 TALENCE cedex,
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 informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
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 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 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 informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More information