Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressure diffusion”

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 271 records · Page 15

In situ investigation of high-pressure hydrogen-induced swelling in elastomers and its correlation with material properties

The resistance of elastomeric materials to high-pressure hydrogen-induced damage is essential for ensuring the reliability of hydrogen infrastructure. Here, in this study, we systematically investigated the swelling behavior and hydrogen transport properties of four elastomer types – EPDM, NBR, FKM, and HNBR – using a custom in-situ view cell system capable of real-time monitoring during decompression from pressures up to 96.5 MPa. Each elastomer was formulated with and without fillers and plasticizers to assess the effects of formulation on swelling response. Thermal desorption analysis (TDA) was employed to determine equilibrium hydrogen content and diffusion coefficients, providing insight into gas uptake and mobility within each material. Correlation analyses using Pearson and Spearman coefficients revealed that the diffusion coefficient showed a stronger relationship with swelling behavior than hydrogen content, highlighting the dominant role of hydrogen mobility. Filled elastomers, particularly those with carbon black, consistently showed reduced swelling due to enhanced stiffness and reduced diffusivity. These results deepen our understanding of diffuso-mechanical interactions in elastomers and support the rational design of sealing materials for high-pressure hydrogen systems.

EPDM↗

Evaluation of Common Thermoplastic Polymers in High-Pressure Cycling Hydrogen Under Ambient and Cold Environments as Applicable to the Hydrogen Infrastructure

Thermoplastic polymers are required to perform unfailingly under stringent conditions of changing pressures (35 MPa to 70 MPa) and temperatures (-40°C to +85°C) in storage and fueling operations. Although diffusivity in glassy polymers appear to be an order of magnitude lower than elastomers, polymer microstructural attributes such as degree of crystallinity and presence of polar and non-polar groups on the main chain with and without branching can play a substantial role in influencing their behaviors in hydrogen environments. In the work described here, the effect of high-pressure hydrogen cycling on PEEK, PTFE, PA11, HDPE, POM (and different commercial grades of these polymers) under ambient and cold (-40°C) conditions is addressed. Ex-situ characterization for polymer changes included density measurements for impact of hydrogen retention, hardness changes using nanoindentation, storage modulus and glass transition changes using dynamic mechanical and thermal analysis (DMTA), degree of crystallization with Differential Scanning Calorimetry (DSC) and tensile testing following ASTM D412. In other evaluations, solid-state nuclear magnetic resonance (NMR), attenuated total reflectance Fourier transform infra-red spectroscopy (ATRFTIR), X-ray diffraction (XRD), and X-ray CT results are presented. The overall impact of these evaluations is to establish a technical basis for the behaviors of common thermoplastics in cycling hydrogen environments under ambient and cold temperatures while establishing the relationship between polymer structure-based properties and hydrogen transport effects.

Menon, Nalini C.↗

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Resistive Switching in SrFeO2.5/Nb:SrTiO3 Heterostructures with Growth-Controlled Film Orientation

Resistive switching, a behavior found in many oxide materials, has the potential to enable emerging computer hardware technologies and architectures. We present resistive switching devices fabricated from epitaxial brownmillerite SrFeO2.5 films with two distinct film orientations, wherein facile oxygen ion diffusion planes are aligned parallel (in-plane) and perpendicular (out-of-plane) with the electrodes. SrFeO2.5 films were grown on (001) oriented Nb:SrTiO3 to enable high-quality interfaces and future integration with Si CMOS technologies. Post-growth vacuum annealing and growth pressure were used to control film orientations, as confirmed by transmission electron microscopy and x-ray diffraction measurements. Films grown with diffusion planes oriented in-plane had oxygen-rich, perovskite-like nanodomains spread throughout the film, and fabricated devices exhibited worse switching consistency and more stochasticity. In contrast, films grown with diffusion planes oriented out-of-plane had a more uniform oxygen-rich perovskite interfacial layer above the bottom electrode, and devices built from this film orientation showed significant statistical improvements in switching voltages and cycling consistency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

Comparison of passivation properties of plasma-assisted ALD and APCVD deposited Al 2 O 3 with SiN x capping

In this paper, we report on the direct comparison of Al 2 O 3 films deposited by plasma-assisted atomic layer deposition (ALD) and atmospheric pressure chemical vapor deposition (APCVD) techniques on the passivation properties of both bare and boron diffused Si surfaces. It is found that ALD and APCVD deposited Al 2 O 3 layers with SiN x capping have very similar chemical composition, hydrogen concentration at the Al 2 O 3 /Si interface, and interface defect density (D it ), but the ALD film has slightly higher negative fixed charge (Q f ) after a high temperature (~750 °C) contact firing cycle typically used for solar cell fabrication. Both films showed excellent surface recombination velocities of <5 cm/s on planar and textured un-diffused Si wafers, but the ALD films gave 1 to 8 fA/cm 2 lower emitter saturation current densities (J 0e ) on implanted boron emitters with sheet resistance in the range of 90–180 Ω/. This is attributed to the observed slightly higher negative charge-induced field-effect passivation in the ALD films. However, SiN x /Al 2 O 3 /p + /n/n + passivated emitter rear totally diffused (PERT) cell structures fabricated with 110 and 150 Ω/ emitter in this study showed no appreciable difference in cell efficiency. Device analysis and simulations showed that even for >23% high-performance solar cells design, with very low bulk and rear contact recombination, the observed small difference in J 0e has negligible effect on cell efficiency. Furthermore, APCVD Al 2 O 3 passivation of B emitter can produce comparable efficiencies to their counterpart ALD Al 2 O 3 passivated n-type Si solar cells.

