pavankonchada123
Subscriber

It is the nucleation of vapor to water molecules when the humidity ratio reaches a certain criteria. UDF will check for the humidity condition in each and every cell and if it reaches , small particles with critical radius will be injected into space. So I have selected volume injection with injection_init macro

 loop(p,I->p_init)  /* 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 =log(supersatratio);
criticalradius= (2*0.072)/(998*461.52*C_T(cell,cthread)*(val));
  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);
 Message("critical  %f\n",nucl);
  
 if(4.15e34*val1*saturatedvapordensity(C_T(cell,cthread))*saturatedvapordensity(C_T(cell,cthread)) >10e12)
 {
PP_RHO(p) = 998;  
  PP_FLOW_RATE(p) = nucl*(PP_RHO(p)*M_PI/6*pow((2*criticalradius),3))/0.001;
  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)
  //Message("Volume %f\n",k); 
   PP_T(p) = C_T(cell,cthread);
  
   PP_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3);
   PP_INIT_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3);
 }