Skip to content

Commit

Permalink
working bsm2p settings
Browse files Browse the repository at this point in the history
  • Loading branch information
@joyxyz1994
joyxyz1994 committed Sep 27, 2024
1 parent 8631e2b commit cde6715
Show file tree
Hide file tree
Showing 3 changed files with 78 additions and 62 deletions.
5 changes: 3 additions & 2 deletions exposan/adm/test_adm1p_interfaces.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -268,5 +268,6 @@

#%%
def run():
sys.simulate(state_reset_hook='reset_cache', t_span=(0, 200), method='BDF')
# run()
sys.simulate(state_reset_hook='reset_cache', t_span=(0, 100), method='BDF')

run()
Binary file modified exposan/bsm2/data/bsm2p_init.xlsx
View file
Binary file not shown.
135 changes: 75 additions & 60 deletions exposan/bsm2/system_P.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,8 @@ def create_system(flowsheet=None, default_init_conds=True):
electron_acceptor_dependent_decay=True,
k_h=2.46, mu_H=4.23, q_fe=2.11, b_H=0.28, mu_PAO=0.82,
q_PP=1.23, q_PHA=2.46, b_PAO=0.14, b_PP=0.14, b_PHA=0.14,
mu_AUT=0.61, b_AUT=0.09)
mu_AUT=0.61, b_AUT=0.09
)
thermo_asm = qs.get_thermo()

# Influent
Expand All @@ -83,56 +84,64 @@ def create_system(flowsheet=None, default_init_conds=True):
)

DO_ID = 'S_O2'
aer1 = aer2 = pc.DiffusedAeration('aer1', DO_ID, KLa=120, DOsat=SOSAT1, V=V_ae)
aer3 = pc.DiffusedAeration('aer3', DO_ID, KLa=60, DOsat=SOSAT1, V=V_ae)
gstrip = True
anae_kwargs = dict(V_max=V_anae, aeration=None, suspended_growth_model=asm, gas_stripping=gstrip)
anox_kwargs = dict(V_max=V_anox, aeration=None, suspended_growth_model=asm, gas_stripping=gstrip)
ae_kwargs = dict(V_max=V_ae, DO_ID=DO_ID, suspended_growth_model=asm, gas_stripping=gstrip)

# aer1 = aer2 = pc.DiffusedAeration('aer1', DO_ID, KLa=120, DOsat=SOSAT1, V=V_ae)
# aer3 = pc.DiffusedAeration('aer3', DO_ID, KLa=60, DOsat=SOSAT1, V=V_ae)
# anae_kwargs = dict(V_max=V_anae, aeration=None, suspended_growth_model=asm, gas_stripping=gstrip)
# anox_kwargs = dict(V_max=V_anox, aeration=None, suspended_growth_model=asm, gas_stripping=gstrip)
# ae_kwargs = dict(V_max=V_ae, DO_ID=DO_ID, suspended_growth_model=asm, gas_stripping=gstrip)
c2_kwargs = dict(
underflow=Q_ras, wastage=Q_was,
surface_area=1500, height=4, N_layer=10, feed_layer=5,
X_threshold=3000, v_max=474, v_max_practical=250,
rh=5.76e-4, rp=2.86e-3, fns=2.28e-3
)

A1 = su.CSTR('A1', ins=[C1-0, 'RAS', carb], **anae_kwargs)
# A1 = su.CSTR('A1', ins=[C1-0, 'RAS'], **anae_kwargs)
A2 = su.CSTR('A2', A1-0, **anae_kwargs)
# A3 = su.CSTR('A3', [A2-0, 'RWW', carb], **anox_kwargs)
A3 = su.CSTR('A3', [A2-0, 'RWW'], **anox_kwargs)
A4 = su.CSTR('A4', A3-0, **anox_kwargs)
O1 = su.CSTR('O1', A4-0, aeration=aer1, **ae_kwargs)
O2 = su.CSTR('O2', O1-0, aeration=aer2, **ae_kwargs)
O3 = su.CSTR('O3', O2-0, [1-A3, 'treated'], split=[Q_intr, Q+Q_ras],
aeration=aer3, **ae_kwargs)
C2 = su.FlatBottomCircularClarifier(
'C2', O3-1, ['effluent', 1-A1, 'WAS'],
**c2_kwargs
rh=5.76e-4, rp=2.86e-3, fns=2.28e-3,
# maximum_nonsettleable_solids=8.0,
)

