8.1
general documentation
cs_air_props.h
Go to the documentation of this file.
1 #ifndef __CS_AIR_PROPS_H__
2 #define __CS_AIR_PROPS_H__
3 
4 /*============================================================================
5  * Specific laws for air properties (temperature, enthalpy)
6  *============================================================================*/
7 
8 /*
9  This file is part of code_saturne, a general-purpose CFD tool.
10 
11  Copyright (C) 1998-2023 EDF S.A.
12 
13  This program is free software; you can redistribute it and/or modify it under
14  the terms of the GNU General Public License as published by the Free Software
15  Foundation; either version 2 of the License, or (at your option) any later
16  version.
17 
18  This program is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
20  FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
21  details.
22 
23  You should have received a copy of the GNU General Public License along with
24  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
25  Street, Fifth Floor, Boston, MA 02110-1301, USA.
26 */
27 
28 /*----------------------------------------------------------------------------*/
29 
30 /*----------------------------------------------------------------------------
31  * Standard C library headers
32  *----------------------------------------------------------------------------*/
33 
34 /*----------------------------------------------------------------------------
35  * Local headers
36  *----------------------------------------------------------------------------*/
37 
38 #include "cs_base.h"
39 
40 /*----------------------------------------------------------------------------*/
41 
42 BEGIN_C_DECLS
43 
44 /*============================================================================
45  * Structure definition
46  *============================================================================*/
47 
48 /* Structure associated to general properties */
49 
50 typedef struct {
51 
52  cs_real_t humidity0; /* Reference humidity */
53  cs_real_t cp_a; /* Specific heat of dry air */
54  cs_real_t cp_v; /* Specific heat of vapor */
55  cs_real_t cp_l; /* Specific heat of liquid water */
56  cs_real_t hv0; /* Chaleur latente */
57  cs_real_t rho_l; /* Masse volumique de l eau */
58  cs_real_t lambda_h; /* Humid air conductivity */
59  cs_real_t lambda_l; /* Water conductivity */
60  cs_real_t droplet_diam; /* Drop diameter for rain zones */
61  cs_real_t molmass_rat; /* Ratio of the molar mass (H2O) over the
62  molar mass (air) */
63  cs_real_t sigma; /* Surface tension between water and air */
64 
65 } cs_air_fluid_props_t;
66 
67 extern cs_air_fluid_props_t *cs_glob_air_props;
68 
69 /*============================================================================
70  * Public function prototypes for Fortran API
71  *============================================================================*/
72 
73 /*============================================================================
74  * Prototypes of public function
75  *============================================================================*/
76 
77 /*----------------------------------------------------------------------------*/
85 /*----------------------------------------------------------------------------*/
86 
87 cs_real_t
88 cs_air_h_l(cs_real_t t_l);
89 
90 /*----------------------------------------------------------------------------*/
98 /*----------------------------------------------------------------------------*/
99 
100 cs_real_t
101 cs_air_hvap(cs_real_t t_vap);
102 
103 /*----------------------------------------------------------------------------*/
112 /*----------------------------------------------------------------------------*/
113 
114 cs_real_t
115 cs_air_dxsath(cs_real_t th,
116  cs_real_t p0);
117 
118 /*----------------------------------------------------------------------------*/
127 /*----------------------------------------------------------------------------*/
128 
129 cs_real_t
130 cs_air_cp_humidair(cs_real_t x,
131  cs_real_t x_s);
132 
133 /*----------------------------------------------------------------------------*/
144 /*----------------------------------------------------------------------------*/
145 
146 cs_real_t
147 cs_air_h_humidair(cs_real_t cp_h,
148  cs_real_t x,
149  cs_real_t x_s,
150  cs_real_t t_h);
151 
152 /*----------------------------------------------------------------------------*/
163 /*----------------------------------------------------------------------------*/
164 
165 cs_real_t
166 cs_air_t_humidair(cs_real_t cp_h,
167  cs_real_t x,
168  cs_real_t x_s,
169  cs_real_t h_h);
170 
171 /*----------------------------------------------------------------------------*/
179 /*----------------------------------------------------------------------------*/
180 
181 cs_real_t
182 cs_liq_h_to_t(cs_real_t h_l);
183 
184 /*----------------------------------------------------------------------------*/
192 /*----------------------------------------------------------------------------*/
193 
194 cs_real_t
195 cs_liq_t_to_h(cs_real_t t_l);
196 
197 /*----------------------------------------------------------------------------*/
198 
199 /*----------------------------------------------------------------------------*/
208 /*----------------------------------------------------------------------------*/
209 
210 cs_real_t
211 cs_air_x_sat(cs_real_t t_c,
212  cs_real_t p);
213 
214 /*----------------------------------------------------------------------------*/
224 /*----------------------------------------------------------------------------*/
225 
226 cs_real_t
227 cs_air_yw_sat(cs_real_t t_c,
228  cs_real_t p);
229 
230 /*----------------------------------------------------------------------------*/
239 /*----------------------------------------------------------------------------*/
240 
241 cs_real_t
242 cs_air_pwv_sat(cs_real_t t_c);
243 
244 /*----------------------------------------------------------------------------*/
253 /*----------------------------------------------------------------------------*/
254 
255 cs_real_t
256 cs_air_x_to_yw(const cs_real_t x);
257 
258 /*----------------------------------------------------------------------------*/
267 /*----------------------------------------------------------------------------*/
268 
269 cs_real_t
270 cs_air_yw_to_x(cs_real_t qw);
271 
272 /*----------------------------------------------------------------------------*/
284 /*----------------------------------------------------------------------------*/
285 
286 void
287 cs_rho_humidair(cs_real_t ywm, //TODO rename yw_h
288  cs_real_t t_liq,
289  cs_real_t p,
290  cs_real_t *yw_liq,
291  cs_real_t *t_h,
292  cs_real_t *rho_h);
293 
294 /*----------------------------------------------------------------------------*/
307 /*----------------------------------------------------------------------------*/
308 
309 cs_real_t
310 cs_air_rho_humidair(cs_real_t x,
311  cs_real_t rho0,
312  cs_real_t p0,
313  cs_real_t t0,
314  cs_real_t molmassrat,
315  cs_real_t t_h);
316 
317 /*----------------------------------------------------------------------------*/
318 
319 END_C_DECLS
320 
321 #endif /* __CS_AIR_PROPERTIES_H__ */
cs_liq_h_to_t
cs_real_t cs_liq_h_to_t(cs_real_t h_l)
Calculation of the temperature of liquid water.
Definition: cs_air_props.c:352
cs_glob_air_props
cs_air_fluid_props_t * cs_glob_air_props
Definition: cs_air_props.c:87
cs_air_yw_sat
cs_real_t cs_air_yw_sat(cs_real_t t_c, cs_real_t p)
Calculation of the air water mass fraction at saturation for a given temperature.
Definition: cs_air_props.c:434
cs_air_x_sat
cs_real_t cs_air_x_sat(cs_real_t t_c, cs_real_t p)
Calculation of the air humidity at saturation for a given temperature.
Definition: cs_air_props.c:394
cs_air_cp_humidair
cs_real_t cs_air_cp_humidair(cs_real_t x, cs_real_t x_s)
Calculation of the Cp of humid air.
Definition: cs_air_props.c:260
cs_liq_t_to_h
cs_real_t cs_liq_t_to_h(cs_real_t t_l)
Calculation of the specific enthalpy of liquid water.
Definition: cs_air_props.c:373
cs_air_x_to_yw
cs_real_t cs_air_x_to_yw(const cs_real_t x)
Convert the absolute humidity of humid air to the air water mass fraction qwt = Ym = mw/mh.
Definition: cs_air_props.c:534
cs_air_h_l
cs_real_t cs_air_h_l(cs_real_t t_l)
Calculation water mass enthalpy.
Definition: cs_air_props.c:125
cs_air_rho_humidair
cs_real_t cs_air_rho_humidair(cs_real_t x, cs_real_t rho0, cs_real_t p0, cs_real_t t0, cs_real_t molmassrat, cs_real_t t_h)
Calculation of the density of humid air.
Definition: cs_air_props.c:644
cs_air_yw_to_x
cs_real_t cs_air_yw_to_x(cs_real_t qw)
Convert the air water mass fraction qwt = Ym = mw/mh to the absolute humidity of humid air.
Definition: cs_air_props.c:555
cs_rho_humidair
void cs_rho_humidair(cs_real_t ywm, cs_real_t t_liq, cs_real_t p, cs_real_t *yw_liq, cs_real_t *t_h, cs_real_t *rho_h)
Calculation of the density of humid air.
Definition: cs_air_props.c:580
cs_air_hvap
cs_real_t cs_air_hvap(cs_real_t t_vap)
Calculation water vapor mass enthalpy.
Definition: cs_air_props.c:143
cs_air_dxsath
cs_real_t cs_air_dxsath(cs_real_t th, cs_real_t p0)
Calculation of the derivate of the absolute humidity at saturation.
Definition: cs_air_props.c:167
cs_air_pwv_sat
cs_real_t cs_air_pwv_sat(cs_real_t t_c)
Computes the saturation water vapor pressure function of the temperature (in Celsius)
Definition: cs_air_props.c:459
cs_air_t_humidair
cs_real_t cs_air_t_humidair(cs_real_t cp_h, cs_real_t x, cs_real_t x_s, cs_real_t h_h)
Calculation of the temperature of humid air.
Definition: cs_air_props.c:324
cs_air_h_humidair
cs_real_t cs_air_h_humidair(cs_real_t cp_h, cs_real_t x, cs_real_t x_s, cs_real_t t_h)
Calculation of the specific enthalpy of humid air.
Definition: cs_air_props.c:291
cs_base.h
BEGIN_C_DECLS
#define BEGIN_C_DECLS
Definition: cs_defs.h:514
cs_real_t
double cs_real_t
Floating-point value.
Definition: cs_defs.h:319
END_C_DECLS
#define END_C_DECLS
Definition: cs_defs.h:515
h_l
@ h_l
Definition: cs_field_pointer.h:202
p
@ p
Definition: cs_field_pointer.h:67
t_l
@ t_l
Definition: cs_field_pointer.h:201
atimbr::qw
double precision, dimension(:,:,:), allocatable qw
Definition: atimbr.f90:109
cstphy::p0
real(c_double), pointer, save p0
reference pressure for the total pressure.
Definition: cstphy.f90:170
cstphy::t0
real(c_double), pointer, save t0
reference temperature.
Definition: cstphy.f90:217
ppthch::th
double precision, dimension(npot), save th
temperature (in K)
Definition: ppthch.f90:119
cs_air_fluid_props_t
Definition: cs_air_props.h:50
cs_air_fluid_props_t::cp_l
cs_real_t cp_l
Definition: cs_air_props.h:55
cs_air_fluid_props_t::cp_v
cs_real_t cp_v
Definition: cs_air_props.h:54
cs_air_fluid_props_t::humidity0
cs_real_t humidity0
Definition: cs_air_props.h:52
cs_air_fluid_props_t::sigma
cs_real_t sigma
Definition: cs_air_props.h:63
cs_air_fluid_props_t::droplet_diam
cs_real_t droplet_diam
Definition: cs_air_props.h:60
cs_air_fluid_props_t::lambda_l
cs_real_t lambda_l
Definition: cs_air_props.h:59
cs_air_fluid_props_t::molmass_rat
cs_real_t molmass_rat
Definition: cs_air_props.h:61
cs_air_fluid_props_t::rho_l
cs_real_t rho_l
Definition: cs_air_props.h:57
cs_air_fluid_props_t::cp_a
cs_real_t cp_a
Definition: cs_air_props.h:53
cs_air_fluid_props_t::lambda_h
cs_real_t lambda_h
Definition: cs_air_props.h:58
cs_air_fluid_props_t::hv0
cs_real_t hv0
Definition: cs_air_props.h:56