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 433 records · Page 24

Real-Time Testbed for Studying Cyberattacks and Defense in DER-integrated Smart Inverter Systems

In this paper, we propose a Hardware-in-the-Loop (HIL) simulation testbed suitable for the implementation and testing of realistic cyberattacks on grid-tied smart inverter systems integrated with Distributed Energy Resources (DER) that use the Distributed Network Protocol-3 (DNP3) protocol for communications between grid components. Specifically, our testbed combines a Real-Time Digital Simulator (RTDS) NovaCor device, outfitted with GNETx2 network interface cards, a gridtied DER topology implemented via the RTDS software package RSCAD, and a custom virtual network that emulates a man in the middle attacker. The Man-in-the-Middle (MITM) attacker captures DNP3 traffic and falsifies telemetry data in DNP3 packets to trigger unwarranted commands from a DNP3 controller that exploit smart inverter grid support functions. We choose DNP3 and implement grid support functions according to the IEEE Std. 1547-2018 mandated for the interconnection and interoperability of DER power systems with associated power components. Furthermore, we develop a protocol payload agnostic attack detection framework that leverages the round-trip time (RTT) anomalies between DNP3 requests and responses and can detect the presence of attacks without having to analyze the payload’s contents, while balancing trade-offs between false alarm counts, missed detections, and time to detection. To facilitate further research, we publicly release benign and attack network traffic exchanged between various sensors, controllers, and actuators in our grid-tied inverter testbed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗

Modeling the impact of extreme weather events and future climate on the radiologically contaminated sites of Enewetak Atoll

Enewetak Atoll underwent 43 historical nuclear tests from 1948 to 1958, including the first hydrogen bomb test, resulting in a substantial nuclear material fallout contaminating the Atoll and the lagoon waters. The radionuclide fallout material deposited in lagoon sediments and land soil will remain for decades to come. With intensifying climate and extreme weather events, the possibility of redistribution of deposited radionuclide material has become a great concern. This study uses a numerical modeling approach to estimate the potential elevated radionuclide concentrations that can be distributed during storm events under current and future climates. We simulated three historical storm scenarios that are most likely to impact Atoll’s environment and remobilize the radionuclide-bound sediments. WRF-ARW was used to reconstruct these storm scenarios under current year (2015) and future year (2090) climates. Storm-induced ocean hydrodynamics conditions were generated using FVCOM. FVCOM-ICM was externally coupled to simulate the fate and transport of radionuclides. Given that the 239 Pu is the largest inventory of the lagoon and Atoll islands, the model results show the highest average 239 Pu concentration that an island may be exposed to is 3.25E-4 Bq/m 3 (becquerel per cubic meters), which is an increase of 84 times the average baseline/existing 239 Pu concentration without the storm conditions. The overall increase in 239 Pu average over all the islands of Atoll is about 20 folds relative to the baseline concentration. Despite the high relative increase ratios, the significantly low activity concentrations may not pose an immediate exposure risk. However, due to the limitations of the study and uncertainties/biases in the historical data used, further research supported by field surveys to better characterize the current contamination level may be needed to make more accurate predictions.

54 ENVIRONMENTAL SCIENCES↗

Quantitative assessment of Ni + and He + ion irradiation damage in a tungsten heavy alloy under the simulated nuclear fusion environment

A 90W-7Ni-3Fe (wt.%) tungsten heavy alloy has been sequentially Ni + and He + ion irradiated at 700 °C to simulate the high temperature irradiation environment of a fusion reactor interior. W/Ni–Fe-W dual-phase alloys have been proposed to serve as plasma facing materials and require detailed investigation of their behavior under fusion relevant conditions to assess their overall applicability. To evaluate material performance under five years of simulated fusion reactor service, microstructural characterization of the nanoscale defect distribution has been performed on both constituent phases, revealing peak swelling in the W phase of approximately 0.03%. The γ-phase (Ni–Fe-W) is found to swell approximately 0.68% under the same irradiation conditions, indicating significant cavity formation and growth. Additionally, a novel multi-projection imaging approach has been applied to determine the extent of damage segregation along the dual-phase W-to-γ interface and exposes that these interfaces act as sink sites for the accumulation of cavities. Interphase boundaries are noted to possess an 11.8% areal coverage of defects along the boundary plane, primarily on the γ-phase side of the boundary. The accumulation of cavities at these interphase boundaries is anticipated to adversely affect overall material toughness, and this work reveals a pressing need for mechanical property testing of irradiated W–Ni-Fe dual-phase alloys.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Porting ATLAS Fast Calorimeter Simulation to GPUs with Performance Portable Programming Models

