CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

why this printStack erroe has occured?

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

Like Tree1Likes
  • 1 Post By adambarfi

Reply
 
LinkBack Thread Tools Display Modes
Old   January 2, 2013, 12:45
Default why this printStack erroe has occured?
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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 is offline   Reply With Quote

Old   January 3, 2013, 11:10
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
any idea is appreciated!urgent help!
immortality is offline   Reply With Quote

Old   January 4, 2013, 06:12
Default
  #3
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
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.
adambarfi is offline   Reply With Quote

Old   January 6, 2013, 08:00
Default
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by immortality View Post
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
__________________
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
gschaider is offline   Reply With Quote

Old   January 9, 2013, 13:43
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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 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<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;
        
    }
}
immortality is offline   Reply With Quote

Old   January 10, 2013, 05:39
Default
  #6
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
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
ngj is offline   Reply With Quote

Old   January 10, 2013, 08:06
Default
  #7
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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?
immortality is offline   Reply With Quote

Old   January 10, 2013, 08:20
Default
  #8
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
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
ngj is offline   Reply With Quote

Old   January 10, 2013, 09:17
Default
  #9
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
may be this error is from the pow function?does it accept floating point power?
immortality is offline   Reply With Quote

Old   January 10, 2013, 10:13
Default
  #10
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by immortality View Post
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
gschaider is offline   Reply With Quote

Old   January 10, 2013, 10:48
Default
  #11
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
thank bernhard.now whats the way to find and see pow values that are calculated?how to do this?
immortality is offline   Reply With Quote

Old   January 10, 2013, 11:35
Default
  #12
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by immortality View Post
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
gschaider is offline   Reply With Quote

Old   January 10, 2013, 13:31
Default
  #13
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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 09:51
Erroe while compling a UDF Kappe Fluent UDF and Scheme Programming 0 August 31, 2009 07:29
Gambit fatal erroe mausmi FLUENT 6 July 28, 2005 07:25
ERROR #004100018 occured Ted CFX 1 April 24, 2005 12:44


All times are GMT -4. The time now is 22:12.