Detailed Description

Definition at line 39 of file DeactivatedSubdomainDirichlet.h.

#include <DeactivatedSubdomainDirichlet.h>

Inheritance diagram for ProcessLib::DeactivatedSubdomainDirichlet:

Collaboration diagram for ProcessLib::DeactivatedSubdomainDirichlet:

Public Member Functions
	DeactivatedSubdomainDirichlet (MeshLib::PropertyVector< unsigned char > const &is_active, MathLib::PiecewiseLinearInterpolation time_interval, ParameterLib::Parameter< double > const &parameter, bool const set_outer_nodes_dirichlet_values, DeactivatedSubdomainMesh const &subdomain, NumLib::LocalToGlobalIndexMap const &dof_table_bulk, int const variable_id, int const component_id)

void	getEssentialBCValues (const double t, GlobalVector const &x, NumLib::IndexValueVector< GlobalIndexType > &bc_values) const override
	Writes the values of essential BCs to `bc_values`.

Public Member Functions inherited from ProcessLib::BoundaryCondition
virtual void	applyNaturalBC (const double, std::vector< GlobalVector * > const &, int const, GlobalMatrix &, GlobalVector &, GlobalMatrix *)

virtual void	preTimestep (const double, std::vector< GlobalVector * > const &, int const)

virtual void	postTimestep (const double, std::vector< GlobalVector * > const &, int const)

virtual	~BoundaryCondition ()=default

Private Member Functions
void	config (NumLib::LocalToGlobalIndexMap const &dof_table_bulk)

Private Attributes
ParameterLib::Parameter< double > const &	_parameter

DeactivatedSubdomainMesh const &	_subdomain

std::unique_ptr< NumLib::LocalToGlobalIndexMap const >	_dof_table_boundary

int const	_variable_id

int const	_component_id

MathLib::PiecewiseLinearInterpolation const	_time_interval

MeshLib::PropertyVector< unsigned char > const &	_is_active

bool const	_set_outer_nodes_dirichlet_values

Constructor & Destructor Documentation

◆ DeactivatedSubdomainDirichlet()

ProcessLib::DeactivatedSubdomainDirichlet::DeactivatedSubdomainDirichlet	(	MeshLib::PropertyVector< unsigned char > const &	is_active,
		MathLib::PiecewiseLinearInterpolation	time_interval,
		ParameterLib::Parameter< double > const &	parameter,
		bool const	set_outer_nodes_dirichlet_values,
		DeactivatedSubdomainMesh const &	subdomain,
		NumLib::LocalToGlobalIndexMap const &	dof_table_bulk,
		int const	variable_id,
		int const	component_id )

Definition at line 20 of file DeactivatedSubdomainDirichlet.cpp.

    : _parameter(parameter),
      _subdomain(subdomain),
      _variable_id(variable_id),
      _component_id(component_id),
      _time_interval(std::move(time_interval)),
      _is_active(is_active),
      _set_outer_nodes_dirichlet_values(set_outer_nodes_dirichlet_values)
{
    config(dof_table_bulk);
}

References config().

Member Function Documentation

◆ config()

void ProcessLib::DeactivatedSubdomainDirichlet::config ( NumLib::LocalToGlobalIndexMap const & dof_table_bulk )

private

Definition at line 39 of file DeactivatedSubdomainDirichlet.cpp.

{
    checkParametersOfDirichletBoundaryCondition(_subdomain.mesh, dof_table_bulk,
                                                _variable_id, _component_id);
 
    std::vector<MeshLib::Node*> const& bc_nodes = _subdomain.mesh.getNodes();
    MeshLib::MeshSubset subdomain_mesh_subset(_subdomain.mesh, bc_nodes);
 
    // Create local DOF table from the BC mesh subset for the given variable
    // and component id.
    _dof_table_boundary = dof_table_bulk.deriveBoundaryConstrainedMap(
        _variable_id, {_component_id}, std::move(subdomain_mesh_subset));
}

References _component_id, _dof_table_boundary, _subdomain, _variable_id, ProcessLib::checkParametersOfDirichletBoundaryCondition(), NumLib::LocalToGlobalIndexMap::deriveBoundaryConstrainedMap(), MeshLib::Mesh::getNodes(), and ProcessLib::DeactivatedSubdomainMesh::mesh.

Referenced by DeactivatedSubdomainDirichlet().

◆ getEssentialBCValues()

void ProcessLib::DeactivatedSubdomainDirichlet::getEssentialBCValues	(	const double	,
		GlobalVector const &	,
		NumLib::IndexValueVector< GlobalIndexType > &	) const

overridevirtual

Writes the values of essential BCs to bc_values.

Reimplemented from ProcessLib::BoundaryCondition.

Definition at line 54 of file DeactivatedSubdomainDirichlet.cpp.

{
    [[maybe_unused]] auto const& bulk_node_ids =
        *MeshLib::bulkNodeIDs(_subdomain.mesh);
    [[maybe_unused]] auto const& bulk_element_ids =
        *MeshLib::bulkElementIDs(_subdomain.mesh);
 
    auto is_inactive_id = [&](std::size_t const bulk_element_id)
    { return _is_active[bulk_element_id] == 0; };
 
    auto is_inactive_element = [&](MeshLib::Element const* const e)
    { return is_inactive_id(bulk_element_ids[e->getID()]); };
 
    std::vector<std::size_t> inactive_nodes_in_bc_mesh;
    std::copy_if(begin(_subdomain.inner_nodes), end(_subdomain.inner_nodes),
                 back_inserter(inactive_nodes_in_bc_mesh),
                 [&](std::size_t const n)
                 {
                     const auto& connected_elements =
                         _subdomain.mesh.getElementsConnectedToNode(n);
 
                     return std::all_of(begin(connected_elements),
                                        end(connected_elements),
                                        is_inactive_element);
                 });
 
    if (_set_outer_nodes_dirichlet_values)
    {
        std::copy_if(begin(_subdomain.outer_nodes), end(_subdomain.outer_nodes),
                     back_inserter(inactive_nodes_in_bc_mesh),
                     [&](std::size_t const n)
                     {
                         const auto& connected_elements =
                             _subdomain.mesh.getElementsConnectedToNode(n);
 
                         return std::all_of(begin(connected_elements),
                                            end(connected_elements),
                                            is_inactive_element);
                     });
    }
    else
    {
        for (std::size_t i = 0; i < _subdomain.outer_nodes.size(); ++i)
        {
            auto const& connected_elements = _subdomain.outer_nodes_elements[i];
            if (std::all_of(begin(connected_elements), end(connected_elements),
                            is_inactive_id))
            {
                inactive_nodes_in_bc_mesh.push_back(_subdomain.outer_nodes[i]);
            }
        }
    }
 
    auto time_interval_contains = [&](double const t)
    {
        return _time_interval.getSupportMin() <= t &&
               t <= _time_interval.getSupportMax();
    };
    if (time_interval_contains(t))
    {
        getEssentialBCValuesLocal(
            _parameter, _subdomain.mesh, inactive_nodes_in_bc_mesh,
            *_dof_table_boundary, _variable_id, _component_id, t, x, bc_values);
        return;
    }
 
    bc_values.ids.clear();
    bc_values.values.clear();
}