Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuser Flows”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A percolating path to green iron

About 1.9 gigatons of steel is produced every year, emitting 8% (3.6 gigatons) of global CO 2 in the process. More than 50% of the CO 2 emissions come from a single step of steel production, known as ironmaking. Hydrogen- based direct reduction (HyDR) of iron oxide to iron has emerged as an emission-free ironmaking alternative. However, multiple physical and chemical phenomena ranging from nanometers to meters inside HyDR reactors alter the microstructure and pore networks in iron oxide pellets, in ways that resist gaseous transport of H 2 /H 2 O, slow reaction rates, and disrupt continuous reactor operation. Using synchrotron nano X-ray computed tomography and percolation theory, we quantify the evolution of pores in iron oxide pellets and demonstrate how nanoscale pore connectivity influences micro- and macroscale flow properties such as permeability, diffusivity, and tortuosity. Our modeling framework connects disparate scales and offers opportunities to accelerate HyDR.

hydrogen↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2025

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics↗
DOE OSTI2025

A Percolating Path to Green Iron

About 1.9 gigatonnes of steel is produced every year emitting 7% (2.7 gigatonnes) of global CO 2 in the process. More than 50% of the CO 2 emissions come from a single step of steelmaking, known as ironmaking. Hydrogen based direct reduction (HyDR) of iron oxide to iron has emerged as an emissions free ironmaking alternative. However multi-scale phenomena ranging from nanometers to meters inside HyDR reactors exhibit detrimental microstructure evolution which resists gaseous transport of H 2 /H 2 O, slows reaction rates and disrupts continuous reactor operation. To resolve the conundrum between atomic and reactor scales, we devise a percolation-theory model to reconcile nanoscale porosity with macroscopic properties relevant to reactor design models. Using synchrotron nano X-ray computed-tomography, we quantify the evolution of pores in iron oxide pellets, and demonstrate how nano-scale pore networks influence micro and macro-scale flow properties such as permeability, diffusivity and tortuosity. Our new modeling framework bridges the gap between scales and offers the criteria to accelerate HyDR by at least 5x via feedstock-reactor synergies based on percolation.

Paul, Subhechchha↗
DOE OSTI2025

Probabilistic Diffusion Models Advance Extreme Flood Forecasting

Extreme floods pose escalating risks in a changing climate, yet forecasting remains challenging due to peak flow underestimation and high uncertainty. We introduce diffusion-based runoff model (DRUM), a probabilistic deep learning (DL) approach that advances extreme flood forecasting across representative basins in the contiguous United States. DRUM outperforms state-of-the-art benchmarks, enhancing nowcasting skill for the top 1‰ of flows in 72.3% of studied basins. Under operational scenarios, DRUM extends reliable lead times by nearly a full day for 20- and 50-year floods. When evaluated with measured precipitation, an ideal condition, recall improves by 0.3–0.4 and the early warning window extends by 2.3 days for 50-year floods. The enhancement potential varies regionally, with precipitation-driven flood zones in the eastern and northwestern US benefiting most, gaining 3–7 days in lead time. These findings highlight the transformative potential of diffusion models as a cutting-edge generative AI technique for advancing hydrology and broader Earth system sciences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2025

Mitigating spectral bias in neural operators via high-frequency scaling for physical systems

Neural operators have emerged as powerful surrogates for modeling complex physical problems. However, they suffer from spectral bias making them oblivious to high-frequency modes, which are present in multiscale physical systems. Therefore, they tend to produce over-smoothed solutions, which is particularly problematic in modeling turbulence and for systems with intricate patterns and sharp gradients such as multi-phase flow systems. In this work, we introduce a new approach named high-frequency scaling (HFS) to mitigate spectral bias in convolutional-based neural operators. By integrating HFS with proper variants of UNet, we demonstrate a higher prediction accuracy by mitigating spectral bias in single and two-phase flow problems. Unlike Fourierbased techniques, HFS is directly applied to the latent space, thus eliminating the computational cost associated with the Fourier transform. Additionally, we investigate alternative spectral bias mitigation through a diffusion model conditioned on neural operators. While the diffusion model integrated with the standard neural operator may still suffer from significant errors, these errors are substantially reduced when the diffusion model is integrated with a HFS-enhanced neural operator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Impact of the Out‐Of‐Plane Flow Shear on Magnetic Reconnection at the Flanks of Earth's Magnetopause

