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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Total thermal neutron cross section measurements of yttrium hydride from 0.0005 - 3 eV

Yttrium hydride serves as a neutron moderator material that enables compact, high temperature nuclear reactors. However, in order to accurately design and simulate a nuclear system relying upon yttrium hydride, the fundamental nuclear data of yttrium hydride must be well understood. Thermal neutron scattering law (TSL) evaluations represent an important aspect of nuclear data as thermal scattering can drastically alter the neutron multiplication factor of a system. Therefore, to support evaluation and validation of thermal neutron scattering for yttrium hydride, researchers at Rensselaer Polytechnic Institute (RPI) performed total thermal neutron cross section measurements for YH 1.68 and YH 1.85 over the energy range of 0.0005 - 3 eV. Further, these measurements represent the first total cross section measurements for yttrium hydride that encompass the entire thermal region. Comparisons were made against the ENDF-B/VIII.0, Zerkle & Holmes and Oak Ridge National Laboratory TSL evaluations, where generally good agreement was found.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Mapping and Characterizing the Visual Impacts of the Existing US Wind Turbine Fleet

Visual impacts of wind turbines have been a persistent concern for wind energy development in the United States (US) for decades and remain a major source of project delays and cancellations. Assessments of visual impacts are frequently performed at a local scale for individual projects, but a comprehensive understanding of broader geographic patterns in visual impacts across the US is lacking. This paper presents a visual impact assessment of the existing land-based wind turbine fleet of the contiguous United States (CONUS). The assessment combines geographic information systems and 3D simulation methods to account for key factors driving the visual magnitude of impacts from the installed turbines. The results indicate that, despite the deployment of approximately 70,000 turbines and over 144 gigawatts of land-based wind in the CONUS, the visual impacts are very small when measured as a proportion of land area, population, and sensitive visual resources. Nonetheless, visual impacts are not evenly distributed: people experience a concentrated share in a small number of natural settings, primarily including plains, prairies, and deserts. Finally, we find that although increased density of wind development consistently leads to visual impacts across a greater proportion of land, it does not always lead to impacts to a greater share of the population. These findings suggest that visual impacts from wind energy are generally well-mitigated across the CONUS to date but also highlight the need for a deeper understanding of landscape sensitivity and individual perceptions of wind turbines in the most heavily impacted natural settings.

3D simulation↗
DOE OSTI2022

Blending and cavern formation within non-Newtonian fluids in stirred tanks: Application to nuclear waste fluid processing

A numerical approach for predicting the time-accurate fluid flow and mixing properties of non-Newtonian fluids is presented. This approach, which is based on the lattice Boltzmann method, is used to characterize the blending in shear-thinning fluids, as well as cavern formation in yield-stress fluids. Predictions compare favorably to measured data and expectations from first-principles theory. Importantly, across the range of fluids and systems here, the simulations require no reparameterization or retuning between scenarios. This generality is exploited to model a nuclear waste processing unit operation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Methodology for physics-informed generation of synthetic neutron time-of-flight measurement data

Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Regional suitability of heating tower heat pumps based on technical and economic performance maps

The heating tower heat pump (HTHP) technology is evolving rapidly due to its high performance in summer and frost-free characteristic in winter. However, lacking comprehensive comparisons between the HTHP and conventional methods, such as air-source heat pump (ASHP), and chiller and boiler, greatly limits its applications. To address this problem, this paper carries out technological and economic performance analyses of the three systems in 353 locations over China by the processes of building load calculation, system sizing, annual simulation, performance comparison, and demonstration. Energy saving principle and potential of the HTHP are presented and explained. Based on that, economic performance comparison is conducted by considering the energy conversion efficiency and prices. The results show that the HTHP has advantage over the ASHP in most locations in China. The improvement in annual coefficient of performance (COP) is about 10% to 20%. For the locations in southwest China like Chongqing, the improvement is also less than 10% due to the large regeneration penalization in winter and small difference between dry-bulb and wet-bulb temperature in summer. The annual COP of the HTHP is 10% to 23% higher than that of the chiller and boiler in hot summer and cold winter zone, and cold zone. However, the HTHP shows no advantage over the chiller and boiler in severe cold zone because of the low efficiency of heat pump in winter, neither in hot summer and warm winter zone due to the small heating load. Based on the results, technical and economic performance maps are developed to guide the application of the HTHP in different climate regions.

heating tower heat pump, air-source heat pump, chi↗
DOE OSTI2026

Quantum computing approach for building surface sunlit in urban-scale energy modeling

