Loading...
Searching...
No Matches
ParticleForce.C
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011-2017 OpenFOAM Foundation
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26\*---------------------------------------------------------------------------*/
27
28#include "ParticleForce.H"
29
30// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
31
32template<class CloudType>
34(
36 const fvMesh& mesh,
37 const dictionary& dict,
38 const word& forceType,
39 const bool readCoeffs
40)
41:
42 owner_(owner),
43 mesh_(mesh),
44 coeffs_(readCoeffs ? dict : dictionary::null)
45{
46 if (readCoeffs && (coeffs_.dictName() != forceType))
47 {
48 FatalIOErrorInFunction(dict)
49 << "Force " << forceType << " must be specified as a dictionary"
50 << exit(FatalIOError);
51 }
52}
53
54
55template<class CloudType>
57:
58 owner_(pf.owner_),
59 mesh_(pf.mesh_),
60 coeffs_(pf.coeffs_)
61{}
62
63
64// * * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * //
65
66template<class CloudType>
69
70
71// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
73template<class CloudType>
75{}
76
77
78template<class CloudType>
80(
81 const typename CloudType::parcelType&,
82 const typename CloudType::parcelType::trackingData& td,
83 const scalar dt,
84 const scalar mass,
85 const scalar Re,
86 const scalar muc
87) const
88{
89 forceSuSp value;
90 value.Su() = Zero;
91 value.Sp() = 0.0;
92
93 return value;
94}
95
96
97template<class CloudType>
99(
100 const typename CloudType::parcelType&,
101 const typename CloudType::parcelType::trackingData& td,
102 const scalar dt,
103 const scalar mass,
104 const scalar Re,
105 const scalar muc
106) const
107{
108 forceSuSp value;
109 value.Su() = Zero;
110 value.Sp() = 0.0;
111
112 return value;
113}
114
115
116template<class CloudType>
118(
119 const typename CloudType::parcelType& p,
120 const typename CloudType::parcelType::trackingData& td,
121 const scalar mass
122) const
123{
124 return 0.0;
125}
126
127
128// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
129
130#include "ParticleForceNew.C"
131
132// ************************************************************************* //
Abstract base class for particle forces.
const fvMesh & mesh() const noexcept
Return the mesh database.
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.
virtual void cacheFields(const bool store)
Cache fields.
virtual ~ParticleForce()
Destructor.
const CloudType & owner() const noexcept
Return const access to the cloud owner.
virtual scalar massAdd(const typename CloudType::parcelType &p, const typename CloudType::parcelType::trackingData &td, const scalar mass) const
Return the added mass.
ParticleForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict, const word &forceType, const bool readCoeffs)
Construct from mesh.
virtual forceSuSp calcCoupled(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 coupled force.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition dictionary.H:133
Helper container for force Su and Sp terms.
Definition forceSuSp.H:63
const vector & Su() const
Return const access to the explicit contribution [kg.m/s2].
Definition forceSuSpI.H:54
scalar Sp() const
Return const access to the implicit coefficient [kg/s].
Definition forceSuSpI.H:60
Mesh data needed to do the Finite Volume discretisation.
Definition fvMesh.H:85
A class for handling words, derived from Foam::string.
Definition word.H:66
volScalarField & p
wallPoints::trackData td(isBlockedFace, regionToBlockSize)
DSMCCloud< dsmcParcel > CloudType
static constexpr const zero Zero
Global zero (0).
Definition zero.H:127
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dictionary dict