# Is there a simple way to avoid carbuncles for FD WENO methods?

I have implemented finite-difference WENO scheme for Euler equations (with some variants - WENO-JS, WENO-Z, WENO-M, different flux splitting). It works well, but have problem with so-called carbuncles (shock instability in the areas where its front is almost aligned with the grid). This problem seems to be well-studied for finite-volume methods (quick search shows that there are some fixes for Riemann solvers), but I can't find almost any information regarding carbuncles in finite-difference methods.

Are there any known fixes for finite-difference methods?

• How about H-correction? Some ideas: doi.org/10.1006/jcph.1998.6047 Nov 2 '17 at 18:21
• What flux splitting are you using ? Can you give more details in your question ? Nov 3 '17 at 3:49
• @Praveen It seems that correct flux splitting is very important. I was using 3 versions (global and local Lax-Friedrich's and Roe-LLF) and was constructing flux splitting after weno-reconstruction of U and F(U) in local characteristic field. Now I weno-reconstruct already splitted flux. There are no more large-scaled carbuncles, but still are some weak shock front instabilities, so the question still stands. Best results so far are obtained by proper Roe-LLF splitting with Roe averaging of Jacobians at $x_{j+1/2}$. Nov 5 '17 at 12:05
• Right, FD-WENO involves first doing a flux splitting and then applying weno on the split fluxes. Can you give me some reference for the Roe-LLF that you used ? Nov 5 '17 at 12:15
• @Praveen It's formula (2.6) from Liu, X. D., Osher, S., & Chan, T. (1994). Weighted essentially non-oscillatory schemes. Journal of computational physics, 115(1), 200-212. Here is the link: paper. In order to choose flux (f(a), f(b) or LLF-splitted), I check eigenvalues over the whole 6-point stencil. Curiously, WENO-Z is more prone to such instability type (probably, because of smaller dissipation near shocks) than WENO-JS. Nov 6 '17 at 12:43