SEARCH · Engineering Papers

Results for “diffusion processes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85

All sources NASA NTRS DOE OSTI

DOE OSTI2013

Characterization of statistical prior image constrained compressed sensing (PICCS): II. Application to dose reduction

Purpose: The ionizing radiation imparted to patients during computed tomography exams is raising concerns. This paper studies the performance of a scheme called dose reduction using prior image constrained compressed sensing (DR-PICCS). The purpose of this study is to characterize the effects of a statistical model of x-ray detection in the DR-PICCS framework and its impact on spatial resolution. Methods: Both numerical simulations with known ground truth and in vivo animal dataset were used in this study. In numerical simulations, a phantom was simulated with Poisson noise and with varying levels of eccentricity. Both the conventional filtered backprojection (FBP) and the PICCS algorithms were used to reconstruct images. In PICCS reconstructions, the prior image was generated using two different denoising methods: a simple Gaussian blur and a more advanced diffusion filter. Due to the lack of shift-invariance in nonlinear image reconstruction such as the one studied in this paper, the concept of local spatial resolution was used to study the sharpness of a reconstructed image. Specifically, a directional metric of image sharpness, the so-called pseudopoint spread function (pseudo-PSF), was employed to investigate local spatial resolution. Results: In the numerical studies, the pseudo-PSF was reduced from twice the voxel width in the prior image down to less than 1.1 times the voxel width in DR-PICCS reconstructions when the statistical model was not included. At the same noise level, when statistical weighting was used, the pseudo-PSF width in DR-PICCS reconstructed images varied between 1.5 and 0.75 times the voxel width depending on the direction along which it was measured. However, this anisotropy was largely eliminated when the prior image was generated using diffusion filtering; the pseudo-PSF width was reduced to below one voxel width in that case. In the in vivo study, a fourfold improvement in CNR was achieved while qualitatively maintaining sharpness; images also had a qualitatively more uniform noise spatial distribution when including a statistical model. Conclusions: DR-PICCS enables to reconstruct CT images with lower noise than FBP and the loss of spatial resolution can be mitigated to a large extent. The introduction of statistical modeling in DR-PICCS may improve some noise characteristics, but it also leads to anisotropic spatial resolution properties. A denoising method, such as the directional diffusion filtering, has been demonstrated to reduce anisotropy in spatial resolution effectively when it was combined with DR-PICCS with statistical modeling.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Thin film studies toward improving the performance of accelerator electron sources

Future electron accelerators require DC high voltage photoguns to operate beyond the present state of the art to conduct new experiments that require ultra-bright electron beams with high average current and higher bunch charge. To meet these demands, the accelerators must demonstrate improvements in a number of photogun areas including vacuum, field emission elimination in high voltage electrodes, and photocathodes. This dissertation illustrates how these improvements can be achieved by the application of suitable thin-films to the photogun structure for producing ultra-bright electron beams. This work is composed of three complementary studies. First, the outgassing rates of three nominally identical 304L stainless steel vacuum chambers were studied to determine the effects of chamber coatings (silicon and titanium nitride) and heat treatments. For an uncoated stainless steel chamber, the diffusion limited outgassing was taken over by the recombination limited process as soon as a low outgassing rate of ~1.79(±0.05) x 10- 13 Torr L s -1 cm -2 was achieved. An amorphous silicon coating on the stainless steel chambers exhibited recombination limited behavior and any heat treatment became ineffective in reducing the outgassing rate. A TiN coated chamber yielded the smallest apparent outgassing rate of all the chambers: 6.44(±0.05) x 10 -13 Torr L s -1 cm -2 following an initial 90 °C bake and 2(±20) x 10 -16 Torr L s -1 cm -2 following the final bake in the series. This perceived low outgassing rate was attributed to the small pumping nature of TiN coating itself. Second, the high voltage performance of three TiN-coated aluminum electrodes, before and after gas conditioning with helium, were compared to that of bare aluminum electrodes and electrodes manufactured from titanium alloy (Ti-6Al-4V). This study suggests that aluminum electrodes, coated with TiN, could simplify the task of implementing photocathode cooling, which is required for future high current electron beam applications. The best performing TiN-coated aluminum electrode demonstrated less than 15 pA of field emission current at -- 175 kV for a 10 mm cathode/anode gap, which corresponds to a field strength of 22.5 MV/m. Third, the effect of antimony thickness on the performance of bialkali-antimonide photocathodes was studied. The high-capacity effusion source enabled us to successfully manufacture photocathodes having a maximum QE around 10% and extended low voltage 1/e lifetime (> 90 days) at 532 nm via the co-deposition method, with relatively thick layers of antimony (≥ 300 nm). We speculate that alkali co-deposition provides optimized stoichiometry for photocathodes manufactured using thick Sb layers, which could serve as a reservoir for the alkali. In summary, this research examined the effectiveness of thin films applied on photogun chamber components to achieve an extremely high vacuum, to eliminate high voltage induced field emission from electrodes, and to generate photocurrent with high quantum yield with an extended operational lifetime. Simultaneous implementation of these findings can meet the challenges of future ultra-bright photoguns.

