Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface properties”

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 1,531 records · Page 85

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sintering time and atmosphere influences on the microstructure and mechanical properties of tungsten heavy alloys

The present paper traces the microstructure and mechanical property variations vs sintering time under conditions where the atmosphere was adjusted to suppress pore generation. Attention is given to tracing the mechanical property links to the grain size under optimized sintering conditions. The optimal cycle is used to sinter two heavy alloy compositions (88 and 95 wt pct W) for times up to 600 minutes at 1480 C. The 88-pct-W samples slumped, but the 95-pct-W samples were fully densified and suitable for tensile testing. At long sintering times, the tungsten grains flattened and the tungsten contiguity decreased, indicating a transition to low-energy configurations for the solid-liquid interfaces. The cube of the mean grain size varied linearly with the isothermal sintering time. This allowed the determination of grain size effects on mechanical properties, showing a decreasing yield strength with increasing time in agreement with the Hall-Petch behavior. The tensile strength and elongation were highest for sintering times from 30 to 90 min, reflecting a minimum in the residual porosity.

German, R. M.↗

In-situ resonant band engineering of solution-processed semiconductors generates high performance n-type thermoelectric nano-inks

Thermoelectric devices possess enormous potential to reshape the global energy landscape by converting waste heat into electricity, yet their commercial implementation has been limited by their high cost to output power ratio. No single “champion” thermoelectric material exists due to a broad range of material-dependent thermal and electrical property optimization challenges. While the advent of nanostructuring provided a general design paradigm for reducing material thermal conductivities, there exists no analogous strategy for homogeneous, precise doping of materials. Here, we demonstrate a nanoscale interface-engineering approach that harnesses the large chemically accessible surface areas of nanomaterials to yield massive, finely-controlled, and stable changes in the Seebeck coefficient, switching a poor nonconventional p- type thermoelectric material, tellurium, into a robust n -type material exhibiting stable properties over months of testing. These remodeled, n -type nanowires display extremely high power factors (~500 µW m -1 K -2 ) that are orders of magnitude higher than their bulk p -type counterparts.

42 ENGINEERING↗

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗

Experimental study of the minority-carrier transport at the polysilicon-monosilicon interface

This paper presents the results of an experimental study designed to explore both qualitatively and quantitatively the mechanism of the improved current gain in bipolar transistors with polysilicon emitter contacts. Polysilicon contacts were deposited and heat treated at different conditions. The electrical properties were measured using p-n junction test structures that are much more sensitive to the contact properties than are bipolar transistors. A simple phenomenological model was used to correlate the structural properties with electrical measurements. Possible transport mechanisms are examined and estimates are made about upper bounds on transport parameters in the principal regions of the devices. The main conclusion of this study is that the minority-carrier transport in the polycrystalline silicon is dominated by a highly disordered layer at the polysilicon-monosilicon interface characterized by very low minority-carrier mobility. The effective recombination velocity at the n(+) polysilicon-n(+) monosilicon interface was found to be a strong function of fabrication conditions. The results indicate that the recombination velocity can be much smaller than 10,000 cm/s.

Neugroschel, A.↗

Engineering metal–metal oxide surfaces for high-performance oxygen reduction on Ag–Mn electrocatalysts

Understanding fundamental material–property relationships in mixed-element catalyst systems is crucial to advancing the viability of renewable electrochemical energy technologies, an important part of creating a more sustainable future. Herein, we report our insight on the nature and dynamics of highly active silver–manganese oxide (Ag–MnOx) catalyst surfaces for the oxygen reduction reaction (ORR) via a combined experimental–theoretical approach. Experimentally, we synthesize well-mixed Ag–Mn co-deposited thin films that are measurably flat and smooth, despite Mn surface migration and oxidation upon air exposure and electrochemical measurements. Cyclic voltammetry in 0.1 M KOH demonstrates up to 10-fold specific activity enhancements over pure Ag at 0.8 V vs. RHE for Ag-rich films (70–95% Ag in bulk). To further understand the Ag–Mn system, separate samples were synthesized with small amounts of Mn sequentially deposited onto the surface of a pure Ag thin film (Mn@Ag), ranging from partial to full surface coverage (down to 0.3 nm Mn $cm^{-2}_{geo}$ ~ 0.2 μg Mn $cm^{-2}_{geo}$). These sequentially deposited Mn@Ag films show analogous performance to their co-deposited counterparts indicating similar enhanced active sites. With density functional theory (DFT), we calculate that this enhancement arises from the tuned d-band of these material surfaces owing to the optimal hybridization of the electronic structures in specific Ag and MnO x geometries. Together, electrochemical measurements, DFT calculations, X-ray absorption spectroscopy, and valence-band X-ray photoelectron spectroscopy suggest synergistic electronic interactions between Ag and MnOx yield enhanced oxygen adsorption, and thus ORR activity, with DFT highlighting the Ag–MnO x interface sites as the most enhanced. This work demonstrates how combined experimental–theoretical methods can help design electrocatalysts with enhanced electrocatalytic properties and understand the nature of complex mixed metal–metal oxide surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermo-Mechanical Properties of SiC/SiC Composites with Hybrid CVI-PIP Matrices

For long term structural service, the upper temperature capability for slurry-cast melt infiltrated (MI) SiC/SiC composites is limited to approx. 1315 C because of silicon reaction with the SiC fibers. For applications requiring material temperatures in excess of 1315 C, alternate methods of manufacturing the SiC matrices without silicon are being investigated, such as a hybrid combination of CVI and PIP. In this study, stacked fabric plies of Sylramic i-BN SiC fibers were coated with a CVI BN interface layer followed by a partial CVI SiC matrix. The remaining porosity in the SiC/SiC preforms was then infiltrated with silicon carbide matrix by PIP. Thermo-mechanical property measurements indicate that these composites are stable to 1700 C in inert environments under no load conditions for 100 h and under load conditions to 1450 C in air for 300 h. The advantages, disadvantages, and potential of this composite system for high temperature applications will be discussed.

Bhatt, R. T.↗

MMIC Amplifiers and Wafer Probes for 350 to 500 GHz

A wireless avionics interface exploits the constrained nature of data networks in flight systems to use a lightweight routing method. This simplified routing means that a processor is not required, and the logic can be implemented as an intellectual property (IP) core in a field-programmable gate array (FPGA). The FPGA can be shared with the flight subsystem application. In addition, the router is aware of redundant subsystems, and can be configured to provide hot standby support as part of the interface. This simplifies implementation of flight applications requiring hot stand - by support. When a valid inbound packet is received from the network, the destination node address is inspected to determine whether the packet is to be processed by this node. Each node has routing tables for the next neighbor node to guide the packet to the destination node. If it is to be processed, the final packet destination is inspected to determine whether the packet is to be forwarded to another node, or routed locally. If the packet is local, it is sent to an Applications Data Interface (ADI), which is attached to a local flight application. Under this scheme, an interface can support many applications in a subsystem supporting a high level of subsystem integration. If the packet is to be forwarded to another node, it is sent to the outbound packet router. The outbound packet router receives packets from an ADI or a packet to be forwarded. It then uses a lookup table to determine the next destination for the packet. Upon detecting a remote subsystem failure, the routing table can be updated to autonomously bypass the failed subsystem.

Samoska, Lorene A.↗

Interface Supports Lightweight Subsystem Routing for Flight Applications

A wireless avionics interface exploits the constrained nature of data networks in flight systems to use a lightweight routing method. This simplified routing means that a processor is not required, and the logic can be implemented as an intellectual property (IP) core in a field-programmable gate array (FPGA). The FPGA can be shared with the flight subsystem application. In addition, the router is aware of redundant subsystems, and can be configured to provide hot standby support as part of the interface. This simplifies implementation of flight applications requiring hot stand - by support. When a valid inbound packet is received from the network, the destination node address is inspected to determine whether the packet is to be processed by this node. Each node has routing tables for the next neighbor node to guide the packet to the destination node. If it is to be processed, the final packet destination is inspected to determine whether the packet is to be forwarded to another node, or routed locally. If the packet is local, it is sent to an Applications Data Interface (ADI), which is attached to a local flight application. Under this scheme, an interface can support many applications in a subsystem supporting a high level of subsystem integration. If the packet is to be forwarded to another node, it is sent to the outbound packet router. The outbound packet router receives packets from an ADI or a packet to be forwarded. It then uses a lookup table to determine the next destination for the packet. Upon detecting a remote subsystem failure, the routing table can be updated to autonomously bypass the failed subsystem.

Lux, James P.↗

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗

Nanoconfined Interfaces for Highly Selective Separation of Critical Rare Earth Elements

Industrial demand for rare earth elements (REEs) has surged over the past three decades due to their unique properties that support sustainable energy and new technologies. Separating individual REEs is challenging and hazardous, typically done through liquid-liquid extraction. There is an urgent need for environmentally friendly and efficient separation technologies for REEs. Porous materials offer promising advances for sustainable REE separation via ion-selective capture. We hypothesize that REE separation can be efficiently achieved in reactive nanopores, such as Zr(IV) and Cr(III) metal-organic frameworks (MOFs), through surface functionalization. By integrating material synthesis, interfacial chemistry experiments, theory, computation, and machine learning, we gained insights into the chemical factors controlling REE speciation and their competitive adsorption on MOFs. Our findings show that these materials’ selectivity can be tuned by surface functionalization. The machine learning component addressed ion-specific diffusion based on MOF topology and chemistry.

36 MATERIALS SCIENCE↗

