Search NASA⌕ Search

SEARCH · Search NASA

Results for “process 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 577 records · Page 32

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Modeling Strategies and Simulations of TOGA-COARE Convective Systems - Part 1

Interactions between deep tropical clouds over the western Pacific warm pool and the larger-scale environment are key to understanding climate change. Cloud models are an extremely useful tool in simulating and providing statistical information on heat and moisture transfer processes between cloud systems and the environment, and can therefore be utilized to substantially improve cloud parameterizations in climate models. In this paper, the Goddard Cumulus Ensemble (GCE) cloud-resolving model is used in multi-day simulations of deep tropical convective activity over the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). Large-scale temperature and moisture advective tendencies, and horizontal momentum from the TOGA-COARE Intensive Flux Array (IFA) region, are applied to the GCE version which incorporates cyclical boundary conditions. Sensitivity experiments show that grid domain size produces the largest response to domain-mean temperature and moisture deviations, as well as cloudiness, when compared to grid horizontal or vertical resolution, and advection scheme. It is found that a minimum grid-domain size of 500 km is needed to adequately resolve the convective cloud features. The control experiment shows that the atmospheric heating and moistening is primarily a response to cloud latent processes of condensation/evaporation, and deposition/sublimation, and to a lesser extent, melting of ice particles. Air-sea exchange of heat and moisture is found to be significant, but of secondary importance, while the radiational response is small. The simulated rainfall and atmospheric heating and moistening, agrees well with observations, and performs favorably to other models simulating this case.

Johnson, Daniel E.↗

On Fast Post-Processing of Global Positioning System Simulator Truth Data and Receiver Measurements and Solutions Data

Post-Processing of data related to a Global Positioning System (GPS) simulation is an important activity in qualification of a GPS receiver for space flight. Because a GPS simulator is a critical resource it is desirable to move off the pertinent simulation data from the simulator as soon as a test is completed. The simulator data files are usually moved to a Personal Computer (PC), where the post-processing of the receiver logged measurements and solutions data and simulated data is performed. Typically post-processing is accomplished using PC-based commercial software languages and tools. Because of commercial software systems generality their general-purpose functions are notoriously slow and more than often are the bottleneck problem even for short duration experiments. For example, it may take 8 hours to post-process data from a 6-hour simulation. There is a need to do post-processing faster, especially in order to use the previous test results as feedback for a next simulation setup. This paper demonstrates that a fast software linear interpolation algorithm is applicable to a large class of engineering problems, like GPS simulation data post-processing, where computational time is a critical resource and is one of the most important considerations. An approach is developed that allows to speed-up post-processing by an order of magnitude. It is based on improving the post-processing bottleneck interpolation algorithm using apriori information that is specific to the GPS simulation application. The presented post-processing scheme was used in support of a few successful space flight missions carrying GPS receivers. A future approach to solving the post-processing performance problem using Field Programmable Gate Array (FPGA) technology is described.

Kizhner, Semion↗

Relationship between Surface Self-Diffusion and Bubble Mobility in Solids: Theory and Atomistic Simulation

Diffusion processes in solids are characterized by complex mechanisms occurring at the atomistic level. The relevant theoretical concepts are sufficiently developed in the case of the diffusion of point defects. At the same time, the diffusion of objects such as clusters of vacancies or nanometer-size cavities in the crystal lattice is described only in terms of continuous-medium approximation. This paper attempts to link the existing theory of diffusion of bubbles in a crystalline matrix with the atomistic description of the process on the basis of the molecular dynamics method. For simplicity, the diffusion of cavities in a two-dimensional Lennard-Jones lattice is considered. A controlled molecular dynamics method is proposed to speed up the direct computation of the diffusion rate of bubbles. Based on the results of simulation, a practical interpretation is given to a parameter of the continuum theory such as the rate of surface self-diffusion. The applicability range of the continuum theory is demonstrated, and principles for the refinement of the theory in the case of minimum-size bubbles are formulated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parallel processing for nonlinear dynamics simulations of structures including rotating bladed-disk assemblies

