Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA conservative and energy stable discontinuous spectral element method for the shifted wave equation in second order form
108
0
0.0
(
0
)
Added by
Kenneth Duru
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we develop a provably energy stable and conservative discontinuous spectral element method for the shifted wave equation in second order form. The proposed method combines the advantages and central ideas of very successful numerical techniques, the summation-by-parts finite difference method, the spectral method and the discontinuous Galerkin method. We prove energy-stability, discrete conservation principle, and derive error estimates in the energy norm for the (1+1)-dimensions shifted wave equation in second order form. The energy-stability results, discrete conservation principle, and the error estimates generalise to multiple dimensions using tensor products of quadrilateral and hexahedral elements. Numerical experiments, in (1+1)-dimensions and (2+1)-dimensions, verify the theoretical results and demonstrate optimal convergence of $L^2$ numerical errors at subsonic, sonic and supersonic regimes.
rate research
Read More
We develop an energy-based finite difference method for the wave equation in second order form. The spatial discretization satisfies a summation-by-parts (SBP) property. With boundary conditions and material interface conditions imposed weakly by the simultaneous-approximation-term (SAT) method, we derive energy estimates for the semi-discretization. In addition, error estimates are derived by the normal mode analysis. The energy-based discretization does not use any mesh-dependent parameter, even in the presence of Dirichlet boundary conditions and material interfaces. Furthermore, similar to upwind discontinuous Galerkin methods, numerical dissipation can be added to the discretization through the boundary conditions. We present numerical experiments that verify convergence and robustness of the proposed method.
In this paper we present a numerical discretization of the coupled elasto-acoustic wave propagation problem based on a Discontinuous Galerkin Spectral Element (DGSE) approach in a three-dimensional setting. The unknowns of the coupled problem are the displacement field and the velocity potential, in the elastic and the acoustic domains, respectively, thereby resulting in a symmetric formulation. After stating the main theoretical results, we assess the performance of the method by convergence tests carried out on both matching and non-matching grids, and we simulate realistic scenarios where elasto-acoustic coupling occurs. In particular, we consider the case of Scholte waves and the scattering of elastic waves by an underground acoustic cavity. Numerical simulations are carried out by means of the code SPEED, available at http://speed.mox.polimi.it.
The locally modified finite element method, which is introduced in [Frei, Richter: SINUM 52(2014), p. 2315-2334] is a simple fitted finite element method that is able to resolve weak discontinuities in interface problems. The method is based on a fixed structured coarse mesh, which is then refined into sub-elements to resolve an interior interface. In this work, we extend the locally modified finite element method to second order using an isoparametric approach in the interface elements. Thereby we need to take care that the resulting curved edges do not lead to degenerate sub-elements. We prove optimal a priori error estimates in the $L^2$-norm and in a modified energy norm, as well as a reduced convergence order of ${cal O}(h^{3/2})$ in the standard $H^1$-norm. Finally, we present numerical examples to substantiate the theoretical findings.
We build a multi-element variant of the smoothness increasing accuracy conserving (SIAC) shock capturing technique proposed for single element spectral methods by Wissink et al. (B.W. Wissink, G.B. Jacobs, J.K. Ryan, W.S. Don, and E.T.A. van der Weide. Shock regularization with smoothness-increasing accuracy-conserving Dirac-delta polynomial kernels. Journal of Scientific Computing, 77:579--596, 2018). In particular, the baseline scheme of our method is the nodal discontinuous Galerkin spectral element method (DGSEM) for approximating the solution of systems of conservation laws. It is well known that high-order methods generate spurious oscillations near discontinuities which can develop in the solution for nonlinear problems, even when the initial data is smooth. We propose a novel multi-element SIAC filtering technique applied to the DGSEM as a shock capturing method. We design the SIAC filtering such that the numerical scheme remains high-order accurate and that the shock capturing is applied adaptively throughout the domain. The shock capturing method is derived for general systems of conservation laws. We apply the novel SIAC filter to the two-dimensional Euler and ideal magnetohydrodynamics (MHD) equations to several standard test problems with a variety of boundary conditions.
We present unconditionally energy stable Runge-Kutta (RK) discontinuous Galerkin (DG) schemes for solving a class of fourth order gradient flows. Our algorithm is geared toward arbitrarily high order approximations in both space and time, while energy dissipation remains preserved without imposing any restriction on time steps and meshes. We achieve this in two steps. First, taking advantage of the penalty free DG method introduced by Liu and Yin [J Sci. Comput. 77:467--501, 2018] for spatial discretization, we reformulate an extended linearized ODE system by the energy quadratization (EQ) approach. Second, we apply an s-stage algebraically stable RK method for temporal discretization. The resulting fully discrete DG schemes are linear and unconditionally energy stable. In addition, we introduce a prediction-correction procedure to improve both the accuracy and stability of the scheme. We illustrate the effectiveness of the proposed schemes by numerical tests with benchmark problems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
