Detailed Description

template<int DisplacementDim>
struct ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >

Definition at line 29 of file LocalAssemblerInterface.h.

#include <LocalAssemblerInterface.h>

Inheritance diagram for ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >:

Collaboration diagram for ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >:

Public Member Functions
	LocalAssemblerInterface (MeshLib::Element const &e, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, TH2MProcessData< DisplacementDim > &process_data)

virtual std::size_t	setIPDataInitialConditions (std::string_view name, double const *values, int const integration_order)=0

virtual std::vector< double > const &	getIntPtDarcyVelocityGas (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtDarcyVelocityLiquid (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtDiffusionVelocityVapourGas (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtDiffusionVelocityGasGas (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtDiffusionVelocitySoluteLiquid (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtDiffusionVelocityLiquidLiquid (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

virtual std::vector< double > const &	getIntPtEnthalpySolid (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0

unsigned	getNumberOfIntegrationPoints () const

int	getMaterialID () const

std::vector< double >	getMaterialStateVariableInternalState (std::function< std::span< double >(typename ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim >::MaterialStateVariables &)> const &get_values_span, int const &n_components) const

ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim >::MaterialStateVariables const &	getMaterialStateVariablesAt (unsigned integration_point) const

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

virtual void	setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, int const process_id)

virtual void	initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)

virtual void	preAssemble (double const, double const, std::vector< double > const &)

virtual void	assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)

virtual void	assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)

virtual void	assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)

virtual void	assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)

virtual void	computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)

virtual void	preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)

virtual void	postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)

void	postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)

virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const

virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double > > const &) const
	Fits to staggered scheme.

Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const =0
	Provides the shape matrix at the given integration point.

virtual	~ExtrapolatableElement ()=default

Static Public Member Functions
static auto	getReflectionDataForOutput ()

Public Attributes
TH2MProcessData< DisplacementDim > &	process_data_

std::vector< typename ConstitutiveRelations::StatefulData< DisplacementDim > >	current_states_

std::vector< typename ConstitutiveRelations::StatefulDataPrev< DisplacementDim > >	prev_states_

std::vector< ConstitutiveRelations::MaterialStateData< DisplacementDim > >	material_states_

std::vector< ConstitutiveRelations::OutputData< DisplacementDim > >	output_data_

NumLib::GenericIntegrationMethod const &	integration_method_

MeshLib::Element const &	element_

bool const	is_axially_symmetric_

ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim > const &	solid_material_

Constructor & Destructor Documentation

◆ LocalAssemblerInterface()

template<int DisplacementDim>

ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::LocalAssemblerInterface	(	MeshLib::Element const &	e,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		TH2MProcessData< DisplacementDim > &	process_data )

inline

Definition at line 32 of file LocalAssemblerInterface.h.

        : process_data_(process_data),
          integration_method_(integration_method),
          element_(e),
          is_axially_symmetric_(is_axially_symmetric),
          solid_material_(MaterialLib::Solids::selectSolidConstitutiveRelation(
              process_data_.solid_materials, process_data_.material_ids,
              element_.getID()))
    {
        unsigned const n_integration_points =
            integration_method_.getNumberOfPoints();
 
        current_states_.resize(n_integration_points);
        prev_states_.resize(n_integration_points);
        output_data_.resize(n_integration_points);
 
        material_states_.reserve(n_integration_points);
        for (unsigned ip = 0; ip < n_integration_points; ++ip)
        {
            material_states_.emplace_back(
                solid_material_.createMaterialStateVariables());
 
            // Set initial stress field to zero. Might be overwritten by
            // integration point data or initial stress.
            current_states_[ip].eff_stress_data.sigma =
                ConstitutiveRelations::KelvinVector<DisplacementDim>::Zero();
        }
    }

References ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::current_states_, NumLib::GenericIntegrationMethod::getNumberOfPoints(), ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::integration_method_, ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::material_states_, ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::output_data_, ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::prev_states_, and ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::solid_material_.

Member Function Documentation

◆ getIntPtDarcyVelocityGas()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDarcyVelocityGas	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtDarcyVelocityLiquid()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDarcyVelocityLiquid	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtDiffusionVelocityGasGas()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDiffusionVelocityGasGas	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtDiffusionVelocityLiquidLiquid()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDiffusionVelocityLiquidLiquid	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtDiffusionVelocitySoluteLiquid()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDiffusionVelocitySoluteLiquid	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtDiffusionVelocityVapourGas()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtDiffusionVelocityVapourGas	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getIntPtEnthalpySolid()

template<int DisplacementDim>

virtual std::vector< double > const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getIntPtEnthalpySolid	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache ) const

pure virtual

Implemented in ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >.

◆ getMaterialID()

template<int DisplacementDim>

int ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getMaterialID ( ) const

inline

Definition at line 117 of file LocalAssemblerInterface.h.

    {
        return process_data_.material_ids == nullptr
                   ? 0
                   : (*process_data_.material_ids)[element_.getID()];
    }

References ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::element_, MeshLib::Element::getID(), and ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::process_data_.

◆ getMaterialStateVariableInternalState()

template<int DisplacementDim>

std::vector< double > ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getMaterialStateVariableInternalState	(	std::function< std::span< double >(typename ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim >::MaterialStateVariables &)> const &	get_values_span,
		int const &	n_components ) const

inline

Definition at line 124 of file LocalAssemblerInterface.h.

    {
        return ProcessLib::getIntegrationPointDataMaterialStateVariables(
            material_states_,
            &ConstitutiveRelations::MaterialStateData<
                DisplacementDim>::material_state_variables,
            get_values_span, n_components);
    }

References ProcessLib::getIntegrationPointDataMaterialStateVariables(), and ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::material_states_.

◆ getMaterialStateVariablesAt()

template<int DisplacementDim>

ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim >::MaterialStateVariables const & ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getMaterialStateVariablesAt ( unsigned integration_point ) const

inline

Definition at line 140 of file LocalAssemblerInterface.h.

    {
        return *material_states_[integration_point].material_state_variables;
    }

References ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::material_states_.

◆ getNumberOfIntegrationPoints()

template<int DisplacementDim>

unsigned ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getNumberOfIntegrationPoints ( ) const

inline

Definition at line 112 of file LocalAssemblerInterface.h.

    {
        return integration_method_.getNumberOfPoints();
    }

References NumLib::GenericIntegrationMethod::getNumberOfPoints(), and ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::integration_method_.

◆ getReflectionDataForOutput()

template<int DisplacementDim>

static auto ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::getReflectionDataForOutput ( )

inlinestatic

Definition at line 145 of file LocalAssemblerInterface.h.

    {
        using Self = LocalAssemblerInterface<DisplacementDim>;
 
        return ProcessLib::Reflection::reflectWithoutName(
            &Self::current_states_, &Self::output_data_);
    }

References ProcessLib::Reflection::reflectWithoutName().

◆ setIPDataInitialConditions()

template<int DisplacementDim>

virtual std::size_t ProcessLib::TH2M::LocalAssemblerInterface< DisplacementDim >::setIPDataInitialConditions	(	std::string_view	name,
		double const *	values,
		int const	integration_order )