Search NASASearch

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 145 records · Page 8

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)

Scheduler Modeling of Distributed Energy Resources for Providing Ancillary Services

Distribution energy resources (DERs) have been integral components of modern power systems, and their capability to provide grid services has been widely studied. To promote the deployment of these resources in providing grid services in real-world utility operations, this paper proposes a day-ahead scheduler model for a distribution system connected DER plant. A certain amount of generation capacity of this DER plant is reserved for frequency services, and some ancillary services for the distribution system-including peak load reduction, voltage regulation, and power factor control-are integrated into the model. The model is tested on a real-world distribution system. From the simulation results, the energy and reserve schedule of the solar photovoltaic (PV) unit and battery energy storage system (BESS) can be determined, and voltage and power factor are well maintained. Additionally, in order to demonstrate the specific characteristics of the co-located and hybrid operation modes for the PV and BESS, these two modes are analyzed both theoretically and through real-time simulation. Simulation results show that most of PV's variability is transferred to the net power in the co-located mode, whereas it is transferred to the BESS in the hybrid mode. This proposed scheduler model and the comparison of co-located and hybrid modes can provide practical guidance for the applications of DER plant in the real-world utility.

14 SOLAR ENERGY

Cross Correlating Cosmological Probes for Stage-IV Surveys

The next few years will be populated with Stage-IV cosmological surveys. This paper forecasts the cosmological information we expect to constrain using probes from the Rubin Observatory Legacy Survey of Space and Time (LSST) and CMB Stage-4 surveys. We explore the effectiveness of different large-scale-structure probes, including the power spectra of galaxy weak lensing, galaxy positions, Cosmic Microwave Background (CMB) lensing, and their cross-correlation. We use correlated lognormal simulations with the expected redshift distribution, number density of galaxies, and noise levels of the LSST survey, as well as CMB-S4 lensing simulations with anticipated noise levels, to obtain the expected error bars of each of the probes. We investigate the constraining power of the cosmological parameters for the individual probes and their combinations in an idealized scenario, where we do not account for astrophysical and observational parameters such as galaxy bias variations and photometric redshift uncertainties.

Gibbins, Grace

The Effects of Turbulence on Heliosheath Ions and Implications for Energetic Neutral Atoms

The distribution of ions in the heliosheath—the region between the heliospheric termination shock and the heliopause—is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distributions were estimated previously based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but these estimates only provide the distributions near the shock. In this work, we use self-consistent hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase in the temperature of thermal solar wind ions. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.

79 ASTRONOMY AND ASTROPHYSICS

Controllable Formation of Threefold-Coordinated Oxygen in Graphene by Low-Energy Ion Implantation

The atomically precise engineering of impurities in graphene and the understanding of their structural and carrier-dependent electronic properties at the nanoscale are crucial for advancing graphene-based nanoelectronics, catalysis, and energy technologies. Here, we demonstrate controllable incorporation of the elusive 3-fold-coordinated O substitutions into graphene using low-energy O + ion implantation under ultrahigh-vacuum conditions. By combining high-resolution scanning tunneling microscopy and spectroscopy (STM/S), bond-resolved noncontact atomic force microscopy techniques, and density functional theory (DFT) calculations, we resolve both the structural and electronic properties of the O-related defects. The STM/S measurements, corroborated by DFT calculations, uncover a characteristic impurity state that is energetically pinned to the Dirac point across different charge-carrier doping regimes. Molecular dynamics simulations further reveal the distribution of implantation-induced configurations and identify the formation of 3-fold-coordinated O dopants. Furthermore, this work provides a viable route to incorporate 3-fold-coordinated O dopants and opens new opportunities for controlled defect engineering in graphene.

3-fold-coordinated oxygen

Mapping the structural–mechanical landscape of amorphous carbon with ReaxFF molecular dynamics

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropy and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.

Dernov, A. (ORCID:0009000004220973)

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

IEEE 123-Bus System A Matrix

System component matrix for IEEE 123 bus distribution system for dynamic simulations.

