Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Wood, William A.↗
NASA NTRS2001

The Oxidation Kinetics of Continuous Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental observations and results suggest two primary regimes as a function of temperature, i.e., diffusion and reaction controlled kinetics. Thermogravimetric analysis of carbon fiber in flowing oxygen gave an activation energy of 64.1 kJ/mol in the temperature range of 500 to 600 C and an apparent activation energy of 7.6 kJ/mol for temperatures from 600 to 1400 C. When C/SiC composite material was unstressed, matrix effects at temperatures from 900 to 1400 C protected the internal fibers. When under stress, self-protection was not observed. Increasing the stress from 10 to 25 ksi caused a 67 to 82 percent reduction in times to failure at temperatures from 750 to 1500 C. Based on experimental results, observation, and theory, a finite difference model was developed, which simulates the diffusion of oxygen into a matrix crack that is bridged by carbon fibers. The model allows the influence of important variables on oxidation kinetics to be studied systematically, i.e., temperature, reaction rate constant, diffusion coefficient, environment, and sample geometry.

Halbig, Michael C.↗
DOE OSTI2026

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗
NASA NTRS1967

A Study to Analyze the Permeation of High Density Gases and Propellant Vapors Through Single Layer Teflon or Teflon Structure Materials and Laminations

This report contains the results of a fifteen month analytical and experimental study of the leakage rate of the pressurant gases (N2, He) and the propellant vapors (N2O4,N2H4) through bladder structures consisting of two layers of Teflon separated by a metallic foil diffusion barrier containing microscopic or larger holes. Results were obtained for the steady state leakage rate through circular holes and long rectangular openings in the barrier for arbitrary thicknesses of the two Teflon layers. The effect of hole shape and relative hole position on the leakage rate were studied. The transient problem was analyzed and it was shown that steady state calculations are adequate for estimating the leakage rate. A computer program entitled "Diffusion Analyzer Program" was developed to calculate the leakage rate, both transient and steady state. Finally, the analytical results were compared to experimentally determined values of the leakage rate through a model laminated bladder structure. The results of the analysis are in good agreement with experiment. The experimental effort (Part II of the Bladder Permeation Program) measured the solubility, diffusion coefficient and permeability of helium, nitrogen and nitrogen tetroxide vapor through Teflon TFE and FEP membranes. Data were obtained in the temperature range of 25 to 100 C at pressures ranging from near vacuum to about 20 atmospheres. Results of the experimental effort were compared with the limited data previously reported. As a verification to the applicability of results to actual bladder systems, counter diffusion tests were performed with a laminated sample containing aluminum foil with a selected group of holes.

Johnson, R. L.↗
DOE OSTI2024

High-entropy Li-rich layered oxide cathode for Li-ion batteries

High-entropy oxides (HEOs) are emerging as promising cathode materials for Li-ion batteries (LIBs) due to their stable solid-state phase and compositional flexibility. Herein, we investigate the structural and electrochemical properties of a novel non-equimolar high-entropy cathode material, termed high-entropy Li-rich layered oxide (HE-LLO, Li 1.15 Na 0.05 Ni 0.19 Mn 0.56 Fe 0.02 Mg 0.02 Al 0.02 O 1.97 F 0.03 ), in comparison to a pristine Li-rich layered oxide (PR-LLO, Li1.2Ni0.2Mn0.6O2). The incorporation of multiple cations (Na + , Al 3+ , Mg 2+ , Fe 3+ ) and anion (F - ) into HE-LLO introduces compositional diversity, enhancing structural stability through the entropy stabilization effect. Theoretical calculations confirm a significantly higher configurational entropy in HE-LLO compared to PRLLO, supporting its high-entropy nature. Electrochemical evaluations demonstrate that HE-LLO exhibits considerable capacity retention, preserving 76.8 % of its discharge capacity at 0.5C after 200 cycles, compared to only 36.2 % for PR-LLO. Even under high-temperature conditions, HE-LLO outperformed PR-LLO, maintaining 76.1 % of its discharge capacity after 100 cycles at 5C, while PR-LLO retained only 12.4 %. These enhancements are attributed to the improved phase reversibility and higher Li + ion diffusion coefficients of HE-LLO, validated by ex-situ characterizations using a synchrotron X-ray technique, along with density functional theory (DFT) calculations. In conclusion, these findings highlight the promise of non-equimolar HEOs as a novel design strategy for highperformance cathode materials.

