SolutionDependentDirichletBoundaryCondition.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <string>
 
#include "BoundaryCondition.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
#include "ParameterLib/MeshNodeParameter.h"
 
namespace ProcessLib
{
struct SolutionDependentDirichletBoundaryConditionConfig
{
    std::string property_name;
    std::string initial_value_parameter_string;
};

class SolutionDependentDirichletBoundaryCondition final
    : public BoundaryCondition
{
public:
    SolutionDependentDirichletBoundaryCondition(
        std::string property_name,
        ParameterLib::Parameter<double> const& parameter,
        MeshLib::Mesh const& bc_mesh,
        NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
        int const variable_id, int const component_id);
 
    void getEssentialBCValues(
        double const t, GlobalVector const& x,
        NumLib::IndexValueVector<GlobalIndexType>& bc_values) const override;

    void postTimestep(double const /*t*/,
                      std::vector<GlobalVector*> const& x,
                      int const process_id) override;
 
private:
    MeshLib::Mesh const& _bc_mesh;
    int const _variable_id;
    int const _component_id;
    std::unique_ptr<NumLib::LocalToGlobalIndexMap const> _dof_table_boundary;
    std::unique_ptr<ParameterLib::MeshNodeParameter<double>> _parameter;
    // avoid repetitively retrieving the desired property from the boundary mesh
    // over time steps
    MeshLib::PropertyVector<double>* _solution_dependent_bc;
};
 
SolutionDependentDirichletBoundaryConditionConfig
parseSolutionDependentDirichletBoundaryCondition(
    BaseLib::ConfigTree const& config);
 
std::unique_ptr<SolutionDependentDirichletBoundaryCondition>
createSolutionDependentDirichletBoundaryCondition(
    SolutionDependentDirichletBoundaryConditionConfig const& config,
    MeshLib::Mesh const& bc_mesh,
    NumLib::LocalToGlobalIndexMap const& dof_table_bulk, int const variable_id,
    int const component_id,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
        parameters);
 
}  // namespace ProcessLib