FastCaloSim is a parameterized simulation of the particle energy response and of the energy distribution in the ATLAS calorimeter. It is a relatively small and self-contained package with massive inherent parallelism and captures the essence of GPU offloading via important operations like data transfer, memory initialization, floating point operations, and reduction. It was identified by the High Energy Physics Center for Computational Excellence project as a good testbed for evaluating the performance and ease of portability of programming models. In this paper, we will discuss the results of our evaluation of the porting process to Kokkos, SYCL, Alpaka, OpenMP and std::par (nvc++), and compare performance on NVIDIA, AMD and Intel GPUs, as well as multicore CPUs.

97 MATHEMATICS AND COMPUTING↗

Isolating solvent–solute hydrogen bonding interactions via 2D IR solvation shell spectroscopy

The solvation shell around a solute is a fundamental feature of liquid-phase solutions, determining the behavior and properties of both the solute and the overall solution. Direct experimental measurements of the solvation shell properties are challenging due to the strong signals generated from the bulk solvent, which overwhelm the small contribution of the solvation shell. Here, we use ultrafast two dimensional infrared (2D IR) spectroscopy and intermolecular cross-peaks to isolate the IR absorption spectrum of methanol molecules in the solvation shell surrounding the solute N-methylacetamide. We demonstrate that the intermolecular coupling between the solvent and solute vibrations is indirectly mediated by a low-frequency hydrogen-bonding mode, suggesting an important mechanism for anharmonic coupling induced by hydrogen bonds. From the relative frequency shifts and cross-peak anisotropy, we find that methanol molecules surrounding N-methylacetamide form stronger and distinctly oriented hydrogen bonds than those in the bulk solvent. Here, we also compare these results with the solvent spectra of the solute N,N-dimethylacetamide to investigate how solute structural changes alter the solvation shell and the contribution of N–H hydrogen bond donation. Our results are supported by molecular dynamics simulations, which provide detailed insights into the hydrogen-bonding distributions. Through these results, we demonstrate 2D solvation shell spectroscopy to be a valuable method for investigating solvation structures and dynamics without interference from the bulk solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calibrating Bayesian generative machine learning for Bayesiamplification

Recently, combinations of generative and Bayesian deep learning have been introduced in particle physics for both fast detector simulation and inference tasks. These neural networks aim to quantify the uncertainty on the generated distribution originating from limited training statistics. The interpretation of a distribution-wide uncertainty however remains ill-defined. We show a clear scheme for quantifying the calibration of Bayesian generative machine learning models. For a Continuous Normalizing Flow applied to a low-dimensional toy example, we evaluate the calibration of Bayesian uncertainties from either a mean-field Gaussian weight posterior, or Monte Carlo sampling network weights, to gauge their behaviour on unsteady distribution edges. Well calibrated uncertainties can then be used to roughly estimate the number of uncorrelated truth samples that are equivalent to the generated sample and clearly indicate data amplification for smooth features of the distribution.

97 MATHEMATICS AND COMPUTING↗

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗

Dehydration-Driven Ion Aggregation and the Onset of Gelation in ZnCl2 Solution

A minimal model of ionic aggregation in concentrated ZnCl2 is developed, guided by molecular dynamics simulations with a machine-learned potential. It explicitly incorporates solvent-site depletion, correlated chloride binding, and allows for loops within Zn–Cl clusters. Dehydration drives two coordination-controlled structural changes: a crossover at Z=2 from predominantly isolated Zn-containing species to Cl-bridged clusters and the onset of gelation near Z≈3. The model reproduces the concentration-dependent trends observed in the molecular dynamics simulations with two fitted parameters, while the cluster-size distribution at the highest concentration is consistent with three-dimensional percolation scaling over the accessible range of cluster sizes.

25 ENERGY STORAGE↗

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗

COSMIC DAWN: Distributed Analysis of Wireless at Nextscale

Distributed Analysis of Wireless at Nextscale (DAWN) is a novel simulation framework for large-scale design-space exploration (DSE) of unmodified software-defined radio (SDR) applications interacting in a scalable, high-fidelity, virtual physics environment. The software-defined nature of the coupled software-physics simulation leverages hardware emulation to permit in-depth examination and modification of not only the electromagnetic environment, including each signal in flight, but also the precise state of system software and components. DAWN supports modular, customizable physics environments allowing realistic propagation effects so that computationally efficient empirical models, reduced order/surrogate models, or large-scale, high-fidelity, site-specific simulations can be used as a propagation medium based on scenario requirements. This paper introduces DAWN’s design and initial implementation, detailing key architectural components, including the Physics Realization Engine (PhyRE), Runtime Infrastructure for Simulation Environments (RISE), and the design space exploration (DSE) suite. It concludes with demonstrations using unmodified 4G/LTE software available from srsRAN on computing resources ranging from a small cluster to ORNL’s Frontier Exascale system.

