Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended finite element method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Fatigue Crack Propagation of Nickel-Based Superalloy: Experiments and Simulations with Extended Finite Element Method

Numerical simulation based on extended finite element method was employed to investigate the fatigue crack propagation of nickel-based superalloy at room temperature. Experimental tests on compact tension specimens have performed to obtain fatigue crack propagation parameters in Paris region. The extended finite element method has presented a new approach to solve the stress intensity factors and can effectively predict crack propagation without re-meshing at crack tip. The simulation results are in good accordance with experimental data in real 3D cases.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling curved interfaces without element‐partitioning in the extended finite element method

Summary In this paper, we model holes and material interfaces (weak discontinuities) in two‐dimensional linear elastic continua using the extended finite element method on higher‐order (spectral) finite element meshes. Arbitrary parametric curves such as rational Bézier curves and cubic Hermite curves are adopted in conjunction with the level set method to represent curved interfaces. Efficient computation of weak form integrals with polynomial integrands is realized via the homogeneous numerical integration scheme—a method that uses Euler's homogeneous function theorem and Stokes' theorem to reduce integration to the boundary of the domain. Numerical integration on cut elements requires the evaluation of a one‐dimensional integral over a parametric curve, and hence, the need to partition curved elements is eliminated. To improve stiffness matrix conditioning, ghost penalty stabilization and the Jacobi preconditioner are used. For material interface problems, we develop an enrichment function that captures weak discontinuities on spectral meshes. Taken together, we show through numerical experiments that these advances deliver optimal algebraic rates of convergence with h ‐refinement ( p =1,2,…,5) and exponential rates of convergence with p ‐refinement ( p =1,2,…,7) for elastostatic problems with holes and material inclusions on Cartesian p th‐order spectral finite element meshes.

Chin, Eric B.↗
DOE OSTI2019

Ceramic nuclear fuel fracture modeling with the extended finite element method

Ceramic fuel pellets used in nuclear light water reactors experience significant fracture due to the high thermal gradients experienced under normal operating conditions. This has important effects on the performance of the fuel system. Because of this, a realistic, physically based fracture modeling capability is essential to predict fuel behavior in a wide variety of normal and off-normal conditions. The extended finite element method (X-FEM) is a powerful method to represent arbitrary propagating discrete cracks in finite element models, and has many characteristics that make it attractive for nuclear fuel performance analysis. This paper describes the implementation of X-FEM in a multiphysics fuel performance code and presents applications of that capability. These applications include several thermal mechanics fracture benchmark problems, which demonstrate the accuracy of this approach. It also includes application of this capability to study nuclear fuel fracture, both on stationary and propagating cracks. The study on stationary cracks shows the effects that varying the number of radial cracks and the length of those cracks have on the energy release rate. Here, the propagating crack case demonstrates random initiation and subsequent propagation of interacting thermally induced cracks during an initial ramp to full power with fresh fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Pellet Cladding Mechanical Interaction Modeling Using the Extended Finite Element Method

As a brittle material, the ceramic UO2 used as light water reactor fuel experiences significant fracturing throughout its life, beginning with the first rise to power of fresh fuel. This has multiple effects on the thermal and mechanical response of the fuel/cladding system. One such effect that is particularly important is that when there is mechanical contact between the fuel and cladding, cracks that extending from the outer surface of the fuel into the volume of the fuel cause elevated stresses in the adjacent cladding, which can potentially lead to cladding failure. Modeling the thermal and mechanical response of the cladding in the vicinity of these surface-breaking cracks in the fuel can provide important insights into this behavior to help avoid operating conditions that could lead to cladding failure. Such modeling has traditionally been done in the context of finite-element-based fuel performance analysis by modifying the fuel mesh to introduce discrete cracks. While this approach is effective in capturing the important behavior at the fuel/cladding interface, there are multiple drawbacks to explicitly incorporating the cracks in the finite element mesh. Because the cracks are incorporated in the original mesh, the mesh must be modified for cracks of specified location and depth, so it is difficult to account for crack propagation and the formation of new cracks at other locations. The extended finite element method (XFEM) has emerged in recent years as a powerful method to represent arbitrary, evolving, discrete discontinuities within the context of the finite element method. Development work is underway by the authors to implement XFEM in the BISON fuel performance code, and this capability has previously been demonstrated in simulations of fracture propagation in ceramic nuclear fuel. These preliminary demonstrations have included only the fuel, and excluded the cladding for simplicity. This paper presents initial results of efforts to apply XFEM to model stress concentrations induced by fuel fractures at the fuel/cladding interface during pellet cladding mechanical interaction (PCMI). This is accomplished by enhancing the thermal and mechanical contact enforcement algorithms employed by BISON to permit their use in conjunction with XFEM. The results from this methodology are demonstrated to be equivalent to those from using meshed discrete cracks. While the results of the two methods are equivalent for the case of a stationary crack, it is demonstrated that XFEM provides the additional flexibility of allowing arbitrary crack initiation and propagation during the analysis, and minimizes model setup effort for cases with stationary cracks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Creep and Creep Fracture Modeling with Surrogate Creep Models and the Extended Finite Element Method

