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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Grizzly Usage and Theory Manual

Grizzly is a multiphysics simulation code for characterizing the behavior of nuclear power plant (NPP) structures, systems and components (SSCs) subjected to a variety of age-related aging mechanisms. Grizzly simulates both the progression of aging processes, as well as the capacity of aged components to safely perform. This initial beta release of Grizzly includes capabilities for engineering-scale thermo-mechanical analysis of reactor pressure vessels (RPVs). Grizzly will ultimately include capabilities for a wide range of components and materials. Grizzly is in a state of constant development, and future releases will broaden the capabilities of this code for RPV analysis, as well as expand it to address degradation in other critical NPP components.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1976

Application of a novel finite difference method to dynamic crack problems

A versatile finite difference method (HEMP and HEMP 3D computer programs) was developed originally for solving dynamic problems in continuum mechanics. It was extended to analyze the stress field around cracks in a solid with finite geometry subjected to dynamic loads and to simulate numerically the dynamic fracture phenomena with success. This method is an explicit finite difference method applied to the Lagrangian formulation of the equations of continuum mechanics in two and three space dimensions and time. The calculational grid moves with the material and in this way it gives a more detailed description of the physics of the problem than the Eulerian formulation.

Chen, Y. M.↗
NASA NTRS1997

Progressive Fracture of Fiber Composite Build-Up Structures

Damage progression and fracture of built-up composite structures is evaluated by using computational simulation. The objective is to examine the behavior and response of a stiffened composite (0/ +/- 45/90)(sub s6) laminate panel by simulating the damage initiation, growth, accumulation, progression and propagation to structural collapse. An integrated computer code, CODSTRAN, was augmented for the simulation of the progressive damage and fracture of built-up composite structures under mechanical loading. Results show that damage initiation and progression have significant effect on the structural response. Influence of the type of loading is investigated on the damage initiation, propagation and final fracture of the build-up composite panel.

Gotsis, Pascal K.↗
NASA NTRS1997

Progressive Fracture of Fiber Composite Build-Up Structures

Damage progression and fracture of built-up composite structures is evaluated by using computational simulation. The objective is to examine the behavior and response of a stiffened composite (0 +/-45/90)(sub s6) laminate panel by simulating the damage initiation, growth, accumulation, progression and propagation to structural collapse. An integrated computer code CODSTRAN was augmented for the simulation of the progressive damage and fracture of built-up composite structures under mechanical loading. Results show that damage initiation and progression to have significant effect on the structural response. Influence of the type of loading is investigated on the damage initiation, propagation and final fracture of the build-up composite panel.

Minnetyan, Levon↗
NASA NTRS2005

Thermal and Environmental Barrier Coatings for Advanced Turbine Engine Applications

Ceramic thermal and environmental barrier coatings (T/EBCs) will play a crucial role in advanced gas turbine engine systems because of their ability to significantly increase engine operating temperatures and reduce cooling requirements, thus help achieve engine low emission and high efficiency goals. Advanced T/EBCs are being developed for the low emission SiC/SiC ceramic matrix composite (CMC) combustor applications by extending the CMC liner and vane temperature capability to 1650 C (3000 F) in oxidizing and water vapor containing combustion environments. Low conductivity thermal barrier coatings (TBCs) are also being developed for metallic turbine airfoil and combustor applications, providing the component temperature capability up to 1650 C (3000 F). In this paper, ceramic coating development considerations and requirements for both the ceramic and metallic components will be described for engine high temperature and high-heat-flux applications. The underlying coating failure mechanisms and life prediction approaches will be discussed based on the simulated engine tests and fracture mechanics modeling results.

Zhu, Dong-Ming↗
DOE OSTI2020

Simulating current flow through a well casing and an induced fracture

A system includes a receiver component that receives a model of a well system, the model comprising a representation of an electrically energized well casing and/or a fracture as a transmission line that leaks electric current as the current traverses the well casing and/or fracture. The receiver component also receives a value that indicates an amount of the current that is applied to the well casing and a location of a source of the current on the well casing. The system further includes an electromagnetic field calculator component that calculates an estimated electromagnetic field at at least one location on the surface of the earth based at least in part upon the representation of the electrically energized well casing and/or fracture, the value that indicates the amount of current that is applied to the well casing, and the location of the source of the current on the well casing.

