Search NASASearch

SEARCH · Search NASA

Results for “Transient Dynamic Analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,052 records · Page 14

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics

Dissipative Phase Transition in the Two-Photon Dicke Model

We explore the dissipative phase transition of the two-photon Dicke model, a topic that has garnered significant attention recently. Our analysis reveals that while single-photon loss does not stabilize the intrinsic instability in the model, the inclusion of two-photon loss restores stability, leading to the emergence of superradiant states, which coexist with the normal vacuum states. Using a second-order cumulant expansion for the photons, we derive an analytical description of the system in the thermodynamic limit, which agrees well with the exact calculation results. Additionally, we present the Wigner function for the system, shedding light on the breaking of the 𝑍4 symmetry inherent in the model. These findings offer valuable insights into stabilization mechanisms in open quantum systems and pave the way for exploring complex nonlinear dynamics in two-photon Dicke models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Atomistic Mechanisms of Stress-Dependent Molten Salt Corrosion in NiCr Alloys

Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.

36 - MATERIALS SCIENCE

The Dark Energy Survey supernova program: a reanalysis of cosmology results and evidence for evolving dark energy with an updated Type Ia supernova calibration

We present improved cosmological constraints from a re-analysis of the Dark Energy Survey (DES) 5-year sample of Type Ia supernovae (DES-SN5YR). This re-analysis includes an improved photometric cross-calibration, recent white dwarf observations to cross-calibrate between DES and low-redshift surveys, retraining the salt3 light-curve model and fixing a numerical approximation in the host-galaxy colour law. Our fully recalibrated sample, which we call DES-Dovekie, comprises ~1600 likely Type Ia SNe from DES and ~200 low-redshift SNe from other surveys. With DES-Dovekie, we obtain Ω m = 0.330 ± 0.015 in flat Lambda-cold dark matter (⁠ΛCDM) which changes Ω m by –0.022 compared to DES-SN5YR. Combining DES-Dovekie with cosmic microwave background data from Planck, Atacama Cosmology Telescope, and South Pole Telescope and the DESI DR2 measurements in a flat CDM cosmology, we find ω 0 = –0.803 ± 0.054 and ω a = –0.72 ± 0.21⁠. Our results hold a significance of 3.2σ, reduced from 4.2σ for DES-SN5YR, to reject the null hypothesis that the data are compatible with the cosmological constant. This significance is equivalent to a Bayesian model preference odds of approximately 5:1 in favour of the flat ω 0 ω a CDM model. Using generally accepted thresholds for model preference, our updated data exhibits only a weak preference for evolving dark energy.

dark energy

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics

Terahertz conductivity of two-dimensional materials: a review

Two-dimensional (2D) van der Waals materials are shaping the landscape of next-generation devices, offering significant technological value thanks to their unique, tunable, and layer-dependent electronic and optoelectronic properties. Time-domain spectroscopic techniques at terahertz (THz) frequencies offer noninvasive, contact-free methods for characterizing the dynamics of carriers in 2D materials. They also pave the path toward the applications of 2D materials in detection, imaging, manufacturing, and communication within the increasingly important THz frequency range. In this paper, we overview the synthesis of 2D materials and the prominent THz spectroscopy techniques: THz time-domain spectroscopy, optical-pump THz-probe technique, and optical pump–probe THz spectroscopy. Through a confluence of experimental findings, numerical simulation, and theoretical analysis, we present the current understanding of the rich ultrafast physics of technologically significant 2D materials: graphene, transition metal dichalcogenides, MXenes, perovskites, topological 2D materials, and 2D heterostructures. Finally, we offer a perspective on the role of THz characterization in guiding future research and in the quest for ideal 2D materials for new applications.

2D materials

From Exploration Flight Test-1 to Artemis II--A NASA Langley's Orion Aerosciences Overview

