OpenLB
1.7
Loading...
Searching...
No Matches
src
particles
subgrid3DLegacyFramework
forces
linearDampingForceForMagDipoleMoment3D.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
26
#ifndef LinearDampingForceForMagDipoleMoment3D_HH
27
#define LinearDampingForceForMagDipoleMoment3D_HH
28
29
#include "
utilities/omath.h
"
30
#include <vector>
31
#include "
linearDampingForceForMagDipoleMoment3D.h
"
32
33
#ifndef M_PI
34
#define M_PI 3.14159265358979323846
35
#endif
36
37
namespace
olb
{
38
39
template
<
typename
T,
template
<
typename
U>
class
PARTICLETYPE,
typename
DESCRIPTOR>
40
LinearDampingForceForMagDipoleMoment3D<T, PARTICLETYPE, DESCRIPTOR>::LinearDampingForceForMagDipoleMoment3D
(T dynVisc, T frictionFac):
41
Force3D
<T, PARTICLETYPE>(), _dynVisc(dynVisc), _frictionFac(frictionFac)
42
{
43
// this->_name = "LinearDampingForceForMagDipoleMoment3D";
44
}
45
46
template
<
typename
T,
template
<
typename
U>
class
PARTICLETYPE,
typename
DESCRIPTOR>
47
void
LinearDampingForceForMagDipoleMoment3D<T, PARTICLETYPE, DESCRIPTOR>::applyForce
(
48
typename
std::deque<PARTICLETYPE<T> >::iterator p,
int
pInt,
49
ParticleSystem3D<T, PARTICLETYPE>
& pSys)
50
{
51
52
Vector<T, 3>
dampingForce = {T(0), T(0), T(0)} ;
53
54
for
(
int
i = 0; i < 3; i++) {
55
56
dampingForce[i] = - 8 *
M_PI
*
util::pow
(p->getRad(), 3.) * _dynVisc * p->getAVel()[i] * _frictionFac ;
57
p->getTorque()[i] += dampingForce[i] ;
58
}
59
}
60
61
}
62
63
#endif
// LinearDampingForceForMagDipoleMoment3D_HH
M_PI
#define M_PI
Definition
analyticalF.hh:43
olb::Force3D
Prototype for all particle forces.
Definition
force3D.h:43
olb::LinearDampingForceForMagDipoleMoment3D::LinearDampingForceForMagDipoleMoment3D
LinearDampingForceForMagDipoleMoment3D(T dynVisc, T frictionFac)
Definition
linearDampingForceForMagDipoleMoment3D.hh:40
olb::LinearDampingForceForMagDipoleMoment3D::applyForce
void applyForce(typename std::deque< PARTICLETYPE< T > >::iterator p, int pInt, ParticleSystem3D< T, PARTICLETYPE > &psSys) override
Definition
linearDampingForceForMagDipoleMoment3D.hh:47
olb::ParticleSystem3D
Definition
particleSystem3D.h:73
olb::Vector
Plain old scalar vector.
Definition
vector.h:47
linearDampingForceForMagDipoleMoment3D.h
olb::util::pow
cpu::simd::Pack< T > pow(cpu::simd::Pack< T > base, cpu::simd::Pack< T > exp)
Definition
pack.h:112
olb
Top level namespace for all of OpenLB.
Definition
boundaryPostProcessors2D.h:34
omath.h
Generated on Tue Mar 5 2024 09:57:52 for OpenLB by
1.10.0