FreeSurfaceOps::OpRhsH< I > Struct Template Reference

Rhs for H (phase-field residual) More...

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

Inheritance diagram for FreeSurfaceOps::OpRhsH< I >:

Collaboration diagram for FreeSurfaceOps::OpRhsH< I >:

Public Member Functions
	OpRhsH (const std::string field_name, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< VectorDouble > dot_h_ptr, boost::shared_ptr< VectorDouble > h_ptr, boost::shared_ptr< MatrixDouble > grad_h_ptr, boost::shared_ptr< MatrixDouble > grad_g_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	uPtr
boost::shared_ptr< VectorDouble >	dotHPtr
boost::shared_ptr< VectorDouble >	hPtr
boost::shared_ptr< MatrixDouble >	gradHPtr
boost::shared_ptr< MatrixDouble >	gradGPtr

Detailed Description

template<bool I>
struct FreeSurfaceOps::OpRhsH< I >

Rhs for H (phase-field residual)

Parameters

field_name	Name of the field associated with the operator
u_ptr	Pointer to the velocity matrix (U)
dot_h_ptr	Pointer to the free surface height time derivative vector (∂H/∂t)
h_ptr	Pointer to the free surface height vector (H)
grad_h_ptr	Pointer to the free surface height gradient matrix (∇H)
grad_g_ptr	Pointer to the chemical potential gradient matrix (∇G)

Returns

MoFEMErrorCode

Maps the conserved phase evolution (Lovric 3.1c): R_h = ∫_Ω v(∂h/∂t + u·∇h) + ∇v·(M(h) ∇g)) dΩ

The template boolean I selects initialization (I=true) vs evolution (I=false).

Examples

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

Definition at line 960 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpRhsH()

template<bool I>

FreeSurfaceOps::OpRhsH< I >::OpRhsH ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  u_ptr, boost::shared_ptr< VectorDouble >  dot_h_ptr, boost::shared_ptr< VectorDouble >  h_ptr, boost::shared_ptr< MatrixDouble >  grad_h_ptr, boost::shared_ptr< MatrixDouble >  grad_g_ptr  )

inline

Definition at line 962 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        uPtr(u_ptr), dotHPtr(dot_h_ptr), hPtr(h_ptr), gradHPtr(grad_h_ptr),
        gradGPtr(grad_g_ptr) {}

Member Function Documentation

iNtegrate()

template<bool I>

MoFEMErrorCode FreeSurfaceOps::OpRhsH< I >::iNtegrate ( EntitiesFieldData::EntData &  data )

inline

Definition at line 971 of file FreeSurfaceOps.hpp.

                                                           {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    auto t_base = data.getFTensor0N();
    auto t_diff_base = data.getFTensor1DiffN<SPACE_DIM>();
 
#ifndef NDEBUG
    if (data.getDiffN().size1() != data.getN().size1())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 1");
    if (data.getDiffN().size2() != data.getN().size2() * SPACE_DIM) {
      MOFEM_LOG("SELF", Sev::error)
          << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
      MOFEM_LOG("SELF", Sev::error) << data.getN();
      MOFEM_LOG("SELF", Sev::error) << data.getDiffN();
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong size 2");
    }
#endif
 
    if constexpr (I) {
      // Find h ∈ H¹(Ω) to minimise ∫(h - h_init)² dΩ + ∫∇v·(M ∇g) dΩ
 
      auto t_h = getFTensor0FromVec(*hPtr);
      auto t_grad_g = getFTensor1FromMat<SPACE_DIM>(*gradGPtr);
 
      for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
        const double r = t_coords(0);
        const double alpha = t_w * vol * cylindrical(r);
 
        const double set_h = init_h(t_coords(0), t_coords(1), t_coords(2));
        const double m = get_M(set_h) * alpha;
 
        int bb = 0;
        for (; bb != nbRows; ++bb) {
          locF[bb] += (t_base * alpha) * (t_h - set_h); // minimize ∫(h - h_init)² dΩ
          locF[bb] += (t_diff_base(i) * m) * t_grad_g(i);
          ++t_base;
          ++t_diff_base;
        }
 
        for (; bb < nbRowBaseFunctions; ++bb) {
          ++t_base;
          ++t_diff_base;
        }
 
        ++t_h;
        ++t_grad_g;
 
        ++t_coords;
        ++t_w;
      }
 
    } else {
 
      auto t_dot_h = getFTensor0FromVec(*dotHPtr);
      auto t_h = getFTensor0FromVec(*hPtr);
      auto t_u = getFTensor1FromMat<U_FIELD_DIM>(*uPtr);
      auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
      auto t_grad_g = getFTensor1FromMat<SPACE_DIM>(*gradGPtr);
 
      for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
        const double r = t_coords(0);
        const double alpha = t_w * vol * cylindrical(r);
 
        const double m = get_M(t_h) * alpha;
 
        int bb = 0;
        for (; bb != nbRows; ++bb) {
          locF[bb] += (t_base * alpha) * (t_dot_h); // time derivative ∂h/∂t
          locF[bb] += (t_base * alpha) * (t_grad_h(i) * t_u(i)); // convection u·∇h
          locF[bb] += (t_diff_base(i) * t_grad_g(i)) * m; // mobility/divergence (M ∇g)
          ++t_base;
          ++t_diff_base;
        }
 
        for (; bb < nbRowBaseFunctions; ++bb) {
          ++t_base;
          ++t_diff_base;
        }
 
        ++t_dot_h;
        ++t_h;
        ++t_grad_g;
        ++t_u;
        ++t_grad_h;
 
        ++t_coords;
        ++t_w;
      }
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dotHPtr

template<bool I>

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpRhsH< I >::dotHPtr

private

Definition at line 1071 of file FreeSurfaceOps.hpp.

gradGPtr

template<bool I>

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsH< I >::gradGPtr

private

Definition at line 1074 of file FreeSurfaceOps.hpp.

gradHPtr

template<bool I>

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsH< I >::gradHPtr

private

Definition at line 1073 of file FreeSurfaceOps.hpp.

hPtr

template<bool I>

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpRhsH< I >::hPtr

private

Definition at line 1072 of file FreeSurfaceOps.hpp.

uPtr

template<bool I>

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsH< I >::uPtr

private

Definition at line 1070 of file FreeSurfaceOps.hpp.

---

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

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