#include <Op_Conv_EF_Stab_PolyMAC_MPFA_Face.h>

Inheritance diagram for Op_Conv_EF_Stab_PolyMAC_MPFA_Face:

Collaboration diagram for Op_Conv_EF_Stab_PolyMAC_MPFA_Face:

Public Member Functions
void	completer () override
	Associe l'operateur au domaine_dis, le domaine_Cl_dis, et a l'inconnue de son equation.
double	calculer_dt_stab () const override
	Computes the stability time step for the convection operator.
void	dimensionner_blocs (matrices_t matrices, const tabs_t &semi_impl={ }) const override
	Dimensions the matrix blocks for the PolyMAC_MPFA Face convection operator with EF stabilization.
void	ajouter_blocs (matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={ }) const override
	Adds convection contributions to the system matrix and right-hand side vector.
Public Member Functions inherited from Op_Conv_EF_Stab_PolyMAC_HFV_Face
int	has_interface_blocs () const override
void	check_multiphase_compatibility () const override
Public Member Functions inherited from Op_Conv_EF_Stab_PolyMAC_CDO_Face
DoubleTab &	ajouter (const DoubleTab &inco, DoubleTab &resu) const override
void	contribuer_a_avec (const DoubleTab &, Matrice_Morse &) const override
	on assemble la matrice.
void	modifier_pour_Cl (Matrice_Morse &, DoubleTab &) const override
	DOES NOTHING - to override in derived classes.
void	dimensionner (Matrice_Morse &mat) const override
	DOES NOTHING - to override in derived classes.
void	set_incompressible (const int flag) override
Public Member Functions inherited from Op_Conv_PolyMAC_CDO_base
void	completer () override
	Associe l'operateur au domaine_dis, le domaine_Cl_dis, et a l'inconnue de son equation.
double	calculer_dt_stab () const override
	Calcul dt_stab.
DoubleTab &	calculer (const DoubleTab &inco, DoubleTab &resu) const override
	calcule la contribution de la convection, la range dans resu renvoie resu
int	impr (Sortie &os) const override
	DOES NOTHING - to override in derived classes.
void	associer_domaine_cl_dis (const Domaine_Cl_dis_base &) override
void	associer (const Domaine_dis_base &, const Domaine_Cl_dis_base &, const Champ_Inc_base &) override
void	associer_vitesse (const Champ_base &) override
Public Member Functions inherited from Operateur_Conv_base
double	dt_stab_conv () const
void	fixer_dt_stab_conv (double dt)
virtual void	associer_norme_vitesse (const Champ_base &norme_vitesse)
virtual void	associer_vitesse_pour_pas_de_temps (const Champ_base &vitesse)
virtual void	associer_champ_temp (const Champ_Inc_base &, bool) const
virtual void	set_transporting_velocity_phase_index (int idx)
Public Member Functions inherited from Operateur_base
virtual void	associer_champ (const Champ_Inc_base &, const std::string &nom_ch)
virtual void	dimensionner_bloc_vitesse (Matrice_Morse &matrice) const
virtual void	contribuer_bloc_vitesse (const DoubleTab &, Matrice_Morse &) const
virtual void	contribuer_au_second_membre (DoubleTab &) const
	DOES NOTHING - to override in derived classes.
void	tester_contribuer_a_avec (const DoubleTab &, const Matrice_Morse &)
virtual void	dimensionner_termes_croises (Matrice_Morse &, const Probleme_base &autre_pb, int nl, int nc) const
virtual void	ajouter_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, DoubleTab &resu) const
virtual void	contribuer_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, Matrice_Morse &matrice) const
virtual void	calculer_dt_local (DoubleTab &) const
virtual void	mettre_a_jour (double temps)
	DOES NOTHING - to override in derived classes.
virtual void	abortTimeStep ()
virtual void	resetTime (double time)
void	associer_eqn (const Equation_base &)
	Associe une equation a l'operateur.
