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

Detailed Description

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

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

Definition at line 163 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 >::AnisotropicCalculator::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 374 of file LiquidFlowLocalAssembler-impl.h.

{
    const double fac = ip_data.integration_weight / mu;
    local_K.noalias() += fac * ip_data.dNdx.transpose() *
                         permeability_with_density_factor * ip_data.dNdx;
 
    if (has_gravity)
    {
        local_b.noalias() += (fac * rho_L) * ip_data.dNdx.transpose() *
                             permeability_with_density_factor *
                             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 >::AnisotropicCalculator::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 397 of file LiquidFlowLocalAssembler-impl.h.

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

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

---

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

• ProcessLib/LiquidFlow/LiquidFlowLocalAssembler.h
• ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h