CFD Online Discussion Forums - why this printStack erroe has occured?

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - why this printStack erroe has occured? (https://www.cfd-online.com/Forums/openfoam-solving/111285-why-printstack-erroe-has-occured.html)

why this printStack erroe has occured?

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$

any idea is appreciated!urgent help!

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

Quote:

Originally Posted by immortality (Post 399786)

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<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$

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

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$

for example p is as below:(don'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<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;

        

    }

}

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

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?

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

may be this error is from the pow function?does it accept floating point power?

Quote:

Originally Posted by immortality (Post 401185)

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.

thank bernhard.now whats the way to find and see pow values that are calculated?how to do this?

Quote:

Originally Posted by immortality (Post 401215)

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

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.