Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMAL STRESSES”

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 379 records · Page 21

Influence of AC Loss on Stress and Strain of Superconducting Coils

The second-generation (2G) high-temperature superconducting (HTS) coated conductors (CC) are increasingly used in power systems recently, especially in large-capacity superconducting magnetic energy storage (SMES). HTSCC in superconducting energy storage coil is subjected to thermal stress which is caused by thermal contraction due to AC loss. The thermal stress will degrade the critical currents and may cause the displacement of the superconducting coil and insulating epoxy resin. In this work, the AC loss of the superconducting coil is firstly obtained by finite element method (FEM) software COMSOL with the method of iterative conductivity. Furthermore, the effect of AC loss on the superconducting coil is thus analyzed by solving the coupling problems which include magnetic field, mechanical behaviors, and heat transfer in solids. The result shows that the AC loss of the superconducting coil is concentrated in the middle of the coil and the maximum AC loss density is 5.96 × 10{sup 5} W/m{sup 3}; the maximum stress laid at the inner radius of the coil and the maximum stress is 0.01 MPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

Examining polymer‐protein biophysical interactions with small‐angle x‐ray scattering and quartz crystal microbalance with dissipation

Abstract Polymer‐protein hybrids can be deployed to improve protein solubility and stability in denaturing environments. While previous work used robotics and active machine learning to inform new designs, further biophysical information is required to ascertain structure–function behavior. Here, we show the value of tandem small‐angle x‐ray scattering (SAXS) and quartz crystal microbalance with dissipation (QCMD) experiments to reveal detailed polymer‐protein interactions with horseradish peroxidase (HRP) as a test case. Of particular interest was the process of polymer‐protein complex formation under thermal stress whereby SAXS monitors formation in solution while QCMD follows these dynamics at an interface. The radius of gyration ( R g ) of the protein as measured by SAXS does not change significantly in the presence of polymer under denaturing conditions, but thickness and dissipation changes were observed in QCMD data. SAXS data with and without thermal stress were utilized to create bead models of the potential complexes and denatured enzyme, and each model fit provided insight into the degree of interactions. Additionally, QCMD data demonstrated that HRP deforms by spreading upon surface adsorption at low concentration as shown by longer adsorption times and smaller frequency shifts. In contrast, thermally stressed and highly inactive HRP had faster adsorption kinetics. The combination of SAXS and QCMD serves as a framework for biophysical characterization of interactions between proteins and polymers which could be useful in designing polymer‐protein hybrids.

60 APPLIED LIFE SCIENCES↗

Thermal shock resistance of ceramic matrix composites

The experimental and analytical investigation of the thermal shock phenomena in ceramic matrix composites is detailed. The composite systems examined were oxide-based, consisting of an aluminosilicate matrix with either polycrystalline aluminosilicate or single crystal alumina fiber reinforcement. The program was divided into three technical tasks; baseline mechanical properties, thermal shock modeling, and thermal shock testing. The analytical investigation focused on the development of simple expressions for transient thermal stresses induced during thermal shock. The effect of various material parameters, including thermal conductivity, elastic modulus, and thermal expansion, were examined analytically for their effect on thermal shock performance. Using a simple maximum stress criteria for each constituent, it was observed that fiber fracture would occur only at the most extreme thermal shock conditions and that matrix fracture, splitting parallel to the reinforcing fiber, was to be expected for most practical cases. Thermal shock resistance for the two material systems was determined experimentally by subjecting plates to sudden changes in temperature on one surface while maintaining the opposite surface at a constant temperature. This temperature change was varied in severity (magnitude) and in number of shocks applied to a given sample. The results showed that for the most severe conditions examined that only surface matrix fracture was present with no observable fiber fracture. The impact of this damage on material performance was limited to the matrix dominated properties only. Specifically, compression strength was observed to decrease by as much as 50 percent from the measured baseline.

Carper, D. M.↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Role of mechanical loads in inducing in-cycle tensile stress in thermally grown oxide

Experimental in situ synchrotron x-ray diffraction results tracking the strain behavior of the various layers during a cycle, under thermo-mechanical conditions are presented in this work. The quantitative strain measurements here show that the thermally grown oxide briefly experiences in-plane tensile stress ({sigma}{sub 22} = +36.4 MPa) with increased mechanical loading during ramp-up in the thermal cycle. These findings are the first in situ experimental observations of these strains under thermo-mechanical conditions, envisaged to serve as a catalyst for crack initiation. The depth resolved measurements of strain taken during applied thermal and mechanical load in this work are a significant step towards achieving realistic testing conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microstructure actuation and gas sensing by the Knudsen thermal force

