CFD中文网

@东岳 老师，DEM这块不是我编写的，工作主要是在CFD部分，据我了解DEM是计算单个颗粒的曳力，不用指定气泡直径大小吧；那我再看一下CFD部分曳力计算吧，谢谢您

@夜阑烟寒 从昨天开始domecross用不了了，你那是不是也是？我用了得一年半了吧，第一次开始这么长时间用不了

@东岳 嗯嗯，确实是尺度问题，谢谢东岳老师！

谢谢老师，我找出来看一下

如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

二维的vorticity只有一个方向 实际上从涡的形成来看 二维似乎并不适合一些涡量的模拟
但是对于一些简单的流场 实际上2D具备非常大的优势 他的处理简单 很多算法很好自编程实现
并且网格可以完全mapped 最大的优势 还是在于后处理的一目了然
这就是为什么大多数理论研究还是集中2D 但是2D可能不是太适用Fluent Fluent主要接轨工程实际问题

@ibelief 请问是视频教程么，可否分享一下，最近在自学这本书，有视频学起来应该效果更好

不能的，目前只能二选一。

@两月三年 如果重装软件都不行，你可以试一下我说的方法，我这方法在其他两个同学那里生效过。

大佬解决这个问题了吗？

网格问题？液滴被设置成壁面了？

@肖艺 STAR CCM+
例子： https://www.cfd-schuck.de/en/cfd/anwendungsbeispiele/simulation-of-gear.html

看错了 你这个是icem吗

https://mirrors.tuna.tsinghua.edu.cn/#

下载速度很快

0_1543372305704_捕获.JPG

https://www.ansys.com/-/media/ansys/corporate/resourcelibrary/article/aa-v2-i2-coupling-momentum-and-continuity.pdf

刚看了下，Fluent那面也只是速度和压力耦合求解。能量抛在了外面。

但是另一个文档速度压力能量放在了一起 ： http://www.bakker.org/dartmouth06/engs150/05-solv.pdf

速度压力耦合求解在高超那面用的多，如果分析N-S方程的数学特征（双曲抛物），你需要把u, v, w, p当成矢向量，然后求特征值特征向量。接下来求解这个矩阵就是耦合算法。我目前只是了解到思路，没亲自植入过。可

大家好！
我想创建两个等腰直角三角形的Y型网格划分，2d不能直接创建Y-BLOCK，只查到了先O型划分再移动节点，间接建立Y Grid，但是这样就会影响我相邻块的关联，有什么办法能实现其划分吗？如下图，红色区域就是两个直角等腰三角形，其相邻的块都是矩形块，我用O型划分了一个，但是周围的块的关联已经很奇怪了。在此先谢过各位前辈了！

cfd2fd2a-2a3a-457e-90da-8343185b0c87-image.png

@blgy 非常的感谢，你这个是在哪看到的啊？你这个图片也看不出来用什么软件，用什么方法划分的啊，仅仅是欣赏一下:ok:

@SnoWolf2018 @玉玺

多谢！！

@东岳 李老师，其实我说的是一个T型管，已知主管道进口速度、出口压力和侧支管道的出口速度。鉴于某些目的，我必须采用可压求解器，主管道进口给质量流量，主管道出口给压力出口边界条件，侧支管道给速度出口，虽然采用了可压求解器，但是流速不超过20m/s，我认为在侧支管道出口依然可以按照不可压处理，采用速度边界条件是可行的

//2.4s label meshSize = this->thermo().rho()().size(); for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < meshSize; icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } } //24s for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < this->thermo().rho()().size(); icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } }

上边两种用法，计算耗时相差 10 倍！
改用 thermo.T() 之后计算时间又变成 2.4 s 了，所以原因应该是 rho 后边跟了两个括号，深层原因是啥呢？

//2.4s label meshSize = this->thermo().rho()().size(); for(label i=0; i < this->nSpecie_; i++) { for(label icell = 0; icell < this->thermo().T().size(); icell++) { label izone = ZoneNumber[icell]; this->RR_[i][icell] = RRCCM[i][izone]*correctCH[icell]; } }

此外，还可以使用 forAll(this->thermo().T(), icell)，效果和 for 循环几乎一样。