Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass diffusivity”

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 127 records · Page 7

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↗

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Impact of binary interactions on the diffuse supernova neutrino background

Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%–20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. As a result, the enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomistic processes of surface-diffusion-induced abnormal softening in nanoscale metallic crystals

Abstract Ultrahigh surface-to-volume ratio in nanoscale materials, could dramatically facilitate mass transport, leading to surface-mediated diffusion similar to Coble-type creep in polycrystalline materials. Unfortunately, the Coble creep is just a conceptual model, and the associated physical mechanisms of mass transport have never been revealed at atomic scale. Akin to the ambiguities in Coble creep, atomic surface diffusion in nanoscale crystals remains largely unclear, especially when mediating yielding and plastic flow. Here, by using in situ nanomechanical testing under high-resolution transmission electron microscope, we find that the diffusion-assisted dislocation nucleation induces the transition from a normal to an inverse Hall-Petch-like relation of the strength-size dependence and the surface-creep leads to the abnormal softening in flow stress with the reduction in size of nanoscale silver, contrary to the classical “alternating dislocation starvation” behavior in nanoscale platinum. This work provides insights into the atomic-scale mechanisms of diffusion-mediated deformation in nanoscale materials, and impact on the design for ultrasmall-sized nanomechanical devices.

36 MATERIALS SCIENCE↗

Single-Molecule Tracking Measurements Reveal the Detailed Mechanisms of Molecular Diffusion in Solvent Mixtures under Nanoconfinement

Understanding mass transport mechanisms in nanopores is important for developing advanced materials for chemical separations, chemical sensing, and energy storage. This paper reports a novel imaging platform that is employed for the first time to investigate the detailed diffusion dynamics of single rhodamine B (RhB) dye molecules confined within solution-filled cylindrical anodic aluminum oxide (AAO) nanopores. The imaging platform relies on illumination of horizontally-oriented AAO nanopores in a highly inclined and laminated optical (HILO) light sheet microscopy geometry. The method was used to investigate the translational and orientational dynamics of single rhodamine B (RhB) molecules within horizontally-oriented 5- and 10-nm diameter AAO nanopores filled with water–ethanol mixtures. The established platform enabled the observation of one-dimensional motion along the pore axis involving occasional short- or long-term immobilization at the single-molecule level. Analysis of cumulative squared-displacement distributions revealed fast (5 – 30 µm²/s), intermediate (1 – 5 µm²/s), and slow (< 1 µm²/s) diffusion components. From the effects of mixture composition and pore size on the contributions of these three components, we inferred that the fast, intermediate, and slow components could be assigned to desorption-mediated hopping, crawling, and wiggling motions, respectively. The platform based on the horizontally-oriented AAO nanopores also permitted single-molecule emission polarization measurements that revealed the negligible steric confinement of individual diffusing RhB molecules within the AAO nanopores. The imaging platform based on AAO membranes and HILO microscopy provided a unique means to investigate how solvation-mediated surface interactions and nanoconfinement govern molecular transport in nanoporous environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tritium Transport Phenomena in Molten-Salt Reactors: Molten Salt Tritium Transport Experiment Design

Tritium is produced from neutron interactions with both lithium and beryllium. Large quantities of tritium are generated in Molten Salt Reactors (MSRs) which use LiF/BeF2 (FLiBe) as the fuel salt. Tritium is unique among the radionuclide hazards as it readily permeates through metal structural materials at high temperatures. All metal surfaces are potential release paths for tritium. For adequate safety analysis and eventual licensing of new reactors, predictive models for tritium transport and release from MSRs must be developed. These models must account for the multiple transport phenomena involved with tritium: fuel salt phase mass transport, dissociation/recombination reactions on metal surfaces, interstitial diffusion through the metal structure, and salt or gas phase mass transport in the downstream fluid. These models also must also be validated with representative experiments. Our previous report outlined tritium transport phenomena involved in MSRs, made suggestions on gaps in the transport dataset, and proposed an experimental test stand to test combined transport effects – tritium transport through pipe walls in a convective salt flow. In this report, we summarize an updated analysis framework for tritium transport in MSRs, report our results on hydrogen and deuterium permeation through Hastelloy N, and describe the final design of the Molten Salt Tritium Transport Experiment (MSTTE, pronounced “misty”). The MSRE provides the only wholistic experimental data set for tritium transport in MSRs and understanding the transport phenomena involved in the MSRE is crucial for future model development. One set of parameters in our analysis framework was unknown for the MSRE—surface reaction rates for tritium on Hastelloy N. This warranted our hydrogen and deuterium permeation campaign to assess the permeability, diffusivity, and solubility of hydrogen isotopes in clean Hastelloy N. Surface reaction rate constants were probed by low pressure measurements, however, no surface effects were observed in the limits of our permeation apparatus. Permeation experiments on oxidized Hastelloy N were not performed for this report but are planned in future work. The experimental test stand, MSTTE, measures combined transport properties of the salt-metal system. MSTTE is a forced convection FLiBe loop with custom designed test section to measure tritium transport through candidate structural materials. We use MSRE relevant dimensionless numbers to design and scale the test section. Hastelloy N is a candidate loop and test section material due to the relevance for the MSRE and related designs, however, other metals are being considered (e.g. 316H SS) which may better align with current vendor concepts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Lithium‐Containing Hybrid SEI Layer Enabling High Mass Loading and Anode‐Less Sodium Metal Batteries

