FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp > Struct Template Reference

#include "tutorials/adv-2/src/FiniteThermalOps.hpp"

Inheritance diagram for FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >:

Collaboration diagram for FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >:

Public Member Functions
	OpCalculateQdotQRhs (const std::string field_name, boost::shared_ptr< MatrixDouble > vec, ScalarFun resistance_function, boost::shared_ptr< MatrixDouble > mat_Grad_Ptr, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	fluxVec
ScalarFun	resistanceFunction
boost::shared_ptr< MatrixDouble >	matGradPtr

Detailed Description

template<int DIM, typename DomainEleOp>
struct FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >

Definition at line 27 of file FiniteThermalOps.hpp.

Constructor & Destructor Documentation

OpCalculateQdotQRhs()

template<int DIM, typename DomainEleOp >

FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >::OpCalculateQdotQRhs ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  vec, ScalarFun  resistance_function, boost::shared_ptr< MatrixDouble >  mat_Grad_Ptr, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 29 of file FiniteThermalOps.hpp.

      : DomainEleOp(field_name, field_name, DomainEleOp::OPROW), fluxVec(vec),
        resistanceFunction(resistance_function), matGradPtr(mat_Grad_Ptr) {}

Member Function Documentation

iNtegrate()

template<int DIM, typename DomainEleOp >

MoFEMErrorCode FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >::iNtegrate ( EntitiesFieldData::EntData &  data )

protected

Definition at line 45 of file FiniteThermalOps.hpp.

                                    {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'J', DIM> J;
  FTensor::Index<'m', DIM> m;
  FTensor::Index<'N', DIM> N;
 
  const auto nb_integration_pts = DomainEleOp::getGaussPts().size2();
 
  FTensor::Tensor2<double, DIM, DIM> t_F;
  FTensor::Tensor1<double, DIM> t_K_inv_Q_over_J;
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
  auto t_grad = getFTensor2FromMat<DIM, DIM>(*(matGradPtr));
 
  // get element volume
  const double vol = DomainEleOp::getMeasure();
  // get integration weights
  auto t_w = DomainEleOp::getFTensor0IntegrationWeight();
  // get base function gradient on rows
  auto t_row_base = data.getFTensor1N<3>();
  // get flux values
  auto t_flux = getFTensor1FromMat<DIM>(*fluxVec);
  // get coordinate at integration points
  auto t_coords = AssemblyDomainEleOp::getFTensor1CoordsAtGaussPts();
 
  // loop over integration points
  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
    // take into account Jacobian
    const double alpha = t_w * vol;
 
    t_F(i, J) = t_grad(i, J) + t_kd(i, J); // Deformation gradient
    auto t_J = determinantTensor(t_F);     // Volume jacobian
 
    t_K_inv_Q_over_J(J) =
        resistanceFunction(t_coords(0), t_coords(1), t_coords(2)) * t_F(m, J) *
        t_F(m, N) * t_flux(N) / t_J; // Value which multiplies delta Q
 
    // loop over rows base functions
    size_t rr = 0;
    for (; rr != DomainEleOp::nbRows; ++rr) {
      DomainEleOp::locF[rr] += t_K_inv_Q_over_J(J) * alpha * (t_row_base(J));
      ++t_row_base;
    }
    for (; rr < DomainEleOp::nbRowBaseFunctions; ++rr)
      ++t_row_base;
    ++t_w; // move to another integration weight
    ++t_coords;
    ++t_flux;
    ++t_grad;
  };
  MoFEMFunctionReturn(0);
}

Member Data Documentation

fluxVec

template<int DIM, typename DomainEleOp >

boost::shared_ptr<MatrixDouble> FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >::fluxVec

protected

Definition at line 38 of file FiniteThermalOps.hpp.

matGradPtr

template<int DIM, typename DomainEleOp >

boost::shared_ptr<MatrixDouble> FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >::matGradPtr

protected

Definition at line 40 of file FiniteThermalOps.hpp.

resistanceFunction

template<int DIM, typename DomainEleOp >

ScalarFun FiniteThermalOps::OpCalculateQdotQRhs< DIM, GAUSS, DomainEleOp >::resistanceFunction

protected

Definition at line 39 of file FiniteThermalOps.hpp.

---

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

• tutorials/adv-2/src/FiniteThermalOps.hpp