Fluent UDF教程
UDF中文教程
目录
第一章. 介绍 (4)
1.1什么是UDF？ (4)
1.2为什么要使用UDF？ (4)
1.3 UDF的局限 (5)
1.4Fluent5到Fluent6UDF的变化 (5)
1.5 UDF基础 (6)
1.6 解释和编译UDF的比较 (8)
1.7一个step-by-stepUDF例子 (9)
第二章.UDF的C语言基础 (16)
2.1引言 (16)
2.2注释你的C代码 (17)
2.3FLUENT的C数据类型 (17)
2.4常数 (17)
2.5变量 (17)
2.6自定义数据类型 (20)
2.7强制转换 (20)
2.8函数 (20)
2.9 数组 (20)
2.10指针 (21)
2.11 控制语句 (22)
2.12常用的C运算符 (24)
2.13 C库函数 (24)
2.14 用#define实现宏置换 (26)
2.15 用#include实现文件包含 (27)
2.16 与FORTRAN 的比较 (27)
UDF 第3章写UDF (27)
3.1概述（Introduction） (28)
3.2写解释式UDF的限制 (28)
3.3 FLUENT求解过程中UDF的先后顺序 (29)
3.4 FLUENT 网格拓扑 (31)
3.5 FLUENT数据类型 (32)
3.6 使用DEFINE Macros定义你的UDF (33)
3.7在你的UDF源文件中包含udf.h文件 (34)
3.8在你的函数中定义变量 (34)
3.9函数体（Functin Body） (35)
3.10 UDF任务（UDF Tasks） (35)
3.11为多相流应用写UDF (41)
3.12在并行下使用你的UDF (50)
第四章DEFINE宏 (51)
4.1 概述 (51)
4.2 通用解算器DEFINE宏 (52)
4.3 模型指定DEFINE宏 (61)
4.4 Multiphase DEFINE Macros (105)
4.5 DPM DEFINE Macros (112)
第五章使用宏存取FLUENT解算器变量 (140)
5.1 介绍 (141)
5.2 单元格宏 (142)
5.3表面宏 (145)
5.4 几何宏 (146)
5.5 节点宏 (149)
5.6多相宏 (149)
5.7 DPM宏 (150)
5.8 Nox宏 (150)
6. Utilities （工具） (151)
6.1简要 (152)
6.2一般目的的循环宏 (152)
6.3多相组分查询宏 (156)
6.4设置面变量 (160)
6.5访问没有赋值的自变量 (161)
6.6访问邻近网格和线的变量 (169)
6.7 用户为网格定义内存 (170)
6.8 矢量工具 (171)
6.9 与非定常数值模拟有关的宏 (174)
6.10其它各种工具宏 (175)
第七章UDF的编译与链接 (178)
7.1 介绍 (178)
7.2 UDF解释 (180)
7.3 编译UDF (184)
第八章在FLUENT中激活你的UDF (197)
8.1 激活通用求解器UDF (197)
8.2 激活模型明确UDF (199)
8.3激活多相UDF (206)
8.4 激活DPM UDF (209)
第九章自定义标量及用法 (214)
9.1 介绍 (214)
9.2 理论 (214)
9.3 UDS的定义，求解，后处理 (216)
第一章.介绍
本章简要地介绍了用户自定义函数(UDF)及其在Fluent中的用法。在1.1到1.6节中我们会介绍一下什么是UDF；如何使用UDF，以及为什么要使用UDF，在1.7中将一步步的演示一个UDF例子。
1.1 什么是UDF？
1.2 为什么要使用UDF？
1.3 UDF的局限
1.4 Fluent5到Fluent6 UDF的变化
1.5 UDF基础
1.6 解释和编译UDF的比较
1.7一个step-by-stepUDF例子
1.1什么是UDF？
用户自定义函数，或UDF，是用户自编的程序，它可以动态的连接到Fluent求解器上来提高求解器性能。用户自定义函数用C语言编写。使用DEFINE宏来定义。UDF中可使用标准C语言的库函数，也可使用Fluent Inc.提供的预定义宏，通过这些预定义宏，可以获得Fluent求解器得到的数据。
UDF使用时可以被当作解释函数或编译函数。解释函数在运行时读入并解释。而编译UDF则在编译时被嵌入共享库中并与Fluent连接。解释UDF用起来简单，但是有源代码和速度方面的限制不足。编译UDF执行起来较快，也没有源代码限制，但设置和使用较为麻烦。
1.2为什么要使用UDF？
一般说来，任何一种软件都不可能满足每一个人的要求，FLUENT也一样，其标准界面及功能并不能满足每个用户的需要。UDF正是为解决这种问题而来，使用它我们可以编写FLUENT代码来满足不同用户的特殊需要。当然，FLUENT的UDF并不是什么问题都可以解决的，在下面的章节中我们就会具体介绍一下FLUENT UDF的具体功能。现在先简要介绍一下UDF的一些功能：
定制边界条件，定义材料属性，定义表面和体积反应率，定义FLUENT输运方程
中的源项，用户自定义标量输运方程（UDS）中的源项扩散率函数等等。
●在每次迭代的基础上调节计算值
●方案的初始化
●（需要时）UDF的异步执行
●后处理功能的改善
●FLUENT模型的改进（例如离散项模型，多项混合物模型，离散发射辐射模型）
由上可以看出FLUENT UDF并不涉及到各种算法的改善，这不能不说是一个遗憾。当然为了源代码的保密我们还是可以理解这样的做法的。其实，如果这些代码能够部分开放，哪怕就一点点，我想FLUENT会像LINUX一样发展更为迅速，使用更为广泛。遗憾的是，从目前来看，这只是一种幻想。什么时候中国人可以出自己的精品？
1.3 UDF的局限
尽管UDF在FLUENT中有着广泛的用途，但是并非所有的情况都可以使用UDF。UDF 并不能访问所有的变量和FLUENT模型。例如，它不能调节比热值；调节该值需要使用求解器的其它功能。如果您不知道是否可以用UDF解决某个特定的问题，您可以求助您的技术支持。
1.4Fluent5到Fluent6UDF的变化
如果你有FLUENT5的UDF编程经验，请注意在FLUENT6种的下列变化：
FLUENT6中加入了大量的通用多相模型。When one of these general multiphase models is enabled, storage must be set aside for the mixture as well as the individual phases. This functionality is manifested in the code through the use of additional thread and domain data structures. Consequently, some predefined macros have been added that allow access to data contained within mixture-level and phase-level domain and thread structures. See Section 3.11 for details on writing UDF for multiphase applications.
If you have a FLUENT 5 UDF with an external domain declaration that you want to use in FLUENT 6, then the extern statement must be replaced by a call to the Get_Domain utility and assignment to a Domain pointer as shown below. The Fluent-provided utility, Get_Domain(1), returns the pointer to the mixture-level domain. See Section 6.5.1for more details on Get_Domain.
Example
extern Domain *domain;
DEFINE_ON_DEMAND(my_udf)
{
...
}
is to be replaced by
DEFINE_ON_DEMAND(my_udf)
{
Domain *domain;
domain = Get_Domain(1);
...
}
The macro C_VOF accesses volume fraction values from the FLUENT solver. C_VOF(c, pt[i]) has two arguments, c and pt[i]. c is the cell identifier. pt[i] is the pointer to the corresponding phase-level thread for the i th phase, where i is the phase_domain_index.
For example, C_VOF(c,pt[i]) can be used to return the volume fraction of the i th phase fluid at cell c. The pointer pt[i] can also be retrieved using THREAD_SUB_THREAD, discussed in Section 6.5.4, using i as an argument.
?For compiled UDF, the makefile called Makefile.udf that was provided in previous FLUENT releases has been renamed to makefile.udf2. See Section 7.3.2for more details.
?For multiphase flow problems, you will need to supply your own user-defined scalar flux function instead of using the default function provided by FLUENT.
?DEFINE_PROPERTY is to be used to define UDF for particle or droplet diameter for the mixture model, previously the Algebraic Slip Mixture Model (ASMM), instead of the DEFINE_DRIFT_DIAM macro.
1.5 UDF基础
? 1.5.1输运方程
? 1.5.2单元（Cells）,面，区域（Zones）和线（Threads）
? 1.5.3操作
? 1.5.4求解器数据
? 1.5.5运行
1.5.1输运方程
FLUENT求解器建立在有限容积法的基础上，这种方法将计算域离散为有限数目的控制体或是单元。网格单元是FLUENT中基本的计算单元，这些单元的守恒特性必须保证。也就是说普通输运方程，例如质量，动量，能量方程的积分形式可以应用到每个单元：
(1.5.1)
此处，
是描述普通输运数量的变量（a general transportable quantity）,根据所求解的输运
方程它可取不同的值。下面是在输运方程中可求解的的子集。
守恒与否需要知道通过单元边界的通量。因此，需计算出单元和面上的属性值（properties）。
1.5.2单元（Cells）,面，区域（Zones）和线（Threads）
单元和单元面被组合为一些区域（zones），这些区域规定了计算域（例如，入口，出口，壁面）的物理组成（physical components）。当用户使用FLUENT中的UDF 时，用户的UDF 可调用流体区域或是边界区域的计算变量（solution variables）。UDF需要获得适当的变量，比如说是区域参考（a zone reference）和单元ID，以便标定各个单元。
区域（A zone）是一群单元或单元面的集合，它可以由模型和区域的物理特征（比如入口，出口，壁面，流体区域）来标定。例如，一些被指定为面域（a face zone）的单元面可以被指定为velocity-inlet 类型，由此，速度也就可指定了。线（A thread）是FLUENT数据结构的内部名称，可被用来指定一个区域。Thread 结构可作为数据储存器来使用，这些数据对于它所表示的单元和面来说是公用的（The Thread structure acts as a container for data that is common to the group of cells or faces that it represents）。
1.5.3操作
多数的UDF任务需要在一个线的所有单元和面上重复执行。比如，定义一个自定义轮廓函数（a custom profile function）则会对一个面线上（in a face thread）的所有单元和面进行循环。为了用户方便，Fluent Inc.向用户提供了一些循环宏工具（looping macro utilities）来执行对单元，面，节点（nodes）和线（threads）的重复操作。例如，单元循环宏（Cell-looping macros ）可以对给定单元线上的所有单元进行循环操作（loop over cells in a given cell thread allowing access to all of the cells）。而面循环宏（Face-looping macros）则可调用所有给定面线（a given face thread）的面。Fluent提供的循环工具请见Chapter 6。
在某些情况下，UDF 需要对某个变量操作，而这个变量恰恰又不能直接被当作变量来传递调用。比如，如果用户使用DEFINE_ADJUST宏来定义UDF，求解器将不会向它传递thread指针。这种情况下，用户函数需要用Fluent提供的宏来调用线指针（thread pointer）。见Chapter 6。
1.5.4求解器数据
通过FLUENT用户界面将C 函数（它已被编译和连接）连接到求解器上可实现调用求解器变量。一旦UDF和求解器正确连接，无论何时，函数都可调用求解器数据。这些数据将会被作为用户变量自动地传递给UDF。注意，所有的求解器变量，不管是求解器传递给UDF的，还是UDF传递给求解器的，都使用SI单位。
1.5.5运行
UDF将会在预定时刻被FLUENT调用。但是，也可对它们进行异步执行，使用DEFINE_ON_DEMAND宏，还可在需要时（on demand）执行。详情请见4.2.3
1.6 解释和编译UDF的比较
编译UDF和FLUENT的构建方式一样。脚本Makefile 被用来调用C编译器来构建一个当地目标代码库（a native object code library）。目标代码库包含高级C语言源代码的机器语言翻译。代码库在FLUENT运行时由“动态加载”（``dynamic loading''）过程连接到FLUENT上。连接后，与共享库的联系（the association with the shared library）将会被保存在用户的case文件中,这样，当FLUENT以后再读入case文件时，此编译库将会与FLUENT 自动连接。这些库是针对计算机的体系结构和一定版本的FLUENT使用的。所以，当FLUENT更新，或计算机操作系统改变，或是在不同类型的机器上运行时，这些库必须重新构建。
而解释UDF则是在运行时，直接从C语言源代码编译和装载（compiled and loaded directly from the C source code）。在FLUENT运行中，源代码被编译为中介的独立于物理结构的使用C预处理程序的机器代码（an intermediate, architecture-independent machine code）。当UDF 被调用时，机器代码由内部仿真器（an internal emulator），或注释器（interpreter）执行。注释器不具备标准C编译器的所有功能；它不支持C语言的某些原理（elements）。所以，在使用interpreted UDF 时，有语言限制（见3.2）。例如，interpreted UDF 不能够通过废弃结构（dereferencing structures）来获得FLUENT数据。要获得数据结构，必须使用由FLUENT提供的预定义宏。另一个例子是FLUENT interpreter 不能识别指针数组。这些功能必须由compiled UDF 来执行。
编译后，用户的C函数名称和内容将会被储存在case文件中。函数将会在读入case文件时被自动编译。独立于物理结构的代码的外层（This extra layer of architecture-independent code）可能会导致执行错误（a performance penalty），但却可使UDF共享不同的物理结构，操作系统，和Fluent版本。如果运行速度较慢，UDF不用被调节就可以编译代码的形式（in compiled mode）运行。FLUENT中的compiled 和interpreted UDF 请见Chapter 7。
选择interpreted UDF或是compiled UDF时，注意以下内容：
?Interpreted UDF
o对其它平台是便捷的（portable）。
o可作为（compiled UDF）来运行。
o不需C编译器。
o比compiled UDF慢。
o需要较多的代码。
o在使用C语言上有限制。
o不能与编译系统或用户库（compiled system or user libraries）连接。
o只能使用预定义宏来获得FLUENT结构中的数据。（见Chapters 5和6）。
?Compiled UDF
o比interpreted UDF运行快。
o在使用C语言上不存在限制。
o可用任何ANSI-compliant C 编译器编译。
o能调用以其他语言编写的函数(specifics are system- and compiler-dependent)。
o机器物理结构需要用户建立FLUENT(2D or 3D) 的每个版本的共享库（a shared library for each version of FLUENT (2D or 3D) needed for your machine
architecture）。
o如果包含有注释器（interpreter）不能处理得C语言元素，则不能作为（interpreted UDF ）运行。
总的来说，当决定使用那种类型的UDF 时：
?使用interpreted UDF作为简单的函数
?使用compiled UDF作为复杂的函数，这些函数
o对CPU有较大要求（例如每次运行时，在每个单元上均须调用的属性UDF （a property UDF）。
o需要使用编译库（require access to a compiled library）。
1.7一个step-by-stepUDF例子
编辑UDF 代码，并且在用户的FLUENT模型中有效使用它，须遵循以下七个基本步骤：
1. 定义用户模型。
2. 编制C 语言源代码。
3. 运行FLUENT，读入，并设置case文件。
4. 编译或注释（Compile or interpret）C语言源代码。
5. 在FLUENT中激活UDF。
6. 开始计算。
7. 分析计算结果，并与期望值比较。
在开始解决问题前，用户必须使用UDF 定义希望解决的问题（Step 1）。例如，加入用户希望使用UDF来定义一个用户化的边界条件（a customized boundary profile ）。用户首先需要定义一系列数学方程来描述这个条件。
接下来用户需要将这些数学方程（概念设计，conceptual design）用C语言写成一个函数（Step 2）。用户可用文本编辑器来完成这一步。以.c为后缀名来把这个文件保存在工作路径下。
写完C语言函数后，用户即可运行FLUENT并且读入或设置case文件（Step 3）。对Ｃ语言源代码进行注释，编译，和调试（interpret, compile, and debug），并在FLUENT中激活用户函数（Step 5）。最后，运行计算（Step 6），分析结果并与期望值比较。（Step 7）。根据用户对结果的分析，可将上述整个过程重复几次。具体如下。
Step 1: 定义用户模型
生成和使用UDF 的第一步是定义用户的模型方程。
如图Figure1.7.1所示的涡轮叶片。模拟叶片周围的流场使用了非结构化网格。计算域由底端的周期性边界（ a periodic boundary on the bottom ）延伸到顶端的相同部分（an identical one on the top），速度入口在左边，压力出口在右边。
Figure 1.7.1: The Grid for the Turbine Vane Example
文中对入口x速度为常数分布和抛物线分布的流场进行了比较。分段线性的分布可由边界场选项得到（the application of a profile using a piecewise-linear profile is available with the boundary profiles option），而多项式分布则只能使用用户自定义函数得到。
进口速度为常数（20 m/s）的结果如图１．７．２和１．７．３所示。当流动沿着涡轮叶片进行时，初始速度场被改变了。
Figure 1.7.2: Velocity Magnitude Contours for a Constant Inlet x Velocity
Figure 1.7.3: Velocity Vectors for a Constant Inlet x Velocity
假定现在要设涡轮叶片入口速度x不是一常数值，其分布如下
变量y在入口中心处为０．０，在入口上部和下部则分别为0.0745 m ，而入口中心处的x速度为20 m/s ，边界上为０。
用户可用UDF 描述这一分布，并将它应运到FLUENT模型中来解决这类问题。
Step 2: 编制C 语言源代码。
选定方程定义UDF 后，用户可用任意文本编辑器来书写Ｃ语言代码。以扩展名.c 保存源代码文件保存到工作路径下。关于UDF 的书写请参考Chapter ３。
下面是一个怎样在UDF中应用方程的例子。UDF的功能由主要的DEFINE 宏（the leading DEFINE macro）来定义。此处，DEFINE_PROFILE 宏用来表示下面的代码旨在给求解器提供边界的轮廓信息。书中将在以后部分讨论其它的DEFINE宏。
/*************************************************************************/
/* udfexample.c */
/* UDF for specifying a steady-state velocity profile boundary condition */
/*************************************************************************/
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, index)
{
real x[ND_ND]; /* this will hold the position vector */
real y;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
y = x[1];
F_PROFILE(f, thread, index) = 20. - y*y/(.0745*.0745)*20.;
}
end_f_loop(f, thread)
}
DEFINE_PROFILE 宏的第一个变量inlet_x_velocity 用来定义速度入口面板中的函数。名称可任意指定。在给定的边界区域上的所有单元面（identified by f in the face loop）上将会使用函数的这个方程。当用户在FLUENT用户界面选定UDF作为边界条件时，将会自动定义线（thread）。下标由begin_f_loop应用程序自动定义。UDF 中，begin_f_loop 被用来形成对边界区域上所有单元面的循环（loop through all cell faces in the boundary zone）。对于每个面，面的质心（the face centroid）的坐标可由F_CENTROID 宏来获得。抛物线方程中用
到了y坐标y ，速度值通过F_PROFILE 宏来返回给面。begin_f_loop宏和F_PROFILE 宏都是FLUENT提供的宏。详情请见Chapter５。
Step 3: 运行FLUENT，读入，并设置case文件
建立UDF后，用户开始设置FLUENT。
1. 在工作路径下启动FLUENT。
2. 读入（或设置）case文件（如果case文件以前设置过，请确认它是否被保存在了工作路径下）。
Step 4: 编译或注释（Compile or interpret）C语言源代码
这部分将例中的源代码作为interpreted UDF 来编译。注意，这个例子不可应用于Windows 的并行网络（Windows parallel networks）。完整的编译和连接UDF 请见Chapter７。
1. 确认UDF 的case 文件（如果以前设置过）和C语言源代码在工作路径下。
2. 用Interpreted UDF面板编译UDF (例如，udfexample.c) 。
Define User-Defined Functions Interpreted...
Figure 1.7.4: The Interpreted UDF Panel
(a) 在Source File Name下键入C语言源代码文件(例如，udfexample.c) 。
!!如果用户源代码不在目前工作路径下，则在编译UDF时，需在Interpreted UDF 面板中间如文件完整的路径。
(b) 在CPP Command Name一栏里，选择C预处理器。
(c) Stack Size缺省设置为10000。如果用户函数的局部变量数目大于10000，将会
导致堆栈溢出。这种情况下，应将Stack Size设置为比局部变量大的数。
(d) 点击Compile编译UDF。
存储case文件时，C语言代码的名称和内容将会储存在case文件中。如果Interpreted UDF面板中的Display Assembly Listing选项被选中，当编译进行时，控制台窗口中将会显示汇编语言代码。另外，此选项将会被存储在case文件中，当以后用户在执行FLUENT任务时，控制台窗口中将会编译时一样显示汇编语言代码。
inlet\_x\_velocity:
.local.pointer thread (r0)
.local.int nv (r1)
0 .local.end
0 save
.local.int f (r3)
1 push.int 0
.local.pointer x (r4)
3 begin.data 8 bytes, 0 bytes initialized:
7 save
. .
. .
156 pre.inc.int f (r3)
158 pop.int
159 b .L3 (22)
.L2:
161 restore
162 restore
163 ret.v
!!注意，如果编译失败，FLUENT将会给出错误信息，请调试程序。详见7.2.3。
(e) 编译结束后，点击Close。
!!此例中的UDF源代码也可作为compiled UDF运行。
Step 5: 在FLUENT中激活UDF
编译连接完UDF后，FLUENT用户界面面板中将会看到UDF。此例中，可在Velocity Inlet面板中选择UDF。
Define Boundary Conditions...
在X-Velocity下拉列表中，选择udf inlet_x_velocity，此名称是由例中的函数给定的。一旦选中，UDF将会替代X-Velocity中的0 值进行运算。点击OK接受新的边界条件，关闭面板。
Step 6: 开始计算
运算方式和以前一样。
Solve Iterate...
Step 7: 分析计算结果，并与期望值比较
计算收敛后，获得一个修正的速度场。Figure1.7.5为入口x速度为抛物线分布的速度等值线，可与（Figure 1.7.2）所示的入口速度为常数20 m/sec 的流场比较。常数条件下，流动在涡轮叶片周围变形（distorted）。入口抛物线分布，在入口中心处为最大值，边缘上为0。
Figure 1.7.5:Velocity Magnitude Contours for a Parabolic Inlet x
Velocity
第二章.UDF的C语言基础
本章介绍了UDF的C语言基础
2.1引言
2.2注释你的C代码
2.3FLUENT中的C数据类型
2.4常数
2.5变量
2.6自定义数据类型
2.7强制转换
2.8函数
2.9数组
2.10指针
2.11声明
2.12常用C操作符
2.13C库函数
2.14用#define实现宏置换
2.15用#include实现文件包含
2.16与FORTRAN比较
2.1引言
本章介绍了C语言的一些基本信息，这些信息对处理FLUENT的UDF很有帮助。本章首先假定你有一些编程经验而不是C语言的初级介绍。本章不会介绍诸如while-do循环，联合，递归，结构以及读写文件的基础知识。如果你对C语言不熟悉可以参阅C语言的相关书籍。
2.2注释你的C代码
熟悉C语言的人都知道，注释在编写程序和调试程序等处理中是很重要的。注释的每一行以“/*”开始，后面的是注释的文本行，然后是“*/”结尾
如：/* This is how I put a comment in my C program */
2.3 FLUENT的C数据类型
FLUENT的UDF解释程序支持下面的C数据类型：
注意：UDF解释函数在单精度算法中定义real类型为float型，在双精度算法宏定义real为double型。因为解释函数自动作如此分配，所以使用在UDF中声明所有的float和double 数据变量时使用real数据类型是很好的编程习惯。
2.4常数
常数是表达式中所使用的绝对值，在C程序中用语句#define来定义。最简单的常数是十进制整数（如：0，1，2）包含小数点或者包含字母e的十进制数被看成浮点常数。按惯例，常数的声明一般都使用大写字母。例如，你可以设定区域的ID或者定义YMIN和YMAX 如下：#define WALL_ID 5
#define YMIN 0.0
#define YMAX 0.4064
2.5变量
变量或者对象保存在可以存储数值的内存中。每一个变量都有类型、名字和值。变量在使用之前必须在C程序中声明。这样，计算机才会提前知道应该如何分配给相应变量的存储类型。
2.5.1声明变量
变量声明的结构如下：首先是数据类型，然后是具有相应类型的一个或多个变量的名字。变量声明时可以给定初值，最后面用分号结尾。变量名的头字母必须是C所允许的合法字符，变量名字中可以有字母，数字和下划线。需要注意的是，在C程序中，字母是区分大小写的。下面是变量声明的例子：
int n; /*声明变量n为整型*/
int i1, i2; /*声明变量i1和i2为整型*/
float tmax = 0.; /* tmax为浮点型实数，初值为0 */
real average_temp = 0.0; /* average_temp为实数，赋初值为0.1*/
2.5.2局部变量
局部变量只用于单一的函数中。当函数调用时，就被创建了，函数返回之后，这个变量就不存在了，局部变量在函数内部（大括号内）声明。在下面的例子中，mu_lam和temp是局部变量。
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu_lam;
real temp = C_T(cell, thread);
if (temp > 288.)
mu_lam = 5.5e-3;
else if (temp > 286.)
mu_lam = 143.2135 - 0.49725 * temp;
else
mu_lam = 1.;
return mu_lam;
}
2.5.3全局变量
全局变量在你的UDF源文件中是对所有的函数都起作用的。（调用一个UDF源文件可能会包括一系列的连接函数。）它们是在单一函数的外部定义的。全局变量一般是在预处理程序之后的文件开始处声明。
2.5.4外部变量
如果全局变量在某一源代码文件中声明，但是另一个源代码的某一文件需要用到它，那么你必须在另一个文件中声明它是外部变量。外部变量的声明很简单，你只需要在变量声明的最前面加上extern即可。如果有几个文件涉及到该变量，最方便的处理方法就是在头文件（.h）中加上extern的定义，然后在所有的.c文件中引用该头文件即可。只有一个.c文件应该包括没有extern关键字的变量声明，如下所示。注意：extern只用于编译过的UDF。
例子：
/* filea.h */
/*包含外部定义的头文件*/
extern real volume;
/* filea.c */
/*调用头文件filea.h中声明的volumn的C函数*/
#include "udf.h"
#include "filea.h"
real volume;
DEFINE_ADJUST(compute_volume, domain)
{
/*计算某些区域volumn的代码*/
volume = ....
}
/* fileb.c */
/*调用头文件filea.h中声明的volumn的另一个C函数*/
#include "udf.h"
#include "filea.h"
DEFINE_SOURCE(heat_source,c,t,ds,eqn)
{
/* 用总数来计算每个单位体积的源项的代码*/
/*fliea.c的compute_volum计算出的volume*/
real total_source = ...;
real source;
source = total_source/volume;
return source;
}
2.5.5静态变量
static声明对于全局变量和局部变量的影响是不一样的。静态局部变量在函数调用返回之后，该变量不会被破坏。静态全局变量则在定义该变量的.c源文件之外对任何函数保持不可见。静态声明也可以用于函数，使该函数只对定义它的.c源文件保持可见。下面是静态全局变量声明的例子。注意：extern只用于编译过的UDF。
例子：
#include "udf.h"
static real abs_coeff = 1.0; /* 吸收系数*/
real source;
DEFINE_SOURCE(energy_source, c, t, dS, eqn)
{
int P1 = ....;
dS[eqn] = -16.* abs_coeff * SIGMA_SBC * pow(C_T(c,t),3.);
source =-abs_coeff *(4.* SIGMA_SBC * pow(C_T(c,t),4.) - C_UDSI(c,t,P1));
return source;
}
DEFINE_SOURCE(p1_source, c, t, dS, eqn)
{
int P1 = ...;
dS[eqn] = -abs_coeff;
source = abs_coeff *(4.* SIGMA_SBC * pow(C_T(c,t),4.) - C_UDSI(c,t,P1));
return source;
}
2.6自定义数据类型
C还允许你用结构和typedef创建自定义数据类型。下面是一个结构列表的定义。注意：typedef 只用于编译过的UDF。
例子：
typedef struct list
{
int a;
real b;
int c;
}
mylist; /* mylist为类型结构列表*/
mylist x,y,z; /* x,y,z为类型结构列表*/
2.7强制转换
你可以通过强制转换将某一数据类型转换为另一种。强制由类型来表示，其中的类型包括int，float等等，如下例所示：
int x = 1;
real y = 3.14159;
int z = x+((int) y); /* z = 4 */
2.8函数
函数是用完成一定任务的一系列语句。在定义该函数的同一源代码中，这些任务可能对其它的函数有用，也可能会被用于完成源文件以外的函数中。每个函数都包括一个函数名以及函数名之后的零行或多行语句，其中有大括号括起来的函数主体可以完成所需要的任务。函数可以返回特定类型的数值。C函数通过数值来传递数据。
函数有很多数据类型，如real，void等，其相应的返回值就是该数据类型，如果函数的类型是void就没有任何返回值。要确定定义UDF时所使用的DEFINE宏的数据类型你可以参阅udf.h文件中关于宏的#define声明一节，也可以参阅附录A的列表。
!! C函数不能改变它们的声明，但是可以改变这些声明所指向的变量。
2.9 数组
数组的定义格式为：名字[数组元素个数]，C数组的下标是从零开始的。变量的数组可以具有不同的数据类型。
例子
int a[10], b[10][10];
real radii[5];