Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPRESSION TESTING”

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 217 records · Page 12

FY2025 Status Report: Model 9975 O Ring Fixture Long-Term Leak Performance

Leak testing experiments to monitor the aging performance of Viton® GLT and GLT-S O-rings used in the model 9975 shipping package have been ongoing since 2004 at Savannah River National Laboratory. Seventy tests using mock-up 9975 primary containment vessels (PCVs) with GLT O-rings were assembled and heated to temperatures ranging from 200 to 450 °F. Due to material substitution, fourteen tests with GLT-S O-rings were initiated in 2008 and heated to temperatures ranging from 200 to 400 °F. The conditioning temperatures are elevated compared to the calculated maximum O-ring temperature in a 9975 package in storage, 158 °F, to accelerate aging and enable observations of O-ring failures in a reasonable time frame. The mock-up PCV fixtures are leak tested periodically, and all GLT O-ring fixtures aged at 350 °F or above have failed to maintain a leak-tight seal. Eight GLT O-ring fixtures aged at 300 °F have failed after 2.8 to 5.7 years at temperature while the remaining fixtures at 300 °F were retired from testing following more than five years of aging without failure. Two of these retired fixtures were returned to testing and heated to 350 °F to evaluate the impact of additional heating at higher temperature for aged O-rings. Fixture #20 failed after 3 months while fixture #18 failed after 9 months at 350 °F. These O-rings demonstrated that aged and in-service O-rings can continue to be used, even at higher temperatures, after being in storage, and their leak performance are consistent with other samples at 350 °F. There has been one GLT O-ring fixture which failed after 13.4 years of aging at 200 °F. However, 20 other GLT O-rings aging at 200 °F have remained leak-tight for over 16.9 years and remain in test. There are two GLT O-ring fixtures at 270 °F; one fixture has failed after 12.9 years while the other fixture remains in test after 12.5 years. All GLT-S O-ring fixtures aged at 300 °F or above have failed their leak test. No failures have yet been observed in GLT-S O-ring fixtures aging at 250 °F for 14.9 years, while one GLT-S O-ring fixture failed after 12.4 years at 200 °F. The leak testing data to date suggest the GLT and GLT-S O-rings aging in the K-Area Complex (KAC) storage at temperatures of 158 °F might maintain a leak-tight seal for up to 59 years. Data from the O-ring fixtures are generally consistent with results from compression stress-relaxation testing and provide confidence in the predictive models based on those results. However, uncertainty exists in extrapolating these elevated temperature results to the lower temperatures of interest for normal storage in KAC. The collective data from these test efforts suggest the minimum O-ring service life at KAC normal storage conditions should be at least 34 years for GLT and GLT-S O-rings. Measurement of compression set in O-rings removed from failed fixtures, compared to that from KAC surveillance O-rings, indicate significant margin remains for O-rings still in service in 9975 packages in KAC. Aging and periodic leak testing will continue for the remaining 24 mock-up PCV fixtures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Low thermal inertia thermal barrier coatings for spark ignition engines: An experimental study

Application of thermal barrier coatings in spark ignition engines have historically been avoided due to the knock penalty associated with higher surface temperatures induced by the ceramic layer. However, advances in low thermal inertia coatings (i.e. temperature swing coatings) that combine low thermal conductivity with low volumetric heat capacity can prevent excessively high surface temperatures during the intake stroke and reduce or avoid knock while improving performance and efficiency. In this study, a novel coating material was tested in a low compression ratio spark ignition engine to evaluate the performance of an advanced low thermal inertia coating during steady and cold-start conditions. A total of four pistons with this novel material was tested, with an additional piston coated with gadolinium zirconate. Spark timing sweeps demonstrated a maximum 0.8% relative thermal efficiency gain with a thin coating of the novel material. Here, the novel material coated above 200-microns showed a deterioration in performance and efficiency due to charge heating increasing knock propensity. Cold-start tests demonstrated that charge heating is beneficial for reducing unburned hydrocarbons and particle matter emissions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Low-carbon fuels for spark-ignited engines: A comparative study of compressed natural gas and liquefied petroleum gas on a CFR engine with exhaust gas recirculation

