Remaillage_FT

Remaillage_FT#

#include <Remaillage_FT.h>

: class Remaillage_FT Cette classe implemente les procedures de remaillage des interfaces pour le Front - Tracking :

How to cite this class in this doc#

RST

:ref:`Remaillage\_FT <class-remaillage-ft>`

Markdown

[Remaillage\_FT](#class-remaillage-ft)

Detailed description#

Inherits from#

public : Objet_U

Inherited by#

public : Remaillage_FT_IJK

List of Public Static Methods#

cree_instance
info
self_cast
self_cast

Complete Member Function Documentation#

a_lisser

int Remaillage_FT::a_lisser(double temps) const

a_remailler

int Remaillage_FT::a_remailler(double temps, const Maillage_FT_Disc &maillage) const

associer_domaine

void Remaillage_FT::associer_domaine(const Domaine_dis_base &domaine_dis)

Cette fonction stocke le domaine_dis dans refdomaine_dis_.

barycentrer_lisser_apres_remaillage

void Remaillage_FT::barycentrer_lisser_apres_remaillage(Maillage_FT_Disc &maillage, ArrOfDouble &var_volume)

idem mais avec le nombre d’iterations de lissage si remaillage

Precondition: pas d’elements virtuels (on doit avoir appele nettoyer_elements_virtuels())

barycentrer_lisser_systematique

void Remaillage_FT::barycentrer_lisser_systematique(double temps, Maillage_FT_Disc &maillage)

applique barycentrage, lissage et correction de volume.

On applique le nombre d’iterations de lissage systematique. Precondition: pas de facettes virtuelles

calculer_variation_volume

double Remaillage_FT::calculer_variation_volume(const Maillage_FT_Disc &maillage, const DoubleTab &position_initiale, ArrOfDouble &varVolume) const

Cette fonction calcule le volume de phase 0 engendre par chaque sommet lors du deplacement de l’interface entre la position initiale et la position actuelle.

du maillage (le volume engendre est positif si l’interface se deplace dans la direction de la normale aux facettes). On definit d’abord le volume engendre par le deplacement de chaque facette : C’est le volume d’un polyedre construit comme reunion de trois tetraedres, engendres par deplacement successifs des trois sommets dans un ordre conventionnel (deux sommets d’une arete sont toujours deplaces dans le meme ordre pour que les deux triangles voisins engendrent des tetraedres conformes, sans laisser de trous et sans se chevaucher) Ensuite, le volume engendre par chaque facette est divise en trois parts egales et reparti sur les trois sommets de la facette.

corriger_volume

void Remaillage_FT::corriger_volume(Maillage_FT_Disc &maillage, ArrOfDouble &var_volume)

deplacement des sommets se sorte a produire la variation de volume prescrite a chaque sommet.

Precondition: pas de facettes virtuelles

corriger_volume_

void Remaillage_FT::corriger_volume_(Maillage_FT_Disc &maillage, ArrOfDouble &var_volume, const int nb_iter_corrections_vol)

duplique

int Remaillage_FT::duplique() const override

Reimplements:

duplique

get_info

const Type_info * Remaillage_FT::get_info() const override

Donne des informations sur le type de l’ Objet_U .

Reimplements:

get_info

get_is_solid_particle

const bool & Remaillage_FT::get_is_solid_particle() const

get_nb_iter_bary_volume_seul

int Remaillage_FT::get_nb_iter_bary_volume_seul()

lisser_dvolume

void Remaillage_FT::lisser_dvolume(const Maillage_FT_Disc &maillage, ArrOfDouble &var_volume, const int nb_iterations) const

Regularise le champ scalaire “var_volume” defini aux sommets du “maillage”.

On regularise en faisant calculant une valeur moyenne par facette, puis en recalculant une moyenne aux sommets en fonction de la moyenne aux faces. Les moyennes sont ponderees par la surface des facettes… Le lissage est repete un nombre donne de fois. Cette methode est utilisee pour lisser les corrections de volume et eviter l’apparition de pointes sur les interfaces.

regulariser_courbure

void Remaillage_FT::regulariser_courbure(Maillage_FT_Disc &maillage, const double coeff, ArrOfDouble &dvolume) const

Regularise le maillage en deplacant les sommets pour reduire les gradients de courbure.

L’equation d’evolution du maillage est une diffusion de masse le long de l’interface proportionnelle au gradient de courbure.

Remaillage_FT

Remaillage_FT::Remaillage_FT()

remaillage_local_interface

void Remaillage_FT::remaillage_local_interface(double temps, Maillage_FT_Disc &maillage)

Verifie les criteres de remaillage locaux (longueur des aretes) et effectue les remaillages locaux si necessaires.

reprendre

int Remaillage_FT::reprendre(Entree &) override

Reprise d’un Objet_U sur un flot d’entree Methode a surcharger.

Reimplements:

reprendre

sauvegarder

int Remaillage_FT::sauvegarder(Sortie &) const override

Sauvegarde d’un Objet_U sur un flot de sortie Methode a surcharger.

Reimplements:

sauvegarder

set_is_solid_particle

void Remaillage_FT::set_is_solid_particle(const bool is_solid_particle)

set_param

void Remaillage_FT::set_param(Param &p)

Methode appelee par readOn.

Declaration des parametres a lire dans le .data

supprimer_facettes_nulles

int Remaillage_FT::supprimer_facettes_nulles(Maillage_FT_Disc &maillage) const

Cette fonction “supprime” les facettes de surface nulle : Elle les reduit a 1 sommet (le sommet 0, pour ne pas changer de processeur)

taille_memoire

unsigned Remaillage_FT::taille_memoire() const override

Reimplements:

taille_memoire

traite_adherence

int Remaillage_FT::traite_adherence(Maillage_FT_Disc &maillage) const

Cette fonction permet de gerer l’adherence de l’interface a la paroi Si une facette possede 3 sommets sur un paroi, elle est supprimee.

traite_decollement

int Remaillage_FT::traite_decollement(Maillage_FT_Disc &maillage, const DoubleTab &deplacement) const

Cette fonction permet de gerer le decollement de l’interface de la paroi Si un sommet de bord n’a pas de sommet voisin de bord.

et que son deplacement tend a le faire rentrer dans le domaine, il est marque comme sommet interne

~Remaillage_FT

Remaillage_FT::~Remaillage_FT()

cree_instance

Objet_U * Remaillage_FT::cree_instance()

info

const Type_info * Remaillage_FT::info()

self_cast

const Remaillage_FT & Remaillage_FT::self_cast(const Objet_U &)

self_cast

Remaillage_FT & Remaillage_FT::self_cast(Objet_U &)

calculer_barycentre_facettes_voisines

int Remaillage_FT::calculer_barycentre_facettes_voisines(const Maillage_FT_Disc &maillage, DoubleTab &barycentres) const

Cette fonction calcule pour chaque sommet le barycentre de l’ensemble des facettes voisines du sommet.

Les “dimension” premieres colonnes contiennent les coordonnees du barycentre, la colonne “dimension” contient la somme des surfaces des facettes voisines du sommet.

calculer_connectivites_sommetFacettes

int Remaillage_FT::calculer_connectivites_sommetFacettes(const Maillage_FT_Disc &maillage, ArrOfInt &fa7VoisinesSom_index, IntTab &fa7VoisinesSom_data) const

Cette fonction calcule les connectivites sommet - >facettes voisines Les facettes voisines des sommets sont stockees dans le tableau fa7VoisinesSom_data.

sous forme de liste chainee. Pour le sommet som, le premier index de la liste est index = fa7VoisinesSom_index[som] la premiere facette est alors fa7VoisinesSom_data(index,1) et l’indice suivant index = fa7VoisinesSom_data(index,0) La chaine est terminee lorsque index = = - 1.

calculer_correction_deplacement

int Remaillage_FT::calculer_correction_deplacement(DoubleTab &deplacement, const ArrOfDouble &varVolume, const DoubleTab &deplacement_varVolume, const ArrOfDouble &norme2_deplacement_varVolume) const

Cette fonction calcule une correction sur un deplacement liee a une variation de volume imposee Utile pour IJK.

calculer_differentielle_volume

int Remaillage_FT::calculer_differentielle_volume(const Maillage_FT_Disc &maillage, DoubleTab &differentielle_volume) const

Calcul de la differentielle du volume de phase 0 par rapport au deplacement de chaque sommet de l’interface.

C’est une forme lineaire qu’on exprime sous la forme d’un vecteur v tel que differentielle(deplacement) = deplacement scalaire v. En un certain sens, v est la normale a l’interface evaluee aux sommets. Si on deplace le sommet dans un plan orthogonal au vecteur, le volume des phases est conserve a l’ordre 1 par rapport a l’amplitude du deplacement.

calculer_longueurIdeale2_arete

double Remaillage_FT::calculer_longueurIdeale2_arete(const Maillage_FT_Disc &maillage, int som0, double x, double y, double z) const

Cette fonction calcule le carre de la longueur ideale d’une arete Dans un premier, cette longueur ideale est la racine cubique (en 3D) de la moyenne.

des volumes des elements euleriens contenant les sommets

calculer_variation_volume_facette_2D

double Remaillage_FT::calculer_variation_volume_facette_2D(int fa7, const Maillage_FT_Disc &maillage, const DoubleTab &position_initiale) const

Cette fonction calcule la difference de volume au niveau d’une facette par rapport a une position initiale des sommets.

version 2D

calculer_variation_volume_facette_3D

double Remaillage_FT::calculer_variation_volume_facette_3D(int fa7, const Maillage_FT_Disc &maillage, const DoubleTab &position_initiale) const

calculer_volume_sommets_supprimes

double Remaillage_FT::calculer_volume_sommets_supprimes(const Maillage_FT_Disc &maillage, const ArrOfInt &tab_somSupp, ArrOfDouble &varVolume) const

Cette fonction calcule la variation de volume liee a la suppression de sommets.

chercher_arete_tab

int Remaillage_FT::chercher_arete_tab(int tmp, const ArrOfInt &tab_index, const IntTab &tab_aretes, int pe0, int numOwner0, int pe1, int numOwner1) const

diviser_grandes_aretes

int Remaillage_FT::diviser_grandes_aretes(Maillage_FT_Disc &maillage) const

Cette fonction divise les grandes aretes en 2.

inserer_tab_aretes

int Remaillage_FT::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

marquer_aretes

int Remaillage_FT::marquer_aretes(Maillage_FT_Disc &maillage, IntTab &tab_aretesMarquees, ArrOfInt &tab_somD, DoubleTab &tab_deplacement_somD, int drap) const

nettoyer_maillage

int Remaillage_FT::nettoyer_maillage(Maillage_FT_Disc &maillage) const

Cette fonction nettoie le maillage : Elle supprime les facettes reduites a 1 sommet.

Elle supprime les sommets non utilises

permuter_aretes

int Remaillage_FT::permuter_aretes(Maillage_FT_Disc &maillage) const

Cette fonction effectue des permutations d’aretes afin d’ameliorer la qualite du maillage.

printOn

Sortie & Remaillage_FT::printOn(Sortie &x) const override

Cette fonction ne doit pas etre appelee.

Reimplements:

printOn

qualiteTriangle

double Remaillage_FT::qualiteTriangle(const FTd_vecteur3 &som0, const FTd_vecteur3 &som1, const FTd_vecteur3 &som2, double &aire) const

Cette methode calcule, pour un triangle donne, sa qualite : celle-ci est comprise dans ]0,1], et vaut 1 pour un triangle equilateral.

