Homec4science
Diffusion SpecMiCP / ReactMiCP (single) Browse

specmicp/src/reactmicp/systems/saturated_diffusion/transport_neutrality_parameters.hpp
e882181496b7single

Last Change
Raw File

transport_neutrality_parameters.hpp
View Options

#ifndef SPECMICP_REACTMICP_SATURATEDDIFFUSIONNEUTRALITY_TRANSPORTPARAMETERS_HPP
#define SPECMICP_REACTMICP_SATURATEDDIFFUSIONNEUTRALITY_TRANSPORTPARAMETERS_HPP
#include "common.hpp"
namespace specmicp {
namespace reactmicp {
namespace systems {
namespace siasaturated {
//! \brief Parameters for the diffusion system
class SaturatedNeutralityDiffusionTransportParameters
{
public:
SaturatedNeutralityDiffusionTransportParameters(
uindex_t nb_nodes,
uindex_t nb_components,
uindex_t nb_secondary
):
m_diffusion_coefficient_component(nb_components),
m_diffusion_coefficient_secondary(nb_secondary),
m_reduction_factor(nb_nodes),
m_porosities(nb_nodes)
{}
SaturatedNeutralityDiffusionTransportParameters(
uindex_t nb_nodes,
uindex_t nb_component,
uindex_t nb_secondary,
scalar_t the_diffusion_coefficient,
scalar_t the_porosity,
scalar_t reduction_factor
):
SaturatedNeutralityDiffusionTransportParameters(nb_nodes, nb_component, nb_secondary)
{
m_diffusion_coefficient_component.setConstant(the_diffusion_coefficient);
m_diffusion_coefficient_secondary.setConstant(the_diffusion_coefficient);
m_reduction_factor.setConstant(reduction_factor);
m_porosities.setConstant(the_porosity);
}
//! \brief Density of water (kg/m^3)
scalar_t density_water() {return 1e3;}
//! Diffusion coefficient of component 'component' at 'node'
scalar_t diffusion_component(index_t node, index_t component) const {
return m_diffusion_coefficient_component(component)*m_reduction_factor(node);
}
//! Diffusion coefficient of secondary species 'secondary' at 'node'
scalar_t diffusion_secondary(index_t node, index_t secondary) const {
return m_diffusion_coefficient_secondary(secondary)*m_reduction_factor(node);
}
//! Reduction factor of node 'node'
scalar_t& reduction_factor(index_t node) {return m_reduction_factor(node);}
//! Diffusion coefficient of component 'component'
scalar_t& diff_coeff_component(index_t component) {
return m_diffusion_coefficient_component(component);
}
//! Diffusion coefficient of secondary species 'secondary'
scalar_t& diff_coeff_secondary(index_t secondary) {
return m_diffusion_coefficient_secondary(secondary);
}
//! \brief Porosity (at a node (function of composition))
scalar_t porosity(index_t node) const {return m_porosities(node);}
//! \brief Porosity at a node
scalar_t& porosity(index_t node) {return m_porosities(node);}
private:
Vector m_diffusion_coefficient_component;
Vector m_diffusion_coefficient_secondary;
Vector m_reduction_factor;
Vector m_porosities;
};
} // end namespace siasaturated
} // end namespace systems
} // end namespace reactmicp
} // end namespace specmicp
#endif // SPECMICP_REACTMICP_SATURATEDDIFFUSIONNEUTRALITY_TRANSPORTPARAMETERS_HPP
c4science · Help