Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied here. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Manufacturing and characterization of continuous fiber-reinforced thermoplastic tape overmolded long fiber thermoplastic

Light-weight construction, design freedom, integration of functions, and compelling cost element are desired by Original Equipment Manufacturers (OEMs) and their suppliers. The emergence of overmolding of continuous-discontinuous reinforcement enables design freedom and ability to tailor stiffness, strength, and damage tolerance for structural applications. In this work, long fiber thermoplastics (LFT) are overmolded with continuous fiber-reinforced thermoplastic (CFRTP) tape are combined using extrusion compression molding process to evaluate the structural performance. The CFRTP tape overmolded LFT samples were characterized using nondestructive and destructive techniques to track fiber alignment, fiber distribution, manufacturing defects, interfacial bonding of the tape – LFT and effect of CFRTP tape on LFT. Mode 1 fracture toughness, G 1c for tape overmolded LFT was higher by 25–30% compared to literature reported G 1c . This response was attributed to excellent fiber distribution, good fiber wetting, and absence of voids at the interface. Three-point bend test indicated that CFRTP tape overmolded LFT composites were better able to resist damage under the bending load compared to constituent LFT composite. Flexural strength of the overmolded composite was higher by 119–142%, and modulus higher by 77–65% compared to constituent LFT composite. Simulated flexural results accurately represents the mechanical behavior of composites. Here, the penetration energy of tape overmolded LFT composites determined by LVI test was in the range of 27.66–30.15 J, which is significantly higher than constituent LFT composite, 7.76 J. CFRTP tape overmolded LFT composite exhibits progressive fiber fracture, matrix cracking, and interfacial debonding failure, whereas, constituent LFT composite showed catastrophic fiber fracture.

42 ENGINEERING↗
DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating the bond strength and fracture mechanisms of cold-sprayed zinc coating on AZ91 magnesium substrate via a combined experimental and computational approach

Magnesium (Mg) alloys are ideal candidates for automotive applications due to their high strength to weight ratio, castability, recyclability etc., however, they lack corrosion and oxidation resistance. Solid-state deposition techniques, such as cold spray, have been demonstrated to enhance their corrosion resistance as it relies on the severe plastic deformation of powder particles upon impact with the substrate to form a metallurgical bond with the substrate and within the coating. At cold sprayed interfaces, a heterogeneous microstructure is formed that includes some porosity, oxides and intermetallics which can significantly affect coating performance. Thus, establishing a direct correlation between the interface microstructure and its properties can aid in designing optimal cold spray parameters. In this study, we investigated the microstructure and mechanical properties of a zinc (Zn) coating deposited on a high pressure die cast (HPDC) AZ91 Mg substrate via high resolution scanning transmission electron microscopy, in situ micro-tensile testing, and finite element method (FEM) modeling. Micro-tensile pillars fabricated using the plasma focused ion beam (PFIB) successfully isolates the coating-substrate interface within the gauge length. The average bond strength of Zn-Mg interface was determined to be ∼140 MPa with failure occurring partially at the interface and mostly into the coatings. A detailed microstructural characterization revealed evidence of a strong metallurgical bonding at the Zn-Mg interface and formation of the C14 MgZn 2 laves phase interlayer resulting in a mixed mode of fracture during the micro-tensile experiments. FEM modeling reveals the stress distribution along the interfaces and suggests that a MgZn 2 layer thickness between 200–400 nm is optimum to increase the bond strength and minimize the triaxiality. Such a site-specific interfacial analysis with correlative computational modeling provides crucial insight into the overall performance of cold spray interfaces.

Bond strength↗
DOE OSTI2021

The mechanical response of micron-sized molecular crystals

Abstract Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mixed imbibition controls the advance of wetting fluid in multiscale geological media

