STARCD学习经验.doc

1、 将列车的六个表面分别设置成wall和与将列车的六个表面设为一个wall所得到的力和力矩的对比2、 match edges功能3、 通过erd得到的力矩数值的来源4、 Icemcfd经验5、 SOLUTION DIVERGES/UNSTABLE6、 TOLERANCE CONSTANT SITOL1 (area tol) FOR MATCHING ASI SET TO 1.00000Dimension 出错。7、 能不能在一个位置同时创建couple和attach8、 在正常无法打开瞬态pstt文件的情况下，如何打开？9、 删除block出现错误时的解决办法10、 出现couples错误的解决办法11、 Symmetry Plane的知识12、 稳态的模型和瞬态的初始模型必须完全一致13、 进行瞬态计算时首先需要考虑到的问题14、 增大TE的松弛因子对于力和力矩的影响15、 瞬态计算中力的波动过大的原因16、 规则列车网格模型的测试总结17、 总结18、 当starcd的设置与prep的设置不相符19、 将icemcfd中划分的网格导入到gambit20、 稳态残差和瞬态残差大小的对应关系21、 瞬态计算的info中某一步出错计算所需时间22、 网格间距的增量23、 “Biexponential”网格划分方式24、 无横风和有横风的对比25、 计算区域的后部长度和上风侧的宽度对瞬态计算横向力波形图的影响26、 瞬态计算时的cells数量问题27、 创建cylindrical坐标系及相关操作28、 列车过曲线中STARCD的VCEX解决的办法29、 STARCD中cylindrical坐标系统操作时的注意事项30、 创建couple出现错误和曲线瞬态计算时出现错误的现象及原因31、 在icemcfd中对模型进行旋转的办法32、 出现“*ERROR - MASTER CELL 1004699 IN COUPLE 1219 IS UNDEFINED”问题的解决办法1、 将列车的六个表面分别设置成wall所得到的力和力矩，在x、y、z方向加起来的总和与将列车的六个表面设为一个wall基本一致。2、 利用pre-mesh paramsmatch edges功能可以将列车斜面的最上部和最下部的edge边界层的间距设置为一致。Reference edge选择最上部的edge，method选择“copy to edge”，“copy absolute”打勾，copy to edge选择最下部的edge，就可以了。3、 力矩的计算是针对每个cell face进行的，而不是将力集中到某一点上来进行计算（假设力的作用点在每个cell face的中心上）。4、 有几次突然打开时在model tree中找不到blocking或geometry模块了，重启了一下电脑就又可以了。加入了挡风墙之后网格的划分变得更加复杂，位于挡风墙和列车之间区域的网格划分需要用不同的原始block来完成。Icemcfd11中创建的模型在icemcfd5.1中打开会出现错误。5、 计算存在错误，错误为SOLUTION DIVERGES/UNSTABLE，Diss出现了NaN的数据。原因在于滑移交界面两侧的网格差异太大。尤其需要注意的是：不光是计算开始时滑移交界面两侧的网格差异不能太大，在计算进行的过程中滑移交界面两侧的网格差异也不能太大。上述错误的原因在于中间部分单独向前移动了2层，导致模型顶部一部分交界处两侧差异太大。原来的模型划分的交界处两侧差异不大，但是移动了之后差异变大，导致计算发散。还有就是模型和原来稳态的也不完全一致了，所以导致无法从pst开始。中间部分的网格一直在向前移动中，而有一部分（中间部分）的网格是很密的，虽然在开始的时候将其交界面两侧的网格调整为差别不大，但是移动之后，交界面两侧的网格差异加大，因此而导致了这种情况。解决办法1、全部而不是局部网格添加/删除；2、加密左右部分的网格，使得在整个移动过程中滑移交界面两侧的网格差异都不是很大。针对第2种方案我考虑了一下，加密网格可能会导致加大不少的计算量，因此不建议采取。6、 默认的，STAR-CD允许an area mismatch为2%，这样的mismatches发生在part of a cell face on one side of the sliding interface cannot be matched to any face on the other side。如果face areas大于默认值就会导致unmatched，计算就会停止并且伴随着一个错误信息。当faces严重扭曲时或者master and slave faces的相对大小差异较大时这种情况可能会发生，或者存在移动网格的情况。matching tolerance可用通过下面的设置来放宽：RCON，12，。将其设置为2则上述的错误消除。7、 经过测试，在一个位置上创建了couple之后，这个位置无法再创建attach了，否则check everything的时候会在attach上出错。8、 、首先在不打开moving mesh的情况下“open post file”，直接载入，然后“add file to list”和“open transient post file”，然后再然后随便选择一个time step，“stor iter/time”，再打开moving mesh，“open post file”， 会弹出一个错误窗口，然后“add file to list”和“open transient post file”载入，这样就不会出错了。见F:work滑移交界面的测试-STARCD-修改6。、步骤和上述略有不同，首先在不打开moving mesh的情况下“open post file”，直接载入，然后“add file to list”和“open transient post file”，然后再然后随便选择一个time step，“stor iter/time”，然后在输出窗口中“ifile prep.inp”，再打开moving mesh，“open post file”， 会弹出一个错误窗口，然后“add file to list”和“open transient post file”载入，这样就不会出错了。见F:work滑移交界面的测试-STARCD5。9、 当删除block出现问题时，或者赋予block给part时（直接退出程序），可以使用blockinginit output blocks，就不会出错了。10、 在filewrite geometry file时出现很多关于couples的错误，随便列举一个如下ERROR - CELL 95811 FACE 1 OF NON-PARTIAL COUPLE 22423 IS 22.26 PERCENT UNCOVERED.修正错误的步骤、通过toolscouple toolcouple check，打开“build new couple set with failed couples”功能，然后利用cnewcouple setboth可以看到错误的couples位于哪些位置。、couple tooldelete couplescpset删除上述那些错误的couples。、然后edit types对couple进行编辑，打开partial boundaries功能。、然后利用create couple功能重新创建couple，注意选择“cell set”。此时再次使用couple check功能进行检查，发现刚才的错误没有了。目前稳态的模型是没有错误的，但是移动一定距离之后有关couple的错误又出现了。不知道是什么原因。经过检查，该模型在xz方向移动没有问题，但是在y方向移动就会出现问题。创建couple时所有的couple都打开partial boundaries功能则没有出现该问题。同时需要经常利用toolscheck toolcrack检查一下，如果打开了partial boundaries功能而crack还存在问题，可以通过增大tolerance解决。11、 很多参考书籍中提到，应该将计算区域的上表面和两侧面设置成对称面。根据流场的尺寸选取的条件, 除“进口”、“出口”外的这些面附近的流体运动已经趋于静止，均定义成对称面。Symmetry Plane Denotes a surface at which the normal velocity and normal gradients of all other variables are zero正向速度为0，所有变量的正向梯度为0。12、 利用稳态开始进行的瞬态计算，稳态的模型和瞬态的初始模型必须完全一致，而不是和瞬态删除掉初始须删除的网格之后的模型一致，否则就会出现下面的错误。* READ RESTART FIELD FROM FILE star.pst *ERROR - READING POST FILE:ND1 / ND2 Error P IN R9DP1 : DIMENSION ERROR13、l 做瞬态计算之前先测试一下稳态计算，如果稳态计算没有问题再进行瞬态的计算，省得做无用功。稳态的模型和瞬态的初始模型必须完全一致(仅仅是指模型的网格，模型的坐标位置可以不一致)，而不是和瞬态删除掉初始须删除的网格之后的模型一致。l 进行瞬态计算之前先检查一下滑移交界面两侧的网格差异，滑移交界面两侧的网格划分差异不要太大。l 如果涉及到全部添加/删除网格，那么应该首先考虑一下cgrid.cgrd、uparm、prep的设置问题。l 建立模型需要考虑到瞬态计算时多出去的那部分网格的处理，在z坐标为0的位置，左右滑移交界面两侧的网格必须在z方向上重合，这样才方便删除位于z坐标为负位置的网格。l 对cgrid.cgrd的调试可以直接通过egrid read cgrid.cgrd命令来进行，这样如果存在问题的话更容易看到问题的所在。l 在瞬态计算的过程中，应该随时监控的状况，如果出现错误（例如MAXIMUM COURANT NUMBER太大），那么应该及时终止。14、 将TE的松弛因子增大加速其收敛对于力和力矩的影响不大，力和力矩的差距小于5。15、l 几乎所有的先不可压缩再可压缩的方式都会在可压缩的第一步出现错误，戚振宕说这是由于不可压缩与可压缩的切换所造成的。l 提到无论是“稳态结果”还是“不可压缩-30步”或是“可压缩-370步”其最终结果都应该差不多。l 首先分析的应该是稳态结果是否符合常理，然后才是瞬态不可压的结果与稳态结果相近，最后再说瞬态可压的结果。l 说明：计算起始时列车尾部距离计算区域尾部为0也无所谓，因为前面的一段时间数据的发展还没有进入稳定期，这一段时间的数据是不准确的，基本不用。l 提到瞬态计算中速度和TE的残差也不能设大了，最好为10-4，由于压力速度耦合方程的原因。16、 有关移动网格中当starcd的设置与prep的设置不相符时，以prep的设置为准。17、 从icemcfd中导出：fileexport meshwrite stl file。得到一个mesh文件，采用原始名字，不要改名字，类型为证书信任列表，拷贝到gambit目录下。导入到gambit：fileimportmeshtype选择stl，即可以打开了。18、 如果瞬态计算的info中某一步出错，错误如下： * WARNING #042 * PROCESS NOT FULLY CONVERGED - EQ.,NIT,RESI: UMOM 1001 0.56488E-02 * WARNING #044 * SOLUTION RESCUED AT RESIDUAL LEVEL 5.573350E-03那么这一步计算所需的时间比没错的某一步计算所需的时间要多的多，大概是好几倍的关系。19、 网格间距的增量大点