Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced 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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING (Presentation)

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Fuel Fabrication and Single Stage Aqueous Process Modeling

The Material Protection, Accounting, and Control Technologies program utilizes modeling and simulation to assess Material Control and Accountability (MC&A) concerns for a variety of nuclear facilities. Single analyst tools allow for rapid design and evaluation of advanced approaches for new and existing nuclear facilities. A low enriched uranium (LEU) fuel conversion and fabrication facility simulator is developed to assist with MC&A for existing facilities. Measurements are added to the model (consistent with current best practices). Material balance calculations and statistical tests are also added to the model. In addition, scoping work is performed for developing a single stage aqueous reprocessing model. Preliminary results are presented and discussed, and next steps outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

JSPEC: A Program for IBS and Electron Cooling Simulation

JSPEC (JLab Simulation Package on Electron Cooling) is an open-source C++ program developed at Jefferson Lab to simulate the evolution of the ion beam under the intrabeam scattering effect and/or the electron cooling effect. JSPEC includes various models of the ion beam, the electron beam, and the friction force, aiming to reflect the latest advances in the field and to provide a useful tool to the community. JSPEC has been benchmarked against other cooling simulation codes and experimental data. It has been used to support the cooler design for JLEIC, an earlier JLab design for the Electron-Ion Collider. A Python wrapper of the C++ code, pyJSPEC, for Python 3.x environment has also been developed and released. It allows users to run JSPEC simulations in a Python environment and makes it possible for JSPEC to collaborate with other accelerator and beam modeling programs, as well as plentiful Python tools in data visualization, optimization, machine learning, etc. In this report, we introduce the features of JSPEC, with a focus on the latest development, and demonstrate how to use JSPEC and pyJSPEC with sample codes and numerical examples.

Accelerator Physics↗
DOE OSTI2022

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1979

Developments in the computation of turbulent boundary layers

Computational techniques applicable to turbulent boundary layers are classified into solutions of Reynolds-averaged equations, in which all the effects of the turbulence are modelled, and solutions of three-dimensional, time dependent Navier-Stokes equations, in which the large eddies are calculated and only the turbulence at scales smaller than the computational mesh spacings has to be modelled. Current computation costs place engineering computations in the first of these categories; large eddy simulations are appropriate currently for special studies of the dynamical processes of turbulence in idealized flow fields. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparisons with experimental data for a variety of flow conditions. An example of a large-eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate-of-strain correlations and dissipation in the Reynolds-averaged equations.

Rubesin, M. W.↗
DOE OSTI2021

Integration of Storage in the DC Link of a Full Converter-Based Distributed Wind Turbine

Energy storage is known to support the dispatchability of variable renewable resources. In this paper, we model a battery energy storage system (BESS) integrated with the DC link of a Type IV full converter-based wind turbine and the necessary controls to achieve efficient dispatch. To support the validation of control methodologies, we build a detailed model of a Type IV research wind turbine at the National Renewable Energy Laboratory (NREL), the Controls Advanced Research Turbine (CART 3), and we integrated a lithium-ion BESS model in grid-following mode into the model. The simulation results illustrate the sizing and control of the DC link-integrated BESS for a given variable wind resource and varying dispatch strategies (i.e., under constant, uncertain, and ramping wind scenarios). The integrated storage can smooth variabilities in distributed wind output, hedge against uncertainties, provide the ramping capability, as well as support stability under voltage and frequency transients. All of these have been illustrated in MATLAB/Simulink.

DC-link voltage↗
DOE OSTI2023

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗
DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
DOE OSTI2016

Influence of impact conditions on plasma generation during hypervelocity impact by aluminum projectile

For describing hypervelocity impact (relative low-speed as related to space debris and much lower than travelling speed of meteoroids) phenomenon associated with plasma generation, a self-developed 3D code was advanced to numerically simulate projectiles impacting on a rigid wall. The numerical results were combined with a new ionization model which was developed in an early study to calculate the ionized materials during the impact. The calculated results of ionization were compared with the empirical formulas concluded by experiments in references and a good agreement was obtained. Then based on the reliable 3D numerical code, a series of impacts with different projectile configurations were simulated to investigate the influence of impact conditions on hypervelocity impact generated plasma. It was found that the form of empirical formula needed to be modified. A new empirical formula with a critical impact velocity was advanced to describe the velocity dependence of plasma generation and the parameters of the modified formula were ensured by the comparison between the numerical predictions and the empirical formulas. For different projectile configurations, the changes of plasma charges with time are different but the integrals of charges on time almost stayed in the same level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ion temperature and rotation fluctuation measurements with ultra-fast charge exchange recombination spectroscopy (UF-CHERS) in the DIII-D tokamak

