Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling and 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A software package for modeling and simulating fault graphs

This report describes a novel fault graph modeling language and a simulation tool for executing models specified in the language. The modeling language has three primary features that distinguish it from similar reliability analysis tools. These are (1) a random variable modeling several distinct outcomes of a single fault; (2) chains of faults in which one fault triggers another; and (3) time to fail sampled from probability distributions including positive normal, exponential, Weibull with a minimum, or immediate. These features are motivated by their use in a historical analysis of centrifuge reliability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Numerical simulation of the hot-tail runaway electron production mechanism using CQL3D and comparison with Smith–Verwichte analytical model

Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.

Physics↗
DOE OSTI2020

Survey of Modeling and Simulation Techniques for Advanced Manufacturing Technologies Volume I – Predicting Initial Microstructures

This report summarizes the current state of modeling and simulation methods for predicting the initial structure and properties of material in components assembled using advanced manufacturing technologies (AMTs). The report is the first volume in a two-volume series. This first volume focuses on predicting initial microstructures of AM material. The second volume discusses predicting material properties given the initial microstructure. The focus is on technologies of particular relevance to the design and manufacture of nuclear reactor components. The purpose of the report is to help develop the technical knowledge base to support regulatory decisions that will be needed to assess nuclear components manufactured with AMTs as these components are installed at nuclear power plants (NPPs). The report develops a list of AMTs of particular interest to the NRC and summarizes the key microstructural features and corresponding processing parameters relevant to each technology. Further sections describe available physically-based and data-driven prediction methods and survey widely available software tools. The report concludes with a summary of gaps that may be of particular interest to the NRC when evaluating modeling and simulation methods for AMTs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A protocol for model intercomparison of impacts of marine cloud brightening climate intervention

A modeling protocol (defined by a series of climate model simulations with specified model output) is introduced. Studies using these simulations are designed to improve the understanding of climate impacts using a strategy for climate intervention (CI) known as marine cloud brightening (MCB) in specific regions; therefore, the protocol is called MCB-REG (where REG stands for region). The model simulations are not intended to assess consequences of a realistic MCB deployment intended to achieve specific climate targets but instead to expose responses to interventions in six regions with pervasive cloud systems that are often considered candidates for such a deployment. A calibration step involving simulations with fixed sea surface temperatures (SSTs) is first used to identify a common forcing, and then coupled simulations with forcing in individual regions and combinations of regions are used to examine climate impacts. Synthetic estimates constructed by superposing responses from simulations with forcing in individual regions are considered a means of approximating the climate impacts produced when MCB interventions are introduced in multiple regions. A few results comparing simulations from three modern climate models (CESM2, E3SMv2, and UKESM1) are used to illustrate the similarities and differences between model behavior and the utility of estimates of MCB climate responses that were synthesized by summing responses introduced in individual regions. Cloud responses to aerosol injections differ substantially between models (CESM2 clouds appear much more susceptible to aerosol emissions than the other models), but patterns in precipitation and surface temperature responses were similar when forcing is imposed with similar amplitudes in the same regions. A previously identified La Niña-like response to forcing introduced in the Southeast Pacific is evident in this study, but the amplitude of the response was shown to markedly differ across the three models. Other common response patterns were also found and are discussed. Forcing in the Southeast Atlantic consistently (across all three models) produces weaker global cooling than that in other regions, and the Southeast Pacific and South Pacific show the strongest cooling. This indicates that the efficiency of a given intervention depends on not only the susceptibility of the clouds to aerosol perturbations, but also the strength of the underlying radiative feedbacks and ocean responses operating within each region. These responses were generally robust across models, but more studies and an examination of responses with ensembles would be beneficial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Discrete-event Simulation Process Model for the Pyrochemical Processing of Plutonium at Los Alamos National Laboratory

