Search NASA⌕ Search

SEARCH · Search NASA

Results for “multi-region 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.

491 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Orbital Processing of High-Quality CdTe Compound Semiconductors

CdZnTe crystals were grown in one-g and in micro-g for comparative analysis. The two micro-g crystals were grown in the Crystal Growth Furnace during the First United States Microgravity Laboratory mission (USML-1). The samples were analyzed for chemical homogeneity, structural perfection, and optoelectronic performance (infrared transmission). Fourier Transform Infrared (FTIR) transmission of both ground and flight materials showed that the infrared transmission was close to theoretical, 63% versus 66%, suggesting that the material was close to the stochiometric composition during both the ground and flight experiments. Infrared microscopy confirmed that the principal precipitates were Te and their size (1-10 microns) and density suggested that the primary flight and ground base samples experienced similar cooling rates. Macrosegregation was predicted, using scaling analysis, to be low even in one-g crystals and this was confirmed experimentally, with nearly diffusion controlled growth achieved even in the partial mixing regime on the ground. Radial segregation was monitored in the flight samples and was found to vary with fraction solidified, but was disturbed due to the asymmetric grvitational and thermal fields experienced by the flight samples. The flight samples, however, were found to be much higher in structural perfection than the ground samples produced in the same furnace under identical growth conditions except for the gravitational level. Rocking curve widths were found to be substantially reduced, from 20/35 (one-g) to 9/20 (micro-g) for the best regions of the crystals. The full width at half maximum (FWHM) of 9 arc seconds is as good as the best reported terrestrially for this material. The ground samples were found to have a fully developed mosaic structure consisting of subgrains, whereas the flight sample dislocations were discrete and no mosaic substructure was evident. The defect density was reduced from 50-100,000 (one-g) to 500-25000 EPD (micro-g). These results were confirmed using rocking curve analysis, synchrotron topography, and etch pit analysis. The low dislocation density is thought to have resulted from the near-absence of hydrostatic pressure which allowed the melt to solidify with minimum or no wall contact, resulting in very low stress being exerted on the crystal during growth or during post-solidification cooling.

D J Larson, Jr.↗
NASA NTRS2026

Developing a Pyrolysis Gas Thermal Blocking Model for Reentry Demise

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These heating coefficients were derived for typical metallic materials that melt and do not have a strong gas-phase contribution to the flow in the boundary layer. Modern satellites typically feature fiber-reinforced polymer (FRP) components, such as solar array booms, facesheets of sandwich panels, or overwraps for composite-overwrapped pressure vessels (COPV). These FRP materials do not behave the same as metals in the reentry environment, but instead will pyrolyze and develop significant volumes of gas into the boundary layer. Accurately predicting the reentry demise of FRP components is critical to assessing the reentry casualty risk for modern spacecraft. Research in recent years has shown that this demisability can depend heavily on how the expulsion of gaseous pyrolysis products through the outer surface of the material affects the heat flux at the surface. The ODPO has been developing a reduced-order model of the effect of pyrolysis gas blowing on the heat flux based on correlations between a blowing factor and a non-dimensional heat flux to be incorporated in the upcoming version 7.3 of the Object Reentry Survivability Analysis Tool (ORSAT). This presentation discusses the progress of this development project and the challenges remaining for generalizing the model across families of FRP materials.

Benton Greene↗
NASA NTRSDate not supplied

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗
NASA NTRS2023

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig↗
NASA NTRS2026

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↗
NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗
NASA NTRS2026

Modeling Study of Hatch Spacing’s Effect on Grain Morphology in Repairing Damaged Metal Parts with Welding and Hatch Spacing’s Potential Use in Lunar Exploration

The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.

Welding for repairing↗
NASA NTRS2019

Physics-Based Modeling and Simulation of Emerging Battery Technologies for Aerospace

Recently there is a growing interest in the aviation sector to reduce air and noise pollution. Electrochemical energy storage devices such as batteries coupled with a distributed electric propulsion system can reduce noise concerns as well as emissions and allow the concept of Urban Air Mobility to come to fruition. The battery performance needs to improve considerably compared from current state-of-art Li-ion battery (about 200Wh/Kg) to realize all-electric passenger jet for short flights (up to 690 miles). NASA is exploring lithium-oxygen battery chemistry to power hybrid (battery powered electrical system) and all-electric aircraft for short distance and long-distance flights. Li-O2 is one of the advanced Li-ion technologies that promise to provide specific energy of more than 750Wh/Kg. For this presentation, we present our work on improving power density of Li-O2 batteries through the use of multiphysics simulations. Next, a path is outlined to port these model to simulate performance for a new battery chemistry for space application, Li-CO2. Li-CO2 uses carbon dioxide as the active material instead of oxygen. Although this technology is in its early development, the offers two benefits: it can be used as a CO2 scrubber, as oxygen is one of the by-products on charging, and as a backup or a standalone battery for various Mars or Venus missions, where the carbon dioxide content in the atmosphere is high and need battery to operate at higher temperatures.

Mehta, Mohit↗
NASA NTRS2000

The NASA Inductrack Model Rocket Launcher at the Lawrence Livermore National Laboratory

