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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2024

Deep Generative Models for Fast Photon Shower Simulation in ATLAS

The need for large-scale production of highly accurate simulated event samples for the extensive physics programme of the ATLAS experiment at the Large Hadron Collider motivates the development of new simulation techniques. Building on the recent success of deep learning algorithms, variational autoencoders and generative adversarial networks are investigated for modelling the response of the central region of the ATLAS electromagnetic calorimeter to photons of various energies. The properties of synthesised showers are compared with showers from a full detector simulation using GEANT4 . Both variational autoencoders and generative adversarial networks are capable of quickly simulating electromagnetic showers with correct total energies and stochasticity, though the modelling of some shower shape distributions requires more refinement. This feasibility study demonstrates the potential of using such algorithms for ATLAS fast calorimeter simulation in the future and shows a possible way to complement current simulation techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Investigating the XENON1T low-energy electronic recoil excess using NEST

The search for dark matter, the missing mass of the Universe, is one of the most active fields of study within particle physics. The XENON1T experiment recently observed a 3.5⁢𝜎 excess potentially consistent with dark matter, or with solar axions. Here, we will use the Noble Element Simulation Technique (NEST) software to simulate the XENON1T detector, reproducing the excess. We utilize different detector efficiency and energy reconstruction models, but they primarily impact sub-keV energies and cannot explain the XENON1T excess. However, using NEST, we can reproduce their excess in multiple, unique ways, most easily via the addition of 31 ± 11 37 Ar decays. Furthermore, this results in new, modified background models, reducing the significance of the excess to ≤2.2⁢𝜎 at least using non-Profile Likelihood Ratio (PLR) methods. This is independent confirmation that the excess is a real effect, but potentially explicable by known physics. Many cross-checks of our 37 Ar hypothesis are presented.

particle dark matter↗
DOE OSTI2023

Numerical Studies of the On-Line Reactivity Monitoring Techniques in Subcritical Systems: Revisit of the Area-Ratio Method, the Neutron Source-Modulation Method, and the Flux-to-Current Ratio Method

In the subcritical systems, the subcriticality level must be monitored continuously to ensure the safe operation of the systems. Techniques such as the neutron source-modulation method and the flux-to-current ratio method have been proposed for on-line reactivity measurements of the subcritical systems. However, the accuracy of these methods is unknown in the actual applications. In this report, the fundamental concepts of these techniques are visited, with numerical simulation techniques developed to validate these on-line reactivity monitoring methodologies. Specifically, the Monte Carlo numerical simulation results show that the neutron source-modulation method is not accurate in determining the reactivity of the subcritical systems. The flux-to-current ratio method can provide on-line reactivity measurements, and its spatial effects have been studied and well understood through the numerical simulations of the fuel loading of the neutron source facility of Ukraine. In addition, a modified version of the area-ratio method is proposed for the on-line measurements of the absolute reactivity of the subcritical systems. Monte Carlo numerical simulations showed that this method can provide the absolute reactivity of the subcritical systems within ~5% error, if the neutron detectors are placed close to the fuel assembly interface with the reflector zone or in the reflector zones away from the neutron source.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Structure and Dynamics of Aqueous Electrolytes Confined in 2D-TiO 2 /Ti 3 C 2 T 2 MXene Heterostructures

The synthesis of heterostructures of different two-dimensional (2D) materials offers an approach to combine advantages of different materials constituting the heterostructure and ultimately enhance their performance for applications such as electrochemical energy storage, achieving high energy, and high-power densities. Understanding the behavior of ions and solvents in confinement between these dissimilar layers is critical to understand their performance and control. Considering aqueous electrolytes, we explore the heterostructure of 2D lepidocrocite-type TiO 2 (2D-TiO 2 ) and hydroxylated or O-terminated Ti 3 C 2 MXene using ReaxFF molecular dynamics simulations and elastic/quasielastic neutron scattering techniques. Simulating a bilayer water intercalation, we find that the extent of interlayer hydration is impacted most by the surface terminations on the MXene and is marginally affected by 2D-TiO 2 . However, the introduction of 2D-TiO 2 decreases the water self-diffusion due to the notch sites (i.e., surface oxygen ridges) entrapping water molecules. Intercalating alkali cations into the heterostructures, we find that Li + is predominantly adsorbed at the 2D-TiO 2 surface instead of the MXenes with the preferential occupation of the notch sites. In contrast, Na+ forms a planar solvation with water, while K + is adsorbed both at the O-terminated MXene and 2D-TiO 2 . This behavior is altered when OH-terminated MXene is involved—the repulsion from the protons on the MXene surface forces the K + ions to be adsorbed exclusively to 2D-TiO 2 , while Na + retains some of its solvation in the water layer due to its smaller size. In OH-terminated MXenes, we see a consistent transfer of protons from the MXene surface toward 2D-TiO 2 , implying a greater capacity to store protons in the heterostructures. Of the three cations simulated, Na + hinders the proton migration the least and Li + the most because of its position near the 2D-TiO 2 surface. Furthermore, 2D-TiO 2 /MXene heterostructures are likely to exhibit a higher energy density but lower power density, especially with Na + intercalation.

