CFD中文网

结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说，结构化和非结构化区分是文件保存格式的区分，结构化是使用索引，非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格，因为我最近也在学openfoam和icem，icem可以通过自动化分成四面体非结构网格，也可以分块化成6面体网格，但是我通常画完6面体网格后，icem有可以转化成非结构网格按钮，将画好的网格转化成非结构体网格（转化后依然是6面体）。
所以如果你要话非结构化网格也可以分块（画了6面体转化一下），如果你说的是画四面体非结构化网格分块的话，这方面知识我就不知道了

@赵一铭 目前想通过现有商业软件做一做···openfoam有点难，项目要求又很急····

@李东岳 好的，谢谢老师

@李东岳 alpha.c4这里的c4就是作为爆源的高温高压气体，好的,那我去看看怎么根据压强温度密度计算这个值，因为看之前的0文件夹中的都是alpha.c4 1，所以我以为是高温高压气体站所处方块的体积分数。非常感谢您的指导！
f106b6c610998150ac697385f731b6c.png
计算前
04f98e3c28e1bfb0ba959238d02940c.png
计算后

udf可以进行

多谢老师解答，问题已经解决了

在帮助文档找到了模型。:xinlei:

@李东岳 这个应该就是buoyantPimpleFoam的处理方法

@李东岳 李老师，在写fvm::laplacian项之前已经通过体积分数alpha计算过nu了。因为我觉得直接修改turbulence::类里面的nu很麻烦。所以我就自己定义了提及标量场nu然后正常使用alpha计算之后带入laplacian，所以获得了与文献一致的结果。

@JokerZhu 是的，但是githubfast已经关停了好像，去github官网搜就行了

流域里是有流体的，球受到的阻力也是随速度增大的。从静止下降1s后对称面上的压力分布如下图所示，球的直径是1m。等有时间了会做个与正常6DOF模拟（就是让网格变形那种的对比）的对比。20200722114110.png

您好，我想在您这个回答的基础上再问一个问题。由于在rigidBodyMotion求解器中，centraOfMass是相对于centraOfRotation，那么sixDoFRigidBodyMotion求解器中的momentOfInertia改成rigidBodyMotion求解器中的interia应该做何种变换呢？

@卡洛 无痛苦ns方程里面有，专门讨论的flux()函数，你看下

测试出来了QAQ：
𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑠𝑜𝑢𝑟𝑐𝑒=(𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑐ℎ𝑎𝑛𝑔𝑒×𝐶𝑢𝑟𝑟𝑒𝑛𝑡 𝑝ℎ𝑎𝑠𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛)/(𝐶𝑒𝑙𝑙 𝑣𝑜𝑙𝑢𝑚𝑒 ×𝑇𝑖𝑚𝑒 𝑠𝑡𝑒𝑝 𝑠𝑖𝑧𝑒)

Eg: Bin 1 bin fraction decrease 2% in 0.1s, cell volume is 0.001 m3, inclusion phase volume fraction is 1%

322c153e-e638-4775-9e9a-4e0939981818-image.png
2 is the number that we give fluent, 0.02 is actual bin fraction change.

@赵一铭 整个程序我用的SIMPLE算法，所以连续性方程改成了压力修正方程来求解压力，不能直接求解出密度，怎么办？

reFineMesh可以加密一个方向，这个需要指定区域，你试试topoSet+refineMesh呢》

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

您好，请问一下，您知道有什么好的方法可以提取重叠网格的位移吗？能否分享一个算例呢？@孤岛

@Prometheus10 谢谢您的回复，我找到原因了，原来我之前没有进入mnt文件夹，所以看不到自已本来的硬盘，进入mnt后就能正常使用了。

用的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 其他的都没动