Alloy components in advanced nuclear reactors will be subjected to environmental conditions that could include high temperatures, irradiation, and exposure to corrosive salts. These conditions could lead to the formation of crack-like defects, which could grow over time in a mechanism known as creep crack growth (CCG). Predicting growth rates of these defects is important for assessing the safe operating life of advanced reactors. This project documents progress toward developing and testing next-generation data-driven constitutive models for deformation creep. It also documents the application of the extended finite element method in conjunction with surrogate creep models to predict CCG parameters under a variety of conditions. These important incremental developments contribute to the longer-term objective of developing microstructure-aware constitutive models that can be used for predicting creep deformation and CCG at the component scale with improved accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Modeling 3D PCMI using the Extended Finite Element Method with higher order elements

This report documents the recent development to enable XFEM to work with higher order elements. It also demonstrates the application of higher order (quadratic) elements to both 2D and 3D models of PCMI problems, where discrete fractures in the fuel are represented using XFEM. The modeling results demonstrate the ability of the higher order XFEM to accurately capture the effects of a crack on the response in the vicinity of the intersecting surfaces of cracked fuel and cladding, as well as represent smooth responses in the regions away from the crack.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗
DOE OSTI2018

Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions

A discrete De Rham complex enables compatible, structure-preserving discretizations for a broad range of partial differential equations problems. Such discretizations can correctly reproduce the physics of interface problems, provided the grid conforms to the interface. However, large deformations, complex geometries, and evolving interfaces makes generation of such grids difficult. We develop and demonstrate two formally equivalent approaches that, for a given background mesh, dynamically construct an interface-conforming discrete De Rham complex. Both approaches start by dividing cut elements into interface-conforming subelements but differ in how they build the finite element basis on these subelements. The first approach discards the existing non-conforming basis of the parent element and replaces it by a dynamic set of degrees of freedom of the same kind. The second approach defines the interface-conforming degrees of freedom on the subelements as superpositions of the basis functions of the parent element. These approaches generalize the Conformal Decomposition Finite Element Method (CDFEM) and the extended finite element method with algebraic constraints (XFEM-AC), respectively, across the De Rham complex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

An arc-length method for controlled cohesive crack propagation using high-order XFEM and Irwin’s crack closure integral

Numerical modeling of cohesive crack growth in quasi-brittle materials is challenging, primarily due to the combination of (i) nonlinearity associated with the fracture process zone (FPZ), (ii) arbitrary directions to which a crack may propagate, and (iii) snap-back or snap-through instabilities encountered in the response of the structure. To address these challenges, here we propose a novel arc-length method that can follow the equilibrium path of cohesive crack propagation. The proposed approach is based on the extended finite element method (XFEM) with scalar high-order enrichment functions and Irwin’s crack closure integral, which allows for direct control of the applied loads necessary to propagate cohesive cracks. This is achieved by augmenting a constraint equation written in terms of stress intensity factors (SIFs), and expressed explicitly in terms of the enriched degrees of freedom, which is an attractive feature achieved with Irwin’s integral, since SIFs can be written in closed-form. Note that singular enrichments are active in an unstable crack propagation state and automatically vanish in stable crack configurations. Furthermore, to propagate cracks in arbitrary directions, we employ a maximum circumferential stress criterion implemented by (i) direct usage of the SIFs, and by (ii) a new stress-based nonlocal implementation of this principle. Various benchmark problems including pure mode I and mixed-mode fracture are solved to demonstrate the predictive capability of the present framework for cohesive crack modeling.

42 ENGINEERING↗
DOE OSTI2022

Efficient co-solution of time step size and independent state in simulations of fluid-driven fracture propagation with embedded meshes

Here we present an efficient time-continuation scheme for fluid-driven fracture propagation problems in the extended finite element method framework. The approach applies a monolithic solution strategy to a fully coupled and implicit approximation of hydro-mechanical systems in conjunction with simultaneous linear elastic propagation of multiple fractures. At the end of each time step, the process ensures that the weakest fracture tip is in an equilibrium propagation regime. Furthermore, the solution process provides an initialization procedure for the newly created fracture spaces and an a priori estimate of the stress intensity factor growth rate, improving simulation robustness, and efficiency. The solution process is validated using the Kristianovich-Geertsma-de Klerk analytical solution under the toughness- and viscosity-dominated regimes. It is also extended to and demonstrated on problems with multiple fractures undergoing simultaneous propagation with stress shadow interactions. Numerical examples demonstrate that the solution process can reduce the required computational cost by one order of magnitude compared to other existing methods.

42 ENGINEERING↗
DOE OSTI2023

Initial Fracture Propagation Modeling of Graphite Components with Grizzly

Graphite has historically been extensively used in power reactor cores and will be used in multiple types of advanced reactors currently under development. These graphite structural components can experience significant stresses due to nonuniform volumetric strains induced by irradiation and thermal expansion, which can lead to fracture. Robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors are important for evaluating component integrity, developing design standards, and interpreting experimental results to characterize graphite performance. The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program has been developing degradation models for other structural components in nuclear reactors within the Grizzly and BlackBear codes. This report documents an effort to develop initial capabilities for modeling graphite fracture within these codes, building on prior efforts to model fracture in other materials. Major elements of this effort include developing a new system for modeling fracture nucleation and growth in two dimensions using the extended finite element method and incorporating a damage and plasticity model. These capabilities are applied here to model a representative graphite component and a splitting disc experiment used to obtain tensile strength.

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