====== RCHIM ======
FIXME THIS PAGE IS FALSE, PLEASE SEE [[laws:rchim|RCHIM]].
===== Description =====
Biochemical law for chemical reaction, whose PREPRO is LRCHIM.F and the law is implicit in the [[elements:react|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ù :
  * $\theta$ est la teneur en eau ;
  * $\theta_{sat} = \theta_{sat}\times n$ est la teneur en eau saturée, obtenue sur base des courbes de rétention d’eau ;
  * $\theta_{sat} = \theta_{sat}\times n$ est la teneur en eau résiduelle, obtenue sur base des courbes de rétention ;
  * $n$ est la porosité ;

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ù : 
  * $\frac{162}{60}$ est un rapport stœchiométrique ;
  * $\theta$ est la teneur en eau ;
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ù : 
  * $k_{AVG}$ est la constante d’inhibition choisie afin que la concentration maximale en AGV corresponde celles observées sur site ;
  * $c\ \left[g/m^3\ \text{d’eau}\right]$ est la concentration en acides gras volatils;
__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ù :
  * $k_0\ \left[s^{-1}\right]$ est le taux de croissance spécifique maximum ;
  * $k_{MC}\ \left[g/m^3\ \text{d’eau}\right]$ la constante de demi-saturation de Monod ;
  * $m\ \left[g/m^3\ \text{d’eau}\right]$ est la concentration en méthane ;

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ù :
  * $Y$ est un facteur de rendement pour contrôler les masses relatives ;
La dégradation du méthane est quant à elle donnée par :
\[r_k = k_2\times m\]

Où : 
  * $k_2$ est le taux de dégradation du méthane

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)^^
|IL|Law 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) ^^
|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|
===== 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 |