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

divide error despite of being everything in comment

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 16, 2013, 10:02
Default divide error despite of being everything in comment
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
OF sends a divide error while i have commented all terms except "type".why it continue showing this error message?
Code:
right
    {
       //type zeroGradient;
        type groovyBC;

        /*variables (

                  "r=.0795775;"
                  "rpm=23860.806298;"
                  "omega=rpm*pi/30;"
                  "v_r=r*omega;"
                  "w_cell=.008333;"
                  "n=1;"
                  
                  "w_w0=n*w_cell;"
                  "w_w4=(91.679714-42.208987)*r*pi/180;"
                  "w_w5=(180-137.372516)*r*pi/180;"
                  "w_p2=(137.372516-91.679714)*r*pi/180;"
                  "w_p4=(42.208987-0)*r*pi/180;"
                
                  "c1=w_p4/v_r;"
                  "c2=(w_p4+w_w4)/v_r;"
                  "c3=(w_p4+w_w4+w_p2)/v_r;"
                  "c4=(w_p4+w_w4+w_p2+w_w5)/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_2=2533125;"
                  "T0_2=1131.102789;"
                  "p0_4=1304314.009315;"
                  "T0_4=1233.460862;"
                  "gamma=1.4;"
                  "R=287.14;"
                  "part1=1+(gamma-1)/2*sqr(mag(internalField(U).x))/(gamma*R*internalField(T));"
                  "part2=1-(gamma-1)*sqr(mag(internalField(U).x))/(2*gamma*R*T0_2);"
                  "part4=1-(gamma-1)*sqr(mag(internalField(U).x))/(2*gamma*R*T0_4);"
   

);*/

        //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)||((t1+c2<time()&&time()<t1+c3)||(t2+c2<time()&&time()<t2+c3)||(t3+c2<time()&&time()<t3+c3)||(t4+c2<time()&&time()<t4+c3)||(t5+c2<time()&&time()<t5+c3)||(t6+c2<time()&&time()<t6+c3)||(t7+c2<time()&&time()<t7+c3)||(t8+c2<time()&&time()<t8+c3)||(t9+c2<time()&&time()<t9+c3)||(t10+c2<time()&&time()<t10+c3)||(t11+c2<time()&&time()<t11+c3))&&(internalField(U).x/sqrt(gamma*R*internalField(T))<1)?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_4*pow(part1,-3.5):p0_2*pow(part1,-3.5)";//p0_4*pow(part4,3)*sqrt(part4)//1167464.555675:2386438.388181
        //gradientExpression "0"; 
        //value uniform 1200000;//1304314.009315
        //debug CommonDriver true;
        //debug true;
    }
