Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation technique”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Third-order photon correlations extract single-nanocrystal multiexciton properties in solution

Colloidal semiconductor nanocrystals are considered promising materials for high-flux optical applications, including lasing, light-emitting diodes, biological imaging, and quantum optics. In high-flux applications, multiexcitons can significantly contribute to emission, influencing its brightness, spectral purity, and kinetics. As a result, understanding and controlling multiexciton emission in colloidal nanocrystal materials is of the utmost importance. In the past, single-nanocrystal photon correlation methods have been applied to understand biexciton and triexciton efficiencies, lifetimes, and spectra. While powerful, such methods suffer from user selection bias and require stable emission from single nanocrystals. To compensate for this shortcoming, second-order correlation methods were developed to extract sample-averaged biexciton properties from a solution of nanocrystals. Until now, however, the analogous third-order solution photon correlation methods remained unexplored. In this work, we present a pair of third-order photon correlation techniques to obtain the sample-averaged single-nanocrystal triexciton quantum yield and lifetime in a solution-phase experiment. These techniques derive from the relationship between the Poisson probability of nanocrystal photon absorption and the intrinsic probability of nanocrystal photon emission. We validate the theoretical background of these techniques by creating a numerical model to simulate the diffusion and emission of many nanocrystals in solution. Our simulations confirm that the average triexciton quantum yield and triexciton lifetime can be extracted from a solution of nanocrystals. These techniques will enable researchers to gain a better understanding of the fundamental multiexciton properties of colloidal nanocrystals.

