Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Internal gravity waves in massive stars: II. Frequency analysis across stellar mass

Stars that are over 1.6 solar masses are generally known to possess convective cores and radiative envelopes, which allows for the propagation of outwardly travelling internal gravity waves (IGWs). Here, we study the generation and propagation of IGWs in such stars using two-dimensional, fully non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for five different stellar masses at the middle of the main sequence: 3 M ⊙ , 5 M ⊙ , 7 M ⊙ , 10 M ⊙ , and 13 M ⊙ . From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compared our simulation results with stellar photometric data from a recent survey and we found that for frequency intervals above 8 μHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars. This indicates that the brightness variations are caused by core-generated IGWs.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2020

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗
DOE OSTI2024

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

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

Non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor instability using the macroscopic forcing method

The importance of non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor Instability (RTI) is investigated. The macroscopic forcing method is utilized to measure spatio-temporal moments of the eddy diffusivity kernel representing passive scalar transport in the ensemble averaged fields. Presented in this work are several studies assessing the importance of the higher-order moments of the eddy diffusivity, which contain information about non-locality, in models for RTI. First, it is demonstrated through a comparison of leading-order models that a purely local eddy diffusivity is insufficient to capture the mean field evolution of the mass fraction in RTI. Therefore, higher-order moments of the eddy diffusivity operator are not negligible. Models are then constructed by utilizing the measured higher-order moments. It is demonstrated that an explicit operator based on the Kramers–Moyal expansion of the eddy diffusivity kernel is insufficient. An implicit operator construction that matches the measured moments is shown to offer improvements relative to the local model in a converging fashion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2023

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2022

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗
DOE OSTI2025

Enhancing water and oxygen transport through electrode engineering for AEM water electrolyzers

Anion-exchange membrane water electrolyzers (AEMWEs) can accelerate the deployment of more efficient and affordable hydrogen production solutions. Here, electrode structure is shown to affect water back-diffusion and oxygen transport, which, in return, governs overpotential behaviors in AEMWEs. Measurements indicate that electrode with copious catalytic sites produces water close to the AEM, creating a higher water gradient and driving water back-diffusion, which improves membrane hydration and mass transport. In situ measurement reveals a high pH gradient near the anode surface, which affects anode kinetics. Operando measurement shows reduced oxygen accumulation when decoupling oxygen production and transport on anode. Catalyst ink rheology and stability are tuned with additives to realize scalable fabrication of electrodes with enhanced transport features, allowing AEMWE to operate at 2 A cm−2 for over 1,000+ h at a 2.3 μV h−1 degradation rate. Analysis during and post-durability provides insights into degradation mechanisms. This work demonstrates an electrode design strategy for efficient and durable AEMWEs.

AEM↗
DOE OSTI2025

Regional-Scale Modeling Parameterizations for Secondary Organic Aerosol Formation from Isoprene Epoxydiols: Experimentally Based Evaluation and Optimization

Isoprene is an abundant volatile organic compound emitted from broadleaf forests. Under low nitric oxide concentrations, isoprene is photochemically oxidized to form gas-phase isoprene epoxydiols (IEPOX). In the presence of acidified sulfate aerosols, IEPOX enhances the secondary organic aerosol (SOA) formation. Predictions of IEPOX-SOA in regional-scale models, e.g., the Community Multiscale Air Quality Model (CMAQ), are uncertain due to homogeneous aerosol assumptions, underpredictions of water uptake (hygroscopicity), and aerosol surface area. Here, we used experimental measurements of IEPOX-SOA tracers, 2-methyltetrols (2-MT) and 2-methyltetrol sulfates (2-MTS), formed at initial IEPOX-to-inorganic sulfate ratios ranging from 1–10.5, at ∼50% relative humidity to constrain key IEPOX-SOA parameters: phase separation, organic shell diffusivity (D org ), acidity, hygroscopic growth, mass accommodation, and kinetics. The base CMAQ parametrization overpredicted experimental IEPOX-SOA with an average normalized mean bias (NMB average ) of 1.63. CMAQ with phase separation underpredicted IEPOX-SOA (NMB average = −0.71). Using the phase-separated model, CMAQ model performance was optimized (NMB average = 0.077) with an increased D org = 2 × 10 –16 m 2 s –1 and increased rate constants (k 2-MT = 1 × 10 –3 M 2 s –1 , k 2-MTS = 8.83 × 10 –3 M 2 s –1 ). The optimized model explicitly accounted for hygroscopic growth by utilizing experimentally derived growth rates, improving aerosol surface area predictions. Our model highlights the importance of the aerosol mixing state (homogeneous versus phase-separated), aerosol size dynamics, and hygroscopic growth in modeling heterogeneous reactive uptake of IEPOX.

