Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FRACTURES”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Variable resolution Poisson-disk sampling for meshing discrete fracture networks

Here, we present the near-Maximal Algorithm for Poisson-disk Sampling (nMAPS) to generate point distributions for variable resolution Delaunay triangular and tetrahedral meshes in two and three-dimensions, respectively. nMAPS consists of two principal stages. In the first stage, an initial point distribution is produced using a cell-based rejection algorithm. In the second stage, holes in the sample are detected using an efficient background grid and filled in to obtain a near-maximal covering. Extensive testing shows that nMAPS generates a variable resolution mesh in linear run time with the number of accepted points. We demonstrate nMAPS capabilities by meshing three-dimensional discrete fracture networks (DFN) and the surrounding volume. The discretized boundaries of the fractures, which are represented as planar polygons, are used as the seed of 2D-nMAPS to produce a conforming Delaunay triangulation. The combined mesh of the DFN is used as the seed for 3D-nMAPS, which produces conforming Delaunay tetrahedra surrounding the network. Under a set of conditions that naturally arise in maximal Poisson-disk samples and are satisfied by nMAPS, the two-dimensional Delaunay triangulations are guaranteed to only have well-behaved triangular faces. While nMAPS does not provide triangulation quality bounds in more than two dimensions, we found that low-quality tetrahedra in 3D are infrequent, can be readily detected and removed, and a high-quality balanced mesh is produced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

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

A non-uniform corrosion model and meso-scale fracture modelling of concrete

Corrosion-induced concrete cracking is a significant durability problem for reinforced concrete structures. Considerable research has been carried out in the last few decades to understand and model the expansion mechanism of the corrosion products around the reinforcing bar and simulate the cracking behaviour of the concrete cover. In this paper, a new corrosion model based on non-uniform corrosion expansion is formulated and validated against experimental data. A meso-scale fracture model, consisting of aggregates, cement paste/mortar and ITZ, is established for the cases of both middle and side reinforcing bars. Under the developed corrosion and concrete fracture model, the cracking phenomena of the concrete cover are accurately simulated. It has been found that the non-uniform corrosion model can be used to express the realistic corrosion rust progression around the reinforcing bar, with the best accuracy. It has also been found that some microcracks occur before they are connected to form the dominating discrete crack which usually appears on the concrete surface. Moreover, the effects of the corrosion variables, as well as other key material and geometric parameters, on surface cracking of concrete are investigated.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep reinforcement learning control of hydraulic fracturing

Hydraulic fracturing is a technique to extract oil and gas from shale formations, and obtaining a uniform proppant concentration along the fracture is key to its productivity. Recently, various model predictive control schemes have been proposed to achieve this objective. But such controllers require an accurate and computationally efficient model which is difficult to obtain given the complexity of the process and uncertainties in the rock formation properties. In this article, we design a model-free data-based reinforcement learning controller which learns an optimal control policy through interactions with the process. Deep reinforcement learning (DRL) controller is based on the Deep Deterministic Policy Gradient algorithm that combines Deep-Q-network with actor-critic framework. In addition, we utilize dimensionality reduction and transfer learning to quicken the learning process. We show that the controller learns an optimal policy to obtain uniform proppant concentration despite the complex nature of the process while satisfying various input constraints.

42 ENGINEERING↗
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 OSTI2016

Modeling of heat extraction from variably fractured porous media in Enhanced Geothermal Systems

Modeling of heat extraction in Enhanced Geothermal Systems is presented. The study builds on recent studies on the use of directional wells to improve heat transfer between doublet injection and production wells. The current study focuses on the influence of fracture orientation on production temperature in deep low permeability geothermal systems, and the effects of directional drilling and separation distance between boreholes on heat extraction. The modeling results indicate that fracture orientation with respect to the well-pair plane has significant influence on reservoir thermal drawdown. As a result, the vertical well doublet is impacted significantly more than the horizontal well doublet

15 GEOTHERMAL ENERGY↗
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 OSTI2012

Radiation-Induced Rib Fractures After Hypofractionated Stereotactic Body Radiation Therapy: Risk Factors and Dose-Volume Relationship