An upgraded detector and several optimizations have significantly improved the Ultra-Fast Charge Exchange Recombination Spectroscopy (UF-CHERS) diagnostic sensitivity to ion temperature and parallel velocity fluctuations at turbulence relevant spatio-temporal scales. Normalized broadband ion temperature and parallel velocity fluctuations down to x̃x∼1% (x = Ti, v∥) and up to ∼450 kHz have been measured in a variety of plasmas. The multi-field nature of the CHERS technique also allows measurements of the cross-phase angles of the fluctuating fields. UF-CHERS is optimized to observe emissions from the electron exchange reaction between intrinsic C6+ and hydrogenic neutral beam injected particles near 529 nm. UF-CHERS consists of two chords separated by ∼1 cm radially, less than the turbulence correlation length in DIII-D plasmas, which enables correlated measurements to suppress incoherent electronic and photon noise. The optical components of the spectrometer include a volume-phase-holographic grating with >90% transmission between 528 and 530 nm and f/2 200-mm lenses, selected to maximize the optical efficiency and photon flux. Diffracted light from each chord is collected in eight spectral bins, each with a bandwidth of ∼0.25 nm, and detected and amplified by chilled avalanche photodiodes and custom high-gain, wide bandwidth low-noise preamplifiers to achieve the optimal signal-to-noise ratio. The resulting signals are digitized at 1 MHz, 103–104× faster than the conventional CHERS diagnostics. Spatial coverage is achieved by repositioning a motorized fiber tray between plasmas. UF-CHERS measurements will advance the understanding of turbulent ion transport and contribute to the validation of transport models and simulations.

Truong, D. D. (ORCID:0000000285732539)↗
DOE OSTI2020

Advances in the simulation of space-charge effects

The nonlinear space-charge effects are an important topic in high intensity accelerator beam dynamics and have been extensively studied using macroparticle tracking simulations. Here, we report on recent advances in the simulation of space-charge effects using a symplectic multiparticle space-charge model. The transverse space-charge limit was explored for a periodic focusing and defocusing lattice. The artificial numerical emittance growth in the macroparticle space-charge simulation was analyzed using a one-dimensional model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Commercial Building Energy Saver, API (CBES API) v2.0

The CBES API provides Application Programming Interface to a suite of functions to improve energy efficiency of buildings, including building energy benchmarking, preliminary retrofit analysis using a pre-simulation database DEEP, and detailed retrofit analysis using energy modeling with the EnergyPlus simulation engine. New capabilities in version 2.0 include renewable technologies (PV, storage), advanced HVAC systems, demand response measures, TDV (Time Dependent Valuation energy) metrics, as well as shading from adjacent buildings, trees or roof obstructions. The CBES API is used to power the LBNL CBES Web App. It can be used by third party developers and vendors.

Hong, Tianzhen↗
DOE OSTI2016

Low Cost, High Capacity Regenerable Sorbent for Carbon Dioxide Capture from Existing Coal-fired Power Plants

In this project TDA Research, Inc (TDA) has developed a new post combustion carbon capture technology based on a vacuum swing adsorption system that uses a steam purge and demonstrated its technical feasibility and economic viability in laboratory-scale tests and tests in actual coal derived flue gas. TDA uses an advanced physical adsorbent to selectively remove CO 2 from the flue gas. The sorbent exhibits a much higher affinity for CO 2 than N 2 , H 2 O or O 2 , enabling effective CO 2 separation from the flue gas. We also carried out a detailed process design and analysis of the new system as part of both sub-critical and super-critical pulverized coal fired power plants. The new technology uses a low cost, high capacity adsorbent that selectively removes CO 2 in the presence of moisture at the flue gas temperature without a need for significant cooling of the flue gas or moisture removal. The sorbent is based on a TDA proprietary mesoporous carbon that consists of surface functionalized groups that remove CO 2 via physical adsorption. The high surface area and favorable porosity of the sorbent also provides a unique platform to introduce additional functionality, such as active groups to remove trace metals (e.g., Hg, As). In collaboration with the Advanced Power and Energy Program of the University of California, Irvine (UCI), TDA developed system simulation models using Aspen PlusTM simulation software to assess the economic viability of TDA’s VSA-based post-combustion carbon capture technology. The levelized cost of electricity including the TS&M costs for CO 2 is calculated as $116.71/MWh and $113.76/MWh for TDA system integrated with sub-critical and super-critical pulverized coal fired power plants; much lower than the $153.03/MWhand $147.44/MWh calculated for the corresponding amine based systems. The cost of CO 2 captured for TDA’s VSA based system is $38.90 and $39.71 per tonne compared to $65.46 and $66.56 per tonne for amine based system on 2011 $ basis, providing 40% lower cost of CO 2 captured. In this analysis we have used a sorbent life of 4 years. If a longer sorbent life can be maintained (which is not unreasonable for fixed bed commercial PSA systems), this would lower the cost of CO 2 captured by $0.05 per tonne (e.g., to $38.85 and $39.66 per tonne at 5 years sorbent replacement). These system analysis results suggest that TDA’s VSA-based post-combustion capture technology can substantially improve the power plant’s thermal performance while achieving near zero emissions, including greater than 90% carbon capture. The higher net plant efficiency and lower capital and operating costs results in a substantial reduction in the cost of carbon capture and cost of electricity for the power plant equipped with TDA’s technology.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

