fluent动网格算例设置.doc

前段时间曾有网友问过该cas的做法，拖的时间太长了，实在不好意思，今天上午就抽空做了一下（这可是工作时间，真对不住给我人民币的老板了!)闲话少说，直插正题：制作该cas需要用到以下软件：gambit（做网格的，网上有的下）fluent （计算用的，网上也有的下）fluent小工具（就是可以在编译udf不用装vc的，我的blog上有的下）计算边界定义说明：1.生成网格，本文采用的是非结构网格，相信这个过程就不需要我赘述了吧，如果不会的话，请参阅gambit生成网格的文献。2.将本文的udf附件（box.c）复制到你的工作文件夹（也就是你放mesh的文件夹）。启动fluent，读入网格。选择define/User-defined/s/compiled。如图：在source files中 add 你的udf（也就是box.c），点击build，然后点击load。3.在Define/models/solver选择非定常求解器unstead.4.在Material面板中选择fluent database,在左侧的菜单中选择water，点击copy。5.在Dynamic Mesh中点击parameters，激活Dynamic Mesh，选中soomthing和remesh,具体参数如下图所示：6.在dynamic mesh选项下选择zones，设置box为rigid body。设置wall为deforming7.在boundary conditons 将in设置为速度入口，速度为2m/s,out设置为outflow或者压力出口都可以。8.初始化。9.将时间步长设为0.05（或者更小），时间步设为200，开始计算。附件：udf：box.txtcas：feijiegou.rar网格变化过程的动画：sequence2.rar后记:由于论坛上传附件格式的限制，麻烦各位下载box.txt后，将后缀txt改为c即可食用。该udf直接copy自该网友的帖子，我觉得该udf有点问题的：它将box所受的压力看成是所受x方向上的合力。我认为最合适的还是直接将6dof的udf稍加改造即可。该教程只是一个简单的制作步骤，只是详细的描述了一下动网格部分的设置，如果各位网友需要其他常规操作的具体步骤描述（比如边界条件的设置，动画的设置，后处理等等），请参考fluent的tutorial等参考资料。box.txt#include udf.hstatic real v_prev = 0.0; DEFINE_CG_MOTION(moving_body, dt, vel, omega, time, dtime)Thread *t;face_t f;real NV_VEC (A);real force, dv;NV_S (vel, =, 0.0);NV_S (omega, =, 0.0);if (!Data_Valid_P () return;t = DT_THREAD (dt);force = 0.0;begin_f_loop (f, t) F_AREA (A, f, t); force = force + F_P (f, t) * A0;end_f_loop