25 ENERGY STORAGE↗
DOE OSTI2020

CsI(Tl) pulse shape discrimination with the Belle II electromagnetic calorimeter as a novel method to improve particle identification at electron–positron colliders

Here we describe the implementation and performance of CsI(Tl) pulse shape discrimination for the Belle II electromagnetic calorimeter, representing the first application of CsI(Tl) pulse shape discrimination for particle identification at an electron–positron collider. The pulse shape characterization algorithms applied by the Belle II calorimeter are described. Control samples of $γ, μ^+, π^±, K^±$ and $ρ/ \bar{ρ}$ are used to demonstrate the significant insight into the secondary particle composition of calorimeter clusters that is provided by CsI(Tl) pulse shape discrimination. Comparisons with simulation are presented and provide further validation for newly developed CsI(Tl) scintillation response simulation techniques, which when incorporated with GEANT4 simulations allow the particle dependent scintillation response of CsI(Tl) to be modelled. Comparisons between data and simulation also demonstrate that pulse shape discrimination can be a new tool to identify sources of improvement in the simulation of hadronic interactions in materials. The $K^0_L$ efficiency and photon-as-hadron fake-rate of a multivariate classifier that is trained to use pulse shape discrimination is presented and comparisons are made to a shower-shape based approach. CsI(Tl) pulse shape discrimination is shown to reduce the photon-as-hadron fake-rate by over a factor of 3 at photon energies of 0.2 GeV and over a factor 10 at photon energies of 1 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Review of Basic Energy Reconstruction Techniques in Liquid Xenon and Argon Detectors for Dark Matter and Neutrino Physics Using NEST

Detectors based upon the noble elements, especially liquid xenon as well as liquid argon, as both single- and dual-phase types, require reconstruction of the energies of interacting particles, both in the field of direct detection of dark matter (weakly interacting massive particles WIMPs, axions, etc.) and in neutrino physics. Experimentalists, as well as theorists who reanalyze/reinterpret experimental data, have used a few different techniques over the past few decades. In this paper, we review techniques based on solely the primary scintillation channel, the ionization or secondary channel available at non-zero drift electric fields, and combined techniques that include a simple linear combination and weighted averages, with a brief discussion of the application of profile likelihood, maximum likelihood, and machine learning. Comparing results for electron recoils (beta and gamma interactions) and nuclear recoils (primarily from neutrons) from the Noble Element Simulation Technique (NEST) simulation to available data, we confirm that combining all available information generates higher-precision means, lower widths (energy resolution), and more symmetric shapes (approximately Gaussian) especially at keV-scale energies, with the symmetry even greater when thresholding is addressed. Near thresholds, bias from upward fluctuations matters. For MeV-GeV scales, if only one channel is utilized, an ionization-only-based energy scale outperforms scintillation; channel combination remains beneficial. We discuss here what major collaborations use.

argon↗
DOE OSTI2022

Cold-atom quantum simulator for string and hadron dynamics in non-Abelian lattice gauge theory

We propose an analog quantum simulator for simulating real-time dynamics of (1+1)-dimensional non-Abelian gauge theory well within the existing capacity of ultracold-atom experiments. The scheme calls for the realization of a two-state ultracold fermionic system in a one-dimensional bipartite lattice, and the observation of subsequent tunneling dynamics. Being based on the loop string hadron formalism of SU(2) lattice gauge theory, this simulation technique is completely SU(2) invariant and simulates accurate dynamics of physical phenomena such as string breaking and/or pair production. The scheme is scalable and particularly effective in simulating the theory in the weak-coupling regime, and also a bulk limit of the theory in the strong-coupling regime up to certain approximations. This paper also presents a numerical benchmark comparison of the exact spectrum and real-time dynamics of lattice gauge theory to that of the atomic Hamiltonian with an experimentally realizable range of parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Integrating Marine Hydrokinetic and Offshore Wind Energy: A Review of Technologies, Deployment, and Challenges

Together, offshore wind (OSW) and marine hydrokinetic (MHK) technologies have vast potential to expand the world’s access to abundant energy resource. With more than 60 GW of offshore wind energy capacity and 527 MW of ocean energy deployed globally by 2023, there is a significant amount of available resources; however, technical and non-technical challenges prevent the combined large-scale deployment of these technologies. There is still a lack of research that provides a parallel review of both MHK and OSW technologies in order to better understand their synergistic working principles. This paper aimed to address that research gap by presenting a comprehensive side-by-side review of the worldwide technological landscape, global deployment trends, integration strategies, and modeling approaches for MHK and OSW. A particular focus has been given on analyzing existing modeling and simulation techniques, assessing integration and control strategies, and comparing technologies based on water depth. Furthermore, this study provides important insights into the readiness levels of both technologies by highlighting ongoing international projects. By addressing these issues, this review will give researchers and industry stakeholders an outline for assessing the maturity of OSW and MHK systems and facilitating their transition to large-scale, sustainable deployment.

16 - TIDAL AND WAVE POWER↗
DOE OSTI2024

Analysis of dislocation configurations in SiC crystals through X-ray topography aided by ray tracing simulations

Silicon carbide as a wide bandgap semiconductor is of great research interest for its widespread deployment in a range of electronic and optoelectronic devices, particularly in power electronics. However, defects in silicon carbide crystals are still major concerns that is hampering the development of high-performance devices. X-ray topography, particularly using the synchrotron beam has been instrumental in characterizing and analyzing defect configurations in silicon carbide crystals to optimize crystal growth as well as understand the effect of defects on device performance. Here, in recent years, the use of ray-tracing simulation technique based on the orientation contrast mechanism to simulate contrast of defects observed on actual X-ray topographs has proven to be an effective approach to investigate the nature of crystallographic defects in various semiconductors. This review discusses the principle of ray-tracing simulation and its application and modifications to incorporate the effects of surface relaxation and photoelectric absorption to better simulate different dislocations observed in 4H–SiC as well as 6H–SiC crystals of various orientations. The adaptation to weak beam topography and plane wave topography is also discussed. The application of ray-tracing simulation in dislocation characterization of silicon carbide of different polytypes is systematically reviewed including different types of dislocations observed in both off-axis wafers and axial-sliced samples through synchrotron X-ray topography under various beam conditions, recording geometries and reflections. The result of ray-tracing simulation is further utilized in other studies including the investigation of effective penetration depth of all types of dislocations lying on the basal plane on grazing-incidence X-ray topography.

Engineering↗
DOE OSTI2020

Confining Liquids inside Carbon Nanotubes: Accelerated Molecular Dynamics with Spliced, Soft-Core Potentials and Simulated Annealing

Understanding emergent phenomena of fluids under physical confinement requires the development of advanced tools for rapid and accurate simulation of their physiochemical properties. Simulating liquid molecules commensurate in size with the nanoscale enclosures that confine them is a key challenge. In this work, we demonstrate an accelerated molecular dynamics simulation technique that combines soft-core potentials (SCP) and simulated annealing (SA) to analyze confined liquids. This integrated SCP/SA method relies on a new spliced soft-core potential (SSCP), which enables tunable accuracy with respect to the target hard-core potential (HCP). SCP/SA enables the packing of enclosures with bulk material in a controlled, thermodynamically consistent manner. The enhanced SSCP accuracy is a critical feature of SCP/SA, enabling a smooth transition between the SCP and the HCP at a desired SCP hardness. We applied SCP/SA to the problem of filling a carbon nanotube (CNT) in periodic boundary conditions with a popular ionic liquid (IL), 1- butyl-3-methylimidazolium hexafluorophosphate $[\text{BMIM}^+][\text{PF}^–_6]$. We performed a series of triplicate simulations on systems with varying CNT diameter and charge to demonstrate SCP/SA’s versatility. Beyond this IL/CNT system, the SCP/SA simulation framework has a broad range of potential applications, not limited to nanoscale enclosures and interfaces, including both solid-state and biological systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Supporting New Advanced Nuclear Technologies for Commercial-Maritime Applications

The ANS summary doesn't require an abstract, but I will produce one for the purpose of LRS: The large demand for maritime nuclear power underscores the need for experimental campaigns and modeling and simulation of new advanced reactors, which offer numerous advantages in terms of safety, efficiency, and compactness. INL, through the work conducted by NRIC and ABS, has addressed some of the technical, regulatory, and economic aspects of potential nuclear commercial maritime applications. However, on the technical side, there remain important physical phenomena, particularly for advanced reactors, that are not yet fully understood. Addressing these knowledge gaps requires a combination of experiments and advanced modeling and simulation techniques. INL possesses significant expertise in Multiphysics modeling and simulation. By collaborating with INL, the maritime nuclear sector can leverage this expertise to advance the development and deployment of innovative nuclear technologies.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Phase behavior of continuous-space systems: A supervised machine learning approach

The phase behavior of complex fluids is a challenging problem for molecular simulations. Supervised machine learning (ML) methods have shown potential for identifying the phase boundaries of lattice models. In this work, we extend these ML methods to continuous-space systems. We propose a convolutional neural network model that utilizes grid-interpolated coordinates of molecules as input data of ML and optimizes the search for phase transitions with different filter sizes. We test the method for the phase diagram of two off-lattice models, namely, the Widom–Rowlinson model and a symmetric freely jointed polymer blend, for which results are available from standard molecular simulations techniques. The ML results show good agreement with results of previous simulation studies with the added advantage that there is no critical slowing down. We find that understanding intermediate structures near a phase transition and including them in the training set is important to obtain the phase boundary near the critical point. The method is quite general and easy to implement and could find wide application to study the phase behavior of complex fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