Lee Rudd

Tema 3: Simulación
estacionaria.
Modelado y simulación en Ingeniería Química. Manuel Rodríguez

2
1.Introducción
1.1 Funcionalidades de los simuladores en estado estacionario.
Problemas de simulación
Dadas unas entradas y parámetros computar la salida
Estudios de sensibilidad
Como la simulación pero permite estudiar la influencia de ciertas variables,
variándolas sistemáticamente.
Problemas de diseño
Estudia diferentes conjuntos de especificaciones. Se puede “quitar” un parámetro y
añadir una nueva especificación. La ejecución dará un valor a ese parámetro.
Más difícil de resolver, y de mal especificar.
Problemas de optimización
Se plantea una función objetivo para mejorar el proceso. Se añaden desigualdades. Se
“quitan” especificaciones (grados de libertad para la optimización) y se resuelve el
problema.
Problemas de síntesis
Se conocen las entradas y las salidas y se quiere obtener el diagrama de flujo, así
como las condiciones de operación y los parámetros de los equipos.

3
Simulación y análisis de sensibilidad

4
Diseño de especificación

5
Optimización

6
Síntesis de procesos

7
1.2 Arquitectura de software de un entorno de modelado y simulacion.

8
Pasos para realizar una simulación.
1 Seleccionar un conjunto de unidades consistente.
2 Especificar los componentes (especies químicas) que participan en el proceso.
3 Seleccionar un paquete de propiedades físicas
4 Definir la estructura del diagrama de flujo. Qué unidades participan y los

flujos (corrientes) de materia y/o energía entre ellas..
5 Seleccionar un modelo para cada unidad de proceso.
6 Satisfacer los grados de libertad del proceso
7 Seleccionar las opciones de resolución numérica
8 Ejecutar el modelo
9 Examinar los resultados

9
1.3 Estrategias de resolución: Secuencial Modular vs. Orientado a ecuaciones
Secuencial Modular Orientado a ecuaciones
Se simula una unidad cada vez Se resuelven todas las unidades a la

vez
Se descompone el diagrama de flujo Ordenación de ecuaciones
Iterar en las corrientes de “rasgado” Se actualizan todas las incógnitas

simultáneamente
Menos flexible, más robusto Más flexible, menos robusto
Inicialización es importante Inicialización es muy importante
No necesita mucho almacenamiento Requisitos de almacenamiento

grandes.
· Muy empleados en la actualidad: ej. Aspen Plus · Arquitectura preferida en nuevos simuladores: ej.
(AspenTech), ChemCAD, PROII (SimSci), Hysim Aspen Custom Modeler (AspenTech), RTO-OPT
(Hyprotech), Hysys (Hyprotech, comprada (AspenTech), NOVA (Nova), gPROMS, ABACUSS.
recientemente por AspenTechnologies)

10
Sistema de ecuaciones no lineales
Orientado a ecuaciones Secuencial Modular

Asignación de variables de entrada Variables de entrada asignadas
(Lee-Christensen-Rudd)
(Dulmage-Mendelsohn decomposition)
Particionado (identificación de ciclos)
Orden de resolución
Output assignment (Norman)
(Steward) (Keham & Sachem)
(Duff) (y detecta singularidad estructural) (Sargent&Westerberg)
Particionado (matriz de incidencia) (Tarjan)
Busca matriz triangular Selección de variables de corte (tear)
Orden de resolución (Barkeley & Motard)
(Tarjan) (Westerberg&Motard)
(Sargent&Westerberg) (Upadyhe&Grens)(Mejor convergencia)
(Christensen&Rudd)
Incialización
Resolución del sistema: 3 niveles

Resolución del sistema. En cada módulo:
Cada sistema (bucle algebraico): Propiedades físicas
Wegstein Ecuaciones del módulo
Aprox. Sucesivas Diagrama de flujo
Newton
Broyden Wegstein
Aprox. Sucesivas
Resolución de LAEs Newton
(métodos Gaussianos) Broyden
Resolución de LAEs
(métodos Gaussianos)

