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 1,027 records · Page 57

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗

Dynamics and control of hypersonic aeropropulsive/aeroelastic vehicles

The guidance and control of hypersonic vehicles is examined by studying the airframe/engine/structural-dynamic interactions of a generic vehicle with scramjet propulsion. The pitch-attitude dynamics are described for the vehicle configuration that can sustain hypersonic flight at near-orbital altitudes. These aerospacecraft have strong airframe/engine/elastic coupling in attitude dynamics and engine responses with static instability with respect to pitch. An integrated airframe-engine control system is presented for the control of the system's strong aeropropulsive/aeroelastic coupling. The control methodology utilizes feedback of measured/synthesized values of angle of attack, blended pitch rate, thrust, and combustor-inlet pressure. The multiinput/multioutput engine controller requires high-bandwidth actuation of the fuel-flow control and the effective diffuser ratio. The proposed control laws do not provide optimized performance in terms of pitch response suggesting that additional control crossfeeds and filtering is needed.

Schmidt, David K.↗

Theory and observations of electromagnetic ion cyclotron waves in Saturn's inner magnetosphere

High-resolution Voyager 1 magnetic field observations of Saturn's inner magnetosphere are examined for the presence of ULF waves. Quasi-circular left-hand polarized transverse oscillations are found in the near-equatorial region of 5-7 Rs with a wave period about 10 s and peak amplitude of about 2 nT. The wave is identified as the electromagnetic oxygen cyclotron mode occurring at a frequency just below the O(+) ion cyclotron frequency. A theoretical model of wave excitation based on gyroresonant coupling through a temperature anisotropy of O(+) pickup ions is developed which accounts for the principal features of the wave spectrum. It is hypothesized that wave-particle interactions provide a level of scattering commensurate with the weak pitch angle diffusion regime but nonetheless one that regulates and maintains a constant thermal anisotropy of ions along the magnetic field. Arguments are also presented that O(+) was the dominant thermal ion of the Dione-Tethys plasma torus at the time of the Pioneer 11 encounter the year previous to the Voyager 1 measurements.

Barbosa, D. D.↗

Effect of Rain-Adjusted Satellite Sea Surface Salinity on ENSO Predictions from the GMAO S2S Forecast System

ENSO has a significant impact on climate variability throughout the world and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of altimetry and in situ data lead to improved short-term forecasts of the coupled system. However, few studies have focused on improving the near-surface density/mixing through assimilation of satellite sea surface salinity (SSS). For expediency, most projects that do assimilate SSS do so as if these data were observed at the top model layer instead of at the surface. In rainy regions where buoyant water sits as a fresh lens at the surface, this assumption is likely invalid. Therefore, we adjust SSS so that it more accurately represents the salinity at the first model layer. The Rain Impact Model (RIM) uses a simple diffusion model to determine the near surface salinity gradient (i.e., 1 cm to 5 m). Satellite SSS data are modified using this near-surface salinity gradient, so the salinity values are now valid at the first model layer (we call this SSS@5m). We assess the impact of satellite SSS observations for near-surface dynamics within ocean reanalyses and how these impact dynamical ENSO forecasts using the NASA GMAO S2S coupled system. For all reanalysis experiments, all available along-track sea level and in situ observations are assimilated for 2011 to 2020 using the LETKF scheme (Penny et al., 2013). One reanalysis assimilates Aquarius /SMAP SSS as before. An additional reanalysis is performed assimilating the SSS@5m data. Validation statistics are compared for experiments that assimilate SSS (sub-optimally as before) versus the SSS@5m. We also compare results of coupled forecasts that are initialized from these reanalyses in spring. For all but the big La Niña in 2017, all NINO3.4 forecasts were improved by using the RIM. We will show that improved SSS estimates upgrades density and near-surface mixing leading to more accurate coupled air/sea interaction and better forecasts.

ENSO↗

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A 2D Kaleidoscope of Electron Heat Fluxes Driven by Auroral Electron Precipitation

