UDF总结

UDF使用技巧1、查找相应的函数的时候，可以现在word里面找到相应的函数名字，然后依次去中文帮助文档、英文帮助文档和网页帮助文档，看看详细解释并找找是否有相应的例子。2、打个比方来说，thread就是公路，连接的cell和face,cell和face就相当于公路上汽车停靠的站点，cell_t这个面向的是单元，而face_t面向的是边或者面（二维或三维）在fluent循环过程中，一般是用thread作线程检索，而cell或者face作检索过程中位置（相当于指示位置的参数）参数的指示3、对于UDF来说，积分就是做加法，把通过面上每个网格的质量流量相加4、cell和face的区别，什么时候用cell，什么时候用face5、1.begin,end_c_loopmacroisusedforloopingoverallthecellsinparticularthreadforserialprocessing.Forparallelprocessing,thecellsinsideapartitioncanbecategorizedasinteriorandexteriorcells.Themacrosbegin,end_c_loop_int;begin,end_c_loop_extandbegin,end_c_loop_allareusedforloopingoverinterior,exteriorandallthecells（inapartition）respectively.。Inparallelsimulations,bothbegin,end_c_loopandbegin,end_c_loop_allmacroswilldothesamejob.Forfacestheloopingmacroinparallelarebegin,end_f_loop_int;begin,end_f_loop_extandbegin,end_f_loopforloopingoverinterior,boundaryandallfacesrespectively.Forallpracticalpurpose,theuserneednotseparatetheinteriorandboundaryfacesofapartition.Hence,begin,end_f_loop_intandbegin,end_f_loop_extmacrosarerarelyused.实际问题1、DEFINE_UDS_UNSTEADY中的apu包括的函数是不是不包括当前时刻的变量，而su包含前一时刻的变量，所以用了C_STORAGE_R存储前一时刻的变量。2、等于零是因为计算源项的时候温度还没有更新，所以两个温度值是相等的。这个时候其实是需要两个UDM，分别保存上一步和再上一步的温度。而且可以考虑全部在源项里完成，不用Adjust宏：source=(C_UDMI(c,t,1)-C_UDMI(c,t,0))/TIME_STEP;C_UDMI(c,t,0)=C_UDMI(c,t,1);C_UDMI(c,t,1)=C_T(c,t);returnsource;当然，在刚开始计算的时候要注意下UDM的初始问题，防止出错。从第三个时间步开始，就是正常的了。这个方法挺好的，只要一开始的用DEFINE_INIT,对C_UDMI（c,t,0）=0.；C_UDMI（c,t,1）=0.；一、数据类型~cell_tc单元格标识符face_tf面标识符Thread*t指向线程的指针Thread**pt象限矩阵线程的指针，Node*node指向节点的指针Domain*d指向定义域的指针二、node宏函数NODE_X(node)《节点在x方向上的坐标（real）NODE_Y(node)节点在y方向上的坐标（real）NODE_Z(node)节点在z方向上的坐标（real）F_NNODES(f,t)面上节点的数目（int）三、cell宏函数1、C_CENTROID(x,c,t)realx[ND_ND]单元格的重心坐标x[]注意：这句话直接使用，结果会输出到求解器中，需要使用的时候，调用x[]；其中，x[0]代表x方向坐标，x[1]代表y方向坐标，x[2]代表z方向坐标。2、C_VOLUME(c,t)单元格的体积（real）C_NNODES(c,t)单元格的节点的数目（int）C_NFACES(c,t)单元格的面的数目（int）3、C_FACE(c,t,i)返回给定单元cell_tc和Thread*t的面face_tf；注意：索引号i可以被传递给c_face_loop。C_FACE_THREAD(c,t,i)返回面face_tf（由C_FACE返回）的Thread*t；注意：索引号i可以被传递给c_face_loop。4、C_R(c,t)返回线程t的单元c上的密度C_P(c,t)压力C_U(c,t)速度分量u)C_V(c,t)速度分量vC_W(c,t)速度分量wC_T(c,t)温度C_H(c,t)焓…C_K(c,t)湍流动能kC_NUT(c,t)Spalart-Allmaras模型中的湍流密度C_D(c,t)湍流动能耗散率（和e像，看小木虫）C_O(c,t)比耗散率w,C_YI(c,t,i)组分质量分数5、C_R_G(c,t)返回线程t的单元c上的密度梯度C_P_G(c,t)压力梯度C_U_G(c,t)}速度分量u梯度C_V_G(c,t)速度分量v梯度C_W_G(c,t)速度分量w梯度C_T_G(c,t)温度梯度C_H_G(c,t)?焓梯度C_K_G(c,t)湍流动能k梯度C_NUT_G(c,t)Spalart-Allmaras模型中的湍流密度梯度C_D_G(c,t)湍流动能耗散率梯度（和e像，看小木虫）C_O_G(c,t)"比耗散率w梯度C_YI_G(c,t,i)组分质量分数梯度注意：①C_T_G(c,t)［0］:温度梯度x方向的分量；［1］：y方向；［2］：z方向C_R_G(c,t)只能用于密度基求解器，C_P_G(c,t)只能用于压力基求解器③C_YI_G(c,t,i)只能用于密度基求解器，用于压力基求解器要设置，一切设置好后，输入rpsetvar’species/save-gradients#t④只有当已经求解出包含这个变量的方程时才能得到梯度变量，需要如下设置solve/set/expert回车然后两个yes。输入q可以返回上级目录6、…C_R_RG(c,t)返回线程t的单元c上的密度改造梯度C_P_RG(c,t)压力改造梯度C_U_RG(c,t)速度分量u改造梯度C_V_RG(c,t)速度分量v改造梯度(C_W_RG(c,t)速度分量w改造梯度C_T_RG(c,t)温度改造梯度C_H_RG(c,t)焓改造梯度C_K_RG(c,t)湍流动能k改造梯度&C_NUT_RG(c,t)Spalart-Allmaras模型中的湍流密度改造梯度C_D_RG(c,t)湍流动能耗散率改造梯度(和e像，看小木虫)C_O_RG(c,t)比耗散率w改造梯度C_YI_RG(c,t,i)组分质量分数改造梯度$注意：①C_T_RG(c,t)［0］：温度梯度x方向的分量；［1］：y方向；［2］：z方向C_R_RG(c,t)只用于密度基求解器，C_P_RG(c,t)只用于压力基求解器C_YI_RG(c,t,i)只能用于密度基求解器只有当已经求解出包含这个变量的方程时才能得到梯度变量，需要如下设置solve/set/expert回车然后两个yes。输入q可以返回上级目录7、C_R_M1(c,t)返回上一个时刻(t-At)的密度C_P_M1(c,t)返回上一个时刻(t-At)的压力C_U_M1(c,t)返回上一个时刻(t-At)的速度分量uC_V_M1(c,t)返回上一个时刻(t-At)的速度分量vC_W_M1(c,t)返回上一个时刻(t-At)的速度分量wC_T_M1(c,t)!返回上一个时刻(t-At)的温度C_YI_M1(c,t,i)返回上一个时刻(t-At)的组分质量分数注意：①datafromC_T_M1isavailableonlyifuser-definedscalarsaredefined.②详细可以看宏DEFINE_UDS_UNSTEADY(name,c,t,i,apu,su)。8、C_R_M2(c,t)(返回上两个时刻(t-2At)的密度C_P_M2(c,t)返回上两个时刻(t-2At)的压力C_U_M2(c,t)返回上两个时刻(t-2At)的速度分量uC_V_M2(c,t)返回上两个时刻(t-2At)的速度分量vC_W_M2(c,t)返回上两个时刻(t-2At)的速度分量wC_T_M2(c,t)返回上两个时刻(t-2At)的温度C_YI_M2(c,t,i)返回上两个时刻(t-2At)的组分质量分数注意：①datafromC_T_M2isavailableonlyifuser-definedscalarsaredefined.9、#C_MU_L(c,t)层流粘度C_MU_T(c,t)湍流粘度C_MU_EFF(c,t)有效粘度C_K_L(c,t)导热系数>C_K_T(c,t,prt)湍流热导率C_K_EFF(c,t,prt)有效热导率C_DIFF_L(c,t,i,j)层流组分扩散率C_DIFF_EFF(c,t,i)有效组分扩散率)注意：prt是湍流普朗特常数。10、C_CP(c,t)比热容C_RGAS(c,t)气体常数/相对分子质量C_NUT(c,t)》Spalart-Allmaras湍流粘度11、C_FMEAN(c,t)primarymeanmixturefractionC_FMEAN2(c,t)secondarymeanmixturefractionC_FVAR(c,t)】primarymixturefractionvarianceC_FVAR2(c,t)secondarymixturefractionvarianceC_PREMIXC(c,t)reactionprogressvariableC_LAM_FLAME_SPEED(c,t)层流火焰速度C_SCAT_COEFF(c,t),散射系数C_ABS_COEFF(c,t)吸收系数C_CRITICAL_STRAIN_RATE(c,t)临界应变率C_LIQF(c,t)cell中的液体分数C_POLLUT(c,t,i)(污染物的质量分数注意:①C_LIQFisavailableonlyinfluidcellsandonlyifsolidificationisturnedONC_POLLUT(c,t,i)：0-MassFractionofNO；MassFractionofHCN；MassFractionofNH3；MassFractionofN2O；SootMassFraction；