CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Foam::error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 15, 2017, 00:21
Default Foam::error
  #1
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
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
yangzhuan is offline   Reply With Quote

Old   July 15, 2017, 02:24
Default
  #2
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 114
Rep Power: 4
saddy is on a distinguished road
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
saddy is offline   Reply With Quote

Old   July 15, 2017, 03:21
Default
  #3
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Quote:
Originally Posted by saddy View Post
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
Attached Files
File Type: zip 0.zip (8.2 KB, 3 views)
File Type: zip system.zip (7.0 KB, 2 views)
File Type: zip constant.zip (3.0 KB, 0 views)
yangzhuan is offline   Reply With Quote

Old   July 17, 2017, 12:18
Default
  #4
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 114
Rep Power: 4
saddy is on a distinguished road
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
saddy is offline   Reply With Quote

Old   July 18, 2017, 04:50
Default
  #5
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Quote:
Originally Posted by saddy View Post
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 is offline   Reply With Quote

Old   July 30, 2017, 07:03
Default
  #6
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Quote:
Originally Posted by saddy View Post
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"
yangzhuan is offline   Reply With Quote

Old   July 30, 2017, 10:45
Default
  #7
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 114
Rep Power: 4
saddy is on a distinguished road
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 is offline   Reply With Quote

Old   August 2, 2017, 01:23
Default
  #8
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 114
Rep Power: 4
saddy is on a distinguished road
i ran your case....i used openfoam v1612+
it was running like a charm without any problem
saddy is offline   Reply With Quote

Old   August 3, 2017, 22:53
Default
  #9
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Quote:
Originally Posted by saddy View Post
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.
yangzhuan is offline   Reply With Quote

Old   August 3, 2017, 23:18
Default
  #10
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Hi,
I forgot to say that I used openfoam 3.0.0. Maybe it's a version problem?
yangzhuan is offline   Reply With Quote

Old   August 4, 2017, 10:20
Default
  #11
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 114
Rep Power: 4
saddy is on a distinguished road
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
saddy is offline   Reply With Quote

Old   August 5, 2017, 00:50
Default
  #12
New Member
 
Alice
Join Date: Jun 2017
Posts: 26
Rep Power: 3
yangzhuan is on a distinguished road
Quote:
Originally Posted by saddy View Post
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.
yangzhuan is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 02:54.