Decades of work on low-carbon fuels have established their potential for substantial emissions reductions; however, their adoption is still limited by infrastructure concerns and engine efficiency deficits. As infrastructures have begun to evolve, research on strategies that maximize engine efficiency through interactions with fuel properties must also now take center place. This paper compares the performance, emissions, and combustion characteristics of two forefront low-carbon fuels: compressed natural gas (CNG) and liquefied petroleum gas (LPG) in a cooperative fuel research (CFR) engine over a range of compression ratios and engine loads. The effects of exhaust gas recirculation (EGR), end-gas auto-ignition, and a novel combustion control tool, the combustion intensity metric (CIM), were also evaluated at different stoichiometric engine operating conditions. In comparison to LPG, CNG operation demonstrated an extended knock-free regime, allowing engine operation at higher engine loads and compression ratios, but LPG operation exhibited enhanced combustion characteristics with higher peak pressures and faster apparent heat release rates (AHRR). LPG operation achieved higher brake thermal efficiencies and lower equivalent CO 2 emissions compared to CNG operation at the tested engine loads and compression ratios. LPG demonstrated significantly higher EGR tolerance limits compared to CNG, with a maximum of 28% EGR rate, compared to 23% for CNG. This improved EGR dilution tolerance was responsible for a 90% reduction in NOx emissions for LPG compared to a maximum of 70% with CNG. EGR dilution also exhibited more effective knock mitigation potential with LPG, suppressing knock intensity values by up to 98% and transitioning the engine operation towards normal combustion from heavy knocking conditions. As a result, the CIM was found to decrease burn durations and improve the quality of combustion by controlling the desired fraction of end-gas auto-ignition.

33 ADVANCED PROPULSION SYSTEMS↗

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗

A Novel Rotating Flexure-Test Technique for Brittle Materials with Circular Geometries

The four-point bend test is one of the simplest and often the preferred flexural strength evaluation method for brittle materials. In this loading mode, fracture often initiates from a critical surface (or subsurface) flaw when subjected to a tensile stress state. However, if the critical flaw exists on the compression side of the test specimen, it may not activate to grow a crack and hence the resulting flexural strength will be higher than the true value. The goal of this study is to measure the true flexural strength of a solid or hollow cylindrical brittle specimen by ensuring that failure occurs at its weakest point by rotating along its longitudinal axis, thereby exposing and activating its critical surface flaw during a four-point bend test. A novel test fixture has been designed and fabricated, and the true flexural strengths of cylindrical brittle and quasi-brittle tubular specimens have been measured and compared to existing experimental data obtained through traditional four-point bend tests. In this work, experimental results showcase the orientation dependance on flexural strength for various materials. Additionally, similarities between experimental findings and those available in literature, including observations of fracture surfaces and relationships between surface roughness and material strength, are discussed.

36 MATERIALS SCIENCE↗

Effect of Transverse Compression on the Dielectric Breakdown of Power Cable with Ethylene Propylene Rubber Insulation

Understanding of the failure mechanisms of submarine dynamic power cables (SDPC) is critical for innovative design to meet the 2035 cost reduction target of U.S. DOE Floating Offshore Wind Shot. This is important because the current design suffers a significant failure rate in the field. This project explored a multiphysics approach and conducted systematic electro-thermo-mechanical (ETM) experimental study on the power cores extracted from a 15 kV power cable. The power cable had three cores of copper conductor and ethylene propylene rubber (EPR) insulation. An ETM testing system with transverse compression capability was developed. Mechanical loading induced dielectric breakdown. Increasing temperature from room temperature (RT, 22 degrees Celsius) to 90 degrees Celsius resulted in the 67% decrease in the failure mechanical load as defined by the dielectric breakdown. With the same temperature rise, the creep mechanical load was reduced by 70% at a given dielectric breakdown time. The post-test measurement revealed an impressive recovery of electrical resistance, although the dielectric breakdown occurred in the ETM test that would not be captured without the use of ETM testing.

Wang, Hong [ORNL] (ORCID:0000000201730545)↗

Flowback Test Analyses at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) Site

Abstract In 2017 and 2019, injection testing was carried out in three zones in a vertical well in granite at the Frontier Observatory for Research in Geothermal Energy site near Milford, Utah, USA. In several injection cycles, flowback was implemented rather than shut-in. The goal was to explore an alternative to prolonged shut-in periods for inferring closure stress, formation compressibility, and formation permeability (permeability thickness product). The flowback procedures involved a cyclic flowback/shut-in, while pressure decreased. The flowback data are presented, and analyses are shown. The inferred closure stress(es) from flowback analyses are lower than for equivalent injection cycles that were strictly shut-in. Relatively high formation compressibility obtained from the flowback analyses indicates an extensive, fractured system. This study also includes numerical simulation of the flowback events. The numerical model shows that the rebound pressure is not necessarily the lower bound of the minimum principal stress. The signature of stiffness change can be identified as the process when the depletion mainly transitions from hydraulic fracture to natural fractures from numerical analysis. Overall, flowback potentially has advantages over shut-in because of the reduced time to closure.

