Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermophysical Property Measurements: Improved Density, Viscosity and Thermal Diffusivity Methods

The technical bases and methods developed at Argonne to measure thermophysical properties of molten salt mixtures are summarized and representative results provided to show the accuracy and precision achieved. Three key measurements addressed in this report are density, viscosity and thermal diffusivity. The densities of two eutectic FLiNaK mixtures were measured at temperatures to approximately 750 °C and the results are compared with literature values. The results are used to show how the effect of salt surface tension is taken into account in the density calculation and how the surface tension for other salt mixtures is being measured. Viscosity measurements of eutectic FLiNaK are shown to illustrate improvements to the measurement method being made to address and minimize the sources of error and uncertainties in the reported values. Thermal diffusivity measurements of FLiNaK are shown to illustrate the measurement method developed at Argonne and the achieved precision. The thermal diffusivity of the graphite used to make the cells that contain the salts during analyses was measured to assess the relative contribution to the response. The improved procedures developed in this work are being applied in property measurements of other salts of interest.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three-dimensional electrostatic simulations to parameterize the perturbed electric field profile for nonideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These nonideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the nonideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern; the distortion of the electron diffusion size and shape; and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread functions and the sublinear variance in flats. The method is implemented in the Photon Simulator and the code is publicly available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

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

Isotope effects on the high pressure viscosity of liquid water measured by differential dynamic microscopy

In this paper, differential dynamic microscopy is performed in diamond anvil cells to measure the viscosity of water along the 24 °C isotherm to high-pressure by determination of the tracer diffusion coefficient of monodisperse silica spheres of known diameter and application of the Stokes-Einstein-Sutherland equation. This technique allows liquid samples to be compressed to greater pressure prior to freezing than with other viscometry methods. The highest-pressure measurement was made at 1.67 GPa, considerably deeper into the supercompressed regime than previously reported. The effect of isotopic composition is investigated with samples of normal water, heavy water, and partially deuterated water. When data below 0.25 GPa are excluded a free volume model fits the observed viscosities well yielding a theoretical glass transition density close to that observed in very-high-density amorphous ice. The improved fit above 0.25 GPa coincides with the loss of other anomalous behaviors in liquid water caused by hydrogen bonding and represents a transition to properties closer to those of a simple liquid.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Sodium doping of solution-processed amine-thiol based CIGS solar cells by thermal evaporation of NaCl

Poor crystallinity, high degree of porosity and rough surfaces are the main drawbacks of solution-processed CIGS absorbers resulting in lower power conversion efficiencies when compared to vacuum-based CIGS solar cells. Therefore, promoting absorber grain growth is key to further improve solution-based solar cell performance. The effect of alkali elements such as Na in CIGS absorbers is generally recognised to have beneficial effects not only on the absorber opto-electronic properties but also on the grain growth. In this work, thermal evaporation of a thin layer of NaCl prior to selenisation resulted in absorbers with significantly larger CIGS grains than previously seen with Na diffusing directly from the from soda-lime glass substrate. NaCl is non-toxic, abundant and readily available compound that has not been typically used as an evaporation source, but rather as an additive into CIGS precursor solution. The effect of Na on these solution-processed CIGS devices was primarily observed in the spectacular morphological changes leading to improved carrier collection and minority carrier lifetimes, but less on the absorber doping. Transmission electron microscopy (TEM) revealed voids forming around large CIGS grains upon NaCl addition and these had a negative effect on inter-grain carrier transport. Nonetheless, the resulting device performance doubled from 5% to 10% with addition of Na using this doping approach; however, a compromise between the optimum grain growth and optimum electronic properties had to be made. This study demonstrates a novel, simple and effective Na-doping strategy for CIGS absorbers and reveals the current limitations of the Na-doping in solution-processed atmospherically deposited cells.

14 SOLAR ENERGY↗
DOE OSTI2020

Electrical control of spin relaxation anisotropy during drift transport in a two-dimensional electron gas

Spin relaxation was studied in a two-dimensional electron gas confined in a wide GaAs quantum well. Recently, the control of the spin relaxation anisotropy by diffusive motion was first shown in D. Iizasa et al., arXiv:2006.08253. Here, we demonstrate electrical control by drift transport in a system with two subbands occupied. The combined effect of in-plane and gate voltages was investigated using time-resolved Kerr rotation. Furthermore, the measured relaxation time presents strong anisotropy with respect to the transport direction. For an in-plane accelerating electric field along [110], the lifetime was strongly suppressed irrespective of the applied gate voltage. Remarkably, for transport along [$1\bara{1}0$], the data shows spin lifetime that was gate dependent and longer than in the [110] direction regardless of the in-plane voltage. In agreement, independent results of anisotropic spin precession frequencies are also presented. Nevertheless, the long spin lifetime, strong anisotropy and drift response seen in the data are beyond the existing models for spin drift and diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Imogolite in water: Simulating the effects of nanotube curvature on structure and dynamics

