src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 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>
  31 Foam::SRFForce<CloudType>::SRFForce
  32 (
  33     CloudType& owner,
  34     const fvMesh& mesh,
  35     const dictionary& dict
  36 )
  37 :
  38     ParticleForce<CloudType>(owner, mesh, dict, typeName, false),
  39     srfPtr_(NULL)
  40 {}
  41
  42
  43 template<class CloudType>
  44 Foam::SRFForce<CloudType>::SRFForce
  45 (
  46     const SRFForce& srff
  47 )
  48 :
  49     ParticleForce<CloudType>(srff),
  50     srfPtr_(NULL)
  51 {}
  52
  53
  54 // * * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * //
  55
  56 template<class CloudType>
  57 Foam::SRFForce<CloudType>::~SRFForce()
  58 {}
  59
  60
  61 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  62
  63 template<class CloudType>
  64 void Foam::SRFForce<CloudType>::cacheFields(const bool store)
  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_ = NULL;
  75     }
  76 }
  77
  78
  79 template<class CloudType>
  80 Foam::forceSuSp Foam::SRFForce<CloudType>::calcNonCoupled
  81 (
  82     const typename CloudType::parcelType& p,
  83     const scalar dt,
  84     const scalar mass,
  85     const scalar Re,
  86     const scalar muc
  87 ) const
  88 {
  89     forceSuSp value(vector::zero, 0.0);
  90
  91     const typename SRF::SRFModel& srf = *srfPtr_;
  92
  93     const vector& omega = srf.omega().value();
  94
  95     const vector& r = p.position();
  96
  97     // Coriolis and centrifugal acceleration terms
  98     value.Su() = mass*(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
  99
 100     return value;
 101 }
 102
 103
 104 // ************************************************************************* //