FreeSurfaceOps::OpWettingAngleRhs Struct Reference

[OpWettingAngleRhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpWettingAngleRhs:

Collaboration diagram for FreeSurfaceOps::OpWettingAngleRhs:

Public Member Functions
	OpWettingAngleRhs (const std::string field_name, boost::shared_ptr< MatrixDouble > grad_h_ptr, boost::shared_ptr< Range > ents_ptr=nullptr, double wetting_angle=0)
MoFEMErrorCode	iNtegrate (DataForcesAndSourcesCore::EntData &row_data)

Private Attributes
boost::shared_ptr< MatrixDouble >	gradHPtr
boost::shared_ptr< Range >	entsPtr
double	wettingAngle

Detailed Description

[OpWettingAngleRhs]

Boundary wetting-angle RHS

Parameters

field_name	Name of the field associated with the operator
grad_h_ptr	Pointer to the gradient of the free surface height matrix
ents_ptr	Pointer to the range of entities (edges) to consider
wetting_angle	Contact wetting angle in degrees

Returns

MoFEMErrorCode

Implements the boundary term from eq. 3.2d:

• ∫_Γ sh ε² ||∇ϕ|| cos(α) dΓ The operator computes rhs_wetting = s * η² * ||∇h|| * cos(angle) and assembles it into the G residual on selected meshset edges.

Examples

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

Definition at line 201 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpWettingAngleRhs()

FreeSurfaceOps::OpWettingAngleRhs::OpWettingAngleRhs ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  grad_h_ptr, boost::shared_ptr< Range >  ents_ptr = nullptr, double  wetting_angle = 0  )

inline

Definition at line 202 of file FreeSurfaceOps.hpp.

      : AssemblyBoundaryEleOp(field_name, field_name,
                              AssemblyBoundaryEleOp::OPROW),
        gradHPtr(grad_h_ptr), entsPtr(ents_ptr), wettingAngle(wetting_angle) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode FreeSurfaceOps::OpWettingAngleRhs::iNtegrate ( DataForcesAndSourcesCore::EntData &  row_data )

inline

Definition at line 210 of file FreeSurfaceOps.hpp.

                                                                      {
    MoFEMFunctionBegin;
    if (entsPtr) {
      if (entsPtr->find(AssemblyBoundaryEleOp::getFEEntityHandle()) ==
          entsPtr->end())
        MoFEMFunctionReturnHot(0);
    }
    const double area = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_row_base = row_data.getFTensor0N();
    auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto s = wetting_angle_sub_stepping(getFEMethod()->ts_step);
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double r = t_coords(0);
      const double alpha = t_w * cylindrical(r) * area;
      const double h_grad_norm = sqrt(t_grad_h(i) * t_grad_h(i) + eps); // ||∇h||
      const double cos_angle = std::cos(M_PI * wettingAngle / 180);
      const double rhs_wetting = s * eta2 * h_grad_norm * cos_angle; // s * η² * ||∇h|| * cos(angle)
 
      // cerr << "pass "
      //      << h_grad_norm <<"\n";
      int bb = 0;
      for (; bb != nbRows; ++bb) {
        locF[bb] += alpha * t_row_base * rhs_wetting;
        ++t_row_base;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb)
        ++t_row_base;
 
      ++t_w;
      ++t_grad_h;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

entsPtr

boost::shared_ptr<Range> FreeSurfaceOps::OpWettingAngleRhs::entsPtr

private

Definition at line 254 of file FreeSurfaceOps.hpp.

gradHPtr

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpWettingAngleRhs::gradHPtr

private

Definition at line 253 of file FreeSurfaceOps.hpp.

wettingAngle

double FreeSurfaceOps::OpWettingAngleRhs::wettingAngle

private

Definition at line 255 of file FreeSurfaceOps.hpp.

---

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

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