SEARCH · Engineering Papers

Results for “Discontinuity factors”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Interface crack in a nonhomogeneous elastic medium

The linear elasticity problem for an interface crack between two bonded half planes is reconsidered. It is assumed that one of the half planes is homogeneous and the second is nonhomogeneous in such a way that the elastic properties are continuous throughout the plane and have discontinuous derivatives along the interface. The problem is formulated in terms of a system of integral equations and the asymptotic behavior of the stress state near the crack tip is determined. The results lead to the conclusion that the singular behavior of stresses in the nonhomogeneous medium is identical to that in a homogeneous material provided the spacial distribution of material properties is continuous near and at the crack tip. The problem is solved for various values of the nonhomogeneity parameter and for four different sets of crack surface tractions, and the corresponding stress intensity factors are tabulated.

Delale, F.↗

DOE OSTI2021

Optimized Carbon Fiber Intermediate Development to Enable High-Volume Manufacturing of Lightweight Automotive Composites

The ongoing pursuit of improved fuel economy and reduced greenhouse gas emissions has resulted in sustained interest for lightweight materials technologies. In this context, carbon fiber composites have captured the imagination of automotive engineers due to the potential to achieve substantial mass reduction when compared to traditional steel construction. That stated, the use of carbon fiber composites in automotive has been limited for the most part to premium supercars and other derivative platforms. In these cases, manufacturing costs are less of an obstacle to implementation, and the performance benefits of carbon fiber have enabled production of structures offering more than 50% weight savings. In practice, translating these low-volume demonstrations onto high-volume vehicle platforms has remained challenging. This can be attributed to several factors, with the absence of suitable high throughput production methodologies being a key impediment. To date, structural, crash critical components have relied upon manufacturing techniques born out of the aerospace industry. This has created a disconnect between automotive production systems that are accustomed to manufacturing multiple parts per minute and the aerospace technologies that have cycle times in the order of hours. Consequently, the focus of this project is the development of manufacturing process technology for carbon fiber composites that can support a mainstream vehicle program at an assumed throughput of 100,000 vehicles per year. In practice, this translates to a part-to-part cycle time of less than 3 minutes. Project participation included contributions from a broad range of academic, industrial and national lab partners. The primary scope of work, being the development of new carbon fiber epoxy compounds that are stable at room temperature and suited to high throughput automated processing. For project management, the work streams were divided into six key areas, with the lead organization in parentheses. • Carbon fiber/epoxy materials formulation development and scale up (Dow) • Simulation of discontinuous near isotropic meso-structure intermediates (Purdue) • Simulation of mechanical performance of compression molded components (Purdue) • Meso-Scale morphological analysis and correlation with structural performance (UTK) • Paint and adhesion durability analysis (MSU) • Demonstrator part design, prototype production, and validation testing (Ford). The primary goal at the commencement of the project was development of a chopped carbon fiber sheet molding compound (SMC) that offered a three times improvement in tensile modulus over a comparable glass-based SMC. In addition to meeting mechanical performance targets, the resin kinetics were modified to achieve a processing cycle time of less than 3 minutes. Other critical-to-quality (CTQ) specifications were also stipulated to account for a broad range of materials and processing characteristics. To achieve the above, staff scientists at Dow Chemical created an extensive series of new epoxy blends for testing and validation. Throughout this development, a key challenge was attaining material performance goals without impacting processing behavior and paintability of finished components. The latter required a new internal mold release system being developed by Dow that was designed to complement the kinetics of the rapid cure epoxy. As a complement to work studies at the industrial partners, the teams from academia executed a series of analytical and experimental studies to investigate potential factors influencing CF-SMC performance. Unit cell models were developed to capture the meso-scale representations of the fiber matrix architecture. Results of this analysis and subsequent morphological investigations led to the design of a novel composite derivative comprising carbon fiber platelets embedded in an epoxy matrix; the platelet size and aspect ratio playing significant role in final composite properties. This approach was a departure from previous research in CF-SMC development whereby bulk filamentization or disassembly of the carbon fiber rovings had been considered the most effective means of achieving both fiber wet through and wet-out. As the course of the academia studies progressed, the aspect ratio of the fiber constituents was further optimized before finalizing material attributes and processing conditions. For the purposes of technology validation, the Ford team led a work stream devoted to the design, fabrication and testing of demonstration components. The carbon fiber SMC material has the potential to displace numerous stampings and castings on an automotive structure but ultimately vehicle closure applications were selected to showcase the abilities of the CF-SMC to achieve both mass reduction and business case for large complex structures. Using target properties established by the Dow staff scientists, the complete closure system for a full-size sedan decklid and a mid-size wagon liftgate were engineered. Prototypes for both applications were fabricated using production representative processing methods to allow for physical testing and performance validation of the CF-SMC structures. Following completion of a testing program that concluded with a FMVSS301 55 mph offset rear crash, the CF-SMC formulation was declared by the Ford team to have met all engineering requirements. To summarize, the joint development activities during this project led to significant technical breakthroughs and achievement of all milestones. The result was the development of a novel, tack-free carbon fiber molding compound that is suited to automated processing. This combined room temperature stability, fast cure kinetics, and internal mold release system facilitates cycle times that are conducive to high-volume production. The VORAFUSE M6400 successfully passed technology validation at Ford and is now eligible for consideration on future production commercial vehicle programs.

