CFD中文网

@CiCi 同是做可燃介质的燃烧爆炸的数值模拟，应用背景放到了安全而不是燃烧利用。我想我猜到师姐在哪个组了。上网不涉密，涉密不上网（手动狗头)。目前这个方向的话大家好像都在往锂电池的热失控排气及诱发的爆炸模拟上转。

@alex 看到你的网格了，感谢分享，感觉划分的非常不错，好厉害！学习了！

如果招聘信息可以提供更全面，待遇可以明盘的话，我可以在CFD界推一下

是的，理想气体，明天我更新下。

另外，温度本身不是守恒变量，因此，CFD一般求解守恒变量，但是温度方程通常比较简单，一些不太成熟的代码会使用温度防尘

老铁是不是毕业了，后来咋样了

:baobao: :baobao:

@bestucan 非常感谢老师指导:xinxin: 我回去试一试

我用系统自带的网格（ 1 75 1）可以算出来，但是我改成（30 2000 1）就算崩溃了，打算算高2m宽15.4mm的圆管截面

请教一个比较初级的问题，查了帮助文档没有找到想得到的答案。
ICEM采用Octree方法划分非结构网格的时候，帮助文档和网上的一些教程中都写到：“建议全局网格参数设置中的Max element值要是2的幂”
这里2的幂是指：2的正整数幂？2的正数幂？2的实数幂？
实际测试过感觉比如填0.05啥的也可以划分，网格也确实要比2的密很多。所以想问下各位高手这里这个数值填写应该有什么经验或者规律呢？尤其是模型尺寸很小的话，有时候确实需要很小的尺寸啊~

各位CFD大佬，小弟最近在用fluent做模拟时遇到一个问题，物理背景就是管内的热水以射流的方式喷出然后被周围的环境冷水稀释，在做小模型模拟（x=40m y=2m z=13.3cm）的时候，模拟结果尚可，对称性很完美，但在做大模型的时候（x=4000m y=2000m z=10m）,模拟结果很不好。
用的稳态模拟，大模型网格数量大概七百多万，用的reliazable k-e湍流模型，进口的湍流参数设置用CFD计算器估算过了，另外在模拟的时候fluent一直有湍流粘度比超限的提示，但小模型就不会有，不知道最后的不对称结果和这个是不是有关系。
1564dc71-8ed2-4905-ba1b-67fbca1936e3-图片.png
热水从圆的四周均匀喷射出来，然后被流动的冷水稀释，冷水流动方向从左至右，可以看到图中并没有达到一个对称效果，等值线都没有对称。
21f863e0-2b3c-4e9f-994a-ae3797889acf-图片.png
热水进口在上图的最底部，垂直向上喷射，中间比较黑的那一块有一块顶盖，作用是让热水向周围扩散开来而不是垂直向上冲，图中看到网格有些瑕疵的地方是tecplot的显示问题，可以忽略。
a652dbd8-6ce0-4c30-9c4a-97e6fee8e965-892fea64760b5fcdfe9299ba371783e.png
从这张图可以看出热水向上喷射，被顶盖挡住，最后向四周扩散上浮的过程。

总结一下几点疑问：
1.模拟时fluent的湍流粘度比超限警告是否和最后模拟结果的不对称有关系？
2.小模型模拟结果对称，而大模型就不行，和网格的大小有关系吗？小模型的最大网格大概也就在厘米级别，而大模型的最大网格一般在米级别（最大有五米左右）。
其实最核心的诉求还是想解决对称性的问题，最好不要用对称边界条件，能用整个模型就用整个模型。

@星星星星晴 可能我得太复杂了，一共200多页，用他那个就出错

我知道了，公式里面用不了，文本里面可以用

好吧，我查了一下，附加一个yhmath这个包 :136:

https://tex.stackexchange.com/questions/100574/really-wide-hat-symbol/101136#101136

@李东岳 网格是用fluentMeshing画的，严格按照周期性边界设置的，导入OpenFOAM之后，当matchTolerance很小的时候，使用createPatch设置周期性边界的时候会报错，

以下是 createPatch 代码

patches ( { name cyclic1; patchInfo { type cyclic; neighbourPatch cyclic2; transform rotational; rotationAxis (0 0 1); rotationCentre (0 0 0); matchTolerance 1e10; } constructFrom patches; patches (cyclic_1); set f0; } { name cyclic2; patchInfo { type cyclic; neighbourPatch cyclic1; transform rotational; rotationAxis (0 0 1); rotationCentre (0 0 0); matchTolerance 1e10; } constructFrom patches; patches (cyclic_2); set f0; } );

以下是报错信息

--> FOAM FATAL ERROR: More than one patch accessing the same transform but not of the same sign. patch:cyclic2 transform:0 sign:1 current transforms:1(1) From function Foam::label Foam::globalIndexAndTransform::addToTransformIndex(Foam::label, Foam::label, bool, Foam::scalar) const in file lnInclude/globalIndexAndTransformI.H at line 210. FOAM exiting

我翻了一下tutorial中的案例，只找到了如下的一段controlDic中的function，是否是对其进行一定修改来达成目的呢？求各位大佬指教
1721204958853.png

@suzuhaz 大师求助帮忙看一下简单的一个模型块的划分0_1532943202163_05.zip

@zhe 在 在interFoam中通过setFields设置水体速度，但结果没有 中说：

@李东岳 不好意思啊，东岳老师，因为mesh还挺多，并且是我改动过的求解器，所以没发上来一起算算。因为我这个求解器是对比MPPICInterFoam加入了粒子模型。我一直没尝试过给案例里通过设置setFields里的速度，这次一试，发现不对。

我返回去用interFoam算了一下，水流没有问题了。我又用OpenFOAM自带的MPPICInterFoam算了一下，一样的设置条件。0时刻，setFields起到了作用，速度场没有问题。在第一个时间布长，速度全部归零了。看了压强，跟interFoam的结果完全不同。这两个求解器没有特别大的不同，只是在MPPICInterFoam里考虑了粒子的体积占比，但即使不注入粒子，也一样速度归零。不知道有没有谁看过MPPICInterFoam，这个在算水流时，有什么大的不同吗？谢谢

你好，请问问题解决了么？我也遇到了相同的情况。
用修改的interFoam设置不同相速度，0时刻没问题，计算第一个结果时速度场也归0😅

通过VOF和solidification模型仿真液滴结冰，有没有大佬实现过体积膨胀？

advection:可以理解为平流，主要由物质浓度存在梯度引起的变化。convection：可以理解为对流，主要是由热量不均匀引起的垂直方向的对流。好像在流体力学中对这两者区别不是很大:xiabanle:

我认为这个是和几何无关的，衡量网格质量的参数中有一个叫aspect ratio，这个值为1时最好，像那种类似薄片或者细长的网格，这个值会很大，这种网格一般质量很差。

好像是分内部点（所有面都不是边界）和边界点，内部点根据cell插值（距离倒数分之一），边界点根据face插值。