Detailed Description

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

Definition at line 49 of file SteadyStateDiffusionFEM.h.

#include <SteadyStateDiffusionFEM.h>

Inheritance diagram for ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >:

[legend]

Collaboration diagram for ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >:

[legend]

Public Member Functions
	LocalAssemblerData (MeshLib::Element const &element, std::size_t const, bool is_axially_symmetric, unsigned const integration_order, SteadyStateDiffusionData const &process_data)

void	assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &, std::vector< double > &, std::vector< double > &local_K_data, std::vector< double > &) override

Eigen::Vector3d	getFlux (MathLib::Point3d const &p_local_coords, double const t, std::vector< double > const &local_x) const override

Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const override
	Provides the shape matrix at the given integration point. More...

std::vector< double > const &	getIntPtDarcyVelocity (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override

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

void	setInitialConditions (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, bool const use_monolithic_scheme, 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	assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_xdot, 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_xdot, const double dxdot_dx, const double dx_dx, 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_xdot, const double dxdot_dx, const double dx_dx, 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_dot, 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, double const t, double const dt)

void	postNonLinearSolver (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, GlobalVector const &xdot, double const t, double const dt, bool const use_monolithic_scheme, int const process_id)

virtual std::vector< double >	interpolateNodalValuesToIntegrationPoints (std::vector< double > const &)

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

Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual	~ExtrapolatableElement ()=default

Private Types
using	ShapeMatricesType = ShapeMatrixPolicyType< ShapeFunction, GlobalDim >

using	ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

using	LocalAssemblerTraits = ProcessLib::LocalAssemblerTraits< ShapeMatricesType, ShapeFunction::NPOINTS, NUM_NODAL_DOF, GlobalDim >

using	NodalMatrixType = typename LocalAssemblerTraits::LocalMatrix

using	NodalVectorType = typename LocalAssemblerTraits::LocalVector

using	GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType

Private Attributes
MeshLib::Element const &	_element

SteadyStateDiffusionData const &	_process_data

IntegrationMethod const	_integration_method

std::vector< ShapeMatrices, Eigen::aligned_allocator< ShapeMatrices > >	_shape_matrices

Member Typedef Documentation

◆ GlobalDimVectorType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType

private

Definition at line 59 of file SteadyStateDiffusionFEM.h.

◆ LocalAssemblerTraits

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::LocalAssemblerTraits = ProcessLib::LocalAssemblerTraits< ShapeMatricesType, ShapeFunction::NPOINTS, NUM_NODAL_DOF, GlobalDim>

private

Definition at line 54 of file SteadyStateDiffusionFEM.h.

◆ NodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::NodalMatrixType = typename LocalAssemblerTraits::LocalMatrix

private

Definition at line 57 of file SteadyStateDiffusionFEM.h.

◆ NodalVectorType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::NodalVectorType = typename LocalAssemblerTraits::LocalVector

private

Definition at line 58 of file SteadyStateDiffusionFEM.h.

◆ ShapeMatrices

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

private

Definition at line 52 of file SteadyStateDiffusionFEM.h.

◆ ShapeMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

using ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

private

Definition at line 51 of file SteadyStateDiffusionFEM.h.

Constructor & Destructor Documentation

◆ LocalAssemblerData()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::LocalAssemblerData	(	MeshLib::Element const &	element,
		std::size_t const	,
		bool	is_axially_symmetric,
		unsigned const	integration_order,
		SteadyStateDiffusionData const &	process_data
	)

inline

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

Definition at line 64 of file SteadyStateDiffusionFEM.h.

         : _element(element),
           _process_data(process_data),
           _integration_method(integration_order),
           _shape_matrices(
               NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                         GlobalDim>(
                   element, is_axially_symmetric, _integration_method))
     {
     }

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

void ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble	(	double const	t,
		double const	dt,
		std::vector< double > const &	local_x,
		std::vector< double > const &	,
		std::vector< double > &	,
		std::vector< double > &	local_K_data,
		std::vector< double > &
	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 79 of file SteadyStateDiffusionFEM.h.

     {
         auto const local_matrix_size = local_x.size();
         // This assertion is valid only if all nodal d.o.f. use the same shape
         // matrices.
         assert(local_matrix_size == ShapeFunction::NPOINTS * NUM_NODAL_DOF);
  
         auto local_K = MathLib::createZeroedMatrix<NodalMatrixType>(
             local_K_data, local_matrix_size, local_matrix_size);
  
         unsigned const n_integration_points =
             _integration_method.getNumberOfPoints();
  
         ParameterLib::SpatialPosition pos;
         pos.setElementID(_element.getID());
  
         auto const& medium =
             *_process_data.media_map->getMedium(_element.getID());
         MaterialPropertyLib::VariableArray vars;
         vars[static_cast<int>(MaterialPropertyLib::Variable::temperature)] =
             medium
                 .property(
                     MaterialPropertyLib::PropertyType::reference_temperature)
                 .template value<double>(vars, pos, t, dt);
  
         for (unsigned ip = 0; ip < n_integration_points; ip++)
         {
             pos.setIntegrationPoint(ip);
             auto const& sm = _shape_matrices[ip];
             auto const& wp = _integration_method.getWeightedPoint(ip);
  
             double p_int_pt = 0.0;
             NumLib::shapeFunctionInterpolate(local_x, sm.N, p_int_pt);
             vars[static_cast<int>(
                 MaterialPropertyLib::Variable::phase_pressure)] = p_int_pt;
             auto const k = MaterialPropertyLib::formEigenTensor<GlobalDim>(
                 medium.property(MaterialPropertyLib::PropertyType::diffusion)
                     .value(vars, pos, t, dt));
  
             local_K.noalias() += sm.dNdx.transpose() * k * sm.dNdx * sm.detJ *
                                  sm.integralMeasure * wp.getWeight();
         }
     }

References ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_element, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_integration_method, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_process_data, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_shape_matrices, MaterialPropertyLib::diffusion, MeshLib::Element::getID(), ProcessLib::SteadyStateDiffusion::SteadyStateDiffusionData::media_map, ProcessLib::SteadyStateDiffusion::NUM_NODAL_DOF, MaterialPropertyLib::phase_pressure, MaterialPropertyLib::reference_temperature, ParameterLib::SpatialPosition::setElementID(), ParameterLib::SpatialPosition::setIntegrationPoint(), NumLib::shapeFunctionInterpolate(), and MaterialPropertyLib::temperature.

◆ getFlux()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

Eigen::Vector3d ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getFlux	(	MathLib::Point3d const &	p_local_coords,
		double const	t,
		std::vector< double > const &	local_x
	)		const

inlineoverridevirtual

Computes the flux in the point p_local_coords that is given in local coordinates using the values from local_x.

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 130 of file SteadyStateDiffusionFEM.h.

     {
         // TODO (tf) Temporary value not used by current material models. Need
         // extension of getFlux interface
         double const dt = std::numeric_limits<double>::quiet_NaN();
  
         // Eval shape matrices at given point
         // Note: Axial symmetry is set to false here, because we only need dNdx
         // here, which is not affected by axial symmetry.
         auto const shape_matrices =
             NumLib::computeShapeMatrices<ShapeFunction, ShapeMatricesType,
                                          GlobalDim>(
                 _element, false /*is_axially_symmetric*/,
                 std::array{p_local_coords})[0];
  
         // fetch hydraulic conductivity
         ParameterLib::SpatialPosition pos;
         pos.setElementID(_element.getID());
         auto const& medium =
             *_process_data.media_map->getMedium(_element.getID());
  
         MaterialPropertyLib::VariableArray vars;
         vars[static_cast<int>(MaterialPropertyLib::Variable::temperature)] =
             medium
                 .property(
                     MaterialPropertyLib::PropertyType::reference_temperature)
                 .template value<double>(vars, pos, t, dt);
         double pressure = 0.0;
         NumLib::shapeFunctionInterpolate(local_x, shape_matrices.N, pressure);
         vars[static_cast<int>(MaterialPropertyLib::Variable::phase_pressure)] =
             pressure;
  
         auto const k = MaterialPropertyLib::formEigenTensor<GlobalDim>(
             medium.property(MaterialPropertyLib::PropertyType::diffusion)
                 .value(vars, pos, t, dt));
  
         Eigen::Vector3d flux(0.0, 0.0, 0.0);
         flux.head<GlobalDim>() =
             -k * shape_matrices.dNdx *
             Eigen::Map<const NodalVectorType>(local_x.data(), local_x.size());
  
         return flux;
     }

References ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_element, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_process_data, NumLib::computeShapeMatrices(), MaterialPropertyLib::diffusion, MeshLib::Element::getID(), ProcessLib::SteadyStateDiffusion::SteadyStateDiffusionData::media_map, MaterialPropertyLib::phase_pressure, MaterialPropertyLib::reference_temperature, ParameterLib::SpatialPosition::setElementID(), NumLib::shapeFunctionInterpolate(), and MaterialPropertyLib::temperature.

◆ getIntPtDarcyVelocity()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

std::vector<double> const& ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getIntPtDarcyVelocity	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_table,
		std::vector< double > &	cache
	)		const

inlineoverridevirtual

Implements ProcessLib::SteadyStateDiffusion::SteadyStateDiffusionLocalAssemblerInterface.

Definition at line 185 of file SteadyStateDiffusionFEM.h.

     {
         // TODO (tf) Temporary value not used by current material models. Need
         // extension of secondary variable interface.
         double const dt = std::numeric_limits<double>::quiet_NaN();
  
         auto const n_integration_points =
             _integration_method.getNumberOfPoints();
  
         int const process_id = 0;  // monolithic scheme
         auto const indices =
             NumLib::getIndices(_element.getID(), *dof_table[process_id]);
         assert(!indices.empty());
         auto const local_x = x[process_id]->get(indices);
         auto const local_x_vec =
             MathLib::toVector<Eigen::Matrix<double, ShapeFunction::NPOINTS, 1>>(
                 local_x, ShapeFunction::NPOINTS);
  
         cache.clear();
         auto cache_mat = MathLib::createZeroedMatrix<
             Eigen::Matrix<double, GlobalDim, Eigen::Dynamic, Eigen::RowMajor>>(
             cache, GlobalDim, n_integration_points);
  
         ParameterLib::SpatialPosition pos;
         pos.setElementID(_element.getID());
  
         auto const& medium =
             *_process_data.media_map->getMedium(_element.getID());
  
         MaterialPropertyLib::VariableArray vars;
         vars[static_cast<int>(MaterialPropertyLib::Variable::temperature)] =
             medium
                 .property(
                     MaterialPropertyLib::PropertyType::reference_temperature)
                 .template value<double>(vars, pos, t, dt);
         double pressure = 0.0;
         for (unsigned i = 0; i < n_integration_points; ++i)
         {
             pos.setIntegrationPoint(i);
             NumLib::shapeFunctionInterpolate(local_x, _shape_matrices[i].N,
                                              pressure);
             vars[static_cast<int>(
                 MaterialPropertyLib::Variable::phase_pressure)] = pressure;
  
             auto const k = MaterialPropertyLib::formEigenTensor<GlobalDim>(
                 medium.property(MaterialPropertyLib::PropertyType::diffusion)
                     .value(vars, pos, t, dt));
             // dimensions: (d x 1) = (d x n) * (n x 1)
             cache_mat.col(i).noalias() =
                 -k * _shape_matrices[i].dNdx * local_x_vec;
         }
  
         return cache;
     }

References ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_element, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_integration_method, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_process_data, ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_shape_matrices, MathLib::createZeroedMatrix(), MaterialPropertyLib::diffusion, MeshLib::Element::getID(), NumLib::getIndices(), ProcessLib::SteadyStateDiffusion::SteadyStateDiffusionData::media_map, MaterialPropertyLib::phase_pressure, MaterialPropertyLib::reference_temperature, ParameterLib::SpatialPosition::setElementID(), NumLib::shapeFunctionInterpolate(), and MaterialPropertyLib::temperature.

◆ getShapeMatrix()

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

Eigen::Map<const Eigen::RowVectorXd> ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getShapeMatrix ( const unsigned integration_point ) const

inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 176 of file SteadyStateDiffusionFEM.h.

     {
         auto const& N = _shape_matrices[integration_point].N;
  
         // assumes N is stored contiguously in memory
         return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
     }

References ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_shape_matrices.

Member Data Documentation

◆ _element

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

MeshLib::Element const& ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_element

private

Definition at line 245 of file SteadyStateDiffusionFEM.h.

Referenced by ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble(), ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getFlux(), and ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getIntPtDarcyVelocity().

◆ _integration_method

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

IntegrationMethod const ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_integration_method

private

Definition at line 248 of file SteadyStateDiffusionFEM.h.

Referenced by ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble(), and ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getIntPtDarcyVelocity().

◆ _process_data

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

SteadyStateDiffusionData const& ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_process_data

private

Definition at line 246 of file SteadyStateDiffusionFEM.h.

Referenced by ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble(), ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getFlux(), and ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getIntPtDarcyVelocity().

◆ _shape_matrices

template<typename ShapeFunction , typename IntegrationMethod , int GlobalDim>

std::vector<ShapeMatrices, Eigen::aligned_allocator<ShapeMatrices> > ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::_shape_matrices

private

Definition at line 250 of file SteadyStateDiffusionFEM.h.

Referenced by ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble(), ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getIntPtDarcyVelocity(), and ProcessLib::SteadyStateDiffusion::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::getShapeMatrix().

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

ProcessLib/SteadyStateDiffusion/SteadyStateDiffusionFEM.h

Detailed Description

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

Public Member Functions

Private Types

Private Attributes

Member Typedef Documentation

◆ GlobalDimVectorType

◆ LocalAssemblerTraits

◆ NodalMatrixType

◆ NodalVectorType

◆ ShapeMatrices

◆ ShapeMatricesType

Constructor & Destructor Documentation

◆ LocalAssemblerData()

Member Function Documentation

◆ assemble()

◆ getFlux()

◆ getIntPtDarcyVelocity()

◆ getShapeMatrix()

Member Data Documentation

◆ _element

◆ _integration_method

◆ _process_data

◆ _shape_matrices

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