Detailed Description

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

Definition at line 33 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

#include <VariableDependentNeumannBoundaryConditionLocalAssembler.h>

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

[legend]

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

[legend]

Public Member Functions
	VariableDependentNeumannBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const local_matrix_size, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, VariableDependentNeumannBoundaryConditionData const &data)
void	assemble (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, double const t, std::vector< GlobalVector * > const &x, int const process_id, 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

Private Types
using	Base
using	NodalVectorType = typename Base::NodalVectorType
using	NodalMatrixType = typename Base::NodalMatrixType

Private Attributes
VariableDependentNeumannBoundaryConditionData const &	_data
std::size_t const	_local_matrix_size

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::VariableDependentNeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::Base

private

Initial value:

GenericNaturalBoundaryConditionLocalAssembler<ShapeFunction, GlobalDim>

ProcessLib::GenericNaturalBoundaryConditionLocalAssembler::GenericNaturalBoundaryConditionLocalAssembler

GenericNaturalBoundaryConditionLocalAssembler(MeshLib::Element const &e, bool is_axially_symmetric, NumLib::GenericIntegrationMethod const &integration_method)

Definition GenericNaturalBoundaryConditionLocalAssembler.h:79

Definition at line 37 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

◆ NodalMatrixType

template<typename ShapeFunction, int GlobalDim>

using ProcessLib::VariableDependentNeumannBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::NodalMatrixType = typename Base::NodalMatrixType

private

Definition at line 40 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunction, int GlobalDim>

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

private

Definition at line 39 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ VariableDependentNeumannBoundaryConditionLocalAssembler()

template<typename ShapeFunction, int GlobalDim>

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

inline

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

Definition at line 45 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

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

References _data, and _local_matrix_size.

Member Function Documentation

◆ assemble()

template<typename ShapeFunction, int GlobalDim>

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

inlineoverridevirtual

Implements ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface.

Definition at line 57 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

    {
        NodalVectorType _local_rhs(_local_matrix_size);
        _local_rhs.setZero();
        // Get element nodes for the interpolation from nodes to
        // integration point.
        NodalVectorType const constant_node_values =
            _data.constant.getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
        NodalVectorType const coefficient_current_variable_node_values =
            _data.coefficient_current_variable
                .getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
        NodalVectorType const coefficient_other_variable_node_values =
            _data.coefficient_other_variable
                .getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
        NodalVectorType const coefficient_mixed_variables_node_values =
            _data.coefficient_mixed_variables
                .getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
        unsigned const n_integration_points =
            Base::_integration_method.getNumberOfPoints();
 
        auto const indices_current_variable =
            NumLib::getIndices(mesh_item_id, dof_table_boundary);
        auto const indices_other_variable = NumLib::getIndices(
            mesh_item_id, *_data.dof_table_boundary_other_variable);
        std::vector<double> const local_current_variable =
            x[process_id]->get(indices_current_variable);
        std::vector<double> const local_other_variable =
            x[process_id]->get(indices_other_variable);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& n_and_weight = Base::_ns_and_weights[ip];
            auto const& N = n_and_weight.N;
            auto const& w = n_and_weight.weight;
 
            double current_variable_int_pt = 0.0;
            double other_variable_int_pt = 0.0;
 
            NumLib::shapeFunctionInterpolate(local_current_variable, N,
                                             current_variable_int_pt);
            NumLib::shapeFunctionInterpolate(local_other_variable, N,
                                             other_variable_int_pt);
            NodalVectorType const neumann_node_values =
                constant_node_values +
                coefficient_current_variable_node_values *
                    current_variable_int_pt +
                coefficient_other_variable_node_values * other_variable_int_pt +
                coefficient_mixed_variables_node_values *
                    current_variable_int_pt * other_variable_int_pt;
            _local_rhs.noalias() += N * neumann_node_values.dot(N) * w;
        }
 
        b.add(indices_current_variable, _local_rhs);
    }