鍔ㄩ噺婧愰」 - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

#include "udf.h" 

#define  Pm 0.03 

#define  Pmx 0.9 

#define  Et 20 

DEFINE_SOURCE(xmom_source,c,t,dS,eqn) 

{ 

  

  

real xc[ND_ND];  

   

   

   

    

  

  

real force; 

  

  

cell_t cell; 

Thread*cell_thread; 

Domain*subdomain; 

real A,B,C,source; 

begin_c_loop_all(cell,cell_thread) 

{ 

  

   

   

   

 C_CENTROID(xc,cell,cell_thread); 

  

   

   

   

   

   

 if(xc[0]>=-25 && xc[0]<=-12 ) 

  

   

   

   

   

   { 

 

  

   

   

   

   

   

 if(xc[1]>=7 && xc[1]<=12) 

  

  

  

  

force=3200*pow(10,-275*fabs(xc[0])-275*fabs(xc[1])); 

  

  

  

  

A=1/(C_R(c,t)*force); 

  

   

   

   

   

   

 B=2*Pm*Pmx*Et/C_R(c,t); 

  

   

   

   

   

   

 source=C_R(c,t)*(sqrt(A*C_U(c,t)*C_U(c,t)+B)-C_U(c,t)); 

  

   

   

   

   

   

 dS[eqn]=A*C*pow((A*C_U(c,t)*C_U(c,t)+B),3/2)-C; 

  

  

  

  

  

  

  

  

return source; 

  

   

   

   

   

   

   

  

} 

end_c_loop_all(cell,cell_thread) 

  

   } 

}