int	get_decal_temps () const
int	set_decal_temps (int)
int	get_nb_ss_pas_de_temps () const
int	set_nb_ss_pas_de_temps (int)
const Matrice &	get_matrice () const
Matrice &	set_matrice ()
const SolveurSys &	get_solveur () const
SolveurSys &	set_solveur ()
Entree &	lire_solveur (Entree &)
virtual int	systeme_invariant () const
virtual void	ajouter_contribution_explicite_au_second_membre (const Champ_Inc_base &inconnue, DoubleTab &derivee) const
const Champ_Inc_base &	mon_inconnue () const
bool	has_champ_inco () const
const std::string &	nom_inconnue () const
void	ouvrir_fichier (SFichier &os, const Nom &, const int flag=1) const
	Ouverture/creation d'un fichier d'impression d'un operateur A surcharger dans les classes derivees.
void	ouvrir_fichier_partage (EcrFicPartage &, const Nom &, const int flag=1) const
	Ouverture/creation d'un fichier d'impression d'un operateur A surcharger dans les classes derivees.
void	set_fichier (const Nom &)
const Nom	fichier () const
void	set_description (const Nom &nom)
const Nom	description () const
DoubleTab &	flux_bords ()
DoubleTab &	flux_bords () const
void	creer_champ (const Motcle &motlu) override
const Champ_base &	get_champ (const Motcle &nom) const override
bool	has_champ (const Motcle &nom, OBS_PTR(Champ_base) &ref_champ) const override
bool	has_champ (const Motcle &nom) const override
void	get_noms_champs_postraitables (Noms &nom, Option opt=NONE) const override
void	calculer_pour_post (Champ_base &espace_stockage, const Nom &option, int comp) const override
Motcle	get_localisation_pour_post (const Nom &option) const override
virtual void	ajouter_flux (const DoubleTab &inconnue, DoubleTab &contribution) const
virtual void	calculer_flux (const DoubleTab &inconnue, DoubleTab &flux) const
virtual void	preparer_calcul ()
bool	has_impr_file () const
Public Member Functions inherited from Objet_U
	~Objet_U () override
	Destructeur, supprime l'objet de la liste d'objets enregistres dans "memoire".
int	numero () const
	Renvoie l'indice de l'objet dans Memoire::data.
virtual int	duplique () const =0
virtual Sortie &	printOn (Sortie &) const
	Ecriture de l'objet sur un flot de sortie Methode a surcharger.
virtual Entree &	readOn (Entree &)
	Lecture d'un Objet_U sur un flot d'entree Methode a surcharger.
virtual unsigned	taille_memoire () const =0
virtual int	est_egal_a (const Objet_U &) const
	Renvoie 1 si l'objet x et *this sont une seule et meme instance (meme adresse en memoire).
virtual const Nom &	le_nom () const
	Donne le nom de l'Objet_U Methode a surcharger : renvoie "neant" dans cette implementation.
virtual void	nommer (const Nom &)
	Donne un nom a l'Objet_U Methode virtuelle a surcharger.
virtual int	reprendre (Entree &)
	Reprise d'un Objet_U sur un flot d'entree Methode a surcharger.
virtual int	sauvegarder (Sortie &) const
	Sauvegarde d'un Objet_U sur un flot de sortie Methode a surcharger.
int	get_object_id () const
	Renvoie l'identifiant unique de l'objet object_id_.
virtual const Type_info *	get_info () const
	Donne des informations sur le type de l'Objet_U.
const Nom &	que_suis_je () const
	renvoie la chaine identifiant la classe.
const char *	le_type () const
	Donne le nom du type de l'Objet_U.
virtual int	change_num (const int *const)
	Change le numero interne de l'Objet_U.
virtual int	lire_motcle_non_standard (const Motcle &motlu, Entree &is)
	Lecture des parametres de type non simple d'un objet_U a partir d'un flot d'entree.
virtual int	associer_ (Objet_U &)
	Associe l'Objet_U a un autre Objet_U Methode virtuelle a surcharger.
