MixTransport::MixTransportElement::OpRhsBcOnValues Struct Reference

calculate \( \int_\mathcal{S} {\boldsymbol\tau} \cdot \mathbf{n}u \textrm{d}\mathcal{S} \) More...

#include <users_modules/tutorials/cor-0to1/src/MixTransportElement.hpp>

Inheritance diagram for MixTransport::MixTransportElement::OpRhsBcOnValues:

Collaboration diagram for MixTransport::MixTransportElement::OpRhsBcOnValues:

Public Member Functions
	OpRhsBcOnValues (MixTransportElement &ctx, const std::string fluxes_name, Vec f)
	Constructor. More...
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	Integrate boundary condition. More...

Public Attributes
MixTransportElement &	cTx
Vec	F
VectorDouble	nF

Detailed Description

calculate \( \int_\mathcal{S} {\boldsymbol\tau} \cdot \mathbf{n}u \textrm{d}\mathcal{S} \)

This terms comes from differentiation by parts. Note that in this Dirichlet boundary conditions are natural.

Definition at line 1055 of file MixTransportElement.hpp.

Constructor & Destructor Documentation

OpRhsBcOnValues()

MixTransport::MixTransportElement::OpRhsBcOnValues::OpRhsBcOnValues	(	MixTransportElement &	ctx,
		const std::string	fluxes_name,
		Vec	f
	)

Constructor.

Definition at line 1064 of file MixTransportElement.hpp.

        : MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator(
              fluxes_name, UserDataOperator::OPROW),
          cTx(ctx), F(f) {}

Member Function Documentation

doWork()

MoFEMErrorCode MixTransport::MixTransportElement::OpRhsBcOnValues::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

Integrate boundary condition.

Parameters

side	local index of entity
type	type of entity
data	data on entity

Returns

error code

Definition at line 1079 of file MixTransportElement.hpp.

                                                          {
      MoFEMFunctionBegin;
      if (data.getFieldData().size() == 0)
        MoFEMFunctionReturnHot(0);
      EntityHandle fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
      nF.resize(data.getIndices().size());
      nF.clear();
      int nb_gauss_pts = data.getN().size1();
      for (int gg = 0; gg < nb_gauss_pts; gg++) {
        const double x = getCoordsAtGaussPts()(gg, 0);
        const double y = getCoordsAtGaussPts()(gg, 1);
        const double z = getCoordsAtGaussPts()(gg, 2);
        
        double value;
        CHKERR cTx.getBcOnValues(fe_ent, gg, x, y, z, value);
        ;
        double w = getGaussPts()(2, gg) * 0.5;
        if (getNormalsAtGaussPts().size1() == (unsigned int)nb_gauss_pts) {
          noalias(nF) +=
              w * prod(data.getVectorN<3>(gg), getNormalsAtGaussPts(gg)) *
              value;
        } else {
          noalias(nF) += w * prod(data.getVectorN<3>(gg), getNormal()) * value;
        }
      }
      Vec f = (F != PETSC_NULL) ? F : getFEMethod()->ts_F;
      CHKERR VecSetValues(f, data.getIndices().size(), &data.getIndices()[0],
                          &nF[0], ADD_VALUES);
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

cTx

MixTransportElement& MixTransport::MixTransportElement::OpRhsBcOnValues::cTx

Definition at line 1058 of file MixTransportElement.hpp.

F

Vec MixTransport::MixTransportElement::OpRhsBcOnValues::F

Definition at line 1059 of file MixTransportElement.hpp.

nF

VectorDouble MixTransport::MixTransportElement::OpRhsBcOnValues::nF

Definition at line 1070 of file MixTransportElement.hpp.

---

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

• users_modules/tutorials/cor-0to1/src/MixTransportElement.hpp