Unstructured Grid-Based Discontinuous Galerkin Method for Broadband Electromagnetic Simulations

A parallel, unstructured, high-order discontinuous Galerkin method is developed for the time-dependent Maxwell's equations, using simple monomial polynomials for spatial discretization and a fourth-order Runge–Kutta scheme for time marching. Scattering results for a number of validation cases are computed employing polynomials of up to third order. Accurate solutions are obtained on coarse meshes and grid convergence is achieved, demonstrating the capabilities of the scheme for time-domain electromagnetic wave scattering simulations.     

Arbitrary High-Order Discontinuous Galerkin Schemes for the Magnetohydrodynamic Equations

In this paper, we propose a discontinuous Galerkin scheme with arbitrary order of accuracy in space and time for the magnetohydrodynamic equations. It is based on the Arbitrary order using DERivatives (ADER) methodology: the high order time approximation is obtained by a Taylor expansion in time. In this expansion all the time derivatives are replaced by space derivatives via the Cauchy-Kovalevskaya procedure. We propose an efficient algorithm of the Cauchy-Kovalevskaya procedure in the case of the three-dimensional magneto-hydrodynamic (MHD) equations. Parallel to the time derivatives of the conservative variables the time derivatives of the fluxes are calculated. This enables the analytic time integration of the volume integral as well as that of the surface integral of the fluxes through the grid cell interfaces which occur in the discrete equations. At the cell interfaces the fluxes and all their derivatives may jump. Following the finite volume ADER approach the break up of all these jumps into the different waves are taken into account to get proper values of the fluxes at the grid cell interfaces. The approach under considerations is directly based on the expansion of the flux in time in which the leading order term may be any numerical flux calculation for the MHD-equation. Numerical convergence results for these equations up to 7th order of accuracy in space and time are shown.     

A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations

In this paper, we develop a class of high order conservative semi-Lagrangian (SL) discontinuous Galerkin methods for solving multi-dimensional linear transport equations. The methods rely on a characteristic Galerkin weak formulation, leading to \(L^2\) stable discretizations for linear problems. Unlike many existing SL methods, the high order accuracy and mass conservation of the proposed methods are realized in a non-splitting manner. Thus, the detrimental splitting error, which is known to significantly contaminate long term transport simulations, will be not incurred. One key ingredient in the scheme formulation, borrowed from CSLAM (Lauritzen et al. in J Comput Phys 229(5):1401–1424, 2010), is the use of Green’s theorem which allows us to convert volume integrals into a set of line integrals. The resulting line integrals are much easier to approximate with high order accuracy, hence facilitating the implementation. Another novel ingredient is the construction of quadratic curves in approximating sides of upstream cell, leading to quadratic-curved quadrilateral upstream cells. Formal third order accuracy is obtained by such a construction. The desired positivity-preserving property is further attained by incorporating a high order bound-preserving filter. To assess the performance of the proposed methods, we test and compare the numerical schemes with a variety of configurations for solving several benchmark transport problems with both smooth and nonsmooth solutions. The efficiency and efficacy are numerically verified.     

Arbitrary High-Order Explicit Hybridizable Discontinuous Galerkin Methods for the Acoustic Wave Equation

We propose a new formulation of explicit time integration for the hybridizable discontinuous Galerkin (HDG) method in the context of the acoustic wave equation based on the arbitrary derivative approach. The method is of arbitrary high order in space and time without restrictions such as the Butcher barrier for Runge–Kutta methods. To maintain the superconvergence property characteristic for HDG spatial discretizations, a special reconstruction step is developed, which is complemented by an adjoint consistency analysis. For a given time step size, this new method is twice as fast compared to a low-storage Runge–Kutta scheme of order four with five stages at polynomial degrees between two and four. Several numerical examples are performed to demonstrate the convergence properties, reveal dispersion and dissipation errors, and show solution behavior in the presence of material discontinuities. Also, we present the combination of local time stepping with h-adaptivity on three-dimensional meshes with curved elements.     

A New Nonsymmetric Discontinuous Galerkin Method for Time Dependent Convection Diffusion Equations

