[Yage]

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

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

[smog]

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

[akidess]

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?

[Yage]

Originally Posted by smog

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

Dear Shashank,

The setFieldsDict are as follow:

regions
(
    boxToCell
    {
        box     (-1 -1 -1.5) (1.55  1 0.02);
        fieldValues
        (
            volScalarFieldValue alpha.water 1
            volVectorFieldValue U (2 0 0)
        );
    }

[Yage]

Originally Posted by akidess

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.

[Yage]

Originally Posted by akidess

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.

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

[venkcy]

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