# why this printStack erroe has occured?

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

 January 2, 2013, 12:45 why this printStack erroe has occured? #1 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 I use grovyBC.is it possible this error has been initiated from that?how? Code: ```diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 2.74848e-07, Final residual = 5.02194e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 3.91135e-07, Final residual = 3.67893e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 3.33139e-06, Final residual = 8.32991e-15, No Iterations 3 ExecutionTime = 184.8 s ClockTime = 283 s Mean and max Courant Numbers = 0.00676626 0.101999 deltaT = 7.48578e-09 Time = 7.94685e-05 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 2.70326e-07, Final residual = 5.06875e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 4.03258e-07, Final residual = 3.73739e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 3.67996e-06, Final residual = 8.35622e-15, No Iterations 3 ExecutionTime = 184.83 s ClockTime = 283 s Mean and max Courant Numbers = 0.00658396 0.139851 deltaT = 5.3149e-09 Time = 7.94738e-05 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 2.00399e-07, Final residual = 4.79827e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 3.47318e-07, Final residual = 3.89382e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 3.73577e-06, Final residual = 7.01848e-15, No Iterations 3 [0] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #2 in "/lib/x86_64-linux-gnu/libc.so.6" [0] #3 in "/lib/x86_64-linux-gnu/libm.so.6" [0] #4 Foam::pow(Foam::Field&, Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #5 Foam::pow(Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #6 parserPatch::PatchValueExpressionParser::parse() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" [0] #7 Foam::PatchValueExpressionDriver::parse(std::string const&) in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" [0] #8 Foam::groovyBCFvPatchField::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so" [0] #9 Foam::mixedFvPatchField::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [0] #10 Foam::GeometricField::GeometricBoundaryField::evaluate() in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [0] #11 [0] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [0] #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [0] #13 [0] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [thesis-X58A-UD7:13825] *** Process received signal *** [thesis-X58A-UD7:13825] Signal: Floating point exception (8) [thesis-X58A-UD7:13825] Signal code: (-6) [thesis-X58A-UD7:13825] Failing at address: 0x3e800003601 [thesis-X58A-UD7:13825] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36420) [0x7f8095a51420] [thesis-X58A-UD7:13825] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f8095a513a5] [thesis-X58A-UD7:13825] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36420) [0x7f8095a51420] [thesis-X58A-UD7:13825] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x60f6) [0x7f8095fd80f6] [thesis-X58A-UD7:13825] [ 4] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UListIdEERKd+0x41) [0x7f8096aeb5e1] [thesis-X58A-UD7:13825] [ 5] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam3powERKNS_5UListIdEERKd+0x62) [0x7f8096aedd02] [thesis-X58A-UD7:13825] [ 6] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so(_ZN11parserPatch26PatchValueExpressionParser5parseEv+0x6ff0) [0x7f8087e08a30] [thesis-X58A-UD7:13825] [ 7] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so(_ZN4Foam26PatchValueExpressionDriver5parseERKSs+0x6b) [0x7f8087e999cb] [thesis-X58A-UD7:13825] [ 8] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so(_ZN4Foam20groovyBCFvPatchFieldIdE12updateCoeffsEv+0x55) [0x7f808f74cbb5] [thesis-X58A-UD7:13825] [ 9] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam17mixedFvPatchFieldIdE8evaluateENS_8UPstream10commsTypesE+0x226) [0x7f8098844976] [thesis-X58A-UD7:13825] [10] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam(_ZN4Foam14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEE22GeometricBoundaryField8evaluateEv+0x1b2) [0x43caf2] [thesis-X58A-UD7:13825] [11] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam() [0x423f45] [thesis-X58A-UD7:13825] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f8095a3c30d] [thesis-X58A-UD7:13825] [13] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam() [0x4287ed] [thesis-X58A-UD7:13825] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 0 with PID 13825 on node thesis-X58A-UD7 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Killing PID 13821 PyFoam WARNING on line 232 of file /usr/local/lib/python2.7/dist-packages/PyFoam/Execution/FoamThread.py : Process 13821 was already dead Getting LinuxMem: [Errno 2] No such file or directory: '/proc/13821/status' thesis@thesis-X58A-UD7:~/Desktop/laminar_gap_.42\$```

 January 3, 2013, 11:10 #2 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 any idea is appreciated!urgent help!

 January 4, 2013, 06:12 #3 Senior Member     Mostafa Mahmoudi Join Date: Jan 2012 Posts: 319 Rep Power: 8 please give more information about your solver, your geometry, if you change the solver, your BC and etc. nevertheless, see the below links: http://openfoamwiki.net/index.php/HowTo_debugging http://openfoamwiki.net/index.php/Ma...:Ostream.26.29 http://www.cfd-online.com/Forums/ope...rintstack.html immortality likes this.

