Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 649 records · Page 36

Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are impractical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The empirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which satisfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗

Uncertainty in counting ice nucleating particles with continuous flow diffusion chambers

Abstract. This study investigates the measurement of ice nucleating particle (INP) concentrations and sizing of crystals using continuous flow diffusion chambers (CFDCs). CFDCs have been deployed for decades to measure the formation of INPs under controlled humidity and temperature conditions in laboratory studies and by ambient aerosol populations. These measurements have, in turn, been used to construct parameterizations for use in models by relating the formation of ice crystals to state variables such as temperature and humidity as well as aerosol particle properties such as composition and number. We show here that assumptions of ideal instrument behavior are not supported by measurements made with a commercially available CFDC, the SPectrometer for Ice Nucleation (SPIN), and the instrument on which it is based, the Zurich Ice Nucleation Chamber (ZINC). Non-ideal instrument behavior, which is likely inherent to varying degrees in all CFDCs, is caused by exposure of particles to different humidities and/or temperatures than predicated from instrument theory of operation. This can result in a systematic, and variable, underestimation of reported INP concentrations. We find here variable correction factors from 1.5 to 9.5, consistent with previous literature values. We use a machine learning approach to show that non-ideality is most likely due to small-scale flow features where the aerosols are combined with sheath flows. Machine learning is also used to minimize the uncertainty in measured INP concentrations. We suggest that detailed measurement, on an instrument-by-instrument basis, be performed to characterize this uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Sodium Pumping via Condensation within a Non-Wetting Porous Structure

The transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is required to complete the thermodynamic cycle of a sodium thermal electrochemical converter (Na-TEC). A unique capillary pump is explored in this work, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > pi/2) stainless steel porous structure. Due to the curvature of the non-wetting liquid-vapor interface, the liquid adjacent to the interface is at a higher pressure than the vapor. This is in contrast to traditional wicks, where the liquid adjacent to the interface has a lower pressure. Conjugate heat transfer techniques are used to model the coupled momentum and thermal transport processes within this non-wetting porous structure. Due to the low mass flowrates in the Na-TEC (< 0.1 mg/s), higher order friction terms are neglected and the classical Darcy’s law is used for the non-isothermal, compressible vapor flow. At lower temperatures, the sodium vapor density becomes sufficiently small and there is a transition from classic Poiseuille flow (Kn < 0.1) to Knudsen flow (Kn > 10). Furthermore, transport due to binary diffusion is also considered for situations where there is a substantial presence of non-condensable argon in the sodium vapor plenum. Initial modeling results showing expected mass flowrates, temperature profiles, and liquid-vapor interface locations are presented for various physical conditions. The main design features of an experimental set-up used to demonstrate this capillary pumping is also described in detail, including the method used to measure the low sodium mass flowrates. Initial experimental results demonstrating sodium pumping for a range of temperatures (800K – 1100K) will be presented and compared against the model predictions. Finally, several of the challenges with this experiment will be discussed, and specific strategies used to address these are presented.

14 SOLAR ENERGY↗

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Equilibration in fermionic systems

Highlights: •Local equilibration in finite fermion systems is based on a nonlinear diffusion equation. •The fermionic diffusion equation is solved analytically in closed form. •The exact solutions agree with the numerical results. •Particle-number and energy conservation are shown to be fulfilled. •The model is applied to hadrons at low energies, and relativistic quarks. -- Abstract: The time evolution of a finite fermion system towards local statistical equilibrium is investigated using analytical solutions of a nonlinear partial differential equation that had been derived earlier from the Boltzmann collision term. The solutions of this fermionic diffusion equation are rederived in closed form, evaluated exactly for simplified initial conditions, and applied to hadron systems at low energies in the MeV-range, as well as to quark systems at relativistic energies in the TeV-range where antiparticle production is abundant. Conservation laws for particle number including created antiparticles, and for the energy are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

An isogeometric variational multiscale method for large-eddy simulation of coupled multi-ion transport in turbulent flow