CFD Modeling of Concentration Polarization in High-Performance Flat-Sheet Modules

Membrane-based gas separation is a promising, regeneration-free route for industrial CO₂ capture, but its performance is often limited by concentration polarization (CP). Classical 1-D models neglect CP and thus overpredict fluxes, especially for today’s high-permeance, high-selectivity membranes. Therefore, we built a 3-D CFD model that resolves local species transport, boundary-layer resistance, permeate-side effects, and pressure drop in flat-sheet modules. Using the Concentration Polarization Coefficient (CPC) and the overall Effective Membrane CP Index (EAC), the model quantifies how geometry and operating conditions shape CP intensity. Simulations show CP becomes pronounced for advanced membranes, with feed-channel height emerging as the dominant geometric lever: shrinking the gap mitigates CP but raises pressure drop, forcing a trade-off between mass-transfer enhancement and energy cost. These insights emphasize the need for CFD-level analysis in next-generation module design, where 1-D treatments can severely underestimate required membrane area and misguide scale-up.

carbon capture↗
DOE OSTI2026

TChem-atm (v2.0.0): scalable performance-portable multiphase atmospheric chemistry

We present TChem-atm, a performance-portable approach that enables efficient simulation of chemically detailed and multiphase atmospheric chemistry on modern heterogeneous computing architectures. Unlike previous efforts that rely on architecture-specific code or focus exclusively on gas-phase chemistry, TChem-atm supports fully coupled gas–aerosol systems with execution across CPUs, NVIDIA GPUs, and AMD GPUs through the Kokkos programming model. It integrates the flexible multiphase capabilities of the Community Atmospheric Model Chemistry Package (CAMP) with the high-performance kinetic routines of TChem, and includes automatic Jacobian construction with support for a range of stiff ODE solvers. In a proof-of-concept integration with the particle-resolved model PartMC, TChem-atm reproduces the existing PartMC–CAMP implementation within solver tolerances and delivers substantial GPU speedups, especially for large particle populations. Performance benchmarks reveal substantial speedups on GPU platforms, particularly for large particle populations, with consistent results across hardware backends. TChem-atm enables performance-portable execution across CPUs and GPUs, though optimal efficiency may require modest architecture-specific tuning (e.g., team and vector sizes), with up to a twofold improvement on the NVIDIA H100. It directly supports sectional and particle-resolved host models, while modal aerosol schemes require minor adaptation to provide particle-scale quantities such as representative diameters. By enabling chemically detailed, multiphase simulations with performance portability and host-model flexibility, TChem-atm facilitates the incorporation of advanced chemistry into atmospheric models.

Díaz-Ibarra, Oscar Homero [Sandia National Laborat↗
NASA NTRS1985

SSME fuel preburner phase 1/phase 2+ comparison

An analytical comparison of the predicted combustion flowfields for the Phase 1 and Phase 2+ fuel preburner configurations using an advanced CFD combustion code is currently underway at Rocketdyne. The Phase 2+ injector has a modified element design. The redesigned element provides for improved atomization through a shift in the injector element/face bleed flow split while maintaining an acceptable pressure drop. Cold-flow tests of the new design have been encouraging. As a complement to those tests, the Advanced Rocket Injector Combustor Code (ARICC) was selected to conduct a detailed combustion flowfield simulation. ARICC models two-dimensional (or axisymmetric), transient, turbulent, two-phase mixing and combustion flowfields. It is derived from a member of the Los Alamos ICEd-ALE family of codes. Unique features of ARICC include explicit representation of the coaxial LOX jet and hydrogen gas flows, distributed atomization processes, droplet breakup, and supercritical vaporation processes. Predictions of temperature and OH concentration profiles for the Phase 1 and Phase 2+ injector flowfields indicate a trend toward more uniform temperature distributions and shorter flame lengths with the Phase 2+ design.

Jensen, R. J.↗