fluent计算三维机翼不收敛问题

IRESS

机翼弦长230mm,展长800mm，Re=20万，来流10m/s,用过k-e模型和S-A,可升力阻力都不收敛。是计算域太小了吗？还是设置有问题，麻烦大家看看，谢谢！

[ 本帖最后由 IRESS 于 2011-3-25 11:07 编辑 ]

xiangchuang

好像一般远场的距离至少为模型弦长的20倍

IRESS

原帖由 xiangchuang 于 2011-3-25 11:57 发表

                               
登录/注册后可看大图

好像一般远场的距离至少为模型弦长的20倍

指的是x方向的，对吧？
那高度方向也是厚度的20倍，宽度是展长的20倍吗？
这样网格量就要很大了

xrs333

不妨按风洞实验条件进行数值计算，也好与实验值对照，验证模型。

IRESS

原帖由 xrs333 于 2011-3-25 13:35 发表

                               
登录/注册后可看大图

不妨按风洞实验条件进行数值计算，也好与实验值对照，验证模型。

那流场的边界就更小了，现在遇到的问题是不知道所用的湍流模型，边界条件对不对。之前用同样的模型算过2维的（用的是压力远场边界），可以得到收敛的结果，到3维时，结果不收敛。是边界的问题吗？fluent中没有自由滑移壁，除了一个速度入口，一个压力出口，一个对称面，其他3个面用对称壁合理吗？

sliu2011

您的这个要算一个低雷诺模拟的的了。上下都用symmetry，应该不合适。如果气流有迎角的话，那不就不对了吗？流动都不能穿过去。你二维的既然有算过，应该知道用velocity inlet就是可以的。
算低雷诺数的问题，以及你选用的湍流模型，你需要检查一下Yplus值，如果这个值和模型设计的时候的值对不上的话，就是不当的使用。SA的话，把网格做到3以下，我一般都做到1一下。kepsilon的话，使用了wall function，应该是30和100之间。
因此你的网格要加密。从图上看，正交性也不太好。你用ICEM做的吧，再改改。

边界条件要对。特别是湍流那项。低雷诺数的仿真对于freestream turbulence intensity敏感的，湍流模型也是敏感的。
保证速度场的模拟是对的，检查一下velocity contour

onebyte

"在流场内的流动及边界形状具有镜像对称性时，可以在计算中设定使用对称边界条件"，你的流场整体不具有对称性，为什么使用对称边界？

stonesss

问一下 sliu2011 你说的3以下或是1以下是说的Y+值吗？算一个三维导弹的升力系数，工况2.5Ma，迎角4， 8等与试验数据差50%，随着计算步数的增加，误差越来越大，残差始终在10e-2处震荡，我用的Y+值为10-70，若进一步减少到1以下，边界层第一层的厚度会仅为5*10e-4，导致边界层内网格被压的很薄，同时长宽又不能调太低，那样的网格数量急剧增加，600万的网格计算时间已经不能接受了，在多的计算机内存就不接受了，还有用外o—BLOCK建的边界层，参考网格划分技术实例详解中设置为21层，可以设置的更多吗?另外看到一些文章里要求10-15层，边界层的第一层高度很重要，层数和整个边界层厚度要求严格吗？还是只要密度合适，整体厚度对计算误差大吗？我看大多数文章强调第一层的厚度，有些文章也对层数和厚度给出了计算方法，现在不知道层数和厚度如何确定，求大虾们指教