Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plasma surface modification coupled with thermal and step-over distance effects on significant fracture improvement of adhesively-bonded metal-$\mathrm{CFRTP}$ dissimilar materials

Here, this work proposes an approach to improve the interfacial bonding between carbon-fiber-reinforced thermoplastic polymers (CFRTP) and thermoset adhesive by using the surface modification with the combination of air plasma, thermal, and step-over effects. Thanks to the proper heating which causes the exposure of carbon fibers on the CFRTP surface, both polymer matrix and these exposed carbon fibers can have improved bonding with adhesive due to air plasma treatment. By conducting the Double Cantilever Beam (DCB) tests performed on adhesively-bonded AA5052/CFRP-PA66 dissimilar joints, it was shown in this work that the average Mode I specific fracture energy can be improved up to about 180% compared to non-treated counterparts. The foregoing improvement can be further increased by considering the step-over distance of the surface treatment, which triggers the tortuous damage path in the failure process of adhesively-bonded structures. By manipulating this treatment parameter, often overlooked in the literature, and other parameters (i.e., treatment speed, nozzle tip-to-surface distance, etc.) investigated in this study, the average Mode I specific fracture energy of adhesively-bonded AA5052/CFRP-PA66 dissimilar joints can be improved reaching up to about 410% compared to non-treatment counterparts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗
DOE OSTI2026

Effect of Specimen Thickness on Fracture Toughness and Plane Stress to Plane Strain Transition in Medium-Density Polyethylene

This study investigates the effect of specimen thickness on fracture toughness and the transition from plane stress to plane strain conditions in Medium-Density Polyethylene (MDPE) using Single Edge Notched Bend (SENB) specimens. Three thickness groups (7.5 mm, 9.0 mm, and 12.1 mm) were tested following ASTM D5045 protocol. Conditional stress intensity factors (KQ) increased from 2.2 MPavm to 2.8 MPavm with increasing thickness, demonstrating significant size dependency. Confocal microscopy revealed a 34% reduction in maximum crack tip opening displacement (CTOD) from 0.478 mm to 0.314 mm as thickness increased, with plastic zone lateral extent decreasing by 32%. This quantitative evidence validates Irwin's theoretical prediction of plastic zone size transition from plane stress r_y˜ (1/2p) (K_I/s_y )^2 to plane strain conditions r_y˜ (1/6p) (K_I/s_y )^2.Scanning electron microscopy confirmed progressive suppression of shear lips and evolution from ductile tearing with extensive polymer chain drawing in thin specimens to localized crazing in thick specimens. Despite all specimens satisfying the ASTM D5045 criterion (B = 2.5(KQ/s_y )^2 ) ,residual shear lips persisted even in the thickest specimens, demonstrating that nominal size requirements alone are insufficient for achieving complete plane strain conditions in highly ductile polymers. The findings emphasize the necessity of integrated analytical and morphological assessment for valid fracture toughness characterization, with critical implications for life assessment and integrity prediction in thick-section polymer components such as natural gas distribution pipelines.

Medium-Density Polyethylene↗
DOE OSTI2021

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2020

Mesoscale informed parameter estimation through machine learning: A case-study in fracture modeling

Scale bridging is a critical need in computational sciences, where the modeling community has developed accurate physics models from first principles, of processes at lower length and time scales that influence the behavior at the higher scales of interest. However, it is not computationally feasible to incorporate all of the lower length scale physics directly into upscaled models. This is an area where machine learning has shown promise in building emulators of the lower length scale models, which incur a mere fraction of the computational cost of the original higher fidelity models. We demonstrate the use of machine learning using an example in materials science estimating continuum scale parameters by emulating, with uncertainties, complicated mesoscale physics. Additionally, we describe a new framework to emulate the fine scale physics, especially in the presence of microstructures, using machine learning, and showcase its usefulness by providing an example from modeling fracture propagation. Our approach can be thought of as a data-driven dimension reduction technique that yields probabilistic emulators. Our results show well-calibrated predictions for the quantities of interests in a low-strain simulation of fracture propagation at the mesoscale level. Furthermore, on average, we achieve ~10% relative errors on time-varying quantities like total damage and maximum stresses. Successfully replicating mesoscale scale physics within the continuum models is a crucial step towards predictive capability in multi-scale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The effect of grain boundary yield and disconnections on grain boundary fracture strength

