
[Sponsors] 
July 15, 2017, 00:21 
Foam::error

#1 
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3 
Hi,
I simulate the watergas 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.42333e67 Time = 0.218749 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 8.31475e06, Final residual = 4.80174e10, No Iterations 1 Phase1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.842399 Min(alpha.water) = 7.82983e22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.31197e06, Final residual = 4.81265e10, No Iterations 1 Phase1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.842399 Min(alpha.water) = 7.836e22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.30939e06, Final residual = 4.82519e10, No Iterations 1 Phase1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.842399 Min(alpha.water) = 7.836e22 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 8.3066e06, Final residual = 4.83629e10, No Iterations 1 Phase1 volume fraction = 0.842399 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.842399 Min(alpha.water) = 7.84216e22 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.19752e07, global = 1.38624e11, cumulative = 5.61794e05 DICPCG: Solving for p_rgh, Initial residual = 0.0180542, Final residual = 7.93864e08, No Iterations 196 time step continuity errors : sum local = 1.10086e11, global = 2.96969e14, cumulative = 5.61794e05 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #2 ? in "/lib64/libc.so.6" #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #4 double Foam::gSumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&, int) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #5 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #8 Foam::fvMatrix<double>::solve() in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" #10 Foam::RASModels::realizableKE<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" #11 ? in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #12 __libc_start_main in "/lib64/libc.so.6" #13 ? in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" 浮点数例外 (core dumped) Thank you 

July 15, 2017, 02:24 

#2 
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 114
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: 26
Rep Power: 3 
Quote:
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. Cell to face addressing 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 (nonclosed singly connected) in2 108 126 ok (nonclosed singly connected) bimian 6766 6800 ok (nonclosed singly connected) out 276 294 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (0 0.0249969 0.025) (1 0.0249969 0.025) Mesh has 3 geometric (nonempty/wedge) directions (1 1 1) Mesh has 3 solution (nonempty) directions (1 1 1) Boundary openness (1.21207e16 7.78735e17 1.49379e16) OK. Max cell openness = 1.42787e15 OK. Max aspect ratio = 45.477 OK. Minimum face area = 3.31589e07. Maximum face area = 2.60875e05. Face area magnitudes OK. Min volume = 1.66628e09. Max volume = 1.16271e07. Total volume = 0.00195001. Cell volumes OK. Mesh nonorthogonality Max: 32.8264 average: 8.01901 Nonorthogonality check OK. Face pyramids OK. Max skewness = 1.31682 OK. Coupled point location match (average 0) OK. Mesh OK. End Best 

July 17, 2017, 12:18 

#4 
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 114
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: 26
Rep Power: 3 
Quote:
Best 

July 30, 2017, 07:03 

#6  
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3 
Quote:
Code:
Courant Number mean: 3.40147e05 max: 0.502766 Interface Courant Number mean: 1.48331e05 max: 0.0517021 deltaT = 4.061e68 Time = 4.71186 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 2.03688e06, Final residual = 2.93968e11, No Iterations 1 Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 2.03667e06, Final residual = 2.93849e11, No Iterations 1 Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 2.03639e06, Final residual = 2.93744e11, No Iterations 1 Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 volume fraction = 0.956421 Min(alpha.water) = 2.21102e28 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 2.03613e06, Final residual = 2.93615e11, No Iterations 1 Phase1 volume fraction = 0.956421 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase1 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.43246e10, global = 1.84738e10, cumulative = 2.56315e06 DICPCG: Solving for p_rgh, Initial residual = 0.00538723, Final residual = 6.52414e08, No Iterations 224 time step continuity errors : sum local = 1.02756e13, global = 6.67559e16, cumulative = 2.56315e06 [3] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" [3] #1 Foam::sigFpe::sigHandler(int) in "/opt/software/OpenFOAM/OpenFOAMv3.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/OpenFOAMv3.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/OpenFOAMv3.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/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" [3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so" [3] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" [3] #8 Foam::fvMatrix<double>::solve() in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" [3] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" [3] #10 Foam::RASModels::realizableKE<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so" [3] #11 ? in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" [3] #12 __libc_start_main in "/lib64/libc.so.6" [3] #13 ? in "/opt/software/OpenFOAM/OpenFOAMv3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" 

July 30, 2017, 10:45 

#7 
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 114
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: 114
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: 26
Rep Power: 3 
Quote:
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: 26
Rep Power: 3 
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: 114
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: 26
Rep Power: 3 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
mpirun interFoam very stable > then blows up  ghadab  OpenFOAM Running, Solving & CFD  3  October 27, 2013 10:34 