cs_thermal_system.h

Go to the documentation of this file.

#ifndef CS_THERMAL_SYSTEM_H
#define CS_THERMAL_SYSTEM_H
 
/*============================================================================
 * Functions to handle cs_thermal_system_t structure
 *============================================================================*/
 
/*
  This file is part of code_saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2026 EDF S.A.
 
  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.
 
  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/
 
#include "base/cs_field.h"
#include "base/cs_time_step.h"
#include "cdo/cs_advection_field.h"
#include "cdo/cs_equation.h"
#include "cdo/cs_property.h"
#include "cdo/cs_source_term.h"
#include "cdo/cs_xdef.h"
#include "mesh/cs_mesh.h"
 
/*============================================================================
 * Macro definitions
 *============================================================================*/
 
/* Generic name given to fields and equations related to this module */
 
#define CS_THERMAL_EQNAME           "thermal_equation"
#define CS_THERMAL_CP_NAME          "thermal_capacity"
#define CS_THERMAL_LAMBDA_NAME      "thermal_conductivity"
 
#define CS_THERMAL_POST_ENTHALPY               (1 << 0) /* 1 */
 
/*=============================================================================
 * Structure and type definitions
 *============================================================================*/
 
typedef cs_flag_t  cs_thermal_model_type_t;
 
typedef enum {
 
  CS_THERMAL_MODEL_STEADY                     = 1<<0,  /* =  1 */
  CS_THERMAL_MODEL_ADVECTION                  = 1<<1,  /* =  2 */
 
  /* Physical model associated with the thermal system */
 
  CS_THERMAL_MODEL_NAVSTO                     = 1<<2,  /* =  4 */
 
  /* Main variable to consider (by default the temperature in Kelvin)
     ---------------------------------------------------------------- */
 
  CS_THERMAL_MODEL_USE_TEMPERATURE            = 1<<3,  /* =  8 */
  CS_THERMAL_MODEL_USE_ENTHALPY               = 1<<4,  /* = 16 */
  CS_THERMAL_MODEL_USE_TOTAL_ENERGY           = 1<<5,  /* = 32 */
 
  /* Treatment of the diffusion term
     ------------------------------- */
 
  CS_THERMAL_MODEL_ANISOTROPIC_CONDUCTIVITY   = 1<<6,  /* = 64 */
 
  /* Additional bit settings
     ----------------------- */
 
  CS_THERMAL_MODEL_IN_CELSIUS                 = 1<<7   /* = 128 */
 
} cs_thermal_model_type_bit_t;
 
/* Set of parameters related to the thermal module */
 
typedef struct {
 
  cs_flag_t   model;                  /* Choice of modelling */
  cs_flag_t   numeric;                /* General numerical options */
  cs_flag_t   post;                   /* Post-processing options */
 
  /* Equation associated to this module */
  /* ---------------------------------- */
 
  cs_equation_t  *thermal_eq;
 
  /* Properties associated to this module */
  /* ------------------------------------ */
 
  cs_property_t  *unsteady_property; /* Property of the unsteady term: rho.cp */
  cs_property_t  *lambda;            /* Thermal conductivity */
  cs_property_t  *cp;                /* Heat capacity */
  cs_property_t  *rho;               /* Mass density (may be shared) */
  cs_property_t  *kappa;             /* lambda/cp may be null*/
 
  /* value of lambda/cp in each cell (allocated if the related property kappa
   * is non-null and defined by array) */
  cs_real_t      *kappa_array;
 
  /* Fields associated to this module */
  /* -------------------------------- */
 
  cs_field_t     *temperature;
  cs_field_t     *enthalpy;
  cs_field_t     *total_energy;
 
  /* Additional members */
  /* ------------------ */
 
  cs_real_t       ref_temperature;
 
  /* N.B.: Other reference values for properties are stored within each
   * property structure */
 
} cs_thermal_system_t;
 
/*============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_thermal_system_get_reference_temperature(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_set_reference_temperature(cs_real_t    ref);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_flag_t
cs_thermal_system_get_model(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
bool
cs_thermal_system_needs_navsto(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
bool
cs_thermal_system_is_activated(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_thermal_system_t *
cs_thermal_system_activate(cs_flag_t         model,
                           cs_flag_t         numeric,
                           cs_flag_t         post);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_destroy(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_equation_t *
cs_thermal_system_get_equation(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t *
cs_thermal_system_get_face_temperature(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_init_setup(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_finalize_setup(const cs_cdo_connect_t     *connect,
                                 const cs_cdo_quantities_t  *quant,
                                 const cs_time_step_t       *time_step);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_compute_steady_state(const cs_mesh_t              *mesh,
                                       const cs_cdo_connect_t       *connect,
                                       const cs_cdo_quantities_t    *quant,
                                       const cs_time_step_t         *time_step);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_compute(bool                          cur2prev,
                          const cs_mesh_t              *mesh,
                          const cs_cdo_connect_t       *connect,
                          const cs_cdo_quantities_t    *quant,
                          const cs_time_step_t         *time_step);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_init_values(const cs_mesh_t             *mesh,
                              const cs_cdo_connect_t      *connect,
                              const cs_cdo_quantities_t   *quant,
                              const cs_time_step_t        *ts);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_update(const cs_mesh_t             *mesh,
                         const cs_cdo_connect_t      *connect,
                         const cs_cdo_quantities_t   *quant,
                         const cs_time_step_t        *ts,
                         bool                         cur2prev);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_extra_op(const cs_cdo_connect_t      *connect,
                           const cs_cdo_quantities_t   *cdoq);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_extra_post(void                      *input,
                             int                        mesh_id,
                             int                        cat_id,
                             int                        ent_flag[5],
                             cs_lnum_t                  n_cells,
                             cs_lnum_t                  n_i_faces,
                             cs_lnum_t                  n_b_faces,
                             const cs_lnum_t            cell_ids[],
                             const cs_lnum_t            i_face_ids[],
                             const cs_lnum_t            b_face_ids[],
                             const cs_time_step_t      *time_step);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_log_setup(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_system_cht_boundary_conditions_setup(void);
 
#endif /* CS_THERMAL_SYSTEM_H */