Anaerobic Digestion Model No. 1 (ADM1)¶
- class qsdsan.processes.ADM1(components=None, path=None, N_xc=0.002686, N_I=0.004286, N_aa=0.007, f_ch_xc=0.2, f_pr_xc=0.2, f_li_xc=0.3, f_xI_xc=0.2, f_fa_li=0.95, f_bu_su=0.13, f_pro_su=0.27, f_ac_su=0.41, f_va_aa=0.23, f_bu_aa=0.26, f_pro_aa=0.05, f_ac_aa=0.4, f_ac_fa=0.7, f_pro_va=0.54, f_ac_va=0.31, f_ac_bu=0.8, f_ac_pro=0.57, Y_su=0.1, Y_aa=0.08, Y_fa=0.06, Y_c4=0.06, Y_pro=0.04, Y_ac=0.05, Y_h2=0.06, q_dis=0.5, q_ch_hyd=10, q_pr_hyd=10, q_li_hyd=10, k_su=30, k_aa=50, k_fa=6, k_c4=20, k_pro=13, k_ac=8, k_h2=35, K_su=0.5, K_aa=0.3, K_fa=0.4, K_c4=0.2, K_pro=0.1, K_ac=0.15, K_h2=7e-06, b_su=0.02, b_aa=0.02, b_fa=0.02, b_c4=0.02, b_pro=0.02, b_ac=0.02, b_h2=0.02, KI_h2_fa=5e-06, KI_h2_c4=1e-05, KI_h2_pro=3.5e-06, KI_nh3=0.0018, KS_IN=0.0001, pH_limits_aa=(4, 5.5), pH_limits_ac=(6, 7), pH_limits_h2=(5, 6), T_base=298.15, pKa_base=[14, 9.25, 6.35, 4.76, 4.88, 4.82, 4.86], Ka_dH=[55900, 51965, 7646, 0, 0, 0, 0], kLa=200, K_H_base=[0.00078, 0.0014, 0.035], K_H_dH=[-4180, -14240, -19410], **kwargs)¶
Anaerobic Digestion Model No.1. [1], [2]
- Parameters:
components (class:CompiledComponents, optional) – Components corresponding to each entry in the stoichiometry array, defaults to thermosteam.settings.chemicals.
path (str, optional) – Alternative file path for the Petersen matrix. The default is None.
N_xc (float, optional) – Nitrogen content of composite materials [kmol N/kg COD]. The default is 2.686e-3.
N_I (float, optional) – Nitrogen content of inert organics [kmol N/kg COD]. The default is 4.286e-3.
N_aa (float, optional) – Nitrogen content of amino acids [kmol N/kg COD]. The default is 7e-3.
f_ch_xc (float, optional) – Fraction of carbohydrates from composite disintegration [kg COD/kg COD]. The default is 0.2.
f_pr_xc (float, optional) – Fraction of proteins from composite disintegration [kg COD/kg COD]. The default is 0.2.
f_li_xc (float, optional) – Fraction of lipids from composite disintegration [kg COD/kg COD]. The default is 0.3.
f_xI_xc (float, optional) – Fraction of inert particulates from composite disintegration [kg COD/kg COD]. The default is 0.2.
f_fa_li (float, optional) – Fraction of long chain fatty acids (LCFAs) from hydrolysis of lipids [kg COD/kg COD]. The default is 0.95.
f_bu_su (float, optional) – Fraction of butyrate from sugars [kg COD/kg COD]. The default is 0.13.
f_pro_su (float, optional) – Fraction of propionate from sugars [kg COD/kg COD]. The default is 0.27.
f_ac_su (float, optional) – Fraction of acetate from sugars [kg COD/kg COD]. The default is 0.41.
f_va_aa (float, optional) – Fraction of valerate from amino acids [kg COD/kg COD]. The default is 0.23.
f_bu_aa (float, optional) – Fraction of butyrate from amino acids [kg COD/kg COD]. The default is 0.26.
f_pro_aa (float, optional) – Fraction of propionate from amino acids [kg COD/kg COD]. The default is 0.05.
f_ac_aa (float, optional) – Fraction of acetate from amino acids [kg COD/kg COD]. The default is 0.4.
f_ac_fa (float, optional) – Fraction of acetate from LCFAs [kg COD/kg COD]. The default is 0.7.
f_pro_va (float, optional) – Fraction of propionate from LCFAs [kg COD/kg COD]. The default is 0.54.
f_ac_va (float, optional) – Fraction of acetate from valerate [kg COD/kg COD]. The default is 0.31.
f_ac_bu (float, optional) – Fraction of acetate from butyrate [kg COD/kg COD]. The default is 0.8.
f_ac_pro (float, optional) – Fraction of acetate from propionate [kg COD/kg COD]. The default is 0.57.
Y_su (float, optional) – Biomass yield of sugar uptake [kg COD/kg COD]. The default is 0.1.
Y_aa (float, optional) – Biomass yield of amino acid uptake [kg COD/kg COD]. The default is 0.08.
Y_fa (float, optional) – Biomass yield of LCFA uptake [kg COD/kg COD]. The default is 0.06.
Y_c4 (float, optional) – Biomass yield of butyrate or valerate uptake [kg COD/kg COD]. The default is 0.06.
Y_pro (float, optional) – Biomass yield of propionate uptake [kg COD/kg COD]. The default is 0.04.
Y_ac (float, optional) – Biomass yield of acetate uptake [kg COD/kg COD]. The default is 0.05.
Y_h2 (float, optional) – Biomass yield of H2 uptake [kg COD/kg COD]. The default is 0.06.
q_dis (float, optional) – Composites disintegration rate constant [d^(-1)]. The default is 0.5.
q_ch_hyd (float, optional) – Carbohydrate hydrolysis rate constant [d^(-1)]. The default is 10.
q_pr_hyd (float, optional) – Protein hydrolysis rate constant [d^(-1)]. The default is 10.
q_li_hyd (float, optional) – Lipid hydrolysis rate constant [d^(-1)]. The default is 10.
k_su (float, optional) – Sugar uptake rate constant [d^(-1)]. The default is 30.
k_aa (float, optional) – Amino acid uptake rate constant [d^(-1)]. The default is 50.
k_fa (float, optional) – LCFA uptake rate constant [d^(-1)]. The default is 6.
k_c4 (float, optional) – Butyrate or valerate uptake rate constant [d^(-1)]. The default is 20.
k_pro (float, optional) – Propionate uptake rate constant [d^(-1)]. The default is 13.
k_ac (float, optional) – Acetate uptake rate constant [d^(-1)]. The default is 8.
k_h2 (float, optional) – H2 uptake rate constant [d^(-1)]. The default is 35.
K_su (float, optional) – Half saturation coefficient of sugar uptake [kg COD/m3]. The default is 0.5.
K_aa (float, optional) – Half saturation coefficient of amino acid uptake [kg COD/m3]. The default is 0.3.
K_fa (float, optional) – Half saturation coefficient of LCFA uptake [kg COD/m3]. The default is 0.4.
K_c4 (float, optional) – Half saturation coefficient of butyrate or valerate uptake [kg COD/m3]. The default is 0.2.
K_pro (float, optional) – Half saturation coefficient of propionate uptake [kg COD/m3]. The default is 0.1.
K_ac (float, optional) – Half saturation coefficient of acetate uptake [kg COD/m3]. The default is 0.15.
K_h2 (float, optional) – Half saturation coefficient of H2 uptake [kg COD/m3]. The default is 7e-6.
b_su (float, optional) – Decay rate constant of sugar-uptaking biomass [d^(-1)]. The default is 0.02.
b_aa (float, optional) – Decay rate constant of amino-acid-uptaking biomass [d^(-1)]. The default is 0.02.
b_fa (float, optional) – Decay rate constant of LCFA-uptaking biomass [d^(-1)]. The default is 0.02.
b_c4 (float, optional) – Decay rate constant of valerate- or butyrate-uptaking biomass [d^(-1)]. The default is 0.02.
b_pro (float, optional) – Decay rate constant of propionate-uptaking biomass [d^(-1)]. The default is 0.02.
b_ac (float, optional) – Decay rate constant of acetate-uptaking biomass [d^(-1)]. The default is 0.02.
b_h2 (float, optional) – Decay rate constant of H2-uptaking biomass [d^(-1)]. The default is 0.02.
KI_h2_fa (float, optional) – H2 inhibition coefficient for LCFA uptake [kg COD/m3]. The default is 5e-6.
KI_h2_c4 (float, optional) – H2 inhibition coefficient for butyrate or valerate uptake [kg COD/m3]. The default is 1e-5.
KI_h2_pro (float, optional) – H2 inhibition coefficient for propionate uptake [kg COD/m3]. The default is 3.5e-6.
KI_nh3 (float, optional) – Free ammonia inhibition coefficient for acetate uptake [M]. The default is 1.8e-3.
KS_IN (float, optional) – Inorganic nitrogen (nutrient) inhibition coefficient for soluble substrate uptake [M]. The default is 1e-4.
pH_limits_aa (2-tuple, optional) – Lower and upper limits of pH inhibition for acidogens and acetogens, unitless. The default is (4,5.5).
pH_limits_ac (2-tuple, optional) – Lower and upper limits of pH inhibition for aceticlastic methanogens, unitless. The default is (6,7).
pH_limits_h2 (2-tuple, optional) – Lower and upper limits of pH inhibition for H2-utilizing methanogens, unitless. The default is (5,6).
T_base (float, optional) – Base temperature for kinetic parameters [K]. The default is 298.15.
pKa_base (iterable[float], optional) – pKa (equilibrium coefficient) values of acid-base pairs at the base temperature, unitless, following the order of ADM1._acid_base_pairs. The default is [14, 9.25, 6.35, 4.76, 4.88, 4.82, 4.86].
Ka_dH (iterable[float], optional) – Heat of reaction of each acid-base pair at base temperature [J/mol], following the order of ADM1._acid_base_pairs. The default is [55900, 51965, 7646, 0, 0, 0, 0].
kLa (float, optional) – Liquid-gas mass transfer coefficient [d^(-1)]. The default is 200.
K_H_base (iterable[float], optional) – Henry’s Law coefficients of biogas species at the base temperature [M dissolved in liquid/bar]. Follows the order of ADM1._biogas_IDs. The default is [7.8e-4, 1.4e-3, 3.5e-2].
K_H_dH (iterable[float], optional) – Heat of reaction of liquid-gas transfer of biogas species [J/mol]. Follows the order of ADM1._biogas_IDs. The default is [-4180, -14240, -19410].
Examples
>>> from qsdsan import processes as pc >>> cmps = pc.create_adm1_cmps() >>> adm1 = pc.ADM1() >>> adm1.show() ADM1([disintegration, hydrolysis_carbs, hydrolysis_proteins, hydrolysis_lipids, uptake_sugars, uptake_amino_acids, uptake_LCFA, uptake_valerate, uptake_butyrate, uptake_propionate, uptake_acetate, uptake_h2, decay_Xsu, decay_Xaa, decay_Xfa, decay_Xc4, decay_Xpro, decay_Xac, decay_Xh2, h2_transfer, ch4_transfer, IC_transfer])
References
- check_stoichiometric_parameters()¶
Check whether product COD fractions sum up to 1 for each process.
- set_KI_nh3(K)¶
Set inhibition coefficient for free ammonia [M].
- set_KS_IN(K)¶
Set inhibition coefficient for inorganic nitrogen as a secondary substrate [M nitrogen].
- set_h2_inhibit_K(KI, process)¶
Set the H2 inhibition coefficient [kg/m3] for a process given its ID.
- set_half_sat_K(K, process)¶
Set the substrate half saturation coefficient [kg/m3] for a process given its ID.
- set_pH_inhibit_bounds(process, lower=None, upper=None)¶
Set the upper and/or lower limit(s) of pH inhibition [unitless] for a process given its ID.
- set_pKas(pKas)¶
Set the pKa values of the acid-base reactions at the base temperature.
- set_parameters(**parameters)¶
Set values to stoichiometric parameters in ADM1._stoichio_params.
- set_rate_constant(k, process)¶
Set the reaction rate constant [d^(-1)] for a process given its ID.