FreeSurfaceOps::OpNormalConstrainRhs Struct Reference

[OpNormalConstrainRhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpNormalConstrainRhs:

Collaboration diagram for FreeSurfaceOps::OpNormalConstrainRhs:

Public Member Functions
	OpNormalConstrainRhs (const std::string field_name, boost::shared_ptr< MatrixDouble > u_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< MatrixDouble >	uPtr

Detailed Description

[OpNormalConstrainRhs]

Boundary normal constraint RHS

Parameters

field_name	Name of the field associated with the operator
u_ptr	Pointer to the velocity matrix on the free surface

Returns

MoFEMErrorCode

Computes the integral ∫_Γ ​w_L​·(n·u) dΓ, where:

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

Examples

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

Definition at line 89 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpNormalConstrainRhs()

FreeSurfaceOps::OpNormalConstrainRhs::OpNormalConstrainRhs ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  u_ptr  )

inline

Definition at line 90 of file FreeSurfaceOps.hpp.

      : AssemblyBoundaryEleOp(field_name, field_name,
                              AssemblyBoundaryEleOp::OPROW),
        uPtr(u_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode FreeSurfaceOps::OpNormalConstrainRhs::iNtegrate ( EntitiesFieldData::EntData &  row_data )

inline

Definition at line 96 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_normal = getFTensor1Normal();
    auto t_u = getFTensor1FromMat<SPACE_DIM>(*uPtr);
    auto t_row_base = row_data.getFTensor0N();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double r = t_coords(0);
      const double alpha = t_w * cylindrical(r);
 
      int bb = 0;
      for (; bb != nbRows; ++bb) {
        locF[bb] += alpha * t_row_base * (t_normal(i) * t_u(i)); // λ * (n · u)
        ++t_row_base;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb)
        ++t_row_base;
 
      ++t_w;
      ++t_u;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

uPtr

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpNormalConstrainRhs::uPtr

private

Definition at line 128 of file FreeSurfaceOps.hpp.

---

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

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