CFD中文网

@韬智tz 谢谢

@buaa_hgj 我也在做柔性变形，前辈您有什么经验可以指点一二吗？

@jinlinna 感谢！

如果我没记错的话，流体中文网上的ICEM版块里有人画过4119桨的结构网格，也传了文件在上面，你可以去上面找一下看还在不在。（耐心往后翻一下）

@tens 明白了，十分感谢！

@硫酸亚铜
看来这几何给你之前刚从土里挖出来

@cfd-china 直接用WIN自带的远程连接就可以了

SA是Spalart-Allmaras，Reynolds-Averaged Navier-Stokes Spalart-Allmaras turbulence和Spalart-Allmaras Detached Eddy
Spalart-Allmaras采用的是标准的SA，不是高雷诺数的。DES也是标准的DES模型

@cfd-china 好的

休斯顿市区内150平米loft卖15万美元：https://www.zillow.com/homes/for_sale/Houston-TX/pmf,pf_pt/39051_rid/globalrelevanceex_sort/29.893043,-95.053368,29.656432,-95.539513_rect/10_zm/

@一朵苍凉 我也不太懂，可以看看控制领域的相频特性，也许有帮助:zoule:

fluent自定义过饱和函数时应该用那个压力数据来计算？我的静压只有不到100pa，太小了

@两月三年 good ，好主意。我试试哈。

glgoo，我一直在用的

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？ 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
3.png

@李东岳 多谢李老师回复。上传的结果是关闭Decompose polyhedra选项后的结果，Tecplot里面也是这种情况。用的是2D的问题，采用的是选择x、y方向加密的方式。

在源码里看到refineMesh调用了cutDirToEdge，其中加密方向是根据单元边与选择方向的关系判断加密的方向（下图的cosAngle，maxEdgeI）。O型网格按照这样处理的方式，处于圆柱上方的网格和后方的网格加密的方向是不一致的（第二张图片是只加密x方向的结果，网格的加密方式不一样了）。不知道是否可以通过设置参数进行控制？

845840da-6fb4-4df4-8d1d-59871760b046-image.png

3a078fc9-11f2-4e2c-bc86-9baa876c65fd-image.png

Library Genesis : http://gen.lib.rus.ec/
这上面恐怕能搜到大部分上述电子书

针对一个高度为3mm，宽度为1mm的矩形微通道，用fluent mesh 画网格后，在入口流速为0.3 0.4 m/s时，计算都可以收敛。但是当入口流速为0.5m/s时，计算收敛不了，这可能是什么原因导致的，尝试过增加网格数量也无效

好的，我看看，文章不长，图很多，简单点。。。