Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 199 records · Page 11

Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow

To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54 turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJkg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilizing the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly under predicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJkg condition, but with up to a factor of three over prediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJkg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.

Earth entry↗

Surface disordering of Pb(110)

A molecular dynamics simulation of Pb(110), using embedded atom method energy functional and two-body potential derived for the purpose, shows that near 400 K the three surface layers begin to disorder well before the bulk (interior) layers. Transfer of atoms from the first subsurface layer to the surface layer occurs. Disordering propagates beyond the top three layers above 550 K, accompanied by formation of an adlayer and vacancies in the top three layers. Behavior of the two-dimensional layer structure factors indicates that disordering is anisotropic. Simulation results are consistent with experimental observations of surface roughening near 400 K and more extensive surface ordering above 525 K. Results are consistent with simulations for Ni and Al.

Tibbits, P.↗

Tokamak Spectroscopy for X-Ray Astronomy

This paper presents the measured x-ray and Extreme Ultraviolet (XUV) spectra of three astrophysically abundant elements (Fe, Ca and Ne) from three different tokamak plasmas. In every case, each spectrum touches on an issue of atomic physics that is important for simulation codes to be used in the analysis of high spectral resolution data from current and future x-ray telescopes. The utility of the tokamak as a laboratory test bed for astrophysical data is demonstrated. Simple models generated with the HULLAC suite of codes demonstrate how the atomic physics issues studied can affect the interpretation of astrophysical data.

Fournier, Kevin B.↗

The measurement of atomic oxygen in the mesosphere and lower thermosphere

A discussion of the laboratory calibration, calibration simulation, and accuracy of a resonance fluorescence instrument for measuring oxygen atoms is presented. The discussion demonstrates that a sensor can be calibrated with good accuracy to obtain in situ measurements of the oxygen density in the mesosphere/lower thermosphere. Oxygen data are reported from six rocket flights. These data represent four midlatitude flights and two auroral flights. There are two night and two day profiles for the midlatitude experiments. The maximum densities found from these sets of data in the 90-100 km regions are between 1.5 x 10 to the 11th/cu cm and 3.5 x 10 to the 11th/cu cm. Agreement with the MSIS-86 model predictions is obtained for altitudes above 120 km for midlatitude flights. The auroral data are less than the model prediction by about a factor of 2. The data show structure that is probably indicative of gravity wave effects of a nominal 5-km vertical wavelength.

Sharp, William E.↗

Materials selection for long life in LEO: A critical evaluation of atomic oxygen testing with thermal atom systems

The use of thermal atom test methods as a materials selection and screening technique for low-Earth orbit (LEO) spacecraft is critically evaluated. The chemistry and physics of thermal atom environments are compared with the LEO environment. The relative reactivities of a number of materials determined to be in thermal atom environments are compared to those observed in LEO and in high quality LEO simulations. Reaction efficiencies measured in a new type of thermal atom apparatus are one-hundredth to one-thousandth those observed in LEO, and many materials showing nearly identical reactivities in LEO show relative reactivities differing by as much as a factor of 8 in thermal atom systems. A simple phenomenological kinetic model for the reaction of oxygen atoms with organic materials can be used to explain the differences in reactivity in different environments. Certain specific thermal test environments can be used as reliable materials screening tools. Using thermal atom methods to predict material lifetime in LEO requires direct calibration of the method against LEO data or high quality simulation data for each material.

Koontz, S. L.↗

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗

Nitric Oxide Concentrations from Solar Maximum to Solar Minimum Years Derived from SABER

Nitric oxide (NO) is of particular importance in the upper atmosphere nominally the region between 100 km and 250 km, as it plays the role in cooling and regulating the energy budget of the system. The NO cooling rates have been routinely derived for the past two decades from the measurements made by Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite (TIMED) satellite since 2002. However, a comprehensive long-term dataset of global thermospheric NO concentration as a function of local time is currently unavailable. Given that solar activity significantly influences NO abundance on timescales ranging from days to decades, analyzing the behavior of NO concentration during solar storm and solar cycle can provide valuable insights into thermospheric energy balance, improving storm-time and climatological knowledge. In this study, we used newly calculated NO cooling rates and simulations of neutral kinetic temperature and atomic oxygen number density from the MSIS 2.1 model to derive 20-year NO concentration data. The global annual mean NO concentration shows a strong correlation with the solar cycle, with high concentrations and upper boundaries up to 250 km during solar maximum years and lower concentrations and upper boundaries up to ~230 km during solar minimum years. We also compared NO concentrations obtained from SABER measurements and MSIS 2.1 simulations during solar maximum (2014) and solar minimum (2009) years and found good agreement between the two during both day and night.

Ningchao Wang↗

Atomic oxygen effects on thin film space coatings studied by spectroscopic ellipsometry, atomic force microscopy, and laser light scattering