14 SOLAR ENERGY↗

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

Temperature effect on turbulent burning velocity of lean premixed hydrogen/air flames

Hydrogen has drawn great attention in recent years as a carbon-free fuel. The turbulent burning velocity (S T ) is an important parameter for the design and modeling of hydrogen-fueled engines given the high propagation speed of hydrogen flames. It has been well documented that S T of hydrogen flames can be dramatically increased by thermo-diffusive effects which are sensitive to thermodynamic conditions. Previous studies have mainly focused on the pressure effect on S T of lean hydrogen flames, while the temperature effect has been largely ignored. Here, in the present study, the turbulent burning velocity for a lean hydrogen/air mixture over a wide range of temperatures (300–641 K) and pressures (1–15 atm) is investigated through direct numerical simulations of statistically planar turbulent premixed flames. Results show that the variation of normalized turbulent burning velocity (S T /S L , where S L is the laminar flame speed) with temperature and pressure is mainly controlled by the variation of the stretching factor I 0 . While S T /S L is only marginally dependent on temperature at the atmospheric pressure, it exhibits a decreasing trend with temperature at an elevated pressure (10 atm). This is associated with different temperature dependencies of flame surface area enlargement at the two different pressures, despite the monotonically decreasing trends of I 0 with temperature at both pressures. In addition, under engine-relevant conditions where the temperature and pressure increase simultaneously, the promotion effect of pressure is found to be largely canceled out by the suppression effect of temperature, leading to only a slight increase in I 0 and S T /S L . The observed trends are further explained through detailed flame dynamic analysis. Furthermore, I 0 at different temperatures and pressures is found to correlate very well with the enhancement of fuel consumption rate in the critically strained laminar flames. The present study elucidates the strong impact of temperature on S T of lean premixed hydrogen/air flames at elevated pressures and provides new insights into the modeling of S T , especially for engine-relevant conditions.

