Search NASASearch

SEARCH · Search NASA

Results for “Human Thermal Modeling”

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,521 records · Page 14

Extrusion-based Additive Manufacturing of Regolith-Filled Shape Memory Vitrimer Composite for Lunar Construction

The National Aeronautics and Space Administration (NASA) is visiting the moon again. This time, the objective is to explore establishing a permanent lunar base. To achieve both longterm human habitation on the moon and future deep space travel, it is crucial to make the most of the in-situ resources and build autonomous systems on the moon to support the construction of a lunar habitat. NASA’s In-situ Resource Utilization (ISRU) program aims to minimize the need to ship heavy prefabricated structures, reducing cost and enhancing sustainability. Here, we developed an economical extrusion method for printing lunar regolith-based composites using shape memory vitrimer as a binder. A rheological study is conducted to determine the extrudability of the composite with different regolith weight percentages. Several characterizations were conducted on the composites. The as-printed composites exhibited compressive and flexural strengths of 73.32 MPa and 156.59 MPa, respectively, and good impact tolerance. The composite maintained 57.92% of its mechanical properties even after the second crack healing cycle. The composites also exhibit shape fixity ratio of 90.02% and shape recovery ratio of 83.46%. The simple synthesis method, sustainability, and good thermomechanical properties make the 3D printed composite an ideal material for lunar construction applications.

Kingsley Yeboah Gyabaah

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling

Coexistence of Synchronization and Stochasticity in Thermally Coupled Mott Oscillators

Synchronization is conventionally regarded as a mechanism for suppressing variability and enforcing order in coupled systems, from pendula and lasers to neurons and electronic oscillators. Here, we show that synchronization can also embed stochasticity at finer scales. We observe this phenomenon in thermally coupled VO 2 neuristors, where robust in-phase synchronization at the microsecond scale coexists with spike onset fluctuations at the nanosecond scale, with no fixed leader. The coexistence of order and disorder originates from stochastic domain-level physics of the insulator–metal and metal–insulator transitions, where local variations in transition temperature drive cycle-to-cycle randomness in nucleation, percolation, and relaxation. A stochastic domain model reproduces this effect by generating synchronized spike trains with random lead–lag jitter, and experimental interspike interval statistics confirm the persistence of fine-scale variability despite macroscopic phase locking. These findings establish that synchronization and stochasticity can coexist within the same physical platform, revealing hidden disorder within collective order. Furthermore, this insight reframes synchronization as not purely deterministic, but as a universal context where microscopic variability can persist, with implications for electronics, cryptography, and the fundamental physics of order–disorder coexistence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

On-Orbit Calibration Assessment of NOAA-21 VIIRS Thermal Emissive Bands and Implications for JPSS-4 VIIRS

The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the NOAA-21 (N21) satellite has been successfully operating for over three years collecting valuable scientific measurements. A large suite of weather forecasting models and scientific research applications are supported with the VIIRS measurements from N21 combined with those acquired from the VIIRS instrument onboard the S-NPP and NOAA-20 spacecraft. Among the 22 VIIRS spectral bands, there are seven thermal emissive bands (TEB) covering the 3.7 to 12.2 µm spectral range at two different spatial resolutions. The VIIRS TEB detectors are calibrated onboard using a blackbody at a controlled temperature and a space view for background signal measurement. In this paper, we discuss the on-orbit performance of the N21 VIIRS TEB using various parameters such as the detector gain, noise, and offset. This on-orbit assessment provides critical insights into instrument behavior and calibration performance that will inform pre-launch testing and post-launch validation strategies for the upcoming JPSS-4 VIIRS mission scheduled for launch in 2027. The methodologies and lessons learned from characterizing the N21 TEB performance will enable more efficient and comprehensive radiometric assessment of the next VIIRS instrument. Specific focus is placed on identifying performance trends and anomaly signatures observed during the N21 commissioning phase that can enhance the JPSS4 instrument checkout procedures and accelerate the transition to routine on-orbit operations. The established performance baselines and uncertainty estimates will serve as an acceptance criterion for the JPSS-4 commissioning activities, ensuring mission readiness and data continuity for operational users.

