8.0
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 
62 } cs_air_fluid_props_t;
63 
64 extern cs_air_fluid_props_t *cs_glob_air_props;
65 
66 /*============================================================================
67  * Public function prototypes for Fortran API
68  *============================================================================*/
69 
70 /*============================================================================
71  * Prototypes of public function
72  *============================================================================*/
73 
74 /*----------------------------------------------------------------------------*/
82 /*----------------------------------------------------------------------------*/
83 
84 cs_real_t
85 cs_air_h_l(cs_real_t t_l);
86 
87 /*----------------------------------------------------------------------------*/
95 /*----------------------------------------------------------------------------*/
96 
97 cs_real_t
98 cs_air_hvap(cs_real_t t_vap);
99 
100 /*----------------------------------------------------------------------------*/
109 /*----------------------------------------------------------------------------*/
110 
111 cs_real_t
112 cs_air_dxsath(cs_real_t th,
113  cs_real_t p0);
114 
115 /*----------------------------------------------------------------------------*/
124 /*----------------------------------------------------------------------------*/
125 
126 cs_real_t
127 cs_air_cp_humidair(cs_real_t x,
128  cs_real_t x_s);
129 
130 /*----------------------------------------------------------------------------*/
141 /*----------------------------------------------------------------------------*/
142 
143 cs_real_t
144 cs_air_h_humidair(cs_real_t cp_h,
145  cs_real_t x,
146  cs_real_t x_s,
147  cs_real_t t_h);
148 
149 /*----------------------------------------------------------------------------*/
160 /*----------------------------------------------------------------------------*/
161 
162 cs_real_t
163 cs_air_t_humidair(cs_real_t cp_h,
164  cs_real_t x,
165  cs_real_t x_s,
166  cs_real_t h_h);
167 
168 /*----------------------------------------------------------------------------*/
176 /*----------------------------------------------------------------------------*/
177 
178 cs_real_t
179 cs_liq_h_to_t(cs_real_t h_l);
180 
181 /*----------------------------------------------------------------------------*/
189 /*----------------------------------------------------------------------------*/
190 
191 cs_real_t
192 cs_liq_t_to_h(cs_real_t t_l);
193 
194 /*----------------------------------------------------------------------------*/
195 
196 /*----------------------------------------------------------------------------*/
205 /*----------------------------------------------------------------------------*/
206 
207 cs_real_t
208 cs_air_x_sat(cs_real_t t_c,
209  cs_real_t p);
210 
211 /*----------------------------------------------------------------------------*/
221 /*----------------------------------------------------------------------------*/
222 
223 cs_real_t
224 cs_air_yw_sat(cs_real_t t_c,
225  cs_real_t p);
226 
227 /*----------------------------------------------------------------------------*/
236 /*----------------------------------------------------------------------------*/
237 
238 cs_real_t
239 cs_air_pwv_sat(cs_real_t t_c);
240 
241 /*----------------------------------------------------------------------------*/
250 /*----------------------------------------------------------------------------*/
251 
252 cs_real_t
253 cs_air_x_to_yw(const cs_real_t x);
254 
255 /*----------------------------------------------------------------------------*/
264 /*----------------------------------------------------------------------------*/
265 
266 cs_real_t
267 cs_air_yw_to_x(cs_real_t qw);
268 
269 /*----------------------------------------------------------------------------*/
281 /*----------------------------------------------------------------------------*/
282 
283 void
284 cs_rho_humidair(cs_real_t ywm, //TODO rename yw_h
285  cs_real_t t_liq,
286  cs_real_t p,
287  cs_real_t *yw_liq,
288  cs_real_t *t_h,
289  cs_real_t *rho_h);
290 
291 /*----------------------------------------------------------------------------*/
304 /*----------------------------------------------------------------------------*/
305 
306 cs_real_t
307 cs_air_rho_humidair(cs_real_t x,
308  cs_real_t rho0,
309  cs_real_t p0,
310  cs_real_t t0,
311  cs_real_t molmassrat,
312  cs_real_t t_h);
313 
314 /*----------------------------------------------------------------------------*/
315 
316 END_C_DECLS
317 
318 #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:371
cs_glob_air_props
cs_air_fluid_props_t * cs_glob_air_props
Definition: cs_air_props.c:85
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:453
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:413
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:279
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:392
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:553
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:144
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:663
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:574
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:599
cs_air_hvap
cs_real_t cs_air_hvap(cs_real_t t_vap)
Calculation water vapor mass enthalpy.
Definition: cs_air_props.c:162
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:186
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:478
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:343
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:310
cs_base.h
BEGIN_C_DECLS
#define BEGIN_C_DECLS
Definition: cs_defs.h:509
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:510
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:167
cstphy::t0
real(c_double), pointer, save t0
reference temperature.
Definition: cstphy.f90:214
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::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::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