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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1950

Isentropic phase changes in dissociating gases and the method of sound dispersion for the investigation of homogeneous gas reactions with very high speed : Conclusion

The analytical results of Part I are also applied to sound dispersion by friction and heat conduction, An irreversible change of momentum, energy, and type of particle corresponding to friction, heat conduction, and diffusion effects can appear both in the direction of the sound field and traverse to it. Longitudinal damping, the coupling of longitudinal damping and that due to chemical and physical changes, and coupling of diffusion and compositional changes are treated for a plane sound wave of infinite extent. The same principles are also applied to sound effects in cylindrical tubes. The limitations of the method are discussed in some detail.

SOUND VELOCITY - GASES↗

NASA NTRS1977

Two-dimensional hygrothermal diffusions into a finite width composite laminate

Two dimensional hygroscopic diffusion into composite laminates with free edges was obtained through the use of a second order finite difference scheme. A computer program was developed for the transient nonlinear coupled hygrothermal diffusion into finite width laminates with varying surface conditions along two edges. The formulation permits the diffusion coefficient to be a function of temperature, moisture concentration, and fiber orientation. The moisture distributions thus obtained are necessary for analysis of the moisture induced interlaminar stresses in the vicinity of free edges. It is shown that large diffusion times tend to eliminate the significance of stacking sequence, but gradients within and between layers are significant for all diffusion times.

Farley, G. L.↗

DOE OSTI2021

Role of a Multivalent Ion–Solvent Interaction on Restricted Mg 2+ Diffusion in Dimethoxyethane Electrolytes

