CFD中文网

确实和概率论有关，谢谢，记录一下

\begin{equation}
\left( \begin{matrix}
k\\
x
\end{matrix}\right )=\frac{k!}{(k-x)!x!}
\end{equation}

@陈宏洋 至少得有目录页才能制作吧

信息太少了

@gooseEast 在 rhoCentralFoam发散，出现Maximum number of iterations exceeded 中说：

在thermoI.H文件第46行标量f是什么

OpenFOAM代码中，这个f是Cp，定压比热容，单位是[J/(kg K)]。这个T函数是通过定压比热容的值迭代求出温度。limit，F，dFdt分别对应的是hConstThermoI.H文件里limit，Cp，dCpdt三个函数。

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

自己做了一个300m/s的来流，然后打在一个子弹形状的壁面上，想观察来流经子弹产生空化现象。但是现在问题是，300m/s的来流经过子弹形状的壁面输出压力云图是全蓝的，不知道是哪里设置错误，有没有大佬帮帮我。

我使用CFDEM中的解析CFD-DEM模型+RANS+k epsilon计算水平圆管中颗粒沉降的时候(管径115mm，入口流速3m/s，粘度0.01Pa·s，密度1000kg/m3；颗粒密度1100kg/m3，直径15mm；雷诺数几万)，发现颗粒越靠近圆管底部，下沉速度越慢，甚至在距离圆管中心0.4mm的时候，颗粒没有下沉，反而向上移动了。想请教一下各位老师，这种情况是不是不对呢？
印象中层流里面，管壁附近重颗粒才会往管道中间移动，不知道我这个认识对不对。

CFD-online 用户访问数量

2017年

替代文字

2016年

替代文字

2015年

替代文字

2014年

替代文字

2013年

替代文字

2012年

替代文字

2011年

替代文字

2010年

替代文字

2009年

替代文字

2008年

替代文字

@wzlr2023 我看你这个算例里只有两相流，正常用 interFoam足够了。这个相关资料还挺多的。这个多相流的求解器正如我前面所说我没有去验证其是否可以正常求解。就你这个算例来讲，里面涉及到了六自由度运动求解。造成两相交界面求解发散的问题也可能是六自由度求解引起的流场发散。建议可以先设置一个简单的动网格算例，如动边界造波来验证多相流求解是否稳定。如发现是多相流求解出问题了，欢迎进行补充修正。

网格质量如何？

fluid@zhoushan:~/$WM_PROJECT_USER_DIR/FluidSolidInteraction/src$ ./Allwmake

wmake libso fluidSolidInteraction
./Allwmake: 4: wmake: not found wmake solvers/fsiFoam
./Allwmake: 6: wmake: not found wmake solvers/ampFsiFoam
./Allwmake: 7: wmake: not found wmake solvers/weakFsiFoam
./Allwmake: 8: wmake: not found wmake solvers/fluidFoam
./Allwmake: 9: wmake: not found wmake solvers/solidFoam
./Allwmake: 10: wmake: not found wmake solvers/thermalSolidFoam
./Allwmake: 11: wmake: not found wmake libso utilities/functionObjects/pointHistory
./Allwmake: 13: wmake: not found wmake libso utilities/functionObjects/patchAvgTractionHistory
./Allwmake: 14: wmake: not found wmake libso utilities/functionObjects/centrifugalBodyForce
./Allwmake: 15: wmake: not found
请问上述的not found怎么解决呀 在网上参考了各种方法 都无法成功

可不可以结合codestream实现旋转呢？

请教各位，该怎么对边界附近的单元赋值，比如下面这个UDF.我想将一个区域的 C_UDMI(c, tt, 0)全部设置为0，然后边界附近的 C_UDMI(c, tt, 0)设置为1，该怎么操作？谢谢！

void set_zone(Thread* t) { cell_t c; real vol; Thread *tt = THREAD_SUPER_THREAD(t); begin_c_loop(c, t) { if (THREAD_ID(t) == 1047) { C_UDMI(c, tt, 0) = 1; } end_c_loop(c, lct) } void set_cells(Thread *t) { Thread *tt = THREAD_SUPER_THREAD(t); face_t f; Thread *t_nbr; cell_t c0, c1; Domain *mix_domain = Get_Domain(1); for (int i = 0; i < num_boundaries; i++) { t_nbr = Lookup_Thread(mix_domain, 1); begin_f_loop(f, t_nbr){ c0 = F_C0(f, t_nbr); c1 = F_C1(f, t_nbr); if (c0 != -1 && THREAD_ID(THREAD_T0(t_nbr)) == 1047) { C_UDMI(c0, tt, 1) = 0.0; } if (c1 != -1 && THREAD_ID(THREAD_T1(t_nbr)) == 1047) { C_UDMI(c1, tt, 1) = 0.0; } } end_f_loop(f, t_nbr) } } DEFINE_ADJUST(aaa, mix_domain) { Thread* mix_thread; int phase_domain_index = 0; thread_loop_c(mix_thread, mix_domain) { Thread* t = THREAD_SUB_THREAD(mix_thread, phase_domain_index); if (t != NULL) { set_capillary_pressure_sources(t); set_boundary_cells(t); } } }

你这一会1250 一会1350 到底多少度的，而且你这个截面是截的多高处的呢，另外得看一下你的残差图

这个问题我也遇到过，比如对于同一个流场，同样的空化数，来流速度10m/s就不会产生负压，来流速度20m/s就有负压产生。增加SS空化模型的系数n的值或许能解决问题。

今天偶尔看到一个文章。好像意思是通过Optuna + OpenFOAM，可以实现最小的计算时间。感觉类似通过Optuna跑好几次，然后寻找出来最合适的松弛因子。不知道有没有哪个大佬感兴趣玩一玩。

https://ss1.xrea.com/penguinitis.g1.xrea.com/study/OpenFOAM/relax_factor_optimization.html

@李东岳 车灯热分析。

@promisegwj 啊不对 是用拉梅系数···

@李东岳 问题已经解决了，是我在录制脚本的时候进行了多余的不必要操作，再次感谢李东岳老师，您的建议大大的减少了我的工作量。

@Samuel-Tu 请问你解决了这个问题了吗？我是在controlDict中添加code进行计算得到的一些自定义的变量，但是结果在写入到文件中时，除了时间步数之外，其他的参数就一直是6位数或更少，小数点后边的就更少了。也尝试改了李老师说的这个，但是好像没用。