Una alternativa numérica en la solución de un sistema que modela la producción de biogás

<!-- metadata commented in wiki content
==Una alternativa numérica en la solución de un sistema que modela la producción de biogás==

'''Maximiliano Machado H 

Facultad de Ciencias Naturales y Matemáticas  

Universidad de Ibagué 

Calle 67, Carrera 22, Ambalá, CP 730001, Ibagué, Colombia,

Tel: +57(8)2709400, Fax: Tel: +57(8)2709401.

e-mail: maximiliano.machado@unibague.edu.co  Laura H Alvarez G 

Facultad de Ingenierías

Programa de Ingeniería Industrial 

Universidad de Ibagué 

Calle 67, Carrera 22, Ambalá, CP 730001, Ibagué, Colombia,

Tel: +57(8)2709400, Fax: Tel: +57(8)2709401.

e-mail: 2320111031@estudiantesunibague.edu.co'''

==Resumen==

Se estudia un modelo simplificado del Anaerobic Digestion Model Number One (ADM1) que representa la Digestión Anaerobia para la producción de biogás y tratamiento de aguas residuales. Este sistema está modelado por sistemas de ecuaciones diferenciales no lineales. Primero se implementa el Análisis Dimensional para obtener un modelo equivalente con variables y parámetros adimensionales.  Luego se soluciona numéricamente tanto el modelo inicial como el modelo adimensional y se realiza un análisis, numérico y gráfico, comparativo entre sus correspondientes soluciones numéricas en donde se puede observar las bondades de trabajar con el modelo adimensional con respecto al modelo original.

''Palabras Clave.'' Análisis dimensional, biogás, digestión anaerobia, reescalamiento, teorema de Pi-Buckingham.

==A Numerical Alternative In The Solution Of A System That Models The Biogas Production    ==
-->

==Abstract==

We study a simplified model of the Anaerobic Digestion Model Number One (ADM1) that represents Anaerobic Digestion for the biogas production and wastewater treatment. This system is modeled by systems of nonlinear differential equations. First, we implemented the Dimensional Analysis to obtain an equivalent model with non-dimensional variables and parameters. Then, numerical solutions are solved both the initial model and the non-dimensional model, and we comparative graph between their corresponding numerical solutions. where we see the advantages of working with the non-dimensional model with respect to the original model.

''Keywords.'' Anaerobic digestion, biogas, dimensional analysis, Pi-Buckingham theorem.

==1 Introducción==