Code:
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
--> FOAM Warning : 
    From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(const fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
    in file groovyBCFvPatchField.C at line 124
    No value defined for p on right therefore using 40{0}
Reading field U

--> FOAM Warning : 
    From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(const fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
    in file groovyBCFvPatchField.C at line 124
    No value defined for U on right therefore using 40{(0 0 0)}
Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 1e-08

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
swak4Foam: Allocating new repository for sampledGlobalVariables
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 9.235371434e-12, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 9.235534354e-12, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#6  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
Floating point exception
immortality is offline   Reply With Quote

Old   February 16, 2013, 16:12
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,920
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by immortality View Post
OF sends a divide error while i have commented all terms except "type".why it continue showing this error message?

Code:
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
--> FOAM Warning : 
    From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(const fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
    in file groovyBCFvPatchField.C at line 124
    No value defined for p on right therefore using 40{0}
Reading field U

--> FOAM Warning : 
    From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(const fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
    in file groovyBCFvPatchField.C at line 124
    No value defined for U on right therefore using 40{(0 0 0)}
Reading/calculating face flux field phiCreating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 1e-08

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
swak4Foam: Allocating new repository for sampledGlobalVariables
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 9.235371434e-12, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 9.235534354e-12, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#6  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
Floating point exception
There is no groovyBC in the stack-trace. But the problem is your commenting out everything. Have a look at the warning
Code:
No value defined for p on right therefore using 40{0}
and think which density the perfect gas equation gives for a p=0 and what might happen if you divide through that density
__________________
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   February 16, 2013, 18:39
Default
  #3
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
you are right.but i get same error when they're not commented.i have not such an error at beginning the run but there was during the run because flow turning to high velocities.
Code:
right
    {
       //type zeroGradient;
        type groovyBC;
        variables (

                  "r=.0795775;"
                  "rpm=23860.806298;"
                  "omega=rpm*pi/30;"
                  "v_r=r*omega;"
                  "w_cell=.008333;"
                  "n=1;"
                  
                  "w_w0=n*w_cell;"
                  "w_w4=(91.679714-42.208987)*r*pi/180;"
                  "w_w5=(180-137.372516)*r*pi/180;"
                  "w_p2=(137.372516-91.679714)*r*pi/180;"
                  "w_p4=(42.208987-0)*r*pi/180;"
                
                  "c1=w_p4/v_r;"
                  "c2=(w_p4+w_w4)/v_r;"
                  "c3=(w_p4+w_w4+w_p2)/v_r;"
                  "c4=(w_p4+w_w4+w_p2+w_w5)/v_r;"
                  "t1=pos().y/v_r;"//(w_w0-pos().y)
                  "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_2=2533125;"
                  "T0_2=1131.102789;"
                  "p0_4=1304314.009315;"
                  "T0_4=1233.460862;"
                  "gamma=1.4;"
                  "R=287.14;"
                  "part1=1+(gamma-1)/2;"//*sqr(mag(internalField(U).x))/(gamma*R*internalField(T))
                  "part2=1-(gamma-1)*sqr(mag(internalField(U).x))/(2*gamma*R*T0_2);"
                  "part4=1-(gamma-1)*sqr(mag(internalField(U).x))/(2*gamma*R*T0_4);"
   

);

        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)||((t1+c2<time()&&time()<t1+c3)||(t2+c2<time()&&time()<t2+c3)||(t3+c2<time()&&time()<t3+c3)||(t4+c2<time()&&time()<t4+c3)||(t5+c2<time()&&time()<t5+c3)||(t6+c2<time()&&time()<t6+c3)||(t7+c2<time()&&time()<t7+c3)||(t8+c2<time()&&time()<t8+c3)||(t9+c2<time()&&time()<t9+c3)||(t10+c2<time()&&time()<t10+c3)||(t11+c2<time()&&time()<t11+c3))&&(internalField(U).x/sqrt(gamma*R*internalField(T))<1)?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_4*pow(part1,-3.5):p0_2*pow(part1,-3.5)";//p0_4*pow(part4,3)*sqrt(part4)//1167464.555675:2386438.388181
        gradientExpression "0"; 
        value uniform 1304314.009315;//1304314.009315
        //debug CommonDriver true;
        //debug true;
    }
the error:
Code:
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 1e-08

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 9.235667626e-12, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9999999922, Final residual = 5.030244961e-12, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#6  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
Floating point exception
immortality is offline   Reply With Quote

Old   February 17, 2013, 07:00
Default
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,920
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by immortality View Post
you are right.but i get same error when they're not commented.i have not such an error at beginning the run but there was during the run because flow turning to high velocities.

the error:
Code:
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 1e-08

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 9.235667626e-12, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9999999922, Final residual = 5.030244961e-12, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#6  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/sonicFoam"
Floating point exception
Error still happens outside swak/groovy. But the involvement of groovy probably is that somewhere the value gets set to 0. But it is not obvious to me where that might have happened.
__________________
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   February 17, 2013, 08:10
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
i already downgraded it to 0.2.0 version and it starts without any error on dividing although it stops after about 500 iterations.what may be the difference related to?
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
dsmcInitialise - dsmcFoam archymedes OpenFOAM Pre-Processing 93 February 11, 2014 02:22
converting Fluent mesh to openfoam standard mesh deepesh OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 17 October 25, 2013 03:59
fluent3DMeshToFoam bego OpenFOAM 27 May 29, 2013 13:08
Mesh conversion MultiphaseFlowsLab OpenFOAM 13 May 2, 2012 12:52
give your comment for this problem rt Main CFD Forum 2 August 3, 2006 00:03


All times are GMT -4. The time now is 20:50.