[boffin5]

Hi,
My focus here is simulations with porous media defined by Darcy Forchheimer coefficients. There is a simpleCar tutorial in v2012; it simulates the airflow through a car radiator. When I run it there it gives reasonable results. Then I adapted it to run in OF8, which I prefer to use. In this case, it doesn't recognize the presence of the radiator. Also, under polymesh there is a 'sets' folder that doesn't appear in v2012. I'm not sure why.
I need to understand why the OF8 version isn't running properly. Attached are images of the paraview outputs. The files that should pertain to the problem are fvOptions and topoSetDict. Here they are for v2012 and OF8, respectively:


Hoping for advice, solutions and/or information!


fvOptionsv2012

porosity1
{
    type            explicitPorositySource;
    active          true;

    explicitPorositySourceCoeffs
    {
        type            DarcyForchheimer;
        selectionMode   cellZone;
        cellZone        porousZone;

        DarcyForchheimerCoeffs
        {
            //d   d [0 -2 0 0 0 0 0] (5e7 -1000 -1000);
            d   d [0 -2 0 0 0 0 0] (5e7 -5000 -5000);
            f   f [0 -1 0 0 0 0 0] (0 0 0);

            coordinateSystem
            {
                origin  (0 0 0);
                e1      (1 0 0);
                e2      (0 1 0);
            }
        }
    }
}

fvOptionsOF8

porosity1
{
    type            explicitPorositySource;
    active          true;

    explicitPorositySourceCoeffs
    {
        //type            DarcyForchheimer;
        selectionMode   cellZone;
        cellZone        porousZone;        
        
        type            DarcyForchheimer;
        
        // D 100;  // Very little blockage
        // D 200;  // Some blockage but steady flow
        // D 500;  // Slight waviness in the far wake
        // D 1000; // Fully shedding behavior
        //DarcyForchheimerCoeffs
        
        
            d   d [0 -2 0 0 0 0 0] (5e7 -5000 -5000);
            f   f [0 -1 0 0 0 0 0] (0 0 0);
            
            //d    (50 -50 -50);
            //f    (0 0 0);

            coordinateSystem
            {
                type    cartesian;    //added
                origin  (0 0 0);      //added
                coordinateRotation  {    //added
                //origin  (0 0 0);
                type    axesRotation;
                e1      (1 0 0);
                e2      (0 1 0);
                }
            }
      }
}

topoSetDictv2012

actions
(
    {
        name    porousCells;
        type    cellSet;
        action  new;
        source  boxToCell;
        box     (2.05 0.4 -1) (2.1 0.85 1);
    }
    {
        name    porousZone;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        set     porousCells;
    }

    {
        name    airIntake;
        type    faceSet;
        action  new;
        source  patchToFace;
        patch   body;
    }
    {
        name    airIntake;
        type    faceSet;
        action  subset;
        source  boxToFace;
        box     (2.6 0.75 0)(2.64 0.8 0.1);
    }


);

topoSetDictOF8

actions
(
    {
        name    porousCells;
        type    cellSet;
        action  new;
        source  boxToCell;
        sourceInfo {
            box     (2.05 0.4 -1) (2.1 0.85 1);
            }
    }
    {
        name    porousZone;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        sourceInfo {
            cellSet     porousCells;
            }
    }

   {
        name    airIntake;
        type    faceSet;
        action  new;
        source  patchToFace;
        sourceInfo {
            //patch   body;
            name   "body";
            }
    }
    {
        name    airIntake;
        type    faceSet;
        action  subset;
        source  boxToFace;
        sourceInfo {
            box     (2.6 0.75 0)(2.64 0.8 0.1);
            }
    }
);