CFD中文网

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

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

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

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

恳请大佬们指点。

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

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

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

@队长别开枪 感谢感谢 😊

@ibelief 在新版本的帮助文档中写了，运动区域的夹角必须是90，120,180这样可以被整除的数。

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。

密度就是单位体积的质量，速度的物质导数就是加速度。
所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？

力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。

物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。

要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。

另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。

/OpenFOAM/OpenFOAM-2.4.0/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/T.water
这是我想修改的温度文件的位置。
初始蒜粒里面wall给的边界条件是zeroGradient.我想做一个壁面散热的效果（壁面以热通量W/m2的形式向外散热）很明显zeroGradient不再适合。
我想求教：这种热通量向外散热的边界条件该如何设置？是采用fixedGradient还是其他形式的壁面函数？
如果是fixedGradient，那么它的具体值的单位是什么？如果采用其他方式，那么有没有能使twoPhaseEulerFoam运行起来的成功的壁面函数？如果有，还望各路大神不吝赐教！！！
最后还是麻烦东岳老师帮我看一下，或者号召一下大家伙，看看我的问题。拜谢拜谢！
@东岳

亲测，可以对齐~如图，word2016，自带的公式编辑器！
0_1520851049441_d3fa154b-70c5-438f-9674-805ee760bf75-image.png

同问该问题～

嗯嗯，矢量分量等值面有了，谢谢！

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

@东岳 谢谢李博士，我感觉好像是不算的，因为在后处理的过程中，MassFlow的进口和出口的所有流体的质量流量值和液相的相等，但是固相的就是固相的质量流量，好像不包括，好像DPM模型中的液相和固相是采用分离流计算的，那我们在设置出口质量流量的时候就只考虑液相的质量流量，不过可能这个是这个模型只用于计算固相体积分数较少的情况

@东岳 谢谢您！

@昝磊 @史浩 百年好合！ :laile:

二变量高斯分布：
\begin{equation}
f(u,v)=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left(-\frac{1}{2(1-\rho^2)}\left[\frac{(u-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(u-\mu_1)(v-\mu_2)}{\sigma_1\sigma_2}+\frac{(v-\mu_2)^2}{\sigma_2^2}\right]\right)
\end{equation}
MGF为：
\begin{equation}
m_{i,j}=\exp\left(i\mu_1+j\mu_2+0.5(\sigma_1^2i^2+2\rho\sigma_1\sigma_2ij+\sigma_2^2j^2)\right)
\end{equation}
纯矩计算方法为
\begin{equation}
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=\mu\\
m_{2,0}&=\mu^2+\sigma^2\\
m_{3,0}&=\mu^3+3\mu\sigma^2\\
\end{split}
\end{equation}
假设$\mu_1=10,\mu_2=20,\sigma_1=\sigma_2=2,\rho=0.5$，有纯矩：
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=10\\
m_{2,0}&=104\\
m_{3,0}&=1120\\
m_{0,1}&=20\\
m_{0,2}&=404\\
m_{0,3}&=8240\\
\end{split}
同时有混合矩
\begin{equation}
m_{i,j}=\exp\left(36\right)
\end{equation}
不行，混合矩计算方法不对

只算过解析解，不知道你用的是什么方法算flowmap?

一个单核代码能出错其实还是不常见的

今天在Peric那本书里面看到说

如果能量守恒的话，速度在计算域内应该是有界的

page 161

@King_RVM 你好，使用starCCM绘制网格，然后导入到openFoam里面跑算例，计算结果会出现质量不好的情况么？请教一下~