FreeSurfaceOps::OpLhsU_dG Struct Reference

Lhs for U dG (Jacobian block ∂R_U/∂G) More...

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

Inheritance diagram for FreeSurfaceOps::OpLhsU_dG:

Collaboration diagram for FreeSurfaceOps::OpLhsU_dG:

Public Member Functions
	OpLhsU_dG (const std::string field_name_u, const std::string field_name_h, boost::shared_ptr< MatrixDouble > grad_h_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Attributes
boost::shared_ptr< MatrixDouble >	gradHPtr

Detailed Description

Lhs for U dG (Jacobian block ∂R_U/∂G)

Parameters

field_name_u	Name of the field associated with the row operator (U)
field_name_h	Name of the field associated with the column operator (H)
grad_h_ptr	Pointer to the free surface height gradient matrix (∇h)

Returns

MoFEMErrorCode

This operator assembles the contribution of the chemical potential G into the momentum Jacobian. It corresponds to linearization of the surface tension term −κ g ∇h wrt g.

Examples

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

Definition at line 884 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpLhsU_dG()

FreeSurfaceOps::OpLhsU_dG::OpLhsU_dG ( const std::string  field_name_u, const std::string  field_name_h, boost::shared_ptr< MatrixDouble >  grad_h_ptr  )

inline

Definition at line 886 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name_u, field_name_h,
                            AssemblyDomainEleOp::OPROWCOL),
        gradHPtr(grad_h_ptr) {
    sYmm = false;
    assembleTranspose = false;
  }

Member Function Documentation

iNtegrate()

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

inline

Definition at line 895 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_row_base = row_data.getFTensor0N();
    auto t_w = getFTensor0IntegrationWeight();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double r = t_coords(0);
      const double alpha = t_w * vol * cylindrical(r);
 
      FTensor::Tensor1<double, SPACE_DIM> t_phase_force_dg;
      t_phase_force_dg(i) = -alpha * kappa * t_grad_h(i); // -κ ∇h
 
      int rr = 0;
      for (; rr != nbRows / U_FIELD_DIM; ++rr) {
        auto t_mat =
            getFTensor1FromMat<U_FIELD_DIM, 1,
                               DataLayoutTraits<DataLayout::CoeffsByGauss>>(
                locMat, rr * U_FIELD_DIM);
        auto t_col_base = col_data.getFTensor0N(gg, 0);
 
        for (int cc = 0; cc != nbCols; ++cc) {
          const double bb = t_row_base * t_col_base;
          t_mat(i) += t_phase_force_dg(i) * bb;
 
          ++t_mat;
          ++t_col_base;
        }
 
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_base;
      }
 
      ++t_grad_h;
      ++t_coords;
      ++t_w;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

gradHPtr

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

private

Definition at line 946 of file FreeSurfaceOps.hpp.

---

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

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