vhlle/src/trancoeff.cpp

0001 #include "trancoeff.h"
0002 #include "eos.h"
0003 #include "inc.h"
0004
0005 TransportCoeff::TransportCoeff(double _etaS, double _zetaS, EoS *_eos)
0006 {
0007  etaS = _etaS ;
0008  zetaS = _zetaS ;
0009  eos = _eos ;
0010 }
0011
0012 void TransportCoeff::getEta(double e, double T, double &_etaS, double &_zetaS)
0013 {
0014  _etaS = etaS ;
0015  _zetaS=zetaS*(1./3.-eos->cs2(e))/(exp((0.16-T)/0.001)+1.) ;
0016 }
0017
0018
0019 void TransportCoeff::getTau(double T, double &_taupi, double &_tauPi)
0020 {
0021     if(T>0.) _taupi=std::max(3.0/5.068*etaS/T,0.003) ; else _taupi=0.1 ;
0022     //if(T>0.) _tauPi=std::max(3.0/5.068*(1./4./C_PI)/T,0.005) ; else _tauPi=0.1 ;
0023   _tauPi = 0.1 ; // for I-S analytical solution
0024 }