Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polymorphous nature of cubic halide perovskites

Many common crystal structures can be described by a single (or very few) repeated structural motif (“monomorphous structures”) such as octahedron in cubic halide perovskites. Interestingly, recent accumulated evidence suggests that electronic structure calculations based on such macroscopically averaged monomorphous cubic ($\textit{Pm} – 3m$) halide perovskites obtained from x-ray diffraction, show intriguing deviations from experiment. These include systematically too small band gaps, dielectric constants dominated by the electronic, negative mixing enthalpy of alloys, and significant deviations from the measured pair distribution function. We show here that a minimization of the systems $\textit{T}$ = 0 internal energy via density functional theory reveals a distribution of different low-symmetry local motifs, including tilting, rotations, and B-atom displacements (“polymorphous networks”). This is found only if one allows for larger-than-minimal cell size that does not geometrically exclude low symmetry motifs. As the (super) cell size increases, the energy is lowered relative to the monomorphous cell, and stabilizes after ~32 formula units (≥160 atoms) are included. Being a result of nonthermal energy minimization in the internal energy without entropy, this correlated set of displacements must represent the intrinsic geometry preferred by the underlying chemical bonding (lone pair bonding), and as such has a different origin than the normal, dynamic thermal disorder modeled by molecular dynamics. Indeed, the polymorphous network, not the monomorphous ansatz, is the kernel structure from which high temperature thermal agitation develops. The emerging physical picture is that the polymorphous network has an average structure with high symmetry, yet the local structural motifs have low symmetries. In this work we find that, compared with monomorphous counterparts, the polymorphous networks have significantly lower predicted total energies, larger band gaps, and ionic dominated dielectric constants, and agree much more closely with the observed pair distribution functions. An analogous polymorphous situation is found in the paraelectric phase of a few cubic oxide perovskites where local polarization takes the role of local displacements in halide perovskites, and in the paramagnetic phases of a few 3 d oxides where the local spin configuration takes that role.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Hyperheavy spherical and toroidal nuclei: The role of shell structure

The properties of toroidal hyperheavy even-even nuclei and the role of toroidal shell structure are extensively studied within covariant density functional theory. The general trends in the evolution of toroidal shapes in the Z ≈ 130–180 region of the nuclear chart are established for the first time. These nuclei are stable with respect to breathing deformations. The most compact fat toroidal nuclei are located in the Z ≈ 136, N ≈ 206 region of the nuclear chart, but thin toroidal nuclei become dominant with increasing proton number and toward proton and neutron drip lines. The roles of toroidal shell structure, its regularity, supershell structure, and shell gaps as well as the role of different groups of the pairs of the orbitals in its formation are investigated in detail. The lowest in energy solutions at axial symmetry are characterized either by large shell gaps or low density of the single-particle states in the vicinity of the Fermi level in at least one of the subsystems (proton or neutron). Related quantum shell effects are expected to act against the instabilities in breathing and sausage deformations for these subsystems. The investigation with a large set of covariant energy density functionals reveals that substantial proton Z = 154 and 186 and neutron N = 228, 308, and 406 spherical shell gaps exist in all functionals. The nuclei in the vicinity of the combination of these particle numbers form the islands of stability of spherical hyperheavy nuclei. Furthermore, the study suggests that the N = 210 toroidal shell gap plays a substantial role in the stabilization of fat toroidal nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Johnson, Praise Noah↗
DOE OSTI2021

Charge Radius of Neutron-Deficient Ni 54 and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii

he nuclear root-mean-square charge radius of Ni 54 was determined with collinear laser spectroscopy to be R ( Ni 54 ) = 3.737 ( 3 ) fm . In conjunction with the known radius of the mirror nucleus Fe 54 , the difference of the charge radii was extracted as Δ R ch = 0.049 ( 4 ) fm . Based on the correlation between Δ R ch and the slope of the symmetry energy at nuclear saturation density ( L ), we deduced 21 ≤ L ≤ 88 MeV . The present result is consistent with the L from the binary neutron star merger GW170817, favoring a soft neutron matter EOS, and barely consistent with the PREX-2 result within 1 σ error bands. Our result indicates the neutron-skin thickness of Ca 48 as 0.15–0.21 fm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nanoscale interplay of native point defects near Sr-deficient SrxTiO 3 /SrTiO 3 interfaces