The Orion Aerosciences program at NASA Langley has played a central role in advancing the aerodynamic and aeroheating prediction capabilities required for the Orion crew vehicle’s return from deep space. This presentation provides a technical overview of aerosciences contributions spanning Exploration Flight Test-1 (EFT-1), Artemis I, and the ongoing post-flight analysis of Artemis II. EFT-1 provided the first high-energy entry dataset for Orion, enabling critical validation of aerodynamic force and moment predictions, static and dynamic stability characteristics, and aeroheating environments at relevant flight Mach and Reynolds numbers. Flight-derived pressure data were used to refine the Flush Air Data System (FADS) methodology for atmospheric density reconstruction and to improve Best Estimated Trajectory (BET) solutions. The EFT-1 data also offered key insights into heat shield performance, including char layer recession, in-depth thermal response, and material retention behavior under flight conditions, informing updates to both thermal response models and uncertainty quantification practices. Building on EFT-1, Artemis I extended the database to true lunar-return conditions. Observations of heat shield performance, including localized char loss, bondline response, and recession variability, provided an unprecedented opportunity to reassess Thermal Protection System (TPS) and aeroheating modeling assumptions. Aerodynamic reconstruction efforts incorporated improved FADS calibration, enhanced atmospheric modeling, and refined force and moment databases to reduce trajectory and load uncertainties. Aeroheating comparisons between pre-flight predictions and flight data enabled targeted model updates, particularly in transitional flow environments and wake heating regions. For Artemis II, these lessons were systematically incorporated into the pre-flight prediction process. Updates included refined aerodynamic databases anchored to flight-validated corrections, improved density estimation and BET methodologies using enhanced database interpolation algorithm and FADS modeling, and revised aeroheating design environments informed by Artemis I material response observations. By the time of the workshop, Artemis II post-flight analysis will be underway, and preliminary findings will be presented where available, including early comparisons of aerodynamic reconstruction, atmospheric density estimation, and thermal protection system performance relative to updated predictions. Collectively, this body of work is a testament to the dedicated and multidisciplinary team whose sustained efforts have contributed to the program’s success and to the progressive maturation of Orion aerosciences modeling through numerical modeling, ground and flight data assimilation. The integrated advancement of aerodynamics, trajectory reconstruction, FADS-based density estimation, and aeroheating analysis has reduced predictive uncertainty and strengthened confidence for future crewed lunar and deep-space missions.

Orion

Robust and tunable oxide nanoscrolls for solar-driven H 2 generation and storage

Hydrogen gas is a promising alternative to fossil fuels due to its high energy output and environmentally safe byproducts. Various morphologies of photocatalytic materials have been explored for high-efficiency H 2 production, for instance, quasi-1D nanoscroll structures that provide a larger surface-to-volume ratio. Recently, we predicted layer-by-layer formation of stable oxide nanoscrolls directly from dichalcogenide precursors, eliminating the need for costly formation of two-dimensional oxides for a roll-up synthesis of nanoscrolls. Here, in this study, we evaluate the suitability of those oxide nanoscroll materials—MoO 3 , WO 3 , PdO 2 , HfO 2 , and GeO 2 —for solar-driven photocatalytic H 2 production and storage. Using excited state theory coupled with Bethe–Salpeter equation simulations, we discern their electronic and optical properties as a function of interlayer scroll spacing and find them to be highly conducive for solar-driven photocatalysis. Additionally, using ab initio molecular dynamics simulations, we show that they are also suitable for H 2 storage as the nanoscrolls exhibit an effective trapping of hydrogen, even in the presence of defects and vacancies in the oxides. This work thus demonstrates the discovery of robust and tunable oxide nanoscrolls as materials for advancing solar-driven hydrogen technologies.

Chemical compounds

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE

Lunar Electrical Conductivity, Permeability and Temperature From Apollo Magnetometer Experiments