January 6, 2013, 08:00
#4
Assistant Moderator

Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,997
Rep Power: 42
Quote:
 Originally Posted by immortality I use grovyBC.is it possible this error has been initiated from that?how? Code: ```smoothSolver: Solving for Uy, Initial residual = 3.47318e-07, Final residual = 3.89382e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 3.73577e-06, Final residual = 7.01848e-15, No Iterations 3 [0] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #2 in "/lib/x86_64-linux-gnu/libc.so.6" [0] #3 in "/lib/x86_64-linux-gnu/libm.so.6" [0] #4 Foam::pow(Foam::Field&, Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #5 Foam::pow(Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #6 parserPatch::PatchValueExpressionParser::parse() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" [0] #7 Foam::PatchValueExpressionDriver::parse(std::string const&) in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" [0] #8 Foam::groovyBCFvPatchField::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so" [0] #9 Foam::mixedFvPatchField::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [0] #10 Foam::GeometricField::GeometricBoundaryField::evaluate() in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [0] #11 [0] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [0] #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [0] #13 [0] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" [thesis-X58A-UD7:13825] *** Process received signal *** [thesis-X58A-UD7:13825] Signal: Floating point exception (8) [thesis-X58A-UD7:13825] Signal code: (-6) [thesis-X58A-UD7:13825] Failing at address: 0x3e800003601 [thesis-X58A-UD7:13825] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36420) [0x7f8095a51420] [thesis-X58A-UD7:13825] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f8095a513a5] [thesis-X58A-UD7:13825] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36420) [0x7f8095a51420] [thesis-X58A-UD7:13825] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x60f6) [0x7f8095fd80f6] [thesis-X58A-UD7:13825] [ 4] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UListIdEERKd+0x41) [0x7f8096aeb5e1] [thesis-X58A-UD7:13825] [ 5] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam3powERKNS_5UListIdEERKd+0x62) [0x7f8096aedd02] [thesis-X58A-UD7:13825] [ 6] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so(_ZN11parserPatch26PatchValueExpressionParser5parseEv+0x6ff0) [0x7f8087e08a30] [thesis-X58A-UD7:13825] [ 7] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so(_ZN4Foam26PatchValueExpressionDriver5parseERKSs+0x6b) [0x7f8087e999cb] [thesis-X58A-UD7:13825] [ 8] /home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so(_ZN4Foam20groovyBCFvPatchFieldIdE12updateCoeffsEv+0x55) [0x7f808f74cbb5] [thesis-X58A-UD7:13825] [ 9] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam17mixedFvPatchFieldIdE8evaluateENS_8UPstream10commsTypesE+0x226) [0x7f8098844976] [thesis-X58A-UD7:13825] [10] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam(_ZN4Foam14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEE22GeometricBoundaryField8evaluateEv+0x1b2) [0x43caf2] [thesis-X58A-UD7:13825] [11] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam() [0x423f45] [thesis-X58A-UD7:13825] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f8095a3c30d] [thesis-X58A-UD7:13825] [13] /opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam() [0x4287ed] [thesis-X58A-UD7:13825] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 0 with PID 13825 on node thesis-X58A-UD7 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Killing PID 13821 PyFoam WARNING on line 232 of file /usr/local/lib/python2.7/dist-packages/PyFoam/Execution/FoamThread.py : Process 13821 was already dead Getting LinuxMem: [Errno 2] No such file or directory: '/proc/13821/status' thesis@thesis-X58A-UD7:~/Desktop/laminar_gap_.42\$```