Electron heat flux is an important value for ionospheric space weather modeling networks. Utilizing the 2D array of Time History of Events and Macroscale Interactions during Substorms all-sky-imager (ASI) observations, Gabrielse et al. (2021, https://doi.org/10.3389/fphy.2021.744298) described a new method that estimates the auroral scale sizes of intense precipitating electron energy fluxes and their mean energies during two substorms on 16 February 2010. These parameters in combination with SuperThermal Electron Transport code were used to develop a new methodology to calculate electron thermal fluxes from data inputs in 2D during one of the substorms at 09:40:00 UT across Canada and Alaska. To test the effect of various precipitation lifetimes on electron heat flux values, boxcar averages ranging from 0 to 900 s were applied to the ASI data. These data are then combined with the newly developed kinetic simulation to determine the thermal fluxes associated with the observed diffuse and discrete precipitation.

Electron precipitation in aurora↗

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Giant interfacial spin Hall angle from Rashba-Edelstein effect revealed by the spin Hall Hanle process

The Rashba-Edelstein effect (REE), which generates interfacial spin polarization and subsequent spin current, is a compelling spin-charge conversion mechanism for spintronic applications, since it’s not limited by the elemental spin-orbit couplings. In this work we demonstrate REE at Pt/ferroelectric interfaces by showing a positive correlation between polarization and effective spin Hall angle in the recently elucidated spin Hall Hanle effects (SHHE), in which a Larmor precession of spin polarization in a diffusion process from the interface manifest as magnetoresistance and Hall effect. We show that REE leads to a large enhancement of the effective spin Hall angle of ferroelectric interface Pt/h-LuFeO 3 compared with that of Pt/Al 2 O 3 , without obvious differences in the spin relaxation time. Modeling using SHHEs involving REE as an additional source of interfacial polarization suggests that REE can lead to an interfacial spin Hall angle (0.4 ± 0.1) in Pt/h-LuFeO 3 that is one order of magnitude larger than the bulk value of Pt. Our results demonstrate that a ferroelectric interface can produce large spin-charge conversion and that SHHEs are a sensitive tool for characterizing interfacial spin-transport properties.

36 MATERIALS SCIENCE↗

Effective field theory of fluctuating wall in open systems: from a kink in Josephson junction to general domain wall

We investigate macroscopic behaviors of fluctuating domain walls in nonequilibrium open systems with the help of the effective field theory based on symmetry. Since the domain wall in open systems breaks the translational symmetry, there appears a gapless excitation identified as the Nambu-Goldstone (NG) mode, which shows the non-propagating diffusive behavior in contrast to those in closed systems. After demonstrating the presence of the diffusive NG mode in the (2+1)-dimensional dissipative Josephson junction, we provide a symmetry-based general analysis for open systems breaking the one-dimensional translational symmetry. A general effective Lagrangian is constructed based on the Schwinger-Keldysh formalism, which supports the presence of the gapless diffusion mode in the fluctuation spectrum in the thin wall regime. Besides, we also identify a term peculiar to the open system, which possibly leads to the instability in the thick-wall regime or the nonlinear Kardar-Parisi-Zhang coupling in the thin-wall regime although it is absent in the Josephson junction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modelling of Laser-Enhanced Chemical Vapor Deposition

Research is directed at development of a detailed model of mass and heat transfer and chemical reaction in the pyrolysis of silane for the growth of thin amorphous silicon substrates incorporating laser heating of the gas phase above the film. The model will be the basis for evaluation of the relative importances of the decomposition of SiH4 in the vapor phase, mass transfer of the intermediate species, e.g., SiH2, and the evolution of hydrogen gas. Plans are also underway for developing a model for homogeneous nucleation of Si in the vapor phase to model the rate limitations observed at high gas-phase temperatures and high partial pressures of silane. Work was concentrated on an almost one-dimensional model for the coupling of the CO2 laser beam for heat transfer of the vapor phase with simple kinetic models for SiH4 decomposition and subsequent absorption of Si vapor on the substrate. Mass transfer in the vapor phase is assumed to be solely by diffusion. The role of convection in the vapor phase caused by the large changes in density in and around the center of the laser beam will be analyzed to evaluate the potential of microgravity experiments for increasing the uniformity of the film and the deposition rate.

Brown, R. A.↗

A computational study of the flowfield surrounding the Aeroassist Flight Experiment vehicle

A symmetric total variation diminishing (STVD) algorithm has been applied to the solution of the three-dimensional hypersonic flowfield surrounding the Aeroassist Flight Experiment (AFE) vehicle. Both perfect-gas and chemical nonequilibrium models have been used. The perfect-gas flows were computed at two different Reynolds numbers, including a flight trajectory point at maximum dynamic pressure, and on two different grids. Procedures for coupling the solution of the species continuity equations with the Navier-Stokes equations in the presence of chemical nonequilibrium are reviewed and tested on the forebody of the AFE and on the complete flowfield assuming noncatalytic wall and no species diffusion. Problems with the STVD algorithm unique to flows with variable thermodynamic properties (real gas) are identified and algorithm modifications are suggested. A potential heating problem caused by strong flow impingement on the nozzle lip in the near wake at 0-deg angle of attack has been identified.

Gnoffo, Peter A.↗

The effect of neutrino transport on the collapse of iron stellar cores

A multigroup flux-limited diffusion approximation to neutrino transport, correct to first order in material velocities, is described. The role of neutrino-electron scattering in determining core deleptonization is studied using a modified Fokker-Planck approximation in which adjacent neutrino energy groups are coupled. The moderate deleptonizaton obtained during infall restricts homologous core masses to less than 0.8 solar mass. With 0.6-0.8 solar mass of outer core to traverse, nuclear dissociation alone stalls the shock waves produced at core bounce. In addition, neutrino energy losses remove at least another 4 x 10 to the 51st ergs. In no case is matter ejected promptly as a direct result of the initial shock wave.

Myra, E. S.↗

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗

Proton-Mediated and Ir-Catalyzed Iron/Iron-Oxide Redox Kinetics for Enhanced Rechargeability and Durability of Solid Oxide Iron–Air Battery

Long duration energy storage (LDES) is an economically attractive approach to accelerating clean renewable energy deployment. The newly emerged solid oxide iron–air battery (SOIAB) is intrinsically suited for LDES applications due to its excellent low-rate performance (high-capacity with high efficiency) and use of low-cost and sustainable materials. However, rechargeability and durability of SOIAB are critically limited by the slow kinetics in iron/iron-oxide redox couples. Here the use of combined proton-conducting BaZr 0.4 Ce 0.4 Y 0.1 Yb 0.1 O 3 (BZC4YYb) and reduction-promoting catalyst Ir to address the kinetic issues, is reported. It is shown that, benefiting from the facilitated H + diffusion and boosted FeOx-reduction kinetics, the battery operated under 550 °C, 50% Fe-utilization and 0.2 C, exhibits a discharge specific energy density of 601.9 Wh kg –1 -Fe with a round-trip efficiency (RTE) of 82.9% for 250 h of a cycle duration of 2.5 h. Under 500 °C, 50% Fe-utilization and 0.2 C, the same battery exhibits 520 Wh kg –1 -Fe discharge energy density with an RTE of 61.8% for 500 h. This level of energy storage performance promises that SOIAB is a strong candidate for LDES applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