The continuous rupturing and rebuilding of unstable solid electrolyte interphase (SEI) layer during cycling would block Na + diffusion and induce Na dendrite formation, ultimately limiting the practical application of high-energy-density sodium metal batteries. Herein, a hybrid SEI layer containing Li-species is dexterously constructed on the surface of sodium metal anode. Li-containing inorganic components (Li 3 N, LiF, and Li 2 CO 3 ) are introduced to stabilize the Na/electrolyte interface and enhance the mechanical and diffusion kinetic properties of the SEI layer, which can reduce the side reactions and gas generation, regulate Na + flux during cycling and promote rapid Na + migration for uniform dendrite-free Na deposition. As a result, the constructed Na symmetric cells achieve low overpotential and long cycle life of 5900, 1800, and 500 h at current densities of 3, 10, and 30 mA cm −2 , respectively. Furthermore, the full cells paired with the Na₃V₂(PO₄)₃ cathode demonstrate high specific capacity and excellent cycle stability, even at an ultra-high cathode loading of 39.3 mg cm −2 and a low N/P ratio (negative/positive electrode capacity ratio of 1.21).

Diffusion kinetic↗

Electrical Analysis of Pulsed Laser Annealed Poly-Si:Ga/SiOx Passivating Contacts

High-efficiency single crystalline silicon (c-Si) solar cells require precise control of dopant diffusion profiles and highly active doping concentrations through thermal annealing. Conventional furnace annealing has been successfully employed for dopant diffusion in aluminum back-surface field (Al-BSF), passivated emitter and rear contact (PERC), and Topcon cells. However, furnace annealing limits some next-generation polycrystalline silicon (poly-Si) on SiOx passivating contacts. This limitation largely affects p-type passivating contacts by not being able to: 1) control the dopant diffusion accurately to prevent B segregation at the SiOx and c-Si interface; and 2) provide highly activated doping concentrations for low contact resistivity. In this contribution, we use a nanosecond excimer laser to examine the passivation quality and electrical performance of both B- and Ga-doped poly-Si/SiOx passivating contacts. The core of this method is to take advantage of the nonequilibrium nature of the anneal through rapid melting and recrystallizing the poly-Si in a short timescale and to achieve doping concentrations above the solid solubility limit. Simulations were performed on polished surfaces to visualize the doping diffusion profiles under different laser conditions, and secondary ion mass spectrometry (SIMS) was used to verify the experimental diffusion profiles post pulsed laser melting (PLM). The results show the dopant diffusion profiles can be tuned precisely through PLM. The electrical analysis using VdP-Hall measurements reveals that the active doping concentration reached 10^21 cm-3 for B and ~2×10^21 cm-3 for Ga in poly-Si, far exceeding their solid solubility limit in Si, with a nearly 100% dopant activation achieved for B, and 20% for Ga. This highly activated dopant profile results in a low contact resistivity of <15 mO·cm^2 for B-doped contacts and <40 mO·cm^2 for Ga-doped contacts. We demonstrate that the passivation quality of the laser annealed contacts shows a linear dependency with the laser energy density and the number of pulses. PLM was also performed on random pyramid textured samples, and we found that preferential melting exists for the tips of pyramids, which can be pursued further for selective pinhole opening. Additionally, we perform thermal stability tests of these hyperdoped passivating contacts and show that under certain laser conditions, the B-doped samples exhibit good thermal stability with near full activation for temperatures up to 800°C. Lastly, we transfer this technique to a larger area by overlapping laser spots and show uniform passivation quality and electrical performance. We are currently implementing this PLM annealed p-type passivating contacts into double-side back-junction poly-Si/SiOx passivating contact devices. Our results show that PLM can not only enable the next generation of high-efficiency c-Si photovoltaics technology but can also benefit the c-Si integrated circuit industry through hyperdoping other semiconductors with atoms that may otherwise have a low solubility.

doping diffusion profile↗