42 ENGINEERING↗

DOE OSTI2019

Microstructural evolution of thin polycrystalline metallic films under extreme conditions

Thin films are materials systems that are widely used in applications ranging from electronics and optical devices to industrial and biomedical ones. However, these systems are unstable against various homogenization processes and aging mechanisms (grain growth, coarsening, surface evolution, and diffusion of species) even at low service temperatures. In this work, we examine the role of various aspects of microstructure (grain boundary types and characters, free surfaces, surface diffusion, and thermal grooves) on the thermal aging of such systems. Existing experimental tools will be leveraged to characterize thin films, i.e., in-situ quantitative thermal annealing, and provide direct comparisons to predications emerging from a recently developed meso-scale model via Precession Electron Diffraction (PED). Parametric studies will be conducted to gain insights on the role of each of the aforementioned aspects of microstructure on the dynamics. Herein, we will focus on “hard” gold thin films (with alloying elements like Co, Ni, or Fe) as they are materials of choice in a wide range of applications at Sandia. Initially, pure gold systems are examined and future studies will focus on microstructure evolution in the presence of alloying elements and second-phase particles.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Integration of Pretreatment With Simultaneous Counter-Current Extraction of Energy Sorghum for High-Titer Mixed Sugar Production

Sorghum (Sorghum bicolor L. Moench) offers substantial potential as a feedstock for the production of sugar-derived biofuels and biochemical products from cell wall polysaccharides (i. e., cellulose and hemicelluloses) and water-extractable sugars (i.e., glucose, fructose, sucrose, and starch). A number of preprocessing schemes can be envisioned that involve processes such as sugar extraction, pretreatment, and densification that could be employed in decentralized, regional-scale biomass processing depots. In this work, an energy sorghum exhibiting a combination of high biomass productivity and high sugar accumulation was evaluated for its potential for integration into several potential biomass preprocessing schemes. This included counter-current extraction of water-soluble sugars followed by mild NaOH or liquid hot water pretreatment of the extracted bagasse. A novel processing scheme was investigated that could integrate with current diffuser-type extraction systems for sugar extraction. In this approach, mild NaOH pretreatment (i.e., <90°C) was performed as a counter-current extraction to yield both an extracted, pretreated bagasse and a high-concentration mixed sugar stream. Following hydrolysis of the bagasse, the combined hydrolysates derived from cellulosic sugars and extractable sugars were demonstrated to be fermentable to high ethanol titers (>8%) at high metabolic yields without detoxification using a Saccharomyces cerevisiae strain metabolically engineered and evolved to ferment xylose.

09 BIOMASS FUELS↗

DOE OSTI2023

Design and Characterization of a Burst Mode 20 Mfps Low Noise CMOS Image Sensor

