SEARCH · Engineering Papers

Results for “fault 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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Application of Hapke photometric model to three geologic surfaces using PARABOLA bidirectional reflection data

The Geologic Remote Sensing Field Experiment (GRSFE) was conducted in July and September of 1989 to collect data with both ground and airborne instrumentation. A major objective of GRSFE was to collect data which could be used to test radiative transfer models for the extraction of composition and textural surface properties from remotely acquired data. Reported here are the initial results from an application of the Hapke photometric model, using data from the Portable Apparatus for Remote Acquisition of Bidirectional Observations of Land and Atmosphere (PARABOLA), a ground based radiometer with three spectral channels. PARABOLA data was collected in the Lunar Crater Volcanic Field in Nevada, specifically from the region of Lunar Lake, a playa. The Hapke model was found to be inadequate for three relatively common geologic surfaces (a clay-rich, hard packed surface with decimeter sized mudcracks; a cobble site, similar to a playa site, but strewn with basaltic cobbles and pebbles; and a surface mantled basalt lava flow). The model is not at fault; rather, the complexity of most geologic surfaces is not accounted for in the initial assumptions.

Shepard, Michael K.↗

NASA NTRS2015

Towards Real-Time, On-Board, Hardware-Supported Sensor and Software Health Management for Unmanned Aerial Systems

For unmanned aerial systems (UAS) to be successfully deployed and integrated within the national airspace, it is imperative that they possess the capability to effectively complete their missions without compromising the safety of other aircraft, as well as persons and property on the ground. This necessity creates a natural requirement for UAS that can respond to uncertain environmental conditions and emergent failures in real-time, with robustness and resilience close enough to those of manned systems. We introduce a system that meets this requirement with the design of a real-time onboard system health management (SHM) capability to continuously monitor sensors, software, and hardware components. This system can detect and diagnose failures and violations of safety or performance rules during the flight of a UAS. Our approach to SHM is three-pronged, providing: (1) real-time monitoring of sensor and software signals; (2) signal analysis, preprocessing, and advanced on-the-fly temporal and Bayesian probabilistic fault diagnosis; and (3) an unobtrusive, lightweight, read-only, low-power realization using Field Programmable Gate Arrays (FPGAs) that avoids overburdening limited computing resources or costly re-certification of flight software. We call this approach rt-R2U2, a name derived from its requirements. Our implementation provides a novel approach of combining modular building blocks, integrating responsive runtime monitoring of temporal logic system safety requirements with model-based diagnosis and Bayesian network-based probabilistic analysis. We demonstrate this approach using actual flight data from the NASA Swift UAS.

Unmanned Aerial System↗

DOE OSTI2024

Structural and Electronic Properties of Indium-Doped n-type Cd-Se-Te Crystals

Here, we present a comprehensive investigation into the potential of n-type indium-doped cadmium selenide telluride (CST:In) as a high-performance candidate for solar cell applications, without the need for resource-intensive post-growth treatments that are required for CdTe:In. We compared undoped CST and CST:In crystals under different growth conditions, analyzing their structural and electronic properties using x-ray diffraction (XRD), electron probe microanalysis (EPMA), current-voltage (IV) and Hall effect measurements, time-resolved photoluminescence (TRPL), optical transmission, and photoluminescence (PL) mapping. The results reveal that as-grown CST:In crystals achieve nearly 100% carrier activation, yielding an electron concentration of 9.5x1018 cm -3 , mobility of 653 cm 2 /V.s and a 5 ns lifetime which approaches the radiative limit. Furthermore, comparison of PL maps from crystal growths having different cooling profiles suggests a strong effect of cooling rate on selenium segregation and cubic/hexagonal/polytype phase distribution. Slower cooling leads to a more homogeneous cubic structure with lower Se segregation, while a faster cooling rate results in increased Se segregation, and twin boundaries and stacking faults with polytypic and hexagonal character.

14 SOLAR ENERGY↗

DOE OSTI2021

Determining quantum phase diagrams of topological Kitaev-inspired models on NISQ quantum hardware

