有关多孔介质求解器中阻滞系数不均匀该怎么处理

unravelkaka

这是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是不是并不支持这种赋值方式，是不是需要修改源代码才能实现呢？

unravelkaka

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

李东岳

你这个应该改源代码了。不过很简单。

unravelkaka

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

李东岳

你这个最简单的是不去修改porous库。而是直接修改求解器顶层。手动写一个多孔介质源项过去。

或者你参考这个，也可以写 https://cfd-china.com/topic/7607/设置vectorcodedsource类型源项-源项量纲和程序量纲不匹配/5

unravelkaka

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

CFD中文网

有关多孔介质求解器中阻滞系数不均匀该怎么处理