TH2MProcessData.h
/**
* \file
* \copyright
* Copyright (c) 2012-2022, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*/
#pragma once
#include <Eigen/Core>
#include <memory>
#include <utility>
#include "MaterialLib/MPL/MaterialSpatialDistributionMap.h"
#include "ParameterLib/Parameter.h"
#include "PhaseTransitionModels/NoPhaseTransition.h"
#include "PhaseTransitionModels/PhaseTransition.h"
#include "PhaseTransitionModels/PhaseTransitionModel.h"
namespace MaterialLib
{
namespace Solids
{
template <int DisplacementDim>
struct MechanicsBase;
}
} // namespace MaterialLib
namespace ProcessLib
{
namespace TH2M
{
template <int DisplacementDim>
struct TH2MProcessData
{
MeshLib::PropertyVector<int> const* const material_ids = nullptr;
std::unique_ptr<MaterialPropertyLib::MaterialSpatialDistributionMap>
media_map = nullptr;
/// The constitutive relation for the mechanical part.
/// \note Linear elasticity is the only supported one in the moment.
std::map<int, std::unique_ptr<
MaterialLib::Solids::MechanicsBase<DisplacementDim>>>
solid_materials;
std::unique_ptr<PhaseTransitionModel> phase_transition_model_ = nullptr;
ParameterLib::Parameter<double> const& reference_temperature;
/// Optional, initial stress field. A symmetric tensor, short vector
/// representation of length 4 or 6, ParameterLib::Parameter<double>.
ParameterLib::Parameter<double> const* const initial_stress;
/// Specific body forces applied to solid and both fluid phases.
/// It is usually used to apply gravitational forces.
/// A vector of displacement dimension's length.
Eigen::Matrix<double, DisplacementDim, 1> const specific_body_force;
bool const apply_mass_lumping;
const bool use_TaylorHood_elements;
MeshLib::PropertyVector<double>* element_saturation = nullptr;
MeshLib::PropertyVector<double>* gas_pressure_interpolated = nullptr;
MeshLib::PropertyVector<double>* capillary_pressure_interpolated = nullptr;
MeshLib::PropertyVector<double>* temperature_interpolated = nullptr;
MeshLib::PropertyVector<double>* liquid_pressure_interpolated = nullptr;
EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
};
} // namespace TH2M
} // namespace ProcessLib