(完整版)naca0012翼型俯仰振荡实例进行讲解动网格.pdfVIP

以 NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程； 首先需要声明的是，这个例子也是来源于网络，原作者不详，在此向他表示感谢。 1、问题描述 本例是想对作简谐振荡运动的 NACA0012翼型的气动特性 （升力系数，阻力系数和力矩系数） 进行数值计算，来流速度为 V，攻角的变化规律为： Alpha （t ）=A/2*sin(omega*t) 其中， A=10 度， omega=10*pi 弧度 / 秒。 2、该例需要使用动网格来实现，首先需要编写刚体运动 UDF实现翼型的俯仰运动，由于在 FLUENT的 UDF中只能指定速度，角速度；所以，需要将攻角对时间求导，得到转动角速度 的规律： D （alpha ）/dt=A*omega/2*cos(omega*t) 编写的 UDF在附件中。 3、由于本例只是为了讲述动网格的实现，至于其他方面的设置及分析就不再讨论；这里详 细讲述下动网格的建立以及动网格的预览的结果。步骤如下： 1）将 mesh 文件读入到 FLUENT中，Grid ：check ，scale …，Smooth/Swap…；Display Grid ； 2 ）定义求解器为， Define ：pressure-based ，2D，unsteady ，Implicit ，Green-Gauss Node Based （因本例使用的是三角形单元） . 3）编译 UDF，Define-User-Defined-Functions-Complied … 此时打开了 Complied UDFs的窗口， Add…在选择 UDF的对话框中找到 NACA0012DM文件夹中 的 airfoil.c 文件，选中， ok ；此时返回到 Complied UDFs的窗口点击 Build ，FLUENT开始 进行编译，可以在 FLUENT窗口看到编译的一些过程提示；等编译完成，点击 Load ；就将已 经编译好的 UDF加载到 FLUENT中了。 4 ）定义动网格参数， Define-Dynamic Mesh （选勾，激活动网格模型） -Parameters … 此时打开了 Dynamic Mesh Parameters 窗口，在 Models 中只选取 Dynamic Mesh，本例的网 格类型为三角形单元，要实现的运动为小幅度的转动，因此选用的动网格更新方法为 Smoothing+Remeshing ；开始依次对这两种更新方法进行参数设定： Smoothing 中的参数设定： Spring Constant Factor （弹簧倔强系数），该值设定为一个较小的值，在 0.01 到 0.1 之 间，本例选取 0.08 ； Boundary Node Relaxation （边界节点松弛），设定为 0.5 ； Convergence Tolerance （收敛判据），保持默认的 0.001 ； Number of Iterations （迭代次数），保持默认的 20 ； Remeshing 中的参数设定： 为了得到较好的网格更新，本例在使用局部网格重新划分方法时，使用尺寸函数，也就是 Remeshing+Must Improve Skewness+Size Function 的策略。 将 Minimum Length Scale 及 Maximum Length Scale 均设置为 0 ，为了使所有的区域都被标 记重新划分； Maximum Cell Skewness （最大单元畸变），参考 Mesh Scale Info …中的参考值 0.51 ，将 其设定为 0.4 ，以保证更新后的单元质量； Size Remesh Interval （依照尺寸标准重新划分的间隔），将这个值设定为 1，在 FLUENT， 不满足最大网格畸变的网格在每个时间步都会被标记， 而后重新划分， 而不满足最小， 最大 及尺寸函数的网格，只有在 Current Time= （