VERIFICATION AND VALIDATION OF A PARALLEL 3D DIRECT SIMULATION MONTE CARLO SOLVER FOR ATMOSPHERIC ENTRY APPLICATIONS
|
|
- Jordan Owens
- 5 years ago
- Views:
Transcription
1 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 Stuttgart, 569 Stuttgart, Germany nizenkov@irs.uni-stuttgart.de ABSTRACT The in-house Direct Simulation Monte Carlo (DSMC) solver is verified and validated, which enables parallel, three-dimensional simulations of rarefied gas flows on high-performance clusters. Theoretical aspects of the DSMC method and newly employed schemes enabling grid-independence are briefly discussed. Cases including theoretical calculations, experimental measurements and other numerical simulations are selected from literature and compared to the code. Simulation results of the hypersonic flow around a cylinder, rarefied gas flow around a blunted-cone as well as multiple points of the re-entry trajectory of the Orion capsule are presented in terms of drag and heat flux. An outlook on future code development and applications is given.. INTRODUCTION Rarefied, non-equilibrium gas flows as encountered during atmospheric aerocapture and entry manoeuvres pose a great challenge for the design of spacecraft and their thermal protection system. Even in the continuum flow regime, regions of transitional and free-molecular flow are present in the wake of the heat shield, where control elements and/or the payload might be located. Moreover, the complex rarefied flow environment with significant disturbances is outside of the capabilities of conventional continuum methods such as Navier-Stokes based computational fluid dynamics (CFD). A microscopic approach is required to enable an accurate description of such flows. It models fluid properties on a molecular level and is mathematically described by the Boltzmann equation. Based on these ideas, the Direct Simulation Monte Carlo (DSMC) method was introduced in the 96s by Graeme A. Bird and relies on statistical particle treatment and phenomenological modelling of collisions, energy relaxation and chemical reactions []. The DSMC method describes the fluid flow in a discrete manner with simulation particles, where every simulation particle represents a certain amount of real particles. The basic assumptions are that particle movement and collisions can be decoupled and only binary collisions have to be considered. While the DSMC method does not numerically solve the Boltzmann equation, it was shown that in the infinite particle number limit, the DSMC method approximates the Boltzmann equation []. Moreover, the method was further developed and validated in the last 5 years by different research groups. An excellent summary of current developments and capabilities of the method is given in [3]. The Institute of Space Systems (IRS) and the Institute of Aerodynamics and Gas Dynamics (IAG) at the University of Stuttgart cooperatively develop a coupled Particle in Cell (PIC) and DSMC solver called, which enables the simulation of reactive plasma flows [4]. Besides atmospheric entry simulation, the code is applied to the simulation of electric propulsion systems and highlyionized plasmas. The focus of the present paper lies solely on the DSMC implementation of the code, which is referred to as in the following.. SIMULATION The general simulation procedure during one time step is depicted in Fig.. A reservoir or free-stream flow is modelled by simulation particles which are inserted at the beginning of the simulation or every time step, respectively. The number of inserted particles is determined by the number density, weighting factor (ratio of real to simulation particles) and insertion volume. Particles are then moved through the computational domain and possible gas-surface interactions are considered. Subsequently, collision pairs are formed in each cell and a statistical decision process based on random numbers for chemical reactions and treatment of the internal degrees of freedom is conducted. During this process, particles are only associated with microscopic information such as individual internal energy and velocity. Finally, macroscopic values are sampled through spatial and temporal averaging of weighted particle properties before the next time step. The code utilizes several phenomenological models to describe a physical flow. Collision partners are found with the nearest neighbour scheme and an octree-sorting algorithm [5], avoiding the introduction of numerical diffusion if the fluid domain is not resolved with the mean free path. The sorting algorithm accelerates the search for the nearest collision partner and guarantees a physically correct determination of the collision probability, independent of the grid resolution. Collisions are treated with the
2 Particle movement Boundary treatment Sampling Δt Nearest neighbour Octree-sorting Particle pairing Collision process Chemical reactions Relaxation Figure. Schematic of a DSMC time step in. variable hard sphere (VHS) model, where the collisional cross-section depends on the relative translational energy and probabilities are calculated by the Natural Sample Size method [6]. Energy transfer between kinetic and internal modes is treated with the Larsen-Borgnakke model []. While rotational energy modes are assumed to be continuous, vibrational modes are quantized with the simple harmonic oscillator. The probabilities for rotational and vibrational relaxation are fixed values of P rot =. and P vib =.4 based on empirical considerations. Chemical reactions are enabled through the use of the extended Arrhenius equation. Surface interactions are modelled by the Maxwell model with surface accommodation coefficients, currently without any catalytic effects. The probabilistic nature of the method introduces statistical fluctuations for steady-state problems. An increase of the sample size (longer sampling duration and/or more simulation particles) can reduce these fluctuations. The quality of simulation results is established by several measures. To avoid a non-physical time step, maximal collision probabilities are monitored (P max < ) and the time step adjusted accordingly. The weighting factor is chosen to achieve a meaningful statistical ensemble. Specifically, a sufficient number of simulation particles per cell in low density regions (with respect to the free-stream) has to be ensured. Although a generally applicable value cannot be given, mean values of N part,cell > 5 have proven to produce acceptable results. To achieve an accurate transient simulation, an increased number of particles has to be present at each time step in lieu of temporal averaging. For steady-state problems, convergence is determined based on the transient development of the variable of interest. 3. VERIFICATION AND VALIDATION This section presents cases collected from literature to verify and validate different capabilities of the code. In Section 3., the chemistry model is verified by the means of comparison to semi-empirical and numerical solutions of reservoir simulations, while the following sections present more complex geometries such as a cylinder, a spherically blunted cone and the Orion capsule. In addition, simulated flow conditions increase in difficulty starting with a monatomic gas advancing to a diatomic gas with internal energy relaxation and finally high-temperature air including chemical reactions. 3.. Chemistry Model The chemistry module of the DSMC code is verified based on simple reservoir simulations to exclude the influence of a complex geometry. In the simulation of entry manoeuvres, dissociation and recombination of molecules is of great importance. The first case focuses on such a process for a single species A A + A. The degree of dissociation can be determined by α = n A n A + n A, () where n A and n A are the number densities of atoms and molecules, respectively. The equilibrium composition of a dissociating gas can be described as a function of density and temperature [] α α = ρ ( d ρ exp T ) d, () T where the subscript d denotes species characteristic values, for example T d = 3 K and ρ d =.3 5 kg m 3 for nitrogen. Although the characteristic density ρ d is dependent on the prevalent temperature, differences are found to be negligible and a constant value is chosen []. For the nitrogen reservoir calculation with, a cube with a side length of 5 m and specular walls is simulated. The viscosity coefficient ω =. and the characteristic temperature T ref = K of the simulated species have to be equal due to the VHS model, resulting in reference diameters of d ref = 3 m and d ref = 4. m for atomic and diatomic nitrogen, respectively. The characteristic vibrational temperature for diatomic nitrogen is set to Θ vib = 3395 K. The density ratio ρ d /ρ is chosen as, resulting in a number density of n =.9 3 m 3. The employed reaction rates are summarized in Tab.. A range of temperatures was simulated including two dissociative and recombination reactions each. After thermal and chemical equilibrium was observed, the degree of dissociation was calculated with Eq. and compared to the equilibrium composition resulting from Eq.. Generally, good agreement between the theoretical and numerical solution, as shown in Fig., is observed. While a slightly earlier onset of dissociation can be observed for the DSMC Table. Arrhenius coefficients for nitrogen dissociation and recombination. A [ m 3 s, m 6 s ] b [ ] E a [K] N + N N + N N + N + N N + N + N. 39.6
3 Degree of dissociation α [ ] Theory Temperature T [K] Figure. Degree of dissociation of diatomic nitrogen at equilibrium, ρ d /ρ =. simulation, differences are within reason considering uncertainties in the input parameters such as the characteristic values for Eq. and the Arrhenius coefficients for the DSMC simulation. The next test case considers a transient reservoir simulation of five-species air (N,O, N, O, NO) with 34 potential chemical reaction paths and an initial composition of 9% N and % O. The case is representative for chemical processes occurring during high-enthalpy atmospheric entry. Three different initial conditions (T = K, K, and 3 K at p =.35 kpa) were simulated with and compared to a numerical solution which can be found by means of simple first-order Euler integration [9], if thermal equilibrium is assumed (equilibrium between different energy modes). Simulation parameters such as number density n, cube side length λ and total initial number of particles are given in Tab.. While the same species parameters are utilized for N and N as in the previous test case, the reference diameter for O, NO and O at T ref = K and ω =. is d ref = 3.96 m, 4. m, and 3. m, respectively. Characteristic vibrational temperatures are Θ vib,o = 39 K and Θ vib,no = K. Arrhenius coefficients are set according to []. To augment thermal equilibrium, the rotational and vibrational relaxation probabilities were set to unity. Since this does not guarantee thermal equilibrium, total temperature of each species was calculated with T tot = 3T tra + T rot + ξ vib T vib 5 + ξ vib, (3) where ξ vib is the corresponding vibrational degree of freedom and subscripts tra, rot and vib denote the transla- Table. Parameters for the reservoir simulation of fivespecies air at p =.35 kpa. T [K] n [ m 3 ] λ [m] N part [ ] tional, rotational, and vibrational energy mode. Total reservoir temperature was determined by the particle-weighted mean of species temperatures. Fig. 3 shows the transient development of the species composition at an initial temperature of 3 K. The number density of the respective species is normalized by the initial number density. Excellent agreement is found between the numerical solution and particle simulation. An early onset of O and N dissociation followed by formation of NO can be observed as expected for high-temperature air. Moreover, the deviation from thermal equilibrium does not seem to have a great effect on the species composition when compared to the numerical solution. Statistical noise can be seen in the O and NO concentration towards chemical equilibrium due to very low particle numbers of the respective species and consequently an insufficient ensemble size. However, the simulation result appears to be oscillating around the numerical solution. ns/n [ ] N 3 4 O O NO N Numerical 9 6 Time [s] Figure 3. Species composition of a five-species air reservoir during chemical relaxation (T = 3 K). The transient species composition of a reservoir with an initial temperature of K is shown in Fig. 4. Dissociation of both diatomic species is considerably slower compared to the previous test case. Consequently, results for the atomic nitrogen are accompanied by statistical noise due to low particle numbers. Nonetheless, the particle simulation agrees very well with the numerical solution. The temperature development is of great interest in order to verify the implementation of the chemistry as well as relaxation module. Fig. 5 illustrates the transient reservoir temperature for the three different initial conditions. The equilibrium temperature predicted by the numerical solution is compared to the total reservoir temperature determined by Eq. 3. The agreement between both computations is excellent for all three initial temperatures. In conclusion, the chemistry module is able to reproduce Arrhenius type reaction rates and the corresponding energy variation in a reservoir simulation. Transient species composition and temperature are in excellent agreement with theoretical results.
4 N ns/n [ ] 3 4 O O N NO Numerical 9 6 Time [s] Figure 4. Species composition of a five-species air reservoir during chemical relaxation (T = K). Temperature Teq,Ttot [K] K K K Numerical 9 6 Time [s] Figure 5. Temperature comparison of a five-species air reservoir during chemical relaxation. 3.. Cylinder Simulation results of a simple, two-dimensional cylinder exposed to a Mach argon flow (u = 634. m s, T = K) are compared to other simulations []. The usage of the monatomic argon eliminates the possibility of internal energy relaxation and simplifies the code-to-code verification process. Results of a Navier- Stokes solver DPLR, developed at NASA s Ames Research Center, and the DSMC code MONACO, developed by Iain D. Boyd and his research group at the University of Michigan, are utilized for the verification. Two different flow regimes, characterized by the Knudsen number, Kn =. and Kn =.5, are investigated. While the first test case is in the continuum regime with a number density of n = 4.4 m 3, the second test case represents a rarefied gas flow at a number density of n = m 3. Species parameters are chosen according to []. The cylinder diameter is.34 m (= in) with complete diffuse reflection as well as thermal accommodation and a constant wall temperature of 5 K. The nomenclature and geometry are given in Fig. 6. M Φ Figure 6. Definition of cylinder test case. At first, the drag of the cylinder given per unit length as a two-dimensional problem is considered. Excellent agreement between all three simulations can be seen in Tab. 3 at the lower Knudsen number Kn =.. Although both DSMC simulations compute a lower total drag per length with respect to the CFD simulation, the difference is negligible. For the rarefied gas flow at Kn =.5, both DSMC simulations compare very well but a strong deviation of the CFD results can be observed due to continuum breakdown. Table 3. Comparison of the total drag/length [N m ]. Kn MONACO DPLR < % 4. < %.5.. < %.64 % In the following, the heat flux along the cylinder angle Φ is compared. As depicted in Fig., a good agreement between all codes as expected due to flow conditions close to the continuum regime is achieved. The peak heat flux is lower by < % and 4% than the other DSMC and CFD code, respectively. The results of the test case in the transitional flow regime at Kn =.5 can be seen in Fig.. Again, good agreement is achieved between the DSMC codes with a deviation of just % at the stagnation point, while the CFD code overpredicts the heat flux due to the breakdown of the continuum assumption. The peak heat flux at the stagnation point is overpredicted by 35%, with an increase of the relative deviation further downstream in accordance with an increasingly rarefied gas flow. In conclusion, the first results of an incident flow suggest a correct implementation of the elastic collision module as well as wall treatment routines in and allow the advancement to more complex problems. 4 3 MONACO DPLR Φ [ ] Figure. Heat flux along cylinder surface, Kn =..
5 6 5 4 MONACO DPLR Φ [ ] Figure. Heat flux along cylinder surface, Kn = Blunted Cone The next validation case is a spherically blunted cone geometry in a Mach molecular nitrogen flow. As opposed to the previous case, relaxation of internal energy modes such as rotation and vibration is now possible. al results [] are available for three different inflow conditions and different angles of attack (α =,,, and 3 ). The experimental test conditions are given in [] and the derived free-stream conditions for the simulation are shown in Tab. 4. The VHS model parameters of nitrogen are d ref = 4. m as the reference diameter at a temperature of T ref = 3 K and a viscosity coefficient of ω =.4. The characteristic vibrational temperature is set to Θ v = 3395 K. Table 4. Free-stream conditions of the nitrogen flow. Case v [ m s ] T [K] n [m 3 ] Set Set Set The model geometry and positions of the thermocouples can be seen in Fig. 9, while details of the experimental setup can be found in []. al measurements were taken on the forebody (points -4), the base (points 5-6) and the sting (-9). A global uncertainty of ±% is given in [] for all experimental measurements of the heat flux except when stated otherwise. The temperature of the thin-walled model used in the experiments fluctuated between 9 K and 3 K, thus a constant wall temperature of 3 K was chosen for the simulation. Full thermal and complete diffuse reflection was assumed. For the simulation without an angle of attack, only / (= 3 slice) of the axisymmetric fluid domain was simulated while for the other cases half of the fluid domain was modelled, taking advantage of the symmetry in the xy-plane. Beginning with the first set of test conditions (n = 3.5 m 3 ), the heat flux along the surface of the model at α = is shown in Fig.. Labels of the experimental data points correspond to the positions of the thermocouples and the path variable S is normalized with the nose radius R n =.5 m. Good qualitative agreement can be seen on the forebody of the blunted cone M [mm] R b 5. R c.5 R j. R n.5 R s 6.5 y S R b 3 R c R n z x Figure 9. blunted cone geometry and positions of thermocouples. R j R s # S/R n as well as on the sting in the wake. While the error at the first three measurement points is within the experimental uncertainty of ±%, the heat flux is overpredicted by 9% at the fourth thermocouple. al results of thermocouples 5 and 6 were given as a maximal value (<. kw m ) since experimental accuracy cannot be guaranteed due to very low heating rates. On the sting, the error increases from 6% to 4% for the last three thermocouples Figure. Comparison of heat flux along surface of blunted cone (Set, α = ). Fig. shows the comparison for the second set of test conditions (n =.5 m 3 ) at α =. Again, good qualitative agreement with the experiment can be observed on the forebody, base and sting of the model. The error compared to the first four thermocouples is below %. Although the measurements on the base are still very low and expected to be equal or less than.4 kw m, the error is below %. Heating rates at thermocouples and deviate by less than 3%. The last measurement cannot be matched by the simulation with almost an error of 5%. Since the correct tendency is predicted for the other flow conditions, especially at the sting, it is not clear whether the deviation is due to inaccu-
6 rate simulation results or uncertainties in the experimental measurements. Except for the last two measurements points, the particle simulation was able to stay within the experimental uncertainty for the second set of flow conditions. The heat flux at the most dense flow conditions (n =. m 3 ) at α = is shown in Fig.. Although, the heat flux on the forebody is overpredicted by % to 4%, good quantitative agreement can be seen on the base with errors below 6%. Errors below % and thereby within the experimental uncertainty are determined on the sting of the model. As noted in [], high wall-temperature gradients along the forebody may have possibly led to local conduction effects resulting in underestimated heating rates with the strongest deviation up to thermocouple 4. This assertion corresponds to the simulation results of and an other DSMC code [3] Figure. Comparison of heat flux along surface of blunted cone (Set, α = ) Figure. Comparison of heat flux along surface of blunted cone (Set 3, α = ). Finally, Fig. 3 illustrates results of full, three-dimensional simulations of Set at angles of attack α = and α = 3. As logarithmic plots can be deceiving and the values for thermocouples 5 and 6 can be excluded from the comparison of Set, regular axes are utilized. Good qualitative agreement can be observed for both angles of attack. Excellent error values below 5% are achieved on the forebody at α = and acceptable agreement with errors below 3% can be seen on the sting. A very different heat flux behaviour is observed at α = 3, with the maximum heat flux shifted to the shoulder of the forebody. While the error stays below 3% for the first four thermocouples, better agreement is achieved at the sting with the last two thermocouples within the experimental uncertainty. Concluding, this test case demonstrated the ability of to model flows of diatomic molecules with internal degrees of freedom around complex, threedimensional geometry (a) α = (b) α = 3 Figure 3. Comparison of heat flux along surface of blunted cone at different angles of attack (Set ) Orion Capsule The last comparison is based on the Orion Crew Module (CM) and the comparison with other DSMC and CFD simulations [4]. The free-stream conditions given in Tab. 5 feature a hypersonic flow with a wide range of flow conditions from the continuum to the rarefied gas regime. Fig. 4 shows the geometric definition of the capsule. The simulation included high-temperature chemistry with five-species air (parameters identical to the reservoir test case in Section 3.) and 34 chemical reactions (Arrhenius coefficients []) while ionization was neglected as it was found to have little influence on aerodynamic coefficients. The initial atmospheric composition for the corresponding altitude can be found in [4]. Full diffuse reflection and
7 Table 5. Free-stream conditions and wall temperature for Orion CM (u = 6 m s, α = 6 ) [4]. Alt. [km] n [m 3 ] T [K] Kn [-] T w [K] [m] R b.546 R n 6.35 R s.55 R as.5 L 3.3 β 3.5 α R b S R s R n β R as Drag coefficient [ ] LAURA DS3V DAC 4 6 Altitude [km] Figure 5. Drag coefficient over altitude of Orion CM (u = 6 m s, α = 6 ). To summarize, the presented DSMC code was applied to a complex geometry such as the Orion re-entry capsule. Results of the drag coefficient showed very good agreement with well-established CFD and DSMC codes, free-stream conditions ranging from the continuum to rarefied gas flow regime. It should be noted that although the DSMC method is applicable in the continuum regime, it comes with high computational effort due to high particle numbers and a low required time step. For example, the total simulation duration of the most dense flow at 5 km altitude with simulation particles and a time step of t = 4 9 s was over hours on 4 Intel Xeon CPU E5-6 v3 processors (96 cores). Figure 4. Geometric definition of the Orion CM. complete thermal accommodation was assumed with an isothermal wall temperature T w presented in Tab. 5. The drag coefficient was calculated based on the base radius R b, free-stream velocity u and density ρ. The study [4] included simulations with two DSMC codes, DS3V and DAC, and the CFD code LAURA. DS3V, which is freely available, is developed by Graeme A. Bird and DAC by G. J. LeBeau from NASA s Johnson Space Center. The CFD code LAURA is developed by NASA s Langley Research Center. Details regarding these codes and respective references are available in [4]. While calculations with the continuum approach were conducted from 5 km up to 5 km altitude, particle-based simulations ranged from 5 km to km altitude. The results of the simulations in terms of drag coefficient over altitude compared to are presented in Fig. 5. Excellent agreement of the drag coefficient is observed between the DSMC codes at all altitudes. The deviation from the DAC and DS3V results is below % and 3%, respectively. Furthermore, excellent agreement with the CFD code is achieved at lower altitudes of 5 km and 95 km, with a deviation of less than %. At higher Knudsen number, the breakdown of continuum is evident and the difference increases to roughly %. 4. CONCLUSION The capabilities of were presented by means of several verification and validation cases. Through reservoir simulations, the chemistry module was compared to semi-empirical and numerical relations, showing excellent agreement. Following that, more complex geometries and flow conditions were investigated. First, simulations of hypersonic flow of argon in the continuum and transitional regime around a cylinder were compared to results of other DSMC as well as CFD codes. Very good agreement with other DSMC simulations is achieved whereas good agreement with CFD results is only observed in the continuum regime. Second, a comparison with experimental measurements of the heat flux on a spherically blunted cone and its sting in the wake was made. The results of the threedimensional simulation of the rarefied, hypersonic nitrogen flow compared well on the forebody and sting of the model. Finally, simulation of multiple points of Orion s reentry trajectory comprised a full, three-dimensional case with high-temperature air chemistry and was compared to other CFD and DSMC simulations. To summarize the capabilities of the presented DSMC code: Full, three-dimensional simulation, Parallel computation on high-performance clusters, Chemistry modelling of diatomic molecules, and Grid-independence enabled through nearest neighbour scheme with octree-sorting algorithm.
8 The code is subject to on-going research, some of which is applicable to atmospheric entry. This includes the extension to polyatomic molecules required for the simulation of Mars and Titan entry. Automatic convergence detection and variable relaxation probabilities are implemented as part of student research projects. Moreover, coupling with the Low Diffusion method, applicable in the continuum regime, will enable the efficient simulation of flows with large density gradients, such as the flow around blunted entry vehicles. In addition, the application of the coupled PIC-DSMC solver to electromagnetic shielding concepts is envisioned. ACKNOWLEDGEMENTS The authors would like to thank Airbus Defence and Space, namely Dr. Martin Konopka and Peter Noeding, for supporting the research in the framework of a research training group. The authors are also grateful for the financial support by the Landesgraduiertenförderung Baden-Württemberg. Computational resources have been provided by the High Performance Computing Center Stuttgart (HLRS) of the University of Stuttgart. REFERENCES [9] Haas, B.L. & McDonald, J.D. (993). Validation of chemistry models employed in a particle simulation method. Journal of Thermophysics and Heat Transfer (), 4 4. [] Haas, B.L. (99). Fundamentals of chemistry modeling applicable to a vectorized particle simulation, in 5th Joint Thermophysics and Heat Transfer Conference, AIAA, Reston, Virigina. [] Lofthouse, A.J., Boyd, I.D. & Wright, M.J. (). Effects of continuum breakdown on hypersonic aerothermodynamics. Physics of Fluids 9(). [] Allègre, J., Bisch, D. & Lengrand, J.C. (99). al Rarefied Heat Transfer at Hypersonic Conditions over -Degree Blunted Cone. Journal of Spacecraft and Rockets 34(6), 4. [3] Moss, J.N., Dogra, V.K., Price, J.M. & Hash, D.B. (995). Comparison of DSMC and experimental results for hypersonic external flows, in 3th AIAA Thermophysics Conference, AIAA, Reston, Virigina. [4] Moss, J.N., Boyles, K. & Greene, F.A. (6). Orion Aerodynamics for Hypersonic Free Molecular to Continuum Conditions, in 4th AIAA/AHI Space Planes and Hypersonic Systems and Technologies Conference, AIAA, Reston, Virigina, pp [] Bird, G.A. (994). Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford. [] Wagner, W. (99). A convergence proof for Bird s direct simulation Monte Carlo method for the Boltzmann equation. Journal of Statistical Physics 66(3-4), 44. [3] Boyd, I.D. (4). Computation of Hypersonic Flows Using the Direct Simulation Monte Carlo Method. Journal of Spacecraft and Rockets (Articles in Advance), 6. [4] Munz, C.D., Auweter-Kurtz, M., Fasoulas, S. et al. (4). Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows. Comptes Rendus Mécanique 34(-), [5] Pfeiffer, M., Mirza, A. & Fasoulas, S. (3). A gridindependent particle pairing strategy for DSMC. Journal of Computational Physics 46, 36. [6] Baganoff, D. & Mcdonald, J.D. (99). A collisionselection rule for a particle simulation method suited to vector computers. Physics of Fluids A: Fluid Dynamics (), [] Borgnakke, C. & Larsen, P.S. (95). Statistical collision model for Monte Carlo simulation of polyatomic gas mixture. Journal of Computational Physics (4), [] Lighthill, M.J. (95). Dynamics of a dissociating gas Part I Equilibrium flow. Journal of Fluid Mechanics (), 3.
arxiv: v1 [physics.comp-ph] 12 Nov 2018
arxiv:1811.04742v1 [physics.comp-ph] 12 Nov 2018 Load balancing strategies for the DSMC simulation of hypersonic flows using HPC T. Binder 1, S. Copplestone 2, A. Mirza 1, P. Nizenkov 1, P. Ortwein 2,
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationLecture1: Characteristics of Hypersonic Atmosphere
Module 1: Hypersonic Atmosphere Lecture1: Characteristics of Hypersonic Atmosphere 1.1 Introduction Hypersonic flight has special traits, some of which are seen in every hypersonic flight. Presence of
More information21st International Symposium on Rareed Gas Dynamics half-angle cone forebodies, while Stardust and Microprobe have forebody cone half-angles of 60 and
Rareed Transitional Bridging of Blunt Body Aerodynamics R. G. Wilmoth, R. C. Blanchard, J. N. Moss NASA Langley Research Center, Hampton, VA, USA 1 Introduction Bridging relations are frequently used to
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 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 informationDSMC solver for an optimized Space Crew reentry Orbital Vehicle
DSMC solver for an optimized Space Crew reentry Orbital Vehicle Vigneshwaran Krishnamurthy 1,a, Ram Arvinth Shanmugam Periasamy 1,b, Vertika Saxena 1,c, Nandita Nurani Hari 1,d and Suriyaprabha Chandrasekar
More informationEffects of Continuum Breakdown on Hypersonic Aerothermodynamics
Effects of Continuum Breakdown on Hypersonic Aerothermodynamics Andrew J. Lofthouse* and lain D. Boydt Department ofaerospace Engineering University of Michigan, Ann Arbor, MI, 48109-2140 Michael J. Wrightl
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 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 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 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 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 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 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 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 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 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 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 informationReconsideration of DSMC Models for Internal Energy Transfer and Chemical Reactions
Reconsideration of Models for Internal Energy Transfer and Chemical Reactions N.E. Gimelshein, S.F. Gimelshein, D.A. Levin, M.S. Ivanov and I.J. Wysong Penn State University, University Park, PA 16802
More informationGas-Surface Interaction Effect on Round Leading Edge Aerothermodynamics
Brazilian Journal of Physics, vol. 37, no. 2A, June, 2007 337 Gas-Surface Interaction Effect on Round Leading Edge Aerothermodynamics Wilson F. N. Santos National Institute for Space Research, Cachoeira
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 informationApplication of the Transition Probability Matrix Method to High Knudsen Number Flow Past a Micro-Plate
Application of the Transition Probability Matrix Method to High Knudsen Number Flow Past a Micro-Plate Andrew J. Christlieb, W. Nicholas G. Hitchon, Quanhua Sun and Iain D. Boyd Department of Aerospace
More informationGas-Surface Interaction Effect on Aerodynamics of Reentry Vehicles
International Journal of Educational Research and Information Science 2018; 5(1): 1-9 http://www.openscienceonline.com/journal/eris Gas-Surface Interaction Effect on Aerodynamics of Reentry Vehicles Zay
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 informationTRAJECTORY AND AEROTHERMODYNAMIC ANALYSIS OF TOWED-BALLUTE AEROCAPTURE USING DIRECT SIMULATION MONTE CARLO
AAS 07-307 TRAJECTORY AND AEROTHERMODYNAMIC ANALYSIS OF TOWED-BALLUTE AEROCAPTURE USING DIRECT SIMULATION MONTE CARLO Kristin L. Gates Medlock, * Alina A. Alexeenko, and James M. Longuski We investigate
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 informationA Hybrid Continuum / Particle Approach for Micro-Scale Gas Flows
A Hybrid Continuum / Particle Approach for Micro-Scale Gas Flows Quanhua Sun *, Iain D. Boyd * and Graham V. Candler * Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 489 Department
More informationDSMC Study of Flowfield and Kinetic Effects on Vibrational Excitations in Jet-Freestream Interactions
DSMC Study of Flowfield and Kinetic Effects on Vibrational Excitations in Jet-Freestream Interactions David H. Campbell a and Ingrid J. Wysong b a ERC, Incorporated Air Force Research Laboratory, Edwards
More informationDirect Molecular Simulation of Hypersonic Flows
Direct Molecular Simulation of Hypersonic Flows Thomas E. Schwartzentruber Aerospace Engineering & Mechanics University of Minnesota UMN Students and Researchers: Savio Poovathingal 1 Paul Norman 3 Chonglin
More informationCopyright 2015 American Institute of Aeronautics and Astronautics, Inc.
nn Scanlon, T. J., White, C., Borg, M. K., Palharini, R. C., Farbar, E., Boyd, I. D., Reese, J. M., and Brown, R. E. (2015) Open souce DSMC chemistry modelling for hypersonic flows. AIAA Journal. Copyright
More informationComputational Modeling of Hypersonic Nonequilibrium Gas and Surface Interactions
Computational Modeling of Hypersonic Nonequilibrium Gas and Surface Interactions Iain D. Boyd, Jae Gang Kim, Abhilasha Anna Nonequilibrium Gas & Plasma Dynamics Laboratory Department of Aerospace Engineering
More informationInvestigation of Drag and Heat Transfer for Martian Dust Particles
Investigation of Drag and Heat Transfer for Martian Dust Particles T. Ozawa, T. Suzuki, H. Takayanagi, and K. Fujita Aerospace Research and Development Directorate Japan Aerospace Exploration Agency, Chofu,
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 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 informationThermochemical and Particulate Interfacing for Hybrid High-Altitude Plume and Control Jet Simulations
Thermochemical and Particulate Interfacing for Hybrid High-Altitude Plume and Control Jet Simulations D.B. VanGilder, R.G. Wilmoth, and J.L. Papp Combustion Research and Flow Technology, Inc., 6210 Keller
More informationStefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran
Stefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran Outlines: Introduction DSMC Collision Schemes Results Conclusion
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 informationDSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations
DSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations J. R. Torczynski and M. A. Gallis Engineering Sciences Center, Sandia National Laboratories, P. O. Box
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 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 informationSIMULATION TECHNIQUES IN HYPERSONIC AEROTHERMODYNAMICS
ICAS CONGRESS SIMULATION TECHNIQUES IN HYPERSONIC AEROTHERMODYNAMICS Vladimir V. Riabov Rivier College, Nashua, New Hampshire 36, USA Keywords: hypersonic non-equilibrium rarefied-gas flows, aero- & thermodynamic
More informationNumerical Simulation of Rarefied-Gas Flows about a Rotating Cylinder
Numerical Simulation of Rarefied-Gas Flows about a Rotating Cylinder Vladimir V. Riabov Department of Computer Science, Rivier College, 42 South Main Street, Nashua, NH 36-86, USA Abstract. Subsonic and
More informationRarefaction Effects in Hypersonic Aerodynamics
Rarefaction Effects in Hypersonic Aerodynamics Vladimir V. Riabov Department of Mathematics and Computer Science, Rivier College, 4 S. Main St., Nashua, NH 6, USA Abstract. The Direct Simulation Monte-Carlo
More informationHYPERSONIC AERODYNAMICS OF TOROIDAL BALLUTES
7 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES HYPERSONIC AERODYNAMICS OF TOROIDAL BALLUTES Vladimir V. Riabov Rivier College, Nashua, New Hampshire 36, USA Keywords: DSMC method, aerodynamic
More informationGeneration of a Hybrid DSMC/CFD Solution for Gas Mixtures with Internal Degrees of Freedom
Generation of a Hybrid DSMC/CFD Solution for Gas Mixtures with Internal Degrees of Freedom Kelly A. Stephani, David B. Goldstein and Philip L. Varghese Department of Aerospace Engineering, The University
More informationNumerical Simulation of Microwave Plasma Thruster Flow
Numerical Simulation of Microwave Plasma Thruster Flow IEPC-2007-211 Presented at the 30 th International Electric Propulsion Conference, Florence, Italy September 17-20, 2007 Mao-lin Chen *, Mao Gen-wang,
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 informationSIMULATION OF GAS FLOW OVER MICRO-SCALE AIRFOILS USING A HYBRID CONTINUUM-PARTICLE APPROACH
33rd AIAA Fluid Dynamics Conference and Exhibit 3-6 June 3, Orlando, Florida AIAA 3-44 33 rd AIAA Fluid Dynamics Conference and Exhibit / Orlando, Florida / 3-6 Jun 3 SIMULATION OF GAS FLOW OVER MICRO-SCALE
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 informationSMILE System for 2D/3D DSMC Computations
SMILE System for 2D/3D DSMC Computations M.S. Ivanov 1, A.V. Kashkovsky 1, S.F. Gimelshein 2, G.N. Markelov 3, A.A. Alexeenko 2, Ye.A. Bondar 1, G.A. Zhukova 1, S.B. Nikiforov 1, P.V. Vaschenkov 1 1 Khristianovich
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 informationDirect Molecular Simulation of Nonequilibrium Flows
Direct Molecular Simulation of Nonequilibrium Flows Thomas E. Schwartzentruber Aerospace Engineering and Mechanics University of Minnesota AFOSR Young Investigator Program Grant # FA9550-10-1-0075 Direct
More informationComparison and Analysis of Chemical Models, methodologies and parameters in Hypersonic Rarefied Flows
University of Naples Federico II Ph.D. thesis Comparison and Analysis of Chemical Models, methodologies and parameters in Hypersonic Rarefied Flows by Luigi Morsa Coordinator: Prof. Antonio Moccia Tutor:
More informationSIMULATION OF TITAN ATMOSPHERE BY AN ARC-HEATED FACILITY
SIMULATION OF TITAN ATMOSPHERE BY AN ARC-HEATED FACILITY V. Carandente (1), V. Caso (2), A. Esposito (3), R. Savino (4), G. Zuppardi (5) University of Naples "Federico II", P.le V. Tecchio, 80 80125 Naples
More informationComputation of Shock Waves in Inert Binary Gas Mixtures in Nonequilibrium Using the Generalized Boltzmann Equation
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0361 Computation of Shock Waves in Inert Binary Gas Mixtures
More informationInstitute of Mechanics Lomonosov Moscow State University Institute for Problems in Mechanics of Russian Academy of Sciences
Institute of Mechanics Lomonosov Moscow State University Institute for Problems in Mechanics of Russian Academy of Sciences HEAT TRANSFER IN UNDER-EXPANDED NON-EQUILIBRIUM CARBON DIOXIDE JETS: AN EXPERIMENT
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 informationFlow Simulation over Re-Entry Bodies at Supersonic & Hypersonic Speeds
International Journal of Engineering Research and Development eissn : 2278-067X, pissn : 2278-800X, www.ijerd.com Volume 2, Issue 4 (July 2012), PP. 29-34 Flow Simulation over Re-Entry Bodies at Supersonic
More informationCFD AEROTHERMAL ANALYSIS OF SPACE CAPSULE RE-ENTRY INTO EARTHS ATMOSPHERE AT DIFFERENT ALTITUDES
CFD AEROTHERMAL ANALYSIS OF SPACE CAPSULE RE-ENTRY INTO EARTHS ATMOSPHERE AT DIFFERENT ALTITUDES V. Jyothi Swaroop 1 & N. Mahesh kumar 2 M.tech 1, Associate Professor 2 Sri Venkateswara College of Engg
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 informationA COMBINED CFD-DSMC METHOD FOR NUMERICAL SIMULATION OF NOZZLE PLUME FLOWS
A COMBINED CFD-DSMC METHOD FOR NUMERICAL SIMULATION OF NOZZLE PLUME FLOWS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements
More informationStati sti cal si mul ati on of reentry capsul e aero dynami cs i n hyp ersoni c near-conti nuum flows
Stati sti cal si mul ati on of reentry capsul e aero dynami cs i n hyp ersoni c near-conti nuum flows Mikhail S. Ivanov Computational Aerodynamics Laboratory, Khristianovich Institute of Theoretical and
More informationNUMERICAL SIMULATION OF KINETIC EFFECTS IN LOW-DENSITY HYPERSONIC AERODYNAMICS
25 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL SIMULATION OF KINETIC EFFECTS IN LOW-DENSITY HYPERSONIC AERODYNAMICS Vladimir V. Riabov Rivier College, Nashua, New Hampshire 36, USA
More informationDetached Eddy Simulation on Hypersonic Base Flow Structure of Reentry-F Vehicle
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (2014) 000 000 www.elsevier.com/locate/procedia APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology,
More informationParticle Simulation of Hall Thruster Plumes in the 12V Vacuum Chamber
Particle Simulation of Hall Thruster Plumes in the 12V Vacuum Chamber IEPC-2005-138 Presented at the 29 th International Electric Propulsion Conference, Princeton University, Iain D. Boyd 1, Quanhua Sun
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 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 informationAn Open-Source Hybrid CFD-DSMC Solver for High Speed Flows
Espinoza, Daniel and Casseau, Vincent and Scanlon, Thomas and Brown, Richard (2016) An open-source hybrid CFD-DSMC solver for high speed flows. In: Proceedings of the 30th International Symposium on Rarefied
More informationA MODULAR PARTICLE-CONTINUUM NUMERICAL ALGORITHM FOR HYPERSONIC NON-EQUILIBRIUM FLOWS
A MODULAR PARTICLE-CONTINUUM NUMERICAL ALGORITHM FOR HYPERSONIC NON-EQUILIBRIUM FLOWS by Thomas Edward Schwartzentruber A dissertation submitted in partial fulfillment of the requirements for the degree
More informationThree fast computational approximation methods in hypersonic aerothermodynamics
819 Three fast computational approximation methods in hypersonic aerothermodynamics V.V. Riabov* Rivier College, Computer Science Department, 420 S. Main Street, Nashua, NH 03060, USA Abstract The applicability
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 informationCollisional energy transfer modeling in non-equilibrium condensing flows
Collisional energy transfer modeling in non-equilibrium condensing flows Natalia Gimelshein, Ingrid Wysong and Sergey Gimelshein ERC, Inc, Edwards AFB, CA 93524 Propulsion Directorate, Edwards AFB, CA
More informationApplication of Rarefied Flow & Plasma Simulation Software
2016/5/18 Application of Rarefied Flow & Plasma Simulation Software Yokohama City in Japan Profile of Wave Front Co., Ltd. Name : Wave Front Co., Ltd. Incorporation : March 1990 Head Office : Yokohama
More informationResearch Article. Slip flow and heat transfer through a rarefied nitrogen gas between two coaxial cylinders
Available online wwwjocprcom Journal of Chemical and Pharmaceutical Research, 216, 8(8):495-51 Research Article ISSN : 975-7384 CODEN(USA) : JCPRC5 Slip flow and heat transfer through a rarefied nitrogen
More informationHypersonic flow and flight
University of Stuttgart, Aerospace Engineering and Geodesy Dept. - Lecture - Hypersonic flow and flight Master Level, Specialization 4 lecture hours per week in WS, 3-6 LPs/ECTS Lecturer: Dr. Markus J.
More informationParticle Photon Monte Carlo Method for Simulating Atomic Radiation in Hypersonic Reentry Flows in DSMC
4th Thermophysics Conference, 23-26 June 28, Seattle, Washington Particle Photon Monte Carlo Method for Simulating Atomic Radiation in Hypersonic Reentry Flows in DSMC T. Ozawa, A. Wang, + D. A. Levin
More informationNumerical Studies of Heat Transfer on Hypersonic Blunt Bodies with Gas Injection
Paper 130 Civil-Comp Press, 2010 Proceedings of the Seventh International Conference on Engineering Computational Technology, B.H.V. Topping, J.M. Adam, F.J. Pallarés, R. Bru and M.L. Romero, (Editors),
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 informationMonte Carlo Simulation of Hypersonic Rarefied Gas Flow using Rotationally and Vibrationally Inelastic Cross Section Models
Monte Carlo Simulation of Hypersonic Rarefied Gas Flow using Rotationally and Vibrationally Inastic Cross Section Mods Hiroaki Matsumoto a a Division of systems research, Yokohama National University 79-5
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 informationSimulation of Homogeneous Condensation of Ethanol in High Pressure Supersonic Nozzle Flows using BGK Condensation Model
Simulation of Homogeneous Condensation of Ethanol in High Pressure Supersonic Nozzle Flows using BGK Condensation Model Rakesh Kumar and D. A. Levin Pennsylvania State University, University Park, PA 1682
More informationGridless DSMC. Spencer Olson. University of Michigan Naval Research Laboratory Now at: Air Force Research Laboratory
Gridless DSMC Spencer Olson University of Michigan Naval Research Laboratory Now at: Air Force Research Laboratory Collaborator: Andrew Christlieb, Michigan State University 8 June 2007 30 June 2009 J
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 information