#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
	Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its 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
	Assemble the matrix.
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
	Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its equation.
double	calculer_dt_stab () const override
	Computes dt_stab.
DoubleTab &	calculer (const DoubleTab &inco, DoubleTab &resu) const override
	Computes the convection contribution and stores it in resu, then returns 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 &)
	Associates an equation with the operator.
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
	Opening/creation of a file for printing an operator. To override in derived classes.
void	ouvrir_fichier_partage (EcrFicPartage &, const Nom &, const int flag=1) const
	Opening/creation of a shared file for printing an operator. To override in derived classes.
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
	Destructor. Removes the object from the list of objects registered in "memory".
int	numero () const
	Returns the index of the object in Memoire::data.
virtual int	duplique () const =0
virtual Sortie &	printOn (Sortie &) const
	Writes the object to an output stream. Virtual method to override.
virtual Entree &	readOn (Entree &)
	Reads an Objet_U from an input stream. Virtual method to override.
virtual unsigned	taille_memoire () const =0
virtual int	est_egal_a (const Objet_U &) const
	Returns 1 if x and *this are the same instance (same memory address).
virtual const Nom &	le_nom () const
	Returns the name of the Objet_U. Virtual method to override: returns "neant" in this implementation.
virtual void	nommer (const Nom &)
	Assigns a name to the Objet_U. Virtual method to override.
virtual int	reprendre (Entree &)
	Restores an Objet_U from an input stream. Virtual method to override.
virtual int	sauvegarder (Sortie &) const
	Saves an Objet_U to an output stream. Virtual method to override.
int	get_object_id () const
	Returns the unique identifier of the object (object_id_).
virtual const Type_info *	get_info () const
	Returns type information for the Objet_U.
const Nom &	que_suis_je () const
	Returns the string identifying the class.
const char *	le_type () const
	Returns the type name of the Objet_U.
virtual int	change_num (const int *const)
	Changes the internal number of the Objet_U.
virtual int	lire_motcle_non_standard (const Motcle &motlu, Entree &is)
	Reads non-simple-type parameters of an Objet_U from an input stream.
virtual int	associer_ (Objet_U &)
	Associates the Objet_U with another Objet_U. Virtual method to override.
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 &)
	Associates an equation with the object.
const Equation_base &	equation () const
	Returns the reference to the equation pointed to by MorEqn::mon_equation.
Equation_base &	equation ()
	Returns the reference to the equation pointed to by 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 ()
	Returns a constant reference to the case name. This method is static.
static Nom &	get_set_nom_du_cas ()
	Returns a non-constant reference to the case name (to allow modification). This method is static.
static const Type_info *	info ()
	Returns type information for the Objet_U (static version).
static const Objet_U &	self_cast (const Objet_U &)
	Method added for casting in Python.
static Objet_U &	self_cast (Objet_U &)
Static Public Member Functions inherited from Process
static int	me ()
	Returns the rank of the local processor in the current communication group. See Comm_Group::rank() and PE_Groups::current_group().
static int	nproc ()
	Returns the number of processors in the current group. See Comm_Group::nproc() and PE_Groups::current_group().
static bool	is_parallel ()
static void	exit (int exit_code=-1)
	Exit routine for TRUST within a Kokkos region.
static double	mp_sum (double)
	Computes the sum of x over all processors in the current group.
static float	mp_sum (float)
static trustIdType	mp_sum (trustIdType)
	Computes the sum of x over all processors in the current group.
static double	mp_max (double)
static double	mp_min (double)
static int	mp_max (int)
	Returns the maximum value of x across all processors in the current group.
static int	mp_min (int)
	Returns the minimum value of x across all processors in the current group.
static double	mp_sum_as_double (int v)
static trustIdType	mppartial_sum (trustIdType i)
	Computes the partial sum of x over processors 0 to me()-1 (returns 0 on processor 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)
	Computes the logical AND of b across all processors in the current group.
static bool	mp_or (bool)
static int	check_int_overflow (trustIdType)
static int	je_suis_maitre ()
	Returns 1 if on the master processor of the current group (i.e. me() == 0), 0 otherwise.
static KOKKOS_INLINE_FUNCTION void	Kokkos_exit (const char *)
	Exit routine for TRUST within a Kokkos region.
static int	node_master ()
	Returns 1 if on the NUMA node master processor, 0 otherwise.
static void	exit (const Nom &message, int exit_code=-1)
static bool	is_sequential ()
static void	barrier ()
	Synchronizes all processors in the current group (waits until all processors have reached the barrier).
static void	abort ()
	Abort routine for TRUST on a fatal error.
static Sortie &	Journal (int message_level=0)
	Returns a static Sortie object used as an event journal.
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 ()
	Default constructor: assigns a unique identifier to the object (object_id_) and registers the object in "memory" by giving it a num_obj number.
	Objet_U (const Objet_U &)
	Copy constructor.
const Objet_U &	operator= (const Objet_U &)
	Assignment operator: does nothing (the number and identifier are preserved).
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

Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its equation.

Exceptions

no	associated equation

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/trust-code/src/PolyMAC_family/Operateurs/Op_Conv/Op_Conv_EF_Stab_PolyMAC_MPFA_Face.h
/home/docs/trust-code/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()