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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2020

Fission Product Diffusion in Silicon Carbide: As Revealed by Computer Simulations and Experimental Measurements

The diffusion and release of fission products through silicon carbide in the tri-structural isotropic (TRISO) fuel particles remain unsolved for decades. The underlying mechanism is quite challenging to be determined. To help unveil the mysterious story, the current work applies molecular dynamics method to show the stability of silver, palladium, ruthenium and iodine as an interstitial and their atomic diffusion along coincident site lattice (CSL) boundary, especially ?3 grain boundary (GB). The major finding presents a much faster diffusion along GB than in bulk for all elements considered. The reasonably close estimate to experiments and simulations where available has confirmed the important role of grain boundary diffusion of Ag and Pd in SiC. However, the discrepancy addressed in Ag with measurements from fuel studies suggest a more complicated mechanism, which might be in correlation with high energy grain boundaries or the presence of crack. The subsequent characterization of Ru and I distribution in SiC-PyC-SiC diffusion couples, which have been ion irradiated at 900?C to 10 dpa and 20 dpa, has performed by secondary ion mass spectrometry (SIMS) analysis. The experiment measurements correlate well with the grain boundary diffusion by simulation, which provide further evidences that the grain boundary diffusion cannot be neglected once the fission products are accessible at the grain boundary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1980

Dynamic properties of the thermosphere inferred from Pioneer Venus mass spectrometer measurements

The data obtained by the Pioneer Venus spectrometer experiments indicate that the day-night temperature contrast on Venus is associated with wind velocities of about 200 m/s which transport oxygen, helium, and hydrogen toward the night side. A mass exchange with the mesosphere commensurate with an eddy diffusion coefficient of 3 x 10 to the 7th is required to buffer the horizontal advection so as to reproduce the observed day time bulge in oxygen and the small diurnal variations in helium. The observed time response and magnitude of the day-night density variations require transport processes to be effective over time periods between five and ten days, implying a superrotation rate or prevailing winds in excess of 50 m/s at the equator. Nonlinear mass transport results in wave steepening and contributes to the amplification of the density extrema in hydrogen and helium.

Mayr, H. G.↗
NASA NTRS1975

Film boiling of mercury droplets

Vaporization times of mercury droplets in Leidenfrost film boiling on a flat horizontal plate are measured in an air atmosphere. Extreme care was used to prevent large amplitude droplet vibrations and surface wetting; therefore, these data can be compared to film boiling theory. Diffusion from the upper surface of the drop appears as a dominant mode of mass transfer from the drop. A closed-form analytical film boiling theory is developed to account for the diffusive evaporation. Reasonable agreement between data and theory is seen.

Baumeister, K. J.↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
NASA NTRS1993

Alveolar ventilation to perfusion heterogeneity and diffusion impairment in a mathematical model of gas exchange

This study describes a two-compartment model of pulmonary gas exchange in which alveolar ventilation to perfusion (VA/Q) heterogeneity and impairment of pulmonary diffusing capacity (D) are simultaneously taken into account. The mathematical model uses as input data measurements usually obtained in the lung function laboratory. It consists of two compartments and an anatomical shunt. Each compartment receives fractions of alveolar ventilation and blood flow. Mass balance equations and integration of Fick's law of diffusion are used to compute alveolar and blood O2 and CO2 values compatible with input O2 uptake and CO2 elimination. Two applications are presented. The first is a method to partition O2 and CO2 alveolar-arterial gradients into VA/Q and D components. The technique is evaluated in data of patients with chronic obstructive pulmonary disease (COPD). The second is a theoretical analysis of the effects of blood flow variation in alveolar and blood O2 partial pressures. The results show the importance of simultaneous consideration of D to estimate VA/Q heterogeneity in patients with diffusion impairment. This factor plays an increasing role in gas alveolar-arterial gradients as severity of COPD increases. Association of VA/Q heterogeneity and D may produce an increase of O2 arterial pressure with decreasing QT which would not be observed if only D were considered. We conclude that the presented computer model is a useful tool for description and interpretation of data from COPD patients and for performing theoretical analysis of variables involved in the gas exchange process.

NASA Discipline Environmental Health↗
DOE OSTI2021

Large-eddy simulation of a multi-injection flame in a diesel engine environment using an unsteady flamelet/progress variable approach

In this work, large-eddy simulations (LESs) are conducted for a multiple-injection flame in a diesel engine environment using an unsteady flamelet/progress variable (UFPV) approach in which differential diffusion is considered. The suitability of the UFPV tabulation approach is first evaluated through a priori analyses using the state-of-the-art direct numerical simulation (DNS) dataset. Both the instantaneous data and the conditional values for the major and minor species' mass fractions are compared between the UFPV and the DNS. The comparisons show that the proposed UFPV tabulation approach can give good predictions for the multiple-injection flame at different injection phases. While the gas temperature and major species mass fractions can be accurately predicted with or without differential diffusion being considered in the UFPV flamelet library, the prediction accuracy for the highly diffusive species (e.g., hydrogen) in the main injection phase can be noticeably improved when differential diffusion is taken into account. The fully coupled LES/UFPV simulations show that the overall structure of the multiple-injection flame can be predicted, and the conditional thermo-chemical values are close to the filtered DNS dataset. The reasons for the remaining discrepancies found in the a priori analyses and the a posteriori simulations using the UFPV approach are analyzed.

Wen, Xu (ORCID:0000000331334292)↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2021

Bulk and Short-Circuit Anion Diffusion in Epitaxial Fe 2 O 3 Films Quantified Using Buried Isotopic Tracer Layers

