Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulation”

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.

At least 433 records · Page 24

The Use of Microgravity Simulators for Space Research

The spaceflight environment is known to influence biological processes ranging from stimulation of cellular metabolism to possible impacts on cellular damage repair, suppression of immune functions, and bone loss in astronauts. Microgravity is one of the most significant stress factors experienced by living organisms during spaceflight, and therefore, understanding cellular responses to altered gravity at the physiological and molecular level is critical for expanding our knowledge of life in space. Since opportunities to conduct experiments in space are scarce, various microgravity simulators and analogues have been widely used in space biology ground studies. Even though simulated microgravity conditions have produced some, but not all of the biological effects observed in the true microgravity environment, they provide test beds that are effective, affordable, and readily available to facilitate microgravity research. Kennedy Space Center (KSC) provides ground microgravity simulator support to offer a variety of microgravity simulators and platforms for Space Biology investigators. Assistance will be provided by both KSC and external experts in molecular biology, microgravity simulation, and engineering. Comparisons between the physical differences in microgravity simulators, examples of experiments using the simulators, and scientific questions regarding the use of microgravity simulators will be discussed.

Simulation↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

The formation of interstellar jets

The formation of interstellar jets by convergence of supersonic conical flows and the further dynamical evolution of these jets are investigated theoretically by means of numerical simulations. The results are presented in extensive graphs and characterized in detail. Strong radiative cooling is shown to result in jets with Mach numbers 2.5-29 propagating to lengths 50-100 times their original widths, with condensation of swept-up interstellar matter at Mach 5 or greater. The characteristics of so-called molecular outflows are well reproduced by the simulations of low-Mach-number and quasi-adiabatic jets.

Tenorio-Tagle, G.↗

Maps of Magnetic Field Strength in the OMC-1 Using HAWC+FIR Polarimetric Data

Far-infrared dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high-quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. In this work, the Davis–Chandrasekhar–Fermi technique is used to make maps of the plane-of-sky (POS) component of the magnetic field in the Orion Molecular Cloud (OMC-1) by combining polarization maps at 53, 89, 154 and 214 μm from HAWC+/SOFIA with maps of density and velocity dispersion. In addition, maps of the local dispersion of polarization angles are used in conjunction with Zeeman measurements to estimate a map of the strength of the line-of-sight (LOS) component of the field. Combining these maps, information about the threedimensional magnetic field configuration (integrated along the LOS) is inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG are observed near the BN/KL object, while the OMC-1 bar shows strengths of up to a few hundred μG. These estimates of the magnetic field components are used to produce maps of the mass-to-magnetic-flux ratio (M/Φ)—a metric for probing the conditions for star formation in molecular clouds— and determine regions of sub- and supercriticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds. Unified Astronomy Thesaurus concepts: Molecular clouds (1072); Giant molecular clouds (653); Interstellar magnetic fields (845); Far infrared astronomy (529)

Jordan A Guerra↗

Laboratory simulation of dust spectra

Laboratory studies of the IR spectra of interstellar dust are reviewed. Studies of the absorption spectra of dense molecular clouds are discussed, including methods to produce interstellar ice analogues, simulations of astronomical spectra, and IR absorption features caused by ices. Comparisons are made between observational and experimental results of interstellar dust studies. Also, the interstellar emission features associated with dusty regions exposed to UV radiation are examined, including bands related to PAHs and PAH-related materials. It is shown that interstellar spectra are more consistant with emission from free PAHs than with emission from particles.

Allamandola, L. J.↗

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗

Using Polarized PCRTM-Solar Model to Validate EMIT Surface Reflectance L2A Product

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer scattering code, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate hyperspectral radiation fields from visible to shortwave infrared (VSWIR), i.e., 250 nm to 2500 nm. Since only limited accurate line-by-line calculation using VLIDORT is needed, this coupled PCRTM-VLIDORT approach is very efficient for simulating hyperspectral VSWIR spectra. These simulations incorporate accurately defined optical properties (gas absorptions, molecular and aerosol/cloud particle scattering) and detailed surface reflectance inputs. In this study we used the current L2A products of the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as the inputs to PCRTM-VLIDORT to simulate EMIT spectra. Comparing the simulated spectra with EMIT measurements provides an efficient method to validate the current EMIT product.

