[Liz1219]

Hi all! The p_rgh appeared "nan" when my case is solved with chtMultiRegionFoam for fluid region in first iteration.I tried many ways, but still couldn't solve it.Now, I am very upset then wishing lovely people can give me some pertinent advice.Vielen Dank!
Well ,these my codes of fluid.Have a look ,please.

mesh


field p

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location 0/domain;
object p;
}
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 101325;
boundaryField
{
domain_INLET
{
type calculated;
value $internalField;
}
domain_WALL
{
type calculated;
value $internalField;
}
domain_OUTLET
{
type calculated;
value $internalField;
}
domain_TOP
{
type calculated;
value $internalField;
}
domain_to_left
{
type cyclicAMI;
value $internalField;
}
domain_to_front
{
type cyclicAMI;
value $internalField;
}
domain_to_right
{
type cyclicAMI;
value $internalField;
}
domain_to_feet
{
type calculated;
value $internalField;
}
domain_to_shoe
{
type calculated;
value $internalField;
}
left_to_domain
{
type cyclicAMI;
value $internalField;
}
left_to_shoe
{
type calculated;
value $internalField;
}
front_to_domain
{
type cyclicAMI;
value $internalField;
}
front_to_shoe
{
type calculated;
value $internalField;
}
right_to_domain
{
type cyclicAMI;
value $internalField;
}
right_to_shoe
{
type calculated;
value $internalField;
}
}

field p_rgh

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location 0/domain;
object p_rgh;
}
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 101325;
boundaryField
{
domain_INLET
{
type fixedFluxPressure;
value $internalField;
}
domain_WALL
{
type fixedFluxPressure;
value $internalField;
}
domain_OUTLET
{
type uniformTotalPressure;
p0 $internalField;
}
domain_TOP
{
type fixedFluxPressure;
value $internalField;
}
domain_to_left
{
type cyclicAMI;
value $internalField;
}
domain_to_front
{
type cyclicAMI;
value $internalField;
}
domain_to_right
{
type cyclicAMI;
value $internalField;
}
domain_to_feet
{
type fixedFluxPressure;
value $internalField;
}
domain_to_shoe
{
type fixedFluxPressure;
value $internalField;
}
left_to_domain
{
type cyclicAMI;
value $internalField;
}
left_to_shoe
{
type fixedFluxPressure;
value $internalField;
}
front_to_domain
{
type cyclicAMI;
value $internalField;
}
front_to_shoe
{
type fixedFluxPressure;
value $internalField;
}
right_to_domain
{
type cyclicAMI;
value $internalField;
}
right_to_shoe
{
type fixedFluxPressure;
value $internalField;
}
}

field T

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location 0/domain;
object T;
}
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
domain_INLET
{
type fixedValue;
value $internalField;
}
domain_WALL
{
type zeroGradient;
}
domain_OUTLET
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}
domain_TOP
{
type zeroGradient;
}
domain_to_left
{
type cyclicAMI;
value $internalField;
}
domain_to_front
{
type cyclicAMI;
value $internalField;
}
domain_to_right
{
type cyclicAMI;
value $internalField;
}
domain_to_feet
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
}
domain_to_shoe
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
}
left_to_domain
{
type cyclicAMI;
value $internalField;
}
left_to_shoe
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
}
front_to_domain
{
type cyclicAMI;
value $internalField;
}
front_to_shoe
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
}
right_to_domain
{
type cyclicAMI;
value $internalField;
}
right_to_shoe
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod fluidThermo;
}
}

field U

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location 0/domain;
object U;
}
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
domain_INLET
{
type surfaceNormalFixedValue;
refValue uniform -3;
}
domain_WALL
{
type noSlip;
}
domain_OUTLET
{
type pressureInletOutletVelocity;
value $internalField;
}
domain_TOP
{
type noSlip;
}
domain_to_left
{
type cyclicAMI;
value $internalField;
}
domain_to_front
{
type cyclicAMI;
value $internalField;
}
domain_to_right
{
type cyclicAMI;
value $internalField;
}
domain_to_feet
{
type noSlip;
}
domain_to_shoe
{
type noSlip;
}
left_to_domain
{
type cyclicAMI;
value $internalField;
}
left_to_shoe
{
type noSlip;
}
front_to_domain
{
type cyclicAMI;
value $internalField;
}
front_to_shoe
{
type noSlip;
}
right_to_domain
{
type cyclicAMI;
value $internalField;
}
right_to_shoe
{
type noSlip;
}
}

field fvSchems

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system/domain";
object fvSchemes;
}
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(div(phi,U)) Gauss linear;
div(phi,U) Gauss linearUpwindV grad(U);
div(phi,h) bounded Gauss upwind;
div(phid,p) bounded Gauss upwind;
div(phi,K) Gauss linear;
div(phi,k) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,R) bounded Gauss upwind;
div((rho*R)) Gauss linear;
div(R) Gauss linear;
div(U) Gauss linear;
div(phi,muTilda) Gauss linear;
div(phi,Ekp) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear uncorrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default uncorrected;
}
fluxRequired
{
default no;
p;
}
wallDist
{
method meshWave;
}

field fvSolution

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system/domain";
object fvSolution;
}
solvers
{
U
{
solver PBiCGStab;
smoother DILU;
preconditioner DILU;
tolerance 1e-6;
relTol 0.01;
nSweeps 2;
}
epsilon
{
solver PBiCGStab;
smoother DILU;
preconditioner DILU;
tolerance 1e-6;
relTol 0.01;
nSweeps 2;
}
h
{
solver PBiCGStab;
smoother DILU;
preconditioner DILU;
tolerance 1e-6;
relTol 0.01;
nSweeps 2;
}
k
{
solver PBiCGStab;
smoother DILU;
preconditioner DILU;
tolerance 1e-6;
relTol 0.01;
nSweeps 2;
}
p_rgh
{
solver GAMG;
smoother GaussSeidel;
preconditioner GaussSeidel;
tolerance 1e-6;
relTol 0.01;
nSweeps 2;
}
rho
{
solver smoothSolver;
smoother GaussSeidel;
preconditioner GaussSeidel;
tolerance 1e-5;
relTol 0.01;
nSweeps 2;
}
UFinal
{
$U;
relTol 0;
}
epsilonFinal
{
$epsilon;
relTol 0;
}
hFinal
{
$h;
relTol 0;
}
kFinal
{
$k;
relTol 0;
}
p_rghFinal
{
$p_rgh;
relTol 0;
}
rhoFinal
{
$rho;
relTol 0;
}
}
PIMPLE
{
momentumPredictor yes;
transonic no;
consistent no;
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
residualControl
{
e 1e-3;
h 1e-3;
p 1e-2;
p_rgh 1e-3;
rho 1e-3;
t 1e-3;
turbulence 1e-3;
u 1e-4;
}
}
relaxationFactors
{
"p.*" 0.3;
p_rgh 0.3;
"rho.*" 0.1;
equations
{
"U.*" 0.7;
"h.*" 0.5;
"e.*" 0.5;
"(k|epsilon|omega|R|nuTilda|pollutant|epsilon) " 0.7;
T 0.5;
}
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
PhiRefCell 0;
PhiRefValue 0;
}
cache
{
grad(U);
}

Thanks all!