Topological protection is employed in fault-tolerant error correction and in developing quantum algorithms with topological qubits. But, topological protection intrinsic to models being simulated , also robustly protects calculations, even on NISQ hardware. We leverage it by simulating Kitaev-inspired models on IBM quantum computers and accurately determining their phase diagrams. This requires constructing conventional quantum circuits for Majorana braiding to prepare the ground states of Kitaev-inspired models. The entanglement entropy is then measured to calculate the quantum phase boundaries. We show how maintaining particle-hole symmetry when sampling through the Brillouin zone is critical to obtaining high accuracy. This work illustrates how topological protection intrinsic to a quantum model can be employed to perform robust calculations on NISQ hardware, when one measures the appropriate protected quantum properties. It opens the door for further simulation of topological quantum models on quantum hardware available today.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2011

Relaxing USOS Solar Array Constraints for Russian Vehicle Undocking

With the retirement of Space Shuttle cargo delivery capability and the ten year life extension of the International Space Station (ISS) more emphasis is being put on preservation of the service life of ISS critical components. Current restrictions on the United States Orbital Segment (USOS) Solar Array (SA) positioning during Russian Vehicle (RV) departure from ISS nadir and zenith ports cause SA to be positioned in the plume field of Service Module thrusters and lead to degradation of SAs as well as potential damage to Sun tracking Beta Gimbal Assemblies (BGA). These restrictions are imposed because of the single fault tolerant RV Motion Control System (MCS), which does not meet ISS Safety requirements for catastrophic hazards and dictates 16 degree Solar Array Rotary Joint position, which ensures that ISS and RV relative motion post separation, does lead to collision. The purpose of this paper is to describe a methodology and the analysis that was performed to determine relative motion trajectories of the ISS and separating RV for nominal and contingency cases. Analysis was performed in three phases that included ISS free drift prior to Visiting Vehicle separation, ISS and Visiting Vehicle relative motion analysis and clearance analysis. First, the ISS free drift analysis determined the worst case attitude and attitude rate excursions prior to RV separation based on a series of different configurations and mass properties. Next, the relative motion analysis calculated the separation trajectories while varying the initial conditions, such as docking mechanism performance, Visiting Vehicle MCS failure, departure port location, ISS attitude and attitude rates at the time of separation, etc. The analysis employed both orbital mechanics and rigid body rotation calculations while accounting for various atmospheric conditions and gravity gradient effects. The resulting relative motion trajectories were then used to determine the worst case separation envelopes during the clearance analysis. Analytical models were developed individually for each stage and the results were used to build initial conditions for the following stages. In addition to the analysis approach, this paper also discusses the analysis results, showing worst case relative motion envelopes, the recommendations for ISS appendage positioning and the suggested approach for future analyses.

Menkin, Evgeny↗

DOE OSTI2024

New Safety Feedback Control Design to Guarantee Adequate Frequency Performance in Microgrids

Safety analysis of power systems is concerned with the system's ability to maintain critical variables within specified limits following a disturbance. Frequency control adequacy has become increasingly important as the system inertia decreases due to the increase in renewable energy penetration. Various controllers for inverters have been proposed to improve the system frequency response and few are capable to ensure the safety of the response. In this article, a diesel-wind energy system is considered and modeled as a switching system between normal, faulted, and post-fault modes. Further, a safety feedback controller is designed as a supplementary signal for a wind turbine generator such that the speed of the diesel generator stays within a permissible range in the presence of a finite energy disturbance. Numerical results on the modified 33-bus microgrid system obtained of the proposed novel approach indicate that the suggested control configuration can guarantee adequate frequency response without excessive conservativeness.

barrier function↗

DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗

DOE OSTI2025

Morphology, Deformations, and Photocatalytic Activity of Thermally Treated Brookite Titanium Dioxide Thin Films