25 ENERGY STORAGE↗
NASA NTRS1979

Two-stream theory of reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption parameters of the radiative transfer equation - the single scattering albedo, the optical depth, and the integrated phase function are obtained from the grain size and density of snow. Analytical expressions for the intensity within the snowpack, the reflectance, and the asymptotic flux extinction coefficient, are given. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance, and the asymptotic flux extinction coefficient. The theory may also be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗
NASA NTRS1984

Problems in simulating the stratocumulus-topped boundary layer with a third-order closure model

The Andre et al. (1976, 1978) third-order closure model, in which the time rate of change terms, the relaxation and rapid effects for pressure-related terms, and the clipping approximation are used along with the quasi-normal closure, is invoked in the study of turbulence in a cloudy layer that is radiatively cooled from above. A spurious oscillation whose greatest amplitude lies near the inversion is shown by analysis to arise from the mean gradient and buoyancy terms of the triple-moment equations. An attempt is made to damp the oscillation through the introduction of diffusion terms into the triple-moment equations. The results obtained are noted to be sensitive to the ad hoc eddy coefficient applied in the third-moment equations.

Moeng, C.-H.↗
DOE OSTI2025

Self-Trapped-Exciton Radiative Recombination in β–Ga 2 O 3 : Impact of Two Concurrent Nonradiative Auger Processes

The peculiarities of radiative and nonradiative processes associated with self-trapped intrinsic eXcitons in the excited β-Ga 2 O 3 crystals are studied via time-resolved techniques of induced absorption, transient grating, and photoluminescence (PL) at room temperature. The excitation above the bandgap is produced by laser pulses with linear light polarization parallel and orthogonal in the (–201) and (001) planes. We elucidate that the nonradiative recombination rate occurring in the eXciton prevails over its radiative emission rate in a wide range of free carrier concentration composed of excited and equilibrium electrons. Hence, the nonradiative recombination has no effect on the strong anisotropy and the shape of the eXciton emission band. However, we find out that the conventional ABC model of electron effective lifetime is insufficient for explanation of the excitation dependences. Inclusion of two nonradiative Auger mechanisms in a modified ABC formula provides excellent agreement of these dependences. We conclude that the trap-assisted Auger process is in proportion to the free electron density with coefficient B = 1.1 × 10 –11 cm 3 /s and appears at low/intermediate excitation, while the triple-particle Auger process is in proportion to Δn 2 with coefficient C = 8 × 10 –30 cm 6 /s and appears at high excitation conditions. The transition between two Auger mechanisms is accompanied by a rise of the eXciton diffusivity in preferred crystallographic directions where the radiative PL intensity is maximal. The diffusion length LD in these directions can reach values ~300 nm, but, at high excitations, LD becomes limited by Auger lifetimes. These findings pave the way for the implementation of self-trapped eXcitons into specific optoelectronic devices.

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

Transient thermal stress problem for a circumferentially cracked hollow cylinder

The transient thermal stress problem for a hollow elasticity cylinder containing an internal circumferential edge crack is considered. It is assumed that the problem is axisymmetric with regard to the crack geometry and the loading, and that the inertia effects are negligible. The problem is solved for a cylinder which is suddenly cooled from inside. First the transient temperature and stress distributions in an uncracked cylinder are calculated. By using the equal and opposite of this thermal stress as the crack surface traction in the isothermal cylinder the crack problem is then solved and the stress intensity factor is calculated. The numerical results are obtained as a function of the Fourier number tD/b(2) representing the time for various inner-to-outer radius ratios and relative crack depths, where D and b are respectively the coefficient of diffusivity and the outer radius of the cylinder.

Nied, H. F.↗
DOE OSTI2022

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗
NASA NTRS1994

Loss of solar He and Ne from IDPS in subducting sediment: Diffusion and the effect of phase changes