This paper presents a novel ultra-high speed, high conversion-gain, low noise CMOS image sensor (CIS) based on charge-sweep transfer gates implemented in a standard 180 nm CIS process. Through the optimization of the photodiode geometry and the utilization of charge-sweep transfer gates, the proposed pixels achieve a charge transfer time of less than 10 ns without requiring any process modifications. Moreover, the gate structure significantly reduces the floating diffusion capacitance, resulting in an increased conversion gain of 183 µV/e−. This advancement enables the image sensor to achieve the lowest reported noise of 5.1 e− rms. To demonstrate the effectiveness of both optimizations, a proof-of-concept CMOS image sensor is designed, taped-out and characterized.

47 OTHER INSTRUMENTATION↗

DOE OSTI2016

Magnetically self-regulated formation of early protoplanetary disks

The formation of protoplanetary disks during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disk radius, r≃18 au (η{sub AD}/0.1s){sup 2/9}(B{sub z}/0.1G){sup −4/9}(M/0.1M{sub ⊙}){sup 1/3}, where M is the total disk plus protostar mass, η{sub AD} is the ambipolar diffusion coefficient, and B {sub z} is the magnetic field in the inner part of the core. This is significantly smaller than the disks that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disk radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disk radii upon the various relevant quantities, suggesting weak variations among class-0 disk sizes. In some cases, we note the onset of large spiral arms beyond this radius.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Distribution of Neutrons in a Cold Lucite Moderator

A series of experiments was performed in 2010 to study the behavior of neutrons in a cold Lucite moderator. The diffusion of neutrons in a cold moderator is described as a two-step process. In the first step, the fast neutrons produced by a neutron source collide with the nuclei of Lucite, losing energy until they reach thermal energies. In the second step, the thermal neutrons continue to diffuse without any significant loss of energy until they are absorbed or leak out of the system. However, because the Lucite moderator is at low temperature, there is less up scattering in the neutron energy compared to a Lucite moderator at room temperature. Experiments were performed to study the behavior of neutrons and their distribution in a cold Lucite moderator and their results compared to simulations. Based on these measurements, an estimation of the absorption cross section at -78.5 °C yielded a value of 0.54 barns ± 0.02 barns.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Simulating the Effect of Fracture Connectivity on Repository Performance with GDSA Framework - 18589

This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2014

Evaluation of Low Power ARM Processors for Monte Carlo Particle Transport: MCNP6 on the ARM Cortex-A8 - Paper 83

In the push towards exascale computing, the limiting factors are power requirements and heat dissipation. Using current technology, the power required to drive an exascale supercomputer would near 400-MW. Clearly, future computing hardware will have to be more energy efficient in order to achieve exascale performance. Most of the current research in advanced computing architectures is in the use of standard CPUs with accelerators. However, these accelerators are not generally supported by the traditional scientific programming languages (C, C++, and Fortran) and the scientific software must be tailored specifically to take advantage of the accelerators. Large scale scientific software (hundreds-of-thousand to millions of lines of code) can not be easily ported to the computing platform du jour. The widespread use of low power ARM based processors in the mobile computing market offer a low power alternative to accelerators that are also supported by the traditional scientific programming languages and programming paradigms such as MPI and OpenMP. Computing architectures for scientific computing are often assessed using standard benchmarks such as LINPACK benchmark or using suites of mini-applications that are greatly simplified stand-ins for scientific software. Often the real applications that will deployed on the scientific supercomputers are not ported or tested until delivery of the machines. In order to evaluate the use of ARM processors for scientific computing the Monte Carlo particle transport code, MCNP version 6 was compiled on the ARM Cortex-A8 and a set of benchmark problems were used to assess performance in comparison to an Intel Xeon E5-2670 and the low power Intel Atom N270. MCNP is a continuous energy Monte Carlo particle transport code in development at Los Alamos National Laboratory since 1977. MCNP is used to transport neutrons, photon, electrons, and charged particles using continuous energy interaction data. It is commonly used in fields of nuclear power, nuclear medicine, radiation protection, and high-energy physics through-out the world. MCNP version 6 is approximately 450 thousand lines of code, written primary in Fortran90, with some C and C++. MCNP supports distributed multiprocessing using MPI and threading with OpenMP. The ARM Cortex-A8 executes MCNP6 using about twice as much power as the Intel Xeon E5-2670 and about the same amount of power as the Intel Atom N270. However, these results are encouraging for the use of low-power ARM chips in future supercomputers. The ARM Cortex-A8 is a 32-bit CPU and it's NEON general-purpose SIMD (single-instruction multiple data) engine does not support double-precision floating-point arithmetic. More recent ARM processors (ARM Cortex-A53 and Cortex-A57) support 64-bit addressing and have NEON co-processors that support double-precision floating-point arithmetic. ARM claims that the A50 series delivers a performance increase of 15 times the Cortex-A8 performance for similar power usage. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2013

