Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FRACTURE PROPERTIES”

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 559 records · Page 31

A parametric study of fracture toughness of fibrous composite materials

Impacts to fibrous composite laminates by objects with low velocities can break fibers giving crack-like damage. The damage may not extend completely through a thick laminate. The tension strength of these damage laminates is reduced much like that of cracked metals. The fracture toughness depends on fiber and matrix properties, fiber orientations, and stacking sequence. Accordingly, a parametric study was made to determine how fiber and matrix properties and fiber orientations affect fracture toughness and notch sensitivity. The values of fracture toughness were predicted from the elastic constants of the laminate and the failing strain of the fibers using a general fracture toughness parameter developed previously. For a variety of laminates, values of fracture toughness from tests of center-cracked specimens and values of residual strength from tests of thick laminates with surface cracks were compared to the predictions to give credibility to the study. In contrast to the usual behavior of metals, it is shown that both ultimate tensile strength and fracture toughness of composites can be increased without increasing notch sensitivity.

Poe, C. C., Jr.↗

Modification of the Local Cubic Law of fracture flow for weak inertia, tortuosity, and roughness

The classical Local Cubic Law (LCL) generally overestimates flow through real fractures. We thus developed and tested a modified LCL (MLCL) which takes into account local tortuosity and roughness, and works across a low range of local Reynolds Numbers. The MLCL is based on (1) modifying the aperture field by orienting it with the flow direction and (2) correcting for local roughness changes associated with local flow expansion/contraction. In order to test the MLCL, we compared it with direct numerical simulations with the Navier-Stokes equations using real and synthetic three-dimensional rough-walled fractures, previous corrected forms of the LCL, and experimental flow tests. The MLCL performed well and the effective errors (δ) in volumetric flow rate range from -6.7% to 13.4% with an arithmetic mean of |δ| (<|δ|>) equal to 3.7%. The MLCL is more accurate than previous modifications of the LCL. We also investigated the error associated with applying the Cubic Law (CL) while utilizing modified aperture field. The δ from the CL ranges from -14.2% to 11.2%, with a slightly higher <|δ|> = 6.1% than the MLCL. The CL with the modified aperture field considering local tortuosity and roughness may also be sufficient for predicting the hydraulic properties of rough fractures.

58 GEOSCIENCES↗

Microstructural and Textural Evolutions during β → α/ω Transformation in Biomedical Ti-15Mo Alloy and Its Effects on Mechanical Performance

This study investigated the precipitation behaviors and mechanical performance of Ti-15Mo alloy subjected to different treatment conditions. The microstructures were analyzed via scanning electron microscopy and transmission electron microscopy. The tensile properties and resultant fracture surfaces were investigated to determine the relationship between the phase transformation and mechanical performance. The Ti-15Mo alloy had a notable elongation due to no observable precipitation during solution treatment plus quenching (STQ). A high ultimate tensile strength (UTS) was obtained with a considerable amount of ω precipitates transformed during air cooling after solution treatment (ST). The coarse precipitates and additional precipitation free zones in the Ti-15Mo alloy treated with furnace cooling (FC) after ST induced a decrease in UTS. Furthermore, the absence of a continuous grain boundary α phase resulted in a lower strength in the specimen annealed at 600 °C after hot rolling. The alloy iso-aged after STQ exhibited good tensile properties with a UTS of 920 MPa and an elongation of 15%. The texture and variant relationship of the precipitates was analyzed with x-ray diffraction and electron backscatter diffraction. The common variants of the β matrix influenced the nucleation of precipitates in the Ti-15Mo alloy annealed after hot rolling, and the precipitates showed an obvious {01-10}α main texture. The variant selection of α precipitates at the grain boundary was caused by the two adjacent β grains in the alloy treated by FC after ST.

36 MATERIALS SCIENCE↗

A probabilistic assessment of geomechanical reservoir integrity during CO 2 sequestration in flood basalt formations