The diffusion behavior of Mg 2+ in electrolytes is not as readily accessible as that from Li + or Na + utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting 25 Mg NMR. In MgTFSI 2 /DME solutions, “bound” DME (coordinating to Mg 2+ ) and “free” DME (bulk) are distinguishable from 1 H NMR. With the exchange rates between them obtained from 2D 1 H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg 2+ ) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. Further, the high activation enthalpy for exhange (65–70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C–H bond length. Comparison of the diffusion behaviors of Mg 2+ , TFSI – , DME, and Li + reveals a relative restriction to Mg 2+ diffusion that is caused by the long-range interaction between Mg 2+ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Metal Removal by a Free Surface Constructed Wetland and Prediction of Metal Bioavailability and Toxicity with Diffusive Gradients in Thin Films (DGT) and Biotic Ligand Model (BLM)

The H-02 constructed wetland is a free water surface wetland to remove copper (Cu) and zinc (Zn) from the industrial wastewater. In this study, we evaluated the performance of the wetland from 2018 to 2019 and coupled the diffusive gradients in thin films (DGTs) and biotic ligand model (BLM) to explore metal speciation and bioavailability in wetland waters. Surface water samples were collected and piston DGTs were deployed in different sites of the wetland. The H-02 wetland functioned well during the sampling period with high removal efficiencies (Cu: 73.8 ± 1.2% and Zn: 75.2 ± 16.0%). In our study, with the assumption that the combination of BLM predicted inorganic metals species, BLM Cu(II) and BLM Zn(II), were the bioavailable and toxic species, DGT-Cu did not correlate to BLM Cu(II) (P = 0.47), but DGT-Zn positively correlated to BLM Zn(II) (R 2 = 0.35, P < 0.001). Compared to the modeling results of BLM, DGT-indicated labile and/or bioavailable Cu included not only free Cu ions and inorganic Cu complexes but also a high percentage of Cu-labile organic matter complexes. DGT-indicated Zn included free Zn ion, inorganic Zn, and only a low percentage of Zn-labile organic matter complexes. Our findings illustrated the appropriate use of passive sampling techniques and geological modeling when biomonitoring could be substituted. Furthermore, the close monitoring of metal concentrations, speciation, and bioavailability helps us understand metal biogeochemistry and metal removal processes and ensure the long-term sustainability of the constructed wetland.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Interfacial Void Formation and Self-Healing in Oxide Scales on Al-containing High-Entropy Alloy

The exceptional high-temperature oxidation resistance of Al-containing high-entropy alloys (HEAs) is often attributed to the formation of a protective α-Al 2 O 3 scale. However, the dynamic, atomic-scale mechanisms governing the stability of this scale—including interfacial void formation and the often-postulated but rarely visualized “self-healing” capacity—remain poorly understood. Herein, we reveal the complex evolution of the triple-layer oxide scale on an Al 10 CoCrFeNi HEA through combined electron microscopy and diffraction study. We show that interfacial voids are an inherent consequence of the scaling process, originating from two distinct mechanisms: the Kirkendall effect at the interface between the γ-Al 2 O 3 /α-Al 2 O 3 and alloy driven by cationic diffusion imbalance and volumetric contraction due to phase transformations at the spinel/Cr 2 O 3 interface. Crucially, we provide microstructural evidence consistent with an intrinsic self-healing response. This process is driven by coupled inward diffusion of oxygen and outward diffusion of metal cations, leading to the in-situ formation of transient θ-Al 2 O 3 and spinel phases that partially fill and seal the voids. Here, these results provide atomic-scale insights into the phase evolution, defect formation, and self-repair of oxide scales in HEAs—highlighting pathways to enhance their oxidation resistance in extreme environments.

High-entropy alloy↗

DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Ion transport and loss in the Earth's quiet ring current. 2: Diffusion and magnetosphere-ionosphere coupling

We have studied the transport and loss of H(+), He(+), and He(++) ions in the Earth's quiet time ring current (1 to 300 keV/e, 3 to 7 R(sub E), Kp less than 2+, absolute value of Dst less than 11, 70 to 110 degs pitchangles, all LT) comparing the standard radial diffusion model developed for the higher-energy radiation belt particles with measurements of the lower energy ring current ions in a previous paper. Large deviations of that model, which fit only 50% of the data to within a factor of 10, suggested that another transport mechanism is operating in the ring current. Here we derive a modified diffusion coefficient corrected for electric field effects on ring current energy ions that fit nearly 80% of the data to within a factor of 2. Thus we infer that electric field fluctuations from the low-latitude to midlatitude ionosphere (ionospheric dynamo) dominated the ring current transport, rather than high-latitude or solar wind fluctuations. Much of the remaining deviation may arise from convective electric field transport of the E less than 30 keV particles. Since convection effects cannot be correctly treated with this azimuthally symmetric model, we defer treatment of the lowest-energy ions to a another paper. We give chi(exp 2) contours for the best fit, showing the dependence of the fit upon the internal/external spectral power of the predicted electric and magnetic field fluctuations.

Sheldon, R. B.↗

NASA NTRS1995

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗

DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗

DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Short efficient ejector systems

A research program was conducted to improve the performance of low pressure ratio ejector systems. The results show that short, efficient ejectors operating at nearly ideal performance are possible through the use of forced mixer lobes. Forced mixer lobes generate large scale axial vorticity which results in rapid mixing and improved diffuser performance. Ejector testing was conducted using both an ejector wind tunnel. Numerous mixer lobe variations were tested to develop lobe design guidelines. The improved performance of mixer-ejectors is presented over a range of operating conditions. Results of mixer lobe angle, penetration, and alignment are presented. Lobe angles of up to 25 deg, coupled with diffuser wall angles over 20 deg, operate without separation, allowing efficient, short ejector systems. Both warm and cold flow test results are presented. Temperature similarity expressions previously developed are further formulated and verified.

Presz, Walter M., Jr.↗

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 OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

NASA NTRS1991

Dynamic magnetic reconnection at an X-type neutral point

The relaxation of a two-dimensional 'X-type' neutral point magnetic field disturbed from equilibrium is considered. Perturbations are shown to possess well-defined azimuthal modes which allow an exact determination of the magnetic annihilation rate. Free magnetic energy is dissipated by oscillatory reconnection which couples resistive diffusion at the neutral point to global advection of the outer field. The decay of azimuthally symmetric (m = 0) modes - the only modes associated with topological reconnection - is limited by the dissipation time scale of the 'fundamental' (n = 0) mode with no radial nodes. This mode decays over typically 100 Alfven times. An analytic treatment shows that the oscillation and decay time scales couple according to a given law.

Craig, I. J. D.↗

DOE OSTI2024

Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations

This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Binder-Free Graphite Anodes for Next-Generation High-Performance Lithium-Ion Batteries

High-energy density anodes are crucial for next-generation lithium-ion batteries (LIBs) particularly for electric vehicle (EV) applications. Sluggish lithium-diffusion kinetics coupled with conventional anode fabrication processes containing polymeric binders hinder fast-charging capabilities and high-energy density of graphite. Herein, we introduce a binder-free graphite anode fabrication strategy using the electrospinning technique that contains ~2.41% carbon nanotubes (CNTs). Our strategy relies on the formation of an interconnecting conductive CNT network coupled with an ultrathin N-doped carbon coating on graphite particles from sacrificial binders. This combination enhances both structural integrity and electrical conductivity and, in turn, improves fast-charging capabilities and high energy density of LIBs. The binder-free graphite anode achieves ~335.0 mAh g–1 capacity at C/3 rate over 400 cycles with capacity retention of >95% and average Coulombic efficiencies >99.95%. These promising results suggest that the binder-free anode fabrication with a multifunctional design approach could elevate the energy-density limits of the graphite anodes, solving high-energy density requirements of EVs, and potentially provides a path forward for the development of economically feasible energy storage systems for various applications.

Ozcan, Muca [ORNL] (ORCID:0000000320020474)↗