|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
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<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so" #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" #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" #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #8 Foam::fvMatrix<double>::solve() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam" #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" #10 Foam::RASModels::realizableKE<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::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) Thank you |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 11 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
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 (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 |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 11 ![]() |
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
|
|
![]() |
![]() |
![]() |
![]() |
#5 | |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
Quote:
Best |
||
![]() |
![]() |
![]() |
![]() |
#6 | |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
Quote:
![]() 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" |
||
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 11 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#8 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 11 ![]() |
i ran your case....i used openfoam v1612+
it was running like a charm without any problem |
|
![]() |
![]() |
![]() |
![]() |
#9 | |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
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. |
||
![]() |
![]() |
![]() |
![]() |
#10 |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
Hi,
I forgot to say that I used openfoam 3.0.0. Maybe it's a version problem? |
|
![]() |
![]() |
![]() |
![]() |
#11 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 11 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#12 |
New Member
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 9 ![]() |
||
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
mpirun interFoam very stable --> then blows up | ghadab | OpenFOAM Running, Solving & CFD | 3 | October 27, 2013 11:34 |