Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic modulus”

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 433 records · Page 24

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗

Exploring Causal Physical Mechanisms via Non-Gaussian Linear Models and Deep Kernel Learning: Applications for Ferroelectric Domain Structures

Rapid emergence of multimodal imaging in scanning probe, electron, and optical microscopies has brought forth the challenge of understanding the information contained in these complex data sets, targeting the intrinsic correlations between different channels, and further exploring the underpinning causal physical mechanisms. Here, we develop such an analysis framework for Piezoresponse Force Microscopy. We argue that under certain conditions, we can bootstrap experimental observations with the prior knowledge of materials structure to get information on certain nonobserved properties, and demonstrate linear causal analysis for PFM observables. We further demonstrate that the strength of individual causal links between complex descriptors can be ascertained using the deep kernel learning (DKL) model. In this DKL analysis, we use the prior information on domain structure within the image to predict the physical properties. This analysis demonstrates the correlative relationships between morphology, piezoresponse, elastic property, etc., at nanoscale. The prediction of morphology and other physical parameters illustrates a mutual interaction between surface condition and physical properties in ferroelectric materials. Overall, this analysis is universal and can be extended to explore the correlative relationships of other multichannel data sets, and allow for high-fidelity reconstruction of underpinning functionalities and physical mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanomechanical mapping and strain rate sensitivity of microcrystalline cellulose

Nanoindentation provides a convenient and high-throughput means for mapping mechanical properties and for measuring the strain rate sensitivity of a material. Here, nanoindentation was applied to the study of microcrystalline cellulose. Here, constant strain rate nanoindentation revealed a depth dependence of nanohardness and modulus, mostly attributed to material densification. Nanomechanical maps of storage modulus and hardness resolved the shape and size of voids present in larger particles. In smaller, denser particles, however, where storage modulus varied little spatially, there was still some spatial dependence of hardness, which can be explained by cellulose’s structural anisotropy. Additionally, hardness changed with the indentation strain rate in strain rate jump tests. The resulting strain rate sensitivity values were found to be in agreement with those obtained by other techniques in the literature.

36 MATERIALS SCIENCE↗

Electrical and Elastic Properties of Individual Single-Layer Nb 4 C 3 T x MXene Flakes

2D carbides and nitrides (MXenes) are widely recognized for their exceptional promise for numerous applications. However, physical property measurements of their individual monolayers remain very limited despite their importance for revealing the intrinsic physical properties of MXenes. The first mechanical and electrical measurements of individual single-layer flakes of Nb 4 C 3 T x MXene, which are prepared via an improved synthetic method are reported. Characterization of field-effect transistor devices based on individual single-layer Nb 4 C 3 T x flakes shows an electrical conductivity of 1024 ± 165 S cm -1 , which is two orders of magnitude higher than the previously reported values for bulk Nb 4 C 3 T x assemblies, and an electron mobility of 0.41 ± 0.27 cm 2 V -1 s -1 . Atomic force microscopy nanoindentation measurements of monolayer Nb 4 C 3 T x membranes yield an effective Young's modulus of 386 ± 13 GPa, assuming a membrane thickness of 1.26 nm. This is the highest value reported for nanoindentation measurements of solution-processable 2D materials, revealing the potential of Nb 4 C 3 T x as a primary component for various mechanical applications. In conclusion, the agreement between the mechanical properties of 2D Nb 4 C 3 T x MXene and cubic NbC suggests that the extensive experimental data on bulk carbides could be useful for identifying new MXenes with improved functional characteristics.

36 MATERIALS SCIENCE↗

Tensile testing and fracture mechanism analysis of polyvinyl alcohol nanofibrous webs

Abstract A tensile properties testing study was conducted to understand the influence of thickness, cross‐head speed (speed of testing), gauge length (GL; specimen test length), and sample shape on important tensile properties of polyvinyl alcohol (PVA) nanofiber webs. The effects of each testing parameter on load at break, extension at break, Young's modulus, and tensile stress–strain curve of PVA nanofiber webs are analyzed. The Welch two sample t ‐tests show the significant difference among tested data. Using interaction plots, two‐way analysis of variance, and margin mean plots, the interaction effects among testing parameters have been analyzed. Of all the factors, cross‐head speed, the interaction among GL, and sample thickness (GL: Thickness) and the interaction among GL, testing speed and sample thickness (GL: Speed: Thickness) have significant influence on the tensile properties of PVA nanofiber webs. Moreover, the hypothesized model of mechanism of tensile strain–stress curve of PVA nanofiber webs has been proposed. Based on the model, the tensile strain–stress curve can be split into three stages: linear elastic, partial break up, and complete breakage. This study will provide a better understanding of tensile testing parameters' effects and their interaction effects on the tensile properties of nanowebs.

