Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Responses of summer mesoscale convective systems to irrigation over the North China Plain based on convection-permitting model simulations

Extensive irrigation activities in the North China Plain (NCP) significantly influence regional weather and climate. However, previous studies focusing on the NCP were primarily based on coarse-resolution models, which are unable to explicitly resolve convection systems, causing large uncertainty in precipitation simulations. In this study, a convection-permitting model coupled with a dynamic irrigation scheme is utilized to investigate the impacts of irrigation on summertime mesoscale convective systems (MCSs) over the NCP. Sensitivity experiments with irrigation off and on are conducted for 5 summers and an MCS identification and tracking algorithm is applied to both satellite observations and model simulations. We find that incorporating irrigation in the model increases MCS precipitation, which agrees more with observations. The probability distributions of MCS lifetime, area, propagation speed, and intensity are all better simulated with irrigation. Irrigation increases the occurrence frequency of MCSs throughout the entire day. The nighttime increase is partly because of more frequent local initiation of MCS developed from isolated deep convection, while the daytime increase is mainly attributed to the changes in MCSs initiating elsewhere and then propagating to the NCP. On average, irrigation induces additional moisture that is more thermodynamically favorable for precipitation, but this effect is partially offset by the weakened ascending air motion primarily caused by irrigation surface cooling. Compared to weak MCS precipitation events, strong MCS precipitation events experience greater enhancement in precipitation intensity when including irrigation because the offset effect from the change in large-scale ascending air motion is insignificant. In addition, irrigation makes the variation of MCS precipitation intensity more correlated with the variation in ascending motion but less correlated with that in atmospheric moisture. Our results suggest the pronounced impacts of irrigation on MCSs over the NCP which should be included in numerical models to improve regional precipitation simulation and prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predicting the Evolution of Shallow Cumulus Clouds With a Lotka‐Volterra Like Model

Abstract In numerical weather prediction and climate models, boundary‐layer clouds are controlled by a wide range of subgrid‐scale processes. However, understanding the nature of these processes and their role in the evolution of the cloud size distribution as a whole has been elusive. To address this issue, we adopt a novel empirical framework from the field of population dynamics to model the evolution of cloud size statistics by using the shallow cumulus properties obtained from a large‐eddy simulation (LES). Our approach involves representing the cloud size distribution and the total cloud area using a revised Lotka‐Volterra model and ridge linear model, respectively. The physical interpretation of the total cloud area and coefficients obtained from the optimization of the models reveals three stages probably interpreted by dominant processes: the formation of new clouds, the growth of single clouds, and a steady state with organized transitions involving the growth and decay of multiple clouds. Furthermore, we showcase the potential of this framework to serve as a component of scale‐aware parameterizations of shallow‐convective clouds in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrogen Bonding in Amorphous Indomethacin

Amorphous Indomethacin has enhanced bioavailability over its crystalline forms, yet amorphous forms can still possess a wide variety of structures. Here, Empirical Potential Structure Refinement (EPSR) has been used to provide accurate molecular models on the structure of five different amorphous Indomethacin samples, that are consistent with their high-energy X-ray diffraction patterns. It is found that the majority of molecules in amorphous Indomethacin are non-bonded or bonded to one neighboring molecule via a single hydrogen bond, in contrast to the doubly bonded dimers found in the crystalline state. The EPSR models further indicate a substantial variation in hydrogen bonding between different amorphous forms, leading to a diversity of chain structures not found in any known crystal structures. The majority of hydrogen bonds are associated with the carboxylic acid group, although a significant number of amide hydrogen bonding interactions are also found in the models. Evidence of some dipole–dipole interactions are also observed in the more structurally ordered models. The results are consistent with a distribution of Z-isomer intramolecular type conformations in the more disordered structures, that distort when stronger intermolecular hydrogen bonding occurs. The findings are supported by 1H and 2H NMR studies of the hydrogen bond dynamics in amorphous Indomethacin.

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

A Power-Hardware-in-the-Loop (PHIL) Evaluation of Service Restoration With Networked Microgrids

