Discrete crack mechanics of disk compression for measurement of low fracture toughness
Not provided.
Materials Science↗
SEARCH · Search NASA
Results for “Discrete fractures”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding
A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.
Andrew C Bergan↗
Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding
A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.
Andrew Bergan↗
Relation of lineaments to sulfide deposits: Bald Eagle Mountain, Centre County, Pennsylvania
The author has identified the following significant results. Discrete areas of finely-fractured and brecciated sandstone float are present along the crest of Bald Mountain and are commonly sites of sulfide mineralization, as evidenced by the presence of barite and limonite gossans. The frequency distributions of the brecciated float as the negative binomial distribution supports the interpretation of a separate population of intensely fractured material. Such zones of concentrated breccia float have an average width of one kilometer with a range from 0.4 to 1.6 kilometers and were observed in a quarry face to have subvertical dips. Direct spatial correlation of the Landsat-derived lineaments to the fractured areas on the ridge is low; however, the mineralized and fracture zones are commonly assymetrical to the lineament positions. Such a systematic dislocation might result from an inherent bias in the float population or could be the product of the relative erosional resistance of the silicified material in the mineralized areas in relation to the erosionally weak material at the stream gaps.
Mcmurtry, G. J.↗
Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation
Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.
Dombard, A. J.↗
A hybrid meshfree discretization to improve the numerical performance of peridynamic models
Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.
42 ENGINEERING↗
In-Situ Visualization of a Growing Brittle Crack in Aluminum Oxynitride Using Synchrotron X-Rays and the Double-Cleavage Drilled Compression Geometry
Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.
Gorske, Sara F.↗
Surface and interfacial properties of carbon fibers
Differences in the adhesion of three carbon fibers (Hercules AS1 and AS4, and Hysol-Grafil XAS) to polycarbonate (PC) have been shown to correlate with the absorptivity of PC on the three fiber types. The absorptivity (energy of absorption) was determined using retention time liquid chromatography and the adhesion was measured using the single embedded filament tensile test. A correlation was also found between adhesion strength and the O/N surface element ratio using XPS analysis. The chemical details for these correlations have not yet been determined. A study of filament fracture statistics has been initiated using single and multiple embedded filament tensile tests. Filament fracture has been measured as a function of strain and for different interfiber distances. Preliminary results indicate that fiber fracture is a discontinuous function of increasing strain and may in fact occur at discrete strain intervals. Fiber-fiber interaction effects on fiber fracture have been found for interfiber distances of up to two to three fiber diameters.
Bascom, W. D.↗
Multiscale fatigue crack initiation in hierarchical additively manufactured alloys
Bioinspired hierarchical microstructures offer a route toward engineered fatigue resistance in additively manufactured alloys. However, it remains unclear how discrete structural constituents independently govern damage accumulation, particularly during the critical fatigue initiation regime where short cracks strongly interact with local microstructure. Here, we investigate multiscale fatigue initiation in a dual-phase, nanolamellar AlCoCrFeNi 2.1 high-entropy alloy. By comparing microscale specimens that isolate the nanolamellar structure against macroscale specimens containing the full melt-pool architecture, we identify size-dependent fatigue initiation mechanisms. We find that failure is dictated by nanolamellar interfaces at the microscale, whereas mesoscale melt pool boundaries serve to initiate fatigue at the macroscale. This mechanistic shift is accompanied by a transition from macroscale quasi-brittle failure to microscale plasticity-driven crack extension. Our results provide a physical framework for understanding how structural hierarchy governs the transition from discrete microstructural deformation to continuum fatigue fracture behavior, informing the design of damage-tolerant, additively manufactured alloys.
36 MATERIALS SCIENCE↗
Molecular dynamics simulation of propagating cracks
Steady state crack propagation is investigated numerically using a model consisting of 236 free atoms in two (010) planes of bcc alpha iron. The continuum region is modeled using the finite element method with 175 nodes and 288 elements. The model shows clear (010) plane fracture to the edge of the discrete region at moderate loads. Analysis of the results obtained indicates that models of this type can provide realistic simulation of steady state crack propagation.
Mullins, M.↗
A parametric analysis of oilfield design factors affecting the detectability and characterization of electrically conductive hydrofracks
Electrical responses in the vicinity of energized steel-cased well sources offer significant potential for monitoring induced fractures. However, the high complexity of well-fracture-host models spanning multiple length scales compels analysts to simplify their numerical models due to enormous computational costs. This consequently limits our understanding regarding monitoring capabilities and the limitations of electrical measurements on realistic hydraulically fracturing systems. In this paper, we use the hierarchical finite element approach to construct geoelectric models in which geometrically complex fractures and steel-cased wells are discretely represented in 3D conducting media without sacrificing the model realism and computation efficiency. We have discovered systematic numerical analyses of the electrical responses to evaluate the influences of borehole material conductivity and the source type as well as the effects of well geometry, conductivity contrast, source location, fracture growth, and fracture propagation. Furthermore, the numerical results indicate that the borehole material property has a strong control on the electrical potentials along the production and monitoring wells. The monopole source located at a steel-cased well results in a current density distribution that decays away from the source location throughout the well length, whereas the dipole source produces a current density that dominates mainly along the dipole length. Moreover, the conductivity contrast between the fractures and host does not change the overall pattern of the electrical potentials but varies its amplitude. The fracture models near different well systems indicate that the well geometry controls the entire distribution of potentials, while the characteristics of the voltage difference profiles along the wells before and after fracturing are insensitive to the well geometry and the well in which the source is located. Further, the hydraulic-fracturing models indicate that the voltage differences along the production well before and after fracturing have strong sensitivity to fracture growth and fracture set propagation.
58 GEOSCIENCES↗
Fabric controls on fracture surface roughness of an architected rock material
Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.
3D printed rock↗
DECOVALEX-2023: Task A Final Report
Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
A multi-resolution approach to hydraulic fracture simulation
Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.
58 GEOSCIENCES↗
The mechanics of brittle granular materials with coevolving grain size and shape
The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.
36 MATERIALS SCIENCE↗
500+ Million Years of Volcanism in Australe
The roughly circular collection of mare deposits in the "Mare Australe" region (38.9°S, 93°E) is ~900 km in diameter. The maria erupted into post-basin impact craters but are generally disconnected from each other. The distribution of maria could reflect what might have been the early stages of other basin filling mare eruptions, but which are preserved in Australe and some farside locations [1-3]. Topographic data fail to indicate a 900-km basin, but gravity and topographic data suggest a 600-km basin in the northern part of Australe [6]. The relationship between the basin and the wider-spread mare deposits is unclear, particularly outside the basin. We utilized images, gravity, topography, and composition to investigate the timing of volcanism in "Mare Australe" to interpret the evolution of magmatism in Australe and how it might relate to other mare sources. Humboldt crater (~3.6 Ga) makes an excellent stratigraphic marker that can be traced across much of the deposits and serves as a stratigraphic calibration for absolute model ages (AMAs) derived from crater size-frequency distributions (CSFDs). We found a range of AMAs from ~3.3 to 3.8 Ga in north Australe, similar to previous studies [4], and there is also evidence for multiple deposits of older cryptomaria. The age range, however, is narrower than for maria in other basins (e.g., Imbrium and South Pole-Aitken [4-5]), suggesting different magmatism in Australe. Most of the mare deposits are >3.5 Ga; younger deposits are co-located in the north basin. The basin may have provided additional pathways for rising magma to the surface over time. We also found evidence of a several interconnected deposits within the north basin. Other, smaller deposits, however, clearly show distinct age differences between units within a single topographic depression (e.g., Lyot crater) and likely represent discrete eruptions. Additionally, there are no floor fractured craters other than Humboldt [6], suggesting that magma generally did not stall out near the surface, also consistent with gravity and (lack of) magnetic anomalies (in contrast to Marius Hills and South Pole-Aitken basin [7-8]).
Julie Dawn Stopar↗
Analysis of fatigue crack propagation
The correlation between fatigue crack propagation and stress intensity factor is analyzed. When determining fatigue crack propagation rate, a crack increment, delta a, and its corresponding increment in load cycles, delta N, are measured. Fatigue crack propagation must be caused by a shear and/or a normal separation mode. Both of these two processes are discrete if one looks at the atomic level. If the average deformation and fracture properties over the crack increments, delta a, can be considered as homogeneous, if the characteristic discrete lengths of sigma a, if the plastic zone size is small, and if a plate is thick enough to insure a plane strain case, da/dN is proportional to delta K squared. Any deviation of empirical data from this relation must be caused by the fact that one or more of these conditions are not satisfied. The effects of plate thickness and material inhomogeneity are discussed in detail. A shear separation mode of fatigue crack propagation is described and is used to illustrate the effects of material inhomogeneity.
Liu, H. W.↗