Lou, Lihua↗

Multiscale Characterization of the Caney Shale — An Emerging Play in Oklahoma

From a hydrocarbon perspective, the Caney Shale has historically been evaluated as a sealing unit, which resulted in limited studies characterizing the rock properties of the Caney Shale and its suitability for hydraulic fracturing. The objective of our research is to help bridge the current knowledge gap through the integration of multiscale laboratory techniques and to characterize the macro- and microscale rock properties of the Caney Shale. We employed an integrated approach for the characterization of the Caney using 200 ft (61 m) of Caney core from a target well in southern Oklahoma. Core observation and petrographic analysis of thin sections were combined to characterize the general rock types and associated fabrics and textures. Mineralogical composition, pore system architecture, and rock fabric were analyzed using x-ray diffraction (XRD), scanning electron microscopy/energy dispersive x-ray spectroscopy (SEM/EDS), and focused ion beam (FIB)-SEM. In addition, rebound hardness and indentation testing were carried out to determine rock hardness (brittleness) and elasticity, respectively. With the integrated multiscale characterization, three mixed carbonate-siliciclastic rock types were identified — mudstone, calcareous siltstone, and silty carbonate — likely representing a spectrum of deposition from low to relatively high energy environments in the distal portions of a ramp system. Silty carbonate contains mostly interparticle pores. The calcareous siltstones and silty mudstones contain a combination of organic matter pores and interparticle pores. Each of the rock types shows unique mineralogical compositions based on XRD. The mudstone lithofacies has the highest clay content and the least carbonate content. Calcareous siltstones show moderate carbonate and clay content. Silty carbonate indicates the highest carbonate content with the least clay content. In an order of mudstone, calcareous siltstone, and silty carbonate, rebound hardness and Young’s modulus show an increasing trend. As a result of rock-fluid interactions, there are potential scaling reactions during completion and production that could ultimately affect permeability and production rates. Overall, the proposed multiscale integration approach is critical for the geologic characterization of most rocks. However, in shale reservoirs dominated by microporosity and microstructure where engineered fractures are expected to provide permeability at a reservoir scale, successful integration is essential. An optimized, integrated geological characterization of the Caney Shale that is well aligned with the engineering designs in drilling, completing, and producing wellbores will ultimately lead to optimal production while providing safe and environmentally responsible operations.

Wang, Yulun↗

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗

Assessment and Non-Destructive Evaluation of the Influence of Residual Solvent on a Two-Part Epoxy-Based Adhesive Using Ultrasonics

Polymers are increasingly being used in higher demanding applications due to their ability to tailor the properties of structures while allowing for a weight and cost reduction. Solvents play an important role in the manufacture of polymeric structures since they allow for a reduction in the polymer’s viscosity or assist with the dispersion of fillers into the polymer matrix. However, the incorrect removal of the solvent affects both the physical and chemical properties of polymeric materials. The presence of residual solvent can also negatively affect the curing kinetics and the final quality of polymers. Destructive testing is mainly performed to characterize the properties of these materials. However, this type of testing involves using lab-type equipment that cannot be taken in-field to perform in situ testing and requires a specific sample preparation. Here, a method is presented to non-destructively evaluate the curing process and final viscoelastic properties of polymeric materials using ultrasonics. In this study, changes in longitudinal sound speed were detected during the curing of an aerospace epoxy adhesive as a result of variations in polymer chemistry. To simulate the presence of residual solvent, samples containing different weight percentages of isopropyl alcohol were manufactured and tested using ultrasonics. Thermogravimetric analysis was used to show changes in the decomposition of the adhesive due to the presence of IPA within the polymer structure. Adding 2, 4, and 6 wt.% of IPA decreased the adhesive’s lap shear strength by 40, 58, and 71%, respectively. Ultrasonics were used to show how the solvent influenced the curing process and the final sound speed of the adhesive. Young’s modulus and Poisson’s ratio were determined using both the longitudinal and shear sound speeds of the adhesive. Using ultrasonics has the potential to non-invasively characterize the quality of polymers in both an in-field and manufacturing settings, ensuring their reliability during use in demanding applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗

