A discontinuous‐Galerkin‐based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user‐defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous‐Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements, boundary locking is avoided and optimal‐order convergence is achieved. This is shown through numerical experiments in reaction–diffusion problems. Copyright © 2008 John Wiley & Sons, Ltd.     

An L∞–Lp mesh‐adaptive method for computing unsteady bi‐fluid flows

This paper discusses the contribution of mesh adaptation to high‐order convergence of unsteady multi‐fluid flow simulations on complex geometries. The mesh adaptation relies on a metric‐based method controlling the L ^p‐norm of the interpolation error and on a mesh generation algorithm based on an anisotropic Delaunay kernel. The mesh‐adaptive time advancing is achieved, thanks to a transient fixed‐point algorithm to predict the solution evolution coupled with a metric intersection in the time procedure. In the time direction, we enforce the equidistribution of the error, i.e. the error minimization in L ^∞ norm. This adaptive approach is applied to an incompressible Navier–Stokes model combined with a level set formulation discretized on triangular and tetrahedral meshes. Applications to interface flows under gravity are performed to evaluate the performance of this method for this class of discontinuous flows. Copyright © 2010 John Wiley & Sons, Ltd.     

A cut‐cell non‐conforming Cartesian mesh method for compressible and incompressible flow

This paper details a multigrid‐accelerated cut‐cell non‐conforming Cartesian mesh methodology for the modelling of inviscid compressible and incompressible flow. This is done via a single equation set that describes sub‐, trans‐, and supersonic flows. Cut‐cell technology is developed to furnish body‐fitted meshes with an overlapping mesh as starting point, and in a manner which is insensitive to surface definition inconsistencies. Spatial discretization is effected via an edge‐based vertex‐centred finite volume method. An alternative dual‐mesh construction strategy, similar to the cell‐centred method, is developed. Incompressibility is dealt with via an artificial compressibility algorithm, and stabilization achieved with artificial dissipation. In compressible flow, shocks are captured via pressure switch‐activated upwinding. The solution process is accelerated with full approximation storage (FAS) multigrid where coarse meshes are generated automatically via a volume agglomeration methodology. This is the first time that the proposed discretization and solution methods are employed to solve a single compressible–incompressible equation set on cut‐cell Cartesian meshes. The developed technology is validated by numerical experiments. The standard discretization and alternative methods were found equivalent in accuracy and computational cost. The multigrid implementation achieved decreases in CPU time of up to one order of magnitude. Copyright © 2007 John Wiley & Sons, Ltd.     

A numerical study of a semi‐algebraic multilevel preconditioner for the local discontinuous Galerkin method

In this work, we consider the local discontinuous Galerkin (LDG) method applied to second‐order elliptic problems arising in the modeling of single‐phase flows in porous media. It has been recently proven that the spectral condition number of the stiffness matrix exhibits an asymptotic behavior of ??(h^?2) on structured and unstructured meshes, where h is the mesh size. Thus, efficient preconditioners are mandatory. We present a semi‐algebraic multilevel preconditioner for the LDG method using local Lagrange‐type interpolatory basis functions. We show, numerically, that its performance does not degrade, or at least the number of iterations increases very slowly, as the number of unknowns augments. The preconditioner is tested on problems with high jumps in the coefficients, which is the typical scenario of problems arising in porous media. Copyright © 2007 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.     

A new stabilized finite element method for reaction–diffusion problems: The source‐stabilized Petrov–Galerkin method

This paper proposes a new stabilized finite element method to solve singular diffusion problems described by the modified Helmholtz operator. The Galerkin method is known to produce spurious oscillations for low diffusion and various alternatives were proposed to improve the accuracy of the solution. The mostly used methods are the well‐known Galerkin least squares and Galerkin gradient least squares (GGLS). The GGLS method yields the exact nodal solution in the one‐dimensional case and for a uniform mesh. However, the behavior of the method deteriorates slightly in the multi‐dimensional case and for non‐uniform meshes. In this work we propose a new stabilized finite element method that leads to improved accuracy for multi‐dimensional problems. For the one‐dimensional case, the new method leads to the same results as the GGLS method and hence provides exact nodal solutions to the problem on uniform meshes. The proposed method is a Galerkin discretization used to solve a modified equation that includes a term depending on the gradient of the original partial differential equation. Copyright © 2008 John Wiley & Sons, Ltd.     

