Search NASASearch

SEARCH · Search NASA

Results for “SPACE 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 73 records · Page 4

Adaptive Computing (AC) [SWR-24-106]

The Adaptive Computing (AC) software stack supports goal-based computing, for which a simulation workload is created on the fly adapting to the results of calculations. Application-specific code defines an objective, which may be to solve an optimization problem or to train a surrogate model with minimal uncertainty. Then, the AC driver decides where in the design parameter space to run simulations to best achieve that objective. This process is iterative and online; as new data is returned from simulations, the AC driver chooses new simulations to run. The AC driver can strategically run simulations on distributed hardware resources (including high performance computing machines, cloud resources, and edge devices) to maximize throughput and obey resource constraints.

Griffin, Kevin [National Renewable Energy Laborato

Adaptive Computing (AC) (Open Source) [SWR-24-106]

The Adaptive Computing (AC) software stack supports goal-based computing, for which a simulation workload is created on the fly, adapting to the results of calculations. Application-specific code defines an objective, which may be to solve an optimization problem or to train a surrogate model with minimal uncertainty. Then, the AC driver decides where in the design parameter space to run simulations to best achieve that objective. This process is iterative and online; as new data is returned from simulations, the AC driver chooses new simulations to run. The AC driver can strategically run simulations on distributed hardware resources (including high performance computing machines, cloud resources, and edge devices) to maximize throughput and obey resource constraints.

Griffin, Kevin [National Laboratory of the Rockies

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics

Deployment and validation of predictive 6-dimensional beam diagnostics through generative reconstruction with standard accelerator elements

Understanding the 6-dimensional phase space distribution of particle beams is essential for optimizing accelerator performance. Conventional diagnostics such as use of transverse deflecting cavities offer detailed characterization but require dedicated hardware and space. Generative phase space reconstruction (GPSR) methods have shown promise in beam diagnostics, yet prior implementations still rely on such components. Here we present the first experimental implementation and validation of the GPSR methodology, realized by the use of standard accelerator elements including accelerating cavities and dipole magnets, to achieve complete 6-dimensional phase space reconstruction. Through simulations and experiments at the Pohang Accelerator Laboratory X-ray Free Electron Laser facility, we successfully reconstruct complex, nonlinear beam structures. Furthermore, we validate the methodology by predicting independent downstream measurements excluded from training, revealing the reconstruction closely resembling ground truth. This advancement establishes a pathway for predictive diagnostics across beamline segments while reducing hardware requirements and expanding applicability to various accelerator facilities.

Kim, Seongyeol [Pohang Univ. of Science and Techno

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Four-dimensional phase-space reconstruction of flat and magnetized beams using neural networks and differentiable simulations

Beams with cross-plane coupling or extreme asymmetries between the two transverse phase spaces are often encountered in particle accelerators. Flat beams with large transverse-emittance ratios are critical for future linear colliders. Similarly, magnetized beams with significant cross-plane coupling are expected to enhance the performance of electron cooling in hadron beams. Preparing these beams requires precise control and characterization of the four-dimensional transverse phase space. In this study, we employ generative phase-space reconstruction techniques to rapidly characterize magnetized and flat-beam phase-space distributions using a conventional quadrupole-scan method. The reconstruction technique is experimentally demonstrated on an electron beam produced at the Argonne Wakefield Accelerator and successfully benchmarked against conventional diagnostics techniques. Specifically, we show that predicted beam parameters from the reconstructed phase-space distributions (e.g., as magnetization and flat-beam emittances) are in excellent agreement with those measured from the conventional diagnostic methods. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS

Coordinate-based simulation of pair distance distribution functions for small and large molecular assemblies: implementation and applications

X-ray scattering has become a major tool in the structural characterization of nanoscale materials. Thanks to the widely available experimental and computational atomic models, coordinate-based X-ray scattering simulation has played a crucial role in data interpretation in the past two decades. However, simulation of real-space pair distance distribution functions (PDDFs) from small- and wide-angle X-ray scattering, SAXS/WAXS, has been relatively less exploited. This study presents a comparison of PDDF simulation methods, which are applied to molecular structures that range in size from β-cyclodextrin [1 kDa molecular weight (MW), 66 non-hydrogen atoms] to the satellite tobacco mosaic virus capsid (1.1 MDa MW, 81 960 non-hydrogen atoms). The results demonstrate the power of interpretation of experimental SAXS/WAXS from the real-space view, particularly by providing a more intuitive method for understanding of partial structure contributions. Furthermore, the computational efficiency of PDDF simulation algorithms makes them attractive as approaches for the analysis of large nanoscale materials and biological assemblies. The simulation methods demonstrated in this article have been implemented in stand-alone software, SolX 3.0 , which is available to download from https://12idb.xray.aps.anl.gov/solx.html.

Zuo, Xiaobing (ORCID:0000000201344804)

Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips

High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.

36 MATERIALS SCIENCE

GP Cosmology Surrogate v1.0

GP Cosmology Surrogate is a Python library for building and training a generalized multi-output Gaussian process (GP) framework of @takhtaganov2021cosmic. In this approach, the surrogate is constructed sequentially, guided by a Bayesian optimization acquisition function that targets reduction of emulation error in the regions most consistent with the observational data. This adaptive design concentrates computational resources where they have the greatest impact on inference accuracy. The library supports efficient training for separable GP kernels, which allows the use of Kronecker algebra to handle high-dimensional input spaces and large numbers of correlated outputs. This makes it well suited for applications such as modeling cosmological power spectra, large-scale physical simulations, and multi-output hyperparameter tuning. By combining scalable multi-output GP modeling with data-driven adaptive sampling, GPsurrogate enables parameter inference and optimization with substantially fewer simulations than conventional space-filling designs.

Lukic, Zarija [Lawrence Berkeley National Laborato

The Impact of the 8–10 March 2012 Geomagnetic Storm on Inner Zone Protons as Measured by Van Allen Probes

The Relativistic Electron Proton Telescope (REPT) instrument on the Van Allen Probes observed a double-peaked inner zone proton population throughout the 7 year lifetime of the mission. M. Hudson et al. (2023) showed that a strong SEP event accompanied by a CME-shock in early March 2012 provided the Solar Energetic Proton (SEP) source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L = 2. The study followed trajectories of SEP protons launched isotropically from a sphere at 7 Re for 2.5 hr in fields calculated by the LFM-RCM global MHD model, which includes electric fields needed to model the transport and trapping of the protons by the shock, and then a radial diffusion simulation was run for 2 years using the result from the test-particle simulation as the initial condition. The simulation result was compared with REPT measurement in November 2013 and showed reasonable agreement. However, the simulation overestimated the Phase Space Density by a factor of four due to lack of field line curvature scattering during the storm in the model. In this study, a test-particle simulation is performed for 2 days following the injection and trapping of protons in March 2012 using TS05 fields to simulate the field line curvature scattering of the trapped SEP due to the buildup of the ring current during the geomagnetic storm. The resulting sample distribution was then weighted using the flux at the end of the two-hour MHD-test particle simulation. A radial diffusion simulation is then run using the initial profile that included the loss effect, with improved comparison with REPT measurements after 2 years.

79 ASTRONOMY AND ASTROPHYSICS

Surface Segregation Studies in Ternary Noble Metal Alloys: Comparing DFT and Machine Learning with Experimental Data

Abstract Surface segregation, whereby the surface composition of an alloy differs systematically from the bulk, has historically been hard to study, because it requires experimental and modeling methods that span alloy composition space. In this work, we study surface segregation in catalytically relevant noble and platinum‐group metal alloys with a focus on three ternary systems: AgAuCu, AuCuPd, and CuPdPt. We develop a data set of 2478 fcc slabs with those compositions including all three low‐index crystallographic orientations relaxed with Density Functional Theory using the PBEsol functional with D3 dispersion corrections. We fine‐tune a machine learning model on this data and use the model in a series of 1800 Monte Carlo simulations spanning ternary composition space for each surface orientation and ternary chemical system. The results of these simulations are validated against prior experimental surface segregation data collected using composition spread alloy films for AgAuCu and AuCuPd. Our findings reveal that simulations conducted using the (110) orientation most closely match experimentally observed surface segregation trends, and while predicted trends qualitatively match observation, biases in the PBEsol functional limit numeric accuracy. This study advances understanding of surface segregation and the utility of computational studies and highlights the need for further improvements in simulation accuracy.

Chemistry

Deforestation Increases Vegetation Vulnerability to Drought Across Biomes

Land use and land cover changes have altered terrestrial ecosystem carbon storage, but their impacts on ecosystem sensitivity to drought and temperature fluctuations have not been evaluated spatially over the globe. We estimate drought and temperature sensitivities of ecosystems using vegetation greenness from satellite observations and vegetation biomass from dynamic global vegetation model (DGVM) simulations. Using a space-for-time substitution with satellite data, we first illustrate the effects of vegetation cover changes on drought and temperature sensitivity and compare them with the effects estimated from DGVMs. We also compare simulations forced by scenarios with and without land cover changes to estimate the historical land cover change effects. Satellite data and vegetation models both show that converting forests to grasslands results in a more negative or decreased positive sensitivity of vegetation greenness or biomass to drought. Significant variability exists among models for other types of land cover transitions. We identify substantial effects of historical land cover changes on drought sensitivity from model simulations with a generally positive direction globally. Deforestation can lead to either an increased negative sensitivity, as drought-tolerant forests are replaced by grasslands based on model ensemble mean, or a decreased negative sensitivity, since forests under current land cover are predicted to exhibit greater drought resistance compared to those under pre-industrial land cover. Overall, our findings emphasize the critical role of forests in maintaining ecosystem stability and resistance to drought and temperature fluctuations, thereby implying their importance in stabilizing the carbon stock under increasingly extreme climate conditions.

Xiao, Chenwei [Max Planck Institute for Biogeochem

Non-propagating structures and propagating waves in solar wind turbulence revealed by simulations and observations

Structures and waves are common features of solar wind turbulence at various scales. The interplay between structures and waves is important for processes such as the turbulent energy cascade, plasma heating, and particle scattering. Our understanding of turbulence has been advanced by not only new space missions and numerical simulations, but also techniques that have been developed to interpret the rapidly growing turbulence data. We review basic models of turbulence with a specific focus on the analysis methods for understanding magnetic structures and waves. MHD and kinetic waves in single-spacecraft time series measurements can be identified through mode decomposition or their characteristic polarization signatures. The structures in this paper are considered as zero-frequency, non-propagating or convected modes embedded in the solar wind. The synergy between observations and simulations is most evident in the application of spatial-temporal analysis to multi-spacecraft observation and turbulence simulations. The spatial-temporal analysis has greatly improved our understanding of structures and waves in turbulence. We conclude by discussing prospects for future research.

79 ASTRONOMY AND ASTROPHYSICS

Kekulé valence bond order in the honeycomb lattice optical Su-Schrieffer-Heeger model and its relevance to graphene

We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Space-time crystals from particle-like topological solitons

Time crystals are unexpected states of matter that spontaneously break time-translation symmetry either in a discrete or continuous manner. However, spatially mesoscale space-time crystals that break both space and time symmetries have not been reported. Here we report a continuous space-time crystal in a nematic liquid crystal driven by ambient-power, constant-intensity unstructured light. Our numerically constructed four-dimensional configurations exhibit good agreement with these experimental findings. Although meeting the established criteria to identify time-crystalline order, both experiments and computer simulations reveal a space-time crystallization phase formed by particle-like topological solitons. The robustness against temporal perturbations and spatiotemporal dislocations shows the stability and rigidity of the studied space-time crystals, which relates to their locally topological nature and many-body interactions between emergent spontaneously twisted, particle-like solitonic building blocks. Their potential technological utility includes optical devices, photonic space-time crystal generators, telecommunications and anti-counterfeiting designs, among others.

Liquid crystals

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY