modelADM1_CSTR

Modelling anaerobic digestion model number 1 equations
Diagram of ADM1_CSTROpen the diagram
▶ Simulate in your browser

Extends from TransiEnt.Basics.Icons.ADM1 (Icon for ADM1 model).

Information

1. Purpose of model

This model contains different equatios for modeling the processes of anaerobic digestion. It can be chosen between the parameters of Bulkowska et al. and the bsm2 parameters.

2. Level of detail, physical effects considered, and physical insight

(Description)

3. Limits of validity

(Description)

4. Interfaces

(none)

5. Nomenclature

(no elements)

6. Governing Equations

(no equations)

7. Remarks for Usage

It is important to check if the parameters of the record and inflow parameters of the adm1 model fit to each other so that only bsm2 parameters or only bulkowska parameters are used.

8. Validation

This model has been validated as part of the master thesis of Philip Jahneke and it was tested in the check model "TestADM1" (look there for further information about the validation). Only the mesophilic operation mode has been validated.

9. References

[1] Bułkowska, K., et al. (2015). .ADM1-based modeling of anaerobic codigestion of maize silage and cattle manure – a feedstock characterisation for model implementation (part I)

[2] Rosén, C., & Jeppsson, U. (2006). Aspects on ADM1 Implementation within the BSM2 Framework

10. Version History

Model created by Philipp Jahneke (philipp.koziol@tuhhl.de), Sept 2018

Model adapted for TransiEnt by Jan Westphal (j.westphal@tuhh.de) in May 2020

Parameters

