[记录贴]参与第二届活动的记录

沈沉

看到白文老师已经放出交流数据，体现了本次活动的开放性。
我也开贴记录参与活动的过程，兼作提醒激励自己用。
上周针对Trap Wing 调试一些参数和计算方法。把Fine网格丢进服务器里跑着了。

coarse网格出了一些结果。放张13度攻角的各个站位的压力系数分布上来。

白文

不知是用的哪套网格？
是否也给出力和力矩系数以便于比对？
打算跟一些我们的计算数据。

沈沉

用的是29基地ICEM网格，没放出系数数据主要是TRAP-wing机翼面积不知道多少。
粗测了一下大约2.07平方米。不敢妄下定论，请白老师指教了。

PS：三阶Quick格式离散 + SST K-Omega 湍流模型

[ 本帖最后由 沈沉 于 2012-7-30 10:44 编辑 ]

白文

参考面积（半模）：22.028ft^2=2046468.16512mm^2
参考长度（平均气动弦长c）：39.634in=39.634*25.4mm=1006.7036mm
半展长：b=85.054in=2160.3716mm

网站上有哦：http://hiliftpw.larc.nasa.gov/Workshop1/geometries.html

此外，那套网格不好算，可能有点问题。

沈沉

----------------------        Mach         Alpha        CL                    CD
Trap-W  Coarse         0.2                13             2.063459                0.3300171  
----------------------        0.2                28             2.912983                0.6688254

放上数据，谢白文老师帮忙

这网格的确不太好算，主要是网格质量不太好，收敛很慢。当时选这套网格主要是看上去觉得比较漂亮。

[ 本帖最后由 沈沉 于 2012-8-1 01:10 编辑 ]

郭承鹏

沈兄，请问，您计算这个花费了多少cpu时！

白文

NASA-TrapWing模型的力矩参考点为：
(34.342in, -0.95in, 0.0in) = (872.2868mm, -24.13mm, 0.0mm)
备注：对应于原始模型，-Y轴为展向。

对于那个TrapWing-Config1-StructureGrid-CARDC-ICEMCFD网格，力矩参考点应为：
(34.342in, 0.0in, 0.0in) = (872.2868mm, 0.0mm, 0.0mm)
备注：这套网格的+Z轴为展向。

[ 本帖最后由 白文 于 2012-8-1 11:04 编辑 ]

白文

# Trap Wing curve fit data for config 1 - tests 478, 506, 513
# For CL, CD, or CM as a function of alpha
VARIABLES = "alpha","CLlower","CLupper","Cmlower","Cmupper","CDlower","CDupper"
ZONE T = "config 1 pre-stall experimental curve fit"
    -4.000       0.255       0.365      -0.252      -0.206       0.108       0.118  
    -3.000       0.387       0.497      -0.290      -0.244       0.101       0.112  
    -2.000       0.497       0.606      -0.316      -0.270       0.100       0.111  
    -1.000       0.633       0.742      -0.353      -0.307       0.105       0.115  
     0.000       0.809       0.918      -0.405      -0.360       0.112       0.123  
     1.000       1.064       1.123      -0.470      -0.440       0.123       0.134  
     2.000       1.160       1.217      -0.479      -0.451       0.136       0.147  
     3.000       1.250       1.306      -0.488      -0.460       0.151       0.161  
     4.000       1.336       1.392      -0.495      -0.468       0.166       0.176  
     5.000       1.418       1.474      -0.501      -0.474       0.182       0.193  
     6.000       1.498       1.554      -0.506      -0.478       0.199       0.209  
     7.000       1.576       1.632      -0.509      -0.482       0.215       0.226  
     8.000       1.652       1.708      -0.511      -0.484       0.233       0.243  
     9.000       1.726       1.782      -0.513      -0.485       0.250       0.261  
   10.000       1.799       1.856      -0.513      -0.486       0.266       0.282  
   11.000       1.872       1.928      -0.514      -0.486       0.285       0.301  
   12.000       1.943       1.999      -0.514      -0.486       0.305       0.321  
   13.000       2.014       2.070      -0.514      -0.487       0.327       0.343  
   14.000       2.083       2.140      -0.514      -0.487       0.349       0.366  
   15.000       2.152       2.208      -0.515      -0.487       0.369       0.393  
   16.000       2.220       2.276      -0.515      -0.487       0.391       0.415  
   17.000       2.286       2.343      -0.515      -0.487       0.414       0.438  
   18.000       2.352       2.408      -0.515      -0.487       0.438       0.462  
   19.000       2.416       2.472      -0.514      -0.487       0.462       0.486  
   20.000       2.478       2.534      -0.513      -0.486       0.482       0.518  
   21.000       2.538       2.595      -0.512      -0.484       0.506       0.542  
   22.000       2.597       2.653      -0.509      -0.482       0.530       0.566  
   23.000       2.653       2.709      -0.506      -0.478       0.554       0.590  
   24.000       2.707       2.763      -0.502      -0.474       0.578       0.613  
   25.000       2.757       2.813      -0.497      -0.469       0.595       0.644  
   26.000       2.804       2.860      -0.491      -0.463       0.620       0.668  
   27.000       2.847       2.904      -0.484      -0.456       0.643       0.691  
   28.000       2.886       2.942      -0.476      -0.448       0.664       0.713  
   29.000       2.919       2.975      -0.467      -0.439       0.685       0.733  
   30.000       2.946       3.003      -0.456      -0.428       0.701       0.751  
   31.000       2.967       3.023      -0.443      -0.415       0.722       0.771  
   32.000       2.979       3.035      -0.427      -0.399       0.738       0.788  
   33.000       2.982       3.038      -0.407      -0.380       0.751       0.800  
   34.000       2.974       3.031      -0.382      -0.354       0.761       0.810  
   35.000       2.919       3.003      -0.348      -0.320       0.772       0.822  
   36.000       2.868       2.954      -0.304      -0.276       0.792       0.841  