readOn

Entree & Remaillage_FT::readOn(Entree &) override

Lecture des parametres de remaillage.

Reimplements:

readOn

redistribuer_sommets

double Remaillage_FT::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

Deplace les sommets du maillage pour les redistribuer de facon homogene.

Le deplacement est la somme d’une composante tangentielle et d’une composante normale a l’interface. La composante tangentielle ramene le sommet vers le barycentre des facettes voisines du sommet. La composante normale est calculee de sorte a produire la variation de volume imposee. Traitement des lignes de contact: le deplacement est decompose en une composante normale a la ligne de contact et une composante tangente.

regulariser_maillage

double Remaillage_FT::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

Algorithme general de lissage du maillage.

Permet de barycentrer, regulariser la courbure, et d’appliquer une correction de volume en deplacant les noeuds dans la direction normale.

supprimer_doublons_facettes

int Remaillage_FT::supprimer_doublons_facettes(Maillage_FT_Disc &maillage) const

Cette fonction marque a supprimer les facettes en double.

Lorsque 2 facettes se trouvent avoir 3 de leur sommets en commun, il faut supprimer les 2 facettes? Marquer a supprimer = condenser les 3 sommets en un seul (le sommet 0) Precondition: Le maillage doit verifier la propriete “proprietaire facette = premier sommet”

