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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
NASA NTRS2007

Prediction of Size Effects in Notched Laminates Using Continuum Damage Mechanics

This paper examines the use of a continuum damage model to predict strength and size effects in notched carbon-epoxy laminates. The effects of size and the development of a fracture process zone before final failure are identified in an experimental program. The continuum damage model is described and the resulting predictions of size effects are compared with alternative approaches: the point stress and the inherent flaw models, the Linear-Elastic Fracture Mechanics approach, and the strength of materials approach. The results indicate that the continuum damage model is the most accurate technique to predict size effects in composites. Furthermore, the continuum damage model does not require any calibration and it is applicable to general geometries and boundary conditions.

Camanho, D. P.↗
DOE OSTI2021

Directly Probing the Fracture Behavior of Ultrathin Polymeric Films

Understanding fracture mechanics of ultrathin polymeric films is crucial for modern technologies, including semiconductor and coating industries. However, up to now, the fracture behavior of sub-100 nm polymeric thin films is rarely explored due to challenges in handling samples and limited testing methods available. In this work, we report a new testing methodology that can not only visualize the evolution of the local stress distribution through wrinkling patterns and crack propagation during the deformation of ultrathin films but also directly measure their fracture energies. Using ultrathin polystyrene films as a model system, we both experimentally and computationally investigate the effect of the film thickness and molecular weight on their fracture behavior, both of which show a ductile-to-brittle transition. Furthermore, we demonstrate the broad applicability of this testing method in semicrystalline semiconducting polymers. We anticipate our methodology described here could provide new ways of studying the fracture behavior of ultrathin films under confinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Skeletal Structural Consequences of Reduced Gravity Environments

The overall goal of this project is to provide structurally meaningful data on bone loss after exposure to reduced gravity environments so that more precise estimates of fracture risk and the effectiveness of countermeasures in reducing fracture risk can be developed. The project has three major components: (1) measure structural changes in the limb bones of rats subjected to complete and partial nonweightbearing, with and without treatment with ibandronate and periodic full weightbearing; (2) measure structural changes in the limb bones of human bedrest subjects, with and without treatment with alendronate and resistive exercise, and Russian cosmonauts flying on the Mir Space Station; and (3) validate and extend the 2-dimensional structural analyses currently possible in the second project component (bedrest and Mir subjects) using 3-dimensional finite element modeling techniques, and determine actual fracture-producing loads on earth and in space.

Ruff, Christropher B.↗
NASA NTRS2006

Computer Simulation of Fracture in Aerogels

Aerogels are of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While the gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. In this work, we investigate the strength and fracture behavior of silica aerogels using a molecular statics-based computer simulation technique. The gels' structure is simulated via a Diffusion Limited Cluster Aggregation (DLCA) algorithm, which produces fractal structures representing experimentally observed aggregates of so-called secondary particles, themselves composed of amorphous silica primary particles an order of magnitude smaller. We have performed multi-length-scale simulations of fracture in silica aerogels, in which the interaction b e e n two secondary particles is assumed to be described by a Morse pair potential parameterized such that the potential range is much smaller than the secondary particle size. These Morse parameters are obtained by atomistic simulation of models of the experimentally-observed amorphous silica "bridges," with the fracture behavior of these bridges modeled via molecular statics using a Morse/Coulomb potential for silica. We consider the energetics of the fracture, and compare qualitative features of low-and high-density gel fracture.

Good, Brian S.↗
DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗
NASA NTRS1994

Establishment of redox conditions during planetary collisions as an origin of chondrites

Collisions between a 'cometlike' body (mixtures of chondritic materials and ice) and a slightly differentiated body were proposed for shock origin of ordinary chondrites. In this model, chondrules were formed with shock melting, and matrices were formed both by fracturing of materials and by recondensation of evaporated materials. This model can explain different redox conditions of chondrite formations by ice evaporation. Although this model was originally proposed for ordinary chondrites, we assume here that the model can be extended to chondrite formation in general. In this paper, redox conditions during chondrite formation by collisions will be discussed in the light of phase diagrams for solid-vapor equilibria.

