教授

之所以压力泊松方程离散出来之后得到的矩阵是对称的且不是严格对角占优，是因为通常泊松方程里面的 laplacian 项是用中心差分来离散的。楼主看一下 src/finiteVolume/finiteVolume/laplacianSchemes/gaussLaplacianScheme里面 fvmLaplacian 函数的定义就明白了。
动量方程为什么对角占优呢？主要来自于

1. 对流项，特定的对流项格式可以保证离散之后对角占优
2. 时间项，时间项的离散得到的系数只会进入系数矩阵的主对角。
3. 隐式松弛，从构建系数矩阵的角度来看，隐式松弛和时间项达到的效果一样，都会增加对角占优
4. 特殊的源项，比如 fvm::Sp 得到的隐式源项，也是只贡献于系数矩阵的对角元素，可以增加对角占优（不正确使用则会减弱对角占优）。

非常感谢！已更新

最后是都忽略了，这样就没问题了吧

感谢分享代码！

但是去debug别人的代码很费事...尤其这种算两步就发散的问题

原来这样。

最近在算好多LES算例，直接把服务器当主机用了。日常用没问题。但是发现最大的问题就是word文件的编辑。因为我这面很多都是别人发来的word文件我在上面继续编辑。但是linux下面用各种word编辑软件都会出现格式问题。有些文档就两页，但是导出pdf的时候格式就不对。有些是商务文件，格式对不齐，感觉不是很专业。

我先试试wps，看看会不会好一些。

我有空试一下。

老铁知道linux这面有什么方便处理word的东西么。默认的libreoffice，稍微复杂一点的文档就格式有问题。office online有点卡，格式有时候也出现问题。linux下开虚拟机开word更卡..

meanDiameter
    {
        type            coded;
        libs            ("libutilityFunctionObjects.so");
        name            error;
    
        codeExecute
        #{
            const volScalarField& d =
                mesh().lookupObject<volScalarField>("d.alpha.oil");

            scalar d32 = d.weightedAverage(mesh().V()).value();
    
            if (Pstream::master())
            {
                std::ofstream file;
                file.open ("d32", std::ofstream::out | std::ofstream::app);
                file << mesh().time().timeName() << " " << d32 << "\n";
                file.close();
            }
        #};
    }

我最近也凑巧要用这个，我用的上面这种方式。

你调用volScalarField p， 然后计算一下fvc::grad(p)不就可以了么

给你个思路，用codedFixedValue可以么

@bestucan master pdf这个linux使用方便么 这个okular不能跟adobe那样进行标注 okular的标注很奇葩 搞出来是这种的：

这是模拟泡泡流么？好像跟我这个文章里面玩的东西是一个东西

https://www.jianguoyun.com/p/DfaaRYkQ9s3ZBhjG5fED