# # AS = su.PFR('AS', ins=[C1-0, 'RAS', carb], outs='treated',
# AS = su.PFR('AS', ins=[C1-0, 'RAS'], outs='treated',
# N_tanks_in_series=7,
# V_tanks=[V_anae]*2+[V_anox]*2+[V_ae]*3,
# # influent_fractions=[[1]+[0]*6]*3,
# influent_fractions=[[1]+[0]*6]*2,
# internal_recycles=[(6,2,Q_intr)],
# kLa=[0]*4+[120,120,60], DO_ID=DO_ID, DO_sat=SOSAT1,
# suspended_growth_model=asm,
# gas_stripping=True)

# A1 = su.CSTR('A1', ins=[C1-0, 'RAS', carb], **anae_kwargs)
# # A1 = su.CSTR('A1', ins=[C1-0, 'RAS'], **anae_kwargs)
# A2 = su.CSTR('A2', A1-0, **anae_kwargs)
# # A3 = su.CSTR('A3', [A2-0, 'RWW', carb], **anox_kwargs)
# A3 = su.CSTR('A3', [A2-0, 'RWW'], **anox_kwargs)
# A4 = su.CSTR('A4', A3-0, **anox_kwargs)
# O1 = su.CSTR('O1', A4-0, aeration=aer1, **ae_kwargs)
# O2 = su.CSTR('O2', O1-0, aeration=aer2, **ae_kwargs)
# O3 = su.CSTR('O3', O2-0, [1-A3, 'treated'], split=[Q_intr, Q+Q_ras],
# aeration=aer3, **ae_kwargs)
# C2 = su.FlatBottomCircularClarifier(
# 'C2', AS-0, ['effluent', 1-AS, 'WAS'],
# 'C2', O3-1, ['effluent', 1-A1, 'WAS'],
# **c2_kwargs
# )

TC1 = su.Thickener('TC1', C2-2, outs=['thickened_sludge', ''],
thickening_perc=7, TSS_removal_perc=98)
M1 = su.Mixer('M1', ins=(C1-1, TC1-0))
AS = su.PFR('AS', ins=[C1-0, 'RAS', carb], outs='treated',
# AS = su.PFR('AS', ins=[C1-0, 'RAS'], outs='treated',
N_tanks_in_series=7,
V_tanks=[V_anae]*2+[V_anox]*2+[V_ae]*3,
influent_fractions=[[1]+[0]*6]*2 + [[0,0,1,0,0,0,0]],
# influent_fractions=[[1]+[0]*6]*2,
internal_recycles=[(6,2,Q_intr*1.1)],
# kLa=[0]*4+[120,120,60],
# DO_setpoints=[0,0,0,0,2.0,2.0,2.0],
kLa=[0]*4+[240,120,60],
DO_ID=DO_ID, DO_sat=SOSAT1,
suspended_growth_model=asm,
gas_stripping=gstrip)

C2 = su.FlatBottomCircularClarifier(
'C2', AS-0, ['effluent', 1-AS, 'WAS'],
**c2_kwargs
)
TC1 = su.IdealClarifier('TC1', C2-2, outs=['', 'thickened_WAS'],
sludge_flow_rate=30.9,
sludge_MLSS=7.0e4,)
M1 = su.Mixer('M1', ins=(C1-1, TC1-1))
# TC1 = su.Thickener('TC1', C2-2, outs=['thickened_sludge', ''],
# thickening_perc=7, TSS_removal_perc=98)
# M1 = su.Mixer('M1', ins=(C1-1, TC1-0))

# Switch to ADM1 components for the anaerobic digester
cmps_adm = pc.create_adm1p_cmps()
pc.create_adm1p_cmps()
thermo_adm = qs.get_thermo()
adm = pc.ADM1p(
f_bu_su=0.1328, f_pro_su=0.2691, f_ac_su=0.4076,
Expand All @@ -143,42 +152,48 @@ def create_system(flowsheet=None, default_init_conds=True):
J1 = su.mASM2dtoADM1p('J1', upstream=M1-0, thermo=thermo_adm, isdynamic=True,
adm1_model=adm, asm2d_model=asm)
AD = su.AnaerobicCSTR('AD', ins=J1.outs[0], outs=('biogas', 'AD_eff'), isdynamic=True,
V_liq=3400, V_gas=300, T=T_ad, model=adm,)
V_liq=3400, V_gas=300, T=T_ad, model=adm,
pH_ctrl=7.0,)
AD.algebraic_h2 = False
J2 = su.ADM1ptomASM2d('J2', upstream=AD-1, thermo=thermo_asm, isdynamic=True,
adm1_model=adm, asm2d_model=asm)
# Switch back to ASM1 components
qs.set_thermo(thermo_asm)

