[matthewshorter]

Hello,
I am working on a shock tube simulation as part of a research project. My current difficulty is getting a non oscillating value for p,T,U, and rho at the contact surface of the fluid domain. In the attached image you can see the small deviation in the plot. I have also attached the fvSchemes file. Currently, I am using gamma 1 interpolation but vanLeer and vanAlbada give similar results, with gauss linear and superbee causing the solver to crash. I have read the original paper on CentralFoam by Greenshields etal. and they briefly mentioned that the contact region had dubious accuracy. I am observing excellent results for shock velocity and strength, but have concerns regarding impact of the oscillations around contact surface affecting rest of simulation.

I am somewhat confused as a limiter is supposed to make sure that the fields are behaving monotonically throughout the domain. Is the numerical error/oscillation caused by the discontinuity inevitable or can I mitigate it by changing the solver settings? This is currently in OF 3.0.

fluxScheme Kurganov;

ddtSchemes
{
default Euler;
}

gradSchemes
{
// default Gauss;// linear;
default leastSquares 1;
}

divSchemes
{
default none;
div(tauMC) Gauss linear;
div(phi,k) Gauss linear;
div(phi,omega) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
reconstruct(rho) Gamma 1;
reconstruct(U) GammaV 1;
reconstruct(T) Gamma 1;
}

snGradSchemes
{
default corrected;
}

In system: maxCo = 0.1