FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true > Struct Template Reference

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

Inheritance diagram for FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true >:

Collaboration diagram for FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true >:

Public Member Functions
	OpCalculateQdotQLhs_dU (const std::string row_field_name, const std::string col_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 &row_data, EntitiesFieldData::EntData &col_data)

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

Detailed Description

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

Definition at line 197 of file FiniteThermalOps.hpp.

Constructor & Destructor Documentation

OpCalculateQdotQLhs_dU()

template<int DIM, typename AssemblyDomainEleOp >

FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true >::OpCalculateQdotQLhs_dU ( const std::string row_field_name, const std::string col_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 200 of file FiniteThermalOps.hpp.

      : AssemblyDomainEleOp(row_field_name, col_field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        fluxVec(vec), resistanceFunction(resistance_function),
        matGradPtr(mat_Grad_Ptr) {
    AssemblyDomainEleOp::sYmm = false;
  }

Member Function Documentation

iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 223 of file FiniteThermalOps.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'J', DIM> J;
  FTensor::Index<'m', DIM> m;
  FTensor::Index<'M', DIM> M;
  FTensor::Index<'l', DIM> l;
  FTensor::Index<'L', DIM> L;
 
  auto t_w = this->getFTensor0IntegrationWeight();
 
  // Always need three components for vectorial basis
  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_F_inv;
  FTensor::Tensor1<double, DIM> t_F_flux;
  FTensor::Tensor1<double, DIM> t_F_row_base;
  FTensor::Tensor1<double, DIM> t_FtF_Q;
  FTensor::Tensor2<double, DIM, DIM> t_FtF_Finv_Q_row_base;
  FTensor::Tensor2<double, DIM, DIM> t_F_Q_row_base;
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
  auto t_grad = getFTensor2FromMat<DIM, DIM>(*(matGradPtr));
 
  // get flux values
  auto t_flux = getFTensor1FromMat<DIM>(*fluxVec);
 
  // for each integration point
  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
    CHKERR invertTensor(t_F, t_J, t_F_inv);
 
    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;
 
    // get the vector of the local matrix for the derivative wrt a vectorial
    // field with a scalar basis
    auto t_vec = getFTensor1FromPtr<DIM>(&this->locMat(0, 0));
 
    // get coordinate at integration points
    auto t_coords = AssemblyDomainEleOp::getFTensor1CoordsAtGaussPts();
 
    auto alpha_resistance_Jinv =
        alpha * resistanceFunction(t_coords(0), t_coords(1), t_coords(2)) / t_J;
    t_F_flux(l) = t_F(l, L) * t_flux(L);
 
    t_FtF_Q(M) = t_right_Cauchy_Green(M, L) * t_flux(L);
 
    // for each row coefficient
    size_t rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
      auto t_col_grad_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
 
      t_FtF_Finv_Q_row_base(J, l) = t_FtF_Q(M) * t_row_base(M) * t_F_inv(J, l);
      t_F_row_base(l) = t_F(l, M) * t_row_base(M);
      t_F_Q_row_base(l, M) = t_F_flux(l) * t_row_base(M);
 
      // for each coupled coefficient in the columns
      // divied by DIM since scalar basis for vectorial field
      for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; ++cc) {
        t_vec(l) += alpha_resistance_Jinv *
                    (-t_FtF_Finv_Q_row_base(J, l) * t_col_grad_base(J) +
                     t_F_Q_row_base(l, M) * t_col_grad_base(M) +
                     t_F_row_base(l) * t_col_grad_base(L) * t_flux(L));
        ++t_col_grad_base;
        ++t_vec;
      }
      ++t_row_base;
    }
    for (; rr < AssemblyDomainEleOp::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 AssemblyDomainEleOp >

boost::shared_ptr<MatrixDouble> FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true >::fluxVec

protected

Definition at line 214 of file FiniteThermalOps.hpp.

matGradPtr

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 216 of file FiniteThermalOps.hpp.

resistanceFunction

template<int DIM, typename AssemblyDomainEleOp >

ScalarFun FiniteThermalOps::OpCalculateQdotQLhs_dU< DIM, GAUSS, AssemblyDomainEleOp, true >::resistanceFunction

protected

Definition at line 215 of file FiniteThermalOps.hpp.

---

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

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