Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation 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.

At least 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Phase-field simulations of GaN growth by selective area epitaxy from complex mask geometries

Three-dimensional phase-field simulations of GaN growth by selective area epitaxy were performed. The model includes a crystallographic-orientation-dependent deposition rate and arbitrarily complex mask geometries. The orientation-dependent deposition rate can be determined from experimental measurements of the relative growth rates of low-index crystallographic facets. Growth on various complex mask geometries was simulated on both c-plane and a-plane template layers. Agreement was observed between simulations and experiment, including complex phenomena occurring at the intersections between facets. The sources of the discrepancies between simulated and experimental morphologies were also investigated. The model provides a route to optimize masks and processing conditions during materials synthesis for solar cells, light-emitting diodes, and other electronic and opto-electronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2001

Microphysics, Radiation and Surface Processes in the Goddard Cumulus Ensemble (GCE) Model

The response of cloud systems to their environment is an important link in a chain of processes responsible for monsoons, frontal depression, El Nino Southern Oscillation (ENSO) episodes and other climate variations (e.g., 30-60 day intra-seasonal oscillations). Numerical models of cloud properties provide essential insights into the interactions of clouds with each other, with their surroundings, and with land and ocean surfaces. Significant advances are currently being made in the modeling of rainfall and rain-related cloud processes, ranging in scales from the very small up to the simulation of an extensive population of raining cumulus clouds in a tropical- or midlatitude-storm environment. The Goddard Cumulus Ensemble (GCE) model is a multi-dimensional nonhydrostatic dynamic/microphysical cloud resolving model. It has been used to simulate many different mesoscale convective systems that occurred in various geographic locations. In this paper, recent GCE model improvements (microphysics, radiation and surface processes) will be described as well as their impact on the development of precipitation events from various geographic locations. The performance of these new physical processes will be examined by comparing the model results with observations. In addition, the explicit interactive processes between cloud, radiation and surface processes will be discussed.

Tao, Wei-Kuo↗
DOE OSTI2021

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
NASA NTRS1989

Approximate simulation model for analysis and optimization in engineering system design

Computational support of the engineering design process routinely requires mathematical models of behavior to inform designers of the system response to external stimuli. However, designers also need to know the effect of the changes in design variable values on the system behavior. For large engineering systems, the conventional way of evaluating these effects by repetitive simulation of behavior for perturbed variables is impractical because of excessive cost and inadequate accuracy. An alternative is described based on recently developed system sensitivity analysis that is combined with extrapolation to form a model of design. This design model is complementary to the model of behavior and capable of direct simulation of the effects of design variable changes.

Sobieszczanski-Sobieski, Jaroslaw↗
DOE OSTI2016

Laser pulse driven terahertz generation via resonant transition radiation in inhomogeneous plasmas

An intense, short laser pulse propagating across a plasma boundary ponderomotively drives THz radiation. Full format PIC simulations and theoretical analysis are conducted to investigate the properties of this radiation. Simulation results which show the THz emission originates in regions of varying density and covers a broad spectrum with maximum frequency close to the maximum plasma frequency. In the case of a sharp vacuum-plasma boundary, the radiation is generated symmetrically at the plasma entrance and exit, and its properties are independent of plasma density when the density exceeds a characteristic value determined by the product of the plasma frequency and the laser pulse duration. For a diffuse vacuum-plasma boundary, the emission from the plasma entrance and exit is asymmetric: increasing and decreasing density ramps enhance and diminish the radiated energy, respectively. Enhancements by a factor of 50 are found and simulations show that a 1.66 J, 50 fs driver pulse can generate ∼400 μJ of THz radiation in a 1.2 mm increasing density ramp. We present a model that attributes this effect to a plasma resonance process in the density ramp. The results from the model match those of the simulations for ramp lengths less than 600 μm. For longer ramps for which simulations are too time consuming, the model shows that the amount of radiation reaches a maximum at a ramp length determined by collisional absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Simulation studies for surfaces and materials strength

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

Halicioglu, T.↗
NASA NTRS1972

Computer system simulation and analysis

