1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
#ifndef __CS_LAGR_ROUGHNESS_H__
#define __CS_LAGR_ROUGHNESS_H__
/*============================================================================
* Surface roughness
*============================================================================*/
/*
This file is part of Code_Saturne, a general-purpose CFD tool.
Copyright (C) 1998-2016 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 "cs_defs.h"
#include "cs_lagr_particle.h"
/*----------------------------------------------------------------------------*/
BEGIN_C_DECLS
/*============================================================================
* Type definitions
*============================================================================*/
typedef struct {
cs_real_t water_permit;
cs_real_t ionic_strength;
cs_real_t phi_p;
cs_real_t phi_s;
cs_real_t *temperature;
cs_real_t valen;
cs_real_t *debye_length;
cs_real_t cstham;
cs_real_t lambda_vdw;
cs_real_t espasg;
cs_real_t denasp;
cs_real_t rayasp;
cs_real_t rayasg;
} cs_lagr_roughness_param_t;
extern cs_lagr_roughness_param_t *cs_lagr_roughness_param;
/*============================================================================
* Public function definitions for Fortran API
*============================================================================*/
/*----------------------------------------------------------------------------
* Clogging initialization:
* - Retrieve various parameters for storing in global structure
* - Compute and store the Debye screening length
*----------------------------------------------------------------------------*/
void
roughness_init (const cs_real_t *water_permit,
const cs_real_t *ionic_strength,
const cs_real_t temperature[],
const cs_real_t *valen,
const cs_real_t *phi_p,
const cs_real_t *phi_s,
const cs_real_t *cstham,
const cs_real_t *lambda_vdw,
const cs_real_t *espasg,
const cs_real_t *denasp,
const cs_real_t *rayasp,
const cs_real_t *rayasg
);
/*=============================================================================
* Function definitions
*============================================================================*/
/*----------------------------------------------------------------------------
* Deallocate the arrays storing temperature and Debye length.
*----------------------------------------------------------------------------*/
void
cs_lagr_roughness_finalize(void);
/*----------------------------------------------------------------------------
* Compute the energy barrier for a rough wall.
*
* parameters:
* particle <-- pointer to particle data
* attr_map <-- pointer to attribute map
* iel <-- id of cell where the particle is
* energy_barrier <-> energy barrier
*----------------------------------------------------------------------------*/
void
cs_lagr_roughness_barrier(const void *particle,
const cs_lagr_attribute_map_t *attr_map,
cs_lnum_t iel,
cs_real_t *energy_barrier);
/*----------------------------------------------------------------------------*/
END_C_DECLS
#endif /* __CS_LAGR_ROUGHNESS_H__ */
|
