SEARCH · Engineering Papers

Results for “complex simulations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Dynamic clay microstructures emerge via ion complexation waves

Clays control carbon, water and nutrient transport in the lithosphere, promote cloud formation5 and lubricate fault slip through interactions among hydrated mineral interfaces. Clay mineral properties are difficult to model because their structures are disordered, curved and dynamic. Consequently, interactions at the clay mineral-aqueous interface have been approximated using electric double layer models based on single crystals of mica and atomistic simulations. We discover that waves of complexation dipoles at dynamically curving interfaces create an emergent long-range force that drives exfoliation and restacking over time- and length-scales that are not captured in existing models. Curvature delocalizes electrostatic interactions in ways that fundamentally differ from planar surfaces, altering the ratio of ions bound to the convex and concave sides of a layer. Multiple-scattering reconstruction of low-dose energy-filtered cryo electron tomography enabled direct imaging of ion complexes and electrolyte distributions at hydrated and curved mineral interfaces with {\aa}ngstrom resolution over micron length scales. Layers exfoliate and restack abruptly and repeatedly over timescales that depend strongly on the counterion identity, demonstrating that the strong coupling between elastic, electrostatic and hydration forces in clays promote collective reorganization previously thought to be a feature only of active matter.

Whittaker, Michael L↗

DOE OSTI2026

Design of Controller Hardware-In-the-Loop Model of Microgrid with Modular Building Blocks and Automated Design Script

The scalability of controller hardware-in-the-loop (CHIL) simulation is critical for validating control coordination and energy management in microgrids with distributed energy resources, especially as these modern systems become more complex and decentralized. This paper presents a CHIL modeling methodology that combines modular building blocks with an automated design script to streamline the development of high-fidelity microgrid models. Standardized subsystem templates for resources, converters, and buses are integrated with a Python-based script that compiles structured JSON configuration files into simulation-ready initialization code. The proposed approach reduces development time, improves model consistency, and enhances simulation fidelity. The methodology is validated on a Typhoon HIL604 platform and is broadly applicable to real-time simulation of complex, networked microgrid systems. This framework establishes a foundation for automated, scalable CHIL validation and accelerates the design of next-generation distributed energy systems.

Kim, Namwon [ORNL] (ORCID:0000000200438489)↗

DOE OSTI2022

Radiation Belt Remediation: A Complex Engineered System (RBR-ACES) (Institutional Computing Report)

We have performed the simulation campaigns in support of LDRD-DR project Radiation Belt Remediation: A Complex Engineered System (RBR-ACES). This project involves modeling wave generation and wave propagation in plasma environment and studying the effect of waves on the distribution of high energy electrons with the focus on electrons appeared after a high-altitude nuclear explosion (HANE). One of the major goal of this research is strong validation of the numerical models with laboratory experiments and space-based experiments and observations.

42 ENGINEERING↗

DOE OSTI2022

Hybrid Analysis of Fusion Data for Online Understanding of Complex Science on Extreme Scale Computers

The current practice for fusion scientists running first principle simulations on high performance computing plat-forms is to either run their simulations and output their data for post-hoc analysis, or to place in situ analytics into their code. In this paper we examine a complex workflow using XGC fusions simulation run on the Oak Ridge Leadership Computing Facility's supercomputer Summit, which also involve three anal-yses as part of the results necessary for scientific discovery. We discuss the challenges faced when implementing these algorithms and present an original hybrid staging technique to help enable the physicists to make discoveries during the execution of the simulation. By creating this infrastructure, we can examine complicated physics results, which may not have been possible without the infrastructure. For example, our work enables the online visualization of turbulent homoclinic tangle around the magnetic X-point, breaking the last confinement surface. This visualization could help fusion scientists to better understand and improve the turbulence spread of plasma exhaust heat, which is crucial toward realizing plasmas beyond the currently accessible physics regimes of present-day tokamak reactors. The physics of turbulent homoclinic tangle will be reported in a future physics publication, by utilizing the original online analysis/visualization framework presented in this paper.

Suchyta, Eric↗

DOE OSTI2025

Benchmarking Near-Surface Winds in the HRRR Analyses Using Multisource Observations over Complex Terrain in the Southeastern United States