IGES transformer and NURBS in grid generation

In the field of Grid Generation and the CAD/CAM, there are numerous geometry output formats which require the designer to spend a great deal of time manipulating geometrical entities in order to achieve a useful sculptured geometrical description for grid generation. Also in this process, there is a danger of losing fidelity of the geometry under consideration. This stresses the importance of a standard geometry definition for the communication link between varying CAD/CAM and grid system. The IGES (Initial Graphics Exchange Specification) file is a widely used communication between CAD/CAM and the analysis tools. The scientists at NASA Research Centers - including NASA Ames, NASA Langley, NASA Lewis, NASA Marshall - have recognized this importance and, therefore, in 1992 they formed the committee of the 'NASA-IGES' which is the subset of the standard IGES. This committee stresses the importance and encourages the CFD community to use the standard IGES file for the interface between the CAD/CAM and CFD analysis. Also, two of the IGES entities -- the NURBS Curve (Entity 126) and NURBS Surface (Entity 128) -- which have many useful geometric properties -- like the convex hull property, local control property and affine invariance, also widely utilized analytical geometries can be accurately represented using NURBS. This is important in today grid generation tools because of the emphasis of the interactive design. To satisfy the geometry transformation between the CAD/CAM system and Grid Generation field, the CAGI (Computer Aided Geometry Design) developed, which include the Geometry Transformation, Geometry Manipulation and Geometry Generation as well as the user interface. This paper will present the successful development IGES file transformer and application of NURBS definition in the grid generation.

Yu, Tzu-Yi↗

Addition of Short Polymer Chains Mechanically Reinforces Glassy Poly(2-vinylpyridine)–Silica Nanoparticle Nanocomposites

The addition of hard fillers to a polymer matrix is a well-known process for achieving mechanical reinforcement. With a decrease in the size of the fillers, the contribution from polymer–particle nanometer-sized interfaces become significant, and these interfaces affect the mechanical performance of polymer nanocomposites (PNCs) beyond the limits established for conventional composites. However, the molecular mechanisms underlying the improvement in the mechanical performance of glassy PNCs remain unresolved, necessitating a deeper understanding of the structure–property relationships in these intrinsically heterogeneous systems. In this effort, by using Brillouin light scattering (BLS) and dynamic mechanical analysis (DMA), we demonstrated that adding shorter chains to a PNC prepared with high molecular weight polymers significantly improved the mechanical properties of the PNC in the glassy state. The strongest enhancement of mechanical properties occurred at an optimum concentration of short chains. This is in contrast to behavior of the glass transition temperature of PNCs which shows a monotonic decrease with an increase in the concentration of shorter chains. Using experimental data and coarse-grained molecular dynamics (MD) simulations, we have identified the molecular mechanism leading to the observed non-monotonic changes in mechanical reinforcement. Here, this mechanism includes changes in the nanoscale organization at the interface, combined with chain stretching amplified by the addition of the short chains. Overall, our approach paves a simple, cost-effective pathway to fabricating glassy PNCs with significantly improved mechanical properties that will fill various practical needs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-Ray Fluorescence Solvent Detection at the Substrate-Adhesive Interface

With environmental regulations limiting the use of volatile organic compounds, low-vapor pressure solvents have replaced traditional degreasing solvents for bond substrate preparation. When used to clean and prepare porous bond substrates such as phenolic composites, low vapor pressure solvents can penetrate deep into substrate pore networks and remain there for extended periods. Trapped solvents can interact with applied adhesives either prior to or during cure, potentially compromising bond properties. Currently, methods for characterizing solvent time-depth profiles in bond substrates are limited to bulk gravimetric or sectioning techniques. While sectioning techniques such as microtome allow construction of solvent depth profiles, their depth resolution and reliability are limited by substrate type. Sectioning techniques are particularly limited near the adhesive-substrate interface where depth resolution is further limited by adhesive-substrate hardness and, in the case of a partially cured adhesive, mechanical properties differences. Additionally, sectioning techniques cannot provide information about lateral solvent diffusion. Cross-section component mapping is an alternative method for measuring solvent migration in porous substrates that eliminates the issues associated with sectioning techniques. With cross-section mapping, the solvent-wiped substrate is sectioned perpendicular rather than parallel to the wiped surface, and the sectioned surface is analyzed for the solvent or solvent components of interest using a two-dimensional mapping or imaging technique. Solvent mapping can be performed using either direct or indirect methods. With a direct method, one or more solvent components are mapped using red or Raman spectroscopy together with a moveable sample stage and/or focal plane array detector. With an indirect method, an elemental "tag" not present in the substrate is added to the solvent before the substrate is wiped. Following cross sectioning, the tag element can then be mapped by its characteristic x-ray emission using either x-ray fluorescence, or electron-beam energy-and wavelength-dispersive x-ray spectrometry. The direct mapping techniques avoid issues of different diffusion or migration rates of solvents and elemental tags, while the indirect techniques avoid spectral resolution issues in cases where solvents and substrates have adjacent or overlapping peaks. In this study, cross-section component indirect mapping is being evaluated as a method for measuring migration of d-limonene based solvents in glass-cloth phenolic composite (GCP) prior to and during subsequent bonding and epoxy adhesive cure.

