Discontinuous Galerkin Finite Element Research Articles

On Error Estimates of a discontinuous Galerkin Method of the Boussinesq System of Equations

Abstract In this paper, we propose and analyze a discontinuous Galerkin finite element method for solving the transient Boussinesq incompressible heat conducting fluid flow equations. This method utilizes an upwind approach to handle the nonlinear convective terms effectively. We discuss new a priori bounds for the semidiscrete discontinuous Galerkin approximations. Furthermore, we establish optimal a priori error estimates for the semidiscrete discontinuous Galerkin velocity approximation in L 2 \mathbf{L}^{2} and energy norms, the temperature approximation in L 2 L^{2} and energy norms and pressure approximation in L 2 L^{2} -norm for t > 0 t>0 . Additionally, under the smallness assumption on the data, we prove uniform in time error estimates. We also consider a backward Euler scheme for full discretization and derive fully discrete error estimates. Finally, we provide numerical examples to support the theoretical conclusions.

Efficient High-Order Space-Angle-Energy Polytopic Discontinuous Galerkin Finite Element Methods for Linear Boltzmann Transport

We introduce an hp-version discontinuous Galerkin finite element method (DGFEM) for the linear Boltzmann transport problem. A key feature of this new method is that, while offering arbitrary order convergence rates, it may be implemented in an almost identical form to standard multigroup discrete ordinates methods, meaning that solutions can be computed efficiently with high accuracy and in parallel within existing software. This method provides a unified discretisation of the space, angle, and energy domains of the underlying integro-differential equation and naturally incorporates both local mesh and local polynomial degree variation within each of these computational domains. Moreover, general polytopic elements can be handled by the method, enabling efficient discretisations of problems posed on complicated spatial geometries. We study the stability and hp-version a priori error analysis of the proposed method, by deriving suitable hp-approximation estimates together with a novel inf-sup bound. Numerical experiments highlighting the performance of the method for both polyenergetic and monoenergetic problems are presented.

A discontinuous Galerkin method for a coupled Stokes–Biot problem

The interaction between fluid and poroelastic structure is a complex problem that couples the Stokes equations with the Biot system. In this work, we rewrite the poroelastic equations using three fields (displacement, fluid pressure, and total pressure) to account for locking phenomenon in fluid–structure interaction and the oscillation of numerical solutions. Then, in order to solve the coupled system, we use the weighted discontinuous Galerkin finite element method for spatial discretisation and propose a semi-discrete scheme. In addition, we propose a fully discrete scheme using the backward Euler’s method for temporal discretisation. In the framework of the Galerkin approximation, the stability and errors of semi- and fully discrete schemes are analysed using the theory of differential algebraic equations and weak compactness arguments. Through numerical experiments, we validate the theoretical rate of convergence, explain the behaviour of the method in modeling the interaction between surface and groundwater hydrological systems, investigate the robustness of the method with respect to physical parameters, and simulate channel filtering. In particular, the numerical results are free of locking features, reducing non-physical numerical oscillations.

Operator-splitting finite element method for solving the radiative transfer equation

AbstractAn operator-splitting finite element scheme for the time-dependent radiative transfer equation is presented in this paper. The streamline upwind Petrov-Galerkin finite element method and discontinuous Galerkin finite element method are used for the spatial-angular discretization of the radiative transfer equation, whereas the backward Euler scheme is used for temporal discretization. Error analysis of the proposed numerical scheme for the fully discrete radiative transfer equation is presented. The stability and convergence estimates for the fully discrete problem are derived. Moreover, an operator-splitting algorithm for the numerical simulation of high-dimensional equations is also presented. The validity of the derived estimates and implementation is illustrated with suitable numerical experiments.

Discontinuous Galerkin finite element method for solving non-linear model of gradient elution chromatography

Discontinuous Galerkin finite element method for solving non-linear model of gradient elution chromatography

Discontinuous Galerkin finite element method for dynamic viscoelasticity models of power‐law type

Discontinuous Galerkin finite element method for dynamic viscoelasticity models of power‐law type

Iterative solution methods for high-order/hp–DGFEM approximation of the linear Boltzmann transport equation

In this article we consider the iterative solution of the linear system of equations arising from the discretisation of the poly-energetic linear Boltzmann transport equation using a high-order/hp–version discontinuous Galerkin finite element approximation in space, angle, and energy. In particular, we develop preconditioned Richardson iterations which may be understood as generalisations of source iteration in the mono-energetic setting, and derive computable a posteriori bounds for the solver error incurred due to inexact linear algebra, measured in a relevant problem-specific norm. We prove that the convergence of the resulting schemes and a posteriori solver error estimates are independent of the mesh size h and polynomial degree p. We also discuss how the poly-energetic Richardson iteration may be employed as a preconditioner for the generalised minimal residual (GMRES) method. Furthermore, we show that standard implementations of GMRES based on minimising the Euclidean norm of the residual vector can be utilized to yield computable a posteriori solver error estimates at each iteration, through judicious selections of left- and right-preconditioners for the original linear system. The effectiveness of poly-energetic source iteration and preconditioned GMRES, as well as their respective a posteriori solver error estimates, is demonstrated through numerical examples arising in the modelling of photon transport. While the convergence of poly-energetic source iteration is independent of h and p, we observe that the number of iterations required to attain convergence when employing GMRES only depends mildly on h and p. Moreover, this latter approach is highly effective in the low energy regime.

