FreeSurfaceOps::OpLhsG_dH< I > Struct Template Reference

Lhs for G dH. More...

#include "tutorials/vec-5_free_surface/src/FreeSurfaceOps.hpp"

Inheritance diagram for FreeSurfaceOps::OpLhsG_dH< I >:

Collaboration diagram for FreeSurfaceOps::OpLhsG_dH< I >:

Public Member Functions
	OpLhsG_dH (const std::string field_name_g, const std::string field_name_h, boost::shared_ptr< VectorDouble > h_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Attributes
boost::shared_ptr< VectorDouble >	hPtr

Detailed Description

template<bool I>
struct FreeSurfaceOps::OpLhsG_dH< I >

Lhs for G dH.

Parameters

field_name_g	Name of the field associated with the row operator (G)
field_name_h	Name of the field associated with the column operator (H)
h_ptr	Pointer to the free surface height vector (H)

Returns

MoFEMErrorCode

(Jacobian block ∂R_G/∂H) Linearization of the chemical potential equation (Lovric 3.1d) w.r.t. H

Examples

mofem/tutorials/vec-5_free_surface/free_surface.cpp.

Definition at line 1457 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpLhsG_dH()

template<bool I>

FreeSurfaceOps::OpLhsG_dH< I >::OpLhsG_dH ( const std::string  field_name_g, const std::string  field_name_h, boost::shared_ptr< VectorDouble >  h_ptr  )

inline

Definition at line 1459 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name_g, field_name_h,
                            AssemblyDomainEleOp::OPROWCOL),
        hPtr(h_ptr) {
    sYmm = false;
    assembleTranspose = false;
  }

Member Function Documentation

iNtegrate()

template<bool I>

MoFEMErrorCode FreeSurfaceOps::OpLhsG_dH< I >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Definition at line 1468 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_h = getFTensor0FromVec(*hPtr);
 
    auto t_row_base = row_data.getFTensor0N();
    auto t_row_diff_base = row_data.getFTensor1DiffN<SPACE_DIM>();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
      const double r = t_coords(0);
      const double alpha = t_w * vol * cylindrical(r);
 
      const double f_dh = get_f_dh(t_h) * alpha; // derivative of f(h)
      const double beta = eta2 * alpha; // coefficient for diffusion term -ε²
 
      int rr = 0;
      for (; rr != nbRows; ++rr) {
 
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        auto t_col_diff_base = col_data.getFTensor1DiffN<SPACE_DIM>(gg, 0);
 
        for (int cc = 0; cc != nbCols; ++cc) {
 
          if constexpr (I == false)
            locMat(rr, cc) -= (t_row_base * t_col_base) * f_dh; // bulk term derivative (g - f(h))
          locMat(rr, cc) -= (t_row_diff_base(i) * beta) * t_col_diff_base(i); // diffusion term derivative -ε²∇h
 
          ++t_col_base;
          ++t_col_diff_base;
        }
 
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      ++t_h;
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

hPtr

template<bool I>

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpLhsG_dH< I >::hPtr

private

Definition at line 1522 of file FreeSurfaceOps.hpp.

---

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

• tutorials/vec-5_free_surface/src/FreeSurfaceOps.hpp