Thermo v. 0.2.26 flash_vl - Raising Exception on PH flash about thermo

[CalebBell]

Pure components should work fine but instead of FlashVL you need to use FlashPureVLS: https://thermo.readthedocs.io/thermo.flash.html#thermo.flash.FlashPureVLS

[CalebBell]

Hi @mikeyj777,

The legacy Stream interface doesn't use the same enthalpy basis as the new Phase/flasher-based stuff, would not recommend trying to use them together. EquilibriumStream is the replacement for Stream. I have put off documenting it & don't have a tutorial; however, to combine to EquilibriumStream objects, you would add their flow rates (.ns attribute) and energy (.energy atribute) and provide them as input to a new EquilibriumStream object, as well as the pressure & the flasher object.

Here's a brief example

from thermo import *
constants, properties = ChemicalConstantsPackage.from_IDs(['methane', 'ethane', 'nitrogen'])
eos_kwargs = {'Pcs': constants.Pcs, 'Tcs': constants.Tcs, 'omegas': constants.omegas}
gas = CEOSGas(PRMIX, eos_kwargs=eos_kwargs, HeatCapacityGases=properties.HeatCapacityGases)
liquid = CEOSLiquid(PRMIX, eos_kwargs=eos_kwargs, HeatCapacityGases=properties.HeatCapacityGases)
flasher = FlashVL(constants, properties, liquid=liquid, gas=gas)
S1 = EquilibriumStream(flasher=flasher, T=300, P=1e5, ns=[1,2,3])
S2 = EquilibriumStream(flasher=flasher, T=305, P=1e5, ns=[.5,.3,.8])

ns_out = [na+nb for na, nb in zip(S1.ns, S2.ns)]
S3 = EquilibriumStream(flasher=flasher, P=1e5, ns=ns_out, energy=S1.energy+S2.energy)
S3.T, S3.ns

(300.9985275567583, [1.5, 2.2999999999999994, 3.7999999999999994])

Sincerely,
Caleb

[mikeyj777]

thank you so much for the great support! I will give it a go and comment back with any issues.

[mikeyj777]

How would you recommend working with a pure component stream? the equilibrium stream failed with a div by zero error in property_package.py, line 455