FreeSurfaceOps::OpRhsG< I > Struct Template Reference

Rhs for G (chemical potential residual) More...

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

Inheritance diagram for FreeSurfaceOps::OpRhsG< I >:

Collaboration diagram for FreeSurfaceOps::OpRhsG< I >:

Public Member Functions
	OpRhsG (const std::string field_name, boost::shared_ptr< VectorDouble > h_ptr, boost::shared_ptr< MatrixDouble > grad_h_ptr, boost::shared_ptr< VectorDouble > g_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< VectorDouble >	hPtr
boost::shared_ptr< MatrixDouble >	gradHPtr
boost::shared_ptr< VectorDouble >	gPtr

Detailed Description

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

Rhs for G (chemical potential residual)

Parameters

field_name	Name of the field associated with the operator
h_ptr	Pointer to the free surface height vector (H)
grad_h_ptr	Pointer to the free surface height gradient matrix (∇H)
g_ptr	Pointer to the chemical potential vector (G)

Returns

MoFEMErrorCode

Implements Lovric 3.1d (excpet for wetting boundary term - see OpWettingAngleRhs/Lhs): R_g = ∫_Ω s(g - f(h)) - ∇s·(ε²∇h) dΩ

The template boolean I selects initialization (I=true) vs evolution (I=false).

Examples

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

Definition at line 1383 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpRhsG()

template<bool I>

FreeSurfaceOps::OpRhsG< I >::OpRhsG ( const std::string  field_name, boost::shared_ptr< VectorDouble >  h_ptr, boost::shared_ptr< MatrixDouble >  grad_h_ptr, boost::shared_ptr< VectorDouble >  g_ptr  )

inline

Definition at line 1385 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        hPtr(h_ptr), gradHPtr(grad_h_ptr), gPtr(g_ptr) {}

Member Function Documentation

iNtegrate()

template<bool I>

MoFEMErrorCode FreeSurfaceOps::OpRhsG< I >::iNtegrate ( EntitiesFieldData::EntData &  data )

inline

Definition at line 1391 of file FreeSurfaceOps.hpp.

                                                           {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_h = getFTensor0FromVec(*hPtr);
    auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
    auto t_g = getFTensor0FromVec(*gPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_base = data.getFTensor0N();
    auto t_diff_base = data.getFTensor1DiffN<SPACE_DIM>();
    auto t_w = getFTensor0IntegrationWeight();
 
    for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
      const double r = t_coords(0);
      const double alpha = t_w * vol * cylindrical(r);
 
      double set_h;
      if constexpr (I)
        set_h = init_h(t_coords(0), t_coords(1), t_coords(2));
      else
        set_h = t_h;
 
      const double f = get_f(set_h);
 
      int bb = 0;
      for (; bb != nbRows; ++bb) {
        locF[bb] += (t_base * alpha) * (t_g - f); // Bulk term: (g - f(h))
        locF[bb] -= (t_diff_base(i) * (eta2 * alpha)) * t_grad_h(i); // Diffusion term: -ε²∇h
        ++t_base;
        ++t_diff_base;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb) {
        ++t_base;
        ++t_diff_base;
      }
 
      ++t_h;
      ++t_grad_h;
      ++t_g;
 
      ++t_coords;
      ++t_w;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

gPtr

template<bool I>

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpRhsG< I >::gPtr

private

Definition at line 1444 of file FreeSurfaceOps.hpp.

gradHPtr

template<bool I>

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsG< I >::gradHPtr

private

Definition at line 1443 of file FreeSurfaceOps.hpp.

hPtr

template<bool I>

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

private

Definition at line 1442 of file FreeSurfaceOps.hpp.

---

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

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