利用FLUENT二次开发(UDF)辅助科学计算.pptVIP

利用FLUENT二次开发(UDF)辅助科学计算 龚 磊 2007.4.26 引言: 随着计算机硬件的发展及CFD软件的成熟,采用CFD商业软件进行工程计算及科研辅助日趋流行.相比传统的编程,CFD软件计算具有建模快,收敛性好,移植性强等优点,因而可以节省出更多时间花在问题本质上,而不是计算实现上. CFD软件使用没有降低对流体基本知识的掌握，只是自己不用离散方程了，数学上的要求降低了。 FLUENT二次开发-UDF UDF(用户自定义程序)是FLUENT提供的二次开发接口,通过大量的宏定义来实现用户和SOLVER中数据的交互的一种方法，来增强FLUENT的功能及模型的应用范围:比如边界条件,修改控制方程,求解自定义方程等等. 采用UDF接口我们可以定义求解UDS(用户自定义标量方程) 举例实现:双电层理论求解 双电层理论的控制方程(P-N-P方程): 解决方案: 采用三个UDS方程,分别对于以上三个方程, 三个标量方程中的 在这里分别代表电位势 ,正离子浓度 , 负离子浓度 . 具体步骤一:方程规范化 将双电层的三个方程分别与标志输运方程对比: 具体步骤一:方程规范化 正负离子方程与标准输运方程对比 具体步骤二:编写UDF文件 取名为uds-flux.c,内容如: #include “udf.h” DEFINE_UDS_FLUX(UDS0_FLUX,f,t,i){…} DEFINE_UDS_FLUX(UDS1_FLUX,f,t,i){…} DEFINE_SOURCE(POSSION_SOURCE,c,t,dS,eqn){…} DEFINE_UDS_UNSTEADY(uds0_unsteady,c,t,i,apu,su){…} DEFINE_UDS_UNSTEADY(uds0_unsteady,c,t,i,apu,su){…} ………… 具体步骤三: FLUENT设置 1,启动FLUENT,读入CAS文件 2,编译编写好的UDF文件uds-flux.c 菜单:defineuser-definedfunctionscomplied 先点击build编译,再点击load读入库. 3:添加三个UDS方程: 菜单:defineuser-definedscalars 默认为0个,翻动上下按钮添加UDS数目到1 UDS indix 指UDS的编号,从0开始,UDS0对应第一个UDS方程,以此类推. 继续翻动number of user-defined scalars右边的上下按钮添加第二个UDS方程(正离子方程) 继续翻动number of user-defined scalars右边的上下按钮添加第三个UDS方程(电位势方程) 点击OK,系统会弹出对话筐显示UDS已经添加成功,并提示到材料面板设置UDS扩散系数D. 再次点击OK确认. 转到材料面板,设置扩散系数D: 菜单:definematerials 弹出菜单如下: 点击UDS Diffusivity右边的Edit按钮,分别修改三个方程对应的扩散系数. 4:设置边界条件:UDS边界条件及添加源项 菜单:defineBoundary Conditions 源项fluid,点击set,弹出左边的面板.因为只有电位势方程有源项,所以点击user scalar 2旁边的edit来添加源项. 添加源项个数为1,选择POSSION_SOURCE 边界条件wall/inlet/outlet…… Fluent内置的UDS边界条件有两种:第一类边界条件:给定标量值;第二类边界条件:给定流通量.如图电位势方程,给定了壁面zeta电位为25.7mv. 5:求解控制 菜单Solvecontrolssolution… 6:初始化,迭代求解. 菜单solveinitializeinitialize solveIterate 至此,整个UDS设置求解流程已经完成,当然不同情况下略有不同,如选择了Unsteady求解器,添加UDS的时候弹出的菜单如下: CFD软件只是实现工具,深入理解所研究问题的本质才是根本. 谢谢大家! * * FLUENT应用范围 航空 航天 船舶 化工 汽车 能源 生物 其它 FLUENT的应用范围“相当广泛” 2.能源 4.船舶 1.航空 3.汽车 UDS理论及实现 标量输运方程： 非定常项 对流项 扩散项 源项 可通过UDF修改 DEFINE_UDS_UNSTE