# BIODG

**BIODG**           input data block introducing parameters needed to simulate the biodegradation processes.

Record **BIODG.1**

&#x20;                       Free format for 7 parameters, or Format (I5,I5,E10.4,10X,4E10.4)

&#x20;                       IMONOD, ICFLAG, BFAC, SW1, SW2, WEA, WSUB

*IMONOD*          selects between multiplicative and minimum Monod model for the substrate degradation rate equation.

&#x20;        0 : multiplicative Monod model (default).

&#x20;        \>0 : minimum Monod model.

*ICFLAG*          selects how to consider the competitive and Haldane inhibition terms in the Monod model.

&#x20;       \=0: competitive and Haldane inhibition factors are applied to all Monod terms in the substrate degradation rate equation, including the electron acceptor and nutrients terms as in the BIOMOC code formulation (default).

&#x20;       \>0: competitive and Haldane inhibition factors are applied only to the substrate Monod term, as in conventional formulations of cometabolic degradation processes modelled using competitive inhibition effects.

*BFAC*               reduction factor criterion for local Newton-Raphson iteration to reduce substrate residual. Nonlinearity of the reaction is handled by local Newton-Raphson iteration for the primary substrate.  *BFAC* is the factor by which the residual must be reduced before convergence.

*SW1*                 lower limit of aqueous phase saturation considered in the saturation inhibition function (if =0, the default value 0.02 will be applied.).

*SW2*                 upper limit of aqueous phase saturation considered in the saturation inhibition function (SW1 < SW2 ≤ 1).

*WEA*                 weighting factor for the linear interpolation of electron acceptor and nutrients concentrations to be used in the substrate degradation rate equation (0 < WEA ≤ 1). Default value is WEA = 0.5. WEA = 1 corresponds to using the concentration evaluated at the end of the time step.

*WSUB*               weighting factor for the linear interpolation of substrate concentration to be used in the substrate degradation rate equation (0 < WSUB ≤ 1). Default value is WSUB = 0.5. WSUB=1 corresponds to using the concentration evaluated at the end of the time step.

Record **BIODG.2**

&#x20;                       Free format for 1 parameter, or Format (I5)

&#x20;                       NPROC

*NPROC*            number of biodegradation processes in the simulation (NPROC ≤ 20).

Record **BIODG.3**

&#x20;                       Free format for 8 parameters, or Format (2(5X, I5), 2E10.4, 3(5X, I5), E10.4)                      &#x20;

NSPROC(IP), IBIO(IP), AMUMAX(IP), YIELD(IP), NCOMP(IP), NNC(IP), NHAL(IP), ENTB(IP)

*NSPROC(IP)*   number of mass components controlling the substrate degradation rate in process IP.

*IBIO(IP)*          index of microbial population involved in process IP.

&#x20;*AMUMAX(IP)*         maximum specific substrate degradation rate in process IP (kg substrate/ (s kg biomass) ).

*YIELD(IP)*        yield coefficient for the growth of biomass due to the degradation of unit mass of substrate in process IP (kg biomass / kg substrate).

*NCOMP(IP)*     number of mass components responsible for competitive inhibition in process IP.

*NNC(IP)*          number of mass components responsible for non-competitive inhibition in process IP.

*NHAL(IP)*          number of mass components responsible for Haldane inhibition in process IP.

*ENTB(IP)*          heat of reaction for the degradation of substrate in process IP (J/kg substrate).

Record **BIODG.4**

&#x20;                       Free format for as more as 12 parameters, or Format ( 6(5X, I5, E10.4) )

&#x20;                       KSPROC(IP,NC), AKS(IP,NC) ; NC=1, NSPROC(IP)

*KSPROC(IP,NC)*     index of mass component # NC controlling the substrate degradation rate in process IP. By convention, the first (NC=1) must refer to the primary substrate of process IP.

*AKS(IP,NC)*     half saturation constant of mass component # NC in the Monod term of substrate degradation rate in process IP (kg substrate / kg aqueous phase). If specified, must be >0.

Record **BIODG.5**

&#x20;                       Free format for as more as 12 parameters, or Format ( 6(5X, I5, E10.4) )

&#x20;                       KSCOMP(IP,NC), AKCOMP(IP,NC) ; NC=1, NCOMP(IP)

*KSCOMP(IP,NC)*     index of mass component # NC responsible for competitive inhibition in process IP.

*AKCOMP(IP,NC)*     competitive inhibition coefficient for mass component # NC in process IP (kg solute / kg aqueous phase).

