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
121
122
123
124
125
126
127
128
129
130
131
132
|
#ifndef __CS_LAGR_CLOGGING_H__
#define __CS_LAGR_CLOGGING_H__
/*============================================================================
* Clogging modeling.
*============================================================================*/
/*
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 jamming_limit;
cs_real_t min_porosity;
cs_real_t valen;
cs_real_t phi_p;
cs_real_t phi_s;
cs_real_t* temperature;
cs_real_t* debye_length;
cs_real_t cstham;
cs_real_t csthpp;
cs_real_t lambda_vdw;
} cs_lagr_clogging_param_t;
/*============================================================================
* Public function definitions
*============================================================================*/
/*----------------------------------------------------------------------------
* Clogging initialization.
*
* - Retrieve various parameters for storing in global structure.
* - Compute and store the Debye screening length
*----------------------------------------------------------------------------*/
void
cloginit(const cs_real_t *water_permit,
const cs_real_t *ionic_strength,
const cs_real_t *jamming_limit,
const cs_real_t *min_porosity,
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 *csthpp,
const cs_real_t *lambda_vdw
);
/*----------------------------------------------------------------------------
* Clogging finalization.
*
* Deallocate the arrays storing temperature and Debye length.
*----------------------------------------------------------------------------*/
void
cs_lagr_clogging_finalize(void);
/*----------------------------------------------------------------------------
* Clogging:
*
* - Compute the number of deposited particles in contact with the depositing
* particle
* - Re-compute the energy barrier if this number is greater than zero
*
* parameters:
* particle <-- pointer to particle data
* attr_map <-- pointer to attribute map
* iel <-- id of cell where the particle is
* face_area <-- area of face
* energy_barrier <-> energy barrier
* surface_coverage <-> surface coverage
* limit <-> jamming limit
* mporos <-> minimum porosity
*
* returns:
* number of deposited particles in contact with the depositing particle
*----------------------------------------------------------------------------*/
int
cs_lagr_clogging_barrier(const void *particle,
const cs_lagr_attribute_map_t *attr_map,
cs_lnum_t iel,
cs_real_t *energy_barrier,
cs_real_t *surface_coverage,
cs_real_t *limit,
cs_real_t *mporos);
/*----------------------------------------------------------------------------*/
END_C_DECLS
#endif /* __CS_LAGR_CLOGGING_H__ */
|
