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

上传人:王** 文档编号:895747 上传时间:2024-02-22 格式:DOCX 页数:4 大小:16.14KB
下载 相关 举报
UDF程序解释-指定自定义标量梯度函数.docx_第1页
第1页 / 共4页
UDF程序解释-指定自定义标量梯度函数.docx_第2页
第2页 / 共4页
UDF程序解释-指定自定义标量梯度函数.docx_第3页
第3页 / 共4页
UDF程序解释-指定自定义标量梯度函数.docx_第4页
第4页 / 共4页
亲,该文档总共4页,全部预览完了,如果喜欢就下载吧!
资源描述

《UDF程序解释-指定自定义标量梯度函数.docx》由会员分享,可在线阅读,更多相关《UDF程序解释-指定自定义标量梯度函数.docx(4页珍藏版)》请在优知文库上搜索。

1、002指定自定义标量梯度函数下面这段程序,作用是使用一个用户自定义的标量作为函数,来表示另一个用户定义标量的梯度的函数,有点绕啊。这个UDF每一次迭代调用一次,执行的时候用汇编的方式。看程序:/*UDFfordefininguser-definedscalarsandtheirgradients*/#includeudf.hDEFINE_ADJUST(adjust_fcn,d)(Thread*t;cell_tc;realK_EL=1.0;*Donothingifgradientisntallocatedyet.*/if(!Data_Valid_P()return;thread_loop_c(t

2、,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)#includedf.h是一个声明,意思是后面告诉计算机,后面这段程序要调用的“udf.h”里面包含的函数。“udf.h”是一个头文件,是每一个UDF程序必须要包含的,形式是includedf.h,出现在每条程序最开始部分。就是一个库函数,包含了DEFlNE宏的定义,程序提供其他的的宏和函数,比如三角函数,鬲指数函数等等。以汇编的形式编入程序,是写UDF必

3、须有的。如果是新手的写的话,直接第一句写这个,肯定没错。DEFINE_ADJUST(adjuSJfCn,d)是一个DEFINE宏。DEFlNE_AD:JUST是一个通用宏,可以用来调整或修改ANSYSFluent中未作为参数传递的变量。简单说一句,所谓宏,就是一系列程序的合集,打包之后可以方便调用,对于用户来讲可以理解为黑箱,只要提供参数,就能得到结果。你不需要去了解宏内部是怎么运行的,只要知道输入需要什么,输出能给你什么就可以了。当然如果你要了解宏的内部,查手册和资料,看代码都是可以的。Thread*t;Thread,线程,是fluent中的一种结构数据类型,用于存储它所表示的单元格或面组的

4、公共数据。简单可以理解为一种数据类型。整句含义是定义了一个指针,(这里需要一些C语言知识,指针的定义),指针的名称是t,存储数据类型是Threadocelljc;定义了一个celLt的变量,名字叫C;cell_t是fluent中的一种数据类型,表示在单元格线程中标识特定单元格的整数数据类型。realK_EL=L0;定义了一个实数变量,名字叫K_EL,并赋值为1。定义好了需要的变量,下面开始程序。先说下这个程序的总体思路。先做一个检测,如果梯度还没有指派,那什么都不做。之后开始在域上循环每个线程,之后在每个线程里的单元进行梯度修改。if(!Data_Valid_P()return;Data_Va

5、IiCLP()是用来检测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)这是在单元线程里循

6、环每个单元,之后进行操作,由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)表示该单元的体积。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 高等教育 > 大学课件

copyright@ 2008-2023 yzwku网站版权所有

经营许可证编号:宁ICP备2022001189号-2

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!