CFD中文网

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

@fubianhanshu 好的，谢谢您，我去试试

@李东岳 WALE不需要vanDriest吧？

@ibelief 可以 发呗

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

数值计算方法

具体是什么，解PDE么，还是什么

如果能实现顶盖驱动MAC算法，simple也没有问题啊~可以先试试MAC，速度边界条件都是第一类边界条件，其余都为0.只有顶盖有速度。

@wsy11 谢谢您~我现在已经可以成功构建网格啦^^

多面体明显质量好过关 过基础计算关 相对正六面体 同样参数控制 数量变小了

1.JPG 捕获.JPG
就是这几项在方程好像没有它们 它们的作用是什么？

你查一下这方面的实验论文吧？
壁面速度低，容易造成管壁聚集

@FluidGao 嗯嗯，感谢。

@程迪 请问linux如何并行计算，在ubuntu下没有问题，但是在linux下并行，显示这样的错误，希望您指教

OpenFOAM-6(23) mpirun -np 2 interFoam parallel -------------------------------------------------------------------------- The value of the MCA parameter "plm_rsh_agent" was set to a path that could not be found: plm_rsh_agent: ssh : rsh Please either unset the parameter, or check that the path is correct

@sr0922msc 我没有尝试过运行4.0之前的版本，即使是4.0版本，我在运行tutorials案例的时候，也遇到了一些问题，有些可以运行，有些不可以

我大部分情况下模拟结果都是这种
上面呢种是因为激波前移叠加了？
不太懂啊
Snipaste_2021-04-02_21-04-11.png

@李东岳 谢谢老师解答，之前一直出现这个问题导致每次画完网格都得要重启，我重新装下系统在试试

已解决。。。先给出ADD的表达式，然后对其进行限制。修改代码为：
屏幕截图 2021-04-07 112950.png

1.你这个不是动网格模型吧？ 也就是船体自身不发生运动
2.如果是1 的话，就和飞机机翼在空气中的模拟类似了，只是你用了VOF
3.这种情况下，垂直船体收到的力是不是本来就应该不是一个稳定数值？

@xpqiu 感谢您的答疑
原来如此，是两个box相当于框定了一个水域区范围间的意思。
zoneID场的值设置明白了，还想问一下zoneID这个场是不是因为重叠网格所以要设置的，具体作用是什么呢？

@wangfei9088 哇，太谢谢了。我也考虑过是单位的问题，但还未在代码中找到数据验证猜想，非常感谢大佬！