Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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

Impact of dynamic specimen shape evolution on the atom probe tomography results of doped epitaxial oxide multilayers: Comparison of experiment and simulation

The experimental atom probe tomography (APT) results from two different specimen orientations (top-down and sideways) of a high oxygen ion conducting Samaria-doped-ceria/Scandia-stabilized-zirconia multilayer thin film solid oxide fuel cell electrolyte was compared with level-set method based field evaporation simulations for the same specimen orientations. This experiment-simulation comparison explains the dynamic specimen shape evolution and ion trajectory aberrations that can induce density artifacts in final reconstruction, leading to inaccurate estimation of interfacial intermixing. This study highlights the importance of comparing experimental results with field evaporation simulations when using APT to study oxide heterostructure interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The structure of liquid water up to 360 MPa from x-ray diffraction measurements using a high Q-range and from molecular simulation

X-ray diffraction measurements of liquid water are reported at pressures up to 360 MPa corresponding to a density of 0.0373 molecules per Å{sup 3}. The measurements were conducted at a spatial resolution corresponding to Q{sub max} = 16 Å{sup −1}. The method of data analysis and measurement in this study follows the earlier benchmark results reported for water under ambient conditions having a density of 0.0333 molecules per Å{sup 3} and Q{sub max} = 20 Å{sup −1} [J. Chem. Phys. 138, 074506 (2013)] and at 70 °C having a density of 0.0327 molecules per Å{sup 3} and Q{sub max} = 20 Å{sup −1} [J. Chem. Phys. 141, 214507 (2014)]. The structure of water is very different at these three different T and P state points and thus they provide the basis for evaluating the fidelity of molecular simulation. Measurements show that at 360 MPa, the 4 waters residing in the region between 2.3 and 3 Å are nearly unchanged: the peak position, shape, and coordination number are nearly identical to their values under ambient conditions. However, in the region above 3 Å, large structural changes occur with the collapse of the well-defined 2nd shell and shifting of higher shells to shorter distances. The measured structure is compared to simulated structure using intermolecular potentials described by both first-principles methods (revPBE-D3) and classical potentials (TIP4P/2005, MB-pol, and mW). The DFT-based, revPBE-D3, method and the many-body empirical potential model, MB-pol, provide the best overall representation of the ambient, high-temperature, and high-pressure data. The revPBE-D3, MB-pol, and the TIP4P/2005 models capture the densification mechanism, whereby the non-bonded 5th nearest neighbor molecule, which partially encroaches the 1st shell at ambient pressure, is pushed further into the local tetrahedral arrangement at higher pressures by the more distant molecules filling the void space in the network between the 1st and 2nd shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A coarse-graining approach for molecular simulation that retains the dynamics of the all-atom reference system by implementing hydrodynamic interactions

We report on a new approach for deriving coarse-grained intermolecular forces that retains the frictional contribution that is often discarded by conventional coarse-graining methods. The approach is tested for water and an aqueous glucose solution, and the results from the new implementation for coarse-grained molecular dynamics simulation show remarkable agreement with the dynamics obtained from reference all-atom simulations. The agreement between the structural properties observed in the coarse-grained and all-atom simulations is also preserved. We discuss how this approach may be applied broadly to any existing coarse-graining method where the coarse-grained models are rigorously derived from all-atom reference systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modulating factors of hydrologic exchanges in a large-scale river reach: Insights from three-dimensional computational fluid dynamics simulations