Stochastic Maximum Principle for Optimal Control of SPDEs

We prove a version of the maximum principle, in the sense of Pontryagin, for the optimal control of a stochastic partial differential equation driven by a finite dimensional Wiener process. The equation is formulated in a semi-abstract form that allows direct applications to a large class of controlled stochastic parabolic equations. We allow for a diffusion coefficient dependent on the control parameter, and the space of control actions is general, so that in particular we need to introduce two adjoint processes. The second adjoint process takes values in a suitable space of operators on L{sup 4}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Extracting the diffusion tensor from molecular dynamics simulation with Milestoning

We propose an algorithm to extract the diffusion tensor from Molecular Dynamics simulations with Milestoning. A Kramers-Moyal expansion of a discrete master equation, which is the Markovian limit of the Milestoning theory, determines the diffusion tensor. To test the algorithm, we analyze overdamped Langevin trajectories and recover a multidimensional Fokker-Planck equation. The recovery process determines the flux through a mesh and estimates local kinetic parameters. Rate coefficients are converted to the derivatives of the potential of mean force and to coordinate dependent diffusion tensor. We illustrate the computation on simple models and on an atomically detailed system—the diffusion along the backbone torsions of a solvated alanine dipeptide.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Fluctuating hydrodynamics of reactive liquid mixtures

Fluctuating hydrodynamics (FHD) provides a framework for modeling microscopic fluctuations in a manner consistent with statistical mechanics and nonequilibrium thermodynamics. This paper presents an FHD formulation for isothermal reactive incompressible liquid mixtures with stochastic chemistry. Fluctuating multispecies mass diffusion is formulated using a Maxwell-Stefan description without assuming a dilute solution, and momentum dynamics is described by a stochastic Navier-Stokes equation for the fluid velocity. We consider a thermodynamically consistent generalization for the law of mass action for non-dilute mixtures and use it in the chemical master equation (CME) to model reactions as a Poisson process. The FHD approach provides remarkable computational efficiency over traditional reaction-diffusion master equation methods when the number of reactive molecules is large, while also retaining accuracy even when there are as few as ten reactive molecules per hydrodynamic cell. We present a numerical algorithm to solve the coupled FHD and CME equations and validate it on both equilibrium and nonequilibrium problems. We simulate a diffusively driven gravitational instability in the presence of an acid-base neutralization reaction, starting from a perfectly flat interface. We demonstrate that the coupling between velocity and concentration fluctuations dominates the initial growth of the instability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Amorphization mechanism of SrIrO 3 electrocatalyst: How oxygen redox initiates ionic diffusion and structural reorganization