Recent field experiments in Iceland and Washington State (USA) show that basalt formations may be favorable targets for carbon capture and sequestration (CCS) because CO 2 mineralization reactions proceed rapidly. These results imply that there is tremendous opportunity for implementing CCS in large igneous provinces. However, the magnitude of this opportunity comprises commensurate levels of uncertainty because basalt reservoirs are characterized by highly heterogeneous, fracture-controlled hydraulic properties. This geologic uncertainty is propagated as parametric uncertainty in quantitative risk models, thus limiting the efficacy of models to predict CCS performance attributes, such as reservoir integrity and storage potential. To overcome these limitations, here we present a stochastic approach for quantifying the geomechanical performance attributes of CCS operations in a highly heterogeneous basalt reservoir. We utilize geostatistical reservoir characterization to develop an ensemble of equally probable permeability distributions in a flood basalt reservoir with characteristics of the Wallula Basalt Pilot Project. We then simulate industrial-scale CO 2 injections within the ensemble and calculate the mean and variance of fluid pressure over a 1-year injection period. These calculations are combined with the state of stress in southeast Washington State to constrain the spatial extent at which shear failure, fracture initiation, and borehole breakdown may occur. Results from this study show that (i) permeability uncertainty alone causes injection pressure to vary over 25 MPa, (ii) shear failure is likely to occur at 7 times greater distances from the injection than the CO 2 migrates, and (iii) joint initiation pressures are localized within the volume comprising the CO 2 plume.

03 NATURAL GAS↗

Atomistically derived cohesive zone model of intergranular fracture in polycrystalline graphene

Pristine single crystal graphene is the strongest known two-dimensional material, and its nonlinear anisotropic mechanical properties are well understood from the atomic length scale up to a continuum description. However, experiments indicate that grain boundaries in the polycrystalline form reduce the mechanical behavior of polycrystalline graphene. Herein, we perform atomistic-scale molecular dynamics simulations of the deformation and fracture of graphene grain boundaries and express the results as continuum cohesive zone models (CZMs) that embed notions of the grain boundary ultimate strength and fracture toughness. To facilitate energy balance, we employ a new methodology that simulates a quasi-static controlled crack propagation which renders the kinetic energy contribution to the total energy negligible. We verify good agreement between Griffith's critical energy release rate and the work of separation of the CZM, and we note that the energy of crack edges and fracture toughness differs by about 35%, which is attributed to the phenomenon of bond trapping. This justifies the implementation of the CZM within the context of the finite element method (FEM). To enhance computational efficiency in the FEM implementation, we discuss the use of scaled traction-separation laws (TSLs) for larger element sizes. As a final result, we have established that the failure characteristics of pristine graphene and high tilt angle bicrystals differ by less than 10%. This result suggests that one could use a unique or a few typical TSLs as a good approximation for the CZMs associated with the mechanical simulations of the polycrystalline graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THERMO-HYDRO-MECHANICAL MODELING OF WORKING FLUID INJECTION AND THERMAL ENERGY EXTRACTION IN EGS FRACTURES AND ROCK MATRIX

Development of enhanced geothermal systems (EGS) will require creation of a reservoir of sufficient volume to enable commercial-scale heat transfer from the reservoir rocks to the working fluid. A key assumption associated with reservoir creation/stimulation is that sufficient rock volumes can be hydraulically fractured via both tensile and shear failure, and more importantly by reactivation of naturally existing fractures (by shearing), to create the reservoir. The advancement of EGS greatly depends on our understanding of the dynamics of the intimately coupled rock-fracture-fluid-heat system and our ability to reliably predict how reservoirs behave under stimulation and production. Reliable performance predictions of EGS reservoirs require accurate and robust modeling for strongly coupled thermal-hydrological-mechanical (THM) processes. Conventionally, these types of problems have been solved using operator-splitting methods, usually by coupling a subsurface flow and heat transport simulators with a solid mechanics simulator via input files. An alternative approach is to solve the system of nonlinear partial differential equations that govern multiphase fluid flow, heat transport, and rock mechanics simultaneously, using a fully coupled, fully implicit solution procedure, in which all solution variables (pressure, enthalpy, and rock displacement fields) are solved simultaneously. This paper describes numerical simulations used to investigate the poro- and thermal- elastic effects of working fluid injection and thermal energy extraction on the properties of the fractures and rock matrix of a hypothetical EGS reservoir, using a novel simulation software FALCON (Podgorney et al., 2011), a finite element based simulator solving fully coupled multiphase fluid flow, heat transport, rock deformation, and fracturing using a global implicit approach. Investigations are also conducted on how these poro- and thermal-elastic effects are related to fracture permeability evolution.

