Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Low-PAPR interleaved single-carrier FDM scheme for optical wireless communications

As a potential complementary access technology for the 5G wireless systems, optical wireless communication (OWC) gains extensive attention for decades owing to its numerous advantages of broad license-free spectrum, immunity to electromagnetic interference and high-level privacy. However, two rigorous obstacles to achieve high data-rate OWC transmission are the severe nonlinear clipping of transmitter and high-frequency attenuation of OWC system. In this paper, we propose the first interleaved single-carrier frequency-division multiplexing (I-SC-FDM) scheme for OWC systems. Compared with orthogonal frequency-division multiplexing (OFDM), I-SC-FDM has a lower peak-to-average power ratio, which makes it more immune to the nonlinearity clipping in OWC systems. Meanwhile, for the bandwidth-limited OWC systems, I-SC-FDM has better performance on resistance to the serious high-frequency distortion. The simulation results indicate that, under the transmitter nonlinearity and optical-wireless diffuse fading channel, the maximum Q{sup 2} of I-SC-FDM is about 3.44 dB and 3.14 dB higher than that of OFDM when 4-ary quadrature amplitude modulation (4-QAM) and 16-QAM are modulated, respectively. The results show the feasibility and advantages of I-SC-FDM for cost-sensitive OWC systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

BioCARS: Synchrotron facility for probing structural dynamics of biological macromolecules

A major goal in biomedical science is to move beyond static images of proteins and other biological macromolecules to the internal dynamics underlying their function. This level of study is necessary to understand how these molecules work and to engineer new functions and modulators of function. Stemming from a visionary commitment to this problem by Keith Moffat decades ago, a community of structural biologists has now enabled a set of x-ray scattering technologies for observing intramolecular dynamics in biological macromolecules at atomic resolution and over the broad range of timescales over which motions are functionally relevant. Many of these techniques are provided by BioCARS, a cutting-edge synchrotron radiation facility built under Moffat leadership and located at the Advanced Photon Source at Argonne National Laboratory. BioCARS enables experimental studies of molecular dynamics with time resolutions spanning from 100 ps to seconds and provides both time-resolved x-ray crystallography and small- and wide-angle x-ray scattering. Structural changes can be initiated by several methods—UV/Vis pumping with tunable picosecond and nanosecond laser pulses, substrate diffusion, and global perturbations, such as electric field and temperature jumps. Studies of dynamics typically involve subtle perturbations to molecular structures, requiring specialized computational techniques for data processing and interpretation. In this review, we present the challenges in experimental macromolecular dynamics and describe the current state of experimental capabilities at this facility. As Moffat imagined years ago, BioCARS is now positioned to catalyze the scientific community to make fundamental advances in understanding proteins and other complex biological macromolecules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Database of Nonaqueous Proton-Conducting Materials

This work presents the assembly of 48 papers, representing 74 different compounds and blends, into a machine-readable database of nonaqueous proton-conducting materials. SMILES was used to encode the chemical structures of the molecules, and we tabulated the reported proton conductivity, proton diffusion coefficient, and material composition for a total of 3152 data points. The data spans a broad range of temperatures ranging from -70 to 260 °C. To explore this landscape of nonaqueous proton conductors, DFT was used to calculate the proton affinity of 18 unique proton carriers. The results were then compared to the activation energy derived from fitting experimental data to the Arrhenius equation. It was found that while the widely recognized positive correlation between the activation energy and proton affinity may hold among closely related molecules, this correlation does not necessarily apply across a broader range of molecules. This work serves as an example of the potential analyses that can be conducted using literature data combined with emerging research tools in computation and data science to address specific materials design problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2024

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Alluvial fan morphology: A self-similar free boundary problem description

In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.

58 GEOSCIENCES↗
DOE OSTI2018

Insights into the Interconnection of the Electrodes and Electrolyte Species in Lithium–Sulfur Batteries Using Spatially Resolved Operando X-ray Absorption Spectroscopy and X-ray Fluorescence Mapping

The lithium–sulfur (Li–S) battery chemistry has attracted great interest in the last decade because of its outstanding theoretical gravimetric energy density compared to the state-of-the-art lithium-ion battery technology. However, practically achieved energy density is still far below the theoretical value, even in small laboratory-scale batteries. The problems seen in laboratory-scale batteries will inevitably increase during scale-up to large application-format cells, as the electrolyte to active material (AM) ratio will need to be reduced in these cells to achieve high gravimetric energy density on cell-level basis. Our paper shows the unique possibility of X-ray fluorescence (XRF) mapping to visualize the spatial distribution of the AM inside operating Li–S batteries in all cell components [working electrode (WE), separator, and counter electrode (CE)]. Through a combination of operando XRF mapping and X-ray absorption spectroscopy, we show that unless self-discharge is efficiently prevented, the AM can completely dissolve and distribute throughout the cell stack within a time frame of 2 h, causing poor capacity retention. Finally, using a polysulfide diffusion barrier between the WE and the CE, we successfully suppress these processes and thereby establish a tool for examining the sealed cathode electrode compartment, enabling sophisticated studies for future optimization of the WE processes.

25 ENERGY STORAGE↗
DOE OSTI2020

