梅花形纵向折流板换热器的流动和换热特性分析研究.docVIP

梅花形纵向折流板换热器的流动和换热特性分析研究 1、相关定义 1.1、缝隙流动研究中用户自定义程序理论 Fluent 中用户自定义函数为二次程序开发,也称 UDF(user defined functions)。 基于标准 C 语言编写程序并动态链接到求解器,程序编写时可以使用 C 语言的 库函数或者使用Fluent提供的内部预定义宏来获得求解器存储器数据。采用UDF 是解决 Fluent 标准界面上不能解决的物理问题,如通过 UDF 定值特殊的边界条 件,定制材料属性,定制输运方程源项等问题。比较复杂的是在 Fluent 的求解过 程,可以通过 UDF 来调节每次迭代的计算值,同时可以初始化求解方案和求解 北京工业大学工学硕士学位论文 20 模型的改进等。但其局限是源代码暗箱操作,不能进行算法的优化。不过 UDF 来定义实现的数值计算功能是很丰富的,用户可以用 UDF 来定义实现以下方面 的程序定义:边界条件、源项、物性定义、表面和体积反应速率、用户自定义标 量输运方程、离散相模型(如体积力,拉力,源项等)、代数滑流混合物模型变 量初始化、壁面热流量以及供用户自定义标量后处理等。 UDF 在使用和编译时可以当作解释性函数和编译函数处理。前者是在求解运 行中读入并解释,后者是编译为共享库并被条用。只要函数成功编译和链接,用 户自定义函数任何时候都可以被求解器调用,同时实现求解器和 UDF 间的数据 交换和传递,编译的 UDF 和 Fluent 都是使用 makefile 脚本调用目标代码库,代 码库能动态加载链接到 Fluent 求解器中,Fluent 本身是由独立的物理结构机器源 代码。 Fluent 的数值求解基础是有限体积法,将求解计算域离散为有限数目的控制 体和单元。作为求解器计算的基本网格单元,必须保证物理守恒。这就需要满足 输运方程中的质量、动量、能量方程在每个单元都能有积分形式: ??t ∫VρφdV + ∫ A ρ φV ? dA = ∫ A Γ? φ ? dA + ∫V S φ dV (2-16) 其中φ 为输运变量,可能为表示下表中任何变量。输运方程的守恒就必需计 算单元边界的通量,那就需要计算单元和面上的相应的值。相应就提出离散出来 的单元信息和计算思想。 在数值求解时需要把物理空间离散为不同模型和流动空间组成的区域 (zones),在把区域离散为单元面(faces)和单元(cells)。对于不同的区域就可 以用 zones 线程(thread)来标定,单元和面单元也可以被 thread 标定。线程(thread) 是 Fluent 内部数据结构名称,Thread 结构以栈的形式被数据存储器使用。其数 据结构栈所表示的单元和面是公用的,在 UDF 用户自定义时就是通过条用区域 或边界区域的标识和单元 ID 来获取和定义变量的。通常在 UDF 中需要循环所有 单元和面上的线程,Fluent 提供了对单元、面、节点和线程执行循环操作的循环 宏工具。如单元循环宏可以对指定单元线上的所有单元进行循环操作,面循环宏 可调用指定面线的面。 Fluent 中提供连接矢量的宏,允许通过宏来连接单元形心以及单元与面心间 的矢量联系和调用。此宏返回单元边界处面通量值,如边界单元的扩散通量。为 了更好的了解宏返回的信息参数,需要通量是怎么计算的。假设存在一个可用的 梯度标量φ ,则面单元上的标量面值为:φf= φ0 + ?φ ? dr ,dr 为指向面上一单 元矢量, φ0 为单元形心处的梯度。面单元 f 上的通量φ ,则其扩散通量D f 就应该 第 2 章 考虑滑移边界的微缝隙流动及数值方法 2 1 为:D f= Γ f ?φ i A ,其中Γ f 是面 f 上扩散系数,如图 5-2 中所示。对于非结构网 格求解的非结构求解器而言,面单元法向梯度可以通过计算相邻单元与面限定方 面的沿平面内方向来解算。具体扩散通量D f 为: f f(1 0 ) f ( s ) s s DA A A e A A = Γφd? sφ A ?? e +Γ ?φ ? ? ?φ ? A ?? e (2-17) 其中第一项代表沿矢量e s 的主梯度,第二项表示扩散项。 A为 C0 单元指向 C1 单元的 f 面法向面积,ds 为相邻单元形心间距,e s 法向单位矢量。? φ 是相邻两 个单元梯度平均。 图 2-5 相邻 C0 和 C1 单元上的矢量和梯度定义 Fig.2-5 Define of vector and gradient of adjacent cells 1.2、流动模型的定义 (l)导入网格模型 利用 import 命令导入在 ICEM CFD 中建立的流体网格 模型; (2)定义流体域模型 单流体域,流体介质为水,流体域模型参考