Las prácticas de hoy han tratado procesos de extracción líquido-líquido y filtración con membranas. También hemos visto cómo emplear MATLAB para obtener resultados similares a los conseguidos con Berkeley Madonna. Consulta la asistencia de hoy en este
enlace.
Planteando cómo modelar el caso de la extracción líquido-líquido para un proceso de N etapas me ha surgido una duda, ya que no veo claro cúal sería el número de variables y, por tanto, los grados de libertad del sistema.
ResponderEliminarEl modelo propuesto sería el siguiente:
d/dt(X[1..N])=(L*(X[i-1])-X[i]-Q[i])/VL[i]
d/dt(Y[1..N])=(G*(Y[i+1]-Y[i])+Q[i])/VG[i]
Q[1..N]=KLA[i]*(X[i]-Xstar[i])*(VL[i]+VG[i])
Xstar[1..N]=Y[i]/M
Se está suponiendo que los volúmenes en cada etapa son constantes, pero no iguales.
El número de ecuaciones es 4*N, pero no estoy segura del número de variables, ¿serían 7*N+4?
Otra duda que tengo, pero ya referida a la práctica de separación mediante membranas, es que no entiendo por qué para el caso particular en que Q=0, definimos P=0.1
Respecto a la práctica de separación por membranas se usaba P=0.1 para evitar el vaciado del tanque (sacamos menos permeado), otra opción era disminuir el stoptime.
ResponderEliminarhttp://rapidshare.com/files/113072988/BLOG.doc.html
ResponderEliminarSe me ha olvidado decir de que trata el último comentario. He tratado la extracción líquido-líquido en multietapas para el funcionamiento en contracorriente y en paralelo.
ResponderEliminarA la pregunta que se nos realizó en la práctica de extracción L-L (Calcular el caudal de disolvente extractor mínimo para rebajar la concentración de soluto en la corriente L a la mitad), nosotros operábamos en Madonna con Parameter Plot para realizar una gráfica de G=f(X1final).
ResponderEliminarUna segunda opción es realizar un sliders con G en la gráfica de X1 e ir viendo dónde nos queda la horizontal de la curva, y cuando esté a la altura de 0.5 (X0/2) leemos el "G" o caudal de extractor
¡Hola Cristina!
ResponderEliminarCreo que tu modelo no está bien planteado, porque tendrías que particularizar X[1] y Y[N] ya que estas variables tienen que definirse teniendo en cuenta la alimentación. Por otra parte, KLA yo creo que puede considerarse constante para un sistema G/L concreto.
Te dejo aquí como lo he planteado yo; aunque he supuesto que VG y LG son iguales en todas las etapas.
{ECUACIONES = 2*N + 2*(N-1) + 2}
X'[1]=(-Q[i]+L*(X0-X[i]))/VL
Y'[N] = (Q[i]+G*(Y0-Y[i]))/VG
X'[2..N]=(-Q[i]+L*(X[i-1]-X[i]))/VL
Y'[1..N-1]=(Q[i]+G*(Y[i+1]-Y[i]))/VG
Q[1..N]=KLA*(X[i]-Xstar[i])*(VL+VG)
Xstar[1..N]=Y[i]/M
{VARIABLES =2*N + 2*(N-1) + 8}
{Y[], X[], Xstar[], Q[], G, L, VL, VG, KLA, M, YN+1, X0}
{GRADOS DE LIBERTAD = 8}
G = 3
L = 1
Y0 = 0
X0 = 1
VG =10
VL = 10
KLA = 2.5
M =0.8
{CONDICIONES INICIALES = 2*N}
INIT Y[1..N] = 0
INIT X[1..N] = 1
A ver si te sirve... Creo que mi modelo está bien, porque lo he comparado con el que cólgó Juan Carlos poniendo 3 etapas y me salía lo mismo.
¡Gracias!
¡Se me ha olvidado decir que el modelo anterior es para una extracción de N etapas en contracorriente!
ResponderEliminar