[y_jiang]

Hi Sufjanst,

Thank you for the suggestions.

I have made it work by change the BCs to the following combo
p:
inlet: totalPressure
outlet: fixedValue
U:
inlet: pressureInletVelocity
outlet: zeroGradient
Then, I gradually increase the inlet pressure to reach the desired pressure difference. It finally gave me reasonable result.
I have not really tried to use potentialFoam to initialize the field. But your suggestion of replacing sensibleInternalEnergy by sensibleEnthalpy do help me. It makes the model converge faster and more stable. Could you explain what is the reason?

Thank you very much!

[Sakun]

Originally Posted by y_jiang

Hello Antoine,


I got the same issue as you had. The calculation stopped at 9th iteration and came up with the error when solving energy equation.

My goal is to solve a pressure driven problem for a tank with 1 inlet and 2 outlet. As you guys discussed here, I knew that a given value of outlet pressure should work fine with a flow rate inlet. So I set my boundary conditions like below. But I still got the error.

Do you have any idea on it?

Any of your input would be appreciated. Thank you!

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    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 600000;

boundaryField
{

    geofull_inlet1 
    {
        type            zeroGradient;
        refValue        uniform 600000;
        refGradient     uniform 0;
    }

    geofull_outlet2 
    {
        type            fixedValue;
        value           uniform 600000;
    }

    geofull_outlet3 
    {
        type            fixedValue;
        value           uniform 600000;
    }

    geofull_needle1
    {
        type            zeroGradient;
    }

    geofull_needle2
    {
        type            zeroGradient;
    }

    geofull_needle3
    {
        type            zeroGradient;
    }

    geofull_tank 
    {
        type            zeroGradient;

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

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

internalField   uniform (0 0 0);

boundaryField
{

    geofull_inlet1 
    {
        type            flowRateInletVelocity;
        volumetricFlowRate    constant 0.5;
        rhoInlet        1;
//        massFlowRate    0.5;
//        rhoInlet        0.5;
    }

    geofull_outlet2 
    {
        type            pressureInletOutletVelocity; 
        value           uniform (0 0 0);
    }

    geofull_outlet3 
    {
        type            pressureInletOutletVelocity; 
        value           uniform (0 0 0);
    }

    geofull_needle1 
    {
        type           noSlip;
    }

    geofull_needle2 
    {
        type           noSlip;
    }

    geofull_needle3 
    {
        type           noSlip;
    }

    geofull_tank 
    {
        type           noSlip;
    }

}

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

sorry for commenting on an old thread,



did you get the solution for your problem ? cause i have use 2 pressure values as well.