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 541 records · Page 30

Current understanding of point defects and diffusion processes in silicon

The effects of oxidation of Si which established that vacancies (V) and Si self interstitials (I) coexist in Si at high temperatures under thermal equilibrium and oxidizing conditions are discussed. Some essential points associated with Au diffusion in Si are then discussed. Analysis of Au diffusion results allowed a determination of the I component and an estimate of the V component of the Si self diffusion coefficient. A discussion of theories on high concentration P diffusion into Si is then presented. Although presently there still is no theory that is completely satisfactory, significant progresses are recently made in treating some essential aspects of this subject.

Tan, T. Y.↗

A Semianalytical Ocean Color Inversion Algorithm with Explicit Water Column Depth and Substrate Reflectance Parameterization

A semianalytical ocean color inversion algorithm was developed for improving retrievals of inherent optical properties (IOPs) in optically shallow waters. In clear, geometrically shallow waters, light reflected off the seafloor can contribute to the water-leaving radiance signal. This can have a confounding effect on ocean color algorithms developed for optically deep waters, leading to an overestimation of IOPs. The algorithm described here, the Shallow Water Inversion Model (SWIM), uses pre-existing knowledge of bathymetry and benthic substrate brightness to account for optically shallow effects. SWIM was incorporated into the NASA Ocean Biology Processing Group's L2GEN code and tested in waters of the Great Barrier Reef, Australia, using the Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua time series (2002-2013). SWIM-derived values of the total non-water absorption coefficient at 443 nm, at(443), the particulate backscattering coefficient at 443 nm, bbp(443), and the diffuse attenuation coefficient at 488 nm, Kd(488), were compared with values derived using the Generalized Inherent Optical Properties algorithm (GIOP) and the Quasi-Analytical Algorithm (QAA). The results indicated that in clear, optically shallow waters SWIM-derived values of at(443), bbp(443), and Kd(443) were realistically lower than values derived using GIOP and QAA, in agreement with radiative transfer modeling. This signified that the benthic reflectance correction was performing as expected. However, in more optically complex waters, SWIM had difficulty converging to a solution, a likely consequence of internal IOP parameterizations. Whilst a comprehensive study of the SWIM algorithm's behavior was conducted, further work is needed to validate the algorithm using in situ data.

Great Barrier Reef↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Studies of extended planetary atmospheres

There was a theoretical study of physical and chemical processes in the stratosphere, later broadened to include the mesosphere. Particular emphasis was laid on testing of proposed height profiles of the eddy diffusion coefficient against observed tracer data. Eventually the effort shifted to study of ozone time series in satellite data, and interpretation in terms of aeronomical processes. Since all this work is computer-intensive, the first year of funding also contributed to the acquisition of a powerful minicomputer system, in collaboration with several other faculty members. This proved to be highly successful and cost effective.

Hunten, D. M.↗

Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure

The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.

Jameson, Antony↗

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 and design of hydrogen/oxygen chemical rocket engines. The complete multi-species, chemically reacting and diffusing Navier-Stokes equations are modelled, finite difference approach that is tailored to be conservative in an axisymmetric coordinate system for both the inviscid and viscous terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and 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 when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent plug-and-spool rocket engine analysis cases performed. Further, the Soret term was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration cooled rocket engine.

Kacynski, Kenneth J.↗

Solvation dynamics of choline fluoride in ethylene glycol – Water mixtures

Due to their favorable physical properties, deep eutectic solvents (DESs) have been of interest in many applications in the energy and industry sectors. In this study, we investigate the new deep eutectic solvent ethalineF (a 1:2 M ratio of choline fluoride:ethylene glycol) and the effect of water addition on the solvation dynamics compared to the better known ethaline (a 1:2 M ratio of choline chloride:ethylene glycol). Femtosecond transient absorption spectroscopy and NMR diffusometry were used to study the solvation dynamics. 19 F NMR spectra were obtained for ethalineF solutions with different water additions. The 19 F NMR spectra show two fluorine signals; one can be assigned to fluoride anions solvated by EG and the other signal to fluoride anions associated with the much heavier and less mobile choline cation. Density, viscosity, conductivity, and E T (30) polarity measurements are also provided for ethalineF. The nuclear spin nature of fluorine atom permits to measure the self-diffusion coefficient of the negative ion in addition of the cation diffusion in a choline halide based deep eutectic solvent system. This allows to calculate the NMR-predicted conductivities σ NMR . The obtained ion dissociation coefficients α D at different water concentrations suggest a partial ion dissociation which was commensurate with the relatively low ion conductivity. Furthermore, the observed modest conductivity of ethalineF is explained in terms of ion association as well as the gradual increase in conductivity caused by the addition of water.

19F NMR↗

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Saatchi, Sasan S.↗

Calculation of Propulsive Nozzle Flowfields in Multidiffusing Chemically Reacting Environments

An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. 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 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Kacynski, Kenneth John↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

Electron Injection via Modified Diffusive Shock Acceleration in High-Mach-number Collisionless Shocks

The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stress corrosion cracking of a titanium alloy in chloride-containing liquid environments

Ti-8 Al-1 Mo-1 V was tested in solvents containing chloride, bromide or iodide ions. It was found that the region II plateau velocity could be correlated with the inverse of the viscosity for several solvents. It seems most likely that this reflects the influence of the viscosity on the diffusion coefficient of some species involved in a mass transfer process which limits the crack growth. Associated tests were performed in which the potential, pH and concentration of Cl- were varied, in some of the same solutions as used for the viscosity tests. It is concluded from these tests that the absence of cathodic protection against cracking in solutions of low pH arises from an ohmic effect which isolates the crack tip from external potential control.

Smyrl, W. H.↗

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Knudsen diffusion - The effect of small pore size and low gas pressure on gaseous transport in soil

The analytical principles and applications of the theory of Knudsen diffusion are reviewed, with emphasis on gas transport in the soils of planetary bodies. Knudsen diffusion occurs when the mean free path of diffusing gas molecules surpasses the size of the pores through which diffusion proceeds. The process is then dominated by collisions with the pore walls. Computational techniques for deriving the Knudsen coefficient for soils with a nonreentrant cross-section shape are reviewed, along with methods of deriving a coefficient for soils which permit both Knudsen and bulk diffusion. Sample calculations for three pore-size distributions are provided to illustrate the decrease in transport efficiency with increasingly smaller soil pore sizes.

Clifford, S. M.↗

Evaluating diffusion and the thermodynamic factor for binary ionic mixtures

Molecular dynamics (MD) simulations are a powerful tool for the calculation of transport properties in mixtures. Not only are MD simulations capable of treating multicomponent systems, they are also applicable over a wide range of temperatures and densities. In plasma physics, this is particularly important for applications such as inertial confinement fusion. While many studies have focused on the effect of plasma coupling on transport properties, here we focus on the effects of mixing. We compute the thermodynamic factor, a measure of ideal/non-ideal mixing, for three binary ionic mixtures. Here, we consider mixtures of hydrogen and carbon, hydrogen and argon, and argon and carbon, each at 500 randomly generated state points in the warm dense matter and plasma regimes. The calculated thermodynamic factors indicate different mixing behavior across phase space, which can significantly affect the corresponding mutual diffusion coefficients. As MD simulations are still computationally expensive, we apply modern data science tools to predict the thermodynamic factor over a large phase space. Further, we propose a more accurate approximation to the mutual diffusion coefficient than the commonly applied Darken relation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