A discontinuous Galerkin finite element method for dynamic and wave propagation problems in non‐linear solids and saturated porous media

A time‐discontinuous Galerkin finite element method (DGFEM) for dynamics and wave propagation in non‐linear solids and saturated porous media is presented. The main distinct characteristic of the proposed DGFEM is that the specific P3–P1 interpolation approximation, which uses piecewise cubic (Hermite's polynomial) and linear interpolations for both displacements and velocities, in the time domain is particularly proposed. Consequently, continuity of the displacement vector at each discrete time instant is exactly ensured, whereas discontinuity of the velocity vector at the discrete time levels still remains. The computational cost is then obviously saved, particularly in the materially non‐linear problems, as compared with that required for the existing DGFEM. Both the implicit and explicit algorithms are developed to solve the derived formulations for linear and materially non‐linear problems. Numerical results illustrate good performance of the present method in eliminating spurious numerical oscillations and in providing much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain. Copyright © 2003 John Wiley & Sons, Ltd.     

A general discontinuous Galerkin method for finite hyperelasticity. Formulation and numerical applications

A discontinuous Galerkin formulation of the boundary value problem of finite‐deformation elasticity is presented. The primary purpose is to establish a discontinuous Galerkin framework for large deformations of solids in the context of statics and simple material behaviour with a view toward further developments involving behaviour or models where the DG concept can show its superiority compared to the continuous formulation. The method is based on a general Hu–Washizu–de Veubeke functional allowing for displacement and stress discontinuities in the domain interior. It is shown that this approach naturally leads to the formulation of average stress fluxes at interelement boundaries in a finite element implementation. The consistency and linearized stability of the method in the non‐linear range as well as its convergence rate are proven. An implementation in three dimensions is developed, showing that the proposed method can be integrated into conventional finite element codes in a straightforward manner. In order to demonstrate the versatility, accuracy and robustness of the method examples of application and convergence studies in three dimensions are provided. Copyright © 2006 John Wiley & Sons, Ltd.     

An adaptive spacetime discontinuous Galerkin method for cohesive models of elastodynamic fracture

This paper describes an adaptive numerical framework for cohesive fracture models based on a spacetime discontinuous Galerkin (SDG) method for elastodynamics with elementwise momentum balance. Discontinuous basis functions and jump conditions written with respect to target traction values simplify the implementation of cohesive traction–separation laws in the SDG framework; no special cohesive elements or other algorithmic devices are required. We use unstructured spacetime grids in a h‐adaptive implementation to adjust simultaneously the spatial and temporal resolutions. Two independent error indicators drive the adaptive refinement. One is a dissipation‐based indicator that controls the accuracy of the solution in the bulk material; the second ensures the accuracy of the discrete rendering of the cohesive law. Applications of the SDG cohesive model to elastodynamic fracture demonstrate the effectiveness of the proposed method and reveal a new solution feature: an unexpected quasi‐singular structure in the velocity response. Numerical examples demonstrate the use of adaptive analysis methods in resolving this structure, as well as its importance in reliable predictions of fracture kinetics. Copyright © 2009 John Wiley & Sons, Ltd.     

A fracture framework for Euler–Bernoulli beams based on a full discontinuous Galerkin formulation/extrinsic cohesive law combination

A new full Discontinuous Galerkin discretization of Euler–Bernoulli beam is presented. The main interest of this framework is its ability to simulate fracture problems by inserting a cohesive zone model in the formulation. With a classical Continuous Galerkin method, the use of the cohesive zone model is difficult because inserting a cohesive element between bulk elements is not straightforward. On one hand if the cohesive element is inserted at the beginning of the simulation, there is a modification of the structure stiffness and on the other hand inserting the cohesive element during the simulation requires modification of the mesh during computation. These drawbacks are avoided with the presented formulation as the structure is discretized in a stable and consistent way with full discontinuous elements and inserting cohesive elements during the simulation becomes straightforward. A new cohesive law based on the resultant stresses (bending moment and membrane) of the thin structure discretization is also presented. This model allows propagating fracture while avoiding through‐the‐thickness integration of the cohesive law. Tests are performed to show that the proposed model releases, during the fracture process, an energy quantity equal to the fracture energy for any combination of tension‐bending loadings. Copyright © 2010 John Wiley & Sons, Ltd.     

On continuous,discontinuous, mixed,and primal hybrid finite element methods for second‐order elliptic problems

Finite element formulations for second‐order elliptic problems, including the classic H¹‐conforming Galerkin method, dual mixed methods, a discontinuous Galerkin method, and two primal hybrid methods, are implemented and numerically compared on accuracy and computational performance. Excepting the discontinuous Galerkin formulation, all the other formulations allow static condensation at the element level, aiming at reducing the size of the global system of equations. For a three‐dimensional test problem with smooth solution, the simulations are performed with h‐refinement, for hexahedral and tetrahedral meshes, and uniform polynomial degree distribution up to four. For a singular two‐dimensional problem, the results are for approximation spaces based on given sets of hp‐refined quadrilateral and triangular meshes adapted to an internal layer. The different formulations are compared in terms of L²‐convergence rates of the approximation errors for the solution and its gradient, number of degrees of freedom, both with and without static condensation. Some insights into the required computational effort for each simulation are also given.     

Numerical analysis of dissolution processes in cementitious materials using discontinuous and continuous Galerkin time integration schemes

The present paper is concerned with the numerical integration of non‐linear reaction–diffusion problems by means of discontinuous and continuous Galerkin methods. The first‐order semidiscrete initial value problem of calcium leaching of cementitious materials, based on a phenomenological dissolution model, an electrolyte diffusion model and the spatial p‐finite element discretization, is used as a highly non‐linear model problem. A p‐finite element method is used for the spatial discretization. In the context of discontinuous Galerkin methods the semidiscrete mass balance and the continuity of the primary variables are weakly formulated within time steps and between time steps, respectively. Continuous Galerkin methods are obtained by the strong enforcement of the continuity condition as special cases. The introduction of a natural time co‐ordinate allows for the application of standard higher order temporal shape functions of the p‐Lagrange type and the well‐known Gauss–Legendre quadrature of associated time integrals. It is shown, that arbitrary order accurate integration schemes can be developed within the framework of the proposed temporal p‐Galerkin methods. Selected benchmark analyses of calcium dissolution demonstrate the robustness of these methods with respect to pronounced changes of the reaction term and non‐smooth changes of Dirichlet boundary conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

A hybrid discontinuous in space and time Galerkin method for wave propagation problems

Medium‐frequency regime and multi‐scale wave propagation problems have been a subject of active research in computational acoustics recently. New techniques have attempted to overcome the limitations of existing discretization methods that tend to suffer from dispersion. One such technique, the discontinuous enrichment method, incorporates features of the governing partial differential equation in the approximation, in particular, the solutions of the homogeneous form of the equation. Here, based on this concept and by extension of a conventional space–time finite element method, a hybrid discontinuous Galerkin method (DGM) for the numerical solution of transient problems governed by the wave equation in two and three spatial dimensions is described. The discontinuous formulation in both space and time enables the use of solutions to the homogeneous wave equation in the approximation. In this contribution, within each finite element, the solutions in the form of polynomial waves are employed. The continuity of these polynomial waves is weakly enforced through suitably chosen Lagrange multipliers. Results for two‐dimensional and three‐dimensional problems, in both low‐frequency and medium‐frequency regimes, show that the proposed DGM outperforms the conventional space–time finite element method. Copyright © 2014 John Wiley & Sons, Ltd.     

Blending in the extended finite element method by discontinuous Galerkin and assumed strain methods

In the extended finite element method (XFEM), errors are caused by parasitic terms in the approximation space of the blending elements at the edge of the enriched subdomain. A discontinuous Galerkin (DG) formulation is developed, which circumvents this source of error. A patch‐based version of the DG formulation is developed, which decomposes the domain into enriched and unenriched subdomains. Continuity between patches is enforced with an internal penalty method. An element‐based form is also developed, where each element is considered a patch. The patch‐based DG is shown to have similar accuracy to the element‐based DG for a given discretization but requires significantly fewer degrees of freedom. The method is applied to material interfaces, cracks and dislocation problems. For the dislocations, a contour integral form of the internal forces that only requires integration over the patch boundaries is developed. A previously developed assumed strain (AS) method is also developed further and compared with the DG method for weak discontinuities and linear elastic cracks. The DG method is shown to be significantly more accurate than the standard XFEM for a given element size and to converge optimally, even where the standard XFEM does not. The accuracy of the DG method is similar to that of the AS method but requires less application‐specific coding. Copyright © 2007 John Wiley & Sons, Ltd.     

粘性不可压流动的全离散非线性Galerkin方法

这篇文章讨论了有限元情形下，粘性不可压流动的全离散非线性Ｇａｌｅｒｋｉｎ方法，此外，给出了用非线性Ｇａｌｅｒｋｉｎ方法进行的数值试验结果，和通常的Ｇａｌｅｒｋｉｎ方法比较，非线性Ｇａｌｅｒｋｉｎ方法是一种稳定性好，计算量少以及逼近精度高的数值方法。     

On the accuracy of finite volume and discontinuous Galerkin discretizations for compressible flow on unstructured grids

This paper presents a comparison between two high‐order methods. The first one is a high‐order finite volume (FV) discretization on unstructured grids that uses a meshfree method (moving least squares (MLS)) in order to construct a piecewise polynomial reconstruction and evaluate the viscous fluxes. The second method is a discontinuous Galerkin (DG) scheme. Numerical examples of inviscid and viscous flows are presented and the solutions are compared. The accuracy of both methods, for the same grid resolution, is similar, although the DG scheme requires a larger number of degrees of freedom than the FV–MLS method. Copyright © 2009 John Wiley & Sons, Ltd.     

Solving the advection–diffusion equation on unstructured meshes with discontinuous/mixed finite elements and a local time stepping procedure

Explicit schemes are known to provide less numerical diffusion in solving the advection–diffusion equation, especially for advection‐dominated problems. Traditional explicit schemes use fixed time steps restricted by the global CFL condition in order to guarantee stability. This is known to slow down the computation especially for heterogeneous domains and/or unstructured meshes. To avoid this problem, local time stepping procedures where the time step is allowed to vary spatially in order to satisfy a local CFL condition have been developed. In this paper, a local time stepping approach is used with a numerical model based on discontinuous Galerkin/mixed finite element methods to solve the advection–diffusion equation. The developments are detailed for general unstructured triangular meshes. Numerical experiments are performed to show the efficiency of the numerical model for the simulation of (i) the transport of a solute on highly unstructured meshes and (ii) density‐driven flow, where the velocity field changes at each time step. The model gives stable results with significant reduction of the computational cost especially for the non‐linear problem. Moreover, numerical diffusion is also reduced for highly advective problems. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of photonic crystals using the hybrid finite‐element/finite‐difference time domain technique based on the discontinuous Galerkin method

Two‐dimensional photonic crystal structures are analyzed by a recently developed hybrid technique combining the finite‐element time‐domain (FETD) method and the finite‐difference time‐domain (FDTD) method. This hybrid FETD/FDTD method uses the discontinuous Galerkin method as framework for domain decomposition. To the best of our knowledge, this is the first hybrid FETD/FDTD method that allows non‐conformal meshes between different FETD and FDTD subdomains. It is also highly parallelizable. These properties are very suitable for the computation of periodic structures with curved surfaces. Numerical examples for the computation of the scattering parameters of two‐dimensional photonic bandgap structures are presented as applications of the hybrid FETD/FDTD method. Numerical results demonstrate the efficiency and accuracy of the proposed hybrid method. Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order immersed boundary discontinuous‐Galerkin method for Poisson's equation with discontinuous coefficients and singular sources

We adopt a numerical method to solve Poisson's equation on a fixed grid with embedded boundary conditions, where we put a special focus on the accurate representation of the normal gradient on the boundary. The lack of accuracy in the gradient evaluation on the boundary is a common issue with low‐order embedded boundary methods. Whereas a direct evaluation of the gradient is preferable, one typically uses post‐processing techniques to improve the quality of the gradient. Here, we adopt a new method based on the discontinuous‐Galerkin (DG) finite element method, inspired by the recent work of [A.J. Lew and G.C. Buscaglia. A discontinuous‐Galerkin‐based immersed boundary method. International Journal for Numerical Methods in Engineering, 76:427‐454, 2008]. The method has been enhanced in two aspects: firstly, we approximate the boundary shape locally by higher‐order geometric primitives. Secondly, we employ higher‐order shape functions within intersected elements. These are derived for the various geometric features of the boundary based on analytical solutions of the underlying partial differential equation. The development includes three basic geometric features in two dimensions for the solution of Poisson's equation: a straight boundary, a circular boundary, and a boundary with a discontinuity. We demonstrate the performance of the method via analytical benchmark examples with a smooth circular boundary as well as in the presence of a singularity due to a re‐entrant corner. Results are compared to a low‐order extended finite element method as well as the DG method of [1]. We report improved accuracy of the gradient on the boundary by one order of magnitude, as well as improved convergence rates in the presence of a singular source. In principle, the method can be extended to three dimensions, more complicated boundary shapes, and other partial differential equations. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparison of implicit and explicit hybridizable discontinuous Galerkin methods for the acoustic wave equation

We describe implicit and explicit formulations of the hybridizable discontinuous Galerkin method for the acoustic wave equation based on state‐of‐the‐art numerical software and quantify their efficiency for realistic application settings. In the explicit scheme, the trace of the acoustic pressure is computed from the solution on the two elements adjacent to the face at the old time step. Tensor product shape functions for quadrilaterals and hexahedra evaluated with sum factorization are used to ensure low operation counts. For applying the inverse mass matrix of Lagrangian shape functions with full Gaussian quadrature, a new tensorial technique is proposed. As time propagators, diagonally implicit and explicit Runge–Kutta methods are used, respectively. We find that the computing time per time step is 25 to 200 times lower for the explicit scheme, with an increasing gap in three spatial dimensions and for higher element degrees. Our experiments on realistic 3D wave propagation with variable material parameters in a photoacoustic imaging setting show an improvement of two orders of magnitude in terms of time to solution, despite stability restrictions on the time step of the explicit scheme. Operation counts and a performance model to predict performance on other computer systems accompany our results. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of continuous and discontinuous finite element methods for parabolic differential equations employing implicit time integration

This article compares the computational cost, stability, and accuracy of continuous and discontinuous Galerkin Finite Element Methods (GFEM) for various parabolic differential equations including the advection–diffusion equation, viscous Burgers’ equation, and Turing pattern formation equation system. The results show that, for implicit time integration, the continuous GFEM is typically 5–20 times less computationally expensive than the discontinuous GFEM using the same finite element mesh and element order. However, the discontinuous GFEM is significantly more stable than the continuous GFEM for advection dominated problems and is able to obtain accurate approximate solutions for cases where the classic, un-stabilized continuous GFEM fails.     

A full‐discontinuous Galerkin formulation of nonlinear Kirchhoff–Love shells: elasto‐plastic finite deformations,parallel computation,and fracture applications

Because of its ability to take into account discontinuities, the discontinuous Galerkin (DG) method presents some advantages for modeling cracks initiation and propagation. This concept has been recently applied to three‐dimensional simulations and to elastic thin bodies. In this last case, the assumption of small elastic deformations before cracks initiation or propagation reduces drastically the applicability of the framework to a reduced number of materials. To remove this limitation, a full‐DG formulation of nonlinear Kirchhoff–Love shells is presented and is used in combination with an elasto‐plastic finite deformations model. The results obtained by this new formulation are in agreement with other continuum elasto‐plastic shell formulations. Then, this full‐DG formulation of Kirchhoff–Love shells is coupled with the cohesive zone model to perform thin body fracture simulations. As this method considers elasto‐plastic constitutive laws in combination with the cohesive model, accurate results compared with the experiments are found. In particular, the crack path and propagation rate of a blasted cylinder are shown to match experimental results. One of the main advantages of this framework is its ability to run in parallel with a high speed‐up factor, allowing the simulation of ultra fine meshes. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of two‐grid methods for reaction‐diffusion equations by expanded mixed finite element methods

We present two efficient methods of two‐grid scheme for the approximation of two‐dimensional semi‐linear reaction‐diffusion equations using an expanded mixed finite element method. To linearize the discretized equations, we use two Newton iterations on the fine grid in our methods. Firstly, we solve an original non‐linear problem on the coarse grid. Then we use twice Newton iterations on the fine grid in our first method, and while in second method we make a correction on the coarse grid between two Newton iterations on the fine grid. These two‐grid ideas are from Xu's work (SIAM J. Sci. Comput. 1994; 15 :231–237; SIAM J. Numer. Anal. 1996; 33 :1759–1777) on standard finite element method. We extend the ideas to the mixed finite element method. Moreover, we obtain the error estimates for two algorithms of two‐grid method. It is showed that coarse space can be extremely coarse and we achieve asymptotically optimal approximation as long as the mesh sizes satisfy H =??(h^¼) in the first algorithm and H =??(h^?) in second algorithm. Copyright © 2006 John Wiley & Sons, Ltd.     

Embedded representation of fracture in concrete with mixed finite elements

As an alternative to the smeared and discrete crack representations, an embedded representation of fracture for finite element analysis of concrete structures is presented. The three-field Hu–Washizu variational statement is extended to bodies with internal discontinuities. The extended variational statement is then utilized for formulating elements with a discontinuous displacement field. The new elements are capable of modelling different deformation modes of an internal discontinuity at the element level. The satisfactory performance of the embedded crack representation is verified by several case studies on concrete fracture.     

A stabilized volume‐averaging finite element method for flow in porous media and binary alloy solidification processes

A stabilized equal‐order velocity–pressure finite element algorithm is presented for the analysis of flow in porous media and in the solidification of binary alloys. The adopted governing macroscopic conservation equations of momentum, energy and species transport are derived from their microscopic counterparts using the volume‐averaging method. The analysis is performed in a single domain with a fixed numerical grid. The fluid flow scheme developed includes SUPG (streamline‐upwind/Petrov–Galerkin), PSPG (pressure stabilizing/Petrov–Galerkin) and DSPG (Darcy stabilizing/Petrov–Galerkin) stabilization terms in a variable porosity medium. For the energy and species equations a classical SUPG‐based finite element method is employed. The developed algorithms were tested extensively with bilinear elements and were shown to perform stably and with nearly quadratic convergence in high Rayleigh number flows in varying porosity media. Examples are shown in natural and double diffusive convection in porous media and in the directional solidification of a binary‐alloy. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupled meshfree and fractal finite element method for mixed mode two‐dimensional crack problems

This paper presents a coupling technique for integrating the element‐free Galerkin method (EFGM) with the fractal finite element method (FFEM) for analyzing homogeneous, isotropic, and two‐dimensional linear‐elastic cracked structures subjected to mixed‐mode (modes I and II) loading conditions. FFEM is adopted for discretization of the domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprise both the EFG and the finite element (FE) shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no special enriched basis functions or no structured mesh with special FEs are necessary and no post‐processing (employing any path independent integrals) is needed to determine fracture parameters, such as stress‐intensity factors (SIFs) and T‐stress. The numerical results show that SIFs and T‐stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also, a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions. A numerical example on mixed‐mode condition is presented to simulate crack propagation. As in the proposed coupled EFGM–FFEM at each increment during the crack propagation, the FFEM mesh (around the crack tip) is shifted as it is to the new updated position of the crack tip (such that FFEM mesh center coincides with the crack tip) and few meshless nodes are sprinkled in the location where the FFEM mesh was lying previously, crack‐propagation analysis can be dramatically simplified. Copyright © 2010 John Wiley & Sons, Ltd.