In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finit...In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.展开更多
A novel high-order three-dimensional(3-D)discontinuous Galerkin time domain(DGTD)method based on a normalized formulation of Maxwell′s equations is developed for modeling and simulating silicon-on-insulator(SOI)thin-...A novel high-order three-dimensional(3-D)discontinuous Galerkin time domain(DGTD)method based on a normalized formulation of Maxwell′s equations is developed for modeling and simulating silicon-on-insulator(SOI)thin-ridge waveguide.The DGTD method employs unstructured meshes and piecewise high-order polynomials for spatial discretization,and Runge-Kutta methods for time integration.It is found that the numerical results of the leakage loss of SOI thin-ridge waveguide agree well with those of analytical solutions,which proves that the proposed method is an ideal tool for the quantitative analysis for SOI thin-ridge waveguide.展开更多
Several major challenges need to be faced for efficient transient multiscale electromagnetic simulations, such as flex- ible and robust geometric modeling schemes, efficient and stable time-stepping algorithms, etc. F...Several major challenges need to be faced for efficient transient multiscale electromagnetic simulations, such as flex- ible and robust geometric modeling schemes, efficient and stable time-stepping algorithms, etc. Fortunately, because of the versatile choices of spatial discretization and temporal integration, a discontinuous Galerkin time-domain (DGTD) method can be a very promising method of solving transient multiscale electromagnetic problems. In this paper, we present the application of a leap-frog DGTD method to the analyzing of the multiscale electromagnetic scattering problems. The uniaxial perfect matching layer (UPML) truncation of the computational domain is discussed and formulated in the leap-frog DGTD context. Numerical validations are performed in the challenging test cases demonstrating the accuracy and effectiveness of the method in solving transient multiscale electromagnetic problems compared with those of other numerical methods.展开更多
A high-order leap-frog based non-dissipative discontinuous Galerkin time-domain method for solving Maxwell's equations is introduced and analyzed.The proposed method combines a centered approximation for the evalu...A high-order leap-frog based non-dissipative discontinuous Galerkin time-domain method for solving Maxwell's equations is introduced and analyzed.The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements,with a Nth-order leap-frog time scheme.Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes.The method is proved to be stable under some CFL-like condition on the time step.The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided.Numerical experiments with high-order elements show the potential of the method.展开更多
A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.Th...A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.The staggered DG scheme defines the discrete pressure on the primal triangular mesh,while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh.In this paper,a new pair of equal-order-interpolation velocity-pressure finite elements is proposed.On the primary triangular mesh(the pressure elements),the basis functions are piecewise polynomials of degree N and are allowed to jump on the boundaries of each triangle.On the dual mesh instead(the velocity elements),the basis functions consist in the union of piecewise polynomials of degree N on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries,while they are continuous inside.In other words,the basis functions on the dual mesh arc built by continuous finite elements on the subtriangles.This choice allows the construction of an efficient,quadrature-free and memory saving algorithm.In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations,the arbitrary high order of accuracy in time is achieved through tire use of time-dependent test and basis functions,in combination with simple and efficient Picard iterations.Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes,not only from a computer memory point of view,but also concerning the computational time.展开更多
For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numeric...For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.展开更多
A Krylov space based time domain method for wave propagation problems is introduced. The proposed method uses the Arnoldi algorithm to obtain broad-band frequency domain solutions. This method is especially advantageo...A Krylov space based time domain method for wave propagation problems is introduced. The proposed method uses the Arnoldi algorithm to obtain broad-band frequency domain solutions. This method is especially advantageous in cases where slow convergence is observed when using traditional time domain methods. The efficiency of the method is examined in several test cases to show its fast convergence in such problems.展开更多
In this paper,we first present the optimal error estimates of the semi-discrete ultra-weak discontinuous Galerkin method for solving one-dimensional linear convection-diffusion equations.Then,coupling with a kind of R...In this paper,we first present the optimal error estimates of the semi-discrete ultra-weak discontinuous Galerkin method for solving one-dimensional linear convection-diffusion equations.Then,coupling with a kind of Runge-Kutta type implicit-explicit time discretization which treats the convection term explicitly and the diffusion term implicitly,we analyze the stability and error estimates of the corresponding fully discrete schemes.The fully discrete schemes are proved to be stable if the time-stepτ≤τ0,whereτ0 is a constant independent of the mesh-size h.Furthermore,by the aid of a special projection and a careful estimate for the convection term,the optimal error estimate is also obtained for the third order fully discrete scheme.Numerical experiments are displayed to verify the theoretical results.展开更多
In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restric...In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restriction of explicit time marching methods,we employ a class of exponential time differencing(ETD)schemes for time integration,which is based on a linear convex splitting principle.Numerical experiments of the accuracy and long time simulations are given to show the efficiency and capability of the proposed numerical schemes.展开更多
时域离散伽辽金法(Discontinuous Galerkin Time Domain,DGTD)同时具有时域有限元算法(Finite Element Time Domain,FETD)非结构网格剖分和时域有限差分算法(Finite Difference Time Domain,FDTD)显式迭代的优点,是一种非常有前途的电...时域离散伽辽金法(Discontinuous Galerkin Time Domain,DGTD)同时具有时域有限元算法(Finite Element Time Domain,FETD)非结构网格剖分和时域有限差分算法(Finite Difference Time Domain,FDTD)显式迭代的优点,是一种非常有前途的电磁计算方法,该文首先描述了基于矢量基函数的时域离散伽辽金法的基本原理。然后,给出了DGTD处理散射问题时平面波入射加入的具体实现方法。最后,给出了金属球、介质球和金属弹头宽带散射的算例,算例结果的比较表明了该文算法的正确性和有效性。该文的研究,为复杂目标雷达散射截面RCS的准确预估打下了坚实的基础。展开更多
Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these m...Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these methods have the flexibilitywhich is not shared by typical finite element methods, such as the allowance of ar-bitrary triangulation with hanging nodes, less restriction in changing the polynomialdegrees in each element independent of that in the neighbors (p adaptivity), and localdata structure and the resulting high parallel efficiency. In this paper, we give a generalreview of the local DG (LDG) methods for solving high-order time-dependent partialdifferential equations (PDEs). The important ingredient of the design of LDG schemes,namely the adequate choice of numerical fluxes, is highlighted. Some of the applica-tions of the LDG methods for high-order time-dependent PDEs are also be discussed.展开更多
One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for differ...One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth (multistep), and second-order total variation diminishing (TVD) Runge-Kutta (multistage) time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic (rectangular, trapezoidal, triangular, and parabolic cross-sections) and non-prismatic (natural river) channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.展开更多
基金supported by the State Key Program of National Natural Science Foundation of China(11931003)the National Natural Science Foundation of China(41974133)。
文摘In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.
基金Supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘A novel high-order three-dimensional(3-D)discontinuous Galerkin time domain(DGTD)method based on a normalized formulation of Maxwell′s equations is developed for modeling and simulating silicon-on-insulator(SOI)thin-ridge waveguide.The DGTD method employs unstructured meshes and piecewise high-order polynomials for spatial discretization,and Runge-Kutta methods for time integration.It is found that the numerical results of the leakage loss of SOI thin-ridge waveguide agree well with those of analytical solutions,which proves that the proposed method is an ideal tool for the quantitative analysis for SOI thin-ridge waveguide.
基金supported by the National Natural Science Foundation of China(Grant Nos.61301056 and 11176007)the Sichuan Provincial Science and Technology Support Program,China(Grant No.2013HH0047)+1 种基金the Fok Ying Tung Education Foundation,China(Grant No.141062)the"111"Project,China(Grant No.B07046)
文摘Several major challenges need to be faced for efficient transient multiscale electromagnetic simulations, such as flex- ible and robust geometric modeling schemes, efficient and stable time-stepping algorithms, etc. Fortunately, because of the versatile choices of spatial discretization and temporal integration, a discontinuous Galerkin time-domain (DGTD) method can be a very promising method of solving transient multiscale electromagnetic problems. In this paper, we present the application of a leap-frog DGTD method to the analyzing of the multiscale electromagnetic scattering problems. The uniaxial perfect matching layer (UPML) truncation of the computational domain is discussed and formulated in the leap-frog DGTD context. Numerical validations are performed in the challenging test cases demonstrating the accuracy and effectiveness of the method in solving transient multiscale electromagnetic problems compared with those of other numerical methods.
基金supported by a grant from the French National Ministry of Education and Research(MENSR,19755-2005)
文摘A high-order leap-frog based non-dissipative discontinuous Galerkin time-domain method for solving Maxwell's equations is introduced and analyzed.The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements,with a Nth-order leap-frog time scheme.Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes.The method is proved to be stable under some CFL-like condition on the time step.The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided.Numerical experiments with high-order elements show the potential of the method.
基金funded by the research project STiMulUs,ERC Grant agreement no.278267Financial support has also been provided by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2022 attributed to DICAM of the University of Trento(Grant L.232/2016)the PRIN2017 project.The authors have also received funding from the University of Trento via the Strategic Initiative Modeling and Simulation.
文摘A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.The staggered DG scheme defines the discrete pressure on the primal triangular mesh,while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh.In this paper,a new pair of equal-order-interpolation velocity-pressure finite elements is proposed.On the primary triangular mesh(the pressure elements),the basis functions are piecewise polynomials of degree N and are allowed to jump on the boundaries of each triangle.On the dual mesh instead(the velocity elements),the basis functions consist in the union of piecewise polynomials of degree N on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries,while they are continuous inside.In other words,the basis functions on the dual mesh arc built by continuous finite elements on the subtriangles.This choice allows the construction of an efficient,quadrature-free and memory saving algorithm.In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations,the arbitrary high order of accuracy in time is achieved through tire use of time-dependent test and basis functions,in combination with simple and efficient Picard iterations.Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes,not only from a computer memory point of view,but also concerning the computational time.
文摘For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.
文摘A Krylov space based time domain method for wave propagation problems is introduced. The proposed method uses the Arnoldi algorithm to obtain broad-band frequency domain solutions. This method is especially advantageous in cases where slow convergence is observed when using traditional time domain methods. The efficiency of the method is examined in several test cases to show its fast convergence in such problems.
基金Research sponsored by NSFC grants 11871428 and 12071214Nature Science Research Program for Colleges and Universities of Jiangsu Province grant 20KJB110011+1 种基金Research is supported in part by NSFC grants U1930402the fellowship of China Postdoctoral Science Foundation(No.2020TQ0030).
文摘In this paper,we first present the optimal error estimates of the semi-discrete ultra-weak discontinuous Galerkin method for solving one-dimensional linear convection-diffusion equations.Then,coupling with a kind of Runge-Kutta type implicit-explicit time discretization which treats the convection term explicitly and the diffusion term implicitly,we analyze the stability and error estimates of the corresponding fully discrete schemes.The fully discrete schemes are proved to be stable if the time-stepτ≤τ0,whereτ0 is a constant independent of the mesh-size h.Furthermore,by the aid of a special projection and a careful estimate for the convection term,the optimal error estimate is also obtained for the third order fully discrete scheme.Numerical experiments are displayed to verify the theoretical results.
基金This work is supported by NSFC grants No.11601490.
文摘In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restriction of explicit time marching methods,we employ a class of exponential time differencing(ETD)schemes for time integration,which is based on a linear convex splitting principle.Numerical experiments of the accuracy and long time simulations are given to show the efficiency and capability of the proposed numerical schemes.
文摘时域离散伽辽金法(Discontinuous Galerkin Time Domain,DGTD)同时具有时域有限元算法(Finite Element Time Domain,FETD)非结构网格剖分和时域有限差分算法(Finite Difference Time Domain,FDTD)显式迭代的优点,是一种非常有前途的电磁计算方法,该文首先描述了基于矢量基函数的时域离散伽辽金法的基本原理。然后,给出了DGTD处理散射问题时平面波入射加入的具体实现方法。最后,给出了金属球、介质球和金属弹头宽带散射的算例,算例结果的比较表明了该文算法的正确性和有效性。该文的研究,为复杂目标雷达散射截面RCS的准确预估打下了坚实的基础。
基金The research of the first author is support by NSFC grant 10601055,FANEDD of CAS and SRF for ROCS SEMThe research of the second author is supported by NSF grant DMS-0809086 and DOE grant DE-FG02-08ER25863.
文摘Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these methods have the flexibilitywhich is not shared by typical finite element methods, such as the allowance of ar-bitrary triangulation with hanging nodes, less restriction in changing the polynomialdegrees in each element independent of that in the neighbors (p adaptivity), and localdata structure and the resulting high parallel efficiency. In this paper, we give a generalreview of the local DG (LDG) methods for solving high-order time-dependent partialdifferential equations (PDEs). The important ingredient of the design of LDG schemes,namely the adequate choice of numerical fluxes, is highlighted. Some of the applica-tions of the LDG methods for high-order time-dependent PDEs are also be discussed.
文摘One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth (multistep), and second-order total variation diminishing (TVD) Runge-Kutta (multistage) time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic (rectangular, trapezoidal, triangular, and parabolic cross-sections) and non-prismatic (natural river) channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.