[gtpunnamarajuaju@gmail.co]

Hi,
I'm trying to run a case with a droplet of water falling in a liquid film of particular thickness with the energy equation. So I've used the compressibleInterFoam solver for this. While executing the code, I got an error, the screenshot of which, is attached.

Attachment 50911

I assume that there is an error in my 'p_rgh' file or my 'fvsolution' file, but I'm unable to find what my error is.
I've attached even my 'p_rgh' and 'fv_solution' files. Could someone help me find out what my error is and help me rectify it.




Thank you so much!


fvsolution file:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
alpha.water
{
nAlphaCorr 1;
nAlphaSubCycles 1;
cAlpha 1;
}

pcorr
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-10;
relTol 0;
smoother DICGaussSeidel;
}
tolerance 1e-05;
relTol 0;
maxIter 100;
}

".*(rho|rhoFinal)"
{
solver diagonal;
}

p_rgh
{
solver GAMG;
tolerance 1e-07;
relTol 0.01;
smoother DIC;
}
/* p_rgh
{
solver GAMG;
preconditioner
{
preconditioner DIC;
tolerance 1e-07;
relTol 0;
nVcycles 2;
smoother DICGaussSeidel;
nPreSweeps 2;
}
tolerance 1e-07;
relTol 0;
maxIter 20;
}*/

p_rghFinal
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-07;
relTol 0;
nVcycles 2;
smoother DICGaussSeidel;
nPreSweeps 2;
}
tolerance 1e-07;
relTol 0;
maxIter 20;
}

U
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0;
nSweeps 1;
}

"(T|k|B|nuTilda).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-08;
relTol 0;
}
}

PIMPLE
{
momentumPredictor no;
transonic no;
nOuterCorrectors 1;
nCorrectors 4;
nNonOrthogonalCorrectors 0;
}


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






p-rgh file:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 1e5;

boundaryField
{
walls
{
type fixedValue;
value uniform 0;
}
bottom
{
type fixedValue;
value uniform 0;
}
atmosphere
{
type totalPressure;
p0 uniform 0;
U U;
T T;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;
}
frontAndBack
{
type fixedValue;
value uniform 0;

}
}


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


If the screen shot isn't clear, the following is the error:

PIMPLE: iteration 1
MULES: Solving for alpha.water
Liquid phase volume fraction = 0.1 Min(alpha.water) = 0 Min(alpha.air) = 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for T, Initial residual = 0.603844, Final residual = 0.365488, No Iterations 1000
min(T) 300
GAMG: Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000
max(U) nan
min(p_rgh) 0
GAMG: Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000


--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /auto/home/vgpunnamaraju/depthCharge3D/system/data.solverPerformance.p_rgh at line 0.

From function Foam::Istream& Foam:perator>>(Foam::Istream&, Foam::doubleScalar&)
in file lnInclude/Scalar.C at line 93.

FOAM exiting