administrators

网格数量太多了。15510 * 16000 * 22个网格？

他那个sigmaDotU是已经离散之后的了，你看他成了面矢量，类似$\nabla\cdot\bfU$，里面$\bfU$是体矢量，但是离散后的形式是$\sum \bfU_f\cdot\bfS_f$，里面$\bfU_f\cdot\bfS_f$就是标量，但是是面标量。

$\nabla\cdot(\tau\cdot\bfU)$中$\tau\cdot\bfU$是体矢量，离散后是$\sum (\tau\cdot\bfU)_f\cdot\bfS_f$，$(\tau\cdot\bfU)_f\cdot\bfS_f$为面标量，这个就是sigmaDotU

τ·u这个是矢量，所以你的sigmaDotU应该是矢量，

看起来可以通过更改源代码的方式，定义一个cellZone，然后再cellZone里面去找。貌似是个路子。

@strive_dut 在 OpenFOAM后处理怎么提取某点随时间变化的数据 中说：

计算过程中添加probe导出的数据存放在了输出的时间步文件夹中

这个不是输出在postProcessing文件夹中么？

应该影响不大。这种颗粒应该不多，不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话，那就会引起收敛性的问题。

@冠竹 在 OpenFOAM setFields 无气泡求助 中说：

相场初始化不成功的原因可能是哪些？

1. 你看下是否真的没有初始化成功，手动打开下alpha文件，如果文件很大就是成功了。paraview有时候不显示

2. 另外就是位置不正确。你把你setFields里面设置的值，去paraview里面画一个sphere，看是不是在你要的位置

目前我遇到的100%是这两个原因。

@hotairballon1997 p.cell()是拉格朗日粒子所在单元的标号，那就是你说的只读取网格中心值(cell方式)

第一次是1没问题。0.1也没问题。这非常正常。不需要优化。

@疏影横斜水清浅 类似这种只要有数学公式，植入都不是问题。所以公式是什么样的？