I'm part of a team trying to generalize a 1D Advection-Diffussion-Reaction code we inherited by extending it to 2D by using dimensional splitting, i.e. solving advection and diffusion for x and y separately with 1D schemes.
Diffusion was (relatively) painless but there are issues with advection.
Originally the code had a flux limited, third order upwind-biased scheme for 1D advection which worked like a charm in the 1D case. However, when the same scheme is extended to 2D, crazy instabilities develop. What supposed to be a radial advection-diffusion profile wrinkles in the transverse direction and shoots a "jet-like" instability:
However when I use a first order upwind scheme I do get a nice, symmetric radial profile. The problem with the first order upwind scheme solution is the excessive false diffusion which the higher order scheme does suppress. But the higher order scheme is even more problematic because of the aforementioned instability.
I was wondering if higher order schemes in general are a bad idea in a multidimensional case? Are multidimensional problems more succeptible to "oscillations" caused by high order schemes? What about flux limiters, should 1D flux limiters prevent this oscillations from happening (in my situation, it seems the flux limiter for the third order case doesn't help)? Should I stick to lower order schemes even if they bring false diffusion? This third order scheme I was using is also flux limited, but the flux limit doesn't stop the instabilities...