ChemicalSolverInterface.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/MPL/VariableType.h"
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
 
namespace MaterialPropertyLib
{
class Medium;
}
 
namespace ParameterLib
{
class SpatialPosition;
}
 
namespace ChemistryLib
{
class ChemicalSolverInterface
{
public:
    ChemicalSolverInterface(GlobalLinearSolver& linear_solver_)
        : linear_solver(linear_solver_)
    {
    }
 
    virtual void initialize() {}
 
    virtual void initializeChemicalSystemConcrete(
        std::vector<double> const& /*concentrations*/,
        GlobalIndexType const& /*chemical_system_id*/,
        MaterialPropertyLib::Medium const& /*medium*/,
        ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/)
    {
    }
 
    virtual void setChemicalSystemConcrete(
        std::vector<double> const& /*concentrations*/,
        GlobalIndexType const& /*chemical_system_id*/,
        MaterialPropertyLib::Medium const* /*medium*/,
        MaterialPropertyLib::VariableArray const& /*vars*/,
        ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/,
        double const /*dt*/)
    {
    }
 
    virtual void setAqueousSolutionsPrevFromDumpFile() {}
 
    virtual void executeSpeciationCalculation(double const dt) = 0;
 
    virtual std::vector<GlobalVector*> getIntPtProcessSolutions() const = 0;
 
    virtual double getConcentration(
        int const /*component_id*/,
        GlobalIndexType const /*chemical_system_id*/) const
    {
        return std::numeric_limits<double>::quiet_NaN();
    }
 
    virtual std::vector<std::string> const getComponentList() const
    {
        return {};
    }
 
    virtual void updateVolumeFractionPostReaction(
        GlobalIndexType const& /*chemical_system_id*/,
        MaterialPropertyLib::Medium const& /*medium*/,
        ParameterLib::SpatialPosition const& /*pos*/, double const /*porosity*/,
        double const /*t*/, double const /*dt*/)
    {
    }
 
    virtual void updatePorosityPostReaction(
        GlobalIndexType const& /*chemical_system_id*/,
        MaterialPropertyLib::Medium const& /*medium*/,
        double& /*porosity*/)
    {
    }
 
    virtual void computeSecondaryVariable(
        std::size_t const /*ele_id*/,
        std::vector<GlobalIndexType> const& /*chemical_system_indices*/)
    {
    }
 
    virtual ~ChemicalSolverInterface() = default;
 
public:
    std::vector<GlobalIndexType> chemical_system_index_map;
    GlobalLinearSolver& linear_solver;
};
}  // namespace ChemistryLib