In this work, we implement and compare two fourth-order compact finite difference (CFD) methods to model the propagation of acoustic waves on a 1-D domain. The first scheme is based on standard implicit CFD constructed on nodal grids and must solve a couple of tridiagonal linear systems at each time iteration. The formulation of the second method uses the novel staggered mimetic CFD that explicitly approximate derivatives at cell center points. Mimetic FD operators inherit this name from the fact that some conservation properties fulfilled by the continuous divergence, gradient, and curl operators, are also satisfied by the discrete mimetic counterparts. Both CFD methods are combined to high-order Runge-Kutta schemes for time integration. In our current preliminary experiments, results show that the mimetic CFD scheme is slightly more accurate and shows similar fourth-order convergence than the nodal scheme. In this paper, we also compare the simulation times spent by the compact schemes to an explicit Leap-frog staggered solver. This comparison reflects the high computational efficiency of the mimetic CFD scheme, and along with its accurate results, suggests the potential of this method for multidimensional wave propagation problems.


Article metrics loading...

Loading full text...

Full text loading...


  1. Abouali, M. and Castillo, J. E.
    [2012] High-order compact Castillo-Grone’s Mimetic Operators. Report of Computational Science Research Center at San Diego State University.CSRCR021-13.
    [Google Scholar]
  2. Carpenter, M., Gottlieb, D., Abarbanel, S.
    [1993] Stable and accurate boundary treatments for compact, high-order finite-difference scheme. Applied Numerical Mathematics, 12, 55–87.
    [Google Scholar]
  3. Castillo, J. E. and Grone, R. D.
    [2003] A matrix analysis approach to higher order approximations for divergence and gradient satisfying a global conservation law. SIAM Journal of Matrix Analysis and Applications, 25, 128–142.
    [Google Scholar]
  4. Lele, S.
    [1992] Compact finite difference schemes with spectral-like resolution. Journal of Computational Physics, 103, 16–42.
    [Google Scholar]
  5. Rojas, O., Day, S., Castillo, J. and Dalguer, L. A.
    [2008] Modelling of rupture propagation using high-order mimetic finite differences. Geophysical. Journal International, 172, 631–650.
    [Google Scholar]
  6. Rojas, O., Otero, B., Castillo, J. E. and Day, S. M.
    [2014] Low dispersive modeling of Rayleigh waves on partly staggered grids. Computational Geosciences, 18(1), 29–43.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error