ProcessLib::RichardsComponentTransport Namespace Reference

Namespaces
	anonymous_namespace{CreateRichardsComponentTransportProcess.cpp}

Classes
struct	IntegrationPointData
class	RichardsComponentTransportLocalAssemblerInterface
class	LocalAssemblerData
class	RichardsComponentTransportProcess
struct	RichardsComponentTransportProcessData

Functions
std::unique_ptr< Process >	createRichardsComponentTransportProcess (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< int, std::shared_ptr< MaterialPropertyLib::Medium >> const &media)

Variables
const unsigned	NUM_NODAL_DOF = 2

Function Documentation

createRichardsComponentTransportProcess()

std::unique_ptr< Process > ProcessLib::RichardsComponentTransport::createRichardsComponentTransportProcess	(	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< int, std::shared_ptr< MaterialPropertyLib::Medium >> const &	media
	)

Input File Parameter:

prj__processes__process__type

Input File Parameter:

prj__processes__process__RichardsComponentTransport__coupling_scheme

Input File Parameter:

prj__processes__process__RichardsComponentTransport__process_variables

Input File Parameter:

prj__processes__process__RichardsComponentTransport__process_variables__concentration

Input File Parameter:

prj__processes__process__RichardsComponentTransport__process_variables__pressure

Input File Parameter:

prj__processes__process__RichardsComponentTransport__specific_body_force

Definition at line 79 of file CreateRichardsComponentTransportProcess.cpp.

{
    config.checkConfigParameter("type", "RichardsComponentTransport");
 
    DBUG("Create RichardsComponentTransportProcess.");
 
    auto const coupling_scheme =
        config.getConfigParameterOptional<std::string>("coupling_scheme");
    const bool use_monolithic_scheme =
        !(coupling_scheme && (*coupling_scheme == "staggered"));
 
    // Process variable.
 
    auto const pv_config = config.getConfigSubtree("process_variables");
 
    std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
        process_variables;
    if (use_monolithic_scheme)  // monolithic scheme.
    {
        auto per_process_variables = findProcessVariables(
            variables, pv_config,
            {
             "concentration",
             "pressure"});
        if (per_process_variables.size() > 2)
        {
            OGS_FATAL(
                "By now RichardsComponentTransport process only supports "
                "single component transport simulation.");
        }
        process_variables.push_back(std::move(per_process_variables));
    }
    else  // staggered scheme.
    {
        OGS_FATAL("The staggered coupling scheme is not implemented.");
    }
 
    // Specific body force parameter.
    Eigen::VectorXd specific_body_force;
    std::vector<double> const b =
        config.getConfigParameter<std::vector<double>>("specific_body_force");
    assert(!b.empty() && b.size() < 4);
    if (b.size() < mesh.getDimension())
    {
        OGS_FATAL(
            "specific body force (gravity vector) has {:d} components, mesh "
            "dimension is {:d}",
            b.size(), mesh.getDimension());
    }
    bool const has_gravity = MathLib::toVector(b).norm() > 0;
    if (has_gravity)
    {
        specific_body_force.resize(b.size());
        std::copy_n(b.data(), b.size(), specific_body_force.data());
    }
 
    auto media_map =
        MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh);
 
    DBUG(
        "Check the media properties of RichardsComponentTransport process ...");
    checkMPLProperties(mesh, *media_map);
    DBUG("Media properties verified.");
 
    RichardsComponentTransportProcessData process_data{
        std::move(media_map), specific_body_force, has_gravity};
 
    SecondaryVariableCollection secondary_variables;
 
    ProcessLib::createSecondaryVariables(config, secondary_variables);
 
    return std::make_unique<RichardsComponentTransportProcess>(
        std::move(name), mesh, std::move(jacobian_assembler), parameters,
        integration_order, std::move(process_variables),
        std::move(process_data), std::move(secondary_variables),
        use_monolithic_scheme);
}

References BaseLib::ConfigTree::checkConfigParameter(), ProcessLib::RichardsComponentTransport::anonymous_namespace{CreateRichardsComponentTransportProcess.cpp}::checkMPLProperties(), MaterialPropertyLib::createMaterialSpatialDistributionMap(), ProcessLib::createSecondaryVariables(), DBUG(), ProcessLib::findProcessVariables(), BaseLib::ConfigTree::getConfigParameter(), BaseLib::ConfigTree::getConfigParameterOptional(), BaseLib::ConfigTree::getConfigSubtree(), MeshLib::Mesh::getDimension(), MaterialPropertyLib::name, OGS_FATAL, and MathLib::toVector().

Referenced by ProjectData::parseProcesses().

Variable Documentation

NUM_NODAL_DOF

const unsigned ProcessLib::RichardsComponentTransport::NUM_NODAL_DOF = 2

Definition at line 55 of file RichardsComponentTransportFEM.h.

Referenced by ProcessLib::RichardsComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::LocalAssemblerData(), and ProcessLib::RichardsComponentTransport::LocalAssemblerData< ShapeFunction, IntegrationMethod, GlobalDim >::assemble().