CFD中文网

@李东岳
老师，高尔夫球和离散气泡模型中的那种效果叫什么名字？在paraview之中要如何实现啊？

@tidedrinker
pimpleDyMFoam是新版本快？

在非压缩的求解器中,我查了4.1之下的创建通量的createPhi文件,发现有两个一个里面写的是fvc::flux(U),还有一个是linearInterpolate(hU) & mesh.Sf().说明里面写了前一个是相对relative.这两者有什么区别呢?我試了非压缩下的标量运输,fvc::flux(U)和linearInterpolate(U) & mesh.Sf()(和上面比少了一个h)两个的结果是一样的,请问下这两者是不是同一种东西两个写法呢?

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

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

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

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

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

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

@李东岳
读取每个cell的值

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

不含有扩散项和源项

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

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

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

@youmengtian 我当时一直找不到答案，都没管了。。。。。谢谢你

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

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


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

成功流动，结果一样


rhog和rho

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

@yhdthu 没处理

@mohui
我是在自己写的空化代码中求 Rayleigh–Plesset 方程，得出气相的数据和气体和液体alpha。把原方程求的rho当成混合流的密度来看，再求出液相的密度。

大家好，我一直是在基于rhoPimpleFOAM去做一个空化的求解器。

基本思路是;
自己写的空化求解代码 -> rhoPimpleFOAM的代码 -> 自己定义的各个量在网格的守恒方程

从第一步的空化求解代码得到四个自定义的值，分别是气泡数，气泡质量，气泡半径和气泡半径速度。

气泡质量 单位kg/m3

气泡数 单位1/m3

气泡半径 实际上是半径乘数量 单位 1/m2

气泡半径速度 实际上是半径速度乘数量 单位 1/m2s

以上四个就是自定义量的守恒方程，用以计算在格子中生成后，跟随流动的变化。
其中自定义通量phit1和密度rhot的定义分别为

在creatfieds处定义了rhot，然后在0文件夹中见了一个rhot文件作如下设置：

在creatfields中定义了phit1，我是仿照通量phi去写的

我自己简单的试了一下好像可以，但是在真的运算中气泡数缺出现了负值如图

不知各位有没有类似定义过守恒方程经验的同学能够提供点意见？

@李东岳
级测试，如您所说改为NO_READ后不再出现警告，谢谢

我自己定义了几个通量，在计算时出现了以下的警告。

我按照提示改成了mustread可是还是出现了相同信息。各位同学有遇过相同的情况吗？

@wwzhao
谢谢赵老师的回答，我尝试下写入不同的值看是否能正常写入文档

我最近在试着一个OF里面的一个简单的功能实现。
就是在任意一个求解器中加入一段这样的代码：

我想对网格上的一个值alphal进行写入1.0，然后看下会不会在0之类的时间文件夹之下的文档输出。
alphal的定义如下：

这是边界条件的文档：

我进行计算后输出的却没有像速度文件之类的输出：

还是像原来一样：

我查了log文件都像这样的写入了，可是为什么没有在时间文件里写入呢？

添加了write的代码还是一样：

@影川风
我在研究的是rhoPimpleFoam，虽然不一样，但是在rhopimpleFoam也出现了psi，我也想对其进行修改，
我之前也是找不到修改的地方所以我是直接thermo的代码给//掉，自己写了一个状态方程，
但在pEqn里面求解时还是有psi，觉得自己修改得不完全。

@李东岳
的确奇怪，我现在都不知道为什么纯粹进行值之间的交换会出现问题，又没有显示命名的重复，
不过我把发生错误的变量名字改了以后，错误就不在出现，虽然解决了，但还是对于错误的原因存疑。

我一直都在试，好像错误总出现在cell的值的交换当中。

你好，请问你这个问题解决了吗？