Table of Contents

RCHIM

FIXME THIS PAGE IS FALSE, PLEASE SEE RCHIM.

Description

Biochemical law for chemical reaction, whose PREPRO is LRCHIM.F and the law is implicit in the REACT element.

The model

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.

Modèles de McDougall (Mcdougall, 2007)

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}\]

Files

Prepro: LRCHIM.F

Availability

Plane stress state /
Plane strain state /
Axisymmetric state /
3D state /
Generalized plane state /

Input file

Parameters defining the type of constitutive law

Line 1 (2I5, 60A1)
ILLaw number
ITYPE 199
COMMENT Any comment (up to 60 characters) that will be reproduced on the output listing

Integer parameters

Line 1 (1I5)
ICOAL= 1
= 2
= 3 to use the biochemical degradation of the organic matter
Line 2 (2G10.0)
FLUXF
FLUXD

Real parameters

If ICOAL = 1

Line 1 (5G10.0)
Bmaximal hydrolysis rate
CK_AGVproduct inhibition factor
CNstructural transformation parameter
Ycell/substrate yield coefficient
CK_2methanogen death rate

If ICOAL = 2

Line 1 (7G10.0)
Bmaximal hydrolysis rate
CK_AGVproduct inhibition factor
CNstructural transformation parameter
Ycell/substrate yield coefficient
CK_2methanogen death rate
CK_0maximum specific growth rate
ORG_INIinitial organic matter content

If ICOAL = 3

Line 1 (7G10.0)
Bmaximal hydrolysis rate
CK_AGVproduct inhibition factor
CNstructural transformation parameter
Ycell/substrate yield coefficient
CK_2methanogen death rate
CK_0maximum specific growth rate
ORG_INIinitial organic matter content
Line 2 (7G10.0)
CK_MChalf saturation constant
POROSITYporosity
SAT_RESresidual saturation
SAT_SATmaximum saturation
CN_VAN_Genuchtenfirst coefficient of Van Genuchten
Cm_Van_Genuchtensecond coefficient of Van Genuchten
Alphapressure normalization factor
Line 3 (4G10.0)
RHO_H2Owater density
E_Nthermal flux per kg of degraded organic matter
CMinitial condition on methanogen biomass concentration
ORGinitial condition on organic matter content

Stresses

Number of stresses

4

Meaning (numerical, not physical)

SIG(1)consumption of water
SIG(2)production of gas
SIG(3)heat flux generated
SIG(4)VFA produced

State variables

Number of state variables

4

List of state variables

Q(1) organic matter content
Q(2) methanogen biomass concentration
Q(3)modified enzymatic hydrolysis rate (VFA accumulation rate)
Q(4) VFA depletion rate