buoyantPimpleFoam相关问题

[yjc95]

@yjc95 东子老师，您有空看看哈

[Afr1yne]

@yjc95 我也是新手试着回答一下，如果有错误请大佬们纠正下
建议你去读下源代码 你的问题都在里面
1.internalField是初始场 type就是个壁面函数 你可以查下源代码看看公式
2.是把rhok赋值给rho value是针对type的跟rho无关
3.层流alphaEff = nu/Pr

[李东岳]

是，楼上说的对。

另外你说的那个内部场都是0的问题，属于openfoam/CFD最初级的内容，这个就是CFD的初始场设置。类似的初级问题我这面回复处于最低优先级。

[yjc95]

@李东岳 东子老师说的对，我基础太差了，我该好好去读您网站上给出的书籍。

[yjc95]

@Afr1yne 大佬，追问一个问题，问题3中的湍流Pr该咋取呀？

[Afr1yne]

@yjc95 我也是新手 不是大佬
案例中一般都取0.7 具体的你可以找找文献看

[yjc95]

@Afr1yne 谢谢您

[yjc95]

@李东岳 东子老师，新年好呀。我带着新问题来向您请教了。我用v2206版本下的heatTransfer-buoyantBoussinesqPimpleFoam-hotRoom作为模板来模拟储热池的放热过程，我把floor设置成低温水进口，把ceiling设置成高温水出口，采用层流，介质为水，算出来的结果温度场很奇怪，中间的斜温层很厚。

我查阅文献，并且用fluent算出来的斜温层应该非常薄才对。
此外，我把所有物性参数扩大和缩小1000倍，温度场没变化。
具体设置如下：

• alphat

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    floor
    {
        type            alphatJayatillekeWallFunction;
        Prt             0.1;
        value           uniform 0;
    }

    ceiling
    {
        type            alphatJayatillekeWallFunction;
        Prt             0.1;
        value           uniform 0;
    }

    fixedWalls
    {
        type            alphatJayatillekeWallFunction;
        Prt             0.1;
        value           uniform 0;
    }
}

• nut

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

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

    ceiling
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    fixedWalls
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
}

• p

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 1e5;

boundaryField
{
    floor
    {
        type            zeroGradient;
    }

    ceiling
    {
        type            zeroGradient;
    }

    fixedWalls
    {
        type            zeroGradient;
    }
}

• p_rgh

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 1e5;

boundaryField
{
    floor
    {
        type            zeroGradient;
        rho             rhok;
    }

    ceiling
    {
        type            zeroGradient;
        rho             rhok;
    }

    fixedWalls
    {
        type            zeroGradient;
        rho             rhok;
    }
}

• T

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 363;

boundaryField
{
    floor
    {
        type            fixedValue;
        value           uniform 283;
    }

    ceiling
    {
        type            zeroGradient;
    }

    fixedWalls
    {
        type            zeroGradient;
    }
}

• U

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      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            fixedValue;
        value           uniform (0 0.1 0);
    }

    fixedWalls
    {
        type            noSlip;
    }
}

• transportProperties

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel Newtonian;

// Laminar viscosity
nu              1e-06;

// Thermal expansion coefficient
beta            3e-04;

// Reference temperature
TRef            273;

// Laminar Prandtl number
Pr              7;

// Turbulent Prandtl number
Prt             0.1;

• turbulenceProperties

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType      laminar;
/*
RAS
{
    RASModel        kEpsilon;

    turbulence      on;

    printCoeffs     on;
}
*/

麻烦您抽空看看，谢谢您！

[李东岳]

把你的openfoam跟fluent的对比图放上来我看看

[yjc95]

@李东岳 东子老师，我又检查了一遍fluent的设置，发现是我把速度入口数值设置错了，重新算了一遍，看上去好像也差不多。耽误您时间了，抱歉。

[yjc95]

@李东岳 东子老师，您在可压+瞬态里面写到"方程(20)构成一个稀疏线性系统。求解方程(20)即可获得速度 "其中的 指的是相邻网格的速度，在这里是用当前时刻已知的速度代入么？同样的，方程(41)中的 应该是要求的控制体体心处的压力吧， 是否是用当前时刻各个面的压力代入，这个压力好像之前并不知道吧，是咋求出的呢？

[yjc95]

@李东岳 东子老师，我悟了，这俩问题我想明白了。我有新的疑问：方程(41)中的 和 有啥联系呢？

[李东岳]

带f是定义在面上的

[yjc95]

@李东岳 您的意思是，体心处的p是面心处p的某种函数，比如体心处的p＝面心处的p的加权平均？

[李东岳]

面心的值是体心的值的插值

[yjc95]

@李东岳 谢谢您