Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface roughness”

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

Thermomechanical Properties of Neutron Irradiated Al 3 Hf-Al Thermal Neutron Absorber Materials

A thermal neutron absorber material composed of Al 3 Hf particles in an aluminum matrix is under development for the Advanced Test Reactor. This metal matrix composite was fabricated via hot pressing of high-purity aluminum and micrometer-size Al 3 Hf powders at volume fractions of 20.0, 28.4, and 36.5%. Room temperature tensile and hardness testing of unirradiated specimens revealed a linear relationship between volume fraction and strength, while the tensile data showed a strong decrease in elongation between the 20 and 36.5% volume fraction materials. Tensile tests conducted at 200 °C on unirradiated material revealed similar trends. Evaluations were then conducted on specimens irradiated at 66 to 75 °C to four dose levels ranging from approximately 1 to 4 dpa. Tensile properties exhibited the typical increase in strength and decrease in ductility with dose that are common for metallic materials irradiated at ≤0.4T m . Hardness also increased with neutron dose. The difference in strength between the three different volume fraction materials was roughly constant as the dose increased. Nanoindentation measurements of Al 3 Hf particles in the 28.4 vol% material showed the expected trend of increased hardness with irradiation dose. Transmission electron microscopy revealed oxygen at the interface between the Al 3 Hf particles and aluminum matrix in the irradiated material. Scanning electron microscopy of the exterior surface of tensile tested specimens revealed that deformation of the material occurs via plastic deformation of the Al matrix, cracking of the Al 3 Hf particles, and to a lesser extent, tearing of the matrix away from the particles. The fracture surface of an irradiated 28.4 vol% specimen showed failure by brittle fracture in the particles and ductile tearing of the aluminum matrix with no loss of cohesion between the particles and matrix. The coefficient of thermal expansion decreased upon irradiation, with a maximum change of –6.3% for the annealed irradiated 36.5 vol% specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bulk Motions in the Black Hole Jet Sheath as a Candidate for the Comptonizing Corona

Abstract Using two-dimensional general relativistic resistive magnetohydrodynamic simulations, we investigate the properties of the sheath separating the black hole jet from the surrounding medium. We find that the electromagnetic power flowing through the jet sheath is comparable to the overall accretion power of the black hole. The sheath is an important site of energy dissipation as revealed by the copious appearance of reconnection layers and plasmoid chains. About 20% of the sheath power is dissipated between 2 and 10 gravitational radii. The plasma in the dissipative sheath moves along a nearly paraboloidal surface with transrelativistic bulk motions dominated by the radial component, whose dimensionless 4-velocity is ∼1.2 ± 0.5. In the frame moving with the mean (radially dependent) velocity, the distribution of stochastic bulk motions resembles a Maxwellian with an “effective bulk temperature” of ∼100 keV. Scaling the global simulation to Cygnus X-1 parameters gives a rough estimate of the Thomson optical depth across the jet sheath, ∼0.01–0.1, and it may increase in future magnetohydrodynamic simulations with self-consistent radiative losses. These properties suggest that the dissipative jet sheath may be a viable “coronal” region, capable of upscattering seed soft photons into a hard, nonthermal tail, as seen during the hard states of X-ray binaries and active galactic nuclei.

Sridhar, Navin (ORCID:0000000255199550)↗
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
DOE OSTI2024

Local Electric Field Effects on Water Dissociation in Bipolar Membranes Studied Using Core–Shell Catalysts

The local electric field strength is thought to affect the rate of water dissociation (WD) in bipolar membranes (BPMs) at the catalyst–nanoparticle surfaces. Here, we study core–shell nanoparticles, where the core is metallic, semiconducting, or insulating, to understand this effect. The nanoparticle cores were coated with a WD catalyst layer (TiO2 or HfO2) via atomic layer deposition (ALD), and the morphology was imaged with transmission electron microscopy. Irrespective of the core material, these core–shell catalysts displayed comparable WD overpotentials at optimal mass loading, despite the hypothesized differences in the electric field strength across the catalyst particle suggested by continuum electrostatic simulations. Substantial atomic interdiffusion between the core and shell was ruled out by X-ray absorption spectroscopy, X-ray photoelectron spectroscopy, and diffuse reflectance optical measurements. However, the optimal mass loading of catalyst was roughly 1 order of magnitude higher for the conductive and high dielectric core materials than for the low dielectric insulating cores. These findings are consistent with the hypothesis that electric field screening within the core material focuses the electric field drop between particles such that larger film thicknesses can be tolerated. Collectively, these data support the idea that it is the local electric field at the molecular level that controls proton-transfer rates and that the metal core/dielectric-shell constructs introduced here modulate that field. Further materials and synthetic design may enable optimization of the electric field strength across the proton-transfer trajectory at the material surface.