aerosols↗
DOE OSTI2021

Enhancing carbon dioxide gas-diffusion electrolysis by creating a hydrophobic catalyst microenvironment

Electroreduction of carbon dioxide (CO 2 ) over copper-based catalysts provides an attractive approach for sustainable fuel production. While efforts are focused on developing catalytic materials, it is also critical to understand and control the microenvironment around catalytic sites, which can mediate the transport of reaction species and influence reaction pathways. Here, we show that a hydrophobic microenvironment can significantly enhance CO 2 gas-diffusion electrolysis. For proof-of-concept, we use commercial copper nanoparticles and disperse hydrophobic polytetrafluoroethylene (PTFE) nanoparticles inside the catalyst layer. Consequently, the PTFE-added electrode achieves a greatly improved activity and Faradaic efficiency for CO 2 reduction, with a partial current density >250 mA cm -2 and a single-pass conversion of 14% at moderate potentials, which are around twice that of a regular electrode without added PTFE. The improvement is attributed to a balanced gas/liquid microenvironment that reduces the diffusion layer thickness, accelerates CO 2 mass transport, and increases CO 2 local concentration for the electrolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gas diffusion electrodes, reactor designs and key metrics of low-temperature CO 2 electrolysers

We report that CO 2 emissions can be recycled via low-temperature CO 2 electrolysis to generate products such as carbon monoxide, ethanol, ethylene, acetic acid, formic acid and propanol. In recent years, progress has been made towards an industrially relevant performance by leveraging the development of gas diffusion electrodes (GDEs), which enhance the mass transport of reactant gases (for example, CO 2 ) to the active electrocatalyst. Innovations in GDE design have thus set new benchmarks for CO 2 conversion activity. In this Review, we discuss GDE-based CO 2 electrolysers, in terms of reactor designs, GDE composition and failure modes, to identify the key advances and remaining shortfalls of the technology. This is combined with an overview of the partial current densities, efficiencies and stabilities currently achieved and an outlook on how phenomena such as carbonate formation could influence the future direction of the field. Our aim is to capture insights that can accelerate the development of industrially relevant CO 2 electrolysers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Defect-mediated diffusion of implanted Mg in GaN: Suppressing dopant redistribution by sequential thermal and microwave annealing

The diffusion behavior of Mg in Mg/N co-implanted GaN is investigated in response to a set of annealing conditions and methodologies, namely, 1000 °C/30 min thermal anneal, by high-temperature pulsed gyrotron microwave annealing at 1420 or 1500 °C, or by thermal and microwave annealing, sequentially. After 1000 °C annealing, the diffusion of Mg in GaN is found to be negligible, as measured by secondary ion mass spectrometry. Annealing by gyrotron microwave annealing alone induces the diffusion of Mg at a rate on the order of 10 −12 cm 2 /s. However, the use of a thermal anneal before microwave gyrotron annealing reduces this rate by an order of magnitude to 10 −13 cm 2 /s. We find that a model that considers Mg diffusion from an inhomogeneous medium that contains a defect-rich implanted region near-surface to a relatively pristine region below the implant range better explains the observed diffusion behavior than a conventional model that assumes a homogeneous medium. By analyzing the diffusion behavior using the Boltzmann–Matano method, we present a discussion of reduction in [V Ga ] by thermal annealing at 1000 °C, leading to a suppressed diffusion coefficient during subsequent high-temperature annealing relative to diffusion after 1420/1500 °C annealing alone. This effect holds potential for improvement in the precision of selectively doped regions for future applications based on the (Al)GaN material system. An improved field profile control in real devices can increase the breakdown and current-handling capabilities in power electronic applications.

