porousmediumflow/sequential/properties.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*****************************************************************************
 *   See the file COPYING for full copying permissions.                      *
 *                                                                           *
 *   This program is free software: you can redistribute it and/or modify    *
 *   it under the terms of the GNU General Public License as published by    *
 *   the Free Software Foundation, either version 3 of the License, or       *
 *   (at your option) any later version.                                     *
 *                                                                           *
 *   This program is distributed in the hope that it will be useful,         *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of          *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the            *
 *   GNU General Public License for more details.                            *
 *                                                                           *
 *   You should have received a copy of the GNU General Public License       *
 *   along with this program.  If not, see <http://www.gnu.org/licenses/>.   *
 *****************************************************************************/
#ifndef DUMUX_SEQUENTIAL_PROPERTIES_HH
#define DUMUX_SEQUENTIAL_PROPERTIES_HH
 
#include <dumux/common/properties.hh>
#include <dumux/common/properties/model.hh>
#include <dumux/common/properties/grid.hh>
#include <dumux/common/defaultmappertraits.hh>
#include <dumux/discretization/method.hh>
#include <dumux/porousmediumflow/sequential/gridadaptproperties.hh>
#include <dumux/porousmediumflow/sequential/gridadaptinitializationindicatordefault.hh>
 

namespace Dumux
{
namespace Properties
{

// Type tags

NEW_TYPE_TAG(SequentialModel, INHERITS_FROM(ModelProperties, GridAdapt, GridProperties));

// Property tags

NEW_PROP_TAG( SolutionTypes);
NEW_PROP_TAG( Indices);
 
// Some properties that have been removed from numeric model
NEW_PROP_TAG( Model ); 
NEW_PROP_TAG( DiscretizationMethod ); 
 
NEW_PROP_TAG( PressureModel ); 
NEW_PROP_TAG( TransportModel ); 
NEW_PROP_TAG( Velocity ); 
NEW_PROP_TAG( NumEq ); 
NEW_PROP_TAG( NumPhases); 
NEW_PROP_TAG( NumComponents); 
NEW_PROP_TAG( Variables); 
NEW_PROP_TAG( CellData );
NEW_PROP_TAG( MaxIntersections ); 
NEW_PROP_TAG( PressureCoefficientMatrix ); 
}
}
 
#include <dune/grid/common/mcmgmapper.hh>
#include <dune/istl/bvector.hh>
 
#include <dumux/common/timemanager.hh>
#include <dumux/common/boundarytypes.hh>
#include<dumux/common/boundaryconditions.hh>
#include<dumux/discretization/method.hh>
 
namespace Dumux
{
 
template<class TypeTag>
class VariableClass;
 
namespace Properties
{

// Properties

SET_TYPE_PROP(SequentialModel, JacobianMatrix, typename GET_PROP_TYPE(TypeTag, PressureCoefficientMatrix));

template<class TypeTag>
struct ModelTraits<TypeTag, TTag::SequentialModel>
{
private:
    struct DummyTraits
    {
        using Indices = typename GET_PROP_TYPE(TypeTag, Indices);
        static constexpr int numEq() { return GET_PROP_VALUE(TypeTag, NumEq); }
    };
public:
    using type = DummyTraits;
};

template<class TypeTag>
struct GridView<TypeTag, TTag::SequentialModel>
{
private:
    using Grid = typename GET_PROP_TYPE(TypeTag, Grid);
 
public:
    using type = typename Grid::LeafGridView;
};

template<class TypeTag>
struct MaxIntersections<TypeTag, TTag::SequentialModel>
{
private:
    using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
    enum
    {
        dim = GridView::dimension
    };
public:
    static const int value = 2*dim;
};
 
// Dumux 3.1 changes the property `FVGridGeometry` to `GridGeometry`.
// For ensuring backward compatibility on the user side, it is necessary to
// stick to the old name for the specializations, see the discussion in MR 1647.
// Use diagnostic pragmas to prevent the emission of a warning message.
// TODO after 3.1: Rename to GridGeometry, remove the pragmas and this comment.
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wdeprecated-declarations"
template<class TypeTag>
struct FVGridGeometry<TypeTag, TTag::SequentialModel>
{
    struct MockFVGridGeometry
    : public DefaultMapperTraits<typename GET_PROP_TYPE(TypeTag, GridView)>
    {
        static constexpr Dumux::DiscretizationMethod discMethod = Dumux::DiscretizationMethod::cctpfa;
        using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
    };
 
public:
    using type = MockFVGridGeometry;
};
#pragma GCC diagnostic pop

template<class TypeTag>
struct SolutionVector<TypeTag, TTag::SequentialModel>
{
public:
    using type = typename GET_PROP(TypeTag, SolutionTypes)::ScalarSolution;
};

template<class TypeTag>
struct SolutionTypes<TypeTag, TTag::SequentialModel>
{
    using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar);
    using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
    using Grid = typename GridView::Grid;
    using Variables = typename GET_PROP_TYPE(TypeTag, Variables);
 
