Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 631 records · Page 35

Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination

Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗

Vapor–Fed Electrolyzers for Carbon Dioxide Reduction Using Tandem Electrocatalysts: Cuprous Oxide Coupled with Nickel–Coordinated Nitrogen–Doped Carbon

We report ethylene is particularly attractive due to its major importance as a feedstock for various applications including the polymer industry. As such, catalyst and electrolyzer developments are crucial to achieve industrially relevant ethylene production and efficiency levels. In this work, a tandem electrocatalyst composed of copper nanocubes and nickel-coordinated nitrogen-doped carbon (NiNC) is presented, which is integrated into gas diffusion electrodes (GDEs) for direct conversion of vapor-fed CO 2 into ethylene. Evaluation of tandem GDEs in the vapor-fed flow electrolyzer shows significantly increased ethylene selectivity in terms of faradaic efficiency and C 2 H 4 /CO ratio compared to a non-tandem copper GDE. The enhancements are attributed to the increased local CO availability near the copper surface via effective CO 2 to CO conversion on neighboring NiNC. The experimental results are validated by 3D resolved continuum simulations, which show increased flux of higher-order prod-ucts with the added CO flux from NiNC. The practical viability of Cu/NiNC catalyst is further evaluated in a membrane electrode assembly electrolyzer, achieving 40% FE toward ethylene at 150 mA cm -2 and 3.2 V. These find-ings highlight the high selectivity and formation rate of ethylene achieved by successful device integration of the Cu/NiNC catalyst, demonstrating the potential for implementation in large-scale sustainable CO 2 electrolyzers.

36 MATERIALS SCIENCE↗

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗

Analytical-based simulation approach for an anion exchange membrane fuel cell

An analytical and empirical-based 1-D, non-isothermal, steady-state model for anion exchange membrane fuel cell capable of capturing two-phase phenomena is presented in this study. Coupled multi-physics including mass and charge transport, electrochemical reactions, heat transfer, and two-phase water transport are considered in the model and the simulated results are compared to experimental data. To better represent actual material properties and localized conditions, the model applies multilayer discretization in the gas diffusion electrode to enhance prediction accuracy. The model successfully predicts the baseline performance at 70 °C, 131 kPa abs., 92% RH with pure H 2 /O 2 gas as well as the limiting current at 10% H 2 . The robust simulation approach allows for simplistic and accurate estimation of cell performance without the complications of applying two-phase parameters and expensive computational need for numerical models. In addition, the results from the sensitivity studies of material properties and operating conditions provide valuable insights on water management strategies and optimal component design for advancing anion exchange membrane fuel cell technology.

1-D model↗

Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data

Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.

MPPD↗

A chemo-mechanical model for describing sorption hysteresis in a glassy polyurethane

Abstract Hysteretic sorption and desorption of water is observed from 0 to 95% relative humidity and 298–333 K on a glassy polyurethane foam. It is postulated that sorption-induced swelling of the glassy polyurethane increases the concentration of accessible hydrogen-bonding adsorption sites for water. The accessibility of sites is kinetically controlled due to the restricted thermal motions of chains in the glassy polymer, causing a difference in accessible site concentrations during sorption and desorption. This discrepancy leads to hysteresis in the sorbed concentrations of water. A coupled chemo-mechanical model relating volumetric strain, adsorption site concentration, and sorbed water concentration is employed to describe water sorption hysteresis in the glassy polyurethane. This model not only describes the final mass uptake for each relative humidity step, but also captures the dynamics of water uptake, which exhibit diffusion and relaxation rate-controlled regimes.

36 MATERIALS SCIENCE↗

Ultraviolet spectroscopy of pressurized and supercritical carbon dioxide

