[Josefina]

hello everybody!

I am running an interFoam case and I realize that in the log file I don't have the values of U.

The example of a time step is given below:

Courant Number mean: 0.0308112 max: 0.577916
Interface Courant Number mean: 0.00290818 max: 0.455244
deltaT = 0.00239286
Time = 0.0162499

MULES: Solving for alpha1
Phase-1 volume fraction = 0.506127 Min(alpha1) = -3.7238e-39 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.506358 Min(alpha1) = -1.08909e-31 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.00222856, Final residual = 0.000102936, No Iterations 3
time step continuity errors : sum local = 0.000328819, global = -3.13238e-05, cumulative = -0.000744183
DICPCG: Solving for p_rgh, Initial residual = 0.000265302, Final residual = 9.37021e-06, No Iterations 9
time step continuity errors : sum local = 3.01083e-05, global = -7.00703e-06, cumulative = -0.00075119
DICPCG: Solving for p_rgh, Initial residual = 3.68985e-05, Final residual = 8.97875e-08, No Iterations 28
time step continuity errors : sum local = 2.86783e-07, global = -1.28258e-07, cumulative = -0.000751318
ExecutionTime = 0.74 s ClockTime = 0 s

Does anyone know why it does not appear?

in every timestep file, I have the U file so OpenFoam calculates it but I would be able to see and plot residuals also.

I thank you very much for your answer

[alexeym]

Hi,

Do you have

PIMPLE
{
...
    momentumPredictor   yes;
...
}

in your fvSolution file?

If you take a look at UEqn.H:

...
    if (pimple.momentumPredictor())
    {
        solve
        (
            UEqn
         ==
            fvc::reconstruct
            (
                (
                    fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
                  - ghf*fvc::snGrad(rho)
                  - fvc::snGrad(p_rgh)
                ) * mesh.magSf()
            )
        );
    }

if you switch off momentum predictor step there will be no solver output about U field.

[Josefina]

Hi!

thank you very much!!
Indeed, the problem was I had put "no" in the "momentumPredictor". As I am quite new with OpenFoam I have yet the good reflexes. :$

So I am very grateful for your answer, now I can see both p and U