CFD中文网

https://www.cfd-online.com/Forums/openfoam-meshing/133850-snappyhexmesh-problems.html

icem网格划分导入fluent时检查失败，显示Error: interface zone 3 has two adjacent cell zones.
4f68eb65-132e-4314-b0f6-6964dc466673-image.png

就是这个扩散通量的正负。一般扩散项到等号左边是付的。

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

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

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

@李东岳 对 有

@李东岳 居然和我方向那么像，流化床和燃烧机理。

@cfd-china 直接用WIN自带的远程连接就可以了

对称几何对称网格对称边界条件，但最后模拟出来的结果不对称的问题。

这个问题是会有发生，就类似圆柱扰流，注定是不对称的。不过如果是一个可以对称的物理，加上一个完全对称的数值设置。发生不对称的结果。就不太好说了。

各位大佬好，模型是一个0.125x0.5x0.6m的光滑通道，上下壁面均为无滑移的300K恒温壁面，进口温度320K，进出口及左右壁面均采用周期性边界条件，模拟步骤是先用RANS模拟，之后使用TUI指令添加扰动然后切换到LES模拟。按理说流体温度应该小于等于320K而且随时间逐渐减小，但是云图画出来最高温度有323K，又模拟了0.5s后最高温度没什么变化其他地方的温度似乎也没什么变化。周期性设置如图（压力梯度最开始为0）搜狗截图20230903172308.png
我想知道是为啥 :136:

@cfd-china 好的

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

多谢，已更新

@李东岳 好的谢谢您，李老师

5右边多写了个$\rho$吧，谢谢！回北京我更新一下

请问一下各位大佬，欧拉-欧拉气固两相流模拟，设置流体域旋转后，wall可以设置成不旋转吗，对结果有没有影响，因为wall设置成旋转之后，就无法设置镜面反射系数了

用google搜索吧。

你现在可以按照提示操作一遍。如果还不行，把你的错误输入到google中去，一般都会找到满意的答案。起码我目前为止是的。

没有google就花钱买v屁恩

我在做gas-liquid compressible flow时，需要water的equation of state 去解液态场的压力。
我知道现在有这样几种方法：
the stiffened-gas equation
the Tait equation
Noble-Abel Stiffened-Gas equation of state
想问一下大家有没有这方面的经验，因为有很多系数，不太清楚怎么选择系数。比如：
0_1520637932983_Screen Shot 2018-03-09 at 5.25.10 PM.png
关于gamma和P_infinite的选择会影响很大

谢谢大家！

网格加密到Y+<1，计算结果变化不大，对流通量用二阶MUSCL, AUSM格式

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

@mohui 如果使用湍流模型，那么就是模化（或部分模化）了流动的脉动。真实流动的脉动越强，湍动能也更大；但是你计算的时候是不会算出来流动的脉动的，只会得到较大的k值。
你在k方程加扰动，对结果影响应该不会很大。