Xiaozhen Xiong↗

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗

Simulations of Moon-Forming Impacts

We have incorporated a newly upgraded version of the equation of state ANEOS that allows for the formation of molecular vapor into a smoothed-particle hydrodynamics code. Results of simulations of Moon-forming impacts with this new EOS are presented. Additional information is contained in the original extended abstract.

Canup, R. M.↗

Synthesis of H2 in dirty ice mantles by fast ion energy loss: New experimental results increase the relevance of this mechanism

Recent experimental results support the importance of H2 production in molecular clouds by cosmic ray bombardment of the mantles of grains. The formation of molecules different from those originally present in the irradiated layer can be explained by the production of molecular fragments induced by the release of energy if the impinging fast particle. One way of considering the process is in terms of a transiently hot cylinder, initially about 50 A in diameter, that exists around the track of an individual fast ion. Since ice has a relatively low thermal conductivity, energy lost by the ion in the ice layers remains confined around the track for time long enough to be thermalized. The hot cylinder increases in diameter and decreases in temperature on a time scale of 10(exp -11) to 10(exp -10) sec. Molecular fragments that are formed in this high temperature region acquire enough mobility to recombine with different partners, forming new molecules. A Monte Carlo simulation of the interaction between cosmic rays and grain mantles, at various depths in the core of a spherical molecular cloud, was performed. The simulation was continued until 40,000 ions had hit each grain of the type and size chosen. During the performed experiments thin icy films made of H2O and CD4 mixed in the gas phase and deposited on a cold finger at 9 K were irradiated with 1.5 MeV helium beams. Among synthesized molecules were found H2, HD, and D2.

Pirranello, Valero↗

Space station molecular sieve development

An essential function of a space environmental control system is the removal of carbon dioxide (CO2) from the atmosphere to control the partial pressure of this gas at levels lower than 3 mm Hg. The use of regenerable solid adsorbents for this purpose was demonstrated effectively during the Skylab mission. Earlier sorbent systems used zeolite molecular sieves. The carbon molecular sieve is a hydrophobic adsorbent with excellent potential for space station application. Although carbon molecular sieves were synthesized and investigated, these sieves were designed to simulate the sieving properties of 5A zeolite and for O2/N2 separation. This program was designed to develop hydrophobic carbon molecular sieves for CO2 removal from a space station crew environment. It is a first phase effort involved in sorbent material development and in demonstrating the utility of such a material for CO2 removal on space stations. The sieve must incorporate the following requirements: it must be hydrophobic; it must have high dynamic capacity for carbon dioxide at the low partial pressure of the space station atmosphere; and it must be chemiclly stable and will not generate contaminants.

Chang, C.↗

Laser supported detonation wave source of atomic oxygen for aerospace material testing

A pulsed high-flux source of nearly monoenergetic atomic oxygen was developed to perform accelerated erosion testing of spacecraft materials in a simulated low-earth orbit (LEO) environment. Molecular oxygen is introduced into an evacuated conical expansion nozzle at several atmospheres pressure through a pulsed molecular beam valve. A laser-induced breakdown is generated in the nozzle throat by a pulsed CO2 TEA laser. The resulting plasma is heated by the ensuing laser-supported detonation wave, and then it rapidly expands and cools. An atomic oxygen beam is generated with fluxes above 10 to the 18th atoms per pulse at 8 + or - 1.6 km/s with an ion content below 1 percent for LEO testing. Materials testing yielded the same surface oxygen enrichment in polyethylene samples as observed on the STS mission, and scanning electron micrographs of the irradiated polymer surfaces showed an erosion morphology similar to that obtained on low earth orbit.

Krech, Robert H.↗

Effects of microgravity or simulated launch on testicular function in rats

Reproductive toxicology and cellular and molecular biology approaches were used to evaluate testicular function in rats from Cosmos 2044. It is found that concentrations of testosterone in testicular tissue or peripheral blood plasma were reduced in flight rates to less than 20 percent of values for simulated-launch or vivarium controls. Spermatogenesis was essentially normal in flight rats, but production of testosterone was severely depressed.

Amann, R. P.↗