OpenLB
1.7
Loading...
Searching...
No Matches
src
particles
subgrid3DLegacyFramework
forces
stokesDragForceForHomVelField3D.hh
Go to the documentation of this file.
1
/* This file is part of the OpenLB library
2
*
3
* Copyright (C) 2017 Marie-Luise Maier, Mathias J. Krause, Sascha Janz
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
25
#ifndef StokesDragForceForHomVelField3D_HH
26
#define StokesDragForceForHomVelField3D_HH
27
28
#include "
utilities/omath.h
"
29
#include "
stokesDragForceForHomVelField3D.h
"
30
31
namespace
olb
{
32
33
template
<
typename
T,
template
<
typename
U>
class
PARTICLETYPE,
typename
DESCRIPTOR>
34
StokesDragForceForHomVelField3D<T, PARTICLETYPE, DESCRIPTOR>::StokesDragForceForHomVelField3D
(T dynVisc,
Vector<T, 3>
fluidVel):
35
Force3D
<T, PARTICLETYPE>(), _dynVisc(dynVisc), _fluidVel(fluidVel)
36
{
37
// this->_name = "StokesDragForceForHomVelField3D";
38
}
39
40
template
<
typename
T,
template
<
typename
U>
class
PARTICLETYPE,
typename
DESCRIPTOR>
41
void
StokesDragForceForHomVelField3D<T, PARTICLETYPE, DESCRIPTOR>::applyForce
(
42
typename
std::deque<PARTICLETYPE<T> >::iterator p,
int
pInt,
43
ParticleSystem3D<T, PARTICLETYPE>
& pSys)
44
{
45
46
Vector<T, 3>
force = {T(0), T(0), T(0)} ;
47
48
for
(
int
i = 0; i < 3; i++) {
49
50
force[i] = -1 * 6 *
M_PI
* p->getRad() * (p->getVel()[i] - _fluidVel[i]) * _dynVisc ;
51
p->getForce()[i] += force[i] ;
52
}
53
54
}
55
}
56
57
#endif
// StokesDragForceForHomVelField3D
M_PI
#define M_PI
Definition
analyticalF.hh:43
olb::Force3D
Prototype for all particle forces.
Definition
force3D.h:43
olb::ParticleSystem3D
Definition
particleSystem3D.h:73
olb::StokesDragForceForHomVelField3D::applyForce
void applyForce(typename std::deque< PARTICLETYPE< T > >::iterator p, int pInt, ParticleSystem3D< T, PARTICLETYPE > &psSys) override
Definition
stokesDragForceForHomVelField3D.hh:41
olb::StokesDragForceForHomVelField3D::StokesDragForceForHomVelField3D
StokesDragForceForHomVelField3D(T dynVisc, Vector< T, 3 > fluidVel)
Definition
stokesDragForceForHomVelField3D.hh:34
olb::Vector
Plain old scalar vector.
Definition
vector.h:47
olb
Top level namespace for all of OpenLB.
Definition
boundaryPostProcessors2D.h:34
omath.h
stokesDragForceForHomVelField3D.h
Generated on Tue Mar 5 2024 09:57:52 for OpenLB by
1.10.0