Solar shadow calculations are needed in building energy modeling and performance simulation of PV systems installed on roofs or facades of buildings. We present a quantum computing approach for calculation of building surface sunlit fractions by recasting solar visibility as a binary optimization problem solved by quantum annealing. Each triangulated surface centroid is encoded as a binary qubit indicating sunlit or shaded status. Geometric visibility constraints are derived from the Möller-Trumbore intersection algorithm and converted into a constrained quadratic binary model compatible with contemporary quantum annealers. The coefficients were embedded to D-Wave quantum computer. To demonstrate feasibility, we conducted a case study in San Francisco for a target building with 52 triangles and roughly 2700 nearby triangles within 50 m evaluated at representative winter and summer solar positions. The results demonstrated that quantum annealing can reliably calculate and distinguish sunlit from shaded surfaces. Quantum samples achieved average accuracy exceeding 92.4 %, with the aggregate surface-level agreement approaching 99.9 %. The outputs of quantum computers agreed closely with classical algorithms, indicating practical feasibility and promising scalability. Finally, the hourly sunlit fractions of building surfaces can be obtained for urban energy modelling. This is the first study to apply quantum computing to the solar shadow and building surface sunlit calculation. It introduces a new paradigm that differs fundamentally from traditional approaches.

Deng, Zhipeng↗
DOE OSTI2021

Model-free control for frequency response support in microgrids utilizing wind turbines

Islanded microgrids integrated with renewable energy resources bring additional challenges to the primary frequency control due to their power-electronic interfaced generations. Developing an efficient control strategy for inertia emulation and frequency response support of microgrids is becoming an important issue. Wind turbine dynamics, however, will degrade the emulated inertial and primary responses with most traditional control strategies. On the other hand, compensation methods often require the knowledge of system models at the design stage and state estimation for feedback once deployed online. To overcome these challenges, this paper explores a model-free control (MFC) strategy to emulate the desired inertia and support the frequency response of a diesel-wind microgrid system. Additionally, the proposed MFC strategy utilizes a model reference concept to obtain guaranteed inertial response and improve the frequency response of microgrids. This work considers different control strategies with different operating points. The proposed controller is implemented and verified using the modified 33-bus full nonlinear model of the three-phase diesel-wind system in Simulink. Simulation results demonstrate the improved frequency response by precisely emulating the desired inertia using the MFC strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

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

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗
DOE OSTI2024

Mechanistic definition and prediction of the mass exchange coefficient between rivers and hyporheic zones: The $α$ of two $Ω$s

Solute transport in interconnected rivers and hyporheic zones is typically modeled through dual-domain models where first-order solute mass transfer between the two domains, Ω R and Ω HZ , is represented by a coefficient α. The transient storage model (TSM) is an example of such an approach. In practice, α is determined by fitting the tails of solute tracer breakthrough curves using a TSM. This approach has led to ambiguity regarding α’s physical meaning and transferability. Here, in this work, we investigated the physical basis for α and tested it with virtual experiments through the fully coupled multiphysics model hyporheicFoam for the Ω R – Ω HZ system. hyporheicFoam explicitly simulated coupled flow and solute transport over a kilometer with centimeter-scale resolution. Model results were analyzed to calculate α following its theoretical definition directly. Using the determined α within a TSM enables accurate reproduction of solute transport, underscoring α’s physical relevance and precision.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effective and asymptotic criticality of structurally disordered magnets

Changes in magnetic critical behaviour of quenched structurally-disordered magnets are usually exemplified in experiments and in MC simulations by diluted systems consisting of magnetic and non-magnetic components. In our study we aim to show that similar effects can be observed not only for diluted magnets with non-magnetic impurities but may be implemented, e.g., by the presence of two (and more) chemically different magnetic components as well. Therefore we consider a model of the structurally-disordered quenched magnet where all lattice sites are occupied by Ising-like spins of different lengths L. In such a random spin length Ising model, the length L of each spin is a random variable governed by the distribution function p (L). We demonstrate that this model belongs to the universality class of the site-diluted Ising model. This proves that both models are described by the same values of asymptotic critical exponents. However, their effective critical behaviour differs. As a case study, we consider a quenched mixture of two different magnets with values of elementary magnetic moments L 1 = 1 and L 2 = s, and of concentration c and 1 - c, correspondingly. We apply field-theoretical renormalization group approach to analyse the renormalization group flow for different initial conditions, triggered by s and c, and to calculate effective critical exponents further away from the fixed points of the renormalization group transformation. We show how the effective exponents are governed by difference in properties of the magnetic components.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Techno-economic comparison of sCO 2 cycles for particle-based CSP at design-point conditions