36 MATERIALS SCIENCE↗

NASA NTRS2008

Automatic Hazard Detection for Landers

Unmanned planetary landers to date have landed 'blind'; that is, without the benefit of onboard landing hazard detection and avoidance systems. This constrains landing site selection to very benign terrain,which in turn constrains the scientific agenda of missions. The state of the art Entry, Descent, and Landing (EDL) technology can land a spacecraft on Mars somewhere within a 20-100km landing ellipse.Landing ellipses are very likely to contain hazards such as craters, discontinuities, steep slopes, and large rocks, than can cause mission-fatal damage. We briefly review sensor options for landing hazard detection and identify a perception approach based on stereo vision and shadow analysis that addresses the broadest set of missions. Our approach fuses stereo vision and monocular shadow-based rock detection to maximize spacecraft safety. We summarize performance models for slope estimation and rock detection within this approach and validate those models experimentally. Instantiating our model of rock detection reliability for Mars predicts that this approach can reduce the probability of failed landing by at least a factor of 4 in any given terrain. We also describe a rock detector/mapper applied to large-high-resolution images from the Mars Reconnaissance Orbiter (MRO) for landing site characterization and selection for Mars missions.

hazard detection↗

NASA NTRS1981

On eigensolutions for discontinuous liners in a duct containing uniform mean flow

Sound attenuation in a rectangular acoustically lined duct containing uniform mean flow is analytically investigated using the generalized Wiener-Hopf technique. Uniqueness of the solution is enforced for lined sections of the finite axial extent by imposing edge conditions at the liner interface. Possible edge conditions are considered, including the Kutta condition, and the causal solution corresponding to edge conditions is considered the best choice. Solution methods such as the mode matching and singularity methods imply differing edge conditions, and results show that power attenuation is insensitive to the imposed edge conditions, although significant differences are observed for the reflection coefficient. The amplitude of the exponentially increasing instability mode in the lined section must be set to zero as a first approximation to the nonlinear situation, and results indicate that measurements of the reflection factor can be used to make a more definite decision about physically appropriate edge conditions.

Koch, W.↗

NASA NTRS1992

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗

NASA NTRS1988

The mode 3 crack problem in bonded materials with a nonhomogeneous interfacial zone

The mode 3 crack problem for two bonded homogeneous half planes was considered. The interfacial zone was modelled by a nonhomogeneous strip in such a way that the shear modulus is a continuous function throughout the composite medium and has discontinuous derivatives along the boundaries of the interfacial zone. The problem was formulated for cracks perpendicular to the nominal interface and was solved for various crack locations in and around the interfacial region. The asymptotic stress field near the tip of a crack terminating at an interface was examined and it was shown that, unlike the corresponding stress field in piecewise homogeneous materials, in this case the stresses have the standard square root singularity and their angular variation was identical to that of a crack in a homogeneous medium. With application to the subcritical crack growth process in mind, the results given include mostly the stress intensity factors for some typical crack geometries and various material combinations.

Erdogan, Fazil↗

NASA NTRS1991

The mode III crack problem in bonded materials with a nonhomogeneous interfacial zone

Erdogan, F.↗

NASA NTRS1988

JPL pulsar timing observations. V - Macro and microjumps in the Vela pulsar 0833-45