Imbibition plays a central role during multiphase flows in geological media; yet, our understanding of imbibition dynamics across the multiscale features of rocks, including fractures, microcracks, and heterogeneous matrices is limited. Here, we deployed a unique combination of micro- and clinical X-ray CT imaging to probe directly the spatiotemporal evolution of imbibition in a natural, multiscale shale sample. The method enabled unraveling of mixed imbibition behavior, with spontaneous imbibition rates in the matrix pores and microcracks driven by capillarity proceeding as $\propto t^{-1/2}$, whereas pressure-driven flows in the macroscopic fractures yield apparent imbibition rates with $\propto t^{-1}$ scaling. Regardless of the aperture, we showed that spontaneous imbibition proceeds uniformly away from the fractures. The effects of microcracks are further investigated by comparing the dynamics of the imbibition fronts with and without microcracks. Finally, these results highlight the paramount role of microcracks on anomalous front roughening and the need to account for microcracks appropriately in the development of predictive theoretical models.

58 GEOSCIENCES↗
DOE OSTI2022

Cracking during pyrolysis of preceramic polymers within glass microtubes

Abstract Cracking of preceramic polymers during pyrolysis under highly‐constrained conditions is examined by X‐ray computed tomography of fine glass microtubes containing the pyrolyzing material. The microtubes represent model geometries that mimic the long channels between fibers during production of ceramic composites by precursor impregnation and pyrolysis. Complementary fracture mechanics analyses of interface cracking and crack kinking are used to glean insights into the conditions under which periodic alternating cracks form. A key finding is that alternating cracks are an inherent feature of constrained pyrolysis . This feature is attributable in large part to the high energy release rates for interface cracks to kink into the pyrolyzing material under the hydrostatic tension developed during pyrolysis. It also requires interfaces with toughness comparable to that of the pyrolyzing material, to prevent large‐scale interface separation. The results further indicate the need for small uniform spaces for pyrolysis within fiber preforms in order to produce networks of fine periodic pyrolysis cracks; these networks in turn facilitate impregnation and pyrolysis in subsequent processing cycles.

Larson, Natalie M.↗
DOE OSTI2023

COMSOL Models and Simulation Results for Closed-Loop Geothermal Systems with Thermal Enhancements

The Closed-Loop Geothermal Working Group and the National Renewable Energy Laboratory (NREL) have studied the performance of disc-shaped and linear thermal enhancements (conductive material introduced around a wellbore) in closed-loop geothermal systems. COMSOL models were developed to investigate the increase in heat production in systems with thermal enhancements compared to those without, across a range of scenarios. This dataset includes two example COMSOL models, a summary of the simulation results, and a link to an accompanying paper. Studied are idealized versions of two thermal enhancement designs: (1) existing or newly created fractures surrounding the wellbore, filled with highly conductive material, and (2) a 'fishbone structure,' where short wellbores are drilled radially outward from the main wellbore and filled with highly conductive material.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Implementation and testing of physics-based pulverization model in BISON

This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multiscale Thermomechanical Damage Model with Internal State Variables for Ceramic Matrix Composites

A length scale-dependent matrix damage model has been developed to simulate the response of ceramic matrix composites (CMCs). This model captures damage initiation and propagation in the brittle matrix material by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress-strain constitutive relationships at each length scale. Material degradation due to matrix cracking is captured by a damage variable determined using fracture mechanics and the self-consistent scheme. An additional state variable is introduced to capture the nucleation and growth of matrix porosity, a key mechanism contributing to matrix nonlinearity. Porosity occurs as a result of material diffusion around grain boundaries and the growth of the porosity state variable is related to the material entropy dissipation and the volumetric strain. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide, unidirectional (UD) carbon fiber/silicon carbide matrix (C/SiC) CMC, and five-harness satin (5HS) woven C/SiC CMC.

Skinner, Travis D.↗
DOE OSTI2021

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗
DOE OSTI2025

Atomic layer deposition of nanofilms on porous polymer substrates: Strategies for success