The pyrochemical metal production operations that occur in the Plutonium Facility at Los Alamos National Laboratory perform plutonium purification with the aim to provide plutonium metal for a variety of defense- and non-defense missions within the National Nuclear Security Administration. The demands and constraints associated with the pyrochemical processing of plutonium are complex, making decision analyses challenging for program managers who require plutonium production for their mission applications. The construction of a discrete-event simulation process model is proposed to measure and report the process capacity, material throughput, equipment requirements, and dose accumulation for operators of the pyrochemical metal production operations. The process model, constructed in the ExtendSim™ software, will represent the cause-and-effect relationships between the pyrochemical processing environment and the process constraints, including criticality limitations, material control and accountability measures, chemical analysis requirements, and equipment availability. An accurate representation of the pyrochemical metal production process capacity through simulation modeling will be helpful to program managers in their efforts to forecast plutonium availability for mission applications. Furthermore, the proposed process model will be vital for future analyses that will measure the interactions between the pyrochemical metal production operations and the aqueous reprocessing operations and their ability to minimize transuranic waste disposal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A scale-invariant bonded particle model for simulating large deformation and failure of continua

In this paper, a scale-invariant bond model for simulating two-dimensional large deformation of continua is presented. An accurate torque form caused by relative rotation between bonded particles is derived based on beam theory. Validations of the proposed bond model are carried out by large deformation simulations of a simply supported beam subjected to central loading and a cantilever beam with end loading. Comparison of simulations with analytical solutions show a very good agreement. Later, Brazilian tests are simulated to demonstrate the ability of the framework to model failure and fracture process of rock.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling and Simulation of Automated Mobility Districts

Increasing interest and investment in connected, automated, and electric vehicles (CAEVs), and mobility-as-a-service (MaaS) concepts are paving the way for the next major shift in transportation through automated and shared mobility. The initial excitement towards rapid deployment and adoption of automated vehicles has subsided, and low-speed automated shuttles are emerging as a more pragmatic pathway for introducing automated mobility in geo-fenced districts. Such shuttles hold the promise to provide a viable alternative for serving short trips in urban districts with high travel densities. As interest in low-speed automated shuttle systems (to improve urban mobility) increases, the need for tools that can inform communities regarding benefits or dis-benefits of automated shuttle deployments is imminent. However, most of the existing transportation planning and simulation tools are not capable of handling emerging shared automated mobility options. This presentation provide the details of a microscopic simulation toolkit that can be used by cities and communities to plan for the deployment of low-speed automated shuttles systems, as well as other shared mobility options. Labeled as Automated Mobility District (AMD) modeling and simulation toolkit, the proposed decision support tool intends to help cities evaluate the mobility and sustainability impacts of deploying shared automated vehicles in geofenced regions.The AMD modeling and simulation task is funded through DOE's SMART Mobility Consortium. This presentation outlines the progress of the AMD modeling and simulation task in FY19, including case studies conducted in Greenville, SC, and Austin, TX.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

History Matching and Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope is ongoing. After more than 3 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is less than that of the same polymer used in a polymer pilot in Argentina.

Wang, Xindan↗
DOE OSTI2019

Wind Energy Modeling and Simulation, Volume 1: Atmosphere and Plant

In order to optimise the yield of wind power from existing and future wind plants, the entire breadth of the system of a plant, from the wind field to the turbine components, needs to be modelled in the design process. The modelling and simulation approaches used in each subsystem as well as the system-wide solution methods to optimize across subsystem boundaries are described in this reference. Chapters are written by technical experts in each field, describing the current state of the art in modelling and simulation for wind plant design. This comprehensive, two-volume research reference will provide long-lasting insight into the methods that will need to be developed for the technology to advance into its next generation. Volume 1 covers the computing challenges in full turbine modelling, then discusses bridging scales in the atmosphere and turbulence modelling, wind forecasting, wind plant flow, and plant level controller design.

17 WIND ENERGY↗
DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗
DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Early Career Researcher Support for Foundations of Molecular Modeling and Simulation Conference

