PAGE\*MERGEFORMAT#一个简单而完整的UDS例子gearboyPAGE\*MERGEFORMAT#标准方程形式:S(默认UJtXjXjFluent中各项的定义如下:名称表达式各项在UDS中需要定义的量在UDS中对应的宏非稳态项t离散形式的非稳态项pArap
要求
对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗
解方程为2tXjXj当t,即如果计算达到稳定时,XjXj2,此方程解析解为xyc,其中c为常数(假定其为0),可由边界条件确定。因此,当解稳定时,xy,等值线为一系列45度的斜直线。将此方程与标准方程对比发现:名称表达式各项在UDS中需要定义的量需要在UDS中定义的宏非稳态项t离散形式的非稳态项VnVn1ttDEFINE_UDS_UNSTEADY对流项Xj诵量1A(单位向量11,1,1)DEFINE_UDS_FLUX扩散项22Xj扩散率1DEFINE_DIFFUSIVITY源项22DEFINESOURCE边界条件c0bXbyb(b代表边界处)DEFINE_PROFILE*注:源项和扩散率为常数,可以不用定义宏,但为了说明问题,仍采用宏来定义。假定如图所示的计算域,其中正方形边长为10,原点在正方形中心,四边均为wall。将如下源码,写入一个UDF文件,例如MyFistUDF.c,与cas文件放在同一目录。#inelude"udf.h"DEFINE_UDS_UNSTEADY(MyUnsteady,c,t,i,apu,su){realphysical_dt,vol,phi_old;physical_dt=RP_Get_Real("physical-time-step");vol=C_VOLUME(c,t);*apu=-vol/physical_dt;/implicitpart*/phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i));*su=vol*phi_old/physical_dt;/*explicitpart*/}DEFINE_UDS_FLUX(MyFlux,f,t,i){realNV_VEC(unit_vec),NV_VEC(A);〃声明矢量变量F_AREA(A,f,t);NV_DS(unit_vec,=,1,1,1,*,1);//单位矢量赋值returnNV_DOT(unit_vec,A);〃矢量点积}DEFINE_DIFFUSIVITY(MyDiff,c,t,i){return1.0;}DEFINE_SOURCE(MySource,c,t,dS,eqn){dS[eqn]=0;return2.0;DEFINE_PROFILE(MyProfile,thread,index){realx[ND_ND];/*thiswillholdthepositionvector*/realxx,yy;face_tf;begin_f_loop(f,thread){F_CENTROID(x,f,thread);xx=x[0];yy=x[i];F_PROFILE(f,thread,index)=(xx+yy);//bxbyb}end_f_loop(f,thread)}编译通过并载入启动非稳态求解器,并在Define-User-Defined-Scalars面板中添加一个UDS,选择自定义的通量和非稳态项在材料面板中,选择自定义扩散率slameMaterialTypeOrderMaairpluid一zl'■*NaimChemicalFormulaFluentFluidMaterials广Chen|airdFlueiMixlureUser-DenonePropertiesDensity(kg/m3]constantEdit」Viscosity(kg^m-sj|1.225User-DefinedFunctionsconstantMyUnsteady::libudfMyFlux::libudfMyDiff;:libudfUDSDiffusivity(Itg/m-s)|1.7894e-85MySource::libudfMyProfile::libudf:CanceHHelpChangc/CrealeDeleteCloseHelp在boundaryconditions的Fluid面板中,选择自定义源项在boundaryconditions的壁面边界条件面板中,选择自定义边界条件在solve-controls-solution中禁止Flow计算,只计算自定义标量-|匚|x|迭代过程中的残差曲线,可以看出正逐渐趋向稳定采用后处理等值线面板,观察自定义标量的等值线未完全稳定时的等值线^FLUTIT[D]FlupidnlXriL33-L^I'卅如L8^->QI1!»-bl<1M-G1l呱心,--•■十4Pa*N-JSLrKW斗葩时C->iiUaiivlUaH.-lmqfeid0EkZ0099^00)EkUDTI£.勺竝玄m.nrthv=K>)最终稳定时的等值线,可以看出是一系列45度的斜直线,与解析解相同。但右上角顶点的值为9.9,而解析解中应为x+y=1O,所以尚存在1%的误差。