CFD中文网

非常感谢

topoSet有一些比如寻找点最近的面，寻找距离stl最近的面之类的，你试试

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

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

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

非矩形网格的歪斜（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点），结构网格的优势是大大的。

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

DPM限制.PNG
如上图，帮助文档中在dpm模型下使用mrf、动网格、sliding mesh会导致颗粒追踪失去意义，请问在dpm模型中使用srf会不会导致颗粒追踪失去意义？或者说在fluent里是不是把srf归为mrf？

用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}

OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform $epsilonInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type kqRWallFunction; value uniform $kInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
nut
internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}

OUTLET { type calculated; value uniform 0; } WALL { type nutUSpaldingWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
p
internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}

#includeEtc "caseDicts/setConstraintTypes"

}
Uinlet (10 0 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type noSlip; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
时间步长试过1 0.1 0.01 其他的都没动

想找一些适合cfd入门学生学习的题目
最好有把cfd的各个环节单列出来的题目（例如 贴体坐标系生成）
以及这些问题用什么工具求解
我是指学习使用哪些工具对于科研实际具有意义（例如使用哪些软件（fluent？turbo？） 还是用哪种语言自己编写代码（例如cpp）？有哪些部分已经出现类库，哪些需要自己编写呢？尤其是科研中是怎样实际操作的）

直接的宏应该是没有的 你可以自己用UDF写个类似于level set的函数计算每个位置到界面的距离

@东岳 我是按这个公式计算的，密度为1，那么u_tau=sqrt(tau_w)，但是我发现，论文里的摩擦速度刚好是我用OF算出来的1.225倍。论文是用fluent算的，fluent里面的空气密度是1.225。我现在想不通 为什么会差这1.225倍。:jingya:

你好，

这个问题我最近研究过。请参考：CFD中的高阶精度是不是扯淡。在这里面已经讨论很详细了，此处不赘述。

大家好，我最近在模拟波浪。在现在interfoam底下的case里，对于波的设置条件里，有一项ramptime，我不清楚它的意思，请大家帮忙解答。谢谢🙏

你是不是想把曲面展开为一个平面？https://blog.csdn.net/mrbaolong/article/details/106949340

@showhand 多谢大神回复！

请教一下你是每一层加密都对三个方向进行加密吗，shm中能能不能只对z方向进行加密呢。@一朵苍凉

@卡洛 请问你是用的哪种方法？

只有在可以设置 flux 的边界条件上，才可以指定梯度边界。
比如，温度的边界条件可以指定 heat flux，这时你就可以用 DEFINE_PROFILEDEFINE_PROFILE 定义一个值，挂载在 flux 的边界条件上。
UDS 也可以指定 flux，同理可以用 UDF 指定梯度边界条件。

大家好

f29cfccc-509a-4c69-bbd0-3eddb612ac78-图片.png

如图所示的二维空间，蓝色圆形的导热系数与白色区域不同，我可以用blockMesh生成规则网格，利用funkySetFieldsDict在已知圆柱中心、半径情况下定义导热系数。

对于三维问题，我利用3dmax生成了不同直径的堆积球，球心、半径的信息我不清楚，如图所示生成stl文件：

d750e2b8-86fd-41c9-aea4-216a88ba6a85-图片.png

接下来我想计算整体的导热问题，这也是一个各向异性导热问题，我仅需要用blockmesh画网格即可，但如何定义球的物性？

用 DEFINE_SOURCE 宏，根据需求挂载在 k 或者 epsilon 的方程后面，试试

具体是啥阀呢，阀内流体是水还是油，我可以试试，工程问题可以用cfd计算一下阀内的噪声，但都是参考，可以通过cfd进行可视化研究，可以看到里面的流场和声场，具体咨询合作可以加微信18647841059

目前想通过编写“.js脚本文件”在DesignModeler中创建三维几何模型，之后再划分网格用Fluent模拟。求教“.js脚本文件”的编写方法，教程资源。感谢！
1.png

多相流的有限体积法跟单相流没什么区别啊