Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical distributions”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

A General Reliability Model for Ni-BaTiO3-Based Multilayer Ceramic Capacitors

The evaluation for potential space project applications of multilayer ceramic capacitors (MLCCs) with Ni electrode and BaTiO3 dielectric material requires an in-depth understanding of the MLCCs reliability. A general reliability model for Ni-BaTiO3 MLCCs is developed and discussed in this paper. The model consists of three parts: a statistical distribution; an acceleration function that describes how a capacitors reliability life responds to external stresses; and an empirical function that defines the contribution of the structural and constructional characteristics of a multilayer capacitor device, such as the number of dielectric layers N, dielectric thickness d, average grain size r, and capacitor chip size A. Application examples are also discussed based on the proposed reliability model for Ni-BaTiO3 MLCCs.

statistical modeling↗
NASA NTRS2014

A Method for Obtaining High Time and Spatial Resolution Convective Cloud Top Data for the TTL

A method for obtaining high time and spatial resolution convective cloud top data for the TTL Leonhard Pfister, Eric Jensen, Rei Ueyama, Eliot Atlas, and Maria Navarro Convective systems in the tropics have a maximum in the cloud top altitude distribution of about 13.5 km. However, there is a significant tail to this distribution -- a few percent reach the cold point tropopause (CPT) at 16.5 km, and there has been clear evidence of convective mass deposited as high as 19 km in the tropics. The region between 13.5 km and the cold point tropopause is transitional, between the free tropical troposphere where convective mixing dominates, and the stratosphere where slow upward ascent dominates. In this region (the Tropical Tropopause Layer), convective injection, slow ascent, and mixing from midlatitudes all have similar time scales. So, even though only a few percent of convective systems reach the CPT, convection is important. Space Based Lidar and cloud radar measurements have yielded information about long term average statistical distributions of cloud altitude as a function of location. However, we also need time-dependent cloud top altitude and cloud top potential temperature information, primarily to understand the water vapor and TTL cloud distributions. This is because the effect of convection depends on the local temperature, and on the subsequent temperature history. Time dependent cloud top information is also needed to understand short-lived tracers because cross-isentropic flow is time and space dependent. This paper presents a method of obtaining time and space dependent convective cloud top theta (and altitude) information using 3-hourly geostationary brightness temperature data, coupled with global 3 -hourly rainfall estimates and temperature analyses. We explore different mixing algorithms to obtain the most reasonable agreement with near-simultaneous observations by cloudsat and calipso. Observations of short-lived tracers from ATTREX, coupled with short-term trajectories are used to test the method's accuracy. An important caveat is the ambiguity of evaluating convective cloud top altitudes under from combined cloudsat and calipso measurements.

cold point tropopause (CPT)↗
DOE OSTI2026

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

The 1990 annual statistics and highlights report

The National Space Science Data Center (NSSDC) has archived over 6 terabytes of space and Earth science data accumulated over nearly 25 years. It now expects these holdings to nearly double every two years. The science user community needs rapid access to this archival data and information about data. The NSSDC has been set on course to provide just that. Five years ago the NSSDC came on line, becoming easily reachable for thousands of scientists around the world through electronic networks it managed and other international electronic networks to which it connected. Since that time, the data center has developed and implemented over 15 interactive systems, operational nearly 24 hours per day, and is reachable through DECnet, TCP/IP, X25, and BITnet communication protocols. The NSSDC is a clearinghouse for the science user to find data needed through the Master Directory system whether it is at the NSSDC or deposited in over 50 other archives and data management facilities around the world. Over 13,000 users accessed the NSSDC electronic systems, during the past year. Thousands of requests for data have been satisfied, resulting in the NSSDC's sending out a volume of data last year that nearly exceeded a quarter of its holdings. This document reports on some of the highlights and distribution statistics for most of the basic NSSDC operational services for fiscal year 1990. It is intended to be the first of a series of annual reports on how well NSSDC is doing in supporting the space and Earth science user communities.

Green, James L.↗
DOE OSTI2020

EMP-Resilient Electric Grid Transformer Analysis

