Deforming Composite Grids for Fluid Structure Interactions
|
|
- Rodney Holt
- 6 years ago
- Views:
Transcription
1 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, CA, USA. 2 Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY, USA. SIAM Conference on Computational Science and Engineering, Reno, Nevada, February 28 - March 4, Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
2 Acknowledgments. 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) Deforming Composite Grids for FSI CSE / 28
3 Outline 1 Background: overlapping grids, Overture and CG. 2 Deforming Composite Grids (DCG) for fluid structure interactions. 3 The elastic-piston problem and the fluid-solid Riemann problem. 4 Stable interface conditions for coupling the Euler equations and the elastic wave equation. 5 Verification problems. 1 super-seismic shock 2 deforming diffuser. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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) Deforming Composite Grids for FSI CSE / 28
5 Previous work: shock hitting rigid cylinders with AMR (Shock hitting rigid cylinders with AMR) randomcylrhowithgrids.mpg Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
6 Current work: shock hitting elastic cylinders (Shock hitting elastic cylinders) shockmultidisk.mpg Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
12 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. Efficient for high-order accurate methods. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
13 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) Deforming Composite Grids for FSI CSE / 28
14 A sample FSI-DCG simulation Gas Solid Mach 2 shock in a gas hitting two elastic cylinders. 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). Adaptive mesh refinement (in progress). Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
15 Fluid and Solid Solvers Fluid solver: we solve the inviscid Euler equations with a second-order extension of Godunov s method (cgcns). WDH., D. W. Schwendeman, Parallel Computation of Three-Dimensional Flows using Overlapping Grids with Adaptive Mesh Refinement, J. Comp. Phys. 227 (2008). WDH., DWS, Moving Overlapping Grids with Adaptive Mesh Refinement for High-Speed Reactive and Nonreactive Flow, J. Comp. Phys. 216 (2005). WDH., DWS, An adaptive numerical scheme for high-speed reactive flow on overlapping grids, J. Comp. Phys. 191 (2003). Solid solver: we solve the elastic wave equation as a first order system with a second-order upwind characteristic scheme (cgsm). Daniel Appelö, JWB, WDH, DWS, "Numerical Methods for Solid Mechanics on Overlapping Grids: Linear Elasticity, LLNL-JRNL , submitted (2010). Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
16 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) Deforming Composite Grids for FSI CSE / 28
17 The elastic piston x = 0 x = F(t) t C + x = c pt 2 3 ū 0 ( x) 4 v 0 ( x) 5 σ 0 ( x) solid (F(τ), τ) (x, t) fluid x = a 0 t 2 4 ρ 3 0 v 0 5 p 0 x The governing equations for the solid and fluid are 8 ū >< t v = 0 v t σ x / ρ = 0, for x < 0, >: σ t ρc p v x 2 = 0 where ρe = p/(γ 1) + ρv 2 /2. The interface conditions are 8 ρ >< t + (ρv) x = 0 (ρv) t + (ρv >: 2 + p) x = 0, for x > F(t), (ρe) t + (ρev + pv) x = 0 v(0, t) = v(f(t), t), σ(0, t) = σ(f(t), t) p(f(t), t) + p e. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
18 The elastic piston x = 0 x = F(t) t C + x = c pt 2 3 ū 0 ( x) 4 v 0 ( x) 5 σ 0 ( x) solid (F(τ), τ) (x, t) fluid x = a 0 t 2 4 ρ 3 0 v 0 5 p 0 x The governing equations for the solid and fluid are 8 ū >< t v = 0 v t σ x / ρ = 0, for x < 0, >: σ t ρc p v x 2 = 0 where ρe = p/(γ 1) + ρv 2 /2. The interface conditions are 8 ρ >< t + (ρv) x = 0 (ρv) t + (ρv >: 2 + p) x = 0, for x > F(t), (ρe) t + (ρev + pv) x = 0 v(0, t) = v(f(t), t), σ(0, t) = σ(f(t), t) p(f(t), t) + p e. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
19 An exact solution to the elastic piston problem For a given x = F(t), and constant ρ 0, p 0, v 0 = 0, the solution in the fluid region F(t) < x < a 0 t is (assuming no shocks) v(x, t) = Ḟ(τ(x, t)), a(x, t) = 1 + γ 1 a 0 2! v(x, t), a 0 p(x, t) p 0 =» x F(τ) = a 0 + γ + 1 Ḟ(τ) (t τ). 2 The general solution for the solid follows from the d Alembert solution,! γ ρ(x, t) = ρ 0! 2γ/(γ 1) a(x, t), a 0 ū( x, t) = f( x c pt) + g( x + c pt), f(ξ) = 1 ˆū0 (ξ) 1 Z ξ v 0 (s)ds for ξ < 0, 2 c p 0 ( 12 ˆū0 (ξ) + 1 R ξ g(ξ) = cp 0 v 0(s)ds for ξ < 0, F(ξ/c p) f( ξ) for ξ > 0. Applying the interface equations gives an ODE for F(t) in terms of the initial conditions, p 0 h1 ρc p 2 + γ 1 i 2γ/(γ 1) Ḟ(t) Ḟ(t) + = ˆū 0 2a 0 c ( cpt) 1 v 0 ( c pt), for t > 0. p c p Alternatively if we choose F(t) = a q tq, we can choose initial conditions in the solid as ū 0 ( x) = p 0 ρ 0 c 2 p Z x to give a smooth solution with the specified interface motion. 0 h1 + γ 1 i 2γ/(γ 1) Ḟ( s/c p) ds, v0 ( x) = Ḟ( x/cp), for x < 0, 2a 0 Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
20 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 2 4 ρ 3 1 v 1 5 p 1 fluid x = (v 1 + a 1 )t x = St 2 4 ρ 3 0 v 0 5 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) Deforming Composite Grids for FSI CSE / 28
21 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) Deforming Composite Grids for FSI CSE / 28
22 The FSI-DCG time-stepping algorithm The FSI-DCG interface approximation is an extension of the scheme of Banks and Sjögreen: J. W. Banks and B. Sjögreen, A Normal Mode Stability Analysis of Numerical Interface Conditions for Fluid/Structure Interaction, Commun. Comput. Phys., The main steps are: 1 The fluid and solid domains are first advanced independently giving provisional interface values. 2 The provisional interface values are projected based on the solution to the fluid-solid Riemann problem. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
23 The FSI-DCG time-stepping algorithm solid fluid... w n 2 w n 1 w n 0 w n 1 x w n 1 w n 0 w n 1 w n 2... x Define discrete approximations, (v n i v(i x i, n t), v n i Solid: w n i = [ūi n, v i n, σ i n ] i= 1, 0, 1,... Fluid: w n i = [ρ n i, vi n, pi n ], i= 1, 0, 1,... v(i x i, n t)) Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
24 The FSI-DCG time-stepping algorithm stage condition type assigns 1a Predict grid and grid velocity extrapolation F p, Ḟ p 1b Advance w n i, wn i, i = 0, 1, 2,... PDE interior, interface 2a Eval v I, σ I, ρ I from FSR projection v I, σ I, ρ I 2b v0 n, v 0 n = v I, p0 n, σn 0 = σ I, ρ n 0 = ρ I projection w n 0, wn 0 2c w n 1 = E(3) +1 wn 0, wn 1 = Ē(3) +1 wn 0, extrapolation wn 1, wn 1 2d Eval : v f = 1 ρ x p,... PDE v f, v s, σ f, σ s 2e Eval v I, σ I from FSR projection v I, σ I 3a 1 ρ x p = v I, ρa 2 x v = σ I compatibility p 1 n, v 1 n 3b ρc p 2 x v = σ I, 1 ρ x σ = v I compatibility v 1 n, σn 1 3c Correct grid and grid velocity projection F n, Ḟ n Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
25 The interface projection step Using the linearized FSR solution, the interface values are an impedance weighted average of the provisional fluid and solid values: 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 Compare: the standard FSI scheme uses the heavy solid limit, z z, velocity-from-solid, stress-from-fluid: v I = v 0 σ I = σ 0 = p + p e The standard scheme is unstable for light solids. Note: for hard problems with shocks hitting the interface, there are advantages to using the full nonlinear solution to the FSR problem. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
26 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) Deforming Composite Grids for FSI CSE / 28
27 Two-dimensional verification problems 1 superseismic shock. 2 deforming diffuser. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
28 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) Deforming Composite Grids for FSI CSE / 28
29 Superseismic shock: grids and computed solution ρ v 2 The red grid for the solid domain is shown adjusted for the displacement. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
30 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) Deforming Composite Grids for FSI CSE / 28
31 The deforming diffuser (0, y b ) supersonic inflow (x d, y d ) supersonic inflow (0, F(0)) fluid domain Ω f interface supersonic outflow slip wall y = F(x) solid domain Ω s (1, F(1)) (0, y a) displacement BC slip wall (x c, y c) an FSI problem with a smooth semi-analytic solution. Fluid: Prandtl-Meyer analytic solution is defined as a function of F(x). Solid: steady elasticity equations are solved on a fine grid. The coupled exact solution and F(x) are determined by iteration. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
32 The deforming diffuser grid and solution Contours of the fluid pressure, [min, max] = [.5, 1.] and norm of the solid stress tensor σ, [min, max] = [0,.054]. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
33 Max-norm errors and convergence rates for the deforming diffuser Solid Grid ū r v r σ r G (2) dd 2.3e-4 8.6e-4 4.0e-2 G (4) dd 5.8e e e G (8) dd 9.9e e e G (16) dd 1.6e e e rate Fluid Grid ρ r v 1 r v 2 r T r G (2) dd 3.6e-2 1.7e-2 2.3e-2 1.2e-2 G (4) dd 8.8e e e e G (8) dd 2.1e e e e G (16) dd 5.0e e e e rate max-norm errors at t = 1 with the Godunov slope-limiter off. solutions are converging to second-order in the max norm. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
34 Summary the deforming composite grid (DCG) approach is being developed to to model fluid-structure interactions (FSI) for compressible fluids and elastic solids. the solution of the fluid-solid Riemann problem can be used to define stable interface approximations for both light and heavy solids. our FSI-DCG approximation was verfied to be second-order accurate in the max-norm on problems with smooth solutions. Future work: support for adaptive mesh refinement (AMR) at the interface. more general solid mechanics models. incompressible fluids. extend to three space dimensions. Henshaw (LLNL) Deforming Composite Grids for FSI CSE / 28
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 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 informationSolving 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. Virginia Tech, Blacksburg
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 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 informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
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 informationAn Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations. Thesis by Cheng Chi
An Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations Thesis by Cheng Chi In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah
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 information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
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 informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More 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 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 informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationA monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
1 / 41 A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion Chen-Yu CHIANG O. Pironneau, T.W.H. Sheu, M. Thiriet Laboratoire Jacques-Louis Lions
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
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 informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
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 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 informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
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 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 informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
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 informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
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 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 informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationLos Alamos National Laboratory Hydrodynamic Methods Applications and Research 1 LA-UR
Rayleigh-Taylor instability is generated when a heavy fluid sits above a lighter fluid in a gravitational field. The flow behavior is described in terms of bubbles of light fluid rising into the heavier
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
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 informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationSELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS
SIAM J. APPL. MATH. Vol. 0, No. 0, pp. 000 000 c XXXX Society for Industrial and Applied Mathematics SELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS ALLEN M. TESDALL, RICHARD SANDERS,
More informationEntry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design
Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
More informationSimulation of sympathetic detonation by a CIP Eulerian code
Computational Ballistics II 107 Simulation of sympathetic detonation by a CIP Eulerian code S. Kubota 1, Z. Liu 2, T. Saburi 1, Y. Ogata 1 & M. Yoshida 1 1 Research Center for Explosion Safety, National
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
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 informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
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 informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
More informationDETONATION DIFFRACTION, DEAD ZONES AND THE IGNITION-AND-GROWTH MODEL
DETONATION DIFFRACTION, DEAD ZONES AND THE IGNITION-AND-GROWTH MODEL G. DeOliveira, A. Kapila and D. W. Schwendeman Rensselaer Polytechnic Institute J. B. Bdzil Los Alamos National Laboratory W. D. Henshaw
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 informationDevelopment of an evaporation boundary condition for DSMC method with application to meteoroid entry
Development of an evaporation boundary condition for DSMC method with application to meteoroid entry F. Bariselli, S. Boccelli, A. Frezzotti, A. Hubin, T. Magin Annual METRO meeting 29th November 2016
More informationAMME2261: Fluid Mechanics 1 Course Notes
Module 1 Introduction and Fluid Properties Introduction Matter can be one of two states: solid or fluid. A fluid is a substance that deforms continuously under the application of a shear stress, no matter
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
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 informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationA Computational Study of High-Speed Droplet Impact
Copyright 2011 Tech Science Press FDMP, vol.7, no.4, pp.329-340, 2011 A Computational Study of High-Speed Droplet Impact T. Sanada 1, K. Ando 2 and T. Colonius 2 Abstract: When a droplet impacts a solid
More informationTHE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS
Journal of Hyperbolic Differential Equations Vol., No. 4 005 909 917 c World Scientific Publishing Company THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS VOLKER ELLING, and TAI-PING LIU, Dept.
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
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 informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
More information1 Introduction to PDE MATH 22C. 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers:
1 Introduction to PDE MATH 22C 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers: y = f(t) Output y 2R Input t 2R Name of function f t=independent
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
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 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 informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
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 informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationMAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring Dr. Jason Roney Mechanical and Aerospace Engineering
MAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring 2003 Dr. Jason Roney Mechanical and Aerospace Engineering Outline Introduction Kinematics Review Conservation of Mass Stream Function
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationParallel Computation of Forced Vibration for A Compressor Cascade
44th AIAA Aerospace Sciences Meeting and Exhibit 9-2 January 26, Reno, Nevada AIAA 26-628 AIAA Paper 26-628 Parallel Computation of Forced Vibration for A Compressor Cascade Zongjun Hu and Gecheng Zha
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationTruncation Error Estimates for Mesh Refinement in Lagrangian Hydrocodes Sarah Lianne Cole
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Truncation Error Estimates for Mesh Refinement in Lagrangian Hydrocodes Sarah Lianne Cole Thesis submitted for the degree of Doctor of Philosophy
More informationAerothermodynamics of high speed flows
Aerothermodynamics of high speed flows AERO 0033 1 Lecture 6: D potential flow, method of characteristics Thierry Magin, Greg Dimitriadis, and Johan Boutet Thierry.Magin@vki.ac.be Aeronautics and Aerospace
More information