Efficiency of high‐order elements for continuous and discontinuous Galerkin methods

To evaluate the computational performance of high‐order elements, a comparison based on operation count is proposed instead of runtime comparisons. More specifically, linear versus high‐order approximations are analyzed for implicit solver under a standard set of hypotheses for the mesh and the solution. Continuous and discontinuous Galerkin methods are considered in two‐dimensional and three‐dimensional domains for simplices and parallelotopes. Moreover, both element‐wise and global operations arising from different Galerkin approaches are studied. The operation count estimates show, that for implicit solvers, high‐order methods are more efficient than linear ones. Copyright © 2013 John Wiley & Sons, Ltd.     

Unstructured,curved elements for the two‐dimensional high order discontinuous control‐volume/finite‐element method

Quadrilateral and triangular elements with curved edges are developed in the framework of spectral, discontinuous, hybrid control‐volume/finite‐element method for elliptic problems. In order to accommodate hybrid meshes, encompassing both triangular and quadrilateral elements, one single mapping is used. The scheme is applied to two‐dimensional problems with discontinuous, anisotropic diffusion coefficients, and the exponential convergence of the method is verified in the presence of curved geometries. Copyright © 2014 John Wiley & Sons, Ltd.     

Weight‐adjusted discontinuous Galerkin methods: Matrix‐valued weights and elastic wave propagation in heterogeneous media

Weight‐adjusted inner products are easily invertible approximations to weighted L² inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time‐domain method for wave propagation which is low storage, energy stable, and high‐order accurate for arbitrary heterogeneous media and curvilinear meshes. In this work, we extend weight‐adjusted DG methods to the case of matrix‐valued weights, with the linear elastic wave equation as an application. We present a DG formulation of the symmetric form of the linear elastic wave equation, with upwind‐like dissipation incorporated through simple penalty fluxes. A semidiscrete convergence analysis is given, and numerical results confirm the stability and high‐order accuracy of weight‐adjusted DG for several problems in elastic wave propagation.     

A high‐order hybridizable discontinuous Galerkin method for elliptic interface problems

We present a high‐order hybridizable discontinuous Galerkin method for solving elliptic interface problems in which the solution and gradient are nonsmooth because of jump conditions across the interface. The hybridizable discontinuous Galerkin method is endowed with several distinct characteristics. First, they reduce the globally coupled unknowns to the approximate trace of the solution on element boundaries, thereby leading to a significant reduction in the global degrees of freedom. Second, they provide, for elliptic problems with polygonal interfaces, approximations of all the variables that converge with the optimal order of k + 1 in the L²(Ω)‐norm where k denotes the polynomial order of the approximation spaces. Third, they possess some superconvergence properties that allow the use of an inexpensive element‐by‐element postprocessing to compute a new approximate solution that converges with order k + 2. However, for elliptic problems with finite jumps in the solution across the curvilinear interface, the approximate solution and gradient do not converge optimally if the elements at the interface are isoparametric. The discrepancy between the exact geometry and the approximate triangulation near the curved interfaces results in lower order convergence. To recover the optimal convergence for the approximate solution and gradient, we propose to use superparametric elements at the interface. Copyright © 2012 John Wiley & Sons, Ltd.     

A face‐centred finite volume method for second‐order elliptic problems

