diff --git a/Cantera/src/transport/DustyGasTransport.cpp b/Cantera/src/transport/DustyGasTransport.cpp index dc0cead8d..9f1a50c25 100644 --- a/Cantera/src/transport/DustyGasTransport.cpp +++ b/Cantera/src/transport/DustyGasTransport.cpp @@ -42,8 +42,6 @@ namespace Cantera { DustyGasTransport::DustyGasTransport(thermo_t* thermo) : Transport(thermo), m_nsp(0), - m_tmin(0.0), - m_tmax(1.0E300), m_mw(0), m_dk(0), m_temp(-1.0), @@ -67,8 +65,6 @@ namespace Cantera { DustyGasTransport::DustyGasTransport(const DustyGasTransport &right) : Transport(), m_nsp(0), - m_tmin(0.0), - m_tmax(1.0E300), m_mw(0), m_dk(0), m_temp(-1.0), @@ -105,8 +101,6 @@ namespace Cantera { Transport::operator=(right); m_nsp = right.m_nsp; - m_tmin = right.m_tmin; - m_tmax = right.m_tmax; m_mw = right.m_mw; m_d = right.m_d; m_x = right.m_x; @@ -183,8 +177,6 @@ namespace Cantera { // constant mixture attributes m_thermo = phase; m_nsp = m_thermo->nSpecies(); - m_tmin = m_thermo->minTemp(); - m_tmax = m_thermo->maxTemp(); m_gastran = gastr; // make a local copy of the molecular weights @@ -332,21 +324,31 @@ namespace Cantera { "invert returned ierr = "+int2str(ierr)); } } - //==================================================================================================================== + //==================================================================================================================== + // Return the Multicomponent diffusion coefficients. Units: [m^2/s]. + /* + * Returns the array of multicomponent diffusion coefficients. + * + * @param ld The dimension of the inner loop of d (usually equal to m_nsp) + * @param d flat vector of diffusion coefficients, fortran ordering. + * d[ld*j+i] is the D_ij diffusion coefficient (the diffusion + * coefficient for species i due to species j). + */ void DustyGasTransport::getMultiDiffCoeffs(const int ld, doublereal* const d) { int i,j; updateMultiDiffCoeffs(); for (i = 0; i < m_nsp; i++) { - for (j = 0; j < m_nsp; j++) { + for (j = 0; j < m_nsp; j++) { d[ld*j + i] = m_multidiff(i,j); } } } - //==================================================================================================================== - /** - * Update temperature-dependent quantities. - */ + // Update temperature-dependent quantities within the object + /* + * The object keeps a value m_temp, which is the temperature at which quantities were last evaluated + * at. If the temperature is changed, update Booleans are set false, triggering recomputation. + */ void DustyGasTransport::updateTransport_T() { if (m_temp == m_thermo->temperature()) return; diff --git a/Cantera/src/transport/DustyGasTransport.h b/Cantera/src/transport/DustyGasTransport.h index dee73bdd3..36129948f 100644 --- a/Cantera/src/transport/DustyGasTransport.h +++ b/Cantera/src/transport/DustyGasTransport.h @@ -77,6 +77,16 @@ namespace Cantera { virtual void setParameters(const int type, const int k, const doublereal* const p); + + //! Return the Multicomponent diffusion coefficients. Units: [m^2/s]. + /*! + * Returns the array of multicomponent diffusion coefficients. + * + * @param ld The dimension of the inner loop of d (usually equal to m_nsp) + * @param d flat vector of diffusion coefficients, fortran ordering. + * d[ld*j+i] is the D_ij diffusion coefficient (the diffusion + * coefficient for species i due to species j). + */ virtual void getMultiDiffCoeffs(const int ld, doublereal* const d); //! Get the molar fluxes [kmol/m^2/s], given the thermodynamic @@ -168,6 +178,12 @@ namespace Cantera { * at. If the temperature is changed, update Booleans are set false, triggering recomputation. */ void updateTransport_T(); + + //! Update concentration-dependent quantities within the object + /*! + * The object keeps a value m_temp, which is the temperature at which quantities were last evaluated + * at. If the temperature is changed, update Booleans are set false, triggering recomputation. + */ void updateTransport_C(); void updateBinaryDiffCoeffs(); @@ -176,12 +192,16 @@ namespace Cantera { void eval_H_matrix(); - // gas attributes + //! Number of species in the gas phase int m_nsp; - doublereal m_tmin, m_tmax; + + //! Local copy of the species molecular weights + /*! + * units kg /kmol + * length = m_nsp; + */ vector_fp m_mw; - //! binary diffusion coefficients DenseMatrix m_d;