/usr/include/palabos/latticeBoltzmann/entropicLbTemplates2D.h is in libplb-dev 1.5~r1+repack1-2build2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 | /* This file is part of the Palabos library.
*
* Copyright (C) 2011-2015 FlowKit Sarl
* Route d'Oron 2
* 1010 Lausanne, Switzerland
* E-mail contact: contact@flowkit.com
*
* The most recent release of Palabos can be downloaded at
* <http://www.palabos.org/>
*
* The library Palabos is free software: you can redistribute it and/or
* modify it under the terms of the GNU Affero General Public License as
* published by the Free Software Foundation, either version 3 of the
* License, or (at your option) any later version.
*
* The library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*/
/* Orestis Malaspinas contributed this code.
*/
/** \file
* A collection of dynamics classes (e.g. BGK) with which a Cell object
* can be instantiated -- header file.
*/
#ifndef ENTROPIC_LB_HELPERS_2D_H
#define ENTROPIC_LB_HELPERS_2D_H
#include "core/globalDefs.h"
#include "latticeBoltzmann/nearestNeighborLattices2D.h"
#include <cmath>
namespace plb {
template<typename T>
struct entropicLbTemplates<T, descriptors::D2Q9Descriptor>
{
typedef descriptors::D2Q9Descriptor<T> Descriptor;
/// Computation of equilibrium distribution with an expansion
/// with respect to a small velocity u
static T equilibrium( plint iPop, T rho, Array<T,2> const& u)
{
T c_u = Descriptor::c[iPop][0]*u[0] + Descriptor::c[iPop][1]*u[1];
T c_u2 = c_u*c_u;
T c_u3 = c_u2*c_u;
T c_u4 = c_u3*c_u;
T c_u5 = c_u4*c_u;
T c_u6 = c_u5*c_u;
T c_u7 = c_u6*c_u;
T uSqr = u[0]*u[0] + u[1]*u[1];
T uSqr2 = uSqr*uSqr;
T uSqr3 = uSqr2*uSqr;
T powUx = u[0]*u[0]*u[0]*u[0]*u[0]; // u_x^5
T powUy = u[1]*u[1]*u[1]*u[1]*u[1]; // u_y^5
T C = Descriptor::c[iPop][0] * powUx + Descriptor::c[iPop][1] * powUy;
powUx *= u[0]; // u_x^6
powUy *= u[1]; // u_y^6
T E = powUx + powUy;
powUx *= u[0]; // u_x^7
powUy *= u[1]; // u_y^7
T F = Descriptor::c[iPop][0] * powUx + Descriptor::c[iPop][1] * powUy;
return Descriptor::t[iPop] * rho *
((T)1
+ c_u*(C*(T)81/(T)20 + uSqr2*(T)27/(T)8 - uSqr*(T)9/(T)2
- E*(T)81/(T)24 - uSqr3*(T)81/(T)48 + (T)3)
+ c_u2*(uSqr2*(T)81/(T)16 - uSqr*(T)27/(T)4
+ C*(T)243/(T)40 + (T)9/(T)2)
+ c_u3*(uSqr2*(T)243/(T)48 - uSqr*(T)81/(T)12 + (T)27/(T)6)
- c_u4*uSqr*(T)243/(T)48
+ c_u4*(T)81/(T)24
- c_u5*uSqr*(T)729/(T)240
+ c_u5*(T)243/(T)120
+ c_u6*(T)729/(T)720
+ c_u7*(T)2187/(T)5040
- C*uSqr*(T)81/(T)40
+ C*(T)27/(T)20 - uSqr3*(T)27/(T)48 - E*(T)27/(T)24
- F*(T)81/(T)56 - uSqr*(T)3/(T)2 + uSqr2*(T)9/(T)8
)
- Descriptor::SkordosFactor() * Descriptor::t[iPop];
}
};
}
#endif
|