From 88dc7a31c958a3c1f2d2b76650cc6af44c672e94 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Thu, 23 Oct 2008 14:55:54 +0000 Subject: [PATCH] Added the voltage term onto vcs_dfe calculation of single phase chemical potentials. It's all currently zero, but I can't figure out why it should be left off. --- Cantera/src/equil/vcs_solve_TP.cpp | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/Cantera/src/equil/vcs_solve_TP.cpp b/Cantera/src/equil/vcs_solve_TP.cpp index 3ad233746..673a94b4d 100644 --- a/Cantera/src/equil/vcs_solve_TP.cpp +++ b/Cantera/src/equil/vcs_solve_TP.cpp @@ -4255,7 +4255,7 @@ namespace VCSnonideal { m_SSfeSpecies[kspec] + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { if (m_SSPhase[kspec]) { - feSpecies[kspec] = m_SSfeSpecies[kspec]; + feSpecies[kspec] = m_SSfeSpecies[kspec] + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { if (molNum[kspec] <= VCS_DELETE_MINORSPECIES_CUTOFF) { iph = m_phaseID[kspec]; @@ -4296,12 +4296,11 @@ namespace VCSnonideal { std::exit(-1); } #endif - feSpecies[kspec] = - m_SSfeSpecies[kspec] + feSpecies[kspec] = m_SSfeSpecies[kspec] + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { if (m_SSPhase[kspec]) { - feSpecies[kspec] = m_SSfeSpecies[kspec]; + feSpecies[kspec] = m_SSfeSpecies[kspec]+ m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { if (molNum[kspec] <= VCS_DELETE_MINORSPECIES_CUTOFF) { iph = m_phaseID[kspec]; @@ -4348,7 +4347,7 @@ namespace VCSnonideal { m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; ; } else { if (m_SSPhase[kspec]) { - feSpecies[kspec] = m_SSfeSpecies[kspec]; + feSpecies[kspec] = m_SSfeSpecies[kspec]+ m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { if (molNum[kspec] <= VCS_DELETE_MINORSPECIES_CUTOFF) { iph = m_phaseID[kspec];