Meyers, V.↗
DOE OSTI2022

Magnetic damping in ferromagnetic/heavy-metal systems: The role of interfaces and the relation to proximity-induced magnetism

Damping and spin transport in spintronic multilayered systems continues to be a topic of active research. The enhancement of damping in ferromagnet (FM)/spacer layer (SL)/heavy-metal (HM) thin-film systems was studied for Co 25 Fe 75 /SL/Pt with a nonmagnetic (NM) SL of either Au or Cu with variable thickness, in order to understand the correlation with proximity-induced magnetism (PIM) in the HM. Structural, PIM and magnetic damping measurements were undertaken on the same samples. Specifically, secondary ion mass spectroscopy, element specific x-ray magnetic reflectivity and x-ray magnetic circular dichroism at the Pt and Au L 3 edges, and ferromagnetic resonance methods were used. With increasing thickness of a Cu or Au SL directly between the FM and the Pt layer, the Pt PIM and the damping both fall rapidly, with a relationship between damping and PIM that depends on the SL material. The PIM observed in the Au layer showed a complex dependence on the layer thickness, suggesting some hybridization with the Pt. The role of the number and location of interfaces on the damping was demonstrated with the addition of a SL within the Pt layer, which showed that the specific details of the NM/HM interface also affects the damping. The insertion of a Cu SL within the Pt showed a measurable increase in the overall enhancement of the damping while the insertion of a Au SL into Pt had almost no effect on the damping. Finally, together these results demonstrate the role of both PIM and of additional interfaces in the enhancement of damping in FM/HM systems, which is not fully accounted for by existing theory

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A chemical kinetics theory for interpreting the non-monotonic temperature dependence of enzymatic reactions

Abstract. One notable observation of enzymatic chemical reactions is that, for a given abundance of enzymes and substrates, temperature increases cause reaction rates to first increase consistent with the Arrhenius relationship, then plateau, and finally fall off quickly to zero at high temperatures. While many mathematical functions have been used to describe this pattern, we here propose a chemical kinetics theory which successfully replicates this observation and provides insights into the processes responsible for these dynamics. The chemical kinetics theory combines the law of mass action, von Smoluchowski's diffusion-limited chemical reaction theory, and Eyring's transition state theory. This new theory reveals that the thermally reversible enzyme denaturation ensured by the ceaseless thermal motion of molecules and ions in an enzyme solution explains the plateau and subsequent decrease in chemical reaction rates with increasing temperature. The temperature-dependent affinity parameter (K) that relates enzymes and substrates through their binding also affects the shape of the emergent temperature response. We demonstrate that with an increase in substrate availability, K shifts the optimal temperature, where reaction rates plateau, towards higher values. Further, we show that the chemical kinetics theory accurately represents 12 sets of published enzyme assay data and includes the popular mechanistic model by Ratkowsky et al. (2005) as a special case. Given its good performance and solid theoretical underpinning, we believe this new theory will facilitate the construction of more mechanistic-based environmental biogeochemical models.

Tang, Jinyun↗
DOE OSTI2023

Aggregating in vitro-grown adipocytes to produce macroscale cell-cultured fat tissue with tunable lipid compositions for food applications

We present a method of producing bulk cell-cultured fat tissue for food applications. Mass transport limitations (nutrients, oxygen, waste diffusion) of macroscale 3D tissue culture are circumvented by initially culturing murine or porcine adipocytes in 2D, after which bulk fat tissue is produced by mechanically harvesting and aggregating the lipid-filled adipocytes into 3D constructs using alginate or transglutaminase binders. The 3D fat tissues were visually similar to fat tissue harvested from animals, with matching textures based on uniaxial compression tests. The mechanical properties of cultured fat tissues were based on binder choice and concentration, and changes in the fatty acid compositions of cellular triacylglyceride and phospholipids were observed after lipid supplementation (soybean oil) during in vitro culture. This approach of aggregating individual adipocytes into a bulk 3D tissue provides a scalable and versatile strategy to produce cultured fat tissue for food-related applications, thereby addressing a key obstacle in cultivated meat production.

Yuen Jr, John Se Kit (ORCID:0000000198541654)↗