OpenLB 1.7
Loading...
Searching...
No Matches
latticePhysViscosity3D.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_VISCOSITY_3D_HH
26#define LATTICE_PHYS_VISCOSITY_3D_HH
27
28#include<vector> // for generic i/o
29#include<cmath> // for lpnorm
30#include<math.h>
31
33#include "superBaseF3D.h"
36#include "dynamics/lbm.h" // for computation of lattice rho and velocity
38#include "blockBaseF3D.h"
41
42namespace olb {
43
44template<typename T, typename DESCRIPTOR>
46 SuperLattice<T, DESCRIPTOR>& sLattice, const UnitConverter<T,DESCRIPTOR>& converter, bool logscale)
47 : SuperLatticePhysF3D<T, DESCRIPTOR>(sLattice, converter, 1)
48{
49 this->getName() = "physViscosity";
50 const int maxC = this->_sLattice.getLoadBalancer().size();
51 this->_blockF.reserve(maxC);
52 for (int iC = 0; iC < maxC; iC++) {
53 this->_blockF.emplace_back(
55 this->_sLattice.getBlock(iC),
56 this->_converter, logscale)
57 );
58 }
59}
60
61template<typename T, typename DESCRIPTOR>
63 BlockLattice<T, DESCRIPTOR>& blockLattice,
64 const UnitConverter<T,DESCRIPTOR>& converter, bool logscale)
65 : BlockLatticePhysF3D<T, DESCRIPTOR>(blockLattice, converter, 1),
66 _logscale(logscale)
67{
68 this->getName() = "physViscosity";
69}
70
71template<typename T, typename DESCRIPTOR>
73{
74 Cell<T,DESCRIPTOR> cell = this->_blockLattice.get(input[0], input[1], input[2]);
75 const auto omegaPL = cell.template getField<descriptors::OMEGA>();
76 if (omegaPL == T()) {
77 output[0] = std::numeric_limits<double>::quiet_NaN();
78 return true;
79 }
80 output[0] = (1./omegaPL - 0.5) / descriptors::invCs2<T,DESCRIPTOR>()
81 * this->_converter.getPhysDeltaX() * this->_converter.getPhysDeltaX() / this->_converter.getPhysDeltaT();
82 if (_logscale) {
83 output[0] = util::log10(output[0]);
84 }
85
86 return true;
87}
88
89}
90#endif
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
functor returns pointwise phys viscosity on local lattices
BlockLatticePhysViscosity3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter, bool logscale)
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
Platform-abstracted block lattice for external access and inter-block interaction.
Highest-level interface to Cell data.
Definition cell.h:148
std::string & getName()
read and write access to name
Definition genericF.hh:51
std::vector< std::unique_ptr< BlockF3D< T > > > _blockF
Super functors may consist of several BlockF3D<W> derived functors.
SuperLattice< T, DESCRIPTOR > & _sLattice
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
SuperLatticePhysViscosity3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter, bool logscale=false)
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
Wrapper functions that simplify the use of MPI.
ADf< T, DIM > log10(const ADf< T, DIM > &a)
Definition aDiff.h:496
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.