Record **BIODG.6**

&#x20;                       Free format for as more as 12 parameters, or Format ( 6(5X, I5, E10.4) )

&#x20;                       KSNC(IP,NC), AKNC(IP,NC) ; NC=1, NNC(IP)

*KSNC(IP,NC)*          index of mass component # NC responsible for non-competitive inhibition in process IP.

*AKNC(IP,NC)*          non-competitive inhibition coefficient for mass component # NC in process IP (kg solute / kg aqueous phase).

Record **BIODG.7**

&#x20;                       Free format for as more as 12 parameters, or Format ( 6(5X, I5, E10.4) )

&#x20;                       KSHAL(IP,NC), AKHAL(IP,NC) ; NC=1, NHAL(IP)

*KSHAL(IP,NC)*       index of mass component # NC responsible for Haldane inhibition in process IP.

*AKHAL(IP,NC)*       Haldane inhibition coefficient for mass component # NC in process IP (kg solute /kg aqueous phase).

Record **BIODG.8**

&#x20;                       Free format for as more as 20 parameters, or Format (8E10.4)

&#x20;                       UPTAKE(IP,L) ; L=1,NumCom

*UPTAKE(IP,L)*       uptake coefficient of component # L in process IP with respect to 1 mole of degraded primary substrate (mole component / mole substrate).

Repeat records BIODG.3 through BIODG.8 for a total of NPROC (up to 20) different biodegradation processes (IP= 1, NPROC).

Record **BIODG.9**

&#x20;                       Free format for 2 parameters, or Format (2I5)

&#x20;                       NPOP, IBFLAG

*NPOP*               number of microbial populations (NPOP ≤ 20).

*IBFLAG*          specifies if initial conditions for biomass concentration are supplied through input blocks INDOM, INCON, or through an INCON file.

0: constant biomass concentrations are given as initial conditions through block BIODG.10.

\>0: variable biomass concentrations are supplied as initial conditions through blocks INDOM, INCON or through an INCON file.

Record **BIODG.10**

&#x20;                       Free format for 5 parameters, or Format (5E10.4)

&#x20;                       BAI(IB), BA0(IB), TMAX(IB), DEATH(IB), AKBIO(IB)

*BAI(IB)*          initial concentration (valid for the entire simulation grid) in the aqueous phase of microbial population # IB (kg biomass / kg aqueous phase).

*BA0(IB)*          minimum concentration in the aqueous phase of microbial population # IB enforced during the simulation (kg biomass / kg aqueous phase).

*TMAX(IB)*       maximum temperature for the calculation of temperature inhibition function in the substrate degradation rate equation (°C).

*DEATH(IB)*      death rate constant, or maintenance constant, for the microbial population # IB (s-1).

*AKBIO(IB)*        inhibition constant for biomass growth of microbial population # IB (kg biomass kg aqueous phase).

Repeat record BIODG.10 for a total of NPOP (up to 20) different microbial populations (IB= 1, NPOP).

**Used in**: All EOS modules

**Example**:

*BIODG*----1----*----2----*----3----*----4----*----5----*----6----*----7----\*----8

*0,  1.e-10,  ,  0.0 ,  .2,  0.9,  0.9*                          //BIODG.1

*2*                                                                              //BIODG.2&#x20;

*2, 1, 1.2732e-5, 1.0, 0, 0, 0, 0.0*                         //BIODG.3 for process 1&#x20;

*4, 2.e-6, 2, 1.e-6*                                                 //BIODG.4 for process 1&#x20;

*0*                                                                            //BIODG.5 for process 1&#x20;

*0*                                                                          //BIODG.6 for process 1&#x20;

*0*                                                                         //BIODG.7 for process 1&#x20;

*-4.0, 9.0, 0.0, 1.0*                                               //BIODG.8 for process 1    &#x20;

*2, 1, 1.1574e-5, 0.8, 0, 0, 0, 0.0*                      //BIODG.3 for process 2&#x20;

*3, 1.e-6, 2, 0.5e-6*                                              //BIODG.4 for process 2&#x20;

*0*                                                                          //BIODG.5 for process 2

*0*                                                                           //BIODG.6 for process 2&#x20;

*0*                                                                           //BIODG.7 for process 2&#x20;

*-3.0, 7.5, 1.0,  0.0*                                              //BIODG.8 for process 2              &#x20;

*1*                                                                           // BIODG.9      &#x20;

*1.530e-6, 1.00e-6, 40.,  2.3148e-7, 0.0*       //BIODG.10
