CFD中文网

@东岳 谢谢东岳老师~
在cfdonline上面也找到了类似的问题 https://www.cfd-online.com/Forums/openfoam-post-processing/125782-how-find-face-area-each-cell-boundary-patch.html
非常感谢

c7a27198-c865-41aa-bc5f-ce28c85f32d3-1728486466827.jpg

请教各位大佬专家，立方体里边，挖掉一个圆柱体，对圆柱顶面进行切割，上半部分为流体入口边界，请问又什么好的网格划分思路吗？

单纯立方体里挖一个圆柱能够实现，但中间多切了一条线，不知道怎么操作比较好。

另外，有没有好的网格划分方面的教程资料推荐？

恳请大佬们指点。

这个人2018年发的贴，2020年更新了一下，转给大家看看

到南京高校不到两年，接连评上了副教授、博导、教授～原本在北京评个副研究员都费劲。大学老师比科研院所，舒服太多了！
过去一年半，是我过得最幸福的一段时光。
对应三重压力都没了。具体来看：
收入原来的累死累活拿到30个，现在轻松30个。原来各种单位和领导的任务，现在没有领导了，只有自己的想干的任务。成果质量数量都有明显提高。
当然也有不爽的地方，比如不如北京有契约精神、红眼病严重、管理水平差的一比等。估计和我到学校不是太好有关系。
但总体来看，还是比北京爽太多，才是我想要的有钱而自由的生活。不过，下次可以考虑去个更好点的高校。

0_1545199915999_1.jpg
不知道为什么不支持gif图。。只能上传一个静图了。

这几天在看有限体积法的书，大佬的论文、《数值传热学》和 《The FVM inCFD OpenFOAM..Matlab》都看了，逛论坛发现都在推荐《Computational Methods for Fluid Dynamics》这本；去SpringerLink一搜正好今年出了新版，网上资源好像还不多，在这里分享一下： 点击下载

可以设计一个量纲的类，里面包含一个整形数组来记录量纲的分量，如[L, M, T]还要定义量纲的运算。然后设计一个标量的类，里面除了浮点数据之外，再加一个量纲，定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0]，速度的量纲是[1, 0, -1]。

@李东岳 感谢李老师，我试试看能不能画出来

@李东岳 好的，谢谢李老师

@mellonie 1. 我推测软件湍流强度 intensity 的计算是根据这个经验公式进行的
I = 0.16(Re)^-0.125
这个公式是从pipe flow 充分发展的案例中得到的经验公式。
请确认一下湍动能和湍动强度。 Fluent 中默认的湍动强度是5% ,所以你说的计算得到的是4%应该是通过一个计算体系得到的。
文献是0.01%，可以追溯一下文献的计算方式。

设置成Intensity and Hydraulic Diameter也是没有问题的

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢

I want to keep my bike close by, so I can fix it

I want to take off (remove) the bike wheel

Allen key

It doesn't fit

Move over, scoot over, can I get passed, please

She likes more and more like me everyday

I rarely (don't often) watch TV

Why did you set the temperature to 20 degrees

Dont rush me

I want to get rid of the extra oil

I can take some, I will finish what's left

She has a loose tooth, and its goin to fall out

She is not fully awake, she is a bit groggy

She is the one who loves me the most, she loves me the most

Put on the striped onesie

Move forwards and backwards

Her ton nail is broken

This is her first time seeing it

This drink is non-alcoholic, you can have a sip

Take a big step forward

A loose thread

Lift my leg and put it down

Stick out your tongue

She is lying on her side

I did it all in one go

Do you mind if I add her to our chat?

Please comment on the sentences above

Look at the lines on your arm. You leaned against something, that left lines on your arm.

Pull the handle to open the wardrobe

She just wants to be held up one your shoulder, she just wants to be held close to your chest

I am afraid of the dark

When she is sleeping, how does she know who is beside her

She does not care about anything, she is not easily bothered by things, she has high tolerance. I have a high pain tolerance.

If someone's ideas can be easily influenced by others,

She does not want to lean back

Give yourself a round of applause

Candy wrapper

I am meshing the pataya？？？

Are u scared of mosquitoes biting you?

Do you hear me calling you?

Another spoonful

How about I weigh you?

并行计算的问题，一般是计算过程出错就会出现这些error文件。
删掉没啥影响，重新打开并行计算就可以了。

weblate ，或者 github。
debian中文手册翻译
但学习成本都挺高。

但是特性好，版本控制啥的，weblate 就是个轻量的版本控制系统。po 文件可以统一术语。后台就是这个软件 gettex

曾经想参与 debian 手册的翻译，没时间入门这些东西。:xinlei:

也有会用软件一个人带着其他人翻，用邮件列表报翻好的片段，然后统一上传。

对这种不是那么重的翻译，用专业的方法成本太大。那种翻译框架搭好了是照着往多种语言翻译的。用土法对个人压力过大，这种压力不只是时间紧迫性和工作量上的。没有好的工具，hold 不住相当碎片化的翻译数据，邮件勉强可以，版本控制器最适合（全历史记录）

可以扣搜扣搜有没有日文版和日文翻译者，取取经:chitang: 一般好多技术文档有国际翻译，八成就有日语翻译。
日本 OpenFOAM 大本营：googlegroups
组织上决定派你去踩点，加油:chigua2:

@一二 都是用matlab做的

感谢东岳老师，我尝试着改小了，问题得到了部分解决。

@东岳 谢谢大佬~

一方面是动网格变形或更替方式不合适，另外我觉得是中间细小的几何部位网格有点粗糙，可以试着加密局部网格，毕竟2D计算，网格资源消耗也不是很大。

导出时候选择了如下这些选项
image.png

说一点我的看法，不一定对，仅供参考：

non-orthogonality

dca809a8-9cd5-4f74-9329-6d8459966264-image.png

non-orthogonality由snappyHexMeshDict.meshQualityControls.maxNonOrtho直接控制，你的snappyHexMeshDict里这个值是45，这个值一般不需要控制的这么小，65是更常见的配置值。

另外，你的snappyHexMeshDict里，relaxed.maxNonOrtho给的是75，这个值通常是和maxNonOrtho = 65相配合的。如果你真的要把maxNonOrtho设为45，按逻辑来说relaxed.maxNonOrtho也应该相应的减少一点。

网格纵横比

这里截取你的背景网格配置：

vertices ( (-340 -302 -3) // Slightly smaller than the STL model bounds (332 -302 -3) (332 341 -3) (-340 341 -3) (-340 -302 95) // Slightly larger than the STL model bounds (332 -302 95) (332 341 95) (-340 341 95) ); blocks ( hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1) // Adjust mesh density as needed );

如上面B老师所说，你网格的大纵横比就是这里导致的。如果不是有意要在Z方向上做加密的话没必要这么分块，我口算一下，大概(70 65 10)就可以。

castellate

从上面可以看到，你背景网格的量级是100 * 100 * 100 = 1,000,000；而你SHM的maxGlobalCells给了2,000,000。这看起来不太对，这几乎没给SHM进一步细分的空间，从结果来看也是这样，网格划分的不够细。

另外，截取一段level的配置：

features ( { file "building.eMesh"; level 3; } ...... ); refinementSurfaces { building { level (3 4); } ...... }

我理解features的level应该不低于Surface的。

胡乱写了一些个人看法，希望能有所帮助。