Remaillage_FT.h

/****************************************************************************
* Copyright (c) 2022, 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 Remaillage_FT_included
#define Remaillage_FT_included
 
#include <TRUSTTabFT_forward.h>
#include <Objet_U.h>
#include <FTd_tools.h>
#include <TRUST_Ref.h>
 
class Domaine_dis_base;
class Maillage_FT_Disc;
class Domaine_VF;
class Param;
 
// POUR DEBUGGER LA CONSERVATION DU VOLUME
#define DEBUG_CONSERV_VOLUME 0
 
/*! @brief : class Remaillage_FT Cette classe implemente les procedures de remaillage des interfaces pour le Front-Tracking :
 *
 *
 *
 * @sa Transport_Interfaces_FT_Disc Maillage_FT_Disc
 */
class Remaillage_FT : public Objet_U
{
  Declare_instanciable_sans_constructeur(Remaillage_FT);
public:
  Remaillage_FT();
  void set_param(Param& p) const override;
  void associer_domaine(const Domaine_dis_base& domaine_dis);
 
  int supprimer_facettes_nulles(Maillage_FT_Disc& maillage) const;
  int traite_adherence(Maillage_FT_Disc& maillage) const;
  int traite_decollement(Maillage_FT_Disc& maillage, const DoubleTab& deplacement) const;
 
  int a_remailler(double temps, const Maillage_FT_Disc& maillage) const;
  int a_lisser(double temps) const;
 
  // Remaillage local avec changement des connectivites:
  void remaillage_local_interface(double temps, Maillage_FT_Disc& maillage);
 
  // Regularisations du maillage (deplacement des sommets sans changer la connectivite)
  void corriger_volume(Maillage_FT_Disc& maillage, ArrOfDouble& var_volume);
  void corriger_volume_(Maillage_FT_Disc& maillage, ArrOfDouble& var_volume, const int nb_iter_corrections_vol);
  void barycentrer_lisser_systematique(double temps, Maillage_FT_Disc& maillage);
  void barycentrer_lisser_apres_remaillage(Maillage_FT_Disc& maillage, ArrOfDouble& var_volume);
 
  int  sauvegarder(Sortie& ) const override;
  int  reprendre(Entree&) override;
 
  double calculer_variation_volume(const Maillage_FT_Disc& maillage,
                                   const DoubleTab& position_initiale,
                                   ArrOfDouble& varVolume) const;
 
  void lisser_dvolume(const Maillage_FT_Disc& maillage,
                      ArrOfDouble& var_volume,
                      const int nb_iterations) const;
 
  virtual void regulariser_courbure(Maillage_FT_Disc& maillage,
                                    const double coeff,
                                    ArrOfDouble& dvolume) const;
 
  int get_nb_iter_bary_volume_seul()
  {
    return nb_iter_bary_volume_seul_;
  };
 
  void set_is_solid_particle(const bool is_solid_particle) {is_solid_particle_=is_solid_particle;}
  const bool& get_is_solid_particle() const {return is_solid_particle_;}
 
#if DEBUG_CONSERV_VOLUME
  double calculer_volume_mesh(const Maillage_FT_Disc& mesh) const;
  double calculer_somme_dvolume(const Maillage_FT_Disc&, const ArrOfDouble&) const;
#endif
 
 
protected:
 
  int tester_a_remailler(const Maillage_FT_Disc& maillage) const;
 
  double regulariser_maillage(Maillage_FT_Disc& maillage,
                              ArrOfDouble& var_volume,
                              const double facteur_barycentrage_tangent,
                              const double coeff_lissage,
                              const int nb_iter_barycentrage,
                              const int nb_iter_lissage,
                              const int max_nb_iter_correction_volume,
                              const double seuil_dvolume) const;
 
  double redistribuer_sommets(Maillage_FT_Disc&   maillage,
                              const double        relaxation_direction_tangente,
                              const double        relaxation_direction_normale,
                              ArrOfDouble& var_volume_impose,
                              ArrOfDouble& var_volume_obtenu) const;
 
  int calculer_barycentre_facettes_voisines(const Maillage_FT_Disc& maillage,
                                            DoubleTab& barycentres) const;
 
  int calculer_connectivites_sommetFacettes(const Maillage_FT_Disc& maillage,ArrOfInt& fa7VoisinesSom_index, IntTab& fa7VoisinesSom_data) const;
  // Used by IJK only:
  int calculer_correction_deplacement(DoubleTab& deplacement,const ArrOfDouble& varVolume,const DoubleTab& deplacement_varVolume, const ArrOfDouble& norme2_deplacement_varVolume) const;
  int calculer_differentielle_volume(const Maillage_FT_Disc& maillage,
                                     DoubleTab& differentielle_volume) const;
 
