
[Sponsors] 
November 10, 2014, 00:52 

#41 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
@RodriguezFatz
Code:
Does the omegabounding vanish, if you use the laplaciansetting I suggested? Code:
Courant Number mean: 0.000113778 max: 1 deltaT = 1.41782e09 Time = 6.59824e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.59824e07 transformation: ((0 0 0) (0.999946 (0.0103643 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.474286, 1.00002, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.443604, 1.05059, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.519849, 1.00249, 0.999845 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00915686, No Iterations 54 GAMG: Solving for pcorr, Initial residual = 0.0243898, Final residual = 0.00464874, No Iterations 1 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.2097 0.276597 smoothSolver: Solving for Ux, Initial residual = 3.16377e05, Final residual = 2.14271e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00139322, Final residual = 1.18808e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00157917, Final residual = 1.27902e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000472102, Final residual = 3.83132e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.76418e05, Final residual = 2.45351e13, No Iterations 1 GAMG: Solving for p, Initial residual = 1.21131e09, Final residual = 1.21131e09, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171506, global = 0.000151276, cumulative = 0.0929682 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.15378e07, Final residual = 1.15378e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.15378e07, Final residual = 1.15378e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171515, global = 0.000151276, cumulative = 0.0931195 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.153e07, Final residual = 1.153e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.153e07, Final residual = 1.153e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171515, global = 0.000151276, cumulative = 0.0932708 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.11641e06, Final residual = 1.37612e10, No Iterations 1 bounding omega, min: 6204.38 max: 6.34422e+09 average: 2.18357e+06 smoothSolver: Solving for k, Initial residual = 2.41327e05, Final residual = 9.59219e10, No Iterations 1 ExecutionTime = 70114.9 s ClockTime = 94366 s regIOobject::readIfModified() : Rereading object fvSchemes from file "/home/eatin/OpenFOAM/eatin2.3.0/run/tutorials/TurboCharger/Trial_run/system/fvSchemes" Courant Number mean: 0.00011376 max: 1 deltaT = 1.41782e09 Time = 6.61242e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.61242e07 transformation: ((0 0 0) (0.999946 (0.0103866 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.474033, 1, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.439444, 1.05157, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.519845, 1.00248, 0.999845 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00996615, No Iterations 57 GAMG: Solving for pcorr, Initial residual = 1.31767e05, Final residual = 1.31767e05, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.02652 0.480629 smoothSolver: Solving for Ux, Initial residual = 3.16636e05, Final residual = 2.11836e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00140312, Final residual = 1.13106e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00159988, Final residual = 1.2634e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000464709, Final residual = 3.32125e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.73208e05, Final residual = 4.49914e13, No Iterations 1 GAMG: Solving for p, Initial residual = 4.48889e13, Final residual = 4.48889e13, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171907, global = 0.000151647, cumulative = 0.0934224 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16407e07, Final residual = 1.16407e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16407e07, Final residual = 1.16407e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171916, global = 0.000151647, cumulative = 0.0935741 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16329e07, Final residual = 1.16329e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16329e07, Final residual = 1.16329e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000171916, global = 0.000151647, cumulative = 0.0937257 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.129e06, Final residual = 2.07432e10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.40993e05, Final residual = 1.02225e09, No Iterations 1 ExecutionTime = 70251.8 s ClockTime = 94547 s Courant Number mean: 0.000113739 max: 1 deltaT = 1.41781e09 Time = 6.6266e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.6266e07 transformation: ((0 0 0) (0.999946 (0.0104089 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.473781, 1, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.435284, 1.05167, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.51984, 1.00247, 0.999844 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00919969, No Iterations 54 GAMG: Solving for pcorr, Initial residual = 1.21672e05, Final residual = 1.21672e05, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.03418 0.471993 smoothSolver: Solving for Ux, Initial residual = 3.16366e05, Final residual = 2.11287e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00140093, Final residual = 1.12678e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00159565, Final residual = 1.25571e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000461355, Final residual = 3.36557e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.71918e05, Final residual = 6.20407e13, No Iterations 1 GAMG: Solving for p, Initial residual = 6.18972e13, Final residual = 6.18972e13, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000172306, global = 0.000152017, cumulative = 0.0938777 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16604e07, Final residual = 1.16604e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16604e07, Final residual = 1.16604e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000172314, global = 0.000152017, cumulative = 0.0940298 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16525e07, Final residual = 1.16525e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16525e07, Final residual = 1.16525e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000172314, global = 0.000152017, cumulative = 0.0941818 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.12444e06, Final residual = 9.55552e11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.40609e05, Final residual = 9.28125e10, No Iterations 1 ExecutionTime = 70374.2 s ClockTime = 94719 s Courant Number mean: 0.000113719 max: 1 deltaT = 1.41781e09 Time = 6.64078e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.64078e07 transformation: ((0 0 0) (0.999946 (0.0104311 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.47353, 1.00001, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.431124, 1.05171, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.519832, 1.00246, 0.999844 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00988888, No Iterations 57 GAMG: Solving for pcorr, Initial residual = 1.30742e05, Final residual = 1.30742e05, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.02687 0.482841 smoothSolver: Solving for Ux, Initial residual = 3.1612e05, Final residual = 2.10315e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00139863, Final residual = 1.12206e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00159138, Final residual = 1.24932e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000457716, Final residual = 3.34466e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.70692e05, Final residual = 5.14521e13, No Iterations 1 GAMG: Solving for p, Initial residual = 5.13329e13, Final residual = 5.13329e13, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000172702, global = 0.000152384, cumulative = 0.0943342 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16259e07, Final residual = 1.16259e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16259e07, Final residual = 1.16259e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00017271, global = 0.000152384, cumulative = 0.0944865 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.16181e07, Final residual = 1.16181e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.16181e07, Final residual = 1.16181e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00017271, global = 0.000152384, cumulative = 0.0946389 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.12412e06, Final residual = 9.48272e11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.40277e05, Final residual = 9.26293e10, No Iterations 1 ExecutionTime = 70511.2 s ClockTime = 94911 s Courant Number mean: 0.000113698 max: 1 deltaT = 1.4178e09 Time = 6.65496e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.65496e07 transformation: ((0 0 0) (0.999945 (0.0104534 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.47328, 1, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.447838, 1.05176, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.519822, 1.00245, 0.99986 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00924126, No Iterations 54 GAMG: Solving for pcorr, Initial residual = 1.2222e05, Final residual = 1.2222e05, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.03417 0.472107 smoothSolver: Solving for Ux, Initial residual = 3.15871e05, Final residual = 2.09785e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00139647, Final residual = 1.11942e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.0015874, Final residual = 1.24392e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000453908, Final residual = 3.31894e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.69419e05, Final residual = 6.18263e13, No Iterations 1 GAMG: Solving for p, Initial residual = 6.16834e13, Final residual = 6.16834e13, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000173096, global = 0.000152751, cumulative = 0.0947917 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.15716e07, Final residual = 1.15716e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.15716e07, Final residual = 1.15716e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000173105, global = 0.000152751, cumulative = 0.0949444 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.15638e07, Final residual = 1.15638e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.15638e07, Final residual = 1.15638e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000173105, global = 0.000152751, cumulative = 0.0950972 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.12485e06, Final residual = 9.43144e11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.39952e05, Final residual = 9.24967e10, No Iterations 1 ExecutionTime = 70649.8 s ClockTime = 95106 s Courant Number mean: 0.000113677 max: 1.00001 deltaT = 1.41779e09 Time = 6.66913e07 solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 6.66913e07 transformation: ((0 0 0) (0.999945 (0.0104757 0 0))) AMI: Creating addressing and weights between 1900 source faces and 32076 target faces AMI: Patch source sum(weights) min/max/average = 0.995413, 1, 0.999762 AMI: Patch target sum(weights) min/max/average = 0.473032, 1.00001, 0.996806 AMI: Creating addressing and weights between 17748 source faces and 5456 target faces AMI: Patch source sum(weights) min/max/average = 0.447836, 1.0518, 1.00015 AMI: Patch target sum(weights) min/max/average = 0.51981, 1.00244, 0.99986 AMI: Creating addressing and weights between 17839 source faces and 1957 target faces AMI: Patch source sum(weights) min/max/average = 0.86998, 1, 0.999108 AMI: Patch target sum(weights) min/max/average = 0.991279, 1, 0.99992 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00988871, No Iterations 57 GAMG: Solving for pcorr, Initial residual = 1.30738e05, Final residual = 1.30738e05, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 rhoEqn max/min : 2.02717 0.48285 smoothSolver: Solving for Ux, Initial residual = 3.15639e05, Final residual = 2.09028e09, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.00139419, Final residual = 1.11562e07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00158337, Final residual = 1.23881e07, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.000450021, Final residual = 3.2803e08, No Iterations 1 GAMG: Solving for p, Initial residual = 6.68107e05, Final residual = 4.95164e13, No Iterations 1 GAMG: Solving for p, Initial residual = 4.94024e13, Final residual = 4.94024e13, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00017349, global = 0.000153117, cumulative = 0.0952503 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.15139e07, Final residual = 1.15139e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.15139e07, Final residual = 1.15139e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000173498, global = 0.000153117, cumulative = 0.0954034 rho max/min : 2 0.5 GAMG: Solving for p, Initial residual = 1.15061e07, Final residual = 1.15061e07, No Iterations 0 GAMG: Solving for p, Initial residual = 1.15061e07, Final residual = 1.15061e07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000173498, global = 0.000153117, cumulative = 0.0955565 rho max/min : 2 0.5 smoothSolver: Solving for omega, Initial residual = 2.12539e06, Final residual = 9.40374e11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.39606e05, Final residual = 9.23655e10, No Iterations 1 ExecutionTime = 70803.7 s ClockTime = 95324 s 

November 10, 2014, 00:53 

#42 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
But Later, after a few iterations my simulation crashes showing this
Code:
[7] #0 Foam::error::printStack(Foam::Ostream&) in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [7] #1 Foam::sigFpe::sigHandler(int) in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [7] #2 in "/lib/x86_64linuxgnu/libc.so.6" [7] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" [7] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" [7] #5 [7] in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam" [7] #6 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6" [7] #7 [7] in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam" [EATStandalone:13533] *** Process received signal *** [EATStandalone:13533] Signal: Floating point exception (8) [EATStandalone:13533] Signal code: (6) [EATStandalone:13533] Failing at address: 0x3e8000034dd [EATStandalone:13533] [ 0] /lib/x86_64linuxgnu/libc.so.6(+0x364a0) [0x7fc667b0b4a0] [EATStandalone:13533] [ 1] /lib/x86_64linuxgnu/libc.so.6(gsignal+0x35) [0x7fc667b0b425] [EATStandalone:13533] [ 2] /lib/x86_64linuxgnu/libc.so.6(+0x364a0) [0x7fc667b0b4a0] [EATStandalone:13533] [ 3] /home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x2ab) [0x7fc66d071f0b] [EATStandalone:13533] [ 4] /home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE7correctEv+0x32) [0x7fc66d07f5f2] [EATStandalone:13533] [ 5] rhoPimpleDyMFoam() [0x41f217] [EATStandalone:13533] [ 6] /lib/x86_64linuxgnu/libc.so.6(__libc_start_main+0xed) [0x7fc667af676d] [EATStandalone:13533] [ 7] rhoPimpleDyMFoam() [0x42660d] [EATStandalone:13533] *** End of error message ***  mpirun noticed that process rank 7 with PID 13533 on node EATStandalone exited on signal 8 (Floating point exception).  [7] #0 Foam::error::printStack(Foam::Ostream&) in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [7] #1 Foam::sigFpe::sigHandler(int) in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [7] #2 in "/lib/x86_64linuxgnu/libc.so.6" [7] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" [7] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" [7] #5 [7] in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam" [7] #6 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6" [7] #7 [7] in "/home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/bin/rhoPimpleDyMFoam" [EATStandalone:13837] *** Process received signal *** [EATStandalone:13837] Signal: Floating point exception (8) [EATStandalone:13837] Signal code: (6) [EATStandalone:13837] Failing at address: 0x3e80000360d [EATStandalone:13837] [ 0] /lib/x86_64linuxgnu/libc.so.6(+0x364a0) [0x7fd517db24a0] [EATStandalone:13837] [ 1] /lib/x86_64linuxgnu/libc.so.6(gsignal+0x35) [0x7fd517db2425] [EATStandalone:13837] [ 2] /lib/x86_64linuxgnu/libc.so.6(+0x364a0) [0x7fd517db24a0] [EATStandalone:13837] [ 3] /home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x2ab) [0x7fd51d318f0b] [EATStandalone:13837] [ 4] /home/eatin/OpenFOAM/OpenFOAM2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE7correctEv+0x32) [0x7fd51d3265f2] [EATStandalone:13837] [ 5] rhoPimpleDyMFoam() [0x41f217] [EATStandalone:13837] [ 6] /lib/x86_64linuxgnu/libc.so.6(__libc_start_main+0xed) [0x7fd517d9d76d] [EATStandalone:13837] [ 7] rhoPimpleDyMFoam() [0x42660d] [EATStandalone:13837] *** End of error message ***  mpirun noticed that process rank 7 with PID 13837 on node EATStandalone exited on signal 8 (Floating point exception).  [1]+ Exit 136 mpirun np 8 rhoPimpleDyMFoam parallel > log I am unable to attach my output log file as it is exceeding max file size even after compressing it. I will mail it to you if you give me your mail id's. Thanks 

November 10, 2014, 02:55 

#43 
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 25 
Did you set "default Gauss linear uncorrected"?
I would try to use the "uncorrected" just for the omega term. Or firstly better try "default Gauss linear limited 0.5".
__________________
The skeleton ran out of shampoo in the shower. 

November 10, 2014, 03:00 

#44 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
ok..
I will change it to Gauss Linear Limited 0.5 Please check post #40 

November 10, 2014, 03:03 

#45 
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 25 
Better ask Tobias about that stuff... I won't be any help.
__________________
The skeleton ran out of shampoo in the shower. 

November 10, 2014, 03:58 

#47 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
Please check attachments


November 10, 2014, 05:02 

#48 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
Thanks,
Try and use the following settings: controlDict Code:
deltaT 1e8; writeInterval 2e08; writePrecision 10; maxCo 2; fvSolutions Code:
"(rhoUh)" { solver PBiCG; preconditioner DILU; tolerance 1e06; relTol 0.1; } "(rhoUh)Final" { $U; relTol 0; } "(kepsilonomega)" { solver PBiCG; preconditioner DILU; tolerance 1e10; relTol 0.1; } "(kepsilonomega)Final" { $k; relTol 0; } PIMPLE { momentumPredictor yes; transonic no; nOuterCorrectors 50; nCorrectors 2; nNonOrthogonalCorrectors 1; turbOnFinalIterOnly false; rhoMin rhoMin [ 1 3 0 0 0 ] 0.1; rhoMax rhoMax [ 1 3 0 0 0 ] 2.5; } relaxationFactors { fields { p 0.3; pFinal 1; } equations { "(Uhkepsilonomega)" 0.4; "(Uhkepsilonomega)Final" 1; } } Code:
gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss upwind; //linearUpwind grad(U); div(phi,h) Gauss upwind; //linearUpwind grad(h); div(phi,K) Gauss linear; div(meshPhi,p) Gauss linear; div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div((muEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default none ; laplacian(rAUf,pcorr) Gauss linear corrected; laplacian(muEff,U) Gauss linear corrected; laplacian(alphaEff,h) Gauss linear corrected; laplacian(rhorAUf,p) Gauss linear corrected; laplacian(DepsilonEff,epsilon) Gauss linear uncorrected; laplacian(DkEff,k) Gauss linear uncorrected; laplacian(DomegaEff,omega) Gauss linear uncorrected; } I feel very sad that you did not change your files like we mentioned. All advices we gave to you are useless if you do not take it into account. That makes me very sad and especially the imagination that some people offers their time to reading threads, checking out the cases and errors to give any advice to solve the problem is demotivating because if I spend time to a problem and give some advice and none of the advice is used and the thread starter always ask the same questions, its like »time wasted« ... especially in old topics like » using the PIMPLE alogithm The old posts contained the following changes (that are now included in the three files):
Again this mean no harm to you but in some cases it makes me very very sad to investigate time for such problems if the hints are not used. At least  now  I can imagine why there are too less support of advanced FOAM users due to the above mentioned things. I hope you will »copy paste« or »modify« your files and give feedback.
__________________
Keep foaming, Tobias Holzmann Last edited by Tobi; November 10, 2014 at 06:07. 

November 10, 2014, 05:10 

#49 
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 21 
Just to confirm, your rotor is rotating with 300000 RPM right??


November 10, 2014, 05:15 

#50 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
For me I thought its 300 RPM due to the comma 300,000 > for me 300 RPM, but maybe your are right. If this would be correct (300.000) then there are a few things to change... but never did simulations in that RPM range.
__________________
Keep foaming, Tobias Holzmann 

November 10, 2014, 05:28 

#51 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
Hi tobi,
I really dint mean to make you sad I was implementing the changes to the best of my knowledge. I have read about the PIMPLE ALGORITHM on your blog and tried to implement relaxation factors for variables but ended up getting this Code:
> FOAM FATAL ERROR: previous iteration field IOobject: volScalarField p "/home/eatin/OpenFOAM/OpenFOAM2.3.0/TurboCharger/Trail_run2/0" not stored. Use field.storePrevIter() at start of iteration. From function GeometricField<Type, PatchField, GeoMesh>::prevIter() const in file /home/eatin/OpenFOAM/OpenFOAM2.3.0/src/OpenFOAM/lnInclude/GeometricField.C at line 863. FOAM aborting Code:
field.storePrevIter( ) Philipp has suggested me to use Gauss linear uncorrected for the omega,As i had no knowledge about the other things in laplacian i changed to default Gauss linear uncorrected for the laplacian schemes. I had posted the output in post #41 then he suggested me to use default Gauss linear limited 0.5. please check post #43 I am using OpenFOAM for the first time and getting to know little about it day by day from experts like you through this forum and reading lot of threads. Please dont think i am not following your instructions, i am trying my level best to first understand your replies and implement them. I really thank you guys for helping me out and really appreaciate your efforts to patiently explain concepts to me. 

November 10, 2014, 05:31 

#52 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
@vasava
Yes 300000 rpm 

November 10, 2014, 06:07 

#53 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
Hi,
as I told you its without mean on harm...
__________________
Keep foaming, Tobias Holzmann 

November 10, 2014, 06:24 

#54 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
Hi tobi,
Ya have modified the files as you suggested, changed the timStep to 1e10 and running the simulation. Lot of pimple iterations per time step.Not even completed the first timeStep. Will update you soon with the output. Can you please check the output i got with the previous settings , they are here http://www.cfdonline.com/Forums/ope...dymfoam3.html Post #41 and #42. Is the error because of my thermoPhysicalProperties ?? please confirm Code:
/** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.3.0   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } mixture { specie { nMoles 1; // No. of moles molWeight 28.96; // Molecular weight of air } thermodynamics { Cp 1.0063e+3; // specific heat Hf 2.544e+06; // heat of fusion calculated using h=2501+1.84t where t=25 0C } transport { // mu 1.844e05; // dynamic viscosity // Pr 0.7; // prandtl number As 1.458e06; //sutherland coeffecient Ts 110.4; //sutherland temperature } } // ************************************************************************* // 

November 10, 2014, 06:42 

#55 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
Few clarifications
1. Should i increase my timeStep to 1e09 after few iterations , this was calculated for courant no. 1? or 1e10 is fine 2. Lot of pimple iterations per timeStep , taking high computational time , is this alright or shall i decrease the nOuterCorrectors ? ( Not even 1 timeStep completed yet ) 3. Yes , I am using rhoPimpleDyMFoam and have implemented this change in the code mentioned here http://www.openfoam.org/mantisbt/view.php?id=909 to include the mesh motion. The term was missing before. 

November 10, 2014, 07:08 

#56 
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
Hi,
due to the output I think its your thermodynamic model. You divide through zero (floating point exeption). Problem is that I dont know when this error occur but I think befor calculating the enthalpy h. Then you have problems with the thermodynamic class but I think in your case its more due to some instability problems of your solution. So just use my mentioned modifications and have a look if everything is working. Use pyFoam for that (if you do not know it, it would be nice to check it out). Especially pyFoamPlotWatcher.py...
__________________
Keep foaming, Tobias Holzmann 

November 10, 2014, 07:23 

#57  
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
Quote:
Quote:
But I forgot something. The fvSolution file should look like that: Code:
PIMPLE { momentumPredictor yes; transonic no; nOuterCorrectors 100; nCorrectors 2; nNonOrthogonalCorrectors 1; turbOnFinalIterOnly false; rhoMin rhoMin [ 1 3 0 0 0 ] 0.1; rhoMax rhoMax [ 1 3 0 0 0 ] 2.5; residualControl { "(Ukomega)" { tolerance 1e5; relTol 0; } p { tolerance 1e4; relTol 0; } } } Now you underrelax the equations or matrixes and be much more stable. You use this relaxation factors till you reach your residual control criterions for each variable you set and the last iteration is with no relaxation factors (factors = 1). Normally the first time step will not converge if you have a complex problem but after a few time steps your Pimple loop can decrease from 200 loops to 12 loops. And you will get bigger timesteps and stable simulation. Thats the reason why PIMPLE is better than PISO if you want to have a look at bigger time steps. If you set nOuterCorrections to 1 then you have PISO mode. Therefor maxCo should be set to 0.1  0.7 (depend on your case) and you are not allowed to underrelax anymore. If you use PISO your dT is something around 1e5 to 1e8 (depend on case and problem). For PIMPLE you can get much bigger timesteps. One iteration is more expensive but you go on in time faster. Hope its clear now.
__________________
Keep foaming, Tobias Holzmann Last edited by Tobi; November 11, 2014 at 03:10. 

November 11, 2014, 03:16 

#58 
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 21 
@Tobi: Are rhoMin and rhoMax the value of minimum and maximum density from rhoCoeffs function respectively??


November 11, 2014, 03:30 

#59 
Member
Abhijit
Join Date: Jul 2014
Posts: 75
Rep Power: 10 
Hi ,
So far the simulation is running well, it has reached to 8.5e08 s. I have included the residual control and the pimple is getting converged for 35 iterations each time step. 

November 11, 2014, 04:36 

#60  
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,670
Blog Entries: 6
Rep Power: 48 
Quote:
Now you have to analyse your output. To fasten up the simulation you have several options now:
Example:
Quote:
No. In the pEqn.H the solver calculates your density with your thermodynamic model. After that you bound your density with the values you insert in your fvSolution. That is helpful if you know that the density range of your fluid should be between 0.8 and 1 (or other values). If the thermodynamic calculates 0.001 and 3 in some cells it will be cutted and set to the values you set. That is very helpful if there are only some cells which give unphysical results. It stabilizes your solution BUT you should always check the output cutting is unphysical too. So if you set rhoMin to 0.5 and the solver calculate rhoMin = 0.43, you should decrease your rhoMin because 0.43 seems to be correct.
__________________
Keep foaming, Tobias Holzmann 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Study of the EEqn.H in rhoPimpleDyMFoam.  Horacio Aguerre  OpenFOAM Programming & Development  11  August 19, 2022 02:47 
Developing a rhoPimpleDyMFoam solver  bvieira  OpenFOAM Programming & Development  20  October 9, 2014 12:12 
rhoPimpleDymFoam  jvd.mechanic  OpenFOAM Running, Solving & CFD  0  June 15, 2014 05:20 
Divergence in rhoPimpleDyMFoam  bvieira  OpenFOAM Running, Solving & CFD  1  July 19, 2012 02:22 