Simulation of Weakly Ionized Hypersonic Flows with Reactive Species Weighting Scheme in the Direct Simulation Monte Carlo Method
Abstract
The conservative species weighting scheme enables the direct simulation Monte Carlo (DSMC) method to resolve trace species such as electrons and ions in weakly ionized flows. The present study extends this scheme to reacting flows, utilizing a split-merge concept for particles with species-dependent weights. Reactive collisions between differently weighted species are processed according to their respective weights, introducing the production of an integer number of product particles and a probabilistic approach applied to both reactant and product particles to satisfy mass conservation. The extended scheme was validated in single-cell dissociation and exchange reaction test cases, verifying mass and total energy conservation and the mole-fraction histories against conventional DSMC. This extension provides the electron number density required for predicting the onset of communication blackout at high altitudes during hypersonic reentry. The method was applied to the Orbital Re-entry Experiment (OREX) flight conditions from 105.0 km to 92.8 km altitude. Translational and internal temperatures showed good agreement with conventional DSMC results, validating the macroscopic flow-field computation. The scheme enabled continuous spatial distributions of trace electron number density, consistent with the communication blackout onset observed in the flight data. The computed ion current reproduced the overall trend of the in-flight electrostatic probe data. The influences of the ion recombination at the vehicle surface and of the uncertainty of the N + O associative-ionization rate were quantified, and the range of the computed currents encompassed the available measurements. The 105.0 km case extended the prediction beyond available measurements, consistent with the overall trend.
Versions
Related
- Reactive Trace Species Modeling with the Simulation Tool PICLas
- Simulation of quantum entangled states by classical probabilistic models
- High-order Magnus expansion for Hamiltonian simulation
- Robust high-order quantum simulation using finite-width pulses
- Spectroscopy of analogue black holes using simulation-based inference
BibTeX
@misc{morimoto2026simulation,
title = {{Simulation of Weakly Ionized Hypersonic Flows with Reactive Species Weighting Scheme in the Direct Simulation Monte Carlo Method}},
author = {Takato Morimoto and Virgile Charton and Eiichiro Yamaoka and Kiyoshi Kinefuchi},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.08476},
eprint = {2610.08476},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.08476},
}
RIS
TY - UNPB
TI - Simulation of Weakly Ionized Hypersonic Flows with Reactive Species Weighting Scheme in the Direct Simulation Monte Carlo Method
AU - Takato Morimoto
AU - Virgile Charton
AU - Eiichiro Yamaoka
AU - Kiyoshi Kinefuchi
PY - 2026
DA - 2026/10/06
PB - arXiv
DO - 10.48550/arxiv.2610.08476
UR - https://arxiv.org/abs/2610.08476
AB - The conservative species weighting scheme enables the direct simulation Monte Carlo (DSMC) method to resolve trace species such as electrons and ions in weakly ionized flows. The present study extends this scheme to reacting flows, utilizing a split-merge concept for particles with species-dependent weights. Reactive collisions between differently weighted species are processed according to their respective weights, introducing the production of an integer number of product particles and a probabilistic approach applied to both reactant and product particles to satisfy mass conservation. The extended scheme was validated in single-cell dissociation and exchange reaction test cases, verifying mass and total energy conservation and the mole-fraction histories against conventional DSMC. This extension provides the electron number density required for predicting the onset of communication blackout at high altitudes during hypersonic reentry. The method was applied to the Orbital Re-entry Experiment (OREX) flight conditions from 105.0 km to 92.8 km altitude. Translational and internal temperatures showed good agreement with conventional DSMC results, validating the macroscopic flow-field computation. The scheme enabled continuous spatial distributions of trace electron number density, consistent with the communication blackout onset observed in the flight data. The computed ion current reproduced the overall trend of the in-flight electrostatic probe data. The influences of the ion recombination at the vehicle surface and of the uncertainty of the N + O associative-ionization rate were quantified, and the range of the computed currents encompassed the available measurements. The 105.0 km case extended the prediction beyond available measurements, consistent with the overall trend.
ER -