Iniciamos recordando que dentro del objetivo de crear una plataforma común y facilitar la simulación de una gran variedad de procesos anaerobios, la Asociación Internacional del agua IWA (2002) publicó el Modelo para Digestión Anaerobia No. 1 (ADM1) con el cual pretende proporcionar una base unificada para la modelización de la digestión anaerobia. De este modelo se derivan otros modelos como los de lodos activados (ASM1, ASM2, ASM2d, ASM3), entre otros <span id='citeF-8'></span>[[#cite-8|[8]]].

En la literatura se encuentra como aplicación de la digestión anaerobia el tratamiento de aguas industriales <span id='citeF-3'></span>[[#cite-3|[3]]], aguas residuales <span id='citeF-11'></span>[[#cite-11|[11]]],<span id='citeF-14'></span>[[#cite-14|[14]]] y en la biodegradación de aguas en el despulpado del café <span id='citeF-10'></span>[[#cite-10|[10]]],<span id='citeF-19'></span>[[#cite-19|[19]]], caña de azúcar <span id='citeF-13'></span>[[#cite-13|[13]]] entre otras. De igual forma se han hecho estudios sobre su modelacion <span id='citeF-16'></span>[[#cite-16|[16]]],<span id='citeF-20'></span>[[#cite-20|[20]]],<span id='citeF-21'></span>[[#cite-21|[21]]],<span id='citeF-22'></span>[[#cite-22|[22]]] y de análisis de bifurcación y estabilidad <span id='citeF-2'></span>[[#cite-2|[2]]],<span id='citeF-9'></span>[[#cite-9|[9]]] y de estrategias de optimización apropiadas <span id='citeF-18'></span>[[#cite-18|[18]]].

La digestión anaerobia es un proceso muy complejo debido a las reacciones bioquímicas que en ella ocurren y por la cantidad y variedad de bacterias involucradas, las cuales difieren en cuanto a su velocidad de crecimiento y a la sensibilidad a cada producto intermedio como los inhibidores (hidrógeno, ácido acético o amoniaco producido en la acidogénesis de aminoácidos). El proceso de digestión anaerobia está dividido en cuatro etapas subsiguientes: hidrólisis, acidogénesis, acetogénesis y metanogénesis, en las cuales se debe vigilar el equilibrio del proceso para evitar la acumulación de compuestos intermedios inhibidores que generan unos bajos niveles de pH y a la vez impiden la producción óptima de biogás <span id='citeF-5'></span>[[#cite-5|[5]]],<span id='citeF-6'></span>[[#cite-6|[6]]], otra de las aplicaciones directas de la digestión anaerobia.

El biogás es una mezcla de gases, considerado una fuente de energía renovable que se obtiene como resultado de la fermentación de materia orgánica con altas concentraciones de humedad, como lo son los desechos de animales, residuos vegetales, basuras domésticas, aguas residuales industriales o domésticas, entre otros; en condiciones anaeróbicas mediante la digestión de un grupo de bacterias, que puede suceder en condiciones naturales o en dispositivos específicos como los reactores. La composición del biogás depende principalmente del material utilizado en la digestión anaerobia y las condiciones en las que se realiza el proceso, sin embargo, está constituido aproximadamente por Metano entre 50% y 70%, Dióxido de Carbono entre 30% y 50% y un porcentaje menor al 5% de otros gases considerados impurezas (hidrógeno, nitrógeno, oxígeno, sulfuro de hidrógeno o ácido sulfhídrico). De este proceso también se obtiene una solución acuosa o lodo en donde se encuentran los microorganismos que se encargan de la degradación de la materia orgánica <span id='citeF-13'></span>[[#cite-13|[13]]].

Luego de la publicación del modelo ADM1 el doctor Olivier Bernard del INRIA de Francia, junto con otros investigadores, publica un artículo en donde plantea un modelo matemático macroscópico simplificado que representa la dinámica del tratamiento de aguas residuales por medio de dos de las reacciones de la digestión anaerobia: la acidogénesis y la metanogénesis; y que además conserva parámetros cinéticos e incorpora equilibrios electroquímicos para poder incluir la alcalinidad, propiedad que desempeña un papel muy importante en relación con el monitoreo y estrategias de control de la planta. De igual manera propone un método de identificación de parámetros para la digestión anaerobia; a este modelo de denominó AMOCO <span id='citeF-3'></span>[[#cite-3|[3]]].

De todos estos modelos simplificados del ADM1 que existen, consideramos el AMOCO debido a que, ha sido uno de los modelos más completos en la inclusión de parámetros cinéticos, después de ADM1.

Este modelo tiene las siguientes características:

* Es un modelo dinámico sobre el desarrollo y la identificación de los parámetros del proceso de digestión anaerobia.
* El balance de masas se ha considerado a dos pasos, acidogénesis y metanogénesis.
* El modelo incorpora los equilibrios electroquímicos a fin de incluir la alcalinidad.
* El modelo trabaja con la medida de demanda química de oxígeno (DQO), la cual determina la cantidad de oxígeno requerido para oxidar la materia orgánica en una muestra de agua residual.

Entre los factores ambientales que afectan el proceso de digestión anaerobia se encuentra el pH y el control de éste, en los sistemas anaerobios, es de vital importancia para el desarrollo del proceso, debido a la gran suceptibilidad de los microorganismos, especialmente las metanobacterias, frente a las especies ácidas. El equilibrio ácido-base que tiene lugar en la operación de los digestores anaerobios es muy importante por la presencia de diversas especies de esta naturaleza que están en el medio y que requieren ser neutralizadas para restituir el <math display="inline">pH</math> del medio, por lo que junto con este parámetro debe considerarse a la alcalinidad que es la forma de control del <math display="inline">pH</math> dice <span id='citeF-6'></span>[[#cite-6|[6]]]. Por tanto, se considera la alcalinidad <math display="inline">A</math> como variable en el modelo.

Como se menciona en <span id='citeF-7'></span>[[#cite-7|[7]]], para la estabilidad del pH es importante el equilibrio <math display="inline">CO_2</math>-bicarbonato. Dado que el total de carbono inorgánico es la suma de las especies de carbono inorgánico en una solución, esto es, dioxido de carbono (<math display="inline">{CO}_2</math>), acido carbónico, bicarbonato y carbonato.

El total de Carbono inorgánico total es un parámetro que se utiliza para valorar la calidad de las aguas de un determinado lugar y debido a que éste se mide con la acidificación de una muestra de agua, lo que lleva a un equilibrio del <math display="inline">{CO}_2</math>, entonces otra de las variables en el modelo será el total de carbono inorgánico <math display="inline">C</math>.

===1.1 El modelo===

Un esquema en el que podemos observar la reacciones en el proceso es:

&nbsp; <math>S_1</math> &nbsp;<math display="inline">\stackrel{r_1}{\longrightarrow }</math> &nbsp; <math>S_2</math> + <math>{CO}_{2}</math>

&nbsp; <math>S_2</math> &nbsp;<math display="inline">\stackrel{r_2}{\longrightarrow }</math> &nbsp; <math>{CH}_{4} + {CO}_{2}</math>

donde:

* <math display="inline">r_{1} = \mu _1(S_1) X_1</math> es la tasa de reacción en la acetogénesis con tasa de crecimiento bacterial <math display="inline">\mu _1(S_1)</math>
* <math display="inline">r_{2} = \mu _2(S_2) X_2</math> es la tasa de reacción en la metanogénesis con tasa de crecimiento bacterial <math display="inline">\mu _2(S_2)</math>.

Las variables en nuestra dinámica son:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>S_{1}:=  \hbox{ Cantidad de sustrato organico}  [g/l]</math>
|-
| style="text-align: center;" | <math> X_{1}:=  \hbox{ Cantidad de biomasa acetogenica}  [g/l]</math>
|-
| style="text-align: center;" | <math> S_{2}:=  \hbox{ Cantidad de sustrato acetogenico}  [g/l]</math>
|-
| style="text-align: center;" | <math> X_{2}:=  \hbox{ Cantidad de biomasa metanogenica}  [g/l]</math>
|-
| style="text-align: center;" | <math> A:=  \hbox{ Cantidad de alcalinidad}  [mmol/l]</math>
|-
| style="text-align: center;" | <math> C:=  \hbox{ Cantidad total de carbono inorganico ( }CO_2\hbox{)}  [mmol/l] </math>
|}
|}

y los parámetros son:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\mu _{1_{max}}:= \hbox{ Velocidad especifica de crecimiento maximo, asociada a }S_1</math>
|-
| style="text-align: center;" | <math> K_{S_1}:= \hbox{ Constante de saturacion media, asociada a }S_1</math>
|-
| style="text-align: center;" | <math> \mu _{2_{max}}:= \hbox{ Velocidad especifica de crecimiento maximo, asociada a }S_2</math>
|-
| style="text-align: center;" | <math> K_{S_2}:= \hbox{ Constante de saturacion media, asociada a }S_2</math>
|-
| style="text-align: center;" | <math> K_{I2}:= \hbox{ Constante de inhibicion, asociada a }S_2</math>
|-
| style="text-align: center;" | <math> k_{1}:=  \hbox{ Coeficiente cinetico asociado al sustrato }S_1</math>
|-
| style="text-align: center;" | <math> k_{2}:=  \hbox{ Coeficiente cinetico asociado al sustrato }S_2</math>
|-
| style="text-align: center;" | <math> k_{3}:=  \hbox{ Coeficiente cinetico asociado al consumo de sustrato }S_2</math>
|-
| style="text-align: center;" | <math> k_{4}:=  \hbox{ Coeficiente cinetico asociado a la produccion de }CO_2\hbox{ por }X_1</math>
|-
| style="text-align: center;" | <math> k_{5}:=  \hbox{ Coeficiente cinetico asociado a la produccion de }CO_2\hbox{ por }X_2</math>
|-
| style="text-align: center;" | <math> K_{L_{a}}:=  \hbox{ Es el coeficiente volumetrico de transferencia de oxigeno}</math>
|-
| style="text-align: center;" | <math> K_{H}:=  \hbox{ Constante de Henry}</math>
|-
| style="text-align: center;" | <math> P_{C}:=  \hbox{ Presion parcial del gas }C </math>
|}
|}



La modelo representa sólo dos pasos de la digestión anaerobia: acidogénesis y metanogénesis; el cual está dado por los siguientes balances ( ver  <span id='citeF-3'></span>[[#cite-3|[3]]] )

<span id="eq-1.a"></span>
<span id="eq-1.b"></span>
<span id="eq-1.c"></span>
<span id="eq-1.d"></span>
<span id="eq-1.e"></span>
<span id="eq-1.f"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\frac{dX_{1}}{dt}  = \left(\mu _{1_{max}} \frac{S_1}{S_1 + K_{S_1}} - \alpha D \right) X_{1}\triangleq F_{1}</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.a)
|-
| style="text-align: center;" | <math> \frac{dS_{1}}{dt}  = D \left(S_{1in} - S_1 \right)- k_1  \mu _{1_{max}} \frac{S_1}{S_1 + K_{S_1}}  X_{1}\triangleq F_{2}</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.b)
|-
| style="text-align: center;" | <math> \frac{dX_{2}}{dt} = \left(\mu _{2_{max}} \frac{S_2}{\frac{S_2^2}{K_{I2}} + S_2 + K_{S_2}} - \alpha D \right)X_{2}\triangleq F_{3} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.c)
|-
| style="text-align: center;" | <math> \frac{dS_{2}}{dt} = D \left(S_{2in} - S_2 \right)+ k_2 \left(\frac{ \mu _{1_{max}} S_1}{S_1 + K_{S_1}} \right)X_{1} - k_3 \left(\frac{\mu _{2_{max}} S_2}{\frac{S_2^2}{K_{I2}} + S_2 + K_{S_2}} \right)X_{2}\triangleq F_{4} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.d)
|-
| style="text-align: center;" | <math> \frac{dA}{dt} = D \left(A_{in} - A \right)\triangleq F_{5}</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.e)
|-
| style="text-align: center;" | <math> \frac{dC}{dt}  = D \left(C_{in} - C \right)+  \left(\frac{k_4 \mu _{1_{max}} S_1 X_{1}}{S_1 + K_{S_1}} \right)</math>
|-
| style="text-align: center;" | <math>  \; \; + \; \left(\frac{ k_5 \mu _{2_{max}} S_2 X_{2}}{\frac{S_2^2}{K_{I2}} + S_2 + K_{S_2}} \right)- K_{L_a}[C + S_2 - A - K_{H}P_{C}]\triangleq F_{6} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.f)
|}
|}

con condiciones iniciales

<span id="eq-1.g"></span>
<span id="eq-1.h"></span>
<span id="eq-1.i"></span>
<span id="eq-1.j"></span>
<span id="eq-1.k"></span>
<span id="eq-1.l"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>X_{1}(0) \; \; = \; \; c_{1},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.g)
|-
| style="text-align: center;" | <math> S_{1}(0) \; \; = \; \; c_{2},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.h)
|-
| style="text-align: center;" | <math> X_{2}(0) \; \; = \; \; c_{3},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.i)
|-
| style="text-align: center;" | <math> S_{2}(0) \; \; = \; \; c_{4},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.j)
|-
| style="text-align: center;" | <math> A(0) \;\; = \; \; c_{5},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.k)
|-
| style="text-align: center;" | <math> C(0) \; \; = \; \; c_{6}, </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1.l)
|}
|}

