Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation framework”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

The Predictive Simulation Capability development team of the STMD Game Changing Development sponsored PSI project is implementing computational simulation capability for the efficient and accurate prediction of Plume-Surface Interaction induced surface erosion and cratering in Martian and Lunar environments. The status of the Focus Area 3 of the PSI project in the generation and efficient application of accurate soil particle composition modeling in the Gas-Granular Flow Solver (GGFS) computational framework is presented. The process of constitutive closure model database generation using DEM particle interaction modeling for capturing the effects of irregular particle shape and poly-disperse mixture distribution effects is outlined. This capability has now been ported to NASA supercomputer assets and NASA engineers successfully demonstrated technology and skillset transfer in model generation for spherical and irregularly shaped, mono-disperse and bi-disperse mixture compositions. Assessment of the computational efficiency and practicality of the academic serially executed DEM tools on NASA supercomputers identified the need to migrate to a DEM framework capable of performing parallel simulations in a simultaneous process orchestrated in an automated setup, execution, database extraction, and dataset delivery ready for application simulations. The LIGGGHTS DEM toolset has been selected as the most suitable tool to migrate the DEM simulations. Once the soil model generation process is implemented, models capturing the shape and poly-dispersity effects will be generated to perform much refined validation simulations against the experiments performed under the PSI project. The application readiness of the soil models currently operational in GGFS was presented for the example of a full scale, 3-D simulation of the plume induced erosion and crater formation of the Apollo LM at an elevation of 5m above ground in a low pressure, near vacuum background.

Peter A Liever↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Sensitivity and Effective Parameterization of a Multi-Scale Model of Proton-Exchange-Membrane Water Electrolysis

The development of proton-exchange-membrane water electrolysis (PEMWE) is crucial to establishing a green-hydrogen infrastructure. Active research in hydrogen technologies is focused on minimizing costs, improving efficiency, and mitigating safety risks. In this study, sensitivity analyses were performed using a physics-based model and used to identify critical parameters capable of reducing overpotentials and safety risks. These parameters include operational (i.e. temperature), structural (i.e. thickness of the components), and material properties (i.e. exchange current density). Simulated results suggest that better engineering of the membrane electrode assemblies (MEA) can offset intrinsic limitations. Membrane thickness and conductivity yielded the largest percent change on performance. Gas crossover was observed to be sensitive to diffusion, solubility, and bubble formation. Current density regimes were identified, where each crossover parameter was most influential. MEA parameterization serves as a framework for theoretical PEMWE optimization simulations and input for design criteria for future research targets.

Dizon, A↗
DOE OSTI2023

System and method for simulation-assisted additive manufacturing involving estimating anisotropic deformation

A simulation-assisted additive manufacturing system and method for generating an anisotropic compensation to account for non-uniform deformation due to additive manufacturing and service loading. The predicted deformation may not be fully defined. The present disclosure provides a system and method for estimating the missing deformation data through regression analysis. The present disclosure also provides an integrated framework where the various simulated-assisted design modules are configured for two-way communication and sharing access to changes to the model.

Kim, Seokpum↗
DOE OSTI2024

Simulation of triggering and evolution of ELM by pellet injection in EAST under BOUT++ framework

A BOUT ++ three-field magnetohydrodynamic model is employed to study the triggering and evolution of edge localized mode (ELM) by Li pellets injected along the outer mid-plane in the EAST configuration. The linear simulation shows that compared with a large deposition on the pedestal top (scenario I), a smaller deposition within the steep-gradient pedestal region (scenario II) can stimulate much larger linear growth rates of all-n peeling-ballooning modes (PBMs). The nonlinear simulation shows that there exists a pellet size threshold for ELM triggering for two deposition locations; the threshold for scenario I predicted in the present study matches the EAST observation well. Comparison of the two scenarios reveals that a smaller deposition is sufficient to trigger an ELM in a much shorter time in scenario II, whose ELM size is comparable to that in scenario I. This conclusion confirms previous DIII-D and ASDEX-Upgrade observations, suggesting that the steep-gradient pedestal region is a favorable deposition location for ELM triggering with minimum pellet size. Simulation analyses also find that the positive radial gradient of the hump-like pressure profile in the outer mid-plane induced by the pellet deposition plays a different role in the two scenarios. In scenario I, the force resulting from the gradient hinders the outflow of core plasmas and in return, the perturbation is suppressed from spreading inwards after ELM crashes. In scenario II, with a sizable deposition, the gradient results in another competitive perturbation growth region during the linear phase, thus dispersing the free energy and reducing the efficiency of destabilizing PBMs by pellet injection. The suppressing effect of saturated zonal flow on other modes, the short ELM fast crash phase, and the restricting transport effect of the positive radial pressure gradient work together to constrain the pedestal energy loss, especially when the pellet deposition amount is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Modeling the Formation and Eruption of Coronal Structures by Linking Data-Driven Magnetofrictional and MHD Simulations for AR 12673

