[anydimenx]

Hi All,
In OpenFOAM version 4.1, in the directory:
tutorials/multiphase/compressibleInterDyMFoam/ras/sloshingTank2D/0

there are both p and p_rgh file
in file p:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / 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 1e6;

......

I don't understand why the object is p_rgh in file p, I think the object name should be consistent with the file name.

in file p_rgh:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / 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 1e6;

......

both files set the initial field 1e6 uniformly, if I run compressibleInterDyMFoam, it will succeed. However, if I change the initial field to 0 uniformly, which I think the dynamic pressure should be 0 initially, then the solver compressibleInterDyMFoam crash, it says the floating point exception, does it mean something is divided by the pressure?

So why there are two p_rgh objects. In the cases of other solvers, there is only p_rgh file, because the static pressure can be calculated out, we only need to specify the dynamic pressure p_rgh.

And in another directory:
tutorials/multiphase/compressibleInterFoam/laminar/depthCharge2D/0.orig:

there are also p and p_rgh files,

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

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

internalField uniform 1e5;

......

Now the object name is consistent with the file name.

in file p_rgh:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / 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;

......

So why should we specify both p and p_rgh? Why p_rgh is not 0 initially? Could you please explain them? I am confusing now.