Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective 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 1,297 records · Page 72

Measurement of gas-concentration-driven permeation for the examination of permeability, solubility, and diffusivity in varying materials

This manuscript describes the development and operation of an apparatus for the measurement of steady-state and transient gas permeation through different types of solid materials with varying geometries. It is capable of operation from 293 K to 673 K and could theoretically be used with any non-corrosive gas or a mix of gases, although only hydrogen isotopes are used in the current study. A quadrupole mass spectrometer is used to measure permeation fluxes as low as 1011 molecules/s. This unique test setup allows for the simultaneous measurement of diffusivity, solubility, and permeability. Furthermore, varying the pressure in the fore-sample volume allows for tests of Sievert’s law and can give information on surface effects.

Young, Katherine T. (ORCID:0000000285845576)↗

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Forecasting fluid-injection induced seismicity to choose the best injection strategy for safety and efficiency

Induced seismicity poses a challenge to the development of Enhanced Geothermal Systems (EGS). Improving monitoring and forecasting techniques is essential to mitigate induced seismicity and thereby fostering a positive perception of EGS projects among local authorities and population. Induced seismicity is the result of complex and coupled thermo-hydro-mechanical-chemical mechanisms. Injection flux and pressure are crucial controlling parameters for both hydraulic stimulation and circulation protocols. We develop a methodology combining a hydro-mechanical model with a seismicity rate model to estimate the magnitude and frequency of mainshocks and aftershocks induced by fluid injection. We apply the methodology to the case of the Basel EGS (2006, Switzerland) to compare the effects of progressive, cyclic and constant injections on the mechanical response of discrete faults. Results from the coupled hydro-mechanical models show that the pore pressure diffusion and consequent enhancement of fault permeability are limited to the vicinity of the injection well during cyclic injection. Additionally, constant injection induces seismicity from the start of the injection but enhances the permeability of most of the faults within a shorter duration, inducing less post-injection seismicity. The methodology can be adapted to any numerical model and allows new projects to be developed by anticipating the safest injection protocol. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.

58 GEOSCIENCES↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗

Defect Kinetics and Control for Module Reliability

Potential induced degradation is currently one of the most important module degradation mechanisms. It has been suggested that stacking faults decorated with sodium from the module glass are responsible for this effect and authors have also shown the reversibility of this effect upon reverse biasing of the module. The importance of sodium in the failure mechanism is clear, however, little is known regarding the factors that control its diffusion into the wafer, making it nearly impossible to predict the performance of a given module and engineer it to be better. Sodium migration from module glass into silicon cells and the resulting module degradation is a clear example of how defect kinetics can determine overall module performance and long-term reliability. To the detriment of the industry and its bankability, no quantitative models yet exist to predict defect-assisted module degradation, limiting the progress in improving reliability. In particular, the understanding of defect behavior under high electric fields, under stresses imparted by encapsulation or temperature, and under real operating conditions over long periods of time is a crucial gap in the current state-of-the-art. In this work we developed a Defect-Device-Degradation model to predict defect behavior and its impact on device performance over the module operational lifetime using experimentally-determined defect parameterizations. The validated model will provide a platform for manufacturing process optimization across input materials and architectures to avoid deleterious defects upstream and enable enhanced module robustness.

36 MATERIALS SCIENCE↗

Optimizing Sub-Sampled STEM Imaging for Beam Sensitive Materials and Dynamic Processes

The goal of this research is to quantify how sub-sampling and inpainting (broadly classified as methods in Artificial Intelligence) can be used to obtain the highest resolution images under any experimental conditions, and to use this to reproducibly control dynamic in-situ processes in gases and liquids on the nanoscale. As atomic diffusion under non-extreme conditions is typically in the sub-ms time domain, the ability to use this sub-sampling approach to reduce the effect of the electron beam and to speed-up the imaging process represents the optimum combination of spatial and temporal resolution achievable in a conventional microscope. This LDRD funding demonstrated the feasibility of imaging dynamics in complex energy systems using this methodology by establishing the optimum experimental sampling conditions (optimizing the number of pixels, speed, contrast and resolution) for each experiment. For in-situ liquid cell experiments this analysis shows that simply by distributing the same dose and dose rate maximally in space and time, the kinetics for radiolysis products can be modified significantly. The final stage of this work, i.e. the experimental demonstration of its use for nucleation and growth, was severely curtailed by issues with COVID19. However, all the experiments are designed and currently underway, with the goal to complete the final demonstration and publication as soon as safety concerns for laboratory access are met.

36 MATERIALS SCIENCE↗

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗

Ballistic-diffusive phonon transport in cellulose nanocrystals by ReaxFF molecular dynamics simulations

