FreeSurfaceOps::OpNormalForceRhs Struct Reference

[OpNormalConstrainRhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpNormalForceRhs:

Collaboration diagram for FreeSurfaceOps::OpNormalForceRhs:

Public Member Functions
	OpNormalForceRhs (const std::string field_name, boost::shared_ptr< VectorDouble > lambda_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< VectorDouble >	lambdaPtr

Detailed Description

[OpNormalConstrainRhs]

Boundary normal force RHS

Parameters

field_name	Name of the field associated with the operator
lambda_ptr	Pointer to the Lagrange multiplier vector on the free surface

Returns

MoFEMErrorCode

Assembles the contribution of the Lagrange multiplier (λ) normal force ie ∫_Γ N_u · (λ n) dΓ into the momentum residual.

Examples

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

Definition at line 141 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpNormalForceRhs()

FreeSurfaceOps::OpNormalForceRhs::OpNormalForceRhs ( const std::string  field_name, boost::shared_ptr< VectorDouble >  lambda_ptr  )

inline

Definition at line 142 of file FreeSurfaceOps.hpp.

      : AssemblyBoundaryEleOp(field_name, field_name,
                              AssemblyDomainEleOp::OPROW),
        lambdaPtr(lambda_ptr) {}

Member Function Documentation

iNtegrate()

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

inline

Definition at line 148 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_normal = getFTensor1Normal();
    auto t_lambda = getFTensor0FromVec(*lambdaPtr);
    auto t_row_base = row_data.getFTensor0N();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      auto t_nf = getFTensor1FromArray<U_FIELD_DIM, U_FIELD_DIM>(locF);
 
      const double r = t_coords(0);
      const double alpha = t_w * cylindrical(r);
 
      int bb = 0;
      for (; bb != nbRows / U_FIELD_DIM; ++bb) {
 
        t_nf(i) += alpha * t_row_base * t_normal(i) * t_lambda; // N_u · (λ n)
        ++t_row_base;
        ++t_nf;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb)
        ++t_row_base;
 
      ++t_w;
      ++t_lambda;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

lambdaPtr

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpNormalForceRhs::lambdaPtr

private

Definition at line 184 of file FreeSurfaceOps.hpp.

---

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

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