[Manta]

Hey everyone,

I've got a problem with negative temeratures using the rhoCentralFoam solver. In my 2d case a moving piston is placed in a cylinder. While the piston is moving, the fluid shoud flow along that piston.
I used the movingConeTopoFvMesh file together with layerAdditionRemoval. After a few iterations, a new layer of cells is added. One iteration later you can already see the temperature diverging and another iteration later I get the error "negative initial temperature".

My boundary conditions are:
-p has always zerogradient to all walls and my symmetryaxis has symmetric
-T exactly like p
-U walls noSlip and the moving piston got movingWallVelocity.
-My front and back bc's are always empty.

My fvSolution looks like this:

solvers
{
    p
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-06;
        relTol          0.01;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "pcorr.*"
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-2;
        relTol          0;
    }

    "rho.*"
    {
        solver          diagonal;
        tolerance       1e-05;
        relTol          0;
    }

    "(U|h|e|R|k|epsilon|omega)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0.1;
    }

    "(U|h|e|R|k|epsilon|omega)Final"
    {
        $U;
        relTol          0;
    }

    cellMotionUx
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-08;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor   yes;
    correctPhi          yes;
    nOuterCorrectors    2;
    nCorrectors         2;
    transonic           yes;
    nNonOrthogonalCorrectors 0;

    rhoMin              0.1;
    rhoMax              100.0;
    temperatureMin	 250;
}

relaxationFactors
{
    equations
    {	
        ".*"        1;
    }
}

My fvSchemes looks like this:

fluxScheme      Kurganov;

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    
}

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;

    laplacian(diffusivity,cellMotionU) Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
    reconstruct(rho) vanLeer;
    reconstruct(U)  vanLeerV;
    reconstruct(T)  vanLeer;
}

snGradSchemes
{
    default         corrected;
}

Does anyone see my mistake?