Classes
struct	IntegrationPointData

class	TwoPhaseFlowWithPrhoLocalAssembler

class	TwoPhaseFlowWithPrhoLocalAssemblerInterface

class	TwoPhaseFlowWithPrhoMaterialProperties

class	TwoPhaseFlowWithPrhoProcess
	A class to simulate the two-phase flow process with P-rho model in porous media. More...

struct	TwoPhaseFlowWithPrhoProcessData

Functions
std::unique_ptr< TwoPhaseFlowWithPrhoMaterialProperties >	createTwoPhaseFlowPrhoMaterialProperties (BaseLib::ConfigTree const &config, MeshLib::PropertyVector< int > const *const material_ids)

void	checkMPLProperties (std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)

std::unique_ptr< Process >	createTwoPhaseFlowWithPrhoProcess (std::string const &name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)

Variables
const unsigned	NUM_NODAL_DOF = 2

Function Documentation

◆ checkMPLProperties()

void ProcessLib::TwoPhaseFlowWithPrho::checkMPLProperties ( std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const & media )

Definition at line 28 of file CreateTwoPhaseFlowWithPrhoProcess.cpp.

{
    std::array const required_medium_properties = {
        MaterialPropertyLib::permeability, MaterialPropertyLib::porosity};
    std::array const required_liquid_properties = {
        MaterialPropertyLib::viscosity, MaterialPropertyLib::density};
    std::array const required_gas_properties = {
        MaterialPropertyLib::viscosity, MaterialPropertyLib::density,
        MaterialPropertyLib::molar_mass};
 
    for (auto const& m : media)
    {
        checkRequiredProperties(*m.second, required_medium_properties);
        checkRequiredProperties(m.second->phase("AqueousLiquid"),
                                required_liquid_properties);
        checkRequiredProperties(m.second->phase("Gas"),
                                required_gas_properties);
    }
}

References MaterialPropertyLib::density, MaterialPropertyLib::molar_mass, MaterialPropertyLib::permeability, MaterialPropertyLib::porosity, and MaterialPropertyLib::viscosity.

Referenced by createTwoPhaseFlowWithPrhoProcess().

◆ createTwoPhaseFlowPrhoMaterialProperties()

std::unique_ptr< TwoPhaseFlowWithPrhoMaterialProperties > ProcessLib::TwoPhaseFlowWithPrho::createTwoPhaseFlowPrhoMaterialProperties	(	BaseLib::ConfigTree const &	config,
		MeshLib::PropertyVector< int > const *const	material_ids )

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__id

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__capillary_pressure

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__relative_permeability

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__relative_permeability__relative_permeability

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__relative_permeability__relative_permeability__id

Definition at line 31 of file CreateTwoPhaseFlowPrhoMaterialProperties.cpp.

{
    DBUG("Reading material properties of two-phase flow process.");
 
    // Get porous properties
    std::vector<int> mat_ids;
    std::vector<int> mat_krel_ids;
    std::vector<
        std::unique_ptr<MaterialLib::PorousMedium::CapillaryPressureSaturation>>
        _capillary_pressure_models;
    std::vector<
        std::unique_ptr<MaterialLib::PorousMedium::RelativePermeability>>
        _relative_permeability_models;
 
    auto const& poro_config = config.getConfigSubtree("porous_medium");
    for (auto const& conf : poro_config.getConfigSubtreeList("porous_medium"))
    {
        auto const id = conf.getConfigAttributeOptional<int>("id");
        mat_ids.push_back(*id);
 
        auto const& capillary_pressure_conf =
            conf.getConfigSubtree("capillary_pressure");
        auto pc = MaterialLib::PorousMedium::createCapillaryPressureModel(
            capillary_pressure_conf);
        _capillary_pressure_models.emplace_back(std::move(pc));
 
        auto const& krel_config =
            conf.getConfigSubtree("relative_permeability");
        for (
            auto const& krel_conf :
            krel_config.getConfigSubtreeList("relative_permeability"))
        {
            auto const krel_id =
                krel_conf.getConfigAttributeOptional<int>("id");
            mat_krel_ids.push_back(*krel_id);
            auto krel_n =
                MaterialLib::PorousMedium::createRelativePermeabilityModel(
                    krel_conf);
            _relative_permeability_models.emplace_back(std::move(krel_n));
        }
        BaseLib::reorderVector(_relative_permeability_models, mat_krel_ids);
    }
 
    return std::make_unique<TwoPhaseFlowWithPrhoMaterialProperties>(
        material_ids,
        std::move(_capillary_pressure_models),
        std::move(_relative_permeability_models));
}

References MaterialLib::PorousMedium::createCapillaryPressureModel(), MaterialLib::PorousMedium::createRelativePermeabilityModel(), DBUG(), BaseLib::ConfigTree::getConfigAttributeOptional(), BaseLib::ConfigTree::getConfigSubtree(), and BaseLib::reorderVector().