A simulation model has been developed as a tool for studying different computer system configurations. The model simulates a computer system in terms of the number and speed of hardware subsystems and the runs processed by the system in terms of equivalent adds for computation and data bits transferred for I/O. The model was written in the GPSS simulation language.

Williams, T. G.↗
NASA NTRS2023

Electrification Within Wintertime Stratiform Regions Sampled During the 2020/2022 NASA IMPACTS Field Campaign

Two nor'easter events—sampled during the NASA Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign—were examined to characterize the microphysics in relation to the underlying electrification processes within wintertime stratiform regions. A theoretical model was developed to determine whether accretion or diffusion growth regimes were preferential during periods of greatest electrification. Model simulation with electrification parameterization was used to provide supplemental context to the physical processes of in-cloud microphysics and electrification. The strongest electric fields (i.e., ∼80 V m −1 at 20 km) during the 2020 NASA IMPACTS deployment was associated with large non-rimed ice crystals colliding with each other. During the 29–30 January 2022 science flight, the NASA P-3 microphysical probe data demonstrated that non-inductive charging was possible off the coastline of Cape Cod, Massachusetts. Later in the science flight, when the NASA P-3 and ER-2 were coordinating with each other, measured electric fields consistently were less than 8 V m −1 and electrification was subdued owing to reduced concentrations of graupel and large ice hydrometeors. Altogether, the in-situ observations provide evidence for the non-riming collisional charging mechanism and demonstrates that graupel and supercooled liquid water may not be necessary for weak electrification within wintertime stratiform regions. Model output from simulation of both events suggested that the main synoptic snowbands were associated with elevated hydrometeor snow charge density and electric fields.

snowfall↗
DOE OSTI2025

Advancing the Understanding of Cloud Microphysical Processes and Aerosol Indirect Effects in High-Latitude Mixed-Phase Clouds by Linking ARM Measurements with Climate Model Simulations (Final Report)

The key objectives of this project were to advance our understanding of cloud microphysical characteristics and aerosol indirect effects on mixed-phase clouds in high latitudes. To improve the representation of ice and mixed-phase clouds in Earth System Models (ESMs), we propose an integrated observation and modeling study of cloud macro- and microphysical properties, including spatial heterogeneities, mass partitioning between ice crystals and supercooled liquid water, effects of ice nucleating particles (INPs), and efficiency of secondary ice production (SIP), etc. Specifically, we took four main approaches in this project: (1) examining macro- and microphysical properties of ice and mixed-phase clouds based on in-situ and ground-based observations from multiple field campaigns funded by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program, including the Mixed-Phase Arctic Cloud Experiment (M-PACE), Indirect and Semi-Direct Aerosol Campaign (ISDAC), Ice Nucleating Particle Sources at Oliktok Point (INPOP), ARM West Antarctic Radiation Experiment (AWARE), Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS), and Macquarie Island Cloud and Radiation Experiment (MICRE); (2) evaluating the DOE Energy Exascale Earth System Model (E3SM) simulations based on observations, particularly for ice and mixed-phase cloud microphysical properties; (3) examining the impacts of INPs on ice and mixed-phase clouds. Specifically, a series of comparisons were conducted using observations over the Arctic, Southern Ocean, and Antarctica, including comparisons between the lower and higher southern latitudes as well as comparisons between the northern and southern hemispheres. In addition, aerosol indirect effects from distinct sources of dust particles were examined; and (4) investigating the impacts of SIP. Ultimately, these results helped to improve cloud microphysics and aerosol-cloud interaction parameterizations in the E3SM model. Overall, the project provided improved understanding regarding various factors, including thermodynamic, dynamic, and aerosol conditions, on the micro- and macrophysical properties of ice and mixed-phase clouds in the high latitudes. Resulting analysis helped to provide an improved physical basis for refining the current cloud microphysics parameterizations related to ice and mixed-phase clouds in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Theoretical detection threshold of the proton-acoustic range verification technique

