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
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
|
/* This file is part of the OpenLB library
*
* Copyright (C) 2006, 2016 Robin Trunk
* E-mail contact: info@openlb.net
* The most recent release of OpenLB can be downloaded at
* <http://www.openlb.net/>
*
* Generic version of the collision, which modifies the particle
* distribution functions, by Orestis Malaspinas.
*
* 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.
*/
#include "core/postProcessing.h"
#include "core/blockLattice3D.h"
namespace olb {
/**
* This class interpolates missing f_i from values
* near the boundary to get a more stable outflow
* condition for the density. It is assumed that the
* next two cells in direction of this f_i have
* viable values.
*/
template<typename T, typename DESCRIPTOR>
class ConvectionBoundaryProcessor3D : public LocalPostProcessor3D<T,DESCRIPTOR> {
public:
ConvectionBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
int z1_, int discreteNormalX_,
int discreteNormalY_, int discreteNormalZ_);
int extent() const override
{
return 0;
}
int extent(int whichDirection) const override
{
return 0;
}
void process(BlockLattice3D<T,DESCRIPTOR>& blockLattice) override;
void processSubDomain ( BlockLattice3D<T,DESCRIPTOR>& blockLattice,
int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
private:
int interpolationPop[DESCRIPTOR::q];
int x0, x1, y0, y1, z0, z1;
};
template<typename T, typename DESCRIPTOR>
class ConvectionBoundaryProcessorGenerator3D : public PostProcessorGenerator3D<T,DESCRIPTOR> {
public:
ConvectionBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_,
int z0_, int z1_, int discreteNormalX_,
int discreteNormalY_, int discreteNormalZ_);
PostProcessor3D<T,DESCRIPTOR>* generate() const override;
PostProcessorGenerator3D<T,DESCRIPTOR>* clone() const override;
private:
int discreteNormalX;
int discreteNormalY;
int discreteNormalZ;
};
/**
* This class copies missing values in the
* external field from the neighbour in normal direction.
* Therefore it is assumed this neighbour is a fluid cell.
*/
template<typename T, typename DESCRIPTOR>
class ExtFieldBoundaryProcessor3D : public LocalPostProcessor3D<T,DESCRIPTOR> {
public:
ExtFieldBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
int z1_, int discreteNormalX_, int discreteNormalY_,
int discreteNormalZ_, int offset_);
int extent() const override
{
return 0;
}
int extent(int whichDirection) const override
{
return 0;
}
void process(BlockLattice3D<T,DESCRIPTOR>& blockLattice) override;
void processSubDomain ( BlockLattice3D<T,DESCRIPTOR>& blockLattice,
int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
private:
int x0, x1, y0, y1, z0, z1;
int discreteNormalX, discreteNormalY, discreteNormalZ;
int offset;
bool par;
};
template<typename T, typename DESCRIPTOR>
class ExtFieldBoundaryProcessorGenerator3D : public PostProcessorGenerator3D<T,DESCRIPTOR> {
public:
ExtFieldBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_, int z0_,
int z1_, int discreteNormalX_, int discreteNormalY_,
int discreteNormalZ_, int offset_);
PostProcessor3D<T,DESCRIPTOR>* generate() const override;
PostProcessorGenerator3D<T,DESCRIPTOR>* clone() const override;
private:
int discreteNormalX, discreteNormalY, discreteNormalZ;
int offset;
};
/**
* This class resets some values of the distribution
* on the boundary that can have arbitrary values
* to be zero and thus ensures a correct computation
* of the density that is about to leave the domain.
*/
template<typename T, typename DESCRIPTOR>
class ZeroDistributionBoundaryProcessor3D : public LocalPostProcessor3D<T,DESCRIPTOR> {
public:
ZeroDistributionBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
int z1_, int discreteNormalX_, int discreteNormalY_,
int discreteNormalZ_);
int extent() const override
{
return 0;
}
int extent(int whichDirection) const override
{
return 0;
}
void process(BlockLattice3D<T,DESCRIPTOR>& blockLattice) override;
void processSubDomain ( BlockLattice3D<T,DESCRIPTOR>& blockLattice,
int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
private:
int resetPop[DESCRIPTOR::q];
int x0, x1, y0, y1, z0, z1;
};
template<typename T, typename DESCRIPTOR>
class ZeroDistributionBoundaryProcessorGenerator3D : public PostProcessorGenerator3D<T,DESCRIPTOR> {
public:
ZeroDistributionBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_,
int z0_, int z1_, int discreteNormalX_,
int discreteNormalY_, int discreteNormalZ_);
PostProcessor3D<T,DESCRIPTOR>* generate() const override;
PostProcessorGenerator3D<T,DESCRIPTOR>* clone() const override;
private:
int discreteNormalX;
int discreteNormalY;
int discreteNormalZ;
};
}
|
