ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::IsotropicCalculator Struct Reference

Detailed Description

template<typename ShapeFunction, int GlobalDim>
struct ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::IsotropicCalculator

Calculator of the Laplacian and the gravity term for isotropic permeability tensor

Definition at line 176 of file LiquidFlowLocalAssembler.h.

Static Public Member Functions
static void	calculateLaplacianAndGravityTerm (Eigen::Map< NodalMatrixType > &local_K, Eigen::Map< NodalVectorType > &local_b, IntegrationPointData< GlobalDimNodalMatrixType > const &ip_data, GlobalDimMatrixType const &permeability_with_density_factor, double const mu, double const rho_L, GlobalDimVectorType const &specific_body_force, bool const has_gravity)
static Eigen::Matrix< double, GlobalDim, 1 >	calculateVelocity (Eigen::Map< const NodalVectorType > const &local_p, IntegrationPointData< GlobalDimNodalMatrixType > const &ip_data, GlobalDimMatrixType const &permeability, double const mu, double const rho_L, GlobalDimVectorType const &specific_body_force, bool const has_gravity)

Member Function Documentation

calculateLaplacianAndGravityTerm()

template<typename ShapeFunction, int GlobalDim>

void ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::IsotropicCalculator::calculateLaplacianAndGravityTerm ( Eigen::Map< NodalMatrixType > & local_K, Eigen::Map< NodalVectorType > & local_b, IntegrationPointData< GlobalDimNodalMatrixType > const & ip_data, GlobalDimMatrixType const & permeability_with_density_factor, double const mu, double const rho_L, GlobalDimVectorType const & specific_body_force, bool const has_gravity )

static

Definition at line 355 of file LiquidFlowLocalAssembler-impl.h.

{
    const double K = permeability_with_density_factor(0, 0) / mu;
    const double fac = K * ip_data.integration_weight;
    local_K.noalias() += fac * ip_data.dNdx.transpose() * ip_data.dNdx;
 
    if (has_gravity)
    {
        local_b.noalias() +=
            (fac * rho_L) * ip_data.dNdx.transpose() * specific_body_force;
    }
}

References ProcessLib::LiquidFlow::IntegrationPointData< GlobalDimNodalMatrixType >::dNdx, and ProcessLib::LiquidFlow::IntegrationPointData< GlobalDimNodalMatrixType >::integration_weight.

calculateVelocity()

template<typename ShapeFunction, int GlobalDim>

Eigen::Matrix< double, GlobalDim, 1 > ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::IsotropicCalculator::calculateVelocity ( Eigen::Map< const NodalVectorType > const & local_p, IntegrationPointData< GlobalDimNodalMatrixType > const & ip_data, GlobalDimMatrixType const & permeability, double const mu, double const rho_L, GlobalDimVectorType const & specific_body_force, bool const has_gravity )

static

Definition at line 377 of file LiquidFlowLocalAssembler-impl.h.

{
    const double K = permeability(0, 0) / mu;
    // Compute the velocity
    GlobalDimVectorType velocity = -K * ip_data.dNdx * local_p;
    // gravity term
    if (has_gravity)
    {
        velocity += (K * rho_L) * specific_body_force;
    }
    return velocity;
}

References ProcessLib::LiquidFlow::IntegrationPointData< GlobalDimNodalMatrixType >::dNdx.

---

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

• LiquidFlowLocalAssembler.h
• LiquidFlowLocalAssembler-impl.h