ProcessLib::TwoPhaseFlowWithPrho Namespace Reference

Classes
struct	IntegrationPointData
class	TwoPhaseFlowWithPrhoLocalAssemblerInterface
class	TwoPhaseFlowWithPrhoLocalAssembler
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, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const &parameters)
std::unique_ptr< Process >	createTwoPhaseFlowWithPrhoProcess (std::string 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< std::string, std::unique_ptr< MathLib::PiecewiseLinearInterpolation >> const &curves)

Variables
const unsigned	NUM_NODAL_DOF = 2

Function Documentation

createTwoPhaseFlowPrhoMaterialProperties()

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

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__fluid

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__liquid_density

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__gas_density

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__liquid_viscosity

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__gas_viscosity

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__permeability

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__porosity

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__material_property__porous_medium__porous_medium__storage

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 34 of file CreateTwoPhaseFlowPrhoMaterialProperties.cpp.

{
    DBUG("Reading material properties of two-phase flow process.");
 
    auto const& fluid_config = config.getConfigSubtree("fluid");
 
    // Get fluid properties
    auto const& rho_conf = fluid_config.getConfigSubtree("liquid_density");
    auto _liquid_density =
        MaterialLib::Fluid::createFluidDensityModel(rho_conf);
    auto const& rho_gas_conf = fluid_config.getConfigSubtree("gas_density");
    auto _gas_density =
        MaterialLib::Fluid::createFluidDensityModel(rho_gas_conf);
    auto const& mu_conf = fluid_config.getConfigSubtree("liquid_viscosity");
    auto _viscosity = MaterialLib::Fluid::createViscosityModel(mu_conf);
    auto const& mu_gas_conf = fluid_config.getConfigSubtree("gas_viscosity");
    auto _gas_viscosity = MaterialLib::Fluid::createViscosityModel(mu_gas_conf);
 
    // Get porous properties
    std::vector<int> mat_ids;
    std::vector<int> mat_krel_ids;
    std::vector<std::unique_ptr<MaterialLib::PorousMedium::Permeability>>
        _intrinsic_permeability_models;
    std::vector<std::unique_ptr<MaterialLib::PorousMedium::Porosity>>
        _porosity_models;
    std::vector<std::unique_ptr<MaterialLib::PorousMedium::Storage>>
        _storage_models;
    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& permeability_conf = conf.getConfigSubtree("permeability");
        _intrinsic_permeability_models.emplace_back(
            MaterialLib::PorousMedium::createPermeabilityModel(
                permeability_conf, parameters));
 
        auto const& porosity_conf = conf.getConfigSubtree("porosity");
        auto n = MaterialLib::PorousMedium::createPorosityModel(porosity_conf,
                                                                parameters);
        _porosity_models.emplace_back(std::move(n));
 
        auto const& storage_conf = conf.getConfigSubtree("storage");
        auto beta = MaterialLib::PorousMedium::createStorageModel(storage_conf);
        _storage_models.emplace_back(std::move(beta));
 
        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);
    }
 
    BaseLib::reorderVector(_intrinsic_permeability_models, mat_ids);
    BaseLib::reorderVector(_porosity_models, mat_ids);
    BaseLib::reorderVector(_storage_models, mat_ids);
 
    return std::make_unique<TwoPhaseFlowWithPrhoMaterialProperties>(
        material_ids, std::move(_liquid_density), std::move(_viscosity),
        std::move(_gas_density), std::move(_gas_viscosity),
        std::move(_intrinsic_permeability_models), std::move(_porosity_models),
        std::move(_storage_models), std::move(_capillary_pressure_models),
        std::move(_relative_permeability_models));
}

References MaterialPropertyLib::beta, MaterialLib::PorousMedium::createCapillaryPressureModel(), MaterialLib::Fluid::createFluidDensityModel(), MaterialLib::PorousMedium::createPermeabilityModel(), MaterialLib::PorousMedium::createPorosityModel(), MaterialLib::PorousMedium::createRelativePermeabilityModel(), MaterialLib::PorousMedium::createStorageModel(), MaterialLib::Fluid::createViscosityModel(), DBUG(), BaseLib::ConfigTree::getConfigAttributeOptional(), BaseLib::ConfigTree::getConfigSubtree(), and BaseLib::reorderVector().

Referenced by createTwoPhaseFlowWithPrhoProcess().

createTwoPhaseFlowWithPrhoProcess()

std::unique_ptr< Process > ProcessLib::TwoPhaseFlowWithPrho::createTwoPhaseFlowWithPrhoProcess	(	std::string	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< std::string, std::unique_ptr< MathLib::PiecewiseLinearInterpolation >> const &	curves
	)

Input File Parameter:

prj__processes__process__type

Input File Parameter:

prj__processes__process__TWOPHASE_FLOW_PRHO__process_variables

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

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 26 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.");
    }
 
    std::unique_ptr<TwoPhaseFlowWithPrhoMaterialProperties> material =
        createTwoPhaseFlowPrhoMaterialProperties(mat_config, material_ids,
                                                 parameters);
 
    TwoPhaseFlowWithPrhoProcessData process_data{
        specific_body_force, has_gravity, mass_lumping,       diff_coeff_b,
        diff_coeff_a,        temperature, std::move(material)};
 
    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), mat_config,
        curves);
}

References BaseLib::ConfigTree::checkConfigParameter(), ProcessLib::createSecondaryVariables(), createTwoPhaseFlowPrhoMaterialProperties(), DBUG(), ProcessLib::findProcessVariables(), BaseLib::ConfigTree::getConfigParameter(), BaseLib::ConfigTree::getConfigSubtree(), INFO(), MeshLib::materialIDs(), MaterialPropertyLib::name, and MathLib::toVector().

Referenced by ProjectData::parseProcesses().

Variable Documentation

NUM_NODAL_DOF

const unsigned ProcessLib::TwoPhaseFlowWithPrho::NUM_NODAL_DOF = 2

Definition at line 59 of file TwoPhaseFlowWithPrhoLocalAssembler.h.

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