15 GEOTHERMAL ENERGY↗

Cantor-derived medium-entropy alloys: bridging the gap between traditional metallic and high-entropy alloys

The year 2004 marked the beginning of a new era in the design of metallic materials, as the concept of multiple principal component alloys, commonly known as High-Entropy Alloys (HEAs), was proposed by Cantor and Yeh. The unexpected single-phase microstructure, instead of the expected brittle intermetallic compounds, was attributed to the large entropy of mixing and immediately caught the attention of the scientific community. Today, HEAs are considered important advanced materials and a broad range of alloys using nominally the same design principle have been investigated. Despite that, the CrMnFeCoNi (Cantor) alloy stands out as the most successful HEA due to its outstanding mechanical properties and microstructure. In this scenario, variants of the Cantor alloy, named medium-entropy alloys (MEAs), are gaining significant interest as they display a better industrial potential than both HEAs and traditional alloys. These variants of the Cantor alloy with only three or four main elements result in 15 possible combinations. The microstructure of these alloys is discussed in terms of advanced characterization as well as thermodynamic parameters and computational simulation. Their phase stability is addressed over a wide range of temperatures and strain rates. The mechanical properties, especially the fracture toughness, of the CrFeCoNi and CrCoNi alloys have been reported to be even superior to those of the Cantor alloy and most modern engineering alloys. This is associated with the formation of a continuous sequence of strengthening mechanisms, including hierarchical twin networks, which serve to prolong the strain hardening. The present article reviews and critically assesses, for the first time, recent advances in these Cantor-derived MEAs.

36 MATERIALS SCIENCE↗

Chemo-Mechanics of α-V 2 O 5 During Lithiation and Implications for Rechargeable Battery Cathodes

Chemo-mechanical degradation of layered oxide electrodes is strongly influenced by crystallographic anisotropy, local stress evolution, and ion insertion, yet the intrinsic mechanical response of layered materials remains incompletely understood. Indeed, most prior studies have focused on polycrystalline materials but single crystals enable direct observation of coupling between anisotropic ion diffusion and mechanical response. This study aims to determine how crystallographic anisotropy and lithiation affect deformation, fracture, and mechanical properties in single-crystal V 2 O 5 , and compares this behavior with polycrystalline counterparts. Polycrystalline V 2 O 5 thin films and single-crystal α-V 2 O 5 were studied using nanoindentation, scanning electron microscopy, focused ion beam cross-sectioning, and Raman spectroscopy. Single crystals were tested in pristine and chemically lithiated states, including experiments in which crystals were first plastically deformed via nanoindentation and subsequently lithiated. Polycrystalline films exhibited significantly higher hardness and elastic modulus than single crystals. Single crystals indented normal to the exposed (001) basal plane exhibited pronounced anisotropic deformation, including directional slip, crystallographically-guided cracking, anisotropic crack propagation, interlayer separation, and shear localization. Lithiation caused substantial softening, reduced hardness and modulus, and suppressed displacement bursts during nanoindentation, while previously indented regions showed crack formation and growth upon lithiation. Mechanical behavior of α-V 2 O 5 is strongly governed by crystallographic anisotropy and further altered by lithiation, with pre-existing deformation serving as a strong driver of fracture during ion insertion. These findings illuminate the coupling among ion insertion, deformation, and fracture in layered oxides and provide a basis for understanding and mitigating mechanical failure in electrochemical energy-storage materials.

Anisotropy↗

Development and Demonstration of Material Properties Database and Software for the Simulation of Flow Properties in Cementitious Materials

This report describes work performed by the Savannah River National Laboratory (SRNL) in fiscal year 2014 to develop a new Cementitious Barriers Project (CBP) software module designated as FLOExcel. FLOExcel incorporates a uniform database to capture material characterization data and a GoldSim model to define flow properties for both intact and fractured cementitious materials and estimate Darcy velocity based on specified hydraulic head gradient and matric tension. The software module includes hydraulic parameters for intact cementitious and granular materials in the database and a standalone GoldSim framework to manipulate the data. The database will be updated with new data as it comes available. The software module will later be integrated into the next release of the CBP Toolbox, Version 3.0. This report documents the development efforts for this software module. The FY14 activities described in this report focused on the following two items that form the FLOExcel package; 1) Development of a uniform database to capture CBP data for cementitious materials. In particular, the inclusion and use of hydraulic properties of the materials are emphasized; and 2) Development of algorithms and a GoldSim User Interface to calculate hydraulic flow properties of degraded and fractured cementitious materials. Hydraulic properties are required in a simulation of flow through cementitious materials such as Saltstone, waste tank fill grout, and concrete barriers. At SRNL these simulations have been performed using the PORFLOW code as part of Performance Assessments for salt waste disposal and waste tank closure.

