CFD中文网

@wsy11 那就是你的blockmesh 设置 把底面 ground 在设置低一些

这二者都可以得到cs=0.065，Ce=1.048，Ck=0.065是标准系数。然后$k_{sgs}$也跟这两个系数的比有关。不过取别的值也能达到这个效果。 @coolhhh 在这方面应该有更深的见解。

5f74eac5-a02d-4d83-afe1-2def61286ea5-图片.png RT，初始流场是零速度，但是发现导出的第一个时间步流体运动位置为图中涡卷起地方，理论上讲下游的各个地方速度场应该为0，但是云图上显示下游已经有速度了，而且随时间发展流体未到达的位置速度还会变。
这种现象是要怎么去理解？是因为边界条件设置导致的吗？

template<class TurbulenceModel, class BasicTurbulenceModel> tmp<fvScalarMatrix> kOmegaSST<TurbulenceModel, BasicTurbulenceModel>::kSource() const { return tmp<fvScalarMatrix> ( new fvScalarMatrix ( k_, dimVolume*this->rho_.dimensions()*k_.dimensions()/dimTime ) ); } tmp<fvScalarMatrix> kEqn ( fvm::ddt(alpha, rho, k_) + fvm::div(alphaRhoPhi, k_) - fvm::laplacian(alpha*rho*DkEff(F1), k_) == alpha()*rho()*Pk(G) - fvm::SuSp((2.0/3.0)*alpha()*rho()*divU, k_) - fvm::Sp(alpha()*rho()*epsilonByk(F1, F23), k_) + kSource() + fvOptions(alpha, rho, k_) );

湍动能输运方程中的kSource()应该是源项，它为什么是k_，请问有没有人明白？

@hachikoi 还没做过多浮体的例子，最近做了一个overset+AMR+PLIC的算例，结果还行。你的算例几何尺寸可以分享一下么，圣诞节放假的时候我可以用我的求解器试着算算。

我是个小白啊，感觉有时候处理这种空间后，反而打不开了，所以还是买硬盘吧

多谢两位老师

@kcol 相当于volScalarField，但是只定义内部网格，在边界网格上没有定义

$$
\begin{array}{l}
\mathbf{F}_y(t)=-\rho\left(g \sin \theta(t)+y \omega(t)^2+z \varepsilon(t)+2 \omega(t) u_y(t)\right) \boldsymbol{j} \\
\mathbf{F}_z(t)=-\rho\left(g \cos \theta(t)+z \omega(t)^2-y \varepsilon(t)-2 \omega(t) u_z(t)\right) \boldsymbol{k}
\end{array}
$$

这个方程里面非粗体的$\omega,\varepsilon$以及$\theta_m,T$你知道是什么么

@yhdthu 好的 谢谢大佬

IMG_20231127_193223.jpg
在计算单相流作为初场，可以收敛，但是湍流粘度比超限，继续设置为多相流空化计算的时候，开始发散，不知道有哪些解决办法，已经重新画过网格。
边界条件是速度入口，压力出口，湍流sst k-w，稳态计算

@李东岳 老师，最近我也遇到了这个问题，单相流计算的时候会提醒，但是结果也可以收敛，但是计算完初场的时候，改为多相流空化计算，马上就发散了，不知道怎么可以解决这个问题，已经重新画了好几次网格了

一个流弊的低雷诺数SA模型在CFD中文网于2023年11月27日诞生！

当气泡尺寸过小时，特别是对于水中自由浮升的微气泡，就下面这篇帖子曾说到的，interFoam算不出来。

大家主要是从壁面效应带来的阻力，解释了微气泡无法浮升。

在远离壁面的时候，也无法上升么？

如果要改湍流强度，你直接用codedFixedValue去重写一下入口InletK会更好操作一些。

@李东岳 李老师你好，是CFDEMcoupling 5.x版本自带案例的periodicchannel

我没有这两个求解器啊老铁 ..

@coolhhh 对，我使用这个命令：alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'时，还没有把anaconda的环境给注释掉，OpenFOAM也可以正常使用，但是并行命令不能使用。
我大概懂大佬你的意思了，如果使用了alias命令，只要没有输入Anaconda3应该就不会激活anaconda环境。可能是因为我这边前人在给anaconda设置环境的时候，有点繁琐（具体我也没弄清楚），所以我可能只使用一条alias命令，不能完全把相应的anaconda环境给包含进去。所以才导致使用alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'命令还不能消除使用并行命令带来的报错。再次感谢您的回答，谢谢！

@hachikoi 就是说目前还没有实现需要自己修改源码是吗

分享几个自己了解的开源软件：

OpenFOAM，不用说了。