11
2. Estrategia secuencial-modular
Grafos
Un grafo G es un conjunto de vértices (V) y un conjunto de arcos que

unen los vértices.
5
4
1 3
2
Digrafo: es un grafo en el cual los arcos tienen una dirección. La

dirección que marcan indica el flujo de una cantidad.
5
4
1 3
2

12
Matriz de adyacencia o estructural vértices

1 2 3 4 5
5 1 1 Matriz
vértices
4 2 1 simétrica
1 3 3 1 1
Indica que vértices
2 4 1 1 1 1 están unidos
5 1 1
v. destino
1 2 3 4 5
1 1
d
v. origen
2 Son valores lógicos.
5
a 4 Booleanos.
3 1 1
1 b c e
3 4 1 1
5
2
arcos
a b c d e
1 -1 Matriz
2 +1 signos
vértices
3 -1 -1 Negativo: salida
4 +1 -1 +1 -1 Positivo: llegada
5 +1 +1

13
Identificación de bucles
Sobre la matriz de adyacencia
5
4
1 3
2
conectados en dos pasos

v. destino
1 2 3 4 5 1 2 3 4 5
1 1 1 1 1 1
v. origen
2 1 bucles
2 1 1
A= 3 1 1 A2= 3 1 1
4 1 4 1 1
5 1 5 1
1 2 3 4 5 1 2 3 4 5
1 1 1 1 1 1 1 1
2 1 1 2 1 1 1
A3= 3 1 1 1
A4= 3 1 1 1 1
bucles
4 1 1 4 1 1 1
5 1 1 5 1 1
conectados en tres pasos conectados en cuatro pasos

14
Poner el digrafo correspondiente
1 2 3 4
6 5

15
Diagrama de flujo
¿Hay bucles en el diagrama?¿Cuántos?¿Se pueden identificar de forma automática?

¿Cuál es el mejor orden para resolverlos?
¿Cómo resolver los subsistemas que forman parte de un bucle?

16
Particionado (y ordenación)
El particionado consiste en la identificación de los subsistemas

que están interrelacionados. Es decir, que constituyen un bucle y
que hay que tratar en su conjunto a la hora de resolver el
diagrama de flujo.
Algoritmo de Sargent y Westerberg
1)  Seleccionar un vértice (de forma arbitraria)

2)  Trazar un camino desde el vértice siguiendo un arco no explorado
(formando un conjunto con los vértices que aparecen) hasta que:
a)  Se encuentra un vértice que no tiene salida (o no va a otra unidad).
Borrar el vértice (y sus arcos incidentes) y ponerlo en cabeza de una
lista. Continuar desde el vértice anterior.
b)  Se encuentra un vértice que ya está en el conjunto de vértices
(bucle). Agrupar todos los vértices del bucle en un único conjunto.
3)  Si queda algún vértice repetir otra vez desde el paso uno hasta que no
quede ningún vértice.
Otros algoritmos de particionado que se basan en la matriz de adyacencias: algoritmo de Norman y el

algoritmo de Keham y Shacham

17
Ejemplo
A
B C D
E F G
Empieza por C
C àDàFàG (sin salida) Lista: G
B C D
E F
Sigue con F (sin salida) Lista: F G

18
B C D
Sigue con D (sin salida) Lista: D F G
B C
Sigue con CàE (sin salida) Lista: E D F G
A
Sigue con C (sin salida) Lista: C E D F G
B C
y ya igual A y B
Lista: A B C E D F G Las listas NO son únicas

19
Ejemplo
A
B C D
E F
Empieza por C
C àDàFà (sin salida) Lista:F
B C D
C àDàAàBàC (agrupo en un BUCLE)

20
CDAB
CDABàE (sin salida)
Lista: E F
CDAB (sin salida) Lista: [CDAB ] E F