Magnetometers have been deployed at four Apollo sites on the Moon to measure remanent and induced lunar magnetic fields. Measurements from this network of instruments have been used to calculate the electrical conductivity, temperature, magnetic permeability, and iron abundance of the lunar interior. The measured lunar remanent fields range from 3 gammas (γ) minimum at the Apollo 15 site to 327 γ maximum at the Apollo 16 site. Simultaneous magnetic field and solar plasma pressure measurements show that the remanent fields at the Apollo 12 and 16 sites interact with, and are compressed by, the solar wind. Remanent fields at Apollo 12 and Apollo 16 are increased 16 γ and 32 γ, respectively, by a solar plasma bulk pressure increase of 1.5 x 10 -7 dynes/cm 2 . Global lunar fields due to eddy currents, induced in the lunar interior by magnetic transients, have been analyzed to calculate an electrical conductivity profile for the Moon. From nightside magnetometer data in the solar wind it has been found that deeper than 170 km into the Moon the conductivity rises from 3 x 10 -4 mhos/m to 10 -2 mhos/m at 1000-km depth. Recent analysis of data obtained in the geomagnetic tail, in regions free of complicating plasma effects, yields results consistent with nightside values. Conductivity profiles have been used to calculate the lunar temperature for an assumed lunar material of olivine. In the outer layer (~ 170 km thick) the temperature rises to 110° C, after which it gradually increases with depth to 1500° C at a depth of ~ 1000 km. Simultaneous measurements by magnetometers on the lunar surface and in orbit around the moon are used to construct a whole-moon hysteresis curve, from which the global lunar magnetic permeability is determined to be μ = 1.012 ± 0.006. The corresponding global induced dipole moment is 2 x 10 18 gauss-cm 3 for typical inducing fields of 10 -4 gauss in the lunar environment. Lunar free iron abundance corresponding to the global permeability is determined to be 2.5 ± 2.0 wt. %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior. For a free iron/orthopyroxene lunar composition the total iron content is calculated to be 12.8 ± 1.0 wt. %; for a free iron/olivine composition, total iron content is 5.5 ± 1.2 wt. %. Other lunar models with an iron core and with a shallow iron-rich layer are also discussed in light of the measured global lunar permeability. Velocities and thicknesses of the Earth's magnetopause and bow shock have been estimated from simultaneous magnetometer measurements. Average speeds are determined to be about 50 km/s for the magnetopause and 70 km/s for the bow shock, although there are large variations in the measurement for any particular boundary crossing. Corresponding measured boundary thicknesses average about 2300 km for the magnetopause and 1400 km for the bow shock.

P Dyal

Model Spread in Tropical Low Cloud Feedback Tied to Overturning Circulation Response to Warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

Atmospheric dynamics

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin

Actuating Liquid Crystals Rapidly and Reversibly by Using Chemical Catalysis

Abstract Microtubules and catalytic motor proteins underlie the microscale actuation of living materials, and they have been used in reconstituted systems to harness chemical energy to drive new states of organization of soft matter (e.g., liquid crystals (LCs)). Such materials, however, are fragile and challenging to translate to technological contexts. Rapid (sub‐second) and reversible changes in the orientations of LCs at room temperature using reactions between gaseous hydrogen and oxygen that are catalyzed by Pd/Au surfaces are reported. Surface chemical analysis and computational chemistry studies confirm that dissociative adsorption of H 2 on the Pd/Au films reduces preadsorbed O and generates 1 ML of adsorbed H, driving nitrile‐containing LCs from a perpendicular to a planar orientation. Subsequent exposure to O 2 leads to oxidation of the adsorbed H, reformation of adsorbed O on the Pd/Au surface, and a return of the LC to its initial orientation. The roles of surface composition and reaction kinetics in determining the LC dynamics are described along with a proof‐of‐concept demonstration of microactuation of beads. These results provide fresh ideas for utilizing chemical energy and catalysis to reversibly actuate functional LCs on the microscale.

Chemistry

Introduction to Radiation Issues for International Space Station Extravehicular Activities: Chapter 1