ZONE T = "config 1 post-stall experimental curve fit"
   35.000       2.293       2.395      -0.254      -0.215       0.772       0.822  
   36.000       2.221       2.321      -0.236      -0.197       0.792       0.841  
   37.000       2.148       2.248      -0.219      -0.180       0.828       0.878  
#
# Date file created by J. Hannon:   04-Jun-2010

备注：
1、数据不再是单点的，给出了上下界，也就是说给出了误差带。
2、风洞实验数据修正后，小数点后只给出了3位有效数字，包括阻力系数。
3、感觉人家做事还是比较靠谱的。

[ 本帖最后由 白文 于 2012-8-2 16:04 编辑 ]

白文

原帖由 沈沉 于 2012-8-1 09:06 发表

                               
登录/注册后可看大图

----------------------        Mach         Alpha        CL                    CD
Trap-W  Coarse         0.2                13             2.063459                0.3300171  
----------------------        0.2                28             2.912983                0.6688254

放上数据，谢白文老师帮忙

这网格的确 ...

整理一下。。。
                                Alpha            CL                        CD                      Cm
SKLA-Str-CG-沈沉   13           2.063459             0.3300171            -0.4860940
SKLA-Str-CG-白文   13           2.016603             0.3258261            -0.4632183
Exp. by J. Hannon  13.000      2.042(+-0.028)    0.335(+-0.008)     -0.500(+-0.014)

SKLA-Str-CG-沈沉   28           2.912983             0.6688254            -0.4317842
SKLA-Str-CG-白文   28           2.839007             0.6533513            -0.3916434
Exp. by J. Hannon   28.000     2.914(+-0.028)    0.688(+-0.025)    -0.462(+-0.014)

备注：
1、其中13度攻角为50000步（5万步）计算结果，28度攻角为100000步（10万步）计算结果，感觉还未完全收敛。
2、其中SKLA是空气动力学国家重点实验室的缩写，CG代表Coarse Grid。

[ 本帖最后由 白文 于 2012-8-2 15:21 编辑 ]

白文

升力曲线、阻力曲线和力矩曲线

[ 本帖最后由 白文 于 2012-8-1 14:30 编辑 ]

沈沉

原帖由 郭承鹏 于 2012-8-1 02:04 发表

                               
登录/注册后可看大图

沈兄，请问，您计算这个花费了多少cpu时！

郭兄，小弟当时没注意看，现在回想估计大概每个算例750CPU小时（以单线程1.0GHz计）（约460万网格，每算例迭代3万次）

[ 本帖最后由 沈沉 于 2012-8-1 11:56 编辑 ]

沈沉

白老师小结了一下，有实验数据，我也对自己的计算有信心了
（各广义气动力计算结果都落在实验数据范围内 ）

