Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion 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 1,657 records · Page 92

Heat transfer, diffusion, and evaporation

Although it has long been known that the differential equations of the heat-transfer and diffusion processes are identical, application to technical problems has only recently been made. In 1916 it was shown that the speed of oxidation of the carbon in iron ore depends upon the speed with which the oxygen of the combustion air diffuses through the core of gas surrounding the carbon surface. The identity previously referred to was then used to calculate the amount of oxygen diffusing to the carbon surface on the basis of the heat transfer between the gas stream and the carbon surface. Then in 1921, H. Thoma reversed that procedure; he used diffusion experiments to determine heat-transfer coefficients. Recently Lohrisch has extended this work by experiment. A technically very important application of the identity of heat transfer and diffusion is that of the cooling tower, since in this case both processes occur simultaneously.

GASES - DIFFUSION↗

Scalar rate correlation at a turbulent liquid free surface - A two-regime correlation for high Schmidt numbers

An experimental correlation is derived for gas absorption at a turbulent, shear-free liquid interface. The correlation is expressed in terms of the liquid-side turbulence intensity, liquid-side macroscale, and the properties of the diffusing gas and solvent. The transfer coefficient increases linearly with rms velocity up to a point where the eddy Reynolds number reaches a critical (Schmidt number dependent) value. At higher velocities, there is a more rapid linear rise. The slope of the lower Reynolds number region is proportional to the square root of the diffusivity; at Reynolds numbers much higher than that of the break point, the slope becomes independent of diffusivity.

Khoo, Boo-Cheong↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical study of solar flare particle propagation in the heliosphere

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation. For low-rigidity particles the radial mean free path so obtained is significantly larger than the mean free path calculated from the scattering coefficient due to the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between the scattering mean free path and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

Mechanical and Thermophysical Properties of 3D-Printed SiC-FY20

In the Transformational Challenge Reactor (TCR), the fuel blocks consist of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which are stacked to form fuel columns. The SiC matrix is manufactured using binderjet 3D printing followed by loading the TRISO fuel particles and the chemical vapor infiltration (CVI) process. Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Mechanical and thermophysical properties were measured from various types of specimens printed for two or three orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 286–306 MPa, which are 22–27% lower than that of the CVD SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build directions of SiC samples, while thermal diffusivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. This report also includes discussions on the uniaxial tensile properties of the as-printed SiC before CVI and on ongoing efforts for irradiation effects studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effectiveness of Micro-Blowing Technique in Adverse Pressure Gradients

The impact of the micro-blowing technique (MBT) on the skin friction and total drag of a strut in a turbulent, strong adverse-pressure-gradient flow is assessed experimentally over a range of subsonic Mach numbers (0.3 less than M less than 0.7) and reduced blowing fractions (0 less than or equal to 2F/C (sub f,o) less than or equal to 1.75). The MBT-treated strut is situated along the centerline of a symmetric 2-D diffuser with a static pressure rise coefficient of 0.6. In agreement with presented theory and earlier experiments in zero-pressure-gradient flows, the effusion of blowing air reduces skin friction significantly (e.g., by 60% at reduced blowing fractions near 1.75). The total drag of the treated strut with blowing is significantly lower than that of the treated strut in the limit of zero-blowing; further, the total drag is reduced below that of the baseline (solid-plate) strut, provided that the reduced blowing fractions are sufficiently high. The micro-blowing air is, however, deficient in streamwise momentum and the blowing leads to increased boundary-layer and wake thicknesses and shape factors. Diffuser performance metrics and wake surveys are used to discuss the impact of various levels of micro-blowing on the aerodynamic blockage and loss.

Welch, Gerard E.↗

Performance Characteristics of Plane-Wall Two-Dimensional Diffusers

Experiments have been made at Stanford University to determine the performance characteristics of plane-wall, two-dimensional diffusers which were so proportioned as to insure reasonable approximation of two-dimensional flow. All of the diffusers had identical entrance cross sections and discharged directly into a large plenum chamber; the test program included wide variations of divergence angle and length. During all tests a dynamic pressure of 60 pounds per square foOt was maintained at the diffuser entrance and the boundary layer there was thin and fully turbulent. The most interesting flow characteristics observed were the occasional appearance of steady, unseparated, asymmetric flow - which was correlated with the boundary-layer coalescence - and the rapid deterioration of flow steadiness - which occurred as soon as the divergence angle for maximum static pressure recovery was exceeded. Pressure efficiency was found to be controlled almost exclusively by divergence angle, whereas static pressure recovery was markedly influenced by area ratio (or length) as well as divergence angle. Volumetric efficiency. diminished as area ratio increased, and at a greater rate with small lengths than with large ones. Large values of the static-pressure-recovery coefficient were attained only with long diffusers of large area ratio; under these conditions pressure efficiency was high and. volumetric efficiency low. Auxiliary tests with asymmetric diffusers demonstrated that longitudinal pressure gradient, rather than wall divergence angle, controlled flow separation. Others showed that the addition of even a short exit duct of uniform section augmented pressure recovery. Finally, it was found that the installation of a thin, central, longitudinal partition suppressed flow separation in short diffusers and thereby improved pressure recovery