Yep (the stack-trace leads to groovy). groovyBC doesn't check if your expression can cause a FPE. What caused the FPE depends on your expression (just some examples: division by zero, sqrt of a negative, overflow etc). Have a look at your expressions and use your mathematical common sense
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

 January 9, 2013, 13:43 #5 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 hi bernhard I have checked the groovyBC and everything seems true but this error is persisting since two weeks by rhoCentralFoam. Code: ```diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 7.57613e-08, Final residual = 2.12151e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 2.04106e-07, Final residual = 2.02439e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 5.56355e-08, Final residual = 6.94597e-13, No Iterations 3 ExecutionTime = 505.77 s ClockTime = 507 s Mean and max Courant Numbers = 0.000956127 0.00999806 deltaT = 3.76016e-10 Time = 2.55856e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 7.58979e-08, Final residual = 1.87862e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 2.041e-07, Final residual = 2.13119e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 5.53532e-08, Final residual = 6.95182e-13, No Iterations 3 ExecutionTime = 505.88 s ClockTime = 507 s Mean and max Courant Numbers = 0.000956127 0.00999906 deltaT = 3.76016e-10 Time = 2.55893e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 7.60312e-08, Final residual = 1.8685e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 2.04094e-07, Final residual = 1.89043e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 5.50738e-08, Final residual = 6.94215e-13, No Iterations 3 ExecutionTime = 505.99 s ClockTime = 507 s Mean and max Courant Numbers = 0.000956127 0.0100001 deltaT = 3.75696e-10 Time = 2.55931e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 7.60964e-08, Final residual = 1.76105e-17, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 2.03913e-07, Final residual = 2.02328e-17, No Iterations 3 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 5.4751e-08, Final residual = 6.93696e-13, No Iterations 3 #0 Foam::error::printStack(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 in "/lib/x86_64-linux-gnu/libm.so.6" #4 Foam::pow(Foam::Field&, Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::pow(Foam::UList const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #6 parserPatch::PatchValueExpressionParser::parse() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" #7 Foam::PatchValueExpressionDriver::parse(std::string const&) in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so" #8 Foam::groovyBCFvPatchField::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so" #9 Foam::mixedFvPatchField::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #10 Foam::GeometricField::GeometricBoundaryField::evaluate() in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" #11 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam" thesis@thesis-X58A-UD7:~/Desktop/WR_1\$``` for example p is as belowdon't afraid of large expressions,they're simple!) Code: ```dimensions [1 -1 -2 0 0 0 0]; internalField uniform 303975; boundaryField { right { type groovyBC; variables ( //"pi=3.1415926535;" "ymax=max(pos().y);" "ymin=min(pos().y);" "r=0.02325;" "rpm=32151;" "omega=rpm*pi/30;" "v_r=r*omega;" "w_cell=.004;" "n=1;" "w_w0=n*w_cell;" "w_w1=(115.63-46)*r*pi/180;" "w_w2=(180-159.28)*r*pi/180;" "w_w3=(21.6-0)*r*pi/180;" "w_w4=(83-61.6)*r*pi/180;" "w_w5=(180-128)*r*pi/180;" "w_p1=(160.28-116.63)*r*pi/180;" "w_p2=(71.1-21.6)*r*pi/180;" "w_p3=(46-0)*r*pi/180;" "w_p4=(129-84)*r*pi/180;" "c1=w_p2/v_r;" "c2=(w_p2+w_w4)/v_r;" "c3=(w_p2+w_w4+w_p4)/v_r;" "c4=(w_p2+w_w4+w_p4+w_w5+w_w3)/v_r;" "t1=(w_w0+w_w3-pos().y)/v_r;" "t2=t1+c4;" "t3=t1+2*c4;" "t4=t1+3*c4;" "t5=t1+4*c4;" "t6=t1+5*c4;" "t7=t1+6*c4;" "t8=t1+7*c4;" "t9=t1+8*c4;" "t10=t1+9*c4;" "t11=t1+10*c4;" "p0_1=303975;" "T0_1=440;" "p0_3=932190;" "T0_3=1248;" "p0_2=1023382.5;" "T0_4=973;" "gamma=1.4;" "R=287.14;" "par1=1-(gamma-1)*magSqr(internalField(U))/(2*gamma*R*T0_1);" "par2=1+(gamma-1)/2*magSqr(internalField(U))/(gamma*R*internalField(T));"//"par2=1-(gamma-1)*magSqr(internalField(U))/(2*gamma*R*T0_2);"// ); fractionExpression "(t1

 January 10, 2013, 05:39 #6 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,740 Rep Power: 29 Hi Ehsan, Are you sure that the expression par2 on the right patch is always positive? You evaluate a power function with a negative power, so if you for some physical/numerical reason have a zero or negative value of par2, then the simulation crashes. This is also what you are told in error message number 4, and further in error message 8 it tells you that it is in your groovyBC boundary condition. (This, however, does not mean that groovyBC has a bug!) Kind regards, Niels

 January 10, 2013, 08:06 #7 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 dear Niels I used only isentropic relations when p0/p=(1+(g-1)/g*M^2)^(g/(g-1)) Then reciprocate it to find p.in par1 also i replaced T from T0/T relation into p0/p. How you found the pow error is related to par2?

 January 10, 2013, 08:20 #8 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,740 Rep Power: 29 Oh, I had not seen that you also use power functions in the other groovyBC boundary condition. So I have to modify the argument: the error arises from any of the power-functions you are using. Why one or more of these parameters turn negative is unknown to me. / Niels

 January 10, 2013, 09:17 #9 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 may be this error is from the pow function?does it accept floating point power?

