FreeSurfaceOps::OpNormalConstrainLhs Struct Reference

[OpWettingAngleRhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpNormalConstrainLhs:

Collaboration diagram for FreeSurfaceOps::OpNormalConstrainLhs:

Public Member Functions
	OpNormalConstrainLhs (const std::string field_name_row, const std::string field_name_col)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Detailed Description

[OpWettingAngleRhs]

[OpNormalConstrainLhs]

Boundary normal constraint LHS (Jacobian block)

Parameters

field_name_row	Name of the field associated with the row operator
field_name_col	Name of the field associated with the column operator

Returns

MoFEMErrorCode

Computes the matrix entries for ∫_Γ ​w_L​⋅(n⋅u) dΓ, and its transpose ∫_Γ ​w_u​⋅(λn) dΓ, where:

• w_L : test function for Lagrange multiplier (λ)
• w_u : test function for velocity (u)
• n : outward normal vector on the boundary.

Examples

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

Definition at line 272 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpNormalConstrainLhs()

FreeSurfaceOps::OpNormalConstrainLhs::OpNormalConstrainLhs ( const std::string  field_name_row, const std::string  field_name_col  )

inline

Definition at line 274 of file FreeSurfaceOps.hpp.

      : AssemblyBoundaryEleOp(field_name_row, field_name_col,
                              AssemblyBoundaryEleOp::OPROWCOL) {
    assembleTranspose = true;
    sYmm = false;
  }

Member Function Documentation

iNtegrate()

MoFEMErrorCode FreeSurfaceOps::OpNormalConstrainLhs::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Definition at line 282 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_normal = getFTensor1Normal();
    auto t_row_base = row_data.getFTensor0N();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
      auto t_mat = getFTensor1FromPtr<U_FIELD_DIM>(&locMat(0, 0));
 
      const double r = t_coords(0);
      const double alpha = t_w * cylindrical(r);
 
      int rr = 0;
      for (; rr != nbRows; ++rr) {
 
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        const double a = alpha * t_row_base;
 
        for (int cc = 0; cc != nbCols / U_FIELD_DIM; ++cc) {
          t_mat(i) += (a * t_col_base) * t_normal(i); // α * N_row * N_col * n
          ++t_col_base;
          ++t_mat;
        }
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr)
        ++t_row_base;
 
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  };

---

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

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