ProcessLib::SurfaceFlux Class Reference

Detailed Description

Definition at line 19 of file SurfaceFlux.h.

#include <SurfaceFlux.h>

Public Member Functions
	SurfaceFlux (MeshLib::Mesh &boundary_mesh, std::size_t bulk_property_number_of_components, unsigned const integration_order)
void	integrate (std::vector< GlobalVector * > const &x, MeshLib::PropertyVector< double > &balance, double const t, MeshLib::Mesh const &bulk_mesh, std::vector< std::size_t > const &active_element_ids, std::function< Eigen::Vector3d(std::size_t const, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &)> const &getFlux)

Private Attributes
std::vector< std::unique_ptr< SurfaceFluxLocalAssemblerInterface > >	_local_assemblers

Constructor & Destructor Documentation

SurfaceFlux()

ProcessLib::SurfaceFlux::SurfaceFlux	(	MeshLib::Mesh &	boundary_mesh,
		std::size_t	bulk_property_number_of_components,
		unsigned const	integration_order )

Constructs an object that is able to compute a balance.

Parameters

boundary_mesh	The integration boundary mesh.
bulk_property_number_of_components	The number of components the variable has.
integration_order	Integration order used in local assembly.

Definition at line 19 of file SurfaceFlux.cpp.

{
    DBUG("Create local balance assemblers.");
    // Populate the vector of local assemblers.
    _local_assemblers.resize(boundary_mesh.getElements().size());
 
    // needed to create dof table
    auto mesh_subset_all_nodes = std::make_unique<MeshLib::MeshSubset>(
        boundary_mesh, boundary_mesh.getNodes());
 
    // Collect the mesh subsets in a vector.
    std::vector<MeshLib::MeshSubset> all_mesh_subsets;
    std::generate_n(std::back_inserter(all_mesh_subsets),
                    bulk_property_number_of_components,
                    [&]() { return *mesh_subset_all_nodes; });
 
    // needed for creation of local assemblers
    auto dof_table = std::make_unique<NumLib::LocalToGlobalIndexMap const>(
        std::move(all_mesh_subsets), NumLib::ComponentOrder::BY_LOCATION);
 
    auto const bulk_element_ids =
        boundary_mesh.getProperties().template getPropertyVector<std::size_t>(
            MeshLib::getBulkIDString(MeshLib::MeshItemType::Cell),
            MeshLib::MeshItemType::Cell, 1);
    auto const bulk_face_ids =
        boundary_mesh.getProperties().template getPropertyVector<std::size_t>(
            MeshLib::getBulkIDString(MeshLib::MeshItemType::Face),
            MeshLib::MeshItemType::Cell, 1);
 
    ProcessLib::createLocalAssemblers<SurfaceFluxLocalAssembler>(
        boundary_mesh.getDimension() + 1,  // or bulk_mesh.getDimension()?
        boundary_mesh.getElements(), *dof_table, _local_assemblers,
        NumLib::IntegrationOrder{integration_order},
        boundary_mesh.isAxiallySymmetric(), *bulk_element_ids, *bulk_face_ids);
}

References _local_assemblers, NumLib::BY_LOCATION, MeshLib::Cell, ProcessLib::createLocalAssemblers(), DBUG(), MeshLib::Face, MeshLib::getBulkIDString(), MeshLib::Mesh::getDimension(), MeshLib::Mesh::getElements(), MeshLib::Mesh::getNodes(), MeshLib::Mesh::getProperties(), and MeshLib::Mesh::isAxiallySymmetric().

Member Function Documentation

integrate()

void ProcessLib::SurfaceFlux::integrate	(	std::vector< GlobalVector * > const &	x,
		MeshLib::PropertyVector< double > &	balance,
		double const	t,
		MeshLib::Mesh const &	bulk_mesh,
		std::vector< std::size_t > const &	active_element_ids,
		std::function< Eigen::Vector3d(std::size_t const, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &)> const &	getFlux )

Executes for each element of the mesh the local integration procedure.

Parameters

x	The global solution the integration values will be fetched of.
balance	The vector the integration results will be stored in.
t	The balance will be computed at the time t.
bulk_mesh	Stores a reference to the bulk mesh that is needed to fetch the information for the integration over the surface mesh.
active_element_ids	The IDs of active elements. It is empty if there is no deactivated subdomain.
getFlux	function that calculates the flux in the integration points of the face elements of the bulk element that belongs to the surface.

Definition at line 57 of file SurfaceFlux.cpp.

{
    DBUG("Integrate SurfaceFlux.");
 
    GlobalExecutor::executeSelectedMemberOnDereferenced(
        &SurfaceFluxLocalAssemblerInterface::integrate, _local_assemblers,
        active_element_ids, x, balance, t, bulk_mesh, getFlux);
}

References _local_assemblers, DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), and ProcessLib::SurfaceFluxLocalAssemblerInterface::integrate().

Member Data Documentation

_local_assemblers

std::vector<std::unique_ptr<SurfaceFluxLocalAssemblerInterface> > ProcessLib::SurfaceFlux::_local_assemblers

private

Definition at line 54 of file SurfaceFlux.h.

Referenced by SurfaceFlux(), and integrate().

---

The documentation for this class was generated from the following files:

• ProcessLib/SurfaceFlux/SurfaceFlux.h
• ProcessLib/SurfaceFlux/SurfaceFlux.cpp