好文档就是一把金锄头!
欢迎来到金锄头文库![会员中心]
电子文档交易市场
安卓APP | ios版本
电子文档交易市场
安卓APP | ios版本

一个简单而完整的UDS例子.doc

6页
  • 卖家[上传人]:枫**
  • 文档编号:417088406
  • 上传时间:2023-02-07
  • 文档格式:DOC
  • 文档大小:233KB
  • / 6 举报 版权申诉 马上下载
  • 文本预览
  • 下载提示
  • 常见问题
    • 一个简单而完整的UDS例子 gearboy标准方程形式: (默认)Fluent中各项的定义如下:名称表达式各项在UDS中需要定义的量在UDS中对应的宏非稳态项离散形式的非稳态项(详细见UDF帮助)DEFINE_UDS_UNSTEADY对流项通量 DEFINE_UDS_FLUX扩散项扩散率 DEFINE_DIFFUSIVITY假定需要求解方程为当,即如果计算达到稳定时,,此方程解析解为,其中c为常数(假定其为0),可由边界条件确定因此,当解稳定时,,等值线为一系列45度的斜直线将此方程与标准方程对比发现:名称表达式各项在UDS中需要定义的量需要在UDS中定义的宏非稳态项离散形式的非稳态项DEFINE_UDS_UNSTEADY对流项通量(单位向量)DEFINE_UDS_FLUX扩散项扩散率 DEFINE_DIFFUSIVITY源项22DEFINE_SOURCE边界条件 (b代表边界处)DEFINE_PROFILE*注:源项和扩散率为常数,可以不用定义宏,但为了说明问题,仍采用宏来定义假定如图所示的计算域,其中正方形边长为10,原点在正方形中心,四边均为wall。

      将如下源码,写入一个UDF文件,例如MyFistUDF.c,与cas文件放在同一目录include "udf.h"DEFINE_UDS_UNSTEADY(MyUnsteady,c,t,i,apu,su){real physical_dt, vol, phi_old;physical_dt = RP_Get_Real("physical-time-step");vol = C_VOLUME(c,t);*apu = -vol / physical_dt; /*implicit part*/phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(i));*su = vol*phi_old/physical_dt; /*explicit part*/}DEFINE_UDS_FLUX(MyFlux,f,t,i){real NV_VEC(unit_vec), NV_VEC(A); //声明矢量变量F_AREA(A, f, t);NV_DS(unit_vec, =, 1, 1, 1, *, 1); //单位矢量赋值return NV_DOT(unit_vec, A); //矢量点积}DEFINE_DIFFUSIVITY(MyDiff,c,t,i){return 1.0;}DEFINE_SOURCE(MySource,c,t,dS,eqn){ dS[eqn]=0; return 2.0;}DEFINE_PROFILE(MyProfile,thread,index){real x[ND_ND]; /* this will hold the position vector */real xx,yy;face_t f;begin_f_loop(f,thread) {F_CENTROID(x,f,thread);xx = x[0];yy=x[1];F_PROFILE(f,thread,index) = (xx+yy); //}end_f_loop(f,thread)}编译通过并载入启动非稳态求解器,并在Define-User-Defined-Scalars面板中添加一个UDS,选择自定义的通量和非稳态项在材料面板中,选择自定义扩散率在boundary conditions的Fluid面板中,选择自定义源项在boundary conditions的壁面边界条件面板中,选择自定义边界条件在solve-controls-solution中禁止Flow计算,只计算自定义标量迭代过程中的残差曲线,可以看出正逐渐趋向稳定采用后处理等值线面板,观察自定义标量的等值线未完全稳定时的等值线最终稳定时的等值线,可以看出是一系列45度的斜直线,与解析解相同。

      但右上角顶点的值为9.9,而解析解中应为x+y=10,所以尚存在1%的误差。

      点击阅读更多内容
      关于金锄头网 - 版权申诉 - 免责声明 - 诚邀英才 - 联系我们
      手机版 | 川公网安备 51140202000112号 | 经营许可证(蜀ICP备13022795号)
      ©2008-2016 by Sichuan Goldhoe Inc. All Rights Reserved.