Purpose: Range verification in proton therapy using the proton-acoustic signal induced in the Bragg peak was investigated for typical clinical scenarios. The signal generation and detection processes were simulated in order to determine the signal-to-noise limits. Methods: An analytical model was used to calculate the dose distribution and local pressure rise (per proton) for beams of different energy (100 and 160 MeV) and spot widths (1, 5, and 10 mm) in a water phantom. In this method, the acoustic waves propagating from the Bragg peak were generated by the general 3D pressure wave equation implemented using a finite element method. Various beam pulse widths (0.1–10 μs) were simulated by convolving the acoustic waves with Gaussian kernels. A realistic PZT ultrasound transducer (5 cm diameter) was simulated with a Butterworth bandpass filter with consideration of random noise based on a model of thermal noise in the transducer. The signal-to-noise ratio on a per-proton basis was calculated, determining the minimum number of protons required to generate a detectable pulse. The maximum spatial resolution of the proton-acoustic imaging modality was also estimated from the signal spectrum. Results: The calculated noise in the transducer was 12–28 mPa, depending on the transducer central frequency (70–380 kHz). The minimum number of protons detectable by the technique was on the order of 3–30 × 10{sup 6} per pulse, with 30–800 mGy dose per pulse at the Bragg peak. Wider pulses produced signal with lower acoustic frequencies, with 10 μs pulses producing signals with frequency less than 100 kHz. Conclusions: The proton-acoustic process was simulated using a realistic model and the minimal detection limit was established for proton-acoustic range validation. These limits correspond to a best case scenario with a single large detector with no losses and detector thermal noise as the sensitivity limiting factor. Our study indicated practical proton-acoustic range verification may be feasible with approximately 5 × 10{sup 6} protons/pulse and beam current.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2004

Integrated Modeling Approach for the Terrestrial Planet Finder Mission

Because of the complexity of the Terrestrial Planet Finder (TPF) design concepts, the project will rely heavily on the use of engineering and science simulations to predict on-orbit performance. Furthermore, current understanding of these missions indicates that the 3m to 8m class optical systems need to be as stable as picometers in wavefront and sub-milli arcsec in pointing. These extremely small requirements impose on the models a level of predictive accuracy heretofore never achieved, especially in the area of microgravity effects, material property accuracy, thermal solution convergence, and all other second order modeling effects typically ignored. New modeling tools and analysis paradigms are developed which emphasize computational accuracy and fully integrated analytical simulations. The process is demonstrated on sample problems using the TPF Coronagraph design concept. The TPF project is also planning a suite of testbeds through which various aspects of the models and simulations will be verified.

terrestrial planet finder↗
DOE OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Large Scale GPU Accelerated PPMLR-MHD Simulations for Space Weather Forecast

PPMLR-MHD is a new magnetohydrodynamics (MHD) model used to simulate the interactions of the solar wind with the magnetosphere, which has been proved to be the key element of the space weather cause-and-effect chain process from the Sun to Earth. Compared to existing MHD methods, PPMLR-MHD achieves the advantage of high order spatial accuracy and low numerical dissipation. However, the accuracy comes at a cost. On one hand, this method requires more intensive computation. On the other hand, more boundary data is subject to be transferred during the process of simulation. In this work, we present a parallel hybrid solution of the PPMLR-MHD model implemented using the computing capabilities of both CPUs and GPUs. We demonstrate that our optimized implementation alleviates the data transfer overhead by using GPU Direct technology and can scale up to 151 processes and achieve significant performance gains by distributing the workload among the CPUs and GPUs on Titan at Oak Ridge National Laboratory. The performance results show that our implementation is fast enough to carry out highly accurate MHD simulations in real time.

Computer Science↗
DOE OSTI2015

Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter

