MixTransport::MixTransportElement::OpRhsBcOnValues Struct Reference

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

#include <users_modules/basic_finite_elements/mix_transport/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, DataForcesAndSourcesCore::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 1094 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 1103 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,
		DataForcesAndSourcesCore::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 1118 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++) {
        double x, y, z;
        if (getNormalsAtGaussPts().size1() == (unsigned int)nb_gauss_pts) {
          x = getHoCoordsAtGaussPts()(gg, 0);
          y = getHoCoordsAtGaussPts()(gg, 1);
          z = getHoCoordsAtGaussPts()(gg, 2);
        } else {
          x = getCoordsAtGaussPts()(gg, 0);
          y = getCoordsAtGaussPts()(gg, 1);
          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 1097 of file MixTransportElement.hpp.

F

Vec MixTransport::MixTransportElement::OpRhsBcOnValues::F

Definition at line 1098 of file MixTransportElement.hpp.

nF

VectorDouble MixTransport::MixTransportElement::OpRhsBcOnValues::nF

Definition at line 1109 of file MixTransportElement.hpp.

---

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

• users_modules/basic_finite_elements/mix_transport/src/MixTransportElement.hpp