const Interprete &	interprete () const
Interprete &	interprete ()
Public Member Functions inherited from Process
virtual	~Process ()
Public Member Functions inherited from MorEqn
void	associer_eqn (const Equation_base &)
	Associe une equation a l'objet.
const Equation_base &	equation () const
	Renvoie la reference sur l'equation pointe par MorEqn::mon_equation.
Equation_base &	equation ()
	Renvoie la reference sur l'equation pointe par MorEqn::mon_equation.
int	mon_equation_non_nul () const
Public Member Functions inherited from Champs_compris_interface
virtual	~Champs_compris_interface ()
virtual bool	has_champ (const Motcle &nom, OBS_PTR(Champ_base)&ref_champ) const =0

Additional Inherited Members
Static Public Member Functions inherited from Objet_U
static const Nom &	nom_du_cas ()
	Renvoie une reference constante vers le nom du cas.
static Nom &	get_set_nom_du_cas ()
	Renvoie une reference non constante vers le nom du cas (pour pouvoir le modifier).
static const Type_info *	info ()
	Donne des informations sur le type de l'Objet_U.
static const Objet_U &	self_cast (const Objet_U &)
	methode ajoutee pour caster en python
static Objet_U &	self_cast (Objet_U &)
Static Public Member Functions inherited from Process
static int	me ()
	renvoie mon rang dans le groupe de communication courant.
static int	nproc ()
	renvoie le nombre de processeurs dans le groupe courant Voir Comm_Group::nproc() et PE_Groups::current_group()
static bool	is_parallel ()
static void	exit (int exit_code=-1)
	Routine de sortie de TRUST dans une region Kokkos.
static double	mp_sum (double)
	Calcule la somme de x sur tous les processeurs du groupe courant.
static float	mp_sum (float)
static trustIdType	mp_sum (trustIdType)
	Calcule la somme de x sur tous les processeurs du groupe courant.
static double	mp_max (double)
static double	mp_min (double)
static int	mp_max (int)
	renvoie le plus grand int i sur l'ensemble des processeurs du groupe courant.
static int	mp_min (int)
	renvoie le plus petit int i sur l'ensemble des processeurs du groupe courant.
static double	mp_sum_as_double (int v)
static trustIdType	mppartial_sum (trustIdType i)
	Calul de la somme partielle de i sur les processeurs 0 a me()-1 (renvoie 0 sur le processeur 0).
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2)
	C++14 compatible mp_sum_for_each: combine multiple mp_sum calls into one collective operation Usage: mp_sum_for_each(a, b); mp_sum_for_each(a, b, c); mp_sum_for_each(a, b, c, d); mp_sum_for_each(a, b, c, d, e); All arguments must be of the same type (double or int) and are modified in place. Supports 2-5 parameters.
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2)
	C++14 compatible mp_max_for_each: combine multiple mp_max calls into one collective operation.
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2)
	C++14 compatible mp_min_for_each: combine multiple mp_min calls into one collective operation.
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename _TYPE_>
static void	mp_sum_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
template<typename _TYPE_>
static void	mp_max_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
template<typename _TYPE_>
static void	mp_min_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
static bool	mp_and (bool)
	Calcule le 'et' logique de b sur tous les processeurs du groupe courant.
