Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic recovery”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Creep of chemically vapour deposited SiC fibres

The creep, thermal expansion, and elastic modulus properties for chemically vapor deposited SiC fibers were measured between 1000 and 1500 C. Creep strain was observed to increase logarithmically with time, monotonically with temperature, and linearly with tensile stress up to 600 MPa. The controlling activation energy was 480 + or - 20 kJ/mole. Thermal pretreatments near 1200 and 1450 C were found to significantly reduce fiber creep. These results coupled with creep recovery observations indicate that below 1400 C fiber creep is anelastic with negligible plastic component. This allowed a simple predictive method to be developed for describing fiber total deformation as a function of time, temperature, and stress. Mechanistic analysis of the property data suggests that fiber creep is the result of beta-SiC grain boundary sliding controlled by a small percent of free silicon in the grain boundaries.

Dicarlo, J. A.↗
DOE OSTI2021

Thermal Evolution of Internal Strain in Doped PbTe

Recent improvements in the efficiency of heat-to-electricity energy conversion in lead chalcogenide thermoelectrics involve reducing the thermal conductivity by incorporating large amounts of internal strain. The extent to which typical lead chalcogenide processing techniques (such as doping, ball milling, and densification) increase internal strain and dislocation density must be quantified to improve materials design. In this study, neutron powder diffraction is leveraged to evaluate the internal strain introduced by ball milling in doped and undoped powders. Doping with Na and/or Eu increases internal strain beyond ball milling alone, with the greatest increase from combining the two dopants. Strain recovery occurs in each powder above 400 K but can be suppressed by co-doping, indicating a strong dopant–dislocation interaction in this system. Therefore, high-temperature processing of PbTe powders should be avoided if high internal strain is desired. Low-temperature densification and/or rapid pressing techniques may be key to maintaining internal strain in the final pressed pellet. The diffraction peak asymmetry and correlated elastic softening measured in pressed PbTe pellets in past studies were not observed in the precursor powders measured here, suggesting that measurements of the Debye temperature on final pressed pellets are required to examine the influence of defect-induced elastic softening on thermal conductivity. This work provides key guidance for defect engineering to maximize internal strain and thermoelectric performance in PbTe thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical response and microstructure evolution from multiaxial stress relaxation in textured zircaloy nuclear cladding

Modern energy markets provide motivation for nuclear reactors to increase operational flexibility and load following capacity of existing reactor designs. The operational power changes result in fuel pellet expansion and contraction resulting in pellet-cladding mechanical interaction (PCMI) and high cladding stresses that are strongly related to stressed corrosion cracking (SCC) which can lead to fuel rod failure. This work aims to quantify the stress fields in cladding with a focus on multiaxial stress relaxation response from an imposed strain and the resulting microstructure changes in textured zircaloy claddings. Stress relaxation tests resulting in multiaxial stress fields are performed to measure the material response of Zircaloy-4 nuclear cladding. Multiple experiments using samples from the same piece of cladding are performed so samples with the same processing history may be removed at strategic points throughout mechanical testing for microstructure characterization. The microstructure evolution is analyzed with mechanical response to evaluate the static recovery impact on crystallographic texture and grain morphology. Experiments are performed at Idaho National Laboratory using an experimental device capable of applying internal pressure and independent axial load to commercial nuclear cladding at elevated temperatures with stress- or strain-controlled experimental capability. X-ray diffraction and electron backscatter diffraction are used for microstructure analysis. Rheological models are used to visualize the multiaxial stress-strain relationships and isolate elastic, plastic, and viscoplastic material properties for analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗
NASA NTRS1986

Estimation of high temperature low cycle fatigue on the basis of inelastic strain and strainrate

