safe

transport_parameters.hpp
View Options

			#ifndef SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERS_HPP
			#define SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERS_HPP

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


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

			SaturatedDiffusionTransportParameters(double diffusion_coefficient, double porosity):
			m_diff(diffusion_coefficient), m_poro(porosity)
			{}

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

			//! \brief Diffusion coefficient (element by element) (m^2/s)
			double diffusion_coefficient(int element) {return m_diff;}

			//! \brief Porosity (at a node (function of composition))
			double porosity(int node) {return m_poro;}

			private:
			double m_diff;
			double m_poro;
			};

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

			#endif // SPECMICP_REACTMICP_SATURATEDDIFFUSION_TRANSPORTPARAMETERS_HPP