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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Failure Strain and Related Triaxiality of Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 Metals, Part I: Experimental Investigation

The objective of this study is to develop failure-limit material models for Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 metals, based on parameters of plastic equivalent strain (failure strain) and stress triaxiality. The research is conducted in two parts. This paper presents Part One of the study. In Part One, custom-designed test specimens undergo controlled uniaxial tension and compression testing at ambient temperature. These tests are performed at quasi-static speeds using Universal Testing Machines (UTMs) in accordance with ASTM E8 and ASTM E9 standards. Experimental data, specifically engineering stress–strain and force–displacement curves, are recorded from the onset of loading until specimen fracture, or in the case of compression tests, until the capacity of the testing machine is reached. In Part Two, the emphasis shifts to the calibration of Finite Element Analysis (FEA) models of the custom-designed test specimens. Plastic equivalent strain and the corresponding stress triaxiality values at failure are extracted from each test specimen for the given metal. These values are then systematically plotted onto a single graph to construct the failure-limit curve, which delineates the boundary conditions for material failure. This approach will facilitate the development of a comprehensive material property definition that correlates plastic equivalent strain with stress triaxiality at failure for Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 metals.

Harwell, Ron↗
DOE OSTI2021

Advances and challenges in CO2 foam technologies for enhanced oil recovery in carbonate reservoirs

Utilization of surfactants for generating foam has the potential for fluid-mobility control during CO2 flooding processes, leading to improved oil recovery. This manuscript is a review of the most important aspects for the design of CO2-foams with mobility control in carbonate reservoirs, including recent advances in novel surfactants and analytical techniques for analysis of surfactant adsorption on carbonate minerals and their thermal stability. Several key parameters and properties regarding foam transport in porous media are reviewed such as, the minimum pressure gradient required for foam generation, the effects of the partition coefficient of a CO2-soluble surfactant on foam transport in rock cores, visualization of foam flow by either microfluidic or PET/CT imaging processes in heterogeneous or fractured systems and the effect of oil on foam transport under reservoir conditions. Also are discussed the most recent advances of two major foam modeling methods, including the semi-empirical STARS foam model and the population balance model, as well as CO2 foam pilot tests and the factors governing their successes.

Jian, Guoqing↗
DOE OSTI2021

Employing Weibull Analysis and Weakest Link Theory to Resolve Crystalline Silicon PV Cell Strength Between Bare Cells and Reduced- and Full-Sized Modules

Weibull analysis and weakest link theory are employed to resolve the probability of crystalline silicon PV cell fracture when measured as bare cells and when stressed in reduced- and full-sized modules. Experimental results indicated that the characteristic cell strength is reduced by ~20% once packaged into the laminate of a one-cell module and loaded in four-point flexure (4PF). This experimental observation was shown consistent with a weakest link theory prediction that the strength limiting flaws reside on the surface of the cell's edge. Here, the analysis is ultimately extended to present the equivalent loading of four-cell modules by uniform pressure and 4PF and a uniformly loaded full-sized module and demonstrates that smaller, representative, modules must be loaded to a much higher level than their parent full-sized modules to achieve an equivalent probability for cell fracture.

14 SOLAR ENERGY↗
DOE OSTI2021

EGS Collab Experiment 1: Circulation Testing Processed data

This submission includes processed and reduced data for circulation testing that was conducted at the 164' fracture on the 4850 ft level of the Sanford Underground Research Facility. The circulation tests were done to test the flow through the 164' fracture in the EGS Collab Experiment 1 testbed on the 4850 ft level of the Sanford Underground Research Facility in order to validate computer models concerned with flow processes and heat exchange processes in EGS. This data was processed from raw data that is available in the GDR submission linked below: "EGS Collab Circulation Testing Raw Data and Documentation." Python scripts used to process and reduce the data are included. The scripts also serve as "documentation" as well as examples for data processing. More information about the test itself as well as the rationale and process of data processing is available at wiki page linked below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Unlocking the Electrochemical–Mechanical Coupling Behaviors of Dendrite Growth and Crack Propagation in All-Solid-State Batteries

