动量方程如何加入加压力梯度源项？

[李东岳]

方案二是正确的，方案一mesh下面不能赋值，可以直接写fvc::grad(p)。

同理，下面那个T，如果直接给gradT....，就不需要写mesh了，直接赋值。

[气固两相流]

@d-benjamin 请问按照方案二加入压力梯度后，在算例中应该如何给压力梯度赋值呢？

[李东岳]

他已经赋值了

pGrad = fvc::grad(p);

[气固两相流]

@东岳

东岳老师，我的错误显示是我没有在0文件夹下建立 pGrad 文件，建立pGrad的文件的话，这个地方应该怎么赋值，而且如果已经在blockMesh中设置周期性边界条件cyclic的话，那么pGrad不就和p一样了依然没有压差么？

[李东岳]

volVectorField pGrad
(
    IOobject
    (
        "pGrad",
        runTime.timeName(),
        mesh,
        //IOobject::MUST_READ,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    Zero
    //vector(1, 0, 0)
);

因为你代码pGrad为必须读取，你做上面的改动。可以把pGrad设置为默认的(0 0 0)。也可以写为vector(1, 0, 0)

[D.Benjamin]

@气固两相流 你好，正如东岳老师所说，只要在0文件夹下有p文件，就会按照给的式子进行赋值了

[气固两相流]

@东岳 非常感谢东岳老师的回答，但是直接对他赋具体值zero或者vector（1 0 0）没有成功，但是fvc::grad( p )是没有问题的

[气固两相流]

@d-benjamin 嗯嗯，谢谢您的回答，直接赋值fvc::grad ( p)是没有问题的

[我是河滩]

你好，你在尝试1中加压力梯度源项，为什么对UEqn.H和pEqn.H都做了修改？
simple算法求NS方程是分两步，UEqn.H是预估步，pEqn.H包括校正步。既然在预估步就加入了源项，为什么在校正步pEqn.H还要添加源项呢？
另外，对于尝试3 fvoptions，这种方法的源项是加在simple求解的哪一步中了？是和尝试2相同吗？

[D.Benjamin]

@我是河滩 你好，他只修改了U的方程，并没有修改p方程，只是定义了一个新的压力梯度场。他上面三个尝试应该都是一样的。有一点我不是很明白，就是尝试3中的volumeMode，两个关键字absolute和specific有什么具体的区别。如果东岳老师可以帮忙解答一下就好了。@东岳

[我是河滩]

@东岳 在 动量方程如何加入加压力梯度源项？ 中说：

多相流的的表面张力项

如何把源项加入到压力方程中，和动量方程加源项有何区别？

[张RP]

@D-Benjamin 请参照 http://www.cfd-china.com/topic/208/fvoptions-rho-yi-与scalarsemiimplicitsource/9

[李东岳]

@我是河滩
加入到压力方程更加无震荡
加入到动量方程更简单

[我是河滩]

谢谢。

[李东岳]

@我是河滩 你做什么方向的啊，

[我是河滩]

运动物体计算，浸入边界。