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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A Consistent Spatially Adaptive Smoothed Particle Hydrodynamics Method for fluid-structure interactions

A new consistent, spatially adaptive, smoothed particle hydrodynamics (SPH) method for Fluid-Structure Interactions (FSI) is presented. The method combines several attributes that have not been simultaneously satisfied by other SPH methods. Specifically, it is second-order convergent; it allows for resolutions spatially adapted with moving (translating and rotating) boundaries of arbitrary geometries; and, it accelerates the FSI solution as the adaptive approach leads to fewer degrees of freedom without sacrificing accuracy. The key ingredients in the method are a consistent discretization of differential operators, a posteriori error estimator/distance-based criterion of adaptivity, and a particle-shifting technique. The method is applied in simulating six different flows or FSI problems. The new method's convergence, accuracy, and efficiency attributes are assessed by comparing the results it produces with analytical, finite element, and consistent SPH uniform high-resolution solutions as well as experimental data.

Hu, Wei↗
DOE OSTI2025

Surrogate modeling and optimization of the leaching process in a rare earth elements recovery plant

Critical minerals (CMs) and Rare Earth Elements (REEs) play a vital role in crucial infrastructure technologies such as renewable energy generation and batteries. Recovering them from waste materials has recently been found to significantly reduce environmental impact and supply chain costs related to these materials. In this work, we investigate surrogate modeling techniques aimed to simplify the modeling, simulation, and optimization of the leaching processes involved in CM and REE recovery flowsheets. As there is currently a lack of systematic studies on this topic, we perform extensive computational testing to ascertain which surrogate models are easier to construct and offer high predictive accuracy. Further, our results suggest that sparse quadratic models balance predictive accuracy and computational efficiency. Training and using these surrogates for global optimization of the leaching process requires two orders of magnitude fewer measurements and is up to four orders of magnitude faster than optimizing the original simulation using equation-oriented optimization or derivative-free optimization.

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

Re-splitting δ f method for electro-magnetic gyrokinetic particle-in-cell (PIC) simulation of tokamak plasmas

A new method based on a re-splitting δ f technique is proposed to mitigate the cancellation problem in electro-magnetic gyrokinetic PIC simulation. This new scheme, together with the pull-back mitigation (PBM) scheme (Mishchenko et al., 2014), have recently been implemented in the GEM code (Chen et al., 2003). These two algorithms are systematically compared with the original GEM algorithm for the simulation of the ion temperature gradient (ITG) mode, the kinetic ballooning mode (KBM), the toroidal Alfvén eigenmode (TAE) and the energetic particle mode (EPM). The three algorithms agree well for all of these modes. Moreover, both the re-splitting method and PBM can mitigate the cancellation problem and improve the computational efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A risk map methodology to assess the spatial and temporal distribution of leakage into groundwater from Geologic Carbon Storage

The risks to potable aquifers due to brine leakage through plugged and abandoned (P&A) wells is highly uncertain and a potentially significant contributor to the risk profile in Geologic Carbon Storage (GCS). This uncertainty stems from the unknown location of wells and the large variance of P&A wellbore permeability, making the spatial assessment of P&A brine leakage risk challenging. A new methodology is presented to generate “risk maps”, or spatial distributions of brine leakage risk to groundwater resources as defined with no-impact or Maximum Contaminant Level (MCL) thresholds. The methodology utilizes probability theory, thereby avoiding the use of computationally expensive Monte Carlo simulations while maintaining flexibility in modeling techniques. These maps provide quantitative probabilities of risk as a function of time to inform site selection and monitoring during and post-injection, conducive to the US EPA's permitting of class-VI wells and the so-called “area of review”, AoR. As a demonstration of the methodology, a numerical model of a hypothetical fully-coupled system spanning from the injection reservoir to the USDW is used to assess the evolution of brine leakage through P&A wells. Risk maps of CO 2 leakage can also be generated with this methodology for a comprehensive assessment of GCS leakage risk.

58 GEOSCIENCES↗
DOE OSTI2023

Bridging microscale to macroscale mechanical property measurements of FeCrAl alloys by crystal plasticity modeling

FeCrAl alloys are candidates for accident tolerant fuel cladding of light water reactors. In this work, a microstructure- and temperature-dependent crystal plasticity model is employed to bridge microscale to macroscale mechanical property measurements of FeCrAl alloys. With the visco-plastic self-consistent (VPSC) polycrystal plasticity framework, a mechanism-based single crystal plasticity (MSCP) model adopts the Arrhenius type rate equation to describe the dependence of the critical resolved shear stress for dislocation slips on their temperature-dependent intrinsic frictional resistance and the microstructure-dependent irradiation hardening. The intrinsic frictional resistance associated with {110}<111> and {112}<111> slip systems were measured by in-situ micromechanical testing on unirradiated/irradiated samples at 25-500 °C. The irradiation hardening is estimated by the Bacon-Kocks-Scattergood (BKS) model with density and size of radiation-induced defects measured from microstructural characterization. Several features associated with thermo-mechanical behavior of unirradiated/irradiated polycrystalline FeCrAl alloys are captured. High density of deformation-induced dislocations and radiation-induced defects results in obvious hardening at room temperature, which is weakened at high temperature, and facilitates damage evolution during deformation. Moreover, both high temperature and radiation-induced defects, which facilitate dislocation multiplication, trigger large hardening rate. Finally, the proposed method together with application of accelerator-based ion irradiation technique is a surrogate approach to simulate neutron damage, improving the efficiency associated with evaluation of mechanical properties of FeCrAl alloys exposed to temperature, stress and radiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A review of helium–hydrogen synergistic effects in radiation damage observed in fusion energy steels and an interaction model to guide future understanding

