|
[Sponsors] |
January 2, 2013, 12:45 |
why this printStack erroe has occured?
|
#1 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
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<double>&, Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #5 Foam::pow(Foam::UList<double> 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<double>::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so" [0] #9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [0] #10 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::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,208
Rep Power: 27 |
any idea is appreciated!urgent help!
|
|
January 4, 2013, 06:12 |
|
#3 |
Senior Member
|
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 |
|
January 6, 2013, 08:00 |
|
#4 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
__________________
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,208
Rep Power: 27 |
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<double>&, Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::pow(Foam::UList<double> 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<double>::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so" #9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #10 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::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$ 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<time()&&time()<t1+c1)||(t2<time()&&time()<t2+c1)||(t3<time()&&time()<t3+c1)||(t4<time()&&time()<t4+c1)||(t5<time()&&time()<t5+c1)||(t6<time()&&time()<t6+c1)||(t7<time()&&time()<t7+c1)||(t8<time()&&time()<t8+c1)||(t9<time()&&time()<t9+c1)||(t10<time()&&time()<t10+c1)||(t11<time()&&time()<t11+c1)?1:0"; valueExpression "p0_2*pow(par2,-3.5)" gradientExpression "0"; value uniform 1023382.5; } left { 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_p3/v_r;" "c2=(w_p3+w_w1)/v_r;" "c3=(w_p3+w_w1+w_p1)/v_r;" "c4=(w_p3+w_w1+w_p1+w_w2)/v_r;" "t1=(w_w0-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);" "par3=1-(gamma-1)*magSqr(internalField(U))/(2*gamma*R*T0_3);"//"par2=1+(gamma-1)/2*magSqr(internalField(U))/(gamma*R*internalField(T));"// ); fractionExpression "(t1<time()&&time()<t1+c1)||(t1+c2<time()&&time()<t1+c3)||(t2<time()&&time()<t2+c1)||(t2+c2<time()&&time()<t2+c3)||(t3<time()&&time()<t3+c1)||(t3+c2<time()&&time()<t3+c3)||(t4<time()&&time()<t4+c1)||(t4+c2<time()&&time()<t4+c3)||(t5<time()&&time()<t5+c1)||(t5+c2<time()&&time()<t5+c3)||(t6<time()&&time()<t6+c1)||(t6+c2<time()&&time()<t6+c3)||(t7<time()&&time()<t7+c1)||(t7+c2<time()&&time()<t7+c3)||(t8<time()&&time()<t8+c1)||(t8+c2<time()&&time()<t8+c3)||(t9<time()&&time()<t9+c1)||(t9+c2<time()&&time()<t9+c3)||(t10<time()&&time()<t10+c1)||(t10+c2<time()&&time()<t10+c3)||(t11<time()&&time()<t11+c1)||(t11+c2<time()&&time()<t11+c3)?1:0"; valueExpression "(t1<time()&&time()<t1+c1)||(t2<time()&&time()<t2+c1)||(t3<time()&&time()<t3+c1)||(t4<time()&&time()<t4+c1)||(t5<time()&&time()<t5+c1)||(t6<time()&&time()<t6+c1)||(t7<time()&&time()<t7+c1)||(t8<time()&&time()<t8+c1)||(t9<time()&&time()<t9+c1)||(t10<time()&&time()<t10+c1)||(t11<time()&&time()<t11+c1)? p0_3*pow(par3,3.5) : p0_1*pow(par1,3.5)"; value uniform 932190; gradientExpression "0"; //type fixedValue; //value uniform 250000; //(p0_1/pow(1+(gamma-1)/2*sqr(internalField(U)/(gamma*R*T0_1-(gamma-1)/2*sqr(internalField(U),3.5)) //(w_w1-pos().y)/22<time() && time()<(w_w1+w_p1-pos().y)/22)? //p0_1*pow(1-(gamma-1)*sqr(internalField(U))/(2*gamma*R*T0_1),gamma/(gamma-1)) //2*gamma*R*T0_1 } walls { type zeroGradient; } empty { type empty; } } |
|
January 10, 2013, 05:39 |
|
#6 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37 |
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,208
Rep Power: 27 |
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: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37 |
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,208
Rep Power: 27 |
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: 4,225
Rep Power: 51 |
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,208
Rep Power: 27 |
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: 4,225
Rep Power: 51 |
Quote:
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,208
Rep Power: 27 |
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. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Error occured when I tried to make mesh | sangje2001 | OpenFOAM | 0 | November 12, 2012 14:14 |
Erroe while startng Geometery(design modeller), in workbench | mwtengr | FLUENT | 1 | April 30, 2010 10:51 |
Erroe while compling a UDF | Kappe | Fluent UDF and Scheme Programming | 0 | August 31, 2009 08:29 |
Gambit fatal erroe | mausmi | FLUENT | 6 | July 28, 2005 08:25 |
ERROR #004100018 occured | Ted | CFX | 1 | April 24, 2005 13:44 |