21
Ejercicio: Particionar el siguiente diagrama de flujo.
IJKLMNL
IJK[LMN]
IJK[LMN]OPK
IJ[KLMNOP]
I [JKLMNOPSQR]
FH
FGCDEABC
FG[CDEAB]
FG[CDEAB]F
[FGCDEAB]

22
1 2 3 4 5 6 7 8 9 10 11 12
1 à2 3 6 9à12 4 5à7 8à10à11

23
React. CloracionàColumna àResto

24
Procedimiento de resolución secuencial modular
Especificar
parámetros Especificar
unidades entradas
Poner en
(unidad)k=0
Resolver
unidad k
SI
Es k>N? Solución
N=num unidades final
NO
Resultados
Poner en
k=k+1
¿Qué ocurre si hay un reciclo?

25
Selección de variables de corte. (Rasgado)
Para resolver un grupo de unidades que están interconectadas con algún

bucle es necesario dar un valor inicial a una corriente desconocida. Esa
corriente es la corriente de rasgado y sobre ella es necesario iterar para
obtener el valor final.
corriente de rasgado
valor estimado valor calculado
¿Cómo seleccionar la mejor corriente de rasgado?

¿Cómo actualizar el valor estimado tras cada
iteración para tener una rápida convergencia?

26
Algoritmo de Barkeley & Motard (número mínimo de corrientes)
1)  Realizar el grafo dual. Las corrientes son nodos y los equipos que las unen
son arcos.
2)  Se recorre el grafo. Si una corriente tiene sólo una entrada es eliminada
del grafo (y sustituida por la entrada).
3)  Cuando se llega a un autobucle (autociclo) ésta es una corriente de corte.
Se elimina del grafo y se vuelve al punto 2) hasta eliminar el grafo.
NOTA: Si se llega a una situación donde el algoritmo no puede avanzar se
escoge una variable de corte y se procede.
Criterios para escoger variables de corte:

•  Conocer de forma aproximada el valor que va a tener en la solución.
•  El número de variables que tenga dicha corriente.
•  El número de equipos a los que entra una corriente
•  etc,...

27
Ejemplo S1
S2 S6
Proceso
2 3 6 9
S3
S4 S5
6
S1 S2 S6 9
Grafo dual de
2
corrientes 2 6 6 9
3 6
S3 S4 S5
Las corrientes S1, S4 y S6 sólo tienen una entrada, luego

son eliminadas del grafo.

28
6
S2 9
2 2 6
6 9
3 6
S3 S5
S5 es un autociclo, se
elimina, junto con las S2
corrientes que entran y
salen de el.
S5 es la primera S3 S5
corriente de rasgado.
S3 S2
Se escoge según algún
criterio una de las dos como
corriente de corte. S2
Conjunto final:
S2 es la primera S2
corriente de rasgado. {S5,S2}

29
Ejercicio. Seleccionar las corrientes de rasgado.

30
Grafo del proceso
8
K 5 O 6
1 4
2 7
L M S
3
1 K 8
L S
M ¿Cuáles son las corrientes de corte?
2 7
K L L M S
3 Conjunto final:
L
{S1,S2}
5 4 6
O O

31
Algoritmo de Upadyhe y Grens
1)  Identificar los bucles presentes

a)  Empezar por una unidad cualquiera hasta que aparezca una repetida
b)  Cuando aparezca una repetida anotar las corrientes que conforman el bucle
hallado
c)  Seguir por la corriente anterior a la corriente repetida hasta que vuelva a repetirse
otra unidad
d)  Continuar hasta que se han recorrido todas las corrientes
2)  Escoger la familia de corrientes de rasgado que rompe los bucles una única vez.
a)  Escoger un conjunto de corrientes que rasga todos los bucles una única vez
b)  Seleccionar una unidad que tenga todas sus salidas en ese conjunto de corrientes
c)  Reemplazar en ese conjunto de corrientes de rasgado las salidas de esa unidad por
sus entradas
d)  Repetir hasta que no quedan conjuntos de rasgado nuevos por salir
e)  Todos los conjuntos de rasgado que han aparecido son equivalentes. Como opción
definitiva escoger un conjunto con un menor número de corrientes.

