This class provides the common infrastructure shared by all DG diffusion operators in TRUST. It handles the association with the DG domain and its boundary conditions, the management of an effective diffusivity field (nu_), and the computation of a stable explicit time step. More...

#include <Op_Diff_DG_base.h>

Inheritance diagram for Op_Diff_DG_base:

[legend]

Collaboration diagram for Op_Diff_DG_base:

[legend]

Public Member Functions
void	associer (const Domaine_dis_base &, const Domaine_Cl_dis_base &, const Champ_Inc_base &) override
	Associates the operator with a DG domain and its boundary conditions.
double	calculer_dt_stab () const override
	Computes the maximum stable explicit time step for the diffusion operator.
void	associer_diffusivite (const Champ_base &diffu) override
void	completer () override
	Finalizes the operator setup after all associations have been made.
const Champ_base &	diffusivite () const override
void	mettre_a_jour (double t) override
	DOES NOTHING - to override in derived classes.
void	update_nu () const
	Updates the cached effective diffusivity field nu_ by combining molecular and turbulent contributions.
double	nu (int i, int compo) const
DoubleTab &	calculer (const DoubleTab &, DoubleTab &) const override
	Computes the diffusion operator applied to inco and stores the result in resu.
int	impr (Sortie &os) const override
	DOES NOTHING - to override in derived classes.
