自己写的一个空化求解器的气泡数，气泡质量，气泡半径和气泡半径速度的守恒方程的一些疑问

[100yearsalone]

@李东岳
你好，李老师
我目前就想试试气泡数量密度的传输方程，我试了几个进行了测试。

你的回答公式无法正式显示，所以我自己查了一下这样定义了。


其中ng2是气泡数量密度，phis2是通量，我感觉应该是没问题的，但是实际运行后结果是
ng2没有发生移动。
我试了一下另外两种定义rg0和rg1，phis0=rhog0U，phis1=rhoU

成功流动，结果一样


rhog和rho

rho和rhog不一样，结果是一样，代表我以后传输方程只用随便定义一个密度是吗？如果不是的话，为什么一开始的式子没发生流动呢？

[李东岳]

你开始的代码

fvm::ddt() + fvc::div(phis2)

里面你只是定义了phis2，并没有更新phis2，因此这个值每个时间步不变所以没有移动。你需要在每个时间步里面计算一下:

while(TIMELOOP)
{
    phis2 = phi*fvc::interpolate(s2);
}

[100yearsalone]

@李东岳
我按照你给我那篇论文里的式子

不含有扩散项和源项

其中mg21就是气泡数量密度，上面的phis2是通量

气泡个数发生的测试条件是这样，我觉得我设置得没有什么问题，但是有一点不太明白

就是在结果里面为什么会出现负值？

[李东岳]

forAll是做什么的？

[100yearsalone]

@李东岳
读取每个cell的值

[李东岳]

结果越界了，用迎风试试

[100yearsalone]

@李东岳
老师，您之前给我写了一个气泡数量密度的运输方程是图中的第一个式子

对于有些标量，例如a我想跟随这个N做移动，我设另外一个标量N乘以a等于Na做同一样的运输方程，如图中式子2
求完N和Na以后，再相除求出a，这种方法在OF里面应该是没问题的吧？
我试了一下，好像不太符合我的预期。
我设了ng0作为第一个式子的N，得出下图：

而rg0作为a，a为常值1，ngrd的结果：

我很不解为什么ngrd的运输方程计算后会出现这种结果？

[李东岳]

例如a我想跟随这个N做移动，我设另外一个标量N乘以a等于Na做同一样的运输方程，如图中式子2

不理解你这句话。另外你上面三个图实际上是完全一样的方程。为什么要求解三个？

ngrd是什么？

[100yearsalone]

@李东岳
ngrd就是两个标量的乘积，第一个式子求了ng，第二个式子求ngrd，最后相除得rd
不过听您这么一说我可能又想错了

这图上面ng0是气泡数，mg0是气泡质量，rg0是气泡半径，rhog0是根据这三个计算的密度。
在运输之前只有一个格子有值，运行运输方程后变成三个格子有值了，可是从密度来看另外两个格子都变成了极大的数，
我就是找个方法比较稳定的传输这几个值，让他们运输以后不偏差密度值太多。

[李东岳]

提示一点：如果求解本身不符合物理的传输方程，结果肯定是不复合物理。

[LiuZhuang]

@100yearsalone 新写的phit1里的U，这个U是从0文件夹下的U调用的么？