Amit Angal

Development of a silver-zinc battery system

Summary report is described of historical documentation and detailed design data for development of silver-zinc battery for use on Surveyor spacecraft. Electrical and physical characteristics of battery models are included, along with data on qualification, acceptance, solar-thermal-vacuum, mission simulation testing, and actual flight performance.

Moses, A. J.

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

05 NUCLEAR FUELS

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles with Interrupted Relaxation Periods of TIMETAL 21S at 650 °C

The predictive performance of a multimechanism GVIPS model with saturating hardening function is assessed against a set of high temperature cyclic experiments under both strain-control and stress-control at varying rates of loading as well as several complex interrupted cyclic responses. The material specimens were composed of the titanium alloy, TIMETAL 21S, and tested at 650 °C, see Lissenden et al., 2007. The model did a very good job of predicting the variety of tests, particularly given the fact that only monotonic loading tests (tensile, creep, relaxation) and a single fully reversed cycle was used for characterization. It was particularly interesting that the model was capable of reasonably capturing the rate of accumulated strain during ratchetting under tensile mean stress. Given complex interrupted cyclic response, the model was able to both qualitatively and even quantitatively predict reasonably well both the cyclic and associated relaxation periods (in all quadrants of the stress-strain space) thus confirming the validity of the chosen functional forms for both hardening and thermal recovery. The quantitative inaccuracy is associated with the fact that the material specimens tested by Lissenden et al., 2007 was significantly “softer” than those tested by Castelli in Saleeb et al., 1994 which were used for characterization of the GVIPS model parameters.

Viscoplaticity

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Prime Time for Model-Predictive Control? Assessing the Technical and Market Readiness of Advanced Controls in Buildings

Despite three decades of extensive research and field testing that have consistently validated the benefits of Model Predictive Control (MPC) in building applications, the technology has seen limited market adoption. This paper evaluates the readiness of MPC for widespread deployment, showcases recent demonstrations and field tests across diverse building types, including residential, small commercial, large commercial, and campus settings. Our results demonstrate that MPC can optimize system operations to achieve load shifting, minimize curtailment of on-site generation, and reduce energy costs by up to 80 %, while maintaining or improving occupant comfort. We also show that MPC can effectively control large assets, such as MW-sized thermal storage systems, and respond to dynamic pricing signals. However, achieving scale remains difficult due to labor-intensive workflows, reliance on a “PhD-in-the-loop” for MPC design and maintenance, susceptibility to fragile data infrastructure, and persistent workforce education and acceptance barriers. To bridge this gap, we outline a transition from bespoke, labor intensive prototypes toward streamlined, segment-targeted deployment strategies that leverage model templates, semantic tools, and generative AI. By automating control configuration and reducing engineering effort, these recommendations provide a pathway for transforming successful research demonstrations into scalable, market ready solutions for MPC-based controls.

Pritoni, Marco

The Alpha-Proton-X-Ray Spectrometer Deployment Mechanism: an Anthropomorphic Approach to Sensor Placement on Martian Rocks and Soil

On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Richard S Blomquist

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan

Characterizing and modeling the mechanical behavior of an anion exchange membrane for carbon capture applications

A new direct air capture (DAC) technology uses a moisture swing (MS) process with anion exchange membranes, potentially offering a more energy-efficient way to remove CO₂ from the air. In this MS process, the membrane absorbs CO₂ as it dries and releases it when water is added. Understanding the mechanical behavior of these membranes is essential for improving the design and efficiency of DAC systems and prolonging sorbent lifetime. This study tested one anion exchange membrane, Fumasep’s FAA-3, under mechanical loading and various temperature and humidity conditions to measure its swelling, stiffness, strength, plastic deformation, and stress relaxation. Experimental results were used to identify a mechanical model for FAA-3 that can be used to predict the material’s nonlinear viscous behavior under various loads and environments. Unlike prior studies that assumed linear elastic behavior, this work incorporates humidity-dependent swelling, thermal expansion, and nonlinear viscoelasticity of FAA-3 in both the experiments and the model.

Direct air capture

Large Language Model for Validation, Optical Calibration, and Learning (VOCAL) Distributed Temperature Sensing Interface