Abstract Carbon dioxide (CO 2 ) is prevalent in planetary atmospheres and sees use in a variety of industrial applications. Despite its ubiquitous nature, its photochemistry remains poorly understood. In this work we explore the density dependence of pressurized and supercritical CO 2 electronic absorption spectra by vacuum ultraviolet spectroscopy over the wavelength range 1455-2000 Å. We show that the lowest absorption band transition energy is unaffected by a density increase up to and beyond the thermodynamic critical point (137 bar, 308 K). However, the diffuse vibrational structure inherent to the spectrum gradually decreases in magnitude. This effect cannot be explained solely by collisional broadening and/or dimerization. We suggest that at high densities close proximity of neighboring CO 2 molecules with a variety of orientations perturbs the multiple monomer electronic state potential energy surfaces, facilitating coupling between binding and dissociative states. We estimate a critical radius of ~4.1 Å necessary to cause such perturbations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weak trimerization in the frustrated two-dimensional triangular Heisenberg antiferromagnet Lu y Y 1-y MnO 3

To understand the 2D triangular Heisenberg antiferromagnetic system, we investigated the magnetic structures and the dynamics of Lu y Y 1-y MnO 3 in detail. The substitutions are adjusted to the Mn atomic position close to x Mn =$\frac{1}{3}$. The neutron powder diffraction data claims that the magnetic structure of Lu y Y 1-y MnO 3 is described as a mixture of Γ 3 (P6' 3 cm') and Γ 4 (P6' 3 c'm) at the x Mn position for y=0.15, 0.30, and 0.45. The ratio of Γ 3 and Γ 4 depends on temperature and composition and the fraction of Γ 3 increases upon cooling, while no clear trimerization was observed at the xMn position. We estimated exchange parameters from the analysis of the low-energy part of the spin waves. The results showed a weak trimerization effect on cooling because the nearest-neighbor exchange interaction is slightly enhanced. The temperature dependence of the spin-wave dispersion around the Γ point shows that the spin gap closes with increasing temperature because the exchange interactions in the nearest Mn-Mn neighbor become smaller. Gapless diffusive magnetic excitation from a Mn triangular lattice has been observed in a wide range in Q and E space of Lu y Y 1-y MnO 3 . We found that Lu 0.3 Y 0.7 MnO 3 could be an ideal case to investigate the trimerization, frustrated magnetism, and magnetoelastic coupling often observed in two-dimensional triangular lattice Heisenberg antiferromagnet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Importance of the Antarctic Slope Current in the Southern Ocean Response to Ice Sheet Melt and Wind Stress Change

Abstract We use two coupled climate models, GFDL‐CM4 and GFDL‐ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre‐industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near‐surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large‐scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean‐state representation and meltwater‐driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.

Beadling, R. L.↗

Non-galvanic mass transport in molten fluoride salt isothermal corrosion cells

Mass transport of corrosion products between Ni-201 test samples and 316 L capsules in a 500-hour FLiNaK corrosion experiment has been investigated. In one experiment, the capsule and samples were electrically connected for galvanic coupling, while in the other they were electrically isolated using boron-nitride. Post-corrosion characterization of Ni-201 from the two experiments did not reveal any remarkable differences. Substantial transport of Fe and Cr from the 316 L to the Ni-201 via the salt was observed, with notable diffusion of these elements in the near-surface regions of the Ni-201 samples, while a Cr-rich film was also observed on their surfaces.

36 MATERIALS SCIENCE↗

I – /I 3 – Redox-Assisted Synthesis and Properties of Low Dimensional, Mixed-Valent Gold Iodide Perovskite Derivatives

