Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture simulation”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Interaction Between Proppant Compaction and Single-/Multiphase Flows in a Hydraulic Fracture

Understanding proppant-transport and -deposition patterns in a hydraulic fracture is vital for effective and economical production of petroleum hydrocarbons. In this research, a numerical-modeling approach, combining the discrete-element method (DEM) with single-/multiphase lattice Boltzmann (LB) simulation, was adopted to advance understanding of the interaction between reservoir depletion, proppant-particle compaction, and single-/multiphase flows in a hydraulic fracture. DEM was used to simulate effective-stress increase and the resultant proppant-particle movement and rearrangement during the process of reservoir depletion caused by hydrocarbon production. Simulated pore structure of the proppant pack was extracted and used as boundary conditions in the LB simulation to calculate the time-dependent permeability of the proppant pack. We first validated the DEM-LB numerical work flow, and the simulated proppant-pack permeabilities as functions of effective stress were in good agreement with laboratory measurements. Furthermore, three proppant packs with the same average particle diameter but different diameter distributions were generated to study the role of proppant-size heterogeneity (variation in particle diameter). Specifically, we used the coefficient of variation (COV) of diameter, defined as the ratio of standard deviation of diameter to mean diameter, to characterize the heterogeneity of particle size. We obtained proppant-pack porosity, permeability, and fracture-width reduction (closure distance) as functions of effective stress. Under the same effective stress, a proppant pack with a higher diameter COV had lower porosity and permeability and larger fracture-width reduction. This was because the high diameter COV gave rise to a wider diameter distribution of proppant particles; smaller particles were compressed into the pore space between larger particles with increasing stress, leading to larger closure distance and lower porosity and permeability. With multiphase LB simulation, relative permeability curves were obtained, which are critical for larger-scale reservoir simulations under various effective stresses. The relative permeability of the oil phase was more sensitive to changes in geometry and stress, compared with the wetting phase. This was because the oil phase occupied larger pores; compaction of the proppant pack affected the structure of the pores, because the pores were farther from the grain contacts and thus more susceptible to collapse. It is also interesting to note that when effective stress increased continuously, oil relative permeability increased first and then decreased. This nonmonotonic behavior was the result of the nonmonotonic development of pore structure and oil connectivity under increasing stress.

Fan, Ming↗
DOE OSTI2022

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
DOE OSTI2023

Economic and operational investigation of CO 2 sequestration through enhanced oil recovery in unconventional reservoirs in Colorado, USA

The ongoing CCUS commercial projects are highly relied on the support of government incentives due to massive capital investment. Here, this study analyzes the economics of carbon capture utilization and sequestration (CCUS) projects, shows a state-wide CCUS deployment exercise, followed by simulation results of enhanced oil recovery (EOR) based CO 2 storage in unconventional reservoirs. The comprehensive economic analysis of capture, transportation, sequestration costs, enhanced 45Q tax credits, and EOR revenue implies the practicality of CO 2 -EOR to offset the high CCUS costs. With the economics understanding, we study the top CO 2 sources, existing CO 2 pipelines, and sequestration sinks in the state of Colorado, USA. This paper next presents results from EOR simulation in one section of the unconventional Denver-Julesburg (DJ) Basin Niobrara and Codell reservoirs. The simulation model is based on a geological static model, incorporated with hydraulic fracture stimulation, history matched to production, and calibrated to the microseismic and time-lapse surface seismic data. The CO 2 -EOR simulation results show that oil production can be increased and more CO 2 stored with: a longer primary production period; the presence of a shut-in period; higher injection rates; and multi-well injectors. The modeling results show that about 7–10 Mscf of CO 2 will be stored when recovering 1 stb of EOR oil. By adding the enhanced oil revenue and the carbon credits together, it is estimated that the most economic case can generate $\$13$ MM when oil price is assumed to be $\$80$/stb, and the EOR oil revenue is 3.4 times greater than that generated from 45Q incentives. It corresponds to the scenario that a five-year primary production is followed by CO 2 injection into four wells with the sequence of injection (4 MMscf/day for 6 months), shut-in (6 months) and production (12 months). The best practices in this study will provide valuable insights for similar CCUS projects in other unconventional fields. Furthermore, this study defines a term named “Carbon Neutrality Index (CNI)” by comparing the amount of CO 2 stored with that burned by EOR oil. The CNI value of 0 indicates the enhanced oil is carbon neutral; a negative CNI value implies there is a net reduction in carbon emission. The 4-year huff-n-puff (HnP) simulation leads to a positive CNI value, indicating that the EOR oil generated in this process is not carbon neutral yet.