The generation of forces and moments on structures immersed in rarefied non-isothermal gas flows has received limited practical implementation since first being discovered over a century ago. The formation of significant thermal stresses requires both large thermal gradients and characteristic dimensions which are comparable to the gas molecular mean free path. For macroscopic geometries, this necessitates impractically high temperatures and very low pressures. At the microscale, however, these conditions are easily achieved, allowing the effects to be exploited, namely, for gas-property sensing and microstructure actuation. In this letter, we introduce and experimentally evaluate performance of a microelectromechanical in-plane Knudsen radiometric actuator, a self-contained device having Knudsen thermal force generation, sensing, and tuning mechanisms integrated onto the same platform. Sensitivity to ambient pressure, temperature gradient, as well as gas composition is demonstrated. Results are presented in terms of a non-dimensional force coefficient, allowing measurements to be directly compared to the previous experimental and computational data on out-of-plane cantilevered configurations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Triiodide Attacks the Organic Cation in Hybrid Lead Halide Perovskites: Mechanism and Suppression

Molecular I 2 can be produced from iodide-based lead perovskites under thermal stress; triiodide, I 3 − , is formed from this I 2 and I − . Triiodide attacks protic cation MA + - or FA + -based lead halide perovskites (MA + , methylammonium; FA + , formamidinium) as explicated through solution-based nuclear magnetic resonance (NMR) studies: triiodide has strong hydrogen-bonding affinity for MA + or FA + , which leads to their deprotonation and perovskite decomposition. Triiodide is a catalyst for this decomposition that can be obviated through perovskite surface treatment with thiol reducing agents. In contrast to methods using thiol incorporation into perovskite precursor solutions, no penetration of the thiol into the bulk perovskite is observed, yet its surface application stabilizes the perovskite against triiodide-mediated thermal stress. Thiol applied to the interface between FAPbI 3 and Spiro-OMeTAD (“Spiro”) prevents oxidized iodine species penetration into Spiro and thus preserves its hole-transport efficacy. Surface-applied thiol affects the perovskite work function; it ameliorates hole injection into the Spiro overlayer, thus improving device performance. It helps to increase interfacial adhesion (“wetting”): fewer voids are observed at the Spiro/perovskite interface if thiols are applied. Perovskite solar cells (PSCs) incorporating interfacial thiol treatment maintain over 80% of their initial power conversion efficiency (PCE) after 300 h of 85 °C thermal stress.

36 MATERIALS SCIENCE↗

Interfacial Shear Strength of Oxide Scale and SS 441 Substrate

Recent developments on decreasing the operating temperature for Solid Oxide Fuel Cells (SOFCs) have enabled the use of high temperature ferritic alloys as interconnect materials. Oxide scale will inevitably grow on the ferritic interconnects in a high temperature oxidation environment of SOFCs. The growth of the oxide scale induces growth stresses in the scale layer and on the scale/substrate interface. These growth stresses combined with the thermal stresses induced upon stacking cooling by the thermal expansion coefficient mismatch between the oxide scale and the substrate may lead to scale delamination/buckling and eventual spallation, which may lead to serious cell performance degradation. Hence the interfacial adhesion strength between the oxide scale and the substrate is crucial to the reliability and durability of the metallic interconnect in SOFC operating environments. In this paper, we applied an integrated experimental/modeling methodology to quantify the interfacial adhesion strength between the oxide scale and the SS 441 metallic interconnect. The predicted interfacial strength is discussed in details.

30 DIRECT ENERGY CONVERSION↗

Laser-zone growth in a Ribbon-To-Ribbon (RTR) process. Silicon sheet growth development for the large area sheet task of the low-cost solar array project

A new calculation of the effects of thermal stresses during growth on silicon ribbon quality is reported. Thermal stress distributions are computed for ribbon growth under a variety of temperature profiles. A growth rate of 55 cu cm/min with a single ribbon was achieved. The growth of RTR ribbon with a fairly uniform parallel dendritic structure was demonstrated. Results with two approaches were obtained for reducing the Mo impurity level in polycrystalline feedstock. Coating the Mo substrate with Si3N4 does not effect thermal shear separation of the polyribbon; this process shows promise of improving cell efficiency and also increasing the useful life of the molybdenum substrate. A number of solar cells were fabricated on RTR silicon grown from CVD feedstock.

Baghdadi, A.↗

Determination of elastic properties of polycrystalline U 3 Si 2 using resonant ultrasound spectroscopy

