Fluent關于UDF實作速度入口随時間變化

Fluent關于UDF實作速度入口随時間變化](這裡寫自定義目錄标題)

fluent編譯過程

小編也是第一次實作UDF在fluent中實作控制速度入口條件随時間變化，主要的學習過程還是網上的一些資料，編譯之前需要設定fluent(我用的是19.1)與VS2013的環境配置，隻有正确的環境配置才能通過編譯。

程式設計語言

程式設計語言用的C，代碼可以自己在TXT檔案中寫好，喜歡用IDE編寫 的話也可去VS中或則其他編譯器。

1. 代碼如下

   #include “udf.h”

   DEFINE_PROFILE(inlet_velocity_O2, thread, position)

   {

   real t,v;

   face_t f;

   begin_f_loop(f, thread)

   {

   t=RP_Get_Real(“flow-time”);

   v =1.736exp(-((t- 29.69)/3.18)((t- 29.69)/3.18)) +31.14exp(-((t-42.86)/22.89)((t-42.86)/22.89));

   F_PROFILE(f,thread,position)=v;

   }

   end_f_loop(f,thread)

   }

   DEFINE_PROFILE(inlet_velocity_CO, thread, position)

   {

   real t,v;

   face_t f;

   begin_f_loop(f, thread)

   {

   t=RP_Get_Real(“flow-time”);

   v=1.034exp(-((t- 29.69)/2.872)((t- 29.69)/2.872)) + 18.59exp(-((t-40.59)/21.4)((t-40.59)/21.4));

   F_PROFILE(f,thread,position)=v;

   }

   end_f_loop(f,thread)

   }

   DEFINE_PROFILE(inlet_velocity_Air, thread, position)

   {

   real t,v;

   face_t f;

   begin_f_loop(f, thread)

   {

   t=RP_Get_Real(“flow-time”);

   v =3.539exp(-((t-29.69)/2.842)((t-29.69)/2.842)) + 63.67exp(-((t-40.37)/21.26)((t-40.37)/21.26));

   F_PROFILE(f,thread,position)=v;

   }

   end_f_loop(f,thread)

   }

   DEFINE_PROFILE(inlet_velocity_CO2, thread, position)

   {

   real t,v;

   face_t f;

   begin_f_loop(f, thread)

   {

   t=RP_Get_Real(“flow-time”);

   v=11.45exp(-((t-33.63)/ 14)((t-33.63)/ 14)) + 10.83exp(-((t-23.19)/30.45)*((t-23.19)/30.45));

   F_PROFILE(f,thread,position)=v;

   }

   end_f_loop(f,thread)

   }

   

這樣就okl!!!