Sahu, Vibhuti [Oak Ridge National Laboratory] (ORC

Cross Correlating Cosmological Probes for LSST & CMB-S4

The upcoming years will be populated with state-of-the-art Stage-IV cosmological surveys. This paper forecasts the cosmological information we expect to constrain using probes from the Rubin Observatory Legacy Survey of Space and Time (LSST) and CMB-S4. We explore the effectiveness of different large-scale-structure probes, including the power spectra of galaxy weak lensing, galaxy position, Cosmic Microwave Background (CMB) lensing, and their cross-correlations. We use correlated lognormal simulations with the expected redshift distributions, galaxy number densities, and noise levels of the LSST survey. For CMB weak lensing, we use CMB-S4 lensing simulations with anticipated noise levels to obtain the expected error bars of each probe. We investigate the constraining power of the cosmological parameters for the individual probes and their combinations in an idealized scenario. We do not take into account astrophysical and observational parameters such as galaxy bias variations and photometric redshift un- certainties. Overall, we find the auto-correlated probes hold stronger constraints on cosmological parameters than the cross-correlated probes due to kernel differences. However, the cross-correlation of galaxy clustering and CMB lensing is very comparable to the auto-correlation of galaxy clustering with only an 8% stronger constraint in S8. It will be important for future surveys to use these auto and cross probes in combination due to the different astrophysical and systematic effects involved in CMB and late-time galaxy data. When systematics are added to this analysis, the cross probes will serve as a check for systematic biases in individual probes. In addition, we find a significant increase in constraint using all four probes in combination. This can be illustrated by the increase in constraint of S8 by 93.2% comparing galaxy clustering auto-correlation to a combination of all probes. Our study is a first step toward forecasting the high-precision cosmological constraints we expect to obtain using the next generation of large-scale structure probes.

Gibbins, Grace

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Impact of atomic substitution on core-hole relaxation dynamics: A study of Br 2 and IBr

Understanding inner-shell decay processes in heavy-element molecules is essential for unraveling x-ray-induced photodynamics and advancing molecular imaging techniques. Here, in this study, we investigate the influence of atomic substitution on core-hole relaxation dynamics and molecular fragmentation in Br 2 and IBr, initiated by x-ray ionization absorption at the Br K-edge. Using a combination of x-ray/ion coincidence measurements and Monte Carlo/molecular dynamics simulations, we track the charge distribution and the kinetic energy release (KER) of fragment ions with a total charge from 2+ to 8+. For both molecules, the simulated KER values show good agreement with experiment across different fragmentation channels. Our comparison reveals that substituting Br with the heavier I atom in IBr has a minimal impact on the inner-shell electronic decay process but significantly influences nuclear motion, leading to slower dissociation and thus a KER close to the Coulomb limit—an effect attributed to the atomic mass. These findings highlight the interplay between electronic and nuclear effects in molecular fragmentation, particularly in heavy-element species, and provide new insights into medical therapies, structural biology, and astrophysics.

Bhat, Nivedita [Argonne National Laboratory (ANL),

Ray-tracing image simulations of transparent objects with complex shape and inhomogeneous refractive index

Optical images of transparent three-dimensional objects can be different from a replica of the object’s cross section in the image plane, due to refraction at the surface or in the body of the object. Simulations of the object’s image are thus needed for the visualization and validation of physical models. We report ray tracing image simulations that achieved high physical fidelity, reproducing optical behaviors and image features not rendered in previous studies. We replicated brightfield microscopy images of drops with complex shapes, and images of pressure and shock waves traveling inside them. For high physical fidelity, the simulations must replicate the spatial and angular distribution of illumination rays, and both the experiment and the simulation must be designed for accurate optical modeling. The simulations are highly sensitive to the properties of the drops and can be used to diagnose and refine fluid dynamics models. The simulated images can also be optimized to extract multiple 3D properties from experimental images. Compared to specialized single-shot 3D imaging methods, this approach has the advantage that it preserves the experimental simplicity, the high resolution, and the visual interpretability characteristic to basic optical imaging. The techniques introduced here are directly applicable to optical microscopy, so they can be used in other fields, such as microfluidics and biology, to expand the type and the accuracy of three-dimensional information that can be extracted from basic optical images.

Cavitation

Equivalent Point Sources for Simulating Nuclear/Radiological Search Operations

Modeling and simulation can be used to estimate the detection of different types and strengths of radioactive sources. Simulations of shielded gamma-ray sources in complex environments, such as searches in urban areas using vehicle-mounted detectors, can demonstrate what source strengths and shielding combinations should be able to be detected. For faster calculations, gamma-emitting sources that could have different shielding layers can be simulated ahead of time, resulting in equivalent point source emission distributions that can be used in larger simulation models. Larger-scale models will then run faster using equivalent point sources because the shielding layers around the sources have already been accounted for. Accompanying this paper are 260 equivalent point source emission files for 72 sources with different shielding layers, formatted for three different Monte Carlo codes. Users wanting to request versions of the SCALE and/or MCNP software can find details here: https://rsicc.ornl.gov/

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Enabling Grid-Forming Control with Fault Ride-Through in Unbalanced Distribution Networks

Distribution networks are often unbalanced, causing oscillatory responses in inverter control designed for balanced conditions. Here, to address this problem, this paper proposes a novel time-domain transformation appropriate for inverter control and enables the decomposition of three-phase unbalanced signals into constant positive and negative components. Relations useful for calculating unbalanced active and reactive power are derived from first principle, providing insight into vector products of unbalanced three-phase signals. Furthermore, a grid-forming control effective under unbalanced conditions is developed, which delivers superior performance while meeting UNIFI1 specifications for grid-forming control under unbalanced conditions. specifications applicable to category 4 inverter-based resource, like setting and regulating frequency/voltage, providing voltage support, sharing active power, injecting negative sequence current, and riding through faults. A current limiter is proposed for safe fault ride-through and integrates with the grid-forming control featuring frequency/voltage droop controllers and current and voltage control loops. The transformation of interconnected inverters is formulated and stability of the proposed control analyzed to support robust parameter selections. The effectiveness of the proposed transformation and grid-forming control is demonstrated through analytical results and real-time simulation of a IEEE 123 distribution network on the Real-Time Digital Simulator. Comparison with existing methods shows that the proposed strategy satisfies the UNIFI specifications with a much better performance.

24 - POWER TRANSMISSION AND DISTRIBUTION

CHARGE-MAP: An integrated framework to study the multicriteria EV charging infrastructure expansion problem

The widespread adoption of electric vehicles (EVs) in recent years has necessitated the development of effective charging infrastructures. However, charging infrastructure expansion is a multifaceted problem that requires careful consideration of the existing infrastructure, spatiotemporal distribution of charging demands, power-grid capacity, and budget constraints. Here, to approach this complex problem, we present CHARGE-MAP, a data-driven simulation-optimization framework, focused on ensuring meaningful charging experience for individual EV owners. CHARGE-MAP integrates three modules: an agent-based simulation module that estimates spatiotemporal distribution of charging demands by modeling EV adopter mobility and charging behavior; an optimization module that determines optimal new charging station/charger locations and capacities, while minimizing expected detour distances and wait-times with a limited number of new stations; and a power module that determines how to connect the stations to the power grid while maintaining its stability. Using the state of Virginia (consisting of 95 counties and 38 independent cities) as a case study, our results show that CHARGE-MAP can meet the demand of ~198,600 predicted EVs with 1,305 new public charging stations and 2,164 new chargers. It reduces average detour distances for charging by 66% and wait-times at stations by 72% compared to the existing infrastructure. Furthermore, transformer capacity requirement analysis reveals that only 1.8% of residential transformers require upgrades, while over 80% of commercial charging locations can be supported with modest transformer infrastructure (25 to 50 kVA). This indicates that targeted investments can facilitate cost-effective EV integration. Consequently, CHARGE-MAP provides policymakers and urban planners with crucial data-driven insights for effective EV charging infrastructure expansion. Sign up for PNAS alerts.

charging infrastructure

Impact of position and density of nanoscale voids on fracture initiation in iron from phase field fracture simulation

Understanding the impact of these bubbles on crack propagation, like that of helium bubble-induced cracking in irradiated materials is incredibly complex. A useful first study towards understanding bubble effects on fracture is to examine how voids impact fracture first. In this work, we used phase-field fracture simulations to examine the influence of voids and their distribution on Mode I fracture in Fe. Assuming brittle fracture, two simulation configurations were considered: (1) nanoscale systems with one or two voids, and (2) nanoscale systems with an experimentally relevant distribution of voids, with up to 20% void area. Results from simulations with one and two voids showed that voids within 10 nm of a crack tip reduce the stress required for crack growth, with the magnitude of reduction depending on void-to-crack orientation. Comparisons with linear elastic fracture mechanics and evaluation of one versus two void systems revealed deviations from linear superposition, implying complex interactions between void and crack tip stress fields. In multi-void simulations, as void sizes increase, the nearest void to the crack tip exerts a greater influence on fracture stress than the overall porosity. Furthermore this study provides valuable insights into the relationship between void size and concentration, and the stress necessary for crack growth, marking a step forward towards understanding He bubble-induced fracture in ferrous materials.

36 MATERIALS SCIENCE