Navier_Stokes_FT_Disc_interne.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 Navier_Stokes_FT_Disc_interne_included
#define Navier_Stokes_FT_Disc_interne_included
 
 
class Navier_Stokes_FT_Disc_interne
{
public:
  Navier_Stokes_FT_Disc_interne() :
    terme_gravite_(GRAVITE_GRAD_I),   // Par defaut terme gravite ft sans courants parasites
    type_interpol_indic_pour_dI_dt_(INTERP_STANDARD), // Default is the historical interpolation
    OutletCorrection_pour_dI_dt_(NO_CORRECTION),   // Default is the historical
    is_penalized(0),                  // Par defaut, pas de penalisation L2 du forcage
    eta(1.0),                         // Par defaut, coefficient de penalisation L2 = 1.
    p_ref_pena(-1.e40),               // Par defaut, pas de penalisation L2 de la pression sinon valeur reference
    is_pfl_flottant(0),               // Traitement local Dirichlet pression si les CL pression sont toutes en Neumann homogene
    x_pfl_imp(-1.e40),                // Par defaut, x, y, z du point de modification de la pression fluide
    y_pfl_imp(-1.e40), z_pfl_imp(-1.e40)
  { }
 
  // Par defaut, pas de prise en compte de la courbure pour corriger le champ etendu delta_vitesse
  int correction_courbure_ordre_ = 0;
 
  // Par defaut, mpoint cree un saut de vitesse
  bool mpoint_inactif = false;
  // Par defaut, mpointv  cree un saut de vitesse
  bool mpointv_inactif = false;
 
  OWN_PTR(Champ_Fonc_base)  second_membre_projection;
  OWN_PTR(Champ_Fonc_base)  second_membre_projection_jump_;
  OWN_PTR(Champ_Fonc_base)  derivee_u_etoile;
  OWN_PTR(Champ_Fonc_base)  gradient_pression;
  OWN_PTR(Champ_Fonc_base)  terme_diffusion;
  OWN_PTR(Champ_Fonc_base)  terme_convection;
  OWN_PTR(Champ_Fonc_base)  terme_source;
  OWN_PTR(Champ_Fonc_base)  terme_source_interfaces;
  OWN_PTR(Champ_Fonc_base)  indicatrice_p1b;
  OWN_PTR(Champ_Fonc_base)  gradient_indicatrice;
  OWN_PTR(Champ_Fonc_base)  potentiel_faces;
  OWN_PTR(Champ_Fonc_base)  potentiel_elements;
  // delta_u_interface = la partie "saut de vitesse" du champ de vitesse a l'interface
  OWN_PTR(Champ_Inc_base) delta_u_interface;
  OWN_PTR(Champ_Fonc_base)  laplacien_d;
  OWN_PTR(Champ_Fonc_base)  mpoint;
  OWN_PTR(Champ_Fonc_base)  mpoint_vap;
  // Variation temporelle indicatrice de phase
  OWN_PTR(Champ_Fonc_base)  derivee_temporelle_indicatrice;
  OWN_PTR(Champ_Fonc_base)  ai; // Eulerian interfacial area.
  OWN_PTR(Champ_Inc_base) vitesse_jump0_; // Extended Velocity of phase 0.
  OWN_PTR(Champ_Inc_base) vitesse_jump1_; // Extended Velocity of phase 1.
 
  OWN_PTR(Champ_Fonc_base) contact_force_source_term; // for fpi module, solid-solid collision source term
 
  LIST(OBS_PTR(Champ_base)) liste_champs_compris;
 
  // Si matrice_pression_invariante != 0,
  //   on ne recalcule pas la matrice de pression a chaque pas de temps.
  // Par defaut, recalculer la matrice pression
  bool matrice_pression_invariante = false;
 
  // Si on veut ajouter une interface a vitesse imposee :
  //  reference a l'equation de transport correspondante :
  VECT(OBS_PTR(Transport_Interfaces_FT_Disc)) ref_eq_interf_vitesse_imposee;
 
  // Si le fluide est diphasique, c'est l'indicatrice de l'equation suivante
  // qui est utilisee pour determiner les proprietes du fluide:
  // (masse volumique, viscosite, tension superficielle, ...)
  OBS_PTR(Transport_Interfaces_FT_Disc) ref_eq_interf_proprietes_fluide;
 
  // Si le fluide est diphasique, la reference au fluide:
  OBS_PTR(Fluide_Diphasique) ref_fluide_diphasique;
 
  OBS_PTR(Convection_Diffusion_Temperature_FT_Disc) ref_equation_mpoint_;
  OBS_PTR(Convection_Diffusion_Temperature_FT_Disc) ref_equation_mpoint_vap_;
 
  // Valeur maximale de courbure autorisee pour calculer le
  // terme source de tension de surface (clipping si valeur superieur)
  // Par defaut, pas de clipping
  double clipping_courbure_interface = 1e40;
 
  enum Terme_Gravite
  {
    GRAVITE_RHO_G, GRAVITE_GRAD_I
  };
  Terme_Gravite terme_gravite_;
  Noms equations_concentration_source_fluide_;
 
  // Si is_explicite != 0,
  //   on calcul vpoint de facon explicite dans l etape de prediction des vitesses.
  // Par defaut, calcul explicite de vpoint etape predicition
  bool is_explicite = true;
 
  // Si is_boussinesq_ != 0, on calcul une force par le modele de boussinesq
  // Par defaut, l'hypothese de Boussinesq n'est pas utilisee pour la flottabilite dans les phases.
  bool is_boussinesq_ = false;
 
  // Flag pour la method du saut de vitesse :
  // Par defaut, on utilise la methode historique pour imposer le saut de vitesse du changement de phase.
  bool new_mass_source_ = false;
 
  // Flag if to shift second member 2 due to phase change:
  // Par defaut, secmem2 is located at interface cells
  bool shift_secmem2_ = false;
 
  enum Type_interpol_indic_pour_dI_dt
  {
    INTERP_STANDARD, INTERP_MODIFIEE, INTERP_AI_BASED, INTERP_STANDARD_UVEXT, INTERP_MODIFIEE_UVEXT, INTERP_AI_BASED_UVEXT, INTERP_STANDARD_UIEXT, INTERP_MODIFIEE_UIEXT, INTERP_AI_BASED_UIEXT
  };
  Type_interpol_indic_pour_dI_dt type_interpol_indic_pour_dI_dt_;
 
  enum OutletCorrection_pour_dI_dt
  {
    NO_CORRECTION, CORRECTION_GHOST_INDIC, ZERO_NET_FLUX_ON_MIXED_CELLS, ZERO_OUT_FLUX_ON_MIXED_CELLS
  };
  OutletCorrection_pour_dI_dt OutletCorrection_pour_dI_dt_;
 
  // Si is_penalized != 0,
  //   on penalise L2 le terme de forcage.
  int is_penalized;
 
  // Valeur de l'inverse du coefficient de penalisation L2 du terme de forcage.
  double eta;
 
  // Valeur pour la penalisation L2 de la pression.
  double p_ref_pena;
 
  // Point de penalisation L2 de la pression du fluide
  int is_pfl_flottant; // Traitement local Dirichlet pression si les CL pression sont toutes en Neumann homogene
  double x_pfl_imp;
  double y_pfl_imp;
  double z_pfl_imp;
 
};
 
#endif /* Navier_Stokes_FT_Disc_interne_included */