SEARCH · Engineering Papers
Results for “diffusion coefficients”
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.
Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores
Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.
The Martian climate and energy balance models with CO2/H2O atmospheres
The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.
Variations in eddy diffusion and associated transport processes
Variation in the eddy diffusion coefficient constitutes a momentum source. This variation near the turbopause affects both the composition and the wind and temperature fields of the thermosphere. Two types of calculations are carried out, for horizontal winds considered self-consistently in a three-dimensional dynamic model and horizontal winds artificially suppressed to simulate a one-dimensional model. Momentum transfer associated with variations in eddy diffusion coefficients represents the source function in the analysis presented. It is shown that for long-term variations in the thermosphere, the compositional dependence on variations in eddy diffusion coefficients can be entirely different in one and three dimensional models with horizontal winds included. Horizontal winds can greatly reduce the effects of eddy diffusion in general and the ratio between He and O amplitudes in particular.
Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo
Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.
Electrochemical measurement and analysis of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 in molten eutectic LiCl-KCl
For this work, electrochemical measurements of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 were obtained in molten LiCl-KCl eutectic at 773 K using a three-electrode configuration. Tungsten rods were used as the working and counter electrodes. Two quasi- and one reference electrode(s) (glassy carbon, nickel wire, and Ni/NiO) were used to collect electrochemical measurements of the analyte salts. Cyclic voltammetry data were analyzed to calculate values of diffusion coefficient, exchange current density, and charge transfer resistance and were determined to be on the order of 10 –5 cm 2 s –1 , 10 –2 A cm –2 , and 10 0 Ω; respectively. In general, the values of diffusion coefficients were found to be consistent with those reported of high-temperature molten salts. Relatively large values of exchange current density corresponded with smaller values of charge transfer resistance. These values were found to be reasonable in comparison to results available in the literature. The measured cyclic voltammograms were normalized with respect to both lithium reduction and chloride ion oxidation potentials. Such a normalization technique is effective for comparing experimentally obtained cyclic voltammetry data to those that have been published. An analysis of experimentally obtained results indicates the relation of electrode reactions (e.g., differences of reduction potentials) are independent of the choice of reference electrode. Additionally, the choice of reference electrode did not affect the electrochemical window, exchange current density, or diffusion coefficient values.
In Situ Electrochemical Study of the Coexistence of Eu 3+ and Eu 2+ in Molten LiCl-KCl by Rotating Disc Electrode
In this study, a graphite rotating disc electrode was successfully designed and applied to conduct electrochemical measurements under flow conditions in molten LiCl-KCl at 773 K. The concentration and diffusion coefficients of Eu 3+ and Eu 2+ in molten salt were quantitatively measured with this electrode by applying in situ cyclic voltammetry and chronoamperometry techniques. The initial concentration ratio of Eu 3+ /Eu 2+ in LiCl-KCl-0.5 wt% EuCl 3 was calculated to be 2.37. The calculated diffusion coefficients of Eu 3+ and Eu 2+ in the salt were 9.31 ± 0.06 × 10 -5 cm 2 s -1 and 9.65 ± 0.6 × 10 -5 cm 2 s -1 , respectively. A decrease of the diffusion coefficient of Eu 3+ and Eu 2+ was observed at a higher concentration of EuCl 3 , which implies the enhanced Eu-Eu ion interaction. This electrode design is expected to be utilized for the concentration measurement of multivalent ions in a flowing molten salt.
Constraints on gravity wave induced diffusion in the middle atmosphere
A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H2O and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is, at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.
Constraints on gravity wave induced diffusion in the middle atmosphere
A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H20 and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.
External field effects on diffusion and solidification derived from the free-volume model
Expressions for the diffusion coefficient and the solidification rate from the free-volume model of liquids developed by Turnbull and Cohen have been used to estimate the effects which microgravity and magnetic fields will have on these quantities. The mathematical formalism describing changes of the diffusion coefficient and the solidification rate is the same for both the microgravity and magnetic field cases, but the difference between the magnitudes of the two effects is quite large. The change in the two parameters is found to be less than .0001% for the microgravity case and on the order of 0.1 to 1.1% for the magnetic field case for four representative materials. The diffusion coefficient and the solidification rate are found to increase under the influence of an applied magnetic field, and this is in agreement with experimental observations.
Analysis of redox additive-based overcharge protection for rechargeable lithium batteries
The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection, has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. Digital simulation of the overcharge experiment leads to numerical representation of the potential transients, and estimate of the influence of diffusion coefficient and interelectrode distance on the transient attainment of the steady state during overcharge. The model has been experimentally verified using 1,1-prime-dimethyl ferrocene as a redox additive. The analysis of the experimental results in terms of the theory allows the calculation of the diffusion coefficient and the formal potential of the redox couple. The model and the theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.
Analysis of Jovian low frequency radio emissions
The density of ions in the Io plasma torus and the scattering of these ions by low frequency electromagnetic emissions detected by Voyager 1 were studied. The ion density profile was investigated using whistler dispersion measurements provided by the Voyager plasma instrument. The scale height and absolute density of H+ ions in the vicinity of the plasma torus were determined by combining the measured plasma densities with the whistler dispersion measurements. A theoretical analysis of the modes of propagation of low frequency electromagnetic emissions in the torus was undertaken. Polarization reversal effects and rough estimates of the ion diffusion coefficient were utilized. Numerical evaluation of the ion diffusion coefficients in the torus were made using the observed Voyager 1 wave intensities. Results show that the observed wave intensities produce significant ion diffusion effects in the ion torus.
Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials
Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.
Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC
In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model
Chemical zoning and diffusion of Ca, Al, Mn, and Cr in olivine of springwater pallasite
The pallasites, consisting mainly of Fe-Ni metal and olivine, are thought to represent the interior of a planetary body which slowly cooled from high temperature. Although the olivines are nearly homogeneous, ion microprobe studies revealed variations of Ca, Ti, Co, Cr, and Ni near grain edges. These variations were thought to represent diffusion in response to falling temperature of the parent body. Pallasite cooling rates have been estimated based on kamacite taenite textures but results differ by x100. In principle elemental profiles in olivine can allow estimates of cooling rate if diffusion coefficients are known; in addition, given a cooling rate, diffusion coefficients could be derived. Data are presented which show that apparent diffusion profiles can be measured for Al, Ca, Cr, and Mn which qualitatively agree with expected diffusion rates and have the potential of providing independent estimates of pallasite cooling rates.
Comparison of GEOS-5 AGCM Planetary Boundary Layer Depths Computed with Various Definitions
Accurate models of planetary boundary layer (PBL) processes are important for forecasting weather and climate. The present study compares seven methods of calculating PBL depth in the GEOS-5 atmospheric general circulation model (AGCM) over land. These methods depend on the eddy diffusion coefficients, bulk and local Richardson numbers, and the turbulent kinetic energy. The computed PBL depths are aggregated to the Koppen climate classes, and some limited comparisons are made using radiosonde profiles. Most methods produce similar midday PBL depths, although in the warm, moist climate classes, the bulk Richardson number method gives midday results that are lower than those given by the eddy diffusion coefficient methods. Additional analysis revealed that methods sensitive to turbulence driven by radiative cooling produce greater PBL depths, this effect being most significant during the evening transition. Nocturnal PBLs based on Richardson number are generally shallower than eddy diffusion coefficient based estimates. The bulk Richardson number estimate is recommended as the PBL height to inform the choice of the turbulent length scale, based on the similarity to other methods during the day, and the improved nighttime behavior.
Thermophysical properties of simple liquid metals: A brief review of theory
In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.
Simulation of pitch angle diffusion of charged particles in a disordered magnetic field
Results are reported for computer simulation experiments in which a statistical ensemble of random magnetic field realizations is generated, orbits of charged particles in the random fields are followed, and a pitch-angle diffusion coefficient is derived from the temporal evolution of the orbits. Diffusion coefficients predicted by three nonlinear theories are compared with the derived coefficients for the standard quasilinear theory of velocity diffusion, and the goals of future simulations are outlined.