The principal objective of this research is to develop, test, and implement coarse-grained, parallel-processing strategies for nonlinear dynamic simulations of practical structural problems. There are contributions to four main areas: finite element modeling and analysis of rotational dynamics, numerical algorithms for parallel nonlinear solutions, automatic partitioning techniques to effect load-balancing among processors, and an integrated parallel analysis system.

Hsieh, Shang-Hsien↗

A Graphics Processing Unit–Based, Industrial Grade Compositional Reservoir Simulator

Summary Recently, graphics processing units (GPUs) have been demonstrated to provide a significant performance benefit for black-oil reservoir simulation, as well as flash calculations that serve an important role in compositional simulation. A comprehensive approach to compositional simulation based on GPUs has yet to emerge, and the question remains as to whether the benefits observed in black-oil simulation persist with a more complex fluid description. We present a positive answer to this question through the extension of a commercial GPU-based black-oil simulator to include a compositional description based on standard cubic equations of state (EOSs). We describe the motivations for the selected nonlinear formulation, including the choice of primary variables and iteration scheme, and support for both fully implicit methods (FIMs) and adaptive implicit methods (AIMs). We then present performance results on an example sector model and simplified synthetic case designed to allow a detailed examination of runtime and memory scaling with respect to the number of hydrocarbon components and model size, as well as the number of processors. We finally show results from two complex asset models (synthetic and real) and examine performance scaling with respect to GPU generation, demonstrating that performance correlates strongly with GPU memory bandwidth. NOTE: This paper is also published as part of the 2021 SPE Reservoir Simulation Conference Special Issue.

Engineering↗

Determination of Chemical Kinetic Rate Constants of a Model for Carbothermal Processing of Lunar Regolith Simulant Using Methane

We have previously developed a chemical conversion model of the carbothermal processing of lunar regolith using methane to predict the rate of production of carbon monoxide. In this carbothermal process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Hydrogen is carried away by the exiting stream of gases and carbon is deposited on the melt surface. The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. In our model, we assume that the flux of carbon deposited is equal to the product of the surface reaction rate constant gamma and the concentration of methane adjacent to the melt surface. Similarly, the rate of consumption of carbon per unit volume in the melt is equal to the product of the melt reaction rate constant k and the concentrations of carbon and metal oxide in the melt. In this paper, we describe our effort to determine gamma and k by comparison of the predictions from our model with test data obtained by ORBITEC (Orbital Technologies Corporation). The concentration of methane adjacent to the melt surface is a necessary input to the model. It is inferred from the test data by a mass balance of methane, adopting the usual assumptions of the continuously-stirred-tank-reactor model, whereby the average concentration of a given gaseous species equals its exit concentration. The reaction rates gamma and k have been determined by a non-linear least-squares fit to the test data for the production of carbon monoxide and the fraction of the incoming methane that is converted. The comparison of test data with our model predictions using the determined chemical kinetic rate constants provides a consistent interpretation of the process over the full range of temperatures, pressures, and methane flow rates used in the tests, thereby increasing our confidence to use the model for scale-up purposes.

Balasubramaniam, R↗

Estimating outflow masses and velocities in merger simulations: Impact of r -process heating and neutrino cooling

Here, the determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fall-back timescale of any bound ejecta. We derive a model for the asymptotic velocity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full three-dimensional fluid simulations of the outflows, as nonlocal effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long three-dimensional simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian stability and force modeling for uncertain machining processes

Accurately simulating machining operations requires knowledge of the cutting force model and system frequency response. However, this data is collected using specialized instruments in an ex-situ manner. Bayesian statistical methods instead learn the system parameters using cutting test data, but to date, these approaches have only considered milling stability. This paper presents a physics-based Bayesian framework which incorporates both spindle power and milling stability. Initial probabilistic descriptions of the system parameters are propagated through a set of physics functions to form probabilistic predictions about the milling process. The system parameters are then updated using automatically selected cutting tests to reduce parameter uncertainty and identify more productive cutting conditions, where spindle power measurements are used to learn the cutting force model. The framework is demonstrated through both numerical and experimental case studies. Results show that the approach accurately identifies both the system natural frequency and cutting force model.

42 ENGINEERING↗

The R -process Alliance: Enrichment of r -process Elements in a Simulated Milky Way–like Galaxy