For this research, the elastic properties of U 3 Si 2 at room temperature have been measured via resonant ultrasound spectroscopy. Results show that the average value of Young's and the bulk modulus for U 3 Si 2 are 130.4±0.5 and 68.3±0.5 GPa, respectively. Further, a numerical model to assess thermal stress in an operating fuel is evaluated. Lastly, the thermal stress evolved in U 3 Si 2 is compared to UO 2 to facilitate an estimation of the probability of crack formation in U 3 Si 2 under representative light water reactor operating conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High temperature thermal energy storage in steel and sand

The technical and economic potential for high temperature (343 C, 650 F) thermal energy storage in hollow steel ingots, pipes embedded in concrete, and for pipes buried in sand was evaluated. Because it was determined that concrete would separate from pipes due to thermal stresses, concrete was replaced by sand, which is free from thermal stresses. Variations of the steel ingot concept were not cost effective compared to the sand-pipe approach, therefore, the sand-pipe thermal storage unit (TSU) was evaluated in depth to assess the approximate tube spacing requirements consistent with different system performance characteristics and also attendant system costs. For large TSUs which do not require fast response times, the sand-pipe approach offers attractive possibilities. A pipe diameter about 9 cm (3.5 in) and pipe spacing of approximately 25 cm (10 in), with sand filling the interspaces, appears appropriate. Such a TSU system designed for 8 hours charge/discharge cycle has an energy unit storage cost (CE) of $2.63/kWhr-t and a power unit storage cost (Cp) of $42/kW-t (in 1977 dollars).

Turner, R. H.↗

LNG Cascading Damage Study Volume I: Fracture Testing Report

As part of the liquefied natural gas (LNG) Cascading Damage Study, a series of structural tests were conducted to investigate the thermal induced fracture of steel plate structures. The thermal stresses were achieved by applying liquid nitrogen (LN 2 ) onto sections of each steel plate. In addition to inducing large thermal stresses, the lowering of the steel temperature simultaneously reduced the fracture toughness. Liquid nitrogen was used as a surrogate for LNG due to safety concerns and since the temperature of LN 2 is similar (-190°C) to LNG (-161°C). The use of LN 2 ensured that the tests could achieve cryogenic temperatures in the range an actual vessel would encounter during a LNG spill. There were four phases to this test series. Phase I was the initial exploratory stage, which was used to develop the testing process. In the Phase II series of tests, larger plates were used and tested until fracture. The plate sizes ranged from 4 ft square pieces to 6 ft square sections with thicknesses from ¼ inches to ¾ inches. This phase investigated the cooling rates on larger plates and the effect of different notch geometries (stress concentrations used to initiate brittle fracture). Phase II was divided into two sections, Phase II-A and Phase II-B. Phase II-A used standard A36 steel, while Phase II-B used marine grade steels. In Phase III, the test structures were significantly larger, in the range of 12 ft by 12 ft by 3 ft high. These structures were designed with more complex geometries to include features similar to those on LNG vessels. The final test phase, Phase IV, investigated differences in the heat transfer (cooling rates) between LNG and LN 2 . All of the tests conducted in this study are used in subsequent parts of the LNG Cascading Damage Study, specifically the computational analyses.

03 NATURAL GAS↗

Thermal drawdown and late‐stage seismic‐slip fault reactivation in enhanced geothermal reservoirs

Abstract Late‐stage seismic slip in geothermal reservoirs has been shown as a potential mechanism for inducing seismic events of magnitudes to ~2.6 as late as two decades into production. We investigate the propagation of fluid pressures and thermal stresses in a prototypical geothermal reservoir containing a centrally located critically stressed fault from a doublet injector and withdrawal well to define the likelihood, timing, and magnitude of events triggered by both fluid pressures and thermal stresses. We define two bounding modes of fluid production from the reservoir. For injection at a given temperature, these bounding modes relate to either low‐ or high‐relative flow rates. At low relative dimensionless flow rates the pressure pulse travels slowly, the pressure‐driven changes in effective stress are muted, but thermal drawdown propagates through the reservoir as a distinct front. This results in the lowest likelihood of pressure‐triggered events but the largest likelihood of late‐stage thermally triggered events. Conversely, at high relative non‐dimensional flow rates the propagating pressure pulse is larger and migrates more quickly through the reservoir but the thermal drawdown is uniform across the reservoir and without the presence of a distinct thermal front, and less capable of triggering late‐stage seismicity. We evaluate the uniformity of thermal drawdown as a function of a dimensionless flow rate Q D that scales with fracture spacing s ( m ), injection rate q ( kg/s ), and the distance between the injector and the target point L ∗ ( Q D ∝ qs 2 / L ∗ ). This parameter enables the reservoir characteristics to be connected with the thermal drawdown response around the fault and from that the corresponding magnitude and timing of seismicity to be determined. These results illustrate that the dimensionless temperature gradient adjacent to the fault dT D / dx D is exclusively controlled by the factor Q D . More significantly, this temperature gradient correlates directly with both the likelihood and severity of triggered events, enabling the direct scaling of likely magnitudes and timing to be determined a priori and directly related to the characteristics of the reservoir. This dimensionless scaling facilitates design for an optimum Q D value to yield both significant heat recovery and longevity of geothermal reservoirs while minimizing associated induced seismicity.