Hyporheic exchange is a critical mechanism shaping hydrological and biogeochemical processes along a river corridor. Recent studies on quantifying the hyporheic exchange were mostly limited to local scales due to field inaccessibility, computational demand, and complexity of geomorphology and subsurface geology. Surface flow conditions and subsurface physical properties are well known factors on modulating the hyporheic exchange, but quantitative understanding of their impacts on the strength and direction of hyporheic exchanges at reach scales is absent. In this study, a high resolution computational fluid dynamics (CFD) model that couples surface and subsurface flow and transport is employed to simulate hyporheic exchanges in a 7-km long reach along the main-stem of the Columbia River. Assuming that the hyporheic exchange does not affect surface water flow conditions due to its negligible magnitude compared to the volume and velocity of river water, we developed a one-way coupled surface and subsurface water flow model using the commercial CFD software STAR-CCM+. The model integrates the Reynolds-averaged Navier-Stokes (RANS) equation solver with a realizable k-e two-layer turbulence model, a two-layer all y^+ wall treatment, and the volume of fluid (VOF) method, and is used to simulate hyporheic exchanges by tracking the free water-air interface as well as flow in the river and the subsurface porous media. The model is validated against measurements from acoustic Doppler current profiler (ADCP) in the stream water and hyporheic fluxes derived from a set of temperature profilers installed across the riverbed. The validated model is then employed to systematically investigate how hyporheic exchanges are influenced by surface water fluid dynamics strongly regulated by upstream dam operations, as well as subsurface structures (e.g. thickness of riverbed and subsurface formation layers) and hydrogeological properties (e.g. permeability). The results suggest that the thickness of riverbed alluvium layer is the dominant factor for reach-scale hyporheic exchanges, followed by the alluvium permeability, the depth of the underlying impermeable layer, and the assumption of hydrostatic pressure.

Bao, Jie↗

Simulation of Heterogeneous Atom Probe Tip Shape Evolution during Field Evaporation Using a Level Set Method and Different Evaporation Models

Accurate reconstruction of spatial positions of field evaporated ions from detector patterns depends on a correct understanding of the tip shape evolution and evaporation laws of component atoms. Reconstruction artifacts introduced in the heterogeneous tip can be contributed by the assumption on homogeneous evaporation behavior. A level set method based specimen shape evolution model is developed in this study to simulate the evaporation of synthetic layered-structured atom probe tip. The simulation results of the shape evolution by the present level set model agree well with the literature data using other numerical approaches, i.e. finite difference method. The asymmetric evolving shape predicted by the level set model demonstrates that the complex evaporation behavior of heterogeneous tips and interface curvature can potentially lead to the tip reconstruction distortion. Compared with other shape evaporation methods the new method provides considerably higher resolution and computational efficiency with the aid of the intrinsic subgrid scale accuracy.

atom probe tomography, field evaporation, tip shap↗

One-dimensional particle simulations of Knudsen-layer effects on D-T fusion

Particle simulations are used to solve the fully nonlinear, collisional kinetic equation describing the interaction of a high-temperature, high-density, deuterium-tritium plasma with absorbing boundaries, a plasma source, and the influence of kinetic effects on fusion reaction rates. Both hydrodynamic and kinetic effects influence the end losses, and the simulations show departures of the ion velocity distributions from Maxwellian due to the reduction of the population of the highest energy ions (Knudsen-layer effects). The particle simulations show that the interplay between sources, plasma dynamics, and end losses results in temperature anisotropy, plasma cooling, and concomitant reductions in the fusion reaction rates. However, for the model problems and parameters considered, particle simulations show that Knudsen-layer modifications do not significantly affect the velocity distribution function for velocities most important in determining the fusion reaction rates, i.e., the thermal fusion reaction rates using the local densities and bulk temperatures give good estimates of the kinetic fusion reaction rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuum stress intensity factors from atomistic fracture simulations

Stress intensity factors (SIFs) are utilized in continuum fracture mechanics to quantify the stress fields surrounding a crack in a homogeneous material in the linear elastic regime. Critical values of the SIFs define an intrinsic measure of the resistance of a material to propagate a crack. At atomic scales, yet, fracture occurs as a series of atomic bonds breaking, differing from the continuum description. As a consequence, a formal analog of the continuum SIFs calculated from atomistic simulations can have spatially localized, microstructural contributions that originate from varying bond configurations. The ability to characterize fracture at the atomic scale in terms of the SIFs offers both an opportunity to probe the effects of chemistry, as well as how the addition of a microstructural component affects the accuracy. We present a novel numerical method to determine SIFs from molecular dynamics (MD) simulations. The accuracy of this approach is first examined for a simple model, and then applied to atomistic simulations of fracture in amorphous silica. MD simulations provide time and spatially dependent SIFs, with results that are shown to be in good agreement with experimental values for fracture toughness in silica glass.