Referenced by createTwoPhaseFlowWithPrhoProcess().

◆ createTwoPhaseFlowWithPrhoProcess()

std::unique_ptr< Process > ProcessLib::TwoPhaseFlowWithPrho::createTwoPhaseFlowWithPrhoProcess	(	std::string const &	name,
		MeshLib::Mesh &	mesh,
		std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&	jacobian_assembler,
		std::vector< ProcessVariable > const &	variables,
		std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &	parameters,
		unsigned const	integration_order,
		BaseLib::ConfigTree const &	config,
		std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &	media )

Input File Parameter: prj__processes__process__type

Process Variables

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__process_variables

Primary process variables as they appear in the global component vector:

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__process_variables__liquid_pressure

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__process_variables__overall_mass_density

Process Parameters

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__specific_body_force

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__mass_lumping

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__diffusion_coeff_component_b

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__diffusion_coeff_component_a

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__temperature

Input File Parameter: prj__processes__process__TWOPHASE_FLOW_PRHO__material_property

Definition at line 49 of file CreateTwoPhaseFlowWithPrhoProcess.cpp.

{
    config.checkConfigParameter("type", "TWOPHASE_FLOW_PRHO");
 
    DBUG("Create TwoPhaseFlowProcess with Prho model.");
 
 
    auto const pv_config = config.getConfigSubtree("process_variables");
 
    auto per_process_variables = findProcessVariables(
        variables, pv_config,
        {
         "liquid_pressure",
         "overall_mass_density"});
    std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
        process_variables;
    process_variables.push_back(std::move(per_process_variables));
 
    SecondaryVariableCollection secondary_variables;
 
    ProcessLib::createSecondaryVariables(config, secondary_variables);
    // Specific body force
    std::vector<double> const b =
        config.getConfigParameter<std::vector<double>>("specific_body_force");
    assert(!b.empty() && b.size() < 4);
    Eigen::VectorXd specific_body_force(b.size());
    bool const has_gravity = MathLib::toVector(b).norm() > 0;
    if (has_gravity)
    {
        std::copy_n(b.data(), b.size(), specific_body_force.data());
    }
 
    auto const mass_lumping = config.getConfigParameter<bool>("mass_lumping");
    // diffusion coeff
    auto const& diff_coeff_b = ParameterLib::findParameter<double>(
        config,
        "diffusion_coeff_component_b", parameters, 1, &mesh);
    auto const& diff_coeff_a = ParameterLib::findParameter<double>(
        config,
        "diffusion_coeff_component_a", parameters, 1, &mesh);
    auto const& temperature = ParameterLib::findParameter<double>(
        config,
        "temperature", parameters, 1, &mesh);
 
    auto const& mat_config = config.getConfigSubtree("material_property");
 
    auto const material_ids = materialIDs(mesh);
    if (material_ids != nullptr)
    {
        INFO("The twophase flow is in heterogeneous porous media.");
    }
    else
    {
        INFO("The twophase flow is in homogeneous porous media.");
    }
 
    auto media_map =
        MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh);
    checkMPLProperties(media);
 
    std::unique_ptr<TwoPhaseFlowWithPrhoMaterialProperties> material =
        createTwoPhaseFlowPrhoMaterialProperties(mat_config, material_ids);
 
    TwoPhaseFlowWithPrhoProcessData process_data{
        specific_body_force, has_gravity,         mass_lumping,
        diff_coeff_b,        diff_coeff_a,        temperature,
        std::move(material), std::move(media_map)};
 
    return std::make_unique<TwoPhaseFlowWithPrhoProcess>(
        std::move(name), mesh, std::move(jacobian_assembler), parameters,
        integration_order, std::move(process_variables),
        std::move(process_data), std::move(secondary_variables));
}

References BaseLib::ConfigTree::checkConfigParameter(), checkMPLProperties(), MaterialPropertyLib::createMaterialSpatialDistributionMap(), ProcessLib::createSecondaryVariables(), createTwoPhaseFlowPrhoMaterialProperties(), DBUG(), ParameterLib::findParameter(), ProcessLib::findProcessVariables(), BaseLib::ConfigTree::getConfigParameter(), BaseLib::ConfigTree::getConfigSubtree(), INFO(), and MathLib::toVector().

Referenced by ProjectData::parseProcesses().

Variable Documentation

◆ NUM_NODAL_DOF

const unsigned ProcessLib::TwoPhaseFlowWithPrho::NUM_NODAL_DOF = 2

Definition at line 60 of file TwoPhaseFlowWithPrhoLocalAssembler.h.

Referenced by ProcessLib::TwoPhaseFlowWithPrho::TwoPhaseFlowWithPrhoLocalAssembler< ShapeFunction, GlobalDim >::assemble().

Classes

Functions

Variables