Analytical Dirichlet boundary conditions. More...

#include <users_modules/basic_finite_elements/src/AnalyticalDirichlet.hpp>

Collaboration diagram for AnalyticalDirichletBC:

Classes
struct	ApproxField
	finite element to approximate analytical solution on surface More...

struct	DirichletBC
	Structure used to enforce analytical boundary conditions. More...

Public Member Functions
	AnalyticalDirichletBC (MoFEM::Interface &m_field)

template<typename FUNEVAL >
MoFEMErrorCode	setApproxOps (MoFEM::Interface &m_field, const std::string field_name, boost::shared_ptr< FUNEVAL > function_evaluator, const int field_number=0, const string nodals_positions="MESH_NODE_POSITIONS")
	Set operators used to calculate the rhs vector and the lhs matrix. More...

MoFEMErrorCode	setFiniteElement (MoFEM::Interface &m_field, string fe, string field, Range &tris, string nodals_positions="MESH_NODE_POSITIONS")
	set finite element More...

MoFEMErrorCode	setUpProblem (MoFEM::Interface &m_field, string problem)
	set problem solver and create matrices and vectors More...

MoFEMErrorCode	solveProblem (MoFEM::Interface &m_field, string problem, string fe, DirichletBC &bc, Range &tris)
	solve boundary problem More...

MoFEMErrorCode	solveProblem (MoFEM::Interface &m_field, string problem, string fe, DirichletBC &bc)
	solve boundary problem More...

MoFEMErrorCode	destroyProblem ()
	Destroy problem. More...

Public Attributes
ApproxField	approxField

Mat	A

Vec	D

Vec	F

KSP	kspSolver

Detailed Description

Analytical Dirichlet boundary conditions.

Definition at line 18 of file AnalyticalDirichlet.hpp.

Constructor & Destructor Documentation

◆ AnalyticalDirichletBC()

AnalyticalDirichletBC::AnalyticalDirichletBC ( MoFEM::Interface & m_field )

Definition at line 188 of file AnalyticalDirichlet.cpp.

189 : approxField(m_field){};

Member Function Documentation

◆ destroyProblem()

MoFEMErrorCode AnalyticalDirichletBC::destroyProblem ( )

Destroy problem.

Destroy matrices and vectors used to solve boundary problem, i.e. finding values of DOFs which approximate analytical solution.

Returns: error code

Definition at line 268 of file AnalyticalDirichlet.cpp.

                                                      {
   MoFEMFunctionBegin;
   CHKERR KSPDestroy(&kspSolver);
   CHKERR MatDestroy(&A);
   CHKERR VecDestroy(&F);
   CHKERR VecDestroy(&D);
   MoFEMFunctionReturn(0);
 }

◆ setApproxOps()

template<typename FUNEVAL >

MoFEMErrorCode AnalyticalDirichletBC::setApproxOps	(	MoFEM::Interface &	m_field,
		const std::string	field_name,
		boost::shared_ptr< FUNEVAL >	function_evaluator,
		const int	field_number = `0`,
		const string	nodals_positions = `"MESH_NODE_POSITIONS"`
	)

inline

Set operators used to calculate the rhs vector and the lhs matrix.

To enforce analytical function on boundary, first function has to be approximated by finite element base functions. This is done by solving system of linear equations, Following function set finite element operators to calculate the left hand side matrix and the left hand side matrix.

Parameters

m_field	interface
field_name	field name
function_evaluator	analytical function to evaluate
field_number	field index
nodals_positions	name of the field for ho-geometry description

Returns: error code

Definition at line 158 of file AnalyticalDirichlet.hpp.

                                                                       {
     MoFEMFunctionBeginHot;
     if (approxField.getLoopFeApprox().getOpPtrVector().empty()) {
       if (m_field.check_field(nodals_positions))
         CHKERR addHOOpsFace3D(nodals_positions, approxField.getLoopFeApprox(),
                               false, false);
       approxField.getLoopFeApprox().getOpPtrVector().push_back(
           new ApproxField::OpLhs(field_name));
     }
     approxField.getLoopFeApprox().getOpPtrVector().push_back(
         new ApproxField::OpRhs<FUNEVAL>(field_name, function_evaluator,
                                         field_number));
     MoFEMFunctionReturnHot(0);
   }

◆ setFiniteElement()

MoFEMErrorCode AnalyticalDirichletBC::setFiniteElement	(	MoFEM::Interface &	m_field,
		string	fe,
		string	field,
		Range &	tris,
		string	nodals_positions = `"MESH_NODE_POSITIONS"`
	)

set finite element

Parameters

m_field	mofem interface
fe	finite element name
field	field name
tris	faces where analytical boundary is given
nodals_positions	field having higher order geometry description

Returns: error code

Definition at line 192 of file AnalyticalDirichlet.cpp.

                                                                  {
   MoFEMFunctionBegin;
  
   CHKERR m_field.add_finite_element(fe, MF_ZERO);
   CHKERR m_field.modify_finite_element_add_field_row(fe, field);
   CHKERR m_field.modify_finite_element_add_field_col(fe, field);
   CHKERR m_field.modify_finite_element_add_field_data(fe, field);
   if (m_field.check_field(nodals_positions)) {
     CHKERR m_field.modify_finite_element_add_field_data(fe, nodals_positions);
   }
   CHKERR m_field.add_ents_to_finite_element_by_type(tris, MBTRI, fe);
   MoFEMFunctionReturn(0);
 }