42 ENGINEERING↗

Dehalogenation of electrochemical processing salt simulants with ammonium phosphates and immobilization of salt cations in an iron phosphate glass waste form

This study demonstrates the conversion of simulated pyroprocessing salt wastes [i.e., ER(SF) and ERV2] into chemically durable iron phosphate glass waste forms through reactions with ammonium hydrogen phosphate precursors [i.e., NH4H2PO4 and (NH4)2HPO4] and added Fe2O3. The reaction of the salt simulants with the phosphate precursors resulted in dehalogenation of the salt through the production of NH4Cl and NH4I (along with water), of which the NH4Cl can be used to produce UCl3 oxidant for use in the electrorefiner during pyroprocessing. These iron phosphate glasses contain a notably higher waste loading than the baseline glass-bonded sodalite ceramic waste form, which requires about 70 vol% sodalite to accommodate Cl from the salt, because they are dehalogenated prior to vitrification. The high waste loadings coupled with higher bulk densities (~3×103 kg/m3) than glass-bonded sodalite (~2.3×103 kg/m3) mean that these glasses have a higher storage volume, i.e., the volume required to immobilize a given mass of waste. The iron phosphate glass waste forms were characterized using a variety of techniques including physical property measurements and chemical durability testing. Following initial testing, a proof-of-concept demonstration was conducted to show that the NH4Cl generated from this process could be used to produce UCl3 using depleted uranium dendrites as a feedstock.

phosphate glass, pyroprocessing, waste forms, salt↗

Translational diffusion of water inside hydrophobic carbon micropores studied by neutron spectroscopy and molecular dynamics simulation

When water molecules are confined to nanoscale spacings, such as in the nanometer-size pores of activated carbon fiber (ACF), their freezing point gets suppressed down to very low temperatures (~150 K), leading to a metastable liquid state with remarkable physical properties. Here we have investigated the ambient pressure diffusive dynamics of water in microporous Kynol ACF-10 (average pore size ~11.6 Å, with primarily slit-like pores) from temperature T = 280 K in its stable liquid state down to T = 230 K into the metastable supercooled phase. The observed characteristic relaxation times and diffusion coefficients are found to be, respectively, higher and lower than those in bulk water, indicating a slowing down of the water mobility with decreasing temperature. The observed temperature-dependent average relaxation time (${{\tau}}$) when compared to previous findings indicate that it is the width of the slit pores-not their curvature-that primarily affects the dynamics of water for pore sizes larger than 10 Å. The experimental observations are compared to complementary molecular dynamics simulations of a model system, in which we studied the diffusion of water within the 11.6 Å gap of two parallel graphene sheets. We find generally a reasonable agreement between the observed and calculated relaxation times at the low momentum transfer Q (Q ≤ 0.9 Å -1) . At high Q, however, where localized dynamics becomes relevant, this ideal system does not satisfactorily reproduce the measurements. Consequently, the simulations are compared to the experiments at low Q, where the two can be best reconciled. The best agreement is obtained for the diffusion parameter D associated with the hydrogen-site when a representative stretched exponential function, rather than the standard bimodal exponential model, is used to parametrize the self-correlation function I (Q,t).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diagnosing Model Errors in Simulation of Solar Radiation on Inclined Surfaces

Transposition models have been widely used in the solar energy industry to simulate solar radiation on inclined PV panels. Following numerous studies comparing the performance of transposition models, this paper aims to understand the quantitative uncertainty in the state-of-the-art transposition models and the sources leading to the uncertainty. Our results show significant differences between two highly used isotropic transposition models with one substantially underestimating the diffuse plane-of-array (POA) irradiances when diffuse radiation is perfectly isotropic. In the empirical transposition models, the selection of empirical coefficients and land surface albedo can both result in uncertainty in the output. This study can be used as a guide for future development of physics-based transposition models.

14 SOLAR ENERGY↗