# problem with InterFOAM interface courant number

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 11, 2016, 10:29 problem with InterFOAM interface courant number #1 Member   Yage Join Date: May 2014 Posts: 42 Rep Power: 3 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```

 January 11, 2016, 17:06 #2 New Member   shashank moghe Join Date: Feb 2015 Posts: 29 Rep Power: 3 Looks like there is no interface. Can you post your setFieldsDict?

 January 12, 2016, 03:07 #3 Senior Member     Anton Kidess Join Date: May 2009 Location: Delft, Netherlands Posts: 1,014 Rep Power: 18 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? __________________ *On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. *Join the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oHPxcPqde7HtA2

January 12, 2016, 04:03
#4
Member

Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Quote:
 Originally Posted by smog 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)
);
}```

January 12, 2016, 04:17
#5
Member

Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Quote:
 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.

January 12, 2016, 04:39
#6
Member

Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Quote:
 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.

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```

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Saturn CFX 45 February 8, 2016 05:42 Saima CFX 45 September 22, 2015 10:53 mpeppels CFX 10 June 16, 2015 15:48 mshahed91 CFX 3 January 10, 2015 12:19 maka OpenFOAM Pre-Processing 6 August 12, 2010 09:01

All times are GMT -4. The time now is 16:35.