Abstract Magnetic reconnection changes the magnetic field topology and facilitates the energy and particle exchange at magnetospheric boundaries such as the Earth's magnetopause. The flow shear perpendicular to the reconnecting plane prevails at the flank magnetopause under southward interplanetary magnetic field conditions. However, the effect of the out‐of‐plane flow shear on asymmetric reconnection is an open question. In this study, we utilize kinetic simulations to investigate the impact of the out‐of‐plane flow shear on asymmetric reconnection. By systematically varying the flow shear strength, we analyze the flow shear effects on the reconnection rate, the diffusion region structure, and the energy conversion rate. We find that the reconnection rate increases with the upstream out‐of‐plane flow shear, and for the same upstream conditions, it is higher at the dusk side than at the dawn side. The diffusion region is squeezed in the outflow direction due to magnetic pressure which is proportional to the square of the Alfvén Mach number of the shear flow. The out‐of‐plane flow shear increases the energy conversion rate , and for the same upstream conditions, the magnitude of is larger at the dusk side than at the dawn side. This study reveals that out‐of‐plane flow shear not only enhances the reconnection rate but also significantly boosts energy conversion, with more pronounced effects on the dusk‐side flank than on the dawn‐side flank. These insights pave the way for better understanding the solar wind‐magnetosphere interactions.

Liang, Haoming↗
DOE OSTI2025

Generative AI for Wind-Wave Inflows

Extreme wind-wave conditions can create dangerous or damaging environmental conditions for sea-based structures and watercrafts. However, the model of consistent, coupled wind-wave flow conditions is computational expensive, even for simple flow scenarios. This work examines the use of diffusion-based generative artificial intelligence tools to accelerate the generation of these two-phase flows and to characterize probabilities of various states to help with the generation of extreme patterns.

17 WIND ENERGY↗
DOE OSTI2026

Cation-Diffusive Carbon Interlayers Stabilize Na Metal and Double the Current in Na-S Redox-Flow Batteries for Grid-Scale Energy Storage

The sodium-sulfur nonaqueous redox-flow batteries (Na-S NARFBs) using earth-abundant elements are highly attractive due to the low material cost and improved energy density for grid-scale energy storage. However, the low current performance, poor Na0/Na+ redox kinetics, and Na dendrite growth pose severe challenges. We introduce cation-diffusive layers (CDLs): thin and Na+ affinitive interlayers at the Na anode that direct Na+ transport and stabilize Na deposition. Benchmarking three archetypal materials—carbon paper (CP), glass microfiber paper (GF), and foam—across Na-Na and Na-Cu, and Na-S cells identifies CP as the optimum. CP reduces symmetric cell overpotential by more than 70%, achieves 98% Na plating-stripping efficiency, and doubles the Na-S cell current density from 0.5 to 1.0 mA cm−2 without sacrificing capacity or efficiency. Ex situ electrochemical and SEM/XPS analysis, combined with molecular dynamics (MD) studies, reveal that electron-rich carbon fibers disperse supporting salt aggregates, enrich near-surface Na+ density, and create ion transport pathways for fast Na0/Na+ exchange while mitigating membrane degradation. Because of the ion-centric mechanism, CDLs can be generalized to other metal-anode designs. Further, this work establishes CDL design rules—cationic affinity and appropriate micro/nanostructure—as a simple, scalable route to high-current, durable metal-anode flow batteries.

Wu, Wenda [ORNL] (ORCID:000900033307687X)↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Spin-caloric resistance of Dirac plasma in a graphene Corbino device

The thermal resistance of a spin-polarized hydrodynamic Dirac plasma in graphene is considered. Here, a mechanism for the coupling of heat and spin flows is discussed, demonstrating that spin diffusion and spin thermocurrent modify viscous dissipation, leading to a significant enhancement of thermal resistance. Practical calculations are then presented for graphene devices in the Corbino geometry.

Graphene↗
DOE OSTI2024

Understanding Selectivity Control in the Electrocatalytic Reduction of CO2 to Liquid Products in Gas-Fed Electrolyzers

The room-temperature electrochemical reduction of carbon dioxide to liquid products is a soaring carbon utilization technology with an energy and environmental impact, offering a pathway to convert renewable energy into valuable C1 (e.g., formic acid) and C2+ products (e.g., ethanol and n-propanol)1. Gas-fed flow electrolyzers, in which a gas diffusion layer is used to transport gaseous CO2 into the electrode, have emerged as promising electrocatalytic reactors for large-scale applications, reaching competitive production costs for carbon monoxide (CO) and formic acid (HCOOH). Despite their increased use in recent years, several factors governing their performance have yet to be understood.

Berch, John El↗