This paper describes the power-hardware-in-the- loop (PHIL) evaluation of the feasibility of service restoration solutions determined by the PowerModelsONM.jl tool. This tool incorporates microgrids and the networking of microgrids into its determination of an optimal service restoration solution. The paper presents PHIL simulation results for a case study based on a real distribution feeder with multiple microgrids, showcas- ing the effectiveness of networked microgrids in aiding system restoration after an outage. The study leverages high-fidelity, real- time electromagnetic transient models to ensure accuracy in the simulation results. This work is the final output from the Resilient Operation of Networked Microgrids (RONM) project funded by the U.S. Department of Energy Office of Electricity Microgrid Program and led by Los Alamos National Laboratory. RONM focused on the application of PowerModelsONM.jl to enhance the resilience of distribution systems.

fault location isolation and service restoration (↗
DOE OSTI2024

Development of a Distribution Optimal Power Flow Federate for Open-Source OEDI-SI Platform

Increasing numbers of distributed generators in the electric power distribution networks require developing a control strategy to optimize solutions in real time. Linearized optimal distribution flow development has seen growth and acceptance in the distribution systems literature for efficiently modeling the \glspl{opf} for distribution systems. This paper examines the implementation and integration procedure for linearized optimal distribution flow federate to \gls{oedisi} platform. Specifically, we discuss i) the usage of the \gls{oedisi} platform, ii) obtaining a tractable solution using developed \gls{opf} federate, and iii) validation of solutions and bench-marking the \gls{oedisi} platform with developed \gls{opf} federate using OpenDSS. In brief, we demonstrate how a general linearized optimal distribution flow federate can be developed and integrated with a co-simulation environment to mimic real-world examples. The efficacy of the proposed method is demonstrated using the IEEE 123-bus test system under different scenarios to obtain a tractable solution and compare its results.

Sadnan, Rabayet↗
DOE OSTI2025

Real-Time Lifetime Prediction of Semiconductor Devices Using Hardware-in-the-Loop

This paper presents a unique approach to enable real-time lifespan prediction of semiconductor power modules using a Hardware-in-the-Loop (HIL) system. By integrating the module's overall loss characteristics-specifically switching and conduction losses-with a thermoelectric model of the thermal management system, this research demonstrates that the model can dynamically estimates the junction temperature profile of the semiconductor devices in response to a changing torque demand profile for the motor drive system. This capability enables continuous monitoring of the module's operational time and cumulative stress induced on the devices to compute accumulated remaining lifetime or time-to-failure (TTF). This study provides an architectural framework for the HIL system with high-fidelity component models of multiple physical domains, allowing simulation of dynamic behaviors of a closely-coupled motor drive system. The advanced real-time computation and measurement functionalities of the HIL system allow for both dynamic lifetime calculations based on simulated data and aggregate lifetime predictions utilizing historical data. Moreover, this paper details an algorithm that not only computes cumulative damage but also synthesizes these data into a comprehensive aggregated lifetime metric. This methodology can enhance the maintenance scheduling strategies and operational reliability of semiconductor devices in critical applications, ultimately extending their service life while optimizing performance.

hardware-in-the-loop (HIL)↗
DOE OSTI2026

Particle Acceleration and Pitch-angle Evolution in Relativistic Turbulence

Synchrotron radiation detected from relativistic astrophysical objects such as pulsar-wind nebulae and jets from active galactic nuclei depends on the magnetic fields and the distribution functions of energetic electrons in these systems. Relativistic magnetically dominated turbulence has been recognized as an efficient mechanism for structure formation and nonthermal particle acceleration in these environments. Recent numerical simulations of relativistic turbulence have provided insights into the energy distribution functions of accelerated electrons. Much less is currently understood about their pitch angle distributions, which are crucial for accurately interpreting the spectra of synchrotron radiation. We perform a detailed case study of the pitch angle distributions formed during the process of turbulent acceleration for B 0 /δB 0 = 10 and $\tilde{σ}$ 0 ~40, where B0 is the uniform component of the magnetic field, δB0 is the fluctuating component, and $\tilde{σ}$ 0 is the plasma magnetization based on the magnetic fluctuations. We find that even minimal numerical noise can cause substantial pitch angle scattering, but we demonstrate techniques for overcoming the numerical challenges associated with the evolution of very small pitch angles. Our numerical results are consistent with the phenomenological considerations found in C. Vega et al. (2024a, 2025).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Tensor network simulations of quasi-GPDs in the massive Schwinger model

Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effectiveness of denoising diffusion probabilistic models for fast and high-fidelity whole-event simulation in high-energy heavy-ion experiments

