research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Equil] remove unused VCS_SOLVE::reportCSV

Browse source
This commit is contained in:
Ray Speth 2017-08-21 18:58:13 -04:00
parent 08726b0904
commit da801f4cbc
2 changed files with 0 additions and 128 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
include/cantera/equil/vcs_prob.h
View file

@ -257,8 +257,6 @@ public:
*/
size_t addOnePhaseSpecies(vcs_VolPhase* volPhase, size_t k, size_t kT);
void reportCSV(const std::string& reportFile);
//! Set the debug level
/*!
* @param vcs_debug_print_lvl input debug level

126
src/equil/vcs_prob.cpp
View file

@ -289,132 +289,6 @@ size_t VCS_PROB::addOnePhaseSpecies(vcs_VolPhase* volPhase, size_t k, size_t kT)
return kT;
}
void VCS_PROB::reportCSV(const std::string& reportFile)
{
FILE* FP = fopen(reportFile.c_str(), "w");
if (!FP) {
throw CanteraError("VCS_PROB::reportCSV", "Failure to open file");
}
vector_fp volPM(nspecies, 0.0);
vector_fp activity(nspecies, 0.0);
vector_fp ac(nspecies, 0.0);
vector_fp mu(nspecies, 0.0);
vector_fp mu0(nspecies, 0.0);
vector_fp molalities(nspecies, 0.0);
double vol = 0.0;
size_t iK = 0;
for (size_t iphase = 0; iphase < NPhase; iphase++) {
size_t istart = iK;
vcs_VolPhase* volP = VPhaseList[iphase];
size_t nSpeciesPhase = volP->nSpecies();
volPM.resize(nSpeciesPhase, 0.0);
volP->sendToVCS_VolPM(&volPM[0]);
double TMolesPhase = volP->totalMoles();
double VolPhaseVolumes = 0.0;
for (size_t k = 0; k < nSpeciesPhase; k++) {
iK++;
VolPhaseVolumes += volPM[istart + k] * mf[istart + k];
}
VolPhaseVolumes *= TMolesPhase;
vol += VolPhaseVolumes;
}
fprintf(FP,"--------------------- VCS_MULTIPHASE_EQUIL FINAL REPORT"
" -----------------------------\n");
fprintf(FP,"Temperature = %11.5g kelvin\n", T);
fprintf(FP,"Pressure = %11.5g Pascal\n", PresPA);
fprintf(FP,"Total Volume = %11.5g m**3\n", vol);
fprintf(FP,"Number Basis optimizations = %d\n", m_NumBasisOptimizations);
fprintf(FP,"Number VCS iterations = %d\n", m_Iterations);
iK = 0;
for (size_t iphase = 0; iphase < NPhase; iphase++) {
size_t istart = iK;
vcs_VolPhase* volP = VPhaseList[iphase];
const ThermoPhase* tp = volP->ptrThermoPhase();
string phaseName = volP->PhaseName;
size_t nSpeciesPhase = volP->nSpecies();
volP->sendToVCS_VolPM(&volPM[0]);
double TMolesPhase = volP->totalMoles();
activity.resize(nSpeciesPhase, 0.0);
ac.resize(nSpeciesPhase, 0.0);
mu0.resize(nSpeciesPhase, 0.0);
mu.resize(nSpeciesPhase, 0.0);
volPM.resize(nSpeciesPhase, 0.0);
molalities.resize(nSpeciesPhase, 0.0);
int actConvention = tp->activityConvention();
tp->getActivities(&activity[0]);
tp->getActivityCoefficients(&ac[0]);
tp->getStandardChemPotentials(&mu0[0]);
tp->getPartialMolarVolumes(&volPM[0]);
tp->getChemPotentials(&mu[0]);
double VolPhaseVolumes = 0.0;
for (size_t k = 0; k < nSpeciesPhase; k++) {
VolPhaseVolumes += volPM[k] * mf[istart + k];
}
VolPhaseVolumes *= TMolesPhase;
vol += VolPhaseVolumes;
if (actConvention == 1) {
const MolalityVPSSTP* mTP = static_cast<const MolalityVPSSTP*>(tp);
tp->getChemPotentials(&mu[0]);
mTP->getMolalities(&molalities[0]);
tp->getChemPotentials(&mu[0]);
if (iphase == 0) {
fprintf(FP," Name, Phase, PhaseMoles, Mole_Fract, "
"Molalities, ActCoeff, Activity,"
"ChemPot_SS0, ChemPot, mole_num, PMVol, Phase_Volume\n");
fprintf(FP," , , (kmol), , "
" , , ,"
" (J/kmol), (J/kmol), (kmol), (m**3/kmol), (m**3)\n");
}
for (size_t k = 0; k < nSpeciesPhase; k++) {
std::string sName = tp->speciesName(k);
fprintf(FP,"%12s, %11s, %11.3e, %11.3e, %11.3e, %11.3e, %11.3e,"
"%11.3e, %11.3e, %11.3e, %11.3e, %11.3e\n",
sName.c_str(),
phaseName.c_str(), TMolesPhase,
mf[istart + k], molalities[k], ac[k], activity[k],
mu0[k]*1.0E-6, mu[k]*1.0E-6,
mf[istart + k] * TMolesPhase,
volPM[k], VolPhaseVolumes);
}
} else {
if (iphase == 0) {
fprintf(FP," Name, Phase, PhaseMoles, Mole_Fract, "
"Molalities, ActCoeff, Activity,"
" ChemPotSS0, ChemPot, mole_num, PMVol, Phase_Volume\n");
fprintf(FP," , , (kmol), , "
" , , ,"
" (J/kmol), (J/kmol), (kmol), (m**3/kmol), (m**3)\n");
}
for (size_t k = 0; k < nSpeciesPhase; k++) {
molalities[k] = 0.0;
}
for (size_t k = 0; k < nSpeciesPhase; k++) {
std::string sName = tp->speciesName(k);
fprintf(FP,"%12s, %11s, %11.3e, %11.3e, %11.3e, %11.3e, %11.3e, "
"%11.3e, %11.3e,% 11.3e, %11.3e, %11.3e\n",
sName.c_str(),
phaseName.c_str(), TMolesPhase,
mf[istart + k], molalities[k], ac[k],
activity[k], mu0[k]*1.0E-6, mu[k]*1.0E-6,
mf[istart + k] * TMolesPhase,
volPM[k], VolPhaseVolumes);
}
}
iK += nSpeciesPhase;
}
fclose(FP);
}
void VCS_PROB::setDebugPrintLvl(int lvl)
{
vcs_debug_print_lvl = lvl;