Wang, Yiqing [Argonne National Laboratory (ANL), A↗

Spatial profile of argon (1s 5 ) metastables in an electron beam generated plasma

Electron beams with an applied magnetic field generate a secondary cold plasma with a selective chemical composition, featuring low-energy ions and metastable species in the discharge periphery, ideal for low-damage plasma treatment of material substrates. In this work, we studied the plasma generated by an e-beam using a 4 kV voltage in a pure argon gas environment under a magnetic field of 150 G and in the pressure range of 25–90 mTorr. We measured the absolute spatial density profile of argon (1 s 5 ) metastables in an electron beam generated plasma by laser-induced fluorescence and found it to be of the order of 10 16 m −3 . The electron temperature and the electron density measured by a Langmuir probe were of the order of 10 16 m −3 and less than an eV respectively. Electron-impact quenching was identified as a significant loss mechanism for the Ar(1s 5 ) state, leading to the saturation of the metastable density at higher pressures. Outside the primary ionization region, the spatial distribution of argon metastables followed a linear diffusion profile, indicating negligible additional production in those regions.

EEDF↗

Kernel Buffer Volume Fraction Margin of the AGR Designed Fuel Particle

Modeling results used to assess the fuel performance of the TRISO-coated fuel particles as a function of kernel/buffer volume fraction include SiC tangential stress, formation of the buffer/IPyC gap, particle temperature profile, internal particle pressure, fission gas released from the kernel, probability of fuel particle failure, and fission product diffusion. These results were evaluated at two burnup levels and irradiation temperatures to bound expected steady-state irradiation conditions. In general, increasing the kernel/buffer volume fraction increases the SiC stress and subsequently the failure probability of a fuel particle when compared to the AGR designed particle. There was little impact on the fission product diffusion through the particle as the kernel/buffer volume fraction increased.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Stability Transition From Stable to Unstable Frictional Slip With Finite Pore Pressure

Pore fluids are ubiquitous throughout the lithosphere and are commonly invoked as the cause of induced seismicity and slow earthquakes. We perform lab experiments to address these questions for drained fault conditions and low pore pressure. We shear simulated faults at effective normal stress (σ' n ) of 20 MPa and pore pressures Pp from 1 to 4 MPa. We document the full range of lab earthquake behaviors from slow slip to elasto-dynamic rupture and show that slow slip can be explained by the slip rate dependence of the critical rheologic stiffness without dilatancy hardening or other fluid effects. Our fault permeabilities ranges from 10 –18 to 10 –17 m 2 with an initial porosity of 0.1 and estimated fluid diffusion time ≈1 s. Slow slip and quasi-dynamic fault motion may arise from high P p at higher pressures but dilatancy strengthening is not a general requirement.

58 GEOSCIENCES↗

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Field-scale fault reactivation experiments by fluid injection highlight aseismic leakage in caprock analogs: Implications for CO 2 sequestration

Observations on tens-of-meter scale experiments of fault activation by fluid injection conducted in shales allow exploring how aseismic and seismic events may jeopardize the integrity of a sealing caprock overlying a CO 2 sequestration reservoir. We contrast the behavior of shale faults with another set of experiments conducted in carbonates. Significant fluid leakage occurs along the initially low-permeability shale faults when rupture is activated. Most of the leakage pathway closes when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed leakage behavior, which is also caused by fault opening favored by the softness of the shale, and by the structure of the fault zone that prevents fluids from diffusing into the adjacent damage zone. Experiments show a large amount of aseismic deformation. Small-magnitude seismicity (M w < -2.5) is observed outside the pressurized leakage patch. Stress transferred from this aseismic deformation patch can build up to stress-criticality and favor seismicity. Thus, in terms of fault activation in caprocks, aseismic fault slip leading to increased permeability and a loss of seal integrity is of great concern.

58 GEOSCIENCES↗

Understanding Oxide–Metal Interactions During Hot Isostatic Pressing to Diffusion Bond Aluminum Alloy 6061 Plates

The interaction between Mg, Si, and Al 2 O 3 during hot isostatic pressing diffusion bonding of aluminum alloy 6061 (AA6061) plates was investigated through thermodynamic calculations and experimental microstructural characterization. Thermodynamic calculations as functions of temperature, pressure, and composition revealed that the interaction among Mg, Si, and Al 2 O 3 yields Mg 2 Si and either MgO + Al or MgAl 2 O 4 + Al, facilitating the reduction of Al 2 O 3 and allowing Al/Al metallic bonds to form. Total pressure variation had a negligible influence on the oxygen partial pressure, and consequently, the reaction product formation. Oxygen partial pressure variation as a function of temperature and initial amount of Al 2 O 3 determined the formation of either MgO or MgAl 2 O 4 . Experimental Hot Isostatic Pressure (HIP) bonding at 723 K and 833 K under a constant pressure of 1017 atm documented the cooling rate-dependent formation of β–Mg 2 Si precipitates. High-resolution transmission electron microscopy imaging and selected area electron diffraction patterns verified the formation of β–Mg 2 Si and MgO at the interface but did not detect MgAl 2 O 4 . In conclusion, findings from this study clarify the role of thermochemical interactions in oxide disruption and bonding mechanisms during HIP diffusion bonding of AA6061 and provide guidance for optimizing joining processes for monolithic nuclear fuel assemblies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Plasma-droplet interaction study to assess transport limitations and the role of • OH, O • ,H • ,O 2 (a 1 Δ g ),O 3 , He(2 3 S) and Ar(1s 5 ) in formate decomposition

Plasmas interacting with liquid microdroplets are gaining momentum due to their ability to significantly enhance the reactivity transfer from the gas phase plasma to the liquid. This is, for example, critically important for efficiently decomposing organic pollutants in water. In this contribution, the role of •OH as well as non-•OH-driven chemistry initiated by the activation of small water microdroplets in a controlled environment by diffuse RF glow discharge in He with different gas admixtures (Ar, O 2 and humidified He) at atmospheric pressure is quantified. The effect of short-lived radicals such as O• and H• atoms, singlet delta oxygen (O 2 (a 1 Δ g )), O 3 and metastable atoms of He and Ar, besides •OH radicals, on the decomposition of formate dissolved in droplets was analyzed using detailed plasma diagnostics, droplet characterization and ex situ chemical analysis of the treated droplets. The formate decomposition increased with increasing droplet residence time in the plasma, with ~70% decomposition occurring within ~15 ms of the plasma treatment time. The formate oxidation in the droplets is shown to be limited by the gas phase •OH flux at lower H 2 O concentrations with a significant enhancement in the formate decomposition at the lowest water concentration, attributed to e – /ion-induced reactions. However, the oxidation is diffusion limited in the liquid phase at higher gaseous •OH concentrations. Here, the formate decomposition in He/O 2 plasma was similar, although with an order of magnitude higher O• radical density than the •OH density in the corresponding He/H 2 O plasma. Using a one-dimensional reaction–diffusion model, we showed that O 2 (a 1 Δ g ) and O 3 did not play a significant role and the decomposition was due to O•, and possibly •OH generated in the vapor containing droplet-plasma boundary layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