28 namespace ThermalTwoPhaseFlowWithPP
30 template <
typename NodalRowVectorType,
typename GlobalDimNodalMatrixType,
31 typename NodalMatrixType>
35 NodalRowVectorType
const& N_, GlobalDimNodalMatrixType
const& dNdx_,
37 double const& integration_weight_, NodalMatrixType
const mass_operator_,
38 NodalMatrixType
const diffusion_operator_)
47 NodalRowVectorType
const N;
48 GlobalDimNodalMatrixType
const dNdx;
65 std::vector<GlobalVector*>
const& x,
66 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& dof_table,
67 std::vector<double>& cache)
const = 0;
71 std::vector<GlobalVector*>
const& x,
72 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& dof_table,
73 std::vector<double>& cache)
const = 0;
76 template <
typename ShapeFunction,
typename IntegrationMethod,
int GlobalDim>
100 bool const is_axially_symmetric,
101 unsigned const integration_order,
111 unsigned const n_integration_points =
113 _ip_data.reserve(n_integration_points);
114 auto const shape_matrices =
116 GlobalDim>(element, is_axially_symmetric,
118 for (
unsigned ip = 0; ip < n_integration_points; ip++)
120 auto const& sm = shape_matrices[ip];
121 const double integration_factor = sm.integralMeasure * sm.detJ;
124 sm.integralMeasure * sm.detJ *
126 sm.N.transpose() * sm.N * integration_factor *
128 sm.dNdx.transpose() * sm.dNdx * integration_factor *
133 void assemble(
double const t,
double const dt,
134 std::vector<double>
const& local_x,
135 std::vector<double>
const& local_xdot,
136 std::vector<double>& local_M_data,
137 std::vector<double>& local_K_data,
138 std::vector<double>& local_b_data)
override;
141 const unsigned integration_point)
const override
143 auto const& N =
_ip_data[integration_point].N;
146 return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
151 std::vector<GlobalVector*>
const& ,
152 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
153 std::vector<double>& )
const override
161 std::vector<GlobalVector*>
const& ,
162 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
163 std::vector<double>& )
const override
virtual std::vector< double > const & getIntPtSaturation(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0
virtual std::vector< double > const & getIntPtWettingPressure(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0
ShapeMatrixPolicyType< ShapeFunction, GlobalDim > ShapeMatricesType
static const int nonwet_pressure_matrix_index
std::vector< double > _pressure_wetting
static const int cap_pressure_matrix_index
std::vector< double > const & getIntPtSaturation(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &) const override
std::vector< IntegrationPointData< NodalRowVectorType, GlobalDimNodalMatrixType, NodalMatrixType >, Eigen::aligned_allocator< IntegrationPointData< NodalRowVectorType, GlobalDimNodalMatrixType, NodalMatrixType > > > _ip_data
typename ShapeMatricesType::NodalVectorType NodalVectorType
typename ShapeMatricesType::GlobalDimVectorType GlobalDimVectorType
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix(const unsigned integration_point) const override
Provides the shape matrix at the given integration point.
void assemble(double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_xdot, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data) override
std::vector< double > const & getIntPtWettingPressure(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &) const override
MeshLib::Element const & _element
typename ShapeMatricesType::NodalRowVectorType NodalRowVectorType
IntegrationMethod const _integration_method
ThermalTwoPhaseFlowWithPPLocalAssembler(MeshLib::Element const &element, std::size_t const, bool const is_axially_symmetric, unsigned const integration_order, ThermalTwoPhaseFlowWithPPProcessData const &process_data)
typename ShapeMatricesType::GlobalDimMatrixType GlobalDimMatrixType
static const int temperature_size
typename ShapeMatricesType::ShapeMatrices ShapeMatrices
typename ShapeMatricesType::NodalMatrixType NodalMatrixType
typename ShapeMatricesType::GlobalDimNodalMatrixType GlobalDimNodalMatrixType
std::vector< double > _saturation
static const int temperature_matrix_index
typename LocalAssemblerTraits::LocalMatrix LocalMatrixType
ThermalTwoPhaseFlowWithPPProcessData const & _process_data
typename LocalAssemblerTraits::LocalVector LocalVectorType
static const int nonwet_pressure_size
static const int cap_pressure_size
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
const unsigned NUM_NODAL_DOF
MatrixType< ShapeFunction::NPOINTS, ShapeFunction::NPOINTS > NodalMatrixType
NumLib::ShapeMatrices< NodalRowVectorType, DimNodalMatrixType, DimMatrixType, GlobalDimNodalMatrixType > ShapeMatrices
MatrixType< GlobalDim, ShapeFunction::NPOINTS > GlobalDimNodalMatrixType
MatrixType< GlobalDim, GlobalDim > GlobalDimMatrixType
VectorType< GlobalDim > GlobalDimVectorType
VectorType< ShapeFunction::NPOINTS > NodalVectorType
RowVectorType< ShapeFunction::NPOINTS > NodalRowVectorType
ThermalTwoPhaseFlowWithPPMaterialProperties const & mat_property
NodalMatrixType const diffusion_operator
double const integration_weight
IntegrationPointData(NodalRowVectorType const &N_, GlobalDimNodalMatrixType const &dNdx_, ThermalTwoPhaseFlowWithPPMaterialProperties const &material_property_, double const &integration_weight_, NodalMatrixType const mass_operator_, NodalMatrixType const diffusion_operator_)
GlobalDimNodalMatrixType const dNdx
NodalMatrixType const mass_operator
NodalRowVectorType const N
EIGEN_MAKE_ALIGNED_OPERATOR_NEW
std::unique_ptr< ThermalTwoPhaseFlowWithPPMaterialProperties > material
Vector< NNodes *NodalDOF > LocalVector
Matrix< NNodes *NodalDOF, NNodes *NodalDOF > LocalMatrix