Gan, Quan↗

Characterization of SiC-SiC accident tolerant fuel cladding after stress application

In this work, SiC-SiC tubes with representative LWR (water cooled reactor) cladding dimensions were subjected to cycles of out-of-pile thermal stress followed by quenching. Samples were held at elevated temperatures and dropped into a water bath to replicate the high thermal stresses generated during reflooding. The samples were then characterized to measure changes in strength, permeability, and structure. Results were evaluated both for quenching into room temperature and boiling water. The results show that SiC-SiC composite tubes quenched from 1000 Celsius degrees in air into a boiling water bath exhibited no increase in He leak rate through the cylinder wall up to a pressure difference of 34.5 kPa compared to the as-fabricated state. SiC-SiC composite tubes quenched from 300-400 C. degrees into room temperature water showed an increase in He leak rate compared to the as-fabricated state. We attribute this to the more severe thermal gradient generated through the wall thickness due to higher heat flux characteristic of nucleate boiling compared to film boiling. A persistent vapor film forms around the samples quenched into boiling water, whereas no vapor film forms for the samples quenched into room temperature water. Another important result is that not loss in strength was noted compared to unquenched samples and that the stresses generated by the quenching were not enough to cause fiber fracture.

36 MATERIALS SCIENCE↗

A historical perspective of the YF-12A thermal loads and structures program

Around 1970, the Y-F-12A loads and structures efforts focused on numerous technological issues that needed defining with regard to aircraft that incorporate hot structures in the design. Laboratory structural heating test technology with infrared systems was largely created during this program. The program demonstrated the ability to duplicate the complex flight temperatures of an advanced supersonic airplane in a ground-based laboratory. The ability to heat and load an advanced operational aircraft in a laboratory at high temperatures and return it to flight status without adverse effects was demonstrated. The technology associated with measuring loads with strain gages on a hot structure was demonstrated with a thermal calibration concept. The results demonstrated that the thermal stresses were significant although the airplane was designed to reduce thermal stresses. Considerable modeling detail was required to predict the heat transfer and the corresponding structural characteristics. The overall YF-12A research effort was particularly productive, and a great deal of flight, laboratory, test and computational data were produced and cross-correlated.

Jenkins, Jerald M.↗

Expansion Compression Contacts for Thermoelectric Legs

In a proposed alternative to previous approaches to making hot-shoe contacts to the legs of thermoelectric devices, one relies on differential thermal expansion to increase contact pressures for the purpose of reducing the electrical resistances of contacts as temperatures increase. The proposed approach is particularly applicable to thermoelectric devices containing p-type (positive-charge-carrier) legs made of a Zintl compound (specifically, Yb14MnSb11) and n-type (negative charge-carrier) legs made of SiGe. This combination of thermoelectric materials has been selected for further development, primarily on the basis of projected thermoelectric performance. However, it is problematic to integrate, into a practical thermoelectric device, legs made of these materials along with a metal or semiconductor hot shoe that is required to be in thermal and electrical contact with the legs. This is partly because of the thermal-expansion mismatch of these materials: The coefficient of thermal expansion (CTE) of SiGe is 4.5 x 10(exp -6) C (exp -1), while the CTE of Yb14MnSb11 is 20 x 10(exp -6) C(exp -1). Simply joining a Yb14MnSb11 and a SiGe leg to a common hot shoe could be expected to result in significant thermal stresses in either or both legs during operation. Heretofore, such thermal stresses have been regarded as disadvantageous. In the proposed approach, stresses resulting from the CTE mismatch would be turned to advantage.

Sakamoto, Jeffrey↗

Determination of Creep Behavior of Thermal Barrier Coatings Under Laser Imposed High Thermal and Stress Gradient Conditions

A laser sintering/creep technique has been established to determine the creep behavior of thermal barrier coatings under steady-state high heat flux conditions. For a plasma sprayed zirconia-8 wt. % yttria coating, a significant primary creep strain and a low apparent creep activation energy were observed. Possible creep mechanisms involved include stress induced mechanical sliding and temperature and stress enhanced cation diffusion through the splat and grain boundaries. The elastic modulus evolution, stress response, and total accumulated creep strain variation across the ceramic coating are simulated using a finite difference approach. The modeled creep response is consistent with experimental observations.

Zhu, Dongming↗