Predicting Asteroid Material Properties from a DART-like Kinetic Impact

Abstract NASA’s Double Asteroid Redirection Test (DART) mission is the first full-scale test of the kinetic impactor method for asteroid deflection, in which a spacecraft intentionally impacts an asteroid to change its trajectory. DART represents an important first step for planetary defense technology demonstration, providing a realistic assessment of the effectiveness of the kinetic impact approach on a near-Earth asteroid. The momentum imparted to the asteroid is transferred from the impacting spacecraft and enhanced by the momentum of material ejected from the impact site. However, the magnitude of the ejecta contribution is dependent on the material properties of the target. These properties, such as strength and shear modulus, are unknown for the DART target asteroid, Dimorphos, as well as most asteroids since such properties are difficult to characterize remotely. This study examines how hydrocode simulations can be used to estimate material properties from information available post-impact, specifically the asteroid size and shape, the velocity and properties of the impacting spacecraft, and the final velocity change imparted to the asteroid. Across >300 three-dimensional simulations varying seven material parameters describing the asteroid, we found many combinations of properties could reproduce a particular asteroid velocity. Additional observations, such as asteroid mass or crater size, are required to further constrain properties like asteroid strength or outcomes like the momentum enhancement provided by impact ejecta. Our results demonstrate the vital importance of having as much knowledge as possible prior to an impact mission, with key material parameters being the asteroid’s mass, porosity, strength, and elastic properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗

Pressure dependence of intermediate-range order and elastic properties of glassy Baltic amber

Amber is a unique example of a fragile glass that has been extensively aged below its glass transition temperature, thus reaching a state that is not accessible under normal experimental conditions. In this paper, we studied the medium-range order of Baltic amber by x-ray diffraction (XRD) at high pressures. The pressure dependences of the low-angle XRD intensity between 0 and 5 Å -1 were measured from 0 to 7.3 GPa by the energy-dispersive XRD. The first diffraction peak at 1.1 Å -1 and ambient pressure has a doublet structure consisting of the first sharp diffraction peak (FSDP) at 1.05 Å -1 and the second feature at 1.40 Å -1 . The peak position and the width of the FSDP increase as the pressure increases, while the intensity of the FSDP decreases. Below P 0 = 2.4 GPa, the rapid increase of the FSDP peak position was observed, while above P 0 , the gradual increase was observed. Below P 0 , voids and holes in a relatively low-density state are suppressed, whereas above P 0 , the suppression becomes mild. Such a change suggests the crossover from the low- to high-density state at P 0 . There is a close correlation between the pressure dependence of XRD and previously reported sound velocity results. The correlation between the mean-square fluctuation of the shear modulus on the nanometer scale and fragility in amber and other glass formers is also discussed.

36 MATERIALS SCIENCE↗

NanoMechanics: Elasticity and Friction in Nano‐Objects

A large scientific and technological effort is underway to investigate the properties of two‐dimensional (2D) materials to become building blocks in integrated nano‐electronic and photonic circuits, composites, coatings, energy harvesting nano‐ systems, nano‐sensors, and nano‐electro‐mechanical systems (NEMS).While several experiments and calculations have revealed exciting novel phenomena in these nanostructures, many scientific and technological questions remain open. In particular, understanding and controlling the structure and mechanical properties of 2D materials at the interface with a solid surface is of key importance to enable the aforementioned nanotechnologies, as well as to fulfill the potential of strain engineering to tailor their properties. The vision of this DoE research program was to investigate novel mechanical properties and phenomena in 1D and 2D materials with the overarching goal of defining a new basic understanding of mechanical behavior in nano and quantum systems. The group of the PI has developed in the last years several studies on the mechanical properties of Carbon nanotubes and oxide nanobelts, more recently the PI has focused her attention on the properties of two‐dimensional materials, such as graphene and MoS 2 , which are a few‐atomic‐layer thick films and hold a great potential for technological and energy applications. The most studied 2D material is graphene, existing as a single layer of graphite or a few‐layer thick epitaxial graphene film. Graphene possesses a large in‐plane Young’s modulus as well as high intrinsic carrier mobility, and high in‐plane thermal conductivity. Besides graphene, also 2D films of graphene oxide (GO), hexagonal Boron Nitride (h‐BN), and transition metal dichalcogenides such as MoS 2 exhibit unique and excellent properties and hold great promise for nanotechnology applications. This research program was aimed at developing a new basic understanding of the mechanical properties of 1D and 2D materials at the interface with solid substrates, with a focus on the role of defects, materials structure, and substrate interaction. Ultimately, we aimed at developing the basic knowledge and experimental tools for manipulating 1D and 2D materials’ structure, stacking, substrate interaction, defects, and number of layers to realize an entirely new class of ultra‐hard, ultra‐thin, and ultra‐light active materials with ad‐ hoc mechanical properties for a variety of applications.