This work proposes a novel finite volume paradigm, ie, the face‐centred finite volume (FCFV) method. Contrary to the popular vertex and cell‐centred finite volume methods, the novel FCFV defines the solution on the mesh faces (edges in two dimensions) to construct locally conservative numerical schemes. The idea of the FCFV method stems from a hybridisable discontinuous Galerkin formulation with constant degree of approximation, and thus inheriting the convergence properties of the classical hybridisable discontinuous Galerkin. The resulting FCFV features a global problem in terms of a piecewise constant function defined on the faces of the mesh. The solution and its gradient in each element are then recovered by solving a set of independent element‐by‐element problems. The mathematical formulation of FCFV for Poisson and Stokes equation is derived, and numerical evidence of optimal convergence in two dimensions and three dimensions is provided. Numerical examples are presented to illustrate the accuracy, efficiency, and robustness of the proposed methodology. The results show that, contrary to other finite volume methods, the accuracy of the FCFV method is not sensitive to mesh distortion and stretching. In addition, the FCFV method shows its better performance, accuracy, and robustness using simplicial elements, facilitating its application to problems involving complex geometries in three dimensions.     

An interface‐enriched generalized FEM for problems with discontinuous gradient fields

A new generalized FEM is introduced for solving problems with discontinuous gradient fields. The method relies on enrichment functions associated with generalized degrees of freedom at the nodes generated from the intersection of the phase interface with element edges. The proposed approach has several advantages over conventional generalized FEM formulations, such as a lower computational cost, easier implementation, and straightforward handling of Dirichlet boundary conditions. A detailed convergence study of the proposed method and a comparison with the standard FEM are presented for heat transfer problems. The method achieves the optimal rate of convergence using meshes that do not conform to the interfaces present in the domain while achieving a level of accuracy comparable to that of the standard FEM with conforming meshes. Various application problems are presented, including the conjugate heat transfer problem encountered in microvascular materials. Copyright © 2011 John Wiley & Sons, Ltd.     

Higher‐order extensions of a discontinuous Galerkin method for mid‐frequency Helmholtz problems

Recently, a discontinuous Galerkin method with plane wave basis functions and Lagrange multiplier degrees of freedom was proposed for the efficient solution of Helmholtz problems in the mid‐frequency regime. In this paper, this method is extended to higher‐order elements. Performance results obtained for various two‐dimensional problems highlight the advantages of these elements over classical higher‐order Galerkin elements such as Q₂ and Q₄ for the discretization of interior and exterior Helmholtz problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A p‐adaptive method for electromagnetic wave propagation

The discontinuous Galerkin FEM is used for the numerical solution of the three‐dimensional Maxwell equations. Control of errors in the numerical level for the divergence‐free constraint of the magnetic field can be obtained through the use of divergence‐free vector bases. In this work, the so‐called perfectly hyperbolic formulation of the Maxwell equations is used to retain both divergence‐free magnetic field and in the presence of charges to satisfy the Gauss constraint for the electric field at the numerical level. For both approaches, it is found that higher‐order approximations have favorable effect on the preservation of the divergence constraints and that the perfectly hyperbolic formulations retains these errors to a lower level. It is shown that high‐order accuracy in space and time is achieved in unstructured meshes using implicit time marching. For nonuniform meshes, local resolution refinement is used using p‐type adaptivity to ensure accurate electromagnetic wave propagation. Thus, the potential of the method to reach the required higher resolution in anisotropic meshes and obtain accurate electromagnetic wave propagation with reduced computational effort is demonstrated.     

Phase‐field modelling of solidification microstructure formation using the discontinuous finite element method

This paper presents a discontinuous Galerkin finite element computational methodology for solving the coupled phase‐field and heat diffusion equations to predict the microstructure evolution during solidification. The phase‐field modelling of microstructure formation is briefly discussed. The discontinuous Galerkin finite element formulation for the phase‐field model systems for solidification microstructure formation is described in detail. Numerical stability is performed using the Neumann method. The accuracy of the discontinuous model is verified by the analytic solution of a simple 1‐D solidification problem. Numerical calculations using the discontinuous finite element phase‐field model have been performed for simulating the complex 2‐D and 3‐D dendrite structures formed in supercooled melts and the results are compared well with those in literature using the finite difference methods. Parallel computing algorithm is presented and results show that the minimization of the intercommunication between microprocessors is the key to increase the effectiveness in parallel computing with the discontinuous finite element phase‐field model for solidification microstructure formation. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