Self-diffusion is a fundamental physical process that, in solid materials, is intimately correlated with both microstructure and functional properties. In this work, a universal approach is presented to precisely characterize self-diffusion in ionic solids by isotopically enriching anions and/or cations at specific locations within an epitaxial film stack, and characterize their redistribution at high spatial resolution with atom probe tomography. Nanoscale anion diffusivity is quantified in epitaxial α-Fe 2 O 3 thin films deposited by molecular beam epitaxy with a thin (10 nm) buried layer highly enriched in 18 O. The isotopic sensitivity of the atom probe allows precise measurement of 18 O distribution across the sharp interfaces between this layer and the surrounding Fe 2 O 3 after annealing. Short-circuit anion diffusion through 1D and 2D structural defects in Fe 2 O 3 are also directly visualized in 3D. This versatile approach to study precisely tailored thin film samples at high spatial and mass fidelity will facilitate a deeper understanding of atomic-scale diffusion phenomena.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Modeling of turbulent supersonic H2-air combustion with a multivariate beta PDF

Recent calculations of turbulent supersonic reacting shear flows using an assumed multivariate beta PDF (probability density function) resulted in reduced production rates and a delay in the onset of combustion. This result is not consistent with available measurements. The present research explores two possible reasons for this behavior: use of PDF's that do not yield Favre averaged quantities, and the gradient diffusion assumption. A new multivariate beta PDF involving species densities is introduced which makes it possible to compute Favre averaged mass fractions. However, using this PDF did not improve comparisons with experiment. A countergradient diffusion model is then introduced. Preliminary calculations suggest this to be the cause of the discrepancy.

Baurle, R. A.↗
DOE OSTI2023

Phosphorus diffusion and deactivation during SiGe oxidation

Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.

Thornton, Chappel S.↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

Corotation lag limit on mass-loss rate from Io

Considering rapid escape of H2O from Io during an early hot evolutionary epoch, an H2O plasma torus is constructed by balancing dissociation and ionization products against centrifugally driven diffusion, including for the first time the effects of corotation lag resulting from mass loading. Two fundamental limits are found as the mass injection rate increases: (1) an 'ignition' limit of 1.1 x 10 to the 6th kg/s, beyond which the torus cannot ionize itself and photoionization dominates; and (2) the ultimate mass loading limit of 1.3 x 10 to the 7th kg/s, which occurs when neutrals newly created by charge exchange and recombination cannot leave the torus, thereby bringing magnetospherically driven transport to a halt. Connecting this limit with the variations of Io's temperature in its early evolution epoch gives an estimate of the upper limit on the total mass loss from Io, about 3.0 x 10 to the 20th kg (for high-opacity nebula) and about 8.9 x 10 to the 20th kg (for low-opacity nebula). These limits correspond to eroding 8 km and 22 km of H2O from the surface. It is concluded that compared to the other Galilean satellites, Io was created basically dry.

Huang, T. S.↗
NASA NTRS1997

Numerical Modeling of HgCdTe Solidification: Effects of Phase Diagram, Double-Diffusion Convection and Microgravity Level

Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.

Bune, Andris V.↗
NASA NTRS2011

An Examination of the Effect of Boundary Layer Ingestion on Turboelectric Distributed Propulsion Systems

A Turboelectric Distributed Propulsion (TeDP) system differs from other propulsion systems by the use of electrical power to transmit power from the turbine to the fan. Electrical power can be efficiently transmitted over longer distances and with complex topologies. Also the use of power inverters allows the generator and motors speeds to be independent of one another. This decoupling allows the aircraft designer to place the core engines and the fans in locations most advantageous for each. The result can be very different installation environments for the different devices. Thus the installation effects on this system can be quite different than conventional turbofans where the fan and core both see the same installed environments. This paper examines a propulsion system consisting of two superconducting generators, each driven by a turboshaft engine located so that their inlets ingest freestream air, superconducting electrical transmission lines, and an array of superconducting motor driven fan positioned across the upper/rear fuselage area of a hybrid wing body aircraft in a continuous nacelle that ingests all of the upper fuselage boundary layer. The effect of ingesting the boundary layer on the design of the system with a range of design pressure ratios is examined. Also the impact of ingesting the boundary layer on off-design performance is examined. The results show that when examining different design fan pressure ratios it is important to recalculate of the boundary layer mass-average Pt and MN up the height for each inlet height during convergence of the design point for each fan design pressure ratio examined. Correct estimation of off-design performance is dependent on the height of the column of air measured from the aircraft surface immediately prior to any external diffusion that will flow through the fan propulsors. The mass-averaged Pt and MN calculated for this column of air determine the Pt and MN seen by the propulsor inlet. Since the height of this column will change as the amount of air passing through the fans change as the propulsion system is throttled, and since the mass-average Pt and MN varies by height, this capture height must be recalculated as the airflow through the propulsor is varied as the off-design performance point is converged.

Felder, James L.↗
NASA NTRS1997

Measurement of Atomic Oxygen in Diffuse Aurora and Ion Density in the E-Region

An ion mass spectrometer (IMS) was refurbished, calibrated and supplied to the University of Colorado payload (Dr. Charles Barth, P.I.) which was launched from White Sands in September of 1993 as NASA 33.062. The nose cone failed to deploy and their were problems with the ACS so the mission was declared a failure. However, the door covering the IMS deployed and the instrument obtained data. The launch occurred shortly after a payload carrying solar x-ray detectors was launched. Thus a small portion of the Colorado payload science was salvaged; namely, the NO(+)/O2(+) ratio to compare with the measured x-ray flux. Figure I shows the NO(+) to O2(+) ratio vs. altitude. The behavior is typical of the E-region.

Sharp, William E.↗