Dendrite growth and crack propagation are two major hurdles on the road towards the large-scale commercialization of lithium metal all-solid-state batteries (ASSBs). Due to the high multiphysics coupled nature of the underlying dendrite growth mechanism, understanding it has been difficult. Herein, for the first time, an electrochemical-mechanical model is established that directly couples dendrite growth and crack propagation from a physics-based perspective at the cell level. Results reveal that overpotential-driven stress propels a crack to penetrate through the solid electrolyte, creating vacancies for dendrite growth, leading to the short circuit of the battery. Thus, high lithiation/charging rate and low conductivity of electrolytes can accelerate the electrochemical failure of the battery. It is further discovered that Young's modulus E LLZO of the electrolyte has competing contributions to the fracture and dendrite growth; specifically, when E LLZO = 40-100 GPa, the short circuit is triggered early. A larger toughness value hinders the crack propagation and mitigates the Li dendrite growth. The developed multiphysics model provides an in-depth understanding of the coupling of crack propagation and dendrite growth within ASSBs and an insightful mechanistic design guidance map for robust and safe ASSB cells.

25 ENERGY STORAGE↗
DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2023

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Fully Bayesian Analysis With Model Inadequacy Correction For Nuclear Graphite Property Models With Hierarchical Variance Structure

Nuclear-grade graphites are extensively utilized in the core designs of various advanced nuclear reactors. Within the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, such as nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity, forecasting component performance over the reactor's lifespan, and developing design standards necessitate robust tools for predicting fracture initiation and propagation in graphite structural components within nuclear reactors. This code enables the Bayesian calibration of properties for nuclear-grade graphites. Using a hierarchical Bayesian approach, multiple experimental data sources are combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due inadequacies in the model and the inherent spread in the experimental data are quantified.

Dhulipala, Som Lakshmi NarasimhaLakshmi Narasimha ↗
DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗
DOE OSTI2021

Experimental Study of Microgel Conformance-Control Treatment for a Polymer-Flooding Reservoir Containing Superpermeable Channels

Summary Polymer flooding has been widely used to improve oil recovery. However, its effectiveness would be diminished when channels (e.g., fractures, fracture-like channels, void-space conduits) are present in a reservoir. In this study, we designed a series of particular sandwich-like channel models and tested the effectiveness and applicable conditions of micrometer-sized preformed particle gels (PPGs, or microgels) in improving the polymer-flooding efficiency. We studied the selective penetration and placement of the microgel particles, and their abilities for fluid diversion and oil-recovery improvement. The results suggest that polymer flooding alone would be inefficient to achieve a satisfactory oil recovery as the heterogeneity of the reservoir becomes more serious (e.g., permeability contrast kc/km > 50). The polymer solution would vainly flow through the channels and leave the majority of oil in the matrices behind. Additional conformance-treatment efforts are required. We tried to inject microgels in an attempt to shut off the channels. After the microgel treatment, impressive improvement of the polymer-flooding performance was observed in some of our experiments. The water cut could be reduced significantly by as high as nearly 40%, and the sweep efficiency and overall oil recovery of the polymer flood were improved. The conditions under which the microgel-treatment strategy was effective were further explored. We observed that the microgels form an external impermeable cake at the very beginning of microgel injection and prevent the gel particles from entering the matrices. Instead, the microgel particles could selectively penetrate and shut off the superpermeable channels under proper conditions. Our results suggest that the 260-µm microgel particles tested in this study are effective to attack the excessive-water-production problem and improve the oil recovery when the channel has a high permeability (>50 darcies). The gels are unlikely to be effective for channels that are less than 30 darcies because of the penetration/transport difficulties. After the gels effectively penetrate and shut off the superpermeable channel, the subsequent polymer solution is diverted to the matrices (i.e., the unswept oil zones) to displace the bypassed oil. Overall, this study provides important insights to help achieve successful polymer-flooding applications in reservoirs with superpermeable channels.

Engineering↗
DOE OSTI2025

Library of Advanced Materials for Engineering (LAMÉ) (V.5.24)

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options an

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plasticity and fracture behavior of Inconel 625 manufactured by laser powder bed fusion: Comparison between as-built and stress relieved conditions

In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.

36 MATERIALS SCIENCE↗
DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Investigating Material Properties of Subsurface Rock Formations Modified by Engineering Mineral Precipitation (Final Scientific and Technical Report)