The response of a high - voltage (HV) transformer to fast rise time voltages, such as that from a electromagnetic pulse (EMP) can result in interruption of power distribution and possibly system failure. To help identify these potential occurrences, it is necessary to develop a transformer model that not only captures the input/output response of the transformer but also the internal behavior. The model constructed should cover the frequency band of interest while capturing the internal physical and electrical characteristics. This broad-band, high-fidelity model would enable the prediction of unwanted effects through simulation. A proposed modeling scheme for a HV transformer is described in Part 1 of this report. Part 2 of this report details assessments of internal voltage and electrical field holdoff testing of transformer insulation dielectric breakdown, including comparison of low frequency (DC/60 Hz) holdoff to rise times characteristic of lightning (1 s) and EMP E1 transients (10 - 30 ns). This initial project is a path toward establishing electrical grid transformer failure criteria for EMP voltage transients. We developed modeling methods and measured breakdown electrical field statistical distributions for direct current, 60 Hz, lightning and EMP characteristic voltage rise times. Methods of nanosecond-scale capacitive discharge unit high voltage source development, suggestions for derating of 60 Hz insulation maximum electrical fields for EMP nanosecond pulse voltage withstand rating, and potential methods for increasing transformer resilience to such fast rise time pulses are described.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Link statistics of dislocation network during strain hardening

Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.