Context .The data-driven and time-dependent modeling of coronal magnetic fields is crucial for understanding solar eruptions. These efforts are complicated by the challenges of finding a balance between physical realism and computing efficiency. One possible technique is to couple two modeling approaches. Aims . Our aim here is to showcase our progress in using time-dependent magnetofrictional model (TMFM) results as input to dynamical magnetohydrodynamic (MHD) simulations. However, due to the different evolution processes in these two models, using TMFM snapshots in an MHD simulation is nontrivial. We address these issues, both physically and numerically, discuss the incompatibility of the TMFM output to serve as the initial condition in MHD simulations, and show our methods of mitigating this. The evolution of the flux systems and the cause of the eruption are investigated. Methods . TMFM is a prevalent approach that has proven to be a very useful tool in the study of the formation of unstable structures in the solar corona. In particular, it is capable of incorporating observational data as initial and boundary conditions and requires shorter computational time compared to MHD simulations. To leverage the efficiency of data-driven TMFM and also to simulate eruptive events in the MHD framework, one can apply TMFM up to a certain time before the expected eruption(s) and then proceed with the simulation in the full or ideal MHD regime in order to more accurately capture the eruption process. Results . We show the results of a benchmark test case with a linked TMFM and MHD simulation to study the evolution of NOAA active region 12673. A rise of a twisted flux bundle through the MHD simulation domain is observed, but we find that the rate of the rise and the altitude reached depends on the time of the TMFM snapshot that was used to initialize the MHD simulation and the helicity injected into the system. The analysis suggested that torus instability and slip-running reconnection could play an important role in the eruption. Conclusions . The results show that the linkage of TMFM and zero-β MHD models can be successfully used to model the eruptive coronal magnetic fields.

F. Daei↗
NASA NTRS2013

Simulation of the 23 July 2012 Extreme Space Weather Event: What if This Extremely Rare CME Was Earth Directed?

Extreme space weather events are known to cause adverse impacts on critical modern day technological infrastructure such as high-voltage electric power transmission grids. On 23 July 2012, NASA's Solar Terrestrial Relations Observatory-Ahead (STEREO-A) spacecraft observed in situ an extremely fast coronal mass ejection (CME) that traveled 0.96 astronomical units (approx. 1 AU) in about 19 h. Here we use the SpaceWeather Modeling Framework (SWMF) to perform a simulation of this rare CME.We consider STEREO-A in situ observations to represent the upstream L1 solar wind boundary conditions. The goal of this study is to examine what would have happened if this Rare-type CME was Earth-bound. Global SWMF-generated ground geomagnetic field perturbations are used to compute the simulated induced geoelectric field at specific ground-based active INTERMAGNET magnetometer sites. Simulation results show that while modeled global SYM-H index, a high-resolution equivalent of the Dst index, was comparable to previously observed severe geomagnetic storms such as the Halloween 2003 storm, the 23 July CME would have produced some of the largest geomagnetically induced electric fields, making it very geoeffective. These results have important practical applications for risk management of electrical power grids.

SPACE WEATHER↗
DOE OSTI2026

Demand response event simulator and risk-aware bidding tool for industrial customers

Incentive Based Demand Response (IBDR) program participation delivers financial benefits to the consumers and resiliency benefits to the electricity grid. Effectively participating in these programs as an industrial consumer requires bidding strategies that balance financial risk with operational constraints. Existing bidding tools tend not to fully incorporate stochastic IBDR event modeling, program specific baseline and payment/penalty calculations, or demand reduction process control schemes that account for the cascading impacts of shutdown in complex facilities. Here, this work presents an IBDR event simulator and risk-aware bidding framework tool integrating three key components: a flexible, parameterized demand response event generator that rigorously accounts for program structures and stochasticity, a demand response operational simulation model that generates explicit control strategies for load reduction, and a Monte Carlo simulator to evaluate financial risk for varied capacity bids. A case study at a wastewater treatment plant participating in PG&E's Capacity Bidding Program demonstrates the framework's utility. In the peak capacity price month of August, optimal bidding by the wastewater treatment plant nets a mean IBDR benefit of $101,000 (67% of the August electricity bill) with 0.4% probability of a financial loss. This framework enables industrial operators to make informed bidding decisions, negotiate better program terms with demand response load aggregators, and analyze energy flexibility investments at their facilities. Ultimately, this work reduces participation barriers in IBDR programs and supports the broader goal of enhancing grid reliability and renewable energy integration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probabilistic Sizing of Energy Storage Systems for Reliability and Frequency Security in Wind-Rich Power Grids

