cs_thermal_model.h

Go to the documentation of this file.

#ifndef CS_THERMAL_MODEL_H
#define CS_THERMAL_MODEL_H
 
/*============================================================================
 * Base thermal model data.
 *============================================================================*/
 
/*
  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_defs.h"
#include "base/cs_field.h"
 
/*=============================================================================
 * Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Thermal model type
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_THERMAL_MODEL_INIT = -999,
  CS_THERMAL_MODEL_NONE = 0,
  CS_THERMAL_MODEL_TEMPERATURE,
  CS_THERMAL_MODEL_ENTHALPY,
  CS_THERMAL_MODEL_TOTAL_ENERGY,
  CS_THERMAL_MODEL_INTERNAL_ENERGY,
  CS_THERMAL_MODEL_N_TYPES
 
} cs_thermal_model_variable_t;
 
typedef enum {
 
  CS_TEMPERATURE_SCALE_NONE = 0,
  CS_TEMPERATURE_SCALE_KELVIN = 1,
  CS_TEMPERATURE_SCALE_CELSIUS = 2
 
} cs_temperature_scale_t;
 
/* thermal model descriptor */
/*--------------------------*/
 
typedef struct {
 
  cs_thermal_model_variable_t  thermal_variable;   /* Thermal variable */
  cs_temperature_scale_t       temperature_scale;  /* Temperature scale */
 
  /* Has kinetic source term correction */
  int           has_kinetic_st;
  int           cflt;           /* compute the thermal cfl condition */
  int           cflp;           /* compute the pressure cfl condition */
  bool          has_pdivu;
  bool          has_dissipation;
 
} cs_thermal_model_t;
 
/*============================================================================
 * Static global variables
 *============================================================================*/
 
/* Pointer to thermal model structure */
 
extern const cs_thermal_model_t  *cs_glob_thermal_model;
 
/*=============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Provide access to cs_glob_thermal_model
 *
 * needed to initialize structure with GUI
 *----------------------------------------------------------------------------*/
 
cs_thermal_model_t *
cs_get_glob_thermal_model(void);
 
/*----------------------------------------------------------------------------
 * Return thermal field (temperature, enthalpy, total energy according to
 * thermal model).
 *
 * returns:
 *   pointer to thermal field
 *----------------------------------------------------------------------------*/
 
cs_field_t *
cs_thermal_model_field(void);
 
/*----------------------------------------------------------------------------
 * Print the thermal model structure to setup.log.
 *----------------------------------------------------------------------------*/
 
void
cs_thermal_model_log_setup(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_init(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_gamma_d_c_square(const cs_real_t  cp[],
                                  const cs_real_t  cv[],
                                  const cs_real_t  temp[],
                                  const cs_real_t  pres[],
                                  const cs_real_t  fracv[],
                                  const cs_real_t  fracm[],
                                  const cs_real_t  frace[],
                                  cs_real_t        gdc2[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
CS_F_HOST_DEVICE cs_real_t
cs_thermal_model_demdt(const cs_real_t  pres,
                       const cs_real_t  temp,
                       const cs_real_t  yw,
                       const cs_real_t  rvsra,
                       const cs_real_t  cva,
                       const cs_real_t  cvv,
                       const cs_real_t  cpl,
                       const cs_real_t  l00);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
CS_F_HOST_DEVICE cs_real_t
cs_thermal_model_demdt_ecsnt(const cs_real_t  pres,
                             const cs_real_t  temp,
                             const cs_real_t  yw,
                             const cs_real_t  rvsra,
                             const cs_real_t  cva,
                             const cs_real_t  cvv,
                             const cs_real_t  cpl,
                             const cs_real_t  l00);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_kinetic_st_prepare(const cs_real_t  imasfl[],
                                    const cs_real_t  bmasfl[],
                                    const cs_real_t  vela[][3],
                                    const cs_real_t  vel[][3]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_kinetic_st_finalize(const cs_real_t  cromk1[],
                                     const cs_real_t  cromk[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_add_kst(cs_real_t  smbrs[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_cflp(const cs_real_t  croma[],
                      const cs_real_t  trav2[][3],
                      const cs_real_t  cvara_pr[],
                      const cs_real_t  imasfl[],
                      cs_real_t        cflp[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_cflt(const cs_real_t  croma[],
                      const cs_real_t  tempk[],
                      const cs_real_t  tempka[],
                      const cs_real_t  xcvv[],
                      const cs_real_t  vel[][3],
                      const cs_real_t  i_massflux[],
                      const cs_real_t  b_massflux[],
                      cs_real_t        cflt[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_eos_predicted_rho(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_cv(cs_real_t  *xcvv);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_dissipation(const cs_real_t  vistot[],
                             const cs_real_t  gradv[][3][3],
                             cs_real_t        smbrs[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_newton_t(int               method,
                          const cs_real_t  *pk1,
                          const cs_real_t   th_scal[],
                          const cs_real_t   cvar_pr[],
                          const cs_real_t   cvara_pr[],
                          const cs_real_t   yw[],
                          cs_real_t         yv[],
                          cs_real_t         temp[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_thermal_model_pdivu(cs_real_t  smbrs[]);
 
/*----------------------------------------------------------------------------*/
 
#endif /* CS_THERMAL_MODEL_H */