We study the formation of stars with varying amounts of heavy elements synthesized by the rapid neutron-capture process (r-process) based on our detailed cosmological zoom-in simulation of a Milky Way–like galaxy with an N-body/smoothed particle hydrodynamics code, asura. Most stars with no overabundance in r-process elements, as well as the strongly r-process-enhanced (RPE) r-II stars ([Eu/Fe] > +0.7), are formed in dwarf galaxies accreted by the Milky Way within the 6 Gyr after the Big Bang. In contrast, over half of the moderately enhanced r-I stars (+0.3 < [Eu/Fe] ≤ +0.7) are formed in the main in situ disk after 6 Gyr. Our results suggest that the fraction of r-I and r-II stars formed in disrupted dwarf galaxies is larger the higher their [Eu/Fe] is. Accordingly, the most strongly enhanced r-III stars ([Eu/Fe] > +2.0) are formed in accreted components. These results suggest that non-r-process-enhanced stars and r-II stars are mainly formed in low-mass dwarf galaxies that hosted either none or a single neutron star merger, while the r-I stars tend to form in the well-mixed in situ disk. We compare our findings with high-resolution spectroscopic observations of RPE metal-poor stars in the halo and dwarf galaxies, including those collected by the R-Process Alliance. We conclude that observed [Eu/Fe] and [Eu/Mg] ratios can be employed in chemical tagging of the Milky Way’s accretion history.

Chemical enrichment↗

Fast Emulation of Expensive Simulations using Approximate Gaussian Processes [Slides]

Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources

97 MATHEMATICS AND COMPUTING↗

The stratospheric sulfate aerosol layer - Processes, models, observations, and simulations

After briefly reviewing the observational data on the stratospheric sulfate aerosol layer, the chemical and physical processes that are likely to fix the properties of the layer are discussed. We present appropriate continuity equations for aerosol particles, and show how to solve the equations on a digital computer. Simulations of the unperturbed aerosol layer by various published models are discussed and the sensitivity of layer characteristics to variations in several aerosol model parameters is studied. We discuss model applications to anthropogenic pollution problems and demonstrate that moderate levels of aerospace activity (supersonic transport and Space Shuttle operations) will probably have only a negligible effect on global climate. Finally, we evaluate the possible climatic effect of a ten-fold increase in the atmospheric abundance of carbonyl sulfide.

Whitten, R. C.↗

Permanganate Oxidation of Defense Waste Processing Facility (DWPF) Recycle Collection Tank (RCT) Simulants – Protocol Runs for Nominal and Chemical Process Cell (CPC) Foamover Conditions

Savannah River National Laboratory (SRNL) has completed a work scope in support of the Savannah River Remediation (SRR) Nitric-Glycolic flowsheet. Glycolic acid will replace formic acid as a reducing agent in the pre-processing of High-Level Waste sludge performed at the Defense Waste Processing Facility (DWPF). Glycolic acid will perform the same function as formic acid, namely reduction of mercury and adjustment of feed rheology and melter oxidation/reduction potential, while significantly reducing the potential for hydrogen generation in DWPF processing. Development testing has demonstrated glycolic acid virtually eliminates hydrogen production in the pre-processing steps. The Nitric-Glycolic flowsheet significantly improves DWPF’s ability to address one of the Savannah River Site’s key challenges - the incorporation of effluent received from the Salt Waste Processing Facility (SWPF). SWPF will deliver significant effluent volume to DWPF, resulting in a concurrent increase in DWPF effluent returned to the Concentration, Storage and Transfer Facilities (CSTF). This work scope demonstrates that glycolate can be destroyed under the Nitric-Glycolic flowsheet prior to returning the DWPF recycle stream to CSTF. To avoid potential flammability issues due to thermolysis of glycolate in CSTF, SRR tasked SRNL to quantify and mitigate glycolate returns via DWPF’s recycle stream. The development of a strategy for glycolic mitigation was initiated with a system’s engineering workshop. Various chemical and/or physical solutions for how and where to destroy glycolate were considered – consistent with DWPF’s operational capabilities and process requirements. The workshop identified chemical oxidation of glycolate within the DWPF Recycle Collection Tank (RCT) as the most promising option. (The RCT collects offgas condensate from pre-processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF.) A series of down-select experimental tests were performed that showed sodium permanganate was the most effective oxidant. This work was performed to identify an appropriate operational protocol for DWPF implementation. Building upon the downselect experiments, the work targeted an improved description of suitable RCT conditions prior to the strike (pH, etc.), permanganate addition rate, and destruction reaction kinetics. Data collected were also used to determine the best expression of reaction stoichiometry and the effective minimum permanganate addition required.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Low-speed impact phenomena and orbital resonances in the moon- and planet-building process

