github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-22 20:38:01 -05:00
Code Issues Packages Projects Releases Wiki Activity

Add some logic for mixtures to use "normal" TP solver when phase is imposed (not the global search)

Browse Source
This commit is contained in:
Ian Bell
2016-08-22 22:18:49 -06:00
parent 43ba3b9071
commit de455acf10
1 changed files with 17 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/Backends/Helmholtz/FlashRoutines.cpp
View File

@@ -55,17 +55,27 @@ void FlashRoutines::PT_flash_mixtures(HelmholtzEOSMixtureBackend &HEOS)
}
}
else{
// Following the strategy of Gernert, 2014
StabilityRoutines::StabilityEvaluationClass stability_tester(HEOS);
if (!stability_tester.is_stable()){
throw CoolProp::ValueError("PT_INPUTS are two-phase, not yet supported");
if (HEOS.imposed_phase_index == iphase_not_imposed){
// Blind flash call
// Following the strategy of Gernert, 2014
StabilityRoutines::StabilityEvaluationClass stability_tester(HEOS);
if (!stability_tester.is_stable()){
throw CoolProp::ValueError("PT_INPUTS are two-phase, not yet supported");
}
else{
// It's single-phase
double rho = HEOS.solver_rho_Tp_global(HEOS.T(), HEOS.p(), 20000);
HEOS.update_DmolarT_direct(rho, HEOS.T());
HEOS._Q = -1;
HEOS._phase = iphase_liquid;
}
}
else{
// It's single-phase
double rho = HEOS.solver_rho_Tp_global(HEOS.T(), HEOS.p(), 20000);
// It's single-phase, and phase is imposed
double rho = HEOS.solver_rho_Tp(HEOS.T(), HEOS.p());
HEOS.update_DmolarT_direct(rho, HEOS.T());
HEOS._Q = -1;
HEOS._phase = iphase_liquid;
HEOS._phase = HEOS.imposed_phase_index;
}
}
}