# Foam::error

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

 July 15, 2017, 00:21 Foam::error #1 New Member   Alice Join Date: Jun 2017 Posts: 25 Rep Power: 2 Hi, I simulate the water-gas two phase using interFoam. But i have a error message when run interFoam. I don't understand the error. Code: ```Courant Number mean: 0.000593557 max: 0.503771 Interface Courant Number mean: 0.00022715 max: 0.296666 deltaT = 4.42333e-67 Time = 0.218749 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 8.31475e-06, Final residual = 4.80174e-10, No Iterations 1 Phase-1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.842399 Min(alpha.water) = -7.82983e-22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.31197e-06, Final residual = 4.81265e-10, No Iterations 1 Phase-1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.842399 Min(alpha.water) = -7.836e-22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.30939e-06, Final residual = 4.82519e-10, No Iterations 1 Phase-1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.842399 Min(alpha.water) = -7.836e-22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.3066e-06, Final residual = 4.83629e-10, No Iterations 1 Phase-1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.842399 Min(alpha.water) = -7.84216e-22 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 0.0571824, Final residual = 0.00157402, No Iterations 2 time step continuity errors : sum local = 2.19752e-07, global = 1.38624e-11, cumulative = -5.61794e-05 DICPCG: Solving for p_rgh, Initial residual = 0.0180542, Final residual = 7.93864e-08, No Iterations 196 time step continuity errors : sum local = 1.10086e-11, global = -2.96969e-14, cumulative = -5.61794e-05 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #2 ? in "/lib64/libc.so.6" #3 double Foam::sumProd(Foam::UList const&, Foam::UList const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #4 double Foam::gSumProd(Foam::UList const&, Foam::UList const&, int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #5 Foam::PBiCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so" #7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #8 Foam::fvMatrix::solve() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #9 Foam::SolverPerformance Foam::solve(Foam::tmp > const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" #10 Foam::RASModels::realizableKE >::correct() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" #11 ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #12 __libc_start_main in "/lib64/libc.so.6" #13 ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" 浮点数例外 (core dumped)``` Could someone please help me understand. Thank you

 July 15, 2017, 02:24 #2 Senior Member   sandy Join Date: Feb 2016 Location: . Posts: 113 Rep Power: 4 clearly your boundary conditions are to be blamed. this occurs usually bcz of 2 reasons 1. compile a solver for openfoam 2.4.0 in another openfoam say 4.0 so it compiles but its still isn't working correctly 2. when u give incorrect unphysical boundary conditions. hence the solution proceeds initially, but later the problem surfaces and solution is difficult to attain hence. your whole calculation is dumped i believe you haven't done first error but you have committed second one. so rectify your bc. thats it AND OH GOD. HOW COME YOUR deltaT is order of 10^-67 ......total unphysical setup.....better upload your case files here

July 15, 2017, 03:21
#3
New Member

Alice
Join Date: Jun 2017
Posts: 25
Rep Power: 2
Quote:
 Originally Posted by saddy clearly your boundary conditions are to be blamed. this occurs usually bcz of 2 reasons 1. compile a solver for openfoam 2.4.0 in another openfoam say 4.0 so it compiles but its still isn't working correctly 2. when u give incorrect unphysical boundary conditions. hence the solution proceeds initially, but later the problem surfaces and solution is difficult to attain hence. your whole calculation is dumped i believe you haven't done first error but you have committed second one. so rectify your bc. thats it AND OH GOD. HOW COME YOUR deltaT is order of 10^-67 ......total unphysical setup.....better upload your case files here
Hi,
Thank you for your reply. I think my bc is ok, and I couldn't see any wrong in my case. I upload my case here. The polyMesh file should have be in constant file. But the file exceeds the forum's limit of 195.3KB. So I can't upload the ployMesh file. I can post the checkMesh here.
Code:
```Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points:           58800
faces:            168431
internal faces:   161113
cells:            54924
faces per cell:   6
boundary patches: 4
point zones:      0
face zones:       0
cell zones:       0
Overall number of cells of each type:
hexahedra:     54924
prisms:        0
wedges:        0
pyramids:      0
tet wedges:    0
tetrahedra:    0
polyhedra:     0
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch               Faces    Points   Surface topology
in1                 168      189      ok (non-closed singly connected)
in2                 108      126      ok (non-closed singly connected)
bimian              6766     6800     ok (non-closed singly connected)
out                 276      294      ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (0 -0.0249969 -0.025) (1 0.0249969 0.025)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.21207e-16 7.78735e-17 -1.49379e-16) OK.
Max cell openness = 1.42787e-15 OK.
Max aspect ratio = 45.477 OK.
Minimum face area = 3.31589e-07. Maximum face area = 2.60875e-05.  Face area magnitudes OK.
Min volume = 1.66628e-09. Max volume = 1.16271e-07.  Total volume = 0.00195001.  Cell volumes OK.
Mesh non-orthogonality Max: 32.8264 average: 8.01901
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.31682 OK.
Coupled point location match (average 0) OK.
Mesh OK.
End```
Best
Attached Files
 0.zip (8.2 KB, 3 views) system.zip (7.0 KB, 2 views) constant.zip (3.0 KB, 0 views)

 July 17, 2017, 12:18 #4 Senior Member   sandy Join Date: Feb 2016 Location: . Posts: 113 Rep Power: 4 sorry for late reply. i have been very busy i haven.t seen all of your files. i saw ur alpha files and i don;t get it why u using inletoutlet just use fixedvalue and use value 0 or 1 as per your phase. try this hopefuly it will work

