MixTransport::MixTransportElement::OpSkeleton Struct Reference

calculate jump on entities More...

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

Inheritance diagram for MixTransport::MixTransportElement::OpSkeleton:

Collaboration diagram for MixTransport::MixTransportElement::OpSkeleton:

Classes
struct	OpVolSide
	calculate values on adjacent tetrahedra to face More...

Public Member Functions
	OpSkeleton (MixTransportElement &ctx, const std::string flux_name)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
VolumeElementForcesAndSourcesCoreOnSide	volSideFe
	volume element to get values from adjacent tets to face More...
map< int, VectorDouble >	valMap
MixTransportElement &	cTx

Detailed Description

calculate jump on entities

Definition at line 1454 of file MixTransportElement.hpp.

Constructor & Destructor Documentation

OpSkeleton()

MixTransport::MixTransportElement::OpSkeleton::OpSkeleton	(	MixTransportElement &	ctx,
		const std::string	flux_name
	)

Definition at line 1494 of file MixTransportElement.hpp.

        : MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator(
              flux_name, UserDataOperator::OPROW),
          volSideFe(ctx.mField), cTx(ctx) {
      volSideFe.getOpPtrVector().push_back(new OpSkeleton::OpVolSide(valMap));
    }

Member Function Documentation

doWork()

MoFEMErrorCode MixTransport::MixTransportElement::OpSkeleton::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

Definition at line 1501 of file MixTransportElement.hpp.

                                                                 {
      MoFEMFunctionBegin;
 
      if (type == MBTRI) {
        EntityHandle fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
 
        double def_val = 0;
        Tag th_error_jump;
        CHKERR cTx.mField.get_moab().tag_get_handle(
            "ERROR_JUMP", 1, MB_TYPE_DOUBLE, th_error_jump,
            MB_TAG_CREAT | MB_TAG_SPARSE, &def_val);
        double *error_jump_ptr;
        CHKERR cTx.mField.get_moab().tag_get_by_ptr(
            th_error_jump, &fe_ent, 1, (const void **)&error_jump_ptr);
        *error_jump_ptr = 0;
 
        // check if this is essential boundary condition
        EntityHandle essential_bc_meshset =
            cTx.mField.get_finite_element_meshset("MIX_BCFLUX");
        if (cTx.mField.get_moab().contains_entities(essential_bc_meshset,
                                                    &fe_ent, 1)) {
          // essential bc, np jump then, exit and go to next face
          MoFEMFunctionReturnHot(0);
        }
 
        // calculate values form adjacent tets
        valMap.clear();
        CHKERR loopSideVolumes("MIX", volSideFe);
        ;
 
        int nb_gauss_pts = data.getN().size1();
 
        // it is only one face, so it has to be bc natural boundary condition
        if (valMap.size() == 1) {
          if (static_cast<int>(valMap.begin()->second.size()) != nb_gauss_pts) {
            SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
                    "wrong number of integration points");
          }
          for (int gg = 0; gg != nb_gauss_pts; gg++) {
            double x, y, z;
            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);
            if (static_cast<int>(getNormalsAtGaussPts().size1()) ==
                nb_gauss_pts) {
              w *= norm_2(getNormalsAtGaussPts(gg)) * 0.5;
            } else {
              w *= getArea();
            }
            *error_jump_ptr += w * pow(value - valMap.begin()->second[gg], 2);
          }
        } else if (valMap.size() == 2) {
          for (int gg = 0; gg != nb_gauss_pts; gg++) {
            double w = getGaussPts()(2, gg);
            if (getNormalsAtGaussPts().size1() == (unsigned int)nb_gauss_pts) {
              w *= norm_2(getNormalsAtGaussPts(gg)) * 0.5;
            } else {
              w *= getArea();
            }
            double delta = valMap.at(1)[gg] - valMap.at(-1)[gg];
            *error_jump_ptr += w * pow(delta, 2);
          }
        } else {
          SETERRQ1(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
                   "data inconsistency, wrong number of neighbors "
                   "valMap.size() = %d",
                   valMap.size());
        }
      }
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

cTx

MixTransportElement& MixTransport::MixTransportElement::OpSkeleton::cTx

Definition at line 1492 of file MixTransportElement.hpp.

valMap

map<int, VectorDouble> MixTransport::MixTransportElement::OpSkeleton::valMap

store values at integration point, key of the map is sense of face in respect to adjacent tetrahedra

Definition at line 1465 of file MixTransportElement.hpp.

volSideFe

VolumeElementForcesAndSourcesCoreOnSide MixTransport::MixTransportElement::OpSkeleton::volSideFe

volume element to get values from adjacent tets to face

Definition at line 1460 of file MixTransportElement.hpp.

---

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

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