Here, we report a set of three new mixed-valent Au I Au III iodides: (ClPy) 3 [AuI 2 ] 2 [AuI 4 ] [1], (BrPy) 3 [AuI 2 ] 2 [AuI 4 ] [2], and (ClPy) 2 [AuI 2 ][AuI 4 ] [3], as well as three new monovalent Au III iodides: (XPy) 2 [AuI 4 ][I 3 ] (Py = 4-X-pyridinium X = Cl, Br, and I) [4–6]. Two of these mixed-valent compounds (1 and 2) incorporate both monovalent Au I ···Au I (aurophilic bonding) and mixed-valent Au I I···Au III I couples (Au–I halogen bonding), to the best of our knowledge an unprecedented structural feature. These same two mixed-valent compounds also exhibit a rare low-dimensional molecular architecture with respect to second sphere Au···I interactions, namely, 1D chains of Au···I interactions, extending along a single crystallographic axis. All compounds were synthesized with the assistance of the I – /I 3 – redox couple and tacit manipulation of the polyiodide content during synthesis. Air-free synthesis was used to influence the redox process of I – /I 3 – , resulting in better selection for mixed-valent products. Compounds 1–2 and 4–5 exhibit a characteristically narrow bandgap (1.04–1.25 eV), as measured via diffuse reflectance spectroscopy (DRS). Computational analyses were used to rationalize the specific assembly modes of [AuI 2 ] – and [AuI 4 ] – species, and they show that the Au I ···I interaction type is favored over the Au III ···I.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quenched lattice fluctuations in optically driven SrTiO3

Abstract Crystal lattice fluctuations, which are known to influence phase transitions of quantum materials in equilibrium, are also expected to determine the dynamics of light-induced phase changes. However, they have only rarely been explored in these dynamical settings. Here we study the time evolution of lattice fluctuations in the quantum paraelectric SrTiO 3 , in which mid-infrared drives have been shown to induce a metastable ferroelectric state. Crucial in these physics is the competition between polar instabilities and antiferrodistortive rotations, which in equilibrium frustrate the formation of long-range ferroelectricity. We make use of high-intensity mid-infrared optical pulses to resonantly drive the Ti–O-stretching mode at 17 THz, and we measure the resulting change in lattice fluctuations using time-resolved X-ray diffuse scattering at a free-electron laser. After a prompt increase, we observe a long-lived quench in R-point antiferrodistortive lattice fluctuations. Their enhancement and reduction are theoretically explained by considering the fourth-order nonlinear phononic interactions to the driven optical phonon and third-order coupling to lattice strain, respectively. These observations provide a number of testable hypotheses for the physics of light-induced ferroelectricity.

36 MATERIALS SCIENCE↗

Resonant interactions involving local vibrational modes in crystals

When an impurity with a light mass is inserted into a crystal, it can undergo a high-frequency oscillation referred to as a local vibrational mode (LVM). A Fermi resonance may occur between the LVM and lower-frequency modes of the defect. The LVM may also interact with phonons or the electromagnetic field. Understanding these interactions can help model and control diffusion, defect reactions, and thermal conductivity. LVMs have been probed in semiconductors using pressure and alloying as experimental parameters, resulting in anticrossing between localized and extended vibrational modes. Furthermore, these types of vibrational interactions could play an important role in the stability and thermal properties of organic–inorganic hybrid semiconductors. The coupling between an LVM and electromagnetic wave yields an “LVM polariton,” an excitation that has significant vibrational and electric-field amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurements of nonequilibrium interatomic forces using time-domain x-ray scattering

Here, we demonstrate an experimental approach to determining the excited-state interatomic forces using femtosecond x-ray pulses from an x-ray free-electron laser. We determine experimentally the excited-state interatomic forces that connect photoexcited carriers to the nonequilibrium lattice dynamics in the prototypical Peierls-distorted material, bismuth. The forces are obtained by a constrained least-squares fit of a pairwise interatomic force model to the excited-state phonon dispersion relation as measured by the time- and momentum-resolved x-ray diffuse scattering. We find that photoexcited carriers weaken predominantly the nearest-neighbor forces, which drives the measured softening of the transverse acoustic modes throughout the Brillouin zone as well as the zone-center A 1g optical mode. This demonstrates a bond-selective approach to measuring electron-phonon coupling relevant to a broad range of photoinduced phase transitions and transient light-driven states in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗

Amorphization and siliconization of silicon carbide as a first wall material

The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