FreeSurfaceOps::OpRhsU Struct Reference

[OpWettingAngleLhs] More...

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

Inheritance diagram for FreeSurfaceOps::OpRhsU:

Collaboration diagram for FreeSurfaceOps::OpRhsU:

Public Member Functions
	OpRhsU (const std::string field_name, boost::shared_ptr< MatrixDouble > dot_u_ptr, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > grad_u_ptr, boost::shared_ptr< VectorDouble > h_ptr, boost::shared_ptr< MatrixDouble > grad_h_ptr, boost::shared_ptr< VectorDouble > g_ptr, boost::shared_ptr< VectorDouble > p_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	dotUPtr
boost::shared_ptr< MatrixDouble >	uPtr
boost::shared_ptr< MatrixDouble >	gradUPtr
boost::shared_ptr< VectorDouble >	hPtr
boost::shared_ptr< MatrixDouble >	gradHPtr
boost::shared_ptr< VectorDouble >	gPtr
boost::shared_ptr< VectorDouble >	pPtr

Detailed Description

[OpWettingAngleLhs]

[OpRhsU]

Rhs for U (momentum residual)

Parameters

field_name	Name of the field associated with the operator
dot_u_ptr	Pointer to the velocity time derivative matrix (∂U/∂t)
u_ptr	Pointer to the velocity matrix (U)
grad_u_ptr	Pointer to the velocity gradient matrix (∇U)
h_ptr	Pointer to the free surface height vector (h)
grad_h_ptr	Pointer to the free surface height gradient matrix (∇h)
g_ptr	Pointer to the chemical potential vector (g)
p_ptr	Pointer to the pressure vector (p)

Returns

MoFEMErrorCode

This operator implements the volume terms of the momentum equation (not Jacobian) (Lovric 3.1a)

R_u = ∫_Ω w·(ρ(h)(∂u/∂t + u·∇u + a_0) - κg∇h) - (∇·w)p + (∇w):(2μ(h)D(u))dΩ

Code ↔ PDE term mapping (paper notation : code variable)

• ρ(ϕ) ∂u/∂t : t_inertia_force(i) = (rho * alpha) * t_dot_u(i)
• ρ(ϕ) (u·∇)u : t_convection(i) = (rho * alpha) * (t_u(j) * t_grad_u(i,j))
• ∇·(2 μ(ϕ) D(u)) : t_stress built from get_D(2*mu) and t_grad_u, then assembled with t_diff_base into locF
• −∇p (pressure grad) : incorporated via t_kd(i,j)*t_p term inside t_stress
• −κ η ∇ϕ (surface) : t_phase_force(i) = -alpha * kappa * t_g * t_grad_h(i)
• body force ρ a0 : t_gravity(SPACE_DIM-1) = alpha * rho * a0
• buoyancy force : (currently commented out) t_buoyancy(SPACE_DIM-1) = -(alpha * rho * a0) * t_h;

Examples

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

Definition at line 475 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpRhsU()

FreeSurfaceOps::OpRhsU::OpRhsU ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  dot_u_ptr, boost::shared_ptr< MatrixDouble >  u_ptr, boost::shared_ptr< MatrixDouble >  grad_u_ptr, boost::shared_ptr< VectorDouble >  h_ptr, boost::shared_ptr< MatrixDouble >  grad_h_ptr, boost::shared_ptr< VectorDouble >  g_ptr, boost::shared_ptr< VectorDouble >  p_ptr  )

inline

