TEA not working in BioSTEAM #174

Uzair1shah · 2023-10-02T18:34:49Z

Uzair1shah
Oct 2, 2023

Hello

I have tried to run TEA for the preprocessing of a biorefinery using the same code provided in the documentation, however it shows the following error:

*TypeError: unsupported operand type(s) for : 'property' and 'float'

I have tried to fix this error, however I am unable to solve. The code is given as follows:

S101 = Stream('S101', biomass=333334.2, units='kg/hr', phase='s', price=price['biomass'])
S102 = Stream('S102', phase='s')
UO101 = TruckScale('UO101', ins=S101, outs=S102)                                     		# Truck scale

# Simulation of UO101 (Truck scale)
UO101.simulate()
UO101.show( flow='kg/hr', composition = False, N=9)
r1 = UO101.results()
UO101.diagram()
print(r1)


# Truckdumper UOP
S103 = Stream('S103', phase='s')
UO102 = TruckDumper('UO102', ins=S102, outs=S103)                                             # Truck dumper

# Simulation of UO102 (Truck dumper)
UO102.simulate()
UO102.show(flow='kg/hr', composition = False, N=9)
r2 = UO102.results()
UO102.diagram()
print(r2)



# Truck tipper receiving hopper
S104 = Stream('S104', phase='s')
UO103 = DumpHopper('UO103', ins=S103, outs=S104)                                     		# Dump hopper

# Simulation of UO102 (Truck tipper)
UO103.simulate()
UO103.show(flow='kg/hr' , composition = False, N=9)
r3 = UO103.results()
UO103.diagram()
print(r3)

# Vibrating Screen
S105 = Stream('S105', phase='s')
S106 = Stream('S106', phase='s')

UO104 = VibratingScreen('UO104', ins=S104, outs=(S105,S106), split=0.9)                     # Vibrating Screen


# Simulation of UO103 (Vibrating Screen)
UO104.simulate()
UO104.show(flow='kg/hr', composition = False, N=9)
r4 = UO104.results()
UO104.diagram()
print(r4)

                                                                                            # Rotary Dryer
S107 = Stream('S107', phase='g') #Air
S108 = Stream('S108', phase='g') #Natural gas
S109 = Stream('S109', phase='s') #dry feedstock/solids
S110 = Stream('S110', phase='g') #hot air
S111 = Stream('S111', phase='g') #emissions
UO105 = units.DrumDryer('UO105', ins=(S105,S107,S108), outs=(S109,S110,S111),moisture_content = 0.01, split = (dict(moisture=0.9)) )

# Simulation of UO105 (Rotary dryer)
UO105.simulate()
UO105.show(flow='kg/hr', composition = False, N=10)
r5 = UO105.results()
UO105.diagram()
print(r5)

S113 = Stream('S113', phase='s') #recycle stream
S114 = Stream('S114', phase='s')
UO106 = units.Mixer('UO106', ins=(S109, S113), outs=(S114))

# Hammermill
S112 = Stream('S112', phase='s') #crushed dry feedstock
S115 = Stream('S115', phase='s') 

UO107 = units.HammerMill('UO107', ins=(S114), outs=(S115))

#Vibrating Screen
S116 = Stream('S116', phase='s')
UO108 = VibratingScreen('UO108', ins=S115, outs=(S113,S116), split=0.1) 

Preprocessing = bst.main_flowsheet.create_system('Preprocessing')
Preprocessing.simulate()
Preprocessing.diagram('cluster', format='png')

UO106.show(flow='kg/hr', composition = False, N=9)
r6 = UO106.results()
print(r6)

UO107.show(flow='kg/hr', composition = False, N=9)
r7 = UO107.results()
print(r7)

UO108.show(flow='kg/hr', composition = False, N=9)
r8 = UO108.results()
print(r8)

class Preprocessing(bst.TEA):

    def __init__(self, system, IRR, duration, depreciation, income_tax,
                 operating_days, lang_factor, construction_schedule, WC_over_FCI,
                 labor_cost, fringe_benefits, property_tax,
                 property_insurance, supplies, maintenance, administration):
 
        super().__init__(system, IRR, duration, depreciation, income_tax,
                         operating_days, lang_factor, construction_schedule,
                         startup_months=0, startup_FOCfrac=0, startup_VOCfrac=0,
                         startup_salesfrac=0, finance_interest=0, finance_years=0,
                         finance_fraction=0, WC_over_FCI=WC_over_FCI)
        self.labor_cost = labor_cost
        self.fringe_benefits = fringe_benefits
        self.property_tax = property_tax
        self.property_insurance = property_insurance
        self.supplies= supplies
        self.maintenance = maintenance
        self.administration = administration


    def _DPI(self, installed_equipment_cost):
        return installed_equipment_cost


    def _TDC(self, DPI):
        return DPI

 
    def _FCI(self, TDC):
        return TDC


    def _FOC(self, FCI):
        return (FCI*(self.property_tax + self.property_insurance
                     + self.maintenance + self.administration)
                + self.labor_cost*(1+self.fringe_benefits+self.supplies))


tea =  Preprocessing(system=Preprocessing,
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
                  property_insurance=0.11,
                  supplies=0.20,
                  maintenance=0.01,
                  administration=0.005)

tea.show()
tea.get_cashflow_table()

Note: This code would not be able to run it because different unit operations are used other than available in the repository. I am new to using BioSTEAM, so I am not sure if there is any error in my code or the TEA file of BioSTEAM. I would really appreciate if you can help me in this regard. Thanks a lot for always helping.

Best,
Uzair

Answered by yoelcortes

Oct 3, 2023

@Uzair1shah, OK! I was able to reproduce the error. The issue is the first argument to a TEA object should be a System object not a TEA subclass:

