pavankonchada123
Subscriber

Definetly I will try with wet steam model, but now again I have stuck at max diameter query. The value is as given in injection-0 tab the minimum diameter is ok and I am always having the maximum no particles with max diameter

loop(p,I->p)  /* Standard ANSYS FLUENT Looping Macro to get particle
                  streams in an Injection */
 {
  cell = PP_CELL(p);
  cthread = PP_CELL_THREAD(p);
satpress = H2O_Saturation_Pressure(C_T(cell,cthread));
parpress = parpressure(cell,cthread);
supersatratio = parpress/satpress;
if((parpress/satpress)>1)
{
val =1.000*(double)log((double)supersatratio);
criticalradius= (2*0.072)/(998*461.52*C_T(cell,cthread)*(val));
//Message("CR %f\n",criticalradius*10e7);
  mix = THREAD_MATERIAL(cthread);
  mixture_species_loop(mix,sp,i)
  {
   mw[i] = MATERIAL_PROP(sp,PROP_mwi);
   yi[i] = C_YI(cell,cthread,i);
  }
  real val1=exp((-0.301*criticalradius*criticalradius)/(1.38e-23*C_T(cell,cthread)));
 
  nucl=4.15e34*0.001*val1*C_VOLUME(cell,cthread)*C_R(cell,cthread)*C_R(cell,cthread);
  if((nucl*(998*M_PI/6*pow((2*criticalradius),3))/0.001) >10e-9)
  {
Message("Val %f\n",(nucl*(998*M_PI/6*pow((2*criticalradius),3))/0.001));
Message("CR %f\n",criticalradius*10e7);
  }
  
if((nucl*(998*M_PI/6*pow((2*criticalradius),3))/0.001) > 10e-12 && criticalradius > 10e-10)
 {
 
PP_RHO(p) = 998;  
  
  PP_INIT_DIAM(p) = 2*criticalradius;
 PP_DIAM(p) = 2*criticalradius;  
 
 PP_VEL(p)[0] = C_U(cell,cthread);
PP_VEL(p)[1] =C_V(cell,cthread);
PP_VEL(p)[2] = C_W(cell,cthread);//C_VOLUME(cell,c_thread)*C_YI(cell,cthread,0)
PP_T(p) = C_T(cell,cthread);
PP_INIT_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3);  
PP_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3); 
PP_FLOW_RATE(p) = nucl*(PP_RHO(p)*M_PI/6*pow((2*criticalradius),3))/0.001;
Message("Entered %f \n",PP_FLOW_RATE(p)*10e35);
real m_air = (C_VOLUME(cell,cthread)*C_YI(cell,cthread,1)*C_R(cell,cthread));
real m_vaprem=(C_VOLUME(cell,cthread)*C_YI(cell,cthread,0)*C_R(cell,cthread))-(nucl*(PP_RHO(p)*M_PI/6*pow((2*criticalradius),3)));
C_YI(cell,cthread,0)=m_vaprem/(m_vaprem+ m_air);
 }
 else
 {
   PP_FLOW_RATE(p) = 0.0;
   PP_INIT_DIAM(p)=0.0;
PP_DIAM(p) = 0;
   PP_INIT_MASS(p)= 0.0;
    PP_MASS(p) = 0;
 
}