FreeSurfaceOps::OpCalculateLift Struct Reference

Calculate lift force on free surface boundary. More...

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

Inheritance diagram for FreeSurfaceOps::OpCalculateLift:

Collaboration diagram for FreeSurfaceOps::OpCalculateLift:

Public Member Functions
	OpCalculateLift (const std::string field_name, boost::shared_ptr< VectorDouble > p_ptr, boost::shared_ptr< VectorDouble > lift_ptr, boost::shared_ptr< Range > ents_ptr)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< VectorDouble >	pPtr
boost::shared_ptr< VectorDouble >	liftPtr
boost::shared_ptr< Range >	entsPtr

Detailed Description

Calculate lift force on free surface boundary.

Parameters

field_name	Name of the field associated with the operator
p_ptr	Pointer to the pressure vector on the free surface
lift_ptr	Pointer to the lift force vector to be computed
ents_ptr	Pointer to the range of entities (edges) to consider

Returns

MoFEMErrorCode

Computes the integral ∫_Γ ​-p n dΓ, where:

• p : pressure on the free surface
• n : outward normal vector on the boundary

Examples

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

Definition at line 31 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpCalculateLift()

FreeSurfaceOps::OpCalculateLift::OpCalculateLift ( const std::string  field_name, boost::shared_ptr< VectorDouble >  p_ptr, boost::shared_ptr< VectorDouble >  lift_ptr, boost::shared_ptr< Range >  ents_ptr  )

inline

Definition at line 32 of file FreeSurfaceOps.hpp.

      : BoundaryEleOp(field_name, field_name, BoundaryEleOp::OPROW),
        pPtr(p_ptr), liftPtr(lift_ptr), entsPtr(ents_ptr) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBEDGE] = true;
  }

Member Function Documentation

doWork()

MoFEMErrorCode FreeSurfaceOps::OpCalculateLift::doWork ( int  row_side, EntityType  row_type, EntitiesFieldData::EntData &  data  )

inline

Definition at line 42 of file FreeSurfaceOps.hpp.

                                                        {
    MoFEMFunctionBegin;
 
    const auto fe_ent = getFEEntityHandle();
    if (entsPtr->find(fe_ent) != entsPtr->end()) {
 
      auto t_w = getFTensor0IntegrationWeight();
      auto t_p = getFTensor0FromVec(*pPtr);
      auto t_normal = getFTensor1Normal();
      auto t_coords = getFTensor1CoordsAtGaussPts();
      auto t_lift = getFTensor1FromArray<SPACE_DIM, SPACE_DIM>(*liftPtr);
 
      const auto nb_int_points = getGaussPts().size2();
 
      for (int gg = 0; gg != nb_int_points; gg++) {
 
        const double r = t_coords(0);
        const double alpha = cylindrical(r) * t_w;
        t_lift(i) -= t_normal(i) * (t_p * alpha); // −p n
 
        ++t_w;
        ++t_p;
        ++t_coords;
      }
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

entsPtr

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

private

Definition at line 75 of file FreeSurfaceOps.hpp.

liftPtr

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpCalculateLift::liftPtr

private

Definition at line 74 of file FreeSurfaceOps.hpp.

pPtr

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpCalculateLift::pPtr

private

Definition at line 73 of file FreeSurfaceOps.hpp.

---

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

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