OpenLB 1.7
Loading...
Searching...
No Matches
latticePhysCorrBoundaryForce2D.hh
Go to the documentation of this file.
1/* This file is part of the OpenLB library
2 *
3 * Copyright (C) 2012 Lukas Baron, Tim Dornieden, Mathias J. Krause,
4 * Albert Mink
5 * E-mail contact: info@openlb.net
6 * The most recent release of OpenLB can be downloaded at
7 * <http://www.openlb.net/>
8 *
9 * This program is free software; you can redistribute it and/or
10 * modify it under the terms of the GNU General Public License
11 * as published by the Free Software Foundation; either version 2
12 * of the License, or (at your option) any later version.
13 *
14 * This program is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 * GNU General Public License for more details.
18 *
19 * You should have received a copy of the GNU General Public
20 * License along with this program; if not, write to the Free
21 * Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
22 * Boston, MA 02110-1301, USA.
23*/
24
25#ifndef LATTICE_PHYS_CORR_BOUNDARY_FORCE_2D_HH
26#define LATTICE_PHYS_CORR_BOUNDARY_FORCE_2D_HH
27
28#include <vector>
29#include "utilities/omath.h"
30#include <limits>
31
33#include "superBaseF2D.h"
36#include "dynamics/lbm.h" // for computation of lattice rho and velocity
38#include "blockBaseF2D.h"
41
42namespace olb {
43
44template<typename T, typename DESCRIPTOR>
47 FunctorPtr<SuperIndicatorF2D<T>>&& indicatorF,
48 const UnitConverter<T,DESCRIPTOR>& converter)
49 : SuperLatticePhysF2D<T, DESCRIPTOR>(sLattice, converter, 2),
50 _indicatorF(std::move(indicatorF))
51{
52 this->getName() = "physCorrBoundaryForce";
53 for (int iC = 0; iC < this->_sLattice.getLoadBalancer().size(); ++iC) {
54 this->_blockF.emplace_back(
56 this->_sLattice.getBlock(iC),
57 _indicatorF->getBlockIndicatorF(iC),
58 this->_converter));
59 }
60}
61
62template<typename T, typename DESCRIPTOR>
65 SuperGeometry<T,2>& superGeometry, const int material,
66 const UnitConverter<T,DESCRIPTOR>& converter)
68 superGeometry.getMaterialIndicator(material),
69 converter)
70{ }
71
72template <typename T, typename DESCRIPTOR>
74 BlockLattice<T,DESCRIPTOR>& blockLattice,
75 BlockIndicatorF2D<T>& indicatorF,
76 const UnitConverter<T,DESCRIPTOR>& converter)
77 : BlockLatticePhysF2D<T,DESCRIPTOR>(blockLattice, converter, 2),
78 _indicatorF(indicatorF),
79 _blockGeometry(indicatorF.getBlockGeometry())
80{
81 this->getName() = "physCorrBoundaryForce";
82}
83
84template<typename T, typename DESCRIPTOR>
86 T output[], const int input[])
87{
88 for (int i = 0; i < this->getTargetDim(); ++i) {
89 output[i] = T();
90 }
91
92 if (_indicatorF(input)) {
93 for (int iPop = 1; iPop < DESCRIPTOR::q; ++iPop) {
94 // Get direction
95 const Vector<int,2> c = descriptors::c<DESCRIPTOR>(iPop);
96 // Get next cell located in the current direction
97 // Check if the next cell is a fluid node
98 if (_blockGeometry.get({input[0] + c[0], input[1] + c[1]}) == 1) {
99 // Get f_q of next fluid cell where l = opposite(q)
100 T f = this->_blockLattice.get(input[0] + c[0], input[1] + c[1])[iPop];
101 // Get f_l of the boundary cell
102 // Add f_q and f_opp
103 f += this->_blockLattice.get(input[0], input[1])[descriptors::opposite<DESCRIPTOR>(iPop)];
104 // Update force
105 for (int i = 0; i < this->getTargetDim(); ++i) {
106 output[i] -= c[i] * (f - 2. * descriptors::t<T,DESCRIPTOR>(iPop));
107 }
108 }
109 }
110 for (int i = 0; i < this->getTargetDim(); ++i) {
111 output[i] = this->_converter.getPhysForce(output[i]);
112 }
113 }
114 return true;
115}
116
117}
118#endif
Base block indicator functor (discrete)
functor returns pointwise phys force acting on a indicated boundary on local lattice see: Caiazzo,...
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysCorrBoundaryForce2D(BlockLattice< T, DESCRIPTOR > &blockLattice, BlockIndicatorF2D< T > &indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Platform-abstracted block lattice for external access and inter-block interaction.
Smart pointer for managing the various ways of passing functors around.
Definition functorPtr.h:60
std::string & getName()
read and write access to name
Definition genericF.hh:51
std::vector< std::unique_ptr< BlockF2D< T > > > _blockF
Super functors may consist of several BlockF2D<W> derived functors.
Representation of a statistic for a parallel 2D geometry.
SuperLattice< T, DESCRIPTOR > & _sLattice
functor to get pointwise phys force acting on a boundary with a given material on local lattice see: ...
SuperLatticePhysCorrBoundaryForce2D(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
Plain old scalar vector.
Definition vector.h:47
Wrapper functions that simplify the use of MPI.
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.