03 NATURAL GAS↗
DOE OSTI2022

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗
DOE OSTI2024

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗
DOE OSTI2014

Deformation behavior of sub-micron and micron sized alumina particles in compression

The ability to integrate ceramics with other materials has been limited due to high temperature (>800°C) ceramic processing. Recently, researchers demonstrated a novel process, aerosol deposition (AD), to fabricate ceramic films at room temperature (RT). In this process, sub-micron sized ceramic particles are accelerated by pressurized gas, impacted on the substrate, plastically deformed, and form a dense film under vacuum. This AD process eliminates high temperature processing thereby enabling new coatings and device integration, in which ceramics can be deposited on metals, plastics, and glass. However, knowledge in fundamental mechanisms for ceramic particles to deform and form a dense ceramic film is still needed and is essential in advancing this novel RT technology. In this work, a combination of experimentation and atomistic simulation was used to determine the deformation behavior of sub-micron sized ceramic particles; this is the first fundamental step needed to explain coating formation in the AD process. High purity, single crystal, alpha alumina particles with nominal sizes of 0.3 µm and 3.0 µm were examined. Particle characterization, using transmission electron microscopy (TEM), showed that the 0.3 µm particles were relatively defect-free single crystals whereas 3.0 µm particles were highly defective single crystals or particles contained low angle grain boundaries. Sub-micron sized Al 2 O 3 particles exhibited ductile failure in compression. In situ compression experiments showed 0.3µm particles deformed plastically, fractured, and became polycrystalline. Moreover, dislocation activity was observed within these particles during compression. These sub-micron sized Al 2 O 3 particles exhibited large accumulated strain (2-3 times those of micron-sized particles) before first fracture. In agreement with the findings from experimentation, atomistic simulations of nano-Al 2 O 3 particles showed dislocation slip and significant plastic deformation during compression. On the other hand, the micron sized Al 2 O 3 particles exhibited brittle fracture in compression. In situ compression experiments showed 3µm Al 2 O 3 particles fractured into pieces without observable plastic deformation in compression. Particle deformation behaviors will be used to inform Al 2 O 3 coating deposition parameters and particle-particle bonding in the consolidated Al 2 O 3 coatings.

Sarobol, Pylin↗
NASA NTRS1985

Elevated temperature crack growth

The purpose is to determine the ability of currently available P-I integrals to correlate fatigue crack propagation under conditions that simulate the turbojet engine combustor liner environment. The utility of advanced fracture mechanics measurements will also be evaluated during the course of the program. To date, an appropriate specimen design, a crack displacement measurement method, and boundary condition simulation in the computational model of the specimen were achieved. Alloy 718 was selected as an analog material based on its ability to simulate high temperature behavior at lower temperatures. Tensile and cyclic tests were run at several strain rates so that an appropriate constitutive model could be developed. Suitable P-I integrals were programmed into a finite element post-processor for eventual comparison with experimental data.

Malik, S. N.↗
NASA NTRS1996

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗
DOE OSTI2016

Role of Defects in Swelling and Creep of Irradiated SiC