TypeNameDefaultDescription
RealPMPetersenMatrix.PM
RealMINModelica.Constants.epsvery small number used to avoid numerical instability of the acidbase equations, (biggest number so that 1.0+eps = 1.0)
RealY_proif operationMode == "thermophilic" then PetersenMatrix.Parameters.Y_pro_therm else PetersenMatrix.Parameters.Y_pro_meso
ADM1_Units.KineticRateConstantk_decif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_dec_therm else PetersenMatrix.Parameters.k_dec_mesorate constant of microbial decay
ADM1_Units.KineticRateConstantk_dec_acif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_dec_therm else PetersenMatrix.Parameters.k_dec_ac_mesodecay rate constant of acetate degraders
ADM1_Units.KineticRateConstantk_disif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_dis_therm else PetersenMatrix.Parameters.k_dis_mesodisintegration rate constant
ADM1_Units.KineticRateConstantk_hyd_chif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_hyd_ch_therm else PetersenMatrix.Parameters.k_hyd_ch_mesohydrolysis rate constant of carbohydrates
ADM1_Units.KineticRateConstantk_hyd_prif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_hyd_pr_therm else PetersenMatrix.Parameters.k_hyd_pr_mesohydrolysis rate constant of proteins
ADM1_Units.KineticRateConstantk_hyd_liif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_hyd_li_therm else PetersenMatrix.Parameters.k_hyd_li_mesohydrolysis rate constant of lipids
ADM1_Units.KineticRateConstantk_m_suif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_su_therm else PetersenMatrix.Parameters.k_m_su_mesoMonod maximum uptake rate constant for sugar degraders
ADM1_Units.KineticRateConstantk_m_aaif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_aa_therm else PetersenMatrix.Parameters.k_m_aa_mesoMonod maximum uptake rate constant for amino-acid degraders
ADM1_Units.KineticRateConstantk_m_faif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_fa_therm else PetersenMatrix.Parameters.k_m_fa_mesoMonod maximum uptake rate constant for LCFA degraders
ADM1_Units.KineticRateConstantk_m_c4if operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_c4_therm else PetersenMatrix.Parameters.k_m_c4_mesoMonod maximum uptake rate constant for valerate and butyrate degraders
ADM1_Units.KineticRateConstantk_m_proif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_pro_therm else PetersenMatrix.Parameters.k_m_pro_mesoMonod maximum uptake rate constant for propionate degraders
ADM1_Units.KineticRateConstantk_m_acif operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_ac_therm else PetersenMatrix.Parameters.k_m_ac_mesoMonod maximum uptake rate for acetate degraders
ADM1_Units.KineticRateConstantk_m_h2if operationMode == "thermophilic" then PetersenMatrix.Parameters.k_m_h2_therm else PetersenMatrix.Parameters.k_m_h2_mesoMonod maximum uptake rate constant for hydrogen degraders
ADM1_Units.ConcentrationCODK_S_suif operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_su_therm else PetersenMatrix.Parameters.K_S_su_mesosugar concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_aaif operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_aa_therm else PetersenMatrix.Parameters.K_S_aa_mesoamino acid concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_faif operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_fa_therm else PetersenMatrix.Parameters.K_S_fa_mesoLCFA concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_c4if operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_c4_therm else PetersenMatrix.Parameters.K_S_c4_mesoValerate/Butyrate concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_proif operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_pro_therm else PetersenMatrix.Parameters.K_S_pro_mesoPropionate concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_acif operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_ac_therm else PetersenMatrix.Parameters.K_S_ac_mesoAcetate concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.ConcentrationCODK_S_h2if operationMode == "thermophilic" then PetersenMatrix.Parameters.K_S_h2_therm else PetersenMatrix.Parameters.K_S_h2_mesoHydrogen concentration below which uptake is reduced to less than half of maximum uptake
ADM1_Units.KineticRateConstantk_m{k_m_su, k_m_aa, k_m_fa, k_m_c4, k_m_pro, k_m_ac, k_m_h2}vector of Monod maximum uptake rates
RealComponent_in{S_su_in, S_aa_in, S_fa_in, S_va_in, S_bu_in, S_pro_in, S_ac_in, S_h2_in, S_ch4_in, S_IC_in, S_IN_in, S_I_in, X_c_in, X_ch_in, X_pr_in, X_li_in, X_su_in, X_aa_in, X_fa_in, X_c4_in, X_pro_in, X_ac_in, X_h2_in, X_I_in}vector of Inflow Concentrations
SI.VolumeV_liquidgeometryCSTR.V_fluidVolume of Substrate in CSTR
SI.VolumeV_gasgeometryCSTR.V_gasGas Volume inside CSTR
SI.Timet_res20*86400mean residence time of substrate
BooleanComponentInputfalse=true, if Input Concentration is defined by input
BooleanuseBSM2truetrue if BSM2 parameters shall be used
StringoperationMode"mesophilic"
SI.ConcentrationCat40
SI.ConcentrationAn20Concentration of Cations and Anions of strong acids or base, whose equilibrium is assumed not to be affected by pH-Value
SI.MassConcentrationS_su_in0.01Concentration of Monosaccharides in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_aa_in0.001Concentration of Amino Acids in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_fa_in0.001Concentration of Long Chain Fatty Acids in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_va_in0.001Concentration of total Valerate in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_bu_in0.001Concentration of total Butyrate in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_pro_in0.001Concentration of total Propionate in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_ac_in0.001Concentration of total Acetate in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_h2_in1e-8Concentration of dissolved Hydrogen gas in inflowing substrate in kgCOD/m^3
SI.MassConcentrationS_ch4_in1e-5Concentration of dissolved Methane gas in inflowing substrate in kgCOD/m^3
SI.ConcentrationS_IC_in40Concentration of Inorganic Carbon in inflowing substrate in mol/m^3
SI.ConcentrationS_IN_in10Concentration of Inorganic Nitrogen in inflowing substrate in mol/m^3
SI.MassConcentrationS_I_in0.02Concentration of Soluble Inerts in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_c_in2Concentration of Composites in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_ch_in5Concentration of Carbohydrates in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_pr_in20Concentration of Proteins in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_li_in5Concentration of Lipids in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_su_in0.00Concentration of Sugar degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_aa_in0.01Concentration of Amino Acid degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_fa_in0.01Concentration of LCFA degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_c4_in0.01Concentration of Valerate and Butyrate degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_pro_in0.01Concentration of Propionate degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_ac_in0.01Concentration of Acetate degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_h2_in0.01Concentration of Hydrogen degraders in inflowing substrate in kgCOD/m^3
SI.MassConcentrationX_I_in25Concentration of Particulate Inerts in inflowing substrate in kgCOD/m^3

Connectors

TypeNameDefaultDescription
Modelica.Blocks.Interfaces.RealOutput[24]Components{S_su, S_aa, S_fa, S_va, S_bu, S_pro, S_ac, S_h2, S_ch4, S_IC, S_IN, S_I, X_c, X_ch, X_pr, X_li, X_su, X_aa, X_fa, X_c4, X_pro, X_ac, X_h2, X_I}output composition of outflow which equals composition in CSTR
Modelica.Blocks.Interfaces.RealInput[24]Components_Ininput composition if set from outside or taken outflow composition of other CSTR as input

