Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

Planar optical waveguide for refractive index determining with high sensitivity and dual-band characteristic for Nano-sensor application

In this paper, we develop an optical waveguide based on metal–insulator-metal as a refractive index sensor with dual-band feature and a high figure of merit (FOM) and sensitivity. In this model, we utilize a Z-shape slot based on two stubs model in both metal layers and this structure is placed over a dielectric layer of SiO{sub 2} as a substrate for planar application in the optical systems. In addition, the parametric studies are used to clarifying the elements influence on controlling the transmission and working frequency to achieve a dual-band characteristic at 1550 and 3100 nm which is suitable for optical fiber application. In fact, the compound stubs technique is suggested to modify the resonances of the waveguide for dual-band attributes based on the transmission line model. We simulate this waveguide using the Finite Integrated Technique (FIT) method. It is supposed that the waveguide is filled by various materials with the different refractive indexes in the range of 1 to 1.5. We obtain the FOM for these materials while the Fano response provides high FOM about the 3800 RIU{sup −1} and the sensitivity is 1790 nm/RIU. Therefore, the compound stubs strategy help to improve the FOM factor in the optical waveguide.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
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 OSTI2020

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗
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↗
DOE OSTI2020

Measurement of space charge effects in the MicroBooNE LArTPC using cosmic muons

Large liquid argon time projection chambers (LArTPCs), especially those operating near the surface, are susceptible to space charge effects. In the context of LArTPCs, the space charge effect is the build-up of slow-moving positive ions in the detector primarily due to ionization from cosmic rays, leading to a distortion of the electric field within the detector. This effect leads to a displacement in the reconstructed position of signal ionization electrons in LArTPC detectors ("spatial distortions"), as well as to variations in the amount of electron-ion recombination experienced by ionization throughout the volume of the TPC. We present techniques that can be used to measure and correct for space charge effects in large LArTPCs by making use of cosmic muons, including the use of track pairs to unambiguously pin down spatial distortions in three dimensions. Finally, the performance of these calibration techniques are studied using both Monte Carlo simulation and MicroBooNE data, utilizing a UV laser system as a means to estimate the systematic bias associated with the calibration methodology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Application of adaptive multilevel splitting on coupled neutron-photon TRIPOLI-4{sup R} Monte Carlo simulations - 25427

In the context of radiation protection simulations, the Adaptive Multilevel Splitting (AMS) algorithm is a challenging variance reduction (VR) technique that has been recently investigated in the field of particle transport simulation. It has been implemented in the forthcoming version 11 of the Monte Carlo code TRIPOLI-4{sup R} and successfully tested in neutron-only and photon-only configurations. This paper addresses the application of the AMS algorithm to coupled simulations, and particularly to neutron-photon Monte Carlo calculations. The branching process occurring during the Monte Carlo coupled transport is taken into account in the new coupled-AMS algorithm and is explained in this paper. Two different neutron-photon configurations are then investigated, leading to a comparison of the coupled-AMS algorithm with the analog simulation on the one hand, and with the Exponential Transform (ET) on the other hand, which is the standard VR technique of TRIPOLI-4. Gains up to 30 are obtained in terms of Figure of Merit relatively to the analog simulation, which is about 4 to 6 times more efficient than the ET method for these configurations. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

The (Un)reasonable Effectiveness of Neural Network in Cherenkov Calorimetry

We report a greater than factor of two improvement in the hadronic energy resolution of a simulated Cherenkov calorimeter by estimating the energy with machine learning over traditional techniques. The prompt signal formation and energy threshold properties of Cherenkov radiation provide identifiable features that machine learning techniques can exploit to produce a superior model for energy reconstruction. We simulated a quartz-fiber calorimeter via the GEANT4 framework to study the reconstruction techniques in single events. We compared the machine learning-based reconstruction performance to the traditional simple sum of signal and dual-readout techniques that use both Cherenkov and scintillation signals. We describe why this game-changing approach to Cherenkov hadron calorimetry excels and our plans for a dedicated beam test to validate these findings with a fast, radiation-hard hadron calorimeter prototype.

Akchurin, Nural↗
DOE OSTI2019

Strain-Induced Structural Deformation Study of 2D Mo x W (1- x ) S 2

The potential ability of tuning properties and its possible applications in the fields of optoelectronics and/or flexible electronics, has increased the demand for 2D alloys in recent times. Understanding the mechanical performance of 2D materials under extreme conditions, such as strain, stress, and fracture is essential for the reliable electronic devices based on these structures. Here, combined molecular dynamics (MD) simulations and in situ Raman spectroscopic techniques are used to study the mechanical performance of a 2D alloy system, MoxW (1- x ) S 2 . It is determined that W substitution in MoS 2 causes solid-solution strengthening and increase in the Young's modulus values. Higher W content decreases failure strain for MoS 2 . Due to spatially resolved Raman spectroscopy and MD simulations results, a detailed model to explain failure mechanisms in MoxW (1- x ) S 2 alloys is proposed.

36 MATERIALS SCIENCE↗