We propose a discontinuous Galerkin finite element method for convection diffusion equations that involves a new methodology handling the diffusion term. Test function derivative numerical flux term is introduced in the scheme formulation to balance the solution derivative numerical flux term. The scheme has a nonsymmetric structure. For general nonlinear diffusion equations, nonlinear stability of the numerical solution is obtained. Optimal kth order error estimate under energy norm is proved for linear diffusion problems with piecewise P ^k polynomial approximations. Numerical examples under one-dimensional and two-dimensional settings are carried out. Optimal (k+1)th order of accuracy with P ^k polynomial approximations is obtained on uniform and nonuniform meshes. Compared to the Baumann-Oden method and the NIPG method, the optimal convergence is recovered for even order P ^k polynomial approximations.     

An hp-Analysis of the Local Discontinuous Galerkin Method for Diffusion Problems

We present an hp-error analysis of the local discontinuous Galerkin method for diffusion problems, considering unstructured meshes with hanging nodes and two- and three-dimensional domains. Our estimates are optimal in the meshsize h and slightly suboptimal in the polynomial approximation order p. Optimality in p is achieved for matching grids and polynomial boundary conditions.     

An A Posteriori Error Estimate for the Local Discontinuous Galerkin Method

In this work an a posteriori global error estimate for the Local Discontinuous Galerkin (LDG) applied to a linear second order elliptic problem is analyzed. Using a mixed formulation, an upper bound of the error in the primal variable is derived from explicit computations. Finally, a local adaptive scheme based on explicit error estimators is studied numerically using one dimensional problems.This revised version was published online in July 2005 with corrected volume and issue numbers.     

Discontinuous Galerkin Method for Investigating Ice Strength

This paper discusses the numerical modeling of various ice-strength measurement experiments, including uniaxial compression and bending, and it also compares the data obtained by field and numerical experiments. Numerical simulation is based on a dynamic system of continuum mechanics equations with ice considered as an elasto-plastic medium with brittle and crushing fracture dynamic criteria. The simulation software developed by the authors is based on the discontinuous Galerkin method and runs on high-performance systems with a distributed memory. Estimating the explicit values used by the mathematical models poses a major problem because some of them cannot be directly measured in field experiments due to the multiple interferences of physical processes. In practice, it is only possible to directly measure their total influence. However, this problem can be solved by comparing the numerical experiment with the field data. As a result of this work, the elasto-plastic ice model is verified and some missing physical properties are obtained by the numerical experiments.     

A Uniform Additive Schwarz Preconditioner for High-Order Discontinuous Galerkin Approximations of Elliptic Problems

In this paper we design and analyze a uniform preconditioner for a class of high-order Discontinuous Galerkin schemes. The preconditioner is based on a space splitting involving the high-order conforming subspace and results from the interpretation of the problem as a nearly-singular problem. We show that the proposed preconditioner exhibits spectral bounds that are uniform with respect to the discretization parameters, i.e., the mesh size, the polynomial degree and the penalization coefficient. The theoretical estimates obtained are supported by numerical tests.     

Discontinuous Galerkin Method for Linear Free-Surface Gravity Waves

In this paper, we discuss a discontinuous Galerkin finite (DG) element method for linear free surface gravity waves. We prove that the algorithm is unconditionally stable and does not require additional smoothing or artificial viscosity terms in the free surface boundary condition to prevent numerical instabilities on a non-uniform mesh. A detailed error analysis of the full time-dependent algorithm is given, showing that the error in the wave height and velocity potential in the L²-norm is in both cases of optimal order and proportional to O(Δt²+h^p+1), without the need for a separate velocity reconstruction, with p the polynomial order, h the mesh size and Δt the time step. The error analysis is confirmed with numerical simulations. In addition, a Fourier analysis of the fully discrete scheme is conducted which shows the dependence of the frequency error and wave dissipation on the time step and mesh size. The algebraic equations for the DG discretization are derived in a way suitable for an unstructured mesh and result in a symmetric positive definite linear system. The algorithm is demonstrated on a number of model problems, including a wave maker, for discretizations with accuracy ranging from second to fourth order.This revised version was published online in July 2005 with corrected volume and issue numbers.     