Wind energy plays a crucial role in sustainable power generation, yet its full potential in the southeast United States (SEUS) remains underexplored. This study advances wind resource assessment in the SEUS using existing regional modeling and multisource observations. We find that terrain complexity, influenced by orography and forest canopy, significantly impacts the accuracy of modeled wind speed. While the High-Resolution Rapid Refresh (HRRR) model effectively simulates wind profiles over flat, nonforested terrain, larger errors are produced above forest canopies, and particularly for those stands in hilly and mountainous terrain where there are no observations aloft that could be used by HRRR’s data assimilation software. Seasonal variations, e.g., changes in the leaf area index, further complicate the relationship between terrain complexity and simulation errors. Here, our findings emphasize the critical need for comprehensive wind profile measurements, extending from below the canopy to height above the canopy and through the lower planetary boundary layer, to fully quantify model performance over and near forests, and for improving wind resource assessment, planning, and management. This study provides valuable insights not only for further model development but also for future field campaign deployment with the goal of further improving our understanding of wind resources in the region.

Complex terrain↗

DOE OSTI2025

Quantum magic and multipartite entanglement in the structure of nuclei

Motivated by the Gottesman-Knill theorem, we present a detailed study of the quantum complexity of p -shell and s d -shell nuclei. Valence-space nuclear shell-model wave functions generated by the code are mapped to qubit registers using the Jordan-Wigner mapping (12 qubits for the p shell and 24 qubits for the s d shell), from which measures of the many-body entanglement ( n -tangles) and magic (nonstabilizerness) are determined. While exact evaluations of these measures are possible for nuclei with a modest number of active nucleons, Monte Carlo simulations are required for the more complex nuclei. The broadly applicable Pauli-string I ̂ Z ̂ exact (PSIZe) Markov chain Monte Carlo (MCMC) technique is introduced to accelerate the evaluation of measures of magic in deformed nuclei (with hierarchical wave functions), by factors of ≈ 8 for some nuclei. Significant multinucleon entanglement is found in the s d shell, dominated by proton-neutron configurations, along with significant measures of magic. This is evident not only for the deformed states, but also for nuclei on the path to instability via regions of shape coexistence and level inversion. These results indicate that quantum-computing resources will accelerate precision simulations of such nuclei and beyond. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Evaporation Induced Entrainment of Contaminants from Evaporating and Burning Liquids.

This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A parametric SPICE model for the simulation of spark gap switches

The use of SPICE-based software for the simulation of pulsed power systems—even large complex systems—has become commonplace in the pulsed power community. This is in contrast to earlier work in the field that relied on specially developed simulation codes such as Sandia’s Screamer or the Navy Research Lab’s Bertha, which natively incorporated models for common pulsed power components such as spark gap switches. Unlike these programs, SPICE programs provide a simple and familiar user interface and wide availability. However, SPICE programs do not include realistic models for key pulsed power circuit devices—including the spark gap switch. While simple switch models do exist in SPICE programs, these can only crudely approximate the behavior of a spark gap. This effort focuses on developing an SPICE circuit model for a gas-filled spark gap switch that is physically realistic while being simple enough to permit simulations to run in reasonable times on typical personal computers. Detailed information is provided for implementation in two common versions of SPICE: LTspice and Orcad PSPICE. Adaptation to other SPICE programs is possible with minimal modification. The model is intended as a design tool that uses physical parameters as inputs to connect it directly to the development of useable pulsed power systems. Data collected from the operation of a high-pressure pulsed-charged switch and a complete 12-stage Marx generator have been used to demonstrate the implementation and accuracy of the model over a wide range of parameters.

Pouncey, J. Cameron (ORCID:0000000153500816)↗

DOE OSTI2024

Multi-Federate Co-Convergence with HELICS

In co-simulation studies, convergence refers to the ability of different simulation tools involved to achieve a consistent and stable solution. Convergence is a critical aspect of co-simulation as it determines the accuracy of the results obtained from the simulation. Convergence in co-simulation studies depends on several factors, this includes system complexity, the accuracy of the models used, and the numerical methods employed by the simulation tools. It is essential to ensure that the coupling interfaces between the different tools are designed to allow for data exchange in a consistent and accurate manner. Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS) is an open-source co-simulation framework developed for the energy domain. This paper explores the convergence performance of a set of co-simulation use cases. We further explore the use of a co-convergence helper federate to help with co-simulation convergence. The convergence efficacy of several algorithms (both gradient-based and gradient-free) is tested against these use cases. Finally, the sensitivity of these algorithms to several factors, such as system scaling and others, is tested and detailed in this paper. Our results show that for a subset of use cases, the co-convergence helper federate is able to improve co-simulation convergence significantly.

