src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/HeatTransferModel/HeatTransferModel.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "HeatTransferModel.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 template<class CloudType>
  32 Foam::HeatTransferModel<CloudType>::HeatTransferModel(CloudType& owner)
  33 :
  34     dict_(dictionary::null),
  35     owner_(owner),
  36     coeffDict_(dictionary::null)
  37 {}
  38
  39
  40 template<class CloudType>
  41 Foam::HeatTransferModel<CloudType>::HeatTransferModel
  42 (
  43     const dictionary& dict,
  44     CloudType& owner,
  45     const word& type
  46 )
  47 :
  48     dict_(dict),
  49     owner_(owner),
  50     coeffDict_(dict.subDict(type + "Coeffs"))
  51 {}
  52
  53
  54 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  55
  56 template<class CloudType>
  57 Foam::HeatTransferModel<CloudType>::~HeatTransferModel()
  58 {}
  59
  60
  61 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  62
  63 template<class CloudType>
  64 const CloudType& Foam::HeatTransferModel<CloudType>::owner() const
  65 {
  66     return owner_;
  67 }
  68
  69
  70 template<class CloudType>
  71 const Foam::dictionary& Foam::HeatTransferModel<CloudType>::dict() const
  72 {
  73     return dict_;
  74 }
  75
  76
  77 template<class CloudType>
  78 const Foam::dictionary& Foam::HeatTransferModel<CloudType>::coeffDict() const
  79 {
  80     return coeffDict_;
  81 }
  82
  83
  84 template<class CloudType>
  85 Foam::scalar Foam::HeatTransferModel<CloudType>::h
  86 (
  87     const scalar dp,
  88     const vector& Ur,
  89     const scalar rhoc,
  90     const scalar rhop,
  91     const scalar cpc,
  92     const scalar cpp,
  93     const scalar muc
  94 ) const
  95 {
  96     const scalar Re = rhoc*mag(Ur)*dp/(muc + ROOTVSMALL);
  97
  98 //    const scalar Pr = muc/alphac;
  99     const scalar Pr = this->Pr();
 100
 101     const scalar Nu = this->Nu(Re, Pr);
 102
 103     const scalar kappa = cpc*muc/Pr;
 104
 105     return Nu*kappa/dp;
 106 }
 107
 108
 109 template<class CloudType>
 110 Foam::scalar Foam::HeatTransferModel<CloudType>::Cu
 111 (
 112     const scalar dp,
 113     const vector& Ur,
 114     const scalar rhoc,
 115     const scalar rhop,
 116     const scalar cpc,
 117     const scalar cpp,
 118     const scalar muc
 119 ) const
 120 {
 121     const scalar Re = rhoc*mag(Ur)*dp/(muc + ROOTVSMALL);
 122
 123 //    const scalar Pr = muc/alphac;
 124     const scalar Pr = this->Pr();
 125
 126     const scalar Nu = this->Nu(Re, Pr);
 127
 128     const scalar kappa = cpc*muc/Pr;
 129
 130     const scalar htc = Nu*kappa/dp;
 131
 132     return 6.0*htc/(dp*rhop*cpp);
 133 }
 134
 135
 136 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 137
 138 #include "NewHeatTransferModel.C"
 139
 140
 141 // ************************************************************************* //
 142