97 MATHEMATICS AND COMPUTING↗

Microstructure and Strain Rate-Dependent Tensile Deformation Behavior of Fiber Laser-Welded Butt Joints of Dual-Phase Steels

The microstructure and tensile deformation behavior of the fiber laser-welded similar and dissimilar dual-phase (DP) steel joints over a wide range of strain rates from 10{sup −3} to 10{sup 3} s{sup −1} were investigated for the further applications on the lightweight design of vehicles. The high strain rate dynamic tensile deformation process and full-field strain distribution of the base metals and welded joints were examined using the digital image correlation method and high-speed photography. The strain rate effects on the stress–strain responses, tensile properties, deformation, and fracture behavior of the investigated materials were analyzed. The yield stress (YS) and ultimate tensile strength (UTS) of the dissimilar DP780/DP980 welded joints were lying in-between those of the DP780 and DP980 base metals, and all materials exhibited positive strain rate dependence on the YS and UTS. Owing to the microstructure heterogeneity, the welded joints showed relatively lower ductility in terms of total elongation (TE) than those of the corresponding base metals. The strain localization started before the maximum load was reached, and the strain localization occurred earlier during the whole deformation process with increasing strain rate. As for the dissimilar welded joint, the strain localization tended to occur in the vicinity of the lowest hardness value across the welded joint, which was in the subcritical HAZ at the DP780 side. As the strain rate increased, the typical ductile failure characteristic of the investigated materials did not change.

36 MATERIALS SCIENCE↗

A point process model for generating biofilms with realistic microstructure and rheology

Biofilms are communities of bacteria that exhibit a multitude of multiscale biomechanical behaviours. Recent experimental advances have led to characterisations of these behaviours in terms of measurements of the viscoelastic moduli of biofilms grown in bioreactors and the fracture and fragmentation properties of biofilms. These properties are macroscale features of biofilms; however, a previous work by our group has shown that heterogeneous microscale features are critical in predicting biofilm rheology. In this paper, we use tools from statistical physics to develop a generative statistical model of the positions of bacteria in biofilms. Specifically, the model is a type of pairwise interaction model (PIM). We show through simulation that the macroscopic mechanical properties of biofilms depend on the choice of microscale spatial model. A key finding is that uniform and non-uniform sets of points lead to differing mechanical properties. This distinction appears not to have been previously considered in mathematical biofilm literature. We also found that realisations of a biologically informed PIM have realistic in silico mechanical properties, and have statistical properties that closely match experimental data. We also note that a Poisson spatial point process of suitable number density also yields realistic mechanical properties, but that the spatial distribution of points does not reflect those occurring in our experimentally observed biofilm.

Mathematics↗

Permeability Evolution and Frictional Stability of Fabricated Fractures With Specified Roughness

Roughness is widely observed on natural fractures, and its impact on the potential for induced seismicity and associated fluid migration in the subsurface remains unclear. Here we perform fracture shearing and fluid flow experiments on artificially fabricated fractures with specified roughness to investigate the role of fracture roughness on frictional properties and permeability evolution. Given the experimental conditions, we observe that rough fractures show high roughness ratio S q /L w and return higher frictional strength due to the presence of cohesive interlocking asperities. Rough fracture surfaces show velocity strengthening behavior in the initial shearing stage, which may evolve to velocity neutral and velocity weakening at greater displacements—suggesting a dynamic weakening that rough fractures become less stable with shearing. The surface roughness exerts a dominant control on permeability evolution over the entire shearing history. Permeability declines monotonically by about 2 orders of magnitude for smooth fractures. For high roughness fractures, the permeabilities evolve episodically due to cycled compaction and dilation during shearing. With a slip distance of 6 to 8 mm, permeability of the rough surface may enhance up to an order of magnitude, but significant permeability reduction may also occur for rough samples when asperities are highly worn with gouge clogging flow paths. However, there is no obvious correlation between permeability evolution and frictional behavior for rough fracture samples when fractures are subject to sudden sliding velocity changes.