Distributed temperature sensing (DTS) using fiber optic sensors (FOS) offers a promising method for temperature measurements in advanced reactors, such as sodium fast reactors and molten salt cooled reactors. To support the calibration and validation of DTS measurements, Argonne National Laboratory developed the Validation, Optical Calibration, and Learning (VOCAL) software package. This report describes the integration of a local large language model (LLM) with a retrieval-augmented generation (RAG) system into the VOCAL interface to serve as an interactive user assistant. The LLM framework enhances the VOCAL platform’s accessibility to users by explaining interface components, clarifying inputs and outputs, and answering user queries dynamically in real-time. The accuracy of the LLM assistant performance was evaluated with 20 queries regarding the interface and its parameters using experimental data from the Thermal Hydraulic Experimental Test Article (THETA) facility. Results demonstrate that the LLM achieved a 95% accuracy rate, with a BERTScore of 0.8816 and SBERT value of 0.7417. Furthermore, validation of the RAG system within the LLM framework showed optimal accuracy with k-values between 1 and 2 using the k-refinement convergence test. The prompt perturbation analysis demonstrated good initial consistency for the RAG system, exhibiting the highest accuracy under punctuation variations and the greatest sensitivity under query reordering. Notably, the model’s errors were limited to data retrieval failures rather than factual hallucinations, reinforcing its baseline reliability. The integration of LLM provides a highly accurate, userfriendly enhancement to the VOCAL platform without disrupting its core computational capabilities for FOS calibration and validation.

Hong, Evan

Controls for Electrified Aircraft Propulsion

Advanced electrified aircraft propulsion (EAP) concepts with integrated power, propulsion, and thermal systems require the development of equally sophisticated controllers to fly safe and efficient missions. Hybrid electric aircraft utilize electric machines mechanically coupled to the engine shafts to extract and insert power for a variety of purposes. Additionally, electric machines may drive propulsive fans pulling from a combination of on-board energy storage devices and engine extracted power. NASA’s investment in hybrid electric aircraft hardware-in-the-loop testing enables controls research on representative engine models using a novel emulation and scaling methodology. High-voltage, high-power electrical powertrain presents several technical risks at altitude. High-power requires high efficiency to minimize losses. Superconducting electric machines and power distribution research at NASA has identified key challenges with potential controls solutions. These risks necessitate the development of controllers robust to model uncertainty, disturbances, and fault conditions. System health management and fault detection schemes play an important role in a multi-layer controls approach that utilizes a supervisor which oversees inner loops for the highly coupled power, propulsion, and thermal systems. Such a supervisory controller enables coordinates energy transfer between engines, energy storage devices, electric machines, propulsors, and heat exchangers. Existing methods have shown an improvement in engine operability using the hybrid electric powertrain.

Halle E Buescher

Characterization of the Mechanical Trash Compactor Bag Removal Force

Managing trash accumulation is a critical challenge in human space explorations since the trash volume and mass can significantly influence launch cost and operational efficiency. This study aids the development of a Mechanical Trash Compactor (Mpactor) by testing the removal force required to remove a “trash bag” from the Mpactor. To do this, a Pneumatic Pressure Mandrel System was designed to simulate the normal force trash would exerts on the trash bag and subsequently the sidewalls of the Mpactor. For various levels of sidewall normal force, the removal force required to remove the simulated trash bag from the Mpactor compaction chamber was measured. This testing is used to identify optimal bag materials. Additionally, testing was done on both the native stainless steel walls of the compaction chamber as well as the walls after being Teflon lined. Split-plot ANOVA, simple main effects analysis, Holm’s post hoc comparisons for interaction, and sensitivity analyses were conducted to assess the influence of trash bag and sidewall materials on the required removal force. It was found that Teflon reduced the force required to remove the trash bag as well as reduced material-dependent variation in removal force. Testing showed that Teflon reduced the required removal force with a Polyvinylidene Fluoride , an unidentified Ethylene acrylate - Low Density Polyethylene copolymer, and a Polyvinyl Chloride all requiring significantly less removal force than the other materials tested.

pressure mandrel

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

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