Always crash when solve a C-D nozzle flow field using rhoCentralFoam
 

dear all:
I have used the rhocentralFoam in OpenFoam to solve an CD nozzle with subsonic inlet and supersonic outlet.
But there is always an error after serval iteration steps even if I have tried varied boundary conditions for the case (you can find them in the notes of the P、U、T files in the 0 folder). Can anyone give me some advice?
thanks!

source code"0"
 

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1013250;

boundaryField
{
symmetry.5
{
type symmetryPlane;

}
wall.2
{
type zeroGradient;
}
/* pressure_inlet.3
{
type totalPressure;
p0 uniform 1.01325e+06;//:p_in
U U;
phi phi;
rho none;
psi none;
gamma 1.4;
value uniform 1.01325e+06;
}
*/
/*
pressure_outlet.4
{
type waveTransmissive;
value uniform 100; //important for correct I/O
field p; //the name of the field that we are working on
gamma 1.4; //the ratio of specific heats
phi phi; //the name of the volumetric flux field (or if you use the mass flux phi, it will be divided by rho)
rho rho; //the name of the density field
psi psi; //the name of the field that is the deriv. of density with respect to pressure
lInf 1; //a measure of how far away the far-field condition should be
fieldInf 100; //the far-field value to be applied to p
}*/
pressure_inlet.3
{
type fixedValue;
value uniform 1.01325e+06;
}
pressure_outlet.4
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 2000;
boundaryField
{
symmetry.5
{
type symmetryPlane;

}
wall.2
{
type zeroGradient;
}
pressure_outlet.4
{
type zeroGradient;
}
pressure_inlet.3
{
type fixedValue;
value uniform 2000;
}

frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (30 0 0);
boundaryField
{

symmetry.5
{
type symmetryPlane;

}
wall.2
{
type fixedValue;
value uniform (0 0 0);
}
pressure_outlet.4
{
type zeroGradient;
}
pressure_inlet.3
{
type pressureInletVelocity;
// type pressureDirectedInletOutletVelocity;
value uniform (30 0 0);
}
frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //

source code"constant"
 

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType ePsiThermo<pureMixture<constTransport<specieThermo <eConstThermo<perfectGas>>>>>;
mixture air 1 28.9 717.5 0 1.82e-05 0.73;
//thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;
//mixture air 1 28.9 1007 0 1.84e-05 0.7;
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType laminar;

// ************************************************** *********************** //

source code"system"
 

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application rhoCentralFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 100;
deltaT 1e-4;
writeControl adjustableRunTime;
writeInterval 1e-1;
cycleWrite 0;
writeFormat ascii;
writePrecision 15;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
adjustTimeStep yes;
maxCo 0.2;
maxDeltaT 1;

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fluxScheme Kurganov;
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(tauMC) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
reconstruct(rho) MUSCL;
reconstruct(U) MUSCLV;
reconstruct(T) MUSCL;
}
snGradSchemes
{
default corrected;
}

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
rho
{
solver diagonal;
}
rhoU
{
solver diagonal;
}
rhoE
{
solver diagonal;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-09;
relTol 0.01;
}
e
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-10 relTol 0;
}

}

// ************************************************** *********************** //