Discontinuous Galerkin method coupled with an integral representation for solving the three-dimensional time-harmonic Maxwell equations

In this paper, we present a mathematical and numerical studies of the three-dimensional time-harmonic Maxwell equations. The problem is solved by a discontinuous Galerkin DG method coupled with an integral representation. This study was completed by some numerical tests to justify the effectiveness of the proposed approach. The numerical simulation was done by an iterative solver implemented in FORTRAN.     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

Locally implicit discontinuous Galerkin method for time domain electromagnetics

In the recent years, there has been an increasing interest in discontinuous Galerkin time domain (DGTD) methods for the solution of the unsteady Maxwell equations modeling electromagnetic wave propagation. One of the main features of DGTD methods is their ability to deal with unstructured meshes which are particularly well suited to the discretization of the geometrical details and heterogeneous media that characterize realistic propagation problems. Such DGTD methods most often rely on explicit time integration schemes and lead to block diagonal mass matrices. However, explicit DGTD methods are also constrained by a stability condition that can be very restrictive on highly refined meshes and when the local approximation relies on high order polynomial interpolation. An implicit time integration scheme is a natural way to obtain a time domain method which is unconditionally stable but at the expense of the inversion of a global linear system at each time step. A more viable approach consists of applying an implicit time integration scheme locally in the refined regions of the mesh while preserving an explicit time scheme in the complementary part, resulting in an hybrid explicit–implicit (or locally implicit) time integration strategy. In this paper, we report on our recent efforts towards the development of such a hybrid explicit–implicit DGTD method for solving the time domain Maxwell equations on unstructured simplicial meshes. Numerical experiments for 3D propagation problems in homogeneous and heterogeneous media illustrate the possibilities of the method for simulations involving locally refined meshes.     

Implicit DG Method for Time Domain Maxwell's Equations Involving Metamaterials

An implicit discontinuous Galerkin method is introduced to solve the time-domain Maxwell's equations in metamaterials. The Maxwell's equations in metamaterials are represented by integral-differential equations. Our scheme is based on discontinuous Galerkin method in spatial domain and Crank-Nicolson method in temporal domain. The fully discrete numerical scheme is proved to be unconditionally stable. When polynomial of degree at most $p$ is used for spatial approximation, our scheme is verified to converge at a rate of $\mathcal{O}(τ^2+h^{p+1/2})$. Numerical results in both 2D and 3D are provided to validate our theoretical prediction.     

A discontinuous Galerkin method for the Vlasov–Poisson system

A discontinuous Galerkin method for approximating the Vlasov–Poisson system of equations describing the time evolution of a collisionless plasma is proposed. The method is mass conservative and, in the case that piecewise constant functions are used as a basis, the method preserves the positivity of the electron distribution function and weakly enforces continuity of the electric field through mesh interfaces and boundary conditions. The performance of the method is investigated by computing several examples and error estimates of the approximation are stated. In particular, computed results are benchmarked against established theoretical results for linear advection and the phenomenon of linear Landau damping for both the Maxwell and Lorentz distributions. Moreover, two nonlinear problems are considered: nonlinear Landau damping and a version of the two-stream instability are computed. For the latter, fine scale details of the resulting long-time BGK-like state are presented. Conservation laws are examined and various comparisons to theory are made. The results obtained demonstrate that the discontinuous Galerkin method is a viable option for integrating the Vlasov–Poisson system.     

高阶间断伽辽金时域有限元方法分析三维谐振腔

从一阶麦克斯韦旋度方程出发,研究一种区域分解时域有限元目的——高阶间断伽辽金时域有限元目的.其中对时间的离散采用Crank-Nicolson差分格式,电场和磁场采用相同阶数的高阶矢量基函数展开.分析三维谐振腔问题,数值结果表明,目的 中时间步长的选取可以摆脱CFL稳定性条件的限制;此外,与基于常用Whitney矢量基函数的目的 相比,采用高阶矢量基函数可以明显地提高计算精度及计算效率.     

Solving coupled nonlinear Schrodinger equations via a direct discontinuous Galerkin method

In this work,we present the direct discontinuous Galerkin(DDG) method for the one-dimensional coupled nonlinear Schrdinger(CNLS) equation.We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system.The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method.Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.     

Smoothed aggregation multigrid solvers for high-order discontinuous Galerkin methods for elliptic problems

We develop a smoothed aggregation-based algebraic multigrid solver for high-order discontinuous Galerkin discretizations of the Poisson problem. Algebraic multigrid is a popular and effective method for solving the sparse linear systems that arise from discretizing partial differential equations. However, high-order discontinuous Galerkin discretizations have proved challenging for algebraic multigrid. The increasing condition number of the matrix and loss of locality in the matrix stencil as p increases, in addition to the effect of weakly enforced Dirichlet boundary conditions all contribute to the challenging algebraic setting.     

Solving coupled nonlinear Schrödinger equations via a direct discontinuous Galerkin method

In this work, we present the direct discontinuous Galerkin (DDG) method for the one-dimensional coupled nonlinear Schrödinger (CNLS) equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.     

Non-conformal domain decomposition method with second-order transmission conditions for time-harmonic electromagnetics

Non-overlapping domain decomposition (DD) methods provide efficient algorithms for solving time-harmonic Maxwell equations. It has been shown that the convergence of DD algorithms can be improved significantly by using high order transmission conditions. In this paper, we extend a newly developed second-order transmission condition (SOTC), which involves two second-order transverse derivatives, to facilitate fast convergence in the non-conformal DD algorithms. However, the non-conformal nature of the DD methods introduces an additional technical difficulty, which results in poor convergence in many real-life applications. To mitigate the difficulty, a corner-edge penalty method is proposed and implemented in conjunction with the SOTC to obtain truly robust solver performance. Numerical results verify the analysis and demonstrate the effectiveness of the proposed methods on a few model problems. Finally, drastically improved convergence, compared to the conventional Robin transmission condition, was observed for an electrically large problem of practical interest.     

Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field

