Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport 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.

At least 703 records · Page 39

A two-dimensional model simulation of the EL Chichon volcanic eruption cloud

Using a two-dimensional model of sulfate photochemistry, transport, and aerosol microphysics, a 2-year period following the eruptions of El Chichon have been simulated. Present calculations suggest that the residence time of the cloud in the stratosphere exceeds 2 years. The model reproduces the observed optical depth, lidar backscatter, and infrared extinction coefficients, if about 10 megatonnes of SO2 are injected. The major deficiency of the model is a somewhat too rapid transport.

Capone, L. A.↗

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗

Role of Eu-Doping in the Electron Transport Behavior in the Zintl Thermoelectric Ca 5-x-y Yb x Eu y Al 2 Sb 6 System

A series of Eu-doped Zintl compounds belonging to theCa 5-x-y Yb x Eu y Al 2 Sb 6 (x = 0, 1.12; 0 ≤ y ≤ 0.63(2)) system have been successfully synthesized by both the arc-melting and the molten Pb-flux methods. All of the five title compounds initially crystallized in the Ca 5 Ga 2 As 6 -type phase (space group Pbam, Z = 2, Pearson code oP26) and maintained their original structure even after the post-heat treatment, unlike the recently reported n-type Zintl analogues in the Ca 5-x-y Yb x RE y Al 2 Sb 6 (RE = Pr, Nd, Sm) systems, which underwent a phase transition from the Ca 5 Ga 2 As 6 -type to the Ca 5 Al 2 Bi 6 -type phase after annealing. This research aimed to understand the origin of the structural preference of the title Ca 5-x-y Yb x Eu y Al 2 Sb 6 system, whether it was affected by the valence electron count or the cationic size. Electrical transport property measurements showed an increase in electrical conductivities and a decrease of Seebeck coefficients for Ca 4.89(1) Eu 0.11 Al 2 Sb 6 , Ca 4.82(1) Eu 0.18 Al 2 Sb 6 , and Ca 4.62(1) Eu 0.38 Al 2 Sb 6 , compared to the parental compound Ca 5 Al 2 Sb 6 . Hole effect measurements proved that these changes should be attributed to the reduced carrier concentration and enhanced carrier mobility. The comprehensive density functional theory calculations including electron density map analysis for the hypothetical model Ca 4.5 Eu 0.5 Al 2 Sb 6 revealed that the polarity between Al and Sb forming the anionic frameworks decreased as the Eu-dopants were introduced, which eventually affected the carrier mobility in the anionic frameworks. Thermal conductivity measurements proved that the Eu-doping successfully lowered the lattice thermal conductivity because of the enhanced atomic disordering. In conclusion, the magnetization measurements for Ca 4.37(2) Eu 0.63 Al 2 Sb 6 showed a typical Curie–Weiss behavior with weak antiferromagnetic nearest-neighbor interactions with θ p = -5.07 K.

36 MATERIALS SCIENCE↗

Encapsulation of Ba in InSb Framework Introduces Chirality in Clathrate-Like BaIn 4 Sb 4

The discovery of Zintl compounds remains a powerful strategy for identifying materials with tunable electronic and thermal transport properties. During a concerted search for new inorganic clathrates with In–Sb frameworks, we discovered BaIn 4 Sb 4 . The composition of this phase deviates from that expected for a type-I clathrate with tetrahedral coordination of all In and Sb atoms (Ba 8 In 31 Sb 15 ). Instead, in the chiral structure of BaIn 4 Sb 4 (space group P 3 1 21, No. 152), a part of the In atoms have a trigonal planar coordination of 1In + 2Sb, forming Sb 2 –In–In–Sb 2 nonplanar fragments isostructural to diborane(4) B 2 H 4 with D 2 d symmetry. Ba atoms are located inside 16-vertex In 8 Sb 8 polyhedra, which share vertices and edges to form a chiral framework around the 3 1 screw axes. The title compound is electron-balanced, [Ba 2+ ][In 2+ ] 2 [In 3+ ] 2 [Sb3–] 4 , which was confirmed by characterization of the charge and heat transport properties. BaIn 4 Sb 4 exhibits a low thermal conductivity and high Seebeck coefficient, suggesting its untapped potential for thermoelectric applications. Finally, density functional theory (DFT) calculations indicate that chemical doping may enhance carrier concentration and improve the originally low electrical conductivity, thus enhancing thermoelectric performance.

electrical properties↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DRACO: An Overview [Slides]

DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Transport of N2O in the stratosphere related to the equatorial semiannual oscillation

