CFD中文网

看到了 GitHub 给每个用户生成了基于 2020 提交历史的 stl 文件，然后就试了试自己的。

我的阻力系数是 4.633370e-01。

计算 case 在这里：https://github.com/randomwangran/randomwangran/tree/main/my-skyline

你的阻力系数是多少？

跑了一个虹吸的模拟，挺有意思的，能模拟一些生活中的流体现象。

https://github.com/randomwangran/capillary-will-rise

以把重力翻过来的话，就和楼主的模拟比较接近了。

Sorry I had hard time to type Chinese.

Thanks for sharing this. Your finding is consistent with what I've
read this in Wilcox's book (P98).

In free shear flow, k-epsilon model is unaffected by the numerical
value of $k$ and $\epsilon$, but k-omega is very sensitive to the
choice of $\omega$.

But the book does not provide any suggestions on how to achieve mesh
independence on k-omega model.

Speaking to mesh independence study. The idea is simple: refining the
mesh at the region with high gradient. In reality, it depends on
researchers to define where are these regions. Even in the simple
context like flow over a circular cylinder. Different people choose
different types of the computational domain, different strategies to
generate the mesh, different parameters to build the block.

There are something I would investigate:

1. continue refine the mesh in the "slot exit"
2. expand the computational domain
3. expand the "slot exit" region, especially where the shear layers
   are.
4. use unsteady solver
5. try LES if you have computational resource

Wilcox, D. C., & others, (1998). Turbulence modeling for cfd. : DCW
industries La Canada, CA.

@波流力 非常感谢,我放在了这个仓库中,会做进一步步的研究,有值得分享的我会回来盖楼的.

求 St 的方法有很多，比如对阻力或者升力系数 fft，或者在wake区域对某个参考点速度检测，然后在垂直于来流方向做fft。需要注意的是对Cd做fft所得的St应该是对Cl的两倍才对.

ref: https://cfd-china.com/topic/774/圆柱扰流计算与物理现象问题思考/2

@hoversoar

还在,:).

时间离散对结果影响大吗，还有Crank-Nicolson后面的系数应该怎么取？

Please see: https://cfd-china.com/topic/3834/分享-openfoam不同离散格式的简单研究

感谢楼主分享,这些数据会让大家少走很多弯路.

我有一个想法,不知楼主能否分享这个工做中 OpenFOAM 的一个算例文件?
后来的人可以基于楼主研究,进一步研究相关问题.比如,在网格质量比较差情的况下,会有怎么样的影响.或者计算效率之类的.

可以先用谷歌代替：

site:https://cfd-china.com/ pisoFoam

深感这三年的付出如同泡沫一般付诸东流，自己无论怎么折腾都只是在瞎蹦跶。

可是，写过的每一行代码，读过的每一本书，推过的每个一个公式，都是你的啊。

国外的世界诚然宽阔，可是不及你心中的世界。

胡思乱想的时候，跑跑步吧，刷 500 圈 400 米跑道；游 500圈 25 米的泳池；骑 500 公里的山路，之后，

我发现，哪怕平平淡淡的生活，都是无比幸福的。

@hoversoar

上次算到step2的中心差分的50.5s处，说可能会阻力又下降。

建议楼主定一个准则.比如计算了 N 个 Vortex Shedding 周.之后又计算了 2N 个周期,然后比较一下所感兴趣的物理量.

那之后我把它继续算到了59.5s，然后速度对流项的linear换到了
limitedLinear。把结果贴上来一下。

不清楚楼主研究的重点是?层主个人来说,是很好奇, linear -> limitedLinear 的区别.

整体看起来的话，三个阶段的定量CdCl差距不是很大。和文献也对的比较上。
初 步认定是速度对流项选择了迎风格式造成的。跟之前的猜想有关，在55s处残
差的Uy震荡明显，造成了Cd上升，Cl振幅也增大. 换成limitedLinear 0.5 之后，
可以发现残差有明显的下降。

不知楼主有没有对这些格式的比较?

UDS 对网格要求更高.

"Peaks or rapid variations in the variables will be smeared out and,
because the rate of error reduction is only first order,
very fine grids are required to obtain accurate solutions."

楼主还可以探索一下这些:

