博研联盟_一个简单而完整的UDS例子

PAGE\*MERGEFORMAT#一个简单而完整的UDS例子gearboy方程形式:型+_^阮—「之、=se（默认侦=^）ct£XjIcXjJFluent中各项的定义如下:名称达式各项在UDS中需要定义的量在UDS中对应的宏非稳态项cP*ct离散形式的非稳态项a屯AtwAt'—.—-'3papuiisu（详细见UDF帮助）DEFINE_UDS_UNSTEADY对流项冰j通量旷DEFINE_UDS_FLUX扩散项12cXj扩散率rDEFINE_DIFFUSIVITY同巾尸J假定需要求解方程为—+—寺=2日跖I次j）当tt七，即如果计算达到稳定时，——*-——=2，此方程解析解为*=x+y+c,欲jI欲j）其中c为常数（假定其为0）,可由边界条件确定。因此，当解稳定时，巾tx+y,等值线为一系列45度的斜直线。将此方程与标准方程对比发现:名称表达式各项在UDS中需要定义的量需要在UDS中定义的宏非稳态项冲离散形式的非稳态项一算—义膏△tAtDEFINE_UDS_UNSTEADY对流项礼通量旷=I.a（单位向量I=G,i,i｝）DEFINE_UDS_FLUX扩散项c2*'2以j扩散率r=1DEFINE_DIFFUSIVITY源项22DEFINE_SOURCE边界条件c=0%=Xb+yb（b代表边界处）DEFINE_PROFILE*注：源项和扩散率为常数，可以不用定义宏，但为了说明问，仍采用宏来定义。假定如图所示的计算域，其中正方形边长为10,原点在正方形中心，四边均为wall。将如下源码，写入一个UDF文件，例如MyFistUDF.c，与cas文件放在同一目录。#include"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[1];F_PROFILE(f,thread,index)=(xx+yy);〃如=xb十yb}end_f_loop(f,thread)}编译通过并载入启动非稳态求解器，并在Define-User-Defined-Scalars面板中添加一个UDS,选择自定义的通量和非稳态项User-DefinedScalar^NumberofUser-DefinedScalarsFluxFunctionMyFlux：：libudtUnsteadyFunctionMyUnsteady::libOK|Cancel|Hkip|在材料面板中，选择自定义扩散率MismeMaterialTypeairChemicalFormulafluid*FluentFluidMaterialsloirzlMixture|nonc」lalerialsOrderMa厚NamLChenFlueiUser-DePropertiesDensity(kg/m3]'.r131constant1.225Viscosity(kg/m-s]'.7131constant1.789Ue-05UDSDiffusivity(kg/m-s]dEdit-,r|U^er-DefinedFunctionsMyUnsteady;:libudfMvDiff；；libudfMySaurce：：libudfMyProfile：：libudf在boundaryconditions的Fluid面板中，选择自定义源项iCa^ellHelpChange/CreateDeleteCloseIHelp|在boundaryconditions的壁面边界条件面板中，选择自定义边界条件在solve-controls-solution中禁止Flow计算，只计算自定义标量迭代过程中的残差曲线，可以看出正逐渐趋向稳定采用后处理等值线面板，观察自定义标量的等值线_|口|XJ未完全稳定时的等值线^FLUENT[0]Flueritfil£心±.19^01SMII沛血11俱㈤[妣flUdta-Wit36c-00§.峰・00■网L7fc*00455<-W-啊皿土加“的4浜*40曷X心莓如皿Hvie,A007E1U516.3(3d.侃qi】即viSLkA,最终稳定时的等值线，可以看出是一系列45度的斜直线，与解析解相同。但右上角顶点的值为9.9,而解析解中应为x+y=10，所以尚存在1%的误差。