Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system 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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Continuum Modeling of Porous Electrodes for Electrochemical Synthesis

Electrochemical synthesis possesses substantial promise to utilize renewable energy sources to power the conversion of abundant feedstocks to value-added commodity chemicals and fuels. Of the potential system architectures for these processes, only systems employing 3-D structured porous electrodes have the capacity to achieve the high rates of conversion necessary for industrial scale. However, the phenomena and environments in these systems are not well understood and are challenging to probe experimentally. Fortunately, continuum modeling is well-suited to rationalize the observed behavior in electrochemical synthesis, as well as to ultimately provide recommendations for guiding the design of next-generation devices and components. In this review, we begin by presenting an historical review of modeling of porous electrode systems, with the aim of showing how past knowledge of macroscale modeling can contribute to the rising challenge of electrochemical synthesis. We then present a detailed overview of the governing physics and assumptions required to simulate porous electrode systems for electrochemical synthesis. Here, leveraging the developed understanding of porous-electrode theory, we survey and discuss the present literature reports on simulating multiscale phenomena in porous electrodes in order to demonstrate their relevance to understanding and improving the performance of devices for electrochemical synthesis. Lastly, we provide our perspectives regarding future directions in the development of models that can most accurately describe and predict the performance of such devices and discuss the best potential applications of future models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Evaluation of a Variable-Aperture Full-Ring SPECT System using Large-Area Pixelated CZT Modules: A Simulation Study for Brain SPECT Applications

Single photon emission computed tomography (SPECT) scanner using cadmium zinc telluride (CZT) offers improved imaging capability over conventional NaI (Tl)-based SPECT scanner. We aim to demonstrate a full-ring SPECT system design with eight large-area CZT detectors that can be used for a broad spectrum of SPECT radiopharmaceuticals and provides higher sensitivity and better spatial resolution than those of conventional NaI (Tl)-based gamma cameras. A newly-designed full-ring SPECT system is composed of 8 large-area CZT cameras (128 mm × 179.2 mm effective area) that can be independently swiveled around their own axes of rotation independently and can have radial motion for varying aperture sizes that can be adapted to different sizes of imaging volume. Extended projection data were generated by conjoining projections of two adjacent detectors to overcome the limited field-of-view (FOV) by each CZT camera. Using Monte Carlo simulations, we evaluated this new system design with digital phantoms including a Derenzo hot-rod phantom and a Zubal brain phantom. Comparison of performance metrics such as spatial resolution, sensitivity, contrast-to-noise ratio (CNR), and contrast-recovery ratio were made between our design and a conventional SPECT scanner. The proposed scanner could result in up to about 3 times faster in acquisition time over conventional scan time at same acquisition time per step. The spatial resolution of our proposed scanner was similar or better to that of the conventional scanner, and there was significant performance improvement over the conventional scanner particularly in sensitivity (approximately 4 times). Overall, we successfully reconstructed the phantom image for both 99 mTc-based perfusion and 123 I-based dopamine transporter (DaT) brain studies simulated for our new design. In particular, the striatal/background contrast-recovery ratio in 3-to-1 reference ratio was over the 0.8 for the 123 I-based DaT study. In conclusion, we demonstrated the potential of our new full-ring CZT SPECT design, showing improved performance metrics such as improved system sensitivity and CNR while maintaining other important imaging parameters such as spatial resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2026

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W↗
DOE OSTI2022

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

2021 Q1 Project Report: Optimized Bifacial PV Systems (Q1 FY2021 Project Report)

This project has four main technical objectives: 1) Develop and improve bifacial performance models by adding the capability to evaluate electrical behavior and performance of bifacial modules and arrays under realistic field conditions including irradiance variability caused by racking, module frame, and position in the array. 2) Instrument and monitor performance of fielded bifacial systems to validate performance models and to measure, analyze and publish on bifacial energy gain. These should include both research and commercial bifacial systems and cover a variety of deployment applications. 3) Evaluate optimal bifacial system designs using simulations leveraging high-performance computing, and also using full sized and miniaturized experimental field deployments. 4) Establish and contribute to international test standards for bifacial system performance, testing, and safety, and work with the community to establish installation and siting best practices.

42 ENGINEERING↗
DOE OSTI2025

A Greening Future Elevates Flash Drought Risk in Northern Mid‐to‐High Latitudes

Flash droughts have become a growing concern, as they can emerge rapidly and increase the risk of crop failure. Although past studies have investigated the meteorological drivers and future changes of flash drought, why flash drought is more frequent over humid and vegetated regions remains underexplored. This study delves further into the mechanism by which vegetation regulates flash drought and its future change using observations from multiple data sets and large ensemble simulations from three Earth system models. On an interannual timescale, both observations and simulations show robust increases in flash drought frequency and a higher flash-to-sub-seasonal drought ratio during spring or antecedent conditions with dense vegetation, supporting the important role of vegetation in flash drought occurrence, especially in the northern mid-to-high latitudes. In the latter regions, the large ensemble simulations show robust increases in flash drought (e.g., 67% and 46% increases in Eastern U.S. and North Asia in 2050–2100 relative to 1950–2000 under the high emission scenario), where the growing season is lengthening. Although greening might suggest reduced drought stress, it drives precipitation-soil moisture-evapotranspiration decoupling by increasing evapotranspiration partitioning to transpiration. As transpiration can access deep soil water through the plant root system, its increased portion can weaken the constraints of concurrent precipitation on evapotranspiration, thus accelerating soil moisture depletion under high evaporative demand, driving a slow-to-rapid drought transition. How vegetation regulates flash drought by regulating surface moisture budget is supported by observations and simulations. Although warming supports early planting, agriculture may increasingly be threatened by surging flash drought risk.

