cs_elec_model.h

Go to the documentation of this file.

#ifndef CS_ELEC_MODEL_H
#define CS_ELEC_MODEL_H
 
/*============================================================================
 * General parameters management.
 *============================================================================*/
 
/*
  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 "cdo/cs_domain.h"
#include "mesh/cs_mesh.h"
#include "mesh/cs_mesh_quantities.h"
 
/*=============================================================================
 * Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Structure to read properties in dp_ELE
 *----------------------------------------------------------------------------*/
 
typedef struct {
 
  int         n_gas;
  int         n_point;
  cs_real_t  *th;
  cs_real_t  *eh_gas;
  cs_real_t  *rhoel;
  cs_real_t  *cpel;
  cs_real_t  *sigel;
  cs_real_t  *visel;
  cs_real_t  *xlabel;
  cs_real_t  *xkabel;
 
} cs_data_elec_t;
 
/*----------------------------------------------------------------------------
 * Structure to read transformer parameters in dp_ELE
 *----------------------------------------------------------------------------*/
 
typedef struct {
 
  int         nbelec;
  int        *ielecc;
  int        *ielect;
  int        *ielecb;
  int         nbtrf;
  int         ntfref;
  int        *ibrpr;
  int        *ibrsec;
  cs_real_t  *tenspr;
  cs_real_t  *rnbs;
  cs_real_t  *zr;
  cs_real_t  *zi;
  cs_real_t  *uroff;
  cs_real_t  *uioff;
 
} cs_data_joule_effect_t;
 
/*----------------------------------------------------------------------------
 * Electrical model options descriptor
 *----------------------------------------------------------------------------*/
 
typedef struct {
 
  int         ixkabe;
  int         ntdcla;
  int         irestrike;
  cs_real_t   restrike_point[3];
  cs_real_t   crit_reca[5];
 
  int         ielcor;
  int         modrec;
  int         idreca;
  int        *izreca;
  cs_real_t   couimp;
  cs_real_t   pot_diff;
  cs_real_t   puisim;
  cs_real_t   coejou;
  cs_real_t   elcou;
  cs_real_t   srrom;
 
} cs_elec_option_t;
 
/*============================================================================
 * Static global variables
 *============================================================================*/
 
/* Pointer to electrical model options structure */
 
extern const cs_elec_option_t        *cs_glob_elec_option;
extern const cs_data_elec_t          *cs_glob_elec_properties;
extern const cs_data_joule_effect_t  *cs_glob_transformer;
 
/* Constant for electrical models */
 
extern const cs_real_t cs_elec_permvi;
extern const cs_real_t cs_elec_epszer;
 
/*=============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Provide access to cs_elec_option
 *----------------------------------------------------------------------------*/
 
cs_elec_option_t *
cs_get_glob_elec_option(void);
 
/*----------------------------------------------------------------------------
 * Provide access to cs_glob_transformer
 *----------------------------------------------------------------------------*/
 
cs_data_joule_effect_t *
cs_get_glob_transformer(void);
 
/*----------------------------------------------------------------------------
 * Initialize structures for electrical model
 *----------------------------------------------------------------------------*/
 
void
cs_electrical_model_initialize(void);
 
/*----------------------------------------------------------------------------
 * Destroy structures for electrical model
 *----------------------------------------------------------------------------*/
 
void
cs_electrical_model_finalize(void);
 
/*----------------------------------------------------------------------------
 * Specific initialization for electric arc
 *----------------------------------------------------------------------------*/
 
void
cs_electrical_model_specific_initialization(void);
 
/*----------------------------------------------------------------------------
 * Read properties file
 *----------------------------------------------------------------------------*/
 
void
cs_electrical_properties_read(void);
 
/*----------------------------------------------------------------------------
 * compute specific electric arc fields
 *----------------------------------------------------------------------------*/
 
void
cs_elec_compute_fields(const cs_mesh_t  *mesh,
                       int               call_id);
 
/*----------------------------------------------------------------------------
 * compute physical properties
 *----------------------------------------------------------------------------*/
 
void
cs_elec_physical_properties(cs_domain_t             *domain);
 
/*----------------------------------------------------------------------------
 * compute source terms for energy
 *----------------------------------------------------------------------------*/
 
void
cs_elec_source_terms(const cs_mesh_t             *mesh,
                     const cs_mesh_quantities_t  *mesh_quantities,
                     int                          f_id,
                     cs_real_t                   *smbrs);
 
/*----------------------------------------------------------------------------
 * compute source terms for vector potential
 *----------------------------------------------------------------------------*/
 
void
cs_elec_source_terms_v(const cs_mesh_t             *mesh,
                       const cs_mesh_quantities_t  *mesh_quantities,
                       int                          f_id,
                       cs_real_3_t                 *smbrv);
 
/*----------------------------------------------------------------------------
 * add variables fields
 *----------------------------------------------------------------------------*/
 
void
cs_elec_add_variable_fields(void);
 
/*----------------------------------------------------------------------------
 * add properties fields
 *----------------------------------------------------------------------------*/
 
void
cs_elec_add_property_fields(void);
 
/*----------------------------------------------------------------------------
 * initialize electric fields
 *----------------------------------------------------------------------------*/
 
void
cs_elec_fields_initialize(const cs_mesh_t  *mesh);
 
/*----------------------------------------------------------------------------
 * scaling electric quantities
 *----------------------------------------------------------------------------*/
 
void
cs_elec_scaling_function(const cs_mesh_t             *mesh,
                         const cs_mesh_quantities_t  *mesh_quantities,
                         cs_real_t                   *dt);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_elec_convert_h_to_t_faces(const cs_real_t  h[],
                             cs_real_t        t[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_elec_convert_h_to_t(const cs_real_t  ym[],
                       cs_real_t        enthal);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_elec_convert_t_to_h_cells(const cs_real_t  t[],
                             cs_real_t        h[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_elec_convert_t_to_h_faces(const cs_lnum_t  n_faces,
                             const cs_lnum_t  face_ids[],
                             const cs_real_t  t[],
                             cs_real_t        h[]);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_elec_convert_t_to_h(const cs_real_t ym[],
                       cs_real_t       temp);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_elec_define_functions(void);
 
/*----------------------------------------------------------------------------*/
 
#endif /* CS_ELEC_MODEL_H */