Detailed Description

template<typename ShapeFunction, typename IntegrationMethod>
class ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >

Definition at line 29 of file HeatTransportBHELocalAssemblerSoil.h.

#include <HeatTransportBHELocalAssemblerSoil.h>

Inheritance diagram for ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >:

[legend]

Collaboration diagram for ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >:

[legend]

Public Types
using	ShapeMatricesType = ShapeMatrixPolicyType< ShapeFunction, 3 >

using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

using	NodalVectorType = typename ShapeMatricesType::NodalVectorType

using	NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType

using	ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

using	GlobalDimNodalMatrixType = typename ShapeMatricesType::GlobalDimNodalMatrixType

Public Member Functions
	HeatTransportBHELocalAssemblerSoil (HeatTransportBHELocalAssemblerSoil const &)=delete

	HeatTransportBHELocalAssemblerSoil (HeatTransportBHELocalAssemblerSoil &&)=delete

	HeatTransportBHELocalAssemblerSoil (MeshLib::Element const &e, bool is_axially_symmetric, unsigned const integration_order, HeatTransportBHEProcessData &process_data)

void	assemble (double const, double const, std::vector< double > const &, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &) override

Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const override
	Provides the shape matrix at the given integration point. More...

Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual	~LocalAssemblerInterface ()=default

void	setInitialConditions (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, bool const use_monolithic_scheme, int const process_id)

virtual void	initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)

virtual void	preAssemble (double const, double const, std::vector< double > const &)

virtual void	assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_xdot, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)

virtual void	assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_xdot, const double dxdot_dx, const double dx_dx, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)

virtual void	assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_xdot, const double dxdot_dx, const double dx_dx, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)

virtual void	computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id)

virtual void	preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)

virtual void	postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double const dt)

void	postNonLinearSolver (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, GlobalVector const &xdot, double const t, double const dt, bool const use_monolithic_scheme, int const process_id)

virtual std::vector< double >	interpolateNodalValuesToIntegrationPoints (std::vector< double > const &)

virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const

virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double >> const &) const
	Fits to staggered scheme. More...

Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual	~ExtrapolatableElement ()=default

Private Attributes
HeatTransportBHEProcessData &	_process_data

std::vector< IntegrationPointDataSoil< NodalRowVectorType, GlobalDimNodalMatrixType >, Eigen::aligned_allocator< IntegrationPointDataSoil< NodalRowVectorType, GlobalDimNodalMatrixType > > >	_ip_data

IntegrationMethod const	_integration_method

std::vector< ShapeMatrices, Eigen::aligned_allocator< ShapeMatrices > >	_shape_matrices

std::size_t const	_element_id

SecondaryData< typename ShapeMatrices::ShapeType >	_secondary_data

Member Typedef Documentation

◆ GlobalDimNodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::GlobalDimNodalMatrixType = typename ShapeMatricesType::GlobalDimNodalMatrixType

Definition at line 40 of file HeatTransportBHELocalAssemblerSoil.h.

◆ NodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

Definition at line 35 of file HeatTransportBHELocalAssemblerSoil.h.

◆ NodalRowVectorType

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType

Definition at line 37 of file HeatTransportBHELocalAssemblerSoil.h.

◆ NodalVectorType

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

Definition at line 36 of file HeatTransportBHELocalAssemblerSoil.h.

◆ ShapeMatrices

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

Definition at line 39 of file HeatTransportBHELocalAssemblerSoil.h.

◆ ShapeMatricesType

template<typename ShapeFunction , typename IntegrationMethod >

using ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, 3 >

Definition at line 33 of file HeatTransportBHELocalAssemblerSoil.h.

Constructor & Destructor Documentation

◆ HeatTransportBHELocalAssemblerSoil() [1/3]

template<typename ShapeFunction , typename IntegrationMethod >

ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::HeatTransportBHELocalAssemblerSoil ( HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod > const & )

delete

◆ HeatTransportBHELocalAssemblerSoil() [2/3]

template<typename ShapeFunction , typename IntegrationMethod >

ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::HeatTransportBHELocalAssemblerSoil ( HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod > && )

delete

◆ HeatTransportBHELocalAssemblerSoil() [3/3]

template<typename ShapeFunction , typename IntegrationMethod >

ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::HeatTransportBHELocalAssemblerSoil	(	MeshLib::Element const &	e,
		bool	is_axially_symmetric,
		unsigned const	integration_order,
		HeatTransportBHEProcessData &	process_data
	)

Definition at line 31 of file HeatTransportBHELocalAssemblerSoil-impl.h.

     : _process_data(process_data),
       _integration_method(integration_order),
       _element_id(e.getID())
 {
     unsigned const n_integration_points =
         _integration_method.getNumberOfPoints();
  
     _ip_data.reserve(n_integration_points);
     _secondary_data.N.resize(n_integration_points);
  
     _shape_matrices =
         NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                   3 /* GlobalDim */>(e, is_axially_symmetric,
                                                      _integration_method);
  
     ParameterLib::SpatialPosition x_position;
     x_position.setElementID(_element_id);
  
     // ip data initialization
     for (unsigned ip = 0; ip < n_integration_points; ip++)
     {
         x_position.setIntegrationPoint(ip);
  
         // create the class IntegrationPointDataBHE in place
         auto const& sm = _shape_matrices[ip];
         double const w = _integration_method.getWeightedPoint(ip).getWeight() *
                          sm.integralMeasure * sm.detJ;
         _ip_data.push_back({sm.N, sm.dNdx, w});
  
         _secondary_data.N[ip] = sm.N;
     }
 }

References ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::_element_id, ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::_integration_method, ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::_ip_data, ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::_secondary_data, ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::_shape_matrices, NumLib::initShapeMatrices(), ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N, and ParameterLib::SpatialPosition::setElementID().

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , typename IntegrationMethod >

void ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::assemble	(	double const	t,
		double const	dt,
		std::vector< double > const &	local_x,
		std::vector< double > const &	,
		std::vector< double > &	local_M_data,
		std::vector< double > &	local_K_data,
		std::vector< double > &
	)

overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 71 of file HeatTransportBHELocalAssemblerSoil-impl.h.

 {
     assert(local_x.size() == ShapeFunction::NPOINTS);
     (void)local_x;  // Avoid unused arg warning.
  
     auto local_M = MathLib::createZeroedMatrix<NodalMatrixType>(
         local_M_data, ShapeFunction::NPOINTS, ShapeFunction::NPOINTS);
     auto local_K = MathLib::createZeroedMatrix<NodalMatrixType>(
         local_K_data, ShapeFunction::NPOINTS, ShapeFunction::NPOINTS);
  
     ParameterLib::SpatialPosition pos;
     pos.setElementID(_element_id);
  
     auto const& medium = *_process_data.media_map->getMedium(_element_id);
     auto const& solid_phase = medium.phase("Solid");
     auto const& liquid_phase = medium.phase("AqueousLiquid");
  
     MaterialPropertyLib::VariableArray vars;
  
     unsigned const n_integration_points =
         _integration_method.getNumberOfPoints();
  
     for (unsigned ip = 0; ip < n_integration_points; ip++)
     {
         pos.setIntegrationPoint(ip);
         auto& ip_data = _ip_data[ip];
         auto const& N = ip_data.N;
         auto const& dNdx = ip_data.dNdx;
         auto const& w = ip_data.integration_weight;
  
         double T_int_pt = 0.0;
         NumLib::shapeFunctionInterpolate(local_x, N, T_int_pt);
  
         vars[static_cast<int>(MaterialPropertyLib::Variable::temperature)] =
             T_int_pt;
  
         // for now only using the solid and liquid phase parameters
         auto const density_s =
             solid_phase.property(MaterialPropertyLib::PropertyType::density)
                 .template value<double>(vars, pos, t, dt);
  
         auto const heat_capacity_s =
             solid_phase
                 .property(
                     MaterialPropertyLib::PropertyType::specific_heat_capacity)
                 .template value<double>(vars, pos, t, dt);
  
         auto const density_f =
             liquid_phase.property(MaterialPropertyLib::PropertyType::density)
                 .template value<double>(vars, pos, t, dt);
  
         auto const heat_capacity_f =
             liquid_phase
                 .property(
                     MaterialPropertyLib::PropertyType::specific_heat_capacity)
                 .template value<double>(vars, pos, t, dt);
  
         auto const porosity =
             medium.property(MaterialPropertyLib::PropertyType::porosity)
                 .template value<double>(vars, pos, t, dt);
  
         auto const velocity =
             liquid_phase
                 .property(MaterialPropertyLib::PropertyType::phase_velocity)
                 .template value<Eigen::Vector3d>(vars, pos, t, dt);
  
         // calculate the hydrodynamic thermodispersion tensor
         auto const thermal_conductivity =
             MaterialPropertyLib::formEigenTensor<3>(
                 medium
                     .property(
                         MaterialPropertyLib::PropertyType::thermal_conductivity)
                     .value(vars, pos, t, dt));
  
         auto thermal_conductivity_dispersivity = thermal_conductivity;
  
         double const velocity_magnitude = velocity.norm();
  
         if (velocity_magnitude >= std::numeric_limits<double>::epsilon())
         {
             auto const thermal_dispersivity_longitudinal =
                 medium
                     .property(MaterialPropertyLib::PropertyType::
                                   thermal_longitudinal_dispersivity)
                     .template value<double>();
             auto const thermal_dispersivity_transversal =
                 medium
                     .property(MaterialPropertyLib::PropertyType::
                                   thermal_transversal_dispersivity)
                     .template value<double>();
  
             auto const thermal_dispersivity =
                 density_f * heat_capacity_f *
                 (thermal_dispersivity_transversal * velocity_magnitude *
                      Eigen::Matrix3d::Identity() +
                  (thermal_dispersivity_longitudinal -
                   thermal_dispersivity_transversal) /
                      velocity_magnitude * velocity * velocity.transpose());
             thermal_conductivity_dispersivity += thermal_dispersivity;
         }
  
         // assemble Conductance matrix
         local_K.noalias() +=
             (dNdx.transpose() * thermal_conductivity_dispersivity * dNdx +
              N.transpose() * velocity.transpose() * dNdx * density_f *
                  heat_capacity_f) *
             w;
  
         // assemble Mass matrix
         local_M.noalias() += N.transpose() * N * w *
                              (density_s * heat_capacity_s * (1 - porosity) +
                               density_f * heat_capacity_f * porosity);
     }
  
     // debugging
     // std::string sep = "\n----------------------------------------\n";
     // Eigen::IOFormat CleanFmt(4, 0, ", ", "\n", "[", "]");
     // std::cout << local_K.format(CleanFmt) << sep;
     // std::cout << local_M.format(CleanFmt) << sep;
 }

◆ getShapeMatrix()

template<typename ShapeFunction , typename IntegrationMethod >

Eigen::Map<const Eigen::RowVectorXd> ProcessLib::HeatTransportBHE::HeatTransportBHELocalAssemblerSoil< ShapeFunction, IntegrationMethod >::getShapeMatrix ( const unsigned integration_point ) const

inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 61 of file HeatTransportBHELocalAssemblerSoil.h.

     {
         auto const& N = _secondary_data.N[integration_point];
  
         // assumes N is stored contiguously in memory
         return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
     }