ComponentTransportProcessData.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 <Eigen/Core>
#include <memory>
 
#include "ChemistryLib/ChemicalSolverInterface.h"
#include "LookupTable.h"
#include "MaterialLib/MPL/MaterialSpatialDistributionMap.h"
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "NumLib/Fem/ShapeMatrixCache.h"
#include "NumLib/NumericalStability/NumericalStabilization.h"
#include "ParameterLib/ConstantParameter.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
class Medium;
}
 
namespace MeshLib
{
template <typename PROP_VAL_TYPE>
class PropertyVector;
}
 
namespace ProcessLib
{
namespace ComponentTransport
{
struct ComponentTransportProcessData
{
    MaterialPropertyLib::MaterialSpatialDistributionMap media_map;
    bool const has_gravity;
    bool const non_advective_form;
    ParameterLib::Parameter<double> const* const temperature;
    bool const chemically_induced_porosity_change;
    ChemistryLib::ChemicalSolverInterface* const chemical_solver_interface;
    std::unique_ptr<LookupTable> lookup_table;
 
    NumLib::NumericalStabilization stabilizer;

    std::vector<Eigen::VectorXd> const projected_specific_body_force_vectors;
 
    int const mesh_space_dimension;

    ParameterLib::Parameter<double> const& aperture_size =
        ParameterLib::ConstantParameter<double>("constant_one", 1.0);
 
    bool const isothermal;

    NumLib::ShapeMatrixCache shape_matrix_cache;
 
    static const int hydraulic_process_id = 0;
    // Thermal process is optional, indicated by -1. If present, it is positive.
    const int thermal_process_id = isothermal ? -1 : 1;
    // TODO (renchao-lu): This variable is used in the calculation of the
    // fluid's density and flux, indicating the transport process id. For now it
    // is assumed that these quantities depend on the first occurring transport
    // process only. The density and flux calculations have to be extended to
    // all processes.
    const int first_transport_process_id = isothermal ? 1 : 2;
 
    MeshLib::PropertyVector<double>* mesh_prop_velocity = nullptr;
    MeshLib::PropertyVector<double>* mesh_prop_porosity = nullptr;
};
 
}  // namespace ComponentTransport
}  // namespace ProcessLib