[整理版]博研联盟_一个简单而完整的uds例子

[整理版]博研联盟_一个简单而完整的uds例子[整理版]博研联盟_一个简单而完整的uds例子一个简单而完整的UDS例子 gearboy标准方程形式: ,,,,,,,,,,,,, (默认) ,,,u，,,,S,,,i,,,t,x,xjj,, Fluent中各项的定义如下: 名称表达式各项在UDS中需要定义的量在UDS中对应的宏离散形式的非稳态项 ,,, 非稳态项 DEFINE_UDS_UNSTEADY ,t (详细见UDF帮助) ,,,,, 通量 ,对流项 DEFINE_UDS_FLUX ,xj 2,, ,,扩散率 ,扩散项 DEFINE...

[整理版]博研联盟_一个简单而完整的uds例子一个简单而完整的UDS例子 gearboy 标准方程形式: ,,,,,,,,,,,,, (默认) ,,,u，,,,S,,,i,,,t,x,xjj,, Fluent中各项的定义如下: 名称表达式各项在UDS中需要定义的量在UDS中对应的宏离散形式的非稳态项 ,,, 非稳态项 DEFINE_UDS_UNSTEADY ,t (详细见UDF帮助) ,,,,, 通量 ,对流项 DEFINE_UDS_FLUX ,xj 2,, ,,扩散率 ,扩散项 DEFINE_DIFFUSIVITY 2,xj ,,,,,,,,,假定需要求解方程为，,,2,,,,t,x,xjj,, ,,,,,,,当，即如果计算达到稳定时，，此方程解析解为，,,2,,x，y，ct,,,,,xx,,jj,, 其中c为常数(假定其为0)，可由边界条件确定。因此，当解稳定时，，等值,,x，y线为一系列45度的斜直线。将此方程与标准方程对比发现: 名称表达式各项在UDS中需要定义的量需要在UDS中定义的宏离散形式的非稳态项 ,,非稳态,V,V DEFINE_UDS_UNSTEADY nn,1,，,, ,t项 ,t,t ,,,通量 ,,I,A,, 对流项 DEFINE_UDS_FLUX ,,xj(单位向量) ,,I,1,1,1 2,, ,,,1扩散率扩散项 DEFINE_DIFFUSIVITY 2,xj 源项 2 2 DEFINE_SOURCE 边界条,,x，y (b代表边界处) c,0 DEFINE_PROFILE bbb件 *注:源项和扩散率为常数，可以不用定义宏，但为了说明问题，仍采用宏来定义。假定如图所示的计算域，其中正方形边长为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); // ,,x，ybbb} 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%的误差。

                    本文档为【[整理版]博研联盟_一个简单而完整的uds例子】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
                    图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。 
 该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。

                    [版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。

                    本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。

                    网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
                

下载需要：免费已有0 人下载

立即下载

[整理版]博研联盟_一个简单而完整的uds例子

你可能还喜欢