Artificial intelligence (AI) generative models, such as generative adversarial networks (GANs), variational autoencoders, and normalizing flows, have been widely used and studied as efficient alternatives for traditional scientific simulations. However, they have several drawbacks, including training instability and inability to cover the entire data distribution, especially for regions where data are rare. This is particularly challenging for whole-event, full-detector simulations in high-energy heavy-ion experiments, such as sPHENIX at the Relativistic Heavy Ion Collider and Large Hadron Collider experiments, where thousands of particles are produced per event and interact with the detector. This work investigates the effectiveness of denoising diffusion probabilistic models (DDPMs) as an AI-based generative surrogate model for the sPHENIX experiment that includes the heavy-ion event generation and response of the entire calorimeter stack. DDPM performance in sPHENIX simulation data is compared with a popular rival, GANs. Results show that both DDPMs and GANs can reproduce the data distribution where the examples are abundant (low-to-medium calorimeter energies). Nonetheless, DDPMs significantly outperform GANs, especially in high-energy regions where data are rare. Additionally, DDPMs exhibit superior stability compared to GANs. The results are consistent between both central and peripheral centrality heavy-ion collision events. Moreover, DDPMs offer a substantial speedup of approximately a factor of 100 compared to the traditional Geant4 simulation method.

42 ENGINEERING↗
DOE OSTI2026

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Enhancing synchrotron radiation micro-CT images using deep learning: an application of Noise2Inverse on bone imaging

In bone-imaging research, in situ synchrotron radiation micro-computed tomography (SRµCT) mechanical tests are used to investigate the mechanical properties of bone in relation to its microstructure. Low-dose computed tomography (CT) is used to preserve bone's mechanical properties from radiation damage, though it increases noise. To reduce this noise, the self-supervised deep learning method Noise2Inverse was used on low-dose SRµCT images where segmentation using traditional thresholding techniques was not possible. Simulated-dose datasets were created by sampling projection data at full, one-half, one-third, one-fourth and one-sixth frequencies of an in situ SRµCT mechanical test. After convolutional neural networks were trained, Noise2Inverse performance on all dose simulations was assessed visually and by analyzing bone microstructural features. Visually, high image quality was recovered for each simulated dose. Lacunae volume, lacunae aspect ratio and mineralization distributions shifted slightly in full, one-half and one-third dose network results, but were distorted in one-fourth and one-sixth dose network results. Following this, new models were trained using a larger dataset to determine differences between full dose and one-third dose simulations. Significant changes were found for all parameters of bone microstructure, indicating that a separate validation scan may be necessary to apply this technique for microstructure quantification. Noise present during data acquisition from the testing setup was determined to be the primary source of concern for Noise2Inverse viability. While these limitations exist, incorporating dose calculations and optimal imaging parameters enables self-supervised deep learning methods such as Noise2Inverse to be integrated into existing experiments to decrease radiation dose.

Obata, Yoshihiro (ORCID:0000000303659129)↗
DOE OSTI2024

Glow Discharge Boronization and Real-Time Boronization Using an Impurity Powder Dropper in LHD

In the Large Helical Device (LHD), boronization is utilized for wall conditioning, employing two distinct methods: glow discharge using diborane gas and the injection of boron powder into the plasma via an impurity powder dropper (IPD). Glow discharge boronization is conducted after maintenance when the vacuum vessel is exposed to air and just before the start of an experimental campaign. In contrast, real-time boronization using the IPD is performed during discharges. The distribution of boron on plasma-facing components was assessed through long-term material probes in the case of glow discharge and through numerical simulations for the IPD method. Results indicated non-uniform boron distributions in both cases. The impact of these boronization techniques on impurity levels and the lifetime of their effects were also examined. Both methods effectively reduced oxygen and iron in the plasma, with the iron impurities originating from the first wall. For glow discharge boronization, the reduction in oxygen levels persisted until the end of the experimental campaign, whereas the reduction in iron levels only lasted a few days. In the case of boron injection during long pulse discharges exceeding 200 s, the suppression of oxygen and iron was observed. The effect on iron quickly diminished within several seconds after the injection ended, whereas the effect on oxygen lasted more than 100 s. The difference in the lifetimes of the boronization effects is further discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in dark sectors that are natural and theoretically motivated extensions to the standard model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the standard model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo's central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host satellite galaxies have cuspier density profiles and smaller stellar half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗
DOE OSTI2024

