Detailed Description

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

Definition at line 32 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

#include <VariableDependentNeumannBoundaryConditionLocalAssembler.h>

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

[legend]

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

[legend]

Public Member Functions
	VariableDependentNeumannBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const local_matrix_size, bool const is_axially_symmetric, unsigned const integration_order, 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, IntegrationMethod, GlobalDim >
	GenericNaturalBoundaryConditionLocalAssembler (MeshLib::Element const &e, bool is_axially_symmetric, unsigned const integration_order)

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

Private Types
using	Base = GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, GlobalDim >

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, IntegrationMethod, GlobalDim >
using	ShapeMatricesType = ShapeMatrixPolicyType< ShapeFunction, GlobalDim >

using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

using	NodalVectorType = typename ShapeMatricesType::NodalVectorType

Protected Attributes inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, GlobalDim >
IntegrationMethod 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 , typename IntegrationMethod , int GlobalDim>

using ProcessLib::VariableDependentNeumannBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, GlobalDim >::Base = GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, GlobalDim>

private

Definition at line 36 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

◆ NodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

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

private

Definition at line 39 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

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

private

Definition at line 38 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ VariableDependentNeumannBoundaryConditionLocalAssembler()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

ProcessLib::VariableDependentNeumannBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, GlobalDim >::VariableDependentNeumannBoundaryConditionLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	local_matrix_size,
		bool const	is_axially_symmetric,
		unsigned const	integration_order,
		VariableDependentNeumannBoundaryConditionData const &	data
	)

inline

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

Definition at line 44 of file VariableDependentNeumannBoundaryConditionLocalAssembler.h.

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

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

void ProcessLib::VariableDependentNeumannBoundaryConditionLocalAssembler< ShapeFunction, IntegrationMethod, 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 56 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);
     }