esi-cfd版块系列教程之二
本帖最后由 nuaawubin 于 2014-12-27 17:23 编辑七、RAE2822跨音速翼型验证
借着第二届流体中文网数值模拟竞赛的东风，将第7个教程做完先。
问题描述：主要用于验证CFD软件对激波的捕捉，激波出现在翼型上面的中间位置，几何模型见附件（是GEOM可读的tab格式）。
计算状态：Ma=0.79，Re=6.5e6，攻角为2.31°。
使用GEOM划分网格，网格总数约为2.5w。
导入网格到FASTRAN中，如图所示。
在PT选项卡下，flow是默认勾选的。
在MO选项卡下，选择flow，选择理想气体，选择求解N-S方程。粘性系数改选为Sutherland公式。湍流模型选择k-eps模型。
在VC（流域设置）选项卡下，全部为fluid域，无需任何改动。
在BC选项卡下设置边界条件。
1、左侧为inlet，上下两侧为inflow/outflow，具体值均按下表给定。右侧为outlet，压力设定为108190pa。翼型上下面为无滑移物面。
Ma 0.79
Re 6.5e6
P 108.19KPa
T 255.56K
α 2.31°
U 233.6m/s
ux 233.42
uy 9.416
k 0.1m^2/s^2
ε 15m^2/s^3
在IC选项卡下初始化流域，与上表保持一致即可。
SC求解设置如图所示，其余保持不变即可。
提交求解计算，处理结果。
马赫数云图：
压力分布云图和压力分布系数如图所示：（激波位置还是不太准，可能与网格和湍流模型有关）
之前也用fluent（SA湍流模型）算了一下，比这次算的要好一些（http://cfluid.com/thread-145187-2-1.html），有兴趣的话大家也可以用Fastran的SA模型算算看。
感谢看帖，如有错漏，请指教。
perfact
楼主，你好！为想请教一下你湍流模型k，e系数是怎么取得？有什么说法吗？
yangjun2851247 发表于 2014-12-27 20:06
楼主，你好！为想请教一下你湍流模型k，e系数是怎么取得？有什么说法吗？
user mannual里有一个估算公式，你可以在湍流模型那一节找到
推送到首页“精华帖子”了。
nuaawubin 发表于 2014-12-27 20:44
user mannual里有一个估算公式，你可以在湍流模型那一节找到
按照帮助文档给的公式算出的k值和你的不一样？？
yangjun2851247 发表于 2014-12-27 23:20
按照帮助文档给的公式算出的k值和你的不一样？？
这个是我的错误，感谢指出。
nuaawubin 发表于 2014-12-29 17:11
这个是我的错误，感谢指出。
:)
狂赞，棒棒的！！
开卷有益，谢谢楼主！
页:
[1]