The Inductrack magnetic levitation system, developed at the Lawrence Livermore National Laboratory, is being studied for its possible use for launching rockets. Under NASA sponsorship, a small model system is being constructed at the Laboratory to pursue key technical aspects of this proposed application. The Inductrack is a passive magnetic levitation system employing special arrays of high-field permanent magnets (Halbach arrays) on the levitating carrier, moving above a "track" consisting of a close-packed array of shorted coils with which are interleaved with special drive coils. Halbach arrays produce a strong spatially periodic magnetic field on the front surface of the arrays, while canceling the field on their back surface. Relative motion between the Halbach arrays and the track coils induces currents in those coils. These currents levitate the carrier cart by interacting with the horizontal component of the magnetic field. Pulsed currents in the drive coils, synchronized with the motion of the carrier, interact with the vertical component of the magnetic field to provide acceleration forces. Motional stability, including resistance to both vertical and lateral aerodynamic forces, is provided by having Halbach arrays that interact with both the upper and the lower sides of the track coils. In its completed form the model system that is under construction will have a track approximately 100 meters in length along which the carrier cart will be propelled up to peak speeds of Mach 0.4 to 0.5 before being decelerated. Preliminary studies of the parameters of a full-scale system have also been made. These studies address the problems of scale-up, including means to simplify the track construction and to reduce the cost of the pulsed-power systems needed for propulsion.

L S Tung↗
NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗
NASA NTRSDate not supplied

Final Thermal Design and Thermal Vacuum Testing of the StarBurst Instrument

The StarBurst Multimessenger Pioneer is a small satellite mission serving as a wide-field gamma-ray observatory designed to capture the initial emissions of short gamma-ray bursts, electromagnetic signatures of neutron star mergers. This paper presents the final thermal design and analysis of the StarBurst Instrument, comprising the bus-to-instrument interface plate, control electronics, and twelve crystal detector units, which form the core of the mission’s science capability. The passive thermal control system design requires consideration of restrictive keep-out zones, unknown orbital parameters, and narrow temperature limits of the detectors. Also summarized is the instrument level thermal vacuum cycle test, correlated model refinements, and updated model results. Following successful completion of the instrument test campaign, the hardware was integrated with the spacecraft bus for spacecraft level testing, including additional thermal vacuum testing. The results from the spacecraft level thermal vacuum test will further inform the instrument thermal model, ensuring accurate flight temperature predictions. StarBurst launches as a secondary payload in 2027 and has a mission duration of at least one year.

StarBurst↗
NASA NTRS1976

NASA's Advanced Control Law Program for the F-8 Digital Fly-By-Wire Aircraft

This paper briefly describes the NASA F-8 Digital Fly-By-Wire (DFBW) program and Langley Research Center's role in investigating and promoting advanced control laws for possible flight experimentation. The primary purpose of the paper is to provide background material for papers describing some of the F-8 DFBW adaptive control law studies. The paper provides a brief description of the Phase II DFBW F-8 aircraft and its control system, discusses some of the advanced control law study objectives and guidelines, and provides some mathematical models which are useful in the control analysis problem.

Optimal Control↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
NASA NTRS2003

Optimal Battery Charging for Damage Mitigation

Our control philosophy is to charge the NiH2 cell in such a way that the damage incurred during the charging period is minimized, thus extending its cycle life. This requires nonlinear dynamic model of NiH2 cell and a damage rate model. We must do this first. This control philosophy is generally considered damage mitigating control or life-extending control. This presentation covers how NiH2 cells function, electrode behavior, an essentialized model, damage mechanisms for NiH2 batteries, battery continuum damage modeling, and battery life models. The presentation includes graphs and a chart illustrating how charging a NiH2 battery with different voltages and currents affects damages the battery and affects its life. The presentation concludes with diagrams of control system architectures for tracking battery recharging.

Hartley, Tom T.↗
NASA NTRSDate not supplied

SR-1 Freedom Thermal Architecture Challenges

The SR-1 Freedom mission aims to demonstrate nuclear electric propulsion (NEP) and deliver the Skyfall helicopter payload to Mars. Repurposing the existing Gateway Power and Propulsion Element (PPE) spacecraft and combining it with a 20 kWe-class nuclear power module (NPM) presents a variety of thermal architecture challenges. These include unprecedented waste heat rejection requirements, integration with a spacecraft bus originally designed for a different mission profile, and novel packaging constraints. This presentation will describe the simplified concept of operations as it relates to on-orbit thermal environments, early-phase architecture trades and supporting analyses, and planned forward work. We will also discuss our strategy for integrated thermal modeling of the complete spacecraft and examples of interface challenges necessitated by this ambitious effort.

Thermal↗
NASA NTRSDate not supplied

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗
NASA NTRSDate not supplied

Predicting Team Functioning in Long Term Space Missions Using Acoustic and Linguistic Measures

Maintaining optimal team functioning is critical for long-duration space exploration missions, yet traditional monitoring methods, such as self-reports and wearable sensors, often impose operational burdens or suffer from bias. This paper investigates a non-intrusive speech-based artificial intelligence (AI) framework to predict degradations in team functioning using data from the Human Exploration Research Analog (HERA) of the U.S. National Aeronautics and Space Administration (NASA). Using acoustic features, linguistic descriptors, and semantic embeddings, we evaluate static non-linear and temporal machine learning models to predict both objective (task accuracy) and subjective (self-reported efficacy and cohesion) team functioning outcomes. Results indicate that temporal models outperform static approaches, with prediction of objective task accuracy in Team Interaction Battery (TIB) improving from near chance to 71%. Self-reported outcomes, including team efficacy and cohesion, are predicted more reliably than task performance, achieving balanced accuracies of up to 85.56% and 78.12%, respectively, and are found to be most strongly associated with acoustic features. In a second interdependent task, the MMSEV–EVA, accuracies of up to 78% are achieved using temporal models with acoustic features. Furthermore, incorporating just 1–2 days of team-specific historical data systematically improved performance, and acoustic markers from informal pre-task interactions provided modest predictive gains. Finally, while automated preprocessing yielded viable accuracy, humancorrected data provided moderate performance gains, though transcription error rates did not significantly correlate with model performance. These findings highlight the potential of speech as a passive, high-fidelity monitoring tool for autonomous habitats.

Temporal modeling↗
NASA NTRSDate not supplied

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↗