viernes, 27 de junio de 2008

TUTORÍAS INDIVIDUALES

Aquellos que deseen consultar dudas de cara al próximo examen, pueden pasarse por el despacho nº 4 del Departamento de Ingeniería Química, de lunes a miércoles a partir de las 18:00.

64 comentarios:

  1. En la práctica Nº 11, para representar el plano de fases, ¿hay que realizar alguna modificación en los ejes, de escalas o algo?

    Un saludo

    ResponderEliminar
  2. Debes poner la siguiente orden: RENAME TIME=TR
    Así, en lugar de que la variable independiente sea T, ahora será TR, que para nosotros es la T del tanque; entonces al escribir
    STARTIME =250, estás indicando que la temperatura de inicio es 250 (porque ahora TIME es TR)
    STOPTIME= 500, la T final será 500.
    Los incrementos ahora son de T no de tiempo.

    No sé si preguntabas esto o te he entendido mal, de todas maneras ánimo!

    ResponderEliminar
  3. Hola! me gustaria saber si los ejercicios que hicimos con matlab y excell en clase podrian caer en el examen o solo lo hicimos a modo de ejemplo.

    ResponderEliminar
  4. Me parece que los ejercicios de matlab y excell son a modo de ejemplo, porque en la convocatoria del examen solo pone que se puede utilizar el Berkeley.

    ResponderEliminar
  5. Sí hay que modificar los ejes. Al localizar los estados estacionarios ves que están en 289 y 441K (estables). Por eso, en el diagrama de fases cambias la escala del eje X (TR) desde 280 hasta 500 (rango en el cual están comprendidos los estados estacionarios). El eje Y (CA) irá desde 0 hasta 2.5, puesto que este es el valor inicial y máximo de concentración.

    Espero que esto resuelva tu duda

    ResponderEliminar
  6. Hola. He planteado una variante a la práctica 5. Como he puesto un dibujillo del proceso el desarrollo se puede ver en el siguiente enlace: http://rapidshare.com/files/126153155/extracci_n_L-L.doc.html

    Si hay problemas comentadlo por aquí. Gracias

    ResponderEliminar
  7. Hola otra vez. He planteado una variación respecto a la práctica 3. El modelo es similar al de clase, he mantenido el reflujo variable y el destilado acumulado, pero lo he hecho con tres componentes, A, B y C, (C no aparece en los balances pues queda definido a partir de los otros dos). La única dificultad viene a la hora de fijar los grados de libertad y las condiciones iniciales.
    No me salen datos coherentes,creo que fallo en algo, jejeje. Por favor si alguien lo ve, que me lo diga. P.D: En las condicionales en lugar de -, es < .Gracias

    {DESTILACIÓN DISCONTÍNUA MEZCLA TERNARIA (MULTICOMPONENTE), 5 PISOS, VELOCIDAD DE VAPORIZACIÓN CONSTANTE, REFLUJO VARIABLE}
    {MÉTODO NUMÉRICO}

    METHOD RK4
    DT=0.01
    STOPTIME=1000

    {ECUACIONES=28}

    {CALDERÍN}
    MB'=L-V
    MBA'=L*XA4-V*YA5
    MBB'=L*XB4-V*YB5
    XA5=MBA/MB
    XB5=MBB/MB
    YA5=AL1*XA5/(AL1*XA5+AL2*XB5+AL3*(1-XA5-XB5))
    YB5=AL2*XB5/(AL1*XA5+AL2*XB5+AL3*(1-XA5-XB5))

    {PISO 4, M=CTE EN LOS PISOS 4, 3 2 1}
    XA4'=(L*(XA3-XA4)+V*(YA5-YA4))/M
    XB4'=(L*(XB3-XB4)+V*(YB5-YB4))/M
    YA4=AL1*XA4/(AL1*XA4+AL2*XB4+AL3*(1-XA4-XB4))
    YB4=AL2*XB4/(AL1*XA4+AL2*XB4+AL3*(1-XA4-XB4))

    {PISO 3}
    XA3'=(L*(XA2-XA3)+V*(YA4-YA3))/M
    XB3'=(L*(XB2-XB3)+V*(YB4-YB3))/M
    YA3=AL1*XA3/(AL1*XA3+AL2*XB3+AL3*(1-XA3-XB3))
    YB3=AL2*XB3/(AL1*XA3+AL2*XB3+AL3*(1-XA3-XB3))

    {PISO 2}
    XA2'=(L*(XA1-XA2)+V*(YA3-YA2))/M
    XB2'=(L*(XB1-XB2)+V*(YB3-YB2))/M
    YA2=AL1*XA2/(AL1*XA2+AL2*XB2+(1-XA2-XB2)*AL3)
    YB2=AL1*XB2/(AL1*XA2+AL2*XB2+(1-XA2-XB2)*AL3)

    {PISO 1, CONDENSADOR TOTAL, REFLUJO VARIABLE}
    L=V-D
    D=V/(R+1)
    R=IF TIME-TIME1 THEN 1E10 ELSE
    IF TIME-TIME2 THEN R1 ELSE
    IF TIME-TIME3 THEN 1E10 ELSE R1
    XA1'=(V*YA2-(L+D)*XA2)/M
    XB1'=(V*YB2-(L+D)*XB2)/M

    {DESTILADO ACUMULADO}
    MD'=D
    MDA'=D*XA1
    MDB'=D*XB1
    XAD=MDA/MD
    XBD=MDB/MD

    {VARIABLES=37}

    {GRADOS DE LIBERTAD=9}
    AL1=7
    AL2=5
    AL3=1
    V=12
    M=20
    R1=8
    TIME1=250
    TIME2=500
    TIME3=900

    {CONDICIONES INICIALES}
    INIT MB=3000
    INIT MBA=1000
    INIT MBB=900
    INIT MD=1E-10
    INIT MDA=0
    INIT MDB=0
    INIT XA4=0
    INIT XB4=0
    INIT XA3=0
    INIT XB3=0
    INIT XA2=0
    INIT XB2=0
    INIT XA1=0

    ResponderEliminar
  8. Fernando, prueba a meter como relaciones de equilibrio las que tienen esta forma: y=alfa*x/(1+(alfa-1)*x)

    ResponderEliminar
  9. Hola, aqui os pongo un ejercicio de examen:
    10 moles de una mezcla equimolar de 2 compuestos (A y D) se someten a destilación diferencial.
    El destilado según se va condensando se va alimentando a un tanque en el que se produce la reacción de equilibrio:
    A + B <-----> C
    Kd = 2 L/mol*h
    Ki = 0.5 1/h
    La velocidad de vaporización es 1 mol/h si el destilado es más rico en el compuesto más volatil y
    nulo en caso contrario.
    El volumen inicial del tanque de reacción es 1 litro y el único compuesto presente inicialmente en el tanque es B,
    en una concentración de 5 mol/l.
    No existe ninguna corriente de salida del tanque:
    a) Escribir en B. Madonna el programa que simula la dinámica del proceso.
    b) Comentar la evolución del sistema, dibujando a mano alzada aquellos gráficos q se consideren relevantes.

    Datos:
    DT = 0.01
    Volatilidades relativas:A= 10 y B=1
    Densidad (g/l): A=1000 y B=1500
    Peso molecular (g/mol): A=20 y B=25

    Suerte con su resolución

    ResponderEliminar
  10. Fernando, le he echado un vistazo a la práctica, lo único que creo que debería cambiar es:

    XA1'=(V*YA2-(L+D)*XA1)/M
    XB1'=(V*YB2-(L+D)*XB1)/M

    (pones XA2 Y XB2)

    Yo lo he hecho con esa modificación y sí sale algo bastante razonable, con MB escalonado según los R.

    La ecuación que usas para el equilibrio es la que usamos para mezclas multicomponentes, la ecuación que dice Juan Carlos es para mezclas binarias y siempre que la volatilidad relativa esté referida al componente más pesado (es decir, el alfa del elemento pesado sea 1)
    En el libro de operaciones básicas
    en la página 14 viene la demostración matemática.

    Ánimo a todos!

    ResponderEliminar
  11. Miguel, los compuestos que se someten a destilacion diferencial son A y D ó A y B? Es que las volatilidades relativas que has puesto son respecto a A y B y los pesos moleculares y la densidad tambien.

    ResponderEliminar
  12. Juan Carlos, la relación de equilibrio que propones es para una mezcla binaria y el problema está planteado para una mezcla de tres componentes.

    Rocio, gracias, he probado lo que me sugeriste y salen datos más coherentes, jajaja.

    Gracias

    ResponderEliminar
  13. Hola a todos.
    Me han pasado un ejercicio de otros años (creo que es interesante porque mezcla dos operaciones):

    Se alimentan a un tanque de reaccion continuo 10 mol/h de una mezcla
    binaria formada por los componentes A y B de composición 1 y 0.1 mol/l
    respectivamente.En el tanque tiene lugar la siguiente reacción :



    A↔B

    KD=1
    KI=0.5 1/h

    El volumen del tanque permanece constante e igual a 1 m³.La corriente
    acuosa que sale del tanque se lleva a un extractor liquido-liquido
    donde se pone en contacto con 3 mol/h de una fase organica hacia la
    que se transferirá un soluto y que solo contiene A.La velocidad de
    transferencia de soluto de la fase acuosa que sale del tanque de
    reaccion a la fase orgánica viene dada por la siguiente ecuación :

    Q=10(CAW-CAO)

    ;siendo CAW la concentracion de A a la salida del extractor de la fase
    acuosa y que en ningun caso podra superar la composición de 0.15 mol/
    l y CAO la concentración de la fase orgánica a la salida del
    extractor.

    El volumen de la fase acuosa y la fase órganica es de 1m³.

    Se pide :

    a) Escribir el programa en Berkeley Madonna que simula la dinámica del
    proceso,incluyendo como comentarias el numero de ecuaciones,variables,
    grados de libertad y condiciones iniciales.
    b) Comentar la evolución del sistema y gráficos ue consideres mas
    importantes.


    >>>>Yo he planteado el modelo, tal y como se dice en el apartado (a). Si ha alguien le sale algo distinto, por favor que lo comente por aquí, Gracias. <<<<<<<<<

    METHOD RK4

    STARTTIME = 0
    STOPTIME=10
    DT = 0.02


    {ECUACIONES=8 }

    CAR'=(W*(CAF-CAR))/VR+RA
    CBR'=(W*(CBF-CBR))/VR+RB
    RA=-K1*CAR+K2*CBR
    RB=-RA
    CAW'=(W*(CAR-CAW)-Q)/VW
    CBW'=(W*(CBR-CBW))/VW
    CAO'=(O*CAO+Q)/VO
    Q=10*(CAW-CAO)

    {VARIABLES= 17}

    {GRADOS DE LIBERTAD=9}

    W=1
    O=3
    CAF=1
    CBF=0.1
    VR=1000
    VW=1000
    VO=1000
    K1=1
    K2=0.5

    {CONDICIONES INICIALES=5}

    INIT CAR=1
    INIT CBR=0.1
    INIT CAW=0
    INIT CBW=0
    INIT CAO=0

    ResponderEliminar
  14. Ahí va mi solución al ejercicio que propone Miguel:

    http://rapidshare.com/files/126380792/ejemplo.doc.html

    Si veis algo mal comentadlo por aquí, por favor.

    ¡Gracias!

    (Fernando de nada a mandar :))

    ResponderEliminar
  15. Beatriz Picazo, el ejercicio que plantea Miguel, se resolvió el 6 de Mayo (Ejercicio Propuesto)

    ResponderEliminar
  16. Lleváis toda la razón, me he confundido. Perdona Fernando.

    ResponderEliminar
  17. Juan Carlos, perdon por nada, jejeje. Gracias a ti por tu atención maquina. Nos vemos.Fer

    ResponderEliminar
  18. AÑADE EL CHAT EN LA COLUMNA DE LA IZQUIERDA Y EL CONTADOR DE VISITAS PARA QUE ESTO COJA AMBIENTE!!!!


    ESPAÑA NECESITA... INGENIERÍA QUÍMICA!!!!!!!!!!!!!

    CHAMPIÑONES!!! CHAMPIÑONES!!!
    OE, OE, OE!!!

    ResponderEliminar
  19. Fer, creo recordar que Antonio comentó que, debido a los malos resultados de otros años, decidieron no mezclar dos operaciones en los ejercicios de examen; de todas formas viene bastante bien para repasar dos cosas a la vez. Gracias a todos por los ejercicios que váis poniendo y ánimo para el jueves!

    ResponderEliminar
  20. Rocío, no puedo ver tu enlace. ¿Puedes copiar el archivo en el blog?
    ¡gracias por adelantado!

    ResponderEliminar
  21. Perdón por el error beatriz,
    todos los datos que he puesto al final del ejercicio son referidos a A y D (no a A y B como está puesto).

    Más tarde pondré la solución

    Saludos

    ResponderEliminar
  22. Anónimo2/7/08, 1:11

    Hola de nuevo,
    Antes he puesto un ejercicio de Examen DESTILADOR-REACTOR aquí teneis el modelo:
    {DESTILACIÓN DIFERENCIAL y REACTOR (SEMICONTINUO)}
    METHOD RK4
    STOPTIME= 10
    DT = 0.01

    {ECUACIONES = 16}
    M' = -V
    MA' = -V*YA
    MD' = -V*YD

    XA = MA/M
    XD = MD/M

    SUMM = ALA*XA + ALD*XD
    YA = ALA*XA/SUMM
    YD = ALD*XD/SUMM

    VR' = V*(YA*PmA/RHOA + YD*PmD/RHOD)

    VRCA' = V*YA + (-Kd*CA*CB +Ki*CC)*VR
    VRCB' = (-Kd*CA*CB +Ki*CC)*VR
    VRCC' = (Kd*CA*CB - Ki*CC)*VR

    CA = VRCA/VR
    CB = VRCB/VR
    CC = VRCC/VR

    V = IF YA>YD THEN 1 ELSE 0

    {VARIABLES = 24}
    {GRADOS DE LIBERTAD = 8}

    ALA = 10
    ALD = 1
    RHOA = 1000
    RHOD = 1500
    PmA = 20
    PmD = 25
    Kd = 2
    Ki = 0.5

    {CONDICIONES INICIALES = 7}
    INIT M 10
    INIT MA = 5
    INIT MD = 5
    INIT VR = 1
    INIT VRCA = 0
    INIT VRCB = 5
    INIT VRCC = 0

    ResponderEliminar
  23. ¡Hola!

    Resulta que en el rapidshare del ejercicio que propuso Miguel hay un fallo,

    escribí:

    VRCD'=0, cuando está entrando constantemente D en el caudal de alimentación, lo correcto sería:

    VRCD'=F0*CD0= YD*V.

    Perdonad... demasiado tiempo de exámenes...

    ResponderEliminar
  24. HOLA A TODOS! EL EJERCICIO PUESTO MAS ABAJO ESTA EN EL BLOG EN LA ENTRADA SIMULACRO DE EXAMEN, ES EL QUE SE PUSO EL AÑO PASADO EN EL SIMULACRO; YO HAY VARIAS COSAS Q NO COMPRENDO:
    - YO LO HE HECHO, LAS ECUACIONES ESTAN BIEN, PERO A MI TIME1 ME SALE 20 HORAS, NO 30.

    - Y DESPUES CUANDO HACE EL 2º IF LE RESTA V2 INICIAL (Q SERIA 0 PERO COMO ESTA DIVIDIENDO LE DAMOS EL VALOR 1E-10) F Y V2 NO TIENE LAS MISMAS UNIDADES, CON LO CUAL NO PODRIA INDICARLO ASI NO?!

    A VER Q OS PARECE.


    Un tanque contiene inicialmente 1000L de una disolucion de A con una
    concentracion de 1mol/L. En dicho tanque tiene lugar las reacciones A->B->C
    k1=k2=0.05 h^-1
    Cuando se alcanza el máximo número de moles de B, el contenido del
    tanque se transpasa a un nuevo tanque (inicialmente vacío) a razón de
    100 L/h. En este nuevo tanque tiene lugar la reacción B->D k3=0.5
    h^-1
    Siendo D el producto deseado, se pide:
    a)Programa en Berkeley Madonna.
    b)Rendimiento del producto deseado.
    c)Tiempo óptimo de transpaso.


    METHOD RK4
    STOPTIME=50
    DT = 0.01
    {ECUACIONES = 17}
    ;REACTOR 1
    V1' = -F
    V1CA1' = -F*CA1 + (-K1*CA1)*V1
    V1CB1' = -F*CB1 + (K1*CA1-K2*CB1)*V1
    V1CC1' = -F*CC1 + (K2*CB1)*V1
    CA1 = V1CA1/V1
    CB1 = V1CB1/V1
    CC1 = V1CC1/V1
    ;REACTOR 2
    V2' = F
    V2CA2' = F*CA1
    V2CB2' = F*CB1 + (-K3*CB2)*V2
    V2CC2' = F*CC1
    V2CD2' = (K3*CB2)*V2
    CA2 = V2CA2/V2
    CB2 = V2CB2/V2
    CC2 = V2CC2/V2
    CD2 = V2CD2/V2
    ;TRASPASO
    F = IF TIME < TIME1 THEN 0 ELSE
    IF TIME < TIME1+10-1E-10 THEN 100 ELSE 0
    {VARIABLES = 21}
    {GRADOS DE LIBERTAD = 4}
    K1 = 0.05
    K2 = 0.05
    K3 = 0.5
    TIME1 = 30
    {CONDICIONES INICIALES = 9}
    INIT V1 = 1000
    INIT V1CA1 = 1000
    INIT V1CB1 = 0
    INIT V1CC1 = 0
    INIT V2 = 1E-10
    INIT V2CA2 = 0
    INIT V2CB2 = 0
    INIT V2CC2 = 0
    INIT V2CD2 = 0

    ResponderEliminar
  25. hola! el el ejercicio del traspaso de los tanques, ¿alguien podria explicar la condicion de traspaso? es que no la entiendo. y por otro lado, en cuanto al ejercicio propuesto por miguel y la resolucion, no aparece en el balance VRCD', sin embargo rocio si lo ha puesto (VRCD'=F0*CD0= YD*V) pero D no interviene en la reaccion, ¿habria que ponerlo o no?
    gracias!

    ResponderEliminar
  26. Maria en el ejercicio del traspaso creo que te refieres al que hicimos en clase el último día en clase, la condicion, es que mientras el caudal F vale cero hasta que la concentracion de B sea máxima, lo hacemos con tiempos mediante un liders podemos ver a que tiempo la concentracion de b se hace maxima, luego cuando el tiempo es menor que otro tiempo vale el caudal que le demos pero a ese segundo tiempo tenemos que sumarle dos que son las horas que tarda el traspaso, esto dos tiempos para hacer la condicional debemos hacerlos a posteriori es decir una ver que hemos ejecutado el modelos y nos corre bien modificamos los tiempos de la condicional correctos un saludo.

    ResponderEliminar
  27. Hola!tengo una duda:imaginad que tenemos un reactor donde ocurren dos reacciones exotérmicas:

    A+B-->C r1=k1*CA*CB
    P*Q-->Z r2=k2*CP*CQ

    El término de generacion de calor en el balance entálpico quedaria:

    V*(r1*H1+r2*H2)/M*CP

    si lo desarrollo poniendo las velocidades de reacción sería:

    V*(-k1*CA*CB*H1-k2*CP*CQ*H2)/M*CP

    Mi duda es si estoy en lo cierto poniendo las 2 velocidades con signo negativo para que al multiplicarlas por calores de reaccion negativos(reacciones exotermicas)den un término positivo, lógico si pensamos que la temperatura del reactor va a ir aumentando.

    Un saludo y gracias!

    ResponderEliminar
  28. María, yo le he hecho balances a todos los elementos que están en el reactor, su concentración sí variará, si te interesa sabre la concentración final de D lo puedes poner, sí es muy importante tenerlo en cuenta de cara al cálculo de F, ya que de ello dependerá el volumen del reactor. Añadirlo o no... sería el balance de un inerte; creo que daría igual; se trata de una variable más y una ecuación más no influye en los grados de libertad...

    ResponderEliminar
  29. Aqui os dejo unejercicio de otros años, os dejo el enunciado y en el enlace resuelto:

    La reacción A->B->C se da en un reactor continuo mezcla perfecta, inicialmente está lleno de disolvente, se uiere obtener B como producto deseado.
    En la corriente de alimentación de alimentación contiene solo A con una concentración de 1 mol/l.
    El volumen de reactor es de 10 l.

    a) Programa en Berkeley Madonna.
    b) Ecuaciones, variaables y grados de libertad.
    c) Volumen óptimo del reactor que hace máximo la concentración de B, a la salida del reactor en estado estacionario.¿Cúal es este valor máximo? Considerar que se pueda disponer de reactivo entre 10-500 l en intervalos de 10 l.

    DATOS: k1=2 h-1 k2=2 h-1 DT=0.1

    Aqui os dejo el enlace:

    http://rapidshare.com/files/126536620/EJERCICIO_EXAMEN_SIMOP2.doc.html

    ResponderEliminar
  30. Aqui os dejo unejercicio de otros años, os dejo el enunciado y en el enlace resuelto:

    La reacción A->B->C se da en un reactor continuo mezcla perfecta, inicialmente está lleno de disolvente, se uiere obtener B como producto deseado.
    En la corriente de alimentación de alimentación contiene solo A con una concentración de 1 mol/l.
    El volumen de reactor es de 10 l.

    a) Programa en Berkeley Madonna.
    b) Ecuaciones, variaables y grados de libertad.
    c) Volumen óptimo del reactor que hace máximo la concentración de B, a la salida del reactor en estado estacionario.¿Cúal es este valor máximo? Considerar que se pueda disponer de reactivo entre 10-500 l en intervalos de 10 l.

    DATOS: k1=2 h-1 k2=2 h-1 DT=0.1

    Aqui os dejo el enlace:

    http://rapidshare.com/files/126536620/EJERCICIO_EXAMEN_SIMOP2.doc.html

    ResponderEliminar
  31. El ejercicio que propuso miguel, se podria resolver tambien asi?

    METHOD RK4

    STARTTIME = 0
    STOPTIME=20
    DT=0.1

    {ECUACIONES=26}
    {RESIDUO}
    M'=-V
    MXA'=-V*YA
    MXD'=-V*YD
    XA=MXA/M
    XD=MXD/M
    YA=(ALA*XA)/SUMM
    YD=(ALA*XD)/SUMM
    SUMM=ALA*XA+ALD*XD

    {DESTILADO}
    D’=V
    DXAD’=V*YA
    DXDD’=V*YD
    XAD=DXAD/D
    XDD=DXDD/D

    {TANQUE}
    VR'=F0
    VRCA1'=F0*CA0+RA*VR
    VRCB1'=RB*VR
    VRCC1'=RC*VR
    RA=-KD*CA1*CB1+Ki*CC1
    RB=-KD*CA1*CB1+Ki*CC1
    RC=KD*CA1*CB1-Ki*CC1
    CA1=VRCA1/VR
    CB1=VRCB1/VR
    CC1=VRCC1/VR

    V=IF YA>YD THEN V0 ELSE 1E-10
    F0=(D*XAD*PMA)/RHOA+(D*XDD*PMD)/RHOD
    CA0=(D*XAD)/F0
    {VARIABLES=35}
    {GRADOS DE LIBERTAD=9}
    V0=1
    ALA=10
    ALD=1
    PMA=20
    RHOA=1000
    PMD=25
    RHOD=1500
    KD=2
    Ki=0.5
    {CONDICIONES INICIALES=10}
    INIT M=10
    INIT MXA=5
    INIT MXD=5
    INIT VR=1
    INIT VRCA1=0
    INIT VRCB1=5
    INIT VRCC1=0

    ResponderEliminar
  32. Hola. Tengo una dudilla. Alguien sabe si entra en el examen la parte de control?

    Gracias

    ResponderEliminar
  33. La parte de control si entra en el examen, un saludo

    ResponderEliminar
  34. Elena, no puedo ver el enlace porque se corta; pongo http://rapidshare.com/files/126536620/EJERCICIO_EXAMEN.html
    y me dice que no lo encuentra o algo así.

    ¿Podrías ponerlo de nuevo?

    Muchas gracias

    ResponderEliminar
  35. He intentado hacer el ejercicio de la mezcla multicomponente en destilación por cargas con array. Tras ver la entrada de Nacho (gracias, gracias y gracias) en la práctica 3 encontré el fallo y me salen datos coherentes, ahí os dejo el programa

    ¡ya mismo acabamos!

    METHOD RK4

    STARTTIME = 0
    STOPTIME=1000
    DT = 0.02

    {NÚMERO DE ECUACIONES=4*NPLATES+6 }

    {CALDERÍN, 3 ECUACIONES}

    MB'=L-V
    D/DT(XA[NPLATES])=(L*(XA[i-1]- XA[i])-V*YA[i])/MB
    D/DT(XB[NPLATES])=(L*(XB[i-1]- XB[i])-V*YB[i])/MB


    {PISOS, 2NPLATES-4}

    D/DT(XA[2..NPLATES-1])= (L*(XA[i-1]-XA[i]) + V*(YA[i+1]-YA[i]))/M
    D/DT(XB[2..NPLATES-1])= (L*(XB[i-1]-XB[i]) + V*(YB[i+1]-YB[i]))/M


    {CONDENSADOR, 4 }

    L= V-D
    D=V/(1+R)

    D/DT(XA[1])=(V*YA[i+1]-(L+D)*XA[i])/M
    D/DT(XB[1])=(V*YB[i+1]-(L+D)*XB[i])/M


    {DESTILADO ACUMULADO, 5}
    MD'=D
    MDA'=D*XA[1]
    MDB'=D*XB[1]
    XAD=MDA/MD
    XBD=MDB/MD

    {EQUILIBRIO, 2NPLATES-2}
    YA[2..NPLATES]= ALA*XA[i]/ (ALA*XA[i]+ALB*XB[i]+ALC*(1-XA[i]-XB[i]))
    YB[2..NPLATES]= ALB*XB[i]/ (ALA*XA[i]+ALB*XB[i]+ALC*(1-XA[i]-XB[i]))


    {VARIABLES=4*NPLATES+19}

    {GRADOS DE LIBERTAD=11}
    ALA=7
    ALB=5
    ALC=1
    V=12
    M=20
    NPLATES = 5
    R=IF TIME"menor"TIME1 THEN 1E20 ELSE
    IF TIME"menor"TIME2 THEN R1 ELSE
    IF TIME"menor"TIME3 THEN 1E20 ELSE R1
    R1=8
    TIME1=250
    TIME2=500
    TIME3=900

    {CONDICIONES INICIALES}
    INIT MB=3000
    INIT XA[NPLATES]=0.333
    INIT XB[NPLATES]=0.3
    INIT MD=1E-10
    INIT MDA= 0
    INIT MDB= 0
    INIT XA[1..NPLATES-1]= 0
    INIT XB[1..NPLATES-1]= 0

    ResponderEliminar
  36. Aqui os dejo otro ejercicio de examen, os pongo el enunciado y en el enlace os lo pongo resuelto:

    80 moles de una mezcla binaria equimolar se somete a destilación diferencial durante 5 horas vaporizamos, la velocidad de vaporización es de 10 mol/h.Se pide:

    a) Programa en Berkeley Madonna
    b) Ecuaciones, variables y grados de libertad.
    c) Calcular la mínima vlatilidad reltiva para que la fracción molar en el residuo, del compuesto más volátil sea menor a 0.15 al final de la operación.

    DATOS:DT=0.01
    ALFA=3
    ALFB=1.2

    Aqui os dejo el enlace:

    http://rapidshare.com/files/126591152/EJERCICIO_EXAMEN_SIMOP2_DOS.doc.html

    ResponderEliminar
  37. os pongo de nuevo los dos enlaces de los ejercicios, es que contan

    ejercicio primero el del reactor:
    http://rapidshare.com/files/
    126592048/EJERCICIO_EXAMEN_SIMOP2.doc.html

    el de la destilacion diferencial:

    http://rapidshare.com/files/
    126591152/EJERCICIO_EXAMEN_SIMOP2_DOS.doc.html

    A

    ResponderEliminar
  38. Hola Talia,
    Esa misma pregunta me la hice yo esta mañana y llegué a la misma conclusión que tú. Da igual que las reacciones sean exotérmicas ó endotérmicas, si se define la velocidad de reacción como positiva, siempre se debe poner el término de generación de calor debido a reacción química negativo para que si tenemos una reacción exotérmica T aumente y si es endotérmica T disminuya.
    Espero haberte aclarado la duda.

    ResponderEliminar
  39. Bea, creo que si tienes en cuenta que el reactor se alimenta a partir del destilado, al plantear los balances al destilado, además de tener en cuenta la entrada que viene de la destilación, habría que tener en cuenta una salida de corriente hacia el reactor. Por lo que, si supones que la acumulación en el destilado es cero, la entrada al destilado es igual a la salida.
    Ánimo con el examen!

    ResponderEliminar
  40. Bueno, yo creo q ya está todo el pescao vendio, no?
    Ya solax hace falta un poco de suerte para todos.
    P.D.: La resolucion correcta del ejercicio que propuse es la que yo puse (valga la redundancia).
    Así fue como la hizo antonio el año pasado.

    Saludos, tranquilidad y buenos alimentos. jejej

    ResponderEliminar
  41. hola de nuevo!
    Elena, no se puede definir YA y YB de ese modo porque, como ya dijo Rocío, la volatilidad relativa debe estar referida al componente más pesado (es decir, ALFB debe ser 1). Hay que poner:
    SUMA=XA*ALFA+XB*ALFB
    YA=(XA*ALFA)/SUMA
    YB=(XB*ALFB)/SUMA
    De este modo, la suma de las fracciones molares tanto en el residuo como es el destilado es siempre igual a 1
    ¡Ya queda menos!

    ResponderEliminar
  42. Hola, en la practica de Control de T cuando pone NORMAL(298,2.5)
    El 2.5 que índica?

    ResponderEliminar
  43. ¡hola por tercera vez!
    miguel, el 2.5 es la desviación estándar de la distribución normal.
    hasta mañana!

    ResponderEliminar
  44. lo que indica 2,5 es el varianza, es una distribucion normal. si te fijas para el ultimo caso es 313, 5. mirando el grafico la "dispersion" es mayor para 5 que para 2,5. perdon por las palabras porque no se si son las adecuadas, pero espero que mas o menos lo entiendas

    ResponderEliminar
  45. nuria lo ha dicho mejor, desviacion estandar.

    ResponderEliminar
  46. hola otra vez, tengo otra duda, en el ejercicio de examen del reactor de elena, ¿por que F=200?

    ResponderEliminar
  47. Talia y Juan Carlos, yo lo que hago es poner el termino como:
    (V*AHA*rA+V*AHP*rP)/(V*RHo*CP)
    siendo rA y rP las velocidades de componentes que serian:
    ra=-k1*cA*cB
    rp=-k2*cP*cQ
    Esto es porque en reactores lo haciamos con el reactivo limitante (A y P)aqui en realidad no te dice cual es el reactivo limitantel, pero weno sera del ke te den la entalpia no? jejeje

    ResponderEliminar
  48. ¡Hola! (ya por última vez... que esto se acaba...) La duda que decía Juan Carlos y Talia, en reactores poníamos el calor de reacción ( es decir -H, que por eso es el cambio de signo) respecto a un reactivo o producto (por mol de A, de B...) y se multiplica por la velocidad de reacción referida también al mismo elemento.
    No sé si con esto os ayudo, de todas maneras yo acabo haciendo lo mismo que decís vosotros.

    ¡¡¡ mañana acabamos!!! ánimo!!!!

    ResponderEliminar
  49. Talia y Juan Carlos, yo tenía la misma duda respecto a los signos, pero claro, para que se cumpla que xa reacción exotérmica T aumenta y si es endotérmica T disminuya, el signo deber ser ese.

    ResponderEliminar
  50. hola gente, ¿alguno puede decirme pq en la práctica 14 la temperatura inicial de salida del calentador (INIT TR) es de 290 K si, a tiempo cero, la entrada se realiza a 298 K? Me entra la duda de si se trata una variable de diseño o no. gracias y suerte a todos

    ResponderEliminar
  51. El ejercicio propuesto por Elena de la destilación diferencial lo he hecho asi:

    METHOD RK4

    STARTTIME = 0
    STOPTIME=5
    DT = 0.01

    M'=V
    MXA'=-V*YA
    MXB'=-V*YB
    XA=MXA/M
    XB=MXB/M
    D'=V
    DXAD'=V*YA
    DXBD'=V*YB
    XAD=DXAD/D
    XBD=DXBD/D
    YA=(ALA*XA)/SUMM
    YB=(ALB*XB)/SUMM
    SUMM=ALA*XA+ALB*XB

    ALA=3
    ALB=1.2
    V=10

    INIT M=80
    INIT MXA=40
    INIT MXB=40
    INIT D=1E-10
    INIT DXAD=0
    INIT DXBD=0

    La alfa minima me sale 0.66 ¿a alguien le sale igual q a mí?

    ResponderEliminar
  52. hola!
    Ana Rosa, te has despistado con el signo, es:
    M'=-V
    Poniéndo así, me sale que el ALFA mínima es 5,82. Para comprobar que no tienes este tipo de fallos tontos (pero que estropean todo el ejercicio) asegúrate que las fracciones molares suman 1 en todo momento.
    Espero haberte servido de ayuda!

    ResponderEliminar
  53. Ana a mi tambien me sale ese resultado

    ResponderEliminar
  54. Gracias Nuria es verdad, que despiste más tonto esperemos q no me pase eso en el examen.

    ResponderEliminar
  55. Bff...ana a mi me ha pasado igual, no te preocupes que estas horas ya son muy malas,saludos!

    ResponderEliminar
  56. Dudas y más dudas:
    El X1star en Extracción Liq-Liq, que representa?

    ResponderEliminar
  57. Miguel, el Xstar es la composición que tendría que haber en el líquido para que el sistema esté en equilibrio con la la otra fase (que tiene concentración Y). Como X>Xstar, el sistema tiende a transferir el compuesto para alcanzar el equilibrio (que se alcanzaría cuando X=Xstar.
    Un saludo!

    ResponderEliminar
  58. En el problema que ha propuesto carmen guerrero, sabeis como se hace el rendimiento y el tiempo optimo de trasbase??
    Es que el rendimiento vimos en reactores que era moles de B transformados en P/moles de B alimentados, que yo lo traduciria en:
    R=VCD(a tiempo final)/VCB(a tiempo 30)
    Ademas el TIME1 me sale 20 y a carmen 30.
    Cuando dice lo del tiempo optimo de trasbase, que se refiere al tiempo que tarda en trasbasar de un reactor a otro o al tiempo que tarda en en producirse el trasbase (TIME1) yo creo ke es lo 1º ke he dicho.
    ESpero no haberos liado mucho!
    Suerte para el examen!

    ResponderEliminar
  59. Nuria, gracias x la duda del Xstar
    suerte para mañana

    ResponderEliminar
  60. hola!
    tengo una propuesta que hacer para otros exámenes con respecto a grapar las hojas desde el comienzo. No me parece mal que se grapen si lo que se quiere es evitar que la gente se pase los folios, pero creo que podrían graparse los 4 folios y dejar el texto con el enunciado libre. De este modo, no tendríamos la molestia de ir leyendo el enunciado en una hoja a la vez que escribimos en otra grapada a ésta.
    ¡Disfrutad del verano!

    ResponderEliminar
  61. YO ESTOY CON NURIA, FUE UN POCO PESADO TENER Q ESTAR CON LAS 4 FOLIOS GRAPADOS, PASANDO DE DELANTE A DETRAS PARA LEER EL ENUNCIADO.

    ResponderEliminar
  62. Nuria, tienes toda la razón en eso, de hecho yo me he pasado por el blog para decir exactamente lo mismo. Yo creo que si lo decimos unos cuantos yo creo que se tendrá en cuenta, además, para estas cosas está el foro, para mejorar algunos aspectos como este...suerte con las notas!!!

    ResponderEliminar
  63. Me uno a vosotros, yo se lo dije incluso antes de que lo grapásemos...

    Así que recomendaciones para próximos años: en el cuaderno de prácticas deberían estar los enunciados escritos, en el examen deberían dejar la hoja de examen (aunque sea) suelta comodidad y más ejercicios en los que se usen condicionales, hacer más ejemplos en prácticas.

    YA ESTAMOS DE VACACIONES!!!

    ResponderEliminar