Definition at line 477 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        dotUPtr(dot_u_ptr), uPtr(u_ptr), gradUPtr(grad_u_ptr), hPtr(h_ptr),
        gradHPtr(grad_h_ptr), gPtr(g_ptr), pPtr(p_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode FreeSurfaceOps::OpRhsU::iNtegrate ( EntitiesFieldData::EntData &  data )

inline

Definition at line 489 of file FreeSurfaceOps.hpp.

                                                           {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_dot_u = getFTensor1FromMat<U_FIELD_DIM>(*dotUPtr);
    auto t_u = getFTensor1FromMat<U_FIELD_DIM>(*uPtr);
    auto t_p = getFTensor0FromVec(*pPtr);
    auto t_grad_u = getFTensor2FromMat<U_FIELD_DIM, SPACE_DIM>(*gradUPtr);
    auto t_h = getFTensor0FromVec(*hPtr);
    auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
    auto t_g = getFTensor0FromVec(*gPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_base = data.getFTensor0N();
    auto t_diff_base = data.getFTensor1DiffN<SPACE_DIM>();
 
    auto t_w = getFTensor0IntegrationWeight();
 
    constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<double>();
    FTensor::Tensor2_symmetric<double, SPACE_DIM> t_stress;
    FTensor::Tensor1<double, U_FIELD_DIM> t_phase_force;
    FTensor::Tensor1<double, U_FIELD_DIM> t_inertia_force;
    FTensor::Tensor1<double, U_FIELD_DIM> t_convection;
    FTensor::Tensor1<double, U_FIELD_DIM> t_buoyancy;
    FTensor::Tensor1<double, U_FIELD_DIM> t_gravity;
    FTensor::Tensor1<double, U_FIELD_DIM> t_forces;
 
    t_buoyancy(i) = 0;
    t_gravity(i) = 0;
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const double r = t_coords(0);
      const double alpha = t_w * vol * cylindrical(r);
 
      const double rho = phase_function(t_h, rho_diff, rho_ave);
      const double mu = phase_function(t_h, mu_diff, mu_ave);
 
      auto t_D = get_D(2 * mu);
 
      t_inertia_force(i) = (rho * alpha) * (t_dot_u(i)); // ρ(h)∂U/∂t
      // t_buoyancy(SPACE_DIM - 1) = -(alpha * rho * a0) * t_h;
      t_gravity(SPACE_DIM - 1) = (alpha * rho * a0); // ρ(h)a_0 (a_0 is gravity)
      t_phase_force(i) = -alpha * kappa * t_g * t_grad_h(i); // -κg∇h
      t_convection(i) = (rho * alpha) * (t_u(j) * t_grad_u(i, j)); // ρ(h)(U·∇)U
 
      t_stress(i, j) =
          alpha * (t_D(i, j, k, l) * t_grad_u(k, l) + t_kd(i, j) * t_p); // ∇·(2μ(h)D(U))
 
      auto t_nf = getFTensor1FromArray<U_FIELD_DIM, U_FIELD_DIM>(locF);
 
      t_forces(i) = t_inertia_force(i) + t_buoyancy(i) + t_gravity(i) +
                    t_convection(i) + t_phase_force(i);
 
      int bb = 0;
      for (; bb != nbRows / U_FIELD_DIM; ++bb) {
 
        // Assemble total force contributions into local residual
        t_nf(i) += t_base * t_forces(i);
        t_nf(i) += t_diff_base(j) * t_stress(i, j); // ∇w : ∇·(2μ(h)D(U)) - ∇p
 
        if (coord_type == CYLINDRICAL) {
          t_nf(0) += (t_base * (alpha / t_coords(0))) * (2 * mu * t_u(0) + t_p);
        }
 
        ++t_base;
        ++t_diff_base;
        ++t_nf;
      }
 
      for (; bb < nbRowBaseFunctions; ++bb) {
        ++t_diff_base;
        ++t_base;
      }
 
      ++t_dot_u;
      ++t_u;
      ++t_grad_u;
      ++t_h;
      ++t_grad_h;
      ++t_g;
      ++t_p;
 
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dotUPtr

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsU::dotUPtr

private

Definition at line 580 of file FreeSurfaceOps.hpp.

gPtr

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpRhsU::gPtr

private

Definition at line 585 of file FreeSurfaceOps.hpp.

gradHPtr

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

private

Definition at line 584 of file FreeSurfaceOps.hpp.

gradUPtr

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpRhsU::gradUPtr

private

Definition at line 582 of file FreeSurfaceOps.hpp.

hPtr

boost::shared_ptr<VectorDouble> FreeSurfaceOps::OpRhsU::hPtr

private

Definition at line 583 of file FreeSurfaceOps.hpp.

pPtr

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

private

Definition at line 586 of file FreeSurfaceOps.hpp.

uPtr

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

private

Definition at line 581 of file FreeSurfaceOps.hpp.

---

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

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