Metastable states, in which the coupling between long-range lattice deformations and electronic properties can be controlled, provide a pathway to tailoring the behavior of photocatalytic materials by directing the flow of photoinduced charge carriers. Brookite is a metastable polymorph of earth-abundant TiO 2 that exhibits photocatalytic function and, due to its high energy relative to the anatase and rutile polymorphs, may serve as a precursor for the formation of transitional metastable structures. In this work, facile thermal annealing is employed to promote the formation of predominantly brookite-phase films, regulate the brookite lattice distortions, and determine the effect of these distortions on charge separation, ultimately directed at enhancing photocatalytic activity. Profile fitting of X-ray diffraction patterns and peak shifts in Raman spectra revealed structural distortions of the brookite lattice. Structural defects, including lattice gliding, dislocations, stacking faults, and twin boundaries, were observed using scanning transmission electron microscopy. First-principles simulations reveal how the lattice distortions associated with stacking faults induce band bending, thus increasing the photocatalytic activity of brookite. In conclusion, this study provides insight into the microstructural tuning of metastable phases to enhance their unique functionalities.

band bending↗

DOE OSTI2021

Fingerprinting shock-induced deformations via diffraction

Abstract During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Helium bubble formation in nickel under in-situ krypton and helium ions dual-beam irradiation

Helium has profound effects on the microstructure and mechanical property of materials used in nuclear power systems. Here to understand the interaction between helium and irradiation-induced defects and its effect on bubble formation, single-beam and dual-beam irradiations on nickel were performed with 16 keV helium ions and 1 MeV krypton ions under in-situ TEM observation. For 1 MeV krypton single-beam irradiation, voids were observed at 600ºC but not at 500ºC and 700ºC. For helium single-beam irradiation at 500ºC, helium moved freely and formed bubble-loop complexes. For dual-beam irradiation at 500ºC, helium was trapped by vacancies and formed bubbles homogeneously. The structure of dislocation loops was also affected by the presence of helium. For krypton single-beam irradiation at 500ºC, the microstructure was {100} periodic wall of faulted and perfect dislocation loops. For dual-beam irradiation at 500ºC, periodic wall was not observed, and dislocation loops were predominantly perfect. For single-beam helium at 500ºC, dislocation loops were produced by athermal process of SIAs ejection or loop punching mechanisms. The interplay between helium and vacancies played an important role on the evolution of helium bubbles and defect microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Exceptional fracture toughness of CrCoNi-based medium- and high-entropy alloys at 20 kelvin

Medium- and high-entropy alloys based on the CrCoNi-system have been shown to display outstanding strength, tensile ductility and fracture toughness (damage-tolerance properties), especially at cryogenic temperatures. Here we examine the J Ic and (back-calculated) K JIc fracture toughness values of the face-centered cubic, equiatomic CrCoNi and CrMnFeCoNi alloys at 20 K. At flow stress values of ~1.5 GPa, crack-initiation K JIc toughnesses were found to be exceptionally high, respectively 235 and 415 MPa√m for CrMnFeCoNi and CrCoNi, with the latter displaying a crack-growth toughness K ss exceeding 540 MPa√m after 2.25 mm of stable cracking, which to our knowledge is the highest such value ever reported. Characterization of the crack-tip regions in CrCoNi by scanning electron and transmission electron microscopy reveal deformation structures at 20 K that are quite distinct from those at higher temperatures and involve heterogeneous nucleation, but restricted growth, of stacking faults and fine nanotwins, together with transformation to the hexagonal closed-packed phase. The coherent interfaces of these features can promote both the arrest and transmission of dislocations to generate respectively strength and ductility which strongly contributes to sustained strain hardening. Indeed, we believe that these nominally single-phase, concentrated solid-solution alloys develop their fracture resistance through a progressive synergy of deformation mechanisms, including dislocation glide, stacking-fault formation, nano-twinning and eventually in situ phase transformation, all of which serve to extend continuous strain hardening which simultaneously elevates strength and ductility (by delaying plastic instability), leading to truly exceptional resistance to fracture.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Observation of a mixed close-packed structure in superionic water

