CFD中文网

另外我在youtube上找到一个较为详细的教程，需要的朋友可以看看。
https://www.youtube.com/watch?v=jJHOeDZJUrU

@李东岳 感谢老师指点，确实是个很巧妙的方法。

@李东岳 老师，实在打扰，是否有时间帮我看看这个小问题？

前期在of v1812版本中开发了某求解器，现在发现该版本的EOS中并没有rPolynominal，只有以下几类（看到的rPolynominal是我后期添加的），

~/OpenFOAM/v1812/OpenFOAM-v1812/src/thermophysicalModels/specie/equationOfState$ ls
Boussinesq       adiabaticPerfectFluid  incompressiblePerfectGas  perfectFluid  rPolynomial
PengRobinsonGas  icoPolynomial          linear                    perfectGas    rhoConst

于是我从v2206中复制了rPolynominal，当我准备编译时，却没有找到Make文件，在上一级文件夹~/OpenFOAM/v1812/OpenFOAM-v1812/src/thermophysicalModels/specie中有Make/file如下，

atomicWeights/atomicWeights.C
specie/specie.C
reaction/reactions/makeReactions.C
reaction/reactions/makeLangmuirHinshelwoodReactions.C

LIB = $(FOAM_LIBBIN)/libspecie

这些好像并不涉及EOS的编译，所以我不明白要添加一个EOS应该怎么操作呢？有没有大佬指教一下

前辈你好，我最经也遇到这个问题，请问有时间交流一下吗？

@李东岳 好的，谢谢老师！

@李东岳 老师，我接触openfoam没多久，理解还没有那么深入。我初步是这样想：首先去读取csv文件，生成一个全域的foamfield场文件；然后在createPorousZones.H文件里加载阻滞系数场，并赋到对应的区块上；最后动量方程里有关多孔介质的源项也改成那个全域的阻滞系数场。老师你看看按这样的思路做有没有问题

@李东岳 老师，是否有时间给点建议？或者openfoam现在有自带的赋值非均匀多孔介质的方法吗？

这是porosityProperties文件和toposet文件有关区块porosity1的内容，

porosity1  //烟丝
{
    type            DarcyForchheimer;
    cellZone        porousz1;
    D1 1.17e9;
    d   ($D1 $D1 $D1);
    f   (0 0 0);
    coordinateSystem
    {
        type     cartesian;
        origin   (0 0 0);
        rotation
        {
            type     axes;
            e1       (1 0 0);
            e2       (0 1 0);
        }
    }
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
    {
        name    porous1;
        type    cellSet;
        action  new;
        source  cylinderToCell;
        sourceInfo
        {
            p1      ( 0 0 0);
            p2      ( 0 0 0.062);
            radius  0.00385;
        }
    }
    {
        name     porousz1;
        type     cellZoneSet;
        action   new;
        source   setToCellZone;
        sourceInfo
        {
            set  porous1;
        }
    }

可见区块1是一个半径0.00385，长0.062的圆柱体，整个区域的阻滞系数都是1.17e9，现在我想将阻滞系数设置为沿轴向z非均匀分布的形式，可以通过外部文件（如txt，csv等）读取阻滞系数参数，然后赋值在轴向的不同位置上。原来的DarcyForchheimer是不是并不支持这种赋值方式，是不是需要修改源代码才能实现呢？