A simulation of collisional and gravitational interaction in the early solar system generates planets approximately 1000 km in diameter from an initial swarm of kilometer sized planetesimals. The model treats collisions according to experimental and theoretical impact results (such as rebound, cratering, and catastrophic fragmentation) for a variety of materials whose parameters span plausible values for early solid objects. The small planets form in approximately 1000 yr, during which time most of the mass of the system continues to reside in particles near the original size. The simulation is terminated when the largest objects' random motion is of smaller dimension than their collision cross-sections. The few 1000 km planets may act as seeds for the subsequent, gradual, accretional growth into full-sized planets.

Chapman, C. R.↗

An evaluation of excitation techniques for time domain based flutter data processing

Numerically simulated, wind tunnel and flight test experimental data were used to evaluate various types of system excitation sources for subcritical flutter testing. Emphasis was placed upon the determination of modal frequencies and damping from relatively short time series records. Parameter identification based upon digital time series models was used to predict modal characteristics from subcritical test data. Various types of excitation were considered. These were: 1) impulsive loading and the resulting transient free vibration response, 2) random forced response, in which the 'mechanical' forcing function could be directly measured, 3) sine sweep forced response, in which the 'mechanical' forcing function could be directly measured, and 4) forced response from unsteady aerodynamic excitation or turbulence, in which the excitation could only be indirectly inferred from other measurements.

Pinkelman, J. K.↗

Parallel point defect identification in molecular dynamics simulations without post-processing: A compute and dump style for LAMMPS

Two classes for the molecular dynamics program LAMMPS – one compute style and one dump style – are presented that are designed to identify, count, and output point defects in cascade damage and related molecular dynamics simulations. The calculations are done in parallel across multiple MPI processes, and can be done at the same time as the original simulation. This drastically reduces storage requirements by eliminating the need to post-process every atom in the system. Non-cubic lattices, free surfaces, and large voids can be eliminated from the output by suitable choices of the reference lattice. The classes are derived from LAMMPS’s Dump and Compute classes, and pose no additional overhead to LAMMPS if they are not used.

42 ENGINEERING↗

A mesh re-zoning technique for finite element simulations of metal forming processes

Based on some fundamental properties of finite element approximations, a mesh re-zoning scheme is proposed for finite element simulations of metal forming problems. It is demonstrated that this technique is indispensable in analyzing many difficult forming processes, especially when there exist corners or very irregular shapes on the boundaries. The algorithm is tested by a backward extrusion process and direct extrusion through a square die.

Cheng, J.-C.↗

Performance-based simulants for Hanford radioactive waste treatment process testing

Radioactive wastes from nuclear fuel processing are stored in large underground storage tanks at the Hanford Nuclear Reservation. Treatment and remediation require that waste feed from the storage tanks be delivered to the Hanford Tank Waste Treatment and Immobilization Plant continually for the duration of the treatment mission. Furthermore, the complex physical and chemical processes required for this mission, the significant scale, and the hazardous nature of the waste necessitate the use of simulants in process testing. A new class of simulants, “performance-based simulants,” has been developed to match the process performance of the simulant to actual waste performance data.

Hanford tank waste treatment↗

Computer simulation of gear tooth manufacturing processes

The use of computer graphics to simulate gear tooth manufacturing procedures is discussed. An analytical basis for the simulation is established for spur gears. The simulation itself, however, is developed not only for spur gears, but for straight bevel gears as well. The applications of the developed procedure extend from the development of finite element models of heretofore intractable geometrical forms, to exploring the fabrication of nonstandard tooth forms.

Mavriplis, Dimitri↗