Port-Hamiltonian discontinuous Galerkin finite element methods

Abstract A port-Hamiltonian (pH) system formulation is a geometrical notion used to formulate conservation laws for various physical systems. The distributed parameter port-Hamiltonian formulation models infinite dimensional Hamiltonian dynamical systems that have a nonzero energy flow through the boundaries. In this paper, we propose a novel framework for discontinuous Galerkin (DG) discretizations of pH-systems. Linking DG methods with pH-systems gives rise to compatible structure preserving semidiscrete finite element discretizations along with flexibility in terms of geometry and function spaces of the variables involved. Moreover, the port-Hamiltonian formulation makes boundary ports explicit, which makes the choice of structure and power preserving numerical fluxes easier. We state the Discontinuous Finite Element Stokes–Dirac structure with a power preserving coupling between elements, which provides the mathematical framework for a large class of pH discontinuous Galerkin discretizations. We also provide an a priori error analysis for the port-Hamiltonian discontinuous Galerkin Finite Element Method (pH-DGFEM). The port-Hamiltonian discontinuous Galerkin finite element method is demonstrated for the scalar wave equation showing optimal rates of convergence.

Quadratic Discontinuous Galerkin Finite Element Methods for the Unilateral Contact Problem

Abstract In this article, we employ discontinuous Galerkin methods for the finite element approximation of the frictionless unilateral contact problem using quadratic finite elements over simplicial triangulation. We first develop a posteriori error estimates in the energy norm wherein, the reliability and efficiency of the proposed a posteriori error estimator is addressed. The suitable construction of the discrete Lagrange multiplier 𝝀 𝒉 {\boldsymbol{\lambda_{h}}} and some intermediate operators play a key role in developing a posteriori error analysis. Further, we establish an optimal a priori error estimates under the appropriate regularity assumption on the exact solution 𝒖 {\boldsymbol{u}} . Numerical results presented on uniform and adaptive meshes illustrate and confirm the theoretical findings.

Capturing Local Temperature Evolution during Additive Manufacturing through Fourier Neural Operators

Abstract High-fidelity, data-driven models that can quickly simulate thermal behavior during additive manufacturing (AM) are crucial for improving the performance of AM technologies in multiple areas, such as part design, process planning, monitoring, and control. However, the complexities of part geometries make it challenging for current models to maintain high accuracy across a wide range of geometries. Additionally, many models report a low mean square error (MSE) across the entire domain of a part. However, in each time step, most areas of the domain do not experience significant changes in temperature, except for the regions near recent depositions. Therefore, the MSE-based fidelity measurement of the models may be overestimated. This paper presents a data-driven model that uses the Fourier Neural Operator to capture the local temperature evolution during the AM process. Beside MSE, the model is also evaluated using the R2 metric, which places great weight on the regions where the temperature changes significantly than MSE. The model was trained and tested on numerical simulations based on the Discontinuous Galerkin Finite Element Method for the Direct Energy Deposition AM process. The results shows that the model maintains 0.983-0.999 R2 over geometries not included in the training data, which is higher than Convolutional Neural Networks and Graph Convolutional Neural Networks we implemented, the two widely used architectures in data-driven predictive modeling.

Discontinuous Galerkin method for the coupled dual-porosity-Brinkman model

In this paper we use the weighted discontinuous Galerkin finite element method to solve the dual-porosity-Brinkman coupling model. The dual-porosity model describes fluid flow in a dual-porosity medium composed of matrix and microfractures, while the Brinkman equation describes fluid flow in hydraulic fractures or conduits containing “debris filling” or “proppant filling”. The two models are coupled by the four physically effective interface conditions, which are constructed based on the fundamental properties of the dual-porosity model and the Brinkman-Darcy model. We provide a weighted discontinuous discrete scheme for the model and derive the optimal error estimates. Through numerical experiments, we verify the properties of the model and the advantages of the numerical method, such as the optimal convergence rate of the numerical solution and the simulation of the flow characteristics around hydraulic fractures and conduits.

Discontinuous Galerkin Method of the First Order Parabolic Differential Equation