donde <math display="inline"> t\in [0,T]=\Omega ,</math> con <math display="inline">T>0 </math> y las partes derechas de las ecuaciones (<math display="inline">F_{i}; i=1,...6,</math>) son funciones clase <math display="inline">C^{0}(\Omega )</math>, es decir funciones contínuas en <math display="inline">\Omega </math>. Con <math display="inline">\mu _{1max}>0</math>, <math display="inline">\mu _{2max}>0</math>, <math display="inline">K_{S_{1}}>0</math>, <math display="inline">K_{S_{2}}>0</math>, <math display="inline">K_{I2}>0</math>, <math display="inline">A_{in}>0</math>, <math display="inline">C_{in}>0</math>, <math display="inline">K_{L_{a}}>0</math>, <math display="inline">S_{1in}>0</math>, <math display="inline">S_{2in}>0</math>, y <math display="inline">k_{i}>0</math>, <math display="inline">i=1,2,3,4,5.</math> constantes positivas. El parámetro <math display="inline">\alpha </math> (<math display="inline">0\leq \alpha \leq 1</math>) refleja la heterogeneidad de este proceso: <math display="inline">\alpha =0</math> corresponde a un reactor ideal de lecho fijo, mientras que <math display="inline">\alpha =1</math> corresponde a un reactor ideal de tanque agitado continuamente(CSTR) <span id='citeF-3'></span>[[#cite-3|[3]]].

La expresión <math display="inline">K_{L}a (C-K_H  P_C)</math>  describe la tasa de flujo molar del carbono inorgánico desde su fase líquida hasta su fase gaseosa y el producto <math display="inline">K_H  P_C = B>0</math> determina la concentración de oxígeno disuelto en <math display="inline">C</math>.

Éste, es un sistema de ecuaciones diferenciales ordinarias no lineales de primer orden y autónomo. Dadas estas características se tiene que toda solución de una ecuación diferencial autónoma o es constante o estrictamente creciente o estrictamente decreciente en todo su intervalo de definición. Otra característica importante es que este tipo de sistemas no posee soluciones oscilatorias.

El sistema [[#1.1 El modelo|1.1]], se ha descompuesto en subsistemas de los cuales el primero conformado por las variables <math display="inline">(X_{1},S_{1})</math> tiene soluciones analíticas semitriviales presentadas en <span id='citeF-15'></span>[[#cite-15|[15]]]. Para calcular soluciones analíticas no triviales se encuentra una relación directa entre la variable <math display="inline">S{1}</math> y la ecuación de Abel de primer tipo y aplicando el método propuesto en <span id='citeF-17'></span>[[#cite-17|[17]]] se logra encontrar soluciones exactas no triviales, para mayores detalles ver <span id='citeF-1'></span>[[#cite-1|[1]]].

===1.2 Sobre el análisis adimensional===

Por modelado matemático se entiende como el proceso por el cual se formula un modelo matemático que representa una realidad, donde se incluyen hipótesis, cantidades importantes y variables que intervienen en el fenómeno. Al resolver una ecuación, que represente un modelo, por alg'un método se hacen las simulaciones o aproximaciones pertinentes por medio de alg'un software, cuyos resultados se comparan con datos reales (validar) y de acuerdo a sus diferencias se puede decir si el modelo es predictivo para otras observaciones similares. En general, el objetivo primordial del modelado matemático es darle sentido al mundo tal y como lo observamos <span id='citeF-12'></span>[[#cite-12|[12]]].

Entonces, como consecuencia de la importancia del modelado es necesario desarrollarlo e interpretarlo por medio de técnicas que permitan esto. Dentro de estas se destacan el análisis dimensional y el reescalamiento. La primera permite comprender la dimensión de las cantidades en las ecuaciones y las implicaciones resultantes de la homogeneidad dimensional y la segunda es una técnica que ayuda a entender la magnitud de los términos que aparecen en el modelo comparando las cantidades con las cantidades de referencia que aparecen naturalmente en el fenómeno físico<span id='citeF-12'></span>[[#cite-12|[12]]].

Las magnitudes se clasifican según su origen en magnitudes fundamentales y derivadas. Una magnitud fundamental es aquella que se define por sí misma y es independiente de las demás. Estas están elegidas adecuadamente y son las que permiten expresar cualquier magnitud física en términos de ellas. Entre ellas se encuentran las contenidas dentro del Sistema Internacional de Unidades (SI) establecidas en la XI Conferencia General de Pesos y Medidas (1960) a las que más adelante en la XIV Conferencia General sobre Pesas y Magnitudes en 1971, se les agregó la unidad más reciente, el mol que representa la cantidad de sustancia (1971). Fourier en su libro “Théorie analytique de la chaleur” (1822), define que “Es necesario hacer notar que cada magnitud, indeterminada o constante, tiene una dimensión que le es propia, y que los términos de una no podrán ser comparados si no tuviesen los mismos exponentes de dimensiones”, lo que nos remite al principio de homogeneidad que debe tener cualquier ecuación. A diferencia de la anterior, una magnitud derivada es toda aquella que se puede formar a partir de expresiones matemáticas que relacionen magnitud {<math display="inline">F_{1},\ldots,F_{m}</math>, (<math display="inline">m<n</math>) de menor cantidad de elementos que permite derivar todas las magnitudes involucradas en un fenómeno específico. Recordemos las 7 magnitudes fundamentales definidas por el Bureau International des Points et Mesures, estas son: longitud (L), masa (M), tiempo (<math display="inline">T</math>), intensidad de corriente eléctrica (I), cantidad de sustancia (N), temperatura termodinánica (T) e intensidad luminosa (K).

Magnitudes adimensionales: En general, si

<span id="eq-2"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>[X_{j}]=F_{1}^{a_{1j}}F_{2}^{a_{2j}}...F_{m}^{a_{mj}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

Con <math display="inline">j=1,\ldots,n</math>. Si <math display="inline">[X_{j}]=1</math>. Entonces se dice que <math display="inline">X_{j}</math> es adimensional. Por lo tanto, las magnitudes adimensionales son un conjunto de <math display="inline">n-</math>magnitudes variables o constantes, donde cada una está definida por <math display="inline">m-</math>magnitudes fundamentales, que generan un nuevo sistema de ecuaciones.

Matriz de dimensión: Dado un conjunto de magnitudes adimensionales <math display="inline">[\mu _{1_{max}},\ldots ,M_n]=[M_j ]_{j=1}^{n}</math> representado por el sistema de magnitudes fundamentales <math display="inline">F_{1},\ldots,F_{m}</math> con (<math display="inline">m<n</math>), entonces las magnitudes adimensionales tienen la forma descrita en la ecuación [[#eq-2|2]]. Entonces se denomina ''matriz de dimensión'' a la matriz que tiene la forma:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> {\mathbf A}=   \begin{bmatrix} a_{11} & a_{12} & ... & a_{1n} \\ a_{21} & a_{22} & ... & a_{2n} \\ \vdots & \ddots & \ddots & \vdots \\ a_{m1} & a_{m2} & ... & a_{mn}   \end{bmatrix} </math>
|}
|}

Ley libre de unidades: Se define que, dada una ley <math display="inline">f(q_1,q_2,\ldots q_n)=0</math> se dice ''libre de unidades'' si para todo real positivo <math display="inline">\lambda _1,\ldots,\lambda _m</math> se cumple que: <math display="inline">\overline{f}(\overline{q}_1,\overline{q}_2,\ldots \overline{q}_n)=0 \Longleftrightarrow f(q_1,q_2,\ldots q_n)=0</math>, con <math display="inline">\overline{q}_1=\lambda _{1}^{b_1}...\lambda _{m}^{b_m} q_{j}</math>.

Ahora debemos mencionar lo que se considera la escencia de la adimensionalización y es el conocido Teorema de Pi-Buckingham. Aunque enunciado por primera vez por Aimé Vaschy en su obra “Sur les lois de similitude en physique” (1982), no fue publicado sino hasta 1914 por Edgar Buckingham. El teorema de Pi es conocido como ''la piedra angular'' del análisis dimensional, pues este indica que si hay una ley física que da una relación entre cierto n'umero de cantidades físicas dimensionadas, entonces hay una ley equivalente que puede expresarse como una relación entre ciertas cantidades adimensionales independientes que pueden formarse a partir de las variables fundamentales, denotadas como <math display="inline">\pi _1, \pi _2, \ldots, \pi _n</math>, de dónde se le ha asignado su nombre <span id='citeF-12'></span>[[#cite-12|[12]]].

<span id='theorem-ley'></span>Teorema 1: [Pi-Buckingham]

Sea

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>f(q_1,q_2,\ldots q_n)=0 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

una ley física libre de unidades que relaciona magnitudes dimensionales <math display="inline">q_1,q_2,\ldots q_n</math>. Y sean <math display="inline">[F_1,\ldots ,F_m]</math> con (<math display="inline">m<n</math>) magnitudes fundamentales con <math display="inline">[q_j]=F_{1}^{a_{1j}}F_{2}^{a_{2j}}...F_{m}^{a_{mj}}</math> con <math display="inline">j=1,\ldots ,n</math> y sea <math display="inline">r</math> el rango de la matriz <math display="inline">\mathbf A</math>, siendo <math display="inline">\mathbf A</math> la matriz de dimensión de <math display="inline">[q_j]</math>. Entonces existen <math display="inline">n-r</math> magnitudes adimensionales independientes <math display="inline">\pi _1, \pi _2, \ldots , \pi _{n-r}</math>, que se pueden formar a partir de <math display="inline">q_1,q_2,\ldots q_n</math>, y además la ecuación <math display="inline">f(\pi _1, \pi _2, \ldots , \pi _{n-r})=0</math>, que está expresada sólo en términos de las cantidades adimensionales, es equivalente a la ley física [[#theorem-ley|1]].  Su prueba se encuentra en <span id='citeF-12'></span>[[#cite-12|[12]]].

==2 Resultados principales==

Se toma la ecuación <span style="text-align: center; font-size: 75%;">

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> f\left(t; S_1, X_1, \mu _{1_{max}}, K_{S_1}, D, S_{1in}, S_2, X_2, \mu _{2_{max}}, K_{S_2}, K_{I2}, S_{2in}, k_2, k_3, A,A_{in},C,C_{in},k_4,k_5,K_{L_a},B\right)= 0 </math>
|}
|}

</span> Para probar que el sistema se puede adimensionalizar tomamos las unidades de cada uno de los parámetros y variables, y comprobamos que estén bien definidos en el sistema. Se revisan las unidades

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> S_1 = g/l \; ; \; X_1 = g/l \; ; \; \mu _{1_{max}} =\frac{1}{dia} \; ; \; K_{S_1} = g/l \; ; \; D= \frac{1}{dia} \; ; \; S_{1in} = g/l \; ; \; S_2 = \frac{m \;mol}{l} \; ; </math>
|}
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> X_2 = g/l \; ; \; \mu _{2_{max}} = \frac{1}{dia} \; ; \; K_{I2} =\frac{m \;mol}{l} \; ; \; K_{S_2} = \frac{m \;mol}{l} \; ; \; S_{2in}= \frac{m \;mol}{l} \; ; k_2 = \frac{m \;mol}{g} \; ; </math>
|}
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> k_3 = \frac{m \;mol}{g} \; ; \; A= \frac{m\;mol}{l}  \; ; \; A_{in}= \frac{m\;mol}{l}  \; ; \; C= \frac{m\;mol}{l}  \; ; \; C_{in}= \frac{m\;mol}{l}  \; ; \; k_4 = \frac{m \;mol}{g} </math>
|}
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> k_5= \frac{m\;mol}{g}  \; ; \; K_{L_a}= \frac{1}{dia}  \; ; \; B = \frac{m\;mol}{l} \; ; \; t = dia </math>
|}
|}

Se sustituyen las unidades en el sistema [[#1.1 El modelo|1.1]]. De aquí deteminamos las magnitudes fundamentales: <math display="inline">T</math>: tiempo,&nbsp;<math display="inline">M</math>: masa, &nbsp;<math display="inline">N</math>: cantidad de sustancia &nbsp;y la magnitud derivada:  &nbsp;<math display="inline">V</math>: volúmen que intervienen en el proceso. El paso siguiente es la contrucción de la matriz de dimensiones <math display="inline">\mathbf{A}</math> con <math display="inline">Rank\left(\mathbf{A} \right)= 3</math> y se determina que: el total de magnitudes (parámetros y variables) en el modelo son 25, dos de los cuales no tienen dimensiones (<math display="inline">\alpha </math> y <math display="inline">k_1</math>). Entonces se tiene que hay 23 magnitudes dimensionales. Ahora se calcula el número de magnitudes adimensiones que quedarán en el modelo adimensional, esto es <math display="inline">n-Rank\left(\mathbf{A} \right)=23-3=20</math>, de donde se tiene que el modelo adimensional es de 20.

Por el Teorema de Pi se tiene que la ecuación

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> F \left(\pi _1, \pi _2, \pi _3, \pi _4,  \pi _5, \pi _6, \pi _7, \pi _8, \pi _9, \pi _{10}, \pi _{11}, \pi _{12}, \pi _{13}, \pi _{14}, \pi _{15}, \pi _{16}, \pi _{17}, \pi _{18}, \pi _{19}, \pi _{20} \right)= 0  </math>
|}
|}

es equivalente a la ley  <span style="text-align: center; font-size: 75%;">

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> f\left(t; S_1, X_1, \mu _{1_{max}}, K_{S_1}, D, S_{1in}, S_2, X_2, \mu _{2_{max}}, K_{S_2}, K_{I2}, S_{2in}, k_2, k_3, A,A_{in},C,C_{in},k_4,k_5,K_{L_a},B\right)= 0. </math>
|}
|}

</span>

===2.1  Adimensionalización del modelo===

Se inicia con el reescalamiento

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\left[\pi _1\right]= S_{1in}^{-1} S_1 \quad ; \quad \left[\pi _2\right]= S_{1in}^{-1} X_1 \quad ; \quad \left[\pi _3\right]= t\;\mu _{1_{max}} \quad ; \quad \left[\pi _4\right]= S_{1in}^{-1} K_{S_1}</math>
|-
| style="text-align: center;" | <math> \left[\pi _5\right]= t \; D \quad ; \quad \left[\pi _6\right]= S_{2in}^{-1} S_2 \quad ; \quad \left[\pi _7\right]= S_{1in}^{-1} X_2 \quad ; \quad \left[\pi _8\right]= t\;\mu _{2_{max}}</math>
|-
| style="text-align: center;" | <math> \left[\pi _9\right]= S_{2in}^{-1} K_{S_2} \quad ; \quad \left[\pi _{10}\right]= S_{2in}^{-1} K_{I2} \quad ; \quad \left[\pi _{11}\right]= S_{1in} S_{2in}^{-1} k_2 \quad ; \quad \left[\pi _{12}\right]= S_{1in} S_{2in}^{-1} k_3</math>
|-
| style="text-align: center;" | <math> \left[\pi _{13}\right]= S_{2in}^{-1} A \quad ; \quad \left[\pi _{14}\right]= S_{2in}^{-1} A_{in} \quad ; \quad \left[\pi _{15}\right]= S_{2in}^{-1} C \quad ; \quad \left[\pi _{16}\right]= S_{2in}^{-1} C_{in}</math>
|-
| style="text-align: center;" | <math> \left[\pi _{17}\right]= S_{1in} S_{2in}^{-1} k_{4} \quad ; \quad \left[\pi _{18}\right]= S_{1in} S_{2in}^{-1} k_{5} \quad ; \quad \left[\pi _{19}\right]= t\;K_{L_a} \quad ; \quad  \left[\pi _{20}\right]= S_{2in}^{-1} B </math>
|}
|}

luego,

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\overline{t} = \pi _5 = t\; D \rightarrow t = \frac{1}{D} \overline{t} \rightarrow  \frac{dt}{d\overline{t}} = \frac{1}{D}</math>
|-
| style="text-align: center;" | <math> \overline{S}_1 = {\pi _1} = \frac{S_1}{S_{1in}} \rightarrow  S_1 = S_{1in} \overline{S}_1 </math>
|-
| style="text-align: center;" | <math> \overline{X}_1 = {\pi _2} =  \frac{X_1}{S_{1in}} \rightarrow  X_1 = S_{1in} \overline{X}_1</math>
|-
| style="text-align: center;" | <math> \vdots \;\;\; \;\;\;  \;\;\; \ldots \;\;\;  \;\;\; \;\;\; \vdots  </math>
|}
|}