July 18, 2017, 04:50
#5
New Member

Alice
Join Date: Jun 2017
Posts: 25
Rep Power: 2
Quote:
 Originally Posted by saddy sorry for late reply. i have been very busy i haven.t seen all of your files. i saw ur alpha files and i don;t get it why u using inletoutlet just use fixedvalue and use value 0 or 1 as per your phase. try this hopefuly it will work
Thank you for your reply. I have solved. The main reason is that gas and liquid speed setting is not reasonable.
Best

July 30, 2017, 07:03
#6
New Member

Alice
Join Date: Jun 2017
Posts: 25
Rep Power: 2
Quote:
 Originally Posted by saddy sorry for late reply. i have been very busy i haven.t seen all of your files. i saw ur alpha files and i don;t get it why u using inletoutlet just use fixedvalue and use value 0 or 1 as per your phase. try this hopefuly it will work
Hi, I have confused. After the mesh is refined, this error occurs again. I tried different gas-liquid speeds and still couldn't solve it. I really don't know how to solve it
Code:
```Courant Number mean: 3.40147e-05 max: 0.502766
Interface Courant Number mean: 1.48331e-05 max: 0.0517021
deltaT = 4.061e-68
Time = 4.71186
PIMPLE: iteration 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03688e-06, Final residual = 2.93968e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03667e-06, Final residual = 2.93849e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03639e-06, Final residual = 2.93744e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = -2.21102e-28  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03613e-06, Final residual = 2.93615e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 0.00709543, Final residual = 0.000345669, No Iterations 8
time step continuity errors : sum local = 5.43246e-10, global = 1.84738e-10, cumulative = 2.56315e-06
DICPCG:  Solving for p_rgh, Initial residual = 0.00538723, Final residual = 6.52414e-08, No Iterations 224
time step continuity errors : sum local = 1.02756e-13, global = -6.67559e-16, cumulative = 2.56315e-06
[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #2  ? in "/lib64/libc.so.6"
[3] #3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #4  double Foam::gSumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&, int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #5  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so"
[3] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #8  Foam::fvMatrix<double>::solve() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so"
[3] #10  Foam::RASModels::realizableKE<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so"
[3] #11  ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #12  __libc_start_main in "/lib64/libc.so.6"
[3] #13  ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"```

 July 30, 2017, 10:45 #7 Senior Member   sandy Join Date: Feb 2016 Location: . Posts: 113 Rep Power: 4 hey i read you post! i am still puzzled with your deta T = 10^-67 ??? Its just not physical plz send me your case files : saddyhasan@gmail.com i'll reply in tomorrow or later.

 August 2, 2017, 01:23 #8 Senior Member   sandy Join Date: Feb 2016 Location: . Posts: 113 Rep Power: 4 i ran your case....i used openfoam v1612+ it was running like a charm without any problem

August 3, 2017, 22:53
#9
New Member

Alice
Join Date: Jun 2017
Posts: 25
Rep Power: 2
Quote:
 Originally Posted by saddy i ran your case....i used openfoam v1612+ it was running like a charm without any problem
Oh, really?
I change the boundary conditions, and try it again and again. It still does not work. There are always the same mistakes(Foam::error...). I'll study it again and find out the mistake.
Best.

 August 3, 2017, 23:18 #10 New Member   Alice Join Date: Jun 2017 Posts: 25 Rep Power: 2 Hi, I forgot to say that I used openfoam 3.0.0. Maybe it's a version problem?

 August 4, 2017, 10:20 #11 Senior Member   sandy Join Date: Feb 2016 Location: . Posts: 113 Rep Power: 4 one thing you can try replace your fvsolution and fvschemes files with the ones in dambreak tutorial located in tutorials/multiphase/interfoam most probably it should work

August 5, 2017, 00:50
#12
New Member

Alice
Join Date: Jun 2017
Posts: 25
Rep Power: 2
Quote:
 Originally Posted by saddy one thing you can try replace your fvsolution and fvschemes files with the ones in dambreak tutorial located in tutorials/multiphase/interfoam most probably it should work
Thank you. I will try it.

 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