 |
 |
 |
|
|
|
[Sponsors] | ||||
divide error despite of being everything in comment |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 16, 2013, 09:02
|
divide error despite of being everything in comment
|
#1 |
|
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26  |
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
|
|
|
|
|
|
February 16, 2013, 15:12
|
|
#2 | |
|
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
Code:
No value defined for p on right therefore using 40{0}
__________________
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 |
||
|
|
|
||
|
February 16, 2013, 17:39
|
|
#3 |
|
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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;
}
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 |
|
|
|
|
|
|
February 17, 2013, 06:00
|
|
#4 | |
|
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
|
|
|
||
|
February 17, 2013, 07:10
|
|
#5 |
|
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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?
|
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| [Commercial meshers] fluent3DMeshToFoam | bego | OpenFOAM Meshing & Mesh Conversion | 31 | August 16, 2023 09:04 |
| [Commercial meshers] converting Fluent mesh to openfoam standard mesh | deepesh | OpenFOAM Meshing & Mesh Conversion | 31 | March 29, 2017 05:59 |
| dsmcInitialise - dsmcFoam | archymedes | OpenFOAM Pre-Processing | 94 | July 15, 2016 16:14 |
| 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 |
Linear Mode