During high-temperature interfacial strain, e.g., creep or sintering, grain boundaries will sample non-equilibrium states associated with the non-equilibrium disconnections that mediate strain. This work tests the hypothesis that those non-equilibrium states should affect the mechanical properties of the grain boundary. High-temperature bicrystal tensile strain experiments performed on Sc-doped ZrO 2 as a function of strain rate enable the observation of fracture before and after grain boundary yield. Yielded boundaries exhibit a median fracture strength ≈ 50 % lower than those that did not yield. Molecular dynamics simulations performed on boundaries with and without interstitial-type disconnection loops indicate that the presence of those loops reduce the failure stress.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mineralogy and reactive fluid chemistry evolution of hydraulically fractured Caney shale of Southern Oklahoma

Here, this study investigates geochemical rock-fluid interactions as a potential cause of rapid loss of permeability and productivity in hydraulically fractured shale reservoirs. It also interrogates the effects of these reactions in transforming depleted shale reservoirs into impermeable carbon storage units. The study employs batch reactor experiments where rock-powder samples are reacted with field fracturing fluid under reservoir temperature (95°C).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Damage evolution and ductile fracture of commercially-pure titanium sheets subjected to simple tension and cyclic bending under tension

This paper describes a study into the evolution of damage in commercial-purity titanium (CP–Ti) sheets subjected to cyclic bending under tension (CBT) and uniaxial tension (simple tension, ST). Sections were taken from sheets that were strained to various levels under both methods and were imaged using X-ray computed tomography (XCT) to reveal insights into the size, density, and distribution of microvoids in the sheets. Digital image correlation (DIC) was also used to observe the surface strain of CBT and ST sheets during testing. These results showed that elongation to failure (ETF) for CBT deformation is about 2.5× higher than for ST, local longitudinal strains in the bulk of the CBT sample are around 1.3× higher than the peak strain in the ST necked region, and 1.7× higher in the localized region of CBT failure. It was found that the volume-density of voids in both sheets followed a similar exponential increase with strain, reaching nearly 45× higher in the CBT than the ST sheets before the onset of failure, mainly due to the higher strain levels achieved. A higher volume density of voids developed in the center of sheets at high levels of strain, for both processes, with the density in the sheet center becoming approximately double that found near the edges. Scanning electron microscopy (SEM) was used to examine the fracture surface of CBT and ST sheets after failure. The observations are presented and discussed highlighting the CBT process as effective to delay ductile fracture of the CP-Ti sheets.

Commercially pure titanium↗
DOE OSTI2022

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular dynamics studies of helium bubble effects on grain boundary fracture vulnerabilities in an Fe 70 Ni 11 Cr 19 –1%H austenitic stainless steel

Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Finally, detailed mechanistic analyses are performed to interpret the simulated results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mapping of fracture and ionic conductivity changes in ion implanted solid electrolytes: Insights from molecular dynamics

Ion implantation emerges as a promising technique to address the persistent challenge of lithium (Li) filament growth in solid-state electrolytes as it can induce compressive stresses inhibiting crack growth and deflect dendrites, de facto mitigating early electrolyte failure. In this study, we examine the potential paradox of ion implantation: while aiming to enhance electrolyte performance, the radiation damage associated with implantation might inadvertently compromise both the ionic conductivity and the intrinsic fracture toughness of the material, rendering the material unsuitable for battery applications. Specifically, we employed molecular dynamics simulations to examine the scope of the downsides of ion implantation, specifically: (i) reduced ionic conductivity (due to radiation-induced amorphization) and (ii) mechanical stability (due to radiation-induced embrittlement) in ion-implanted Li 7 La 3 Zr 2 O 12 (LLZO) solid-state electrolytes. We explore how radiation damage impacts LLZO’s crystalline structure, Li-ion diffusion, and fracture properties at various temperatures and radiation damage levels. The study aims to provide insights into the competing effects of ion implantation and suggest potential engineering strategies for developing more robust solid-state electrolytes with improved conductivity and dendrite resistance.