15 GEOTHERMAL ENERGY↗

Initial Design and Experimental Results of a Novel Near-Isothermal Compressor for Heat Pump Applications

In efforts to increase the efficiency of residential and commercial air conditioners and heat pumps, it is found that the compressor has the highest electrical energy usage of the system. Therefore, it is appropriate to try to increase the efficiency of this component to reduce its energy usage. Another challenge with heat pump design is that some compressor types have drawbacks that make modulation difficult. To answer these challenges, we are developing an isothermal liquid compressor. The compressor uses propylene glycol to compress carbon dioxide. In the compression chamber the propylene glycol can enter either from the bottom to create a liquid piston for compression or it can enter through a spray nozzle at the top of the chamber. In the latter case, heat transfer from the gas to be compressed to the liquid droplets is high. This allows near isothermal operation of the compressor, which increases the efficiency by 17% to 30%, compared to adiabatic compression. In addition, the isothermal liquid compressor enables very efficient and simple part load modulation.Experimental results demonstrating the operation of the liquid compressor are presented. Initial data demonstrated a temperature rise of 7 K at pressure ratios of almost 4. For comparison at the same initial pressure, temperature, and pressure ratio, adiabatic compression would result in a temperature increase of approximately 70 K. Plotting data on a P-h diagram demonstrates that the compression started at superheated state and ended in supercritical state. Testing was later performed with repeated compressions in the superheated region of the P-h diagram at liquid flow rates of 2 x 10-3 m3/min and 2 x 10-3 m3/min to understand the limitations of the prototype for use with an actual heat pump system. This work demonstrated a novel cycle on a T-s diagram. Results from this work will be used to develop a second-generation prototype where more rapid cycling is possible.

Kowalski, Steve↗

Experimental investigation of the influence of drainage condition on change in saturation of sheared sand

Here, this paper investigates the validity of the interpretation of results from testing saturated axisymmetric triaxial compression (ATC) sand specimens utilising three-dimensional (3D) synchrotron micro-computed tomography (SMT) to probe localised events that are completely missed or misinterpreted when analysing ATC measurements based on global standard measurements. Drained and undrained experiments were conducted at low and high back-pressures (BPs) coupled with multiple in situ 3D SMT to acquire high-resolution scans of the specimens at different axial strains. Specimens tested under low BP exhibited a large pore air volume change, which was not detected by the pump system that represents standard volume measurement. The increase in air volume caused a significant reduction in the degree of saturation leading to a possible transition from saturated to partially saturated constitutive behaviour. Undrained experiments exhibited a significant volume change contrary to the assumption of negligible volumetric strain for saturated undrained experiments. Air bubbles within the shear band for drained and undrained low-BP specimens showed opposite capillary pressure responses, increased for drained and decreased for undrained cases, due to the variation in the mechanism by which each of the two experiments predominantly counters the volume expansion within the shear band.

capillary pressure↗

Comprehensive Studies on Hot Compaction and Vibration-Assisted Compaction Tests of Aluminum Powder

Abstract Aluminum powder compaction was studied using both test and simulation. Cold compaction, hot compaction, and vibration-assisted (cold) compaction tests were conducted to achieve different density ratios. First, the hot compaction test (at 300 °C, compression pressure 140 MPa) improved about 6% compared with cold compaction under the same compression pressure. Second, although the relative density ratio does not obviously improve at a vibration-assisted (cold) compaction, the strength of the specimens made under vibration loading is much better than those of cold compaction. Additionally, finite element models with well-calibrated Drucker–Prager Cap (DPC) material constitutive model were built in abaqus/standard to simulate the powder compaction process. The results of the finite element model have very good correlations with test results up to the tested range, and this finite element model further predicts the loading conditions needed to achieve the higher density ratios. Two exponential equations of the predicted density ratio were obtained by combining the test data and the simulation results. A new analytical solution was developed to predict the axial pressure versus the density ratio for the powder compaction according to DPC material model. The results between the analytical solution and the simulation model have a very good match.

Engineering↗

Neural-network learning of SPOD latent dynamics

Here, we aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower -dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

97 MATHEMATICS AND COMPUTING↗

Anisotropy evolution of elastomeric foams during uniaxial compression measured via in-situ X-ray computed tomography

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m 3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

36 MATERIALS SCIENCE↗

Combustion characteristics of diisopropoxymethane, a low-reactivity oxymethylene ether

