Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material balance time”

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

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of Premoderated Pulsed-Neutron Die-Away Experiments for Validating Beryllium Thermal Neutron Scattering Laws

Accurate thermal scattering laws (TSLs) are essential for reliable neutron transport simulations, particularly for systems with low-energy neutrons. Pulsed-neutron die-away (PNDA) experiments provide a highly sensitive platform for TSL validation but face challenges from fast neutron leakage and material constraints in nonhydrogenous targets such as beryllium. This work investigates the use of high-density polyethylene (HDPE) as a premoderator to reduce leakage, improve thermalization, and shorten experimental run times while preserving sensitivity to the Be TSL. Here, three premoderated configurations, encapsulation, slab, and interstitial, were evaluated using MCNP6.3® simulations with ENDF/B-VIII.0 nuclear data. The effects of the Be TSL was quantified by comparing decay constants from simulations with and without 𝑆⁡(𝛼,𝛽) treatments for both beryllium and HDPE. The results show that premoderation enables the use of as little as 5.20% of the beryllium volume required for an unreflected geometry, achieving uncertainties below 0.5% and reducing run time by up to 72.2% in the studied configurations. The interstitial configuration achieved the best balance between low statistical uncertainty and high TSL sensitivity, outperforming the slab and encapsulation designs. These findings demonstrate that carefully optimized premoderation can significantly enhance the efficiency, feasibility, and precision of PNDA experiments for Be TSL validation.

Integral Experiments↗

Use of unbalanced laminates as a screening method for microcracking

State-of-the-art, high temperature polyimide matrix composites, reinforced with continuous graphite fibers are known to be susceptible to intraply cracking when thermally cycled over their useful service temperature range. It is believed that the transply cracking, in part, results from residual stresses caused by differences in coefficients of thermal expansion (CTE) between the polymer matrix and the reinforcement. Thermal cycling tests to investigate this phenomenon involve expensive time and energy consuming programs which are not economically feasible for use as a part of a materials screening process. As an alternative to thermal cycling studies, a study of unbalanced crossply graphite fiber reinforcement composites was conducted to assess the effect of the composite ply layup and surface condition on the residual stresses that remain after the processing of these materials. The residual stresses were assessed by measuring the radii of curvature of the types of laminates that were studied. The temperature at which stress-free conditions existed were determined and a dye penetrant method was used to observe surface damage resulting from excessive residual stress buildup. These results are compared with some published results of thermal cycling tests that were previously conducted on balanced polyimide composites.

Papadopoulos, Demetrios S.↗

Use of unbalanced laminates as a screening method for microcracking

State-of-the-art, high temperature polyimide matrix composites, reinforced with continuous graphite fibers are known to be susceptible to intraply cracking when thermally cycled over their useful service temperature range. It is believed that the transply cracking, in part, results from residual stresses caused by differences in coefficients of thermal expansion (CTE) between the polymer matrix and the reinforcement. Thermal cycling tests to investigate this phenomenon involve expensive time and energy consuming programs which are not economically feasible for use as a part of a materials screening process. As an alternative to thermal cycling studies, a study of unbalanced crossply graphite fiber reinforcement composites was conducted to assess the effect of the composite ply layup and surface condition on the residual stresses that remain after the processing of these materials. The residual stresses were assessed by measuring the radii of curvature of the types of laminates that were studied. The temperature at which stress-free conditions existed were determined and a dye penetrant method was used to observe surface damage resulting from excessive residual stress buildup. These results are compared with some published results of thermal cycling tests that were previously conducted on balanced polyimide composites.

Papadopoulos, Demetrios S.↗

A real-time distributed solid oxide electrolysis cell (SOEC) model for cyber-physical simulation

System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A boundary-based approach to the multiscale microstructural characterization of a W-Ni-Fe tungsten heavy alloy