The results of the diffusion experiment for solar He and Ne in IDP's in a magnetic separate from Pacific Ocean sediment suggest that solar He and Ne would be easily released from IDP grains and hence lost from subducting slabs at shallow depths. However, since the diffusion experiments was conducted under high vacuum, there may be a possibility that magnetite grains, which are supposedly the main constituent of the magnetic fraction, might be partly reduced to form a metal phase due to low oxygen fugacity in the experimental condition. If this is the case, such a phase change might affect the gas release and hence the results of the diffusion coefficients. In order to examine whether or not such a phase change really occurred in the condition of the diffusion experiment, I conducted a heating experiment for a magnetic separate from Pacific Ocean sediment. In the same condition as in the diffusion experiment, and the run products were examined with an x ray diffraction method. Three samples were prepared: they were wrapped with platinum foil, put in a vacuum line, and heated in a molybdenum crucible for two hours at 500 C, 800 C,and 950 C, respectively. After cooling the furnace, the samples were taken out from the crucible and analyzed with an x ray diffraction method.

Hiyagon, H.↗
NASA NTRS1982

Ion distribution effects of turbulence on a kinetic auroral arc model

An inverted-V auroral arc structure plasma-kinetic model is extended to phenomenologically include the effects of electrostatic turbulence, with k-parallel/k-perpendicular being much less than unity. It is shown that, unless plasma sheet ions are very much more energetic than the electrons, anomalous resistivity is not a large contributor to parallel electrostatic potential drops, since the support of the observed potential drop requires a greater dissipation of energy than can be provided by the plasma sheet. Wave turbulence can, however, be present, with the ion cyclotron turbulence levels suggested by the ion resonance broadening saturation mechanism of Dum and Dupree (1970) being comparable to those observed on auroral field lines. The diffusion coefficient and net growth rate are much smaller than estimates based solely on local plasma properties.

Cornwall, J. M.↗
NASA NTRS1992

Phase-field model for isothermal phase transitions in binary alloys

A new phase field model is described which models isothermal phase transitions between ideal binary alloy solution phases. Equations are developed for the temporal and spatial variation of the phase field, which describes the identity of the phase, and of the composition. An asymptotic analysis, as the gradient energy coefficient of the phase field becomes small, was conducted. From the analysis, it is shown that the model recovers classical sharp interface models of this situation when the interfacial layers are thin, and they show how to relate the parameters appearing in the phase field model to material and growth parameters in real systems. Further, three stages of temporal evolution are identified: the first corresponding to interfacial genesis which occurs very rapidly; the second to interfacial motion controlled by the local energy difference across the interface and diffusion; the last taking place on a long time scale in which curvature effects are important and which correspond to Ostwald ripening. The results of the numerical calculations are presented.

Wheeler, A. A.↗
NASA NTRSDate not supplied

TPSAS-NF1676L-34796-DND

The recently completed 2015-2017 EVS-2 mission NAAMES (The North Atlantic Aerosols and Marine Ecosystems Study) provides a unique set of airborne remote sensing and ship-based in situ measurements in a remote and under-sampled region of the Earth?s ocean. The NASA SABOR (Ship-Aircraft Bio-Optical Research) campaign was conducted during the summer of 2014, also in the Atlantic Ocean, out of the Chesapeake Bay and in the eastern coastal region of the United States. For both campaigns, the NASA GISS Research Scanning Polarimeter, a multi-angle, multi-spectral polarimeter measured the upwelling polarized radiances from aircraft. And for both campaigns, ship-based in situ measurements of the ocean were collected jointly with the aircraft measurements. We present results from the new RSP-MAPP (1) that is based on optimal estimation and that can retrieve simultaneous aerosol microphysical properties (including effective radius, single-scattering albedo, and real refractive index) and ocean color products using accurate radiative transfer, Mie calculations, and the DP (Detritus Plankton) series of accurate, numerically efficient ocean bio-optical models. The RSP-MAPP algorithm was applied to data collected during SABOR and NAAMES to retrieve aerosol microphysics and ocean products for all Aerosols-Above-Ocean (AAO) scenes. The RSP-MAPP products are compared against i) collocated aerosol extinction and backscatter profiles collected by the NASA LaRC airborne High Spectral Resolution Lidar (HSRL-1), including lidar depth profiles of the ocean diffuse attenuation coefficient and the hemispherical backscatter coefficient, and ii) collocated ship-based in situ measurements of the ocean.

Snorre Stamnes↗
NASA NTRS1980