多物理数值模拟中一种有效的并行耦合方法

在实现多物理并行数值耦合模拟中,需要处理多个物理过程之间网格、并行区域分解的差异.针对该同题,该文基于三维流体力学与激光传播耦合的并行数值模拟,提出了一种实用的并行耦合方法:引入辅助状态将本地插值与通信相分离;构建并行耦合图并定义主导属性,以确定过程间传输的最小数据集合;提供并行数据重分配算法来完成通信.并行数值结果表明:该方法是有效的,在64台处理机上使整体程序获得50.07的加速比.     

Dispersion and Dissipation Error in High-Order Runge-Kutta Discontinuous Galerkin Discretisations of the Maxwell Equations

Different time-stepping methods for a nodal high-order discontinuous Galerkin discretisation of the Maxwell equations are discussed. A comparison between the most popular choices of Runge-Kutta (RK) methods is made from the point of view of accuracy and computational work. By choosing the strong-stability-preserving Runge-Kutta (SSP-RK) time-integration method of order consistent with the polynomial order of the spatial discretisation, better accuracy can be attained compared with fixed-order schemes. Moreover, this comes without a significant increase in the computational work. A numerical Fourier analysis is performed for this Runge-Kutta discontinuous Galerkin (RKDG) discretisation to gain insight into the dispersion and dissipation properties of the fully discrete scheme. The analysis is carried out on both the one-dimensional and the two-dimensional fully discrete schemes and, in the latter case, on uniform as well as on non-uniform meshes. It also provides practical information on the convergence of the dissipation and dispersion error up to polynomial order 10 for the one-dimensional fully discrete scheme.     

A Discontinuous Galerkin Method for Three-Dimensional Shallow Water Equations

We describe the application of a local discontinuous Galerkin method to the numerical solution of the three-dimensional shallow water equations. The shallow water equations are used to model surface water flows where the hydrostatic pressure assumption is valid. The authors recently developed a DG\linebreak method for the depth-integrated shallow water equations. The method described here is an extension of these ideas to non-depth-integrated models. The method and its implementation are discussed, followed by numerical examples on several test problems.This revised version was published online in July 2005 with corrected volume and issue numbers.     

A Hybrid Staggered Discontinuous Galerkin Method for KdV Equations

A hybrid staggered discontinuous Galerkin method is developed for the Korteweg–de Vries equation. The equation is written into a system of first order equations by introducing auxiliary variables. Two sets of finite element functions are introduced to approximate the solution and the auxiliary variables. The staggered continuity of the two finite element function spaces gives a natural flux condition and trace value on the element boundaries in the derivation of Galerkin approximation. On the other hand, to deal with the third order derivative term an hybridization idea is used and additional flux unknowns are introduced. The auxiliary variables can be eliminated in each element and the resulting algebraic system on the solution and the additional flux unknowns is solved. Stability of the semi discrete form is proven for various boundary conditions. Numerical results present the optimal order of \(L^2\)-errors of the proposed method for a given polynomial order.     

A Discontinuous Galerkin Method for the Subjective Surfaces Problem

The work formulates and evaluates the local discontinuous Galerkin method for the subjective surfaces problem based on the curvature driven level set equation. A new mixed formulation simplifying the treatment of nonlinearities is proposed. The numerical algorithm is evaluated using several artificial and realistic test cases.     

Moving Mesh Discontinuous Galerkin Method for Hyperbolic Conservation Laws

In this paper, a moving mesh discontinuous Galerkin (DG) method is developed to solve the nonlinear conservation laws. In the mesh adaptation part, two issues have received much attention. One is about the construction of the monitor function which is used to guide the mesh redistribution. In this study, a heuristic posteriori error estimator is used in constructing the monitor function. The second issue is concerned with the solution interpolation which is used to interpolates the numerical solution from the old mesh to the updated mesh. This is done by using a scheme that mimics the DG method for linear conservation laws. Appropriate limiters are used on seriously distorted meshes generated by the moving mesh approach to suppress the numerical oscillations. Numerical results are provided to show the efficiency of the proposed moving mesh DG method.