Cellulose nanocrystals (CNCs), a ubiquitous nano-sized natural biopolymer, have drawn considerable attention in recent years due to their excellent mechanical and chemical properties. However, thermal properties of CNCs, which are critical for their potential applications as structural or functional materials, are less investigated. In this work, their thermal properties were systematically studied by molecular dynamics (MD) simulations in terms of the effects of polymorphs, size, temperature, and strain. It was found that the thermal conductivities depend on the polymorphs of CNCs, but the thermal conductance of unit chain of each polymorph is similar. Moreover, the strong dependence of the thermal conductivities of CNCs on their length and cross-sectional area indicates the existence of remarkable ballistic-diffusive phonon transport. Single chain CNCs can achieve a thermal conductivity of ~6 W/(m∙K), implying their potential as high thermal conductive polymer. The mechanical strain can enhance the alignment of cellulose skeleton, which reduces phonon scattering and thus increases thermal conductivity. Possible underlying mechanism of the varying thermal conductivity was further discussed and attributed to the shift of phonon frequencies and varied degree of orientation. This study reveals the sophisticated ballistic-diffusive phonon transport in CNCs, and pave path to future design of CNC-based materials with desirable thermal properties.

42 ENGINEERING↗

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗

Experimental and modeling study of the effect of confinement on the thermal decomposition of organic materials

Simultaneous Thermal Analysis is widely used to study thermal decomposition of a variety of organic materials, and the choice of experimental conditions is very important to obtain reliable mass loss and heat flow curves. Here in this work we investigate the effect of confinement on the thermal decomposition of organic materials. As an example, we explore high-density polyethylene (HDPE) and one of its pyrolysis products, namely eicosane (C 20 H 42 ), through experiments and modeling in open crucibles and crucibles closed with a pierced lid. A new model is developed for evaporation and sublimation from open and pinhole thermal-analysis pans. It considers the kinetic resistance of evaporation at the liquid-vapor interface, inhibition of evaporation by buildup of vapor inside the pan, thermal expansion of vapor inside the pan, and diffusion of vapor from the pan. The model is validated using simultaneous heat-flow and mass-loss measurements of n-eicosane evaporation for various pinhole sizes, down to 50 µm. The importance of product evaporation inhibition for measuring polymer decomposition was demonstrated using HDPE. This study sheds light on the effect of confinement on the mass loss rate of organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Look At Dissolved Organic Carbon In Streams and Its Effect On DGT Hg Collection

Savannah River Site (SRS) and Oak Ridge National Laboratory (ORNL) were both major material production sites for the U.S.'s nuclear program during the 20. century. Both sites discharged mercury (Hg) to streams. Oak Ridge streams have accumulated significant levels of Hg from early 'historical' discharges. Diffusive Gradient in Thin Films (DGT) samplers may be a labor and material saving method to monitor the Hg contamination of streams at both sites. Humic and Fulvic acids are released during the degradation of plant and animal residue. They are part of the Dissolved Organic Matter that colors the black water streams of the southern U.S.. Humic acid and Fulvic acids have a tendency to sequester heavy metals, like Hg, and may interfere with the reactions needed to draw the Hg through the thin films of the DGT samplers. Project Objective: Use a UV-Vis spectrophotometer to estimate the levels of Dissolved Organic Carbon (DOC) in the streams selected to test the effect of DOC on DGT samplers. Conclusion: UV-Vis is a convenient way to estimate the amount of DOC in streams and may be used to provide more context in DGT Hg sampling. efforts. Future Work: Analyze data from SRS and ORNL DGT samples; Compare the results of the DGT samplers with more traditional speciating methods from ORNL; Use a TOC Spectrophotometer to gain a more accurate picture of stream DOC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nitrogen effects on radiation response in 12Cr ferritic/martensitic alloys

Ferritic/martensitic steels are one of the best candidates for structural materials for high dose applications in next generation nuclear reactors. The composition of structural materials must be optimized for reliable service during irradiation. This study reports the effect of interstitial nitrogen on radiation response in 12Cr ferritic/martensitic HT9 steels having a controlled nitrogen concentration. Results show that a high amount of ‘free’ nitrogen in the matrix stabilizes the interstitial clusters which leads to (i) larger loop sizes (ii) lower loop density and (iii) slightly reduced radiation induced hardening. It also affects diffusion mechanism of Ni and formation of Ni/Si-rich precipitates.

36 MATERIALS SCIENCE↗

Collective Summary of sCO2 Materials Development (Supercritical Transformational Electric Power Generation (STEP) Level 2 Milestone Report) (Parts I - II)

