FreeSurfaceOps::OpWettingAngleLhs Struct Reference

#include <tutorials/vec-5/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

Examples

free_surface.cpp.

Definition at line 246 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 248 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 259 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);
 
        int rr = 0;
        for (; rr != nb_rows; ++rr) {
          const double delta = beta * t_row_base;
 
          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 349 of file FreeSurfaceOps.hpp.

colIndicesPtr

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

private

Definition at line 348 of file FreeSurfaceOps.hpp.

entsPtr

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

private

Definition at line 346 of file FreeSurfaceOps.hpp.

gradHPtr

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

private

Definition at line 345 of file FreeSurfaceOps.hpp.

locMat

MatrixDouble FreeSurfaceOps::OpWettingAngleLhs::locMat

private

Definition at line 343 of file FreeSurfaceOps.hpp.

wettingAngle

double FreeSurfaceOps::OpWettingAngleLhs::wettingAngle

private

Definition at line 347 of file FreeSurfaceOps.hpp.

---

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

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