# Dewatering
C3 = su.Centrifuge(ID='C3', ins=J2-0, outs=['digested_sludge', ''],
thickening_perc=28, TSS_removal_perc=98)

M2 = su.Mixer('M2', ins=(TC1-1, C3-1), outs=1-C1)
# C3 = su.Centrifuge(ID='C3', ins=J2-0, outs=['digested_sludge', ''],
# thickening_perc=28, TSS_removal_perc=98)
# M2 = su.Mixer('M2', ins=(TC1-1, C3-1))

C3 = su.IdealClarifier('C3', J2-0, outs=['', 'digested_sludge'],
sludge_flow_rate=9.6,
sludge_MLSS=2.8e5,)
M2 = su.Mixer('M2', ins=(TC1-0, C3-0))
HD = su.HydraulicDelay('HD', ins=M2-0, outs=1-C1)

if default_init_conds:
C1.set_init_conc(**c1init)
# AS.set_init_conc(concentrations=asinit)
asdct = asinit.to_dict('index')
for i in (A1, A2, A3, A4, O1, O2, O3):
i.set_init_conc(**asdct[i.ID])
AS.set_init_conc(concentrations=asinit)
# asdct = asinit.to_dict('index')
# for i in (A1, A2, A3, A4, O1, O2, O3):
# i.set_init_conc(**asdct[i.ID])
C2.set_init_solubles(**c2init['s'])
C2.set_init_sludge_solids(**c2init['x'])
C2.set_init_TSS(c2init['tss'])
AD.set_init_conc(**adinit)

sys = qs.System('bsm2p',
path=(C1, A1, A2, A3, A4, O1, O2, O3, C2,
TC1, M1, J1, AD, J2, C3, M2),
recycle=(O3-0, C2-1, M2-0)
# path=(C1, AS, C2, TC1, M1, J1, AD, J2, C3, M2),
# recycle=(C2-1, M2-0)
# path=(C1, A1, A2, A3, A4, O1, O2, O3, C2,
# TC1, M1, J1, AD, J2, C3, M2, HD),
# recycle=(O3-0, C2-1, HD-0)
path=(C1, AS, C2, TC1, M1, J1, AD, J2, C3, M2, HD),
recycle=(C2-1, HD-0)
)
sys.set_tolerance(mol=1e-5, rmol=1e-5)
sys.maxiter = 5000
sys.set_dynamic_tracker(A1, O3, AD, C2-0)
# sys.set_tolerance(mol=1e-5, rmol=1e-5)
# sys.maxiter = 5000
# sys.set_dynamic_tracker(A1, O3, AD, C2-0)
# sys.set_dynamic_tracker(C1, AS, C2, J1, AD, J2, C3)
# sys.set_dynamic_tracker(AS, AD, C2-0)
sys.set_dynamic_tracker(AD, C2-0)

return sys

Expand Down Expand Up @@ -274,7 +289,7 @@ def run(sys, t, t_step, method=None, **kwargs):
print(f'Time span 0-{t}d \n')

sys.simulate(
state_reset_hook='reset_cache',
# state_reset_hook='reset_cache',
t_span=(0,t),
# t_eval=np.arange(0, t+t_step, t_step),
method=method,
Expand All @@ -291,13 +306,13 @@ def run(sys, t, t_step, method=None, **kwargs):
dct = globals()
dct.update(sys.flowsheet.to_dict())

t = 8
t = 500
t_step = 0.1
# method = 'RK45'
method = 'RK23'
# method = 'RK23'
# method = 'DOP853'
# method = 'Radau'
# method = 'BDF'
method = 'BDF'
# method = 'LSODA'

run(sys, t, t_step, method=method)
Expand Down

0 comments on commit cde6715

Please sign in to comment.