THIS PAGE IS FALSE, PLEASE SEE RCHIM.
Biochemical law for chemical reaction, whose PREPRO is LRCHIM.F and the law is implicit in the REACT element.
This constitutive law is used to take into account the degradation of organic matter and the following production of Volatil Fatty Acid (VFA). Those VFA are then consumed to produce methanogen biomass.
This law is based on McDougall’s model.
Le modèle de McDougall est modèle hydro-bio-mécanique de comportement des déchets.
La partie hydraulique du modèle est basé sur l’équation d’écoulement (« équation de Richard ») en milieu non saturée.
La partie biochimique décrit la phase anaérobique de la digestion de la matière organique. L’intégralité du phénomène de biodégradation est régit par le système d’équation suivant :
Hydrolyse et acidogénèse
Les acides gras volatils (AGV) sont générés en utilisant l’équation de génération suivante :
\[r_g = b \times \theta_e \times \phi \times P\]
Le premier est le taux maximal de génération $b$.
Le second est le la teneur en eau effective :
\[\theta_e = \frac{\theta - \theta_{res}}{\theta_{sat} - \theta_{res}}\]
Où :
Le troisième, est le terme de « digestibilité » qui dénote de la capacité du milieu à être biodégradé : \[\phi = 1 - \left[ \frac{Org_0 - Org}{Org_0}\right]^n\] La variation de la teneur organique en fonction du temps est donnée par : \[Org^{t+\Delta t} = Org^t - \theta \times \frac{162}{60}\times r_g\times \Delta t\] Où :
Le dernier terme est le terme d’inhibition qui tient compte de l’acidification du milieu et de l’accumulation des produits : \[P = e^{-k_{AVG}\times c}\] Où :
Acétogénèse et méthanogénèse
La génération du méthane est décrite par une cinétique de Monod :
\[r_j = \frac{k_0\times c}{k_{MC} + c}\times m\]
Où :
La génération du méthane est directement reliée à la consommation des AGV. Ainsi, le taux de dégradation des AGV est donné par : \[r_h = \frac{r_j}{Y}\]
Où :
La dégradation du méthane est quant à elle donnée par : \[r_k = k_2\times m\]
Où :
La génération, la dégradation et le transport des AGV sont compilés dans l’équation de bilan suivante : \[D_c \frac{\partial^2 c}{\partial x^2} + D_c \frac{\partial^2 c}{\partial z^2} - \frac{q_x}{\theta} \frac{\partial c}{\partial x} - \frac{q_z}{\theta} \frac{\partial c}{\partial z} + \left[ r_g - r_h\right] = \frac{\partial c}{\partial t}\] La génération, la dégradation et le transport du méthane sont compilés dans l’équation de bilan suivante : \[D_m \frac{\partial^2 m}{\partial x^2} + D_m \frac{\partial^2 m}{\partial z^2} - \frac{q_x}{\theta} \frac{\partial m}{\partial x} - \frac{q_z}{\theta} \frac{\partial m}{\partial z} + \left[ r_g - r_h\right] = \frac{\partial m}{\partial t}\]
Prepro: LRCHIM.F
Plane stress state | / |
Plane strain state | / |
Axisymmetric state | / |
3D state | / |
Generalized plane state | / |
Line 1 (2I5, 60A1) | |
---|---|
IL | Law number |
ITYPE | 199 |
COMMENT | Any comment (up to 60 characters) that will be reproduced on the output listing |
Line 1 (1I5) | |
---|---|
ICOAL | = 1 |
= 2 | |
= 3 to use the biochemical degradation of the organic matter | |
Line 2 (2G10.0) | |
FLUXF | |
FLUXD |
If ICOAL = 1
Line 1 (5G10.0) | |
---|---|
B | maximal hydrolysis rate |
CK_AGV | product inhibition factor |
CN | structural transformation parameter |
Y | cell/substrate yield coefficient |
CK_2 | methanogen death rate |
If ICOAL = 2
Line 1 (7G10.0) | |
---|---|
B | maximal hydrolysis rate |
CK_AGV | product inhibition factor |
CN | structural transformation parameter |
Y | cell/substrate yield coefficient |
CK_2 | methanogen death rate |
CK_0 | maximum specific growth rate |
ORG_INI | initial organic matter content |
If ICOAL = 3
Line 1 (7G10.0) | |
---|---|
B | maximal hydrolysis rate |
CK_AGV | product inhibition factor |
CN | structural transformation parameter |
Y | cell/substrate yield coefficient |
CK_2 | methanogen death rate |
CK_0 | maximum specific growth rate |
ORG_INI | initial organic matter content |
Line 2 (7G10.0) | |
CK_MC | half saturation constant |
POROSITY | porosity |
SAT_RES | residual saturation |
SAT_SAT | maximum saturation |
CN_VAN_Genuchten | first coefficient of Van Genuchten |
Cm_Van_Genuchten | second coefficient of Van Genuchten |
Alpha | pressure normalization factor |
Line 3 (4G10.0) | |
RHO_H2O | water density |
E_N | thermal flux per kg of degraded organic matter |
CM | initial condition on methanogen biomass concentration |
ORG | initial condition on organic matter content |
4
SIG(1) | consumption of water |
SIG(2) | production of gas |
SIG(3) | heat flux generated |
SIG(4) | VFA produced |
4
Q(1) | organic matter content |
Q(2) | methanogen biomass concentration |
Q(3) | modified enzymatic hydrolysis rate (VFA accumulation rate) |
Q(4) | VFA depletion rate |