In this paper, central discontinuous Galerkin methods are developed for solving ideal magnetohydrodynamic (MHD) equations. The methods are based on the original central discontinuous Galerkin methods designed for hyperbolic conservation laws on overlapping meshes, and use different discretization for magnetic induction equations. The resulting schemes carry many features of standard central discontinuous Galerkin methods such as high order accuracy and being free of exact or approximate Riemann solvers. And more importantly, the numerical magnetic field is exactly divergence-free. Such property, desired in reliable simulations of MHD equations, is achieved by first approximating the normal component of the magnetic field through discretizing induction equations on the mesh skeleton, namely, the element interfaces. And then it is followed by an element-by-element divergence-free reconstruction with the matching accuracy. Numerical examples are presented to demonstrate the high order accuracy and the robustness of the schemes.     

Direct discontinuous Galerkin method for the generalized Burgers—Fisher equation

In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge-Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability of the method.     

Numerical upscaling for the eddy-current model with stochastic magnetic materials

This paper deals with the upscaling of the time-harmonic Maxwell equations for heterogeneous media. We analyze the eddy-current approximation of Maxwell’s equations to describe the electric field for heterogeneous, isotropic magnetic materials. The magnetic permeability of the materials is assumed to have random heterogeneities described by a Gaussian random field. We apply the so-called Coarse Graining method to develop a numerical upscaling of the eddy-current model. The upscaling uses filtering and averaging procedures in Fourier space which results in a formulation of the eddy-current model on coarser resolution scales where the influence of sub-scale fluctuations is modeled by effective scale- and space-dependent reluctivity tensors. The effective reluctivity tensors can be obtained by solving local partial differential equations which contain a Laplacian as well as a curl–curl operator. We present a computational method how the equation of the combined operators can be discretized and solved numerically using an extended variational formulation compared to standard discretizations. We compare the results of the numerical upscaling of the eddy-current model with theoretical results of Eberhard [J.P. Eberhard, Upscaling for the time-harmonic Maxwell equations with heterogeneous magnetic materials, Physical Review E 72 (3), (2005)] and obtain a very good agreement.     

New numerical method for solving time-harmonic Maxwell equations with strong singularity

In this paper we introduce a notion of R_ν-generalized solution to time-harmonic Maxwell equations with strong singularity in a 2D nonconvex polygonal domain. We develop a new weighted edge FEM. Results of numerical experiments prove the efficiency of this method.     

Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程

应用自适应LWDG方法求解三维双曲守恒律方程组,与传统的二阶RKDG方法相比,该方法具有计算量小和精度高的特点.给出一种自适应策略,其中均衡折中策略适用于非相容四面体网格.将二维情形下的后验误差指示子推广到三维双曲守恒律方程组中,数值实验证明了方法的有效性.     

A Discontinuous Galerkin Method Based on a BGK Scheme for the Navier-Stokes Equations on Arbitrary Grids

A discontinuous Galerkin Method based on a Bhatnagar-Gross-Krook (BGK) formulation is presented for the solution of the compressible Navier-Stokes equations on arbitrary grids. The idea behind this approach is to combine the robustness of the BGK scheme with the accuracy of the DG methods in an effort to develop a more accurate, efficient, and robust method for numerical simulations of viscous flows in a wide range of flow regimes. Unlike the traditional discontinuous Galerkin methods, where a Local Discontinuous Galerkin (LDG) formulation is usually used to discretize the viscous fluxes in the Navier-Stokes equations, this DG method uses a BGK scheme to compute the fluxes which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous representation in the flux evaluation at a cell interface through a simple hybrid gas distribution function. The developed method is used to compute a variety of viscous flow problems on arbitrary grids. The numerical results obtained by this BGKDG method are extremely promising and encouraging in terms of both accuracy and robustness, indicating its ability and potential to become not just a competitive but simply a superior approach than the current available numerical methods.     

Local discontinuous Galerkin method for solving Burgers and coupled Burgers equations

In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.     

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.     

聚合物充填过程中裹气现象的数值模拟

针对聚合物充填过程中的裹气现象,采用一种有限元（FEM）-间断有限元（DG）耦合算法对其进行数值模拟。对于自由运动界面,采用水平集（Level Set）方法进行捕捉;用XPP（eXtended Pom-Pom）本构模型来描述黏弹性流体的流变行为。采用有限元-间断有限元耦合算法求解统一的流场方程,并采用隐式间断有限元求解XPP本构方程、Level Set及其重新初始化方程。数值结果与文献中的实验结果及模拟结果吻合较好,验证了数值方法的稳定性及准确性。分析带有非规则嵌件型腔内,注射速度及浇口尺寸对裹气现象的影响,裹气容易出现在较高注射速度及较小浇口的情形。     

Discontinuous Galerkin methods in nanophotonics

This article reviews the state of the recently developed discontinuous Galerkin finite element method for the efficient numerical treatment of nanophotonic systems. This approach combines the accurate and flexible spatial discretisation of classical finite elements with efficient time stepping capabilities. We describe in detail the underlying principles of the discontinuous Galerkin technique and its application to the simulation of complex nanophotonic structures. In addition, formulations for both time‐ and frequency‐domain solvers are provided and specific advantages and limitations of the technique are discussed. The potential of the discontinuous Galerkin approach is illustrated by modelling and simulating several experimentally relevant systems.