The penetration of wind energy has increased significantly in the power grid in recent times. Although wind is abundant, environment-friendly, and cheap, it is variable in nature and does not contribute to system inertia as much as conventional synchronous generators. Coupled with the low inertia contribution, the generation intermittency of wind power leads to reliability and stability issues in the power system. Energy storage systems (ESSs) are among the most prominent alternatives to alleviate these concerns associated with high wind penetration. This paper proposes a planning strategy to size ESS for the reliability and frequency security of wind-rich power grids. A probabilistic methodology for ESS sizing is developed utilizing a composite reliability-based framework with sequential Monte Carlo simulation (MCS). The MCS generates composite reliability indices for the power system, which are employed to obtain the capacity for a reliability energy storage system (RESS). Simultaneously, the MCS-derived probability of synchronization of conventional generators is integrated into an analytical approach for sizing a frequency support energy storage system (FESS). The effect of wind farm dispersion across geographical regions is incorporated in the framework to study possible reductions in the ESS size while maintaining the system reliability and frequency security. Furthermore, the efficacy of the proposed strategy is demonstrated on the RTS-GMLC test system.

25 ENERGY STORAGE↗
DOE OSTI2024

Model output from simulations of manganese-carbon interactions in temperate forest soil profiles

This archive contains model output, code, and scripts for simulations of coupled manganese-carbon cycling in temperate forest soil profiles. These model results were generated as part of a study investigating how manganese availability influences soil organic carbon stocks and demonstrating a new model framework for coupling carbon and manganese cycling. The simulations were in support of a manuscript: "Modeling interactive effects of manganese bioavailability, nitrogen deposition, and warming on soil carbon storage." The study addresses the research questions: How does Mn bioavailability, as driven by subsurface mineral properties, pH, and redox status, affect temperate forest soil organic carbon and litter carbon stocks?How is the relationship between Mn bioavailability and carbon cycling affected by changes in temperature and nitrogen deposition?"Model simulations were conducted in a reactive transport modeling framework using PFLOTRAN coupled to python. Multiple model simulations testing different Mn-bearing mineral solubilities, hydrological patterns, nitrogen deposition rates, and temperatures are included. Soil properties, including total and exchangeable Mn concentrations, are based on values reported for the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), a temperate forested watershed in central Pennsylvania, U.S.A where Mn cycling through vegetation has been documented.File formats include netCDF (.nc), python script (.py), shell script (.sh), plain text PFLOTRAN input file (.in), and plain text PFLOTRAN database file (.dat), and gzipped tar archive (tar.gz).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A framework for detection and attribution of regional precipitation change: Application to the United States historical record

Abstract Despite the emerging influence of anthropogenic climate change on the global water cycle, at regional scales the combination of observational uncertainty, large internal variability, and modeling uncertainty undermine robust statements regarding the human influence on precipitation. Here, we use output from global climate models in a perfect-data sense to develop a framework for conducting regional detection and attribution (D&A) for precipitation, starting with the contiguous United States (CONUS) where observational uncertainty is lower than in other regions. Our unified approach can simultaneously detect systematic trends in mean and extreme precipitation, attribute trends to anthropogenic forcings, compute the effects of forcings as a function of time, and map the effects of individual forcings. Model output is used to conduct a set of tests that yield a parsimonious representation for characterizing seasonal precipitation over the CONUS for the historical record (1900 to present day), which ensures our D&A is insensitive to structural uncertainty. Our framework is developed using synthetic data in a Pearl-causal perspective wherein causality can be identified using intervention-based simulations. While the hypothesis-based framework and accompanying generalized D&A formula we develop should be widely applicable, we include a strong caution that the hypothesis-guided simplification of the formula for the historical climatic record of CONUS as described in this paper will likely fail to hold in other geographic regions and under future warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A review of methods for the reconstruction and analysis of integrated genome-scale models of metabolism and regulation

The current survey aims to describe the main methodologies for extending the reconstruction and analysis of genome-scale metabolic models and phenotype simulation with Flux Balance Analysis mathematical frameworks, via the integration of Transcriptional Regulatory Networks and/or gene expression data. Although the surveyed methods are aimed at improving phenotype simulations obtained from these models, the perspective of reconstructing integrated genome-scale models of metabolism and gene expression for diverse prokaryotes is still an open challenge.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Rapid Rebuilding of the Outer Radiation Belt