The University of Nebraska is currently evaluating Low Earth Orbit (LEO) simulation techniques as well as a variety of thin film protective coatings to withstand atomic oxygen (AO) degradation. Both oxygen plasma ashers and an electron cyclotron resonance (ECR) source are being used for LEO simulation. Thin film coatings are characterized by optical techniques including Variable Angle Spectroscopic Ellipsometry, Optical spectrophotometry, and laser light scatterometry. Atomic Force Microscopy (AFM) is also used to characterize surface morphology. Results on diamondlike carbon (DLC) films show that DLC degrades with simulated AO exposure at a rate comparable to Kapton polyimide. Since DLC is not as susceptible to environmental factors such as moisture absorption, it could potentially provide more accurate measurements of AO fluence on short space flights.

Synowicki, R. A.↗

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗

Image Processing

A new spinoff product was derived from Geospectra Corporation's expertise in processing LANDSAT data in a software package. Called ATOM (for Automatic Topographic Mapping), it's capable of digitally extracting elevation information from stereo photos taken by spaceborne cameras. ATOM offers a new dimension of realism in applications involving terrain simulations, producing extremely precise maps of an area's elevations at a lower cost than traditional methods. ATOM has a number of applications involving defense training simulations and offers utility in architecture, urban planning, forestry, petroleum and mineral exploration.

Source record↗

Monte Carlo Simulation of Alloy Design Techniques: Fracture and Welding Studied Using the BFS Method for Alloys

Large-scale simulations of dynamic processes at the atomic level have developed into one of the main areas of work in computational materials science. Until recently, severe computational restrictions, as well as the lack of accurate methods for calculating the energetics, resulted in slower growth in the area than that required by current alloy design programs. The Computational Materials Group at the NASA Lewis Research Center is devoted to the development of powerful, accurate, economical tools to aid in alloy design. These include the BFS (Bozzolo, Ferrante, and Smith) method for alloys (ref. 1) and the development of dedicated software for large-scale simulations based on Monte Carlo- Metropolis numerical techniques, as well as state-of-the-art visualization methods. Our previous effort linking theoretical and computational modeling resulted in the successful prediction of the microstructure of a five-element intermetallic alloy, in excellent agreement with experimental results (refs. 2 and 3). This effort also produced a complete description of the role of alloying additions in intermetallic binary, ternary, and higher order alloys (ref. 4).

Bozzolo, Guillermo H.↗

Soft X-Ray Imaging of Earth's Dayside Magnetosheath and Cusps Using Hybrid Simulations

Interactions between solar wind ions and neutral hydrogen atoms in Earth's exosphere can lead to the emission of soft X-rays. Upcoming missions such as SMILE and LEXI aim to use soft X-ray imaging to study the global structure of the magnetosphere. Although the magnetosheath and dayside magnetopause can often be driven by kinetic physics, it has typically been omitted from fluid simulations used to predict X-ray emissions. We study the possible results of soft X-ray imaging using hybrid simulations under quasi-radial interplanetary magnetic fields, where ion-ion instabilities drive ultra-low frequency foreshock waves, leading to turbulence in the magnetosheath, affecting the dynamics of the cusp and magnetopause. We simulate soft X-ray emission to determine what may be seen by missions such as LEXI, and evaluate the possibility of identifying kinetic structures. While kinetic structures are visible in high-cadence imaging, current instruments may not have the time resolution to discern kinetic signals.

J Ng↗

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗

Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification

A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.

DSMC↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Measurements in liquid fuel sprays

A ground test facility is being established at NASA Lewis Research Center to simulate the environmental and flight conditions needed to study adverse weather effects. One of the most important components is the water spray system which consists of many nozzles fitted on spray bars. Water is injected through air-assisted atomizers to generate uniform size drops to simulate icing in clouds. The primary objective is to provide experimental data on drop size distribution over a wide range of operating conditions. Correlation equations for mean drop size and initial injection parameters are being determined to assist in the design and modification of the Altitude Wind Tunnel. Special emphasis is being placed on the study of the aerodynamic structure of the air-assisted atomizer sprays. Detailed measurements of the variation of drop size distribution and velocity as a function of time and space are being made. Accurate initial and boundary conditions are being provided for computer model evaluation.

Chigier, N.↗

Measurements of Polyatomic Molecule Formation on an Icy Grain Analog Using Fast Atoms

Carbon dioxide has been produced from the impact of a monoenergetic O(P-3) beam upon a surface cooled to 4.8 K and covered with a CO ice. Using temperature-programmed desorption and mass spectrometer detection, we have detected increasing amounts of CO2 formation with O(P-3) energies of 2, 5, 10, and 14 eV. This is the first measurement of polyatomic molecule formation on a surface with superthermal atoms. The goal of this work is to detect other polyatomic species, such as CH3OH, which can be formed under conditions that simulate the grain temperature, surface coverage, and superthermal atoms present in shock-heated circumstellar and interstellar regions.

Chutjian, A.↗