32
S1
Ejemplo
S2 S6
2 3 6 9
S3
S4 S5
1)  Identificar los bucles presentes
S3 S4 S5 S1
2 3 6 9 2 Bucle 1: {S3,S4,S5,S1}
S3 S4 S5 S1
2 3 6 9 2
S6 Bucle 2: {S5,S6}
6
S3 S4 S5 S1
2 3 6 9 2
Bucle 3: {S3,S4, S2}
S6
6
S2
2

33
Corriente S1 S2 S3 S4 S5 S6
Bucle
1 X X X X
2 X X
3 X X X
2) Escoger la familia de corrientes de rasgado que rompe los bucles una única vez.
2 9 6 3
{S3,S6} {S1,S2,S6} {S2,S5} {S6,S4} {S6,S3}
El objetivo es rasgar todos los bucles para poder integrar, la idea es que si se escoge un
conjunto que rasga más de una vez un bucle (ej: {S1,S2,S5}), esto implica que la
convergencia va a ser más lenta.

34
Ejercicio. Seleccionar las corrientes de rasgado.
Grafo del proceso
K 5 O 6
1 4
2 7
L M S
3

8 35
K 5 O 6
1 4
2 7
L M S
3
2 7 8 1
Bucle 1: {2,7,8,1}
L M S K L Bucle 2: {1,4,5}
3 4
5
Bucle 3: {8,1,4,6}
L O K
Bucle 4: {2,3}
6
S
Corriente 1 2 3 4 5 6 7 8
Bucle
1 X X X X
2 X X X
3 X X X X
3 X X
K S M O
{1,3} {5,8,3} {5,6,7,3} {5,6,2} {4,2}

36
Realizar el particionado y seleccionar las corrientes de corte mediante el

algoritmo de Barkeley del siguiente digrafo.

37
Seleccionar las corrientes de corte mediante el algoritmo de Barkeley del

siguiente digrafo.
S6
5 4
S1 S3
S5
S7 S2
S8 2 3
1 S4

38
Procedimiento de resolución secuencial modular (con reciclo)
Especificar
parámetros Especificar
unidades entradas
Poner en
(unidad)k=0 ¿Qué ocurre si hay un diseño de especificación?
Dar estimaciones
iniciales al reciclo
Resultados
Poner en Poner Dar nuevas
k=k+1 k=0 estimaciones al reciclo
NO
SI Comprobar
Es k>N?
N=num unidades error<tolerancia
NO SI
Resolver Solución
unidad k final

39
Diseño de especificación Controladores computacionales.
Diseño de especificación
dado
¿x1? y1
Unidad 1 Unidad 2 Unidad 3
¿x2? Reciclo
x1
¿Alguna relación?

40
Diseño de especificación Set Point
CC
medido
¿x1? y1
Iteración implícita
Reciclo
x2,guess
x2,calc
x1
Iteración explícita
¿Pueden resolverse por el método de sustitución sucesiva?

41
Diseño de especificación Set Point
CC
dado
¿x1? y1
Iteración implícita
Hay que iterar sobre una corriente determinada
y1-ySP
h(x)=0; h(x)=
x2,guess- x2,calc
Reciclo
x2,guess
Conv
x2,calc
x1
Iteración explícita
Se puede iterar sobre cualquier corriente

42
“Realizar la simulación del sistema de forma gradual”.

“Utilizar resultados de simulaciones anteriores (más sencillas) para las nuevas simulación”

43
3. Estrategia orientada a ecuaciones

Matriz de incidencia
Matriz (booleana) que indica qué variables aparecen en qué
ecuaciones de un sistema.
Sistema de ecuaciones