Fatigue life at elevated temperature can be predicted by introducing parametric values obtained from monotonic constitutive behavior into the Universal-Slopes Equation. For directionally solidified MAR-M200+HF at 975 C, these parameters are the maximum stress achievable under entirely plastic (time-independent) and purely creep (time-dependent) conditions and the corresponding inelastic strains, as well as the elastic modulus. For materials which exhibit plasticity/creep interaction, two more pairs of monotonic parameters must be evaluated for fatigue life prediction. This life-prediction method based on the Universal-Slopes Equation, resulted from a constitutive model characterizing monotonic and cyclic data as inelastic strainrate as a function of inelastic strain. Characterizing monotonic data is this way, permitted distinction between different material responses such as strain-hardening, strain-softening, and dynamic recovery effects. Understanding and defining the region of influence of each of these effects facilitated formulation of the constitutive model in relation to the mechanical and microstructural processes occurring in the material under cyclic loading.

Berkovits, A.↗
DOE OSTI2023

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗
DOE OSTI2021

Structural response of 3D-printed rubber lattice structures under compressive fatigue

In this work, a new direct UV curing, polyurethane-like, highly stretchable rubber, Elastic ToughRubber 90™ (ETR 90), is demonstrated for 3D printing via photocuring with Texas Instruments DLP® optical systems. An array of test specimens, comprising a sandwich structure with a body-centered cubic lattice core, are 3D printed from ETR 90. The compressive behavior of these specimens is characterized up to 10 k and up to 100 k fatigue-loading cycles. The internal, visually obstructed, structures under loads are visualized using in situ X-ray tomography. All lattice structures remain intact and full recovery is observed, indicating the cyclic durability of the 3D-printed rubber.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Robust multivariable controller design for flexible spacecraft

Large, flexible spacecraft are typically characterized by a large number of significant elastic modes with very small inherent damping, low, closely spaced natural frequencies, and the lack of accurate knowledge of the structural parameters. Summarized here is some recent research on the design of robust controllers for such spacecraft, which will maintain stability, and possible performance, despite these problems. Two types of controllers are considered, the first being the linear-quadratic-Gaussian-(LQG)-type. The second type utilizes output feedback using collocated sensors and actuators. The problem of designing robust LQG-type controllers using the frequency domain loop transfer recovery (LTR) method is considered, and the method is applied to a large antenna model. Analytical results regarding the regions of stability for LQG-type controllers in the presence of actuator nonlinearities are also presented. The results obtained for the large antenna indicate that the LQG/LTR method is a promising approach for control system design for flexible spacecraft. For the second type of controllers (collocated controllers), it is proved that the stability is maintained in the presence of certain commonly encountered nonlinearities and first-order actuator dynamics. These results indicate that collocated controllers are good candidates for robust control in situations where model errors are large.

Joshi, Suresh M.↗
DOE OSTI2021

Dynamic flow stress of pure polycrystalline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based modeling to large strains

The dynamic thermo-mechanical behavior of pure aluminum has attracted renewed interest lately due to experimental observations of an anomalous increase in Hugoniot Elastic Limit (HEL) at incipient plasticity and elevated temperatures in polycrystalline pure metals. In context of current dislocation-mediated plasticity models for metals, this increase in dynamic strength is indicative of a transition in the rate controlling mechanism for dislocation glide from being thermally assisted to phonon-drag restricted due to increase in phonon viscosity at elevated temperatures. Though these studies have helped to shed light on these important mechanisms operative in FCC metals at incipient plasticity, the extent to which they contribute to flow stress, particularly at larger strains, remains unclear. Here, we address these questions through a combined experimental and modeling effort focused on investigating the evolution of dynamic flow stress in polycrystalline aluminum using experimental data gathered from a series of combined pressure-and-shear plate impact (PSPI) experiments designed to reveal the flow stress of pure aluminum at strain rates ~ 10 5 /s, plastic strains of up to 40% and temperatures ranging from room to 866 K. In all cases, the flow stress of aluminum, as inferred from the measured transverse particle velocity histories at the free surface of an elastic tungsten carbide target plate, reveals saturation with increasing plastic strains at stress levels that decrease with increasing test temperatures. Numerical simulations are performed to correlate the experimentally observed temperature and strain rate dependence of flow stress at small and large plastic strains using the Austin-McDowell dislocation-mediated plasticity model parametrized to normal plate impact experiments conducted in an earlier study by Zaretsky and Kanel (2012). Extensions to the model are made to better represent dynamic behavior of pure aluminum at larger plastic strains as observed in elevated temperature split Hopkinson Pressure bar (SPHB) experiments of Samanta (1971) and Lindholm and Yeakly (1965), and the combined pressure-and-shear plate-impact experiments conducted in the present study. The main theoretical extension to the Austin-McDowell model made in this paper is the introduction of a new rate- and temperature- dependent dynamic recovery function, which can potentially allow for an effective reduction in the rate of accumulation of dislocations at large plastic strains. The numerical predictions of the revised and re-calibrated plasticity model are brought into agreement with the experimental observations, correlating sufficiently well with dynamic yield stress at incipient plasticity and the flow stress levels at larger plastic strains, plastic strain rates in the range 10 3 – 10 6 /s, and elevated temperatures up to near melt. The model, in concert with the experimental measurements, suggests that at incipient plasticity the rate governing mechanism for plastic flow is phonon-drag restricted dislocation glide, whereas, at higher magnitudes of plastic strain, it transitions to stress-assisted thermally activated glide. The transition between these two rate governing mechanisms is controlled by the evolution of dislocations throughout the deformation process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Predictive analytics to direct clinical attention to complex patients with elevated suicide risk: enhancement of the Veterans Health Administration REACH VET model

