@李东岳 在 icoFoam中压力泊松方程的求解 中说:
@AppleKiller 我更新后告诉你
好的,老师,我看了看源代码,感觉fvc::div(phiHbyA)
和fvc::div(HbyA)
没什么区别(可能我理解不深,没找到其他的区别),phiHbyA就是$HbyA_f\cdot\pmb{S}_f$,如果说为什么OpenFOAM使用fvc::div(phiHbyA),我觉得使用这个省掉了一些计算(HbyA插值到面 + 和$\pmb{S}_f$做点积的计算)
@李东岳 在 icoFoam中压力泊松方程的求解 中说:
@AppleKiller 我更新后告诉你
好的,老师,我看了看源代码,感觉fvc::div(phiHbyA)
和fvc::div(HbyA)
没什么区别(可能我理解不深,没找到其他的区别),phiHbyA就是$HbyA_f\cdot\pmb{S}_f$,如果说为什么OpenFOAM使用fvc::div(phiHbyA),我觉得使用这个省掉了一些计算(HbyA插值到面 + 和$\pmb{S}_f$做点积的计算)
@李东岳
那么这两种形式有区别吗?有的话是什么区别呢?
fvc::div(phiHbyA)
fvc::div(HbyA)
@AppleKiller 在 在OpenFOAM中如何判断一个单元是否为边界单元? 中说:
各位道友好!
关于这个问题,我最初的想法是遍历该单元的所有面,然后判断这些面是否在边界上。请问有没有其他的方法呢?(任何方法都行)
这是这个方法的代码,有没有其他比较快速的方法呢?
forAll (mesh().cells()[cellID], faceI)
{
const label patchi = mesh().boundaryMesh().which(mesh().cells()[cellID][faceI]);
if
(
patchi > 0
&& !isA<processorPolyPatch>(mesh().boundaryMesh()[patchi])
)
{
Info << "Find it" << endl;
}
}
各位道友好!
关于这个问题,我最初的想法是遍历该单元的所有面,然后判断这些面是否在边界上。请问有没有其他的方法呢?(任何方法都行)
好的,谢谢李老师!! 这个问题我纠结了好几天
”我跑出来的气泡数量比较少“ 是这个意思:在设置中气泡注入的数量是一定的,但是气泡在运动中发生碰撞融合,导致气泡数量变少。
但我感觉我的气泡数量和文章相比看起来少了一些,所以想问问老师我这个结果有问题吗?从哪些地方来判定我这个算例的结果好不好呢?
@李东岳 在 鼓泡床参数设置的一个问题 中说:
0.25m/s可以 差不多就是0.3 0.4左右
李老师您好,之前关于气泡注入的问题解决了,目前也使用自己的E-L模型来跑这个算例(第一条帖子),但是我跑出来的气泡数量比较少,我的图如下(截的20s~30s之间的三张图,原来要跑150s)
我不知道我跑出来的记过对不对?因为我的气泡数量和文章中的看起来少了很多,气泡数量是判别算例跑的好不好的一个标准吗?如何判断这个算例跑的好不好呢?
谢谢大佬,我的Courant数设置的是0.5,时间步长使用的自适应时间步长。我的想法是可不可能是因为alpha.water场某个单元计算出来比较小,导致得到该单元的U特别大,然后求解压力方程后在该单元附近出现负压?
@bestucan 在 不可压缩多相流计算出现负压 中说:
可以从oepnfoam自带的dambreak算例开始,把里面的模型、速度等等逐步改成你的状况。看看改到哪一步崩了。
谢谢老哥的回复,是不是可能因为界面处的网格不够细?
各位大佬好!在MPPICInterFoam的源代码中没看到关于动网格的源文件,所以MPPICInterFoam求解器支持自适应网格技术吗?
各位大佬好,本人最近在算多相流的算例,算到某一时刻在气液界面处出现负压(使用的是绝对压强),通过输出的数据,大概是某单元的速度场过大(这个可能是alpha.water场在该单元的值过小导致的),导致后面计算的压强过大,从而出现负压。请问各位大佬有没有遇到过这种情况?这个问题应该怎么解决呢?
以下几张图是alpha.water,U,p场的截图(图片贴出的是有错误的部分)
这个模型能用于Euler-Lagrange模型吗?
@李东岳 在 MULES::correct() 和 MULES::explicitSolve() 有什么区别? 中说:
MULES要看基金会的版本,这是基金会那面提出的。我记得CMULS已经被删了。我在无痛苦ns方程笔记里面说的很齐全。你看看参考一下
好的,谢谢李老师
各位大佬,我看的是OpenFoam-v2206版本的MULES方法,粗看代码,感觉MULES::correct(...)函数(源代码位于CMULESTemplates.C)和MULES::explicitSolve(...)(源代码位于MULESTemplates.C)的实现很像?想请问各位这两个函数有什么区别吗?(稳定性、精度、实现格式等各个方面都行)
@李东岳 在 bubble colume算例一个 中说:
相关描述请参考 http://dyfluid.com/openfoam.html
李老师,您好,这个算例是运行到一定时间后才会出现周期性现象吗?
@液固两相流 在 openFOAM中颗粒压力梯度力于虚拟质量力的计算问题 中说:
压力梯度力的计算公式如下图1所示:
图1
虚拟质量力的计算公式如下图2所示:
图2
在openFOAM中的计算代码如下所示:
压力梯度力:
value.Su() = mass*td.rhoc()/p.rho()*DUcDt;
虚拟质量里:
value.Su() = Cvm*PressureGradient
OpenFOAM里面的计算公式与实际的计算公式有一定的差别,请教一下大家这种处理是基于什么考虑呢?
大佬,您最后是怎么处理的呢?我也意识到这个问题,如果要计算虚拟质量力的话,我发现粒子速度的全导数不易计算(我模拟的是气泡在水中的上浮)
@pengdi 前辈还有个小问题想问问您
您当时求解RP方程后(计算出当前气泡的半径)
是不是需要重新求解平衡状态下的R0和P0?还是直接使用第一次计算得到R0和P0(求解方程如下)
因为方程(1.12)中的PB需要用到R0,所以我比较疑惑,想问问您
表观气速和柱体截面积之积是不是就是体积流量呢?
老师您一般气泡的初始速度设置多少呢?
@李东岳 在 鼓泡床参数设置的一个问题 中说:
是的,你算的是正确的
老师,还有个小问题,对于气泡直径为4mm,液体为水,气体为空气,在这种情况下气泡的速度大概是多少呢?(我模拟气泡大概是0.25m/s左右,我不知道对不对)