FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp > Struct Template Reference

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

Inheritance diagram for FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >:

Collaboration diagram for FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >:

Public Member Functions
	OpCalculateQdotQLhs_dQ (const std::string row_field_name, const std::string col_field_name, ScalarFun resistance_function, boost::shared_ptr< MatrixDouble > mat_Grad_Ptr, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

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

Detailed Description

template<int DIM, typename AssemblyDomainEleOp>
struct FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 107 of file FiniteThermalOps.hpp.

Constructor & Destructor Documentation

OpCalculateQdotQLhs_dQ()

template<int DIM, typename AssemblyDomainEleOp >

FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateQdotQLhs_dQ ( const std::string row_field_name, const std::string col_field_name, ScalarFun resistance_function, boost::shared_ptr< MatrixDouble > mat_Grad_Ptr, boost::shared_ptr< Range > ents_ptr = nullptr )

inline

Definition at line 110 of file FiniteThermalOps.hpp.

      : AssemblyDomainEleOp(row_field_name, col_field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        resistanceFunction(resistance_function), matGradPtr(mat_Grad_Ptr) {}

Member Function Documentation

iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >

MoFEMErrorCode FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData & row_data, EntitiesFieldData::EntData & col_data )

protected

Definition at line 128 of file FiniteThermalOps.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'J', DIM> J;
  FTensor::Index<'M', DIM> M;
  FTensor::Index<'L', DIM> L;
 
  auto t_w = this->getFTensor0IntegrationWeight();
 
  auto t_row_base = row_data.getFTensor1N<3>();
 
  FTensor::Tensor2<double, DIM, DIM> t_F;
  FTensor::Tensor2<double, DIM, DIM> t_right_Cauchy_Green;
  FTensor::Tensor2<double, DIM, DIM> t_Kinv_over_J;
  FTensor::Tensor1<double, DIM> t_Kinv_over_J_row_base;
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
  auto t_grad = getFTensor2FromMat<DIM, DIM>(*(matGradPtr));
 
  // get coordinate at integration points
  auto t_coords = AssemblyDomainEleOp::getFTensor1CoordsAtGaussPts();
 
  for (size_t gg = 0; gg != AssemblyDomainEleOp::nbIntegrationPts; ++gg) {
 
    t_F(i, J) = t_grad(i, J) + t_kd(i, J); // Deformation gradient
    auto t_J = determinantTensor(t_F);     // Volume jacobian
 
    t_right_Cauchy_Green(M, L) =
        t_F(i, M) *
        t_F(i, L); // The inverse of the material thermal conductivity tensor
 
    const double alpha = this->getMeasure() * t_w;
 
    t_Kinv_over_J(M, L) =
        (alpha * resistanceFunction(t_coords(0), t_coords(1), t_coords(2)) /
         t_J) *
        t_right_Cauchy_Green(M, L);
 
    size_t rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
      t_Kinv_over_J_row_base(L) = t_Kinv_over_J(M, L) * t_row_base(M);
      auto t_col_base = col_data.getFTensor1N<3>(gg, 0);
      for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; ++cc) {
        this->locMat(rr, cc) += t_Kinv_over_J_row_base(L) * t_col_base(L);
        ++t_col_base;
      }
      ++t_row_base;
    }
    for (; rr < AssemblyDomainEleOp::nbRowBaseFunctions; ++rr)
      ++t_row_base;
 
    ++t_w; // move to another integration weight
    ++t_coords;
    ++t_grad;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

matGradPtr

template<int DIM, typename AssemblyDomainEleOp >

boost::shared_ptr<MatrixDouble> FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >::matGradPtr

protected

Definition at line 121 of file FiniteThermalOps.hpp.

resistanceFunction

template<int DIM, typename AssemblyDomainEleOp >

ScalarFun FiniteThermalOps::OpCalculateQdotQLhs_dQ< DIM, GAUSS, AssemblyDomainEleOp >::resistanceFunction

protected

Definition at line 120 of file FiniteThermalOps.hpp.

---

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

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