Suicide is a major public health concern, particularly among Veterans. The U.S. Department of Veterans Affairs Veterans Health Administration (VHA) employs the Recovery Engagement and Coordination for Health–Veterans Enhanced Treatment (REACH VET) model to prioritise high-risk patients for targeted clinical attention. REACH VET 1.0 (RV 1.0) was developed on 2008–2011 data. To reflect changes in clinical practice and populations, VHA updated it to REACH VET 2.0 (RV 2.0). This study describes its development and validation. RV 2.0 used longitudinal data from 7,248,170 VHA patients (4,967 suicide deaths) in 2018–2019, with 650 time-varying demographic, clinical and area-level predictors derived from a 2-year lookback (2016–2019). An ensemble of Elastic-Net logistic regression models was trained on 2018 data and evaluated monthly at the population level in 2019, focusing on the top 0.1% intervention risk tier. Analyses assessed model discrimination, suicide detection, risk concentration, subgroup consistency (sex, age and race/ethnicity) and performance relative to RV 1.0 using the same percentile-based risk strata. RV 2.0 outperformed RV 1.0 across all risk strata, with better discrimination (C-statistic 0.76 vs 0.69) and consistent performance across demographic subgroups. Within the top 0.1% of predicted risk, RV 2.0 identified more deaths, higher suicide rates and greater mortality risk concentration both when averaged across the 12 monthly 2019 test sets (5.6 vs 3.6; 83.6 vs 53.7 per 100,000 person-years; 21.0 vs 14.1) and when annualised for 2019 (67 vs 43; 2.7% vs 1.7%; 1,003 vs 644 per 100,000 person-years; 26.7 vs 17.1). RV 2.0 improves suicide risk stratification among Veterans, demonstrating better performance and consistent prediction across subgroups and highlighting the need for regular model updates and evaluation.

