FreeSurfaceOps::OpLhsU_dH Struct Reference

Lhs for U dH (Jacobian block ∂R_U/∂H) More...

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

Inheritance diagram for FreeSurfaceOps::OpLhsU_dH:

Collaboration diagram for FreeSurfaceOps::OpLhsU_dH:

Public Member Functions
	OpLhsU_dH (const std::string field_name_u, const std::string field_name_h, 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< VectorDouble > g_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_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< VectorDouble >	gPtr

Detailed Description

Lhs for U dH (Jacobian block ∂R_U/∂H)

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)
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)
g_ptr	Pointer to the chemical potential vector (g)

Returns

MoFEMErrorCode

This operator builds the sensitivity of the momentum residual to changes in the phase field h (paper 3.2a cross-coupling). It implements:

In code: rho_dh = d_phase_function_h(t_h, rho_diff), mu_dh = d_phase_function_h(t_h, mu_diff)

Examples

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

Definition at line 754 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpLhsU_dH()

FreeSurfaceOps::OpLhsU_dH::OpLhsU_dH ( const std::string  field_name_u, const std::string  field_name_h, 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< VectorDouble >  g_ptr  )

inline

Definition at line 756 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name_u, field_name_h,
                            AssemblyDomainEleOp::OPROWCOL),
        dotUPtr(dot_u_ptr), uPtr(u_ptr), gradUPtr(grad_u_ptr), hPtr(h_ptr),
        gPtr(g_ptr) {
    sYmm = false;
    assembleTranspose = false;
  }

Member Function Documentation

iNtegrate()

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

inline

Definition at line 770 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_grad_u = getFTensor2FromMat<U_FIELD_DIM, SPACE_DIM>(*gradUPtr);
    auto t_h = getFTensor0FromVec(*hPtr);
    auto t_g = getFTensor0FromVec(*gPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_row_base = row_data.getFTensor0N();
    auto t_row_diff_base = row_data.getFTensor1DiffN<SPACE_DIM>();
 
    auto t_w = getFTensor0IntegrationWeight();
 
    FTensor::Tensor2_symmetric<double, SPACE_DIM> t_stress_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_phase_force_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_inertia_force_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_convection_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_buoyancy_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_gravity_dh;
    FTensor::Tensor1<double, U_FIELD_DIM> t_forces_dh;
 
    t_buoyancy_dh(i) = 0;
    t_gravity_dh(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_dh = d_phase_function_h(t_h, rho_diff); // dρ/dh
      const double mu_dh = d_phase_function_h(t_h, mu_diff); // dμ/dh
 
      auto t_D_dh = get_D(2 * mu_dh);
 
      t_inertia_force_dh(i) = (alpha * rho_dh) * t_dot_u(i); // dρ/dh * ∂U/∂t
      // t_buoyancy_dh(SPACE_DIM - 1) = -(alpha * a0) * (rho + rho_dh * t_h);
      t_gravity_dh(SPACE_DIM - 1) = (alpha * rho_dh * a0); // dρ/dh * a_0
      t_convection_dh(i) = (rho_dh * alpha) * (t_u(j) * t_grad_u(i, j)); // dρ/dh * (U·∇)U
      const double t_phase_force_g_dh = -alpha * kappa * t_g; // -κ g
      t_forces_dh(i) = t_inertia_force_dh(i) + t_buoyancy_dh(i) +
                       t_gravity_dh(i) + t_convection_dh(i); // Total dR_U/dh force
 
      t_stress_dh(i, j) = alpha * (t_D_dh(i, j, k, l) * t_grad_u(k, l)); // d(∇·(2μD(U)))/dh
 
      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);
        auto t_col_diff_base = col_data.getFTensor1DiffN<SPACE_DIM>(gg, 0);
 
        for (int cc = 0; cc != nbCols; ++cc) {
 
          const double bb = t_row_base * t_col_base;
          t_mat(i) += t_forces_dh(i) * bb; // Assemble total dR_U/dh force contribution
          t_mat(i) += (t_phase_force_g_dh * t_row_base) * t_col_diff_base(i); // -κ g ∇h term
          t_mat(i) += (t_row_diff_base(j) * t_col_base) * t_stress_dh(i, j); // Viscous term contribution via dμ/dh
 
          if (coord_type == CYLINDRICAL) {
            t_mat(0) += (bb * (alpha / t_coords(0))) * (2 * mu_dh * t_u(0));
          }
 
          ++t_mat;
          ++t_col_base;
          ++t_col_diff_base;
        }
 
        ++t_row_base;
        ++t_row_diff_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_diff_base;
        ++t_row_base;
      }
 
      ++t_dot_u;
      ++t_u;
      ++t_grad_u;
      ++t_h;
      ++t_g;
      ++t_coords;
      ++t_w;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dotUPtr

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

private

Definition at line 866 of file FreeSurfaceOps.hpp.

gPtr

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

private

Definition at line 870 of file FreeSurfaceOps.hpp.

gradUPtr

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

private

Definition at line 868 of file FreeSurfaceOps.hpp.

hPtr

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

private

Definition at line 869 of file FreeSurfaceOps.hpp.

uPtr

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

private

Definition at line 867 of file FreeSurfaceOps.hpp.

---

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

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