On the implementation of the eXtended Finite Element Method (XFEM) for interface problems
Authors: Thomas Carraro, Sven Wetterauer
Abstract
The eXtended Finite Element Method (XFEM) is used to solve interface problems with an unfitted mesh. We present an implementation of the XFEM in the FEM-library deal.II. The main parts of the implementation are (i) the appropriate quadrature rule; (ii) the shape functions for the extended part of the finite element formulation; (iii) the boundary and interface conditions. We show how to handle the XFEM formulation providing a code that demonstrates the solution of two exemplary interface problems for a strong and a weak discontinuity respectively. In the weak discontinuity case, the loss of conformity due to the blending effect and its remedy are discussed. Furthermore, the optimal convergence of the presented unfitted method is numerically verified.
The papers around it, by date and closeness
Higher is closer to this paper; left is older. Lines join papers close to each other: groups of lines are sub-topics. Select a point to see why it is there. Back to the summary of this paper
Papers around this one
Ranked by how close their title and abstract are to this paper’s. Select some to read them together.
Closest papers
-
Thirty years of the finite volume method for solid mechanicsPhilip Cardiff, Ismet Demirdžić · 2018 · Numerical Analysis · AI summary available
Abstract
Since early publications in the late 1980s and early 1990s, the finite volume method has been shown suitable for solid mechanics analyses. At present, there are several flavours of the method, which can be classified in a variety of ways, such as grid arrangement (cell-centred vs staggered vs vertex-centred), solution algorithm (implicit vs explicit), and stabilisation strategy (Rhie-Chow vs Jameson-Schmidt-Turkel vs Godunov upwinding). This article gives an overview, historical perspective, comparison and critical analysis of the different approaches where a close comparison with the de facto standard for computational solid mechanics, the finite element method, is given. The article finish…
-
Compatible finite element methods for geophysical fluid dynamicsColin J. Cotter · 2023 · Numerical Analysis · AI summary available
Abstract
This article surveys research on the application of compatible finite element methods to large scale atmosphere and ocean simulation. Compatible finite element methods extend Arakawa's C-grid finite difference scheme to the finite element world. They are constructed from a discrete de Rham complex, which is a sequence of finite element spaces which are linked by the operators of differential calculus. The use of discrete de Rham complexes to solve partial differential equations is well established, but in this article we focus on the specifics of dynamical cores for simulating weather, oceans and climate. The most important consequence of the discrete de Rham complex is the Hodge-Helmholtz d…
-
A robust monolithic solver for phase-field fracture integrated with fracture energy based arc-length method and under-relaxationRitukesh Bharali, Somdatta Goswami, Cosmin Anitescu et al. · 2021 · Numerical Analysis · AI summary available
Abstract
The phase-field fracture free-energy functional is non-convex with respect to the displacement and the phase field. This results in a poor performance of the conventional monolithic solvers like the Newton-Raphson method. In order to circumvent this issue, researchers opt for the alternate minimization (staggered) solvers. Staggered solvers are robust for the phase-field based fracture simulations as the displacement and the phase-field sub-problems are convex in nature. Nevertheless, the staggered solver requires very large number of iterations (of the order of thousands) to converge. In this work, a robust monolithic solver is presented for the phase-field fracture problem. The solver adop…
-
An embedded corrector problem for homogenization. Part II: Algorithms and discretizationEric Cancès, Virginie Ehrlacher, Frederic Legoll et al. · 2018 · Numerical Analysis · Not summarized yet
Abstract
This contribution is the numerically oriented companion article of the work [E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131]. We focus here on the numerical resolution of the embedded corrector problem introduced in [E. Canc\`es, V. Ehrlacher, F. Legoll and B. Stamm, CRAS 2015; E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131] in the context of homogenization of diffusion equations. Our approach consists in considering a corrector-type problem, posed on the whole space, but with a diffusion matrix which is constant outside some bounded domain. In [E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv p…
-
A robust solution strategy for the Cahn-Larché equationsErlend Storvik, Jakub Wiktor Both, Jan Martin Nordbotten et al. · 2022 · Numerical Analysis · AI summary available
Abstract
In this paper we propose a solution strategy for the Cahn-Larch\'e equations, which is a model for linearized elasticity in a medium with two elastic phases that evolve subject to a Ginzburg-Landau type energy functional. The system can be seen as a combination of the Cahn-Hilliard regularized interface equation and linearized elasticity, and is non-linearly coupled, has a fourth order term that comes from the Cahn-Hilliard subsystem, and is non-convex and nonlinear in both the phase-field and displacement variables. We propose a novel semi-implicit discretization in time that uses a standard convex-concave splitting method of the nonlinear double-well potential, as well as special treatment…
-
Modeling martensitic transformation in shape memory alloys using multi-phase-field elasticity models based on partial rank-one energy relaxation on pairwise interfacesMohammad Sarhil, Oleg Shchyglo, Dominik Brands et al. · 2023 · Numerical Analysis · Not summarized yet
Abstract
To model the mechanically-driven phase transformations, e.g. martensitic transformation, using the phase-field theory, suitable models are needed for describing the mechanical fields of the individual non-vanishing phase-fields in the interface regions in order to obtain the mechanical driving forces of phase-field motion. Quantitative modeling requires satisfying the interfacial static equilibrium and kinematic compatibility conditions which have already been achieved in the literature for dual-phase-field materials by using the rank-one relaxation (or convexification) of the energy density. A direct generalization to the multi-phase-field case is not applicable without breaking these condi…
-
A Computationally Efficient Vectorized Implementation of Localizing Gradient Damage Method in MATLABSubrato Sarkar · 2023 · Numerical Analysis · AI summary available
Abstract
In this work, a recently developed fracture modeling method called localizing gradient damage method (LGDM) is implemented in MATLAB. MATLAB is well-known in the computational research community for its simple and easy-to-learn coding interface. As a result, MATLAB is generally preferred for the initial development (prototyping) of computational models by researchers. However, MATLAB-developed codes are seldom used for large-scale simulations (after initial development is complete) due to their computational inefficiency. Hence, a computationally efficient implementation of LGDM using MATLAB vectorization is presented in this work. The choice of LGDM (as the fracture modeling method) is base…
-
DelayDiffEq: Generating Delay Differential Equation Solvers via Recursive Embedding of Ordinary Differential Equation SolversDavid Widmann, Chris Rackauckas · 2022 · Numerical Analysis · AI summary available
Abstract
Traditional solvers for delay differential equations (DDEs) are designed around only a single method and do not effectively use the infrastructure of their more-developed ordinary differential equation (ODE) counterparts. In this work we present DelayDiffEq, a Julia package for numerically solving delay differential equations (DDEs) which leverages the multitude of numerical algorithms in OrdinaryDiffEq for solving both stiff and non-stiff ODEs, and manages to solve challenging stiff DDEs. We describe how compiling the ODE integrator within itself, and accounting for discontinuity propagation, leads to a design that is effective for DDEs while using all of the ODE internals. We highlight som…
Go further
How are these papers found?
We compare the title and abstract of this paper with those of every paper summarized on this site, and with the papers Semantic Scholar recommends, those it cites and those citing it (for a much-cited paper, its most recent citations and, from OpenAlex, the most cited ones). Among the prior work and the follow-ups, the papers that use each other in depth, a method it uses and the much-cited follow-ups come a little higher. The levels (very close, close, related) come from the similarity of the texts, not from a judgment of quality; * marks a paper without an abstract, compared by its title. The sentence under a paper is written by an AI from the two abstracts: check it in the papers. A paper marked not relevant by three readers leaves the list.
For now only among the papers summarized on this site: search again to add the prior work, the follow-ups and Semantic Scholar’s recommendations.