Semana #8 - Introducción al análisis de sensibilidad

La entrada de esta semana se basará en el documento "An Introduction to Sensitivity Analysis" producido en MIT bajo supervision de Dr. Jay Forrester y escrito por Lucia Breierova y Mark Choudhari en 1996.

El objetivo de esta semana es conocer que es un análisis de sensibilidad y para que sirve, después de examinar la inconsistencia dimensional de la semana #7 es de gran importancia hacer análisis de sensibilidad de nuestros modelos y aumentar la consistencia y precisión de estos.

¿Que es un análisis de sensibilidad? ¿Para que sirve?

El análisis de sensibilidad se utiliza para determinar cómo de "sensible" es un modelo a los cambios en el valor de los parámetros del modelo y a los cambios en la estructura este. En este trabajo, nos centramos en la sensibilidad de parámetros, que se realiza generalmente como una serie de pruebas en las que el modelador establece diferentes valores de los parámetros para ver cómo un cambio en el parámetro provoca un cambio en el comportamiento dinámico de las poblaciones.Al mostrar cómo el comportamiento del modelo responde a los cambios en los valores de parámetros, el análisis de sensibilidad es una herramienta útil en la construcción de modelos, así como en la evaluación del mismo. 

El análisis de sensibilidad ayuda a construir la confianza en el modelo mediante el estudio de las incertidumbres que se asocian a menudo con los parámetros de los modelos. Muchos de los parámetros en los modelos de dinámica de sistemas representan cantidades que son muy difíciles o incluso imposibles de medir con mucha precisión en el mundo real. Además, algunos valores de los parámetros cambian en el mundo real.

Las pruebas de sensibilidad ayudan al modelador para entender la dinámica de un sistema. Experimentar con una amplia gama de valores puede ofrecer conocimientos sobre el comportamiento de un sistema en 

situaciones extremas.

Ejercicio de aplicación: ¿El infierno es exotérmico o endotérmico?

He preparado un ejemplo muy interesante basado en una historia difundida en Internet. Veamos el infierno desde un punto de vista físico:

Se puede decir que el infierno es un lugar al que van

almas condenadas, es decir, es un contenedor. Si asumimos que las almas están en estado gaseoso se

puede hacer un análisis de la situación del mismo usando herramientas como la Ley de Boyle-Mariotte (un gas se enfría cuando se expande y se calienta cuando se comprime) y la Ley de los gases ideales.

En primer lugar, necesitamos saber en qué medida la masa del Infierno varía con el tiempo.Para ello hemos de saber a qué ritmo entran las almas en el Infierno y a qué ritmo salen. Tengo sin embargo entendido que, una vez dentro del Infierno, las almas ya no salen de él, por lo tanto, no se producen salidas.

En lo referente a cuántas almas entran, veamos lo que dicen las diferentes religiones. La mayoría de ellas declaran que si no perteneces a ellas, irás al Infierno. Dado que hay más de una religión que así se expresa y dado que la gente no pertenece a más de una, podemos concluir que todas las almas van al Infierno. Con las tasas de nacimientos y muertes existentes, podemos deducir que el número de almas en el Infierno crece de forma exponencial. 

Veamos ahora cómo varía el volumen del Infierno. Según la Ley de Boyle, para que la temperatura y la presión del Infierno se mantengan estables, el volumen debe expandirse en proporción a la entrada de almas

Hay dos posibilidades:

1. Si el Infierno se expande a una velocidad menor que la de entrada de almas, la temperatura y la presión en el Infierno se incrementarán hasta que éste se desintegre.
2. Si el Infierno se expande a una velocidad mayor que la de la entrada de almas, la temperatura y la presión disminuirán hasta que el Infierno se congele.

Teniendo en cuenta lo anterior planteo el siguiente modelo:

Como puede observarse en el modelo escogí como la variable de nivel el volumen disponible en el infierno el cual se determina la por la expansión que se produce y la entrada de almas a ese espacio. A su vez la expansión depende de la presión, el flujo de almas y el volumen actual. El flujo de almas depende únicamente de la tasa de mortalidad.

Por otro lado, la presión es una de las variables importantes del modelo (de hecho este modelo puede ser representado con dos niveles) y depende de la temperatura, el volumen, los moles del gas y la constante universal , mediante la ecuación de gases ideales P = NRT/V.  Para conocer el valor de la temperatura solo basta despejar la "T" de la ecuación. Sin embargo para calcular la velocidad de expansión que al final será la expansión del infierno, se usa la ecuación de Boyle V1=(P2*V2)/P1 y V2=(P1*V1)/P2.

Juguemos un poco con los parámetros del modelo: Si aumentamos la velocidad de expansión de tal manera que sea superior a la tasa de mortalidad->flujo de almas entonces sucederá lo expuesto en la segunda posibilidad: el infierno se congelará. Afectará al volumen disponible de tal forma que la presión y la temperatura se decrecerán:

En el caso contrario, es decir si la velocidad de expansión es menor al flujo de almas, entonces la presión y la temperatura aumentaran hasta explotar (imagine inflar un globo hasta explotar).

En cualquiera de los dos casos el modelo se mantiene, ahora veamos que pasa si modificamos otras variables como el numero de moles del gas. No sabemos la composición de un alma humana por lo tanto no podemos determinar con exactitud las moles de este, sin embargo este es un valor fijo y lo único que cambiaría seria que aumentaría y disminuiría la presión y la temperatura según el caso. Si el valor es elevado, se necesitarán menos almas para llegar a la desintegración en caso de que la expansión sea menor a la entrada de almas, y si el valor es mínimo se necesitarán mas almas para que ocurra esto y lo mismo aplica en el caso de la congelación.

Podemos concluir que el modelo es sensible a los cambios y  se mantiene, sin embargo podemos hacer un modelo mas complejo si tomamos en cuenta otros factores como el clima dentro del infierno o las guerras en la tierra y nos tocaría hacer un nuevo análisis de sensibilidad.

Referencias:

• Análisis de sensibilidad - http://ocw.mit.edu/courses/sloan-school-of-management/15-988-system-dynamics-self-study-fall-1998-spring-1999/readings/sensitivityanalysis.pdf
• Is Hell exothermic or endothermic?.Support your answer with a proof. - http://www.pinetree.net/humor/thermodynamics.html