UDF程序解释-指定自定义标量梯度函数.docx

002指定自定义标量梯度函数下面这段程序，作用是使用一个用户自定义的标量作为函数，来表示另一个用户定义标量的梯度的函数，有点绕啊。这个UDF每一次迭代调用一次，执行的时候用汇编的方式。看程序：/*UDFfordefininguser-definedscalarsandtheirgradients*/#include"udf.h"DEFINE_ADJUST(adjust_fcn,d)(Thread*t;cell_tc;realK_EL=1.0;*Donothingifgradientisn'tallocatedyet.*/if(!Data_Valid_P()return;thread_loop_c(t,d)if(FLUID_THREAD_P(t)begin_coop_all(c,t)(CJJDSI(C,t,l)+=K_EL*NV_MAG2(C_UDSI_G(c,t,0)*C_VOLUME(c,t)；)end_c_loop_all(c,t)#include"df.h"是一个声明，意思是后面告诉计算机，后面这段程序要调用的“udf.h”里面包含的函数。“udf.h”是一个头文件，是每一个UDF程序必须要包含的，形式是include"df.h",出现在每条程序最开始部分。就是一个库函数，包含了DEFlNE宏的定义，程序提供其他的的宏和函数，比如三角函数，鬲指数函数等等。以汇编的形式编入程序，是写UDF必须有的。如果是新手的写的话，直接第一句写这个，肯定没错。DEFINE_ADJUST(adjuSJfCn,d)是一个DEFINE宏。DEFlNE_AD：JUST是一个通用宏，可以用来调整或修改ANSYSFluent中未作为参数传递的变量。简单说一句，所谓宏，就是一系列程序的合集，打包之后可以方便调用，对于用户来讲可以理解为黑箱，只要提供参数，就能得到结果。你不需要去了解宏内部是怎么运行的，只要知道输入需要什么，输出能给你什么就可以了。当然如果你要了解宏的内部，查手册和资料，看代码都是可以的。Thread*t;Thread,线程，是fluent中的一种结构数据类型，用于存储它所表示的单元格或面组的公共数据。简单可以理解为一种数据类型。整句含义是定义了一个指针，(这里需要一些C语言知识，指针的定义)，指针的名称是t,存储数据类型是Threadocelljc;定义了一个celLt的变量，名字叫C；cell_t是fluent中的一种数据类型，表示在单元格线程中标识特定单元格的整数数据类型。realK_EL=L0;定义了一个实数变量，名字叫K_EL,并赋值为1。定义好了需要的变量，下面开始程序。先说下这个程序的总体思路。先做一个检测，如果梯度还没有指派，那什么都不做。之后开始在域上循环每个线程,之后在每个线程里的单元进行梯度修改。if(!Data_Valid_P()return;Data_VaIiCLP()是用来检测UDF里的变量是否可以的。定义在id.h的头文件中，id.h在Udfh中定义。DataJZaIieLP()的返回值,如果数据传递有效,返回值为L如果无效，返回值为0。！表示非运算，即1取0,0取1。整句的意思是，如果没有传递的数据，则返回。threadOOP_c(t,d)在域上循环单元线程，即在给定域内循环单元线程。域(domain)可以理解成很多数据的集合,其中包括单元(Cell),也可能有面(face),这个宏只在域(domain)里面循环找单元线程(CeIlthread)。括号里的t是前面定义的，d是由程序给定的。begin-cjoop(c1t)这是在单元线程里循环每个单元，之后进行操作，由ed.cjoop(c1t)结束循环。if(FLUID_THREAD_P(t)；if判断语句，if括号后面为真，则执行之后的语句，为假(0)则执行else后面的语句。FLUIDJHREAD_P(t)又是一个宏，用来检测线程里是否是流体单元，如果是返回L不是流体单元返回0。C.DSI(c,t,l)+=K_EL*NV_MAG2(C_UDSI_G(c,t,0)*C_VOLUME(c,t)；+二表示和自己加，写全就是：CJJDSI(C,t,l)=CJJDSI(Ctl)+K_EL*NV_MAG2(C_UDSI_G(c,t.O)*C_VOLUME(c,t)；K_EL是前面定义的实数，初始值为1。C_UDSI(c,t,l)是一个宏，表示访问为用户定义的标量传输方程计算的单元变量，返回值是UDS单元变量。NV_MAG2()是一个宏，求矢量大小的平方值。C-UDSLG(c,t,0)是一个宏，表示访问为用户定义的标量传输方程计算的单元变量，返回值是UDS梯度。C_VOLUME(c,t)是一个个宏，区域可以看到，是t线程里面的c,就是单元线程里的每个单元。JVOLUME(c,t)表示该单元的体积。