Recent observations by the radiation monitor (RDM) on the spacecraft Akebono have shown several cases of greater than 2.5 MeV radiation belt electron enhancements occurring on timescales of less than a few hours. Similar enhancements are also seen in detectors on board the NOAA/POES and TWINS 1 satellites. These intervals are shorter than typical radial diffusion or wave-particle interactions can account for. We choose two so-called "rapid rebuilding" events that occur during high speed streams (4 September 2008 and 22 July 2009) and simulated them with the Space Weather Modeling Framework configured with global magnetosphere, radiation belt, ring current, and ionosphere electrodynamics model. Our simulations produce a weaker and delayed dipolarization as compared to observations, but the associated inductive electric field in the simulations is still strong enough to rapidly transport and accelerate MeV electrons resulting in an energetic electron flux enhancement that is somewhat weaker than is observed. Nevertheless, the calculated flux enhancement and dipolarization is found to be qualitatively consistent with the observations. Taken together, the modeling results and observations support the conclusion that storm-time dipolarization events in the magnetospheric magnetic field result in strong radial transport and energization of radiation belt electrons.

Glocer, A.↗
DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
DOE OSTI2022

Simulation and Optimization of Volatile Fatty Acid Upgrading Strategies for Sustainable Transportation Fuel Production

Underutilized wet waste is a promising feedstock for production of carbon-neutral or carbon-negative liquid transportation fuels. Arrested methanogenesis of wet wastes by microbes, which produces mixtures of C2-C8 carboxylic acids (volatile fatty acids, VFAs), is practiced industrially at a pilot scale and the VFA products can be catalytically upgraded to molecules suitable for use as transportation fuels (alcohols and alkanes) via sequential ketonization and hydrogenation, both steps of which have been demonstrated with high (>90%) yield at the lab scale. We present in this work a simulation-based decision-making framework which evaluates upgrading strategies of VFAs to alcohol- and hydrocarbon-based neat and blended diesel, jet, and gasoline fuels. The processes utilize sequential upgrading steps (ketonization and hydrogenation) accompanied by one distillation step (either separating VFAs before ketonization or alcohols after ketone hydrogenation) to form two liquid streams: a light alcohol stream best suited as a light-duty fuel and a heavy alcohol or hydrocarbon stream best suited as a diesel or jet fuel. Suitability of the liquid products as transportation fuels and maximum blending levels with petrofuels comes from simulation of critical fuel properties (boiling point, flash point, lower heating value, viscosity, melting point, water solubility, and cetane/octane number). Catalytic upgrading steps and fuel property mixing models used by the simulation are experimentally validated. We demonstrate the flexibility of the decision-making algorithm by evaluating processing scenarios for experimentally-observed VFA profiles with varied carbon chain length distributions. Processing scenarios are optimized for either maximum total renewable carbon utilization or production of heavy-duty (jet or diesel) fuels. We evaluate the economic viability and CO2 emissions of each proposed upgrading scenario using techno-economic and life cycle analyses. The decision-making framework developed in this work can also be used to down-select promising strategies for upgrading of other bio-based feedstocks to sustainable fuels and chemicals.

aviation fuel↗
DOE OSTI2020

Understanding the deformation mechanisms in Ni-based superalloys with using crystal plasticity finite element method

Nl-based superalloy is considered as a good candidate due to its excellent resistance to elevated temperature deformation for long term period application. Understanding the deformation and failure mechanisms of Ni-Based superalloys is very helpful for providing design guidelines for processing Ni-based superalloys. Experimental characterization indicates that the deformation mechanisms of Ni based superalloy is strongly microstructure dependent. Besides, damage transform from the void nucleation to the macro cracks by voids growth leading to the failure of the Ni-based superalloys are also showing strong microstructure sensitivity. Therefore, this work focuses on the prediction and comprehension of the deformation and void growth behavior in Ni based superalloy at different working conditions via crystal plasticity finite element modeling and simulation. Physically based crystal plasticity frameworks were developed for newly Ni- based superalloy Haynes 282. It was found that dislocation shearing through the precipitates were acting as the main contributor to the strength of Haynes 282 at room temperature and 815°C. Our analysis of the creeping behavior of Haynes 282 exhibited that resistance of general climb replaced by the resistance induced by the deposited climb dislocation density. In addition, in the study of void growth behavior, our simulation results demonstrated that as the main loading axis perpendicular to the grain boundary (GB), voids grow more slowly on tilt GBs in bicrystals than those in single and bicrystals with twist GBs. And tilt GBs would promote the void grow into irregular shape.

36 MATERIALS SCIENCE↗