[ali.m.1]

Hello

I'm running a submerged jet impingement test using DPMFoam (3.0.x), and although in the log file it says that the pressure is solving (GAMG: Solving for p, Initial residual = 0.006891, Final residual = 5.435e-05, No Iterations 3), and that it converges, when I open a time directory, the whole field is a uniform 1e5 which is my outlet BC.

Here is my p file in my 0 directory:

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

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

internalField   uniform 1e5;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    walls
    {
        type            zeroGradient;
    }
    pOutTop
    {
        type            totalPressure;
        p0              1e+05;
        U               U.water;
        phi             phi.water;
        rho             none;
        psi             none;
        gamma           1;
        value           uniform 1e+05; 
       /* type            fixedValue;
        value           100000;*/
    }
    bottomOutlet
    {
        type            totalPressure;
        p0              1e+05;
        U               U.water;
        phi             phi.water;
        rho             none;
        psi             none;
        gamma           1;
        value           uniform 1e+05; 
      /*  type            fixedValue;
        value           100000;*/
    }
    impingement
    {
        type            zeroGradient;
    }

    symmetry
    {
        type            symmetry;
    }
}

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

Does anyone have any ideas? I'm running a similar case on reactingTwoPhaseEulerFoam with the same BCs and this gives me expected pressure plots.

Ali

[ali.m.1]

Hi All

I realised that if I change the initial p file then the solver runs fine. The problem was that I was confusing gauge pressure and absolute.

Here is the p file that works, in case it will help anyone else in the future:

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

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

internalField   uniform 0;//1e5;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    walls
    {
        type            zeroGradient;
    }
    pOutTop
    {
        type            totalPressure;
        p0              0;//1e+05;
        U               U.water;
        phi             phi.water;
        rho             none;
        psi             none;
        gamma           1;
        value           uniform 0;//1e+05; 
       /* type            fixedValue;
        value           100000;*/
    }
    bottomOutlet
    {
        type            totalPressure;
        p0              0;//1e+05;
        U               U.water;
        phi             phi.water;
        rho             none;
        psi             none;
        gamma           1;
        value           uniform 0;//1e+05; 
      /*  type            fixedValue;
        value           100000;*/
    }
    impingement
    {
        type            zeroGradient;
    }

    symmetry
    {
        type            symmetry;
    }
}

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