scholarly journals An adaptive finite‐difference method for traveltimes and amplitudes

Geophysics ◽  
2002 ◽  
Vol 67 (1) ◽  
pp. 167-176 ◽  
Author(s):  
Jianliang Qian ◽  
William W. Symes

The point‐source traveltime field has an upwind singularity at the source point. Consequently, all formally high‐order, finite‐difference eikonal solvers exhibit first‐order convergence and relatively large errors. Adaptive upwind finite‐difference methods based on high‐order Weighted Essentially NonOscillatory (WENO) Runge‐Kutta difference schemes for the paraxial eikonal equation overcome this difficulty. The method controls error by automatic grid refinement and coarsening based on a posteriori error estimation. It achieves prescribed accuracy at a far lower cost than does the fixed‐grid method. Moreover, the achieved high accuracy of traveltimes yields reliable estimates of auxiliary quantities such as take‐off angles and geometric spreading factors.

2015 ◽  
Vol 17 (2) ◽  
pp. 317-336 ◽  
Author(s):  
Wei Wang ◽  
Chi-Wang Shu ◽  
H.C. Yee ◽  
Dmitry V. Kotov ◽  
Björn Sjögreen

AbstractIn this paper, we extend the high order finite-difference method with subcell resolution (SR) in [34] for two-species stiff one-reaction models to multispecies and multireaction inviscid chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.


Sign in / Sign up

Export Citation Format

Share Document