impl_mass_transfer.hpp

#ifndef __IMPL_HYDRO_MASS_TRANSFER_HPP__
#define __IMPL_HYDRO_MASS_TRANSFER_HPP__
 
#ifndef NDEBUG
#  pragma GCC diagnostic push
#  pragma GCC diagnostic ignored "-Wunused-but-set-variable"
#  pragma GCC diagnostic ignored "-Wnan-infinity-disabled"
#endif
#include <Eigen/Core>
#include <Eigen/Dense>
#ifndef NDEBUG
#  pragma GCC diagnostic pop
#endif
#include <cma_utils/alias.hpp>
#include <eigen_kokkos.hpp>
#include <mc/domain.hpp>
#include <simulation/mass_transfer.hpp>
 
namespace Simulation
{
  class ScalarSimulation;
};
 
namespace Simulation::MassTransfer
{
  struct MassTransferProxy
  {
    ColMajorMatrixtype<double> mtr;
    Eigen::ArrayXXd kla;
    Eigen::ArrayXd Henry;
    double db;
  };
 
  namespace Impl
  {
 
    void flowmap_gas_liquid_mass_transfer(
        MassTransferProxy& mtr,
        const Eigen::ArrayXXd& liquid_concentration,
        const Eigen::ArrayXXd& gas_concentration,
        const Eigen::MatrixXd& liquid_volume,
        const CmaUtils::IterationStatePtrType& state);
 
    void fixed_kla_gas_liquid_mass_transfer(
        MassTransferProxy& mtr,
        const Eigen::ArrayXXd& liquid_concentration,
        const Eigen::ArrayXXd& gas_concentration,
        const Eigen::MatrixXd& liquid_volume,
        const CmaUtils::IterationStatePtrType& state);
 
    void flowmap_kla_gas_liquid_mass_transfer(
        MassTransferProxy& mtr,
        const Eigen::ArrayXXd& liquid_concentration,
        const Eigen::ArrayXXd& gas_concentration,
        const Eigen::MatrixXd& liquid_volume,
        const CmaUtils::IterationStatePtrType& state);
 
  }; // namespace Impl
 
  struct KlFonctor
  {
    double operator()(double kl);
    Eigen::ArrayXXd operator()(const Eigen::ArrayXd& energy_dissipation_array,
                               const Eigen::ArrayXd& diffusion);
    double operator()(double kla, bool _a);
  };
 
  struct InterfacialAreaFonctor
  {
    double
    operator()()
    {
      return 1;
    }
    double operator()(double a);
    Eigen::ArrayXXd operator()(const Eigen::ArrayXd& a);
    Eigen::ArrayXXd operator()(double diameter, const Eigen::ArrayXd& a);
  };
 
  struct EquilibriumFunctor
  {
    double operator()(double cstar);
    Eigen::ArrayXXd operator()(const Eigen::ArrayXXd& cstar);
    Eigen::ArrayXXd operator()(const Eigen::ArrayXXd& henry, bool _henry);
  };
 
} // namespace Simulation::MassTransfer
 
#endif //__SIMULATION_HYDRO_MASS_TRANSFER_HPP__ __SIM_HYDRO_MASS_TRANSFER__