CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   why this printStack erroe has occured? (http://www.cfd-online.com/Forums/openfoam-solving/111285-why-printstack-erroe-has-occured.html)

immortality January 2, 2013 12:45

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$


immortality January 3, 2013 11:10

any idea is appreciated!urgent help!

adambarfi January 4, 2013 06:12

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

gschaider January 6, 2013 08:00

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

immortality January 9, 2013 13:43

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;
       
    }
}


ngj January 10, 2013 05:39

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

immortality January 10, 2013 08:06

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?

ngj January 10, 2013 08:20

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

immortality January 10, 2013 09:17

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

gschaider January 10, 2013 10:13

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.

immortality January 10, 2013 10:48

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

gschaider January 10, 2013 11:35

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

immortality January 10, 2013 13:31

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.


All times are GMT -4. The time now is 11:05.