SEARCH · Search NASA

Results for “Slip localization”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Photon-Instanton Collider Implemented by a Superconducting Circuit

Instantons, spacetime-localized quantum field tunneling events, are ubiquitous in correlated condensed matter and high-energy systems. However, their direct observation through collisions with conventional particles has not been considered possible. We show how recent advances in circuit quantum electrodynamics, specifically, the realization of galvanic coupling of a transmon qubit to a high-impedance transmission line, allows the observation of inelastic collisions of single microwave photons with instantons (phase slips). We develop a formalism for calculating the photon-instanton cross section, which should be useful in other quantum field theoretical contexts. In particular, we show that the inelastic scattering probability can significantly exceed the effect of conventional Josephson quartic anharmonicity and reach order-unity values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗

NASA NTRS1991

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗

DOE OSTI2024

A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282

Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.

36 MATERIALS SCIENCE↗

NASA NTRS1982

A miniature electro-optical air flow sensor

Miniature sensors are needed for rapid and uncomplicated installation on light aircraft engaged in stability research programs. One particularly difficult sensor to miniaturize to the required degree has been a flow angle and velocity sensor for measuring the local flow ahead of a wing. However, by using an electrooptical technique it was possible to overcome the encountered difficulties and to design a sensor satisfying the requirements. The developed sensor for measuring angle-of-attack, yaw, and airspeed was shown to be suitable for rapid instrumentation of research aircraft because of its small size. The size reduction was accomplished by a design feature which eliminates the need for slip rings and wiring within the movable components of the sensor.

Kershner, D. D.↗

DOE OSTI2022

Frictional properties of Opalinus Clay: influence of humidity, normal stress and grain size on frictional stability

SUMMARY The Opalinus Clay (OPA) is a clay-rich formation considered as a potential host rock for radioactive waste repositories and as a caprock for carbon storage in Switzerland. Its very low permeability (10−19 to 10−21 m2) makes it a potential sealing horizon, however the presence of faults that may be activated during the lifetime of a repository project can compromise the long-term hydrological confinement, and lead to mechanical instability. Here, we have performed laboratory experiments to test the effect of relative humidity (RH), grain size (g.s.) and normal stress on rate-and-state frictional properties and stability of fault laboratory analogues corresponding to powders of OPA shaly facies. The sifted host rock powders at different grain size fractions (<63 μm and 63 < g.s. < 125 μm), at room (∼25 per cent) and 100 per cent humidity, were slid in double-direct shear configuration, under different normal stresses (5–70 MPa). We observe that peak friction, μpeak and steady-state friction, μss, depend on water vapour content and applied normal stress. Increasing relative humidity from ∼25 per cent RH (room humidity) to 100 per cent RH causes a decrease of frictional coefficient from 0.41 to 0.35. The analysis of velocity-steps in the light of rate-and-state friction framework shows that the stability parameter (a–b) is always positive (velocity-strengthening), and it increases with increasing sliding velocity and humidity. The dependence of (a–b) on slip rate is lost as normal stress increases, for each humidity condition. By monitoring the variations of the layer thickness during the velocity steps, we observe that dilation (Δh) is directly proportional to the sliding velocity, decreases with normal stress and is unaffected by humidity. Microstructural analysis shows that most of the deformation is accommodated within B-shear zones, and the increase of normal stress (σn) promotes the transition from strain localization and grain size reduction to distributed deformation on a well-developed phyllosilicate network. These results suggest that: (1) the progressive loss of velocity dependence of frictional stability parameter (a–b) at σn > 35 MPa is dictated by a transition from localized to distributed deformation and (2) water vapour content does not affect the deformation mechanisms and dilation, whereas it decreases steady-state friction (μss), and enhances fault stability.

58 GEOSCIENCES↗

DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

NASA NTRS1987

Investigation of hypersonic rarefied flow on a spherical nose of the AOTV