January 10, 2013, 10:13
#10
Assistant Moderator

Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,997
Rep Power: 42
Quote:
 Originally Posted by immortality may be this error is from the pow function?does it accept floating point power?
Yep. But for instance a fractional power is something like a square root so negative values might blow up. And negative powers might get become infinite for small numbers. And positive powers might become infinite for large numbers.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

 January 10, 2013, 10:48 #11 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 thank bernhard.now whats the way to find and see pow values that are calculated?how to do this?

January 10, 2013, 11:35
#12
Assistant Moderator

Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,997
Rep Power: 42
Quote:
 Originally Posted by immortality thank bernhard.now whats the way to find and see pow values that are calculated?how to do this?
Best way ist to look at the expression and think "What values are possible"

The other could be setting in the patch that has the problem the silent variable

debugCommonDriver true;

(this will produce MASSIVE outputs). Define a variable values that holds the values you use the pow-function on and replace your pow-function with

"pow(values,1.5)"

(1.5 is just an example). Now when it fails go back in the output. Somewhere you find the setting of "value" with the values

But the best way is still thinking "what could have gone wrong". Popular are for instance stuff like "pow(T-T0,0.5)" where T went below T0
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request

 January 10, 2013, 13:31 #13 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 19 thank you dear bernhard.you mean that i just add debugCommonDriver true; and take a run? The variables i have defined (par1,par2) are sufficient?because they are as you said at the form of pow(value,1.5) in fact i have pow(par1,3.5 or -3.5) I cant find the problem in expressions because they're just isentropic relations.thanks.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post sangje2001 OpenFOAM 0 November 12, 2012 14:14 mwtengr FLUENT 1 April 30, 2010 09:51 Kappe Fluent UDF and Scheme Programming 0 August 31, 2009 07:29 mausmi FLUENT 6 July 28, 2005 07:25 Ted CFX 1 April 24, 2005 12:44

All times are GMT -4. The time now is 21:10.