2D Materials↗

Data for: Ultrasonic characterization of material heterogeneities in stainless steel parts fabricated by powder bed fusion

These are the raw ultrasonic waveform files and nanoindentation measurements for three additively manufactured 316L stainless steel components with different fabrication parameters. Each part's length is divided into four segmented regions where their fabrication energy densities change. Part V+ begins at 33 J/mm^3 and increases in energy density by 3 J/mm^3 with each segment. Part C remains at a constant 33 J/mm^3 throughout the part. Part V- begins at 33 J/mm^3 and decreases in energy density by 3 J/mm^3 with each segment. We have found that the ultrasonic waves in regions with higher energy densities will traverse more quickly, as shown by the shorter time of flight, and vice versa. Similarly, nanoindentation measurements of reduced modulus and hardness follow this trend. We have measured the direction in which the parts were fabricated (build direction) and the segmented regions that are perpendicular to the build direction (transverse directions). The build direction longitudinal wave velocity measurements are considerably lower than their transverse direction counterparts, indicating anisotropy between the two directions. In Part C, we have removed material from one of its surfaces and recorded its ultrasonic and nanoindentation measurements with each material removal iteration. Changes in the material properties are more prominent by nanoindentation, suggesting material heterogeneity between the part's surface and interior. The README.txt file has information on how to navigate the files and process the data.

Anisotropy↗

Spin Transitions and Compressibility of ε-Fe 7 N 3 and γ'-Fe 4 N: Implications for Iron Alloys in Terrestrial Planet Cores

Iron nitrides are possible constituents of the cores of Earth and other terrestrial planets. Pressure-induced magnetic changes in iron nitrides and effects on compressibility remain poorly understood. In this work, we report synchrotron X-ray emission spectroscopy (XES) and X-ray diffraction (XRD) results for ε-Fe 7 N 3 and γ'-Fe 4 N up to 60 GPa at 300 K. The XES spectra reveal completion of high- to low-spin transition in ε-Fe 7 N 3 and γ'-Fe 4 N at 43 and 34 GPa, respectively. The completion of the spin transition induces stiffening in bulk modulus of ε-Fe 7 N 3 by 22% at ~40 GPa, but has no resolvable effect on the compression behavior of γ'-Fe 4 N. Fitting pressure-volume data to the Birch-Murnaghan equation of state yields V 0 = 83.29 ± 0.03 (Å 3 ), K 0 = 232 ± 9 GPa, K 0 ' = 4.1 ± 0.5 for nonmagnetic ε-Fe 7 N 3 above the spin transition completion pressure, and V 0 = 54.82 ± 0.02 (Å 3 ), K 0 = 152 ± 2 GPa, K 0 ' = 4.0 ± 0.1 for γ'-Fe 4 N over the studied pressure range. By reexamining evidence for spin transition and effects on compressibility of other candidate components of terrestrial planet cores, Fe 3 S, Fe 3 P, Fe 7 C 3 , and Fe 3 C based on previous XES and XRD measurements, we located the completion of high- to low-spin transition at ~67, 38, 50, and 30 GPa at 300 K, respectively. The completion of spin transitions of Fe 3 S, Fe 3 P, and Fe 3 C induces elastic stiffening, whereas that of Fe 7 C 3 induces elastic softening. Changes in compressibility at completion of spin transitions in iron-light element alloys may influence the properties of Earth's and planetary cores.

58 GEOSCIENCES↗

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