Detailed Description

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

Definition at line 31 of file NeumannBoundaryConditionLocalAssembler.h.

#include <NeumannBoundaryConditionLocalAssembler.h>

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

[legend]

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

[legend]

Public Member Functions
	NeumannBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const local_matrix_size, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, NeumannBoundaryConditionData const &data)

void	assemble (std::size_t const id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, double const t, std::vector< GlobalVector * > const &, int const, GlobalMatrix &, GlobalVector &b, GlobalMatrix *) override

Public Member Functions inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
	GenericNaturalBoundaryConditionLocalAssembler (MeshLib::Element const &e, bool is_axially_symmetric, NumLib::GenericIntegrationMethod const &integration_method)

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

Public Attributes
	EIGEN_MAKE_ALIGNED_OPERATOR_NEW

Private Types
using	Base

using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

using	NodalVectorType = typename Base::NodalVectorType

Private Attributes
NeumannBoundaryConditionData const &	_data

NodalVectorType	_local_rhs

Additional Inherited Members
Protected Types inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

using	NodalVectorType = typename ShapeMatricesType::NodalVectorType

Protected Attributes inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
NumLib::GenericIntegrationMethod const &	_integration_method

std::vector< NAndWeight, Eigen::aligned_allocator< NAndWeight > > const	_ns_and_weights

MeshLib::Element const &	_element

Member Typedef Documentation

◆ Base

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::NeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::Base

private

Initial value:

GenericNaturalBoundaryConditionLocalAssembler<ShapeFunction, GlobalDim>

Definition at line 35 of file NeumannBoundaryConditionLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::NeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::NodalVectorType = typename Base::NodalVectorType

private

Definition at line 38 of file NeumannBoundaryConditionLocalAssembler.h.

◆ ShapeMatricesType

template<typename ShapeFunction , int GlobalDim>

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

private

Definition at line 37 of file NeumannBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ NeumannBoundaryConditionLocalAssembler()

template<typename ShapeFunction , int GlobalDim>

ProcessLib::NeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::NeumannBoundaryConditionLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	local_matrix_size,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		NeumannBoundaryConditionData const &	data )

inline

The neumann_bc_value factor is directly integrated into the local element matrix.

Definition at line 43 of file NeumannBoundaryConditionLocalAssembler.h.

        : Base(e, is_axially_symmetric, integration_method),
          _data(data),
          _local_rhs(local_matrix_size)
    {
    }

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , int GlobalDim>

void ProcessLib::NeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::assemble	(	std::size_t const	id,
		NumLib::LocalToGlobalIndexMap const &	dof_table_boundary,
		double const	t,
		std::vector< GlobalVector * > const &	,
		int const	,
		GlobalMatrix &	,
		GlobalVector &	b,
		GlobalMatrix *	)

inlineoverridevirtual

Implements ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface.

Definition at line 55 of file NeumannBoundaryConditionLocalAssembler.h.

    {
        _local_rhs.setZero();
 
        unsigned const n_integration_points =
            Base::_integration_method.getNumberOfPoints();
 
        NodalVectorType parameter_node_values;
        if (typeid(ParameterLib::MeshNodeParameter<double> const) ==
            typeid(_data.neumann_bc_parameter))
        {
            // Get element nodes for the interpolation from nodes to integration
            // point.
            parameter_node_values =
                _data.neumann_bc_parameter
                    .getNodalValuesOnElement(Base::_element, t)
                    .template topRows<
                        ShapeFunction::MeshElement::n_all_nodes>();
        }
 
        double integral_measure = 1.0;
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& ip_data = Base::_ns_and_weights[ip];
            auto const& N = ip_data.N;
            auto const& w = ip_data.weight;
 
            ParameterLib::SpatialPosition const position{
                std::nullopt, Base::_element.getID(), ip,
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(
                        Base::_element, N))};
 
            if (_data.integral_measure)
            {
                integral_measure = (*_data.integral_measure)(t, position)[0];
            }
            if (typeid(ParameterLib::MeshNodeParameter<double> const) ==
                typeid(_data.neumann_bc_parameter))
            {
                _local_rhs.noalias() +=
                    N * parameter_node_values.dot(N) * w * integral_measure;
            }
            else
            {
                auto const value = _data.neumann_bc_parameter(t, position)[0];
                _local_rhs.noalias() += N * value * w * integral_measure;
            }
        }
 
        auto const indices = NumLib::getIndices(id, dof_table_boundary);
        b.add(indices, _local_rhs);
    }