Silicon carbide is a promising cladding material because of its high strength and relatively good corrosion resistance. However, SiC is brittle and therefore SiC-based components need to be carefully designed to avoid cracking and failure by fracture. In design of SiC-based composites for nuclear reactor applications it is essential to take into account how mechanical properties are affected by radiation and temperature, or in other words, what strains and stresses develop in this material due to environmental conditions. While thermal strains in SiC can be predicted using classical theories, radiation-induced strains are much less understood. In particular, it is critical to correctly account for radiation swelling and radiation creep, which contribute significantly to dimensional instability of SiC under radiation. Swelling typically increases logarithmically with radiation dose and saturates at relatively low doses (damage levels of a few dpa). Consequently, swelling-induced stresses are likely to develop within a few months of operation of a reactor. Radiation-induced volume swelling in SiC can be as high as 2%, which is significantly higher than the cracking strain of 0.1% in SiC. Swelling-induced strains will lead to enormous stresses and fracture, unless these stresses can be relaxed via some other mechanism. An effective way to achieve stress relaxation is via radiation creep. Although it has been hypothesized that both radiation swelling and radiation creep are driven by formation of defect clusters, existing models for swelling and creep in SiC are limited by the lack of understanding of specific defects that form due to radiation in the range of temperatures relevant to fuel cladding in light water reactors (LWRs) (<1000°C). For example, defects that can be detected with traditional transmission electron microscopy (TEM) techniques account only for 10-45% of the swelling measured in irradiated SiC. Here, we have undertaken an integrated experimental and modeling effort to discover the previously invisible defects in irradiated SiC and to determine the contributions of these defects to radiation swelling. Knowledge of the most stable defect structures and the rate controlling processes during defect evolution is essential for development of predictive models for swelling and creep as a function of temperature and radiation dose. This research has been enabled by state-of-the-art imaging techniques, such as the aberration corrected scanning transmission electron microscopy (STEM) (FEI TITAN) closely coupled with multi-scale models of stable defect clusters and their evolution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Study on grain boundary character and strain distribution of intergranular cracking in the CGHAZ of T23 steel

Intergranular reheat cracking in the coarse-grained heat-affected zone of T23 steel was produced by strain to fracture tests on a Gleeble 3500 thermal–mechanical simulator. Then the grain boundary character, as well as the strain distribution after reheat crack propagation, was studied by electron backscatter diffraction technique. The results showed that incoherent Σ3 boundaries were seldom found on the prior austenite grain boundaries. Therefore, only the type of random high-angle boundaries played a crucial role in the intergranular cracking. Microstructurally cavities and small cracks were preferentially initiated from high-angle grain boundaries. Low-angle grain boundaries and high-angle ones with misorientation angles less than 15° were more resistant to the cracking. More importantly, the fraction of high-angle grain boundaries increased with the plastic strain induced by both temperature gradient and stress in the coarse-grained heat-affected zone, which contributed to the crack initiation and propagation. Furthermore, the strain distributions in the vicinity of cavities and cracks revealed the accommodation processes of plastic deformation during stress relaxation. It also reflected the strength differences between grain interior and grain boundary at different heat-treated temperatures, which had a large influence on the cracking mechanism. - Highlights: • The coincidence site lattice boundaries play little role in the reheat cracking. • Cavity and crack occur at high-angle grain boundaries rather than low-angle ones. • The strain leads low-angle grain boundaries to transform to high-angle ones. • Strain distribution differs for cavity and crack zones at different temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Influence of Porosity and Pore Distributions on Strength Properties of Porous Alumina

In this work, strength properties of porous alumina with three different porosities as well as dense one were investigated to clarify porosity effects on strength properties. It was revealed that the strength had a general tendency to decrease drastically as increasing porosity. Spatial and size distributions of pores were characterized by observing specimen surfaces through a laser scanning microscope. Based on the observed pore characteristics, a procedure to estimate porosity effects on strength was newly developed in the framework of fracture mechanics, where a stress intensity factor-based criterion was applied for fracture of porous ceramics. In the procedure, it is assumed that porous ceramics have the same characteristics of crack distribution as those of pore distribution, and pores are presumed to be surrounded by virtual cracks. Monte Carlo simulations by using the procedure were carried out to estimate strength for four respective materials. Simulated strength coincided with experimentally observed one within a scatter band of factor of 2{sup 1/2}. Consequently, the developed procedure was confirmed to be efficient in evaluating effects of bulk porosity and pore distributions on strength of porous alumina.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗
DOE OSTI2014

Edge Fracture Prediction ofTraditional and Advanced Trimming Processes for AA6111-T4 Sheets