58 GEOSCIENCES↗

Age-related changes in the plasticity and toughness of human cortical bone at multiple length-scales

The structure of human cortical bone evolves over multiple length-scales from its basic constituents of collagen and hydroxyapatite at the nanoscale to osteonal structures at nearmillimeter dimensions, which all provide the basis for its mechanical properties. To resist fracture, bone’s toughness is derived intrinsically through plasticity (e.g., fibrillar sliding) at structural-scales typically below a micron and extrinsically (i.e., during crack growth) through mechanisms (e.g., crack deflection/bridging) generated at larger structural-scales. Biological factors such as aging lead to a markedly increased fracture risk, which is often associated with an age-related loss in bone mass (bone quantity). However, we find that age-related structural changes can significantly degrade the fracture resistance (bone quality) over multiple lengthscales. Using in situ small-/wide-angle x-ray scattering/diffraction to characterize sub-micron structural changes and synchrotron x-ray computed tomography and in situ fracture-toughness measurements in the scanning electron microscope to characterize effects at micron-scales, we show how these age-related structural changes at differing size-scales degrade both the intrinsic and extrinsic toughness of bone. Specifically, we attribute the loss in toughness to increased non-enzymatic collagen cross-linking which suppresses plasticity at nanoscale dimensions and to an increased osteonal density which limits the potency of crack-bridging mechanisms at micron-scales. The link between these processes is that the increased stiffness of the cross-linked collagen requires energy to be absorbed by “plastic” deformation at higher structural levels, which occurs by the process of microcracking.

59 BASIC BIOLOGICAL SCIENCES↗

Universality and Stability Phase Diagram of Two-Dimensional Brittle Fracture

The two-dimensional oscillatory crack instability, experimentally observed in a class of brittle materials under strongly dynamic conditions, has been recently reproduced by a nonlinear phase-field fracture theory. Here, we highlight the universal character of this instability by showing that it is present in materials exhibiting widely different near crack tip elastic nonlinearity, and by demonstrating that the oscillations wavelength follows a universal master curve in terms of dissipation-related and nonlinear elastic intrinsic length scales. Moreover, we show that upon increasing the driving force for fracture, a high-velocity tip-splitting instability emerges, as experimentally demonstrated. The analysis culminates in a comprehensive stability phase diagram of two-dimensional brittle fracture, whose salient properties and topology are independent of the form of near tip nonlinearity.

36 MATERIALS SCIENCE↗

Machine learning reveals cyclic changes in seismic source spectra in Geysers geothermal field

The earthquake rupture process comprises complex interactions of stress, fracture, and frictional properties. New machine learning methods demonstrate great potential to reveal patterns in time-dependent spectral properties of seismic signals and enable identification of changes in faulting processes. Clustering of 46,000 earthquakes of 0.3 < ML < 1.5 from the Geysers geothermal field (CA) yields groupings that have no reservoir-scale spatial patterns but clear temporal patterns. Events with similar spectral properties repeat on annual cycles within each cluster and track changes in the water injection rates into the Geysers reservoir, indicating that changes in acoustic properties and faulting processes accompany changes in thermomechanical state. The methods open new means to identify and characterize subtle changes in seismic source properties, with applications to tectonic and geothermal seismicity.

15 GEOTHERMAL ENERGY↗

Corrosion Assessment of Candidate Materials for the SHINE Subcritical Assembly Vessel and Components FY15 Report

In the previous report of this series, a literature review was performed to assess the potential for substantial corrosion issues associated with the proposed SHINE process conditions to produce 99 Mo. Following the initial review, substantial laboratory corrosion testing was performed emphasizing immersion and vapor-phase exposure of candidate alloys in a wide variety of solution chemistries and temperatures representative of potential exposure conditions. Stress corrosion cracking was not identified in any of the exposures up to 10 days at 80°C and 10 additional days at 93°C. Mechanical properties and specimen fracture face features resulting from slow-strain rate tests further supported a lack of sensitivity of these alloys to stress corrosion cracking. Fluid velocity was found not to be an important variable (0 to ~3 m/s) in the corrosion of candidate alloys at room temperature and 50°C. Uranium in solution was not found to adversely influence potential erosion-corrosion. Potentially intense radiolysis conditions slightly accelerated the general corrosion of candidate alloys, but no materials were observed to exhibit an annualized rate above 10 μm/y.

