transport_parameters_base.hpp
745d2292fc54reactivetransport

transport_parameters_base.hpp
View Options

			#ifndef SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERSBASE_HPP
			#define SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERSBASE_HPP

			#include "common.hpp"

			namespace specmicp {
			namespace reactmicp {
			namespace systems {
			namespace siasaturated {


			//! \brief Parameters for the diffusion system
			class SaturatedDiffusionTransportParametersBase
			{
			public:

			SaturatedDiffusionTransportParametersBase(uindex_t nbnodes):
			m_porosities(nbnodes)
			{}

			SaturatedDiffusionTransportParametersBase(
			index_t nbnodes,
			scalar_t the_porosity
			):
			SaturatedDiffusionTransportParametersBase(nbnodes)
			{
			m_porosities.setConstant(the_porosity);
			}

			//! \brief Density of water (kg/m^3)
			scalar_t density_water() {return 1e3;}

			//! \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_porosities;
			};

			} // end namespace siasaturated
			} // end namespace systems
			} // end namespace reactmicp
			} // end namespace specmicp

			#endif // SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERSBASE_HPP