Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA high order finite difference method for the elastic wave equations in bounded domains with nonconforming interfaces
141
0
0.0
(
0
)
Added by
Lu Zhang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We develop a stable finite difference method for the elastic wave equations in bounded media, where the material properties can be discontinuous at curved interfaces. The governing equations are discretized in second order form by a fourth or sixth order accurate summation-by-parts operator. The mesh size is determined by the velocity structure of the material, resulting in nonconforming grid interfaces with hanging nodes. We use order-preserving interpolation and the ghost point technique to couple adjacent mesh blocks in an energy-conserving manner, which is supported by a fully discrete stability analysis. In numerical experiments, we demonstrate that the convergence rate is optimal, and is the same as when a globally uniform mesh is used in a single domain. In addition, with a predictor-corrector time integration method, we obtain time stepping stability with stepsize almost the same as given by the usual Courant Friedrichs Lewy condition.
rate research
Read More
We develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.
This article is an account of the NABUCO project achieved during the summer camp CEMRACS 2019 devoted to geophysical fluids and gravity flows. The goal is to construct finite difference approximations of the transport equation with nonzero incoming boundary data that achieve the best possible convergence rate in the maximum norm. We construct, implement and analyze the so-called inverse Lax-Wendroff procedure at the incoming boundary. Optimal convergence rates are obtained by combining sharp stability estimates for extrapolation boundary conditions with numerical boundary layer expansions. We illustrate the results with the Lax-Wendroff and O3 schemes.
In this paper, we combine the nonlinear HWENO reconstruction in cite{newhwenozq} and the fixed-point iteration with Gauss-Seidel fast sweeping strategy, to solve the static Hamilton-Jacobi equations in a novel HWENO framework recently developed in cite{mehweno1}. The proposed HWENO frameworks enjoys several advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function $phi$. They are now computed from the current value of $phi$ and the previous spatial derivatives of $phi$. This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO scheme. In addition, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed-point fast sweeping method is used to update the numerical approximation. It is an explicit method and does not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton-Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some known issues, including that the iterative number is very sensitive to the parameter $varepsilon$ used in the nonlinear weights, as observed in previous studies. Finally, in order to further reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two-dimensional problems, which demonstrate the good performance of the proposed fixed-point fast sweeping HWENO methods.
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 propose a novel Hermite weighted essentially non-oscillatory (HWENO) fast sweeping method to solve the static Hamilton-Jacobi equations efficiently. During the HWENO reconstruction procedure, the proposed method is built upon a new finite difference fifth order HWENO scheme involving one big stencil and two small stencils. However, one major novelty and difference from the traditional HWENO framework lies in the fact that, we do not need to introduce and solve any additional equations to update the derivatives of the unknown function $phi$. Instead, we use the current $phi$ and the old spatial derivative of $phi$ to update them. The traditional HWENO fast sweeping method is also introduced in this paper for comparison, where additional equations governing the spatial derivatives of $phi$ are introduced. The novel HWENO fast sweeping methods are shown to yield great savings in both computational time and storage, which improves the computational efficiency of the traditional HWENO scheme. In addition, a hybrid strategy is also introduced to further reduce computational costs. Extensive numerical experiments are provided to validate the accuracy and efficiency of the proposed approaches.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
