CFD Online Discussion Forums - problem with InterFOAM interface courant number

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - problem with InterFOAM interface courant number (https://www.cfd-online.com/Forums/openfoam/165184-problem-interfoam-interface-courant-number.html)

problem with InterFOAM interface courant number

Dear Foamers,

I got a problem with interfoam running.

I am trying to simulate two phase flow with a submerged fixed body.
I use setField utility to set the two phase region and the water phase has a constant current speed 2m/s

But I dont know why the interface courant number is still 0 after several time step

Any suggestion about how to address the problem?

Regards,
Yage

Code:

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.727558664  Min(alpha1) = 0  Max(alpha1) = 1.00000004

GAMG:  Solving for p_rgh, Initial residual = 1.4702655e-07, Final residual = 2.15068659e-08, No Iterations 1

GAMG:  Solving for p_rgh, Initial residual = 2.19367882e-08, Final residual = 2.19367882e-08, No Iterations 0

time step continuity errors : sum local = 1.23804762e-11, global = 1.5117775e-13, cumulative = 1.51087321e-11

smoothSolver:  Solving for omega, Initial residual = 0.11466698, Final residual = 5.25033964e-09, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.469281825, Final residual = 2.45700553e-08, No Iterations 2

ExecutionTime = 1609.04 s  ClockTime = 1614 s



Courant Number mean: 0.0204447859 max: 0.499108417

Interface Courant Number mean: 0 max: 0

deltaT = 0.00011661004

Time = 0.00172344337

Looks like there is no interface. Can you post your setFieldsDict?

Shashank, there must be an interface because the phase fraction is 0.7 and min and max are 0 and 1.

Yage, did you look at the field at this point with Paraview? Does it make sense what you see?

Quote:

Originally Posted by smog (Post 580506)

Looks like there is no interface. Can you post your setFieldsDict?

Dear Shashank,

The setFieldsDict are as follow:

Code:

regions

(

    boxToCell

    {

        box     (-1 -1 -1.5) (1.55  1 0.02);

        fieldValues

        (

            volScalarFieldValue alpha.water 1

            volVectorFieldValue U (2 0 0)

        );

    }

Quote:

Originally Posted by akidess (Post 580559)

Shashank, there must be an interface because the phase fraction is 0.7 and min and max are 0 and 1.

Yage, did you look at the field at this point with Paraview? Does it make sense what you see?

Dear Anton,

Yes, I look at the domain with paraview and there is a freesurface.

Quote:

Originally Posted by akidess (Post 580559)

Shashank, there must be an interface because the phase fraction is 0.7 and min and max are 0 and 1.

Yage, did you look at the field at this point with Paraview? Does it make sense what you see?

I checked the free surface at 1.25s, There should be a wave form generated above the submerged body but there is not.

and the outcome after a longer running is still weird. the courant number is very small at the free surface.

Code:

PIMPLE: iteration 3

smoothSolver:  Solving for alpha.water, Initial residual = 1.19031101e-06, Final residual = 2.16770716e-13, No Iterations 2

Phase-1 volume fraction = 0.727573605  Min(alpha1) = 0  Max(alpha1) = 1.00000453

Applying the previous iteration compression flux

MULES: Correcting alpha.water

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.727573605  Min(alpha1) = 0  Max(alpha1) = 1.00000453

smoothSolver:  Solving for alpha.water, Initial residual = 9.37249465e-07, Final residual = 1.43841549e-13, No Iterations 2

Phase-1 volume fraction = 0.727573605  Min(alpha1) = 0  Max(alpha1) = 1.00000453

Applying the previous iteration compression flux

MULES: Correcting alpha.water

MULES: Correcting alpha.water

MULES: Correcting alpha.water

Phase-1 volume fraction = 0.727573605  Min(alpha1) = 0  Max(alpha1) = 1.00000453

GAMG:  Solving for p_rgh, Initial residual = 0.000180148101, Final residual = 1.65743017e-07, No Iterations 8

GAMG:  Solving for p_rgh, Initial residual = 1.42465044e-05, Final residual = 5.41884116e-07, No Iterations 1

time step continuity errors : sum local = 2.37251841e-11, global = -4.43079157e-13, cumulative = -1.07886437e-09

smoothSolver:  Solving for omega, Initial residual = 6.67536969e-06, Final residual = 1.77979817e-08, No Iterations 1

smoothSolver:  Solving for k, Initial residual = 4.46870005e-06, Final residual = 4.20645118e-08, No Iterations 1

ExecutionTime = 35840.95 s  ClockTime = 36006 s



Courant Number mean: 0.0159346502 max: 0.497184392

Interface Courant Number mean: 0.000190676158 max: 0.050893189

deltaT = 7.89779808e-05

Time = 1.26461786

setFields was buggy if I use in the quadrant other than first. Only in x and y axis both positive, setFields was working.