◆ setUpProblem()

MoFEMErrorCode AnalyticalDirichletBC::setUpProblem	(	MoFEM::Interface &	m_field,
		string	problem
	)

set problem solver and create matrices and vectors

Parameters

m_field	mofem interface
problem	problem name

Returns: error code

Definition at line 208 of file AnalyticalDirichlet.cpp.

                                                                    {
   MoFEMFunctionBegin;
   CHKERR m_field.getInterface<VecManager>()->vecCreateGhost(problem, ROW, &F);
   CHKERR m_field.getInterface<VecManager>()->vecCreateGhost(problem, COL, &D);
   CHKERR m_field.getInterface<MatrixManager>()
       ->createMPIAIJWithArrays<PetscGlobalIdx_mi_tag>(problem, &A);
  
   CHKERR KSPCreate(PETSC_COMM_WORLD, &kspSolver);
   CHKERR KSPSetOperators(kspSolver, A, A);
   CHKERR KSPSetFromOptions(kspSolver);
   MoFEMFunctionReturn(0);
 }

◆ solveProblem() [1/2]

MoFEMErrorCode AnalyticalDirichletBC::solveProblem	(	MoFEM::Interface &	m_field,
		string	problem,
		string	fe,
		DirichletBC &	bc
	)

solve boundary problem

This functions solve for DOFs on boundary where analytical solution is given.

Parameters

m_field	mofem interface
problem	problem name
fe	finite element name
bc	Driblet boundary structure used to apply boundary conditions

Returns: [description]

Definition at line 257 of file AnalyticalDirichlet.cpp.

                                                                     {
   MoFEMFunctionBegin;
   EntityHandle fe_meshset = m_field.get_finite_element_meshset("BC_FE");
   Range bc_tris;
   CHKERR m_field.get_moab().get_entities_by_type(fe_meshset, MBTRI, bc_tris);
   return solveProblem(m_field, problem, fe, bc, bc_tris);
   MoFEMFunctionReturn(0);
 }

◆ solveProblem() [2/2]

MoFEMErrorCode AnalyticalDirichletBC::solveProblem	(	MoFEM::Interface &	m_field,
		string	problem,
		string	fe,
		DirichletBC &	bc,
		Range &	tris
	)

solve boundary problem

This functions solve for DOFs on boundary where analytical solution is given. i.e. finding values of DOFs which approximate analytical solution.

Parameters

m_field	mofem interface
problem	problem name
fe	finite element name
bc	Driblet boundary structure used to apply boundary conditions
tris	triangles on boundary

Returns: error code

Definition at line 222 of file AnalyticalDirichlet.cpp.

                                                                 {
   MoFEMFunctionBegin;
  
   CHKERR VecZeroEntries(F);
   CHKERR MatZeroEntries(A);
  
   approxField.getLoopFeApprox().snes_B = A;
   approxField.getLoopFeApprox().snes_f = F;
   CHKERR m_field.loop_finite_elements(problem, fe,
                                       approxField.getLoopFeApprox());
   CHKERR VecGhostUpdateBegin(F, ADD_VALUES, SCATTER_REVERSE);
   CHKERR VecGhostUpdateEnd(F, ADD_VALUES, SCATTER_REVERSE);
   CHKERR VecAssemblyBegin(F);
   CHKERR VecAssemblyEnd(F);
   CHKERR MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);
   CHKERR MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY);
  
   CHKERR KSPSolve(kspSolver, F, D);
   CHKERR VecGhostUpdateBegin(D, INSERT_VALUES, SCATTER_FORWARD);
   CHKERR VecGhostUpdateEnd(D, INSERT_VALUES, SCATTER_FORWARD);
  
   CHKERR m_field.getInterface<VecManager>()->setGlobalGhostVector(
       problem, ROW, D, INSERT_VALUES, SCATTER_REVERSE);
  
   bc.trisPtr = boost::shared_ptr<Range>(new Range(tris));
   bc.mapZeroRows.clear();
   bc.dofsIndices.clear();
   bc.dofsValues.clear();
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ A

Mat AnalyticalDirichletBC::A

Definition at line 219 of file AnalyticalDirichlet.hpp.

◆ approxField

ApproxField AnalyticalDirichletBC::approxField

Definition at line 137 of file AnalyticalDirichlet.hpp.

◆ D

Vec AnalyticalDirichletBC::D

Definition at line 220 of file AnalyticalDirichlet.hpp.

◆ F

Vec AnalyticalDirichletBC::F

Definition at line 220 of file AnalyticalDirichlet.hpp.

◆ kspSolver

KSP AnalyticalDirichletBC::kspSolver

Definition at line 221 of file AnalyticalDirichlet.hpp.

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

users_modules/basic_finite_elements/src/AnalyticalDirichlet.hpp
users_modules/basic_finite_elements/src/impl/AnalyticalDirichlet.cpp

Classes

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ AnalyticalDirichletBC()

Member Function Documentation

◆ destroyProblem()

◆ setApproxOps()

◆ setFiniteElement()

◆ setUpProblem()

◆ solveProblem() [1/2]

◆ solveProblem() [2/2]

Member Data Documentation

◆ A

◆ approxField

◆ D

◆ F

◆ kspSolver