Detailed Description

template<typename ShapeFunction, int GlobalDim>
class ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >

Definition at line 57 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

#include <ConstraintDirichletBoundaryConditionLocalAssembler.h>

Inheritance diagram for ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >:

[legend]

Collaboration diagram for ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >:

[legend]

Public Member Functions
	ConstraintDirichletBoundaryConditionLocalAssembler (MeshLib::Element const &surface_element, std::size_t, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, MeshLib::Mesh const &bulk_mesh, std::vector< std::pair< std::size_t, unsigned > > bulk_ids)

double	integrate (std::vector< GlobalVector * > const &x, double const t, std::function< Eigen::Vector3d(std::size_t const, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &)> const &getFlux) override

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

Protected Types
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

using	NodalVectorType = typename ShapeMatricesType::NodalVectorType

Private Attributes
MeshLib::Element const &	_surface_element

std::vector< IntegrationPointData >	_ip_data

NumLib::GenericIntegrationMethod const &	_integration_method

std::size_t const	_bulk_element_id

Eigen::Vector3d const	_surface_element_normal

Member Typedef Documentation

◆ NodalMatrixType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

protected

Definition at line 62 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

protected

Definition at line 63 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

◆ ShapeMatricesType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

protected

Definition at line 61 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ ConstraintDirichletBoundaryConditionLocalAssembler()

template<typename ShapeFunction , int GlobalDim>

ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::ConstraintDirichletBoundaryConditionLocalAssembler	(	MeshLib::Element const &	surface_element,
		std::size_t	,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		MeshLib::Mesh const &	bulk_mesh,
		std::vector< std::pair< std::size_t, unsigned > >	bulk_ids )

inline

Precomputes the shape matrices for a given surface element.

Parameters

surface_element	The surface element used for precomputing the
is_axially_symmetric	Corrects integration measure for cylinder coordinates.
integration_method	The integration method used.
bulk_mesh	The bulk mesh the process is defined on.
bulk_ids	Pairs of bulk element ids and bulk element face ids.

Definition at line 73 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

        : _surface_element(surface_element),
          _integration_method(integration_method),
          _bulk_element_id(bulk_ids[_surface_element.getID()].first),
          _surface_element_normal(MeshLib::calculateNormalizedSurfaceNormal(
              _surface_element, *(bulk_mesh.getElements()[_bulk_element_id])))
    {
        auto const shape_matrices =
            NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                      GlobalDim, NumLib::ShapeMatrixType::N_J>(
                _surface_element, is_axially_symmetric, _integration_method);
 
        auto const bulk_face_id = bulk_ids[_surface_element.getID()].second;
        auto const& bulk_element = *bulk_mesh.getElement(_bulk_element_id);
 
        auto const n_integration_points =
            _integration_method.getNumberOfPoints();
        _ip_data.reserve(n_integration_points);
 
        for (unsigned ip = 0; ip < n_integration_points; ++ip)
        {
            auto const& wp = _integration_method.getWeightedPoint(ip);
            auto bulk_element_point =
                MeshLib::getBulkElementPoint(bulk_element, bulk_face_id, wp);
            _ip_data.emplace_back(shape_matrices[ip].detJ,
                                  shape_matrices[ip].integralMeasure,
                                  wp.getWeight(),
                                  std::move(bulk_element_point));
        }
    }

References ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::_bulk_element_id, ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::_integration_method, ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::_ip_data, ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::_surface_element, MeshLib::getBulkElementPoint(), MeshLib::Mesh::getElement(), MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), and NumLib::N_J.

Member Function Documentation

◆ integrate()

template<typename ShapeFunction , int GlobalDim>

double ProcessLib::ConstraintDirichletBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::integrate	(	std::vector< GlobalVector * > const &	x,
		double const	t,
		std::function< Eigen::Vector3d(std::size_t const, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &)> const &	getFlux )

inlineoverridevirtual

Integration for the element with the id element_id.

Parameters

x	The global vector containing the values for numerical integration.
t	The point in time the the integration will be performed.
getFlux	The function of the constraining process used to calculate the flux.

Implements ProcessLib::ConstraintDirichletBoundaryConditionLocalAssemblerInterface.

Definition at line 115 of file ConstraintDirichletBoundaryConditionLocalAssembler.h.

    {
        auto const n_integration_points =
            _integration_method.getNumberOfPoints();
        // integrated_value +=
        //   int_{\Gamma_e} \alpha \cdot flux \cdot normal \d \Gamma
        double integrated_value = 0;
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const bulk_flux = getFlux(
                _bulk_element_id, _ip_data[ip].bulk_element_point, t, x);
 
            // TODO find solution for 2d case
            double const bulk_grad_times_normal(
                Eigen::Map<Eigen::RowVectorXd const>(bulk_flux.data(),
                                                     bulk_flux.size())
                    .dot(_surface_element_normal));
 
            integrated_value +=
                bulk_grad_times_normal *
                _ip_data[ip].detJ_times_integralMeasure_times_weight;
        }
        return integrated_value;
    }