Analytical and experimental study of mean flow and turbulence characteristics inside the passages of an axial flow inducer

The inviscid and viscid effects existing within the passages of a three bladed axial flow inducer operating at a flow coefficient of 0.065 are investigated. The blade static pressure and blade limiting streamline angle distributions were determined and the three components of mean velocity, turbulence intensities, and turbulence stresses were measured at locations inside the inducer blade passage utilizing a rotating three sensor hotwire probe. Applicable equations were derived for the hotwire data reduction analysis and solved numerically to obtain the appropriate flow parameters. The three dimensional inviscid flow in the inducer was predicted by numerically solving the exact equations of motion, and the three dimensional viscid flow was predicted by incorporating the dominant viscous terms into the exact equations. The analytical results are compared with the experimental measurements and design values where appropriate. Radial velocities are found to be of the same order as axial velocities within the inducer passage, confirming the highly three dimensional characteristic of inducer flow. Total relative velocity distribution indicate a substantial velocity deficiency near the tip at mid-passage which expands significantly as the flow proceeds toward the inducer trailing edge. High turbulence intensities and turbulence stresses are concentrated within this core region. Considerable wake diffusion occurs immediately downstream of the inducer trailing edge to decay this loss core. Evidence of boundary layer interactions, blade blockage effects, radially inward flows, annulus wall effects, and backflows are all found to exist within the long, narrow passages of the inducer.

Gorton, C. A.↗
DOE OSTI2026

Multiscale Modeling and Experimental Insights into High-Temperature Soil Biodegradation Dynamics of Semi-Crystalline Poly(Lactic Acid) Nonwoven Fabrics

This study investigates the biodegradation of semi-crystalline poly(lactic acid) (PLA) nonwovens (NWs) in soil at 58 °C using both experimental and mathematical modeling approaches. The model utilizes a system of parabolic diffusion-reaction partial differential equations (PDEs) to elucidate chemical transformations over time and in space. It accounts for phenomena such as the diffusion of water and lactic acid monomers through the polymer matrix and into the surrounding soil, along with their microbial breakdown. It also accounts for the initial PLA crystallinity and predicts its evolution in time. The model is solved numerically for a single filament, and the results were used to shed light on PLA NW transformations observed in soil over a 180-day incubation period. Various characterization techniques, including scanning electron microscopy (SEM), differential scanning calorimetry (DSC), and Raman spectroscopy, were employed to assess morphological changes, crystallinity, and molecular changes in the NWs throughout the experiment. By comparing the experimental data with the model predictions, the hydrolysis rate coefficient was found to be 3.37 × 10 -7 s -1 , while the rate of microbial degradation of lactic acid monomers was faster, of the order of 9.63 × 10 -7 s -1 . The findings highlight the significant role of crystallinity in the biodegradation process. The PLA degradation ceases when no amorphous material remains, and the crystallinity reaches 0.8, as observed in the experiments by day 120. Furthermore, this research contributes to a deeper understanding of PLA biodegradation dynamics and offers insights for effectively managing biodegradable materials in environmental settings.

Diffusion−reaction modeling↗
DOE OSTI2023

A dynamic variational multiscale method on unstructured meshes for stationary transport problems

Here, this paper presents a variational multiscale (VMS) based finite element method where the stabilization parameter is computed dynamically. The current dynamic procedure takes in a general structure/form of the stabilization parameter with unknown coefficients and computes them dynamically in a local fashion resulting in a dynamic VMS-based finite element method. Thus, a static stabilization parameter with pre-defined coefficients is not needed. A variational Germano identity (VGI) based local procedure suitable for unstructured meshes is developed to perform the dynamic computation in a local fashion. The local VGI based procedure is applied for each interior vertex in the mesh and unknown coefficients are first determined locally at each vertex, and subsequently, for each element a maximum value is taken over the vertices of the element. To make the current procedure practical, a coarser secondary solution is constructed from the primary coarse-scale solution, which is done locally over a patch of elements around each interior vertex. Further, averaging steps are employed to make the local dynamic procedure robust. Currently, the new dynamic VMS formulation is applied to steady problems governed by the advection-diffusion and incompressible Navier-Stokes equations in both 1D and 2D to demonstrate its efficacy and effectiveness.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