Hydrogen Diffusion Through Stainless Steel

SAND2025-14430O Hydrogen Diffusion Through Stainless Steel is a tool that employs a finite difference method algorithm to solve Fick's law and other mass transport models pertinent to the diffusion of hydrogen through stainless steel. Additionally, the software uses particle swarm optimization to determine physical parameters based on experimental data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Mason, Tyler [Sandia National Lab. (SNL-CA), Liver↗

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗

Sub-100 mA/cm 2 CO 2 -to-CO Reduction Current Densities in Hierarchical Porous Gold Electrocatalysts Made by Direct Ink Writing and Dealloying

While most research efforts on CO 2 -to-CO reduction electrocatalysts focus on boosting their selectivity, the reduction rate, directly proportional to the reduction current density, is another critical parameter to be considered in practical applications. This is because mass transport associated with the diffusion of reactant/product species becomes a major concern at a high reduction rate. Nanostructured Au is a promising CO 2 -to-CO reduction electrocatalyst for its very high selectivity. However, the CO 2 -to-CO reduction current density commonly achieved in conventional nanostructured Au electrocatalysts is relatively low (in the range of 1–10 mA/cm 2 ) for practical applications. In this work, we combine direct ink writing-based additive manufacturing and dealloying to design a robust hierarchical porous Au electrocatalyst to improve the mass transport and achieve high CO 2 -to-CO reduction current densities on the order of 64.9 mA/cm 2 with CO partial current density of 33.8 mA/cm 2 at 0.55 V overpotential using an H-cell configuration. Although the current density achieved in our robust hierarchical porous Au electrocatalyst is one order of magnitude higher than the one achieved in conventional nanostructured electrocatalysts, we found that the selectivity of our system is relatively low, namely 52%, which suggests that mass transport remains a critical issue despite the hierarchical porous architecture. We further show that the bulk dimension of our electrocatalyst is a critical parameter governing the interplay between selectivity and reduction rate. In conclusion, the insights gained in this work shed new light on the design of electrocatalysts toward scale-up CO 2 reduction and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpenetrated Structures for Enhancing Ion Diffusion Kinetics in Electrochemical Energy Storage Devices

The architectural design of electrodes offers new opportunities for next-generation electrochemical energy storage devices (EESDs) by increasing surface area, thickness, and active materials mass loading while maintaining good ion diffusion through optimized electrode tortuosity. However, conventional thick electrodes increase ion diffusion length and cause larger ion concentration gradients, limiting reaction kinetics. We demonstrate a strategy for building interpenetrated structures that shortens ion diffusion length and reduces ion concentration inhomogeneity. This free-standing device structure also avoids short-circuiting without needing a separator. The feature size and number of interpenetrated units can be adjusted during printing to balance surface area and ion diffusion. Starting with a 3D-printed interpenetrated polymer substrate, we metallize it to make it conductive. This substrate has two individually addressable electrodes, allowing selective electrodeposition of energy storage materials. Using a Zn//MnO 2 battery as a model system, the interpenetrated device outperforms conventional separate electrode configurations, improving volumetric energy density by 221% and exhibiting a higher capacity retention rate of 49% compared to 35% at temperatures from 20 to 0 °C. Our study introduces a new EESD architecture applicable to Li-ion, Na-ion batteries, supercapacitors, etc.

25 ENERGY STORAGE↗

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T g ), oxygen to carbon (O:C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core–shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10 2 –10 12 Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity <10 2 Pa s. Phase separation while in a liquid phase state, i.e., liquid–liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM 2.5 ) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Properties of Sodium Borosilicate Glasses as a Function of Sulfur Content

SO3 additions, up to 3.0 wt%, were systematically investigated for effects on the physical properties of sodium borosilicate glass melted in air, with a sulfur-free composition of 50SiO2-10Al2O3-12B2O3-21Na2O-7CaO (wt%). Solubility measurements, using electron microscopy chemical analysis, determined the maximum loading to be ~1.5 wt% SO3. It was found that sulfur (here as sulfate) in this glass increased the glass transition temperature, thermal diffusivity, heat capacity, and thermal conductivity, and decreased the mass density. Structural analysis, performed with Raman spectroscopy, indicated that the borosilicate network polymerized with sulfur additions, presumably due to Na2O being required to charge compensate the ionic SO42- additions, thus becoming unavailable to form non-bridging oxygen in the silicate network. It is postulated that this increased crosslinking of the borosilicate backbone lead to a structure with higher dimensionality and average bond energy. This increased the mean free paths and vibration frequency of the phonons, which resulted in the observed increase in thermal properties.

Thermal properties, Borosilicate glass, Sulfur sol↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