Sarma, Prasad V↗
DOE OSTI2021

The Scale and Activity of Symmetric Instability Estimated from a Global Submesoscale-Permitting Ocean Model

Abstract Symmetric instability (SI) extracts kinetic energy from fronts in the surface mixed layer (SML), potentially affecting the SML structure and dynamics. Here, a global submesoscale-permitting ocean model named MITgcm LLC4320 simulation is used to examine the Stone linear prediction of the maximum SI scale to estimate grid spacings needed to begin resolving SI. Furthermore, potential effects of SI on the usable wind work are estimated roughly: this estimate of SI “activity” is useful for assessing if these modes should be resolved or parameterized. The maximum SI scale varies by latitude with median values from 568 to 23 m. Strong seasonality is observed in the SI scale and activity. The median scale in winter is 188 m globally, 2.5 times of that of summer (75 m). SI is more active in winter: 15% of the time compared with 6% in summer. The strongest SI activity is found in the western Pacific, western Atlantic, and Southern Oceans. The required grid spacings for a global model to begin resolving SI eddies in the SML are 24 m (50% of regions resolved) and 7.9 m (90%) in winter, decreasing to 9.4 m (50%) and 3.6 m (90%) in summer. It is also estimated that SI may reduce usable wind work by an upper bound of 0.83 mW m −2 globally, or 5% of the global magnitude. The sensitivity of these estimates to empirical thresholds is provided in the text.

Dong, Jihai↗
DOE OSTI2024

Effects of Copolymer Structure on Enzyme-Catalyzed Polyester Recycling

Polyesters are an omnipresent material used for a variety of applications (e.g., bottles, packaging, textile, windshield), roughly comprising 8% of plastics produced worldwide. Enzymatic recycling is an emerging solution to deal with the increasingly diverse polyesters that are not suitable for mechanical recycling. However, enzyme activity and efficiency are still the limiting factors impeding enzymatic recycling for different plastic waste forms. The effects of thermal and structural properties (e.g., glass transition temperature, crystallinity, specific surface area), which are determined by chemical composition and preprocessing, directly influence enzyme recycling efficiency. This work investigates two extrusion methods (single screw and twin-screw extrusion) to pretreat a range of copolyesters (RPET, PETG, Ecozen, Tritan and PBT) and modify their properties (i.e., glass transition temperature (T g ), crystallinity (%), and molecular weight (M n )). A PET-specific enzyme, leaf-branch compost cutinase (LCC ICCG ), produced from a fed-batch fermentation of Escherichia coli BL21(DE3), was used for the enzymatic depolymerization of different polyesters. Several copolyesters showed improved depolymerization after pretreatment, as measured by rate and amount of monomers produced. Furthermore, those that did not depolymerize were found to have exceptionally high glass transition temperature or percent crystallinity, highlighting the importance of these physical parameters on conversion efficiency.

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

Evaluation of cloud height, optical thickness, and phase retrievals from the CHROMA algorithm applied to Sentinel-3 OLCI data

We previously developed the Cloud Height Retrieval from O 2 Molecular Absorption (CHROMA) algorithm for the Ocean Color Instrument (OCI) on the new NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission. Here, we apply CHROMA to observations from the Ocean Land Colour Instrument (OLCI) to guide expectations for PACE, as it will take some time to obtain large-scale validation data for OCI. We use cloud top height (CTH), phase, and (for liquid clouds) cloud optical thickness (COT) data from the ground-based Atmospheric Radiation Measurement (ARM) network to evaluate the OLCI retrievals. We found that OLCI and Moderate Resolution Imaging Spectroradiometer (MODIS) CTH compare similarly well to the ARM reference. OLCI has a tendency to underestimate CTH as CTH increases, and algorithm assumptions about cloud geometric thickness may contribute to this. ARM COT from multifilter shadowband radiometers (MFRSR) and Sun photometers are well-correlated with one another, albeit with a roughly 30 % offset on average; OLCI and MODIS COT agree more closely with the MFRSR data. OLCI retrieval uncertainty estimates show skill at telling low-uncertainty cases from high-uncertainty ones, although CTH uncertainties are underestimated. Additionally, we compare the OLCI data to satellite retrievals based on thermal infrared measurements from MODIS and Sea and Land Surface Temperature Radiometer (SLSTR) data. Differences are broadly consistent with physical expectations based on the A-band vs. thermal techniques, although one key challenge in such aggregated comparisons is different cloud masking sensitivities and algorithm failure rates meaning additional sampling differences are introduced. We conclude by discussing the transition to and possible enhancements for PACE OCI.