Drought↗
DOE OSTI2023

Final Cooling with Thick Wedges for a Muon Collider

A muon collider is a particle accelerator that collides muons rather than protons or electrons. An important component of such a collider is cooling, in which the emittance of the muon beam is reduced to make it suitable for use in a collider. In the last stage of this process, called final 4D cooling, emittance in the transverse axes is reduced while emittance in the longitudinal axis is allowed to grow. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on single thick wedges. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We characterized the effects of various parameters on the performance of this system, including the initial emittances and Twiss parameters, the length and angle of the wedges, the length of the drift channel, and the length, frequency, gradient, and phase of the RF cavities. We optimized the design parameters of the cooling channel and produced two conceptual designs (corresponding to two possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ≈ 3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

High-fidelity dimer excitations using quantum hardware

The quantum simulation of entangled spin systems can play a central role in quantum magnetic materials discovery. Additionally, the simulation of spectroscopic signatures, such as the dynamical structure factor accessed in inelastic neutron scattering (INS), necessitates a long timescale for circuit evolution. This is because the energy resolution is directly related to the time over which the circuit could be meaningfully evolved. However, canonical Trotterization requires deep circuits precluding such long-time evolution—even for a small number of qubits. Here, in this study, we demonstrate “direct” resource efficient fast-forwarding (REFF) measurements with short-depth circuits that can be used to capture longer time dynamics of spin Hamiltonians. We showcase the results of the dynamics of a quantum spin dimer, the basic quantum unit of emergent many-body spin systems, whose density of states we simulate accurately. The long temporal evolution and measurement of the two-spin correlation functions enable the calculation of the dynamical structure factor S⁡(Q = 0, ω) measured in the neutron scattering cross-section. We exhibit the clarity of the triplet gap and the triplet splitting of the quantum dimer with class-leading fidelity that enables comparison to experimental neutron data. Our results on current circuit hardware outline an important workflow to predict and benchmark against the outputs of INS experiments of quantum magnets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Minimal implicit-solvent coarse-grained simulation of Pluronic block copolymers with ionic liquids

Pluronic block copolymers, composed of poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) in a triblock structure (PEO–PPO–PEO), are well known for their amphiphilic character and ability to self‐assemble into micelles in aqueous solution. The addition of ionic liquids (ILs) can further modulate the core–shell structures of these copolymers, influencing their stability, critical micellization temperature, and size. However, fully atomistic simulations often become prohibitively expensive due to the size and complexity of these systems. In this work, coarse‐grained simulations using a minimal implicit‐solvent model were performed to examine how two classes of ILs, namely, 1‐alkyl‐3‐methylimidazolium ([C n C 1 im]) and 1‐alkyl‐3‐methylpyrrolidinium ([C n C 1 pyrr]), change the micellization of Pluronic block copolymers in aqueous solution. The effects of IL concentration and alkyl group length were investigated, and the model greatly improved the efficiency of simulating large‐scale micelle systems. Furthermore, the numerical simulations are qualitatively compared with experimental investigations. Our results show that adding ILs expands the micelle core by embedding IL tails among the PPO blocks, thereby increasing overall micelle size. Less polar ILs generally induce more pronounced micellar growth. However, the effect of IL tail length on conformation and micellar packing is non‐monotonic. Up to moderate chain lengths (around C8–C10), the IL tails can extend sufficiently to increase local separation within the micelle; at longer tail lengths, enhanced hydrophobic clustering and steric hindrance cause the tails to bend or fold, capping further expansion. In addition, although block copolymer chains tend to pack more closely in the presence of longer‐tailed ILs, the random coil size of an individual polymer chain does not necessarily shrink. Meanwhile, these insights provide a deeper understanding of how Pluronic/IL systems interact, informing applications in drug delivery, cosmetics, food, and environmental engineering. Finally, our minimal implicit‐solvent model can be applied to larger systems and longer timescales, substantially reducing computational cost while reproducing key structural trends observed experimentally.