Highlights: •We present a computational method for coupled multi-ion transport in turbulent flow. •The underlying formulation is a variational multiscale finite element method. •It is combined with the isogeometric concept for electrochemical systems. •Coupled multi-ion transport in fully turbulent Taylor–Couette flow is simulated. •This example is an important model problem for rotating cylinder electrodes. -- Abstract: Electrochemical processes, such as electroplating of large items in galvanic baths, are often coupled to turbulent flow. In this study, we propose an isogeometric residual-based variational multiscale finite element method for multi-ion transport in dilute electrolyte solutions under turbulent flow conditions. In other words, this means that the concepts of isogeometric discretization and variational multiscale methods are successfully combined for developing a method capable of simulating the challenging problem of coupled multi-ion transport in turbulent flow. We present a comprehensive three-dimensional computational method taking into account, among others, coupled convection–diffusion-migration equations subject to an electroneutrality constraint in combination with phenomenological electrode-kinetics modeling. The electrochemical subproblem is one-way coupled to turbulent incompressible flow via convection. Ionic mass transfer in turbulent Taylor–Couette flow is investigated, representing an important model problem for rotating-cylinder-electrode configurations. Multi-ion transport as considered here is an example for mass transport at high Schmidt number (Sc=1389). An isogeometric discretization is especially advantageous for the present problem, since (i) curved boundaries can be represented exactly, and (ii) it has been proven to provide very accurate solutions for flow quantities when being applied in combination with residual-based variational multiscale modeling. We demonstrate that the method is robust and provides results which are in good agreement with direct numerical simulation results as well as empirical mass-transfer correlations reported in literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive Sampling of Large Deviations

We introduce and test an algorithm that adaptively estimates large deviation functions characterizing the fluctuations of additive functionals of Markov processes in the long-time limit. These functions play an important role for predicting the probability and pathways of rare events in stochastic processes, as well as for understanding the physics of nonequilibrium systems driven in steady states by external forces and reservoirs. The algorithm uses methods from risk-sensitive and feedback control to estimate from a single trajectory a new process, called the driven process, known to be efficient for importance sampling. Its advantages compared to other simulation techniques, such as splitting or cloning, are discussed and illustrated with simple equilibrium and nonequilibrium diffusion models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wall loss of atomic nitrogen determined by ionization threshold mass spectrometry

In the afterglow of an inductively coupled N{sub 2} plasma, relative N atom densities are measured by ionization threshold mass spectrometry as a function of time in order to determine the wall loss time t{sub wN} from the exponential decay curves. The procedure is performed with two mass spectrometers on different positions in the plasma chamber. t{sub wN} is determined for various pressures, i.e., for 3.0, 5.0, 7.5, and 10 Pa. For this conditions also the internal plasma parameters electron density n{sub e} and electron temperature T{sub e} are determined with the Langmuir probe and the rotational temperature T{sub rot}{sup N{sub 2}} of N{sub 2} is determined with the optical emission spectroscopy. For T{sub rot}{sup N{sub 2}}, a procedure is presented to evaluate the spectrum of the transition υ{sup ′}=0→υ{sup ″}=2 of the second positive system (C{sup 3}Π{sub u}→B{sup 3}Π{sub g}) of N{sub 2}. With this method, a gas temperature of 610 K is determined. For both mass spectrometers, an increase of the wall loss times of atomic nitrogen with increasing pressure is observed. The wall loss time measured with the first mass spectrometer in the radial center of the cylindrical plasma vessel increases linearly from 0.31 ms for 3 Pa to 0.82 ms for 10 Pa. The wall loss time measured with the second mass spectrometer (further away from the discharge) is about 4 times higher. A model is applied to describe the measured t{sub wN.} The main loss mechanism of atomic nitrogen for the considered pressure is diffusion to the wall. The surface loss probability β{sub N} of atomic nitrogen on stainless steel was derived from t{sub wN} and is found to be 1 for the present conditions. The difference in wall loss times measured with the mass spectrometers on different positions in the plasma chamber is attributed to the different diffusion lengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Leaching behavior of copper from waste printed circuit boards with Brønsted acidic ionic liquid