sys = bst.main_flowsheet.create_system('sys')
tea =  Preprocessing(system=sys, # Not Preprocessing
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
…

View full answer

yoelcortes · 2023-10-03T13:16:07Z

yoelcortes
Oct 3, 2023
Maintainer

@Uzair1shah,

The code you posted does not run for many reasons (e.g., biosteam is not imported, there is no thermodynamic property package). Please include all the code required to reproduce the error.

Can you also post the following:

BioSTEAM version (pip show biosteam).
Python version
Full traceback of the error

It is helpful to remove any irrelevant code too. For example, if the error happens at an earlier line, remove any lines that comes next.

Thanks!

1 reply

Uzair1shah Oct 3, 2023
Author

@yoelcortes thanks for the response. The whole code for reproducing the error is as follows:

import biosteam as bst
import thermosteam as tmo
import graphviz
from biosteam import main_flowsheet as F, units, Stream, settings
import math as math

from biosteam.units import TruckScale, TruckDumper, DumpHopper, VibratingScreen, ScreenerFeeder, StackerReclaimer, ConveyingBelt, HammerMill, VibratingScreen, Mixer, Splitter

import functools

species = {}		#create an empty collection to hold your chemical species
  
#Biomass cellulose
species['cellulose'] = bst.Chemical('cellulose', phase='s', search_db=False, Cp=1.364, rho=1540, formula="C6H10O5", Hf=-1019000, MW =162.1424)            			#Cellulose, verified
species['cellA'] = bst.Chemical('cellA', search_db=False, default=True, phase='s', formula ="C6H10O5", Hf=-1019000, MW=162.1424)                                   	#Activated cellulose
  
# Biomass hemicellulose
species['gmsw'] = bst.Chemical('gmsw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['xyhw'] = bst.Chemical('xyhw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['hce1'] = bst.Chemical('hce1', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified
species['hce2'] = bst.Chemical('hce2', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified

# Biomass lignin
species['ligc'] = bst.Chemical('ligc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                    					#Verified
species['ligh'] = bst.Chemical('ligh', search_db=False, phase='s', formula="C22H28O9",MW=436.45892, Hf=-1722700)                                   					#Verified
species['ligO'] = bst.Chemical('ligO', search_db=False, phase='s', formula="C20H22O10",MW=422.38868, Hf=-1847500)                                  					#Verified
species['lig'] = bst.Chemical('lig', search_db=False, phase='s', formula="C11H12O4",MW=208.21388, Hf=-729310)                                     					#Secondary lignin intermediate, verified
species['ligcc'] = bst.Chemical('ligcc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                   				#Verified
species['ligoh'] = bst.Chemical('ligoh', search_db=False, phase='s', formula="C19H22O8",MW=378.37888, Hf=-1429200)                                  				#Verified

# Biomass extractives
species['tannin'] = bst.Chemical('tannin', search_db=False, phase='s', formula="C15H12O7", MW=304.25608,Hf=-1079700)                                 				#Tannic acid formula, verified
species['itann'] = bst.Chemical('itann', search_db=False, phase='s', formula="C8H4O4",MW=164.11736, Hf=-616980)                                     				#Intermediate tannin, verified
species['triglycerides'] = bst.Chemical('triglycerides', search_db=False, phase='s', formula="C76H52O46", MW=897.4168, Hf=-1546300)                          		#Triglycerides formula, update the triglycerides Hf

# Biomass inorganics
# This could be defined more explictly later
species['ash'] = bst.Chemical('CaO', phase='s', formula="CaO", Hf=5382027.245)                                         												#Ash, verify HF

# Pyrolysis solid products

species['char'] = bst.Chemical('carbon', phase='s', formula="C")        																							#Carbon, verified, properties from refs 1,2
species['hmwl'] = bst.Chemical('hmwl', search_db=False, phase='l', formula="C24H28O4", MW=380.48392, Hf=-958260)                                   					#Verified, High molecular weight lignin

# Fluidized bed media
species['sand'] = bst.Chemical('SiO2', formula="SiO2", Hf=-910860)                                            														#Sand, inert, properties from ref2

#metaplastic species
species['gco2'] = bst.Chemical('gco2',search_db=False, phase='s', formula="CO2", MW=44.01, Hf=-393510)                                            					#Trapped carbon dioxide, properties from ref2
species['gco'] = bst.Chemical('gco',search_db=False, phase='s', formula="CO",MW=28.01, Hf=-110520)                                               					#Trapped carbon monoxide, properties from ref2
species['gmeoh'] = bst.Chemical('gmeoh',search_db=False, phase='s', formula="CH3OH", MW=32.042, Hf=-238640)                                         				#Trapped methanol, verified, properties from ref2
species['gcoh2_l'] =bst.Chemical('gcoh2_l', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)                                        				#Trapped formaldehyde, properties from ref2, loose
species['gc2h6'] = bst.Chemical('gc2h6', search_db=False, phase='s', formula="C2H6", MW=30.07, Hf=-84667)                                           				#Trapped ethane, properties from ref2
species['gch4'] = bst.Chemical('gch4', search_db=False, phase='s', formula="CH4",MW=16.041, Hf=-74840)                                             					#Trapped methane, properties from ref2
species['gc2h4'] = bst.Chemical('gc2h4', search_db=False, phase='s', formula="C2H4", MW=28.052, Hf=52283)                                            				#Trapped ethylene, properties from ref2
species['gc6h5oh'] = bst.Chemical('gc6h5oh', search_db=False, phase='s', formula="C6H6O", MW=94.113, Hf=-96400)                                       				#Trapped phenol, properties from ref3
species['gcoh2_s'] = bst.Chemical('gcoh2_s', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)	                                      				#Trapped formaldehyde, properties from ref2, stiff
species['gh2'] = bst.Chemical('gh2', search_db=False, phase='s', formula="H2",MW=2.016, Hf=0)		                                              					#Trapped hydrogen, properties from ref2

#non-biomass species
species['argon'] = bst.Chemical('argon', formula="Ar")  
species['nitrogen'] = bst.Chemical('nitrogen', formula="N2")  
species['oxygen'] = bst.Chemical('oxygen', formula="O2")  

#Pyrolysis products
## Gas phase species
species['ethylene'] = bst.Chemical('ethylene', formula="C2H4", search_ID='ethylene', phase='g')                  
species['methane'] = bst.Chemical('methane', formula="CH4", search_ID='methane',  phase='g')                                   										#ref2
species['hydrogen'] = bst.Chemical('hydrogen', formula="H2", search_ID='hydrogen',  phase='g')                                     									#ref2
species['cacrbonmonoxide'] = bst.Chemical('carbon monoxide', formula="CO", search_ID='carbon monoxide', phase='g')                               					#ref2
species['carbondioxide'] = bst.Chemical('carbon dioxide', formula="CO2", search_ID='carbon dioxide', phase='g')                            							#ref2
species['formaldehyde'] = bst.Chemical('formaldehyde', formula="CH2O", search_ID='formaldehyde', phase='g')                            								#ref2

## Liquid phase species
species['acrolein'] = bst.Chemical('acrolein', formula="C3H4O", search_ID='acrolein', phase='l')      										#https://atct.anl.gov/Thermochemical%20Data/version%201.122e/species/?species_number=1209

species['sinapylaldehyde'] = bst.Chemical('sinapyaldehyde', search_db=False, phase='l', formula="C11H12O4", MW=208.21388, Hf=-483800)
species['methyllinoleate'] = bst.Chemical('methyllinoleate', phase='l', search_db=False, formula="C19H34O2", MW=294.4721, Hf=-604880)
species['glyoxal'] = bst.Chemical('glyoxal', search_ID='glyoxal', phase='l')							    														#NIST webbook
species['acetaldehyde'] = bst.Chemical('acetaldehyde', search_ID='acetaldehyde', phase='l')                       													#Acetaldehyde, properties from ref2 
#species['aceticacid'] = bst.Chemical('aceticacid', search_db=False, default=True, formula="CH3COOH", MW=60.05, Hf=-484930, Hvap=0, Psat=0, Tb=0)                   #Acetic acid, properties from ref2 
species['glycolaldehyde'] = bst.Chemical('glycolaldehyde', search_ID='glycolaldehyde', phase='l')                                    								#Source Pyrolysis of the Simplest Carbohydrate, Glycolaldehyde (CHO-CH2OH), and Glyoxal in a Heated Microreactor. Find additional sources
species['ethanol'] = bst.Chemical('ethanol', search_ID='ethanol', phase='l')                             															#Ethanol, properties from ref2 
species['propionaldehyde'] = bst.Chemical('propionaldehyde', search_ID='propionaldehyde', phase='l')                    							#Source: NIST webbook https://webbook.nist.gov/cgi/cbook.cgi?ID=C123386&Mask=2 
species['hydroxypropanal'] = bst.Chemical('hydroxypropanal', search_ID='hydroxypropanal',phase='l')  												#https://www.chemeo.com/cid/20-360-9/2%2C2-Dimethyl-3-hydroxypropionaldehyde
species['furfural'] = bst.Chemical('furfural', search_ID='furfural', phase='l')
species['xylosan'] = bst.Chemical('xylosan', phase='l', search_db=False, formula="C5H8O4", MW=132.115, Hf=-647000)                             				#University of Pittsburg Pitt Quantum Repository, DOI 10.17614/Q4SQ8QT7Q
species['levoglucosan'] = bst.Chemical('levoglucosan', search_ID='levoglucosan', phase='l')                      											#Gas phase Hf value, ref4.
species['phenol'] = bst.Chemical('phenol', search_ID='phenol', phase='l')                                 													#phenol, properties from ref5
species['hydroxymethylfurfural'] = bst.Chemical('hydroxymethylfurfural', search_db=False, formula="C6H6O3", MW=126.1100, Hf=-437400, phase='l')           	#Hydromethyl furfural DOI:10.1021/ie901012g
species['anisole'] = bst.Chemical('anisole', search_ID='anisole', phase='l')                                  												#DOI: 10.1021/jp507267f
species['pcoumarylalcohol'] = bst.Chemical('pcoumarylalcohol', search_db=False, formula="C9H10O2", MW=150.1772, Hf=-193500, phase='l')                 		#Unverified
species['formicacid'] = bst.Chemical('formicacid', search_db=False, formula="CH2O2", MW=46.0254, Hf=-431650, phase='l')                           			#ref: Heat of formation of liquid. Should it be of gas?
species['methanol'] = bst.Chemical('methanol', search_ID='methanol', phase='l')                               												#ref2
species['water'] = bst.Chemical('water', search_ID='water', phase='l')                                   													#ref2

species['propaniocacid'] = bst.Chemical('propaniocacid', search_db=False, formula="C3H6O2", Hf=-199481.24962000002, phase='l')
species['linalylpropioinate'] = bst.Chemical('linalylpropioinate', search_db=False, formula="C19H34O2", Hf=-194296.922694, phase='l')
species['vanillin'] = bst.Chemical('vanillin', search_ID='vanillin', phase='l')
species['cresol'] = bst.Chemical('cresol', search_db=False, formula="C7H8O", Hf=-34212.624394, phase='l')

# Create an empty Chemicals object and add each species to it
chemicals = tmo.Chemicals([])

for key, value in species.items():
    #print(key)
    value.default() # call the default method to populate all other properties not defined explicitly. Todo: some properties default to water's properties. Check individual chemical properties for accuracy once model is developed. 
    chemicals.append(value) # push a chemical object to back of the chemicals collection object

chemicals.compile()								# compile the loaded chemicals. Must do.
chemicals.set_synonym('water', 'moisture')		# set synonyms or common name for select species
chemicals.set_synonym('CaO', 'ash')

settings.set_thermo(chemicals, cache=False)				# Todo: read more set_thermo method.


# set raw material price (USD/kg)
price = {'biomass': 0.08, # @10% moisture content wet basis
         'water': 0.000353,
         'steam': 0.017,
         'waste': -0.33,
         'gasoline': 0.756} # 2 USD/gal

#setting electructy price
settings.electricity_price = 0.1           # set electricity price to your desired rate

# set the name of your flowsheet
bst.main_flowsheet.set_flowsheet('fast_pyrolysis_biorefinery')

# define a chemical group with its composition in mass fraction here
chemicals.define_group('biomass',IDs = ['cellulose','xyhw','ligc','ligh','ligO','triglycerides','tannin','moisture','ash'], composition = [0.2347,0.1857,0.0915,0.0761,0.1290,0.0386,0.0176,0.10,0.0915], wt=True)

# Truckscale UOP
# S101 - Biomass inlet into the truck scale
# S102 - Biomass outlet out of the truck scale
S101 = Stream('S101', biomass=333334.2, units='kg/hr', phase='s', price=price['biomass'])
S102 = Stream('S102', phase='s')
UO101 = TruckScale('UO101', ins=S101, outs=S102)                                     		# Truck scale

# Simulation of UO101 (Truck scale)
UO101.simulate()
UO101.show( flow='kg/hr', composition = False, N=9)
r1 = UO101.results()
print(r1)



# Vibrating Screen
S105 = Stream('S105', phase='s')
S106 = Stream('S106', phase='s')

UO104 = VibratingScreen('UO104', ins=S102, outs=(S105,S106), split=0.9)                     # Vibrating Screen


# Simulation of UO103 (Vibrating Screen)
UO104.simulate()
UO104.show(flow='kg/hr', composition = False, N=9)
r4 = UO104.results()
print(r4)
                                                                                           # Rotary Dryer
S107 = Stream('S107', phase='g') #Air
S108 = Stream('S108', phase='g') #Natural gas
S109 = Stream('S109', phase='s') #dry feedstock/solids
S110 = Stream('S110', phase='g') #hot air
S111 = Stream('S111', phase='g') #emissions
UO105 = units.DrumDryer('UO105', ins=(S105,S107,S108), outs=(S109,S110,S111),moisture_content = 0.01, split = (dict(moisture=0.9)) )
# Simulation of UO105 (Rotary dryer)
UO105.simulate()
UO105.show(flow='kg/hr', composition = False, N=10)
r5 = UO105.results()
print(r5)

S113 = Stream('S113', phase='s') #recycle stream
S114 = Stream('S114', phase='s')
UO106 = units.Mixer('UO106', ins=(S109, S113), outs=(S114))

# Hammermill
S112 = Stream('S112', phase='s') #crushed dry feedstock
S115 = Stream('S115', phase='s') 

UO107 = units.HammerMill('UO107', ins=(S114), outs=(S115))

#Vibrating Screen
S116 = Stream('S116', phase='s')
UO108 = VibratingScreen('UO108', ins=S115, outs=(S113,S116), split=0.1) 

Preprocessing = bst.main_flowsheet.create_system('Preprocessing')
Preprocessing.simulate()
Preprocessing.diagram('cluster', format='png')

UO106.show(flow='kg/hr', composition = False, N=9)
r6 = UO106.results()
print(r6)

UO107.show(flow='kg/hr', composition = False, N=9)
r7 = UO107.results()
print(r7)

UO108.show(flow='kg/hr', composition = False, N=9)
r8 = UO108.results()
print(r8)

class Preprocessing(bst.TEA):

    def __init__(self, system, IRR, duration, depreciation, income_tax,
                 operating_days, lang_factor, construction_schedule, WC_over_FCI,
                 labor_cost, fringe_benefits, property_tax,
                 property_insurance, supplies, maintenance, administration):
        # Huang et. al. does not take into account financing or startup
        # so these parameters are 0 by default
        super().__init__(system, IRR, duration, depreciation, income_tax,
                         operating_days, lang_factor, construction_schedule,
                         startup_months=0, startup_FOCfrac=0, startup_VOCfrac=0,
                         startup_salesfrac=0, finance_interest=0, finance_years=0,
                         finance_fraction=0, WC_over_FCI=WC_over_FCI)
        self.labor_cost = labor_cost
        self.fringe_benefits = fringe_benefits
        self.property_tax = property_tax
        self.property_insurance = property_insurance
        self.supplies= supplies
        self.maintenance = maintenance
        self.administration = administration


    def _DPI(self, installed_equipment_cost):
        return installed_equipment_cost


    def _TDC(self, DPI):
        return DPI

 
    def _FCI(self, TDC):
        return TDC

    
    def _FOC(self, FCI):
            return (FCI*(self.property_tax + self.property_insurance
                     + self.maintenance + self.administration)
                + self.labor_cost*(1+self.fringe_benefits+self.supplies))



tea =  Preprocessing(system=Preprocessing,
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
                  property_insurance=0.11,
                  supplies=0.20,
                  maintenance=0.01,
                  administration=0.005)

tea.show()
tea.get_cashflow_table()

In addition, the biosteam version is as follows:

Name: biosteam
Version: 2.36.4

The python version is as follows:
Python 3.10.9
IPython 8.12.0

The full traceback of error is as follows:

Traceback (most recent call last):

  File ~/opt/anaconda3/envs/env2/lib/python3.10/site-packages/spyder_kernels/py3compat.py:356 in compat_exec
    exec(code, globals, locals)

  File ~/Library/CloudStorage/[email protected]/My Drive/untitled1.py:269
    tea.show()

  File ~/opt/anaconda3/envs/env2/lib/python3.10/site-packages/biosteam/_tea.py:987 in show
    print(self._info())

  File ~/opt/anaconda3/envs/env2/lib/python3.10/site-packages/biosteam/_tea.py:983 in _info
    f' NPV: {self.NPV:,.0f} USD at {self.IRR:.1%} IRR')

  File ~/opt/anaconda3/envs/env2/lib/python3.10/site-packages/biosteam/_tea.py:777 in NPV
    taxable_cashflow, nontaxable_cashflow, depreciation = self._taxable_nontaxable_depreciation_cashflows()

  File ~/opt/anaconda3/envs/env2/lib/python3.10/site-packages/biosteam/_tea.py:804 in _taxable_nontaxable_depreciation_cashflows
    FOC = self._FOC(FCI)

  File ~/Library/CloudStorage/[email protected]/My Drive/untitled1.py:246 in _FOC
    return (FCI*(self.property_tax + self.property_insurance

TypeError: unsupported operand type(s) for *: 'property' and 'float'

Unable to load extension: pydevd_plugins.extensions.types.pydevd_plugin_pandas_types

yoelcortes · 2023-10-03T18:07:51Z

yoelcortes
Oct 3, 2023
Maintainer

I'm still unable to reproduce your error. I get the following error running your code:

ImportError: cannot import name 'TruckScale' from 'biosteam.units'...

Are you running any other scripts before running this one?

Although the error may be unrelated, you may find it helpful to upgrade to BioSTEAM 2.38.8. The documentation always shows the latest version.

1 reply

Uzair1shah Oct 3, 2023
Author

Please refer to this code instead, I have removed the truck scale unit operation. I thought Truck Scale was already in the biosteam unit operations repository, and you might be able to import it. However, it seems like I renamed/created that unit operations. The following code should work for you.

import biosteam as bst
import thermosteam as tmo
import graphviz
from biosteam import main_flowsheet as F, units, Stream, settings
import math as math

from biosteam.units import TruckScale, TruckDumper, DumpHopper, VibratingScreen, ScreenerFeeder, StackerReclaimer, ConveyingBelt, HammerMill, VibratingScreen, Mixer, Splitter

import functools

species = {}		#create an empty collection to hold your chemical species
  
#Biomass cellulose
species['cellulose'] = bst.Chemical('cellulose', phase='s', search_db=False, Cp=1.364, rho=1540, formula="C6H10O5", Hf=-1019000, MW =162.1424)            			#Cellulose, verified
species['cellA'] = bst.Chemical('cellA', search_db=False, default=True, phase='s', formula ="C6H10O5", Hf=-1019000, MW=162.1424)                                   	#Activated cellulose
  
# Biomass hemicellulose
species['gmsw'] = bst.Chemical('gmsw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['xyhw'] = bst.Chemical('xyhw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['hce1'] = bst.Chemical('hce1', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified
species['hce2'] = bst.Chemical('hce2', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified

# Biomass lignin
species['ligc'] = bst.Chemical('ligc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                    					#Verified
species['ligh'] = bst.Chemical('ligh', search_db=False, phase='s', formula="C22H28O9",MW=436.45892, Hf=-1722700)                                   					#Verified
species['ligO'] = bst.Chemical('ligO', search_db=False, phase='s', formula="C20H22O10",MW=422.38868, Hf=-1847500)                                  					#Verified
species['lig'] = bst.Chemical('lig', search_db=False, phase='s', formula="C11H12O4",MW=208.21388, Hf=-729310)                                     					#Secondary lignin intermediate, verified
species['ligcc'] = bst.Chemical('ligcc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                   				#Verified
species['ligoh'] = bst.Chemical('ligoh', search_db=False, phase='s', formula="C19H22O8",MW=378.37888, Hf=-1429200)                                  				#Verified

# Biomass extractives
species['tannin'] = bst.Chemical('tannin', search_db=False, phase='s', formula="C15H12O7", MW=304.25608,Hf=-1079700)                                 				#Tannic acid formula, verified
species['itann'] = bst.Chemical('itann', search_db=False, phase='s', formula="C8H4O4",MW=164.11736, Hf=-616980)                                     				#Intermediate tannin, verified
species['triglycerides'] = bst.Chemical('triglycerides', search_db=False, phase='s', formula="C76H52O46", MW=897.4168, Hf=-1546300)                          		#Triglycerides formula, update the triglycerides Hf

# Biomass inorganics
# This could be defined more explictly later
species['ash'] = bst.Chemical('CaO', phase='s', formula="CaO", Hf=5382027.245)                                         												#Ash, verify HF

# Pyrolysis solid products

species['char'] = bst.Chemical('carbon', phase='s', formula="C")        																							#Carbon, verified, properties from refs 1,2
species['hmwl'] = bst.Chemical('hmwl', search_db=False, phase='l', formula="C24H28O4", MW=380.48392, Hf=-958260)                                   					#Verified, High molecular weight lignin

# Fluidized bed media
species['sand'] = bst.Chemical('SiO2', formula="SiO2", Hf=-910860)                                            														#Sand, inert, properties from ref2

#metaplastic species
species['gco2'] = bst.Chemical('gco2',search_db=False, phase='s', formula="CO2", MW=44.01, Hf=-393510)                                            					#Trapped carbon dioxide, properties from ref2
species['gco'] = bst.Chemical('gco',search_db=False, phase='s', formula="CO",MW=28.01, Hf=-110520)                                               					#Trapped carbon monoxide, properties from ref2
species['gmeoh'] = bst.Chemical('gmeoh',search_db=False, phase='s', formula="CH3OH", MW=32.042, Hf=-238640)                                         				#Trapped methanol, verified, properties from ref2
species['gcoh2_l'] =bst.Chemical('gcoh2_l', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)                                        				#Trapped formaldehyde, properties from ref2, loose
species['gc2h6'] = bst.Chemical('gc2h6', search_db=False, phase='s', formula="C2H6", MW=30.07, Hf=-84667)                                           				#Trapped ethane, properties from ref2
species['gch4'] = bst.Chemical('gch4', search_db=False, phase='s', formula="CH4",MW=16.041, Hf=-74840)                                             					#Trapped methane, properties from ref2
species['gc2h4'] = bst.Chemical('gc2h4', search_db=False, phase='s', formula="C2H4", MW=28.052, Hf=52283)                                            				#Trapped ethylene, properties from ref2
species['gc6h5oh'] = bst.Chemical('gc6h5oh', search_db=False, phase='s', formula="C6H6O", MW=94.113, Hf=-96400)                                       				#Trapped phenol, properties from ref3
species['gcoh2_s'] = bst.Chemical('gcoh2_s', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)	                                      				#Trapped formaldehyde, properties from ref2, stiff
species['gh2'] = bst.Chemical('gh2', search_db=False, phase='s', formula="H2",MW=2.016, Hf=0)		                                              					#Trapped hydrogen, properties from ref2

#non-biomass species
species['argon'] = bst.Chemical('argon', formula="Ar")  
species['nitrogen'] = bst.Chemical('nitrogen', formula="N2")  
species['oxygen'] = bst.Chemical('oxygen', formula="O2")  

#Pyrolysis products
## Gas phase species
species['ethylene'] = bst.Chemical('ethylene', formula="C2H4", search_ID='ethylene', phase='g')                  
species['methane'] = bst.Chemical('methane', formula="CH4", search_ID='methane',  phase='g')                                   										#ref2
species['hydrogen'] = bst.Chemical('hydrogen', formula="H2", search_ID='hydrogen',  phase='g')                                     									#ref2
species['cacrbonmonoxide'] = bst.Chemical('carbon monoxide', formula="CO", search_ID='carbon monoxide', phase='g')                               					#ref2
species['carbondioxide'] = bst.Chemical('carbon dioxide', formula="CO2", search_ID='carbon dioxide', phase='g')                            							#ref2
species['formaldehyde'] = bst.Chemical('formaldehyde', formula="CH2O", search_ID='formaldehyde', phase='g')                            								#ref2

## Liquid phase species
species['acrolein'] = bst.Chemical('acrolein', formula="C3H4O", search_ID='acrolein', phase='l')      										#https://atct.anl.gov/Thermochemical%20Data/version%201.122e/species/?species_number=1209

species['sinapylaldehyde'] = bst.Chemical('sinapyaldehyde', search_db=False, phase='l', formula="C11H12O4", MW=208.21388, Hf=-483800)
species['methyllinoleate'] = bst.Chemical('methyllinoleate', phase='l', search_db=False, formula="C19H34O2", MW=294.4721, Hf=-604880)
species['glyoxal'] = bst.Chemical('glyoxal', search_ID='glyoxal', phase='l')							    														#NIST webbook
species['acetaldehyde'] = bst.Chemical('acetaldehyde', search_ID='acetaldehyde', phase='l')                       													#Acetaldehyde, properties from ref2 
#species['aceticacid'] = bst.Chemical('aceticacid', search_db=False, default=True, formula="CH3COOH", MW=60.05, Hf=-484930, Hvap=0, Psat=0, Tb=0)                   #Acetic acid, properties from ref2 
species['glycolaldehyde'] = bst.Chemical('glycolaldehyde', search_ID='glycolaldehyde', phase='l')                                    								#Source Pyrolysis of the Simplest Carbohydrate, Glycolaldehyde (CHO-CH2OH), and Glyoxal in a Heated Microreactor. Find additional sources
species['ethanol'] = bst.Chemical('ethanol', search_ID='ethanol', phase='l')                             															#Ethanol, properties from ref2 
species['propionaldehyde'] = bst.Chemical('propionaldehyde', search_ID='propionaldehyde', phase='l')                    							#Source: NIST webbook https://webbook.nist.gov/cgi/cbook.cgi?ID=C123386&Mask=2 
species['hydroxypropanal'] = bst.Chemical('hydroxypropanal', search_ID='hydroxypropanal',phase='l')  												#https://www.chemeo.com/cid/20-360-9/2%2C2-Dimethyl-3-hydroxypropionaldehyde
species['furfural'] = bst.Chemical('furfural', search_ID='furfural', phase='l')
species['xylosan'] = bst.Chemical('xylosan', phase='l', search_db=False, formula="C5H8O4", MW=132.115, Hf=-647000)                             				#University of Pittsburg Pitt Quantum Repository, DOI 10.17614/Q4SQ8QT7Q
species['levoglucosan'] = bst.Chemical('levoglucosan', search_ID='levoglucosan', phase='l')                      											#Gas phase Hf value, ref4.
species['phenol'] = bst.Chemical('phenol', search_ID='phenol', phase='l')                                 													#phenol, properties from ref5
species['hydroxymethylfurfural'] = bst.Chemical('hydroxymethylfurfural', search_db=False, formula="C6H6O3", MW=126.1100, Hf=-437400, phase='l')           	#Hydromethyl furfural DOI:10.1021/ie901012g
species['anisole'] = bst.Chemical('anisole', search_ID='anisole', phase='l')                                  												#DOI: 10.1021/jp507267f
species['pcoumarylalcohol'] = bst.Chemical('pcoumarylalcohol', search_db=False, formula="C9H10O2", MW=150.1772, Hf=-193500, phase='l')                 		#Unverified
species['formicacid'] = bst.Chemical('formicacid', search_db=False, formula="CH2O2", MW=46.0254, Hf=-431650, phase='l')                           			#ref: Heat of formation of liquid. Should it be of gas?
species['methanol'] = bst.Chemical('methanol', search_ID='methanol', phase='l')                               												#ref2
species['water'] = bst.Chemical('water', search_ID='water', phase='l')                                   													#ref2

species['propaniocacid'] = bst.Chemical('propaniocacid', search_db=False, formula="C3H6O2", Hf=-199481.24962000002, phase='l')
species['linalylpropioinate'] = bst.Chemical('linalylpropioinate', search_db=False, formula="C19H34O2", Hf=-194296.922694, phase='l')
species['vanillin'] = bst.Chemical('vanillin', search_ID='vanillin', phase='l')
species['cresol'] = bst.Chemical('cresol', search_db=False, formula="C7H8O", Hf=-34212.624394, phase='l')

# Create an empty Chemicals object and add each species to it
chemicals = tmo.Chemicals([])

for key, value in species.items():
    #print(key)
    value.default() # call the default method to populate all other properties not defined explicitly. Todo: some properties default to water's properties. Check individual chemical properties for accuracy once model is developed. 
    chemicals.append(value) # push a chemical object to back of the chemicals collection object

chemicals.compile()								# compile the loaded chemicals. Must do.
chemicals.set_synonym('water', 'moisture')		# set synonyms or common name for select species
chemicals.set_synonym('CaO', 'ash')

settings.set_thermo(chemicals, cache=False)				# Todo: read more set_thermo method.


# set raw material price (USD/kg)
price = {'biomass': 0.08, # @10% moisture content wet basis
         'water': 0.000353,
         'steam': 0.017,
         'waste': -0.33,
         'gasoline': 0.756} # 2 USD/gal

#setting electructy price
settings.electricity_price = 0.1           # set electricity price to your desired rate

# set the name of your flowsheet
bst.main_flowsheet.set_flowsheet('fast_pyrolysis_biorefinery')

# define a chemical group with its composition in mass fraction here
chemicals.define_group('biomass',IDs = ['cellulose','xyhw','ligc','ligh','ligO','triglycerides','tannin','moisture','ash'], composition = [0.2347,0.1857,0.0915,0.0761,0.1290,0.0386,0.0176,0.10,0.0915], wt=True)

# Truckscale UOP
# S101 - Biomass inlet into the truck scale
# S102 - Biomass outlet out of the truck scale
S101 = Stream('S101', biomass=333334.2, units='kg/hr', phase='s', price=price['biomass'])



# Vibrating Screen
S105 = Stream('S105', phase='s')
S106 = Stream('S106', phase='s')

UO104 = VibratingScreen('UO104', ins=S101, outs=(S105,S106), split=0.9)                     # Vibrating Screen


# Simulation of UO103 (Vibrating Screen)
UO104.simulate()
UO104.show(flow='kg/hr', composition = False, N=9)
r4 = UO104.results()
print(r4)
                                                                                           # Rotary Dryer
S107 = Stream('S107', phase='g') #Air
S108 = Stream('S108', phase='g') #Natural gas
S109 = Stream('S109', phase='s') #dry feedstock/solids
S110 = Stream('S110', phase='g') #hot air
S111 = Stream('S111', phase='g') #emissions
UO105 = units.DrumDryer('UO105', ins=(S105,S107,S108), outs=(S109,S110,S111),moisture_content = 0.01, split = (dict(moisture=0.9)) )
# Simulation of UO105 (Rotary dryer)
UO105.simulate()
UO105.show(flow='kg/hr', composition = False, N=10)
r5 = UO105.results()
print(r5)

S113 = Stream('S113', phase='s') #recycle stream
S114 = Stream('S114', phase='s')
UO106 = units.Mixer('UO106', ins=(S109, S113), outs=(S114))

# Hammermill
S112 = Stream('S112', phase='s') #crushed dry feedstock
S115 = Stream('S115', phase='s') 

UO107 = units.HammerMill('UO107', ins=(S114), outs=(S115))

#Vibrating Screen
S116 = Stream('S116', phase='s')
UO108 = VibratingScreen('UO108', ins=S115, outs=(S113,S116), split=0.1) 

Preprocessing = bst.main_flowsheet.create_system('Preprocessing')
Preprocessing.simulate()
Preprocessing.diagram('cluster', format='png')

UO106.show(flow='kg/hr', composition = False, N=9)
r6 = UO106.results()
print(r6)

UO107.show(flow='kg/hr', composition = False, N=9)
r7 = UO107.results()
print(r7)

UO108.show(flow='kg/hr', composition = False, N=9)
r8 = UO108.results()
print(r8)

class Preprocessing(bst.TEA):

    def __init__(self, system, IRR, duration, depreciation, income_tax,
                 operating_days, lang_factor, construction_schedule, WC_over_FCI,
                 labor_cost, fringe_benefits, property_tax,
                 property_insurance, supplies, maintenance, administration):
        # Huang et. al. does not take into account financing or startup
        # so these parameters are 0 by default
        super().__init__(system, IRR, duration, depreciation, income_tax,
                         operating_days, lang_factor, construction_schedule,
                         startup_months=0, startup_FOCfrac=0, startup_VOCfrac=0,
                         startup_salesfrac=0, finance_interest=0, finance_years=0,
                         finance_fraction=0, WC_over_FCI=WC_over_FCI)
        self.labor_cost = labor_cost
        self.fringe_benefits = fringe_benefits
        self.property_tax = property_tax
        self.property_insurance = property_insurance
        self.supplies= supplies
        self.maintenance = maintenance
        self.administration = administration


    def _DPI(self, installed_equipment_cost):
        return installed_equipment_cost


    def _TDC(self, DPI):
        return DPI

 
    def _FCI(self, TDC):
        return TDC

    
    def _FOC(self, FCI):
            return (FCI*(self.property_tax + self.property_insurance
                     + self.maintenance + self.administration)
                + self.labor_cost*(1+self.fringe_benefits+self.supplies))



tea =  Preprocessing(system=Preprocessing,
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
                  property_insurance=0.11,
                  supplies=0.20,
                  maintenance=0.01,
                  administration=0.005)

tea.show()
tea.get_cashflow_table()

Also, I will try updating the biosteam version and see if that helps with the error.

yoelcortes · 2023-10-03T19:21:02Z

yoelcortes
Oct 3, 2023
Maintainer

@Uzair1shah, there are many units listed in your code that are not in BioSTEAM. These units were never in BioSTEAM previously (nor have there been any renaming of unit operations in the last year). When you create new unit subclasses, BioSTEAM automatically registers them in the biosteam.units module, so the reason why your code might work on your side is because those units are being created in another script.

I get the same error with the most recent code you posted:

ImportError: cannot import name 'TruckScale' from 'biosteam.units'

It may be better if you try to minimize your code such that it still reproduces the error you are seeing.

1 reply

Uzair1shah Oct 3, 2023
Author

I have removed the unit operations that might be on my side from the import list:

from biosteam.units import VibratingScreen, HammerMill, Mixer, Splitter

These unit operations were already in the biosteam so it should work for you. If it still does not work for you, I can attach scripts for these unit operations as well to reproduce the same error. However, the main code is as follows:

import biosteam as bst
import thermosteam as tmo
import graphviz
from biosteam import main_flowsheet as F, units, Stream, settings
import math as math

from biosteam.units import VibratingScreen, HammerMill, Mixer, Splitter

import functools

species = {}		#create an empty collection to hold your chemical species
  
#Biomass cellulose
species['cellulose'] = bst.Chemical('cellulose', phase='s', search_db=False, Cp=1.364, rho=1540, formula="C6H10O5", Hf=-1019000, MW =162.1424)            			#Cellulose, verified
species['cellA'] = bst.Chemical('cellA', search_db=False, default=True, phase='s', formula ="C6H10O5", Hf=-1019000, MW=162.1424)                                   	#Activated cellulose
  
# Biomass hemicellulose
species['gmsw'] = bst.Chemical('gmsw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['xyhw'] = bst.Chemical('xyhw', search_db=False, default=True, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      	#Verified
species['hce1'] = bst.Chemical('hce1', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified
species['hce2'] = bst.Chemical('hce2', search_db=False, phase='s', formula="C5H8O4",MW = 132.11612, Hf=-759200)                                      				#Verified

# Biomass lignin
species['ligc'] = bst.Chemical('ligc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                    					#Verified
species['ligh'] = bst.Chemical('ligh', search_db=False, phase='s', formula="C22H28O9",MW=436.45892, Hf=-1722700)                                   					#Verified
species['ligO'] = bst.Chemical('ligO', search_db=False, phase='s', formula="C20H22O10",MW=422.38868, Hf=-1847500)                                  					#Verified
species['lig'] = bst.Chemical('lig', search_db=False, phase='s', formula="C11H12O4",MW=208.21388, Hf=-729310)                                     					#Secondary lignin intermediate, verified
species['ligcc'] = bst.Chemical('ligcc', search_db=False, phase='s', formula="C15H14O4",MW=258.27376, Hf=-759390)                                   				#Verified
species['ligoh'] = bst.Chemical('ligoh', search_db=False, phase='s', formula="C19H22O8",MW=378.37888, Hf=-1429200)                                  				#Verified

# Biomass extractives
species['tannin'] = bst.Chemical('tannin', search_db=False, phase='s', formula="C15H12O7", MW=304.25608,Hf=-1079700)                                 				#Tannic acid formula, verified
species['itann'] = bst.Chemical('itann', search_db=False, phase='s', formula="C8H4O4",MW=164.11736, Hf=-616980)                                     				#Intermediate tannin, verified
species['triglycerides'] = bst.Chemical('triglycerides', search_db=False, phase='s', formula="C76H52O46", MW=897.4168, Hf=-1546300)                          		#Triglycerides formula, update the triglycerides Hf

# Biomass inorganics
# This could be defined more explictly later
species['ash'] = bst.Chemical('CaO', phase='s', formula="CaO", Hf=5382027.245)                                         												#Ash, verify HF

# Pyrolysis solid products

species['char'] = bst.Chemical('carbon', phase='s', formula="C")        																							#Carbon, verified, properties from refs 1,2
species['hmwl'] = bst.Chemical('hmwl', search_db=False, phase='l', formula="C24H28O4", MW=380.48392, Hf=-958260)                                   					#Verified, High molecular weight lignin

# Fluidized bed media
species['sand'] = bst.Chemical('SiO2', formula="SiO2", Hf=-910860)                                            														#Sand, inert, properties from ref2

#metaplastic species
species['gco2'] = bst.Chemical('gco2',search_db=False, phase='s', formula="CO2", MW=44.01, Hf=-393510)                                            					#Trapped carbon dioxide, properties from ref2
species['gco'] = bst.Chemical('gco',search_db=False, phase='s', formula="CO",MW=28.01, Hf=-110520)                                               					#Trapped carbon monoxide, properties from ref2
species['gmeoh'] = bst.Chemical('gmeoh',search_db=False, phase='s', formula="CH3OH", MW=32.042, Hf=-238640)                                         				#Trapped methanol, verified, properties from ref2
species['gcoh2_l'] =bst.Chemical('gcoh2_l', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)                                        				#Trapped formaldehyde, properties from ref2, loose
species['gc2h6'] = bst.Chemical('gc2h6', search_db=False, phase='s', formula="C2H6", MW=30.07, Hf=-84667)                                           				#Trapped ethane, properties from ref2
species['gch4'] = bst.Chemical('gch4', search_db=False, phase='s', formula="CH4",MW=16.041, Hf=-74840)                                             					#Trapped methane, properties from ref2
species['gc2h4'] = bst.Chemical('gc2h4', search_db=False, phase='s', formula="C2H4", MW=28.052, Hf=52283)                                            				#Trapped ethylene, properties from ref2
species['gc6h5oh'] = bst.Chemical('gc6h5oh', search_db=False, phase='s', formula="C6H6O", MW=94.113, Hf=-96400)                                       				#Trapped phenol, properties from ref3
species['gcoh2_s'] = bst.Chemical('gcoh2_s', search_db=False, phase='s', formula="CH2O",MW=30.026, Hf=-115890)	                                      				#Trapped formaldehyde, properties from ref2, stiff
species['gh2'] = bst.Chemical('gh2', search_db=False, phase='s', formula="H2",MW=2.016, Hf=0)		                                              					#Trapped hydrogen, properties from ref2

#non-biomass species
species['argon'] = bst.Chemical('argon', formula="Ar")  
species['nitrogen'] = bst.Chemical('nitrogen', formula="N2")  
species['oxygen'] = bst.Chemical('oxygen', formula="O2")  

#Pyrolysis products
## Gas phase species
species['ethylene'] = bst.Chemical('ethylene', formula="C2H4", search_ID='ethylene', phase='g')                  
species['methane'] = bst.Chemical('methane', formula="CH4", search_ID='methane',  phase='g')                                   										#ref2
species['hydrogen'] = bst.Chemical('hydrogen', formula="H2", search_ID='hydrogen',  phase='g')                                     									#ref2
species['cacrbonmonoxide'] = bst.Chemical('carbon monoxide', formula="CO", search_ID='carbon monoxide', phase='g')                               					#ref2
species['carbondioxide'] = bst.Chemical('carbon dioxide', formula="CO2", search_ID='carbon dioxide', phase='g')                            							#ref2
species['formaldehyde'] = bst.Chemical('formaldehyde', formula="CH2O", search_ID='formaldehyde', phase='g')                            								#ref2

## Liquid phase species
species['acrolein'] = bst.Chemical('acrolein', formula="C3H4O", search_ID='acrolein', phase='l')      										#https://atct.anl.gov/Thermochemical%20Data/version%201.122e/species/?species_number=1209

species['sinapylaldehyde'] = bst.Chemical('sinapyaldehyde', search_db=False, phase='l', formula="C11H12O4", MW=208.21388, Hf=-483800)
species['methyllinoleate'] = bst.Chemical('methyllinoleate', phase='l', search_db=False, formula="C19H34O2", MW=294.4721, Hf=-604880)
species['glyoxal'] = bst.Chemical('glyoxal', search_ID='glyoxal', phase='l')							    														#NIST webbook
species['acetaldehyde'] = bst.Chemical('acetaldehyde', search_ID='acetaldehyde', phase='l')                       													#Acetaldehyde, properties from ref2 
#species['aceticacid'] = bst.Chemical('aceticacid', search_db=False, default=True, formula="CH3COOH", MW=60.05, Hf=-484930, Hvap=0, Psat=0, Tb=0)                   #Acetic acid, properties from ref2 
species['glycolaldehyde'] = bst.Chemical('glycolaldehyde', search_ID='glycolaldehyde', phase='l')                                    								#Source Pyrolysis of the Simplest Carbohydrate, Glycolaldehyde (CHO-CH2OH), and Glyoxal in a Heated Microreactor. Find additional sources
species['ethanol'] = bst.Chemical('ethanol', search_ID='ethanol', phase='l')                             															#Ethanol, properties from ref2 
species['propionaldehyde'] = bst.Chemical('propionaldehyde', search_ID='propionaldehyde', phase='l')                    							#Source: NIST webbook https://webbook.nist.gov/cgi/cbook.cgi?ID=C123386&Mask=2 
species['hydroxypropanal'] = bst.Chemical('hydroxypropanal', search_ID='hydroxypropanal',phase='l')  												#https://www.chemeo.com/cid/20-360-9/2%2C2-Dimethyl-3-hydroxypropionaldehyde
species['furfural'] = bst.Chemical('furfural', search_ID='furfural', phase='l')
species['xylosan'] = bst.Chemical('xylosan', phase='l', search_db=False, formula="C5H8O4", MW=132.115, Hf=-647000)                             				#University of Pittsburg Pitt Quantum Repository, DOI 10.17614/Q4SQ8QT7Q
species['levoglucosan'] = bst.Chemical('levoglucosan', search_ID='levoglucosan', phase='l')                      											#Gas phase Hf value, ref4.
species['phenol'] = bst.Chemical('phenol', search_ID='phenol', phase='l')                                 													#phenol, properties from ref5
species['hydroxymethylfurfural'] = bst.Chemical('hydroxymethylfurfural', search_db=False, formula="C6H6O3", MW=126.1100, Hf=-437400, phase='l')           	#Hydromethyl furfural DOI:10.1021/ie901012g
species['anisole'] = bst.Chemical('anisole', search_ID='anisole', phase='l')                                  												#DOI: 10.1021/jp507267f
species['pcoumarylalcohol'] = bst.Chemical('pcoumarylalcohol', search_db=False, formula="C9H10O2", MW=150.1772, Hf=-193500, phase='l')                 		#Unverified
species['formicacid'] = bst.Chemical('formicacid', search_db=False, formula="CH2O2", MW=46.0254, Hf=-431650, phase='l')                           			#ref: Heat of formation of liquid. Should it be of gas?
species['methanol'] = bst.Chemical('methanol', search_ID='methanol', phase='l')                               												#ref2
species['water'] = bst.Chemical('water', search_ID='water', phase='l')                                   													#ref2

species['propaniocacid'] = bst.Chemical('propaniocacid', search_db=False, formula="C3H6O2", Hf=-199481.24962000002, phase='l')
species['linalylpropioinate'] = bst.Chemical('linalylpropioinate', search_db=False, formula="C19H34O2", Hf=-194296.922694, phase='l')
species['vanillin'] = bst.Chemical('vanillin', search_ID='vanillin', phase='l')
species['cresol'] = bst.Chemical('cresol', search_db=False, formula="C7H8O", Hf=-34212.624394, phase='l')

# Create an empty Chemicals object and add each species to it
chemicals = tmo.Chemicals([])

for key, value in species.items():
    #print(key)
    value.default() # call the default method to populate all other properties not defined explicitly. Todo: some properties default to water's properties. Check individual chemical properties for accuracy once model is developed. 
    chemicals.append(value) # push a chemical object to back of the chemicals collection object

chemicals.compile()								# compile the loaded chemicals. Must do.
chemicals.set_synonym('water', 'moisture')		# set synonyms or common name for select species
chemicals.set_synonym('CaO', 'ash')

settings.set_thermo(chemicals, cache=False)				# Todo: read more set_thermo method.


# set raw material price (USD/kg)
price = {'biomass': 0.08, # @10% moisture content wet basis
         'water': 0.000353,
         'steam': 0.017,
         'waste': -0.33,
         'gasoline': 0.756} # 2 USD/gal

#setting electructy price
settings.electricity_price = 0.1           # set electricity price to your desired rate

# set the name of your flowsheet
bst.main_flowsheet.set_flowsheet('fast_pyrolysis_biorefinery')

# define a chemical group with its composition in mass fraction here
chemicals.define_group('biomass',IDs = ['cellulose','xyhw','ligc','ligh','ligO','triglycerides','tannin','moisture','ash'], composition = [0.2347,0.1857,0.0915,0.0761,0.1290,0.0386,0.0176,0.10,0.0915], wt=True)


S101 = Stream('S101', biomass=333334.2, units='kg/hr', phase='s', price=price['biomass'])



# Vibrating Screen
S105 = Stream('S105', phase='s')
S106 = Stream('S106', phase='s')

UO104 = VibratingScreen('UO104', ins=S101, outs=(S105,S106), split=0.9)                     # Vibrating Screen


# Simulation of UO103 (Vibrating Screen)
UO104.simulate()
UO104.show(flow='kg/hr', composition = False, N=9)
r4 = UO104.results()
print(r4)
                                                                                           # Rotary Dryer
S107 = Stream('S107', phase='g') #Air
S108 = Stream('S108', phase='g') #Natural gas
S109 = Stream('S109', phase='s') #dry feedstock/solids
S110 = Stream('S110', phase='g') #hot air
S111 = Stream('S111', phase='g') #emissions
UO105 = units.DrumDryer('UO105', ins=(S105,S107,S108), outs=(S109,S110,S111),moisture_content = 0.01, split = (dict(moisture=0.9)) )
# Simulation of UO105 (Rotary dryer)
UO105.simulate()
UO105.show(flow='kg/hr', composition = False, N=10)
r5 = UO105.results()
print(r5)

S113 = Stream('S113', phase='s') #recycle stream
S114 = Stream('S114', phase='s')
UO106 = units.Mixer('UO106', ins=(S109, S113), outs=(S114))

# Hammermill
S112 = Stream('S112', phase='s') #crushed dry feedstock
S115 = Stream('S115', phase='s') 

UO107 = units.HammerMill('UO107', ins=(S114), outs=(S115))

#Vibrating Screen
S116 = Stream('S116', phase='s')
UO108 = VibratingScreen('UO108', ins=S115, outs=(S113,S116), split=0.1) 

Preprocessing = bst.main_flowsheet.create_system('Preprocessing')
Preprocessing.simulate()
Preprocessing.diagram('cluster', format='png')

UO106.show(flow='kg/hr', composition = False, N=9)
r6 = UO106.results()
print(r6)

UO107.show(flow='kg/hr', composition = False, N=9)
r7 = UO107.results()
print(r7)

UO108.show(flow='kg/hr', composition = False, N=9)
r8 = UO108.results()
print(r8)

class Preprocessing(bst.TEA):

    def __init__(self, system, IRR, duration, depreciation, income_tax,
                 operating_days, lang_factor, construction_schedule, WC_over_FCI,
                 labor_cost, fringe_benefits, property_tax,
                 property_insurance, supplies, maintenance, administration):
        # Huang et. al. does not take into account financing or startup
        # so these parameters are 0 by default
        super().__init__(system, IRR, duration, depreciation, income_tax,
                         operating_days, lang_factor, construction_schedule,
                         startup_months=0, startup_FOCfrac=0, startup_VOCfrac=0,
                         startup_salesfrac=0, finance_interest=0, finance_years=0,
                         finance_fraction=0, WC_over_FCI=WC_over_FCI)
        self.labor_cost = labor_cost
        self.fringe_benefits = fringe_benefits
        self.property_tax = property_tax
        self.property_insurance = property_insurance
        self.supplies= supplies
        self.maintenance = maintenance
        self.administration = administration


    def _DPI(self, installed_equipment_cost):
        return installed_equipment_cost


    def _TDC(self, DPI):
        return DPI

 
    def _FCI(self, TDC):
        return TDC

    
    def _FOC(self, FCI):
            return (FCI*(self.property_tax + self.property_insurance
                     + self.maintenance + self.administration)
                + self.labor_cost*(1+self.fringe_benefits+self.supplies))



tea =  Preprocessing(system=Preprocessing,
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
                  property_insurance=0.11,
                  supplies=0.20,
                  maintenance=0.01,
                  administration=0.005)

tea.show()
tea.get_cashflow_table()

yoelcortes · 2023-10-03T23:27:49Z

yoelcortes
Oct 3, 2023
Maintainer

@Uzair1shah, OK! I was able to reproduce the error. The issue is the first argument to a TEA object should be a System object not a TEA subclass:

sys = bst.main_flowsheet.create_system('sys')
tea =  Preprocessing(system=sys, # Not Preprocessing
                  IRR=0.15,
                  duration=(2018, 2038),
                  depreciation='MACRS7',
                  income_tax=0.21, # Previously 35% in published study
                  operating_days=200,
                  lang_factor=3,
                  construction_schedule=(0.4, 0.6),
                  WC_over_FCI=0.05,
                  labor_cost=2.5e6,
                  fringe_benefits=0.4,
                  property_tax=0.11,
                  property_insurance=0.11,
                  supplies=0.20,
                  maintenance=0.01,
                  administration=0.005)

Thanks!

1 reply

Uzair1shah Oct 4, 2023
Author

That works. Thanks a lot.

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

TEA not working in BioSTEAM #174

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Replies: 4 comments 4 replies

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

Select a reply

TEA not working in BioSTEAM #174

Uzair1shah Oct 2, 2023

Replies: 4 comments · 4 replies

yoelcortes Oct 3, 2023 Maintainer

Uzair1shah Oct 3, 2023 Author

yoelcortes Oct 3, 2023 Maintainer

Uzair1shah Oct 3, 2023 Author

yoelcortes Oct 3, 2023 Maintainer

Uzair1shah Oct 3, 2023 Author

yoelcortes Oct 3, 2023 Maintainer

Uzair1shah Oct 4, 2023 Author

Uzair1shah
Oct 2, 2023

Replies: 4 comments 4 replies

yoelcortes
Oct 3, 2023
Maintainer

Uzair1shah Oct 3, 2023
Author

yoelcortes
Oct 3, 2023
Maintainer

Uzair1shah Oct 3, 2023
Author

yoelcortes
Oct 3, 2023
Maintainer

Uzair1shah Oct 3, 2023
Author

yoelcortes
Oct 3, 2023
Maintainer

Uzair1shah Oct 4, 2023
Author