Polymers such as PTFE (polytetrafluorethylene or Teflon), EPDM (ethylene propylene diene monomer) rubber, FKM fluoroelastomer (Viton), Nylon 11, Nitrile butadiene (NBR) rubber, hydrogenated nitrile rubber (HNBR) and perfluoroelastomers (FF_202) are commonly employed in super critical CO 2 (sCO2) energy conversion systems. O-rings and gaskets made from these polymers face stringent performance conditions such as elevated temperatures, high pressures, pollutants, and corrosive humid environments. In FY 2019, we conducted experiments at high temperatures (100°C and 120°C) under isobaric conditions (20 MPa). Findings showed that elevated temperatures accelerated degradation of polymers in sCO2, and that certain polymer microstructures are more susceptible to degradation over others. In FY 2020, the focus was to understand the effect of sCO2 on polymers at low (10 MPa) and high pressures (40 MPa) under isothermal conditions (100°C). It was clear that the same selectivity was observed in these experiments wherein certain polymeric functionalities showed more propensity to failure over others. Fast diffusion, supported by higher pressures and long exposure times (1000 hours) at the test temperature, caused increased damage in sCO2 environments to even the most robust polymers. We also looked at polymers under compression in sCO2 at 100°C and 20 MPa pressure to imitate actual sealing performance required of these materials in sCO2 systems. Compression worsened the physical damage that resulted from chemical attack of the polymers under these test conditions. In FY 2021, the effect of cycling temperature (from 50°C to 150°C to 50°C) for polymers under a steady sCO2 pressure of 20 MPa was studied. The aim was to understand the influence of cycling temperatures of sCO2 for typical polymers under isobaric (20 MPa) conditions. Thermoplastic polymers (Nylon, and PTFE) and elastomers (EPDM, Viton, Buna N, Neoprene, FF202, and HNBR) were subjected to 20 MPa sCO2 pressure for 50 cycles and 100 cycles in separate experiments. Samples were extracted for ex-situ characterization at 50 cycles and upon the completion of 100 cycles. Each cycle constituted of 175 minutes of cycling from 50°C to 150°C. The polymer samples were examined for physical and chemical changes by Dynamic Mechanical and Thermal Analysis (DMTA), Fourier Transform Infrared (FTIR) spectroscopy, and compression set. Density and mass changes immediately after removal from test were measured for degree of swell comparisons. Optical microscopy techniques and micro computer tomography (micro CT) images were collected on select specimens. Evaluations conducted showed that exposures to super-critical CO2 environments resulted in combinations of physical and/or chemical changes. For each polymer, the dominance of cycling temperatures under sCO2 pressures, were evaluated. Attempts were made to qualitatively link the permanent sCO2 effects to polymer micro- structure, free volume, backbone substitutions, presence of polar groups, and degree of crystallinity differences. This study has established that soft polymeric materials are conducive to failure in sCO2 through mechanisms of failure that are dependent on polymer microstructure and chemistry. Polar pendant groups, large atom substitutions on the backbone are some of the factors that are influential structural factors.

36 MATERIALS SCIENCE↗

Galactic Bar Resonances with Diffusion: An Analytic Model with Implications for Bar–Dark Matter Halo Dynamical Friction

Abstract The secular evolution of disk galaxies is largely driven by resonances between the orbits of “particles” (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari, and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies—effects that have, with few exceptions, been previously explored only with complex N -body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles’ slow actions. We solve this equation for various values of the dimensionless diffusion strength Δ, and then apply our theory to the calculation of bar–halo dynamical friction. For Δ = 0, we recover the classic result of Tremaine and Weinberg that friction ultimately vanishes, owing to the phase mixing of resonant orbits. However, for Δ > 0, we find that diffusion suppresses phase mixing, leading to a finite torque. Our results suggest that stochasticity—be it physical or numerical—tends to increase bar–halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Structure and Filler on Mechanical Properties of Direct Ink Write Pads

Direct ink write (DIW) is an extrusion-based additive manufacturing (AM) technique that enables superior repeatability, property tuning, and material customization compared to conventional manufacturing techniques. In this research, DIW is used to print silicone foams that must accommodate the thermal expansion of other parts. This expansion requires the material to compress and subsequently recover. Thus, a key performance metric is reducing “compression set” or deformation retained after the material is compressed. For this study, a silicone resin was compounded with a getter material, 1,4 bis(phenylethynyl)benzene (DEB) mixed with carbon-supported palladium (Pd) and printed in a variety of lattice structures. The effects of the filler and print structure on compression set were evaluated. In these experiments, 0.5-inch diameter coupons were compressed to 75% their original height for 22 hours at 70 °C, and the change in height calculated at intervals between 30 minutes and 24 hours post-decompression. The unfilled pad with skin layers had the best compression set at 0.7% ± 0.1%, followed by the unfilled pad with no skins at 1.5% ± 0.1%. As the number of filled layers increased, compression set performance degraded. The pad with filled skin layers and all filled layers had a compression set of 7.0% ± 0.3%. Additionally, the reaction extent of DEB was found to increase compression set, possibly due to the diffusion of DEB to the surface at elevated temperatures. Imaging reveals extensive crystallization at the surface. These investigations show that a non-reinforcing filler will increase compression set but that print structure can mitigate this effect.

36 MATERIALS SCIENCE↗

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗

Evaluating diffusion and the thermodynamic factor for binary ionic mixtures

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