In this paper, 14.5 yr of JPL timing data are analyzed to study the spindown of the Vela pulsar. Discontinuities in the spin frequency and its derivatives have bimodal amplitude distributions. Macrojumps are characterized by relative changes in spin rate of about 10 to the -6th with a derivative of 0.01, while microjumps have relative changes in spin rate of less than about 10 to the -9th with derivative of less than about 0.0001. Macrojumps display a (+,-) signature with a well-defined ratio of the amplitudes. Microjumps, however, show all possible signs of events and a larger range in signature. Six macrojumps and their associated transients are quantified using a 10-parameter model for each jump. Amplitude and decay time parameters vary by a factor of about five over the set of jumps. The results are discussed in terms of superfluid vortex models. The macrojumps are consistent with a model in which angular momentum is suddenly transferred to the crust from a more rapidly rotating superfluid core, followed by relaxation of the spin-down rate.

Cordes, J. M.↗

NASA NTRS1991

Computation of viscous blast wave flowfields

A method to determine unsteady solutions of the Navier-Stokes equations was developed and applied. The structural finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using MacCormack's modified Steger-Warming flux or Roe flux difference splittings with total variation diminishing limiters, while the viscous flux is computed using central differences. The use of implicit boundary conditions in conjunction with a telescoping in time and space method permitted solutions to this strongly unsteady class of problems. Comparisons of numerical, analytical, and experimental results were made in two and three dimensions. These comparisons revealed accurate wave speed resolution with nonoscillatory discontinuity capturing. The purpose of this effort was to address the three-dimensional, viscous blast-wave problem. Test cases were undertaken to reveal these methods' weaknesses in three regimes: (1) viscous-dominated flow; (2) complex unsteady flow; and (3) three-dimensional flow. Comparisons of these computations to analytic and experimental results provided initial validation of the resultant code. Addition details on the numerical method and on the validation can be found in the appendix. Presently, the code is capable of single zone computations with selection of any permutation of solid wall or flow-through boundaries.

Atwood, Christopher A.↗

NASA NTRS1992

Energetic and magnetosheath energy particle signatures of the low-latitude boundary layer at low altitudes near noon

The low-latitude boundary layer (LBL) and its separation from the cusp have previously been identified using observations of particle precipitation at magnetosheath energies. Using S3-3 satellite observations, we have determined that these identifications can also be made from energetic particle observations on polar-orbiting satellites. It is found that the equatorward boundary of the LBL is identifiable as an approximately discontinuous decrease in 33-keV electron fluxes from low to high latitudes. Both the energetic ion and electron fluxes decrease discontinuously at the boundary between the LBL and the cusp or polar cap. A distinct LBL is nearly always identifiable in energetic particle measurements in the 10-14 MLT region when counting rates are statistically significant. The identifications obtained using the energetic particle measurements have been compared to those obtained using criteria developed by Newell and Meng (1988, 1989) for magnetosheath energy particle precipitation. In this way, we have evaluated the accuracy of both techniques and used the energetic particle measurements to supplement the identifications obtained using the Newell and Meng criteria. We propose that the Newell and Meng threshold on ion energy flux can be reduced by a factor of 6. This modification provides identification of the LBL for lower ion intensity levels than has previously been thought possible. Source, acceleration, and scattering processes have also been studied within and in the vicinity of the LBL. Observed trapped pitch angle distributions of energetic electrons imply that the LBL is at least partially on closed field lines. Strong scattering of energetic protons is found within and equatorward of the LBL and thus must occur at least partially along closed field lines. Field-aligned electron acceleration by parallel electric fields can be discerned within and poleward of the LBL, but a more detailed analysis is necessary for a statistical study. Conical ion acceleration was seen relatively frequently within the LBL and about half as often poleward of the LBL. Neither acceleration process could be identified anywhere equatorward of the LBL.

Roeder, J. L.↗

NASA NTRS1998

A New Synoptic Scale Feature of the Auroral Oval: The Nightside Gap

Th Polar Ultraviolet Imager (UVI) has recently identified a new morphological feature of the nightside aurora in which a well defined segment of the pre-midnight auroral oval, typically spanning 1.0- 1.5 hours in local time (710-1100 km), exhibits a significant reduction in luminosity relative to the active regions surrounding it. This feature appears in roughly 7% of the substorm period observed during our initial study spanning December 1996 through February, 1997. The energy flux (erg/sq cm/s) is shown to be lower by as much as a factor of four in this nightside gaps in comparison to the adjacent auroral activity. Substorm intensifications are often observed on either side of the nightside gap, which remains relatively inactive during these substorm periods. We demonstrate that the nightside gap is sometimes terminated by the onset of a substorm directly within the same local time boundaries. The nightside gap has strong implications for the dynamics of magnetosphere-ionosphere coupling and current closure in the midnight sector auroral oval. Specifically, we conjecture that the nightside gap occurs when transverse ionospheric currents do not feed into upward field-aligned currents in the midnight sector, perhaps due to the Harang discontinuity.