44
Grafo bipartito
Grafo en el que hay dos conjuntos de vértices y todos los arcos van
de un vértice del primer conjunto a un vértice del segundo
conjunto.
x1 ec1
x2 ec2
x3 ec3
Se emplea para identificar el conjunto de variables a especificar.

45
Asignación de variables de salida.

Qué variable se va a calcular con qué ecuación.
Dado un sistema cuadrado y un conjunto de salida asignado se

puede obtener un grafo dirigido.
ec1 ec2 ec3
¿Qué implica la presencia de bucles?

Se emplea para realizar el particionado del sistema.

46
Grados de libertad e identificación de un conjunto

de variables a especificar (variables de entrada).
Algoritmo de Christensen-Rudd
1)  Se emparejan las variables que sólo aparecen en una ecuación o las ecuaciones que sólo
tienen una variable. Se eliminan los nodos (y sus arcos) del grafo.
2)  Se procede así sucesivamente hasta que quedan nodos sin arcos, estos son los grados de
libertad.
3)  Si se llega a una situación donde todos los nodos tienen 2 ó más arcos, es necesario
asignar (de forma arbitraria) una variable a una función y proceder a continuación
como en el punto uno.
NOTA: Este algoritmo proporciona una secuencia de resolución acíclica, excepto en el caso del
punto 3) que indica la presencia de bucles algebraicos. En ese caso las variables asignadas
arbitrariamente serían variables de rasgado que permitirían resolver el sistema de forma
acíclica. El uso de variables de rasgado en simulaciones orientados a ecuaciones está hoy
en desuso (NO en simulación secuencial-modular!)

47
Ejemplo
x1 f1
x2 f2
x3 f3
x4 f4
x5 f5
x6 La variable 7 sólo aparece en la ec 5

La variable 5 sólo aparece en la ec 4
x7
La variable 3 sólo aparece en la ec 3
x8

48
x1 f1 La variable 8 sólo aparece en la ec 2
x2 f2 La variable 6 sólo aparece en la ec 1
Quedan como grados a especificar:
x4
x1 x2 x4
x6
Nótese que cuando hay varias posibilidades el
criterio ha sido arbitrario. Es mejor apoyarse en
el conocimiento heurístico del modelo.
x8

49
Descomposición de Dulmage-Mendelsohn.
En lugar de actuar sobre el grafo bipartito actúa sobre la matriz de

incidencia.
Descompone la matriz de incidencia en submatrices:

• Una parte sobredeterminada: indica que hay más ecuaciones que incógnitas, puede ser
alguna redundante,...
• Una parte cuadrada: tiene forma de matriz triangular inferior. El subsistema está
correctamente especificado y por tanto no aporta grados de libertad.
• Una parte infradeterminada: indica que hay más variables que ecuaciones y que por tanto
hay un número de grados de libertad que hay que especificar. La forma de proceder es
seleccionar uno, volver a realizar la descomposición de D-M y escoger otra. Así hasta que se
obtiene únicamente una matriz cuadrada triangular inferior por bloques.

50
Algoritmo de asignación de variables de salida.

Algoritmo de Steward
• Selecciona la fila (columna) con menor número de incidencia.

• En esa fila (columna) selecciona la variable que pertenece a la columna (fila) con menor
número de incidencia.
• Asigna esa variable (fila) a la ecuación (columna) y se eliminan de la matriz de incidencia.
Asigna la variable 2 con la ecuación 1


51
Singularidad estructural
Cuando no se puede obtener un conjunto de salida la matriz es
estructuralmente singular.
Implica que hay un subconjunto de ecuaciones que tiene menos
variables que el número de ecuaciones del subconjunto.
La especificación de las variables de entrada puede hacer que un

sistema de ecuaciones sea o no estructuralmente singular.

52
Particionado del sistema de ecuaciones y orden
de resolución.
Algoritmo de Tarjan
Más eficiente que el de Sargent y Westerberg
(también aplicable). Obtiene una matriz
triangular inferior por bloques.
Cada bloque es un sistema de ecuaciones que
hay que resolver de forma simultánea.
ec1 ec2 ec3

