CFD中文网

@yupeng https://openfoam.org/download/source/
你按照这个网站的流程试试呢

@李东岳
嗯嗯 我试试 感谢李老师

@yupeng
建议foamSystemCheck一下，如果是Ok的那就删除所有OpenFOAM相关的文件重新按照流程走一遍再编译

@李东岳
但是采用RANS-simple迭代了10000步也没有显涡

@cccrrryyy
所有云图都是颗粒固定的情况

@李东岳
1、没有涡显示
2、处于超临界状态下的水，密度80.54，粘度小
3、湍流模型为：KOmegaSST
4、网格如下：

这是剖面，实际为全区域三维建模

颗粒局部：

5、虽然速度小，但是粘度低，对应的雷诺数为5715，属于湍流状态；所以采用湍流模型

@李东岳
感谢李老师回复
1、Y+=2.5左右，位于粘性底层，所以没有使用壁面函数
2、网格使用blockMesh生成的，可能密度不够
3、湍流粘度如下图所示：
确实很大，我设定的湍流粘度比为1；nut是计算出来的，应该如何修改其他的相关参数呢。

之前设定K与Omega如下图：

各位老师，最近我在看湍流模型相关的知识
有几个问题想要跟大家请教一下：
一：对于湍流模型都是处理湍流充分发展的流场，在OpenFOAM中k与Omega文件的初始条件与边界条件足以体现。
问题如下，对于静水沉降问题，达到稳定沉静速度对应的湍流，能否使用湍流模型？如若使用湍流模型初始场设定相关的k与Omega参数，与本质的物理问题相矛盾吧？
是不是只能采用DNS计算或者在OpenFOAM中不附加湍流的quasi-DNS计算。

二：在看湍流模型时，对流体绕流固定颗粒，采用湍流模型进行了一些尝试计算：
Laminar-瞬态-pimpleFoam
RANS-noWallFunction-KOmegaSST-simpleFoam
URANS-noWallFunction-KomegaSST-pimpleFoam

分别进行了上述计算且结果如下：

（1）其中湍流模型KOmegaSST：
1、确定湍流参数：Uinlet=1.08 （雷诺数Re=5715） 湍流强度FSTI=1%（low） Eddy viscosity ratio=1
2、初始文件参数-nearwalltreatment：
K：fixedValue 0 初始值：0.00017496
nut：nutLowReWallFunction
Omega： OmegaWallFunction
3、湍流参数初始化如下图：

（2）经计算：以上三个计算Y+均为2.5左右小于5，处于粘性底层。--关心壁面处的受力所以不采用壁面函数处理，使用LowRe壁面处理

（3）结果分析：
① Laminar-瞬态-pimpleFoam：
速度云图如下：
曳力系数如下：
②RANS-noWallFunction-KOmegaSST-simpleFoam：
速度云图如下：
曳力系数如下：
③URANS-noWallFunction-KOmegaSST-pimpleFoam：
速度云图如下：
曳力系数如下：

可以看到使用RAS湍流模型之后，并没有体现涡的结构以及瞬态性，曳力系数并无周期波动。
问题如下：不使用湍流模型反而可以观察涡的结构以及波动；虽然RAS是进行了平均处理。但是大尺度的涡应该在速度云图有所体现吧。

请各位大佬答疑解惑，感谢感谢。

@袁宝强

haha~ 跟论文原文是可以对起来的~感谢！

下图是Opanfoam中KOmegaSST湍流模型中一些常数的取值：

第二张图片是资料《The Finite Volume Method in Computational Fluid Dynamics
An Advanced Introduction with OpenFOAM®and Matlab》-17.5.4章节（P705）中介绍该湍流模型的参数：

为什么两组常数存在部分不对应：（红线所示）

还请各位老师解答，拜托拜托。
OF版本： OpenFOAM-V2106

@fu 在使用laminar（Quasi-DNS）计算的时候，壁面处的网格Y+<1 ,且整个区域网格也要相对较小；当使用湍流模型的时候，壁面处理跟laminar应该不同，使用壁面函数之类的；如果题主使用一套网格应用于两种流动模型，是否会导致计算结果的偏差呢？还请大佬指点

@wsxfyy 修改controlDict之后，在下一步计算中，就会采用新的controDict中的数据，是可以动态更改的。

@wsxfyy 直接修改文件的内容就可以了，比如你想改库朗数，保存间隔之类的直接修改数字就可以了。

@dingcy
https://www.cfd-online.com/Tools/turbulence.php

@cccrrryyy 老师好，我还没用湍流模型，使用层流模型计算过，最终沉降速度（对应颗粒雷诺数5000）大约有30%的误差.

@李东岳 计划用overInterDyMFoam 求解器

@李东岳 就是单颗粒使用动网格+重叠网格，颗粒的自由沉降过程

@zhouxu颗粒周围的流体状态会影响颗粒的受力呀，会影响颗粒下落过程。

@李东岳 麻烦李老师指点

各位老师，大家好
如果实现一个静止颗粒在静止流体中自由沉降，速度由0增加至自由沉降速度，自由沉降速度对应的颗粒雷诺数为5000多，在这个过程中雷诺数逐渐增大，流动状态也由层流过渡为湍流。
如果关注颗粒的加速过程，应该选用什么湍流模型呢？
是充分发展的湍流模型： kOmegaSST
还是转捩模型：kOmegaSSTLM
还请老师们指点