Montana State University’s (MSU) Energy Research Institute (ERI), in collaboration with the Center for Biofilm Engineering (CBE) and the Department of Civil Engineering (CE), has conducted a long‐term research program aimed at developing a novel cementing agent to address wellbore integrity and reduce the unwanted upward migration of fluids and greenhouse gases from the subsurface. The primary technology developed through this research program is known as ureolysis‐induced calcite precipitation (UICP), which harnesses bio‐chemical processes to precipitate calcium carbonate (CaCO 3 ). The same general process can also be called microbially-induced calcium carbonate precipitation (MICP) when microbes provide the process-catalyzing urease enzyme. Both terms are used in this report. Results have conclusively demonstrated that, if properly controlled, UICP can successfully seal fractures, high permeability zones, and compromised cement in the vicinity of wellbores and in nearby caprock. This technology has been successfully deployed to mitigate annular leakage in two test wells and over sixty commercial wells with a 100% success rate. This success in downhole deployment generates consideration of other subsurface applications where UICP could provide benefit to the energy sector, such as shale property modification for unconventional oil and gas recovery. The focus of this research project was to investigate fundamental material and mechanical properties of select shale cores and analyze how these properties change due to engineered mineral precipitation with the intent to control these properties to achieve a range of engineering objectives. Ultimately, the project aim was to identify valuable new areas where application of UICP might contribute to national energy security and environmental protection. The research workplan coupled UICP treatment of core samples, nuclear magnetic resonance (NMR) characterization, and mechanical strength testing at MSU with advanced X‐Ray micro-computed tomography (μCT) imaging and numerical modeling performed by collaborators at two national laboratories, the National Energy Technology Laboratory (NETL) and Lawrence Berkeley National Laboratory (LBNL). Experimental results are useful to inform geo-mechanical models which could be applied to predict mineralized rock formation behavior at field scale. Our findings suggest that NMR and μCT methods to detect and quantify biomineral formation in shale fractures are complementary and consistent with each other. Either could be used to estimate the volume of new mineral formed by UICP in shale fractures. The use of surfactants and guar gum to enhance biomineral precipitation in shale fractures merits further research. UICP can, under some conditions, increase the tensile strength of sealed shale fractures beyond that of the intact shale. These findings demonstrate that continued research in this area may be valuable to understanding and improving shale resource recovery techniques.

58 GEOSCIENCES↗
DOE OSTI2025

A novel peridynamics-based approach to predict pharmaceutical tablet robustness

The pharmaceutical drug product development process can be greatly accelerated through the use of modeling and simulation techniques to predict the manufacturability and performance of a given formulation. The anticipation and possible mitigation of tablet damage due to manufacturing stresses represents a specific area of interest in the pharmaceutical industry for predicting formulation and tableting performance. While the finite element method (FEM) has been extensively used for predicting the mechanical behavior of powder material in the compaction processes, a shortcoming of the approach is the inherent difficulty to predict discontinuities (e.g., damage or cracking) within a tablet as FEM is a continuum-based approach. In this work, we propose a novel method utilizing peridynamics (PD), a numerical method that can capture discontinuities such as tablet fracture, to predict the evolution of damage and breakage in pharmaceutical tablets. The approach links (1) the finite element method – to elucidate the behavior of powders during die compaction – with (2) the peridynamics modeling technique – to model the discontinuous nature of damage and predict tablet breakage during the critical stages of unloading and ejection from the compression die. This short communication presents a proof of concept including a workflow to calibrate the linked FEM-PD simulation models. Further, it demonstrates promising results from a preliminary experimental validation of the approach. Following further development, this approach could be used to guide the optimization of compression processes through targeted changes to formulation material properties, compression process conditions, and/or tooling geometries to deliver improved process efficiency and tablet robustness.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deterministic and Statistical Modeling of a New Thermal Breakout Technology for Measuring the Maximum Horizontal In-Situ Stress

The current state-of-the-art technology for in-situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., “mini-fracs”). This methodology has achieved wide application in the oil and gas industry but has several limitations that often prevent successful in-situ stress measurements. One major limitation is that breakouts do not appear in all boreholes and are generally only a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to directly measure the maximum horizontal in-situ stress magnitude, the absence of borehole breakouts presents a major data gap for in-situ stress measurements. In response to this data gap, a new US Department of Energy (US DOE)-sponsored thermal breakout technology that will provide a method for thermally inducing borehole breakouts and allow the consistent measurement of the maximum horizontal stress magnitude is currently in development. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which can be directly correlated to the maximum horizontal stress magnitude. The first step in this project was an analytical modeling study of the thermal breakout process. Based on the Kirsch solution (Kirsch 1898), a deterministic and statistical analysis was performed on the pertinent parameters that influence the maximum horizontal stress calculation. As a result of the analysis, the findings indicate that the thermal breakout technology is feasible and provides improved accuracy and/or an enhanced ability to measure the maximum horizontal stress. Future work as part of this US DOE-sponsored project includes additional validation through more detailed numerical modeling, laboratory testing, and field testing of the thermal breakout technology.

58 GEOSCIENCES↗