02 PETROLEUM↗
NASA NTRS2007

Progressive Fracture of [0/90/ + or - Theta]s Composite Structure Under Uniform Pressure Load

S-Glass/epoxy [0/90/plus or minus theta]s for theta =45 deg., 60 deg., and 75 deg. laminated fiber-reinforced composite stiffened plate was simulated to investigated for damage and fracture progression under uniform pressure. An integrated computer code was augmented for the simulation of the damage initiation, growth, accumulation, and propagation to fracture and to structural collapse. Results show in detail the damage progression sequence and structural fracture resistance during different degradation stages. Damage through the thickness of the laminate initiated first at [0/90/plus or minus 45]s at 15.168 MPa (2200 psi), followed by [0/90/plus or minus 60]s at 16.96 MPa (2460 psi) and finally by [0/90/plus or minus 75]s at 19.3 MPa (2800 psi). After damage initiation happened the cracks propagate rapidly to structural fracture.

Gotsis, Pascalis K.↗
DOE OSTI2019

Characterizing the impact of particle behavior at fracture intersections in three-dimensional discrete fracture networks

We characterize the influence of different intersection mixing rules for particle tracking simulations on transport properties through three-dimensional discrete fracture networks. It is too computationally burdensome to explicitly resolve all fluid dynamics within a large three-dimensional fracture network. In discrete fracture network (DFN) models, mass transport at fracture intersections is modeled as a subgrid scale process based on a local Péclet number. The two most common mass transfer mixing rules are (1) complete mixing, where diffusion dominates mass transfer, and (2) streamline routing, where mass follows pathlines through an intersection. Although it is accepted that mixing rules impact local mass transfer through single intersections, the effect of the mixing rule on transport at the fracture network scale is still unresolved. Through the use of explicit particle tracking simulations, we study transport through a quasi-two-dimensional lattice network and a three-dimensional network whose fracture radii follow a truncated power-law distribution. We find that the impact of the mixing rule is a function of the initial particle injection condition, the heterogeneity of the velocity field, and the geometry of the network. Furthermore, our particle tracking simulations show that the mixing rule can particularly impact concentrations on secondary flow pathways. As a result, we relate these local differences in concentration to reactive transport and show that streamline routing increases the average mixing rate in DFN simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Strain-dependent partial slip on rock fractures under seismic-frequency torsion: Seismic-Frequency Fracture Partial Slip

Measurements of nonlinear modulus and attenuation of fractures provide the opportunity to probe their mechanical state. We have adapted a low-frequency torsional apparatus to explore the seismic signature of fractures under low normal stress, simulating low effective stress environments such as shallow or high pore pressure reservoirs. We report strain-dependent modulus and attenuation for fractured samples of Duperow dolomite (a carbon sequestration target reservoir in Montana), Blue Canyon Dome rhyolite (a geothermal analog reservoir in New Mexico), and Montello granite (a deep basement disposal analog from Wisconsin). We use a simple single effective asperity partial slip model to fit our measured stress-strain curves, and solve for the friction coefficient, contact radius, and full slip condition. These observations have the potential to develop into new field techniques for measuring differences in frictional properties during reservoir engineering manipulations and estimate the stress conditions where reservoir fractures and faults begin to fully slip.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗
DOE OSTI2022

Grain-boundary fracture mechanisms in Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolytes: When phase transformation acts as a temperature-dependent toughening mechanism

Garnet-type, solid electrolytes, such as Li 7 La 3 Zr 2 O 12 (LLZO), are a promising alternative to liquid electrolytes for lithium-metal batteries. However, such solid-electrolyte materials frequently exhibit undesirable lithium (Li) metal plating and fracture along grain boundaries. In this work, we employ atomistic simulations to investigate the mechanisms and key fracture properties associated with intergranular fracture along one such boundary. Our results show that, in the case of a Σ5 (310) grain boundary, this boundary exhibits brittle fracture behavior, i.e. the absence of dislocation activity ahead of the propagating crack tip, accompanied with a decrease in work of separation, peak stress, and maximum stress intensity factor as the temperature increases from 300 K to 1500 K. As the crack propagates, we predict two temperature-dependent Li clustering regimes. For temperatures at or below 900 K, Li tends to cluster in the bulk region away from the crack plane driven by a void-coalescence mechanism concomitant a simultaneous cubic-to-tetragonal phase transition. The tetragonalization of LLZO in this temperature regime acts as an emerging toughening mechanism. At higher temperatures, this phase transition mechanism is suppressed leading to a more uniform distribution of Li throughout the grain-boundary system and lower fracture properties as compared to lower temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