Public Member Functions inherited from Operateur_Diff_base
virtual void	associer_diffusivite_pour_pas_de_temps (const Champ_base &)
	Associates the true diffusivity in m^2/s (in QC for example, the operator is applied to rho*u, and the dynamic viscosity is then associated to.
virtual void	associer_diffusivite_volumique (const Champ_base &)
virtual void	calculer_borne_locale (DoubleVect &, double, double) const
virtual void	init_op_ext () const
virtual bool	is_turb () const
virtual const Correlation_base *	correlation_viscosite_turbulente () const
virtual void	calculer_von_mises (const DoubleTab &deplacement, DoubleTab &deformation, DoubleTab &contraintes, DoubleTab &von_mises) const
Public Member Functions inherited from Operateur_base
virtual DoubleTab &	ajouter (const DoubleTab &, DoubleTab &) const
virtual void	associer_champ (const Champ_Inc_base &, const std::string &nom_ch)
virtual void	associer_domaine_cl_dis (const Domaine_Cl_dis_base &)
virtual void	dimensionner (Matrice_Morse &) const
	DOES NOTHING - to override in derived classes.
virtual void	dimensionner_bloc_vitesse (Matrice_Morse &matrice) const
virtual void	modifier_pour_Cl (Matrice_Morse &, DoubleTab &) const
	DOES NOTHING - to override in derived classes.
virtual void	contribuer_a_avec (const DoubleTab &, Matrice_Morse &) const
	DOES NOTHING - to override in derived classes.
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 int	has_interface_blocs () const
virtual void	dimensionner_blocs (matrices_t matrices, const tabs_t &semi_impl={ }) const
virtual void	ajouter_blocs (matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={ }) const
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	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
virtual void	check_multiphase_compatibility () 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
Public Member Functions inherited from Support_Champ_Masse_Volumique
	Support_Champ_Masse_Volumique ()
	Constructor of the class.
virtual	~Support_Champ_Masse_Volumique ()
	Virtual destructor (to avoid warnings).
virtual void	associer_champ_masse_volumique (const Champ_base &)
	Method to be called during problem preparation to ask the object to take into account the density field passed as parameter.
virtual void	dissocier_champ_masse_volumique ()
	Cancels the reference to the density field.
virtual int	has_champ_masse_volumique () const
	Returns 1 if the density field has been associated, 0 otherwise.
virtual const Champ_base &	get_champ_masse_volumique () const
	Returns the density field.
const int &	support_ok () const
Public Member Functions inherited from Op_Diff_Turbulent_base
virtual	~Op_Diff_Turbulent_base ()
void	associer_diffusivite_turbulente (const Champ_Fonc_base &)
const Champ_Fonc_base &	diffusivite_turbulente () const
bool	has_diffusivite_turbulente () const

Protected Member Functions
	OBS_PTR (Domaine_DG) le_dom_dg_
	OBS_PTR (Domaine_Cl_DG) la_zcl_dg_
	OBS_PTR (Champ_base) diffusivite_
Protected Member Functions inherited from Operateur_Diff_base
virtual const Champ_base &	diffusivite_pour_pas_de_temps () const
	Returns the field corresponding to the true diffusivity of the medium used for the time step computation.
	OBS_PTR (Champ_base) diffusivite_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 Member Functions inherited from Support_Champ_Masse_Volumique
virtual void	declare_support_masse_volumique (int ok)
	The constructor of a derived class that uses the density field must call this function with the value 1.

Protected Attributes
SFichier	Flux
SFichier	Flux_moment
SFichier	Flux_sum
int	nu_a_jour_ = 0
DoubleTab	nu_
bool	is_var_
bool	is_aniso_
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_

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_Diff_base
std::vector< const Operateur_Diff_base * >	op_ext
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

Detailed Description

This class provides the common infrastructure shared by all DG diffusion operators in TRUST. It handles the association with the DG domain and its boundary conditions, the management of an effective diffusivity field (nu_), and the computation of a stable explicit time step.

The stable time step calculation (calculer_dt_stab) follows a diffusive CFL criterion: dt ~ h^2 / (2 * dim * alpha_max) with special handling for:

Variable-density flows: rho * h^2 / (2 * dim * nu), with VDF-like and general mesh branches.
Robin / external-heat-exchange boundary conditions (Echange_externe_impose): the effective diffusivity is scaled by the Biot number when Bi > 1, making the criterion more conservative.

Derived classes are responsible for implementing the actual flux assembly (ajouter()) according to their specific DG formulation.

Definition at line 42 of file Op_Diff_DG_base.h.

Member Function Documentation

◆ associer()

void Op_Diff_DG_base::associer	(	const Domaine_dis_base &	domaine_dis,
		const Domaine_Cl_dis_base &	zcl,
		const Champ_Inc_base &	)

overridevirtual

Associates the operator with a DG domain and its boundary conditions.

Implements Operateur_base.

Definition at line 177 of file Op_Diff_DG_base.cpp.

◆ associer_diffusivite()

void Op_Diff_DG_base::associer_diffusivite ( const Champ_base & diffu )

inlineoverridevirtual

Implements Operateur_Diff_base.

Definition at line 50 of file Op_Diff_DG_base.h.

◆ calculer()

DoubleTab & Op_Diff_DG_base::calculer	(	const DoubleTab &	inco,
		DoubleTab &	resu ) const

overridevirtual

Computes the diffusion operator applied to inco and stores the result in resu.

Initializes resu to zero then delegates to ajouter(), which adds the diffusive contribution. This follows the standard TRUST operator pattern where calculer() resets the result before calling ajouter().

Parameters

inco	The input field (e.g., temperature, velocity component).
resu	The output field, zeroed then filled with the diffusion contribution.

Returns: A reference to resu.

Reimplemented from Operateur_base.

Definition at line 192 of file Op_Diff_DG_base.cpp.

◆ calculer_dt_stab()

double Op_Diff_DG_base::calculer_dt_stab ( ) const

overridevirtual

Computes the maximum stable explicit time step for the diffusion operator.

Two branches are handled depending on whether a variable density field is present:

Constant density (standard branch): The criterion is dt = h_min^2 / (2 * dim * alpha_max), where h_min is the minimum mesh size and alpha_max is the maximum diffusivity over all elements. This estimate is conservative: if the maximum diffusivity and the minimum mesh size are not co-located, the time step will be underestimated. Additionally, if a Robin boundary condition (Echange_externe_impose) is present, the Biot number Bi = h_imp_max * sqrt(h_min^2) / lambda_max is computed. If Bi > 1, alpha_max is scaled by the same factor, further tightening the criterion.
Variable density (rho field present): A per-element time step is computed as dt = 0.5 * rho / (nu * h), for VDF-like elements (exactly 2*dim faces), dt = rho * h^2 / (2 * dim * nu), for general elements, and the global minimum is taken across all elements and MPI processes.

Returns: The minimum stable time step across all elements and MPI processes.

Reimplemented from Operateur_base.

Definition at line 60 of file Op_Diff_DG_base.cpp.

◆ completer()

void Op_Diff_DG_base::completer ( )

overridevirtual

Finalizes the operator setup after all associations have been made.

Calls the parent completer(), then allocates and initializes the effective diffusivity table nu_. The number of components is set to 2 for the Transport_K_Epsilon equation (one per turbulent quantity), or to the line size of the diffusivity field for all other equations. The nu_a_jour_ flag is reset to 0 to force a recomputation on the first use.

Reimplemented from Operateur_base.

Reimplemented in Op_Diff_DG_Elem.

Definition at line 206 of file Op_Diff_DG_base.cpp.

◆ diffusivite()

const Champ_base & Op_Diff_DG_base::diffusivite ( ) const

inlineoverridevirtual

Implements Operateur_Diff_base.

Definition at line 58 of file Op_Diff_DG_base.h.

◆ impr()

int Op_Diff_DG_base::impr ( Sortie & os ) const

overridevirtual

DOES NOTHING - to override in derived classes.

Prints the operator to an output stream.

Parameters

(Sortie& os)

Returns: (int) always returns 1

Reimplemented from Operateur_base.

Definition at line 169 of file Op_Diff_DG_base.cpp.

◆ mettre_a_jour()

void Op_Diff_DG_base::mettre_a_jour ( double temps )

inlineoverridevirtual

DOES NOTHING - to override in derived classes.

Update of the operator

Parameters

(double) a time step

Reimplemented from Operateur_base.

Definition at line 59 of file Op_Diff_DG_base.h.

◆ nu()

double Op_Diff_DG_base::nu	(	int	i,
		int	compo ) const

inline

Definition at line 66 of file Op_Diff_DG_base.h.

◆ OBS_PTR() [1/3]

Op_Diff_DG_base::OBS_PTR ( Champ_base )

protected

◆ OBS_PTR() [2/3]

Op_Diff_DG_base::OBS_PTR ( Domaine_Cl_DG )

protected

◆ OBS_PTR() [3/3]

Op_Diff_DG_base::OBS_PTR ( Domaine_DG )

protected

◆ update_nu()

void Op_Diff_DG_base::update_nu ( ) const

Updates the cached effective diffusivity field nu_ by combining molecular and turbulent contributions.

This method is a no-op if nu_a_jour_ is already set. Otherwise, it fills nu_ as follows:

Molecular diffusivity (all equations except Transport_K_Epsilon):
- If the diffusivity has no metadata vector (i.e., it is spatially uniform), nu_ is filled by broadcasting the uniform value to all elements, treating both arrays as flat 1D vectors.
- Otherwise (spatially variable diffusivity), nu_ is directly injected from the diffusivity array, after asserting that both share the same parallel metadata vector.
Not usable for now: Turbulent diffusivity (if has_diffusivite_turbulente() is true):
Warning
This branch is currently unreachable: it starts with a throw statement, indicating that turbulent diffusivity is not yet supported for DG operators.
- For Transport_K_Epsilon, nu_(i,j) = nu_molecular(i) + nu_turb(i,j) for j in {0,1}.
- For other equations, nu_turb is added to nu_ element-wise, handling the uniform case (no metadata vector) and the variable case separately.
After the update, nu_a_jour_ is set to 1 to prevent redundant recomputations.

Definition at line 234 of file Op_Diff_DG_base.cpp.

Member Data Documentation

◆ Flux

SFichier Op_Diff_DG_base::Flux

mutableprotected

Definition at line 74 of file Op_Diff_DG_base.h.

◆ Flux_moment

SFichier Op_Diff_DG_base::Flux_moment

protected

Definition at line 74 of file Op_Diff_DG_base.h.

◆ Flux_sum

SFichier Op_Diff_DG_base::Flux_sum

protected

Definition at line 74 of file Op_Diff_DG_base.h.

◆ is_aniso_

bool Op_Diff_DG_base::is_aniso_

protected

Definition at line 81 of file Op_Diff_DG_base.h.

◆ is_var_

bool Op_Diff_DG_base::is_var_

protected

Definition at line 80 of file Op_Diff_DG_base.h.

◆ nu_

DoubleTab Op_Diff_DG_base::nu_

mutableprotected

Definition at line 78 of file Op_Diff_DG_base.h.

◆ nu_a_jour_

int Op_Diff_DG_base::nu_a_jour_ = 0

mutableprotected

Definition at line 77 of file Op_Diff_DG_base.h.

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

/home/docs/trust-code/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_base.h
/home/docs/trust-code/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_base.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ associer()

◆ associer_diffusivite()

◆ calculer()

◆ calculer_dt_stab()

◆ completer()

◆ diffusivite()

◆ impr()

◆ mettre_a_jour()

◆ nu()

◆ OBS_PTR() [1/3]

◆ OBS_PTR() [2/3]

◆ OBS_PTR() [3/3]

◆ update_nu()

Member Data Documentation

◆ Flux

◆ Flux_moment

◆ Flux_sum

◆ is_aniso_

◆ is_var_

◆ nu_

◆ nu_a_jour_