co-convergence↗

DOE OSTI2025

Data driven methods to recognize patterns in EIC weak-strong simulation

Beam-Beam simulations are currently being studied in preparation for future EIC experiments to study beam-beam effects and, in turn, maximize luminosity. Weak-strong methods are studied for single-particle dynamics during collision. 1 million macro-particles for 1 million turns are typically tracked, corresponding to only 10 seconds in the EIC. The goal of this study is to predict beam properties over the scale of hours. A potential solution focuses on using data-driven methods such as machine learning methods to analyze and extend the insights of the beam properties such as long-term nonlinear effects. This would aid in long-term predictions where results would be more efficiently acquired than a typical tracking simulation. Some limitations such as inaccurate predictions and spatial complexity are also discussed. These methods can then be applied to strong-strong simulations in the future studies.

Accelerator Physics↗

DOE OSTI2022

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Finite-Difference Simulation for Infrasound Generated by Finite-Extent Ground Motions

Underground explosions can produce infrasound in the atmosphere, and the wavefield characteristics are often governed by the ground surface motions. Finite-difference methods are popular for infrasound simulation as their generality and robustness allow for complex atmospheric structures and surface topography. A simple point-source approximation is often used because infrasound wavelengths tend to be large relative to the source dimensions. However, this assumption may not be able to capture the complexity of explosion-induced ground motions if the surface area is not compact, and appropriate source models must be incorporated into the finite-difference simulations for accurate infrasound prediction. In this study, we develop a point source representation of the complex ground motions for infrasound sources. Instead of a single point source, we use a series of point sources distributed over the source area. These distributed point sources can be equivalent to air volume changes produced by the ground motions in the atmosphere. Further, we apply the distributed point-source method to a series of buried chemical explosions conducted during the Source Physics Experiment Phase I. Epicentral ground-motion measurements during the experiments provide a way to calculate accurate distributed point sources. We validate and evaluate the accuracy of distributed point source approach for infrasound simulations by direct comparison with acoustic observations in the field experiment.

58 GEOSCIENCES↗

DOE OSTI2021

Tidal Stream Energy Resource Characterization in the Salish Sea

Tidal stream energy holds great potential as a renewable energy source in regions of the world where tidal currents are strong and predictable. The tidal stream energy resource of a region is strongly controlled by its tidal wave characteristics, the local bathymetry, and coastal geometric features. The Salish Sea, a large estuary on the Pacific Northwest coast, represents a great tidal stream resource because of its strong tidal currents in many tidal channels. However, the tidal resource in the Salish Sea has not been systematically characterized, primarily because of its large area and complex bathymetry and coastlines. This paper presents a modeling study conducted to characterize the tidal energy resource of the Salish Sea based on a high-resolution three-dimensional tidal hydrodynamic model of the Salish Sea, which was extensively validated using data derived from 10 tidal gauge stations and 132 Acoustic Doppler Current Profiler stations. Model validation results indicated the Salish Sea hydrodynamic model is skillful in simulating tidal wave propagation and velocity distributions in the Salish Sea. Based on model results, a total of 16 tidal channels with strong tidal currents, 9 in the San Juan Islands and Rosario Strait regions and 7 in the Puget Sound, were identified as hotspots for potential tidal energy development. Velocity probability distributions and exceedance curves of cross-channel average velocity magnitudes were calculated at all 16 channels based on the recommendations of the International Electrotechnical Commission Technical Specifications for tidal energy resource characterization. The tidal energy resource at the 16 hotspots was also characterized using power density distributions and kinetic energy fluxes. The ranking of the kinetic energy fluxes suggested that Admiralty Inlet, Rosario Strait, and Middle Channel are the top three hotspot sites in the Salish Sea for tidal stream energy development. The study demonstrated the need for a high-resolution 3-D modeling framework for accurate simulation of tidal currents in large complex estuarine systems, in the context of tidal energy resource characterization.

tidal stream energy, resource characterization, Sa↗

DOE OSTI2025

Asynchronous-many-task systems: Challenges and opportunities - Scaling an AMR astrophysics code on exascale machines using Kokkos and HPX

Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