y así sucesivamente con todas las magnitudes, de donde se obtiene el sistema adimensional <span style="text-align: center; font-size: 75%;">

<span id="eq-4.a"></span>
<span id="eq-4.b"></span>
<span id="eq-4.c"></span>
<span id="eq-4.d"></span>
<span id="eq-4.e"></span>
<span id="eq-4.f"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{l} \frac{d \overline{X}_1}{d \overline{t}} &= \left(\frac{\overline{\mu }_{1_{max}}\overline{S}_1}{\overline{S}_1 + \overline{K}_{S_1} } - \alpha \right)\overline{X}_1 \\ \frac{d \overline{S}_1}{d \overline{t}} &= 1 - \overline{S}_1 - k_1 \frac{\overline{\mu }_{1_{max}}\overline{S}_1}{\overline{S}_1 + \overline{K}_{S_1} } \overline{X}_1 \\ \frac{d \overline{X}_2}{d \overline{t}} &= \left(\frac{\overline{\mu }_{2_{max}} \overline{K}_{I2}\overline{S}_2 }{ \overline{S}_2^2 + \overline{K}_{I2}\;\overline{S}_2 + \overline{K}_{I2}\;\overline{K}_{S_2}   } - \alpha \right) \\ \frac{d \overline{S}_2}{d \overline{t}} &= 1 - \overline{S}_2 + \overline{k}_2 \; \frac{\overline{\mu }_{1_{max}}\overline{S}_1}{\overline{S}_1 + \overline{K}_{S_1} } \overline{X}_1 - \overline{k}_3 \frac{\overline{\mu }_{2_{max}} \overline{K}_{I2}\overline{S}_2 }{ \overline{S}_2^2 + \overline{K}_{I2}\;\overline{S}_2 + \overline{K}_{I2}\;\overline{K}_{S_2}   }\overline{X}_2 \\ \frac{d \overline{A} }{d \overline{t} } &= \; \overline{A}_{in} - \overline{A} \\ \frac{d \overline{C}}{d \overline{t}} & = \overline{C}_{in} - \overline{C} +  \overline{k}_4 \; \frac{\overline{\mu }_{1_{max}}\overline{S}_1}{\overline{S}_1 + \overline{K}_{S_1} } \overline{X}_1\\ & \; \; + \; \overline{k}_5 \frac{\overline{\mu }_{2_{max}} \overline{K}_{I2}\overline{S}_2 }{ \overline{S}_2^2 + \overline{K}_{I2}\;\overline{S}_2 + \overline{K}_{I2}\;\overline{K}_{S_2}} \overline{X}_2 - K_{L_a}[C + S_2 - A - B] \end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4.f)
|}

