Hi,
I am trying to use Chebyshev pseudospectral method to simulate 3D wave propagation problems. The solver works but is very unstable. I applied mapping functions to make the boundary mesh grid larger but still instability is very apparent.
I read staggered grid tends is helpful in improving the stability of a finite difference solver. I am wondering if that can applied to pseudospectral method. How can it be used in Chebyshev method, in particular?
Any suggestions will be appreciated. Thanks.
Bo