A performant energy-conserving particle reweighting method for Particle-in-Cell simulations

A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Emergence of Diverse Failure Patterns in Weathering‐Induced Landslides: Insights From Particle Finite Element Simulations

Weathering is a fundamental driver of landslide evolution over geological timescales. Despite its ubiquity and importance, quantifying how weathering drives the progressive destabilization of rock slopes remains challenging. In this work, we develop a unified computational framework based on the particle finite element method to investigate the evolution of weathering‐induced landslides, from long‐term weathering to short‐term slope failure and runout dynamics. The framework integrates key processes, including weathering front propagation, time‐dependent strength degradation, rupture surface development, and post‐failure runout dynamics. Through numerical simulation experiments, we elucidate how interactions among weathering characteristics (type, intensity, and rate law), bedrock strength, fracture distribution, and slope geometry govern the failure modes and kinematics of weathering‐induced landslides. Simulations show that matrix‐dominated weathering leads to shallow translational failures, whereas fracture‐dominated weathering produces deep‐seated rotational and compound landslides. Pre‐existing fractures and slope morphology also strongly influence the movement of destabilized landmasses, affecting the failure pattern (e.g., kinematic mode and rupture surface geometry) and post‐failure behavior (e.g., runout velocity). We further demonstrate that the failure time and volume of weathered slopes are governed by the competition between gravitational driving forces and cohesive resisting forces during progressive destabilization. These findings provide new insights into the fundamental mechanisms that drive the emergence of diverse failure patterns of weathering‐induced landslides with important implications for landslide hazard assessment.

Wang, Liang [Eidgenoessische Technische Hochschule↗
DOE OSTI2025

Tropical Cyclone Super Resolution using conditional diffusion denoising probabilistic model from mesoscale simulation to LES

Accurate modeling of tropical cyclone wind fields is essential for the design, risk assessment, and operational planning of offshore energy infrastructure. While mesoscale simulations are widely used thanks to their computational efficiency, they lack the necessary resolution to capture key features such as wind shear and veer profiles as well as the distribution turbulent kinetic energy (TKE). High-fidelity large-eddy simulation (LES) models on the other hand, can resolve turbulent structures and provide a more accurate representation of the complex wind field, albeit at a higher computational cost. To address this modeling gap, we introduce a two-part generative framework to enhance the resolution and physics-capturing ability of mesoscale simulations. First, a reduced-order model based on Karhunen–Loève (KL) decomposition is used to extract dominant spatial modes from one-dimensional mean wind profiles. A multilayer perceptron (MLP) is trained to map mesoscale mode weights to their LES counterparts, enabling accurate reconstruction of vertical velocity profiles. Second, a conditional Diffusion Denoising Probabilistic Model (DDPM) is developed to super-resolve coarse and low-fidelity mesoscale velocity fields, recovering fine-scale turbulence structures and stress distributions. The framework is evaluated across different tropical cyclone intensity categories defined by the Saffir–Simpson scale and demonstrates strong performance in both interpolation and extrapolation tasks. The generated fields accurately reproduce spatial coherence, stress distributions, and spectral energy characteristics observed in LES data. By bridging the fidelity gap between mesoscale and LES outputs, this approach offers a scalable, data-driven solution for enhancing the representation of tropical cyclone wind fields, enabling more robust offshore energy infrastructure systems design in tropical-cyclone-prone areas.

17 WIND ENERGY↗
DOE OSTI2026

The Role of Magnetic Reconnection in Energizing Protons and Heavier Ions at the Heliospheric Current Sheet

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions, indicating possible energization by magnetic reconnection up to 10 s–100 s keV nucleon −1 . Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff E max consistent with in situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as E max ∝ (Q/M) α with α ∼ 0.8–1.1, approximately consistent with PSP measurements in the broader range α ∼ 0.6–1.7. In the limit where ions are injected at the same energy per nucleon, α can be as low as ∼0.3. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

Murtas, Giulia [West Virginia Univ., Morgantown, W↗