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
|
/* This file is part of the OpenLB library
*
* Copyright (C) 2014 Thomas Henn, Mathias J. Krause
* E-mail contact: info@openlb.net
* The most recent release of OpenLB can be downloaded at
* <http://www.openlb.net/>
*
* 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.
*/
#ifndef MATERIALBOUNDARY3D_H
#define MATERIALBOUNDARY3D_H
#include <set>
#include "geometry/superGeometry3D.h"
#include "particles/particleSystem3D.h"
#include "boundary3D.h"
namespace olb {
template<typename T, template<typename U> class PARTICLETYPE>
class ParticleSystem3D;
/*
* Particle boundary based on the fluids material number.
* If a particle moves in the vincinity of a lattice node with
* specified material number it is set inactive and its velocity is set to 0.
**/
template<typename T, template<typename U> class PARTICLETYPE>
class MaterialBoundary3D : public Boundary3D<T, PARTICLETYPE> {
public:
/// Constructor
MaterialBoundary3D(SuperGeometry3D<T>& sg);
/// Copy constructor
MaterialBoundary3D(MaterialBoundary3D<T, PARTICLETYPE>& f);
/// Constructor with set of material numbers
MaterialBoundary3D(SuperGeometry3D<T>& sg,
std::set<int> material);
~MaterialBoundary3D() override
{
}
/// Add a single material number
void addMaterial(int mat)
{
_materials.insert(mat);
}
/// Add several material numbers
void addMaterial(std::vector<int> mats)
{
for (unsigned i=0; i< mats.size(); ++i) {
_materials.insert(mats[i]);
}
}
/// Apply the boundary condition
void applyBoundary(
typename std::deque<PARTICLETYPE<T> >::iterator& p,
ParticleSystem3D<T, PARTICLETYPE>& psSys) override;
private:
SuperGeometry3D<T>& _sg;
std::set<int> _materials;
std::set<int>::iterator _matIter;
// int x0, y0, z0;
};
}
#endif /* MATERIALBOUNDARY3D_H */
|