Atomistic simulations↗
DOE OSTI2020

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Large eddy simulation of a pressurized reactive fluid–particle system: single-particle dynamics analysis

Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new combustion technologies. Because of the high pressure, the volume fraction of particles in a pressurized oxy-combustor could be higher than that in the conventional combustor, therefore, the particle dynamics in the pressurized vessel need to be evaluated accordingly. The present computational work is to predict the influence of the particle size and particle injection location on the particle trajectory. The Stokes number based on turbulence fluctuation, St, is adopted to evaluate the impact of the gas phase on the particulate phase. A large eddy simulation (LES) with the environment pressure of 15 bar and the thermal input of the fuel particle of 100 kWth has been performed by means of the commercial CFD package. It is shown that that in this pressurized environment St << 1 for the 25 µm particles, St ~ 1 for the 100 µm particles, while for the particles exceeding 200 µm St >> 1. The trajectories of the sampling particles from the LES results support this conclusion: the particles less than 100 µm tend to be uniformly distributed in the space while the particles exceeding 200 µm tend to concentrate in the combustor.

Li, Lei↗
DOE OSTI2021

The IDAES process modeling framework and model library—Flexibility for process simulation and optimization

Abstract Energy systems and manufacturing processes of the 21st century are becoming increasingly dynamic and interconnected, which require new capabilities to effectively model and optimize their design and operations. Such next generation computational tools must leverage state‐of‐the‐art techniques in optimization and be able to rapidly incorporate new advances. To address these requirements, we have developed the Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform, which builds on the strengths of both process simulators (model libraries) and algebraic modeling languages (advanced solvers). This paper specifically presents the IDAES Core Modeling Framework (IDAES‐CMF), along with a case study demonstrating the application of the framework to solve process optimization problems. Capabilities provided by this framework include a flexible, modifiable, open‐source platform for optimization of process flowsheets utilizing state‐of‐the‐art solvers and solution techniques, fully open and extensible libraries of dynamic unit operations models and thermophysical property models, and integrated support for superstructure‐based conceptual design and optimization under uncertainty.

42 ENGINEERING↗
DOE OSTI2026

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗
DOE OSTI2024

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulation and Sensitivity Analysis of a Cryogenic Carbon Capture System for the Cement Industry

The talk will describe modeling and analysis of cryogenic systems for carbon capture from industrial point sources, thermodynamic and process simulation requirements to predict capture from typical flue gases, quantification of temperature and flow profiles for a range of gas feed compositions and capture levels, and impacts of process conditions and heat transfer on energy requirements for heat integration.

Paul, Brandon↗
DOE OSTI2021

Growth of Mesoscale Convective Systems in Observations and a Seasonal Convection-Permitting Simulation over Argentina

Abstract A 6.5-month, convection-permitting simulation is conducted over Argentina covering the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) field campaign and is compared with observations to evaluate mesoscale convective system (MCS) growth prediction. Observed and simulated MCSs are consistently identified, tracked, and separated into growth, mature, and decay stages using top-of-the-atmosphere infrared brightness temperature and surface rainfall. Simulated MCS number, lifetime, seasonal and diurnal cycles, and various cloud-shield characteristics including growth rate are similar to those observed. However, the simulation produces smaller rainfall areas, greater proportions of heavy rainfall, and faster system propagations. Rainfall area is significantly underestimated for long-lived MCSs but not for shorter-lived MCSs, and rain rates are always overestimated. These differences result from a combination of model and satellite retrieval biases, in which simulated MCS rain rates are shifted from light to heavy, while satellite-retrieved rainfall is too frequent relative to rain gauge estimates. However, the simulation reproduces satellite-retrieved MCS cloud-shield evolution well, supporting its usage to examine environmental controls on MCS growth. MCS initiation locations are associated with removal of convective inhibition more than maximized low-level moisture convergence or instability. Rapid growth is associated with a stronger upper-level jet (ULJ) and a deeper northwestern Argentinean low that causes a stronger northerly low-level jet (LLJ), increasing heat and moisture fluxes, low-level vertical wind shear, baroclinicity, and instability. Sustained growth corresponds to similar LLJ, baroclinicity, and instability conditions but is less sensitive to the ULJ, large-scale vertical motion, or low-level shear. Growth sustenance controls MCS maximum extent more than growth rate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A framework for calibrating and validating an air loop dynamic model in an HVAC system in Modelica

The use of Modelica for simulating the dynamic behaviors of building heating, ventilation, and air conditioning (HVAC) systems has gained popularity. Calibration of a model that represents large and complex HVAC systems in Modelica involves the determination of hundreds of parameters using real-world operation data. Considering the coupling effects among various components, such a calibration process is complex and time-consuming. In this study, we propose a systematic framework to efficiently calibrate and validate a complex HVAC system model in Modelica. The framework includes strategies to solicit real-system operation data and to decouple the model. The goal of this calibration framework is to accurately and efficiently determine a set of Modelica model parameters that provide a good match between the simulated and real system behaviors. To demonstrate the validity of the framework, it was applied to calibrate the parameters of a Modelica model of the air loop subsystem of a real AHU-VAV system. The results show that the calibrated model can generate simulated results, including VAV air flow rate, outdoor air flow, and fan power consumption, that match well with the operational data of the real system. The coefficient of the variation of the root mean square error, CV(RMSE), of the air flow rate and power consumption are 14.6 % and 11.5 %, respectively. The results prove that the framework is valid and effective, and it can be used to calibrate other complex HVAC system Modelica models in the future.

Li, Yicheng↗