• Linear Interpolation (CDS)
• Quadratic Upwind Interpolation (QUICK)
• Higher-Order Schemes

@random_ran 是否可以回答前面的问题，虽然P的量级并没有发生变化。但是Uz
的曲线震荡明显（相比中心差分），猜想可能是由于TVD 里头那个0.5造成的，
混入了一部分迎风格式，造成耗散？

我没法准确回答你,因为我自己没有太多探索过.

我这样中途切换差分格式是否ok？

为什么不可以?如果是研究格式的问题,我可能也会这么做,省很多时间.但是如果是分析流体机理,最好从头开始用一种固定格式.

今天才把case文件整理了一下,下一步找合适的网格.

@东岳 岳兄过奖了.
@hoversoar 今天才整理了一下你的 case.

乃有新发现的话留言在这,大家都能看到,或者帮助你,我只能尽力,还是要靠氖自己的.

@hoversoar 台兄过奖了。我只是好奇罢了。自己玩的很杂，可是遇到能一起玩的，实在不易，所以慢慢来。我在找 of 的脚本，自动生成。最好简单可控，小 mesh 大家可以互相交流，大 mesh 大家自己留着做研究就好。

楼主 mesh 太大了,我把设置文件拷贝下来,打算用很小的脚本来生成圆柱网格.

这样的话,别人可以更方便地帮助你.

我们就可以一起讨论了.

仓库地址: https://github.com/randomwangran/chigua

Uz太抢眼了。

不知道坐标系是如何定义？

我看了下之前算的几个例子，发现p的残差几乎都在0.001这个量级。

非常有意思的问题，我之前没怎么考虑过 压力 残差的取值问题。我也好好思考一下。

不知道怎么才能继续降下去？

减小时间步长？提高收敛准则？或者优化网格形态？...

我一开始尝试增加到10个nOuterCorrectors, 发现第4~5个外循环开始残差下降不明显甚至开始上升，No Iteration也几乎为0。 为了节省时间所以设定到5。

我之前也一直在用 piso， 不过是默认值，没有动，我得好好看看 piso 算法再来回答你的问题。

现在这边的例子是光滑圆柱，如果要算看看很高兴分享。（how）

• 我记得论坛里有人分享过 *.zip 包含算例和运行脚本，简单点的办法就是做相同的事情

• 稍微难一点的是可以把楼主的 case 文件，在层主的 compute 仓库下 caseIdoNotUnderstand

然后我们一起来 hack.

前排小板凳，

附上现在1w雷诺数的光滑圆柱绕流的例子的某个时间步的截图，用of里头的monitor监控的residual，我也一直很困惑Uz的曲线不像稳定的样子。

为什么关注点再 U_z 不在 残差最大的 p 上投入精力？

我现在把2个nOuterCorrectors由2提高到5（从31秒开始），算个十秒左右看看结果。

这样做的原因是？ 为什么是从2到5？而不是其他的数字，比如7?

目前的时间步库朗数是控制在1.5以下。

这个十秒计算过之后，我会把速度对流格式换成中心差分继续跑看看结果是否有改善，到时候再过来更新！

这个是 case 是机密么？如果能简单运行，我也想试试你的 case.

Re 估算是 10，不知道高雷诺数效果怎么样。

一篇文章，一直对 Deep Learning 很好奇，这篇文章讨论了一些实现办法：https://arxiv.org/pdf/2002.10558.pdf

欢迎大家交流。

The dynamic viscosity is µ = 5×10−3kg/(m · s) and the density is ρ = 1kg/m

U_max equals to 0.5 m

Diameter = 0.1 m

看到有人讨论圆柱绕流就很激动。

可是我已经挖了太多坑还没填。关注中，加油。

@我是河滩 在 应力张量如何理解？ 中说：

为什么应力张量有9个分量，张量是对一点来说的还是面？对流体微元或控制体而言，不同面上的张量一样吗？

这种问题,或许会让人不开心.我很感兴趣,记录在了我的笔记本上,愿意与你讨论: https://github.com/randomwangran/computing/blob/master/m/talk.org

A video to show flow around circular cylinder in real life:
https://youtu.be/0ThQ_nD97hY?t=306