Reid, Elliott G↗

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗

High‐Throughput Evaluation of Hardening Coefficients of Eight Alloying Elements in Magnesium

Liquid–solid diffusion couples (LSDCs) are employed to generate a composition gradient in the single‐phase hexagonal closed‐packed (hcp) solid solution with compositions up to the solubility limit of various solutes in Mg. Nanoindentation scanning across the composition gradient in LSDCs allows effective evaluation of composition‐dependent hardness of eight alloying elements (Al, Ca, Ce, Gd, Li, Sn, Y, and Zn) in the hcp Mg phase. The hardening coefficients, an indicator of the potency of solid‐solution hardening, are evaluated from the measured composition‐hardness data and correlated with various materials properties such as atomic radius, shear modulus, and elastic modulus of the solutes. The rank of hardening potency of Al, Gd, Sn, Y, and Zn measured by nanoindentation is in good agreement with that measured by microindentation reported in the literature. The hardening coefficient (potency) from the strongest to the weakest is Ce > Ca > Y ≈ Gd > Zn > Al ≈ Sn > Li in Mg‐based hcp binary solid solutions. The hardening coefficient is found to be closely correlated with the strengthening potency.

36 MATERIALS SCIENCE↗

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fractional crystallization of iron meteorites, an experimental study

Measurements of the distribution coefficients of various elements (Ni, Co, P, Pt, Ir, Ge, Cr, and Au) between solid and liquid phases were made for FeNi alloys. A high thermal gradient crystal growing furnace was employed in order to promote plane front solidification, rapid mixing in the liquid and a minimum of diffusion in the solid. Eight measurements of the nickel distribution coefficient gave an average value of 0.91 plus or minus 0.03 and this value was independent of the Ni concentration in the range 5.5 to 8.75 wt. % Ni. Equilibrium distribution coefficients were also calculated by correcting for the diffusion boundary layer present in the liquid. These values are within 20% of the measured distribution coefficient values. Except for Co and Ni, all the measured distribution coefficient values are closer to 1.0 than those calculated by Scott (1972) for the IIIAB irons. Using the experimental values and the calculated equilibrium values for ternary additions to FeNi melts, it appears that fractional crystallization processes will not explain the observed chemical variations in the major iron meteorite groups. Alternative explanations are suggested which also recognize the fact that cooling rates vary widely in several chemical groups.

Goldstein, J. I.↗

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stable Rate of Slip Along the Karakax Section of the Altyn Tagh Fault from Observation of Interglacial and Postglacial Offset Morphology and Surface Dating

Digital elevation maps obtained using TanDEM-X and Pleiades data combined with newly obtained surface age estimates using cosmogenic radionuclide (CRN) and optically simulated luminescence (OSL) methods are used to quantify the slip rate along the western section of the Altyn Tagh Fault in southern Xinjiang. The reconstruction of the conical shape of massive alluvial fans inferred to be from the Eemian (115±7 ka) from CRN dating shows consistent left-lateral offsets of 300±20 m, yielding a slip rate of 2.6±0.3 mm/year. Successive episodes of incision have left cut terraces inset in wide canyons, 10–25 m below the fans' surface. The incision was followed by the deposition of a broad terrace of early Holocene age, which is reincised by modern stream channels. Near the village of Shanxili, a 200 m wide valley is partially dammed by a shutter ridge displaced by the fault. A fill terrace deposited upstream from the ridge has an OSL age of 8.8±0.6 ka. The 23±2 m offset of the riser incising the terrace indicates a minimum postdepositional movement on the fault, yielding a Holocene rate of 2.6±0.5 mm/year, consistent with the 115-ka average slip rate. Furthermore, scarp degradation analysis using mass diffusion reveals a nonlinear relationship between fault displacement and degradation coefficient along the progressively exposed fault scarp, a pattern suggesting either seismic clustering or variable diffusion rate since the Eemian. Together with the Gozha Co-Longmu Co fault to the south, the Karakax section of the Altyn Tagh Fault contributes to the eastward movement of the western corner of Tibet.

