关于多相流的电势方程（Laplace方程）的求解问题

李东岳

@cfd_lilili 在关于多相流的电势方程（Laplace方程）的求解问题中说：

fe4.0

这是啥。。。

cfd_lilili

@东岳 openFoam-extend4.0，不好意思，我以为这个是通用的简写

李东岳

@cfd_lilili 在关于多相流的电势方程（Laplace方程）的求解问题中说：

smoothSolver: Solving for phiE, Initial residual = 1, Final residual = 0.00666667, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.969654, Final residual = 0.955898, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999906, Final residual = 0.999906, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999943, Final residual = 0.999943, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999959, Final residual = 0.999959, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.99997, Final residual = 0.99997, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999977, Final residual = 0.999977, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999982, Final residual = 0.999982, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999986, Final residual = 0.999986, No Iterations 1000
smoothSolver: Solving for phiE, Initial residual = 0.999989, Final residual = 0.999989, No Iterations 1000

你在求解一个方程求解多次。另外，迭代次数1000，没有收敛，

这是一个非常简单的方程，你需要更细心的debug

cfd_lilili

@东岳谢谢老师
我把代码换了写法，问题就解决了
这是原来的写法

        fvScalarMatrix phiEEqn
        (
            fvm::laplacian(  phiE )
        );

        phiEEqn.relax();
        solve(phiEEqn);

把它简单地写为

    solve( fvm::laplacian(sigmaE, phiE) );

问题勉强算解决了，但还不知道之前的问题出在哪。

李东岳

后来这个解决了吗 :-)

gzs20

老师，请教一下，不同相的交界面电势的边界条件采用的那种了

李东岳

但其电导系数（sigma）在不同的相中具有不同的值。

我去年在笔记的附录中整理了这部分内容。这部分采取多计算域的方法处理比较好。

gzs20

@李东岳东岳老师，您好，请问您说的笔记的附录在哪里可以找到，告知一下，目前卡在了电势的求解这一步

李东岳

http://dyfluid.com/index.html 主页有下载

gzs20

@李东岳谢谢东岳老师

CFD中文网

关于多相流的电势方程（Laplace方程）的求解问题