T-Flows，也是一款基于FVM的求解器，没咋用过，听说是小型OF，没有OF那么复杂的类和求解器类型。链接文本

REEF3D，一款船海领域的开源CFD软件。链接文本

@疏影横斜水清浅 多谢老铁，将在下一版本更新

@知行 老铁我这之前忘了回你了，才看到，我处理后将在下一版本更新

@疏影横斜水清浅Yang Y -2017- Consistent inflow boundary conditions for modelling the neutral equilibrium atmospheric boundary layer for the SST k-ω model[J].pdf

@Prometheus10 老铁，你这个帖子我移动到这里了，https://cfd-china.com/topic/6754， 工作日我看看

老哥们，你们是如何解决这个问题
我在使用两套overset做浮体自衰减的时候也遇到类似的问题

@LuoS 在 snappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说：

教程中似乎有个建筑风的案例：
/Tutorials/inco

好的谢谢 我试一下

@xpqiu 感谢大佬回复，问题已解决 请问大佬有遇到过多重overset计算错误的问题吗

厉害啊老铁

@zhangxy 是一个网格的尺度

@xpqiu 谢谢您！我大致明白了，去试一试:xinxin:

@coolhhh 发现问题了，因为我的是两相流，需要设定field 为U.water

@李东岳
参考文献：
Knacke, T. (2013).
Potential effects of Rhie & Chow type interpolations in airframe
noise simulations. In: Schram, C., Dénos, R., Lecomte E. (ed):
Accurate and efficient aeroacoustic prediction approaches for
airframe noise, VKI LS 2013-03.

@Chen_hao
当这个模型用于多相流时，IOobject::groupName("omega", alphaRhoPhi.group()) 会变成类似 omega.air 这样的名字。

用 topoSet，把这些区域内的网格提取到 cellSet，提取过程会告诉你网格数量。

还可以用 topoSet 结合 createPatch 来实现，不过这样得到的面不一定是精确地 1/2，因为createPatch 只是从已经画好的网格上选一些面来生成一个新边界。

@soulx7 参考这个算例：
tutorials/incompressible/pisoFoam/LES/pitzDailyMapped

@wg0632
Gauss LUST grad(U) 当中的 grad(U) 对应一个梯度计算格式，会从 gradSchemes 部分取查找。如果 gradSchemes 部分定义了 grad(U) ，则使用该部分的 grad(U) 对应的格式。如果gradSchemes 部分没有定义grad(U) ，则会使用 gradSchemes 部分的 default 对应的格式。如果 default 也没有则报错。

所以，Gauss LUST grad(U) 和Gauss LUST unlimitedGrad(U) 的区别取决于 gradSchemes 里面如何定义的 grad(U) 和 unlimitedGrad(U)。

参考multiphase\overInterDyMFoam\rigidBodyHull

老哥学明白了吗？

@尚善若水 在 基于reactingFoam的管道内预混氢气点火爆炸 中说：

应该都没有特别考虑这个功能，这个也许可以 https://github.com/JieSun-pku/detonationFoam

对，这个应该可以。

会不会是实验的冷却水产生的实际效果跟你模拟设置的边壁14℃是不一样的？

https://cfd-china.com/topic/6644

你看看这个，我觉得可以用symmetry边界

应该采用laplace 光顺而不是弹簧光顺

@hachikoi 注意到了，似乎在螺旋桨与船连接处存在了一个disk薄片，

有大佬知道如何在V2206版本编译moody系泊插件吗？

@tidedrinker 好嘞好嘞谢谢了，我去安装个ESI版本试试

@Joann 您好，我也在尝试结合 waves2foam 和 overinterdymfoam 。但是当我运行一个案例时，入口处 alpha.water 和U的边界条件需要设置为零梯度才能跑。跑算例的过程会出现计算发散，不知道是不是和inter Foam自带的波浪库冲突了。

@tyhero123 好的谢谢

如果每个网格在同个时刻的平均速度是一样的，更简单的方式是按照C程序写法直接读取这条曲线的离散数据，插值出当前时刻的平均速度，并加到decayingTurbulenceInflowGenerator中的U[faceI]

@wjxhit 这不是网格的问题，是那个位置的粒子正好被注射在了两个cores的中间应该。不知道你具体用的什么注射方式，但只有移除那个粒子才能接着run。

@Jasper-0 你现在能堆积，是因为把youngsModulus系数增大从而增大了碰撞力。这会引起无法堆积的，你可以看你的结果，应该会一直蹦，因为排斥力过大。这个系数因材料而定，不能随便给。比如沙子大概e+06好像。碰撞力的大小需要在计算时调节需要的参数，不能改变粒子的物理属性。