TRUST 1.9.8
HPC thermohydraulic platform
Loading...
Searching...
No Matches
Op_Diff_DG_base Class Reference

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_basediffusivite () 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 &)
 Associe la vraie diffusivite en m^2/s (en QC par exemple, l'operateur est applique a rho*u, et on associe alors la viscosite dynamique a.
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_basecorrelation_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 &)
 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 Matriceget_matrice () const
Matriceset_matrice ()
const SolveurSysget_solveur () const
SolveurSysset_solveur ()
Entreelire_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_basemon_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_baseget_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 SortieprintOn (Sortie &) const
 Ecriture de l'objet sur un flot de sortie Methode a surcharger.
virtual EntreereadOn (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 Nomle_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_infoget_info () const
 Donne des informations sur le type de l'Objet_U.
const Nomque_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 Interpreteinterprete () const
Interpreteinterprete ()
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_baseequation () const
 Renvoie la reference sur l'equation pointe par MorEqn::mon_equation.
Equation_baseequation ()
 Renvoie la reference sur l'equation pointe par 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 ()
 Constructeur de la classe.
virtual ~Support_Champ_Masse_Volumique ()
 Destructeur virtuel (pour eviter les warnings).
virtual void associer_champ_masse_volumique (const Champ_base &)
 Methode a appeler dans la preparation du probleme pour demander a l'objet de tenir compte du champ de masse volumique en parametre.
virtual void dissocier_champ_masse_volumique ()
 Annule la reference a la masse volumique.
virtual int has_champ_masse_volumique () const
 Renvoie 1 si la masse volumique a ete associee, 0 sinon.
virtual const Champ_baseget_champ_masse_volumique () const
 Renvoie le champ de masse volumique.
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_basediffusivite_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_basediffusivite_pour_pas_de_temps () const
 Renvoie le champ_don correspondant a la vraie diffusivite du milieu qui sert pour le calcul du pas de temps.
 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 ()
 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_Uoperator= (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 Member Functions inherited from Support_Champ_Masse_Volumique
virtual void declare_support_masse_volumique (int ok)
 Le constructeur d'une classe derivee qui se sert de la masse volumique doit appeler cette fonction avec la valeur 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 Nomnom_du_cas ()
 Renvoie une reference constante vers le nom du cas.
static Nomget_set_nom_du_cas ()
 Renvoie une reference non constante vers le nom du cas (pour pouvoir le modifier).
static const Type_infoinfo ()
 Donne des informations sur le type de l'Objet_U.
static const Objet_Uself_cast (const Objet_U &)
 methode ajoutee pour caster en python
static Objet_Uself_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 SortieJournal (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_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
incoThe input field (e.g., temperature, velocity component).
resuThe 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.

Imprime l'operateur sur un flot de sortie.

Parameters
(Sortie&os)
Returns
(int) renvoie toujours 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.

Mise a jour de l'operateur

Parameters
(double)un pas de temps

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:

  1. 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.
  2. 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/checkouts/readthedocs.org/user_builds/cea-trust-platform/checkouts/v1.9.8/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_base.h
  • /home/docs/checkouts/readthedocs.org/user_builds/cea-trust-platform/checkouts/v1.9.8/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_base.cpp