The use of renewable electricity to prepare materials and fuels from abundant molecules offers a tantalizing opportunity to address concerns over energy and materials sustainability. The oxygen evolution reaction (OER) is integral to nearly all material and fuel electrosyntheses. However, very little is known about the structural evolution of the OER electrocatalyst, especially the amorphous layer that forms from the crystalline structure. Here, we investigate the interfacial transformation of the SrIrO 3 OER electrocatalyst. The SrIrO 3 amorphization is initiated by the lattice oxygen redox, a step that allows Sr 2+ to diffuse and O 2− to reorganize the SrIrO 3 structure. This activation turns SrIrO 3 into a highly disordered Ir octahedral network with Ir square-planar motif. The final Sr y IrO x exhibits a greater degree of disorder than IrO x made from other processing methods. Our results demonstrate that the structural reorganization facilitated by coupled ionic diffusions is essential to the disordered structure of the SrIrO 3 electrocatalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Precipitation of calcium carbonate and calcium phosphate under diffusion controlled mixing

Multi-component mineral precipitation in porous, subsurface environments is challenging to simulate or engineer when in situ reactant mixing is controlled by diffusion. In contrast to well-mixed systems, the conditions that favor mineral precipitation in porous media are distributed along chemical gradients, which evolve spatially due to concurrent mineral precipitation and modification of solute transport in the media. The resulting physical and chemical characteristics of a mixing/precipitation zone are a consequence of coupling between transport and chemical processes, and the distinctive properties of individual chemical systems. We examined the spatial distribution of precipitates formed in “double diffusion” columns for two chemical systems, calcium carbonate and calcium phosphate. Polyacrylamide hydrogel was used as a low permeability, high porosity medium to maximize diffusive mixing and minimize pressure- and density-driven flow between reactant solutions. In the calcium phosphate system, multiple, visually dense and narrow bands of precipitates were observed that were reminiscent of previously reported Liesegang patterns. In the calcium carbonate system, wider precipitation zones characterized by more sparse distributions of precipitates and a more open channel structure were observed. In both cases, formation of precipitates inhibited, but did not necessarily eliminate, continued transport and mixing of the reactants. A reactive transport model with fully implicit coupling between diffusion, chemical speciation and precipitation kinetics, but where explicit details of nucleation processes were neglected, was able to qualitatively simulate properties of the precipitation zones. The results help to illustrate how changes in the physical properties of a precipitation zone depend on coupling between diffusion-controlled reactant mixing and chemistry-specific details of precipitation kinetics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

An analysis of hydrated proton diffusion in ab initio molecular dynamics

A detailed understanding of the inherently multiscale proton transport process raises a number of scientifically challenging questions. For example, there remain many (partially addressed) questions on the molecular mechanism for long-range proton migration and the potential for the formation of long-lived traps giving rise to burst-and-rest proton dynamics. Using results from a sizeable collection of ab initio molecular dynamics (AIMD) simulations (totaling ∼2.7 ns) with various density functional approximations (Becke-Lee-Yang-Parr (BLYP), BLYP–D3, Hamprecht-Cohen-Tozer-Handy, B3LYP) and temperatures (300–330 K), equilibrium and dynamical properties of one excess proton and 128 water molecules are studied. Two features in particular (concerted hops and weak hydrogen-bond donors) are investigated to identify modes in the system that are strongly correlated with the onset of periods of burst-and-rest dynamics. The question of concerted hops seeks to identify those time scales over which long-range proton transport can be classified as a series of sequential water hopping events or as a near-simultaneous concerted process along compressed water wires. The coupling of the observed burst-and-rest dynamics with motions of a fourth neighboring water molecule (a weak hydrogen-bond donor) solvating the protonated water molecule is also investigated. The presence (absence) of hydrogen bonds involving this fourth water molecule before and after successful proton hopping events is found to be strongly correlated with periods of burst (rest) dynamics (and consistent with pre-solvation concepts). By analyzing several realizations of the AIMD trajectories on the 100-ps time scale, convergence of statistics can be assessed. For instance, it was observed that the probability for a fourth water molecule to approach the hydronium, if not already proximal at the beginning of the lifetime of the hydronium, is very low, indicative of the formation of stable void regions. Furthermore, the correlations of the neighboring water atoms are identified as the fourth water approaches the hydronium. Finally, the temperature effects on structural and dynamical properties are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

