Development of an evaporation boundary condition for DSMC method with application to meteoroid entry
|
|
- Domenic Underwood
- 5 years ago
- Views:
Transcription
1 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
2 Physics of a meteoroid entering the atmosphere Coupled physico-chemical phenomena occurring during meteoroid ablation Velocity: from 12 km/s to 72 km/s Size: from micro size grains to meters 2 / 26
3 Focus on the gas-surface phenomena VKI Plasmatron experiment: basaltic sample, 3 MW/m2 heat flux at the stagnation point Schematic of gas-surface interactions 2 Ordinary chondritic sample, 1 MW/m heat flux at the stagnation point 3 / 26
4 Objectives Goals Asses the importance in the competition among the different phenomena Develop an evaporation boundary condition capable of dealing with magma compositions (mixture of oxides) Couple the flow with the material response Introduce a method to compute chemistry of ablated species 4 / 26
5 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work 4 / 26
6 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Governing equations and numerical method 4 / 26
7 Breakdown of the continuum regime Knudsen number = λ D = mesoscopic scale macroscopic scale Governing equations and numerical method Boltzmann equation 5 / 26
8 Boltzmann equation f(x, ξ, t) = one-particle velocity distribution function f t + ξ f X + F f ξ = { }}{ 4π ( f f 1 ff ) 1 gσdωdξ1 0 From microscopic to macroscopic world: ρ = ρv = ρe = f(x, ξ, t) dξ ξ f(x, ξ, t) dξ 1 ξ ξ f(x, ξ, t) dξ 2 q = ϕ = f(x, ξ, t) [ξ v w] n dξ 1 ξ ξ f(x, ξ, t) [ξ vw] n dξ 2 Governing equations and numerical method Boltzmann equation 6 / 26
9 Direct Simulation Monte Carlo method DSMC algorithm Each simulated particle represents a large number of real particles Free motion decoupled from collisions (for t < ν 1 coll ) Grid cells used to choose collisions partners and sample averages DSMC is not Molecular Dynamics Governing equations and numerical method Direct Simulation Monte Carlo method 7 / 26
10 SPARTA numerical tool Developed by Plimpton and Gallis at Sandia National Labs Open source software ( Object-oriented philosophy enables extensions Parallel implementation through domain decomposition Flow around MIR space station Apollo re-entering the atmosphere Governing equations and numerical method SPARTA numerical tool 8 / 26
11 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Ablative wall 8 / 26
12 Kinetic boundary conditions for Boltzmann equation [ξ v w] n<0 flux emerging from the wall {}}{ flux due to evaporation by the wall {}}{ f(ξ)[ξ v w ] n = g(ξ)[ξ v w ] n + K B (ξ ξ)f(ξ )[ξ v w ] n dξ } {{ } flux due to reflection by the wall K B (ξ ξ) = (1 α c ) [ξ v w] n 2π(RT w ) 2 g(ξ) = α e /α c evaporation/condensation coefficients (usually α e = α c = 1) ρ eq w equilibrium vapor density No need to assume equilibrium of the gas with the wall [ξ v w ] n > 0 ( exp ξ v w 2 ) 2RT w α e ρ eq ( w (2πRT w ) exp ξ v w 2 ) 3/2 2RT w Ablative wall Kinetic boundary conditions 9 / 26
13 MAGMA chemical multi-phase equilibrium solver How do we obtain the equilibrium properties? = ρ eq w i (material composition, T w ) ρ eq w i Developed by Fegley and Cameron, 1987 Mass balance, mass action algorithm Only stoichiometric reactions for change of phase Extensively validated vs. experimental data Already used to model vaporization in silicate lavas (under thermodynamic equilibrium assumption) Ablative wall Melt-vapor equilibrium 10 / 26
14 Extension to mixtures of oxides How do we choose the coefficients? α ei, α ci i {O, Si, SiO 2, Mg, MgO,...} Hp. 0: M x O y (l) x M (g) + y O (g) α ci = 0 i M, O Hp. 1: α cm = α em = α M = α Hp. 2: α co = α eo = α O Thermodynamic equilibrium has to be retrieved by kinetic approach ϕ eq i ϕ eq e M ϕ eq e O y k ϕ eq e Mk = ϕ eq i (α i, ρ eq w i, T w ) = ϕ eq c M M = ϕ eq c O k R ϕeq e Ok = x k ϕ eq e Ok k R = ϕ eq c O Ablative wall Extension to mixtures of oxides 11 / 26
15 Electron concentration around a non-ablating meteoroid Kn = 0.1 D = 1 cm H = 80 km V = 72 km/s non-ablating meteoroid m 3 typical e density for meteors Ambipolar diffusion assumption Gas phase chemistry frozen above 90 km electrons from metal species Ablative wall Extension to mixtures of oxides 12 / 26
16 Ablating meteoroid at 95 km altitude Pure magnesium: T wall = 925 K (melting temperature) v wall = 0 m/s α c = α e = 1 D = 1 cm H = 95 km V = 72 km/s Cooling effect of evaporation at the wall Adiabatic T 2.5M K >> T wall Geometric temperature: particles not characterized by this temperature in the thermal sense Translational temperature Ablative wall Extension to mixtures of oxides 13 / 26
17 Ablating meteoroid at 95 km altitude Delay in the excitation of the internal dofs Strong shielding effects of Mg vapor Rotational temperature Ablation products molar fraction Ablative wall Extension to mixtures of oxides 14 / 26
18 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Material response 14 / 26
19 Material response T Randomly and rapidly rotating sphere: t = k 2 T ρc p r + 2k 2 Ablating wall: moving mesh (fixed reference frame) Re-mapping procedure at each time step Finite differences, explicit time integration 1 ρc p r T r Derived variable [K m] Time [s]: Numerical solution Temperature [K] Time [s]: Numerical solution Distance from the surface of the sphere [m] Verification of the spherical coordinates: unsteady solution Distance from the initial surface position [m] Verification of the moving wall (re-mapping): steady solution Material response Material code 15 / 26
20 Tracking of the melting front CP Liquid (C S P + C L P )/2 + L/(2 T) Solid L C P S C P Position of the liquid-solid interface [m] Latent heat of melting [J/Kg]: 4e+05 4e+06 2e+07 Numerical solution 2 T Heat capacity T Time [s] Verification of the apparent heat capacity method: position of the melting front Apparent heat capacity method No need to deform the mesh to track the solid-liquid interface Position of the melting front obtained a posteriori Material response Apparent heat capacity method 16 / 26
21 Interface flow-material ϕ e ϕ c = ρ w v w ϕ e computed theoretically ϕ c directly from DSMC simulation No contribution of reflection (no net flux) Surface mass balance q c + q r q e + ɛσt 4 = ρ w h w v w + k T r + ɛσtw 4 w ϕ e computed theoretically q c, q r directly from DSMC simulation Surface energy balance Material response Interface flow-material 17 / 26
22 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Trajectory code 17 / 26
23 Trajectory code Python implementation (cython to improve performances) Interface with MAGMA for wall equilibrium properties Interface with NASA atmospheric model for free-stream properties Sub time stepping for material response Flow update resolution can be fixed through input file Trajectory code 18 / 26
24 Trajectory-material response coupling No DSMC coupling: evaporation into vacuum (Knudsen-Langmuir) q flow = Λρ (H)V 3 8 Temperature [K] V = 15 km/s, D = m Core temperature Surface temperature Position from the core [m] 10-3 V = 15 km/s, D = m 1 Evaporating front 0.9 Melting front Molten thickness Altitude [km] Altitude [km] D = 2 mm V = 15 km/s Trajectory code 19 / 26
25 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Gas phase chemistry 19 / 26
26 Lagrangian solver LARSEN LAgrangian Reactor for StrEams in Nonequilibrium Atmospheric entry ows many chemical products (air, ablated species) thermal non-equilibrium (T, Tr, Tv) 2D / 3D geometries TOO EXPENSIVE! Result: only simple models currently employed HOWEVER... METEOROID TRAIL: detailled ionization mechanisms required... IDEA: introduce chemistry a posteriori! Gas phase chemistry Lagrangian solver 20 / 26
27 Chemistry of ablated species computed a posteriori 1) Simple simulation 2) Extract streamlines Y i 4) Refined results u(s) ρ(s) + ICs 3) Lagrangian solver s Velocity and density fields from baseline simulation assumed good enough Gas phase chemistry Lagrangian solver 21 / 26
28 LARSEN formulation From baseline simulation: u, ρ are given We still need: Species mass eq. s y i = ω i J i ρu Total enthalpy eq. s H = Q ext Species internal energy s e in i = Din i + Ω in i h in i ω i ρy i u J i ω i Q u = System of ODEs Gas phase chemistry Lagrangian solver 22 / 26
29 Hypersonic inert flow around a cylinder (b) (b) (b) Axisymmetric DSMC simulation Pure argon M = 10 Kn = 0.05 T correctly reproduced Temperature [K] (a) (b) (a) (a) (a) DSMC SPA LARSENLAR curvilinear abscissa [m] Gas phase chemistry Lagrangian solver 23 / 26
30 Relaxation behind a shock wave - refining chemistry Inviscid flow Fire-II capsule free-stream conditions air5: N 2 O 2 N O NO air11: air5 + N + O + NO + N + 2 O+ 2 e Temperatu Temperature [K] Mass fractions N 2 T well predicted ions concentrations well predicted neutral species concentrations improved distance from shock [m] distance from shock [m] O 2 NO SHOCKING - air5 LARSEN - air5 to air11 SHOCKING - air11 N Gas phase chemistry Lagrangian solver 24 / 26 N O O +
31 1 Governing equations and numerical method Boltzmann equation Direct Simulation Monte Carlo method SPARTA numerical tool 2 Ablative wall Kinetic boundary conditions Melt-vapor equilibrium Extension to mixtures of oxides 3 Material response Material code Apparent heat capacity method Interface flow-material 4 Trajectory code 5 Gas phase chemistry Lagrangian solver 6 Conclusions and future work Conclusions and future work 24 / 26
32 Conclusions and future work Conclusions Develop a DSMC evaporation boundary condition for silicates Couple flow-material-trajectory Develop a method to implement detailed chemistry a posteriori Future work Include dynamic of molten layer to gas surface interactions Add ionization of metallic species Assess recombination time of free electrons in the trail Conclusions and future work 25 / 26
33 Thanks for your attention Thanks to: B. Dias for the useful discussion G. Bellas for the help with LARSEN L. Zavalan, B. Helber, P. Collins for the experiments Conclusions and future work 26 / 26
Lagrangian diffusive reactor for detailed thermochemical computations of plasma flows
Lagrangian diffusive reactor for detailed thermochemical computations of plasma flows Stefano Boccelli 1,2, a), Federico Bariselli 1,2,3, b), Bruno Dias 1,4, c) 1, d) and Thierry E. Magin 1 von Karman
More informationApplication of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows
Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows Timothy R. Deschenes and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor,
More informationAll-Particle Multiscale Computation of Hypersonic Rarefied Flow
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 2010, Orlando, Florida AIAA 2010-822 All-Particle Multiscale Computation of Hypersonic Rarefied
More informationModeling of the Electric Field in a Hypersonic Rarefied Flow
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 2010, Orlando, Florida AIAA 2010-635 Modeling of the Electric Field in a Hypersonic Rarefied Flow
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 INVESTIGATION OF AMPLITUDE DEPENDENCE ON NONLINEAR ACOUSTICS USING THE DIRECT SIMULATION MONTE CARLO METHOD PACS: 43.5.Ed Hanford, Amanda
More informationAtmospheric Entry. Technology, Mathematical Model and Simulation
Atmospheric Entry Technology, Mathematical Model and Simulation Julian Köllermeier RWTH Aachen, August 26th 2016 Outline 1. Introduction to Atmospheric Reentry 2. Rarefied Gases: From Science Fiction to
More informationTAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF
TAU Extensions for High Enthalpy Flows Sebastian Karl AS-RF Contents Motivation Extensions available in the current release: Numerical schemes for super- and hypersonic flow fields Models for gas mixtures,
More informationCALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE
J. Comput. Fluids Eng. Vol.17, No.4, pp.69-74, 2012. 12 / 69 CALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE IN NONEQUILIBRIUM HYPERSONIC FLOW M. Ahn Furudate * Dept. of Mechatronics Engineering,
More informationSimulation of the Interaction Between Two Counterflowing Rarefied Jets
Simulation of the Interaction Between Two Counterflowing Rarefied Jets Cyril Galitzine and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 Abstract. A preliminary
More informationVibrational degrees of freedom in the Total Collision Energy DSMC chemistry model
Vibrational degrees of freedom in the Total Collision Energy DSMC chemistry model Mark Goldsworthy, Michael Macrossan Centre for Hypersonics, School of Engineering, University of Queensland, Brisbane,
More informationVERIFICATION AND VALIDATION OF A PARALLEL 3D DIRECT SIMULATION MONTE CARLO SOLVER FOR ATMOSPHERIC ENTRY APPLICATIONS
VERIFICATION AND VALIDATION OF A PARALLEL 3D DIRECT SIMULATION MONTE CARLO SOLVER FOR ATMOSPHERIC ENTRY APPLICATIONS Paul Nizenkov and Stefanos Fasoulas Institut für Raumfahrtsysteme (IRS), Universität
More informationApplication of SMILE++ Computational Tool to High-enthalpy Ionized Flows
Application of SMILE++ Computational Tool to High-enthalpy Ionized Flows A. Shevyrin, P. Vashchenkov, A. Kashkovsky, Ye. Bondar Khristianovich Institute of Theoretical and Applied Mechanics Novosibirsk
More informationHypersonic Blunt Body Thermophysics Using a Unified Kinetic/Continuum Solver
41st AIAA Thermophysics Conference 22-25 June 2009, San Antonio, Texas AIAA 2009-3838 Hypersonic Blunt Body Thermophysics Using a Unified Kinetic/Continuum Solver Andrew J. Lofthouse U.S. Air Force Institute
More informationKinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods
Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering
More informationProgress Report on Chamber Dynamics and Clearing
Progress Report on Chamber Dynamics and Clearing Farrokh Najmabadi, Rene Raffray, Mark S. Tillack, John Pulsifer, Zoran Dragovlovic (UCSD) Ahmed Hassanein (ANL) Laser-IFE Program Workshop May31-June 1,
More informationParticle Simulation of Nonequilibrium Hypersonic Flows
Particle Simulation of Nonequilibrium Hypersonic Flows Thomas E. Schwartzentruber Assistant Professor Department of Aerospace Engineering and Mechanics (AEM) University of Minnesota [1 /36] Molecular Simulation
More informationKinetic Effects in Spherical Expanding Flows of Binary-Gas Mixtures
Kinetic Effects in Spherical Expanding Flows of Binary-Gas Mixtures Vladimir V. Riabov Rivier College, Nashua, New Hampshire, USA Abstract. Diffusion effects in the spherical expanding flows of argon-helium
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationDSMC Simulation of Entry Vehicle Flowfields Using a Collision-Based Chemical Kinetics Approach
DSMC Simulation of Entry Vehicle Flowfields Using a Collision-Based Chemical Kinetics Approach R.G. Wilmoth a, D.B. VanGilder a, and J.L. Papp a a Combustion Research and Flow Technology, Inc., 6210 Keller
More informationA hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models
A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models Alessandro Alaia, Gabriella Puppo May 31, 2011 Abstract In this work we present a non stationary domain
More informationHigh Altitude Rocket Plume and Thermal Radiation Analysis
High Altitude Rocket Plume and Thermal Radiation Analysis [ Woo Jin Jeon, Seung Wook Baek, Jae Hyun Park and Dong Sung Ha ] Abstract In this study, rocket plume behavior at various altitudes and radiative
More informationHYPERSONIC vehicles generate multi-scale gas flows as a result of both large variations in mean-free-path and
39th AIAA Thermophysics Conference 25-28 June 27, Miami, FL AIAA 27-3892 Hybrid Particle-Continuum Simulations of Low Knudsen Number Hypersonic Flows Thomas E. Schwartzentruber, Leonardo C. Scalabrin and
More informationLECTURE 5 SUMMARY OF KEY IDEAS
LECTURE 5 SUMMARY OF KEY IDEAS Etching is a processing step following lithography: it transfers a circuit image from the photoresist to materials form which devices are made or to hard masking or sacrificial
More informationLecture 6 Gas Kinetic Theory and Boltzmann Equation
GIAN Course on Rarefied & Microscale Gases and Viscoelastic Fluids: a Unified Framework Lecture 6 Gas Kinetic Theory and Boltzmann Equation Feb. 23 rd ~ March 2 nd, 2017 R. S. Myong Gyeongsang National
More informationTarget Thermal Response to Gas Interactions
University of California, San Diego UCSD-ENG-092 Target Thermal Response to Gas Interactions A. R. Raffray, J. Pulsifer and M. S. Tillack June 24, 2002 Fusion Division Center for Energy Research University
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More information(2015) 53 (6) ISSN
Scanlon, Thomas J. and White, Craig and Borg, Matthew K. and Palharini, Rodrigo C. and Farbar, Erin and Boyd, Iain D. and Reese, Jason M. and Brown, Richard E. (2015) Open source Direct Simulation Monte
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationEntry Modeling for Asteroid Threat Assessment
Presenter: Eric Stern Entry Modeling for Asteroid Threat Assessment Asteroid Threat Assessment Project (ATAP), Entry Modeling Lead Team: Parul Agrawal, Dinesh Prabhu, Francesco Panerai, Justin Haskins,
More informationScaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09
Scaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09 Michael Macrossan, Centre for Hypersonics, University of Queensland
More informationFollow this and additional works at:
Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 12-19-2016 The Effects
More informationProblems of Flows Through Short Channels Studied by Means of the Boltzmann Equation
Problems of Flows Through Short Channels Studied by Means of the Boltzmann Equation Aristov V.V., Frolova A.A. and Zabelok S.A. Dorodnicyn Computing Centre of Russian Academy of Sciences Kolobov V.I. and
More informationDirect Monte Carlo Simulations of a Binary Gas Free-Jet Flow Over a Flat Plate
Direct Monte Carlo Simulations of a Binary Gas Free-Jet Flow Over a Flat Plate Eric Cecil Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA 2293, USA Abstract.
More information8th European Symposium on ATD for Space Vehicle, Lisbon, 2-6 March 2015
1 FAST a Pre-Design Tool for the Atmospheric Re-entry of Space Vehicles and Debris F. Sourgen (*) Y. Prévereaud J-L. Vérant E. Laroche J-M. Moschetta (*) frederic.sourgen@onera.fr Onera Midi Pyrénées Center
More informationTOPICAL PROBLEMS OF FLUID MECHANICS PERFORMANCE OF SIMPLE CONDENSATION MODEL IN HIGH-PRESSURES
TOPICAL PROBLEMS OF FLUID MECHANICS 9 DOI: http://dx.doi.org/.43/tpfm.6.9 PERFORMANCE OF SIMPLE CONDENSATION MODEL IN HIGH-PRESSURES V. Hric, J. Halama Department of Technical Mathematics, Faculty of Mechanical
More informationSimulation of a typical reentry vehicle TPS local flow features and material response
Simulation of a typical reentry vehicle TPS local flow features and material response E. V. Titov, and D. A. Levin Pennsylvania State University, University Park, PA 1682 Abstract. Statistical BGK/DSMC
More informationPredicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions
Predicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 Abstract. The breakdown of the
More informationA First Course on Kinetics and Reaction Engineering Unit 2. Reaction Thermochemistry
Unit 2. Reaction Thermochemistry Overview This course is divided into four parts, I through IV. Part I reviews some topics related to chemical reactions that most students will have encountered in previous
More informationA Zooming Approach to Investigate Heat Transfer in Liquid Rocket Engines with ESPSS Propulsion Simulation Tool
A Zooming Approach to Investigate Heat Transfer in Liquid Rocket Engines with ESPSS Propulsion Simulation Tool M. Leonardi, F. Di Matteo, J. Steelant, B. Betti, M. Pizzarelli, F. Nasuti, M. Onofri 8th
More informationHypersonic Flow of Rarefied Gas Near the Brazilian Satellite During its Reentry into Atmosphere
398 Brazilian Journal of Physics, vol. 33, no. 2, June, 2003 Hypersonic Flow of Rarefied Gas Near the Brazilian Satellite During its Reentry into Atmosphere Felix Sharipov Departamento de Física, Universidade
More informationBy ablation we mean the recession of a surface due to heating, usually by a hot gas. It is the key process for
16.50 Lecture 15 Subject: Ablative cooling By ablation we mean the recession of a surface due to heating, usually by a hot gas. It is the key process for a) Re-entry heat shields b) Solid propellant nozzles
More informationSimulation of Condensing Compressible Flows
Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,
More informationModeling of Plasma Formation in Rarefied Hypersonic Entry Flows
45th AIAA Aerospace Sciences Meeting and Exhibit 8-11 January 2007, Reno, Nevada AIAA 2007-206 Modeling of Plasma Formation in Rarefied Hypersonic Entry Flows Iain D. Boyd University of Michigan, Ann Arbor,
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationNumerical studies of hypersonic binary gas-mixture flows near a sphere
Numerical studies of hypersonic binary gas-mixture flows near a sphere V.V. Riabov Diffusive Effects in Binary Gas-Mixture Flows near a Sphere Diffusion processes have a significant effect on the structure
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationNon-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture
Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture E.V. Kustova, E.A. Nagnibeda, Yu.D. Shevelev and N.G. Syzranova Department of Mathematics and Mechanics,
More informationValidation of the Chemistry Module for the Euler Solver in Unified Flow Solver
Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-22-2012 Validation of the Chemistry Module for the Euler Solver in Unified Flow Solver William C. Humphrey Jr. Follow this and
More informationA Hybrid CFD-DSMC Method of Modeling Continuum-Rarefied Flows
42nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 24, Reno, Nevada AIAA 24-8 42nd AIAA Aerospace Sciences Meeting and Exhibit Reno, NV January 5 8, 24 A - Method of Modeling Continuum-Rarefied
More informationANALYSIS OF HEAT TRANSFER IN HYPERSONIC FLOW OVER RE-ENTRY CONFIGURATIONS
IJATSE Volume 3 Number 1 January-June 2012, pp. 1-10 ANALYSIS OF HEAT TRANSFER IN HYPERSONIC FLOW OVER RE-ENTRY CONFIGURATIONS R. Balu 1 and L. Prince Raj *2 ABSTRACT: Inter planetary space mission involves
More informationA New Energy Flux Model in the DSMC-IP Method for Nonequilibrium Flows
36th AIAA Thermophysics Conference 23-26 June 23, Orlando, Florida AIAA 23-3774 A New Energy Flux Model in the -IP Method for Nonequilibrium Flows Wen-Lan Wang and Iain D. Boyd Department of Aerospace
More informationPlanar Laser-Induced Iodine Fluorescence Technique for Flow Visualization and Quantitative Measurements in Rarefied Flows
Planar Laser-Induced Iodine Fluorescence Technique for Flow Visualization and Quantitative Measurements in Rarefied Flows Professor James McDaniel*, Eric Cecil*, Erin Reed* and Josh Codoni* Professor Iain
More informationDevelopment of a General Purpose DSMC Flow Solver on Unstructured Meshes. Code
36th AIAA Thermophysics Conference 23-26 June 2003, Orlando, Florida AIAA 2003-3776 Development of a General Purpose DSMC Flow Solver on Unstructured Meshes V. Tenishev,M.Combi Space Physics Research Laboratory,
More informationCFD ANALYSIS OF THE CORK-PHENOLIC HEAT SHIELD OF A REENTRY CUBESAT IN ARC-JET CONDITIONS INCLUDING ABLATION AND PYROLYSIS
CFD ANALYSIS OF THE CORK-PHENOLIC HEAT SHIELD OF A REENTRY CUBESAT IN ARC-JET CONDITIONS INCLUDING ABLATION AND PYROLYSIS 15 th International Planetary Probe Workshop 14/06/2018 Ata Onur Başkaya 1 A. Turchi
More informationAnalysis of Bridging Formulae in Transitional Regime
Analysis of Bridging Formulae in Transitional Regime Luigi Morsa *, Gennaro Zuppardi *, Antonio Schettino ** and Raffaele Votta ** * Department of Aerospace Engineering University of Naples Federico II,
More informationWUFI Workshop at NTNU /SINTEF Fundamentals
WUFI Workshop at NTNU /SINTEF 2008 Fundamentals Contents: From steady-state to transient Heat storage and -transport Moisture storage and -transport Calculation of coupled transport Model limitations 2
More informationA First Course on Kinetics and Reaction Engineering Unit 19. Analysis of Batch Reactors
Unit 19. Analysis of Batch Reactors Overview As noted in Unit 18, batch reactor processing often follows an operational protocol that involves sequential steps much like a cooking recipe. In general, each
More informationEffective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body
Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body I.G. Brykina a, B.V. Rogov b, I.L. Semenov c, and G.A. Tirskiy a a Institute of Mechanics, Moscow
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationNumerical Prediction of Nonequilibrium Hypersonic Flow Around Brazilian Satellite SARA
48 Brazilian Journal of Physics, vol. 35, no., March, 5 Numerical Prediction of Nonequilibrium Hypersonic Flow Around Brazilian Satellite SARA Ghislain Tchuen, Institut Universitaire de Technologie Fotso
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationOn the Coefficients in Meteor Physics Equations
On the Coefficients in Meteor Physics Equations aria Yu. Khanukaeva epartment of Applied Mathematics, Moscow Institute of Physics and Technology (State University), Institutsky 9, 141700 olgoprudny, Moscow
More informationCompressible Flow - TME085
Compressible Flow - TME085 Lecture 13 Niklas Andersson Chalmers University of Technology Department of Mechanics and Maritime Sciences Division of Fluid Mechanics Gothenburg, Sweden niklas.andersson@chalmers.se
More informationDIRECT SIMULATION CALCULATIONS OF THE RAREFIED HYPERSONIC FLOW PAST A BACKWARD-FACING STEP
9 Brazilian Symposium on Aerospace Eng. & Applications Copyright c 9 by AAB September 14-16, 9, S. J. Campos, SP, Brazil DIRECT SIMULATION CALCULATIONS OF THE RAREFIED HPERSONIC FLOW PAST A BACKWARD-FACING
More informationChemical kinetic and radiating species studies of Titan aerocapture entry
16th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Chemical kinetic and radiating species studies of Titan aerocapture entry Pénélope Leyland 1, Raffaello
More informationarxiv: v1 [physics.comp-ph] 7 Mar 2018
A BGK model for high temperature rarefied gas flows C. Baranger 1, G. Marois 1, J. Mathé 1, J. Mathiaud 1, L. Mieussens arxiv:1803.0617v1 [physics.comp-ph] 7 Mar 018 1 CEA-CESTA 15 avenue des Sablières
More informationABSTRACT. Nomenclature
ABSTRACT The behavior of two different models of gas-surface interactions is studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic
More informationHypersonic Flow over a Wedge with a Particle Flux Method
Hypersonic Flow over a Wedge with a Particle Flux Method M. N. Macrossan, M. V. Metchnik and P. A. Pinto Centre for Hypersonics, School of Engineering, University of Queensland, St Lucia, 47, Australia
More informationOn Notation Thermodynamics of Glaciers. Types of Glaciers. Why we care. McCarthy Summer School
On Notation Thermodynamics of Glaciers McCarthy Summer School Andy Aschwanden Geophysical nstitute University of Alaska Fairbanks, USA (hopefully) consistent with Continuum Mechanics (Truffer) with lots
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationTHE design of hypersonic vehicles requires accurate prediction
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 22, No. 1, January March 2008 Velocity Slip and Temperature Jump in Hypersonic Aerothermodynamics Andrew J. Lofthouse, Leonardo C. Scalabrin, and Iain D.
More informationPROGRESS IN PROBABILISTIC MODELLING OF ATOMIC SPONTANEOUS EMISSION PROCESSES IN DSMC
PROGRESS IN PROBABILISTIC MODELLING OF ATOMIC SPONTANEOUS EMISSION PROCESSES IN DSMC Dejan Petkow (1), Maria Zaretskaya (1), Alexandre Kling (1), Georg Herdrich (2) (1) Advanced Concepts Team (ACT), ESA-ESTEC,
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 informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
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 informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationMultiphysics Coupling: Hypersonic Flow with Radiation and Ablation
Multiphysics Coupling: Hypersonic Flow with Radiation and Ablation Current Results and Future Strategies Paul T. Bauman, Roy H. Stogner The University of Texas at Austin February 20, 2009 Paul T. Bauman,
More informationAtmospheric Reentry of a Hydrazine Tank
Atmospheric Reentry of a Hydrazine Tank Robert L. Kelley, William C. Rochelle, ESCG Houston, TX Introduction The purpose of this white paper is to describe the methods used to predict the survivability
More informationThermodynamic evolution of phase explosion during high-power nanosecond laser ablation
Thermodynamic evolution of phase explosion during high-power nanosecond laser ablation Quanming Lu* School of Earth and Space Sciences, University of Science and Technology of China, Hefei, 230026, China
More information17-1 Ideal Gases. Gases are the easiest state of matter to describe - All ideal gases exhibit similar behavior.
17-1 Ideal Gases Gases are the easiest state of matter to describe - All ideal gases exhibit similar behavior. An ideal gas is one that is thin enough, that the interactions between molecules can be ignored.
More informationAn improved unified gas-kinetic scheme and the study of shock structures
IMA Journal of Applied Mathematics (2011) 76, 698 711 doi:10.1093/imamat/hxr002 Advance Access publication on March 16, 2011 An improved unified gas-kinetic scheme and the study of shock structures KUN
More informationCH1101 Physical Chemistry Tutorial 1. Prof. Mike Lyons.
CH111 Physical Chemistry Tutorial 1. Prof. Mike Lyons. CH111 Section A Annual 1 Internal Energy Units: Joules J Internal Energy (U) : total kinetic & potential energy of system. e.g. Gas in container with
More informationDSMC simulations of thermal escape
DSMC simulations of thermal escape Alexey N. Volkov, R.E. Johnson, O.J. Tucker, J.T. Erwin Materials Science & Engineering University of Virginia, USA Financial support is provided by NASA through Planetary
More informationUQ in Reacting Flows
UQ in Reacting Flows Planetary Entry Simulations High-Temperature Reactive Flow During descent in the atmosphere vehicles experience extreme heating loads The design of the thermal protection system (TPS)
More informationProspects for High-Speed Flow Simulations
Prospects for High-Speed Flow Simulations Graham V. Candler Aerospace Engineering & Mechanics University of Minnesota Support from AFOSR and ASDR&E Future Directions in CFD Research: A Modeling & Simulation
More informationA Numerical Investigation of Laser Heating Including the Phase Change Process in Relation to Laser Drilling
A Numerical Investigation of Laser Heating Including the Phase Change Process in Relation to Laser Drilling I.Z. Naqavi, E. Savory & R.J. Martinuzzi Advanced Fluid Mechanics Research Group Department of
More informationThe deposition efficiency and spatial thickness distribution of films created by Directed
Chapter 8 Vapor Transport Model Development The deposition efficiency and spatial thickness distribution of films created by Directed Vapor Deposition synthesis have been shown to be sensitive functions
More informationWhy Should We Be Interested in Hydrodynamics?
Why Should We Be Interested in Hydrodynamics? Li-Shi Luo Department of Mathematics and Statistics Center for Computational Sciences Old Dominion University, Norfolk, Virginia 23529, USA Email: lluo@odu.edu
More informationRadiation Heat Transfer
Heat Lectures 0- CM30 /5/06 CM30 ransport I Part II: Heat ransfer Radiation Heat ransfer In Unit Operations Heat Shields Professor Faith Morrison Department of Chemical Engineering Michigan echnological
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationThermal Protection System (TPS) Design and the Relationship to Atmospheric Entry Environments
Thermal Protection System (TPS) Design and the Relationship to Atmospheric Entry Environments By Bernard Laub Dr. Michael J. Wright Dr. Ethiraj Venkatapathy NASA Ames Research Center Moffett Field, CA
More information2σ e s (r,t) = e s (T)exp( rr v ρ l T ) = exp( ) 2σ R v ρ l Tln(e/e s (T)) e s (f H2 O,r,T) = f H2 O
Formulas/Constants, Physics/Oceanography 4510/5510 B Atmospheric Physics II N A = 6.02 10 23 molecules/mole (Avogadro s number) 1 mb = 100 Pa 1 Pa = 1 N/m 2 Γ d = 9.8 o C/km (dry adiabatic lapse rate)
More informationdf dz = dp dt Essentially, this is just a statement of the first law in one of the forms we derived earlier (expressed here in W m 3 ) dq p dt dp
A problem with using entropy as a variable is that it is not a particularly intuitive concept. The mechanics of using entropy for evaluating system evolution is well developed, but it sometimes feels a
More informationMODELING AND SIMULATION OF RADIATION FROM HYPERSONIC FLOWS WITH MONTE CARLO METHODS
The Pennsylvania State University The Graduate School Department of Aerospace Engineering MODELING AND SIMULATION OF RADIATION FROM HYPERSONIC FLOWS WITH MONTE CARLO METHODS A Dissertation in Aerospace
More informationC ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127
C ONTENTS Preface xviii Nomenclature xxvi CHAPTER ONE BASICS OF HEAT TRANSFER 1 1-1 Thermodynamics and Heat Transfer 2 Application Areas of Heat Transfer 3 Historical Background 3 1-2 Engineering Heat
More informationNeeds work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations
Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,
More informationA Comparative Study Between Cubic and Ellipsoidal Fokker-Planck Kinetic Models
A Comparative Study Between Cubic and Ellipsoidal Fokker-Planck Kinetic Models Eunji Jun German Aerospace Center (DLR), 3773 Göttingen, Germany M. Hossein Gorji Computational Mathematics and Simulation
More informationVelocity Slip and Temperature Jump in Hypersonic Aerothermodynamics
5th AIAA Aerospace Sciences Meeting and Exhibit - January 7, Reno, Nevada AIAA 7- Velocity Slip and Temperature Jump in Hypersonic Aerothermodynamics Andrew J. Lofthouse, Leonard C. Scalabrin and Iain
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 information