Here, under fusion reactor conditions, large quantities of irradiation defects and transmutation gases are produced per unit time by neutrons, resulting in accelerated degradation of structural candidate ferritic (F) and ferritic/martensitic (F/M) steels. Due to the lack of a suitable fusion neutron testing facility, we must rely on high-dose-rate ion-beam experiments and present-day crude modeling estimates. Of particular interest is the possibility of synergistic (positive feedback) effects on materials properties due to the simultaneous action of He, H, and displacement damage (dpa) during operation. In this paper we discuss the state-of-the-art in terms of the experimental understanding of synergistic effects and carry out simulations of triple-species irradiation under ion-beam conditions using first-of-its-kind modeling techniques. Although, state-of-the-art modeling and simulation is not sufficiently well developed to shed light on the experimental uncertainties, we are able to conclude that it is not clear whether synergistic effects, the evidence of which is still not conclusive, will ultimately play a critical role in material performance under fusion energy conditions. We review here some of the evidence for the synergistic effects of hydrogen in the presence of helium and displacement damage, and also include some recent data from our research. While the experimental results to date suggest possible mechanisms for the observed synergistic effects, it is only with more advanced modeling that we can hope to understand the details underlying the experimental observations. By employing modeling and simulation we propose an interaction model that is qualitatively consistent with experimental observations of dpa/He/H irradiation behavior. Our modeling, the results of which should be helpful to researchers going forward, points to gaps and voids in the current understanding of triple ion-beam irradiation effects (displacement damage produced simultaneously with helium and hydrogen implantation) and the synergistic effects of hydrogen.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Investigation of Compton scattering for gamma beam intensity measurements and perspectives at ELI-NP

Compton γ-ray sources have been used for over 30 years with new facilities being under construction or proposed. The gamma beam system under implementation at the Extreme Light Infrastructure - Nuclear Physics facility in Romania will deliver brilliant γ-ray beams with energies up to 19.5 MeV. Several instruments for measuring the parameters of the γ-ray beam are under development at ELI-NP. One of these instruments based on a High Purity Germanium detector is typically used for beam energy measurements at other facilities. Here we investigate the use of a High Purity Germanium detector to continuously monitor the intensity of the ELI-NP gamma beam by measuring the inelastic scattering of photons. This technique relies on both experimental and simulated data and it has been successfully tested during a recent experiment at the High Intensity γ-ray Source facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Hydroxyl radical yields in the heavy ion radiolysis of water

The yields of hydroxyl radicals in the radiolysis of water with protons and carbon ions have been examined using experimental scavenging techniques coupled with Monte Carlo track simulations. Combined with previous results using helium ions, this set of data gives valuable information on the potential for radiation induced damage to biological systems for light ions with very different track structures. Carbon dioxide production from aerated formic acid solutions in concentrations ranging from 10 -3 to 1 M was used as a probe of hydroxyl radical yields from about 8 ns to 8 µs. Numerical interpolation of the results at slightly different ion energies in combination with data from gamma radiolysis allows for a systematic analysis of both track average and track segment yields. As expected, considerable track chemistry is found to occur on the nanosecond to microsecond time scales. Monte Carlo track simulations employing stochastic diffusion-kinetic calculations of product yields are found to reproduce experimental observations satisfactorily. The track simulations are used to extract hydroxyl radical kinetics in pure water at neutral conditions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Immobilizing Pertechnetate in Ettringite via Sulfate Substitution

Until recently, a nuclear waste form did not exist for the immobilization and long-term storage of the anionic form of radionuclide technetium-99 (99Tc), pertechnetate (TcO 4 - ). Across international nuclear waste sites, the inability to capture TcO 4 - poses a threat to the local environment due to the anion’s long half-life and environmental mobility under environmentally oxic conditions. To combat this challenge, most waste forms would capture 99Tc using materials that reduce TcO 4 - to Tc(IV) species, e.g. 99Tc oxides or 99Tc sulfides. However, this is a short-lived solution unless reducing capacities can be regenerated to last over the long life expectancy of some nuclear waste forms (~10,000 years). In light of recent experimental observations gained from the development of cementitious waste forms and unconfirmed hypotheses from the literature, this work explores how ettringite, a common mineral formed in cementitious materials with demonstrated success in incorporating oxyanions, may be used to immobilize TcO 4 - directly and over the life-time of the waste form. Using batch precipitation experiments, solid phase characterization techniques, and ab initio molecular dynamics simulations, we demonstrate herein successful incorporation of TcO 4 - into ettringite via substitution for SO 4 2- within the ettringite crystal structure. Here, these results are first of their kind and confirm the hypothesis that ettringite acts as an immobilizing mineral within cementitious waste forms for the highly mobile TcO 4 - oxyanion, enhancing the overall performance of these waste forms over time.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Interphase Structures and Dynamics near Nanofiller Surfaces in Polymer Solutions

