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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Evaluation of transport conditions during physical vapor transport growth of opto-electronic crystals

Transport conditions were evaluated during the vapor phase growth of mercurous chloride crystals in a closed tube. Experimentally observed growth rates were much smaller than those calculated by the Hertz-Knudsen (H-K) equation. The Arrhenius behavior of growth rate with the temperature was used to derive the sticking coefficient. A one-dimensional diffusion model was used to calculate the total mass flux and was compared with the condensing flux. It was predicted that growth occurred in the convecto-diffusive range.

Singh, N. B.↗
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↗
NASA NTRS1989

Gravimetric measurements of materials outgassing applied to graphite-epoxy laminates

The outgassing rates of two graphite-epoxy laminates, American Cyanamide 985B-626 and HST-7B-112, were obtained using a gravimetric method. The rates as a function of time and temperature were derived from the measurements of their mass losses at temperatures varying from 25 to 150 C and for a time span of up to 400 hours in a vacuum. The data from those measurements were reduced to obtain the outgassing activation energies, the mass losses per unit mass or area, and the corresponding outgassing rates. The rates are expressed in closed-form equations and are directly usable for medling computations. The procedures to obtain these parameters are shown and may be used for the evaluation of other materials. The results of the tests show that the activation energies of the two materials are: 4630 cal/mole for the 985B-626 materials and 4791 cal/mole for the HST-7B-112 sample no. 10 Graphite Exoxy. The outgassing rates of these materials are in the 10E-5 g/sq cm/hr range and they decay according to a power of time of 0.60 at 25 C, indicating that the outgassing process is mainly a diffusion at that temperature. The normalized mass losses versus time obtained from these tests were compared to the discrete results obtained from the ASTM-E595 tests. The comparison provides general indications on the effects of temperature and time in relation to the ASTM test values obtained at 125 C for a 24-hour test duration.

Scialdone, John J.↗
NASA NTRS1990

Gravimetric measurements of materials outgassing applied to graphite-epoxy laminates

The outgassing rates of two graphite-epoxy laminates, American Cyanamide 985B-626 and HST-7B-112, were obtained using a gravimetric method. The rates as a function of time and temperature were derived from the measurements of their mass losses at temperatures varying from 25 to 150 C and for a time span of up to 400 hours in a vacuum. The data from those measurements were reduced to obtain the outgassing activation energies, the mass losses per unit mass or area, and the corresponding outgassing rates. The rates are expressed in closed-form equations and are directly usable for modeling computations. The procedures to obtain these parameters are shown and may be used for the evaluation of other materials. The results of the tests show that the activation energies of the two materials are: 4630 cal/mole for the 985B-626 materials and 4791 cal/mole for the HST-7B-112 sample no. 10 Graphite Epoxy. The outgassing rates of these materials are in the 1OE-5 g sq cm/hr range and they decay according to a power of time of 0.60 at 25 C, indicating that the outgassing process is mainly a diffusion at that temperature. The normalized mass losses versus time obtained from these tests were compared to the discrete results obtained from the ASTM-E595 tests. The comparison provides general indications on the effects of temperature and time in relation to the ASTM test values obtained at 125 C for a 24-hour test duration.

Scialdone, John J.↗
NASA NTRS1991

Plasma transport in the Io torus - The importance of microscopic diffusion

This paper considers the question of whether the distribution of mass in the Io plasma torus is consistent with the concept of interchange eddy transport. Specifically, the flux tube content exhibits a gradual decrease with increasing radial distance from the source near Io without any evidence for substantial density irregularity associated with the plasma source or loss. Using a simple one-dimensional numerical model to simulate macroscopic interchange eddy transport, it is demonstrated that this smooth equilibrium distribution of mass can occur but only with the inclusion of a minimal level of small scale microscopic mixing at a rate approaching Bohm diffusion. Otherwise, the system exhibits a chaotic appearance which never approaches an equilibrium distribution. Various physical mechanisms for the microscopic diffusion process which is required to provide a sufficiently rapid mixing of material between the macroscopic eddies are discussed.

Mei, YI↗
NASA NTRS1979

Interfacial geometry and D-variation effects in two-phase systems

Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.

Tenney, D. R.↗
NASA NTRS1973

Investigation of the environmental change pattern of Japan

The author has identified the following significant results. ERTS-1 imagery clearly identifies the relationships between the status of erosion, effluent patterns affected by the coastal current, and the cultural construction activities. Simple photographic techniques can be used for detecting water mass distribution separately from cloud cover and also noise caused by reflected sunlight from wave surfaces. Polluted water does not diffuse continuously into the oceanic water, but forms masses in the water in the Kuroshio area. The polluted or turbid water in the area just north of the Tomogashima Channel, the south outlet of the Osaka Bay, shows that the northward tidal current runs in a clockwise eddy at the tidal period when the imagery was taken. Such an eddy-like pattern of tidal current had never been revealed by conventional oceanographic data. A front between an oceanic water mass and a polluted water mass runs in a NW-SE direction in the central part of the Osaka Bay. The patterns of turbid water discharged from the Kii River and Yoshino River show a northward tidal current in the North Kii Straits. The pattern of lighter turbid or polluted water located in the northwest region of the North Kii straits suggests the existence of a clockwise eddy in the straits.

Maruyasu, T.↗
NASA NTRS1994

The Prediction of Nozzle Performance and Heat Transfer in Hydrogen/Oxygen Rocket Engines with Transpiration Cooling, Film Cooling, and High Area Ratios

An advanced engineering computational model has been developed to aid in the analysis of chemical rocket engines. The complete multispecies, chemically reacting and diffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and Dufour energy transfer terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film-cooled nozzle, and a transpiration-cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film-cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. It is demonstrated that thermal diffusion has a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle and was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration-cooled rocket engine.

Kacynski, Kenneth J.↗
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↗
NASA NTRS1988

Constraints on gravity wave induced diffusion in the middle atmosphere

A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H2O and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is, at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.

Strobel, Darrell F.↗
NASA NTRS1989

Constraints on gravity wave induced diffusion in the middle atmosphere

A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H20 and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.

Strobel, Darrell F.↗
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↗