FreeSurfaceOps::OpWettingAngleLhs Struct Reference

[OpNormalConstrainLhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpWettingAngleLhs:

Collaboration diagram for FreeSurfaceOps::OpWettingAngleLhs:

Public Member Functions
	OpWettingAngleLhs (const std::string row_field_name, boost::shared_ptr< MatrixDouble > grad_h_ptr, boost::shared_ptr< std::vector< VectorInt > > col_ind_ptr, boost::shared_ptr< std::vector< MatrixDouble > > col_diff_base_ptr, boost::shared_ptr< Range > ents_ptr=nullptr, double wetting_angle=0)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Private Attributes
MatrixDouble	locMat
boost::shared_ptr< MatrixDouble >	gradHPtr
boost::shared_ptr< Range >	entsPtr
double	wettingAngle
boost::shared_ptr< std::vector< VectorInt > >	colIndicesPtr
boost::shared_ptr< std::vector< MatrixDouble > >	colDiffBaseFunctionsPtr

Detailed Description

[OpNormalConstrainLhs]

[OpWettingAngleLhs]

Boundary wetting-angle LHS (linearization)

Parameters

field_name	Name of the field associated with the operator
grad_h_ptr	Pointer to the gradient of the free surface height matrix
col_ind_ptr	Pointer to the column indices vector
col_diff_base_ptr	Pointer to the column difference base functions vector
ents_ptr	Pointer to the range of entities (edges) to consider
wetting_angle	Contact wetting angle in degrees

Returns

MoFEMErrorCode

Linearization of the wetting-angle boundary term (3.1d - Gamma part) used to assemble the corresponding Jacobian contributions for G (and coupling into H).

δ R_g / δ h = ∫_Γ s ε² (∇h / ||∇h||) · ∇(δ h) cos(α) dΓ

Examples

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

Definition at line 340 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpWettingAngleLhs()

FreeSurfaceOps::OpWettingAngleLhs::OpWettingAngleLhs ( const std::string  row_field_name, boost::shared_ptr< MatrixDouble >  grad_h_ptr, boost::shared_ptr< std::vector< VectorInt > >  col_ind_ptr, boost::shared_ptr< std::vector< MatrixDouble > >  col_diff_base_ptr, boost::shared_ptr< Range >  ents_ptr = nullptr, double  wetting_angle = 0  )

inline

Definition at line 342 of file FreeSurfaceOps.hpp.

      : BoundaryEleOp(row_field_name, BoundaryEleOp::OPROW),
        gradHPtr(grad_h_ptr), colIndicesPtr(col_ind_ptr),
        colDiffBaseFunctionsPtr(col_diff_base_ptr), entsPtr(ents_ptr),
        wettingAngle(wetting_angle) {}

Member Function Documentation

doWork()

MoFEMErrorCode FreeSurfaceOps::OpWettingAngleLhs::doWork ( int  side, EntityType  type, DataForcesAndSourcesCore::EntData &  data  )

inline

Definition at line 353 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
    if (entsPtr) {
      if (entsPtr->find(BoundaryEleOp::getFEEntityHandle()) == entsPtr->end())
        MoFEMFunctionReturnHot(0);
    }
    const double area = getMeasure();
 
    const auto row_size = data.getIndices().size();
    if (row_size == 0)
      MoFEMFunctionReturnHot(0);
 
    auto integrate = [&](auto col_indicies, auto &col_diff_base_functions) {
      MoFEMFunctionBegin;
 
      const auto col_size = col_indicies.size();
 
      locMat.resize(row_size, col_size, false);
      locMat.clear();
      int nb_gp = getGaussPts().size2();
      int nb_rows = data.getIndices().size();
 
      auto t_w = getFTensor0IntegrationWeight();
      auto t_coords = getFTensor1CoordsAtGaussPts();
      auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
      auto t_row_base = data.getFTensor0N();
      int nb_row_base_functions = data.getN().size2();
 
      auto s = wetting_angle_sub_stepping(getFEMethod()->ts_step);
 
      for (int gg = 0; gg != nb_gp; ++gg) {
 
        const double r = t_coords(0);
        const double alpha = t_w * area * cylindrical(r);
        const double h_grad_norm = sqrt(t_grad_h(i) * t_grad_h(i) + eps);
        const double one_over_h_grad_norm = 1. / h_grad_norm;
        const double beta = s * alpha * eta2 * one_over_h_grad_norm *
                            std::cos(M_PI * wettingAngle / 180); // s * α * η² * (1/||∇h||) * cos(angle)
 
        int rr = 0;
        for (; rr != nb_rows; ++rr) {
          const double delta = beta * t_row_base; // β * N_row
 
          auto ptr = &col_diff_base_functions(gg, 0);
          auto t_col_diff_base = getFTensor1FromPtr<SPACE_DIM>(ptr);
 
          for (int cc = 0; cc != col_size; ++cc) {
            locMat(rr, cc) += t_col_diff_base(i) * (delta * t_grad_h(i));
            ++t_col_diff_base;
          }
          ++t_row_base;
        }
 
        for (; rr < nb_row_base_functions; ++rr) {
          ++t_row_base;
        }
 
        ++t_grad_h;
        ++t_w;
        ++t_coords;
      }
 
      MoFEMFunctionReturn(0);
    };
 
    for (auto c = 0; c != colIndicesPtr->size(); ++c) {
 
      auto &col_ind = (*colIndicesPtr)[c];
      if (col_ind.size()) {
        auto &diff_base = (*colDiffBaseFunctionsPtr)[c];
 
        CHKERR integrate(col_ind, diff_base);
 
        CHKERR MatSetValues(getKSPB(), data.getIndices().size(),
                            &*data.getIndices().begin(), col_ind.size(),
                            &*col_ind.begin(), &locMat(0, 0), ADD_VALUES);
      }
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

colDiffBaseFunctionsPtr

boost::shared_ptr<std::vector<MatrixDouble> > FreeSurfaceOps::OpWettingAngleLhs::colDiffBaseFunctionsPtr

private

Definition at line 443 of file FreeSurfaceOps.hpp.

colIndicesPtr

boost::shared_ptr<std::vector<VectorInt> > FreeSurfaceOps::OpWettingAngleLhs::colIndicesPtr

private

Definition at line 442 of file FreeSurfaceOps.hpp.

entsPtr

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

private

Definition at line 440 of file FreeSurfaceOps.hpp.

gradHPtr

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

private

Definition at line 439 of file FreeSurfaceOps.hpp.

locMat

MatrixDouble FreeSurfaceOps::OpWettingAngleLhs::locMat

private

Definition at line 437 of file FreeSurfaceOps.hpp.

wettingAngle

double FreeSurfaceOps::OpWettingAngleLhs::wettingAngle

private

Definition at line 441 of file FreeSurfaceOps.hpp.

---

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

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