diff --git a/Cantera/src/transport/LiquidTransport.cpp b/Cantera/src/transport/LiquidTransport.cpp index d581b3215..ce9b75e9b 100644 --- a/Cantera/src/transport/LiquidTransport.cpp +++ b/Cantera/src/transport/LiquidTransport.cpp @@ -974,6 +974,14 @@ namespace Cantera { * Updates the array of pure species viscosities, and the * weighting functions in the viscosity mixture rule. * The flag m_visc_ok is set to true. + * + * Note that for viscosity, a positive activation energy + * corresponds to the typical case of a positive argument + * to the exponential so that the Arrhenius expression is + * + * \f[ + * \mu = A T^n \exp( + E / R T ) + * \f] */ void LiquidTransport::updateViscosity_T() { int k; @@ -991,7 +999,7 @@ namespace Cantera { //m_coeffVisc_Ns[k][2] holds Tact //m_coeffVisc_Ns[k][3] holds log(A) m_logViscSpecies[k] = coeffk[3] + coeffk[1] * m_logt - - coeffk[2] / m_temp ; + + coeffk[2] / m_temp ; m_viscSpecies[k] = exp( m_logViscSpecies[k] ); } else if ( m_viscTempDepType_Ns[k] == LTR_MODEL_POLY ) { diff --git a/Cantera/src/transport/TransportFactory.cpp b/Cantera/src/transport/TransportFactory.cpp index 1d1ec5ffb..d3d96f53c 100644 --- a/Cantera/src/transport/TransportFactory.cpp +++ b/Cantera/src/transport/TransportFactory.cpp @@ -1126,10 +1126,25 @@ namespace Cantera { if ( transportNode.hasChild("viscosity")) { XML_Node& viscosityNode = transportNode.child("viscosity"); - string viscosityModel = viscosityNode.attrib("model"); - if (viscosityModel == "") { - throw CanteraError("LiquidTransport::initLiquid", + if ( viscosityNode.hasChild("compositionDependence") ) { + string viscosityModel = viscosityNode.child("compositionDependence").attrib("model"); + if (viscosityModel == "") { + throw CanteraError("LiquidTransport::getLiquidInteractionsTransportData", "transport::visosity XML node doesn't have a model string"); + } else if ( viscosityModel == "none"){ + ; + } else if ( viscosityModel == "solvent"){ + throw CanteraError("LiquidTransport::getLiquidInteractionsTransportData", + "solvent interactions not implemented"); + } else if ( viscosityModel == "moleFractions"){ + ; + } else if ( viscosityModel == "massFractions"){ + ; + } else if ( viscosityModel == "logMoleFractions"){ + ; + } else if ( viscosityModel == "pairwiseInteractionEnergy"){ + ; + } } } }