This paper is concerned with the mathematical study of the discontinuous Galerkin (DG) finite element method for the parabolic differential equation. The DG method is a vital numerical method with much mass compensation and more flexible meshing than other methods. In this study, we give a general introduction and discuss about the discontinuous Galerkin Method of first order parabolic problem. The parabolic problem satisfies the condition of the existence and uniqueness of DG solution. The error analysis of this problem is also established. The main goal of this study is to theoretically explore the convergence of the solution of the above methods and show the validity of the results.

Robust convergence result of discontinuous Galerkin stabilization method for two-dimensional reaction–diffusion equation with discontinuous source term

A reaction–diffusion problem with discontinuous source term and Dirichlet's boundary conditions on the unit square is considered in this paper. The proposed problem has been discretized using a combination of standard Galerkin finite element method (FEM) and non-symmetric discontinuous Galerkin finite element method with an interior penalty (NIPG) with bilinear elements. Layer adapted mesh of Shishkin type has been utilized to discretize the domain. Standard Galerkin FEM is applied on the layer part of the domain where the domain is dense enough and NIPG is applied to the outside layer part. By means of special choice of discontinuity-penalization parameters, the scheme is proved to be uniformly convergent of order O ( ε 1 / 4 N − 1 + N − 1 ln ⁡ N ) . Numerical tests are carried out in support of theoretical findings.

An Analytical Solution for Variable Viscosity Flow in Fractured Media: Development and Comparative Analysis With Numerical Simulations

AbstractExplicit fracture models often use analytical solutions for predicting flow in fractured media, usually assuming uniform fluid viscosity for simplicity. This assumption, however, can be inaccurate as fluid viscosity varies due to factors like composition, temperature, and dissolved substances. Our study, recognizing these discrepancies, abandons this uniform viscosity assumption for a more realistic model of variable viscosity flow, focusing on viscous displacement scenarios. This includes instances like injecting viscous surfactants for hydrocarbon recovery in fractured reservoirs or soil decontamination. This presents a significant challenge, enhancing our understanding of transport within fractures, mainly governed by advection. Our study centers on a low‐permeability rock matrix intersected by two fractures with variable apertures. We employ two methods: an analytical approach with a new solution and numerical simulations with two distinct in‐house codes, discretizing both the rock matrix and fractures with two‐dimensional triangular elements. The first code uses a Discontinuous Galerkin finite element method, while the second utilizes a finite‐volume method, allowing a comprehensive comparison of solutions. Additionally, we investigate parameter identifiability, like fracture apertures and viscosity ratios, using breakthrough curves from our analytical solution, applying the Markov Chain Monte Carlo technique.

Finite element modeling of infrasound resonance in an underground tunnel structure

Recent observations of acoustic/infrasound signals emanating from an underground mine in central Utah show a spectral signature consisting of several dominant harmonics whose amplitude decay rate increases with frequency. This spectral character is consistent with a model in which the observed signals are generated from acoustic resonance within the mine tunnels. I present models of these signals computed using a discontinuous-Galerkin finite element formulation. The mesh geometry of these models is based on the known geometry of the mine tunnels obtained using traditional survey methods. Explosion source locations in the mine are known from previous analysis of seismic signals and operator logs and are associated with observations of infrasound resonance at multiple far-field stations. I explore the sensitivity of the modeled signals to changes in the source spectrum, boundary conditions within in the mine and perturbations to tunnel geometry. I discuss the possible application of my results to volcanic monitoring, industrial monitoring of underground structures, and the possible determination of underground tunnel geometries using far-field observations of acoustic resonance.

A fully implicit edge/face centered discontinuous Galerkin / mixed finite element scheme for the advection-dispersion equation

Mixed Finite Element (MFE) method is a robust numerical technique for solving elliptic and parabolic partial differential equations (PDEs). However, MFE can generate solutions with strong unphysical oscillations and/or large numerical diffusion for hyperbolic type PDEs. For its part, Discontinuous Galerkin (DG) finite element method is well adapted to solve hyperbolic systems and can accurately reproduce solutions involving sharp fronts. Therefore, the combination of DG and MFE is a good strategy for solving hyperbolic/parabolic problems such as advection – diffusion/dispersion equations. The classical formulation of the two methods is based on operator and time splitting allowing for separate solutions to advection with an explicit scheme and to dispersion with an implicit scheme. However, this kind of approach has the following drawbacks: (i) it lacks efficiency, as two systems with different unknowns are solved at each time step, (ii) it induces errors generated by the splitting, (iii) it can be CPU wise-expensive because of the CFL constraint, and (iv) it cannot be employed for steady-state transport simulations.To overcome these difficulties, we develop in this work a fully implicit edge/face centered DG-MFE formulation where the two methods share the same unknowns. In this formulation, the DG method is developed on lumping regions associated with the mesh edges/faces instead of mesh elements. Thus, the traces of concentration at mesh edges/faces, which are the Degrees Of Freedom (DOF) of the hybrid-MFE, are also part of the DOFs of the DG. The temporal discretization is based on the Crank–Nicolson method for both advection and dispersion. Numerical tests are performed to validate the new scheme by comparison against an analytical solution and to show its ability to handle steady-state transport simulations.The procedure is developed for 2D triangular meshes but can easily be extended to other 2D and 3D shape elements.