This project enabled 14 early career researchers, including graduate students and postdocs, to attend the 8th triennial conference on Foundations of Molecular Modeling and Simulation (FOMMS 2022), which took place in Delavan, WI (Lake Geneva area), from July 17 to 21, 2022. Molecular-level modeling and simulation, which serve as critical tools for advancing science and engineering in a myriad of applications, were at the forefront of discussions. These ranged from the design of new materials for various energy applications to insights into the interaction between drug molecules and protein receptors within the body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Global tropical cyclone precipitation scaling with sea surface temperature

Understanding the relationship between tropical cyclone (TC) precipitation and sea surface temperature (SST) is essential for both TC hazard forecasting and projecting how these hazards will change in the future due to climate change. This work untangles how global TC precipitation is impacted by present-day SST variability (known as apparent scaling) and by long-term changes in SST caused by climate change (known as climate scaling). A variety of datasets are used including precipitation and SST observations, realistic climate model simulations, and idealized climate model simulations. The apparent scaling rates depend on precipitation metric; examples shown here have ranges of 6.1 to 9.5% per K versus 5.9 to 9.8% per K for two different metrics. The climate scaling is estimated at about 5% per K, which is slightly less than the atmospheric moisture scaling based on thermodynamic principles of about 7% per K (i.e., the Clausius–Clapeyron scaling). The apparent scaling is greater than the climate scaling, which implies that the relationship between TC precipitation and present-day SST variability should not be used to project the long-term response of TC precipitation to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Extension of Virtual Test Bed Advanced Modeling and Simulation Capabilities for Fusion Energy

The National Reactor Innovation Center (NRIC) was established to accelerate the deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing reactor experiments. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test beds. It is a collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the mission to accelerate the deployment of advanced reactors by facilitating the adoption of advanced modeling and simulation (M&S) tools developed by the DOE. This mission has been carried out by the VTB since 2020 by hosting and featuring dozens of advanced fission nuclear reactor models developed by national laboratories and academia. The charter of the NRIC’s definition of advanced reactors also includes fusion nuclear reactors. As the tools developed by the NEAMS program are increasingly used for modeling fusion energy devices, there is an increasing need to host fusion reactor models on the VTB repository. The VTB will be extended in 2024 to support fusion energy modeling and simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗
DOE OSTI2024

Capturing the run-in of a pebble-bed reactor by using thermal feedback and high-fidelity neutronics simulations

Modeling the run-in of a pebble-bed reactor (PBR) can be challenging as a result of changes in the power, fuel type, and temperatures that occur throughout the run-in period. Previous work utilized high-fidelity neutronics simulations or lower-fidelity coupled neutronics/thermal-hydraulics models to capture the general characteristics of the run-in process. Here, the present work employs high-fidelity neutronics simulations (using Serpent) coupled with thermal-hydraulics simulations (using Griffin–Pronghorn) to capture the thermal feedback present during the run-in and approach to equilibrium for a PBR. Incorporating thermal feedback enables important distinctions to be made about conditions occurring inside the core, as the power distribution, discharge burnup, and isotopic compositions are all affected by the temperature distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Analysis of permeability change in dissolving rough fractures using depth-averaged flow and reactive transport models

Many geophysical problems and engineering applications involve fluid flow through evolving faults and fractures. These situations are typical in CO 2 injection scenarios where supercritical CO 2 or CO 2 –saturated fluids are able to dissolve the solid matrix. Dissolution or precipitation within a fracture could lead to self-reinforcing or self-limiting behavior due to feedbacks from changes in permeability (k). There is limited predictive understanding of this phenomenon. Here, we investigated the case of fracture dissolution using a two-dimensional depth-averaged reactive transport model. The model simulated an acidic fluid dissolving a three-dimensional rough fracture. We observed different dissolution patterns across a wide range of geometric properties and different Peclet and Damkohler numbers (960 cases total). Despite this variation, a persistent linear increase in k with time (t) roughly manifested once the dissolution front advanced through the fracture. The time rate of change of k (k/t) is roughly constant for all simulations. The modeled rate is consistent with a theoretical value based on the parallel plates model and with experimentally derived values. Finally, the results confirm the common presence of linear time-dependence of fracture k due to self-reinforcement.

58 GEOSCIENCES↗
DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