static bool	mp_or (bool)
static int	check_int_overflow (trustIdType)
static int	je_suis_maitre ()
	renvoie 1 si on est sur le processeur maitre du groupe courant (c'est a dire me() == 0), 0 sinon.
static KOKKOS_INLINE_FUNCTION void	Kokkos_exit (const char *)
	Routine de sortie de TRUST dans une region Kokkos.
static int	node_master ()
	renvoie 1 si on est sur le processeur maitre du noeud numa, 0 sinon.
static void	exit (const Nom &message, int exit_code=-1)
static bool	is_sequential ()
static void	barrier ()
	Synchronise tous les processeurs du groupe courant (attend que tous les processeurs soient arrives a la barriere).
static void	abort ()
	Routine de sortie de Trio-U sur une erreur abort().
static Sortie &	Journal (int message_level=0)
	Renvoie un objet statique de type Sortie qui sert de journal d'evenements.
static double	ram_processeur ()
static void	imprimer_ram_totale (int all_process=0)
static bool	force_single_file (const int ranks, const Nom &filename)
Public Attributes inherited from Operateur_base
int	col_width_
Static Public Attributes inherited from Objet_U
static double	precision_geom = 1e-10
static constexpr bool	HAS_POINTER = false
static int	dimension =0
static int	format_precision_geom =11
static int	axi =0
static int	bidim_axi =0
static int	DEACTIVATE_SIGINT_CATCH =0
static Type_info	info_obj
static bool	disable_TU =false
	Flag to disable or not the writing of the .TU files.
static bool	stat_per_proc_perf_log =false
	Flag to enable the writing of the statistics detailed per processor in _csv.TU file.
Static Public Attributes inherited from Process
static int	exception_sur_exit =0
static int	multiple_files =5120
Protected Member Functions inherited from Op_Conv_PolyMAC_CDO_base
	OBS_PTR (Domaine_PolyMAC_CDO) le_dom_poly_
	OBS_PTR (Domaine_Cl_PolyMAC_family) la_zcl_poly_
	OBS_PTR (Champ_base) vitesse_
Protected Member Functions inherited from Operateur_Conv_base
	OBS_PTR (Champ_base) la_norme_vitesse
	OBS_PTR (Champ_base) vitesse_pour_pas_de_temps_
Protected Member Functions inherited from Operateur_base
	OBS_PTR (Champ_Inc_base) le_champ_inco
Protected Member Functions inherited from Objet_U
	Objet_U ()
	Constructeur par defaut : attribue un numero d'identifiant unique a l'objet (object_id_), et enregistre l'objet en "memoire".
	Objet_U (const Objet_U &)
	Constructeur par copie.
const Objet_U &	operator= (const Objet_U &)
	Operateur= : ne fait rien (on conserve le numero et l'identifiant).
virtual void	set_param (Param &) const
Protected Member Functions inherited from MorEqn
	OBS_PTR (Equation_base) mon_equation
virtual	~MorEqn ()
Protected Attributes inherited from Op_Conv_EF_Stab_PolyMAC_CDO_Face
double	alpha_ = 1
DoubleVect	porosite_f
DoubleVect	porosite_e
Protected Attributes inherited from Op_Conv_PolyMAC_CDO_base
SFichier	Flux
SFichier	Flux_moment
SFichier	Flux_sum
Protected Attributes inherited from Operateur_Conv_base
double	dt_stab_conv_
int	incompressible_ = 1
Protected Attributes inherited from Operateur_base
int	decal_temps
int	nb_ss_pas_de_temps
SolveurSys	solveur
Matrice	matrice_
Nom	out_
Nom	description_
DoubleTab	flux_bords_
Champs_compris	champs_compris_
std::string	nom_inco_

Detailed Description

Definition at line 21 of file Op_Conv_EF_Stab_PolyMAC_MPFA_Face.h.

Member Function Documentation

◆ ajouter_blocs()

void Op_Conv_EF_Stab_PolyMAC_MPFA_Face::ajouter_blocs	(	matrices_t	matrices,
		DoubleTab &	secmem,
		const tabs_t &	semi_impl = { } ) const

overridevirtual

Adds convection contributions to the system matrix and right-hand side vector.

This method implements the convection operator with EF stabilization for PolyMAC_MPFA Face discretization. It computes face-based convection fluxes using upwind stabilization and assembles the corresponding matrix coefficients and source terms. The method handles both face-face and element-element contributions, accounting for face equivalence relationships and boundary conditions.

Parameters

matrices	Map containing the system matrices indexed by unknown field names
secmem	Right-hand side vector to be updated with convection contributions
semi_impl	Map of semi-implicit field values indexed by field names

Note: The convection flux computation uses a blended upwind scheme with stabilization parameter alpha_.; For multiphase flows, added mass correlations are included in the mass matrix contributions.; The method distinguishes between compressible and incompressible flow formulations.; Dirichlet boundary conditions are enforced through direct substitution in the source term.; Face equivalence relationships are handled to maintain consistency across mesh interfaces.; Performance statistics are automatically tracked during execution.

Reimplemented from Op_Conv_EF_Stab_PolyMAC_HFV_Face.

Definition at line 227 of file Op_Conv_EF_Stab_PolyMAC_MPFA_Face.cpp.

◆ calculer_dt_stab()

double Op_Conv_EF_Stab_PolyMAC_MPFA_Face::calculer_dt_stab ( ) const

overridevirtual

Computes the stability time step for the convection operator.

This method calculates the maximum allowable time step based on the CFL stability condition for the convection operator. The computation considers incoming fluxes through each face of every element and determines the most restrictive time step constraint across the entire domain. For multiphase flows, elements with negligible phase presence are excluded from the time step calculation.

Returns: The global minimum stability time step across all MPI processes

Note: Only incoming convective fluxes are considered in the stability analysis.; For multiphase problems, elements with phase fraction below 1e-3 are ignored.; In axisymmetric calculations, boundary faces with specific conditions may be excluded.; The final time step is synchronized across all MPI processes using the minimum value.; Flux contributions account for face porosity and surface area weighting.

Reimplemented from Op_Conv_EF_Stab_PolyMAC_HFV_Face.

Definition at line 79 of file Op_Conv_EF_Stab_PolyMAC_MPFA_Face.cpp.

◆ completer()

void Op_Conv_EF_Stab_PolyMAC_MPFA_Face::completer ( )

overridevirtual

Associe l'operateur au domaine_dis, le domaine_Cl_dis, et a l'inconnue de son equation.

Exceptions

pas	d'equation associee

Reimplemented from Op_Conv_EF_Stab_PolyMAC_CDO_Face.

Definition at line 57 of file Op_Conv_EF_Stab_PolyMAC_MPFA_Face.cpp.

◆ dimensionner_blocs()

void Op_Conv_EF_Stab_PolyMAC_MPFA_Face::dimensionner_blocs	(	matrices_t	matrices,
		const tabs_t &	semi_impl = { } ) const

overridevirtual

Dimensions the matrix blocks for the PolyMAC_MPFA Face convection operator with EF stabilization.

This method constructs the sparsity pattern and allocates memory for the system matrix by analyzing face-element connectivity and establishing the stencil relationships between degrees of freedom. The resulting matrix structure accounts for face-face and element-element contributions based on the PolyMAC_MPFA Face discretization scheme.

Parameters

matrices	Map containing the system matrices indexed by unknown field names
semi_impl	Map of semi-implicit terms indexed by unknown field names

Note: The method returns early if no diagonal block exists or if semi-implicit treatment is requested.; For multiphase problems with added mass correlation, cross-coupling terms between phases are included.; The stencil construction handles face equivalence relationships and boundary conditions appropriately.; Duplicate entries in the stencil are automatically removed before matrix allocation.

Reimplemented from Op_Conv_EF_Stab_PolyMAC_HFV_Face.

Definition at line 130 of file Op_Conv_EF_Stab_PolyMAC_MPFA_Face.cpp.

The documentation for this class was generated from the following files:

/home/docs/checkouts/readthedocs.org/user_builds/cea-trust-platform/checkouts/master/src/PolyMAC_family/Operateurs/Op_Conv/Op_Conv_EF_Stab_PolyMAC_MPFA_Face.h
/home/docs/checkouts/readthedocs.org/user_builds/cea-trust-platform/checkouts/master/src/PolyMAC_family/Operateurs/Op_Conv/Op_Conv_EF_Stab_PolyMAC_MPFA_Face.cpp

Public Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ ajouter_blocs()

◆ calculer_dt_stab()

◆ completer()

◆ dimensionner_blocs()