CFD中文网

@李东岳 这个案例我用的商软，发现空隙率0.01和0.5差别不大，影响最大的是阻力系数，我试了这两个空隙率，水都能很顺畅的流出来

同问，目前在用PBM处理气液两相流，用discrete方法求解PBM，也想试试qmom方法，上面几个矩代表什么意义？要模拟气液流中的气泡，该怎样设置粒径，数量还是。。。？

@wwq 在 请教关于k和epsilon壁面函数的问题 中说：

https://www.openfoam.com/documentation/guides/latest/doc/guide-bcs-wall-turbulence-epsilonWallFunction.html

非常感谢:146:

@李东岳 谢谢李老师！A fifth-order shock capturing scheme.pdf 李老师，我是想把这个文献里的第二部分模块数值算法植入进去，我昨天上传的是第二部分前面一块，后面那一块主要是进行重构，去计算q_L和q_R，然后通过前面的公式去进行计算，前面那块我看和您发的rhocentralfoam解析中很像，但后面的公式不知道应该怎么植入，没有思路，想请教一下老师:chouchou:

@李东岳 解决了，要用全局代理

对对对，这样是对的，我之前把应变率和湍动能搞混了，谢谢李老师了！

@东岳 老师，我重新画网格后，decomposePar经常会报错，但是不影响后面计算，没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

https://stackoverflow.com/questions/30304794/gnuplot-command-with-multiple-with-arguments

@袁宝强 我记得可压缩求解器里面有个可压缩系数，psi，不知道是否有用

@李东岳 因为换个地方继续算，重新分块的时候就根据最后时刻的网格来划分，所以用scotch 就是每个块基本同样网格数（上下几十浮动），而用simple划就会是那种XYZ方向的拉升

@李东岳 谢谢李哥

能，模型是否适用和这个没太大关系，和状态有关系。
如果流动特别慢都没有气泡，就不适用了。

lz解决了吗？我也遇到同样的问题了，高压气流冲击真空环境，使用rhoCentralFoam，算了一段时间之后就出现negative initial temperature t0，改了好久，可以推迟出现的时间，但总会出现，不知道怎么弄了

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。

密度就是单位体积的质量，速度的物质导数就是加速度。
所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？

力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。

物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。

要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。

另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。

@Hungryandfool
Operator 翻译过来叫做 ”算子“，也就是说，微分方程的空间项可以看作是不同的算子，比如 “对流”（散度），“扩散”（梯度）；
Splitting 的意义在于不同项可以相对独立的去求解，甚至采用完全不同的求解器，这样有利于提高计算效率，减小计算复杂性

@东岳 内网格，尝试着接着划吧

如果我没记错的话，流体中文网上的ICEM版块里有人画过4119桨的结构网格，也传了文件在上面，你可以去上面找一下看还在不在。（耐心往后翻一下）

@李东岳 LaTeX党遇到这个真的心累。只好把它留到最后一口气吃掉。

@李东岳 老师，内部也在动的，截图可能看不出来，只是这个边界更明显。
75376c62-5086-4608-879f-b7823f7ebaa8-1639391378.png
7a6ca8ef-4d95-4e82-b6e1-77c1203adfab-1639391406(1).png