Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slip”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Improved ceramic slip casting technique

A primary concern in modern fluid dynamics research is the experimental verification of computational aerothermodynamic codes. This research requires high precision and detail in the test model employed. Ceramic materials are used for these models because of their low heat conductivity and their survivability at high temperatures. To fabricate such models, slip casting techniques were developed to provide net-form, precision casting capability for high-purity ceramic materials in aqueous solutions. In previous slip casting techniques, block, or flask molds made of plaster-of-paris were used to draw liquid from the slip material. Upon setting, parts were removed from the flask mold and cured in a kiln at high temperatures. Casting detail was usually limited with this technique -- detailed parts were frequently damaged upon separation from the flask mold, as the molded parts are extremely delicate in the uncured state, and the flask mold is inflexible. Ceramic surfaces were also marred by 'parting lines' caused by mold separation. This adversely affected the aerodynamic surface quality of the model as well. (Parting lines are invariably necessary on or near the leading edges of wings, nosetips, and fins for mold separation. These areas are also critical for flow boundary layer control.) Parting agents used in the casting process also affected surface quality. These agents eventually soaked into the mold, the model, or flaked off when releasing the case model. Different materials were tried, such as oils, paraffin, and even an algae. The algae released best, but some of it remained on the model and imparted an uneven texture and discoloration on the model surface when cured. According to the present invention, a wax pattern for a shell mold is provided, and an aqueous mixture of a calcium sulfate-bonded investment material is applied as a coating to the wax pattern. The coated wax pattern is then dried, followed by curing to vaporize the wax pattern and leave a shell mold of the calcium sulfate-bonded investment material. The shell mold is cooled to room temperature, and a ceramic slip is poured therein. After a ceramic shell of desired thickness has set up in the shell mold, excess ceramic slip is poured out. While still wet, the shell mold is peeled from the ceramic shell to expose any delicate or detailed parts, after which the ceramic shell is cured to provide a complete, detailed, precision ceramic article without parting lines.

Buck, Gregory M.↗
NASA NTRS2000

Distribution of Slip at the Northern Sumatran Fault System

We model spatial variations in horizontal displacements of 117 geodetic sites measured during annual surveys in 1989-1996 with the Global Positioning System (GPS) as elastic strain across a locked strike-slip fault to infer the contemporary slip rate, locking depth, and location of the Sumatran fault (SF) in northern Sumatra (1 S-3 N). GPS-derived slip rate estimates increase slightly northward from 23 plus or minus 3 mm/yr at 0.8 deg S to 26 plus or minus 2mm/yr at 2.7 N. They agree with geologic estimates north of the Equator, but at 0.5 S they are about 10 mm/yr higher. Strain appears to be distributed asymmetrically about the fault. South of 2 N, about 5 mm/yr of shear is required within the offshore forearc, west of the fault, to achieve a closer agreement of fault locations inferred from GPS velocities with geologically identified traces of the SF. Locking depth estimates are on the order of 10-20 km. The western branch of the major fault bifurcation near 1 N slips at a rate five times higher than the eastern branch. The two main strands of the fault at the northwestern tip of Sumatra (5.5 N) appear to be nearly free of horizontal strain; significant slip must occur away from the two strands, probably further east at two other geologically active branches. The Banda Aceh embayment is extruded to the northwest at a rate of 5 plus or minus 2 mm/yr. Within the estimated velocity uncertainties of several mm/yr, fault-normal deformation along the SF is insignificant. Almost strain free, the northern part of the back-arc basin is part of a rigid Sunda shelf, while the northern forearc is subjected to 8 plus or minus 5 x 10 (exp -8)/yr of extension nearly parallel to the arc.

Genrich, J. F.↗
NASA NTRS2002

Three-Dimensional Stress Fields and Slip Systems for Single Crystal Superalloy Notched Specimens

