Modele_turbulence_hyd_K_Omega.h

/****************************************************************************
* Copyright (c) 2023, 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 Modele_turbulence_hyd_K_Omega_included
#define Modele_turbulence_hyd_K_Omega_included
 
#include <TRUSTTabs_forward.h>
#include <Transport_K_Omega.h>
#include <Modele_turbulence_hyd_RANS_K_Omega_base.h>
 
/*! @brief Classe Modele_turbulence_hyd_K_Omega Cette classe represente le modele de turbulence (k, omega) pour les
 *
 *     equations de Navier-Stokes.
 *
 * @sa Modele_turbulence_hyd_base Modele_turbulence_hyd_LES_base
 */
class Modele_turbulence_hyd_K_Omega: public Modele_turbulence_hyd_RANS_K_Omega_base, public Modele_turbulence_hyd_RANS_Gen<Modele_turbulence_hyd_K_Omega>
{
  Declare_instanciable(Modele_turbulence_hyd_K_Omega);
public:
  void set_param(Param& param) const override;
  int preparer_calcul() override;
  void completer() override;
  bool initTimeStep(double dt) override;
  void mettre_a_jour(double) override;
  virtual inline Champ_Inc_base& get_set_eq_K_Omega();
  virtual inline const Champ_Inc_base& get_eq_K_Omega() const;
  inline Transport_K_Omega_base& get_set_eq_transport() override;
  inline const Transport_K_Omega_base& get_eq_transport() const override;
  void get_noms_champs_postraitables(Noms& nom,Option opt=NONE) const override;
  void discretiser() override;
 
  inline const DoubleTab& get_tabF1() const;
  inline DoubleTab& get_tabF1();
  inline const DoubleTab& get_tabF2() const;
  inline DoubleTab& get_tabF2();
  inline const DoubleTab& get_enstrophy() const;
  inline DoubleTab& get_enstrophy();
  inline const DoubleTab& get_strain_invariant() const;
  inline DoubleTab& get_strain_invariant();
 
  void controler() override { get_set_eq_transport().controler_K_Omega(); }
  virtual Champ_Fonc_base& calculer_viscosite_turbulente(double );
 
  inline bool is_SST() const { return is_SST_ ;}
  inline bool is_SST_or_BSL() const { return is_BSL_ ;} // if SST is a specialization of BSL
 
protected:
  void fill_turbulent_viscosity_tab(const DoubleTab& K_Omega, DoubleTab& turbulent_viscosity);
  OWN_PTR(Champ_Fonc_base) tabF1_; // Blending field for SST model
  OWN_PTR(Champ_Fonc_base) tabF2_; // for the turbulent viscosity in the SST model
  OWN_PTR(Champ_Fonc_base) enstrophy_; // for the turbulent viscosity in the SST model
  OWN_PTR(Champ_Fonc_base) strain_invariant_; // for the turbulent viscosity in the SST model
  OWN_PTR(Champ_Inc_base) visco_turb_au_format_K_Omega_;
  bool is_SST_ = false; // check if model variant is SST
  bool is_BSL_ = false; // check if model variant is BSL
};
 
/*! @brief Renvoie le champ inconnue du modele de turbulence i.
 *
 * e. : (K, Omega). Cette inconnue est portee
 *     par l'equation de transport K-Omega porte par le modele.
 *     (version const)
 *
 * @return (Champ_Inc_base&) le champ inconnue (K, Omega)
 */
inline const Champ_Inc_base& Modele_turbulence_hyd_K_Omega::get_eq_K_Omega() const
{
  const Champ_Inc_base& inc=get_eq_transport().inconnue();
  return inc;
}
 
/*! @brief Renvoie le champ inconnue du modele de turbulence i.
 *
 * e. : (K, Omega). Cette inconnue est portee
 *     par l'equation de transport K-Omega porte par le modele.
 *
 * @return (Champ_Inc_base&) le champ inconnue (K, Omega)
 */
inline Champ_Inc_base& Modele_turbulence_hyd_K_Omega::get_set_eq_K_Omega()
{
  return get_set_eq_transport().inconnue();
}
 
inline const Transport_K_Omega_base& Modele_turbulence_hyd_K_Omega::get_eq_transport() const
{
  const Transport_K_Omega_base& eq_transport=ref_cast(Transport_K_Omega_base,ptr_eq_transport_.valeur());
  return eq_transport;
}
 
inline Transport_K_Omega_base& Modele_turbulence_hyd_K_Omega::get_set_eq_transport()
{
  Transport_K_Omega_base& eq_transport=ref_cast(Transport_K_Omega_base,ptr_eq_transport_.valeur());
  return eq_transport;
}
 
/*! @brief Returns the blending table F1 for SST.
 *
 *     (version const)
 *
 * @return (DoubleTab&) tabF1_
 */
inline const DoubleTab& Modele_turbulence_hyd_K_Omega::get_tabF1() const
{
  return tabF1_->valeurs();
}
 
/*! @brief Returns the blending table F1 for SST.
 *
 * @return (DoubleTab&) tabF1_
 */
inline DoubleTab& Modele_turbulence_hyd_K_Omega::get_tabF1()
{
  return tabF1_->valeurs();
}
 
/*! @brief Returns the table F2 for SST.
 *
 *     (version const)
 *
 * @return (DoubleTab&) tabF2_
 */
inline const DoubleTab& Modele_turbulence_hyd_K_Omega::get_tabF2() const
{
  return tabF2_->valeurs();
}
 
/*! @brief Returns the table F2 for SST.
 *
 * @return (DoubleTab&) tabF2_
 */
inline DoubleTab& Modele_turbulence_hyd_K_Omega::get_tabF2()
{
  return tabF2_->valeurs();
}
 
/*! @brief Returns the field enstrophy.
 *
 *     (version const)
 *
 * @return (DoubleTab&) enstrophy
 */
inline const DoubleTab& Modele_turbulence_hyd_K_Omega::get_enstrophy() const
{
  return enstrophy_->valeurs();
}
 
/*! @brief Returns the field F2.
 *
 * @return (DoubleTab&) tabF2
 */
inline DoubleTab& Modele_turbulence_hyd_K_Omega::get_enstrophy()
{
  return enstrophy_->valeurs();
}
 
/*! @brief Returns the field strain_invariant.
 *
 *     (version const)
 *
 * @return (DoubleTab&) enstrophy
 */
inline const DoubleTab& Modele_turbulence_hyd_K_Omega::get_strain_invariant() const
{
  return strain_invariant_->valeurs();
}
 
/*! @brief Returns the field F2.
 *
 * @return (DoubleTab&) tabF2
 */
inline DoubleTab& Modele_turbulence_hyd_K_Omega::get_strain_invariant()
{
  return strain_invariant_->valeurs();
}
#endif /* Modele_turbulence_hyd_K_Omega_included */