The Navier-Stokes (NS) equations were integrated numerically for investigating the flow characteristics on the forepart of the spherical nose of a space vehicle such as the AOTV or AFE by a modified Accelerated Successive Replacement (ASR) scheme under hypersonic rarefied conditions. Technical feasibility of the mathematical approach was demonstrated by computing the flowfield on a spherical nose under conditions that the AFE encounters at times t = 15 and 20 seconds after its reentry into the atmosphere. Local similar solutions for the merged layer flow along the stagnation line of the sphere were developed. These are correct to the same degree of accuracy as the NS equations. These solutions provided stagnation line boundary conditions for the domain of integration on the spherical noise. Also, a parametric study of the stagnation line solution was made with a view to understand the flow characteristics in tunnels with different ambient fluids. Analytical expressions for surface slip temperature, jump conditions, and concentration level in the presence of the real gas effects at the top of the Knudsen layer were derived and used to calculate the stagnation line flowfield with nonequilibrium dissociation and ionization. A number of graphics were drawn to illustrate the basic physics of the flowfields. The present analysis can be extended to include real gas effects and to bodies of arbitrary shapes. It can further provide boundary conditions for integrating the NS equations in the near wake region.

Jain, Amolak C.↗

NASA NTRS1974

Fault pattern at the northern end of the Death Valley - Furnace Creek fault zone, California and Nevada

The author has identified the following significant results. The pattern of faulting associated with the termination of the Death Valley-Furnace Creek Fault Zone in northern Fish Lake Valley, Nevada was studied in ERTS-1 MSS color composite imagery and color IR U-2 photography. Imagery analysis was supported by field reconnaissance and low altitude aerial photography. The northwest-trending right-lateral Death Valley-Furnace Creek Fault Zone changes northward to a complex pattern of discontinuous dip slip and strike slip faults. This fault pattern terminates to the north against an east-northeast trending zone herein called the Montgomery Fault Zone. No evidence for continuation of the Death Valley-Furnace Creek Fault Zone is recognized north of the Montgomery Fault Zone. Penecontemporaneous displacement in the Death Valley-Furnace Creek Fault Zone, the complex transitional zone, and the Montgomery Fault Zone suggests that the systems are genetically related. Mercury mineralization appears to have been localized along faults recognizable in ERTS-1 imagery within the transitional zone and the Montgomery Fault Zone.

Liggett, M. A.↗

NASA NTRS1988

Particle-laden weakly swirling free jets: Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns with a standard deviation of 15 microns, were used. A single loading ratio of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0.0 to 0.3. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were compared with predictions from three types of multiphase flow analysis: locally homogeneous flow (LHF); deterministic separated flow (DSF); and stochastic separated flow (SSF). For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, due to the neglect of particle inertia. The LHF model predictions of particle mass flux also showed poor agreement with measurements due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux, because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent dispersion of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗

NASA NTRS1988

Particle-laden weakly swirling free jets - Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns with a standard deviation of 15 microns, were used. A single loading ratio of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0.0 to 0.3. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were compared with predictions from three types of multiphase flow analysis: locally homogeneous flow (LHF); deterministic separated flow (DSF); and stochastic separated flow (SSF). For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, due to the neglect of particle inertia. The LHF model predictions of particle mass flux also showed poor agreement with measurements due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux, because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent disperison of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗

NASA NTRS1994

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗

DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A mean field homogenization model for the mechanical response of ceramic matrix composites

SiC/SiC composites offer exceptional mechanical stability at high temperatures and under irradiation. These ceramic matrix composites are therefore strong candidate materials for future nuclear energy applications. Their mechanical response, which exhibits pseudo-plasticity, is mediated by matrix cracking, fiber debonding, and fiber pull-out due to slip. Here, this study introduces a mechanistic model for the behavior of unidirectionally reinforced SiC/SiC composites. Specifically a mean field homogenization approach is proposed to account for all deformation and degradation modes during mechanical deformation. The homogenization scheme relies on a Mori Tanaka method that is extended to consider the effects of the coating’s elasto-plastic response on the development of micromechanical fields. Further, the model proposed introduces a method to effectively account for the role of localized damage (i.e., cracks) on mechanical fields within both the fiber and the matrix. Upon validating the model against experimental data, the roles of interface sliding, coating dimensions and intrinsic elastic response, as well as of microstructure (e.g. porosity, fiber volume fraction) are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