Fluide_generique_TPPI_base.h

/****************************************************************************
* Copyright (c) 2024, CEA
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
* IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
* OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*****************************************************************************/
 
#ifndef Fluide_generique_TPPI_base_included
#define Fluide_generique_TPPI_base_included
 
#include <Fluide_reel_base.h>
#include <TPPI.h>
 
class Fluide_generique_TPPI_base : public Fluide_reel_base
{
  Declare_base( Fluide_generique_TPPI_base ) ;
public :
 
  MRange unknown_range() const override
  {
    if (tmax_ < -100. )
      return Fluide_reel_base::unknown_range();
 
    return { { "temperature", { tmin_ - 273.15, tmax_ - 273.15 } }, { "pression", { pmin_, pmax_ } } };
  }
 
  MRange unknown_range_h() const override
  {
    if (hmax_ < -100. )
      return Fluide_reel_base::unknown_range_h();
 
    return { { "enthalpie", { hmin_ , hmax_ } }, { "pression", { pmin_, pmax_ } } };
  }
 
  // lois en T
  void rho_(const SpanD T, const SpanD P, SpanD R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_pT(P, Tk_(T), R, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void dP_rho_(const SpanD T, const SpanD P, SpanD dP_R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_dp_pT(P, Tk_(T), dP_R, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void dT_rho_(const SpanD T, const SpanD P, SpanD dT_R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_dT_pT(P, Tk_(T), dT_R, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void h_(const SpanD T, const SpanD P, SpanD H, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_h_pT(P, Tk_(T), H, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void dP_h_(const SpanD T, const SpanD P, SpanD dP_H, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_h_dp_pT(P, Tk_(T), dP_H, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void dT_h_(const SpanD T, const SpanD P, SpanD dT_H, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_h_dT_pT(P, Tk_(T), dT_H, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void cp_(const SpanD T, const SpanD P, SpanD CP, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_cp_pT(P, Tk_(T), CP, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void beta_(const SpanD T, const SpanD P, SpanD B, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_beta_pT(P, Tk_(T), B, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void mu_(const SpanD T, const SpanD P, SpanD M, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_mu_pT(P, Tk_(T), M, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void lambda_(const SpanD T, const SpanD P, SpanD L, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_lambda_pT(P, Tk_(T), L, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void compute_CPMLB_pb_multiphase_(const MSpanD input, MLoiSpanD prop, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_CPMLB_pb_multiphase_pT(input, prop, ncomp, id);
  }
 
  void compute_all_pb_multiphase_(const MSpanD input, MLoiSpanD inter, MLoiSpanD bord, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_all_pb_multiphase_pT(input,inter, bord, ncomp, id);
  }
 
  // lois en h
  void rho_h_(const SpanD h, const SpanD P, SpanD R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_ph(P, h, R, ncomp, id);
  }
 
  void dP_rho_h_(const SpanD h, const SpanD P, SpanD dP_R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_dp_ph(P, h, dP_R, ncomp, id);
  }
 
  void dh_rho_h_(const SpanD h, const SpanD P, SpanD dT_R, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_rho_dh_ph(P, h, dT_R, ncomp, id);
  }
 
  void T_(const SpanD h, const SpanD P, SpanD T, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_T_ph(P, h, T, ncomp, id);
    Tc_(T); /* put back T in C */
  }
 
  void dP_T_(const SpanD h, const SpanD P, SpanD dP_T, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_T_dp_ph(P, h, dP_T, ncomp, id);
  }
 
  void dh_T_(const SpanD h, const SpanD P, SpanD dh_T, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_T_dh_ph(P, h, dh_T, ncomp, id);
  }
 
  void cp_h_(const SpanD h, const SpanD P, SpanD CP, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_cp_ph(P, h, CP, ncomp, id);
  }
 
  void beta_h_(const SpanD h, const SpanD P, SpanD B, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_beta_ph(P, h, B, ncomp, id);
  }
 
  void mu_h_(const SpanD h, const SpanD P, SpanD M, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_mu_ph(P, h, M, ncomp, id);
  }
 
  void lambda_h_(const SpanD h, const SpanD P, SpanD L, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_lambda_ph(P, h, L, ncomp, id);
  }
 
  void compute_CPMLB_pb_multiphase_h_(const MSpanD input, MLoiSpanD_h prop, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_CPMLB_pb_multiphase_ph(input, prop, ncomp, id);
  }
 
  void compute_all_pb_multiphase_h_(const MSpanD input, MLoiSpanD_h inter, MLoiSpanD_h bord, int ncomp = 1, int id = 0) const override
  {
    TPPI_->tppi_get_all_pb_multiphase_ph(input, inter, bord, ncomp, id);
  }
 
protected:
  std::shared_ptr<TPPI> TPPI_ = nullptr;
  Motcle model_name_, fluid_name_;
  double tmin_ = -123., tmax_ = -123., pmin_ = -123., pmax_ = -123., hmin_ = -123., hmax_ = -123.;
};
 
#endif /* Fluide_generique_TPPI_base_included */