supprimer_facettes_bord

int Remaillage_FT::supprimer_facettes_bord(Maillage_FT_Disc &maillage) const

Cette fonction marque a supprimer les facettes ayant leurs 3 sommets de bord Marquer a supprimer = condenser les 3 sommets en un seul (le sommet 0)

MODIF GB 21/07/2015 : Pour le cas sloshing par exemple, si le domaine a des coins, l’interface peut contenir un element dont les 3 sommets sont sur le bord sans pour autant que la facette soit integralement sur le bord. Elle peut avoir son 3eme cote dans le domaine. Dans ce cas, on la conserve.

supprimer_petites_aretes

int Remaillage_FT::supprimer_petites_aretes(Maillage_FT_Disc &maillage, ArrOfDouble &varVolume) const

A l’aide de “marquer_aretes”, on determine les aretes trop petites du maillage.

On les supprime en remplacant un sommet d’une extremite de l’arete par le sommet de l’autre extremite. Pour cela, on cree des sommets virtuels, des facettes nulles (deux ou trois sommets confondus) et on change le volume des phases. La variation de volume engendree est ajoutee a varVolume.

tester_a_remailler

int Remaillage_FT::tester_a_remailler(const Maillage_FT_Disc &maillage) const

Attributes Documentation#

info_obj (public)