Here, a combination electron backscattered diffraction and transmission electron microscopy based approach has been implemented to study the effects of purposefully introduced anisotropy in a tungsten heavy alloy (WHA) through hot-rolling. Particular attention has been paid to changes in number and proportion of various boundary types from a quantitative standpoint; incorporating qualitative behavioral observations from prior analyses to generate experimentally-validated bases for the examination and application of a microstructure which exhibits an optimal balance of strength and ductility. It is asserted that a combination of increased temperatures during rolling and additional isothermal hold time for the post-rolling annealing steps may lead to a reduction in unfavorable textural components due to rolling in the W-phase and a decrease in premature fracture due to W-W microcracking respectively. This is expected to further increase the proportion of interphase boundaries and improve the ductility of these rolled structures, producing a superior rolled WHA microstructure.

36 MATERIALS SCIENCE↗

Mass and energy balance in the 1973 August 9 flare

The mass and energy balance of the thermal plasma during the decay phase of the solar flare of August 9, 1973, are studied. The analysis is based on observationally determined values for the differential emission measure, density, turbulent and bulk velocities, and physical dimensions. The total particle content and total thermal energy content of the flare plasmas with temperatures above 100,000 K and their variation with time are calculated. The particle loss and the energy losses through radiation, conduction, and convection are evaluated. The decrease in total particle content can be accounted for by the convective losses through the loop footprints at 100,000 K. Radiation is the dominant energy loss mechanism although convective losses at 100,000 K can be important. Conductive losses at 100,000 K into cooler chromospheric material appear to be negligible. The decrease in the total energy content during the decay phase is equal to the sum of the energy losses over the period of observation. No requirement is found for continued heating during the decay phase.

Dere, K. P.↗

An Edge Alignment-Based Orientation Selection Method for Neutron Tomography

Neutron computed tomography (nCT) is a 3D char-acterization technique used to image the internal morphology or chemical composition of samples in biology and materials sciences. A typical workflow involves placing the sample in the path of a neutron beam, acquiring projection data at a predefined set of orientations, and processing the resulting data using an analytic reconstruction algorithm. Typical nCT scans require hours to days to complete and are then processed using conventional filtered back-projection (FBP), which performs poorly with sparse views or noisy data. Hence, the main methods in order to reduce overall acquisition time are the use of an improved sampling strategy combined with the use of advanced reconstruction methods such as model-based iterative reconstruction (MBIR). In this paper, we propose an adaptive orientation selection method in which an MBIR reconstruction on previously-acquired measurements is used to define an objective function on orientations that balances a data-fitting term promoting edge alignment and a regularization term promoting orientation diversity. Using simulated and experimental data, we demonstrate that our method produces high-quality reconstructions using significantly fewer total measurements than the conventional approach.

Yang, Diyu↗

Autoclave processing for composite material fabrication. 1: An analysis of resin flows and fiber compactions for thin laminate

High quality long fiber reinforced composites, such as those used in aerospace and industrial applications, are commonly processed in autoclaves. An adequate resin flow model for the entire system (laminate/bleeder/breather), which provides a description of the time-dependent laminate consolidation process, is useful in predicting the loss of resin, heat transfer characteristics, fiber volume fraction and part dimension, etc., under a specified set of processing conditions. This could be accomplished by properly analyzing the flow patterns and pressure profiles inside the laminate during processing. A newly formulated resin flow model for composite prepreg lamination process is reported. This model considers viscous resin flows in both directions perpendicular and parallel to the composite plane. In the horizontal direction, a squeezing flow between two nonporous parallel plates is analyzed, while in the vertical direction, a poiseuille type pressure flow through porous media is assumed. Proper force and mass balances have been made and solved for the whole system. The effects of fiber-fiber interactions during lamination are included as well. The unique features of this analysis are: (1) the pressure gradient inside the laminate is assumed to be generated from squeezing action between two adjacent approaching fiber layers, and (2) the behavior of fiber bundles is simulated by a Finitely Extendable Nonlinear Elastic (FENE) spring.

Hou, T. H.↗

Defining the Operational Envelope for Air Flows in the Miniature Arc-Jet Research Chamber (mARC II)

