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
|