Single crystal superalloys have become increasingly popular for turbine blade and vane applications due to their high strength, and creep and fatigue resistance at elevated temperatures. The crystallographic orientation of a single crystal material greatly affects its material properties, including elastic modulus, shear modulus, and ductility. These directional properties, along with the type of loading and temperature, dictate an anisotropic response in the yield strength, creep resistance, creep rupture ductility, fatigue resistance, etc. A significant amount of research has been conducted to determine the material properties in the <001> orientation, yet the material properties deviating from the <001> orientation have not been assessed for all cases. Based on the desired application and design criteria, a crystal orientation is selected to yield the maximum properties. Currently, single crystal manufacturing is able to control the primary crystallographic orientation within 15 of the target orientation, which is an acceptable deviation to meet both performance and cost guidelines; the secondary orientation is rarely specified. A common experiment is the standard load-controlled tensile test, in which specimens with different orientations can be loaded to observe the material response. The deformation behavior of single-crystal materials under tension and compression is known to be a function of not only material orientation, but also of varying microdeformation (i.e. dislocation) mechanisms. The underlying dislocation motion causes deformation via slip, and affects the activation of specific slip systems based on load and orientation. The slip can be analyzed by observing the visible traces left on the surface of the specimen from the slip activity within the single crystal material. The goal of this thesis was to predict the slip systems activated in three-dimensional stress fields of a notched tensile specimen, as a function of crystal orientation, using finite element analysis without addressing microstructural deformation mechanisms that govern their activation. Out of three orientations tested, the specimen with a [110] load orientation and a [001] growth direction had the lowest maximum resolved shear stress; this specimen orientation appears to be the best design candidate for a tensile application.

Magnan, Shannon M.↗
NASA NTRS2002

[Evaporating Droplet and Imaging Slip Flows]

In this report, we summarize work on Evaporating Droplet and Imaging Slip Flows. The work was primarily performed by post-doc Hue Hu, and partially by grad students Lei Li and Danish Chopra. The work includes studies on droplet evaporation and its effects on temperature and velocity fields in an evaporating droplet, new 3-D microscopic particle image velocimetry and direct visualization on wall slip in a surfactant solution. With the exception of the slip measurements, these projects were those proposed in the grant application. Instead of slip flow, the original grant proposed imaging electro-osmotic flows. However, shortly after the grant was issued, the PI became aware of work on electro-osmotic flows by the group of Saville in Princeton that was similar to that proposed, and we therefore elected to carry out work on imaging slip flows rather than electro-osmotic flows.

Larson, R. G.↗
NASA NTRS1998

Minimization of Retinal Slip Cannot Explain Human Smooth-Pursuit Eye Movements