    enum
    {
        dim = GridView::dimension,
        numEq = GET_PROP_VALUE(TypeTag, NumEq),
        numPhases = GET_PROP_VALUE(TypeTag, NumPhases),
        numComponents = GET_PROP_VALUE(TypeTag, NumComponents),
        maxIntersections = GET_PROP_VALUE(TypeTag, MaxIntersections)
    };
 
public:
    using VertexMapper = Dune::MultipleCodimMultipleGeomTypeMapper<GridView>;

    using ElementMapper = Dune::MultipleCodimMultipleGeomTypeMapper<GridView>;

    using PrimaryVariables = Dune::FieldVector<Scalar, numEq>;
    using ScalarSolution = Dune::BlockVector<Dune::FieldVector<Scalar, 1> >;
    using ComponentProperty = Dune::FieldVector<Dune::BlockVector<Dune::FieldVector<Scalar,1> >, numComponents>;
    using PhaseProperty = Dune::FieldVector<Dune::BlockVector<Dune::FieldVector<Scalar,1> >, numPhases>;
    using FluidProperty = Dune::FieldVector<Dune::BlockVector<Dune::FieldVector<Scalar,1> >, numPhases>;
    using PhasePropertyElemFace = Dune::BlockVector<Dune::FieldVector<Dune::FieldVector<Scalar, numPhases>, maxIntersections> >;
    using DimVecElemFace = Dune::BlockVector<Dune::FieldVector<Dune::FieldVector<Scalar, dim>, maxIntersections> >;
};
 
SET_TYPE_PROP(SequentialModel,  Variables, VariableClass<TypeTag>);
 
SET_TYPE_PROP(SequentialModel,  PrimaryVariables, typename GET_PROP(TypeTag, SolutionTypes)::PrimaryVariables);

SET_TYPE_PROP(SequentialModel, TimeManager, Dumux::TimeManager<TypeTag>);

template<class TypeTag>
struct BoundaryTypes<TypeTag, TTag::SequentialModel>
{ private:
    enum { numEq = GET_PROP_VALUE(TypeTag, NumEq) };
public:
    using type = Dumux::BoundaryTypes<numEq>;
};

template<class TypeTag>
struct ModelDefaultParameters<TypeTag, TTag::SequentialModel>
{
    static void defaultParams(Dune::ParameterTree& params, const std::string& group = "")
    {
        params["GridAdapt.CoarsenTolerance"] = "0.001"; 
        params["GridAdapt.EnableInitializationIndicator"] = "false"; 
        params["GridAdapt.EnableMultiPointFluxApproximation"] = "true"; 
        params["GridAdapt.MaxLevel"] = "1"; 
        params["GridAdapt.MaxInteractionVolumes"] = "4"; 
        params["GridAdapt.MinLevel"] = "0"; 
        params["GridAdapt.RefineAtDirichletBC"] = "false"; 
        params["GridAdapt.RefineAtFluxBC"] = "false"; 
        params["GridAdapt.RefineAtSource"] = "false"; 
        params["GridAdapt.RefineTolerance"] = "0.05"; 
 
        params["Impet.CFLFactor"] = "1.0"; 
        params["Impet.EnableVolumeIntegral"] = "true"; 
        params["Impet.ErrorTermFactor"] = "0.5"; 
        params["Impet.ErrorTermLowerBound"] = "0.1"; 
        params["Impet.ErrorTermUpperBound"] = "0.9"; 
        params["Impet.PorosityThreshold"] = "1e-6"; 
        params["Impet.SubCFLFactor"] = "1.0"; 
        params["Impet.SwitchNormals"] = "false"; 
 
        params["MPFA.CalcVelocityInTransport"] = "false"; 
 
        params["TimeManager.SubTimestepVerbosity"] = "0"; 
 
        params["Vtk.OutputLevel"] = "0"; 
    }
};
 
}
}
 
#endif