53
ALGORITMO DE TARJAN
•  Se va recorriendo el grafo (se empieza por cualquier vértice) y se

ponen los nodos que aparecen por primera vez en una camino
(path) y en una pila (stack)
•  Según van apareciendo se les asigna un puntero
•  Si de un vértice (por ejemplo f) se va a otro que ya está en el
camino (por ejemplo a) entonces se apunta el puntero del vértice f
al vértice a, y lo mismo se hace con todos los vértices que haya en
la pila entre f y a.
•  Si un vértice no tiene más salidas y tiene un puntero que le
apunta, se quita del camino y se pone en cabeza de una lista.
•  Si el vértice que no tiene más salidas tiene más vértices que le
apuntan en la pila, entonces todos los vértices que el apuntan
junto con él mismo forman un bucle y se ponen todos en cabeza
de una lista
•  Si un vértice no tiene más salidas y no tiene un puntero que le
apunta, se quita del camino y se sigue por el anterior nodo del
camino.

54
Matriz con variables de salida asignadas. Grafo dirigido
f1 f2 f3 f4 f5
Ejemplo del algoritmo de particionado de Tarjan.

path stack lowlink
f4
1) f1 f1 f1 f2 f3 f5 f4 f1 f2 f2
5) f5 8)
f1
f2 f3
2) f1 f2 f1 f2
f1 9) f1 f2
f1
3) f1 f2 f3 f3
f4
f2 Selecciona f1 y f2
6) f1 f2 f3 f5 f5
f1
f3
Selecciona f4 y f5 f2
f4
4) f1 f2 f3 f5 f4 f5 f1
f3 Particionado y orden de resolución

f2 f3
f1
7) f1 f2 f3 f2 [f2 f1] f3 [f4 f5]
f1
Selecciona f3

55
Realizar el particionado del siguiente sistema mediante el algoritmo de Tarjan
f1 f2 f3 f4 f5 f6
f5
f4
f3
f1 f2 f3 f4 f5 f2
f1
f6
f3
[f1 f2 f3 f6] f4 f5
f1 f2 f3 f6 f2
f1
f6
f1 f2 f3 f6 f3
f2
f1

56
4 5
1 2 3 6 7
3 5
2 7
1 2 3 3
1 6
1 2 3 6 7 3 1 2 2
3 2 1
1 2 3 2 1
1 5 1 1
7
5
6
7 1 2 3 6 3
6 2
1 2 3 6 7 5
3
1 1 [2 3] 4 [6 7 5]
2
3
1
1 2 2
5
1
7
4
6
1 2 3 6 7 5 3
3
1 2 4 2
2
1
1

57
Ejemplo: cambiador contracorriente.

58
Pasos:
1. Asignar variables de entrada (alg. Christensen-Rudd)
2. Asignar conjunto de salida (alg. Duff)
3. Particionar y definir orden de resolución (alg. Tarjan)

59
1.  Asignar variables de entrada (alg. Christensen-Rudd)
k k
k
Q Q
Q
A
U
U
F1 F1
F1
f1 f1
F2 f1 F2
f2 F2 f2
F3 f2 F3
f3 F3 f3
F4 f3 F4
f4 F4 f4
f4
f5 T1 T1
T1
f6 T2 f6 T2
f6 T2
f7 T3 T3
f7 T3
T4 T4
T4
dT dT
dT
A,f5 U,f7 F4,f2;F2,f1;dt,f6

60
k
k k
Q
Q
F1
F1
F1
F3
f3 F3
F3
f3
f4
T1
T1
T1
T2
T2
T2
T3
T3
T3
T4
T4,f4 Q,f3

61

62
2. Asignar conjunto de salida (alg. Steward)
Q A U F2 F4 T4 d
T
f1 x
f2 x
f3 x
f4 x x
f5 x x x x
f6 x x
f7 x x x x