Existing models assume that pursuit attempts a direct minimization of retinal image motion or "slip" (e.g. Robinson et al., 1986; Krauzlis & Weisberger, 1989). Using occluded line-figure stimuli, we have previously shown that humans can accurately pursue stimuli for which perfect tracking does not zero retinal slip (Neurologic ARCO). These findings are inconsistent with the standard control strategy of matching eye motion to a target-motion signal reconstructed by adding retinal slip and eye motion, but consistent with a visual front-end which estimates target motion via a global spatio-temporal integration for pursuit and perception. Another possible explanation is that pursuit simply attempts to minimize slip perpendicular to the segments (and neglects parallel "sliding" motion). To resolve this, 4 observers (3 naive) were asked to pursue the center of 2 types of stimuli with identical velocity-space descriptions and matched motion energy. The line-figure "diamond" stimulus was viewed through 2 invisible 3 deg-wide vertical apertures (38 cd/m2 equal to background) such that only the sinusoidal motion of 4 oblique line segments (44 cd/m2 was visible. The "cross" was identical except that the segments exchanged positions. Two trajectories (8's and infinity's) with 4 possible initial directions were randomly interleaved (1.25 cycles, 2.5s period, Ax = Ay = 1.4 deg). In 91% of trials, the diamond appeared rigid. Correspondingly, pursuit was vigorous (mean Again: 0.74) with a V/H aspect ratio approx. 1 (mean: 0.9). Despite a valid rigid solution, the cross however appeared rigid in 8% of trials. Correspondingly, pursuit was weaker (mean Hgain: 0.38) with an incorrect aspect ratio (mean: 1.5). If pursuit were just minimizing perpendicular slip, performance would be the same in both conditions.

Stone, Leland S.↗
NASA NTRS2005

Path Following with Slip Compensation for a Mars Rover

A software system for autonomous operation of a Mars rover is composed of several key algorithms that enable the rover to accurately follow a designated path, compensate for slippage of its wheels on terrain, and reach intended goals. The techniques implemented by the algorithms are visual odometry, full vehicle kinematics, a Kalman filter, and path following with slip compensation. The visual-odometry algorithm tracks distinctive scene features in stereo imagery to estimate rover motion between successively acquired stereo image pairs, by use of a maximum-likelihood motion-estimation algorithm. The full-vehicle kinematics algorithm estimates motion, with a no-slip assumption, from measured wheel rates, steering angles, and angles of rockers and bogies in the rover suspension system. The Kalman filter merges data from an inertial measurement unit (IMU) and the visual-odometry algorithm. The merged estimate is then compared to the kinematic estimate to determine whether and how much slippage has occurred. The kinematic estimate is used to complement the Kalman-filter estimate if no statistically significant slippage has occurred. If slippage has occurred, then a slip vector is calculated by subtracting the current Kalman filter estimate from the kinematic estimate. This slip vector is then used, in conjunction with the inverse kinematics, to determine the wheel velocities and steering angles needed to compensate for slip and follow the desired path.

Helmick, Daniel↗
NASA NTRS2011

Strike-Slip Fault Patterns on Europa: Obliquity or Polar Wander?

Variations in diurnal tidal stress due to Europa's eccentric orbit have been considered as the driver of strike-slip motion along pre-existing faults, but obliquity and physical libration have not been taken into account. The first objective of this work is to examine the effects of obliquity on the predicted global pattern of fault slip directions based on a tidal-tectonic formation model. Our second objective is to test the hypothesis that incorporating obliquity can reconcile theory and observations without requiring polar wander, which was previously invoked to explain the mismatch found between the slip directions of 192 faults on Europa and the global pattern predicted using the eccentricity-only model. We compute predictions for individual, observed faults at their current latitude, longitude, and azimuth with four different tidal models: eccentricity only, eccentricity plus obliquity, eccentricity plus physical libration, and a combination of all three effects. We then determine whether longitude migration, presumably due to non-synchronous rotation, is indicated in observed faults by repeating the comparisons with and without obliquity, this time also allowing longitude translation. We find that a tidal model including an obliquity of 1.2?, along with longitude migration, can predict the slip directions of all observed features in the survey. However, all but four faults can be fit with only 1? of obliquity so the value we find may represent the maximum departure from a lower time-averaged obliquity value. Adding physical libration to the obliquity model improves the accuracy of predictions at the current locations of the faults, but fails to predict the slip directions of six faults and requires additional degrees of freedom. The obliquity model with longitude migration is therefore our preferred model. Although the polar wander interpretation cannot be ruled out from these results alone, the obliquity model accounts for all observations with a value consistent with theoretical expectations and cycloid modeling.

Rhoden, Alyssa Rose↗
NASA NTRS2004

Path following using visual odometry for Mars Rover in high-slip environments.

An architecture for autonomous operation of Mars rovers in high slip environments has been designed, implemented, and tested. This architecture is composed of several key technologies that enable the rover to accurately follow a designated path, compensate for slippage, and reach intended goals independent of the terrain over which it is traversing (within the mechanical constraints of the mobility system). These technologies include: visual odometry, full vehicle kinematics, a Kalman filter, and a slop compensation/path follower. Visual odometry tracks distinctive scene features in stereo imagery to estimate rover motion between successively acquired stereo image pairs using a maximum likelihood motion estimation algorithm. The full vehicle kinematics for a rocker-bogie suspension system estimates motion, with a no-slip assumption, by measuring wheel, rates, and rocker, bogie, and steering angles. The Kalman filter merges data from an Inertial Measurement Unit (IMU) and visual odometry. This merged estimate is then compared to the kinematic estimate to determining (tracking into account estimate uncertainties) if and how much slippage has occurred then a slip vector is calculated by differencing the current Kalman filter estimate from the kinematic estimate. This slip vector is then used, in conjunction with the inverse kinematics, to determine the necessary wheel velocities and steering angles to compensate for slip and follow the desired path.

Roumeliotis, Stergios I.↗
DOE OSTI2022

Loading direction dependence of asymmetric response of pyramidal slip in rolled AZ31 magnesium alloy

Textured magnesium alloys usually exhibit anisotropic mechanical behavior due to the asymmetric activation of different twinning and slipping modes. This work focuses on the pyramidal slip responses of rolled AZ31 magnesium alloy under two loading conditions, compressive and tensile loading along the normal direction. Under the condition where the compressive loading direction is closely parallel to the c-axis of the unit cell, tensile twinning and basal slips are prohibited, dislocations then active and tend to accumulate at grain boundaries and form dislocation walls. Meanwhile, these dislocations exhibit zigzag morphologies, which result from the cross-slip from {$10\bar{1}1$} first-order pyramidal plane to {$11\bar{2}2$} second-order pyramidal plane, then back to {$10\bar{1}1$} first-order pyramidal plane. Under the condition where tensile twins are prevalent, {$10\bar{1}1$} first-order and {$11\bar{2}2$} second-order pyramidal dislocations are favorable to be activated. Both types of dislocations behave climb-like dissociations onto the basal plane, forming zigzag dislocations. first-order pyramidal plane to second-order pyramidal plane, then back to first-order pyramidal plane. Under the condition where tensile twins are prevalent, first-order and second-order pyramidal dislocations are favorable to be activated. Both types of dislocations behave climb-like dissociations onto the basal plane, forming zigzag dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain boundary slip – twin transmission in titanium

Here, using a combination of mechanical testing, scanning electron microscopy, and a unified crystal plasticity framework for discrete intragranular shear localization, we investigate intense, localized slip bands on prismatic planes and {10$\bar{1}$2} $\langle$$\bar{1}$011$\rangle$ tensile twins, and their transmission across the grain boundaries (GBs) in commercially pure titanium. The analyses show that the orientation and curvature of the GB influence the local stress fields at the GBs, and consequently, the slip/twin transmission across the boundary. In addition to host grain properties, neighboring grain properties, such as active slip systems and instances of heterogeneity like slip bands and twins, heavily affect the deformation mechanisms in each grain. Finally, the applicability of geometric factors to predict the transmission in the experimentally observed co-located pairs is discussed. The local stress field calculated by the discrete slip and twin band model is shown to be capable of determining whether a transmission has occurred in an observed co-located pair, and also the direction of the transmission.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling lattice rotation fields from discrete crystallographic slip bands in superalloys

Here in this work, we investigate the relationship between an intense slip band (ISB) and the zone of large lattice rotations that forms ahead of the tip of the ISB. We develop a crystal plasticity finite element model of a discrete ISB lying within an oligocrystalline assembly and calculate the local crystalline stress and lattice rotation fields generated by the ISB. The calculations demonstrate that, first, a region of severe lattice rotations, commonly referred to as a microvolume, does not form without the ISB, and second, large amounts of accumulated slip in the ISB are required to enlarge the microvolume to sizes and rotation magnitudes observed experimentally. Ahead of the ISB tip, the quintessential plastic zone always forms, but the atypical microvolume forms when non-concentrated and spatially diffuse slip is activated by the ISB-induced stress field. This result suggests that the detrimental ISB/microvolume pair will likely appear in pairs of crystals in which transmission of the slip from the ISB is severely blocked by the grain boundary, a hypothesis that we verify with a few target cases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

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↗
DOE OSTI2022

Modeling Sequences of Earthquakes and Aseismic Slip (SEAS) in Elasto-Plastic Fault Zones With a Hybrid Finite Element Spectral Boundary Integral Scheme

We present a coupled finite element spectral boundary integral framework for modeling sequences of earthquakes and aseismic slip on a 2-D planar rate-and-state fault with off-fault visco-plastic response in the plane strain approximation. The model resolves both slow aseismic deformation and inertia effects during rapid slip. As an application, we perform two sets of simulations with different choices of cohesion to explore the co-evolution of fault slip, bulk plasticity and local stress fields. The first set implements a relatively large value of the cohesion parameter, which results in limiting inelastic strain accumulation to dynamic rupture phases. The second set implements a smaller cohesion, allowing for plastic strain to accumulate in both seismic and aseismic phases. For the first model, our results indicate that the extent and distribution of plastic strain depend on the angle of maximum compressive principal stress. At larger angles, inelastic strain accumulates on the extensional side of a dynamically propagating rupture. At smaller angles, the extent of plasticity is limited to the compressional side of the domain. At smaller cohesion values, off-fault plasticity may occur during aseismic slip, which alters the nucleation characteristics and earthquake sequence pattern. Furthermore, our results at lower cohesion values indicate that plastic strain accumulation may occur in both the extensional and compressional sides of the off-fault bulk even at higher angles of maximum compression. In conclusion, this produces damage patterns that deviate from the traditional off-fault fan-like distribution observed in dynamic rupture simulations and emphasizes the significance of long-term deformation in interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗
DOE OSTI2023

The Stability Transition From Stable to Unstable Frictional Slip With Finite Pore Pressure

Pore fluids are ubiquitous throughout the lithosphere and are commonly invoked as the cause of induced seismicity and slow earthquakes. We perform lab experiments to address these questions for drained fault conditions and low pore pressure. We shear simulated faults at effective normal stress (σ' n ) of 20 MPa and pore pressures Pp from 1 to 4 MPa. We document the full range of lab earthquake behaviors from slow slip to elasto-dynamic rupture and show that slow slip can be explained by the slip rate dependence of the critical rheologic stiffness without dilatancy hardening or other fluid effects. Our fault permeabilities ranges from 10 –18 to 10 –17 m 2 with an initial porosity of 0.1 and estimated fluid diffusion time ≈1 s. Slow slip and quasi-dynamic fault motion may arise from high P p at higher pressures but dilatancy strengthening is not a general requirement.

58 GEOSCIENCES↗
DOE OSTI2023

Tailoring planar slip to achieve pure metal-like ductility in body-centred-cubic multi-principal element alloys

Uniform tensile ductility (UTD) is crucial for the forming/machining capabilities of structural materials. Normally, planar-slip induced narrow deformation bands localize the plastic strains and hence hamper UTD, particularly in body-centred-cubic (bcc) multi-principal element high-entropy alloys (HEAs), which generally exhibit early necking (UTD < 5%). Here we demonstrate a strategy to tailor the planar-slip bands in a Ti-Zr-V-Nb-Al bcc HEA, achieving a 25% UTD together with nearly 50% elongation-to-failure (approaching a ductile elemental metal), while offering gigapascal yield strength. The HEA composition is designed not only to enhance the B2-like local chemical order (LCO), seeding sites to disperse planar slip, but also to generate excess lattice distortion upon deformation-induced LCO destruction, which promotes elastic strains and dislocation debris to cause dynamic hardening. Further, this encourages second-generation planar-slip bands to branch out from first-generation bands, effectively spreading the plastic flow to permeate the sample volume. Moreover, the profuse bands frequently intersect to sustain adequate work-hardening rate (WHR) to large strains. Our strategy showcases the tuning of plastic flow dynamics that turns an otherwise-undesirable deformation mode to our advantage, enabling an unusual synergy of yield strength and UTD for bcc HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stress-dependent activation entropy in thermally activated cross-slip of dislocations

Cross-slip of screw dislocations in crystalline solids is a stress-driven thermally activated process essential to many phenomena during plastic deformation, including dislocation pattern formation, strain hardening, and dynamic recovery. Molecular dynamics (MD) simulation has played an important role in determining the microscopic mechanisms of cross-slip. However, due to its limited timescale, MD can only predict cross-slip rates in high-stress or high-temperature conditions. The transition state theory can predict the cross-slip rate over a broad range of stress and temperature conditions, but its predictions have been found to be several orders of magnitude too low in comparison to MD results. This discrepancy can be expressed as an anomalously large activation entropy whose physical origin remains unclear. Here, we resolve this discrepancy by showing that the large activation entropy results from anharmonic effects, including thermal softening, thermal expansion, and soft vibrational modes of the dislocation. We expect these anharmonic effects to be significant in a wide range of stress-driven thermally activated processes in solids.

Science & Technology - Other Topics↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