Peluso, Alina [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2022

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗
DOE OSTI2022

Effect of loading path on grain misorientation and geometrically necessary dislocation density in polycrystalline aluminum under reciprocating shear

Solid phase processing (SPP) is a promising alloy fabrication technique to produce fine and homogeneous grain structures for high-performance alloys. However, there is very limited modeling capability to understand and predict the grain refinement during SPP. In this work, the crystal plasticity theory was used to study elastic-plastic deformation in polycrystalline aluminums under large shear deformation. Two approaches, kernel averaged misorientation (KAM) and grain reference orientation deviation (GROD), were used to assess the grain misorientations. The geometrically necessary dislocation (GND) density was computed with the plastic strain rate. The deformation simulations were carried out under two loading conditions to investigate the effect of loading paths on the evolutions of grain misorientation and GND density. The results show that the regions with high misorientation and GND density first appear near grain boundaries. These regions then extend toward interior grains. The loading path affects dislocation system activation and dislocation recovery, hence dislocation evolution and misorientation. In conclusion, both two- and three-dimensional simulations showed that the spatial and temporal evolutions of GROD, KAM, and GND density in are closely correlated, which indicates they all can be used as criteria of grain refinement or recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale Assessment of Caprock Integrity for Geologic Carbon Storage in the Pennsylvanian Farnsworth Unit, Texas, USA

Leakage pathways through caprock lithologies for underground storage of CO2 and/or enhanced oil recovery (EOR) include intrusion into nano-pore mudstones, flow within fractures and faults, and larger-scale sedimentary heterogeneity (e.g., stacked channel deposits). To assess multiscale sealing integrity of the caprock system that overlies the Morrow B sandstone reservoir, Farnsworth Unit (FWU), Texas, USA, we combine pore-to-core observations, laboratory testing, well logging results, and noble gas analysis. A cluster analysis combining gamma ray, compressional slowness, and other logs was combined with caliper responses and triaxial rock mechanics testing to define eleven lithologic classes across the upper Morrow shale and Thirteen Finger limestone caprock units, with estimations of dynamic elastic moduli and fracture breakdown pressures (minimum horizontal stress gradients) for each class. Mercury porosimetry determinations of CO2 column heights in sealing formations yield values exceeding reservoir height. Noble gas profiles provide a “geologic time-integrated” assessment of fluid flow across the reservoir-caprock system, with Morrow B reservoir measurements consistent with decades-long EOR water-flooding, and upper Morrow shale and lower Thirteen Finger limestone values being consistent with long-term geohydrologic isolation. Together, these data suggest an excellent sealing capacity for the FWU and provide limits for injection pressure increases accompanying carbon storage activities.

58 GEOSCIENCES↗
NASA NTRS2015

Cloud-Based Numerical Weather Prediction for Near Real-Time Forecasting and Disaster Response

The use of cloud computing resources continues to grow within the public and private sector components of the weather enterprise as users become more familiar with cloud‐computing concepts, and competition among service providers continues to reduce costs and other barriers to entry. Cloud resources can also provide capabilities similar to high‐performance computing environments, supporting multi‐node systems required for near real‐time, regional weather predictions. Referred to as "Infrastructure as a Service", or IaaS, the use of cloud-based computing hardware in an on‐demand payment system allows for rapid deployment of a modeling system in environments lacking access to a large, supercomputing infrastructure. Use of IaaS capabilities to support regional weather prediction may be of particular interest to developing countries that have not yet established large supercomputing resources, but would otherwise benefit from a regional weather forecasting capability. Recently, collaborators from NASA Marshall Space Flight Center and Ames Research Center have developed a scripted, on‐demand capability for launching the NOAA/NWS Science and Training Resource Center (STRC) Environmental Modeling System (EMS), which includes pre‐compiled binaries of the latest version of the Weather Research and Forecasting (WRF) model. The WRF‐EMS provides scripting for downloading appropriate initial and boundary conditions from global models, along with higher‐resolution vegetation, land surface, and sea surface temperature data sets provided by the NASA Short‐term Prediction Research and Transition (SPoRT) Center. This presentation will provide an overview of the modeling system capabilities and benchmarks performed on the Amazon Elastic Compute Cloud (EC2) environment. In addition, the presentation will discuss future opportunities to deploy the system in support of weather prediction in developing countries supported by NASA's SERVIR Project, which provides capacity building activities in environmental monitoring and prediction across a growing number of regional hubs throughout the world. Capacity‐building applications that extend numerical weather prediction to developing countries are intended to provide near real‐time applications to benefit public health, safety, and economic interests, but may have a greater impact during disaster events by providing a source for local predictions of weather‐related hazards, or impacts that local weather events may have during the recovery phase.

Molthan, Andrew↗
NASA NTRS1984

Shock experiments on maskelynite-bearing anorthosite

A series of shock recovery experiments over 9.9-60.4 GPa have been carried out on naturally shocked anorthosite from the Mistastin impact structure in Labrador consisting primarily of diaplectic plagioclase glass or maskelynite, An(50), and pyroxene. Petrographic observations of the experimental products indicate that the component minerals and diaplectic glasses generally retained their initial character throughout, the only exception being the increase in fracturing which occurred in the 9.9 GPa shot. Reshocking at pressures higher than the initial shock tends to lower the refractive index of maskelynite. The increase in refractive index of maskelynite reshocked to pressures lower than the initial pressure is interpreted as due to shock densification of the diaplectic glass above the Hugoniot elastic limit and below the mixed phase regime. The data suggest that the low-high-low density transition of maskelynite occurs about 8 GPa below that of the crystal of corresponding composition.

Lambert, P.↗
DOE OSTI2021

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗
DOE OSTI2022

Using post-processing heat treatments to elucidate precipitate strengthening of additively manufactured superalloy 718

The poor machinability and extensive work hardening of Ni-based superalloys makes additive manufacturing an attractive option for producing geometrically complex components with distinct microstructures. Although previous studies show recovery of high strength at room temperature, very few studies demonstrate successful properties at elevated temperatures required for industrial applications. Here, the objectives of this study are to present a post-build heat treatment for high strength across a wide temperature range, determine the strength contribution of nanoscale precipitating phases to the overall mechanical properties of superalloy 718, and from these, provide a comprehensive microstructure-property relationship for wrought and AM 718 to guide efforts to simulate the properties of AM components. Laser powder bed fusion–produced superalloy 718 was characterized at multiple length scales using scanning electron microscopy and transmission electron microscopy in the as-built condition and with multiple heat treatments designed to form combinations of γ', γ", and δ precipitates. Uniaxial tensile tests performed from room temperature to 600 °C on subsize specimens determined the yield strength, elastic modulus, ultimate tensile strength, fracture stress, and uniform elongation. Precipitates in this work proved to be weak barriers to dislocation motion through a dispersed barrier model, but they provided strength to the alloy through their consistent high density. The relative contribution to the yield strength from γ" remained consistent between 48% and 57% of the total strength up to 600 °C, the primary influence on the high temperature strength of superalloy 718. The strength factors for γ" and δ precipitates were found to trend inversely with tensile test temperature and may be attributable to the differences in precipitate coherency. A post-build heat treatment is recommended to maintain high strength at elevated temperatures. A quantitative microstructure-property relationship, dependent on precipitate size, density, and morphology, was derived and can estimate the yield strength across a wide temperature range applicable to the operational regimes for superalloy 718.

36 MATERIALS SCIENCE↗
NASA NTRS2007

Models of Active Glacial Isostasy Roofing Warm Subduction: Case of the South Patagonian Ice Field

Modern geodetic techniques such as precise Global Positioning System (GPS) and high-resolution space gravity mapping (Gravity Recovery and Climate Experiment, GRACE) make it possible to measure the present-day rate of viscoelastic gravitational Earth response to present and past glacier mass changes. The Andes of Patagonia contain glacial environments of dramatic mass change. These mass load changes occur near a tectonically active boundary between the Antarctic and South American plates. The mechanical strength of the continental side of this boundary is influenced by Neogene ridge subduction and by the subduction of a youthful oceanic slab. A ridge of young volcanos parallels the Pacific coastline. Release of volatiles (such as water) at depth along this ridge creates a unique rheological environment. To assess the influence of this rheological ridge structure on the observational land uplift rate, we apply a two dimensional viscoelastic Earth model. A numerical study is presented which examines the sensitivity of the glacial loading-unloading response to the complex structure at depth related to the subducting slab, the viscous wedge between slab and continental lithosphere, and the increase of elastic thickness from oceanic to continental lithosphere. A key feature revealed by our numerical experiments is a continuum flow wherein the slab subdues the material transport toward oceanic mantle and crust. The restricted flow is sensitive to the details of slab mechanical strength and penetration into the upper mantle. The reduced viscosity within the mantle wedge, however, enhances the load-induced material transport everywhere within the asthenosphere.

crustal mass changes↗