ConvectiveMassElement::OpEnergy Struct Reference

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

Inheritance diagram for ConvectiveMassElement::OpEnergy:

Collaboration diagram for ConvectiveMassElement::OpEnergy:

Public Member Functions
	OpEnergy (const std::string field_name, BlockData &data, CommonData &common_data, SmartPetscObj< Vec > v)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Public Attributes
BlockData &	dAta
CommonData &	commonData
SmartPetscObj< Vec >	V
bool &	lInear
MatrixDouble3by3	h
MatrixDouble3by3	H
MatrixDouble3by3	invH
MatrixDouble3by3	F
VectorDouble3	v

Detailed Description

Definition at line 261 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpEnergy()

ConvectiveMassElement::OpEnergy::OpEnergy	(	const std::string	field_name,
		BlockData &	data,
		CommonData &	common_data,
		SmartPetscObj< Vec >	v
	)

Definition at line 727 of file ConvectiveMassElement.cpp.

    : VolumeElementForcesAndSourcesCore::UserDataOperator(
          field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
      dAta(data), commonData(common_data), V(v, true),
      lInear(commonData.lInear) {}

Member Function Documentation

doWork()

MoFEMErrorCode ConvectiveMassElement::OpEnergy::doWork	(	int	row_side,
		EntityType	row_type,
		EntitiesFieldData::EntData &	row_data
	)

Definition at line 737 of file ConvectiveMassElement.cpp.

                                                                            {
  MoFEMFunctionBeginHot;
 
  if (row_type != MBVERTEX) {
    MoFEMFunctionReturnHot(0);
  }
  if (dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) ==
      dAta.tEts.end()) {
    MoFEMFunctionReturnHot(0);
  }
 
  {
    double energy = 0;
    for (unsigned int gg = 0; gg < row_data.getN().size1(); gg++) {
      double val = getVolume() * getGaussPts()(3, gg);
      double rho0 = dAta.rho0;
      double rho;
      if (lInear) {
        rho = rho0;
      } else {
        h.resize(3, 3);
        noalias(h) =
            (commonData.gradAtGaussPts[commonData.spatialPositions][gg]);
        if (commonData.dataAtGaussPts["DOT_" + commonData.meshPositions]
                .size() > 0) {
          H.resize(3, 3);
          noalias(H) =
              (commonData.gradAtGaussPts[commonData.meshPositions][gg]);
          auto detH = determinantTensor3by3(H);
          invH.resize(3, 3);
          CHKERR invertTensor3by3(H, detH, invH);
          F.resize(3, 3);
          noalias(F) = prod(h, invH);
        } else {
          F.resize(3, 3);
          noalias(F) = h;
        }
        double detF = determinantTensor3by3(F);
        rho = detF * rho0;
      }
      v.resize(3);
      noalias(v) = commonData.dataAtGaussPts[commonData.spatialVelocities][gg];
      energy += 0.5 * (rho * val) * inner_prod(v, v);
    }
    CHKERR VecSetValue(V, 0, energy, ADD_VALUES);
  }
 
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

commonData

CommonData& ConvectiveMassElement::OpEnergy::commonData

Definition at line 265 of file ConvectiveMassElement.hpp.

dAta

BlockData& ConvectiveMassElement::OpEnergy::dAta

Definition at line 264 of file ConvectiveMassElement.hpp.

F

MatrixDouble3by3 ConvectiveMassElement::OpEnergy::F

Definition at line 272 of file ConvectiveMassElement.hpp.

h

MatrixDouble3by3 ConvectiveMassElement::OpEnergy::h

Definition at line 272 of file ConvectiveMassElement.hpp.

H

MatrixDouble3by3 ConvectiveMassElement::OpEnergy::H

Definition at line 272 of file ConvectiveMassElement.hpp.

invH

MatrixDouble3by3 ConvectiveMassElement::OpEnergy::invH

Definition at line 272 of file ConvectiveMassElement.hpp.

lInear

bool& ConvectiveMassElement::OpEnergy::lInear

Definition at line 267 of file ConvectiveMassElement.hpp.

V

SmartPetscObj<Vec> ConvectiveMassElement::OpEnergy::V

Definition at line 266 of file ConvectiveMassElement.hpp.

v

VectorDouble3 ConvectiveMassElement::OpEnergy::v

Definition at line 273 of file ConvectiveMassElement.hpp.

---

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

• users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp
• users_modules/basic_finite_elements/src/impl/ConvectiveMassElement.cpp