SrTiO 3 has many applications involving interfaces with other materials, but defects that affect the properties of SrTiO 3 films can also play a significant role at its heterointerfaces and even at junctions with nonstoichiometric SrTiO 3 . Depth-resolved cathodoluminescence spectroscopy (DRCLS) combined with systematic cation Sr content reduction in SrxTiO 3 ultrathin films grown on SrTiO 3 showed an interplay of native point defects and electronic structure within the Sr-deficient film and how interplay extends tens of nanometers into the substrate below. Understanding how defects form and affect interface electronic structure during epitaxial growth is central to improving complex oxide devices. Controlling the balance of oxygen vacancy defects with strontium vacancies and other acceptor-like defects can improve control over free carrier densities. Likewise, control over nanoscale defect charge distributions can advance new device features such as two-dimensional hole gases and the performance of existing devices such as ferroelectric tunnel junctions. This study shows how DRCLS directly measures the relative densities and spatial distributions of multiple native defects within and extending away from nanoscale SrxTiO 3 /SrTiO 3 junctions and how their interplay varies with controlled epitaxial growth. In turn, these growth-dependent defect variations can help guide SrTiO 3 epitaxial growth with other complex oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bulk density, grain size, carbon, and nitrogen content of floodplain fine sediment along the East River, Colorado collected from 2014 to 2017

The dataset is comprised of grain size, bulk density, carbon and nitrogen content of soils sampled on the floodplain of the East River 4 to 11 km downstream of Gothic, CO near Crested Butte. This dataset was collected to characterize the sedimentological properties and carbon inventory of the East River floodplain. The data help quantify the flux and storage of sediment and carbon between the East River and the alluvial floodplain bounding the river. The samples were collected between 2014 and 2017. The soil samples come from a total of 156 locations and were sampled in the upper 40 to 60 cm of the floodplain and in abandoned channel segments or oxbows. Not all samples were analyzed for all properties, but 26 samples were analyzed for bulk density, full particle size distribution, carbon, and nitrogen content. An additional 48 samples were analyzed for particle size fraction greater than and less than 2 mm and total organic content determined by loss on ignition. The dataset is comprised of a total of six separate data files based on the sample type and analysis and another file lists the sample identifications and locations of samples with a complete suite of analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, $ΔR$$^{mir}_{ch}$, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability α D . In particular, this quantity correlates poorly with both the neutron skin and L. Here, we demonstrate that $ΔR$$^{mir}_{ch}$ is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on $ΔR$$^{mir}_{ch}$, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Microscopic origin of reflection-asymmetric nuclear shapes

The presence of nuclear ground states with stable reflection-asymmetric shapes is supported by rich experimental evidence. Theoretical surveys of odd-multipolarity deformations predict the existence of pear-shaped isotopes in several fairly localized regions of the nuclear landscape in the vicinity of near-lying single-particle shells with Δℓ = Δj = 3. We analyze the role of isoscalar, isovector, neutron-proton, neutron-neutron, and proton-proton multipole interaction energies in inducing the onset of reflection-asymmetric ground-state deformations. The calculations are performed in the framework of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory using two Skyrme energy-density functionals and density-dependent pairing force. We show that reflection-asymmetric ground-state shapes of atomic nuclei are driven by the odd-multipolarity neutron-proton (or isoscalar) part of the nuclear interaction energy. This result is consistent with the particle-vibration picture, in which the main driver of octupole instability is the isoscalar octupole-octupole interaction giving rise to large E3 polarizability. The necessary condition for the appearance of localized regions of pear-shaped nuclei in the nuclear landscape is the presence of parity doublets involving Δℓ = Δj = 3 proton or neutron single-particle shells. This condition alone is, however, not sufficient to determine whether pear shapes actually appear, and, if so, what the corresponding reflection-asymmetric deformation energies are. The predicted small reflection-asymmetric deformation energies result from dramatic cancellations between even- and odd-multipolarity components of the nuclear binding energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bulk density, grain size, carbon content and isotopes, and nitrogen content of floodplain sediment along the Koyukuk River, Alaska collected in Summer 2018.

