MATHEMATICAL MODELING OF THE LONG TIME EVOLUTION OF THE PULSATING DETONATION WAVE IN THE SHOCK ATTACHED FRAME
|
|
- Austin Summers
- 5 years ago
- Views:
Transcription
1 International Conference Quasilinear equations, Inverse problems and their applications MATHEMATICAL MODELING OF THE LONG TIME EVOLUTION OF THE PULSATING DETONATION WAVE IN THE SHOCK ATTACHED FRAME Utkin P.S., Lopato A.I. Institute for Computer Aided Design Russian Academy of Science Moscow Institute of Physics and Technology The reported study was funded by RFBR according to the research project mol_а. Moscow Institute of Physics and Technology, Dolgoprudny, September 2016
2 What is detonation? Detonation is a hydrodynamic wave process of the supersonic propagation of an exothermic reaction through a substance. The dt detonation ti wave (DW) is a self sustained ti shock wave (SW) discontinuity it behind bhid the front of which a chemical reaction is continuously initiated due to heating caused by adiabatic compression. The detonation wave velocities in gaseous mixtures under normal conditions are about 1 3 km/s, front pressures atmospheres. Profile produced by gas motion behind the ideal detonation started at the closed end of the tube Sheperd J.E. // Proc. Comb. Inst Zeldovich von Neumann Doering (ZND) solution for the steady state detonation 2
3 Detonation theory: state of the art DW propagation is characterized by a complicated nonlinear oscillatory process 1D. Pulsations of parameters behind the DW front. 2D. Transverse compression waves that interact with the DW in two dimensional computations and experiments on the DWs propagation ation in narrow gaps. Detonation cells. 3D. Transverse wave propagating in a spiral spin. Leung C. et al. // Physics of Fluids Taylor B.D. et al. // Proc. Comb. Inst Levin V.A. et al. // Comp. Math. and Math. Phys
4 Detonation analogies and simplified models Hydraulic jump and traffic jam are analogous to self sustained DW (ZND like front structures) Kasimov A.R. // J. FluidMech Flynn M.R. et al. // Phys. Rev. E D asymptotic equations: Faria L.M. et al. // J. Fluid. Mech
5 Some problems in detonation calculations Possibility of DW failure in the simulations of the DW long time propagation Higgins i AJ A.J. Approaching dt detonation ti dynamics as an ensemble of interacting ti waves // Proc. 25th ICDERS. 2nd 7th August Leeds, UK. Paper PL3. He L., Lee J.H.S. The dynamical limit of one dimensional detonations // Physics of Fluids Semenov I. et al. Mathematical modeling of detonation initiation via flow cumulation effects // Progress in Propulsion Physics. 8. Proc. EUCASS July Munich, Germany. Why does the detonation wave fail in the numerical computations? Assumptions: Mistakes Mathematical models Computational method 5
6 Aims The aim of the work numerical investigation of weakly unstable and irregular regimes of pulsating DW propagation in two statements the modeling of DW in the laboratory frame () with detonation initiation near the closed end of the channel and modeling in the shock attached frame (),and quantitative i comparison of results using Fourier analysis of the pulsations. 6
7 Governing system of equations in the laboratory frame () Z mass fraction of the reacting mixture component Q heat release A pre exponent exponent factor E activation energy 1D Euler equations + one stage chemical reaction model 7
8 Governing system of equations in the shock attached frame () Z mass fraction of the reacting mixture component Q heat release A pre exponent factor E activation i energy 1D Euler equations + one stage chemical reaction model + shock velocity evolution equation 8
9 Brief review (shock attached frame) Kasimov A.R., Stewart D.S. On the dynamics of the self sustained onedimensional detonations: A numerical study in the shock attached frame // Physics of Fluids (10): Numerical scheme: first approximation order scheme. Shock speed evolution equation: integration of the governing equation along the C + characteristic nearthe shock. Applicability: test with a shock overtaking another shock demonstrates the stability of the algorithm to detonation waves numerical calculations. Results: main regimes of detonation wave propagation p are obtained. The shock dynamics is shown to be determined entirely by the finite region between the shock and the sonic locus. Henrick A.K., Aslam T.D., Powers J.M. Simulations of pulsating one dimensional detonations with true fifth order accuracy // Journal of Computational Physics V P Numerical scheme: fifth approximation order WENO scheme + fifth order Runge Kutta scheme. Shock speed evolution equation: connecting with the momentum flux gradient. Applicability: restrictions strongly unstable detonation waves. Results: the approximation order is shown to be equal to five (steady detonation). For an unstable regime a stable periodic limit cycle is obtained. The phase portrait confirms the unstable regime. The bifurcation diagram of different activation energies (different regimes) is constructed. 9
10 Dimensionless procedure Characteristic scales parameters in front of the DW and a half reaction length : l 12 pa paμ ρa pa ua = Ta = ρa ρar 1 uznd ( Z) DCJ = 12AZ exp EρZND Z pznd Z ( ( ) ( )) dz Dimensionless ρ p u = = = ˆ T pˆ = = ˆ E = ˆ Q ˆ ˆ ˆ = ˆ D = ˆ l ρ p u T E Q D l = 2 2 ρ p u T ρˆ μu u u l a a a a a a a 1 variables: ˆ ( ) ˆ ˆ A uznd Z DCJ A= = u l 12 12Zexp Eρ ˆˆ ( Z) pˆ ( Z) ( ) a ZND ZND dz 1/2 Dimensionless ZND profiles: Dˆ ˆ + ˆ CJ γ DCJ 1 DCJ γ uˆ 0( Z) = ( 1 + Z) pˆ 0( Z) = 1 + Z 2 γ+ 1 Dˆ + 1 ˆ CJ γ DCJ + 1 ˆ 2 ˆ 2 ( ) γ DCJ + 1 DCJ γ = ρˆ 0 Z 1 Z Dˆ = + ( ) ˆ + ( 2 ) ˆ 2 CJ γ γ 1 Q γ 1 Q γ + 1 Dˆ ( ˆ + ) CJ 2 2 γ DCJ 1 1 Erpenbeck J.J. Stability of steady state equilibrium detonations // Physics of Fluids V. 5. P Semenko R. et al. Set valued solutions for non ideal detonation // Shock waves V. 26. P
11 Computational algorithm Physical processes splitting technique Finite volume method Second approximation order ENO reconstruction Courant Isaacson Rees numerical flux Second order Runge Kutta explicit scheme for time stepping Euler implicit method for solving equations of energy and chemical reactions Parallelization (MPI) Algorithm of integration ti of the LSW velocity evolution () ( { } ) { } ( ) A i+ 12 { } { } CIR n n 1 n n n f i + = 12 f u i + f u i+ 1 + u + i ui { ( ) } + { ( ) } n = Ω n Λ n Ω n Ω n + 1 Λ n + 1 Ω n A i i i i i i i+ 1 2 Lopato A.I., Utkin P.S. // Comp. Research Model
12 Algorithm for shock speed calculation (1).. dξ = u+ c D dτ dp du + ρc ( γ 1) Qρω = 0 dτ dτ n ξ, a= a ξ, ξ, b= b ξ, ξ, c= c ξ, ξ. n n n n n n ( * * ) ( * * ) ( * * ) n 1 * ξ * 1 1 ( ) ( ) ( ) ( ) n n n n n n n n 2 a ξ*, ξ* τ + b ξ*, ξ* = u* + c* D, n n 1 n 1 n n 1 n 1 n 1 n 1 2 a ξ*, ξ * τ + b ξ*, ξ * = u * + c * D. The system is solved numerically with Newton iterations. 12
13 Algorithm for shock speed calculation (2).. dp du + ρc ( γ 1) Qρω = 0 dτ dτ New LSW speed D n+1 is determined dt dfrom the second equation. Three point one sided approximation of the derivatives: ( 0 * * ) ( 0 * * ) ( ) n+ 1 n n 1 n+ 1 n+ 1 n+ 1 n n 1 n+ 1 n+ 1 αp + λp + δp + ρ c αv + λv + δv γ 1 Qρ ω = 0 n n Δ τ Δ τ α= +, λ= + = 1 1, δ. n n n n n n n 1 Δτ Δ τ +Δτ Δτ Δτ Δ τ +Δτ Calculation of the SW velocity D n+1 = M n+1 c n+1 on (n+1) th time layer asaa solution of nonlinear equation. Here, Rankine Hugoniot conditions are applied n+ 1 n+ 1 ( γ )( M ) ( M ) n+ 1 n+ + ( γ )( M ) 2 2 n+ 1 n n n n+ 1 ( M ) Z p γ γ ρ v =, =, =, = 1. P γ γ R 2 1 C γ M n n 1 n n 1 Variables p*, p*, v*, v* are calculated using quadratic interpolation procedure. Lopato A.I., Utkin P.S. Detailed simulation of the pulsating detonation wave in the shock attached frame // Computational Mathematics and Mathematical Physics
14 Shock overtaking another shock test τ 14
15 Stable regime E = 25, Q = 50, γ = 1.2, Δx = calc1 N1= 4000, Δ 1= p front pvn.. calc k ln( Δ2 Δ1) pfront pvn 10. ~ Ch k 1.72 N = 2000, Δ = p p ln 2 calc front vn
16 Weakly unstable regime (1) E = 26, Q = 50, γ = 1.2, Δx = Δp ~ 10 4 T = T LIN = T = Lee H.I., HI Stewart t D.S. DS Cl Calculation lti of linear detonation dt ti instability: one dimensional i instability of plane detonation dt ti // Journal of Fluid Mechanics V P
17 Weakly unstable regime (2) E = 26, Q = 50, γ = 1.2, Δx = t τ Long term numerical research 17
18 Weakly unstable regime (3) E = 26, Q = 50, γ = 1.2, Δx = Fourier spectral analysis (5000 units of time) p = Ak cos 2 vn k= 0 p ( πν t ) k A A ν ν 18
19 Weakly unstable regime (4) Fourier spectra of two statements demonstrate similar values of peak frequencies A ν Peak ν() ν ()
20 Irregular regime (1) E = 28, Q = 50, γ = 1.2, Δx = Δp ~ 10 4 pmax 1.65 = pmin 0.76 p p max 1.68 = min
21 Irregular regime (2) E = 28, Q = 50, γ = 1.2, Δx = t t Long term numerical research 21
22 Irregular regime (3) E = 28, Q = 50, γ = 1.2, Δx = Fourier spectral analysis (5000units of time) p = Ak cos ( 2 πνkt ) p vn k= 0 A A ν ν 22
23 Irregular regime (4) Fourier spectra of two statements demonstrate similar values of peak frequencies Peak ν() ν () A ν 23
24 Strongly unstable regime (1) E = 35, Q = 50, γ = 1.2, Δx =
25 Strongly unstable regime (2) E = 35, Q = 50, γ = 1.2, Δx =
26 Conclusion 1. The numerical investigation of pulsating DW propagation using two statements e ts laboratory ato frame and dshock attached frame is carried out. The second statement includes the special numerical algorithm of the integration of shock speed evolution equation with the second approximation order based on the grid characteristic method. 2. The calculation of the weak unstable mode gives a DW that tends to the stable limit cycle. Both statements give similar results that correlate with the known analytical solution. The spectral Furrier analysis demonstrates that main frequency peaks coincide although the regions between the peaks somewhat differ. 3. Numerical investigation of the irregular mode demonstrates the tendency of the front pressure reducing for laboratory statement. The spectral Furrier analysis confirms chaospresence. 4. Further investigations long term calculation of strongly unstable mode for both statements. 26
27 Publications 1. Лопато, А.И., Уткин, П.С. Математическое моделирование пульсирующей волны детонации с использованием ENO схем различных порядков аппроксимации // Компьютерные исследования и моделирование Т. 6, 5. С Лопато, А.И., Уткин, П.С. Особенности расчета детонационной волны с использованием схем различных порядков аппроксимации // Математическое моделирование Т. 27, 7. С Лопато, А.И., Уткин, П.С. О двух подходах к математическому моделированию детонационной волны // Математическое моделирование Т. 28, 2. С Lopato, A.I., Utkin, PS P.S. Two approaches to the mathematical modeling of detonation wave // Mathematical Models and Computer Simulations V. 8, No. 5. P Лопато, А.И., Уткин, П.С. Детальное математическое моделирование пульсирующей детонационной волны в системе координат, связанной с лидирующим скачком // Журнал вычислительной математики и математической физики Т. 56, 5. С Lopato, A.I., Utkin, P.S. Detailed simulation of the pulsating detonation wave in the shockattached frame // Computational Mathematics and Mathematical Physics V. 56, No. 5. P Lopato, A.I., Utkin, P.S. Towards second order algorithm for the pulsating detonation wave modeling in the shock attached frame // Combustion Science and Technology V. 188, No. 11 (to appear). 27
Verified Calculation of Nonlinear Dynamics of Viscous Detonation
Verified Calculation of Nonlinear Dynamics of Viscous Detonation Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and Joseph
More informationThe Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing
The Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More informationAn analytical model for direct initiation of gaseous detonations
Issw21 An analytical model for direct initiation of gaseous detonations C.A. Eckett, J.J. Quirk, J.E. Shepherd Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena CA 91125,
More informationHigher Order DSD Calibration of Ammonium Nitrate/Fuel Oil
25 th ICDERS August 2 7, 2015 Leeds, UK Higher Order DSD Calibration of Ammonium Nitrate/Fuel Oil Carlos Chiquete, Mark Short, Scott I. Jackson and John B. Bdzil Shock and Detonation Physics Group Los
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 informationDetonation Diffraction
Detonation Diffraction E. Schultz, J. Shepherd Detonation Physics Laboratory Pasadena, CA 91125 MURI Mid-Year Pulse Detonation Engine Review Meeting February 10-11, 2000 Super-critical Detonation Diffraction
More informationInitiation of stabilized detonations by projectiles
Initiation of stabilized detonations by projectiles P. Hung and J. E. Shepherd, Graduate Aeronautical Laboratory, California Institute of Technology, Pasadena, CA 91125 USA Abstract. A high-speed projectile
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationRevised Manuscript. Applying Nonlinear Dynamic Theory to One-Dimensional Pulsating Detonations
Combustion Theory and Modelling Vol. 00, No. 00, July 10, 1 13 Revised Manuscript Applying Nonlinear Dynamic Theory to One-Dimensional Pulsating Detonations Hamid Ait Abderrahmane, Frederick Paquet and
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationNear limit behavior of the detonation velocity
Near limit behavior of the detonation velocity John H.S. Lee 1, Anne Jesuthasan 1 and Hoi Dick Ng 2 1 McGill University Department of Mechanical Engineering Montreal, QC, Canada 2 Concordia University
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationApplying the Staggered Mesh Godunov (SMG) Method to Reactive Flows
New Models and Hydrocodes, Lisbon, Portugal,19-23 May 2008 Applying the Staggered Mesh Godunov (SMG) Method to Reactive Flows Gabi Luttwak 1 and Joseph Falcovitz 2 1 Rafael, P.O. Box 2250, Haifa 31021,
More informationAA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow
AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
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 informationEVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY
EVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY Borzou, B. 1 and Radulescu, M.I. 2 1 Mechanical Engineering Department, University of Ottawa, 161 Louis Pasteur,
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 informationThe evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave
Combustion and Flame 151 (2007) 573 580 www.elsevier.com/locate/combustflame The evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave Barbara B. Botros
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
More informationAnalysis for Steady Propagation of a Generic Ram Accelerator/ Oblique Detonation Wave Engine Configuration
Analysis for Steady Propagation of a Generic Ram Accelerator/ Oblique Detonation Wave Engine Configuration J.M. Powers1, D.R. Fulton2, K.A. Gonthier3, and M.J. Grismer4 Department of Aerospace and Mechanical
More informationReaction Rate Closure for Turbulent Detonation Propagation through CLEM-LES
5 th ICDERS August 7, 05 Leeds, UK through CLEM-LES Brian Maxwell, Matei Radulescu Department of Mechanical Engineering, University of Ottawa 6 Louis Pasteur, Ottawa, KN 6N5, Canada Sam Falle School of
More informationOn thermodynamic cycles for detonation engines
On thermodynamic cycles for detonation engines R. Vutthivithayarak, E.M. Braun, and F.K. Lu 1 Introduction Detonation engines are considered to potentially yield better performance than existing turbo-engines
More informationAstrophysical Combustion: From a Laboratory Flame to a Thermonuclear Supernova
25 th ICDERS August 2 7, 2015 Leeds, UK : From a Laboratory Flame to a Thermonuclear Supernova Alexei Y. Poludnenko Naval Research Laboratory Washington, D.C., USA 1 Introduction Exothermic processes associated
More informationCombustion Behind Shock Waves
Paper 3F-29 Fall 23 Western States Section/Combustion Institute 1 Abstract Combustion Behind Shock Waves Sandeep Singh, Daniel Lieberman, and Joseph E. Shepherd 1 Graduate Aeronautical Laboratories, California
More informationNumerical Simulations of a Stratified Oceanic Bottom Boundary Layer. John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD
Numerical Simulations of a Stratified Oceanic Bottom Boundary Layer John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD Motivation Objective I: Assess and improve parameterizations of the bottom boundary
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationEffects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube
5 th ICDERS August 7, 015 Leeds, UK Effects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube Sung Hwan Yoon and Osamu Fujita Division
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationInfluence of the Reaction Heat on the Linear Stability Spectra of Steady Detonation in the Kinetic Frame
Influence of the Reaction Heat on the Linear Stability Spectra of Steady Detonation in the Kinetic Frame Filipe Carvalho and Ana Jacinta Soares Departamento de Ciências Eactas, Viana do Castelo, Portugal
More informationPeriodic oscillation and fine structure of wedge-induced oblique detonation waves
Acta Mech. Sin. (211) 27(6):922 928 DOI 1.17/s149-11-58-y RESEARCH PAPER Periodic oscillation and fine structure of wedge-induced oblique detonation waves Ming-Yue Gui Bao-Chun Fan Gang Dong Received:
More informationUNIVERSITY OF CALGARY. Stability of Detonation Wave under Chain Branching Kinetics: Role of the Initiation Step. Michael Anthony Levy A THESIS
UNIVERSITY OF CALGARY Stability of Detonation Wave under Chain Branching Kinetics: Role of the Initiation Step by Michael Anthony Levy A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT
More informationPeriod Doubling Cascade in Diffusion Flames
Period Doubling Cascade in Diffusion Flames Milan Miklavčič Department of Mathematics, Michigan State University, East Lansing, MI 48824, USA Combustion Theory and Modelling 11 No 1 (2007), 103-112 Abstract
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationarxiv: v1 [physics.flu-dyn] 14 Jul 2018
A minimal hyperbolic system for unstable shock waves Dmitry I. Kabanov a, Aslan R. Kasimov b,c, a Division of Computer, Electrical and Mathematical Sciences and Engineering King Abdullah University of
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationOn Analogy of 2D and 3D Combustible Mixture Flows
25 th ICDERS August 2 7, 2015 Leeds, UK On Analogy of 2D and 3D Combustible Mixture Flows Vladimir A. Levin 1, Ivan S. Manuylovich 1, Vladimir V. Markov 1 1 Institute of mechanics of the Moscow State University,
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Central concepts: Phase velocity: velocity with which surfaces of constant phase move Group velocity: velocity with which slow
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationUnsteady Propagation Process of Oblique Detonation Waves Initiated by Hypersonic Spherical Projectiles
Trans. JSASS Aerospace Tech. Japan Vol. 10, No. ists28, pp. Pe_1-Pe_6, 2012 Original Paper Unsteady Propagation Process of Oblique Detonation Waves Initiated by Hypersonic Spherical Projectiles By Shinichi
More informationDetonation of Gas Particle Flow
2 Detonation of Gas Particle Flow F. Zhang 2.1 Introduction Fine organic or metallic particles suspended in an oxidizing or combustible gas form a reactive particle gas mixture. Explosion pressures in
More informationAsymptotic theory of evolution and failure of self-sustained detonations
J. Fluid Mech. (5), vol. 55, pp. 161 19. c 5 Cambridge University Press DOI: 1.117/S114599 Printed in the United Kingdom 161 Asymptotic theory of evolution and failure of self-sustained detonations By
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 informationSolving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation
Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as
More informationChemical inhibiting of hydrogen-air detonations Sergey M. Frolov
Chemical inhibiting of hydrogen-air detonations Sergey M. Frolov Semenov Institute of Chemical Physics Moscow, Russia Outline Introduction Theoretical studies Classical 1D approach (ZND-model) Detailed
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationEQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE
EQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE D. Scott Stewart and Sunhee Yoo Department of Theoretical and Applied Mechanics University of Illinois, Urbana, IL 61801, USA William C. Davis
More informationDirect Simulation of Ultrafast Detonations in Mixtures
Direct Simulation of Ultrafast Detonations in Mixtures Patrick D. O Connor *, Lyle N. Long * and James B. Anderson * Department of Aerospace Engineering, The Pennsylvania State University, University Park,
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 informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationMild Ignition Phenomena in Rapid Compression Machines
25 th ICDERS August 2 7, 2015 Leeds, UK Kevin P. Grogan a, S. Scott Goldsborough b, Matthias Ihme a a Stanford University, Stanford, CA 94305 b Argonne National Laboratory, Argonne, IL 60439 1 Introduction
More informationPhysics of Explosions
Physics of Explosions Instructor: Dr. Henry Tan Pariser/B4 MACE - Explosion Engineering School of Mechanical, Aerospace and Civil Engineering The University of Manchester Introduction Equation of state
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 informationFundamentals of Rotating Detonation. Toshi Fujiwara (Nagoya University)
Fundamentals of Rotating Detonation Toshi Fujiwara (Nagoya University) New experimental results Cylindical channel D=140/150mm Hydrogen air; p o =1.0bar Professor Piotr Wolanski P [bar] 10 9 8 7 6 5 4
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
More informationEffect of Initial Disturbance on The Detonation Front Structure of a Narrow Duct. Hua-Shu Dou *, Boo Cheong Khoo,
Published in: Shock Waves, Vol.20, No.2, March 2010, pp.163-173. Effect of Initial Disturbance on The Detonation Front Structure of a Narrow Duct Hua-Shu Dou *, Boo Cheong Khoo, Department of Mechanical
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
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 informationBIFURCATING MACH SHOCK REFLECTIONS WITH APPLICATION TO DETONATION STRUCTURE
BIFURCATING MACH SHOCK REFLECTIONS WITH APPLICATION TO DETONATION STRUCTURE Philip Mach A thesis submitted to the Faculty of Graduate and Postdoctoral Studies in partial fulfillment of the requirements
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationAn artificial nonlinear diffusivity method for shock-capturing in supersonic reacting flows
Center for Turbulence Research Annual Research Briefs 005 57 An artificial nonlinear diffusivity method for shock-capturing in supersonic reacting flows By B. Fiorina AND S. K. Lele. Motivation and objectives
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationCALCULATION OF PRESSURE FIELD IN THE PROBLEM OF SONIC BOOM FROM VARIOUS THIN AXISYMMETRIC BODIES
CALCULATION OF PRESSURE FIELD IN THE PROBLEM OF SONIC BOOM FROM VARIOUS THIN AXISYMMETRIC BODIES А.V. Potapkin, D.Yu. Moskvichev Khristianovich Institute of Theoretical and Applied Mechanics SB RAS 630090,
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More 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 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 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 informationFormation and Evolution of Distorted Tulip Flames
25 th ICDERS August 2 7, 2015 Leeds, UK Formation and Evolution of Distorted Tulip Flames Huahua Xiao 1, 2, Ryan W. Houim 1, Elaine S. Oran 1, and Jinhua Sun 2 1 University of Maryland College Park, Maryland,
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 informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More informationThis paper analyzes the role of mixture reactivity on the dynamics of the Richtmyer-Meshkov instability
39th AIAA Fluid Dynamics Conference - 5 June 9, San Antonio, Texas AIAA 9-357 39th AIAA Fluid Dynamics Conference - 5 Jun 9, Grand Hyatt Hotel, San Antonio, Texas Richtmyer-Meshkov instability in reactive
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationSteady Deflagration Structure in Two-Phase Granular Propellants
12th ICDERS, Ann Arbor, Michigan July 23-28, 1989 Steady Deflagration Structure in Two-Phase Granular Propellants Joseph M. Powers 1, Mark E. Miller2, D. Scott Stewart3, and Herman Krier4 July 23-28, 1989
More informationTurbulent spectra generated by singularities
Turbulent spectra generated by singularities p. Turbulent spectra generated by singularities E.A. Kuznetsov Landau Institute for Theoretical Physics, Moscow Russia In collaboration with: V. Naulin A.H.
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationNumerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry
Numerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry V. S. Ajaev and J. Tausch 2 Southern Methodist University ajaev@smu.edu 2 Southern Methodist University tausch@smu.edu
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationResonant solitons in a polydisperse bubble medium
Journal of Physics: Conference Series PAPER OPEN ACCESS Resonant solitons in a polydisperse bubble medium To cite this article: I A Ogorodnikov 018 J. Phys.: Conf. Ser. 1105 01088 View the article online
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationStructural model for the dynamic buckling of a column under constant rate compression
Structural model for the dynamic buckling of a column under constant rate compression Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg Polytechnical University arxiv:1506.0047v
More informationLecture 7 Detonation Waves
Lecture 7 etonation Waves p strong detonation weak detonation weak deflagration strong deflagration / 0 v =/ University of Illinois at Urbana- Champaign eflagrations produce heat Thermal di usivity th
More informationInternational Conference on Methods of Aerophysical Research, ICMAR 2008
International Conference on Methods of Aerophysical Research, ICMAR 8 EXPERIMENTAL STUDY OF UNSTEADY EFFECTS IN SHOCK WAVE / TURBULENT BOUNDARY LAYER INTERACTION P.A. Polivanov, А.А. Sidorenko, A.A. Maslov
More informationAnalytical and Numerical Investigation of Transient Gas Blow down
Analytical and Numerical Investigation of Transient Gas Blow down Hameed Balassim Mahood, Farhan Lafta Rasheed & Abdulsattar K. Abbas Ministry of Sciences and Technology, Baghdad, Iraq Received: August
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 informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationINVESTIGATION OF INSTABILITIES AFFECTING DETONATIONS: IMPROVING THE RESOLUTION USING BLOCK-STRUCTURED ADAPTIVE MESH REFINEMENT PRASHAANTH RAVINDRAN
INVESTIGATION OF INSTABILITIES AFFECTING DETONATIONS: IMPROVING THE RESOLUTION USING BLOCK-STRUCTURED ADAPTIVE MESH REFINEMENT by PRASHAANTH RAVINDRAN Presented to the Faculty of the Graduate School of
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationA second order scheme on staggered grids for detonation in gases
A second order scheme on staggered grids for detonation in gases Chady Zaza IRSN PSN-RES/SA2I/LIE August 1, 2016 HYP 2016, RWTH Aachen Joint work with Raphaèle Herbin, Jean-Claude Latché and Nicolas Therme
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationAdvanced Multi-Scale Methods for Hypersonic Propulsion
Advanced Multi-Scale Methods for Hypersonic Propulsion By J. M. Powers, D. Wirasaet and S. Paolucci Aerospace and Mechanical Engineering University of Notre Dame, Indiana 46556 NASA Hypersonics NRA Review,
More information