58 GEOSCIENCES↗
NASA NTRS2015

Probabilistic Risk Assessment for Astronaut Post Flight Bone Fracture

Introduction: Space flight potentially reduces the loading that bone can resist before fracture. This reduction in bone integrity may result from a combination of factors, the most common reported as reduction in astronaut BMD. Although evaluating the condition of bones continues to be a critical aspect of understanding space flight fracture risk, defining the loading regime, whether on earth, in microgravity, or in reduced gravity on a planetary surface, remains a significant component of estimating the fracture risks to astronauts. This presentation summarizes the concepts, development, and application of NASA's Bone Fracture Risk Module (BFxRM) to understanding pre-, post, and in mission astronaut bone fracture risk. The overview includes an assessment of contributing factors utilized in the BFxRM and illustrates how new information, such as biomechanics of space suit design or better understanding of post flight activities may influence astronaut fracture risk. Opportunities for the bone mineral research community to contribute to future model development are also discussed. Methods: To investigate the conditions in which spaceflight induced changes to bone plays a critical role in post-flight fracture probability, we implement a modified version of the NASA Bone Fracture Risk Model (BFxRM). Modifications included incorporation of variations in physiological characteristics, post-flight recovery rate, and variations in lateral fall conditions within the probabilistic simulation parameter space. The modeled fracture probability estimates for different loading scenarios at preflight and at 0 and 365 days post-flight time periods are compared. Results: For simple lateral side falls, mean post-flight fracture probability is elevated over mean preflight fracture probability due to spaceflight induced BMD loss and is not fully recovered at 365 days post-flight. In the case of more energetic falls, such as from elevated heights or with the addition of lateral movement, the contribution of space flight quality changes is much less clear, indicating more granular assessments, such as Finite Element modeling, may be needed to further assess the risks in these scenarios.

Risk Analysis↗
NASA NTRS1978

Equilibrium and kinetic aspects of brittle fracture

The lattice statics method for simulating a brittle crack will be employed to determine: (1) the energy-based criterion for crack stability, (2) at a given stress intensity in the stable region, the kinetic barrier for crack motion, and (3) the assumptions that break down in the Griffith theory when generalized to the lattice model. A simple model is used in the analysis, but this model is general enough to provide a conceptual framework for more realistic models.

Esterling, D. M.↗
DOE OSTI2025

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗
DOE OSTI2017

Modeling of internal mechanical failure of all-solid-state batteries during electrochemical cycling, and implications for battery design

This study is the first quantitative analysis of mechanical reliability of all-solid state batteries. Mechanical degradation of the solid electrolyte (SE) is caused by intercalation-induced expansion of the electrode particles, within the constrains of a dense microstructure. A coupled electro-chemo-mechanical model was implemented to quantify the material properties that cause an SE to fracture. The treatment of microstructural details is essential to the understanding of stress-localization phenomena and fracture. A cohesive zone model is employed to simulate the evolution of damage. In the numerical tests, fracture is prevented when electrode-particle's expansion is lower than 7.5% (typical for most Li-intercalating compounds) and the solid-electrolyte's fracture energy higher than G c = 4 J m -2 . Perhaps counter-intuitively, the analyses show that compliant solid electrolytes (with Young's modulus in the order of E SE = 15 GPa) are more prone to micro-cracking. This result, captured by our non-linear kinematics model, contradicts the speculation that sulfide SEs are more suitable for the design of bulk-type batteries than oxide SEs. Mechanical degradation is linked to the battery power-density. Fracture in solid Li-ion conductors represents a barrier for Li transport, and accelerates the decay of rate performance.

25 ENERGY STORAGE↗