The study of superionic (SI) water has been a highly active research area since its theoretical prediction. Despite significant experimental and computational efforts, its melting curve and the stability of different oxygen lattices remain debated, impacting our understanding of SI ice’s peculiar transport properties. Experimental results at lower pressures show disagreement, whereas data at higher pressures are scarce due to the extreme challenges of such experiments. In this work, we present ultrafast X-ray diffraction results of water compressed by multiple shocks to pressures up to ~ 180 GPa. At pressures exceeding 150 GPa and temperatures around 2500 K, our diffraction patterns challenge the pure FCC-SI phase model, providing experimental evidence of the mixed close-packed superionic phase predicted by advanced ab initio calculations. At lower pressures, we observe simultaneous signatures of BCC and FCC structures within a pressure-temperature range consistent with some static-compression experiments, helping to resolve contradictory results in literature. These insights offer new constraints on the stability domains of SI phases and reveal detailed structural features, such as stacking faults. Our results advance the structural understanding of high-pressure SI ice to a level approaching that of ice I polymorphs, with potential implications for water-rich interiors of giant planets.

Andriambariarijaona, Leon [Centre National de la R↗

DOE OSTI2021

Optical vortex manipulation for topological quantum computation

Topological quantum computation based on Majorana bound states may enable new paths to fault-tolerant quantum computing. Several recent experiments have suggested that the vortex cores of topological superconductors, such as iron-based superconductors, may host Majorana bound states at zero energy. However, quantum computation with these zero-energy vortex bound states requires precise and fast manipulation of individual vortices, which is difficult to do in a scalable manner. To address this issue, in this study we propose a control scheme based on local heating via, for example, scanning optical microscopy to braid vortex-bound Majorana zero modes in a two-dimensional topological superconductor. First, we derive the conditions required for transporting a single vortex between two defects in the superconducting material by trapping it with a hot spot generated by local optical heating. Equipped with critical conditions for the vortex motion, we then establish the ideal material properties for vortex braiding and describe how transition errors resulting from finite speed and/or temperature can be minimized. Our work paves the way toward optical or microscopic control of zero-energy vortex bound states in two-dimensional topological superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

The complexity of { 11 2 ¯ 1 } twin interface structure and energetics in HCP materials

Here, the quest for a deeper understanding of the mechanical properties of HCP alloys has been an ongoing effort. Due to the lack of a sufficient number of slip systems, twin nucleation and migration originate as a substantial plastic deformation mechanism in HCP metals. The paper uncovers the energetics and atomistic mechanism of the {$11\overline{2}1$} < $11\overline{2}6$ > which is one of the most prominent twin modes in HCP materials that has not been fully understood. In view of corrugated {$11\overline{2}1$} interfaces, the positioning of lattice/motif atoms needs to be clarified. First, we establish the initial twin boundary (TB) structure unambiguously via the lattice offsets resulting in volume invariancy and correct relative positioning of the twin variants. Secondly, the TB migration mechanism is determined from a crystallographic (zero net shuffle, $Σ\vec{s}$ = 0) and energetic standpoints (minimum energy path). The Generalized Planar Fault Energy (GPFE) curves are calculated for Ti, Hf, Mg, and Zr. Such GPFE curves for this twin mode are calculated for the first time, and point to unusually high unstable energies (300-850 mJ/m 2 ), and low migration barriers (<60 mJ/m 2 ) for twinning accompanied with large shuffle magnitudes, |$\vec{s}$|, comparable to the Burger’s vector (|$\vec{s}$| ≅ $b_e$ > 0.44Å)

36 MATERIALS SCIENCE↗

DOE OSTI2023

Deformation mechanisms in an additively manufactured dual-phase eutectic high-entropy alloy

Nanostructured metals and alloys often exhibit high strengths but at the expense of reduced ductility. Through harnessing the far-from-equilibrium processing conditions of laser powder-bed fusion (L-PBF) additive manufacturing, we develop a dual-phase nanolamellar structure comprised of FCC/L1 2 and BCC/B2 phases in a Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 eutectic high-entropy alloy (EHEA), which exhibits a combination of ultrahigh yield strength (>1.4 GPa) and large tensile ductility (~17%). The deformation mechanisms of the additively manufactured EHEA are studied via in-situ synchrotron X-ray diffraction and high-resolution transmission electron microscopy. The high yield strength mainly results from effective blockage of dislocation motion by the high density of lamellar interfaces. The refined nanolamellar structures and low stacking fault energy (SFE) promote stacking fault (SF)-mediated deformation in FCC/L1 2 nanolamellae. The accumulation of abundant dislocations and SFs at lamellar interfaces can effectively elevate local stresses to promote dislocation multiplication and martensitic transformation in BCC/B2 nanolamellae. The cooperative deformation of the dual phases, assisted by the semi-coherent lamellar interfaces, gives rise to the large ductility of the as-printed EHEA. In addition, here we also demonstrate that post-printing heat treatment allows us to tune the non-equilibrium microstructures and deformation mechanisms. After annealing, the significantly reduced SFE and thicknesses of the FCC nanolamellae facilitate the formation of massive SFs. The dissolution of nano-precipitates in the BCC/B2 nanolamellae reduces spatial confinement and further promotes martensitic transformation to enhance work hardening. Our work provides fundamental insights into the rich variety of deformation mechanisms underlying the exceptional mechanical properties of the additively manufactured dual-phase nanolamellar EHEAs.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

NASA NTRS1992

What can we learn about impact mechanics from large craters on Venus?

More than 50 unequivocal peak-ring craters and multiringed impact basins have been identified on Venus from Earth-based Arecibo, Venera 15/16, and Magellan radar images. These ringed craters are relatively pristine, and so serve as an important new dataset that will further understanding of the structural and rheological properties of the venusian surface and of impact mechanics in general. They are also the most direct analogues for craters formed on the Earth in Phanerozoic time. Finite-element simulations of basin collapse and ring formation were undertaken in collaboration with V. J. Hillgren (University of Arizona). These calculations used an axisymmetric version of the viscoelastic finite element code TECTON, modeled structures on the scale of Klenova or Meitner, and demonstrated two major points. First, viscous flow and ring formation are possible on the timescale of crater collapse for the sizes of multiringed basins seen on Venus and heat flows appropriate to the plant. Second, an elastic lithosphere overlying a Newtonian viscous asthenosphere results mainly in uplift beneath the crater. Inward asthenospheric flow mainly occurs at deeper levels. Lithospheric response is dominantly vertical and flexural. Tensional stress maxima occur and ring formation by normal faulting is predicted in some cases, but these predicted rings occur too far out to explain observed ring spacings on Venus (or on the Moon). Overall, these estimates and models suggest that multiringed basin formation is indeed possible at the scales observed on Venus. Furthermore, due to the strong inverse dependence of solid-state viscosity on stress, the absence of Cordilleran-style ring faulting in craters smaller than Meitner or Klenova makes sense. The apparent increase in viscosity of shock-fluidized rock with crater diameter, greater interior temperatures accessed by larger, deeper craters, and decreased non-Newtonian viscosity associated with larger craters may conspire to make the transition with diameter from peak-ring crater to Orientale-type multiringed basin rather abrupt.

Mckinnon, William B.↗

NASA NTRS1993

Aspects of modelling the tectonics of large volcanoes on the terrestrial planets

Analytic solutions for the responses of planetary lithospheres to volcanic loads have been used to model faulting and infer elastic plate thicknesses. Predictions of the distribution of faulting around volcanic loads, based on the application of Anderson's criteria for faulting to the results of the models, do not agree well with observations. Such models do not give the stress state in the load itself, but only suggest a state of horizontal compressive stress there. Further, these models have considered only the effect of an instantaneously emplaced load. They do not address the time evolution of stresses, nor do they consider the effect of a load which grows. A finite element approach allows us to assign elements to the load itself, and thus permits calculation of the stress state and stress history within the edifice. The effects of episodic load growth can also be treated. When these effects are included, models give much better agreement with observations. We use the finite element code TECTON to construct axisymmetric models of volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have implemented time-dependent material properties in order to simulate incremental volcano growth. The viscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has no significant influence on the results. The code first calculates elastic deformations and stresses and then determines the time-dependent viscous deformations and stresses. Time in the model scales as the Maxwell time tau(m) in the asthenosphere. We consider a volcano 25 km in height and 200 km in radius on an elastic lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano consists of three load increments applied at intervals of 1000 tau(m). Contours of maximum deviatoric stress in the fully-grown edifice at the conclusion of flexure (t = 3000 tau(m)) are shown.

Mcgovern, Patrick J.↗