Atomic layer deposition (ALD) is a versatile technique for engineering the surfaces of porous polymers, imbuing the flexible, high-surface-area substrates with inorganic and hybrid material properties. Previously reported enhancements include fouling resistance, electrical conductance, thermal stability, photocatalytic activity, hydrophilicity, and oleophilicity. However, there are many poorly understood phenomena that introduce challenges in applying ALD to porous polymers. In this paper, we address five common challenges and ways to overcome them: (1) entrapped precursor, (2) embrittlement, (3) film fracture, (4) deformation, and (5) pore collapse. These challenges are often interrelated and can exacerbate one another. To investigate these phenomena, we applied various ALD chemistries to porous polymers including polyethersulfone, polysulfone, polyvinylidene fluoride, and polycarbonate track-etched membranes. Reaction-diffusion modeling revealed why certain precursors and processing conditions result in embrittling subsurface material growth, entrapment of unreacted precursors, and nongrowth. We quantify the limits of ALD processing temperatures that are dictated by thermal expansion mismatch and can lead to fractured ALD films. The results herein allow us to make recommendations to avoid, mitigate, or overcome the difficulties encountered when performing ALD and plasma-enhanced ALD on porous polymers. We intend this article to serve as a “lessons learned” guide informed by previous experience to provide a better understanding of the difficulties and limitations of ALD on porous polymers and knowledge-based guidelines for successful depositions. This knowledge can accelerate future research and help experimentalists navigate and troubleshoot as they expose porous polymers to reactive precursor vapors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Osteoporosis Treatments Affect Bone Matrix Maturation in a Rat Model of Induced Cortical Remodeling

ABSTRACT To test how osteoporosis drugs affect bone matrix maturation during cortical bone remodeling, 72 pregnant rats were switched from a 0.4% to a 0.01% calcium diet at parturition for a 23‐day lactation period. At weaning, eight dams were sacrificed to establish baseline values, while the remaining dams were returned to 0.4% calcium and treated with vehicle (saline), sodium fluoride (NaF), zoledronic acid (ZA), or sclerostin antibody (Scl‐Ab) for either 7 or 28 days (eight animals per group per time point). Femora were examined by μCT, dynamic histomorphometry, Fourier transform infrared imaging, and three‐point bending of notched specimens. Cortical porosity decreased in all groups from baseline to day 28. Intracortical mineralizing surface (MS/BS) and mineral apposition rate (MAR), as well as the mineral‐to‐matrix ratio were unaffected by treatment, but intracortical crystallinity was increased in the ZA group at day 10 compared with vehicle. Cortical area increased in all groups over 28 days mainly because of an addition of bone at the endocortical surface. Endocortical MS/BS did not vary among the groups, but endocortical MAR was suppressed in the NaF group at day 2 and elevated in the Scl‐Ab group at day 4 compared with vehicle. Endocortical mineral‐to‐matrix ratio was increased at days 5 and 10 following NaF treatment and endocortical crystallinity was increased at day 5 following ZA treatment compared with vehicle. Fracture toughness did not differ among the groups. Thus, the treatments affected matrix maturation more strongly at the endocortical then intracortical envelope. In this model of induced remodeling, the bone formation phase is synchronized at multiple sites, facilitating study of the effects of drugs or other bone‐targeting agents on matrix maturation independent of their effects on the initiation of remodeling. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Geological activity shapes the microbiome in deep-subsurface aquifers by advection

Subsurface environments host diverse microorganisms in fluid-filled fractures; however, little is known about how geological and hydrological processes shape the subterranean biosphere. Here, we sampled three flowing boreholes weekly for 10 mo in a 1478-m-deep fractured rock aquifer to study the role of fracture activity (defined as seismically or aseismically induced fracture aperture change) and advection on fluid-associated microbial community composition. We found that despite a largely stable deep-subsurface fluid microbiome, drastic community-level shifts occurred after events signifying physical changes in the permeable fracture network. The community-level shifts include the emergence of microbial families from undetected to over 50% relative abundance, as well as the replacement of the community in one borehole by the earlier community from a different borehole. Null-model analysis indicates that the observed spatial and temporal community turnover was primarily driven by stochastic processes (as opposed to deterministic processes). We, therefore, conclude that the observed community-level shifts resulted from the physical transport of distinct microbial communities from other fracture(s) that outpaced environmental selection. Given that geological activity is a major cause of fracture activity and that geological activity is ubiquitous across space and time on Earth, our findings suggest that advection induced by geological activity is a general mechanism shaping the microbial biogeography and diversity in deep-subsurface habitats across the globe.

59 BASIC BIOLOGICAL SCIENCES↗