Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #include "HeatTransferModel.H"
00027
00028
00029
00030 template<class CloudType>
00031 Foam::HeatTransferModel<CloudType>::HeatTransferModel(CloudType& owner)
00032 :
00033 dict_(dictionary::null),
00034 owner_(owner),
00035 coeffDict_(dictionary::null),
00036 BirdCorrection_(false)
00037 {}
00038
00039
00040 template<class CloudType>
00041 Foam::HeatTransferModel<CloudType>::HeatTransferModel
00042 (
00043 const dictionary& dict,
00044 CloudType& owner,
00045 const word& type
00046 )
00047 :
00048 dict_(dict),
00049 owner_(owner),
00050 coeffDict_(dict.subDict(type + "Coeffs")),
00051 BirdCorrection_(coeffDict_.lookup("BirdCorrection"))
00052 {}
00053
00054
00055
00056
00057 template<class CloudType>
00058 Foam::HeatTransferModel<CloudType>::~HeatTransferModel()
00059 {}
00060
00061
00062
00063
00064 template<class CloudType>
00065 const CloudType& Foam::HeatTransferModel<CloudType>::owner() const
00066 {
00067 return owner_;
00068 }
00069
00070
00071 template<class CloudType>
00072 const Foam::dictionary& Foam::HeatTransferModel<CloudType>::dict() const
00073 {
00074 return dict_;
00075 }
00076
00077
00078 template<class CloudType>
00079 const Foam::dictionary& Foam::HeatTransferModel<CloudType>::coeffDict() const
00080 {
00081 return coeffDict_;
00082 }
00083
00084
00085 template<class CloudType>
00086 const Foam::Switch& Foam::HeatTransferModel<CloudType>::BirdCorrection() const
00087 {
00088 return BirdCorrection_;
00089 }
00090
00091
00092 template<class CloudType>
00093 Foam::scalar Foam::HeatTransferModel<CloudType>::htc
00094 (
00095 const scalar dp,
00096 const scalar Re,
00097 const scalar Pr,
00098 const scalar kappa,
00099 const scalar NCpW
00100 ) const
00101 {
00102 const scalar Nu = this->Nu(Re, Pr);
00103
00104 scalar htc = Nu*kappa/dp;
00105
00106 if (BirdCorrection_ && (mag(htc) > ROOTVSMALL) && (mag(NCpW) > ROOTVSMALL))
00107 {
00108 const scalar phit = min(NCpW/htc, 50);
00109 if (phit > 0.001)
00110 {
00111 htc *= phit/(exp(phit) - 1.0);
00112 }
00113 }
00114
00115 return htc;
00116 }
00117
00118
00119
00120
00121 #include "NewHeatTransferModel.C"
00122
00123
00124
00125