Detailed Description

Definition at line 62 of file ComponentTransportFEM.h.

#include <ComponentTransportFEM.h>

Inheritance diagram for ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface:

Collaboration diagram for ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface:

Public Member Functions
	ComponentTransportLocalAssemblerInterface ()=default

void	setStaggeredCoupledSolutions (std::size_t const, CoupledSolutionsForStaggeredScheme *const coupling_term)

virtual void	setChemicalSystemID (std::size_t const)=0

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

void	setChemicalSystem (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	assembleReactionEquation (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, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, int const process_id)

virtual void	postSpeciationCalculation (std::size_t const ele_id, double const t, double const dt)=0

virtual 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 =0

virtual std::vector< double > const &	getInterpolatedLocalSolution (const double, std::vector< GlobalVector * > const &int_pt_x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const =0

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	assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_xdot, 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_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< double > const &) 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 Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const =0
	Provides the shape matrix at the given integration point. More...

virtual	~ExtrapolatableElement ()=default

Protected Attributes
CoupledSolutionsForStaggeredScheme *	_coupled_solutions {nullptr}

Private Member Functions
virtual void	initializeChemicalSystemConcrete (Eigen::VectorXd const &, double const)=0

virtual void	setChemicalSystemConcrete (Eigen::VectorXd const &, double const, double const)=0

virtual void	assembleReactionEquationConcrete (double const t, double const dt, Eigen::VectorXd const &local_x, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, int const transport_process_id)=0

Constructor & Destructor Documentation

◆ ComponentTransportLocalAssemblerInterface()

ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::ComponentTransportLocalAssemblerInterface ( )

default

Member Function Documentation

◆ assembleReactionEquation()

void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::assembleReactionEquation	(	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,
		GlobalMatrix &	M,
		GlobalMatrix &	K,
		GlobalVector &	b,
		int const	process_id
	)

inline

Definition at line 124 of file ComponentTransportFEM.h.

     {
         std::vector<double> local_x_vec;
  
         auto const n_processes = x.size();
         for (std::size_t pcs_id = 0; pcs_id < n_processes; ++pcs_id)
         {
             auto const indices =
                 NumLib::getIndices(mesh_item_id, *dof_tables[pcs_id]);
             assert(!indices.empty());
             auto const local_solution = x[pcs_id]->get(indices);
             local_x_vec.insert(std::end(local_x_vec),
                                std::begin(local_solution),
                                std::end(local_solution));
         }
         auto const local_x = MathLib::toVector(local_x_vec);
  
         auto const indices =
             NumLib::getIndices(mesh_item_id, *dof_tables[process_id]);
         auto const num_r_c = indices.size();
  
         std::vector<double> local_M_data;
         local_M_data.reserve(num_r_c * num_r_c);
         std::vector<double> local_K_data;
         local_K_data.reserve(num_r_c * num_r_c);
         std::vector<double> local_b_data;
         local_b_data.reserve(num_r_c);
  
         assembleReactionEquationConcrete(t, dt, local_x, local_M_data,
                                          local_K_data, local_b_data,
                                          process_id);
  
         auto const r_c_indices =
             NumLib::LocalToGlobalIndexMap::RowColumnIndices(indices, indices);
         if (!local_M_data.empty())
         {
             auto const local_M =
                 MathLib::toMatrix(local_M_data, num_r_c, num_r_c);
             M.add(r_c_indices, local_M);
         }
         if (!local_K_data.empty())
         {
             auto const local_K =
                 MathLib::toMatrix(local_K_data, num_r_c, num_r_c);
             K.add(r_c_indices, local_K);
         }
         if (!local_b_data.empty())
         {
             b.add(indices, local_b_data);
         }
     }

References MathLib::EigenMatrix::add(), MathLib::EigenVector::add(), assembleReactionEquationConcrete(), NumLib::getIndices(), MathLib::toMatrix(), and MathLib::toVector().

Referenced by ProcessLib::ComponentTransport::ComponentTransportProcess::solveReactionEquation().

◆ assembleReactionEquationConcrete()

virtual void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::assembleReactionEquationConcrete	(	double const	t,
		double const	dt,
		Eigen::VectorXd const &	local_x,
		std::vector< double > &	local_M_data,
		std::vector< double > &	local_K_data,
		std::vector< double > &	local_b_data,
		int const	transport_process_id
	)

privatepure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

Referenced by assembleReactionEquation().

◆ getInterpolatedLocalSolution()

virtual std::vector<double> const& ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::getInterpolatedLocalSolution	(	const double	,
		std::vector< GlobalVector * > const &	int_pt_x,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	,
		std::vector< double > &	cache
	)		const

pure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

◆ getIntPtDarcyVelocity()

virtual std::vector<double> const& ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::getIntPtDarcyVelocity	(	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::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

Referenced by ProcessLib::ComponentTransport::ComponentTransportProcess::initializeConcreteProcess().

◆ initializeChemicalSystem()

void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::initializeChemicalSystem	(	std::size_t const	mesh_item_id,
		std::vector< NumLib::LocalToGlobalIndexMap const * > const &	dof_tables,
		std::vector< GlobalVector * > const &	x,
		double const	t
	)

inline

Definition at line 78 of file ComponentTransportFEM.h.

     {
         std::vector<double> local_x_vec;
  
         auto const n_processes = x.size();
         for (std::size_t process_id = 0; process_id < n_processes; ++process_id)
         {
             auto const indices =
                 NumLib::getIndices(mesh_item_id, *dof_tables[process_id]);
             assert(!indices.empty());
             auto const local_solution = x[process_id]->get(indices);
             local_x_vec.insert(std::end(local_x_vec),
                                std::begin(local_solution),
                                std::end(local_solution));
         }
         auto const local_x = MathLib::toVector(local_x_vec);
  
         initializeChemicalSystemConcrete(local_x, t);
     }

References NumLib::getIndices(), initializeChemicalSystemConcrete(), and MathLib::toVector().

Referenced by ProcessLib::ComponentTransport::ComponentTransportProcess::setInitialConditionsConcreteProcess().

◆ initializeChemicalSystemConcrete()

virtual void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::initializeChemicalSystemConcrete	(	Eigen::VectorXd const &	,
		double const
	)

privatepure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

Referenced by initializeChemicalSystem().

◆ postSpeciationCalculation()

virtual void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::postSpeciationCalculation	(	std::size_t const	ele_id,
		double const	t,
		double const	dt
	)

pure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

◆ setChemicalSystem()

void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::setChemicalSystem	(	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
	)

inline

Definition at line 101 of file ComponentTransportFEM.h.

     {
         std::vector<double> local_x_vec;
  
         auto const n_processes = x.size();
         for (std::size_t process_id = 0; process_id < n_processes; ++process_id)
         {
             auto const indices =
                 NumLib::getIndices(mesh_item_id, *dof_tables[process_id]);
             assert(!indices.empty());
             auto const local_solution = x[process_id]->get(indices);
             local_x_vec.insert(std::end(local_x_vec),
                                std::begin(local_solution),
                                std::end(local_solution));
         }
         auto const local_x = MathLib::toVector(local_x_vec);
  
         setChemicalSystemConcrete(local_x, t, dt);
     }

References NumLib::getIndices(), setChemicalSystemConcrete(), and MathLib::toVector().

Referenced by ProcessLib::ComponentTransport::ComponentTransportProcess::solveReactionEquation().

◆ setChemicalSystemConcrete()

virtual void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::setChemicalSystemConcrete	(	Eigen::VectorXd const &	,
		double const	,
		double const
	)

privatepure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

Referenced by setChemicalSystem().

◆ setChemicalSystemID()

virtual void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::setChemicalSystemID ( std::size_t const )

pure virtual

Implemented in ProcessLib::ComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >.

Referenced by ProcessLib::ComponentTransport::ComponentTransportProcess::initializeConcreteProcess().

◆ setStaggeredCoupledSolutions()

void ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::setStaggeredCoupledSolutions	(	std::size_t const	,
		CoupledSolutionsForStaggeredScheme *const	coupling_term
	)

inline

Definition at line 69 of file ComponentTransportFEM.h.

     {
         _coupled_solutions = coupling_term;
     }

References _coupled_solutions.

Member Data Documentation

◆ _coupled_solutions

CoupledSolutionsForStaggeredScheme* ProcessLib::ComponentTransport::ComponentTransportLocalAssemblerInterface::_coupled_solutions {nullptr}

protected

Definition at line 209 of file ComponentTransportFEM.h.

Referenced by setStaggeredCoupledSolutions().

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

ProcessLib/ComponentTransport/ComponentTransportFEM.h

Detailed Description

Public Member Functions

Protected Attributes

Private Member Functions

Constructor & Destructor Documentation

◆ ComponentTransportLocalAssemblerInterface()

Member Function Documentation

◆ assembleReactionEquation()

◆ assembleReactionEquationConcrete()

◆ getInterpolatedLocalSolution()

◆ getIntPtDarcyVelocity()

◆ initializeChemicalSystem()

◆ initializeChemicalSystemConcrete()

◆ postSpeciationCalculation()

◆ setChemicalSystem()

◆ setChemicalSystemConcrete()

◆ setChemicalSystemID()

◆ setStaggeredCoupledSolutions()

Member Data Documentation

◆ _coupled_solutions