OpenLB 1.7
Loading...
Searching...
No Matches
setInterpolatedConvectionBoundary2D.hh
Go to the documentation of this file.
1/* This file is part of the OpenLB library
2 *
3 * Copyright (C) 2020 Alexander Schulz
4 * E-mail contact: info@openlb.net
5 * The most recent release of OpenLB can be downloaded at
6 * <http://www.openlb.net/>
7 *
8 * This program is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU General Public License
10 * as published by the Free Software Foundation; either version 2
11 * of the License, or (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the GNU General Public
19 * License along with this program; if not, write to the Free
20 * Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
21 * Boston, MA 02110-1301, USA.
22*/
23
24//This file contains the Interpolated Convection Boundary
25//This is a new version of the Boundary, which only contains free floating functions
26#ifndef SET_INTERPOLATED_CONVECTION_BOUNDARY_2D_HH
27#define SET_INTERPOLATED_CONVECTION_BOUNDARY_2D_HH
28
30
31namespace olb {
32
34template<typename T, typename DESCRIPTOR>
35void setInterpolatedConvectionBoundary(SuperLattice<T, DESCRIPTOR>& sLattice,T omega, SuperGeometry<T,2>& superGeometry, int material, T* uAv)
36{
37 setInterpolatedConvectionBoundary<T,DESCRIPTOR>(sLattice, omega, superGeometry.getMaterialIndicator(material), uAv);
38}
40template<typename T, typename DESCRIPTOR>
42{
43 OstreamManager clout(std::cout, "setInterpolatedConvectionBoundary");
44 int _overlap = 1;
45 bool includeOuterCells = false;
46 if (indicator->getSuperGeometry().getOverlap() == 1) {
47 includeOuterCells = true;
48 clout << "WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
49 }
50 for (int iCloc = 0; iCloc < sLattice.getLoadBalancer().size(); ++iCloc) {
51 setInterpolatedConvectionBoundary<T,DESCRIPTOR>(sLattice.getBlock(iCloc), omega,
52 indicator->getBlockIndicatorF(iCloc), uAv, includeOuterCells);
53 }
55 addPoints2CommBC<T, DESCRIPTOR>(sLattice, std::forward<decltype(indicator)>(indicator), _overlap);
56}
57
59template<typename T, typename DESCRIPTOR>
61 T* uAv, bool includeOuterCells)
62{
63 auto& blockGeometryStructure = indicator.getBlockGeometry();
64 const int margin = includeOuterCells ? 0 : 1;
65 OstreamManager clout(std::cout, "setInterpolatedConvectionBoundary");
66 bool _output = false;
67 std::vector<int> discreteNormal(3, 0);
68 blockGeometryStructure.forSpatialLocations([&](auto iX, auto iY) {
69 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY}) >= margin
70 && indicator(iX, iY)) {
71 PostProcessorGenerator2D<T, DESCRIPTOR>* postProcessor = nullptr;
72 discreteNormal = indicator.getBlockGeometry().getStatistics().getType(iX, iY);
73 if (discreteNormal[0] == 0) {
74 if (discreteNormal[1] == -1) {
75 if (_output) {
76 clout << "setInterpolatedConvectionBoundary<" << 0 << ","<< -1 << ">(" << iX << ", "<< iY << ", " << omega << " )" << std::endl;
77 }
79 < T,DESCRIPTOR,0,-1 > (iX, iX, iY, iY, uAv);
80 }
81 else if (discreteNormal[1] == 1) {
82 if (_output) {
83 clout << "setInterpolatedConvectionBoundary<" << 0 << "," << 1 << ">(" << iX << ", "<< iY << ", " << omega << " )" << std::endl;
84 }
86 < T,DESCRIPTOR,0,1 > (iX, iX, iY, iY, uAv);
87 }
88 else if (discreteNormal[2] == -1) {
89 if (_output) {
90 clout << "setInterpolatedConvectionBoundary<" << 1 << "," << -1 << ">(" << iX << ", "<< iY << ", " << omega << " )" << std::endl;
91 }
93 < T,DESCRIPTOR,1,-1 > (iX, iX, iY, iY, uAv);
94 }
95 else if (discreteNormal[2] == 1) {
96 if (_output) {
97 clout << "setInterpolatedConvectionBoundary<" << 1 << "," << 1 << ">(" << iX << ", "<< iY << ", " << omega << " )" << std::endl;
98 }
100 < T,DESCRIPTOR,1,1 > (iX, iX, iY, iY, uAv);
101 }
102 if (postProcessor) {
103 block.addPostProcessor(*postProcessor);
104 }
105 }
106 }
107 });
108}
109
110
111
112
113}//namespace olb
114
115
116#endif
Base block indicator functor (discrete)
BlockGeometry< T, 2 > & getBlockGeometry()
Get underlying block geometry structure.
Platform-abstracted block lattice for external access and inter-block interaction.
virtual void addPostProcessor(std::type_index stage, LatticeR< DESCRIPTOR::d > latticeR, PostProcessorPromise< T, DESCRIPTOR > &&promise)=0
Schedule post processor for application to latticeR in stage.
Smart pointer for managing the various ways of passing functors around.
Definition functorPtr.h:60
class for marking output with some text
Representation of a statistic for a parallel 2D geometry.
std::unique_ptr< SuperIndicatorF< T, D > > getMaterialIndicator(std::vector< int > &&materials)
Returns a material indicator using the given vector of materials.
Super class maintaining block lattices for a cuboid decomposition.
BlockLattice< T, DESCRIPTOR > & getBlock(int locC)
Return BlockLattice with local index locC.
LoadBalancer< T > & getLoadBalancer()
Read and write access to the load balancer.
Top level namespace for all of OpenLB.
void setInterpolatedConvectionBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material, T *uAv)
Initialising the InterpolatedConvectionBoundary on the superLattice domain This is an interpolated bo...