Revision 987c1d8c0cb1df40354ad8acca6c72c229242d3f authored by Thomas Fischer on 01 September 2021, 12:37:54 UTC, committed by Dmitry Yu. Naumov on 05 September 2021, 20:09:22 UTC
1 parent a33bf91
ComponentTransportProcessData.h
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*/
#pragma once
#include <memory>
#include "ChemistryLib/ChemicalSolverInterface.h"
#include "LookupTable.h"
#include "MaterialLib/MPL/CreateMaterialSpatialDistributionMap.h"
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "ParameterLib/Parameter.h"
namespace MaterialPropertyLib
{
class Medium;
}
namespace MeshLib
{
template <typename PROP_VAL_TYPE>
class PropertyVector;
}
namespace ProcessLib
{
namespace ComponentTransport
{
struct ComponentTransportProcessData
{
std::unique_ptr<MaterialPropertyLib::MaterialSpatialDistributionMap>
media_map;
Eigen::VectorXd const specific_body_force;
bool const has_gravity;
bool const non_advective_form;
/// This optional tag provides a simple means of considering the temperature
/// effect on the solute transport process.
ParameterLib::Parameter<double> const* const temperature;
/**
* When this optional tag is on, the feedback of chemical reactions on the
* porosity will be counted. The change of porosity equals to the summation
* over the changes in the volume fractions of solid constituents. The
* change of the volume fraction, in terms of a solid constituent, results
* from chemical reactions.
*
* \note In order to use this optional tag, the amount of solid constituents
* should be given as volume fraction instead of molality. In addition,
* an appropriate molar volume is required for each solid. The relationship
* to calculate volume fractions of m solids from molalities is as follows:
* \f[
* b_i = \frac{n_i}{m^l} = \frac{\phi_i}{\rho^l \phi V_{m,i}}, i=1,...,m
* \f]
* where \f$b_i\f$ is the molality in mol/kg of water,
* \f$n_i\f$ is the amount of solid in mol,
* \f$m^l\f$ is the mass of water in kg,
* \f$\phi_i\f$ is the volume fraction of solid i,
* \f$\rho^l\f$ is the density of water in kg/m\f$^3\f$,
* \f$\phi\f$ is the porosity,
* \f$V_{m,i}\f$ is the molar volume of solid i in m\f$^3\f$/mol.
*/
bool const chemically_induced_porosity_change;
ChemistryLib::ChemicalSolverInterface* const chemical_solver_interface;
std::unique_ptr<LookupTable> lookup_table;
const int hydraulic_process_id;
// 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;
MeshLib::PropertyVector<double>* mesh_prop_velocity = nullptr;
};
} // namespace ComponentTransport
} // namespace ProcessLib
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...