Purpose: The purpose of this study was to clarify the incidence, the clinical risk factors, and the dose-volume relationship of radiation-induced rib fracture (RIRF) after hypofractionated stereotactic body radiation therapy (SBRT). Methods and Materials: One hundred sixteen patients treated with SBRT for primary or metastatic lung cancer at our institution, with at least 6 months of follow-up and no previous overlapping radiation exposure, were included in this study. To determine the clinical risk factors associated with RIRF, correlations between the incidence of RIRF and the variables, including age, sex, diagnosis, gross tumor volume diameter, rib-tumor distance, and use of steroid administration, were analyzed. Dose-volume histogram analysis was also conducted. Regarding the maximum dose, V10, V20, V30, and V40 of the rib, and the incidences of RIRF were compared between the two groups divided by the cutoff value determined by the receiver operating characteristic curves. Results: One hundred sixteen patients and 374 ribs met the inclusion criteria. Among the 116 patients, 28 patients (46 ribs) experienced RIRF. The estimated incidence of rib fracture was 37.7% at 3 years. Limited distance from the rib to the tumor (<2.0 cm) was the only significant risk factor for RIRF (p = 0.0001). Among the dosimetric parameters used for receiver operating characteristic analysis, the maximum dose showed the highest area under the curve. The 3-year estimated risk of RIRF and the determined cutoff value were 45.8% vs. 1.4% (maximum dose, {>=}42.4 Gy or less), 51.6% vs. 2.0% (V40, {>=}0.29 cm{sup 3} or less), 45.8% vs. 2.2% (V30, {>=}1.35 cm{sup 3} or less), 42.0% vs. 8.5% (V20, {>=}3.62 cm{sup 3} or less), or 25.9% vs. 10.5% (V10, {>=}5.03 cm{sup 3} or less). Conclusions: The incidence of RIRF after hypofractionated SBRT is relatively high. The maximum dose and high-dose volume are strongly correlated with RIRF.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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 OSTI2015

Groundwater flow characterization in a fractured bedrock aquifer using active DTS tests in sealed boreholes

In recent years, wireline temperature profiling methods have evolved to offer new insight into fractured rock hydrogeology. Important advances in wireline temperature logging in boreholes make use of active line source heating alone and then in combination with temporary borehole sealing with flexible impervious fabric liners to eliminate the effects of borehole cross-connection and recreate natural flow conditions. Here, a characterization technique was developed based on combining fiber optic distributed temperature sensing (DTS) with active heating within boreholes sealed with flexible borehole liners. DTS systems provide a temperature profiling method that offers significantly enhanced temporal resolution when compared with conventional wireline trolling-based techniques that obtain a temperature–depth profile every few hours. The ability to rapidly and continuously collect temperature profiles can better our understanding of transient processes, allowing for improved identification of hydraulically active fractures and determination of relative rates of groundwater flow. Here, the advantage of a sealed borehole environment for DTS-based investigations is demonstrated through a comparison of DTS data from open and lined conditions for the same borehole. Evidence for many depth-discrete active groundwater flow features under natural gradient conditions using active DTS heat pulse testing is presented along with high resolution geologic and geophysical logging and hydraulic datasets. Implications for field implementation are discussed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2012

Fracture strength and elastic modulus as a function of porosity for hydroxyapatite and other brittle materials Part I

Porous brittle materials are used as filters, catalyst supports, solid oxide fuel cells and biomedical materials. However the literature on the Weibull modulus, m, versus volume fraction porosity, P, is extremely limited despite the importance of m as a gauge of mechanical reliability. In Part I of this study, m is determined for 441 sintered hydroxyapatite (HA) specimens fractured in biaxial flexure for 0.08 P 0.62. In this study, we analyze a combined data set collected from the literature that represents work from a total of 17 different research groups (including the present authors), eight different materials and more than 1560 oxide and non-oxide specimens, the m versus P plot is Ushaped with a wide band of m values for P < 0.1 (Region I) and P > 0.55 (Region III), and a narrower band of m values in the intermediate porosity region of 0.1 < P < 0.55 (Region II). The limited range of m ( 4 < m < 11) in Region II has important implications since Region II includes the P range for the majority of the applications of porous brittle materials. Part II of this study focuses on the P dependence of the mean fracture strength,< f >, and the Young s modulus E for the HA specimens tested in Part I along with literature data for other brittle materials. Both < f > and E are power law functions of the degree of densification, , where = 1 P/PG and PG is the green (unfired) porosity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2012

Fracture strength and elastic modulus as a function of porosity for hydroxyapatite and other brittle materials, Part II

Part I of this paper discussed the Weibull modulus m, versus porosity P behavior of brittle materials, including HA. While the Weibull modulus m deals with the scatter in fracture strength data, this paper (Part II) focuses on two additional key mechanical properties of porous materials, namely the average fracture strength f , and Young s modulus E, for P in the interval from P zero to P PG (the porosity of the unfired compacts). The f versus P data for HA from this study and the literature data for alumina, yttria stabilized zirconia (YSZ) and silicon nitride are describedwell by functions of , where = 1 P/PG = the degree of densification. A similar function of applies to the E versus P behavior of HA from this study and data from the literature for alumina, titanium and YSZ. All of the data analyzed in this study (Part II) are based on partially and fully sintered powder compacts (excluding green powder compacts), thus the f / 0 versus and E /E0 versus relationships may apply only to such specimens.

36 MATERIALS SCIENCE↗