scholarly journals Staggered-Grid Finite Difference Method with Variable-Order Accuracy for Porous Media

2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Jinghuai Gao ◽  
Yijie Zhang
Keyword(s):  
Porous Media ◽  
Finite Difference ◽  
Finite Difference Methods ◽  
Computation Time ◽  
High Order ◽  
Staggered Grid ◽  
Difference Operators ◽  
Variable Order ◽  
Low Velocity ◽  
Order Accuracy

The numerical modeling of wave field in porous media generally requires more computation time than that of acoustic or elastic media. Usually used finite difference methods adopt finite difference operators with fixed-order accuracy to calculate space derivatives for a heterogeneous medium. A finite difference scheme with variable-order accuracy for acoustic wave equation has been proposed to reduce the computation time. In this paper, we develop this scheme for wave equations in porous media based on dispersion relation with high-order staggered-grid finite difference (SFD) method. High-order finite difference operators are adopted for low-velocity regions, and low-order finite difference operators are adopted for high-velocity regions. Dispersion analysis and modeling results demonstrate that the proposed SFD method can decrease computational costs without reducing accuracy.

Temporal high-order staggered-grid finite-difference schemes for elastic wave propagation

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. T207-T224 ◽  
Author(s):  
Zhiming Ren ◽  
Zhen Chun Li
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Finite Difference ◽  
Elastic Wave ◽  
High Order ◽  
Elastic Wave Propagation ◽  
Staggered Grid ◽  
S Wave ◽  
Order Accuracy ◽  
Spatial Derivatives

The traditional high-order finite-difference (FD) methods approximate the spatial derivatives to arbitrary even-order accuracy, whereas the time discretization is still of second-order accuracy. Temporal high-order FD methods can improve the accuracy in time greatly. However, the present methods are designed mainly based on the acoustic wave equation instead of elastic approximation. We have developed two temporal high-order staggered-grid FD (SFD) schemes for modeling elastic wave propagation. A new stencil containing the points on the axis and a few off-axial points is introduced to approximate the spatial derivatives. We derive the dispersion relations of the elastic wave equation based on the new stencil, and we estimate FD coefficients by the Taylor series expansion (TE). The TE-based scheme can achieve ([Formula: see text])th-order spatial and ([Formula: see text])th-order temporal accuracy ([Formula: see text]). We further optimize the coefficients of FD operators using a combination of TE and least squares (LS). The FD coefficients at the off-axial and axial points are computed by TE and LS, respectively. To obtain accurate P-, S-, and converted waves, we extend the wavefield decomposition into the temporal high-order SFD schemes. In our modeling, P- and S-wave separation is implemented and P- and S-wavefields are propagated by P- and S-wave dispersion-relation-based FD operators, respectively. We compare our schemes with the conventional SFD method. Numerical examples demonstrate that our TE-based and TE + LS-based schemes have greater accuracy in time and better stability than the conventional method. Moreover, the TE + LS-based scheme is superior to the TE-based scheme in suppressing the spatial dispersion. Owing to the high accuracy in the time and space domains, our new SFD schemes allow for larger time steps and shorter operator lengths, which can improve the computational efficiency.

A High-Order Flux Reconstruction Method for 2-D Vorticity Transport

2021 ◽  
Author(s):  
Adrin Gharakhani
Keyword(s):  
Finite Difference ◽  
Poisson Equation ◽  
Finite Difference Methods ◽  
Neumann Boundary ◽  
Normal Derivative ◽  
High Order ◽  
Benchmark Problems ◽  
Flux Reconstruction ◽  
Reconstruction Method ◽  
Vorticity Transport

Abstract A compact high-order finite difference method on unstructured meshes is developed for discretization of the unsteady vorticity transport equations (VTE) for 2-D incompressible flow. The algorithm is based on the Flux Reconstruction Method of Huynh [1, 2], extended to evaluate a Poisson equation for the streamfunction to enforce the kinematic relationship between the velocity and vorticity fields while satisfying the continuity equation. Unlike other finite difference methods for the VTE, where the wall vorticity is approximated by finite differencing the second wall-normal derivative of the streamfunction, the new method applies a Neumann boundary condition for the diffusion of vorticity such that it cancels the slip velocity resulting from the solution of the Poisson equation for the streamfunction. This yields a wall vorticity with order of accuracy consistent with that of the overall solution. In this paper, the high-order VTE solver is formulated and results presented to demonstrate the accuracy and convergence rate of the Poisson solution, as well as the VTE solver using benchmark problems of 2-D flow in lid-driven cavity and backward facing step channel at various Reynolds numbers.

Optimized staggered-grid finite-difference operators using window functions

2018 ◽  
Vol 15 (2) ◽  
pp. 253-260 ◽  
Author(s):  
Ying-Jun Ren ◽  
Jian-Ping Huang ◽  
Peng Yong ◽  
Meng-Li Liu ◽  
Chao Cui ◽  
...  
Keyword(s):  
Finite Difference ◽  
Staggered Grid ◽  
Difference Operators ◽  
Window Functions
Sign in / Sign up

Export Citation Format

Share Document