Horowitz, Jonah R. [Massachusetts Institute of Tec↗
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Adamantine 1.0: A Thermomechanical Simulator for Additive Manufacturing

Adamantine is a thermomechanical simulation code that is written in C++ and built on top of deal.II (Arndt et al., 2023), p4est (Burstedde et al., 2011), ArborX (Lebrun-Grandié et al., 2020), Trilinos (The Trilinos Project Team, 2020), and Kokkos (Trott et al., 2022). Adamantine was developed with additive manufacturing in mind and it is particularly well adapted to simulate fused filament fabrication, directed energy deposition, and powder bed fusion. Adamantine employs the finite element method with adaptive mesh refinement to solve a nonlinear anisotropic heat equation, enabling support for various additive manufacturing processes. It can also perform elastoplastic and thermoelastoplastic simulations. It can handle materials in three distinct phases (solid, liquid, and powder) to accurately reflect the physical state during different stages of the manufacturing process. To enhance simulation accuracy, adamantine incorporates data assimilation techniques (Asch et al., 2016). This allows it to integrate experimental data from sensors like thermocouples and infrared (IR) cameras. This combined approach helps account for errors arising from input parameters, material properties, models, and numerical calculations, leading to more realistic simulations that reflect what occurs in a particular print.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

At the beginning of the last decade, Petascale supercomputers (i.e., computers capable of more than 1 petaFLOP) emerged. Now, at the dawn of exascale supercomputing, we provide a review of recent landmark simulations of portions of reactor components with turbulence-resolving techniques that this computational power has made possible. In fact, these simulations have provided invaluable insight into flow dynamics, which is difficult or often impossible to obtain with experiments alone. We focus on simulations performed with the spectral element method, as this method has emerged as a powerful tool to deliver massively parallel calculations at high fidelity by using large eddy simulation or direct numerical simulation. We also limit this paper to constant-property incompressible flow of a Newtonian fluid in the absence of other body or external forces, although the method is by no means limited to this class of flows. We briefly review the fundamentals of the method and the reasons it is compelling for the simulation of nuclear engineering flows. We review in detail a series of Petascale simulations, including the simulations of helical coil steam generators, fuel assemblies, and pebble beds. Even with Petascale computing, however, limitations for nuclear modeling and simulation tools remain. In particular, the size and scope of turbulence-resolving simulations are still limited by computing power and resolution requirements, which scale with the Reynolds number. In the final part of this paper, we discuss the future of the field, including recent advancements in emerging architectures such as GPUbased supercomputers, which are expected to power the next generation of high-performance computers.

computational fluid dynamics↗
DOE OSTI2020

Imaging phonon dynamics with ultrafast electron microscopy: Kinematical and dynamical simulations

Ultrafast x-ray and electron scattering techniques have proven to be useful for probing the transient elastic lattice deformations associated with photoexcited coherent acoustic phonons. Indeed, femtosecond electron imaging using an ultrafast electron microscope (UEM) has been used to directly image the influence of nanoscale structural and morphological discontinuities on the emergence, propagation, dispersion, and decay behaviors in a variety of materials. Here, we describe our progress toward the development of methods ultimately aimed at quantifying acoustic-phonon properties from real-space UEM images via conventional image simulation methods extended to the associated strain-wave lattice deformation symmetries and extents. Using a model system consisting of pristine single-crystal Ge and a single, symmetric Lamb-type guided-wave mode, we calculate the transient strain profiles excited in a wedge specimen and then apply both kinematical- and dynamical-scattering methods to simulate the resulting UEM bright-field images. While measurable contrast strengths arising from the phonon wavetrains are found for optimally oriented specimens using both approaches, incorporation of dynamical scattering effects via a multi-slice method returns better qualitative agreement with experimental observations. Contrast strengths arising solely from phonon-induced local lattice deformations are increased by nearly an order of magnitude when incorporating multiple electron scattering effects. We also explicitly demonstrate the effects of changes in global specimen orientation on the observed contrast strength, and we discuss the implications for increasing the sophistication of the model with respect to quantification of phonon properties from UEM images.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Non-intrusive nonlinear model reduction via machine learning approximations to low-dimensional operators

Abstract Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters; this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $$10^3{\times }$$ 10 3 × in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs. The domain of applications include both parabolic and hyperbolic PDEs, regardless of the dimension of full-order models (FOMs).

42 ENGINEERING↗
DOE OSTI2025

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of Large Scale Extrusion Deposition for Structural Applications

Large Scale Extrusion Deposition (LSED) is an evolving additive manufacturing (AM) technology that research, such as that taking place at Oak Ridge National Laboratory (ORNL) and companies like Local Motors, are continuing to utilize and develop new applications. A major LSED application of interest has been molds and tooling as it allows for much shorter production time and lower cost. However, interest has been expanding into using LSED for more structural applications due to more frequent use of high performing polymer composites as feedstock. The use of LSED for structural applications is of particular interest to Local Motors as it is currently being used to create commercially viable, energy efficient electric vehicles. LSED offers a unique opportunity when compared to traditional manufacturing methods as it can significantly reduce the number of components necessary, while decreasing embodied energy and carbon emissions. Even with these benefits however, it is important that LSED is properly understood from a structural aspect as this field has not been as heavily researched as tooling. To ensure a high level of safety and repeatability is an essential responsibility of a manufacturer whose products are structural in nature. As such it is important to understand the materials and LSED process to make structural objects that the manufacturer can be confident in. The main goals of this project were to: develop and investigate materials that are of interest for structural LSED applications, further develop and understand the current machines used in LSED and develop simulations tools of LSED and the mechanical properties of the created structure. Material development focused on composite materials that have high mechanical properties and are stable in a variety of environments. The materials were tested to determine their as printed mechanical and thermal properties as these are necessary for simulations. Once the materials were investigated thoroughly, it allowed for simulations to be performed to compare to experimental data with simulations. Materials were also vetted to determine candidates for multi-material printing. The machine development focused around the areas of process monitoring, non-destructive evaluation, and quality control. Finally, the goal of the simulations was to develop a realistic model of printed structures, including in-process simulation and dynamic simulation. The routes to get to some of these goals and the depth in which they were investigated changed throughout the project due to personnel changes and the COVID-19 pandemic. This project resulted in many valuable results such as the development of an nondestructive evaluation (NDE) technique for interlayer defects, proof of simulation for warpage in simple parts, the development and utilization of a profilometer to monitor a print for defects or inconsistencies, thorough investigation of a material used commercially for structural LSED applications and valuable experimental data on the applicability and advantages of multi-material crush structures versus their neat counterparts by creation and testing of samples.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗
DOE OSTI2020

Improving inverse Compton sources by avoiding non-linearities

We present a new, more nuanced understanding of non-linear effects in inverse Compton sources. Deleterious non-linear effects can arise even at low laser intensities, a regime previously viewed as linear. After laying out a survey of non-linear phenomena which degrade the effectiveness of inverse Compton sources, we discuss two powerful techniques designed to avoid these non-linearities. Starting with the known technique of non-linear longitudinal chirping of the laser pulse in the high laser field regime, we show that the simple stretching of the laser pulse, while keeping the energy constant, can significantly increase the spectral density of the scattered radiation in many operating regimes. Our numerical simulations show that combining these two techniques avoids detrimental non-linearities and improves the performance of inverse Compton sources over an order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis

The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey Year 3 BAO sample over its full photometric redshift range 0.6 < zphoto < 1.1. The mocks are based upon 488 ICE-COLA fast N-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution – halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function w(θ), angular power spectrum Cℓ, and projected clustering ξp(r⊥) to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.Key words: catalogs / large-scale structure of Universe / galaxies: distances and redshifts / Galaxy: halo / methods: numerical

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

¬¬Integration of Quantification of Margins and Uncertainties Methodology into Parallel Discrete Event Simulator Framework

Parallel Discrete Events Simulation (PDES) is becoming increasingly important to lab efforts in security and intelligence. It is used to model complex asynchronous systems such as computer networks, satellite systems, vehicular traffic, and human performance. Uncertainty Quantification (UQ) techniques have become a mainstay of Verification and Validation (V&V) efforts on physics simulations. However, PDES models are very different from traditional physics simulations, and research into UQ techniques for PDES is in its infancy. It is not clear which traditional UQ techniques can be applied to PDES, or what new techniques will need to be developed. The goal of this project was to identify existing UQ techniques that can be applied to PDES, and to develop new techniques as necessary. The project implemented or developed techniques to handle issues that do not appear in traditional physics simulations, but are common among PDES, including sampling techniques for high-dimensional homogenous inputs, response surfaces for high-variance heteroskedastic output, and characterization of skewed output distributions. These are foundational UQ techniques that must be used in any complete UQ analysis (Tong 2018).

97 MATHEMATICS AND COMPUTING↗