Highlights: • A Brønsted acidic ILs was used to leach Cu from WPCBs for the first time. • The particle size of WPCBs has significant influence on Cu leaching rate. • Cu leaching rate was higher than 99% under the optimum leaching conditions. • The leaching process can be modeled with shrinking core model, and the E{sub a} was 25.36 kJ/mol. - Abstract: In this work, a Brønsted acidic ionic liquid, 1-butyl-3-methyl-imidazolium hydrogen sulfate ([bmim]HSO{sub 4}), was used to leach copper from waste printed circuit boards (WPCBs, mounted with electronic components) for the first time, and the leaching behavior of copper was discussed in detail. The results showed that after the pre-treatment, the metal distributions were different with the particle size: Cu, Zn and Al increased with the increasing particle size; while Ni, Sn and Pb were in the contrary. And the particle size has significant influence on copper leaching rate. Copper leaching rate was higher than 99%, almost 100%, when 1 g WPCBs powder was leached under the optimum conditions: particle size of 0.1–0.25 mm, 25 mL 80% (v/v) ionic liquid, 10 mL 30% hydrogen peroxide, solid/liquid ratio of 1/25, 70 °C and 2 h. Copper leaching by [bmim]HSO{sub 4} can be modeled with the shrinking core model, controlled by diffusion through a solid product layer, and the kinetic apparent activation energy has been calculated to be 25.36 kJ/mol.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Factors controlling surface oxygen exchange in oxides

Decreasing the working temperature of solid oxide fuel cells is critical to their increased commercialization but is inhibited by the slow oxygen exchange kinetics at the cathode, which limits the overall rate of the oxygen reduction reaction. We use ab initio methods to develop a quantitative elementary reaction model of oxygen exchange in a representative cathode material, La 0.5 Sr 0.5 CoO 3-δ , and predict that under operating conditions the rate-limiting step for oxygen incorporation from O 2 gas on the stable, (001)-SrO surface is lateral (surface) diffusion of O-adatoms and oxygen surface vacancies. We predict that a high vacancy concentration on the metastable CoO 2 termination enables a vacancy-assisted O-2 dissociation that is 10 2 -10 3 times faster than the rate limiting step on the Sr-rich (La,Sr)O termination. This result implies that dramatically enhanced oxygen exchange performance could potentially be obtained by suppressing the (La,Sr)O termination and stabilizing highly active CoO 2 termination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stratification of drying particle suspensions: Comparison of implicit and explicit solvent simulations

We report that large scale molecular dynamics simulations are used to study drying suspensions of a binary mixture of large and small particles in explicit and implicit solvents. The solvent is first modeled explicitly and then mapped to a uniform viscous medium by matching the diffusion coefficients and the pair correlation functions of the particles. “Small-on-top” stratification of the particles, with an enrichment of the smaller ones at the receding liquid-vapor interface during drying, is observed in both models under the same drying conditions. With the implicit solvent model, we are able to model much thicker films and study the effect of the initial film thickness on the final distribution of particles in the dry film. Laslty, our results show that the degree of stratification is controlled by the Péclet number defined using the initial film thickness as the characteristic length scale. When the Péclet numbers of large and small particles are much larger than 1, the degree of “small-on-top” stratification is first enhanced and then weakens as the Péclet numbers are increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗

Empirical and mechanistic transient fission gas release model for high-burnup LOCA conditions