The second-generation 30 kW miniature Arc-jet Research Chamber (mARC II) at NASA Ames Research Center produces high enthalpy flows relevant for entry systems ground testing. The mARC II facility has recently undergone upgrades, including the installation of a new vacuum system to address the issues preventing it from maintaining underexpanded flow under test conditions. In this work, we present data obtained from an Integrated Systems Testing campaign and provide an initial assessment of arc-jet performance following the upgrades. Air is used as the working gas for the standard mARC II arc-heater configuration with two constrictor disks. Seven runs were investigated for five test conditions to assess the lowest achievable stagnation point heat fluxes for air flow rates of 0.15 or 0.25 g s −1 . The heat flux was measured using a water-cooled Gardon gauge (Ø4.76 mm, 3/16" hemispherical) at 70 mm from the nozzle exit plane. The new vacuum system produced test box pressures in the medium (fine) vacuum range (~0.03 torr, 4 Pa) prior to gas addition and successfully maintained underexpanded flow after gas addition. The upgrade yielded a ~4X reduction in heat flux relative to the previous system for the same set test conditions. We report the lowest heat fluxes measured in mARC II to date, ranging from 26 to 81 W cm −2 , for sonic flow enthalpies of 4–14 MJ kg −1 . Bulk enthalpies estimated using an energy balance method (EB 2 ) are reported for the first time using mARC II.Initial data suggests EB 2 generally estimates lower enthalpies than sonic flow methods for mARC II. Lastly, laminar axisymmetric Navier–Stokes simulations were performed using the NASA DPLR code. Numerical heat flux results show good agreement with experiments at low arc powers (3–15% difference at the minimum set arc current), but discrepancy increases with arc power (49% difference at the maximum set arc current).

Pitot tube↗

Oxidation of SiC Fiber-Reinforced SiC Matrix Composites with a BN Interphase

SiC-fiber reinforced SiC matrix composites with a BN interphase were oxidized in reduced oxygen partial pressures of oxygen to simulate the environment for hypersonic vehicle leading edge applications. The constituent fibers as well as composite coupons were oxidized in oxygen partial pressures ranging from 1000 ppm O2 to 5% O2 balance argon. Exposure temperatures ranged from 816 C to 1353 C (1500 F to 2450 F). The oxidation kinetics of the coated fibers were monitored by thermogravimetric analysis (TGA). An initial rapid transient weight gain was observed followed by parabolic kinetics. Possible mechanisms for the transient oxidation are discussed. One edge of the composite coupon seal coat was ground off to simulate damage to the composite which allowed oxygen ingress to the interior of the composite. Oxidation kinetics of the coupons were characterized by scanning electron microscopy since the weight changes were minimal. It was found that sealing of the coupon edge by silica formation occurred. Differences in the amount and morphology of the sealing silica as a function of time, temperature and oxygen partial pressure are discussed. Implications for use of these materials for hypersonic vehicle leading edge materials are summarized.

Opila, Elizabeth↗

Analyzing Responses of Chemical Sensor Arrays

NASA is developing a third-generation electronic nose (ENose) capable of continuous monitoring of the International Space Station s cabin atmosphere for specific, harmful airborne contaminants. Previous generations of the ENose have been described in prior NASA Tech Briefs issues. Sensor selection is critical in both (prefabrication) sensor material selection and (post-fabrication) data analysis of the ENose, which detects several analytes that are difficult to detect, or that are at very low concentration ranges. Existing sensor selection approaches usually include limited statistical measures, where selectivity is more important but reliability and sensitivity are not of concern. When reliability and sensitivity can be major limiting factors in detecting target compounds reliably, the existing approach is not able to provide meaningful selection that will actually improve data analysis results. The approach and software reported here consider more statistical measures (factors) than existing approaches for a similar purpose. The result is a more balanced and robust sensor selection from a less than ideal sensor array. The software offers quick, flexible, optimal sensor selection and weighting for a variety of purposes without a time-consuming, iterative search by performing sensor calibrations to a known linear or nonlinear model, evaluating the individual sensor s statistics, scoring the individual sensor s overall performance, finding the best sensor array size to maximize class separation, finding optimal weights for the remaining sensor array, estimating limits of detection for the target compounds, evaluating fingerprint distance between group pairs, and finding the best event-detecting sensors.

Zhou, Hanying↗

Composition, Activity, and Stability of IrO x Oxygen Evolution Reaction Electrocatalysts

The oxygen evolution reaction (OER) is integral to several electrochemical energy conversion and storage technologies, including carbon dioxide reduction to value added fuels, nitrogen reduction to ammonia, reversible fuel cells, rechargeable metal−air batteries, and water electrolysis to produce hydrogen. Iridium oxide (IrO x ) is widely recognized as the benchmark OER catalyst for acidic environments. Despite widespread use of IrO x catalysts, most notably in proton-exchange membrane water electrolyzers (PEMWEs), a comprehensive understanding of the physicochemical properties of commercial catalysts and the impact of these properties on both the activity and stability of these catalysts is lacking. Here, we study commercial IrO x catalysts with different physicochemical properties, three nominally considered amorphous and three rutile, to elucidate how structural and compositional variations affect OER activity and stability. Utilizing standardized aqueous electrochemical protocols, time-resolved dissolution quantification using inductively-coupled plasma mass spectrometry, and physicochemical characterization, including multiple synchrotron X-ray techniques, we systematically correlate catalyst properties with OER performance and degradation behavior aided by principal component analysis (PCA). Our results demonstrate the general trend of amorphous IrO x having higher intrinsic activity but limited stability and crystalline rutile IrO 2 having lower activity but enhanced stability against dissolution. The trends within the amorphous and rutile catalyst groups correlate with inherent material properties, including phase composition and structure, crystallinity, particle size, surface area, and surface structure/chemistry. Notably, we identify a rutile catalyst with the largest crystallite/ domain sizes, moderate surface area, a small fraction of hydrous phase, and a favorable pore structure (trimodal distributions of pore sizes ranging from 2−5 nm) that exhibits the best balance between activity and stability among the six catalysts studied here. These findings illustrate a fundamental structure-governed trade-off between activity and stability and highlight the critical role of surface chemistry modification and structure engineering in IrO x catalyst optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatially varying colloidal phase behavior on multi-dimensional energy landscapes

A method is reported to determine equilibrium concentration profiles and local phase behavior for colloids on multi-dimensional energy landscapes. A general expression is derived based on local particle concentration and osmotic pressure differences that are balanced by forces on colloids due to energy landscape gradients. This analysis is applied to colloidal particles in high frequency AC electric fields within octupolar electrodes, where the energy landscape can be shaped in two dimensions. These results are also directly applicable to any particles having induced dipoles in spatially non-uniform electromagnetic fields. Predictions based on modelling colloids with an effective hard disk equation of state indicate inhomogeneous solid and fluid states coexisting on different shaped energy landscapes including multiple minima. Model predictions show excellent agreement with time-averaged Brownian Dynamic simulations at equilibrium. Findings demonstrate a general approach to understand colloidal phase behavior on energy landscapes due to external fields, which could enable control of colloidal microstructure on morphing energy landscapes and the inverse design of fields to assemble hierarchically structured colloidal materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biological removal of gaseous ammonia in biofilters: space travel and earth-based applications

Gaseous NH3 removal was studied in laboratory-scale biofilters (14-L reactor volume) containing perlite inoculated with a nitrifying enrichment culture. These biofilters received 6 L/min of airflow with inlet NH3 concentrations of 20 or 50 ppm, and removed more than 99.99% of the NH3 for the period of operation (101, 102 days). Comparison between an active reactor and an autoclaved control indicated that NH3 removal resulted from nitrification directly, as well as from enhanced absorption resulting from acidity produced by nitrification. Spatial distribution studies (20 ppm only) after 8 days of operation showed that nearly 95% of the NH3 could be accounted for in the lower 25% of the biofilter matrix, proximate to the port of entry. Periodic analysis of the biofilter material (20 and 50 ppm) showed accumulation of the nitrification product NO3- early in the operation, but later both NO2- and NO3- accumulated. Additionally, the N-mass balance accountability dropped from near 100% early in the experiments to approximately 95 and 75% for the 20- and 50-ppm biofilters, respectively. A partial contributing factor to this drop in mass balance accountability was the production of NO and N2O, which were detected in the biofilter exhaust.

NASA Discipline Life Support Systems↗

Exploring filamentous fungi depolymerization of corn stover in the context bioenergy queuing operations

Recalcitrance of lignocellulosic feedstocks to depolymerization is a significant barrier for bioenergy production approaches that require conversion of monomeric carbohydrates to renewable energy sources. This study assesses how low-cost modifications in the feedstock supply chain can be transformed into targeted pretreatments in the context of the entire bioenergy supply chain. The aim of this research is to overcome the physiochemical barriers in corn stover that necessitate increased severity in conversion in terms of chemical loading, temperature, and residence time. Corn stover samples were inoculated with a selective (Ceriporiopsis subvermispora) and non-selective (Phaenarochaete chrysosporium) lignin degrading filamentous fungal strains, then stored aerobically to determine the working envelope for fungal pretreatment to achieve lignin degradation. Dry matter loss and gross chemical makeup of corn stover varied by the length of treatment (2 and 4 weeks) and by the moisture content of the treated corn stover samples (40 and 60%, wet basis). Dry matter loss in P. chrysosporium inoculated biomass was elevated compared to the C. subvermispora inoculated biomass; however, treatment also induced additional chemical composition changes suggestive of depolymerization. Scanning electron microscope images reveal hyphae attached within cell lumen and suggest structural changes within P. chrysosporium treated corn stover after 60% moisture storage. These results highlight that fungal treatment approaches must balance loss of convertible material with the potential for reduction in recalcitrance. Techno-economic assessment (TEA) of fungal pretreatment in a short-term queuing system indicated the viability of this approach compared to conventional queuing operations. The total queuing system cost was estimated at $\$$1.65/tonne of biomass stored. After applying the credit of $\$$1.48/tonne from energy savings in the conversion phase using fungal pretreated biomass, the total system cost was $0.80 lower than traditional biomass queueing approach. While the TEA results suggested that treating biomass with C. subvermispora is the most economically viable storage method in the designed fungal-assisted queuing system, future research should focus on additional fungal depolymerization such as those observed in the P. chrysosporium inoculated biomass.

09 BIOMASS FUELS↗

Capacitance of Carbon Nanotube/Graphene Composite Electrodes with [BMIM + ][BF$_4^–$]/Acetonitrile: Fixed Voltage Molecular Dynamics Simulations

The capacitance of nanoporous carbon electrode materials is dictated by the physical interactions with electrolyte ions and molecules at the accessible interior electrode surface. While significant progress has been made in designing and synthesizing carbon materials with well-defined and relatively homogeneous nanoporosity, the majority of materials remain heterogeneous, and the resulting properties reflect an average over this distribution. In this regard, computer simulations can be a valuable tool to predict or validate structure/property relationships by systematically investigating well-defined, model electrode morphologies. In this work, we utilize fixed-voltage molecular dynamics simulations to predict structure/capacitance relationships for five different model morphologies of carbon nanotube/graphene (CNT/G) composite electrodes with 1-butyl-3-methylimidazolium tetrafluoroborate/acetonitrile electrolyte. The CNT/G electrode models are inspired by experimental “layer-by-layer” syntheses and strike a balance between realism and computational tractability. Comparison between different model electrode architectures elucidates important structure/property relationships. We find that CNT/graphene contact points serve as “hot spots” with significantly enhanced charge separation relative to the rest of the electrode. Furthermore, we demonstrate a specific nanoconfinement motif that provides substantial 3–4$\times$ enhancement of local capacitance, resulting in a ~40% increase of the total electrode differential capacitance. Because the accessible surface area of the model CNT/G electrodes is precisely determined, a comparison of per-area capacitance across systems is unambiguous. Our results thus complement a prior computational demonstration of capacitance enhancement in nanoconfinement while elucidating additional interaction motifs at CNT/G electrochemical interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