A review of rare earth elements and yttrium in coal ash: Content, modes of occurrences, combustion behavior, and extraction methods

Rare earth elements and yttrium (REY) have attracted considerable attention over the last decade because of their vital roles in clean energy, consumer product, national defense and security applications, among other uses. Due to the retention of REY during coal burning, coal combustion ash is considered as potential alternative sources for REY. Understanding the content, speciation, retention and/or transformation behavior of REY during coal combustion not only expands our knowledge of the combustion behavior of the trace elements in coal, but also provides basis for modeling REY partitioning during coal combustion and for developing economically viable REY recovery technologies. This review makes a critical summary of recent progress in the study of REY in coal ash. The contents and the extraction potentials of REY in coal ash derived from 15 major coal-producing countries worldwide were summarized and evaluated. Various analytical methods for determining REY bulk contents and speciation, together with the solid sample pretreatment, analytical accuracy and precision, advantages and disadvantages were summarized and compared. Modern analytical approaches combined indirect methods (e.g., sequential extraction) shed light on the physical distribution, mineralogy, and the chemical state of REY in coal ash. Three types of REY occurrences in coal ash, including Si-Al glassy association, discrete minerals or compounds, and organic association (bound with unburned carbon) were defined in the review. The glassy association can be further divided into REY minerals closely bound to glass phases and dispersed throughout the glassy structure. REY partitioning in various emission streams, the size distribution, and their enrichment behavior in coal ash were discussed. Additionally, thermal behavior and transformation of various REY forms in coal during combustion process, including organic-associated REY, REY phosphates, REY carbonates, clay-bound REY and among others were summarized. Two possible retention mechanisms of REY by aluminosilicate glass at boiler temperature were proposed: the incorporation of the individual REY phases into the glass as inclusions and the diffusion of REY phases throughout Si-Al glass structures in the melting process. Feed coal mineral types, mineral-mineral associations, boiler conditions, and other factors control the retention process. After coal combustion, the speciation of REY in fly ash may be modified by the reactions of REY phases with flue gas components. Further, an overview of REY transformation mechanisms during coal combustion was deeply discussed. Finally, current extraction techniques for REY recovery from coal combustion ash were introduced. Future outlooks and research problems were also identified.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Verification and validation of a rapid heat transfer calculation methodology for transient melt pool solidification conditions in powder bed metal additive manufacturing

A fundamental understanding of spatial and temporal thermal distributions is crucial for predicting solidification and solid-state microstructural development in parts made by additive manufacturing. While sophisticated numerical techniques that are based on finite element or finite volume methods are useful for gaining insight into these phenomena at the length scale of the melt pool (100 - 500 µm), they are ill-suited for predicting engineering trends over full part cross-sections (> 10 x 10 cm) or many layers over long process times (> many days) due to the necessity of fully resolving the heat source characteristics. On the other hand, it is extremely difficult to resolve the highly dynamic nature of the process using purely in-situ characterization techniques. This article proposes a pragmatic alternative based on a semi-analytical approach to predicting the transient heat conduction during powder bed metal additive manufacturing process. The model calculations were theoretically verified for selective laser melting of AlSi10Mg and electron beam melting of IN718 powders for simple cross-sectional geometries and the transient results are compared to steady state predictions from the Rosenthal equation. It is shown that the transient effects of the scan strategy create significant variations in the melt pool geometry and solid-liquid interface velocity, especially as the thermal diffusivity of the material decreases and the pre-heat of the process increases. With positive verification of the strategy, the model was then experimentally validated to simulate two point-melt scan strategies during electron beam melting of IN718, one intended to produce a columnar and one an equiaxed grain structure. Lastly, through comparison of the solidification conditions (i.e. transient and spatial variations of thermal gradient and liquid-solid interface velocity) predicted by the model to phenomenological CET theory, the model accurately predicted the experimental grain structures.

36 MATERIALS SCIENCE↗