Abstract Hydraulic stimulation of low‐permeability rocks in enhanced geothermal systems, shale resources, and CO 2 storage aquifers can trigger microseismic events, also known as microearthquakes (MEQs). The distribution of microseismic source locations in the reservoir may reveal important information about the distribution of hydraulic and geomechanical rock properties. In this paper, we present a framework for conditioning heterogeneous rock permeability and geomechanical property distributions on microseismic data. To simulate the multiphysics processes in these systems, we combine a fully coupled flow and geomechanical model with the Mohr‐Coulomb type rock failure criterion. The resulting multiphysics simulation constitutes the forecast model that relates microseismic source locations to reservoir rock properties. We adopt this forward model in an ensemble Kalman filter (EnKF) data assimilation framework to jointly estimate reservoir permeability and geomechanical property distributions from injection‐induced microseismic response measurements. We show that integration of a large number of spatially correlated microseismic data with practical ensemble sizes can lead to severe underestimation of ensemble spread, and eventually ensemble collapse. To mitigate the variance underestimation issue, two low‐rank data representation schemes are presented and discussed. In the first approach, microseismic data are projected onto a low‐dimensional subspace defined by the left singular vectors of the perturbed observation matrix. The second method uses a coarser grid for representing the microseismic data. A series of numerical experiments is presented to evaluate the performance of the proposed methods and to illustrate their applicability for assimilating microseismic data into coupled flow and geomechanical forward models to estimate multiphysics rock properties.

Tarrahi, Mohammadali↗
DOE OSTI2019

Simulation of Gas Production from Multilayered Hydrate-Bearing Media with Fully Coupled Flow, Thermal, Chemical and Geomechanical Processes Using TOUGH + Millstone. Part 1: Numerical Modeling of Hydrates

TOUGH + Millstone has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH4-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH + HYDRATE simulator, V2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the behavior of hydrate-bearing systems undergoing changes and includes the most recent advances in the description of the system properties, coupled seamlessly with (b) Millstone V1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of TOUGH + Millstone are demonstrated in the simulation and analysis of the system flow, thermal and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the first paper of this series, we discuss the physics underlying the T + H hydrate simulator, the constitutive relationships describing the physical, chemical (equilibrium and kinetic) and thermal processes, the states of the CH 4 + H 2 O system and the sources of critically important data, as well as the mathematical approaches used for the development of the of mass and energy balance equations and their solution. Finally, we provide verification examples of the hydrate code against numerical results from the simulation of laboratory and field experiments.

58 GEOSCIENCES↗
DOE OSTI2024

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Watershed response of the Feather River Basin, California, United States of America to future climate changes

Increasing mean annual temperatures under climate change are expected to reduce seasonal snowpack, increase evapotranspiration (ET), and alter summer baseflow in headwater watersheds worldwide. Strong regional variability in hydrologic responses highlights the need for catchment-scale, physically based models to assess future flood risk and water availability. This study examines climate-driven changes in the hydrologic response of the Upper Feather River watershed in the Sierra Nevada Mountains, California. Four priority climate models and two representative concentration pathways (RCP4.5 and RCP8.5) are used to evaluate future hydrologic responses of the watershed. Hydrologic processes are simulated using an objectively calibrated Soil and Water Assessment Tool Plus model for a historical baseline (1986–2005) and a future period 2070–2099), driven by observed and projected precipitation and temperature. Results indicate a declining contribution of snowfall to annual precipitation, with peak snowfall and water yield shifting 1–3 months earlier. Long-term annual maximum flows are projected to increase considerably, whereas low-flow responses are mixed, with both increases and decreases projected by the end of the century. These findings highlight the need for adaptive watershed management to enhance flood protection, water storage, and drought resilience. Future water resource planning should also account for one-to-three-month shifts in peak water yield and surface runoff due to changes in snowmelt timing and a lower snowfall-to-rainfall ratio under climate change. The modeling framework and insights are transferable to other snow-dominated headwater watersheds experiencing climate-driven change.

Tigabu, T [UC Davis]↗
NASA NTRS1985

Performance evaluation of a simulated data-flow computer with low-resolution actors

Basic problems related to the exploitation of parallelism in a program include sequencing of the instructions and communication of the data. It is pointed out that the data-flow approach offers an elegant solution to the sequencing problem, since all data dependencies are automatically handled and only instructions with ready input sets are activated. It is shown that a change in the level of subcomputations (actors) affects communications costs. The concept of variable resolution is discussed, and the testbed environment is examined. Attention is given to the architecture of the processing elements, the communication network, and the simulators. A description of the analytical model is also provided. Simulation and results are discussed, taking into account test programs and allocation, the variation of the number of processing elements, the variation of the resolution in directed acyclic graphs, performance in processing loops, and array handling.

Gaudiot, J. L.↗