CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Foam::error (https://www.cfd-online.com/Forums/openfoam/190481-foam-error.html)

yangzhuan July 15, 2017 00:21
Foam::error
 
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)
Could someone please help me understand.
Thank you

saddy July 15, 2017 02:24
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

yangzhuan July 15, 2017 03:21
3 Attachment(s)
Quote:
Originally Posted by saddy (Post 657159)
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.
    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
Please help me see where is wrong and how to modify.
Best

saddy July 17, 2017 12:18
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

yangzhuan July 18, 2017 04:50
Quote:
Originally Posted by saddy (Post 657407)
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

yangzhuan July 30, 2017 07:03
Quote:
Originally Posted by saddy (Post 657407)
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:confused:
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"

saddy July 30, 2017 10:45
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.

saddy August 2, 2017 01:23
i ran your case....i used openfoam v1612+
it was running like a charm without any problem

yangzhuan August 3, 2017 22:53
Quote:
Originally Posted by saddy (Post 659199)
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:rolleyes:.
Best.

yangzhuan August 3, 2017 23:18
Hi,
I forgot to say that I used openfoam 3.0.0. Maybe it's a version problem?

saddy August 4, 2017 10:20
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

yangzhuan August 5, 2017 00:50
Quote:
Originally Posted by saddy (Post 659505)
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.:D


All times are GMT -4. The time now is 06:13.