[vainilreb]

Hi,

I'm trying to set up a case like seen in the following sketch:


The solver is chtMultiRegionFoam because I've also implemented a thermalMass zone around the tube. Buoyancy is turned off by setting g to zero.
I've tried a lot of things but I always end up with large vortices in the porous zone. Everything worked fine without multi-region in rhoPimpleFoam,
so I guess my BC settings are to blame.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 101300;

boundaryField
{
    Front_Plate
    {
        type            zeroGradient;
    }
    End_Plate
    {
        type            zeroGradient;
    }
    Propene_Inlet
    {
        type            zeroGradient;
    }
    Mid_Cylinder
    {
        type            zeroGradient;
    }
    Inner_Quartz_Tube
    {
        type            zeroGradient;
    }
    Probe_Porous_Outlet
    {
        type            fixedMean;
    meanValue    101300;
    value        uniform 101300;
    }
    Probe_Walls
    {
        type            zeroGradient;
    }
    Gas_Mixture_Inlet
    {
        type            zeroGradient;
    }
    flowRegion_to_thermalMass
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 101300;

boundaryField
{
    Front_Plate
    {
        type  calculated;
    value uniform 101300;
    }
    End_Plate
    {
        type  calculated;
    value uniform 101300;
    }
    Propene_Inlet
    {
        type  calculated;
    value uniform 101300;
    }
    Mid_Cylinder
    {
        type  calculated;
    value uniform 101300;
    }
    Inner_Quartz_Tube
    {
        type  calculated;
    value uniform 101300;
    }
    Probe_Porous_Outlet
    {
        type  calculated;
    value uniform 101300;
    }
    Probe_Walls
    {
        type  calculated;
    value uniform 101300;
    }
    Gas_Mixture_Inlet
    {
        type  calculated;
    value uniform 101300;
    }
    flowRegion_to_thermalMass
    {
        type  calculated;
    value uniform 101300;
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
    Front_Plate
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    End_Plate
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Propene_Inlet
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Mid_Cylinder
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Inner_Quartz_Tube
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Probe_Porous_Outlet
    {
        type            inletOutlet;
        value           uniform (0 0 0);
        inletValue      uniform ( 0 0 0 );
    }
    Probe_Walls
    {
        type            fixedValue;
        value           uniform (0 0 0);

    }
    Gas_Mixture_Inlet
    {
        type            uniformFixedValue;
        uniformValue    table
    (
    (0 (0 0 0))
    (1 (2.2366 0 0))
    );
    }
    flowRegion_to_thermalMass
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 293.15;

boundaryField
{
    Front_Plate
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    End_Plate
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    Propene_Inlet
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    Mid_Cylinder
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    Inner_Quartz_Tube
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    Probe_Porous_Outlet
    {
        type            inletOutlet;
        value           uniform 293.15;
        inletValue      $internalField;
    }
    Probe_Walls
    {
        type            zeroGradient;
    value        uniform 293.15;
    }
    Gas_Mixture_Inlet
    {
        type            uniformFixedValue;
        uniformValue    table
    (    
    (0 293.15)
    (1 443.15)
    );
    }
    flowRegion_to_thermalMass
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        value           uniform 293.15;
        neighbourFieldName T;
        kappa           fluidThermo;
        kappaName       none;
    }
}


// ************************************************************************* //

I'd really appreciate your help. Thanks in advance!

Robert