Sayer, Andrew M. [Univ. of Maryland Baltimore Coun↗
DOE OSTI2021

Evaluation of the Passive Cooling Potential of Thermal Mass Inherent in Medium to Large Commercial Buildings

Typical commercial construction inherently contains large surface areas of relatively thin mass located in the floor slabs. This research asks how effective this mass could be for passive cooling, relative to other mass depths, if it were to be exposed and used as thermal mass. Although much of the rich literature on thermal mass has been conducted with far greater mass depths, a more limited amount of research suggests that because the surface heat transfer rate is limited, the focus should be on the surface area of exposed mass, not its depth. This presents an opportunity for more buildings to behave more thermally massively, if the mass inherent in typical floor slabs contains utility, particularly over diurnal cycles. In all climates analyzed, it was found that there is a pronounced shoulder where energy savings from passive cooling of increasing mass depths was steep until roughly 7.5–10 cm, and beyond this point the achieved energy savings diminished rapidly. When considering the embodied energy of concrete, the incremental benefit of added mass beyond a typical topping slab of 10 cm does not justify the incremental embodied energy cost from a total energy standpoint alone.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A novel fiber-fretting test for tribological characterization of the fiber/matrix interface

Ceramic matrix composites exhibit excellent high temperature properties and are candidate materials for structural applications from jet engines to nuclear fuel cladding. The strength and deformation behavior of these composites are intrinsically tied to the properties of the fiber/matrix interface. This research introduces a novel technique to evaluate the fundamental friction and wear characteristics of these interfaces. A case study on SiC f /PyC/SiC m composites applies in situ fiber fretting to investigate the kinetic friction coefficient and cyclic wear characteristics of the interphase. The experimental methodology is presented in explicit detail to encourage other researchers to explore the capabilities of this technique. Additionally, testing was carried out as a function of fiber roughness, PyC thickness (≈10, 500, 1200 nm), cycle count (10, 100, 1000), and frequency (1, 5, and 10 Hz). The friction coefficient and wear behavior changed with cycle length due to a transition from adhesive to abrasive mechanisms. SEM fractography provided insights to tribological evolution from adhesive sliding to three-body abrasive wear. The structural disorder at the tribo-surface relating to graphitic content and dangling bond density was investigated using Raman spectroscopy. Properties are contextualized with comparison to uniaxial tensile behavior of the parent composite.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Open Circuit Potential of a Au Catalyst during Selective Oxidation of Glycerol

It has been of great interest to understand relationships between the potential of a metal catalyst and its thermochemical catalytic activity, especially in an aqueous phase environment. In the literature, there are correlations of open circuit potential with reaction rates or surface concentrations of reaction intermediates. In this study, we measured the open circuit potential (OCP) of a Au gauze during selective oxidation of glycerol to glyceric acid and hydrogen peroxide in a basic solution, as a function of hydroxide and alcoholate concentrations. It is found that applying a potential to the Au catalyst has no influence on the reaction rate. Although a rough correlation appears to exist between OCP and reaction rate, the data are better fit to an equation which assumes that the potential of the metal (i.e. the OCP) is in equilibrium with the electrochemical potential of the solution, which is defined by the thermodynamic activities of the oxidizing and reducing species. The equation: OCP = constant + Σ ox [(RT/FZ ox )ln[Ox]] - Σ red [(RT/FZ red )ln[Red]] implicitly assumes that the Au metal functions as a probing electrode. It is further found that this equation also applies to the Au-catalyzed H 2 O 2 oxidation/decomposition reaction in a basic medium, and possibly to formic acid dehydrogenation. Here, we postulate that the apparent correlation between OCP and reaction rate is due to the fact that the reaction rate is proportional to the concentrations of the reaction products and/or reactants that define the electrochemical potential of the solution.

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

Impact of n > 1 neoclassical tearing modes on fast ions, plasma rotation, and the onset of disruptive tearing modes in DIII-D ITER baseline scenario

We present integrated TRANSP-kick analyses quantitatively characterizing how n > 1 magnetic islands redistribute fast ions (FIs) and modify momentum balance in low-torque DIII-D H-mode scenarios operated with the ITER normalized parameter set and shape (ITER baseline scenario, or IBS). In this plasma scenario, disruptive neoclassical tearing modes (NTMs) are most commonly seeded by nonlinear 3-wave coupling when the differential rotation between the q = 1 and q = 2 surfaces ($δf_{1, 2}$) approaches zero. Analysis of a DIII-D IBS database of shots unstable to the 2/1 NTM demonstrates that the flattening rotation profile is correlated with n > 1 NTM amplitude growth, with a strong drop in $δf_{1, 2}$ occurring when the radial magnetic perturbation amplitude at the rational surface ($δB$ r s ) is estimated to be about 20 G. The interpretive TRANSP-kick simulations show that the experimentally observed reduction in core rotation is due to the fast-ion kicks. Nonlinearity in momentum losses versus NTM amplitude is found to be correlated with overlap of FI island chains, occurring when $δB$ r s is about 35 G, in rough agreement with the experimental results. This agreement suggests the primary mechanism for the conductive momentum loss in the core is due to chaotic FI orbits and FI losses developing from overlapping resonant FI island chains.

ITER baseline scenario↗
DOE OSTI2023

Emergence of liquid following laser melting of gold thin films

X-ray structural science is undergoing a revolution driven by the emergence of X-ray Free-electron Laser (XFEL) facilities. The structures of crystalline solids can now be studied on the picosecond time scale relevant to phonons, atomic vibrations which travel at acoustic velocities. In the work presented here, X-ray diffuse scattering is employed to characterize the time dependence of the liquid phase emerging from femtosecond laser-induced melting of polycrystalline gold thin films using an XFEL. In a previous analysis of Bragg peak profiles, we showed the supersonic disappearance of the solid phase and presented a model of pumped hot electrons carrying energy from the gold surface to scatter at internal grain boundaries. This generates melt fronts propagating relatively slowly into the crystal grains. By conversion of diffuse scattering to a partial X-ray pair distribution function, we demonstrate that it has the characteristic shape obtained by Fourier transformation of the measured F ( Q ). The diffuse signal fraction increases with a characteristic rise-time of 13 ps, roughly independent of the incident pump fluence and consequent final liquid fraction. This suggests the role of further melt-front nucleation processes beyond grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theoretical and Experimental Advances in High-Pressure Behaviors of Nanoparticles

Using compressive mechanical forces, such as pressure, to induce crystallographic phase transitions and mesostructural changes while modulating material properties in nanoparticles (NPs) is a unique way to discover new phase behaviors, create novel nanostructures, and study emerging properties that are difficult to achieve under conventional conditions. In recent decades, NPs of a plethora of chemical compositions, sizes, shapes, surface ligands, and self-assembled mesostructures have been studied under pressure by in-situ scattering and/or spectroscopy techniques. As a result, the fundamental knowledge of pressure–structure–property relationships has been significantly improved, leading to a better understanding of the design guidelines for nanomaterial synthesis. In the present review, we discuss experimental progress in NP high-pressure research conducted primarily over roughly the past four years on semiconductor NPs, metal and metal oxide NPs, and perovskite NPs. We focus on the pressure-induced behaviors of NPs at both the atomic- and mesoscales, inorganic NP property changes upon compression, and the structural and property transitions of perovskite NPs under pressure. We further discuss in depth progress on molecular modeling, including simulations of ligand behavior, phase-change chalcogenides, layered transition metal dichalcogenides, boron nitride, and inorganic and hybrid organic–inorganic perovskites NPs. These models now provide both mechanistic explanations of experimental observations and predictive guidelines for future experimental design. Here, we conclude with a summary and our insights on future directions for exploration of nanomaterial phase transition, coupling, growth, and nanoelectronic and photonic properties.

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

Reversible Interfacial Hydride Transfer as a Complementary Tool To Measure Molecular Hydricity

Hydride transfer is an essential elementary reaction across the chemical value chain, but there are limited methods available for quantifying thermodynamic hydricity (ΔG H- ), particularly among main group reagents. Herein, we exploit facile H 2 activation and reversible hydride transfer from a metal surface to a molecular reagent, the net hydrogen reduction reaction (HRR), to develop a potentiometric method for quantifying ΔG H- of main group reagents recalcitrant to conventional methods. HRR potentiometry is first validated with a benzimidazole-based hydride donor and then applied to uncover the impact of the reaction environment on hydricity. For a benzimidazole-based hydride donor, HRR equilibrium potentials are roughly invariant across solvents, indicating that the solvent dependence of its hydricity largely reflects the differential solvation of H - across media. For formate, HRR potentials and corresponding hydricities depend strongly on water content. Additionally, for borohydrides, HRR potentiometry reveals that effective hydricity values are strongly influenced by Lewis acid-base adduct formation with the hydride acceptor but are minimally influenced by the countercation. Together with these studies, the advantages, limitations, and practical considerations of the HRR potentiometry method are discussed, highlighting the power of this methodology as a complementary tool to measure molecular hydricity.

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

Downwind control of oceanic air by land: the land wake and its sensitivity to CO 2

Oceans are well-known to be directly altered by global climate forcings such as greenhouse gas changes, but how oceans are indirectly influenced by land and its response to such forcings remains less explored. Here, we assess the present-day and projected future state of a little-explored feature of the climate system—a 'land wake' in relative humidity downwind of the east coast of North America, consisting of low-humidity continental air extending roughly 1000 km over the Atlantic ocean. The wake exists throughout the year, but is supported by high continental temperatures in summer and low continental moisture in winter. The wake is well represented in an ensemble of global climate models (GCMs), qualitatively matching reanalysis data. Under increasing atmospheric CO 2 , the land wake intensifies in GCM simulations through two pathways: the radiative effects of CO 2 on surface temperatures, and the biogeochemical effect of CO 2 on terrestrial vegetation. Vegetation responses to increased CO 2 alter the summer wake from Florida to Newfoundland, and both the radiative and biogeochemical effects of CO 2 drive reductions in coastal cloud cover. These changes illustrate the potential of rapidly changing terrestrial climate to influence coastal regions and the ocean environment downwind of continents through both light conditions and the energy balance of the surface ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Extracting information from X-ray diffraction patterns containing Laue oscillations

The presence of Laue oscillations in a film grown on a solid surface is broadly taken as indicating a high quality, crystallographically aligned film of the targeted compound. In this paper we briefly review the origins of both Laue oscillations and Kiessig fringes and show how they can be used together to determine if extra thickness exists above or below the coherently diffracting domains. The differences between experimental and “ideal” films are discussed and the effect of structural features (roughness, different thickness coherently diffracting domains and thickness in addition to the coherently diffracting domains) are illustrated with experimental and simulated data for metal and mixed-metal chalcogenide films of titanium, bismuth, vanadium/iron, and bismuth/molybdenum. Examples are given showing how quantitative information can be extracted from experimental diffraction patterns.

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

Mid- to Far-Infrared Anisotropic Dielectric Function of HfS 2 and HfSe 2

The far-infrared (far-IR) remains a relatively underexplored region of the electromagnetic spectrum extending roughly from 20 to 100 µm in free-space wavelength. Research within this range has been restricted due to a lack of optical materials that can be optimized to reduce losses and increase sensitivity, as well as by the long free-space wavelengths associated with this spectral region. Here the exceptionally broad Reststrahlen bands of two Hf-based transition metal dichalcogenides (TMDs) that can support surface phonon polaritons (SPhPs) within the mid-infrared (mid-IR) into the terahertz (THz) are reported. In this vein, the IR transmission and reflectance spectra of hafnium disulfide (HfS 2 ) and hafnium diselenide (HfSe 2 ) flakes are measured and their corresponding dielectric functions are extracted. These exceptionally broad Reststrahlen bands (HfS 2 : 165 cm -1 ; HfSe 2 : 95 cm -1 ) dramatically exceed that of the more commonly explored molybdenum- (Mo) and tungsten- (W) based TMDs (≈5–10 cm -1 ), which results from the over sevenfold increase in the Born effective charge of the Hf-containing compounds. Further, this work therefore identifies a class of materials for nanophotonic and sensing applications in the mid- to far-IR, such as deeply sub-diffractional hyperbolic and polaritonic optical antennas, as is predicted via electromagnetic simulations using the extracted dielectric function.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