CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Problem with udf radiation model in 2d combustion model..????

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 28, 2009, 10:44
Question Problem with udf radiation model in 2d combustion model..????
  #1
New Member
 
siva kumar
Join Date: Jul 2009
Posts: 13
Rep Power: 17
shshiva17 is on a distinguished road
Hi,,,This is siva kumar from IIT MADRAS....while am interpret my udf the error is showing like.....
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
this error is coming in my fluent model when i interpret the source code as shown in below....fuel is methane
#include "udf.h"
DEFINE_SOURCE(energy,c,t,ds,eqn)
{
/*FILE *fp;*/
real source,tep,sigma,Pamb,tamb;
real MW[18];
real Xh2o,Xco,Xco2,Xch4,ph2o,pco,pco2,pch4;
real MWmix,sum=0;
real kapa,absh2o,absco,absco2,absch4;
real a0,a1,a2,a3,a4,a5;
real b0,b1,b2,b3,b4,b5;
real c0,c1,c2,c3,c4;
real d0,d1,d2,d3,d4;
int i;
/*fp=fopen("properties.txt","w");*/
tep=C_T(c,t);
MW[0]=1;
MW[1]=32;
MW[2]=17;
MW[3]=16;
MW[4]=2;
MW[5]=18;
MW[6]=33;
MW[7]=28;
MW[8]=44;
MW[9]=16;
MW[10]=15;
MW[11]=14;
MW[12]=30;
MW[13]=29;
MW[14]=31;
MW[15]=34;
MW[16]=13;
MW[17]=28;
sigma=0.000000056697;
Pamb=1.0;
tamb=298.0;
a0=18.741;
a1=-121.310;
a2=273.500;
a3=-194.05;
a4=56.31;
a5=-5.8169;
b0=-0.23093;
b1=1.12393;
b2=9.4153;
b3=-2.99880;
b4=0.51382;
b5=-.000018684;
c0=6.6334;
c1=-0.0035686;
c2=0.000000016682;
c3=0.00000000025611;
c4=-0.000000000000026558;
if(tep<=750)
{
d0=4.7869;
d1=-0.06953;
d2=0.000295775;
d3=-0.000000425732;
d4=0.000000000202894;
}
else
if(tep>750)
{
d0=10.09;
d1=-0.01183;
d2=0.0000047753;
d3=-0.000000000587209;
d4=-0.000000000000025334;
}
for(i=0;i<18;i++)
{
sum=sum+C_YI(c,t,i)/MW[i];
}
MWmix=1/sum;
Xh2o=C_YI(c,t,5)*MWmix/18;
Xco=C_YI(c,t,7)*MWmix/28;
Xco2=C_YI(c,t,8)*MWmix/44;
Xch4=C_YI(c,t,9)*MWmix/16;
ph2o=Xh2o*Pamb;
pco=Xco*Pamb;
pco2=Xco2*Pamb;
pch4=Xch4*Pamb;
absch4=c0+(c1*tep)+(c2*pow(tep,2))+(c3*pow(tep,3)) +(c4*pow(tep,4)
);
absco2=a0+(a1*(1000/tep))+(a2*pow(1000/tep,2))+(a3*pow(1000/tep,3
))+(a4*pow(1000/tep,4))+(a5*pow(1000/tep,5));
absco=d0+(d1*tep)+(d2*pow(tep,2))+(d3*pow(tep,3))+ (d4*pow(tep,4))
;
absh2o=b0+(b1*(1000/tep))+(b2*pow(1000/tep,2))+(b3*pow(1000/tep,3
))+(b4*pow(1000/tep,4))+(b5*pow(1000/tep,5));
kapa=(absch4*pch4)+(absco*pco)+(absco2*pco2)+(absh 2o*ph2o);
source=-(4*sigma*kapa*(pow(tep,4)-pow(tamb,4)));
/*fprintf(fp,"%f%f%f\n",tep,source,MWmix);*/
ds[eqn]=0;
return source;
}

why the error is coming i don't know please help me??????
shshiva17 is offline   Reply With Quote

Reply


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Modified k-e turbulence model UDF Travis Fluent UDF and Scheme Programming 7 November 11, 2018 21:21
DO radiation method, non-gray band model, UDF wouter Fluent UDF and Scheme Programming 5 March 23, 2017 23:03
Combustion model for CAI/HCCI in AVL FIRE stoodent Main CFD Forum 0 July 19, 2009 19:37
problem about EDC combustion model Ma Shihu FLUENT 1 January 19, 2005 04:25
Fluent incident radiation problem Michael Schwarz Main CFD Forum 0 October 21, 1999 06:56


All times are GMT -4. The time now is 02:59.