Tsuchiyama, A.↗
NASA NTRS2007

A New Sea Ice Model: Pancakes to Metrics

We are developing, testing, and validating a new sea ice dynamics model that treats the ice cover as an elastic/decohesive material in the permanent pack and includes the correct frazil/pancake behavior in the marginal zone. Two salient features of present ice dynamics models are that they do not: 1) reproduce the oriented fracture patterns of openings and closings in the pack ice, and 2) accurately model the effects of frazil/pancake ice formation in the ice margin. These poorly modeled areas account for a substantial portion of the ice growth, turbulent heat flux to the atmosphere, salt flux to the ocean, and energy dissipation due to slippage, ridging, and rafting, in the Arctic. Existing sea ice models have shown limited success in predicting the degree to which a lead will open for prescribed or observed forcing conditions. An important aspect of the new model we are developing is that the existence of cracks, along with their orientation, opening, and closing, is predicted. To put this effort in perspective a short history of ice dynamics modeling and data collection is presented. The RGPS data set is used to validate the model. As part of the testing and validation of the model, we are working on a new metric for comparing linear features (leads and ridges) in the data and model to be used in data assimilation for this model. The model framework is presented as well as some results showing the creation and development of leads in a simulation of ice dynamics in the Beaufort Sea. Other presentations by the authors will show other results from this effort.

Coon, M.↗
DOE OSTI2021

Representative Modules for Accelerated Thermal Cycling and Static Load Testing

In this work, we explore the influence of module size on the rate of interconnect solder bond thermomechanical fatigue (TMF) damage and the probability of cell fracture. For the solder bond TMF damage evaluation, structural mechanics models of crystalline silicon PV models are created to solve with the Finite Element Method. For the probability of cell fracture evaluation, 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. Results conclusively demonstrate that the rate of solder bond TMF damage is independent of module size, interconnect location across the cell and cell location across the module and that smaller, representative, modules must be loaded to a much higher level than their parent full-sized modules to achieve an equivalent driving force for cell fracture.

accelerated testing↗
NASA NTRS1994

Investigation of the effects of manufacturing variations and materials on fatigue crack detection methods in gear teeth

The fatigue life of a gear tooth can be thought of as the sum of the number of cycles required to initiate a crack, N(sub i), plus the number of cycles required to propagate the crack to such a length that fracture occurs, N(sub p). The factors that govern crack initiation are thought to be related to localized stress or strain at a point, while propagation of a fatigue crack is a function of the crack tip parameters such as crack shape, stress state, and stress intensity factor. During a test there is no clear transition between initiation and propagation. The mechanisms of initiation and propagation are quite different and modeling them separately produces a higher degree of accuracy, but then the question that continually arises is 'what is a crack?' The total life prediction in a fracture mechanics model presently hinges on the assumption of an initial crack length, and this length can significantly affect the total life prediction. The size of the initial crack is generally taken to be in the range of 0.01 in. to 0.2 in. Several researchers have used various techniques to determine the beginning of the crack propagation stage. Barhorst showed the relationship between dynamic stiffness changes and crack propagation. Acoustic emissions, which are stress waves produced by the sudden movement of stressed materials, have also been successfully used to monitor the growth of cracks in tensile and fatigue specimens. The purpose of this research is to determine whether acoustic emissions can be used to define the beginning of crack propagation in a gear using a single-tooth bending fatigue test.

Wheitner, Jeffrey A.↗
DOE OSTI2020

Visco-inertial gas flow through wellbore cement fractures

Understanding the nature of fluid flow through fractured wellbore cement is fundamental for evaluating the leakage potential and risk assessments of leaky wellbores. In this study, the conditions that require considering visco-inertial flow for describing the gas flow through wellbore cement fractures were investigated. Nitrogen gas flow tests were conducted on fractured cement samples under varying pressure conditions and flow rates, covering both viscous and visco-inertial flow regimes. The data substantially deviated from Darcy's law at higher flowrates and were well-fit to Forchheimer's equation for visco-inertial flow. The inertial coefficient and critical Reynolds number were expressed as a function of the hydraulic aperture. Finally, the empirical function obtained from the experiments was used as an input to numerical simulations which showed the significant role of visco-inertial flow in wellhead pressure build-up and leakage rates, and demonstrated the importance of visco-inertial flow when modeling gas flow through wellbore cement fractures.

