OpenLB 1.7
Loading...
Searching...
No Matches
latticePhysTemperature2D.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_PHYS_TEMPERATURE_2D_HH
25#define LATTICE_PHYS_TEMPERATURE_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, typename TDESCRIPTOR>
47 : SuperLatticeThermalPhysF2D<T,DESCRIPTOR,TDESCRIPTOR>(sLattice, converter, 1)
48{
49 this->getName() = "physTemperature";
50 int maxC = this->_sLattice.getLoadBalancer().size();
51 this->_blockF.reserve(maxC);
52 for (int iC = 0; iC < maxC; iC++) {
53 this->_blockF.emplace_back(new BlockLatticePhysTemperature2D<T,DESCRIPTOR,TDESCRIPTOR>(this->_sLattice.getBlock(iC), this->_converter));
54 }
55}
56
57template <typename T, typename DESCRIPTOR, typename TDESCRIPTOR>
60 : BlockLatticeThermalPhysF2D<T,DESCRIPTOR,TDESCRIPTOR>(blockLattice,converter,1)
61{
62 this->getName() = "physTemperature";
63}
64
65
66template <typename T, typename DESCRIPTOR, typename TDESCRIPTOR>
68{
69 T latticeTemperature = this->_blockLattice.get( input[0], input[1] ).computeRho();
70 output[0] = this->_converter.getPhysTemperature(latticeTemperature);
71
72 return true;
73}
74
75}
76#endif
BlockLatticePhysTemperature2D returns pointwise phys temperature from rho on local lattices.
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysTemperature2D(BlockLattice< T, TDESCRIPTOR > &blockLattice, ThermalUnitConverter< T, DESCRIPTOR, TDESCRIPTOR > const &converter)
represents all thermal functors that operate on a DESCRIPTOR with output in Phys, e....
Platform-abstracted block lattice for external access and inter-block interaction.
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.
SuperLattice< T, TDESCRIPTOR > & _sLattice
SuperLatticePhysTemperature2D(SuperLattice< T, TDESCRIPTOR > &sLattice, ThermalUnitConverter< T, DESCRIPTOR, TDESCRIPTOR > const &converter)
represents all thermal functors that operate on a DESCRIPTOR with output in Phys, e....
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.
Conversion between physical and lattice units, as well as discretization specialized for thermal appl...
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.