Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

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 721 records · Page 40

Characterizing Physical Properties of Streambed Interface Sediments Using In Situ Complex Electrical Conductivity Measurements

Streambed sediment physical properties such as surface area, are difficult to quantify in situ but exert a high-level control on a wide range of biogeochemical processes and sorption of contaminants. We introduce the use of complex electrical conductivity (CC) methods (also known as spectral-induced polarization (SIP)) that measure both real and imaginary conductivity to noninvasively and efficiently characterize shallow streambed sediments. We explore the method through synthetic modeling, laboratory, and field measurements to demonstrate the sensitivity of imaginary conductivity to sediment surface area, controlled in part by fine-grained iron oxides produced by anoxic groundwater discharge. Laboratory measurements verify expected relationships between CC parameters and sediment properties. Synthetic modeling using a 1D analytical model illustrates the influence of water layer depth and conductivity on the field CC measurements made at the streambed-stream water interface. Specifically, the inverted sediment imaginary conductivity is less impacted by uncertainty in the water layer depth and conductivity relative to the real conductivity and phase shift. Field CC measurements along a landfill-impacted river reveal discrete streambed zones with enhanced bulk surface area generally corresponding to anoxic groundwater discharge zones with high concentrations of fine-grained iron oxide precipitates.

54 ENVIRONMENTAL SCIENCES↗

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗

Assembly of JAZ–JAZ and JAZ–NINJA complexes in jasmonate signaling

Jasmonates (JAs) are plant hormones with crucial roles in development and stress resilience. They activate MYC transcription factors by mediating the proteolysis of MYC inhibitors called JAZ proteins. In the absence of JA, JAZ proteins bind and inhibit MYC through the assembly of MYC–JAZ–Novel Interactor of JAZ (NINJA)–TPL repressor complexes. However, JAZ and NINJA are predicted to be largely intrinsically unstructured, which has precluded their experimental structure determination. Through a combination of biochemical, mutational, and biophysical analyses and AlphaFold-derived ColabFold modeling, we characterized JAZ–JAZ and JAZ–NINJA interactions and generated models with detailed, high-confidence domain interfaces. We demonstrate that JAZ, NINJA, and MYC interface domains are dynamic in isolation and become stabilized in a stepwise order upon complex assembly. By contrast, most JAZ and NINJA regions outside of the interfaces remain highly dynamic and cannot be modeled in a single conformation. Our data indicate that the small JAZ Zinc finger expressed in Inflorescence Meristem (ZIM) motif mediates JAZ–JAZ and JAZ–NINJA interactions through separate surfaces, and our data further suggest that NINJA modulates JAZ dimerization. This study advances our understanding of JA signaling by providing insights into the dynamics, interactions, and structure of the JAZ–NINJA core of the JA repressor complex.

59 BASIC BIOLOGICAL SCIENCES↗

The Long-Term Evolution of the Atmosphere of Venus: Processes and Feedback Mechanisms

This work reviews the long-term evolution of the atmosphere of Venus, and modulation of its composition by interior/exterior cycling. The formation and evolution of Venus’s atmosphere, leading to contemporary surface conditions, remain hotly debated topics, and involve questions that tie into many disciplines. We explore these various inter-related mechanisms which shaped the evolution of the atmosphere, starting with the volatile sources and sinks. Going from the deep interior to the top of the atmosphere, we describe volcanic out-gassing, surface-atmosphere interactions, and atmosphere escape. Furthermore, we address more complex aspects of the history of Venus, including the role of Late Accretion impacts, how magnetic field generation is tied into long-term evolution, and the implications of geochemical and geodynamical feedback cycles for atmospheric evolution. We highlight plausible end-member evolutionary pathways that Venus could have followed, from accretion to its present-day state, based on modeling and observations. In a first scenario, the planet was desiccated by atmospheric escape during the magma ocean phase. In a second scenario, Venus could have harbored surface liquid water for long periods of time, until its temperate climate was destabilized and it entered a runaway greenhouse phase. In a third scenario, Venus’s inefficient outgassing could have kept water inside the planet, where hy- drogen was trapped in the core and the mantle was oxidized. We discuss existing evidence and future observations/missions required to refine our understanding of the planet’s history and of the complex feedback cycles between the interior, surface, and atmosphere that have been operating in the past, present or future of Venus.

Venus↗

Using Response Surface Methods to Correlate the Modal Test of an Inflatable Test Article

