ANALISIS COMPARTIMENTAL

ANALISIS COMPARTIMENTAL Estudiaremos modelos básicos de análisis compartimental. Se trata de describir mediante una fun

Views 99 Downloads 2 File size 411KB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

ANALISIS COMPARTIMENTAL

Estudiaremos modelos básicos de análisis compartimental. Se trata de describir mediante una función x(t) la cantidad de una sustancia que esté presente en un compartimento en el instante de tiempo t. El “compartimento” puede ser de cualquier tipo: un lago, un tanque de mezclas, etc. La idea básica de estos modelos está en una ley de conservación evidente: la tasa de cambio de la sustancia en el compartimento

𝑑𝑥 𝑑𝑡

será

igual a la velocidad de entrada de la sustancia en el compartimento en el instante t menos la velocidad de salida de la misma:

𝑑𝑥 = 𝑣 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎) − 𝑣 (𝑠𝑎𝑙𝑖𝑑𝑎) 𝑑𝑡 Dado que las velocidades de entrada y salida de la sustancia en el compartimento dependen del proceso en cuestión, poco más podemos decir de manera general sobre estos modelos. Pasamos por tanto directamente a analizar un problema concreto:

EJEMPLO: El agua del Lago Magdalena se está viendo sujeta a un proceso contaminante debido a la concentración de plaguicidas consecuencia de la fumigación de los naranjales cercanos. Por otro lado, el río Aguadulce, que desemboca en el lago, fluye hacia este a razón de 200 l/m portando una concentración de plaguicidas de 5 partes por millón. Si se suspende la fumigación en los alrededores del lago en el momento en el que la concentración de plaguicidas había alcanzado el valor de 40 partes por millón y se supone que en dicho instante el volumen del lago es de 100 millones de litros, calcular el tiempo que transcurrirá hasta que la concentración sea inferior a 20 partes por millón.

¿Qué volumen tendrá el lago en ese instante? Nota: Suponer que el lago pierde agua a razón de 300 l/m.

SOLUCION: Denominaremos x(t) a la cantidad de plaguicidas presentes en el lago a tiempo t. La concentración de plaguicidas en el lago a tiempo t, que denominaremos c(t), será obviamente: c(t) = x(t) vol (t), donde vol (t) denota el volumen del lago en cada instante de tiempo. Es importante, en este tipo de problemas, distinguir con precisión las velocidades de entrada y salida de la disolución y las correspondientes al soluto, es decir a la sustancia en cuestión (en este caso los plaguicidas). Desde este punto de vista la velocidad de entrada de disolución en el compartimento (lago) es v(entrada) dis. = 200 l/m y la de salida v (salida) dis. = 300 l/m, donde l denota litro de disolución. Las respectivas velocidades del soluto se obtendrán multiplicando las de la disolución por la concentración correspondiente. De esta manera, al tener la disolución entrante una concentración de 5 partes por millón (es decir 5 mg/l) tendremos: v (entrada) = 200 l/m 5 mg/l = 1000 mg/m.

Por otro lado, el volumen vol(t) será evidentemente: vol(t) = vol (0) + (v(entrada) dis. – v(salida) dis.) t. Como: vol (0) = 100 ∗ 106 = 108 , se tiene:

𝑣𝑜𝑙(𝑡) = 108 − 100𝑡 = 100(106 − 𝑡) ≫ 𝑐(𝑡) =

𝑥(𝑡) 100(106 − 𝑡)

Podemos calcular finalmente la velocidad de salida:

𝑣 𝑠𝑎𝑙𝑖𝑑𝑎 = 300

1 𝑥(𝑡) 𝑚𝑔 3𝑥(𝑡) = 𝑚𝑔/𝑚 𝑚 100(106 − 𝑡) 1 106 − 𝑡

y así escribir la ecuación diferencial: 𝑑𝑥 3𝑥 = 1000 − 𝑑𝑡 106 − 𝑡 La ecuación es fácilmente identificable como una ecuación lineal de primer orden en las variables (t, x). La ecuación lineal homogénea asociada:

𝑑𝑥 3𝑥 + 6 = 0 ≫ 𝑥 = 𝑘(106 − 𝑡)3 𝑑𝑡 10 − 𝑡 y en consecuencia el cambio de variable adecuado para resolver la ecuación completa será 𝑥 = 𝑢(106 − 𝑡)3 . Este cambio convierte la ecuación en: 𝑑𝑢 1000 = 𝑑𝑡 (106 − 𝑡)3 cuya integración nos conduce a la solución general (tras deshacer el cambio de variable):

∫ 𝑑𝑢 = ∫

1000 (106 − 𝑡)

3

𝑑𝑡

𝑎 = (106 − 𝑡) 𝑑𝑎 = −𝑑𝑡 Reemplazando:

∫ 𝑑𝑢 = − ∫

𝑢=

1000 (𝑎)3

𝑑𝑎

1000 +𝐶 2(106 − 𝑡)2

Luego:

𝑥 = 𝑢(106 − 𝑡)3 ≫≫ 𝑢 =

𝑥 (106 −𝑡)3

𝑥(𝑡) = 500(106 − 𝑡) + 𝐶(106 − 𝑡)3

𝑥(𝑡) =

𝑥(𝑡) 𝑐 ≫≫ 𝑐(𝑡) = 5 + (106 − 𝑡)2 100(106 − 𝑡) 100

Para la cantidad de plaguicidas y para la concentración de plaguicidas en el lago, respectivamente. Si se tiene en cuenta ahora la condición inicial: c (0) = 40, es fácil calcular la constante C que determina la solución particular que estamos buscando: 𝐶 = 35 ∗ 10−10. Podemos escribir finalmente la solución particular: 𝑥(𝑡) = 500(106 − 𝑡) + 35 ∗ 10−10 (106 − 𝑡)3 𝑐(𝑡) = 5 + 35 ∗ 10−12 (106 − 𝑡)2

GRAFICA de X(t) (x en kilogramos, t en días) C(t) (x en mg/ml, t en días)