Wise, Mike [ORNL] (ORCID:0000000266120641)↗

Data Management in the Continuum: Cross-facility Object-based Data Transfers

Scientific workflows are evolving from relying on a monolithic storage subsystem at a single High-Performance Computing (HPC) facility to using geographically distributed file systems, repositories, and cloud storage. As a result, storing, accessing, transferring, and managing scientific data have become highly complex and prone to performance inefficiencies. This paper delves into these challenges by exploring an optimized end-to-end interface designed to seamlessly connect various local and remote storage systems, enabling efficient data movement of objects across HPC–Cloud and HPC–HPC environments. We showcase this capability through an object-focused data management runtime system, discuss the effects of relaxed consistency semantics in distributed object scenarios, and illustrate its application in an earthquake simulation workflow. Besides reducing the amount of data by selectively transferring regions of interest, our facility-local results achieved a speedup of 45 × over an optimized HDF5 usage and 15 × over the HDF5 with caching by using the new interface in PDC-XF.

Bez, Jean Luca↗

MINERvA Experiment Open Data Release

This is an open data product that contains all the neutrino and antineutrino data from the MINERvA experiment, packaged in a way that you can use. The ultimate package will contain both neutrino and antineutrino data, and will contain both our Low energy and Medium Energy data, where the neutrino energy distributions peaked around 3GeV and 6GeV respectively. We also provide simulated data and a way to access our uncertainties on that simulation, including flux, neutrino interaction, and detector uncertainties. The data (and simulated data) has already been pre-selected to contain either a muon candidate or an electron candidate.

Collaboration, MINERvA [Fermi National Accelerator↗

DEM simulation of the compression of crushable sand: does the initial particle shape matter?

Advances in DEM modeling, combined with high-resolution X-ray tomography, opened the way for computer models based on virtual replicas of the particles which preserve nearly all facets of their geometry. This leads to simulation advantages, but also high computational costs. Here we tackle a question stemming from this trend: how accurate should particle models be to ensure accuracy? We address this question for the case of the compression of crushable sand. LS-DEM was used to generate three models of Ottawa sand (exact replicas, ellipsoids, and spheres) from digital images of its grains. Compression-induced crushing was simulated for all sets by tracking evolving size and shape distribution. The results confirm that exact replicas provide the closest match of the measurements. However, intermediate degrees of rendering (e.g. ellipsoids preserving volume and aspect ratio of the real grains) led to satisfactory results only marginally different from those of exact replicas. In conclusion, these findings provide an example of the protocols that may be followed to identify the optimal degree of particle approximation which should be regarded as mandatory to achieve a conscious, sustainable use of computational resources.

58 GEOSCIENCES↗

Final cooling with thick wedges for a muon collider

In the final cooling stages for a muon collider, the transverse emittances are reduced while the longitudinal emittance is allowed to increase. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on using dense wedge absorbers. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We optimized the design parameters of the cooling channel and produced conceptual designs (corresponding to possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ~3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗

Reconstruction of 6D Phase Space

A ‘Phase Space’ is a mathematical representation of all possible physical states of matter of a physical volume or system. Particle beams at Fermilab exist in a six-dimensional phase space: three position dimensions, and three momentum dimensions. To ‘reconstruct’ that space means to take measurement data from virtual detectors along those beams and map out where particles are and how they’re moving, using computational methods. In this experiment, we simulated a particle beam (Fig 1) and used the position distribution of each particle at each detector to map out the magnitude of the beam’s initial momentum.

Shirlee, Jermaine [DuPage Coll.]↗

Distributions of Particles Accelerated by Strong Alfvénic Turbulence

This work presents a model for generating nonthermal power-law tails of particles’ energy probability density functions in turbulent collisionless plasmas, applicable to both nonrelativistic and relativistic scenarios. We propose that strong Alfvénic turbulence energizes plasma particles through curvature acceleration, particularly for particles with Larmor radii comparable to the scales of turbulence. When the energy density of the energized particles increases, the efficiency of the energy exchange process diminishes. As a result, the acceleration process saturates, leading to power-law distributions of particle momentum and energy. In the nonrelativistic case, the momentum probability density function scales as f(p)dp ∝ p −3 dp, while in the ultrarelativistic case, the energy probability density function scales as f(γ)dγ ∝ γ −3 dγ, where γ is the Lorentz factor. This model provides a unified framework for understanding particle acceleration in both energy regimes, complementing existing analytical approaches. The predicted scalings are consistent with available observations of energetic ion distributions in the heliosphere and with the findings from numerical simulations of ultrarelativistic particle acceleration in magnetically dominated plasma turbulence.

Alfven waves↗