Cliua, D.↗

DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

DOE OSTI2021

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗

DOE OSTI2026

Unveiling Swift Heavy Ion Track Morphology in Sr-Based High-Entropy Perovskites

The incorporation of multiple cations on a single lattice site in the high-entropy oxides is considered the key driving factor for modifying the known atomic-level response to energetic ion irradiation due to the presence of structural disorder; however, these effects are not well-understood yet. In this work, we present atomic-level insight into irradiation-induced nanoscale phase transformations in a perovskite-structured high-entropy oxide, Sr(Zr 0.2 Sn 0.2 Ti 0.2 Hf 0.2 Nb 0.2 )O 3 (Sr(HE)O 3 ), subjected to 774 MeV swift Xe heavy ions, where damage is dominated by inelastic ion−lattice interactions. While these ions generally are known to create nanoscale disordered channels, “ion tracks”, along the penetration direction in the material, this study shows the formation of discontinuous and partially recrystallized ion tracks in Sr(HE)O 3 . Compared to SrTiO 3 irradiated under identical energy loss conditions, the ion tracks in Sr(HE)O 3 exhibit significantly reduced diameters and a markedly different interfacial structure. Notably, the crystalline−amorphous interface in Sr(HE)O 3 shows minimal lattice distortion, confined to approximately 2−3 monolayers, in contrast to the extended disordered shell commonly observed in SrTiO 3 . Using in situ atomic-resolution electron microscopy, we further demonstrate that the amorphous/disordered regions within Sr(HE)O 3 ion tracks remain highly stable under electron irradiation, whereas tracks in SrTiO 3 readily recrystallize. This enhanced stability is attributed to the dominance of structural and chemical complexity arising from multiple B-site cations, which suppress defect migration and templated recrystallization driven by electronic excitations and local heating. Overall, this study highlights how high-entropy oxide chemistry fundamentally reshapes irradiation damage evolution, offering insights into defect formation and phase stability under extreme conditions.

36 MATERIALS SCIENCE↗

NASA NTRS1990

High order filtering methods for approximating hyberbolic systems of conservation laws

In the computation of discontinuous solutions of hyperbolic systems of conservation laws, the recently developed essentially non-oscillatory (ENO) schemes appear to be very useful. However, they are computationally costly compared to simple central difference methods. A filtering method which is developed uses simple central differencing of arbitrarily high order accuracy, except when a novel local test indicates the development of spurious oscillations. At these points, the full ENO apparatus is used, maintaining the high order of accuracy, but removing spurious oscillations. Numerical results indicate the success of the method. High order of accuracy was obtained in regions of smooth flow without spurious oscillations for a wide range of problems and a significant speed up of generally a factor of almost three over the full ENO method.

Lafon, F.↗

NASA NTRS1995

Nutritional value of protein hydrolysis products (oligopeptides and free amino acids) as a consequence of absorption and metabolism kinetics

When pigs were submitted to duodenal infusion of solutions containing a large percentage of small peptides (PEP) or free amino acids with the same pattern (AAL) amino acids appear in the portal blood more rapidly and more uniformly after infusion of PEP then after infusion of AAL, with the notable exception of methionine for which the opposite was true. These differences were lowered when a carbohydrate (maltose dextrin) was present in the solution, but nevertheless remained significant for the first hour after the infusion. The long-term (8-hour) uptake of free amino acids into the liver and the peripheral tissues differed in profile according to the nature of the duodenal infusion. Peripheral uptake was appreciably less well balanced after infusion of free amino acids (deficiency of threonine and phenylalanine) than after infusion of small peptides (deficiency of methionine). Accordingly, in the rat, under conditions of discontinuous enteral nutrition the mixture of small peptides was of greater nutritive value than the mixture of free amino acids. It thus appears that the absorption kinetics which results in important variations in the temporal distribution of free amino acids in the tissues may be at the origin of transitory imbalances in tissue amino acid uptake, and as a result of a lower nutritive value.

Intestinal Absorption/drug effects↗

NASA NTRS1995

The development of concentration gradients in a suspension of chemotactic bacteria

When a suspension of bacterial cells of the species Bacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e., they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development when N and C depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

NASA Discipline Cell Biology↗

DOE OSTI2022

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