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/boundarytypes.hh>
#include <dumux/porousmediumflow/sequential/gridadaptproperties.hh>
#include <dumux/porousmediumflow/sequential/gridadaptinitializationindicatordefault.hh>
 

namespace Dumux
{
namespace Properties
{

// Type tags

// Create new type tags
namespace TTag {
struct SequentialModel { using InheritsFrom = std::tuple<GridAdapt, GridProperties, ModelProperties>; };
} // end namespace TTag

// Property tags

template<class TypeTag, class MyTypeTag>
struct  SolutionTypes { using type = UndefinedProperty; };
template<class TypeTag, class MyTypeTag>
struct  Indices { using type = UndefinedProperty; };
template<class TypeTag, class MyTypeTag>
struct SequentialBoundaryTypes { using type = UndefinedProperty; };
 
// Some properties that have been removed from numeric model
template<class TypeTag, class MyTypeTag>
struct  Model  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  DiscretizationMethod  { using type = UndefinedProperty; }; 
 
template<class TypeTag, class MyTypeTag>
struct  PressureModel  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  TransportModel  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  Velocity  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  NumEq  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  NumPhases { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  NumComponents { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  Variables { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  CellData  { using type = UndefinedProperty; };
template<class TypeTag, class MyTypeTag>
struct  MaxIntersections  { using type = UndefinedProperty; }; 
template<class TypeTag, class MyTypeTag>
struct  PressureCoefficientMatrix  { using type = UndefinedProperty; }; 
}
}
 
#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

template<class TypeTag>
struct JacobianMatrix<TypeTag, TTag::SequentialModel> { using type = GetPropType<TypeTag, Properties::PressureCoefficientMatrix>; };

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

template<class TypeTag>
struct MaxIntersections<TypeTag, TTag::SequentialModel>
{
private:
    using GridView = typename GetPropType<TypeTag, Properties::GridGeometry>::GridView;
public:
    static constexpr int value = 2*GridView::dimension;
};

template<class TypeTag>
struct GridGeometry<TypeTag, TTag::SequentialModel>
{
    using GV = typename GetPropType<TypeTag, Properties::Grid>::LeafGridView;
    struct MockFVGridGeometry
    : public DefaultMapperTraits<GV>
    {
        static constexpr Dumux::DiscretizationMethod discMethod = Dumux::DiscretizationMethod::cctpfa;
        using GridView = GV;
    };
public:
    using type = MockFVGridGeometry;
};

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

template<class TypeTag>
struct SolutionTypes<TypeTag, TTag::SequentialModel>
{
    using Scalar = GetPropType<TypeTag, Properties::Scalar>;
    using GridView = typename GetPropType<TypeTag, Properties::GridGeometry>::GridView;
    using Variables = GetPropType<TypeTag, Properties::Variables>;
 
    enum
    {
        dim = GridView::dimension,
        numEq = getPropValue<TypeTag, Properties::NumEq>(),
        numPhases = getPropValue<TypeTag, Properties::NumPhases>(),
        numComponents = getPropValue<TypeTag, Properties::NumComponents>(),
        maxIntersections = getPropValue<TypeTag, Properties::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> >;
};
 
template<class TypeTag>
struct Variables<TypeTag, TTag::SequentialModel> { using type = VariableClass<TypeTag>; };
 
template<class TypeTag>
struct PrimaryVariables<TypeTag, TTag::SequentialModel> { using type = typename GetProp<TypeTag, SolutionTypes>::PrimaryVariables; };

template<class TypeTag>
struct TimeManager<TypeTag, TTag::SequentialModel> { using type = Dumux::TimeManager<TypeTag>; };

template<class TypeTag>
struct SequentialBoundaryTypes<TypeTag, TTag::SequentialModel>
{ private:
    enum { numEq = getPropValue<TypeTag, Properties::NumEq>() };
public:
    using type = Dumux::SequentialBoundaryTypes<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