OpenLB 1.7
Loading...
Searching...
No Matches
setInterpolatedPressureBoundary2D.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 Pressure Boundary
25//This is a new version of the Boundary, which only contains free floating functions
26
27#ifndef SET_INTERPOLATED_PRESSURE_BOUNDARY_2D_HH
28#define SET_INTERPOLATED_PRESSURE_BOUNDARY_2D_HH
29
31
32namespace olb {
34template<typename T, typename DESCRIPTOR, class MixinDynamics>
35void setInterpolatedPressureBoundary(SuperLattice<T, DESCRIPTOR>& sLattice, T omega, SuperGeometry<T,2>& superGeometry, int material)
36{
37 setInterpolatedPressureBoundary<T,DESCRIPTOR, MixinDynamics>(sLattice,omega, superGeometry.getMaterialIndicator(material));
38}
39
41template<typename T, typename DESCRIPTOR, class MixinDynamics>
43{
44 OstreamManager clout(std::cout, "setLocalPressureBoundary2D");
45 bool includeOuterCells = false;
46 int _overlap = 1;
47 if (indicator->getSuperGeometry().getOverlap() == 1) {
48 includeOuterCells = true;
49 clout << "WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
50 }
51 for (int iCloc = 0; iCloc < sLattice.getLoadBalancer().size(); ++iCloc) {
52 setInterpolatedPressureBoundary<T,DESCRIPTOR, MixinDynamics>(sLattice.getBlock(iCloc),omega,
53 indicator->getBlockIndicatorF(iCloc), includeOuterCells);
54 }
56 addPoints2CommBC(sLattice, std::forward<decltype(indicator)>(indicator), _overlap);
57}
58
59
61template<typename T, typename DESCRIPTOR, class MixinDynamics>
62void setInterpolatedPressureBoundary(BlockLattice<T,DESCRIPTOR>& block,T omega, BlockIndicatorF2D<T>& indicator, bool includeOuterCells)
63{
64 using namespace boundaryhelper;
65 auto& blockGeometryStructure = indicator.getBlockGeometry();
66 const int margin = includeOuterCells ? 0 : 1;
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 Dynamics<T, DESCRIPTOR>* dynamics = nullptr;
72 discreteNormal = indicator.getBlockGeometry().getStatistics().getType(iX, iY);
73 if (discreteNormal[0] == 0) {
74 dynamics = block.getDynamics(MixinDynamicsExchangeDirectionOrientationMomenta<T,DESCRIPTOR,
76 >::construct(Vector<int,2>(discreteNormal.data() + 1)));
77 block.addPostProcessor(
78 typeid(stage::PostStream), {iX, iY},
79 promisePostProcessorForDirectionOrientation<T,DESCRIPTOR,StraightFdBoundaryProcessor2D>(
80 Vector<int,2>(discreteNormal.data() + 1)));
81 setBoundary(block,iX,iY, dynamics);
82 }
83 }
84 });
85}
86
87
88}//namespace olb
89
90#endif
91
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.
virtual Dynamics< T, DESCRIPTOR > * getDynamics(DynamicsPromise< T, DESCRIPTOR > &&)=0
Return pointer to dynamics yielded by promise.
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.
Plain old scalar vector.
Definition vector.h:47
Top level namespace for all of OpenLB.
void setInterpolatedPressureBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material)
Initialising the setInterpolatedPressureBoundary function on the superLattice domain Interpolated Bou...
void addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.
void setBoundary(BlockLattice< T, DESCRIPTOR > &block, int iX, int iY, Dynamics< T, DESCRIPTOR > *dynamics, PostProcessorGenerator2D< T, DESCRIPTOR > *postProcessor)
Interface for per-cell dynamics.
Definition interface.h:54
Communication after propagation.
Definition stages.h:36