一个简单UDF程序的解释-汇总湍流耗散并在每次迭代之后在控制台显示.docx

《一个简单UDF程序的解释-汇总湍流耗散并在每次迭代之后在控制台显示.docx》由会员分享，可在线阅读，更多相关《一个简单UDF程序的解释-汇总湍流耗散并在每次迭代之后在控制台显示.docx（4页珍藏版）》请在优知文库上搜索。

1、004汇总湍流耗散并在每次迭代之后在控制台显示DEFINE_EXECUTE_AT_END是一个通用宏，其含义是，稳态计算时，在每一次迭代完成之后执行，瞬态分析时，在每个时间载荷步执行。当想计算特定时间的流量时，可以使用该宏。程序会根据你的设置，稳态或瞬态，自动判定是在程序执行的时机。看程序：/*UDFforintegratingturbulentdissipationanddisplayingitintheconsoleattheendofthecurrentiterationortimestepW*/#includeudf.hDEFINE_EXECUTE_AT_END(execute_at_

2、end)(Domain*d;Thread*t;*Integratedissipation.*/realsum_diss=0;cell_tc;d=Get_Domain(l);*mixturedomainifmultiphase*/thread_loop_c(t,d)if(FLUID_THREAD_P(t)(begin_c_loOP(C,t)sm-diss+=C_D(c,t)*C_V0LUME(c,t)；end_c_loop(c,t)printf(,Volmeintegralofturbulentdissipation:gn,sm-diss);fflsh(stdot);)#includedf.h是

3、一个声明，意思是后面告诉计算机，后面这段程序要调用的“udf.h”里面包含的函数。Idf卜”是一个头文件是每一个UDF程序必须要包含的,调用形式是的ncludedf.h,出现在每条程序最开始部分。是一个库函数，包含了DEFINE宏的定义，程序提供其他的的宏和函数，比如三角函数，鬲指数函数等等。以汇编的形式编入程序，是写UDF必须有的。如果是新手的写的话，直接第一句写这个，肯定没错。DEFINE_EXECUTE_AT_END是一个宏，前面已经解释过了，execute_at_end是这个UDF的名字，用户自己定义，主要是用来标识。符合规范的取名方式都是可以的。Domain*d;定义了域类型的指针，

4、指针名叫d。域的定义，是一群节点，面，单元线程的集合。Thread*t;Thread,线程，是fluent中的一种结构数据类型，用于存储它所表示的单元格或面组的公共数据。简单可以理解为一种数据类型。整句含义是定义了一个指针，(这里需要一些C语言知识，指针的定义)，指针的名称是t,存储数据类型是Threadorealsum_diss=O;定义了一个实数变量，名字叫SUm_diss,并赋值0。这个就是最终湍流耗散的汇总值，都存在这个变量里。celljc;定义了一个CelLt的变量，名字叫c；CelLt是fluent中的一种数据类型，表示在单元格线程中标识特定单元格的整数数据类型。定义好了数据，开始

5、执行程序。d二GeLDomain(I);Get_Domain()是一个宏，其含义是，在域指针没有作为参数传递给UDF时，用GejDomain来检索。在单向流时，里面的参数选1,返回流体域的指针；多相流时，1代表混合相，其他的数字代表其他的每项，具体值可以查看每项的ID值。threaLpf(t,d)在域上循环单元线程，即在给定域内循环单元线程。域(domain)可以理解成很多数据的集合,其中包括单元(CelI),也可能有面(face),这个宏只在域(domain)里面循环找单元线程(CeIlthread)。括号里的t是前面定义的，d是前面定义，并且用GejDomain宏赋值。if(FLUID_T

6、HREAD_P(t)；if判断语句，if括号后面为真，则执行之后的语句，为假(O)则执行else后面的语句。FLUIDJHREAD_P(t)又是一个宏，用来检测线程里是否是流体单元，如果是返回1,不是流体单元返回0。begin-cjoop(c,t)这是在单元线程里循环每个单元，之后进行操作，由ed.cjoop(c,t)结束循环。c,t由前面定义给出。sum_diss+=C-D(c,t)*C_VOLUME(C,t);+=表示和自己加，sm-diss是前面定义的实数。写全就是：sr-diss=SUrrLdiss+C-D(c,t)*C_VOLUME(c,t)；C_D(Ct)和CLVOLUME(C,t)又是两个宏，区域可以看到，是t线程里面的单元C1就是单元线程里的每个单元。CLD(C,t)表示该单元的湍流动能耗散率；CLVoLUME(C,t)表示该单元的体积。printf(,Volmeintegralofturbulentdissipation:%gn,sum_diss);打印并显示sr-diss的值。fflsh(stdot);清空标准输出的缓冲区。因为每次迭代或者载荷步更新较快,可能没有显示完，后面的数字就来了，所以清空缓冲区，保证显示。