Components

TypeNameDefaultDescription
TransiEnt.Producer.Gas.BiogasPlant.Base.GeometryCSTRgeometryCSTR
SI.Pressurep_atm101300pressure of gas after leaving the CSTR through port
SI.TemperatureTTemperatur inside reactor
SI.DensityrhoDensity of substrate
SI.VolumeFlowRateq_gasgas volume flow at inside pressure
SI.VolumeFlowRateq_gas_outq_gas*p/p_atmgas volume flow at nominal pressure
SI.VolumeFlowRateq_liquid_inV_liquid/t_resSubstrate Volume flow flowing into the CSTR
SI.VolumeFlowRateq_liquid_outV_liquid/t_resSubstrate Volume flow leaving the CSTR
BooleangasValveOpen
SI.MassConcentrationS_suConcentration of Monosaccharides in kgCOD/m^3
SI.MassConcentrationS_aaConcentration of Amino Acids in kgCOD/m^3
SI.MassConcentrationS_faConcentration of Long Chain Fatty Acids in kgCOD/m^3
SI.MassConcentrationS_vaConcentration of total Valerate in kgCOD/m^3
SI.MassConcentrationS_buConcentration of total Butyrate in kgCOD/m^3
SI.MassConcentrationS_proConcentration of total Propionate in kgCOD/m^3
SI.MassConcentrationS_acConcentration of total Acetate in kgCOD/m^3
SI.MassConcentrationS_h2Concentration of dissolved Hydrogen gas in kgCOD/m^3
SI.MassConcentrationS_ch4Concentration of dissolved Methane gas in kgCOD/m^3
SI.ConcentrationS_ICConcentration of Inorganic Carbon in mol/m^3
SI.ConcentrationS_INConcentration of Inorganic Nitrogen in mol/m^3
SI.MassConcentrationS_IConcentration of Soluble Inerts in kgCOD/m^3
SI.MassConcentrationX_cConcentration of Composites in kgCOD/m^3
SI.MassConcentrationX_chConcentration of Carbohydrates in kgCOD/m^3
SI.MassConcentrationX_prConcentration of Proteins in kgCOD/m^3
SI.MassConcentrationX_liConcentration of Lipids in kgCOD/m^3
SI.MassConcentrationX_suConcentration of Sugar degraders in kgCOD/m^3
SI.MassConcentrationX_aaConcentration of Amino Acid degraders in kgCOD/m^3
SI.MassConcentrationX_faConcentration of LCFA degraders in kgCOD/m^3
SI.MassConcentrationX_c4Concentration of Valerate and Butyrate degraders in kgCOD/m^3
SI.MassConcentrationX_proConcentration of Propionate degraders in kgCOD/m^3
SI.MassConcentrationX_acConcentration of Acetate degraders in kgCOD/m^3
SI.MassConcentrationX_h2Concentration of Hydrogen degraders in kgCOD/m^3
SI.MassConcentrationX_IConcentration of Particulate Inerts in kgCOD/m^3
SI.MassConcentrationS_gas_h2Concentration of Hydrogen in Gas-Phase in kgCOD/m^3
SI.MassConcentrationS_gas_ch4Concentration of Methane in Gas-Phase in kgCOD/m^3
SI.ConcentrationS_gas_co2Concentration of Carbondioxide in Gas-Phase in mol/m^3
SI.PartialPressurep_h2Partial Pressure of Hydrogen in gas phase
SI.PartialPressurep_ch4Partial Pressure of Methane in gas phase
SI.PartialPressurep_co2Partial Pressure of Carbondioxide in gas phase
SI.PartialPressurep_h2oPartial Pressure of Water Vapor in gas phase
SI.PressurepAbsolute pressure of gas inside Reactor
SI.ConcentrationS_cat_pConcentration of additional Cations not related to ADM1 components in solution
SI.ConcentrationS_an_nConcentration of additional Anions not relted to ADM1 components in solution
SI.ConcentrationS_Va_nConcentration of Valerate Anions in solution
SI.ConcentrationS_Bu_nConcentration of Butyrate Anions in solution
SI.ConcentrationS_Pro_nConcentration of Propionate Anions in solution
SI.ConcentrationS_Ac_nConcentration of Acetate Anions in solution
SI.ConcentrationS_HCO3_nConcentration of Hydrogencarbonate Anions in solution
SI.ConcentrationS_NH3_aqConcentration of Ammonia in solution
SI.ConcentrationS_NH4_pConcentration of Ammonium Cations in solution
SI.ConcentrationS_CO2_aqConcentration of Carbondioxide in solution
SI.ConcentrationS_H_pConcentration of Protons in solution
SI.ConcentrationPHICharge Balance without Protons and Hydroxide Anions
RealpH3 - log10(S_H_p)pH Wert if using mol per cubic meter instead of mol per litre
RealKwTransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.Kw_ref, PetersenMatrix.Parameters.deltH0_w, T)Autoprotolysis Constant of Water
SI.ConcentrationKa_CO2TransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.Ka_CO2_ref, PetersenMatrix.Parameters.deltH0_CO2, T)Acid dissociation constant of CO2 in water
SI.ConcentrationKa_NH4TransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.Ka_NH4_ref, PetersenMatrix.Parameters.deltH0_NH4, T)Acid dissociation constant of Ammonium Ions
SI.ConcentrationKa_HVaPetersenMatrix.Parameters.Ka_HVa_refAcid dissociation constant of Valyric Acid
SI.ConcentrationKa_HBuPetersenMatrix.Parameters.Ka_HBu_refAcid dissociation constant of Butyric Acid
SI.ConcentrationKa_HProPetersenMatrix.Parameters.Ka_HPro_refAcid dissociation constant of Propionic Acid
SI.ConcentrationKa_HAcPetersenMatrix.Parameters.Ka_HAc_refAcid dissociation constant of Acetic Acid
RealKH_h2TransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.KH_h2_ref, PetersenMatrix.Parameters.deltH0_H_h2, T)Henry Coefficient for liquid gas transfer of Hydrogen
RealKH_ch4TransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.KH_ch4_ref, PetersenMatrix.Parameters.deltH0_H_ch4, T)Henry Coefficient for liquid gas transfer of Methane
RealKH_co2TransiEnt.Producer.Gas.BiogasPlant.Base.ADM1.ADM1_Functions.VantHoffEquation(PetersenMatrix.Parameters.KH_co2_ref, PetersenMatrix.Parameters.deltH0_H_co2, T)Henry Coefficient for liquid gas transfer of Carbondioxide
RealI_suCoefficient describing inhibition of sugar degraders
RealI_aaCoefficient describing inhibition of amino acid degraders
RealI_faCoefficient describing inhibition of LCFA degraders
RealI_c4Coefficient describing inhibition of Butyrate and Valerate degraders
RealI_proCoefficient describing inhibition of Propionate degraders
RealI_acCoefficient describing inhibition of Acetate degraders
RealI_h2Coefficient describing inhibition of hydrogen degraders
RealI{I_su, I_aa, I_fa, I_c4, I_pro, I_ac, I_h2}vector of Inhibition constants
RealComponent{S_su, S_aa, S_fa, S_va, S_bu, S_pro, S_ac, S_h2, S_ch4, S_IC, S_IN, S_I, X_c, X_ch, X_pr, X_li, X_su, X_aa, X_fa, X_c4, X_pro, X_ac, X_h2, X_I}Compostion vector of substrate inside CSTR
RealDecayRatecat(1, k_dec*{X_su, X_aa, X_fa, X_c4, X_pro}, {k_dec_ac*X_ac, k_dec*X_h2})Decay Rates of microorganisms
RealProcessRate{k_dis, k_hyd_ch, k_hyd_pr, k_hyd_li}.*{X_c, X_ch, X_pr, X_li}Process Rates of disintegration and hydrolysis
RealMonodRateRates of metabolism of microorganisms
RealTransferRateliquid to gas transfer rates
RealReactionsvector containing all Conversion Reactions of ADM1
SI.EnergyFlowRateQ_flow_REnergy Flow released by metabolism reactions
Modelica.Units.SI.MassConcentrationVSorganic fraction of solids
Modelica.Units.SI.MassConcentrationVS_inorganic fraction of solids in inflow
ADM1_PetersenMatrixPetersenMatrixPetersen Matrix of ADM1