36 MATERIALS SCIENCE↗

(Project 14-6770) An Investigation to Establish Multiphysical Property Dataset of Nuclear Materials Based on in-situ Observations and Measurements

In-core nuclear materials including fuel pins and cladding materials fail due to issues including corrosion, mechanical wear, and pellet cladding interaction. In most such scenario microstructure dependent and corrosioninduced chemistry dependent property changes significantly affect performance of cladding, pellet, and housing. Emphasis of this work was on replace conventional pellet-cladding material models with a new straingradient viscoplasticity model that is informed by transmission electron microscopy (TEM) based measurements and by nanomechanical Raman spectroscopy (NMRS) based measurements. The TEM measurements are quantitative in nature and therefore reveal stress-strain relations with simultaneous insights into mechanisms of deformation at nanoscale. The NMRS measurements reveal the similar information at mesoscale along with additional information on relating local microstructural stresses with applied stresses. The resulting information is used to fit constants in the strain gradient viscoplasticity model as well as to validate one. During TEM measurements, a micro-electro-mechanical system based setup was developed with mechanical actuation, sensing, heating, and electrical loading. Contrary to post-mortem analysis or qualitative visualization, this setup combines direct visualization of the mechanisms behind deformation with measurement of stress, strain, thermal and electrical properties. The unique research philosophy of visualizing the microstructure at high resolution while measuring the properties led to fundamental understanding in grain size and temperature effects on measured mechanical properties such as fracture toughness. A key contribution is the role of mechanical loading boundary conditions to deconvolute the insitu TEM based nanoscale and NMRS based mesoscale data to bulk behavior. First the literature based pellet cladding mechanical interaction model based on the work of Retel’s and Williamson’s in literature work to predict tempurature and stress distribution in cladding and pellet at normal operating condition was analyzed. Later the data was fitted to find constants for a viscoplastic strain gradient model. The developed model still needs to be refined and calibrated using various experimental results. That remains the focus of future work. Overall, a major thrust of the work was therefore on active control of the microstructure (grain size, defect density and types) exploiting the multi-physics coupling in materials. In particular, using experiments the synergy of current density, mechanical stress and temperature were studied to annihilate defects and recrystallize grains. The developed model is being examined for implementation in BISON. Multiple invited talks, international journal publications, and conference publications were performed by students supported on this work. Another output is support multiple PhD and masters thesis students who will be an important asset for future basic nuclear research. Future Work Recommendations: A nuclear reactor operates under significant variations of thermal loads due to energy cycling and mechanical loads due to constraint effects. Significant thermal and chemical diffusion takes place at the pallet-cladding level. While the proposed work established new experimental approach and new dataset for Zircaloy-4, the irradiation level was in the range of 1-2 dpa. Samples with higher dpa need to be examined. Therefore, a continual of support of the performed work is essential. Currently, these are the only experiments that can measure the produced data. The work also needs to be extended to different fuel types and cladding types such as SiC and FeCrAl based claddings. A combination of datasets for these materials can then be used to analyze accurately predict behavior of critical pellet cladding systems in accident scenario with high heat flux and high thermal loads. This is a BIG unknown as if now.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrated Multiphysics Modeling of Environmentally Assisted Brittle Fracture

Brittle materials, such as cement, compose major portions of built infrastructure and are vulnerable to degradation and fracture from chemo-mechanical effects. Currently, methods of modeling infrastructure do not account for the presence of a reactive environment, such as water, on the acceleration of failure. Here, we have developed methodologies and models of concrete and cement fracture that account for varying material properties, such as strength, shrinkage, and fracture toughness due to degradation or hydration. The models have been incorporated into peridynamics, non-local continuum mechanics methodology, that can model intersecting and branching brittle fracture that occurs in multicomponent brittle materials, such as concrete. Through development of new peridynamic capabilities, decalcification of cement and differential shrinkage in clay-cement composites have been evaluated, along with exemplar problems in nuclear waste cannisters and wellbores. We have developed methods to simulate multiphase phenomena in cement and cement-composite materials for energy and infrastructure applications.

36 MATERIALS SCIENCE↗