buoyantPimpleFoam相关问题

[李东岳]

win7太老了，没计划支持win7了老铁

[yjc95]

@李东岳 东子老师，这几个代码的疑问您抽空看下哈

[李东岳]

7，更新边界条件
6，你说得对
5，更新通量
3、4、调整边界条件，确保边界条件设置的没问题
1、就是HbyA，constrain用于保证HbyA边界条件正确

[yjc95]

@李东岳 红框②中phiHbyA是为了更方便的计算$\sum HbyA\ast S_f$才这样弄的吗？

[李东岳]

fvc::flux(HbyA)就是phiHbyA

[yjc95]

@李东岳 那fvc::flux(HbyA)后面加上fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)是啥作用呀？

[李东岳]

ddtCorr这个是Weller自己提出的防止小时间步长影响时间结果的方法，我明天看看更新到笔记

[yjc95]

@李东岳 老师，我看完icoFoam之后，花了一个周末又重新研究了下buoyantPimpleFoam，又有一些新的疑问。

问题1

方程(20)、(24)、(26)中$U_p^t$前面的系数应该是${\rho }_p^{\ast }/\mathrm{\Delta }t$吧？

问题2

这一部分应该是速度方程吧？

问题3

前文未出现$U_f^{\ast \ast }$和$p^{\ast }$，我觉得应该改成$U_p^{{}^{\prime }}$和$p_{rgh,f}^{{}^{\prime }}$比较好理解吧？

问题4

这里应该是(37)吧?

问题5

红框这俩符号应该是一个物理量吧？

问题6

红框部分应该是$p_{rgh}^{\ast }$吧？

问题7

2->3->4->5->6->7->2应该就是一个时间步内的PISO循环吧？

问题8

步骤5中，方程(10)中只更新$U_f^t$为$U_f^{\ast \ast }$，不更新${\rho }_f^t$为${\rho }_f^{\ast \ast }$吗？

问题9

您在笔记中提到的phiHbyA=$\sum U_f\cdot S_f$，我总觉得这里不对，我觉得应该是phiHbyA=$\sum Hby{A}_f\cdot S_f$。因为在不可压+瞬态这篇算法里面，您说fvc::div(phiHbyA)表示的是$▽\cdot Hby{A}^{t+\mathrm{\Delta }t}$，不然就感觉对不上了。

问题10

pEqn.H文件中const volScalarField psip0(psi*p)，这句代码里面的psip0 和 psi啥意思呀？p我知道是压力。因为后面出现了很多语句都有 psi，不知道的话影响读代码了。

问题11

笔记中这三个红框没看懂，前俩为啥是≈，而第三个却是=了？我觉得三个都应该是=呀。

[yjc95]

@李东岳 东子老师，我真的不是在审稿，我是在认真问问题呀

[李东岳]

今天忙 事多没顾上 明天回你

[李东岳]

问题1 是可压瞬态算法么？没理解

问题2 已更新

问题4 已更新

问题5 是一个物理量

问题6 是修正压力

问题7 是的

问题8 不更新

问题9 笔记里面的截图我看一下

问题10 psi是可压缩性，psi0是更新前的可压缩性

问题11 这三个应该都是约等号，我更新一下

[yjc95]

@李东岳 问题1是可压瞬态算法里面的
问题9的图如下

[yjc95]

@李东岳 东子老师，下面这两张截图是来自可压+瞬态算法，其中：

psi*correction(fvm::ddt(p_rgh)) 对应红框1

fvc::ddt(rho) 对应红框2

那么fvc::div(phi) 对应方程(39)哪一项呢？您在可压+稳态算法中提到，可压求解器中通量phi定义为$\rho U\cdot S_f$，那代入fvc::div(phi) 这句代码中，不应该就是$\rho U\cdot S_f$吗？但方程(39)中跟它很像的只有红框3，只是把$U$替换成了$HbyA$？这样替换的原因是啥呀？

fvm::laplacian(rhorAUf, p_rgh) 我猜是对应红框4，那红框5是被省略了还是合并了？

我发现，我用of8版本，压力泊松方程和您给的不一样，出现了phiHbyA，这也是我在上一个问题问您的。

[李东岳]

方程的3跟5合起来是fvc::div(phiHbyA)

有关fvc::div(phiHbyA)是U还是HbyA，看我那个笔记，昨天更新了一下

[yjc95]

@李东岳 那fvc::div(phi)和of8中的fvc::div(phiHbyA)是代表一个意思，只不过是of版本不同吗？

[李东岳]

是的

[yjc95]

@李东岳 东子老师，根据公式(36)，红色箭头部分应该缺了$1/V_p$吧。

[李东岳]

你说的对，已加上体积，已更新

[yjc95]

@李东岳 老师，这俩没加上体积

[yjc95]

@李东岳 东子老师，我把可压+瞬态算法全部看完了，代码也看完了，打算用它来模拟储热过程。我在of8里面找到了buoyantPimpleFoa下面的hotRoom案例，感觉和我想要的储热过程很相似，所以我就把floor设置成了低温水入口，ceiling设置成了高温水出口，具体的设置参数是这样的：

0文件夹下有5个文件，其中：
nut

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 1e-6;

boundaryField
{
    floor
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
    ceiling
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
    fixedWalls
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
}

p

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1e5;

boundaryField
{
    floor
    {
        type            zeroGradient;
    }

    ceiling
    {
        type            fixedValue;
        value           uniform 0;
    }

    fixedWalls
    {
        type            zeroGradient;
    }
}

p_rgh

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    floor
    {
        type            zeroGradient;
    }

    ceiling
    {
        type            zeroGradient;
    }

    fixedWalls
    {
        type            zeroGradient;
    }
}

T

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 363;

boundaryField
{
    floor
    {
        type            fixedValue;
        value           nonuniform 283;
    }
    ceiling
    {
        type            zeroGradient;
    }
    fixedWalls
    {
        type            zeroGradient;
    }
}

U

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    floor
    {
        type            fixedValue;
        value           uniform(0 0.1 0);
    }

    ceiling
    {
        type            zeroGradient;
    }

    fixedWalls
    {
        type            noSlip;
    }
}

constant文件夹下面修改了momentum Transport，pRef和thermophysicalProperties：
momentum Transport

simulationType  laminar;

pRef

dimensions      [1 -1 -2 0 0 0 0];
value           1e5;

thermophysicalProperties

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState Boussinesq;
    specie          specie;
    energy          sensibleEnthalpy;
}

pRef            100000;

mixture
{
    specie
    {
        molWeight       18;
    }
    equationOfState
    {
        rho0          1000;
        T0             273;
        beta          3e-05;
    
    }
    thermodynamics
    {
        Cp             4182;
        Hf              0;
    }
    transport
    {
        mu              1e-03;
        Pr              7;
    }
}

system文件夹下就修改了fvSolution文件，增加了pRefCell，pRefPoint，pRefValue：

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors 1;
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    pRefCell       0;
    pRefPoint      0;
    pRefValue      0;
}

做完上述修改之后，输入buoyantPimpleFoa，报错了：

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0



--> FOAM FATAL IO ERROR: 
keyword version is undefined in dictionary "/home/dyfluid/OpenFOAM/dyfluid-11/run/hotRoom/constant/polyMesh/points"

file: /home/dyfluid/OpenFOAM/dyfluid-11/run/hotRoom/constant/polyMesh/points from line 10 to line 13.

    From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 799.

FOAM exiting

我觉得可能是边界条件没设置好，检查再三感觉也没啥问题，麻烦您抽空看看，谢谢您！