In this work, we compare the techno-economic performance of supercritical carbon dioxide power cycles integrated in a particle CSP system. We model four core cycle configurations: simple (with optional bypass), recompression (with optional bypass), partial cooling, and turbine split flow, which each demonstrate different benefits in a CSP system, such as high efficiency, low cost, or large HTF temperature differences. We parametrically sweep cycle design variables for each configuration. The set of power cycle performance results are then combined with a design point particle CSP system model which calculates the system specific cost. The simple cycle and turbine split flow cycles have the best performance in the baseline results, with system specific costs of 5,912 and 5,899 $\$$/kWe respectively. In addition to the baseline set of results, we also vary key parameters and costs in a sensitivity study. The cycle designs with the best system performance limit their efficiency to ~45 %, despite demonstrating higher maximum efficiencies, due to the rapid increase in cost of recuperation as efficiencies rise. The simple cycle has strong performance in the analysis and is on average only 1.4 % worse than the optimal configuration. Lowered turbine inlet temperatures from the sensitivity study improve performance by reducing the PHX and turbine cost. Decreasing the inlet temperature from 700 to 625°C results in an >8 % decrease in system specific cost. Future work should expand sensitivity analyses to colder turbine inlet temperatures and calculate system performance by simulating annual performance with off-design solar and cycle component models.

14 SOLAR ENERGY↗
DOE OSTI2023

Zirconium nanoparticle coating development for FCCI diffusion barrier in nuclear cladding

Fuel clad chemical interaction (FCCI) is a pervasive issue for metallic nuclear fuels. FCCI can damage cladding, reduce thermal efficiency, and compromise the structural integrity of the fuel rod. In this study, we evaluated the FCCI mitigation capabilities of electrophoretically-deposited Zr nanoparticle coatings on ferritic-martensitic stainless-steel cladding (HT9). Zr was deposited on the internal surface of 6 mm outer diameter cladding using electrophoretic deposition (EPD). Thickness varied among samples with similar EPD conditions. A uniformly coated sample was selected and sent to the Argonne Tandem Linac Accelerator System (ATLAS) to simulate the coating's behavior under irradiation conditions. Furthermore, the coating prevented FCCI in diffusion couple studies. Results of this study confirm that Zr could potentially make an effective FCCI barrier in metallic fuel rods but the deposition process needs improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

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

A High-Efficiency Delayed Update Algorithm for Evaluating Slater Determinants in Quantum Monte Carlo

For quantum Monte Carlo simulations of molecular systems or supercells with thousands of electrons, matrix operations related to Slater determinants lead the computational cost. McDaniel et al. [J. Chem. Phys. 2017, 147, 174107] proposed a delayed update algorithm to increase computational efficiency by using matrix–matrix multiplication when updating the inverse matrices of Slater determinants. However, preparing intermediate matrices for applying the Sherman–Morrison–Woodbury formula remained a bottleneck. Here, in this work, we introduce an improved algorithm for CPUs and GPUs that (1) reduces this bottleneck by iteratively updating the intermediate matrices and (2) is efficient at any acceptance ratio, with no cost for rejected moves on CPUs and minimal cost on GPUs. We show the full scheme of integrating the delayed update algorithm into a single-electron move. The high efficiency of our algorithm is demonstrated on CPUs and GPUs for a 512 atom/6144 valence electron calculation, with 12× and 2× overall speed-up compared to traditional rank-1 update schemes in diffusion quantum Monte Carlo, respectively.

Luo, Ye [Argonne National Laboratory (ANL), Argonn↗
DOE OSTI2023

Dynamics of Aqueous Electrolyte Solutions: Challenges for Simulations

This Perspective article focuses on recent simulation work on the dynamics of aqueous electrolytes. It is well-established that full-charge, nonpolarizable models for water and ions generally predict solution dynamics that are too slow in comparison to experiments. Models with reduced (scaled) charges do better for solution diffusivities and viscosities but encounter issues describing other dynamic phenomena such as nucleation rates of crystals from solution. Polarizable models show promise, especially when appropriately parametrized, but may still miss important physical effects such as charge transfer. First-principles calculations are starting to emerge for these properties that are in principle able to capture polarization, charge transfer, and chemical transformations in solution. Finally, while direct ab initio simulations are still too slow for simulations of large systems over long time scales, machine-learning models trained on appropriate first-principles data show significant promise for accurate and transferable modeling of electrolyte solution dynamics.

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

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