FLUENT个人学习资料

用户自定义函数UDF编写1.边界条件UDF1・1边界条件与坐标有关v=20-20X(-^-)2x 0.075进口速度为抛物线的源代码#include"udfh" /*C语言头文件边界条件，自动获取DEFINE_PROFILE(inlet_x_velocity,thread,index){realx[ND_ND];realy;face_tf;/*variablestatement,positionvector数组*//*variablestatement,positionvalue*//*variablestatement,facethreadforloop*/begin_f_loop(f,thread)/*thread获取Fluent中用户定义的边界线*/{F_CENTROID(x,f,thread); /*获取面质心,视维度将坐标返回给x*/y=x[1]; /*x[1]表示y方向坐标值*/F_PROFILE(f,thread,index)=20.-y*y/(.0745*.0745)*20.;厂end_f_loop(f,thread)/*thread建立的线索引进行face遍历*/}一1.2边界条件与温度有关T=T=300+100Xsin(兀J)X 0.005边界温度为正弦变化的函数#include"udf.h"/*C语言头文件*/#definePI3.141592654 .十…丁边界条件，自动获取DEFINE_PROFILE(inlet_x_velocity,thread,index)/*variablestatement,positionvector数组/*variablestatement,positionvector数组*//*variablestatement,positionvalue*//*variablestatement,facethreadforloop*/realr[ND_ND];realx;face_tf;begin_f_loop(f,thread)/*thread获取Fluent中用户定义的边界线*/{—F_CENTROID(r,f,thread); /*获取面质心,视维度将坐标返回给x*/x=r[0]; /*r[0]表示x方向坐标值*/F_PROFILE(f,thread,index)=300.+100.*sin(PI*x/0.005);厂end_f_loop(f,thread)/*thread建立的线索引进行face遍历*/}一1.3边界条件与时间有关#include"udf.h" /*C语言头文件*/#definePI3.141592654DEFINE_PROFILE(unsteady,thread,index){——face_tf; /*variablestatement,facethreadforloop*/begin_f_loop(f,thread)/*thread获取Fluent中用户定义的边界线*/{—realt=RP_Get_Real(“flow-time”)；/*获取计算过程中当前时间*/F_PROFILE(f,thread,index)=20.+5.0*sin(10.*t);/*修改边界条件*/厂end_f_loop(f,thread) /*thread建立的线索引进行face遍历*/}—2.内部物性UDF2・1流体粘度与温度有关#include"udf.h"DEFINE_PROPERTY(cell_viscosity,cell,thread)/*修改物性*/{ /*用户不选择边界，函数自动遍历整个计算区域*/realmu_lam; /*变量声明，存储粘度值*/realtemp=C_T(cell,thread); /*对计算区域cell遍历，取温度*/if(temp>288.)mu_lam=5.5e-3;elseif(temp>286.)mu_lam=143.2135-0.49725*temp;elsemu_lam=1.;returnmu_lam; /*返回给物性参数，流体粘度选项*/} — /*物性参数修改为返回型函数*/3・源项编写UDF3.1多孔介质区域的x维度源项方程式source=-0.5xC2py\v\v（其中C2是常数，取100）2 XX 2源项source中p可时时获取，对定区域内是常量。同理，质心y是常量。此例中，源项的变量为耳，求解出dS/dq=dS[eqn]。XX•••source=-A\v\v=dS/du=-2^\v\XXXX源项UDF必须给出源项对所求变量的直接导数，源项的导数用于对源项的线性化。下面是根据上述方程编写的UDF的C源程序。

include"udf.h"defineC2100.0DEFINE_SOURCE(xmom_source,cell,thread,dS,eqn){ /*全局cell遍历，x方向源项添加*/realx[ND_ND]; /*存储质心坐标，源项中y值*/realcon,source;C_CENTROID(x,cell,thread); /*获取质心坐标*/con=C2*0.5*C_R(cell,thread)*x[1]; /*得源项中的密度或A值*/source=-con*fabs(C_U(cell,thread))*C_U(cell,thread);/*源项式*/dS[eqn]=-2.*con*fabs(C_U(cell,thread)); /*源项倒数*/returnsource; /*返回源项*/} /*CU、CV、CW三维度下速度*/3.2管内流动源项的编写温度为290K液态金属进入管道，管道上壁280K对液态金属冷却。当液态金属温度低于288K时，就添加动量方程的源项。由于添加了源项，液态金属速度会逐渐降为零。在本例中，仅以液态金属的速度大小来表征是否凝固，如果速度等于零，则认为液态金属已经凝固，而不是求解能量方程。x方向动量源项方程为：(式中，C=20)S=S=—Cv=込=—C=dS[eqn]dvx#include"udf.h"#defineCON20.0DEFINE_SOURCE(cell_x_source,cell,thread,dS,eqn){"""realsource;if(C_T(cell