CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

problem with InterFOAM interface courant number

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 11, 2016, 10:29
Default problem with InterFOAM interface courant number
  #1
Member
 
Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Yage is on a distinguished road
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
Yage is offline   Reply With Quote

Old   January 11, 2016, 17:06
Default
  #2
New Member
 
shashank moghe
Join Date: Feb 2015
Posts: 29
Rep Power: 3
smog is on a distinguished road
Looks like there is no interface. Can you post your setFieldsDict?
smog is offline   Reply With Quote

Old   January 12, 2016, 03:07
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,014
Rep Power: 18
akidess will become famous soon enough
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
akidess is offline   Reply With Quote

Old   January 12, 2016, 04:03
Default
  #4
Member
 
Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Yage is on a distinguished road
Quote:
Originally Posted by smog View Post
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)
        );
    }
Yage is offline   Reply With Quote

Old   January 12, 2016, 04:17
Default
  #5
Member
 
Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Yage is on a distinguished road
Quote:
Originally Posted by akidess View Post
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 is offline   Reply With Quote

Old   January 12, 2016, 04:39
Default
  #6
Member
 
Yage
Join Date: May 2014
Posts: 42
Rep Power: 3
Yage is on a distinguished road
Quote:
Originally Posted by akidess View Post
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
Yage is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 45 February 8, 2016 05:42
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 10:53
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 10 June 16, 2015 15:48
Waterwheel shaped turbine inside a pipe simulation problem mshahed91 CFX 3 January 10, 2015 12:19
DecomposePar unequal number of shared faces maka OpenFOAM Pre-Processing 6 August 12, 2010 09:01


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