We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the...We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the resulting HHO-Newmark scheme is energy-conservative,and this scheme is also amenable to static condensation at each time step.We then consider the formulation of the acoustic wave equation as a frst-order system together with singly-diagonally implicit and explicit Runge-Kutta(SDIRK and ERK)schemes.HHO-SDIRK schemes are amenable to static condensation at each time step.For HHO-ERK schemes,the use of the mixed-order formulation,where the polynomial degree of the cell unknowns is one order higher than that of the face unknowns,is key to beneft from the explicit structure of the scheme.Numerical results on test cases with analytical solutions show that the methods can deliver optimal convergence rates for smooth solutions of order O(hk+1)in the H1-norm and of order O(h^(k+2))in the L^(2)-norm.Moreover,test cases on wave propagation in heterogeneous media indicate the benefts of using high-order methods.展开更多
We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusio...We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms： a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.展开更多
基金The authors would like to thank L.Guillot(CEA/DAM)for insightful discussions and CEA/DAM for partial fnancial support.EB was partially supported by the EPSRC grants EP/P01576X/1 and EP/P012434/1.
文摘We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the resulting HHO-Newmark scheme is energy-conservative,and this scheme is also amenable to static condensation at each time step.We then consider the formulation of the acoustic wave equation as a frst-order system together with singly-diagonally implicit and explicit Runge-Kutta(SDIRK and ERK)schemes.HHO-SDIRK schemes are amenable to static condensation at each time step.For HHO-ERK schemes,the use of the mixed-order formulation,where the polynomial degree of the cell unknowns is one order higher than that of the face unknowns,is key to beneft from the explicit structure of the scheme.Numerical results on test cases with analytical solutions show that the methods can deliver optimal convergence rates for smooth solutions of order O(hk+1)in the H1-norm and of order O(h^(k+2))in the L^(2)-norm.Moreover,test cases on wave propagation in heterogeneous media indicate the benefts of using high-order methods.
文摘We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms： a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.
