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

[李东岳]

也就是说OpenFOAM怎么知道pGrad就是▽p呢？

不知道，在声明pGrad之后，需要赋值，比如：

pGrad = fvc::grad(p);

这样才可以。楼主那个应该没写这部分

U ((10 0 0) 0);

这不是简单的速度矢量了，这是给injectionRateSuSp关键词赋值

[D.Benjamin]

@东岳 老师您好，谢谢你的回复，我还想追问一下，赋值的这部分代码应该写到哪个地方？下面我写了两个方案，哪一个方案是正确的？还有，我发现在laplacianFoam的解析中，有两种定义场的方式，如下面的图片所示，定义场"T"的时候是mesh，而定义场"gradTx"的时候，mesh被替换成了gradT.component(vector::x)，正是因为对此疑惑，我才有上面两种方案是否合理的疑惑，望老师给与答复，感激不尽！

[李东岳]

方案二是正确的，方案一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

[李东岳]

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

[我是河滩]

谢谢。

[李东岳]

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

[我是河滩]

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