CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Compressible, buoyant BC for porous flow (https://www.cfd-online.com/Forums/openfoam-pre-processing/122914-compressible-buoyant-bc-porous-flow.html)

vainilreb August 31, 2013 10:25

Compressible, buoyant BC for porous flow
 
Hi,

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

http://knallrot.org/skizze.png
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.

Code:

/*--------------------------------*- 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;
    }
}


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

Code:

/*--------------------------------*- 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;
    }
}


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

Code:

/*--------------------------------*- 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);
    }
}


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

Code:

/*--------------------------------*- 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


All times are GMT -4. The time now is 17:44.