The double-peak structure observed in the distribution of N2O during the spring equinoctial condition is studied by means of a 2D time-dependent tracer model. The meridional transport circulation is determined from computed diabatic heating rates, and the horizontal eddy diffusion coefficients are computed as residuals by potential vorticity conservation, using observed mean winds to compute the budget of zonal mean potential vorticity. This model is used to predict the 30-d evolution of the observed N2O field. The integrated results are generally in good agreement with observed tracer fields up to the midstratosphere. It is proposed that standard climatologies underestimate the strength of the equatorial semiannual oscillation in the zonal mean wind and temperature fields.

Choi, Woo Kap↗

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Two refinements to the quasi-linear theory of ion radial diffusion are proposed and examined analytically with simulations of particle trajectories. The resonance-broadening correction by Dungey (1965) is applied to the quasi-linear diffusion theory by Faelthammar (1965) for an individual model storm. Quasi-linear theory is then applied to the mean diffusion coefficients resulting from simulations of particle trajectories in 20 model storms. The correction for drift-resonance broadening results in quasi-linear diffusion coefficients with discrepancies from the corresponding simulated values that are reduced by a factor of about 3. Further reductions in the discrepancies are noted following the averaging of the quasi-linear diffusion coefficients, the simulated coefficients, and the resonance-broadened coefficients for the 20 storms. Quasi-linear theory provides good descriptions of particle transport for a single storm but performs even better in conjunction with the present ensemble-averaging.

Chen, Margaret W.↗

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

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

Mayr, H. G.↗

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of gravity wave generated advection and diffusion in transport of tracers in the mesosphere

The roles of advection and turbulent diffusion in controlling tracer transport in the meridional plane of the mesosphere are compared using a version of the global primitive equation model for the middle atmosphere developed by Holton and Wehrbein (1980). The results of model sensitivity tests demonstrated that transport in the mesosphere for long-lived species with large-scale heights, such as H2O, CO, and NO(x), is primarily an advective process; on a global scale, small-scale turbulent diffusion produced by gravity wave breaking plays a secondary role. It is concluded that, for such tracers, the eddy diffusion coefficient used in one-dimensional chemical models of the mesosphere should be regarded as a parameterization of the transport by the meridional circulation, not as a parameterization of mixing by small-scale diffusion.

Holton, James R.↗

Membrane Transport Phenomena (MTP)

The activities during the fourth semi-annual period of the MTP project have involved the completion of the Science Concept Review (SCR) presentation and peer review, continuation of analyses for the mass transfer coefficients measured from MTA experiment data, and development of the second generation (MTP-II) instrument. The SCR panel members were generated several recommendations for the MTP project recommendations are : Table 1 Summary of Primary SCR Panel Recommendations (1) Continue and refine development of mass transfer coefficient analyses (2) Refine and upgrade analytical modeling associated with the MTP experiment. (3) Increase resolution of measurements in proximity of the membrane interface. (4) Shift emphasis to measurement of coupled transport effects (i.e., development of MTP phase II experiment concept).

Mason, Larry W.↗

Studies of the inner shelf and coastal sedimentation environment of the Beaufort Sea from ERTS-1

The author has identified the following significant results. The particulate transport processes involved in the movement of surficial waters were examined using secchi disc readings, light attenuation coefficients, and particulate weights from filtration. Observations gathered during the summers of 1971 and 1972 indicate a remarkable difference in particulate matter and turbidity between the two years. ERTS-1 imagery during August 1972 showed turbid water along the northern Alaska coast. The uniformity of distribution of the turbid water and the fact that the river discharge is low at this period suggest that the turbidity is related to causes other than river effluent. Studies indicate that wave action is a more significant factor influencing particulate transport than believed heretofore. The boundary between the essentially immobile shorefast ice and the moving pack ice has been plotted from several ERTS-1 images and found to occur fairly consistently along the 20 meter contour. Considering the vast difference in the amount of ice movement shoreward and seaward of this boundary, ice-bottom action should also be different on either side of this boundary and for that matter at the shear zone that develops along the boundary.

Reimnitz, E.↗

Dynamics of the middle atmosphere in winter (Dynamics). Interrelation between the different variations of turbulent diffusion and ionospheric absorption originating in the middle atmosphere

The turbulent diffusion coefficient was computed from the parameters of sporadic E layers using the wind shear theory of midlatitude sporadic E and models of the ionosphere as well as that of the neutral upper atmosphere. The turbulent diffusion coefficient obtained for the period of circulation disturbances associated with stratospheric warmings and for the intervals of the winter anomaly indicate changes similar to the ionospheric absorption of radio waves, in the former case decreased, in the latter case increased values. This may hit at the role of turbulent transport in the formation of these anomalies. On the basis of these findings, a seasonal variation of the turbulent diffusion coefficient having a minimum in summer and an increase of this parameter with increasing geomagnetic activity are anticipated.

Bencze, P.↗

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