一个简单而完整的UDS例子

1、一个简单而完整的UDS例子一个简单而完整的UDS例子gearboy假定需要求解方程为乎+ dt标准方程形式:即如果计算达到稳定时，？dx/dxj)=2,此名 称表 达 式各项在UDS 中需要定义 的量在UDS中对应的宏离散形式的非稳态 项非稳态项dp(/) dT一apusuDEFINE_UDS_UNSTE ADY（详细见UDF帮助）对 流 项a f 叫时通量VDEFINE_UDS_FLUX扩散项-宵扩散率rDEFINE_DIFFUSIVITY=2dp(b d_ 厂 7- _jdtdxj)瓦IFluent中各项的定义如下:d L g（b-I方程解析解为心C,其中C为常数（假定其为0）,可由边界条

2、件确定。因此，当解稳定时，Otx+v，等值线为一系列45度的斜直线。 将此方程与标准方程对比发现：名 称表 达 式各项在 UDS中需 要定义的 量需要在UDS中定义的宏非稳态项g dt离散形式的非稳态 项ttDEFINE_UDS_UNSTEADY对 流 项J通量诉3.冒（单位向量 7 = m）DEFINE_UDS_FLUX扩散项冷扩散率r = iDEFINEDIFFUSIVITY源 项22DEFINESOURCE边界条件c = 0S, =xZh（b代表边 界处）DEFINE PROFILE*注：源项和扩散率为常数，可以不用定义宏，但为了说明问题，仍采用宏来定义。假定如图所示的计算域，其中正方形

3、边长为10, 原点在正方形中心，四边均为wall。将如下源码，写入一个UDF文件，例如MyFistUDF.c，与cas文件放在同一目录。#include udf.hDEFINE_UDS_UNSTEADY(MyUnsteady皿i,a pu,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_oldC_STORAGE_R(c,t,SV_UDSI_M 1(

4、i);*su =/explicit part*/vol*phi_old/physical_dt;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);积单位矢l=J矢量点l=JDEFINE_DIFFUSIVITY(MyDiff,c,t,i)return 1.0;DEFINE_SOURCE(MySource,c,t,dS,eqn)dSeqn=0;retur

5、n 2.0;DEFINE_PROFILE(MyProfile,thread,index)real xND_ND; /* this will hold the positionvector */real xx,yy;face_t f;begin_f_loop(f,thread)F_CENTROID(x,f,thread);xx = x0;yy=xi;F_PROFILE(f,thread,index) = (xx+yy);II4 = x + yend_f_loop(f,thread)编译通过并载入启动非稳态求解器，并在Define-User-Defined-Scalars 面板中添加一个UDS，选

6、择自定义的通量和非稳态项在材料面板中，选择自定义扩散率在boundary conditions的Fluid面板中，选择自定义源项在boundary conditions的壁面边界条件面板中，选择自定义边界条件在 solve-controls-solution 中禁止 Flow 计算，只计算自定义标量迭代过程中的残差曲线，可以看出正逐渐趋向稳采用后处理等值线面板，观察自定义标量的等值ContoursOptions匚 FilledP Node Values 匠 Global Range 口 Auto Range 厂 Clip to RangeDraw Profiles Draw Grid_evels SetupSurface Name PatternContours ofUser Defined Scalars.User Scalar 0dMinMax|-9.9|9.9Surfaces如default-interior wallSurface Typesaxisdclip-surf exhaust-fanfanX:Displ