CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Initialisation problem with interFoam (https://www.cfd-online.com/Forums/openfoam-programming-development/127040-initialisation-problem-interfoam.html)

belkadi December 2, 2013 05:52

Initialisation problem with interFoam
 
1 Attachment(s)
Dear Foamer,

I'm trying to use the interFoam to simulate the mixing of two fluids.
I've uploaded the test case:
1. When running the blockMesh there is no errors.
2. After that for runing the setFields (no errors) : please do not forget to do the cp0/alph1.org 0/alpha1 before runing the setFields
3. runing interFoam : I think that there is some issues with the initialisation or may be the mesh, but I could not figure out what it is
exaclty : the error message is :

----------------------------

Create mesh for time = 0
 
PIMPLE: Operating solver in PISO mode
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading g
Calculating field g.h
No finite volume options present
time step continuity errors : sum local = 1.47059e+292, global = 1.47059e+292, cumulative = 1.47059e+292
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 5.19085e-11, No Iterations 66
time step continuity errors : sum local = 7.63557e+281, global = -5.03459e+279, cumulative = 1.47059e+292
Courant Number mean: 2.23052e+293 max: 1.58306e+294
Starting time loop
Courant Number mean: 0.0704495 max: 0.5
Interface Courant Number mean: 0 max: 0
deltaT = 3.15844e-298
Time = 3.15844e-298
MULES: Solving for alpha1
Phase-1 volume fraction = 0.939827 Min(alpha1) = -0.114728 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.938477 Min(alpha1) = -0.229443 Max(alpha1) = 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 Uninterpreted:
 
---------------------------------------
As you can see the sum local is too high, I think taht in the beginig of calculation some dividing by zero occured
Honestly, I don't know how to fix it, do you have any ideas
it would be really appreciated !

Many thanks

Krimo


All times are GMT -4. The time now is 23:59.