The International Space Station (ISS) provides significant challenges for radiation protection of the crew due to a combination of circumstances including: the extended duration of missions for many crewmembers, the exceptionally dynamic nature of the radiation environment in ISS orbit, and the necessity for numerous planned extravehicular activities (EVA) for station construction and maintenance. Radiation protection requires accurate radiation dose measurements and precise risk modeling of the transmission of high fluxes of energetic electrons and protons through the relatively thin shielding provided by the space suits worn during EVA. Experiments and analyses have been performed due to the necessity to assure complete radiation safety for the EVA crew and thereby ensure mission success. The detailed characterization described of the material and topological properties of the ISS space suits can be used as a basis for design of space suits used in future exploration missions. In radiation protection practices, risk from exposure to ionizing radiation is determined analytically by the level of exposure, the detrimental quality of the radiation field, the inherent radiosensitivity of the tissues or organs irradiated, and the age and gender of the person at the time of exposure. During low Earth orbit (LEO) EVA, the relatively high fluxes of low-energy electrons and protons lead to large variations in exposure of the skin, lens of the eye, and tissues in other shallow anatomical locations. The technical papers in this publication describe a number of ground-based experiments that precisely measure the thickness of the NASA extravehicular mobility unit (EMU) and Russian Zvezda Orlan-M suits using medical computerized tomography (CT) X-ray analysis, and particle accelerator experiments that measure the minimum kinetic energy required by electrons and photons to penetrate major components of the suits. These studies provide information necessary for improving the understanding of the current ISS space suits and provide insights into improved approaches for the design of future suits. This chapter begins with a summary of the dynamic ionizing radiation environment in LEO space and introduces the concepts and quantities used to quantify exposure to space radiation in LEO. The space suits used for EVA and the experimental partial human phantom are described. Subsequent chapters report results from measured charged particle fields before and after incident protons and secondary particles are transported through the space suits and into organs and tissues.

M R Shavers

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen

Photophysical and Time‐resolved Infrared Properties of Long‐Lived Rhenium(I) 4,5‐Diazafluorene Tricarbonyl Chromophores

This report investigates the synthesis, structural characterization, fundamental molecular photophysics, electrochemistry, UV-Vis spectroelectrochemistry, and time-resolved infrared spectroscopic properties of eight [fac-Re(dafR)(CO) 3 L] 0/+ complexes, where R=ethyl [(dedaf); 1, 3, 5, 7] or H [(dafH); 2, 4, 6, 8] and L=Cl− (1, 2), imidazole [(Im); 3, 4], 4-ethylpyridine [(4-Etpy); 5, 6], or pyridine [(py); 7, 8]. Universally, 1–8 yield higher energy photoluminescence (PL) emission bands and higher PL quantum yields (up to 53 %) than the classic 2,2’-bipyridine (bpy) and 1,10-phenanthroline (phen) ligated Re(I) tricarbonyl complexes. The excited state lifetimes of 1–8 lie between those corresponding to the bpy and phen derivatives, ranging from 120 and 1300 ns at room temperature. Combinations of reductive UV-Vis spectroelectrochemistry, transient absorption spectroscopy, and time-resolved infrared spectroscopy consistently assigned the lowest excited states in 1–8 being of metal-to-ligand charge transfer (MLCT) character. These new ReI MLCT chromophores follow classic energy gap law behavior and possess the characteristics necessary for serving as valuable photosensitizers suitable to energize excited state electron and energy transfer photochemistry.

Charge transfer

Properties of a Ni19.5Pd30Ti50.5 High-Temperature Shape Memory Alloy in Tension and Compression

Potential applications involving high-temperature shape memory alloys have been growing in recent years. Even in those cases where promising new alloys have been identified, the knowledge base for such materials contains gaps crucial to their maturation and implementation in actuator and other applications. We begin to address this issue by characterizing the mechanical behavior of a Ni 19.5 Pd 30 Ti 50.5 high-temperature shape memory alloy in both uniaxial tension and compression at various temperatures. Differences in the isothermal uniaxial deformation behavior were most notable at test temperatures below the martensite finish temperature. The elastic modulus of the material was very dependent on strain level; therefore, dynamic Young’s Modulus was determined as a function of temperature by an impulse excitation technique. More importantly, the performance of a thermally activated actuator material is dependent on the work output of the alloy. Consequently, the strain-temperature response of the Ni 19.5 Pd 30 Ti 50.5 alloy under various loads was determined in both tension and compression and the specific work output calculated and compared in both loading conditions. It was found that the transformation strain and thus, the specific work output were similar regardless of the loading condition. Also, in both tension and compression, the strain-temperature loops determined under constant load conditions did not close due to the fact that the transformation strain during cooling was always larger than the transformation strain during heating. This was apparently the result of permanent plastic deformation of the martensite phase with each cycle. Consequently, before this alloy can be used under cyclic actuation conditions, modification of the microstructure or composition would be required to increase the resistance of the alloy to plastic deformation by slip.

High-Temperature Shape Memory Alloy

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics