Modélisation D'un Bioréacteur
Modélisation D'un Bioréacteur
Modélisation D'un Bioréacteur
Description du processus
Volume: 50 l Temprature: 37C
Paramtres Type Matriau Longueur (m) Porosit(nm) Module membranaire tubulaire PVDF 1 30
8
1.2
Effluent traiter: Plusieurs transformations possibles. Hautement charg en matire organique Chaque transformation peut tre accomplie par une des diffrentes voies mtaboliques possibles Plusieurs essais de modlisation sont possibles
2
acidognes mthanognes
mthane et gaz carbonique permet de simuler les surcharges organiques qui se traduisent par une accumulation des AGV et une inhibition plus ou moins importante de la mthanogense.
3
Modlisation de linhibition par lhydrogne (Mosey, 1983) 4 populations bactriennes et 7 ractions: Population acidogne Population actogne Population mthanogne actotrophe Population mthanogne hydrognotrophe Rendre compte de linhibition de certaines populations bactriennes par lhydrogne Modle gnrique ADM1 (IWA Task Group for Mathematical Modelling of anaerobic Digestion Processes, 2002) 7 espces bactriennes et de nombreux substrats Malgr lexhaustivit des phnomnes quil peut dcrire, il est trs complexe: 26 quations et plus de 80 paramtres Frein lanalyse mathmatique du modle et son utilisation pour des stratgies de contrle
compromis entre la complexit du modle permettant de reprsenter tous les phnomnes connus et son adquation avec les informations exprimentales disponibles. Deux modles de la littrature ont t choisis et adapts un procd de bioracteur membrane:
Le modle un seul tage: 1 seule substrat et une seule population bactrienne Le modle 3 tages (Hill and Barth (1977), Simeonov (1996)): 3 substrats et 2 populations bactriennes
5
= 2
Avec
2 + +
Modle 3 tages
0 = 0 0 1 0 + 0 1 = (1 1 )1 1 = 1 1 + 1 0 1 1 1 2 = (2 2 )2 2 = 2 2 + 1 1 2 2 2
= 2 2
Avec
1 1 = 1
1 1 + 1
et
2 2 = 2
2 2 2 2 + 2 +
Rsolution
Modle un seul tage
Paramtre S0i Ax, As D K1 K2 Ks KI max j
-1
Dimension Interprtation g/l concentration de la matire organique dans linfluent coefficient de rtention des bactries et du substrat Taux de dilution taux de consommation du substrat coefficient de rendement de biogaz constante de demi-saturation des bactries constante dinhibition des bactries taux de croissance spcifique maximal des bactries
Hypothses
Racteur parfaitement mlang Pas de formation de gteau sur la membrane de filtration Filtration idale coefficients de slectivit de la membrane vis- vis des diffrents substances constants
Utilisant le FORTRAN et la mthode de Runge Kutta du 4me ordre Si on a La mthode RK4 est donne par lquation
h yn 1 yn (k1 2k2 2k3 k4 ) 6
10
Effet de Ks1
Sur S0
14 12 'd:\bi_cks1.dat'using 1:2 'd:\bi_cks2.dat'using 1:2 'd:\bi_cks3.dat'using 1:2 'd:\bi_cks4.dat'using 1:2 'd:\bi_cks5.dat'using 1:2 'd:\bi_cks6.dat'using 1:2 'd:\bi_cks7.dat'using 1:2
10
Ks1=200 Ks1=150
Ks1=110
K =50 Ks1=10Ks1=30 s1
0 100 200
Ks1=80
300 400 500 600 700 800
11
Effet de Ks1
Sur Q
90 80 70 60 50 40 30 20 10 0 0 100 200 300 400 500 600 700 800 'd:\bi_cks1.dat'using 1:7 'd:\bi_cks2.dat'using 1:7 'd:\bi_cks3.dat'using 1:7 'd:\bi_cks4.dat'using 1:7 'd:\bi_cks5.dat'using 1:7 'd:\bi_cks6.dat'using 1:7 'd:\bi_cks7.dat'using 1:7
12
Effet de KI
Sur S2
0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 0 100 200 300 400 500 600 700 800 'd:\bi_cki1.dat'using 1:4 'd:\bi_cki2.dat'using 1:4 'd:\bi_cki3.dat'using 1:4 'd:\bi_cki4.dat'using 1:4 'd:\bi_cki5.dat'using 1:4 'd:\bi_cki6.dat'using 1:4 'd:\bi_cki7.dat'using 1:4
13
Effet de KI
Sur Q
35 30 'd:\bi_cki1.dat'using 1:7 'd:\bi_cki2.dat'using 1:7 'd:\bi_cki3.dat'using 1:7 'd:\bi_cki4.dat'using 1:7 'd:\bi_cki5.dat'using 1:7 'd:\bi_cki6.dat'using 1:7 'd:\bi_cki7.dat'using 1:7
25
20
15
10
0 0 100 200 300 400 500 600 700 800 900 1000
14
Effet de
Sur Q
35 30 'd:\bi_beta0.dat'using 1:7 'd:\bi_beta1.dat'using 1:7 'd:\bi_beta2.dat'using 1:7 'd:\bi_beta3.dat'using 1:7
25
20
15
10
0 0 100 200 300 400 500 600 700 800 900 1000
15
Effet de Yp
Sur S2
0.14 0.12 'd:\bi_yp1.dat'using 1:4 'd:\bi_yp2.dat'using 1:4 'd:\bi_yp3.dat'using 1:4 'd:\bi_yp4.dat'using 1:4 'd:\bi_yp5.dat'using 1:4
0.1
0.08
0.06
0.04
0.02
0 0 100 200 300 400 500 600 700 800 900 1000
16
Effet de Yp
Sur Q
35 30 'd:\bi_yp1.dat'using 1:7 'd:\bi_yp2.dat'using 1:7 'd:\bi_yp4.dat'using 1:7 'd:\bi_yp3.dat'using 1:7 'd:\bi_yp5.dat'using 1:7
25
20
15
10
0 0 100 200 300 400 500 600 700 800 900 1000
17
Perspectives
Aprs simulation de ces deux modles, les donnes exprimentales vont permettre:
Loptimisation des paramtres les plus significatifs Validation des modles.
18