const Type_info Remaillage_FT::info_obj

critere_arete_ (protected)

double Remaillage_FT::critere_arete_ = 0.35

dt_lissage_ (protected)

double Remaillage_FT::dt_lissage_ = - 1.

dt_remaillage_ (protected)

double Remaillage_FT::dt_remaillage_ = - 1.

equilateral_ (protected)

int Remaillage_FT::equilateral_ = 0

facteur_longueur_ideale_ (protected)

double Remaillage_FT::facteur_longueur_ideale_ = - 1.

impr_ (protected)

int Remaillage_FT::impr_ = - 1

is_solid_particle_ (protected)

bool Remaillage_FT::is_solid_particle_ = false

lissage_courbure_coeff_ (protected)

double Remaillage_FT::lissage_courbure_coeff_ = - 0.05

lissage_courbure_iterations_old_ (protected)

int Remaillage_FT::lissage_courbure_iterations_old_ = - 1

lissage_courbure_iterations_si_remaillage_ (protected)

int Remaillage_FT::lissage_courbure_iterations_si_remaillage_ = 0

lissage_courbure_iterations_systematique_ (protected)

int Remaillage_FT::lissage_courbure_iterations_systematique_ = 0

lissage_critere_ (protected)

double Remaillage_FT::lissage_critere_ = 0.

nb_iter_bary_volume_seul_ (protected)

int Remaillage_FT::nb_iter_bary_volume_seul_ = 3

nb_iter_barycentrage_ (protected)

int Remaillage_FT::nb_iter_barycentrage_ = 1

nb_iter_remaillage_ (protected)

int Remaillage_FT::nb_iter_remaillage_ = 1

refdomaine_VF_ (protected)

TRUST_Ref<Domaine_VF *> Remaillage_FT::refdomaine_VF_

relax_barycentrage_ (protected)

double Remaillage_FT::relax_barycentrage_ = 1.

seuil_dvolume_residuel_ (protected)

double Remaillage_FT::seuil_dvolume_residuel_ = 0.

surface_interface_ (protected)

double Remaillage_FT::surface_interface_ = 0.

temps_ (protected)

double Remaillage_FT::temps_ = 0

temps_dernier_lissage_ (protected)

double Remaillage_FT::temps_dernier_lissage_ = - 1.e40

temps_dernier_remaillage_ (protected)

double Remaillage_FT::temps_dernier_remaillage_ = - 1.e40

valeur_longueur_fixe_ (protected)

double Remaillage_FT::valeur_longueur_fixe_ = - 1.

variation_volume_ (protected)

double Remaillage_FT::variation_volume_ = 0.

Remaillage_FT

Contents

Remaillage_FT#

How to cite this class in this doc#

Detailed description#

Inherits from#

Inherited by#

List of Public Methods#

List of Public Static Methods#

List of Protected Methods#

Complete Member Function Documentation#

Attributes Documentation#