Crystal plasticity↗
DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗
DOE OSTI2024

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z = (z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y(0)=(y sub 1(0),...,y sub D(0)) knowledge of the statistical distribution of the random errors in y(0). The data space Y containing y(0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x (e.g., energy or dissipation rate), Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. CSI is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface. Neither the heat flow nor the energy bound is strong enough to permit estimation of B(r) at single points on the CMB, but the heat flow bound permits estimation of uniform averages of B(r) over discs on the CMB, and both bounds permit weighted disc-averages with continous weighting kernels. Both bounds also permit estimation of low-degree Gauss coefficients at the CMB. The heat flow bound resolves them up to degree 8 if the crustal field at satellite altitudes must be treated as a systematic error, but can resolve to degree 11 under the most favorable statistical treatment of the crust. These two limits produce circles of confusion on the CMB with diameters of 25 deg and 19 deg respectively.

Backus, George E.↗
NASA NTRS2011

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2022

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Synergy of tensile strength and high cycle fatigue properties in a novel additively manufactured Al-Ni-Ti-Zr alloy with a heterogeneous microstructure

Alloy design strategies in additive manufacturing (AM) to achieve grain refinement and terminal eutectic solidification have been introduced to engineer Al alloys having microstructural hierarchy and heterogeneity. Such alloy design strategies enable crack-free builds with an expanded AM processing window and pushed the strength limit in Al alloys. However, fatigue performance of Al alloys made by AM is restricted by the presence of process induced defects and its stochasticity. In this work, tensile and high cycle fatigue (HCF) behavior of a novel Al-Ni-Ti-Zr alloy with a heterogeneous microstructure is studied in the as-built condition, supplemented by detailed microstructural and mechanical characterization. Excellent strength-ductility synergy of 342 MPa and 16% failure strain achieved in the alloy was associated with the microstructural attributes that pertain to the novel alloy. Additionally, the alloy showed excellent HCF performance with a fatigue endurance limit to ultimate tensile strength ratio of 0.29 in flexural fatigue mode. The study revealed the existence of multiple crack retardation mechanisms and favorable crack propagation pathways through the fine-grained regions which enabled good fatigue performance to the alloy. Further, a probabilistic model has been used to estimate the fatigue life of the alloy as a function of the stochastic microstructure by utilizing the statistical distribution of pores, solid-state inclusions, and grains in the AM Al alloy. Finally, the model parametric trends are consistent with the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mitigation of spatial nonstationarity with vision transformers

Spatial nonstationarity, the location variance of features’ statistical distributions, is ubiquitous in many natural settings. For example, in geological reservoirs rock matrix porosity varies vertically due to geomechanical compaction trends, in mineral deposits grades vary due to sedimentation and concentration processes, in hydrology rainfall varies due to the atmosphere and topography interactions, and in metallurgy crystalline structures vary due to differential cooling. Conventional geostatistical modeling workflows rely on the assumption of stationarity to be able to model spatial features for geostatistical inference. Nevertheless, this is often not a realistic assumption when dealing with nonstationary spatial data and this has motivated a variety of nonstationary spatial modeling workflows such as trend and residual decomposition, cosimulation with secondary features, and spatial segmentation and independent modeling over stationary subdomains. The advent of deep learning technologies has enabled new workflows for modeling spatial relationships. However, there is a paucity of demonstrated best practice and general guidance on mitigation of spatial nonstationarity with deep learning in the geospatial context. We demonstrate the impact of two common types of geostatistical spatial nonstationarity on deep learning model prediction performance and propose the mitigation of such impacts using self-attention (vision transformer) models. We demonstrate the utility of vision transformers for the mitigation of nonstationarity with relative errors as low as 10%, exceeding the performance of alternative deep learning methods such as convolutional neural networks. We establish best practice by demonstrating the ability of self-attention networks for modeling large-scale spatial relationships in the presence of commonly observed geospatial nonstationarity.

58 GEOSCIENCES↗
DOE OSTI2020

Inverse design of two-dimensional graphene/h-BN hybrids by a regressional and conditional GAN

Design of materials with desired properties is currently laborious and heavily relies on intuition of researchers through a trial-and-error process. To tackle this challenge, in this work we propose a novel regressional and conditional generative adversarial network (RCGAN) for inverse design of representative two-dimensional materials, the graphene and boron-nitride (BN) hybrids. RCGAN incorporates a supervised regressor network, thus overcoming the common technical barrier in the traditional unsupervised GANs, which cannot generate data when fed with continuous and quantitative labels. RCGAN can autonomously generate graphene/BN hybrids given any target bandgap values. These structures are distinguished from the ones used for training and exhibit high diversity for a given bandgap. Moreover, they exhibit high fidelity, yielding bandgaps within ~10% MAE F of the desired bandgaps as validated by density functional theory (DFT) calculations. Analysis by the principle component analysis (PCA) and modified locally linear embedding (MLLE) reveals that the generator has successfully generated structures following the statistical distribution of the real structures. It implies the possibility of the RCGAN in recognizing physical rules hidden in the high-dimensional data. The novel strategy for designing regressional GAN architecture together with the successful application to inverse design of materials would inspire further exploration in research fields beyond materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Onion-Like nanoparticles of the metal–organic framework UiO-66 synthesized by sequential spike crystal growth

Reproducible size and morphology of metal–organic framework (MOF) particles are essential for tuning these materials for applications like heterogeneous catalysis. Particle size measurements from electron micrographs provide information about the impact of changing synthetic parameters on MOF growth, but reported data are usually limited to averages and statistical distributions. Further elucidating synthetic control parameters could facilitate studies of facet-dependent and non-classical growth mechanisms. Here, in this work, phase-pure nanoparticles of the MOF UiO-66 were synthesized using alternating spikes of zirconium and hafnium-based nodes, forming an onion-like structure with alternating layers of Zr-UiO-66 and Hf-UiO-66. In conventional, bright-field transmission electron micrographs, Hf-UiO-66 and Zr-UiO-66 layers appeared darker and lighter, respectively. Alternating layers of Hf-UiO-66 and Zr-UiO-66 were also apparent in elemental maps obtained using scanning transmission electron microscopy with energy-dispersive X-ray spectroscopy. Average particle sizes increased monotonically with each spike, and no Zr-only or Hf-only particles were observed after the first spike, indicating that new material was consistently incorporated onto existing seeds. Surface area of the onion-like particles and their uptake of Ni 2+ from solution were both between the values measured for Zr-UiO-66 and Hf-UiO-66. Results demonstrate the use of sequential spikes of nodes in monitoring crystal growth in MOFs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bulk and position-specific isotope geochemistry of natural gases from the Late Cretaceous Eagle Ford Shale, south Texas

Here, a dataset of bulk and position-specific isotope compositions of shale gases from the Late Cretaceous Eagle Ford Shale, south Texas is reported. The chemical and bulk isotopic compositions of the seven samples in this study and those available from the literature show that across the play, produced natural gas from the Eagle Ford Shale preserves a wide range of gas compositions (wetness, ~5 to >90%) and δ 13 C values of C 1 –C 3 . The depth profiles of the isotopic compositions suggest that the western region of the Eagle Ford Shale deposit had experienced as much as 700–800 m of uplift-erosion or gas migration. The timing of and lost-gas fractions by gas expulsion events very likely affected their gas compositions and bulk isotope compositions. Some deep (>3000–3500 m), matured (%Ro > 1.5) gases with heavy δ 13 C(C 1 –C 3 ) values indicate significant loss (>50%) of the early-stage gases. For the position-specific isotope deviations of propane, the ΔC 2-1 values of the five samples show small decrease with well depths, while the ΔH 2-1 values have a general increasing trend. Although some samples fall very close to the equilibrium model trajectories with reasonable calculated temperatures (138–148 °C), the position-specific isotope compositions of propane from the Eagle Ford Shale gases are likely results of thermal cracking of various organic molecules within the source rocks, which have different activation energy for cracking and non-statistical distributions of C/H isotope within them.

58 GEOSCIENCES↗
DOE OSTI2025

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