A distortion measure to validate and generate curved high‐order meshes on CAD surfaces with independence of parameterization

A framework to validate and generate curved nodal high‐order meshes on Computer‐Aided Design (CAD) surfaces is presented. The proposed framework is of major interest to generate meshes suitable for thin‐shell and 3D finite element analysis with unstructured high‐order methods. First, we define a distortion (quality) measure for high‐order meshes on parameterized surfaces that we prove to be independent of the surface parameterization. Second, we derive a smoothing and untangling procedure based on the minimization of a regularization of the proposed distortion measure. The minimization is performed in terms of the parametric coordinates of the nodes to enforce that the nodes slide on the surfaces. Moreover, the proposed algorithm repairs invalid curved meshes (untangling), deals with arbitrary polynomial degrees (high‐order), and handles with low‐quality CAD parameterizations (independence of parameterization). Third, we use the optimization procedure to generate curved nodal high‐order surface meshes by means of an a posteriori approach. Given a linear mesh, we increase the polynomial degree of the elements, curve them to match the geometry, and optimize the location of the nodes to ensure mesh validity. Finally, we present several examples to demonstrate the features of the optimization procedure, and to illustrate the surface mesh generation process. Copyright © 2015 John Wiley & Sons, Ltd.     

The nsBETI method: an extension of the FETI method to non‐symmetrical BEM‐FEM coupled problems

The finite element tearing and interconnecting (FETI) method is recognized as an effective domain decomposition tool to achieve scalability in the solution of partitioned second‐order elasticity problems. In the boundary element tearing and interconnecting (BETI) method, a direct extension of the FETI algorithm to the BEM, the symmetric Galerkin BEM formulation, is used to obtain symmetric system matrices, making possible to apply the same FETI conjugate gradient solver. In this work, we propose a new BETI variant labeled nsBETI that allows to couple substructures modeled with the FEM and/or non‐symmetrical BEM formulations. The method connects non‐matching BEM and FEM subdomains using localized Lagrange multipliers and solves the associated non‐symmetrical flexibility equations with a Bi‐CGstab iterative algorithm. Scalability issues of nsBETI in BEM–BEM and combined BEM–FEM coupled problems are also investigated. Copyright © 2012 John Wiley & Sons, Ltd.     

eXtended hybridizable discontinuous Galerkin for incompressible flow problems with unfitted meshes and interfaces

The eXtended hybridizable discontinuous Galerkin (X-HDG) method is developed for the solution of Stokes problems with void or material interfaces. X-HDG is a novel method that combines the hybridizable discontinuous Galerkin (HDG) method with an eXtended finite element strategy, resulting in a high-order, unfitted, superconvergent method, with an explicit definition of the interface geometry by means of a level-set function. For elements not cut by the interface, the standard HDG formulation is applied, whereas a modified weak form for the local problem is proposed for cut elements. Heaviside enrichment is considered on cut faces and in cut elements in the case of bimaterial problems. Two-dimensional numerical examples demonstrate that the applicability, accuracy, and superconvergence properties of HDG are inherited in X-HDG, with the freedom of computational meshes that do not fit the interfaces     

Three‐dimensional discontinuous Galerkin elements with plane waves and Lagrange multipliers for the solution of mid‐frequency Helmholtz problems

Recently, a discontinuous Galerkin finite element method with plane wave basis functions and Lagrange multiplier degrees of freedom was proposed for the efficient solution in two dimensions of Helmholtz problems in the mid‐frequency regime. In this paper, this method is extended to three dimensions and several new elements are proposed. Computational results obtained for several wave guide and acoustic scattering model problems demonstrate one to two orders of magnitude solution time improvement over the higher‐order Galerkin method. Copyright © 2005 John Wiley & Sons, Ltd.