Wurth, Laura↗

Entropy Stable h/p-Nonconforming Discretization with the Summation-by-Parts Property for the Compressible Euler and Navier–Stokes Equations

In this paper, we extend the entropy conservative/stable algorithms presented by Del Rey Fernandez and coauthors for the compressible Euler and Navier-Stokes equations on nonconforming p-refined/coarsened curvilinear grids to h/p refinement/coarsening. The main difficulty in developing nonconforming algorithms is the construction of appropriate coupling procedures across nonconforming interfaces. Here, we utilize a computationally simple and efficient approach based upon using decoupled interpolation operators. The resulting scheme is entropy conservative/stable and elementwise conservative. Numerical simulations of the isentropic vortex and viscous shock propagation con firm the entropy conservation/stability and accuracy properties of the method (achieving ~ p + 1 convergence), which are comparable to those of the original conforming scheme. Simulations of the Taylor{Green vortex at R(e) = 1,600 and turbulent flow past a sphere at R(e(infinity)) = 2,000 show the robustness and stability properties of the overall spatial discretization for unstructured grids. Finally, to demonstrate the entropy conservation property of a fully-discrete explicit entropy stable algorithm with h=p refinement/coarsening, we present the time evolution of the entropy function obtained by simulating the propagation of the isentropic vortex using a relaxation Runge-Kutta scheme.

Nonconforming interfaces↗

A Self‐Healing System for Polydicyclopentadiene Thermosets

Self-healing offers promise for addressing structural failures, increasing lifespan, and improving durability in polymeric materials. Implementing self-healing in thermoset polymers faces significant manufacturing challenges, especially due to the elevated temperature requirements of thermoset processing. To introduce self-healing into structural thermosets, the self-healing system must be thermally stable and compatible with the thermoset chemistry. This article demonstrates a self-healing microcapsule-based system stable to frontal polymerization (FP), a rapid and energy-efficient manufacturing process with a self-propagating exothermic reaction (≈200 °C). A thermally latent Grubbs-type complex bearing two N-heterocyclic carbene ligands addresses limitations in conventional G2-based self-healing approaches. Under FP's elevated temperatures, the catalyst remains dormant until activated by a Cu(I) co-reagent, ensuring efficient polymerization of the dicyclopentadiene (DCPD) upon damage to the polyDCPD matrix. The two-part microcapsule system consists of one capsule containing the thermally latent Grubbs-type catalyst dissolved in the solvent, and another capsule containing a Cu(I) coagent blended with liquid DCPD monomer. Using the same chemistry for both matrix fabrication and healing results in strong interfaces as demonstrated by lap-shear tests. In an optimized system, the self-healing system restores the mechanical properties of the tough polyDCPD thermoset. Self-healing efficiencies greater than 90% via tapered double cantilever beam tests are observed.

36 MATERIALS SCIENCE↗

The Role of Interfacial Interactions and Oxygen Vacancies in Tuning Magnetic Anisotropy in LaCrO 3 /LaMnO 3 Heterostructures

The interplay of lattice, electronic, and spin degrees of freedom at epitaxial complex oxide interfaces provides a route to tune their magnetic ground states. Unraveling the competing contributions is critical for tuning their functional properties. The relationship between magnetic ordering and magnetic anisotropy and the lattice symmetry, oxygen content, and film thickness in compressively strained LaMnO 3 (LMO)/LaCrO 3 (LCO) superlattices is investigated. Mn–O–Cr antiferromagnetic superexchange interactions across the heterointerface result in a net ferrimagnetic magnetic structure. Bulk magnetometry measurements reveal isotropic in-plane magnetism for as-grown oxygen-deficient thin samples due to equal fractions of orthorhombic a+a-c-, and a-a+c- twin domains. As the superlattice thickness is increased, in-plane magnetic anisotropy emerges as the fraction of the a+a-c- domain increases. On annealing in oxygen, the suppression of oxygen vacancies results in a contraction of the lattice volume, and an orthorhombic to rhombohedral transition leads to isotropic magnetism independent of the film thickness. The complex interactions are investigated using high-resolution synchrotron diffraction and X-ray absorption spectroscopy. These results highlight the role of the evolution of structural domains with film thickness, interfacial spin interactions, and oxygen-vacancy-induced structural phase transitions in tuning the magnetic properties of complex oxide heterostructures.

36 MATERIALS SCIENCE↗

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