A discontinuous Galerkin finite element method for the Oldroyd model of order one

In this work, we analyze a discontinuous Galerkin finite element method for the equations of motion that arise in the 2D Oldroyd model of order one. We investigate the existence and uniqueness of semidiscrete discontinuous solutions, as well as the consistency of the scheme. We derive new a priori and regularity results for the discrete solution and establish optimal error estimates in ‐norm in time and energy norm in space for the velocity and ‐norm in both time and space for the pressure. Uniform estimates are derived for sufficiently small data. We next apply the backward Euler method to the semidiscrete formulation and establish optimal fully discrete error estimates. At the end, we conduct numerical experiments to support our theoretical results and analyze the findings.

A Discontinuous Galerkin–Finite Element Method for the Nonlinear Unsteady Burning Rate Responses of Solid Propellants

The unsteady combustion of solid propellants under oscillating environments is the key to understanding the combustion instability inside solid rocket motors. The discontinuous Galerkin–finite element method (DG-FEM) is introduced to provide an efficient yet flexible numerical platform to investigate the combustion dynamics of solid propellants. The algorithm is developed for the classical unsteady model, the Zel’dovich–Novozhilov model. It is then validated based on a special analytical solution. The DG-FEM algorithm is then compared with the classical spectral method based on Laguerre polynomials. It is shown that the DG-FEM works more efficiently than the traditional spectral method, providing a more accurate solution with a lower computational cost.

A Stokes–Darcy–Darcy model and its discontinuous Galerkin method on polytopic grids

In this paper, we propose and numerically solve a model considering confined flow in fractured porous media coupled with free flow in conduits. Such a situation arises, for example, for fluid flows in fractured vuggy carbonate reservoirs and caves/conduits, where a matrix and fractures are the main spaces for fluid storage and flow. The flow in caves and conduits is governed by the Stokes equation; the flow in fractured porous media, which consists of both a matrix and fractures, is described by a coupled Darcy–Darcy (bulk–fracture) model. Then, the two models are coupled through suitable (physically consistent) interface conditions on the interface between fractured porous media and caves/conduits, playing a key role in a physically faithful simulation that achieves high accuracy. The construction of all the interface conditions of the coupled model is based on the fundamental properties of the traditional bulk–fracture model and the well-known Stokes–Darcy model. The weak formulation is derived for the proposed model, and the well-posedness of the model is theoretically analyzed. We also propose a polyhedral discontinuous Galerkin finite element approximation for the Stokes flow in conduits and Darcy flow in the porous matrix, which is coupled with a conforming finite element scheme for the flow in the fracture. Optimal error estimates in a suitable energy norm are obtained. Two-dimensional numerical experiments are conducted to validate the proposed model and to demonstrate the features of both the model and the numerical method, such as the optimal convergence rate of the numerical solution; the detail flow characteristics around fractures, the matrix, and conduits; and the applicability to real-world problems.

Numerical simulation of CO2-enhanced oil recovery in fractured shale reservoirs using discontinuous and continuous Galerkin finite element methods

Introduction: This study explores the potential of enhancing shale oil recovery and reducing CO2 emissions through CO2 injection in fractured shale reservoirs. The importance of this approach lies in its dual benefit: improving oil extraction efficiency and addressing environmental concerns associated with CO2 emissions.Method: We employed a discrete fracture-matrix model to simulate CO2 flooding in fractured shale reservoirs, utilizing both discontinuous Galerkin (DG) and continuous Galerkin (CG) finite element methods. The DG-CG FEM’s accuracy was validated against the McWhorter problem, ensuring the reliability of the simulation results. Our model also considered various factors, including reservoir heterogeneity, fracture permeability, CO2 injection volume, and gas injection patterns, to analyze their impact on shale oil recovery.Result: Our simulations revealed that fractured reservoirs significantly enhance shale oil production efficiency compared to homogeneous reservoirs, with an approximate 48.9% increase in production. A notable increase in shale oil production, by 15.8%, was observed when fracture permeability was increased by two orders of magnitude. Additionally, a fourfold increase in CO2 injection rate resulted in a 31.5% rise in shale oil production. Implementing a step-by-step reduction in injection volume while maintaining the total CO2 injection constant proved to be more effective than constant-rate injections.Discussion: The study demonstrates the effectiveness of CO2 flooding in fractured shale reservoirs for enhancing shale oil recovery.