This work examines the traditional and advanced trimming of AA6111-T4 aluminum sheets with finite element simulations. The Rice-Tracy damage model is used for the simulation with damage parameters estimated from experimental observation of grain aspect ratio near the fracture surface of trimmed parts. Fine meshes at the shearing zone, adaptive meshing, and adaptive contact techniques are used to accurately capture the contact interactions between the sharp corner of the trimming tools and the blank to be trimmed. To the knowledge of the authors, these are the first trimming simulations that can predict the effects of shearing clearance on burr heights with quantitative accuracy for AA6111-T4 aluminum sheets. In addition, the models have also accurately reproduced the crack initiation site as well as burr and sliver formation mechanisms observed experimentally.

Adaptive Contact.↗
DOE OSTI2017

Deformed Materials: Towards a Theory of Materials Morphology Dynamics

This grant supported work on the response of crystals to external stress. Our primary work described how disordered structural materials break in two (statistical models of fracture in disordered materials), studied models of deformation bursts (avalanches) that mediate deformation on the microscale, and developed continuum dislocation dynamics models for plastic deformation (as when scooping ice cream bends a spoon, Fig. 9). Glass is brittle -- it breaks with almost atomically smooth fracture surfaces. Many metals are ductile -- when they break, the fracture surface is locally sheared and stretched, and it is this damage that makes them hard to break. Bone and seashells are made of brittle material, but they are strong because they are disordered -- lots of little cracks form as they are sheared and near the fracture surface, diluting the external force. We have studied materials like bone and seashells using simulations, mathematical tools, and statistical mechanics models from physics. In particular, we studied the extreme values of fracture strengths (how likely will a beam in a bridge break far below its design strength), and found that the traditional engineering tools could be improved greatly. We also studied fascinating crackling-noise precursors -- systems which formed microcracks of a broad range of sizes before they broke. Ductile metals under stress undergo irreversible plastic deformation -- the planes of atoms must slide across one another (through the motion of dislocations) to change the overall shape in response to the external force. Microscopically, the dislocations in crystals move in bursts of a broad range of sizes (termed 'avalanches' in the statistical mechanics community, whose motion is deemed 'crackling noise'). In this grant period, we resolved a longstanding mystery about the average shape of avalanches of fixed duration (using tools related to an emergent scale invariance), we developed the fundamental theory describing the shapes of avalanches and how they are affected by the edges of the microscope viewing window, we found that slow creep of dislocations can trigger an oscillating response explaining recent experiments, we explained avalanches under external voltage, and we have studied how avalanches in experiments on the microscale relate to deformation of large samples. Inside the crystals forming the metal, the dislocations arrange into mysterious cellular structures, usually ignored in theories of plasticity. Writing a natural continuum theory for dislocation dynamics, we found that it spontaneously formed walls -- much like models of traffic jams and sonic booms. These walls formed rather realistic cellular structures, which we examined in great detail -- our walls formed fractal structures with fascinating scaling properties, related to those found in turbulent fluids. We found, however, that the numerical and mathematical tools available to solve our equations were not flexible enough to incorporate materials-specific information, and our models did not show the dislocation avalanches seen experimentally. In the last year of this grant, we wrote an invited review article, explaining how plastic flow in metals shares features with other stressed materials, and how tools of statistical physics used in these other systems might be crucial for understanding plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Asteroid diversion considerations and comparisons of diversion techniques

The threat of asteroid impacts on Earth poses a low-probability but high consequence risk, with possible outcomes ranging from regional to global catastrophe. However, unique amongst such global threats we have the capability of averting such disasters. Diversion approaches by either kinetic impactor or nuclear energy deposition are the two most practical technologies for mitigating hazardous near Earth asteroids. One of the greatest challenges in understanding our options is the uncertain response of asteroids to such impulsive techniques, due both to our lack of knowledge of the composition and structure of these objects as well as their highly varied nature. Predicting whether we will simply divert or break up a given object is a crucial: the weak self-gravity and inferred weak structure of typical asteroids present the strong possibility the body will fragment for modest impulses. Predictive modeling of failure and fragmentation is one important tool for such studies. In this paper we apply advances in modeling failure and fracture using Adaptive Smoothed Particle Hydrodynamics (ASPH) to understand mega-cratering on asteroids as a validation exercise, and show examples of diverting the near Earth asteroid Bennu using both a kinetic impactor and ablative blow-off due to nuclear energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