The dataset is comprised of grain size, bulk density, carbon and nitrogen content, stable carbon isotopes, and radiocarbon content of soils sampled on the floodplain of the Koyukuk River near Huslia, AK. This dataset was collected to characterize sediment properties and floodplain organic carbon stocks in the discontinuous permafrost floodplain of the Koyukuk River. The data help quantify the flux and storage of sediment and carbon between the Koyukuk River and the surrounding floodplain as the river migrates laterally. The samples were collected in June 27 - July 9, 2018. The soil samples come from a total of 24 locations and were sampled in the upper 490 cm of the river cutbanks, point bars, and floodplain. Samples were analyzed for bulk density, full particle size distribution, nitrogen content and carbon content and isotopes, though not all samples were included in each analysis. An additional 10 samples of floodplain woody debris were measured for radiocarbon content. Data, detailed methods, analyses, and interpretations were previously published at: Douglas, M. M., Lingappa, U. F., Lamb, M. P., Rowland, J. C., West, A. J., Li, G., et al. (2021). Impact of river channel lateral migration on microbial communities across a discontinuous permafrost floodplain. Applied and Environmental Microbiology, 87(20), AEM.01339-21. https://doi.org/10.1128/AEM.01339-21 Douglas, M. M., Li, G. K., Fischer, W. W., Rowland, J. C., Kemeny, P. C., West, A. J., et al. (2022). Organic carbon burial by river meandering partially offsets bank erosion carbon fluxes in a discontinuous permafrost floodplain. Earth Surface Dynamics, 10(3), 421–435. https://doi.org/10.5194/esurf-10-421-2022 The dataset is comprised of a total of three separate data files based on the sample type and analysis, including a file that lists the sample locations, permafrost presence, and geomorphic landform.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Time-resolved in situ visualization of the structural response of zeolites during catalysis

Zeolites form an important class of materials that are widely used in various fields. Their unique properties come primarily from their size, shape, and connectivity of their sub-nm crystalline voids, the Si/Al ratio of the anionic framework, and the number and the nature of charge-balancing cations.The inhomogeneous distribution of the Si/Al ratio within crystals, and hence the cation density, the presence of organic residues, and their inhomogeneous distribution within crystals sensitively affect their performances because these factors influencethe intra-crystalline diffusion rates of the reactants and products. However, the information regarding these inhomogeneities in zeolite crystals cannot be obtained by conventional analytical tools because they only provide information regardingthe average properties of zeolite crystals. Here we employin situX-ray free electron laser-based time-resolved coherent X-ray diffraction imaging to investigate the internal deformations originating from the inhomogeneous Cu ion distributions in Cu-exchanged ZSM-5 zeolite crystals during the deoxygenation of nitrogen oxides with propene. We show that the interactions between the reactants and the active sites lead to an unusual strain distribution,which is confirmed by density functional theory simulations. Moreover, this elucidates the factors responsible for distinct active regions,which are essential to maximizing some applications with enhanced performances. These observations provide insights into the role of structural inhomogeneity in zeolites during a catalytic processand will assist the future design of zeolites for their applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗
DOE OSTI2021

Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes

Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

Valdivia, Pedro↗
DOE OSTI2023

Creep Properties of Shale and Predicted Impact on Proppant Embedment for the Caney Shale, Oklahoma

The Caney shale is an emerging hydrocarbon play located in southwest Oklahoma, USA. Within the Caney shale exist facies which were initially dubbed “reservoir” and “ductile” based on evaluation of well logging data. While past work has shown the distinction of “brittle” and “ductile” is not mechanically justifiable according to formal definitions, here the current work shows some important differences between nominally ductile and reservoir zones. First, the “ductile” zones are more clay rich and have textural differences which can be expected to lead to differences in mechanical properties. One important impact of these differences is observed in triaxial creep experiments showing the “ductile” zones are more prone to creep deformation. Numerical simulations predict the “reservoir” zones will experience very little proppant embedment due to creep deformation of hydraulic fractures around proppant particles. On the other hand, “ductile” zones can be expected to undergo creep-driven proppant embedment leading to loss of fracture aperture ranging up to 100% loss, depending upon the spatial density of the proppant distribution. Hence, this research shows the identification of nominally “ductile” zones from well logs, while a misnomer, can be useful in finding clay-rich, creep-prone zones which will be the most prone to proppant embedment and hence vulnerable to greater production decline over time.

58 GEOSCIENCES↗