Here, we report the in situ structures and dynamics of hydrogenated polybutadiene (PB) chains bound to carbon black nanoparticle surfaces in polymer solutions composed of deuterated PB and deuterated toluene using small-angle neutron scattering and neutron spin-echo techniques together with molecular dynamics (MD) simulations. The experimental results showed that the swollen bound polymer chains exhibit the collective dynamics (the so-called breathing mode) at polymer concentrations ( c ) below and above the overlap polymer concentration ( c *) (i.e., 0.61 < c/c * < 1.83), where the concentration profiles of the bound polymer remained unchanged with the different c values. Interestingly, the collective dynamics slowed down by a factor of 2 compared to that in pure d-toluene when the chain lengths of the bound polymer and matrix polymer were equal. However, when the free polymer chains were longer than the bound polymer chains, the decrease in collective dynamics was not as significant. MD simulations were performed to explore the interfacial event as a whole. As a result, we found that the matrix polymer chains, whose length is equal to that of the bound polymer, can be accommodated in the bound polymer layer effectively and are “strangulated” by the bound polymer chains, while the longer matrix polymer chains only partly penetrate into the bound chains and the diffusion behavior was hardly affected compared to that in bulk.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A Combined Experimental and Computational Study on the Stability of Nanofluids Containing Metal Organic Frameworks

Computational studies on nanofluids composed of metal organic frameworks (MOFs) were performed using molecular modeling techniques. Grand Canonical Monte Carlo (GCMC) simulations were used to study adsorption behavior of 1,1,1,3,3-pentafluoropropane (R-245fa) in a MIL-101 MOF at various temperatures. To understand the stability of the nanofluid composed of MIL-101 particles, we performed molecular dynamics simulations to compute potentials of mean force between hypothetical MIL-101 fragments terminated with two different kinds of modulators in R-245fa and water. Our computed potentials of mean force results indicate that the MOF particles tend to disperse better in water than in R-245fa. The reasons for this observation were analyzed and discussed. Our results agree with experimental results indicating that the employed potential models and modeling approaches provide good description of molecular interactions and the reliabilities. Work performed by LXD was supported by the U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. Work performed by HVRA, SKN, RKM, and PBM was supported by the Office of Energy Efficiency and Renewable Energy, Geothermal Technologies Program. Pacific Northwest National Laboratory is a multiprogram national laboratory operated for DOE by Battelle.

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

Efficient linear phase contrast in scanning transmission electron microscopy with matched illumination and detector interferometry

The ability to image light elements in soft matter at atomic resolution enables unprecedented insight into the structure and properties of molecular heterostructures and beam-sensitive nanomaterials. In this study, we introduce a scanning transmission electron microscopy technique combining a pre-specimen phase plate designed to produce a probe with structured phase with a high-speed direct electron detector to generate nearly linear contrast images with high efficiency. We demonstrate this method by using both experiment and simulation to simultaneously image the atomic-scale structure of weakly scattering amorphous carbon and strongly scattering gold nanoparticles. Our method demonstrates strong contrast for both materials, making it a promising candidate for structural determination of heterogeneous soft/hard matter samples even at low electron doses comparable to traditional phase-contrast transmission electron microscopy. Ultimately, simulated images demonstrate the extension of this technique to the challenging problem of structural determination of biological material at the surface of inorganic crystals.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Planck 2013 results. XXIV. Constraints on primordial non-Gaussianity

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result f NL local = 2.7 ± 5.8, f NL equil = -42 ± 75, and f NL orth = -25 ± 39 (68% CL statistical). Non-Gaussianity is detected in the data; using skew-C &ell; statistics we find a nonzero bispectrum from residual point sources, and the integrated-Sachs-Wolfe-lensing bispectrum at a level expected in the ΛCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C &ell; , wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, three-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c s ≥ 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r D ≥ 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the four-point function in the local model τ NL < 2800 (95% CL). Taken together, these constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Shear-flow trapped-ion-mode interaction revisited. I. Influence of low-frequency zonal flow on ion-temperature-gradient driven turbulence

Collisionless trapped ion modes (CTIMs) turbulence exhibits a rich variety of zonal flow physics. The coupling of CTIMs with shear flow driven by the Kelvin-Helmholtz (KH) instability has been investigated. The work explores the parametric excitation of zonal flow modified by wave-particle interactions leading to a new type of resonant low-frequency zonal flow. The KH-CTIM interaction on zonal flow growth and its feedback on turbulence is investigated using semi-Lagrangian gyrokinetic Vlasov simulations based on a Hamiltonian reduction technique, where both fast scales (cyclotron plus bounce motions) are gyro-averaged.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