63
2. Asignar conjunto de salida (alg. Steward)
Q A U F2 F4 T4 d
T
f1 x
f2 x
f3 x
f4 x x
f5 x x x x
f6 x x
f7 x x x x

64
Q A U F2 F4 T4 d
T
f3 x
f5 x x x x
f7 x x x x
f1 x
f2 x
f4 x x
f6 x x

65
3. Particionar y definir orden de resolución (alg. Tarjan)
f3 x
f5 x x x x
f7 x x x x
f1 x
f2 x
f4 x x
f6 x x
f1 f2 f3 f4 f5 f6 f7
f3 f4 f6 f2 f1 f7 f5

66
Singularidad funcional
Singularidad funcional implica que el Jacobiano es singular para cualquier
valor de las variables, x. Ejemplo:
La causa de singularidad funcional suele ser un modelo mal construido o

con un conjunto de especificaciones mal seleccionado:
• O es redundante lo que implica un número infinito de soluciones.
• O es inconsistente lo que implica que no existe una solución.
En el ejemplo anterior
si k=-9 el sistema es redundante
si k<>-9 el sistema es inconsistente
Singularidad numérica (o local)
El Jacobiano es singular sólo para algunos valores de las variables.

67

Lee Rudd

Cargado por

Copyright:

Formatos disponibles

Lee Rudd

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Lee Rudd

Cargado por

Copyright:

Formatos disponibles

Tema 3: Simulación

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Simulación y análisis de sensibilidad

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

1.2 Arquitectura de software de un entorno de modelado y simulacion.

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Pasos para realizar una simulación.

1 Seleccionar un conjunto de unidades consistente.

2 Especificar los componentes (especies químicas) que participan en el proceso.

3 Seleccionar un paquete de propiedades físicas

4 Definir la estructura del diagrama de flujo. Qué unidades participan y los

5 Seleccionar un modelo para cada unidad de proceso.

6 Satisfacer los grados de libertad del proceso

7 Seleccionar las opciones de resolución numérica

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

1.3 Estrategias de resolución: Secuencial Modular vs. Orientado a ecuaciones

Secuencial Modular Orientado a ecuaciones

Se simula una unidad cada vez Se resuelven todas las unidades a la

Iterar en las corrientes de “rasgado” Se actualizan todas las incógnitas

Inicialización es importante Inicialización es muy importante

No necesita mucho almacenamiento Requisitos de almacenamiento

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Orientado a ecuaciones Secuencial Modular

Resolución del sistema: 3 niveles

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Un grafo G es un conjunto de vértices (V) y un conjunto de arcos que

Digrafo: es un grafo en el cual los arcos tienen una dirección. La

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Matriz de adyacencia o estructural vértices

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

conectados en dos pasos

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Poner el digrafo correspondiente

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

¿Hay bucles en el diagrama?¿Cuántos?¿Se pueden identificar de forma automática?

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

El particionado consiste en la identificación de los subsistemas

1) Seleccionar un vértice (de forma arbitraria)

Otros algoritmos de particionado que se basan en la matriz de adyacencias: algoritmo de Norman y el

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Sigue con F (sin salida) Lista: F G

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Lista: A B C E D F G Las listas NO son únicas

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

C àDàAàBàC (agrupo en un BUCLE)

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

CDAB (sin salida) Lista: [CDAB ] E F

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Ejercicio: Particionar el siguiente diagrama de flujo.

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Ejercicio: Particionar el siguiente diagrama de flujo.

1 à2 3 6 9à12 4 5à7 8à10à11

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

React. CloracionàColumna àResto

Modelado y simulación en Ingeniería Química. Manuel Rodríguez

Procedimiento de resolución secuencial modular

1)  Seleccionar un vértice (de forma arbitraria)

1)  Identificar los bucles presentes

1)  Identificar los bucles presentes

• Selecciona la fila (columna) con menor número de incidencia.

•  Se va recorriendo el grafo (se empieza por cualquier vértice) y se