songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

porosityModel: The "<porosityModel>Coeffs" sub-dictionary is now optional

Browse Source 
For example the porosity coefficients may now be specified thus:

porosity1
{
    type            DarcyForchheimer;

    cellZone        porosity;

    d   (5e7 -1000 -1000);
    f   (0 0 0);

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

rather than

porosity1
{
    type            DarcyForchheimer;
    active          yes;
    cellZone        porosity;

    DarcyForchheimerCoeffs
    {
        d   (5e7 -1000 -1000);
        f   (0 0 0);

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

support for which is maintained for backward compatibility.
This commit is contained in:
Henry Weller 2017-04-13 14:00:00 +01:00
parent 1805b3c98f
commit af2810149e
4 changed files with 52 additions and 47 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
src/finiteVolume/cfdTools/general/porosityModel/porosityModel/porosityModel.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2012-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -97,7 +97,12 @@ Foam::porosityModel::porosityModel
name_(name),
mesh_(mesh),
dict_(dict),
coeffs_(dict.subDict(modelType + "Coeffs")),
coeffs_
(
dict.found(modelType + "Coeffs")
? dict.subDict(modelType + "Coeffs")
: dict
),
active_(true),
zoneName_(cellZoneName),
cellZoneIDs_(),
@ -105,7 +110,7 @@ Foam::porosityModel::porosityModel
{
if (zoneName_ == word::null)
{
dict.lookup("active") >> active_;
dict.readIfPresent("active", active_);
dict_.lookup("cellZone") >> zoneName_;
}
@ -227,8 +232,17 @@ bool Foam::porosityModel::writeData(Ostream& os) const
bool Foam::porosityModel::read(const dictionary& dict)
{
active_ = readBool(dict.lookup("active"));
dict.readIfPresent("active", active_);
if (dict.found(type() + "Coeffs"))
{
coeffs_ = dict.subDict(type() + "Coeffs");
}
else
{
coeffs_ = dict;
}
dict.lookup("cellZone") >> zoneName_;
cellZoneIDs_ = mesh_.cellZones().findIndices(zoneName_);

5
tutorials/compressible/rhoPorousSimpleFoam/angledDuctImplicit/constant/porosityProperties
View File

@ -18,11 +18,9 @@ FoamFile
porosity1
{
type DarcyForchheimer;
active yes;
cellZone porosity;
DarcyForchheimerCoeffs
{
d (5e7 -1000 -1000);
f (0 0 0);
@ -37,7 +35,6 @@ porosity1
e2 (0 0 1);
}
}
}
}

5
tutorials/incompressible/porousSimpleFoam/angledDuctImplicit/constant/porosityProperties
View File

@ -18,11 +18,9 @@ FoamFile
porosity1
{
type DarcyForchheimer;
active yes;
cellZone porosity;
DarcyForchheimerCoeffs
{
d (5e7 -1000 -1000);
f (0 0 0);
@ -37,7 +35,6 @@ porosity1
e2 (0 0 1);
}
}
}
}

5
tutorials/incompressible/porousSimpleFoam/straightDuctImplicit/constant/porosityProperties
View File

@ -18,11 +18,9 @@ FoamFile
porosity1
{
type DarcyForchheimer;
active yes;
cellZone porosity;
DarcyForchheimerCoeffs
{
d (5e7 -1000 -1000);
f (0 0 0);
@ -37,7 +35,6 @@ porosity1
e2 (0 0 1);
}
}
}
}