  double calculer_variation_volume_facette_2D(int fa7, const Maillage_FT_Disc& maillage,const DoubleTab& position_initiale) const;
  double calculer_variation_volume_facette_3D(int fa7, const Maillage_FT_Disc& maillage,const DoubleTab& position_initiale) const;
 
  virtual double calculer_longueurIdeale2_arete(const Maillage_FT_Disc& maillage, int som0, double x, double y, double z) const;
  //  double calculer_longueurIdeale2_arete(const Maillage_FT_Disc& maillage, int som0, int som1) const;
  int supprimer_petites_aretes(Maillage_FT_Disc& maillage, ArrOfDouble& varVolume) const;
  int diviser_grandes_aretes(Maillage_FT_Disc& maillage) const;
  //int marquer_sommets_petites_aretes(Maillage_FT_Disc& maillage, ArrOfInt & tab_somSupp) const;
  int marquer_aretes(Maillage_FT_Disc& maillage, IntTab& tab_aretesMarquees, ArrOfInt& tab_somD, DoubleTab& tab_deplacement_somD, int drap) const;
  int inserer_tab_aretes(int& nb_tab_aretes,IntTab& tab_aretes,DoubleTab& tab_criteres,int pe0, int numOwner0, int pe1, int numOwner1, int face_bord1,int peRequete, int fa7_peR, int iarete_fa7_peR) const;
  int chercher_arete_tab(int tmp, const ArrOfInt& tab_index, const IntTab& tab_aretes,int pe0, int numOwner0, int pe1, int numOwner1) const;
  double calculer_volume_sommets_supprimes(const Maillage_FT_Disc& maillage, const ArrOfInt& tab_somSupp,ArrOfDouble& varVolume) const;
 
  int supprimer_facettes_bord(Maillage_FT_Disc& maillage) const;
  int supprimer_doublons_facettes(Maillage_FT_Disc& maillage) const;
  int permuter_aretes(Maillage_FT_Disc& maillage) const;
  double qualiteTriangle(const FTd_vecteur3& som0, const FTd_vecteur3& som1, const FTd_vecteur3& som2, double& aire) const;
 
  int nettoyer_maillage(Maillage_FT_Disc& maillage) const;
 
  OBS_PTR(Domaine_VF) refdomaine_VF_;
 
  // Values initialized:
  double temps_ = 0;
  double temps_dernier_remaillage_ = -1.e40;
  double temps_dernier_lissage_ = -1.e40;
 
  // Values initialized but which can be set in the data file:
  double dt_remaillage_ = -1.;
  double dt_lissage_ = -1.;
  int nb_iter_remaillage_ = 1;
  int nb_iter_barycentrage_ = 1;
  int nb_iter_bary_volume_seul_ = 3;
  double seuil_dvolume_residuel_ = 0.;
  double relax_barycentrage_ = 1.;
 
  // Criteres de longueur ideale des aretes
  double critere_arete_ = 0.35;
  int impr_ = -1;
  // L'un de ces deux vaut -1, l'autre doit etre positif:
  double valeur_longueur_fixe_ = -1.;
  double facteur_longueur_ideale_ = -1.;
  int equilateral_ = 0;
  // L'orientation de la facette pour ajuster sa taille a l'element eulerien.
  // alors que equilateral=1 utilise la diagonal de l'element comme longueur de reference.
  // Avec equilateral_ = 0, les facettes sont etirees comme le maillage et facteur_longueur vaut
  // 1. si l'arete traverse tout l'element dans la direction donnee par l'arrete.
 
  double variation_volume_ = 0.;
  double surface_interface_ = 0.;
 
 
  // Coefficient de "diffusion" (en quelque sorte la CFL du schema
  // de lissage).
  double lissage_courbure_coeff_ = -0.05; // valeur typique pour stabilite
  // Nombre d'iterations de lissage a faire lors des operations de lissage (dt_lissage_)
  int lissage_courbure_iterations_systematique_ = 0;
  // Nombre d'iterations de lissage a faire en cas de remaillage local ou global
  //  (declanche uniquement si le maillage est effectivement modifie)
  int lissage_courbure_iterations_si_remaillage_ = 0;
  // Ancien parametre de lissage (voir readOn)
  int lissage_courbure_iterations_old_ = -1;
  // Critere local de declenchement du lissage:
  double lissage_critere_ = 0.; // Default value to 0, when lissage is applied, it is for the whole mesh
 
  bool is_solid_particle_=false;  // for fpi module, pointer to NS_FT_Disc::is_solid_particle_
};
 
 
#endif