UDF的宏用法及相关算例.docx
《UDF的宏用法及相关算例.docx》由会员分享，可在线阅读，更多相关《UDF的宏用法及相关算例.docx（102页珍藏版）》请在冰豆网上搜索。
UDF的宏用法及相关算例
7自定义函数（UDF）
7.1，概述
用户自定义函数（User—DefinedFunctions,即UDFs）可以提高FLUENT程序的标准计算功能。
它是用C语言书写的，有两种执行方式：
interpreted型和compiled型。
Interpreted型比较容易使用，但是可使用代码（C语言的函数等）和运行速度有限制。
Compiled型运行速度快，而且也没有代码使用范围的限制，但使用略为繁琐。
我们可以用UDFs来定义：
a）边界条件
b）源项
c）物性定义（除了比热外）
d）表面和体积反应速率
e）用户自定义标量输运方程
f）离散相模型（例如体积力，拉力，源项等）
g）代数滑流（algebraicslip）混合物模型（滑流速度和微粒尺寸）
h）变量初始化
i）壁面热流量
j）使用用户自定义标量后处理
边界条件UDFs能够产生依赖于时间，位移和流场变量相关的边界条件。
例如，我们可以定义依赖于流动时间的
x方向的速度入口，或定义依赖于位置的温度边界。
边界条件剖面UDFs用宏DEFINE_PROFILE定义。
有关例子可以
在5.1和6.1中找到。
源项UDFs可以定义除了DO辐射模型之外的任意输运方程的源项。
它用宏DEFINE_SOURCE定义。
有关例子在5.2和6.2中可以找到。
物性UDFs可用来定义物质的物理性质，除了比热之外，其它物性参数都可以定义。
例如，我们可以定义依赖于温度的粘性系数。
它用宏DEFINE_PROPERTY定义，相关例子在6.3中。
反应速
率UDFs用来定义表面或体积反应的反应速率，分别用宏DEFINE_SR_RATE和DEFINE_VR_RATE定义，例子见6.4。
离散相模型用宏DEFINE_DPM定义相关参数，见5.4。
UDFs还可以对任意用户自定义标量的输运方程进行初始化，定义壁面热流量，或计算存贮变量值（用用户自定义标量或用户自定义内存量）使之用于后处理。
相关的应用见于5.3,
5.5，5.6和5.7。
UDFs有着广泛的应用，本文并不能一一叙述。
如果在使用中遇到问题，可以联系FLUENT技术支部门要求帮助。
在此推荐一个网站www.cfd-，上面有FLUENT论坛，可进行相关询问和讨论。
7.1.1书写UDFs的基本步骤
在使用UDFs处理FLUENT模型的过程中，我们一般按照下面五步进行：
1.概念上函数设计
2.使用C语言书写
3.编译调试C程序
4.执行UDF
5.分析与比较结果
第一步分析我们所处理的模型，目的是得到我们要书写的UDF的数学表达式。
第二步将数学表达式转化成C语
言源代码。
第三步编译调试C语言源代码。
第四步在FLUENT中执行UDF。
最后一步，将所得到的结果与我们
要求的进行比较，如果不满足要求，则需要重复上面的步骤，直到与我们期望的吻合为止。
7.1.2Interpreted型与Compiled型比较
CompiledUDFs执行的是机器语言，这和FLUENT本身运行的方式是一样
的。
一个叫做Makefile的过程能够激活C编辑器，编译我们的C语言代码，从而建立一个目标代码库，目标代码库中包含有高级C语言的低级机器语言诠释。
在运行的时候，一个叫做“dynamicloading”的过程将目标代码库与FLUENT
连接。
一旦连接之后，连接关系就会在case文件中与目标代码库一起保存，所以读入case文件时，FLUENT就会自动
加载与目标代码库的连接。
这些库的建立是基于特定计算机和特定FLUENT版本的，所以升级FLUENT版本后，就必
须重新建立相应的库。
相反，InterpretedUDFs是在运行的时候直接装载编译C语言代码的。
在这种情况下，生成的机器代码不依赖于计
算机和FLUENT版本。
编译后，函数信息将会保存在case文件中，所以读入case文件时，FLUENT也会自动加载相
应的函数。
InterpretedUDFs具有较强的可移植性，而且编译比较简单。
对于简单的UDFs，如果对运行速度要求不高，
一般就采用Interpreted型的。
下面列出的是两种UDFs的一些特性：
1.InterpretedUDFs
——独立于计算机结构；
——能够完全当作CompiledUDFs使用；——不能与其它编译系统或用户库连接；
――只支持部分C语言，不能包含：
goto语句，非ANSI—C语法，结构，联合，函数指针，函数数组等。
!
InterpretedUDFs能够使用FLUENT提供的宏，间接引用存贮于FLUENT的变量，详见2.10。
2．CompiledUDFs
――运行速度比InterpretedUDFs快；
――能够完全于C语言结合；
――能够用任何兼容ANSI—C的编辑器编译；
――对不同FLUENT版本（2D或3D）需要建立不同的共享库；
――不一定能够当作InterpretedUDFs使用。
我们在使用中首先要根据具体情况，明确使用哪种UDFs，然后才能进一步去实现，在使用中要注意上面叙述的事项。
第二节书写UDFs
7.2.1概述
书写Interpreted型和Compiled型用户自定义函数的过程和书写格式是一样的。
主要的区别在于与C语言的结合程
度，Compiled型能够完全使用C语言的语法，而Interpreted型只能使用其中一小部分，这在前面有过论述。
7.2.2UDF格式
通用的UDF格式由三部分组成：
1.定义恒定常数和包含库文件，分别由#DEFINE和#INCLUDE陈述（见2.3）;
2．宏DEFINE_*定义UDF函数（见2.4）；
3.函数体部分（见2.9）
包含库有udf.h,sg.h,mem.h,prop.h,dpm.h等，其中udf.h是必不可少的，书写格式为#includeUdf.h”；所有数值都应采用SI单位制；函数体部分字母采用小写，Interpreted型只能够包含FLUENT支持的C语言语法和函数。
7.2.3包含库udf.h
库文件udf.h必须C函数开头包含。
7.2.4定义函数
7.2.4.1简介
Fluent公司提供了一套宏，来帮助我们定义函数。
这些宏都以DEFINE_开始，对它们的解释包含在udf.h文件中，所以我们必需要包含库udf.h。
为了方便使用，我们把对udf.h文件中解释宏列在附录A中。
UDF使用宏DEFINE_定义，括号列表中第一个参数代表函数名。
例如
DEFINE_PROFILE（inlet_x_velocity，thread，position）
定义了一个名为inlet_x_velocity的函数。
!
所有函数名必须小写
紧接着函数名的是函数的输入参数，如上函数inlet_x_velocity有两个输入参数：
thread和position,thread是一个
指针，指向数据类型Thread，position是个整数（见2.4.3）。
UDF编译和连接之后，函数名就会出现在FLUENT相应的下拉列表内。
如上述函数，编译连接之后，就能在相应的边界条件面板内找到一个名为inlet_x_velocity的函数，选定之后就可以使用。
7.2.4.2udf.h文件中对宏DEFINE_的解释
在udf.h文件中，对附录A的宏作了解释，例如：
#defineDEFINE_PROFILE（name，t，I）voidname（Thread*t，inti）通用的宏解释格式为
#definemacroreplacement-text
在编译前，C预处理器（即cpp）先进行宏替代。
例如
DEFINE_PROFILE（inlet_x_velocity，thread，position）
替代为
voidinlet_x_velocity（Thread*thread，intposition）替代后的函数返回实型值或不返回任何值。
如上述函数由于是void型的，所以不返回任何值。
7.2.4.3宏DEFINE
宏DEFINE是用来定义UDFs的，可以分为三类：
通用的，离散相的和多相的。
从宏DEFINE下划线的后缀，我们可以看出该宏是用来定义哪种类型函数的。
如DEFINE_SOURCE定义的函数用来修改输运方程源项，DEFINE_PROPERTY定义的函数用来定义物质的物理性质。
通用的宏在2.5详述，离散相和多相的分别在2.6中详述。
下面是附录A的简列。
通用类型：
1.DEFINE_ADJUST
2.DEFINE_DIFFUSIVITY
3.DEFINE_HEAT_FLUX
4.DEFINE_INIT
5.DEFINE_ON_DEMAND
6.DEFINE_PROFILE
7.DEFINE_PROPERTY
8.DEFINE_RW_FILE
9.DEFINE_SCAT_PHASE_FUNC
10.DEFINE_SOURCE
11.DEFINE_SR_RATE
12.DEFINE_UDS_FLUX
13.DEFINE_UDS_UNSTEADY
14.DEFINE_VR_RATE
离散相模型：
1.DEFINE_DPM_BODY_FORCE
2.DEFINE_DPM_DRAG
3.DEFINE_DPM_EROSION
4.DEFINE_DPM_INJECTION_INIT
5.DEFINE_DPM_LAW
6.DEFINE_DPM_OUTPUT
7.DEFINE_DPM_PROPERTY
8.DEFINE_DPM_SCALAR_UPDATE
9.DEFINE_DPM_SOURCE
10.DEFINE_DPM_SWITCH
多相模型：
1.DEFINE_DRIFT_DIAMETER
2.DEFINE_SLIP_VELOCITY
7.2.4.2数据类型的定义
作为对C语言数据类型的补充，FLUENT定义了几种特殊的数据类型，最常用的是：
Thread，cell_t,face_t，Node
和Domain。
Thread是相应边界或网格区域的结构类型数据；cell_t表示单独一个控制体体积元，用来定义源项或物性；face_t
对应于网格面，用来定义入口边界条件等；Node表示相应的网格节点；Domain是一种结构，其中包含所有的threads,
cells，faces禾口nodes。
!
Thread，cell_t，face_t，Node和Domain要区分大小写。
7.2.5通用宏及其定义的函数
宏DEFINE用来定义UDFs，下面是通用宏的具体解释。
7.2.5.1DEFINEADJUST
Name
Arguments
ArgumentsType
ReturnType
DEFINEADJUST
domain
Domain*domain
void
该函数在每一步迭代开始前，即在求解输运方程前执行。
可以用来修改调节流场变量，计算积分或微分等。
参数domain在执行时，传递给处理器，通知处理器该函数作用于整个流场的网格区域。
如何激活该函数请参见4.6,具体求解例子见5.3，5.6和5.7。
7.2.5.2DEFINEDIFFUSIVITY
Name
Arguments
ArgumentsType
ReturnType
DEFINE_DIFFUSIVITY
c，t，i
cell_tc，Thread*t，inti
real
该函数定义的是组分扩散系数或者用户自定义标量输运方程的扩散系数，c代表网格，t是指向网格线的指针，i
表示第几种组分或第几个用户自定义标量（传递给处理器）。
函数返回的是实型数据。
例子见5.3。
725.3DEFINE_HEAT_FLUX
Name
Arguments
ArgumentsType
ReturnType
DEFINE_HEAT_FLUX
f,t,c0,t0,cid,cir
face_tf,Thread*t,cell_tc,Thread*t0,realcid[],realcir[]
void
该函数定义的是网格与邻近壁面之间扩散和辐射热流量。
f表示壁面，t指向壁面线，cO表示邻近壁面的网格，tO
指向网格线。
函数中需要给出热扩散系数（cid）和辐射系数（cir），才能求出扩散热流量（qid）和辐射热流量（qir）。
在计算时，FLUENT按照下面的公式求解：
qid=cid[0]+cid[1]xC_T（cO,tO）—cid[2]xF_T（f,t）—cid[3]xpow（F_T（f,t）,4）
qir=cir[0]+cir[1]xC_T（c0,t0）—cir[2]xF_T（f,t）—cir[3]xpow（F_T（f,t）,4）
该函数无返回值。
如何激活函数见4.7,例子见5.6。
7.2.6DEFINEINIT
Name
Arguments
ArgumentsType
ReturnType
DEFINEINIT
domain
Domain*domain
void
该函数用于初始化流场变量，它在FLUENT默认的初始化之后执行。
作用区域是全场，无返回值。
函数的激活见
4.5,例子见5.4.1和5.5。
7.2.7DEFINE_ON_DEMAND
Name
Arguments
ArgumentsType
ReturnType
DEFINEONDEMAND
void
该函数不是在计算中由FLUENT自动调用，而是根据需要手工调节运行。
如何执行见4.12，例子见5.8。
7.2.8DEFINEPROFILE
Name
Arguments
ArgumentsType
ReturnType
DEFINEPROFILE
t,i
celltc,Thread*t
real
该函数定义边界条件。
t指向定义边界条件的网格线，i用来表示边界的位置。
函数在执行时，需要循环扫遍所有的边界网格线，值存贮在F_PROFILE（f,t,i）中，无返回值。
选择使用本函数见4.1，例子见5.1.1,5.1.2,5.1.3,5.3,6.1.1和6.1.2。
7.2.9DEFINEPROPERTY
Name
Arguments
ArgumentsType
ReturnType
DEFINEPROPERTY
c,t
celltc,Thread*t
real
该函数用来定义物质物性参数。
c表示网格，t表示网格线，返回实型值。
使用见4.3,例子见6.3.1。
7.2.10DEFINERWFILE
Name
Arguments
ArgumentsType
ReturnType
DEFINERWFILE
fp
FILE*fp
void
该函数用于读写case和data文件。
fp是指向所读写文件的指针。
使用见4.11，例子见2.9.8。
7.2.11DEFINESCATPHASEFUNC
Name
Arguments
ArgumentsType
ReturnType
DEFINE_SCAT
PHASEFUNC
c,f
realc,real*f
real
该函数定义DO（DiscreteOrdinate）辐射模型中的散射相函数（radiationscatteringphasefunction）。
计算两个变量：
从i向到j向散射的辐射能量分数和前向散射因子（forwardscatteringfactor）。
c表示的是i和j向夹角的余弦值，散射
的能量分数由函数返回，前向散射因子存贮在指针f所指的变量中。
处理器对每种物质，都会调用此函数，分别建立
各物质的散射矩阵。
7.2.12DEFINESOURCE
Name
Arguments
ArgumentsType
ReturnType
DEFINE_SOURCE
c,t,dS,i
cell_tc,Thread*t,realdS[],inti
real
该函数定义，除了DO辐射模型之外，输运方程的源项。
在计算中，函数需要扫描全场网格。
c表示网格，t表示
网格线，dS表示源项对所求输运方程的标量的偏导数，用于对源项的线性化；i标志所定义源项对应于哪个输运方程。
使用见4.2，例子见5.2.1，522，5.3，6.2.1。
7.2.13DEFINE_SR_RATE
Name
Arguments
ArgumentsType
ReturnType
DEFINE_SR_RATE
f，t，r，mw，yi，rr
face_tf，Thread*t，
Reaction*r，real
*mw，real*yi，real*rr
void
该函数定义表面化学反应速率。
f表示面，t表示面的线，r是结构指针，表示化学反应；mw和yi是个实型指针
数组，mw存贮物质的分子量，yi存贮物质的质量分数，rr设置函数的一个相关参数。
函数无返回值，使用见4.8。
7.2.14DEFINEUDSFLUX
Name
Arguments
ArgumentsType
ReturnType
DEFINE_UDS_FLUX
f，t，i
face_tf，Thread*t，inti
real
该函数定义用户自定义标量输运方程（user-definedscalartransportequations）的对流通量。
f，t分别表示所求通量
的面和面的线，i表示第几个输运方程（有处理器传递给本函数）。
7.2.15DEFINEUDSUNSTEADY
Name
Arguments
ArgumentsType
ReturnType
DEFINE_UDS_UNSTEADY
c，t，i，apu，su
cell_tc，Thread*t，inti，real*apu，real*su
void
该函数定义用户自定义标量输运方程的非稳态项。
c表示网格，t表示网格线，i表示第几个输运方程。
在FLUENT
中，非稳态项移到RHS中，并以下面的方式离散：
unsteady_term=-J=（凶dV
't
［（叭_（卄］
曲—迪V
］加一
二t二t
方程右边第一项为apu,第二项为sU。
本函数无返回值。
7.2.16DEFINEVRRATE
Name
Arguments
ArgumentsType
ReturnType
DEFINE_VR_RATE
c，t，r，mw，yi，rr，rr_t
cell_tc，Thread*t，Reaction*r，real*mw，real*yi，real*rr，real*rrt
void
该函数定义体积化学反应速率。
c表示网格，t表示网格线，r表示结构指针，表示化学反应过程，mw指针数组
指向存贮物质分子量的变量，yi指向物质的质量分数；rr和rr_t分别设置层流和湍流时函数相关参数。
函数无返回值，
使用见4.8，例子见641。
7.2.6离散相模型宏及其定义的函数
离散模型（DPM）的宏定义的函数与通用宏所定义的函数书写格式是一样的。
对于离散相需要强调结构指针p，
可以用它得到颗粒的性质和相关信息。
下面是具体的宏定义。
7.2.6.1DEFINEDPMBODYFORCE
Name
Arguments
ArgumentType
ReturnType
DEFINE_DPM
BODYFORCE
p，i
Tracked_Particle*p，inti
real
该函数用于定义除了重力和拉力之外的所有体积力。
p为结构指针，i可取0,1,3分别表示三个方向的体积力。
函
数返回的是加速度。
使用见4.4,例子见542。
7.2.6.2DEFINE_DPM_DRAG
Name
Arguments
ArgumentType
ReturnType
DEFINE_DPM_DRAG
Re,p
Tracked_Particle*p,realRe
real
该函数定义拉力系数Cd,Re为ReynoIds数，与颗粒直径和相对于液相速度有关。
拉力定义为:
_18」
D2
CdRe
24
Ip
函数返回的值是18*CD*Re/24。
使用见4.4，例子见543。
7.263DEFINEDPMEROSION
Name
Arguments
ArgumentType
ReturnType
DEFINE_DPM_EROSION
p,t,f,normal,alpha,
Vmag,mdot
Tracked_Particle*p,
Thread*t,face_tf,realalpha
realnormal,realVmag,realmdot
void
该函数定义颗粒撞击壁面湮灭或产生速率。
t为撞击面的线，f为撞击面；数组normal存贮撞击面的单位法向量;
alpha中存贮颗粒轨道与撞击面的夹角；Vmag存贮颗粒速度大小，mdot存贮颗粒与壁面撞击率。
函数无返回值，颗粒
湮灭或产生的计算结果存贮在面变量F_STORAGE_R（f,t,SV_DPMS_EROSION）和F_STORAGE_R（f,t,
SV_DPMS_ACCRETION）中。
使用见4.4。
7.264DEFINEDPMINJECTIONINIT
Name
Arguments
ArgumentType
ReturnType
DEFINE_DPMinjection」NIT
i
Injection*I
void
该函数用于定义颗粒注入轨道时的物理性质。
I是指针，指向颗粒产生时的轨道。
对每一次注入，该函数需要在第
一步DPM迭代前调用两次，在随后颗粒进入区域前每一次迭代中再调用一次。
颗粒的初始化，诸如位置，直径和速度可以通过该函数设定。
函数无返回值。
7.265DEFINEDPMLAW
Name
Arguments
ArgumentType
ReturnType
DEFINE_DPM_LAW
p,ci
Tracked_Particle
升级会员 