IDENTIFICATION OF ADSORBA TION OF ADSORBATE FT-IR BANDS USING IN-SI -IR BANDS USING IN-SITU TECHNIQUES: Pd SPECIATION AND ADSORPTION CHEMISTRY OF Pd-ZEOLITES FOR PASSIVE NOx ADSORPTION

To meet increasingly stringent automotive emissions standards, further improvements in catalytic converter design are necessary. Current automotive catalyst systems are effective at eliminating emission of nitrogen oxides (NOx) once the catalyst reaches operational temperature (~200 °C). NOx emitted at lower catalyst temperatures now comprises most of the NOx released during a typical test cycle. Referred to as “the cold start problem” this issue has come to the forefront of automotive catalyst development, as mitigating these emissions is necessary to further reduce automotive emissions. Passive NOx adsorbers present an appealing solution to the cold start problem, these being a class of materials that chemisorb exhaust components such as NOx, carbon monoxide (CO) and hydrocarbons at near-ambient temperatures, and then desorb these compounds once the downstream catalyst has reached operational temperature. An effective passive NOx adsorber must have several properties: high NOx adsorption at near-ambient temperatures, near-complete NOx desorption at temperatures within the operational range, high thermal stability, and resistance to automotive exhaust components at high temperatures. The potential environmental impact of such a system is substantial, as NOx emissions currently result in the formation of millions of tons of smog and acid rain each year. Pd-exchanged zeolites have shown promise for deployment as Passive NOx adsorbers, though much remains to be understood about their adsorption chemistry and deactivation. In-situ Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS) provides a convenient probe of adsorbed species, most automotive exhaust components possessing IR-active chemical bonds. By examining the evolution of IR bands under various pretreatments and adsorbates, the overall Pd-speciation and adsorptive zeolite sites of each material can be characterized, and the identities of IR bands can be deduced. In this work, microreactor-MS analysis of the adsorption and desorption behavior of these materials was also examined, these results being coupled with in-situ DRIFTS temperature programmed desorption (TPD) to correlate desorption events with specific adsorbed species. A pair of zeolite frameworks of similar Si/Al ratio but differing pore size were examined, Beta zeolite (BEA) and Chabazite (CHA) representing a medium- and small-pore framework, respectively. The effect of Pd-loading on BEA was examined, as well as the various deactivation pathways and active sites of each material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Gradient nanostructuring via compositional means

Nanocrystalline metals are inherently unstable against thermal and mechanical stimuli, commonly resulting in significant grain growth. Also, while these metals exhibit substantial Hall-Petch strengthening, they tend to suffer from low ductility and fracture toughness. With regard to the grain growth problem, alloying elements have been employed to stabilize the microstructure through kinetic and/or thermodynamic mechanisms. And to address the ductility challenge, spatially-graded grain size distributions have been developed to facilitate heterogeneous deformation modes: high-strength at the surface and plastic deformation in the bulk. In the present work, we combine these two strategies and present a new methodology for the fabrication of gradient nanostructured metals via compositional means. Here, we have demonstrated that annealing a compositionally stepwise Pt-Au film with a homogenous microstructure results in a film with a spatial microstructural gradient, exhibiting grains which can be twice as wide in the bulk compared to the outer surfaces. Additionally, phase-field modeling was employed for the comparison with experimental results and for further investigation of the competing mechanisms of Au diffusion and thermally induced grain growth. This fabrication method offers an alternative approach for developing the next generation of microstructurally stable gradient nanostructured films.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Simulating the smallest ring world of Chariklo

A ring system consisting of two dense narrow rings has been discovered around Centaur Chariklo. The existence of these rings around a small object poses various questions about their origin, stability, and lifetime. In order to understand the nature of Chariklo’s rings, we perform global N-body simulations of the self-gravitating collisional particle rings for the first time. We find that Chariklo should be denser than the ring material in order to avoid the rapid diffusion of the rings. If Chariklo is denser than the ring material, fine spiral structures called self-gravity wakes occur in the inner ring. These wakes accelerate the viscous spreading of the ring significantly and typically occur on timescales of about 100 years for m-sized ring particles, which is considerably shorter than the timescales suggested in previous studies. The existence of these narrow rings implies smaller ring particles or the existence of shepherding satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An Extended Trajectory Mechanics Approach for Calculating the Path of a Pressure Transient: Derivation and Illustration

Following an approach used in quantum dynamics, an exponential representation of the hydraulic head transforms the diffusion equation governing pressure propagation into an equivalent set of ordinary differential equations. Using a reservoir simulator to determine one set of dependent variables leaves a reduced set of equations for the path of a pressure transient. Unlike the current approach for computing the path of a transient, based on a high-frequency asymptotic solution, the trajectories resulting from this new formulation are valid for arbitrary spatial variations in aquifer properties. For a medium containing interfaces and layers with sharp boundaries, the trajectory mechanics approach produces paths that are compatible with travel time fields produced by a numerical simulator, while the asymptotic solution produces paths that bend too strongly into high permeability regions. The breakdown of the conventional asymptotic solution, due to the presence of sharp boundaries, has implications for model parameter sensitivity calculations and the solution of the inverse problem. For example, near an abrupt boundary, trajectories based on the asymptotic approach deviate significantly from regions of high sensitivity observed in numerical computations. In contrast, paths based on the new trajectory mechanics approach coincide with regions of maximum sensitivity to permeability changes.

42 ENGINEERING↗