1、我和白老师计算结果差异主要与采用不同的离散格式与湍流模型有关。
2、高阶精度离散算得的气动力数据更大也是容易理解的，估计主要差在粘性部分。
3、大攻角时S-A可能差得多些，这与流动分离有关。
4、收敛性确实不好。

这里补上力矩系数
alpha 13           Cm   0.4860940
alpha 28           Cm   0.4317842

实验数据很规范，比起有些风洞的结果要好多了 。
Medium网格的计算结果(算例2a)估计要过段时间了。

[ 本帖最后由 沈沉 于 2012-8-2 06:36 编辑 ]

白文

这个TrapWing-Config1-StructureGrid-CARDC-ICEMCFD网格（准确称号是TrapWing-Config1-StructureGrid-SKLA-ICEMCFD网格，以后简称SKLA网格）是套结构网格，将其转换为非结构网格，采用AVICFD-X软件进行了计算，多重网格不收敛，只得死算。现将收敛情况报告如下，谨供参考和交流讨论。

                                     Alpha          CL                        CD                      Cm
SKLA-CG-050000steps   13       2.016603             0.3258261            -0.4632183
SKLA-CG-100000steps   13       2.032511             0.3296754            -0.4708202   
SKLA-CG-110000steps   13       2.033000             0.3297700            -0.4710500
SKLA-CG-120000steps   13       2.033300             0.3298100            -0.4712700
SKLA-CG-130000steps   13       2.033500             0.3298300            -0.4714500
SKLA-CG-140000steps   13       -----------              -------------             --------------
SKLA-CG-150000steps   13       -----------              -------------             --------------   
Exp. by J. Hannon            13       2.042                   0.335                    -0.500
                                               (+-0.028)             (+-0.008)                (+-0.014)

SKLA-CG-015000steps   28      2.659208            0.6142657            -0.3204359
SKLA-CG-050000steps   28      2.775500            0.6316400            -0.3605500
SKLA-CG-100000steps   28      2.839007            0.6533513            -0.3916434
SKLA-CG-150000steps   28      2.854905            0.6585423            -0.4012937
SKLA-CG-160000steps   28      2.855700            0.6588100            -0.4018900
SKLA-CG-170000steps   28      2.856300            0.6589900            -0.4024000
SKLA-CG-180000steps   28      2.856600            0.6590700            -0.4027300
SKLA-CG-190000steps   28      ------------            -------------             --------------
SKLA-CG-200000steps   28      ------------            -------------             --------------
Exp. by J. Hannon            28      2.914                  0.688                    -0.462
                                              (+-0.028)            (+-0.025)                (+-0.014)
SKLA-MG-015000steps   28      2.677437           0.6114225            -0.3214453
SKLA-MG-050000steps   28      2.791744           0.6353854            -0.3639362
SKLA-MG-088000steps   28      2.838800           0.6525700            -0.3875200
SKLA-MG-100000steps   28      2.845100           0.6538044            -0.3906060
SKLA-MG-110000steps   28      2.851800           0.6562500            -0.3940800
SKLA-MG-120000steps   28      ------------           -------------             --------------
SKLA-MG-130000steps   28      ------------           -------------             --------------
SKLA-MG-140000steps   28      ------------           -------------             --------------
SKLA-MG-150000steps   28      ------------           -------------             --------------

粗网格（CG）的网格点数为：4885867，即大约489万。
计算资源：4节点*2CPU*4核（Intel Xeon E5430 2.66GHz），每节点32GB内存。
粗网格每时间步大约3秒，迭代50000步共用时大约41.5小时。

中等密度网格（MG）的网格点数为：12840777，即大约1280万。
计算资源：128核。
中等密度网格每时间步大约4秒。

先到此为止，为了引起注意 ，另开帖讨论了。

[ 本帖最后由 白文 于 2012-8-10 15:48 编辑 ]

周华

GAMBIT还是无法导入CGNS，可能是header标识有问题，不知道GAMBIT的版本是不是太老了。今天也用ICEMCFD转换了一个非结构的粗网格例子，先试算一下，开会回来再逐个尝试。

周华

原帖由 白文 于 2012-8-1 12:18 发表

                               
登录/注册后可看大图

整理一下。。。
                                Alpha            CL                        CD                      Cm
SKLA-Str-CG-沈沉   13           2.063459             0.3300171            -0. ...

Hannon的试验条件是怎么取的，没有特殊的流动条件吧？