Conductivity↗
DOE OSTI2025

An FFT-based micromechanical model for gradient enhanced brittle fracture

Damage models incorporated within FFT-based micromechanical methods have received much attention recently because of the need to better understand and predict brittle and ductile fracture. An important aspect of a damage model is non-local regularization, which removes the mesh dependence of the predictions that otherwise become physically unacceptable upon grid refinement. In this work, the Helmholtz-type equation for non-local gradient regularization of a damage model on a distorted grid is solved using an FFT-based approach. Further, the resulting system of equations is solved using the Jacobi iterative method. The model is applied to simulate brittle fracture of an intermetallic. The influence of the time and space discretization, the length-scale parameter, and intermetallic crystallographic orientation on crack evolution is studied.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Implications of post-processing induced microstructural changes on the deformation and fracture response of additively manufactured Ti–6Al–4V

Ti–6Al–4V parts fabricated via electron beam melting (EBM) powder bed fusion are often subjected to hot isostatic pressing below the β-transus temperature to mitigate defects. During which the pressure aids in pore closure and the thermal exposure results in coarsening of α phase while retaining the columnar prior-β phase grain morphology present in the as-fabricated condition. The same post-processing treatment can also be carried out above β-transus temperature which is an effective way to modify both the grain morphology and the associated α textures. Here, the objective of this work is to correlate thermally induced microstructural changes to the deformation and fracture response of the EBM processed Ti–6Al–4V. To this end, we have carried out both sub-transus and super-transus heat-treatments of the as-fabricated material. The mechanical response of the as-fabricated and all heat-treated materials are characterized by in-situ tensile tests under a high-resolution digital optical microscope. This enabled us to capture large-scale panoramic images of the deforming microstructure, and overcome the trade-off between the image resolution and the field of view during in-situ experiments. The series of images captured throughout the imposed deformation are subsequently used to perform microstructure-based digital image correlation to measure microstructural-scale strains over a large area. The results of the in-situ tests together with detailed fractographic analyses are then used to elucidate how the heterogeneous deformation spanning over multiple microstructural length-scales, for example, at the scale of lamellae, colonies of lamellae and grain boundaries, affect the overall deformation and fracture response of the post-processed materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evaluating Fracture Behavior of Bioinspired Alumina-YSZ Composites through Static and Dynamic Mechanical Testing

Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.

Bioinspired materials↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2022

The mechanical behavior of metal-halide perovskites: Elasticity, plasticity, fracture, and creep

Wide ranging mechanical properties — elasticity, plasticity, fracture, and creep — most relevant to the mechanical reliability of perovskite solar cells (PSCs) are systematically investigated. High quality bulk single-crystals of the commonly studied metal-halide perovskites (MHPs) relevant to PSCs are fabricated and studied: CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and CH 3 NH 3 PbI 3 (MAPbI 3 ). The first direct measurement of MHP Young's modulus (E) using uniaxial compression reveals E <100> of 13.1 ± 1.3 and 10.6 ± 1.0 GPa for MAPbBr 3 and MAPbI 3 , respectively. Further, the Vickers micro-hardness H (100) of MAPbBr 3 and MAPbI 3 is 0.54±0.02 GPa and 0.76±0.05 GPa, respectively. The Vickers micro-indentation fracture toughness K IC of MAPbBr 3 and MAPbI 3 is estimated at 0.20±0.03 MPa∙m 0.5 and 0.18±0.03 MPa∙m 0.5 , respectively. The stress-exponent, n, extracted from nanoindentation creep data is ~8 and ~10 for MAPbBr 3 and MAPbI 3 , respectively. The trends in these properties are discussed. These properties are best estimates and are recommended for use in future mechanical behavior and reliability analyses of MHPs and PSCs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Seismic Imaging of a Shale Landscape Under Compression Shows Limited Influence of Topography–Induced Fracturing

We used seismic refraction to image the P-wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill-valley-hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. Here, this contradiction could be explained by the greater importance of infiltration-driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

58 GEOSCIENCES↗