</span> con condiciones iniciales

<span style="text-align: center; font-size: 75%;">

<span id="eq-5.a"></span>
<span id="eq-5.b"></span>
<span id="eq-5.c"></span>
<span id="eq-5.d"></span>
<span id="eq-5.e"></span>
<span id="eq-5.f"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> \overline{X}_{1}(0) \; \; = \; \; \overline{c}_{1},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.a)
|-
| style="text-align: center;" | <math> \overline{S}_{1}(0) \; \; = \; \; \overline{c}_{2},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.b)
|-
| style="text-align: center;" | <math> \overline{X}_{2}(0) \; \; = \; \; \overline{c}_{3},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.c)
|-
| style="text-align: center;" | <math> \overline{S}_{2}(0) \; \; = \; \; \overline{c}_{4},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.d)
|-
| style="text-align: center;" | <math> \overline{A}(0) \;\; = \; \; \overline{c}_{5},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.e)
|-
| style="text-align: center;" | <math> \overline{C}(0) \; \; = \; \; \overline{c}_{6}, </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5.f)
|}
|}

</span>

===2.2  Aproximaciones numéricas===

La solución numérica de una Ecuación Diferencial Ordinara (EDO) se utiliza cuando no es posible obtener una solución analítica mediante integración de la función <math display="inline">f(x,y)</math> y se obtiene mediante algunos métodos como Runge-Kutta, métodos multipaso y los métodos de extrapolación. No sobra mencionar que cuando se halla la solución numérica, no se obtiene una expresión analítica sino un conjunto de valores numéricos.