This paper presents a practical application of response surface methods (RSM) to correlate a finite element model of a structural modal test. The test article is a quasi‐cylindrical inflatable structure which primarily consists of a fabric weave, with an internal bladder and metallic bulkheads on either end. To mitigate model size, the fabric weave was simplified by representing it with shell elements. The task at hand is to represent the material behavior of the weave. The success of the model correlation is measured by comparing the four major modal frequencies of the analysis model to the four major modal frequencies of the test article. Given that only individual strap material properties were provided and material properties of the overall weave were not available, defining the material properties of the finite element model became very complex. First it was necessary to determine which material properties (modulus of elasticity in the hoop and longitudinal directions, shear modulus, Poisson's ratio, etc.) affected the modal frequencies. Then a Latin Hypercube of the parameter space was created to form an efficiently distributed finite case set. Each case was then analyzed with the results input into RSM. In the resulting response surface it was possible to see how each material parameter affected the modal frequencies of the analysis model. If the modal frequencies of the analysis model and its corresponding parameters match the test with acceptable accuracy, it can be said that the model correlation is successful.

Gupta, Anju↗

Toward an Integrated Solution to Mitigate the Impact of Volcanic Ash to Aviation

The science community is making a concerted effort to improve the reliability of dispersion models for the forecasting of volcanic ash plumes. Toward this end, it has been observed that the assimilation of diverse, accurate and frequent surface, airborne and satellite observations of the source and distal ash plumes may hold the key. Various international research organizations and operational agencies make these observations using a variety of active and passive remote sensing systems and use them to initialize atmospheric trajectory and dispersion models. These observation systems range from surface LIDAR and ceilometers, to airborne radiometers and nephelometers, to satellite radiometers, multi-spectral imagers, LIDAR and UV-photometers. None of these systems alone is a panacea, however, their synergistic application holds great promise toward solving this complex problem. Additionally, the aeronautical and science communities are working to better understand the quantitative thresholds and tolerances of aviation systems to volcanic ash to better inform scientists of the accuracy requirements for dispersion model forecasts. A number of the most recent and promising efforts in all of these area are discussed in this presentation.

Murray, John J.↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗

Structural models of the MscL gating mechanism

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

NASA Discipline Cell Biology↗

Deciphering the capricious precipitation response: irrigation impact in the North China Plain

Intensive irrigation in the North China Plain (NCP) raises environmental concerns, yet its climate impacts remain inconsistent among current modeling studies. These inconsistencies may arise from deficiencies in existing models. To better capture irrigation impact, we employed a recently improved irrigation model tailored for the NCP, which significantly enhances simulations of crop growth and irrigation application. Results reveal that irrigation exerts competing effects on precipitation: enhancing it in May and June by increasing rainfall frequency and intensity, while suppressing it in July and August by reducing intensity. Spatial precipitation changes correlate with upper-tropospheric humidity, while inter-annual changes link to geopotential height and wind anomalies. This suggests that irrigation-induced surface cooling and moisturizing modulate atmospheric thermodynamic and dynamic structures, influencing convective precipitation unevenly across regions and time, leading to complex precipitation changes. These findings highlight the complexity of irrigation impacts and emphasize the need for improved land representation in climate models.

54 ENVIRONMENTAL SCIENCES↗

Discrete Element Modeling of Triboelectrically Charged Particles

Tribocharging of particles is common in many processes including fine powder handling and mixing, printer toner transport and dust extraction. In a lunar environment with its high vacuum and lack of water, electrostatic forces are an important factor to consider when designing and operating equipment. Dust mitigation and management is critical to safe and predictable performance of people and equipment. The extreme nature of lunar conditions makes it difficult and costly to carry out experiments on earth which are necessary to better understand how particles gather and transfer charge between each other and with equipment surfaces. DEM (Discrete Element Modeling) provides an excellent virtual laboratory for studying tribocharging of particles as well as for design of devices for dust mitigation and for other purposes related to handling and processing of lunar regolith. Theoretical and experimental work has been performed pursuant to incorporating screened Coulombic electrostatic forces into EDEM, a commercial DEM software package. The DEM software is used to model the trajectories of large numbers of particles for industrial particulate handling and processing applications and can be coupled with other solvers and numerical models to calculate particle interaction with surrounding media and force fields. While simple Coulombic force between two particles is well understood, its operation in an ensemble of particles is more complex. When the tribocharging of particles and surfaces due to frictional contact is also considered, it is necessary to consider longer range of interaction of particles in response to electrostatic charging. The standard DEM algorithm accounts for particle mechanical properties and inertia as a function of particle shape and mass. If fluid drag is neglected, then particle dynamics are governed by contact between particles, between particles and equipment surfaces and gravity forces. Consideration of particle charge and any tribocharging and electric field effects requires calculation of the forces due to these effects.

Hogue, Michael D.↗

Organic-rich fluids in Pluto’s interior—A source of cryovolcanic activity

The highly variegated surface of Pluto revealed by the New Horizons spacecraft is composed of both volatile ices (N2, CH4, CO) and H2O ice, all of which are intrinsically colorless. The colors on Pluto reveal a non-ice component generally acknowledged to be a refractory complex organic material (tholins) produced by photolysis in the planet’s atmosphere [1], and by photolysis and radiolysis of the surface ices [2-4]. These processes have produced a range of colors from pale yellow to red to brown, and result in a variation of a factor of 10in albedo across the surface. Tholins are disordered polymer-like materials made of repeating chains of linked subunits and complex combinations of functional groups containing carbon. Nitrogen, oxygen, and other elements can be substituted in both the aliphatic and aromatic subunits. The detection of the spectral absorption of an ammoniated compound associated with surface exposures of red-colored H2O ice (herein, RAW) in the LEISA spectral mapping data from New Horizons suggests that a third source of a colored material is ejected from Pluto’s interior as a water-rich fluid(cryomagma) that rapidly freezes at the low temperature (~40K) of the surface [5,6]. The ammoniated compound is thought to be an ammonia hydrate or an ammoniated salt; the spectral data cover a limited wavelength range, resulting in ambiguity in the identification of the form of the ammonia. The distribution of RAW in the vicinity of tectonically stressed and deformed geological structures, particularly graben (fossae), is seen as the result of both fluid effusion and fountaining ejection of cryoclastic materials in Pluto’s Virgil Fossae region [5,6]. RAW is seen elsewhere on Pluto, both in association with graben complexes and on broad expanses of old surface, but primarily in zones of tectonic stress. The nature of the putative subsurface fluid is unknown, but it must have existed at relatively shallow depths such that the graben faulting could reach it. Geochemical modeling of fluids in the interiors of small planetary bodies in the outer Solar System[7] shows that ammonia-bearing H2O in contact with rock and metal components produce a rich chemical soup, including gases that can, in principle, help propel the fluid from shallow reservoirs through crustal fractures to reach the surface. An organic component of the fluid cryomagma is predicted to originate from thermal processing of the components [ref] and may be an example of one-pot synthesis of such complexes as amino acid precursors [ref].

Dale P Cruikshank↗

Evaluating the Impact of Chemical Complexity and Horizontal Resolution on Tropospheric Ozone Over the Conterminous US With a Global Variable Resolution Chemistry Model

A new configuration of the Community Earth System Model (CESM)/Community Atmosphere Model with full chemistry (CAM-chem) supporting the capability of horizontal mesh refinement through the use of the spectral element (SE) dynamical core is developed and called CESM/CAM-chem-SE. Horizontal mesh refinement in CESM/CAM-chem-SE is unique and novel in that pollutants such as ozone are accurately represented at human exposure relevant scales while also directly including global feedbacks. CESM/CAM-chem-SE with mesh refinement down to ∼14 km over the conterminous US (CONUS) is the beginning of the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICAv0). Here MUSICAv0 is evaluated and used to better understand how horizontal resolution and chemical complexity impact ozone and ozone precursors over CONUS as compared to measurements from five aircraft campaigns, which occurred in 2013. This field campaign analysis demonstrates the importance of using finer horizontal resolution to accurately simulate ozone precursors such as nitrogen oxides and carbon monoxide. In general, the impact of using more complex chemistry on ozone and other oxidation products is more pronounced when using finer horizontal resolution where a larger number of chemical regimes are resolved. Large model biases for ozone near the surface remain in the Southeast US as compared to the aircraft observations even with updated chemistry and finer horizontal resolution. This suggests a need for adding the capability of replacing sections of global emission inventories with regional inventories, increasing the vertical resolution in the planetary boundary layer, and reducing model biases in meteorological variables such as temperature and clouds.

MUSICA↗

On-Surface Reactions of Electronically Active Self-Assembled Monolayers for Electrode Work Function Tuning

Self-assembled monolayers (SAMs) help improve the performance of organic electronic devices through interface passivation and enhanced carrier transport. Yet, there is limited information regarding the chemical structure of the SAMs upon functionalization and subsequent thermal treatment. Here, we studied the on-surface reaction of carbazole-derived SAMs on model gold electrodes, focusing on the chemical structure changes induced by thermal treatments. Furthermore, we correlate the microscopic changes with their impact on the electrode’s work function. The carbazole-based SAMs first transform into organometallic complexes. At higher annealing temperatures, SAMs convert to oligomeric complexes. The observed chemical reactions significantly reduce the electrode work function and facilitate electron injection in n-type organic thin-film transistors. Our results highlight the on-surface synthesis of electronically active SAMs as an alternative approach for modifying the work function of electrodes for organic electronics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States

Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.

Geochemistry & Geophysics↗

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