Oxymethylene ethers (OMEs) have been studied for use as low-sooting diesel fuel additives or substitutes; very little literature discusses OMEs as spark-ignition (SI) fuels due to their typically high cetane numbers. In this work, a lower-reactivity, branched OME, diisopropoxymethane (DIPM), is evaluated to determine its effectiveness as a spark-ignition fuel, as it is the lowest-reactivity OME (as determined by Indicated Cetane Number) thus far evaluated in the literature. DIPM is synthesized in-house via acetalization from isopropanol (iPrOH) and trioxane using standard OME production practices. DIPM was then tested in a rapid compression machine (RCM) for autoignition and spark ignition characteristics, and in a modified CFR engine to determine effective octane numbers. In the RCM, an autoignition temperature sweep was performed at stoichiometric conditions from 1000/T = 1.7 - 1.0, at 5:1 inert ratio (comparable to approximately 25% EGR), where it was found that DIPM has ignition delay times 5–10x faster than isooctane and displays NTC ignition behavior. Blends with iPrOH indicate that reactivity can be matched with isooctane with low blend ratios of iPrOH in DIPM. Flame speeds were tested with a laser spark for ignition in the RCM, where the flame speed of DIPM and isooctane is determined to be comparable at engine relevant conditions. In the CFR engine, effective RON and MON based on pressure trace frequency domain measurements were determined for DIPM and a 15 vol% iPrOH in DIPM blend. Neat DIPM has (R+M)/2 = 59 and a negative sensitivity of S = -18, consistent with its NTC behavior and higher reactivity. Furthermore the DIPM/iPrOH blend has positive sensitivity and a pump-gasoline range (R+M)/2 = 89.3. DIPM on its own is unlikely to be an effective SI fuel, however, when blended with iPrOH as an ON booster, it may be a promising SI candidate fuel.

09 BIOMASS FUELS↗

A Comparative Study of Compressive Effects on the Morphology and Performance of Carbon Paper and Cloth Electrodes in Redox Flow Batteries

Compressing porous carbon electrodes is a common approach to improve flow battery performance, but the resulting impact on electrode structure, fluid dynamics, and cell performance is not well understood. Herein, microtomographic imaging, load cell testing, and flow cell diagnostics are employed to characterize how compression‐induced changes impact pressure drop, polarization, and mass‐transfer scaling. Five different compressions are tested, spanning ranges typically used in literature, for AvCarb 1071 cloth (0%, 9%, 20%, 25%, 32%) and Freudenberg H23 paper (0%, 8%, 12%, 17%, 22%). It is found that the two electrode structures have distinct responses to compression, resulting in differing optimal conditions identified for each material; specifically, the Freudenberg H23 exhibits lower combined ohmic, charge‐transfer, and mass‐transport values at 8% compression, resulting in improved electrochemical performance across all compressive values, as compared to the optimal AvCarb 1071 compression (20%). Overall, Freudenberg H23 exhibits a greater sensitivity to compression with peak electrochemical activity correlating with increased permeability, whereas AvCarb 1071 is insensitive to compressive loads but produces lower electrochemical performance. Herein, the trade‐offs of mechanical robustness on fluid‐dynamic and electrochemical performance between the two electrodes are demonstrated by the aforementioned findings, suggesting each could be used for specific operating environments.

25 ENERGY STORAGE↗

A high-order WENO-limited finite-volume algorithm for atmospheric flow using the ADER-differential transform time discretization

A high-order-accurate weighted essentially non-oscillatory (WENO) limited upwind finite-volume scheme is detailed for the compressible, nonhydrostatic, inviscid Euler equations using an arbitrary derivatives (ADER) time-stepping scheme based on differential transforms (DTs). A second-order-accurate alternating Strang dimensional splitting is compared against multidimensional simulation with 2D transport using solid body rotation of various data. The two were found to give nearly identical accuracy in orthogonal, Cartesian coordinates. Orders of convergence are demonstrated at up to ninth-order accuracy with 2D transport. 1D transport is used to confirm that error decreases monotonically with increasing order of accuracy with WENO limiting even for discontinuous data. Further, WENO limiting always decreased the error compared with simulation without limiting in the L 1 norm. A series of standard 2D compressible nonhydrostatic Euler equation test cases were validated against previous results from literature. Finally, it was demonstrated that increasing the order of accuracy led to better resolved features and increased power for kinetic energy at small wavelengths.

54 ENVIRONMENTAL SCIENCES↗

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗

Simulation of hardened cement degradation and estimation of uncertainty in predicted failure times with peridynamics

Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, here we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.

36 MATERIALS SCIENCE↗