OpenLB 1.7
Loading...
Searching...
No Matches
latticeVolumeFractionApproximation2D.hh
Go to the documentation of this file.
1/* This file is part of the OpenLB library
2 *
3 * Copyright (C) 2019 Albert Mink, Mathias J. Krause, Lukas Baron
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#ifndef LATTICE_VOLUME_FRACTION_APPROXIMATION_2D_HH
25#define LATTICE_VOLUME_FRACTION_APPROXIMATION_2D_HH
26
27#include <vector>
28#include "utilities/omath.h"
29#include <limits>
30
32#include "dynamics/lbm.h" // for computation of lattice rho and velocity
35#include "blockBaseF2D.h"
36#include "functors/genericF.h"
40
41
42namespace olb {
43
44template<typename T, typename DESCRIPTOR>
46 SuperLattice<T,DESCRIPTOR>& sLattice, SuperGeometry<T,2>& superGeometry,
47 IndicatorF2D<T>& indicator, int refinementLevel,
48 const UnitConverter<T,DESCRIPTOR>& converter, bool insideOut)
49 : SuperLatticeF2D<T, DESCRIPTOR>(sLattice, 1)
50{
51 this->getName() = "volumeFractionApproximation";
52 int maxC = this->_sLattice.getLoadBalancer().size();
53 this->_blockF.reserve(maxC);
54 for (int iC = 0; iC < maxC; iC++) {
55 this->_blockF.emplace_back(new BlockLatticeVolumeFractionApproximation2D<T, DESCRIPTOR>(this->_sLattice.getBlock(iC),
56 superGeometry.getBlockGeometry(iC),
57 indicator, refinementLevel,
58 converter, insideOut));
59 }
60}
61
62template<typename T, typename DESCRIPTOR>
64 BlockGeometry<T,2>& blockGeometry,
65 IndicatorF2D<T>& indicator,
66 int refinementLevel,
67 const UnitConverter<T,DESCRIPTOR>& converter,
68 bool insideOut)
69 : BlockLatticeF2D<T, DESCRIPTOR>(blockLattice, 1),
70 _blockGeometry(blockGeometry), _indicator(indicator), _refinementLevel(refinementLevel), _converter(converter), _insideOut(insideOut),
71 _physSubGridMinPhysRshift((1./ T(_refinementLevel * 2.) - 0.5) * _converter.getPhysDeltaX()),
72 _physSubGridDeltaX(1. / T(_refinementLevel) * _converter.getPhysDeltaX()),
73 _latticeSubGridVolume(1. / (_refinementLevel * _refinementLevel))
74{
75 this->getName() = "volumeFractionApproximation";
76}
77
78template<typename T, typename DESCRIPTOR>
80{
81 output[0] = 0.;
82 T physR[2];
83 bool inside[1];
84 _blockGeometry.getPhysR(physR, {input[0], input[1]});
85
86 T subGridMinPhysR[2];
87 subGridMinPhysR[0] = physR[0] + _physSubGridMinPhysRshift;
88 subGridMinPhysR[1] = physR[1] + _physSubGridMinPhysRshift;
89
90 for (int x = 0; x < _refinementLevel; x++) {
91 for (int y = 0; y < _refinementLevel; y++) {
92 physR[0] = subGridMinPhysR[0] + x * _physSubGridDeltaX;
93 physR[1] = subGridMinPhysR[1] + y * _physSubGridDeltaX;
94 _indicator(inside, physR);
95 if (!_insideOut) {
96 if (inside[0]) {
97 output[0] += _latticeSubGridVolume;
98 }
99 }
100 else {
101 if (!inside[0]) {
102 output[0] += _latticeSubGridVolume;
103 }
104 }
105 }
106 }
107 return true;
108}
109
110}
111#endif
Representation of a block geometry.
represents all functors that operate on a DESCRIPTOR in general, e.g. getVelocity(),...
functor returns pointwise an approximation for the volume fraction
BlockLatticeVolumeFractionApproximation2D(BlockLattice< T, DESCRIPTOR > &blockLattice, BlockGeometry< T, 2 > &blockGeometry, IndicatorF2D< T > &indicator, int refinementLevel, const UnitConverter< T, DESCRIPTOR > &converter, bool insideOut)
bool operator()(T output[], const int input[])
has to be implemented for 'every' derived class
Platform-abstracted block lattice for external access and inter-block interaction.
std::string & getName()
read and write access to name
Definition genericF.hh:51
IndicatorF2D is an application from .
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.
BlockGeometry< T, D > & getBlockGeometry(int locIC)
Read and write access to a single block geometry.
represents all functors that operate on a SuperLattice in general, e.g. getVelocity(),...
SuperLattice< T, DESCRIPTOR > & _sLattice
SuperLatticeVolumeFractionApproximation2D(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 2 > &superGeometry, IndicatorF2D< T > &indicator, int refinementLevel, const UnitConverter< T, DESCRIPTOR > &converter, bool insideOut=false)
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
The description of a generic interface for all functor classes – header file.
This file contains indicator functions.
Wrapper functions that simplify the use of MPI.
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.