Imogolite is a fascinating inorganic nanotube that is found in nature or synthesized in a laboratory. The synthesis process is carried out in liquid media, and leads to the formation of almost monodisperse diameter nanotubes. Here we investigate, employing classical molecular dynamics simulations, the interaction of water and imogolite for nanotubes of several radii. We established that water penetrates the pores of N = 9 and larger nanotubes, and adopts a coaxial arrangement in it. Also, while water molecules can diffuse along the center of the nanotube, the molecules next to the inner imogolite walls have very low mobility. At the outer nanotube wall, an increase of water density is observed, this effect extends up to 1 nm, beyond which water properties are bulk-like. Both phenomena are affected by the imogolite curvature.

58 GEOSCIENCES↗
DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
DOE OSTI2021

Deuterium transport and retention properties of representative fusion blanket structural materials

Reduced activation ferritic-martensitic (RAFM) steels have been developed for decades for use as fusion blanket structural materials, and have advantages in both mechanical properties and irradiation resistance following careful engineering of the microstructure. However, the hydrogen isotope behavior in these proposed fusion structural materials is not well understood, but is important to assess since it impacts the fusion reactor safety and self-sufficient tritium fuel cycle. Here, we investigated deuterium transport and retention in representative advanced RAFM steels, including castable nanostructured alloys (CNAs), and oxide-dispersion-strengthened (ODS) steels. A gas-driven permeation (GDP) system was used to measure the permeability, diffusivity and solubility of the studied materials, covering the temperature range from 623 K to 873 K, and the loading pressures from 1.8 x 10 4 to 1.0 x 10 5 Pa. The results indicated that the deuterium permeability has little material dependence. In contrast, the deuterium diffusivity of the studied materials showed significant variation. The deuterium diffusivity in ODS steels is one order of magnitude lower than that in RAFM steels and CNAs, and correspondingly, have an effective solubility that is 2–10 times larger than RAFM steels and CNAs. In addition, thermal desorption spectroscopy (TDS) measurements were performed to assess the deuterium retention and desorption of these materials following a static thermal deuterium charging at 723 Kfor 1 hour under the deuterium pressure of 1.0 x 10 5 Pa. It was found that ODS steels exhibit the highest deuterium retention and have broader desorption peaks. Microstructural features contributing to deuterium retention and impacting deuterium transport are discussed to rationalize the observed deuterium behavior in the studied RAFM steels.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Constitutional under-potential plating (CUP) – New insights for predicting the morphological stability of deposited lithium anodes in lithium metal batteries

Lithium (Li) cycling in lithium metal batteries (LMBs) inevitably witnesses a non-planar Li/electrolyte interface (cellular/dendritic microstructure) affecting the performance, safety and reliability of LMBs. Herein, the “constitutional under-potential plating” (CUP) concept is presented for the first time to understand the evolution, formation and predict the morphological stability/instability of the Li/electrolyte interface. Here, the CUP condition, a measure of the driving force for interface perturbation, occurs when the true potential gradient at the interface is smaller than the equilibrium electrodeposition potential gradient at the Li/electrolyte interface resulting in the electrolyte ahead of the Li-electrolyte interface to be in an under-potential state. Detailed CUP criterion calculations conducted herein establishes the relationship of battery operating conditions (e.g. current density, potential gradient) with the electrolyte properties (e.g. Li-ion transference number and diffusion coefficient) predicting the morphological stability condition at the interface. In addition to the driving force for interface perturbation, the capillarity effect or the Gibbs Thomson effect, a measure of the resistance to perturbation of the perturbed interface is discussed for understanding the influence of interfacial energy/adhesion energy between the current collector and the deposited Li on the overall morphological stability of the Li/electrolyte interface, a critical criterion for the unwanted dendritic morphology.

25 ENERGY STORAGE↗
DOE OSTI2021

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗
DOE OSTI2023

Learning effective SDEs from Brownian dynamic simulations of colloidal particles

We construct a reduced, data-driven, parameter dependent effective stochastic differential equation (eSDE) for electric-field mediated colloidal crystallization using data obtained from Brownian dynamics simulations. We use diffusion maps (a manifold learning algorithm) to identify a set of useful latent observables. In this latent space we identify an eSDE using a deep learning architecture inspired by numerical stochastic integrators and compare it with the traditional Kramers–Moyal expansion estimation. We show that the obtained variables and the learned dynamics accurately encode the physics of the Brownian dynamic simulations. We further illustrate that our reduced model captures the dynamics of corresponding experimental data. Further, our dimension reduction/reduced model identification approach can be easily ported to a broad class of particle systems dynamics experiments/models.

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

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

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

Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++

Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spin Hall effect of vorticity

Using mapping between topological defects in an easy-plane magnet and electrical charges, we study interplay between vorticity and spin currents. Here, we demonstrate that the flow of vorticity is accompanied by the transverse spin current generation; an effect which can be termed as the spin Hall effect of vorticity. We study this effect across the BKT transition and establish the role of dissipation and spin nonconservation in the crossover from spin superfluidity to diffusive spin transport. Our results pave the way for low power computing devices relying on vorticity and spin flows that can propagate over long distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