SRFForce.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "SRFForce.H"
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class CloudType>
32 (
33  CloudType& owner,
34  const fvMesh& mesh,
35  const dictionary& dict
36 )
37 :
38  ParticleForce<CloudType>(owner, mesh, dict, typeName, false),
39  srfPtr_(nullptr)
40 {}
41 
42 
43 template<class CloudType>
45 (
46  const SRFForce& srff
47 )
48 :
50  srfPtr_(nullptr)
51 {}
52 
53 
54 // * * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * //
55 
56 template<class CloudType>
58 {}
59 
60 
61 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
62 
63 template<class CloudType>
65 {
66  if (store)
67  {
68  const typename SRF::SRFModel& model = this->mesh().template
69  lookupObject<SRF::SRFModel>("SRFProperties");
70  srfPtr_ = &model;
71  }
72  else
73  {
74  srfPtr_ = nullptr;
75  }
76 }
77 
78 
79 template<class CloudType>
81 (
82  const typename CloudType::parcelType& p,
83  const typename CloudType::parcelType::trackingData& td,
84  const scalar dt,
85  const scalar mass,
86  const scalar Re,
87  const scalar muc
88 ) const
89 {
90  forceSuSp value(Zero, 0.0);
91 
92  const typename SRF::SRFModel& srf = *srfPtr_;
93 
94  const vector& omega = srf.omega().value();
95 
96  const vector& r = p.position();
97 
98  // Coriolis and centrifugal acceleration terms
99  value.Su() =
100  mass*(1.0 - td.rhoc()/p.rho())
101  *(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
102 
103  return value;
104 }
105 
106 
107 // ************************************************************************* //
dictionary dict
Calculates particle SRF reference frame force.
Definition: SRFForce.H:52
const vector & Su() const
Return const access to the explicit contribution [kg m/s^2].
Definition: forceSuSpI.H:56
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
virtual ~SRFForce()
Destructor.
Definition: SRFForce.C:57
Abstract base class for particle forces.
Definition: ParticleForce.H:53
virtual void cacheFields(const bool store)
Cache fields.
Definition: SRFForce.C:64
Helper container for force Su and Sp terms.
Definition: forceSuSp.H:61
dynamicFvMesh & mesh
const Type & value() const
Return const reference to value.
const dimensionedVector & omega() const
Return the angular velocity field [rad/s].
Definition: SRFModel.C:115
static const zero Zero
Definition: zero.H:97
SRFForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
Definition: SRFForce.C:32
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Definition: DSMCCloud.H:215
virtual forceSuSp calcNonCoupled(const typename CloudType::parcelType &p, const typename CloudType::parcelType::trackingData &td, const scalar dt, const scalar mass, const scalar Re, const scalar muc) const
Calculate the non-coupled force.
Definition: SRFForce.C:81
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:69
Top level model for single rotating frame.
Definition: SRFModel.H:62