CFD中文网

最新进展：我将0文件中的pd边界条件进行了修改，可以实现中断后继续运算，但原问题依旧不知该如何解决。

代码如下：

/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | foam-extend: Open Source CFD | | \\ / O peration | Version: 4.0 | | \\ / A nd | Web: http://www.foam-extend.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object pd; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { out { //type zeroGradient; type fixedValue; value $internalField; } "(sym1|sym2)" { type symmetryPlane; } in { //type fixedFluxPressure; //value $interField; type zeroGradient; } cylinder { //type fixedFluxPressure; //value $interField; type zeroGradient; } "(front|back)" { type empty; } } // ************************************************************************* //

Hi，

我遇到的使用三菱柱的情况就是四面体的边界层网格。当然需要使用三棱柱来做边界层，有利于使用三棱柱的上下面把误差抵消掉。

计算效率是指？
目前出了纯结构网格，其他的网格在做索引的时候编程都要比结构网格复杂。
如果值得是精度，三棱柱在边界层肯定要比4面体好，这个在教材中都有提及。

没想到图越来越多。手机看的注意了。我这个图多。

结构网格的拥护者简单说一下为什么我喜欢结构网格。以下只针对简单几何以及不是太复杂的几何。复杂到不行的网格我选择混合网格。我自己做的算例不用纯四面体网格。

六面体网格具有以下优点：

非矩形网格的歪斜（skewness）严重
这需要做面法相梯度计算修正。稍微增加计算步骤。著名的CFD非矩形修正因此而来

截面重构不光顺
举例，一个正方盒子，一半是水，我们用矩形网格重构的气液截面是这样的：
alt text
黄色为一相，白色为一相。即使网格很糙，我们的截面也是平的。下图为三角形非结构网格：
alt text
可见糙网格下界面不平。这导致后处理需要耗费精力。

高度非线性物理问题影响矩阵对角占优
这也是非线性动力学，例如金属撞击、生物力学分析中结构网格肆虐的原因。（这个我没有测试，之前在网上看到的结论，来自链接文本

非结构网格数量更多增加计算资源
网格单元的面越多网格越少。六面体网格并不是数量最少的。数量最少的是多面体网格（star推荐的？）如下图
alt text
再看这个情况对比：只有第一个是结构网格，第二个网格多。第三个第四个和均为非正交且和流线不贴合
alt text
alt text
因此目前自动生成的非结构网格一般数量是即为巨大的。计算的时间会超过你花一星期做结构网格的时间。另外，万一中途发散了怎么办。

费结构网格高度不可控
几何稍微一点点的变化会导致全场四面体网格做位移。结构网格（ICEMblock切出来的）具有即为即为高度可控性。可以对基本所有的网格单元进行控制。俩者区别天差地别。

我了解的所有的离散格式对结构网格适应性更强
也就是说结构网格对流场具有强的多的顺应性。例如旋转机械，结构网格和速度流线高度贴合。

矩形网格对边界更好的捕获
为何非结构网格要添加边界层网格？矩形网格对壁面贴合更好。

复杂CFD模型在结构网格上收敛更强
个人经验。比如双流体模型，LTS局部时间步框架下的模型，激波捕获等。当然那种不可压缩单相流不算

太多就不写了...

引用从业CFD27年OpenFOAM创始人Henry的一段话：

hex and polyhedral cells are generally superior to tet cells in finite-volume and the difference in both accuracy and convergence behaviour can be large. However, it also depends on the quality of the hex and tet cells and many of the problems with tet meshes stem from the poor quality of some tets produced by many automatic tet mesh generators. In order to really appreciate the difference it might be a good idea for you to try both for your problems. For VOF calculations I would expect the solution on hex/poly cell meshes to be dramatically better.

大意为：

在有限体积中，六面体和多面体网格相对于四面体网格优势明显。主要在精度和收敛上。对于VOF（我上文提及的第2点），结构网格的优势是大大的。

最后，需要耗费大量时间生成的六面体网格至今没有消失，必然有存在的道理。

通过阅读东岳老师以前的一篇文章：拉格朗日中的湍流分散力模型，我的理解是这个随机游走模型不应该是在粒子上添加一个随机游走速度吗？但是代码的最后有关Uc+UTurb，如果我理解没错，这不是连续相的速度加上一个湍流脉动速度么？难道是通过随机更改流体速度，间接再去更改颗粒受力，从而达到颗粒的随机？

你的理解是正确的。他的随机，主要体现在随机数的生成上（不是一个固定的数）。

楼主您好，我也是做流化床传热传质fluent模拟的，我做的生物质流化床气化模拟，您模拟做的怎么样了，能给一下您的联系方式（QQ）吗？有太多问题想咨询了，可以有偿咨询哈! 感谢感谢！

我看看

flowRatePatch自动适配，phi是质量通量就计算的质量通量，phi是体积通量就计算的体积通量

好的 东岳老师。

我已经找到原因了，是因为材料粘度引起的沿程损失低于重力引起位压时就会出现这种出口比进口压力大的现象。

问了一下别人，好像是使用snappyHexMesh命令画2D网格就是会出现面贴合效果不好的情况

谢谢您的意见，提取的肯定就是面上各个点的速度数据了，那如果这个时候求液相平均速度，是否就是根据液相体积分数做速度的加权平均呢？

0_1460680809776_捕获.JPG

苹果retina屏那个本也不错。就是能不能装win系统？

http://item.jd.com/2819597.html

@小八 ICER，文章已经找不到了。化工方向的:+1:

按照教程编译结束了老师，但还是提示这个错误。

@xpqiu
.................
25T硬盘确实花了不少钱。:sunglasses:

大家好，我是一名钻探的研究生。我想研究这样的一个情况：一个直径40mm,高度1000mm的圆柱，在圆柱底部每秒产生一定质量的岩屑，岩屑直径1mm，这是岩屑的产生速度，岩屑的实际运动速度为0，同时在圆柱底部有速度已知的空气，现在的问题是想研究空气的速度能不能把不断产生的岩屑都吹出来。我现在发现dpm可以实现这个问题，但是dpm的后处理出的图太少，写论文不够用，欧拉模型出的图挺多，但是不能解决我的问题？大家可以给我出点主意吗？还有别的模型可以解决我的问题吗？

@showhand 多谢大神回复！

请问有前辈知道gmesh如何能划分如图所示得网格吗，图1是图2的block

block-tuya.jpg model-tuya.jpg

很遗憾，木有啊:sad: