|
[Sponsors] |
June 18, 2011, 10:00 |
Foam::error::PrintStack
|
#1 |
New Member
Join Date: May 2011
Location: Germany/Bosnia
Posts: 9
Rep Power: 15 |
hi,
i have following errormessage in OpenFoam, as solver I use BuoyantSimpleFoam. I donīt understand that error. Maybe someone can help me? almir@ubuntu:~/OpenFOAM/zylinder$ buoyantSimpleFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 1.7.x | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 1.7.x-3776603e4c6c Exec : buoyantSimpleFoam Date : Jun 15 2011 Time : 12:26:48 Host : ubuntu PID : 5430 Case : /home/almir/OpenFOAM/zylinder nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading g Reading thermophysical properties Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>> Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting RAS turbulence model kOmegaSST kOmegaSSTCoeffs { alphaK1 0.85034; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.85616; Prt 1; gamma1 0.5532; gamma2 0.4403; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; c1 10; } Calculating field g.h Reading field p_rgh Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00987294, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0157, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00987846, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.0105804, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.899378, Final residual = 0.00305647, No Iterations 4 time step continuity errors : sum local = 18.9686, global = -1.40909e-15, cumulative = -1.40909e-15 rho max/min : 1.22108 1.13449 DILUPBiCG: Solving for omega, Initial residual = 0.999913, Final residual = 0.0105502, No Iterations 2 bounding omega, min: -902.694 max: 24331.5 average: 741.476 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0973743, No Iterations 1 bounding k, min: -0.000335494 max: 0.0031317 average: 0.000848312 ExecutionTime = 0.11 s ClockTime = 0 s Time = 2 DILUPBiCG: Solving for Ux, Initial residual = 0.11425, Final residual = 0.00308422, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.0715711, Final residual = 2.95363e-05, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.109402, Final residual = 0.00196119, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 0.177633, Final residual = 0.00377685, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.997056, Final residual = 0.00735437, No Iterations 4 time step continuity errors : sum local = 11.1294, global = 4.01181e-15, cumulative = 2.60272e-15 rho max/min : 309747 -321318 DILUPBiCG: Solving for omega, Initial residual = 0.594095, Final residual = 0.0340323, No Iterations 1 bounding omega, min: -7.04999e+17 max: 1.99134e+09 average: -1.78979e+14 DILUPBiCG: Solving for k, Initial residual = 0.999984, Final residual = 0.0558935, No Iterations 2 bounding k, min: -7.30668e+08 max: 1.51545e+08 average: -1.226e+06 ExecutionTime = 0.15 s ClockTime = 0 s Time = 3 DILUPBiCG: Solving for Ux, Initial residual = 0.125933, Final residual = 0.0012217, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.097672, Final residual = 0.000856571, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.130559, Final residual = 0.00122656, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 0.467341, Final residual = 0.00933279, No Iterations 1 #0 Foam::error::PrintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #2 in "/lib/libc.so.6" #3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so" #8 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/buoyantSimpleFoam" #9 __libc_start_main in "/lib/libc.so.6" #10 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/buoyantSimpleFoam" Gleitkomma-Ausnahme almir@ubuntu:~/OpenFOAM/zylinder$ greets almir |
|
June 18, 2011, 16:53 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Greetings Almir,
At the risk of sending you off in the wrong direction, you can try this answer: My program stops with an output that starts with #0 Foam::error:: PrintStack(Foam::Ostream&) But in an attempt to send you in the right direction:
Best regards, Bruno
__________________
Last edited by wyldckat; June 18, 2011 at 16:56. Reason: improve legibility ... |
|
February 9, 2012, 07:22 |
|
#4 |
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 16 |
Hi Foamers,
I am running twoPhaseEulerFoam and i have increased the mesh size 6000(which was in tutorial bed2) to 24000 and I am getting following error. I tried for 12000 again same.in blockMeshDict it was (30 200 1) first I have changed it to (30 400 1) then (60 400 2) an so on.Another problem is that I have to change the file 0/alpha everytime.is there any other practical solution for that? Courant Number mean: 0.263255 max: 12.2832 Max Ur Courant Number = 3.77181e+06 Time = 0.071 DILUPBiCG: Solving for alpha, Initial residual = 1.1014e-05, Final residual = 6.17658e-11, No Iterations 33 Dispersed phase volume fraction = 0.3 Min(alpha) = -1.92847 Max(alpha) = 2.81043 DILUPBiCG: Solving for alpha, Initial residual = 0.00010103, Final residual = 5.2889e-11, No Iterations 8 Dispersed phase volume fraction = 0.3 Min(alpha) = -0.247369 Max(alpha) = 1.92082 kinTheory: max(Theta) = 1000 kinTheory: min(nua) = 1.3774e-12, max(nua) = 0.0231854 kinTheory: min(pa) = -9295.94, max(pa) = 1.14803e+09 GAMG: Solving for p, Initial residual = 0.996287, Final residual = 0.0429939, No Iterations 1 time step continuity errors : sum local = 201567, global = 28.235, cumulative = 28.235 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam" Floating point exception Please help me...sorry for stupid questions thanksss a lot!!! |
|
February 9, 2012, 09:17 |
|
#5 |
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 16 |
Hi all,
thanks I did it alone... Thanks... |
|
February 17, 2012, 21:33 |
|
#6 |
Member
Elh. A2. BAH
Join Date: Jan 2012
Posts: 64
Rep Power: 14 |
Hello all,
I am having a similar error in using pimpleDyMFoam. Below is the error output: ------------------------------------ Courant Number mean: 0.00997899 max: 0.807965 deltaT = 1.32295e-104 --> FOAM Warning : From function Time:perator++() in file db/Time/Time.C at line 982 Increased the timePrecision from 267 to 268 to distinguish between timeNames at time 1.97982e-05 Time = 1.979823486337861903608045799352055382769322022795 67718505859375e-05 solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 1.97982e-05 transformation: ((0 0 0) (1 (0 0 0.000103663))) AMI: Creating addressing and weights between 16 source faces and 16 target faces AMI: Patch source weights min/max/average = 1, 1.0007, 1.00035 AMI: Patch target weights min/max/average = 0.986951, 0.987248, 0.987099 smoothSolver: Solving for Ux, Initial residual = 0.140328, Final residual = 5.75579e-08, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 0.140962, Final residual = 5.70666e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 0.814891, Final residual = 0.00591061, No Iterations 3 time step continuity errors : sum local = 0.00031744, global = 5.68686e-06, cumulative = 0.00082973 #0 Foam::error:rintStack(Foam::Ostream&) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/pimpleDyMFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/pimpleDyMFoam" -------------------------- Can you tell how you went along in solving the issue? Thank you for your help. |
|
February 18, 2012, 05:56 |
|
#7 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Greetings ebah6,
Quote:
Bruno
__________________
|
||
February 18, 2012, 16:54 |
|
#8 |
Member
Elh. A2. BAH
Join Date: Jan 2012
Posts: 64
Rep Power: 14 |
Thank you Bruno.
I appreciate. Let me go through this and to see how I can correct my mistakes. I will probably get back to you for more help. My best regards. |
|
April 3, 2012, 00:20 |
|
#9 |
Member
Elh. A2. BAH
Join Date: Jan 2012
Posts: 64
Rep Power: 14 |
Hello Bruno and everyone else,
Allow that I follow up on this thread for I am experience similar issues as those for which the thread was initiated. My log file is as follows: PHP Code:
Could you please have a look at this issue. Thanks in advance. |
|
April 3, 2012, 03:32 |
|
#10 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Look at the courant number: it's increasing a lot. I think you should reduce the `maxCo' value in the system/controlDict file.
The solution will be slower, but I think it'll works, |
|
April 3, 2012, 05:45 |
|
#11 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
To add to samiam1000's answer:
Best regards, Bruno
__________________
Last edited by wyldckat; April 3, 2012 at 05:52. Reason: see "edit:" |
|
April 3, 2012, 05:55 |
|
#12 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Bruno, Dear All,
thanks for the links that you added. I think they are very useful. Also, I do have a problem with buoyantPimpleFoam. I am trying to impose either temerature or velocity in some cells. I did the same with buoyanSimpleFoam, but I have problems with the steady state. Could you check the folder I modified, please? I am attaching the latest version of my solver, here. Thanks a lot, Samuele |
|
April 4, 2012, 19:14 |
|
#13 | |
Member
Elh. A2. BAH
Join Date: Jan 2012
Posts: 64
Rep Power: 14 |
Quote:
I did some dummy test cases: 1) a square box that rotates with the AMI; structured mesh 2) same thing with unstructured mesh. Both these cases don't seem to give any error. 3) I did my learning case with unstructured mesh; it consists on two cylindrical rotors (Darrieus). But here I run into trouble with the problem described above. Attached are some pictures to see how the meshes look like. For the latter case, you can see the velocity field is messing in the outer domain where no rotation is happening. Another question I had is how to export hybrid mesh from pointwise to openfoam? By hybrid I mean unstructured in x-y place and we extrude in the z-direction which will then be structured. I tried that but only the structured boundary faces are exported not the unstructured ones. Thank you for your and my best regards. |
||
April 5, 2012, 09:03 |
|
#14 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Hi ebah6,
Mmm... I'm not an expert on this subject, but this is what I can see that might be the source of the problems:
The other theory is that the thickness of the paddles is having a very big effect on the development of vortexes... and if these are not properly solved, it's only natural that some seriously crazy "fluid pressure shocks" (not a very technical term) will occur. Another issue might be the speed at which the rotor is running. Proper field initialization might be required to induce the solver to start with good starting values; otherwise, you probably will have to simulate starting with the rotation speed at 0 RPM. I'm not very familiar with these solvers, but my guess is that if you only want to have an "averaging" result, then one of the LTS solvers might come in handy... although you would have to create one that would LTS with AMI... Best regards, Bruno
__________________
|
|
April 10, 2012, 18:02 |
|
#15 |
Member
Elh. A2. BAH
Join Date: Jan 2012
Posts: 64
Rep Power: 14 |
Hello,
Yes Bruno, some of the possible issues are as you mentioned; thanks for your insight. In particular, as the body becomes thinner, I run into problems. However, I am only encountering problems when using a turbulence model: the laminar case runs fine. For the cases with a turbulent flow, I refined the mesh again and again but it still crashes with a skyrocketing Courant Number. My pressing issue is that I need to deal with thin bodies, so I need to a work-around. Also, you suggested the STL snappyHexMesh. I did that in a recent past but the sliding interfaces show step like shape dispite the refinement. Thanks for your help. |
|
April 11, 2012, 04:38 |
|
#16 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Hi ebah6,
Mmm, if it's not the mesh, then you've got to start tuning the "fvSchemes" file and possibly the "fvSolution" one as well. Unfortunately I don't know much how to configure them properly for each scenario, so I suggest that you check all of the relevant tutorials in OpenFOAM, as well as the User Guide. Good luck! Bruno
__________________
|
|
May 31, 2012, 12:11 |
Similar Quandary
|
#17 |
Member
Edward Leonard
Join Date: May 2012
Location: Calumet, MI
Posts: 31
Rep Power: 14 |
Good Afternoon, Everyone!
I come to this place with a similar issue, and upon reading the above comments and filtering through the User Guide for more information about initial conditions for k-epsilon and about the Courant number, I'm still having a heck of a time performing a run. Let me explain the situation to you (I can't post the 0/ files for various reasons): A 7m long blunt object is situated in a 10m/s wind-tunnel, with the floor of the tunnel moving with it (so we're in the blunt object's reference frame). The "ground" is of species 1 (alpha1), the wind-tunnel (or atmosphere) is of species 2 (alpha2) and the blunt object is spewing species 3 out of its' side at 40m/s (alpha3). So, in my back-of-the-envelope calculations (inspired by the User Guide), I set my initial value of k=2.5 and epsilon=0.25. I also set the initial time-step to 0.005sec, and for the sake of early testing I turned off "adjustTimeStep." With all of that said and done, it only completes one iteration of calculation, the output for which is here: Code:
Courant Number mean: 0.256274 max: 3.73333 PIMPLE: Operating solver in PISO mode time step continuity errors : sum local = 0.000307692, global = -1.74165e-05, cumulative = -1.74165e-05 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.01496e-11, No Iterations 589 time step continuity errors : sum local = 4.60233e-10, global = -4.49286e-16, cumulative = -1.74165e-05 Starting time loop Courant Number mean: 0.27796 max: 9.09506 Interface Courant Number mean: 0 max: 0 Time = 0.005 diagonal: Solving for alpha1, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for alpha2, Initial residual = 1, Final residual = 4.83778e-07, No Iterations 1 Air phase volume fraction = 0 Min(alpha1) = 0 Max(alpha1) = 1 Liquid phase volume fraction = 1 Min(alpha2) = 1 Max(alpha2) = 1 diagonal: Solving for alpha1, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for alpha2, Initial residual = 0.354522, Final residual = 1.43335e-07, No Iterations 1 Air phase volume fraction = 0 Min(alpha1) = 0 Max(alpha1) = 1 Liquid phase volume fraction = 1 Min(alpha2) = 1 Max(alpha2) = 1 DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0479761, No Iterations 305 time step continuity errors : sum local = 2.33578, global = 2.22882, cumulative = 2.2288 [kaleva:14255] 5 more processes have sent help message help-mpi-common-sm.txt / mmap on nfs [kaleva:14255] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages DICPCG: Solving for p_rgh, Initial residual = 0.683683, Final residual = 0.0337928, No Iterations 320 time step continuity errors : sum local = 3.34528, global = 2.22882, cumulative = 4.45762 DICPCG: Solving for p_rgh, Initial residual = 0.473874, Final residual = 9.50935e-08, No Iterations 467 time step continuity errors : sum local = 2.35844, global = 2.22882, cumulative = 6.68644 DILUPBiCG: Solving for epsilon, Initial residual = 0.999999, Final residual = 25.758, No Iterations 1001 bounding epsilon, min: -7.75298e+11 max: 1.02524e+12 average: -256836 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 89.7085, No Iterations 1001 bounding k, min: -7.39749e+11 max: 6.80932e+11 average: -9.18358e+06 time step continuity errors : sum local = 2.35844, global = 2.22882, cumulative = 8.91526 ExecutionTime = 22.39 s ClockTime = 23 s Courant Number mean: 611.356 max: 4.60701e+06 Interface Courant Number mean: 0 max: 0 Time = 0.01 diagonal: Solving for alpha1, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for alpha2, Initial residual = 1, Final residual = 1.1291e-10, No Iterations 1 Air phase volume fraction = -1.11442 Min(alpha1) = -339.921 Max(alpha1) = 1 Liquid phase volume fraction = -0.114296 Min(alpha2) = -4.60701e+06 Max(alpha2) = 614.595 [0] #0 Foam::error::printStack(Foam::Ostream&) in "/home/leonard/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #1 Foam::sigFpe::sigHandler(int) in "/home/leonard/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #2 __restore_rt at sigaction.c:0 [0] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double,Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/home/leonard/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interMixingFoam" I would much appreciate any input anyone has on this matter. Thanks! |
|
May 31, 2012, 14:32 |
|
#18 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Greetings Edward,
OK, if you've read about the Courant number, then you should know that you should check the smallest cell size you've got:
Oh, and if checkMesh tells you that you've got bad cells or faces, then that's another source of your problems Best regards, Bruno
__________________
|
|
June 1, 2012, 09:40 |
|
#19 | |
Member
Edward Leonard
Join Date: May 2012
Location: Calumet, MI
Posts: 31
Rep Power: 14 |
Quote:
However, I've since decided to go a different route because this takes a painful amount of time to process. I had interFoam running on this large mesh (see checkMesh output below) on 8 CPUs, left it over-night and it had only gotten to 0.08sec by the following morning. Since my goal is a steady-state solution, I think what I want to try is to add the phase mixing of interFoam to the SIMPLE solver of simpleFoam. I took a quick look at it yesterday, and it seems like it will be a formidable task. Any insight on the matter before I hit the ground running? Thanks! ~Ed |
||
June 1, 2012, 09:50 |
|
#20 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Hi Edward,
Mmm, you forgot to attach your checkMesh log. Anyway, if you want the steady state solution with interFoam, then probably this is what you want: http://www.openfoam.org/version2.0.0/steady-vof.php Best regards, Bruno
__________________
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
FoamerrorprintStack | mayank | OpenFOAM Running, Solving & CFD | 38 | November 25, 2011 22:58 |