Unable to run simulation for Lactic Acid from biorefineries

[phiparth]

The following code :
from biorefineries import lactic as la la.load(kind='SSCF', print_results=True)

gives error :

Exception has occurred: ValueError
alias 'Lime' already in use by Chemical('Ash')

After several trial errors, the issue still persists or gives some other errors in biorefineries source code's system.py

Any solution or help would be greatly appreciated

[yalinli2]

@yoelcortes did you add this Lime alias when initializing biosteam for Ash? I'm not sure if this is good... Yeah we use CaO as the default ash, but adding this alias just doesn't seem right.

[phiparth]

@yoelcortes did you add this Lime alias when initializing biosteam for Ash? I'm not sure if this is good... Yeah we use CaO as the default ash, but adding this alias just doesn't seem right.

I used the default code from the folder at https://github.com/BioSTEAMDevelopmentGroup/Bioindustrial-Park/tree/f0579f602a13c4b33d21133c5fbe0bd8aaf9e030/biorefineries/lactic for initialization without changing anything within it.
Though if in the following line within _chemicals.py within biorefineries\lactic :
add_chemical('Ash', search_ID='CaO', phase='s', Hf=-151688*_cal2joule, HHV=0, LHV=0)
If the search_ID part is removed, the error alias 'Lime' already in use by Chemical('Ash') vanishes.

But the next error shows up :

type object 'Reactor' has no attribute '_BM'

After trying some edits within units.py in class Reactor, the error still persists.

[yoelcortes]

@phiparth, thank you for all the details and debugging efforts (very helpful!). We will take care of it hopefully soon.

@yalinli2,

Thanks for verifying this with me. BioSTEAM now automatically adds aliases to chemicals when they are compiled based on the 'common_name' and 'formula' attributes of a chemical. This introduces some issues in some of the old biorefineries, but they are not hard to fix. Instead of this:

from biosteam import Chemical, Chemicals
_cal2joule = 4.184
# This is effectively what happens in the lactic acid chemicals.py module
chemicals = Chemicals([])
chemicals .append(
    Chemical('Ash', search_ID='CaO', phase='s', Hf=-151688*_cal2joule, HHV=0, LHV=0, CAS='Ash')
)
...
chemicals .compile() 

You can do this:

from biosteam import Chemical, Chemicals
_cal2joule = 4.184
chemicals = Chemicals([])
dummy = Chemical('CaO', phase='s', Hf=-151688*_cal2joule, HHV=0, LHV=0)
chemicals .append(CaO.copy(ID='Ash')) # Use copy instead
...
chemicals.compile() 

Would you like me to help update the lactic biorefinery to get it running (not just for this bug)?

Thanks,

[phiparth]

@yoelcortes Thank you for the troubleshoot with the problem.

After replacing the code with the specified snippet, another error pops up

Exception has occurred: UndefinedChemicalAlias
'Ash'
KeyError: ('Glucan', 'Xylan', 'Lignin', 'Ash', 'Acetate', 'Protein', 'Arabinan', 'Galactan', 'Mannan', 'Sucrose', 'Extractives', 'SuccinicAcid')

During handling of the above exception, another exception occurred:

KeyError: 'Ash'

Also there are other errors like in the Reactor class.

It would be a great help if you could help update the lactic biorefinery to get it running

[yoelcortes]

We'll get it running by the end of the week and notify you as soon as we do.

Thanks,

[yalinli2]

thanks @yoelcortes for the suggestions, it was an easy fix (with your help).

@phiparth if you pull the commit BioSTEAMDevelopmentGroup/Bioindustrial-Park@2edf265, you should be able to run the module, thanks for reporting.

[phiparth]

@yalinli2 Thankyou for the help in troubleshooting, but the change seems to create yet another error :

Exception has occurred: UndefinedChemicalAlias
'Ash'
ModuleNotFoundError: No module named 'biorefineries.lactic.load'

During handling of the above exception, another exception occurred:

KeyError: ('Glucan', 'Xylan', 'Lignin', 'Ash', 'Acetate', 'Protein', 'Arabinan', 'Galactan', 'Mannan', 'Sucrose', 'Extractives', 'SuccinicAcid')

During handling of the above exception, another exception occurred:

KeyError: 'Ash'

After trying several changes within _chemical.py, the error still persists.

[yalinli2]

huh interesting, @yoelcortes can you see if the following work on your end? I thought it worked on my end but not sure if any changes weren't pushed to remote. Thanks!

from biorefineries import lactic as la
la.load(kind='SSCF', print_results=True)

[yoelcortes]

@yalinli2, @phiparth,

It is working on my side. After pulling the github repo and running the simulation, I get the following:

---------- Simulation Results ----------
MPSP is $1.390/kg
GWP is 4.296 kg CO2-eq/kg lactic acid
FEC is 53.52 MJ/kg lactic acid
------------------------------------------

@phiparth, please use the latest biosteam versions on PyPI and the biorefineries version github. Let me know if it still does not work.

Thanks,

[phiparth]

@yoelcortes, The 'SSCF' configuration works fine giving the exact simulation results but on switching to 'SHF' configuration, the error occurs :

Exception has occurred: AssertionError
�[38;2;229;63;229m<HeatExchangerNetwork: HXN>�[0m
Not equal to tolerance rtol=0.05, atol=0.001

Mismatched elements: 1 / 1 (100%)
Max absolute difference: 26.932
Max relative difference: 0.071
x: array(352.06808)
y: array(378.999999)

Also the versions are updated along I have used the updated version of biorefineries from github.

Additionally if I try to run the 1_uncertainities.py under the analyses folder for SSCF configuration, another error occurs :

in create_model
else: lactic_sys = flowsheet.system.lactic_sys
AttributeError: 'str' object has no attribute 'system'

[yoelcortes]

@yalinli2,

The SHF biorefinery was not working be cause the pump outlet pressure does not default to 101325 anymore. The low pressures was messing with the heat exchanger network. I took the liberty to add P=101325 to pumps in the lactic acid code and it works now:

---------- Simulation Results ----------
MPSP is $1.864/kg
GWP is 7.323 kg CO2-eq/kg lactic acid
FEC is 100.39 MJ/kg lactic acid
------------------------------------------

@phiparth, thanks for submitting the issue. It should work now.

Thanks!

[phiparth]

@yoelcortes, thank you for the help with SCF configuration, but the problem with analyses part remains.
When I run the code for 1_uncertainities.py for the SSCF configuration, the error occurs :

in create_model
else: lactic_sys = flowsheet.system.lactic_sys
AttributeError: 'str' object has no attribute 'system'

Whereas for the running it for SCF configuration gives :

ValueError: <System: separation_sys> <System: SYS2> <Esterification: R402> assignment destination is read-only

[yoelcortes]

@yalinli2, @phiparth

I made a minor update for compatibility with new and old BioSTEAM versions (the TicToc class was renamed to Timer; to not confuse any new users jaja).

I will be out on business travel until Friday. But I will get back to helping afterwards.

Thanks,

[phiparth]

@yoelcortes
Thanks for the help, but the issue in 1_uncertainities.py under analyses folder remains, possibly due to some attribute error

For SSCF:

in create_model
else: lactic_sys = flowsheet.system.lactic_sys
AttributeError: 'str' object has no attribute 'system'

For SHF :

ValueError: <System: separation_sys> <System: SYS2> <Esterification: R402> assignment destination is read-only