58 GEOSCIENCES↗

Small to Full-Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

Crystalline silicotitanate (CST) ion exchanger is planned to be used to remove cesium (137Cs) from the aqueous phases of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. Column scale up testing was conducted to evaluate performance of Cs exchange onto the CST. Batch contact testing was conducted to assess exchange kinetics of Cs exchange at four different CST PSDs. Column testing was conducted at three column sizes, small (2.5% full height), medium (12% full height), and full height, to assess Cs load performance behavior. Testing at the small scale was compared to actual Hanford tank waste testing to verify the validity of simulant tests to accurately represent full height column performance. A change in CST particle size was essential in scaling the small column dynamics up to full scale. Based on these results, a determination of intraparticle and film diffusion impacts on overall mass transfer coefficients will allow future modeling of the breakthrough performance at a range of process conditions.

Westesen, Amy M.↗

Measurements of He-Ar, He-N 2 , and He-air thermal creep slip and accommodation coefficients at high temperatures

Understanding gas flows in capillaries has many applications in modeling the transport of gases in nano-structured, porous, or fractured media. A network of capillaries can often approximate these media, and also information on gas-surface interactions obtained from capillary experiments can be used for modeling flows in these media. Experimental data on flows of different mixtures of He-Ar, He-N 2 , and He-air in the slip regime and in stainless steel capillaries at high temperatures were obtained by using a two-bulb apparatus. An accurate expression for the thermal creep slip coefficients with the Lennard-Jones potential parameters and diffuse-specular reflection gas-surface interaction conditions were then used to obtain the accommodation coefficients for the different gases. The experimental data are best described with values of accommodation coefficients in the range of 0.1–0.3 for He and 0.5–1.0 for Ar, N 2 , and air. Further, the use of values in this range is suggested for modeling gaseous flows in capillaries and nano-structured, porous, or fractured media if other direct measured values for a particular medium are unavailable.

42 ENGINEERING↗

Mechanical and Thermophysical Properties of 3D-Printed SiC before and after Neutron Irradiation – FY21

This report presents the mechanical and thermophysical properties of 3D-printed SiC before and after neutron irradiation that have been evaluated to assess the fuel matrix material for the Transformational Challenge Reactor (TCR). The TCR fuel form consists of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which is manufactured through a newly developed processing route combining binderjet 3D printing, TRISO fuel particle loading, and chemical vapor infiltration (CVI). Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Baseline mechanical and thermophysical properties were measured from the disk specimens printed for different and sizes orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 280–310 MPa, which are slightly lower than that of the chemical vapor deposition (CVD) SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build direction of SiC samples, while thermal conductivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. Neutron irradiation tests were carried out on the 3D-printed 6-mm diameter SiC disk specimens. Irradiation was carried to 2.3 dpa over a temperature range of 360–880°C. No significant degradation in strength was observed in SiC after irradiations in various conditions and with different orientations. Anisotropy that had been observed in the thermal conductivity of 3D-printed SiC prior to irradiation vanished after irradiation as the irradiation defect thermal resistivity accumulated in the material. Electron microscopy of the microstructure after neutron irradiation showed distinct defect morphologies in the heterogenous material, but no evidence for irradiation-induced cracking or degradation in the microstructure was observed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Creating Novel Gas Separation Constructs by Manipulating Free-Volume Distributions in Polymer-Grafted Nanoparticles

Polymer membranes are critical to sustainability applications, and in this work we focused on one key application, the efficient separation of gas mixtures. Despite their widespread use, important challenges remain in the use of polymer membranes in this context – i.e., finding membrane materials with selectively improved transport of the desired mixture component(s) while possessing enhanced strength and improved aging behavior relative to the best current benchmarks. The important separation figures of merit are the gas flux, which is proportional to the permeability, P i = D i × S i (Di is the gas diffusivity and S i its solubility coefficient) and selectivity (i.e., gas purity,α ij = P i /P j ). The goal is to simultaneously maximize P i and α ij . Most research to date has empirically targeted the development of new glassy polymers with improved separation performance. These include thermally rearranged (TR) polymers and polymers of intrinsic microporosity (PIMs).

36 MATERIALS SCIENCE↗