In response to the nuclear industry desire to extend burnup beyond current licensing practices, the US Nuclear Regulatory Commission (NRC) released a research information letter (RIL) that provides a basis for analyzing fuel fragmentation, relocation, and dispersal (FFRD) in light-water reactors. Of the five elements discussed, the most ambiguous is the significance of transient fission gas (FG) release (FGR) (tFGR) and its effects on fuel performance under loss-of-coolant accident conditions. In fresh fuel, FG migration and eventual release is primarily governed by diffusion-based mechanisms at higher temperatures (>1,000 °C). However, the mechanisms governing FGR changes as burnup increases. More recent research indicates that FGR increases as burnup increases, specifically under temperature transient conditions, and this release occurs at lower temperatures with a new release mechanism. This behavior has been attributed to microcracking in the fuel and is likely related to microstructure embrittlement with the presence of over pressurized FG bubbles. The NRC RIL outlines the complexity of the phenomenon and a need for a deeper understanding to adequately address FFRD for regulatory application. Here, therefore, this manuscript intends to summarize the publicly available tFGR data and discuss the observed dependencies (e.g., burnup, heating rate, sample geometry, terminal temperature). An empirical model has been developed and benchmarked against recently published experimental data. However, this empirical model is limited to conditions for which fitting data exist and, therefore, a high-level discussion is included to provide a roadmap for atomistically-informed multiscale modeling in conjunction with experimental data collection to develop a mechanistic tFGR model widely applicable to a broad range of nuclear fuel conditions at high burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling CO 2 flow in support of a shallow subsurface controlled leakage field test

Here, controlled release of carbon dioxide (CO 2 ) into the soil and atmosphere is performed to test detection and monitoring tools, for which several field laboratories were established by a number of institutions worldwide. Numerical simulations of CO 2 behavior in the shallow subsurface region are other forms of validation and verification of the leakage pathways and destinations. These reports aim to improve monitoring and verification of CO 2 in case of unexpected leakages for public assurance. In this work, we introduce the results of a numerical modeling study conducted to simulate the injection of CO 2 as carried out during a field test in Viamão, southern Brazil, where 20 kg day –1 of CO 2 was pumped for 30 days through a vertical well 3 m below ground in an altered granitic soil. Multiphase flow simulations were performed with the TOUGH2/EOS7CA software for unsaturated porous media, using field data and injection parameters, including sensitivity tests to permeability direction, diffusivity, and boundary conditions. Results with increased horizontal permeabilities are in better agreement with the field observations. In this condition, mass balance calculations indicate approximately 90% of injected CO2 (20 kg day –1 during 30 days) remains in the soil after 180 days from injection start, consistent with the measured flow through the soil–atmosphere interface.

03 NATURAL GAS↗

Comparative study of solute trapping and Gibbs free energy changes at the phase interface during alloy solidification under local nonequilibrium conditions

An analytical model has been developed to describe the influence of solute trapping during rapid alloy solidification on the components of the Gibbs free energy change at the phase interface with emphasis on the solute drag energy. For relatively low interface velocity V < V{sub D}, where V{sub D} is the characteristic diffusion velocity, all the components, namely mixing part, local nonequilibrium part, and solute drag, significantly depend on solute diffusion and partitioning. When V ≥ V{sub D}, the local nonequilibrium effects lead to a sharp transition to diffusionless solidification. The transition is accompanied by complete solute trapping and vanishing solute drag energy, i.e. partitionless and “dragless” solidification.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Efficacy of Repeat Voltage Holds for Conditioning and Calendar Life Testing of Graphite and Silicon Cells

Voltage-hold (V-hold) protocols have shown promise toward calendar lifetime analysis of cells with graphite (Gr) and silicon (Si) anodes. In this work, repeat V-holds are performed on Gr and Si cells paired with lithium iron phosphate cathodes to delineate their beneficial role in formation and conditioning. We find that V-hold at the top of charge supplements constant current cycling in conditioning the cell to higher capacities for both Gr and Si cells after the first V-hold. A reduced order model provides the irreversible capacity proportions of each V-hold. With each repeat V-hold, parasitic loss of lithium to the solid electrolyte interphase (SEI) decreases on both Gr and Si cells. Gr cells show the square-root-of-time capacity loss behavior within 200 h of V-hold, indicative of its fast relaxation and low impact of reference performance test cycles on the SEI growth. Lifetime estimates from repeat V-holds on Gr can reach years. Si exhibits longer transition times from kinetic to diffusion-limited SEI growth, evidenced by the 400 h and 200 h holds showing square-root-of-time and linear behavior, respectively. Lifetime predictions from repeat V-holds on Si only reach 1–2 months, highlighting its limitations. Recommended duration of V-holds for Si cells should be ≥400 h.

25 ENERGY STORAGE↗