Los métodos numéricos para resolver EDO se pueden agrupar en dos conjuntos:

* Métodos de un paso: Euler (simple) y Euler modificado, Taylor y Runge-Kutta (R-K).
* Métodos multipaso: Adams-Bashfort y Adams-Moulton.

Para determinar qué método numérico aplicar a un sistema de EDO, lo primero es determinar si dicho sistema es Stiff (rígido) o No.

La concepción de un sistema Stiff aún no tiene una definición como tal en la solución numérica de EDOs. En la literatura se encuentran dos aproximaciones <span id='citeF-4'></span>[[#cite-4|[4]]]:

* Un sistema de EDOs es Stiff cuando los métodos Backward Differentiation Formulae, normalmente funcionan mejor que los métodos explícitos.
* Un sistema de EDOs es Stiff cuando para algunos métodos se necesita una disminución del tamaño de paso para evitar inestabilidad del sistema.

Una manera practica para identificar(presunción) si un sistema de EDOs es Stiff es apoyarnos en la estabilidad de sus equilibrios, esto es (i) si un equilibrio es estable y todos sus autovalores asociados tienen parte real distinta, se podría asumir que el sistema es Stiff. (ii) Si el equilibrio es inestable o algún par de sus autovalores tienen parte real igual, se prodría asumir que el sistema es No-Stiff. 

Dicho esto, para sistemas Stiff se recomienda aplicar métodos implícitos y para los no-Stiff los métodos explícitos ver <span id='citeF-4'></span>[[#cite-4|[4]]].

Tomando el sistema original [[#1.1 El modelo|1.1]], lo primero es determinar si el sistema es estable, y para ello linealizamos el sistema en uno de los equilibrios más representativos para el proceso (punto de operación). Después de varias pruebas con distintos valores se observa:

<ol>

<li>En todos los casos se presentan uno o dos autovalores con parte real positiva, esto nos dice que el sistema en el equilibrio es inestable.   </li>
<li>De igual forma en todos los casos se presenta un autovalor con parte real negativa pero de multiplicidad algebraica 3, esto es, que aparece tres veces. </li>

</ol>

Con estos resultados se puede asumir que el sistema es No-Stiff (no rígido) y por esto se utiliza un Runge-Kutta de orden 2-3 o de orden 4-5. De igual manera después de varias pruebas se obseva que para este modelo AMOCO [[#1.1 El modelo|1.1]] se obtienen mejores salidas con un R-K de 4o orden como se puede ver en la Figura [[#img-1|1]] y la Figura [[#img-2|2]]:

<div id='img-1'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig1.png|337px|Graficas de X₁ y S₁ comparando los RK]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 1:''' Graficas de <math>X_1</math> y <math>S_1</math> comparando los RK
|}

<div id='img-2'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig2.png|337px|Graficas de X₂, S₂, A y C comparando los RK]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 2:''' Graficas de <math>X_2</math>, <math>S_2</math>, <math>A</math> y <math>C</math> comparando los RK
|}

A simple vista parece que no se encontraran mayores diferencias entre los RK2 y RK4. Pero haciendo una ampliación en una de ellas, en particular en la grafica de <math display="inline">C</math>, como se muestra en la Figura [[#img-3|3]] se puede detallar que el trazo del RK4 en más contíuo que el del RK2.

<div id='img-3'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig3.png|337px|Ampliación de la gráfica de C]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 3:''' Ampliación de la gráfica de <math>C</math>
|}

0.85 
{|  class="floating_tableSCP wikitable" style="text-align: left; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-1'></span>Table. 1 Los parámetros para las gráficas fueron tomados de <span id='citeF-3'></span>[[#cite-3|[3]]].
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" |   '''Parámetro''' 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  '''Valor''' 
| style="border-left: 2px solid;border-right: 2px solid;" |  '''Unidades''' 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  '''Desviación estandar'''
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>S_{1}(0)</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  2 
| style="border-left: 2px solid;border-right: 2px solid;" |  [g/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  2.5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>D</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  0.395 
| style="border-left: 2px solid;border-right: 2px solid;" |  [d<math>^{-1}</math>] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  0.135 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>S_{1in}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  10 
| style="border-left: 2px solid;border-right: 2px solid;" |  [g/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>K_{S_{1}}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  7.1 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>S_{2}(0)</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  10 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol /l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  10 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>S_{2in}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  87.5 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol /l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  12.5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>K_{I2}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  256 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  320 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>K_{S_{2}}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  9.28 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  13.7 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>A(0)</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  87.5 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  12.5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>A_{in}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  87.5 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  12.5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>C(0)</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  87.5 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  12.5 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>C_{in}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  9 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  1 
|- style="border-top: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>K_{L_a}</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  19.8 
| style="border-left: 2px solid;border-right: 2px solid;" |  [d<math>^{-1}</math>] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  3.5 
|- style="border-top: 2px solid;border-bottom: 2px solid;"
| style="border-left: 2px solid;border-right: 2px solid;" | <math>B</math> 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |  1777.6 
| style="border-left: 2px solid;border-right: 2px solid;" |  [mmol/l ] 
| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" |   

|}

Ahora graficamos las soluciones numéricas del modelo adimensional [[#2.1 Adimensionalización del modelo|2.1]]

<div id='img-4'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig4.png|337px|Gráficas de ̅X₁ y ̅S₁ comparando los RK]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 4:''' Gráficas de <math>\overline{X}_1</math> y <math>\overline{S}_1</math> comparando los RK
|}

<div id='img-5'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig5.png|337px|Graficas de ̅X₂, ̅S₂, ̅A y ̅C comparando los RK]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 5:''' Graficas de <math>\overline{X}_2</math>, <math>\overline{S}_2</math>, <math>\overline{A}</math> y <math>\overline{C}</math> comparando los RK
|}

Aquí se presenta una situación similar a la del modelo original y es que en forma aparente no hay mayores diferencias entre los RK2 y RK4. Pero haciendo una ampliación en la gráfica de <math display="inline">\overline{S}_2</math> (Figura [[#img-6|6]]) se obserba, también, que el trazo del RK4 en más contíuo que el del RK2.

<div id='img-6'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig6.png|337px|Ampliación de la gráfica de ̅S₂]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 6:''' Ampliación de la gráfica de <math>\overline{S}_2</math>
|}

===2.3 Comparación de ambos modelos===

Si observamos las salidas de las graficas del RK4 (o RK2) del modelo inicial y lo comparamos con el RK4 (o RK2) del modelo adimensional, se observa que el modelo adimensional tiende a estabilizarse en menos tiempo que el modelo inicial como se puede ver en las Figuras [[#img-7|7]] y [[#img-8|8]].

<div id='img-7'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig7.png|337px|Comparación de las salidas de ambos modelos para X₁]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 7:''' Comparación de las salidas de ambos modelos para <math>X_{1}</math>
|}

<div id='img-8'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_Machado_497328595-fig8.png|337px|Comparación de las salidas de ambos modelos para S₂]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 8:''' Comparación de las salidas de ambos modelos para <math>S_{2}</math>
|}

==Conclusiones==

* Comparando las gráficas de los dos sistemas podemos observar que las salidas de las aproximaciones del sistema adimensional se estabilizan en menos tiempo que el sistema original, esto es un resultado relevante en el proceso dado que para este tipo de modelos al llevarlos a escala industrial el ahorro del tiempo minimiza costos y puede optimizar resultados.
* Por lo anterior se puede afirmar que además de reducir parámetros en el modelo AMOCO, la adimensionalización del sistema, permite conservar todas las propiedades, principios y leyes que están presentes en el proceso de Digestión Anaerobia para la producción de biogás.
* Se presentó una metodología práctica para el estudio de sistemas de ecuaciones diferenciales ordinarias no lineales que modela la Digestión Anaerobia y que puede ser extendida a otros bioprocesos representados por este tipo de ecuaciones.

===BIBLIOGRAFÍA===

<div id="cite-1"></div>
'''[[#citeF-1|[1]]]'''  P.B. Acosta-Humánez, M. Machado-Higuera and A.V. Sinitsyn, ''A Model Of Anaerobic Digestion For Biogas Production Using Abel Equations'', Far East Journal of Mathematical Sciences (FJMS), 101(6), 1295&#8211;1311, (2017). http://dx.doi.org/10.17654/MS101061295.

<div id="cite-2"></div>
'''[[#citeF-2|[2]]]'''   Abdelhamid Ajbar, Khalid Alhumaizi. (2012). ''Dynamics of the Chemostat: A Bifurcation Theory Approach''; CRC Press Taylor & Francis Group. http://www.crcpress.com/Dynamics-of-the-Chemostat-A-Bifurcation-Theory-Approach/Ajbar-Alhumaizi/p/book/9781138112780.

<div id="cite-3"></div>
'''[[#citeF-3|[3]]]'''  O. Bernard, Z. Hadj-Sadok, D. Dochain, A. Genovesi and J-P. Steyer, ''Dynamical model development and parameter identification for an anaerobic wastewater treatment process'', Biotechnology and Bioengineering, 75(4), 424&#8211;438, (2001). https://www.ncbi.nlm.nih.gov/pubmed/11668442.

<div id="cite-4"></div>
'''[[#citeF-4|[4]]]'''  C.F. Curtiss and J.O. Hirschfelder, ''Integration of Stiff equations'', Proceedings of the National Academy of Sciences, 38(3), 235&#8211;243, (1952). https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1063538/.

<div id="cite-5"></div>
'''[[#citeF-5|[5]]]'''  María Teresa Varnero Moreno. (2011). ''Manual de Biogás. Santiago de Chile: FAO''; Organización de las Naciones Unidas para la Alimentación y la Agricultura (FAO). http://www.fao.org/family-farming/detail/es/c/342734/.

<div id="cite-6"></div>
'''[[#citeF-6|[6]]]'''  G.M. Mejía Sánchez. (1996). ''Digestión anaerobia''; Universidad Autónoma de Yucatán (UADY).  http://www.libreria.uady.mx/viewlib.php?i=398.

<div id="cite-7"></div>
'''[[#citeF-7|[7]]]'''  Besel, S. A. (Departamento de Energía). (2007). ''Biomasa: Digestores anaerobios''; Instituto para la Diversificación y Ahorro de la Energía.  http://www.idae.es/publicaciones/biomasa-digestores-anaerobios.

<div id="cite-8"></div>
'''[[#citeF-8|[8]]]'''   D.J. Batstone, J. Keller, I. Angelidaki, S.V. Kalyuzhnyi, S.G. Pavlostathis, A. Rozzi, W.T.M. Sanders, H. Siegrist and V.A. Vavilin, ''The IWA Anaerobic Digestion Model No.1 (ADM1)'', Water Science and Technology, 45(10), 65&#8211;73, (2002).  http://wst.iwaponline.com/content/45/10/65.

<div id="cite-9"></div>
'''[[#citeF-9|[9]]]'''   B. Benyahia, T. Sari, B. Cherky, J. Harmand, ''Bifurcation and stability analysis of a two step model for monitoring anaerobic digestion processes'',  Elsevier Ltd., Journal of Process Control 22(6), 1008&#8211;1019, (2012). http://www.sciencedirect.com/science/article/pii/S0959152412000972.

<div id="cite-10"></div>
'''[[#citeF-10|[10]]]'''  J. Del Real Olvera and J. Islas Gutiérrez, ''Biodegradación anaerobia de las aguas generadas en el despulpado del café'', Revista Colombiana de Biotecnología, 12(2), 230&#8211;239, (2010).  https://revistas.unal.edu.co/index.php/biotecnologia/article/view/18773.

<div id="cite-11"></div>
'''[[#citeF-11|[11]]]'''  W.J. Jewell, ''Anaerobic Sewage Treatment. Part 6'',  Environmental Science and Technology, 21(1), 14&#8211;21, (1987).  http://pubs.acs.org/doi/abs/10.1021/es00155a002.

<div id="cite-12"></div>
'''[[#citeF-12|[12]]]'''  J. David Logan. (2006). ''Applied Mathematics''; Wiley Interscience. https://books.google.com.co/books?id=nUk_AQAAIAAJhl=essource=gbs_book_other_versions.

<div id="cite-13"></div>
'''[[#citeF-13|[13]]]'''  Y. Lorenzo Acosta y M.C. Obaya Abreu, ''La Digestión Anaerobia. Aspectos teóricos. Parte I'',  ICIDCA. Sobre los Derivados de la Caña de Azúcar, XXXIX(1), 35&#8211;48, (2005).  http://www.redalyc.org/articulo.oa?id=223120659006.

<div id="cite-14"></div>
'''[[#citeF-14|[14]]]'''  M. Hernández y L. Delgadillo, ''Aplicación del modelo Adm1 en la digestión anaerobia de aguas residuales y desechos sólidos'',  Revista Tumbaga, 1(6), 29&#8211;42, (2011).  http://revistas.ut.edu.co/index.php/tumbaga/article/view/46.

<div id="cite-15"></div>
'''[[#citeF-15|[15]]]'''  M. Machado-Higuera and A.V. Sinitsyn, ''Existence of lower and upper solutions in reverse order with respect to a variable in a model of acidogenesis to anaerobic digestion'',   Bulletin of the South Ural State University. Ser. Mathematical Modelling, Programming & Computer Software, 8(2), 55&#8211;68, (2015).  https://www.scilit.net/article/00cd41ee5077298166fe0d2ff968ca97.

<div id="cite-16"></div>
'''[[#citeF-16|[16]]]'''  F.&nbsp;Mariet, O.&nbsp;Bernard, M.&nbsp;Ras, L.&nbsp;Lardon and J.&nbsp;Steyer, ''Modeling anaerobic digestion of microalgae using ADM1'',  Bioresource Technology 102(13), 6823&#8211;6829, (2011).  doi:10.1016/j.biortech.2011.04.015.

<div id="cite-17"></div>
'''[[#citeF-17|[17]]]'''  E.&nbsp;Salinas, R.&nbsp;Muñoz, J.&nbsp;Sosa, B.&nbsp;López, ''Analysis to the solutions of Abel's Differential Equations of the first kind under transformation <math>y = u(x)z(x)+v(x)</math>'', Applied Mathematical Sciences 7(42), 2075&#8211;2092, (2013).  http://www.m-hikari.com/ams/ams-2013/ams-41-44-2013/index.html.

<div id="cite-18"></div>
'''[[#citeF-18|[18]]]'''  M. Sbarciog, M. Loccufier and E. Noldus, ''Determination of appropriate operating strategies for anaerobic digestion systems'', Elsevier Ltd., Biochemical Engineering Journal 51(3), 180&#8211;188, (2010).  https://www.sciencedirect.com/science/article/pii/S1369703X10001865

<div id="cite-19"></div>
'''[[#citeF-19|[19]]]'''  B.&nbsp;Sialve, N.&nbsp;Bernet and O.&nbsp;Bernard,  ''Anaerobic digestion of microalgae as a necessary step to make microalgal biodiesel sustainable'', Biotechnology Advances 27(4), 409&#8211;416, (2009).    doi:10.1016/j.biotechadv.2009.03.001.

<div id="cite-20"></div>
'''[[#citeF-20|[20]]]'''  V.M. Trejos, J. Fontalvo y M.A. Gómez, ''Descripción Matemática y Análisis de Estabilidad de Procesos Fermentativos'',  Dyna 76(158), 111&#8211;121, (2009). https://revistas.unal.edu.co/index.php/dyna/article/view/10251.

<div id="cite-21"></div>
'''[[#citeF-21|[21]]]'''  A.O. Wagner, C. Malin, P. Lins, G. Gstraunthaler and P. Illmer, ''Reactor performance of a 750 m<math>^3</math> anaerobic digestion plant: Varied substrate input conditions impacting methanogenic community'',  Anaerobe 29, 29&#8211;33, (2014). https://www.sciencedirect.com/science/article/pii/S1075996414000328.

<div id="cite-22"></div>
'''[[#citeF-22|[22]]]'''  M. Weedermann, G. Seo and G.S. Wolkowicz, ''Mathematical model of anaerobic digestion in a chemostat: effects of syntrophy and inhibition'',  Journal of Biological Dynamics 7(1), 59&#8211;85, (2013). [ ]http://www.tandfonline.com/doi/abs/10.1080/17513758.2012.755573.