03 NATURAL GAS↗
NASA NTRS1993

The phase composition of Triton's polar caps

Triton's polar caps are modeled as permanent nitrogen deposits hundreds of meters thick. Complex temperature variations on Triton's surface induce reversible transitions between the cubic and hexagonal phases of solid nitrogen, often with two coexisting propagating transition fronts. Subsurface temperature distributions are calculated using a two-dimensional thermal model with phase changes. The phase changes fracture the upper nitrogen layer, increasing its reflectivity and thus offering an explanation for the surprisingly high southern polar cap albedo (approximately 0.8) seen during the Voyager 2 flyby. The model has other implications for the phase transition phenomena on Triton, such as a plausible mechanism for the origin of geyser-like plume vent areas and a mechanism of energy transport toward them.

Duxbury, N. S.↗
NASA NTRS1988

Analytical and experimental investigation of the notched strength of thick laminates with surface notches

Experimental and analytical investigations of the fracture of a graphite/epoxy laminate with part-through semielliptic surface notches have been conducted. Experimental values of notched strength were determined for a variety of notch aspect ratios and notch depths. Analytical predictions of laminate failure were obtained from three different models. The three models were equally successful in predicting the fracture of laminates with deep notches. The experimental results of laminates with shallow notches indicated more notch sensitivity (that is greater strength reduction) than did the analytical models.

Harris, Charles E.↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
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↗
NASA NTRS1994

Microwave signatures and surface properties of Ovda Regio and surroundings, Venus

Magellan observations of horizontally and vertically polarized emissivity and radar specific cross sections are jointly modeled to separate dielectric constants from textural effect for lobate plains, fractured plains, fracture belts, festoon flows, and tessera in Ovda Regio and surroundings. The model assumes surface emission and scattering dominate and that signatures are controlled by a combination of quasi-specular and diffuse-scale mechanisms. Textural parameters are found to correlate with geology whereas dielectric constants are found to depend on elevation. A gradual increase in dielectric constant from rock-like values (4 to 5) to a value of approximately 50 is evident over the interval from 6054 to 6056 km, where elevations are derived from Magellan altimetry observations. Stereo radargrammetric analyses demonstrate that a return to rock-like values occurs over approximately 0.5 km above the terrain with highest dielectric constants. The abrupt return to rock-like values occurs independently of geology and reinforces the hypothesis of elevation control on dielectric constant values. Variations in dielectric constant with elevation are modeled using an 1000-ppm concentration of a ferroelectric mineral with a Curie temperature of 707 K. The model reproduces the abrupt change at highest elevations, as the Curie temperature is reached, and the gradual change as lower elevations are encountered. To our knowledge no other model explains the trends with such simplicity.

Arvidson, R. E.↗
DOE OSTI2021

Lithium vanadium oxide (Li 1.1 V 3 O 8 ) thick porous electrodes with high rate capacity: utilization and evolution upon extended cycling elucidated via operando energy dispersive X-ray diffraction and continuum simulation

The phase distribution of lithiated LVO in thick (~500 μm) porous electrodes (TPEs) designed to facilitate both ion and electron transport was determined using synchrotron-based operando energy dispersive X-ray diffraction (EDXRD). Probing 3 positions in the TPE while cycling at a 1C rate revealed a homogeneous phase transition across the thickness of the electrode at the 1st and 95th cycles. Additionally, continuum modelling indicated uniform lithiation across the TPE in agreement with the EDXRD results and ascribed decreasing accessible active material to be the cause of loss in delivered capacity between the 1st and 95th cycles. The model was supported by the observation of significant particle fracture by SEM consistent with loss of electrical contact. Overall, the combination of operando EDXRD, continuum modeling, and ex situ measurements enabled a deeper understanding of lithium vanadium oxide transport properties under high rate extended cycling within a thick highly porous electrode architecture.

25 ENERGY STORAGE↗