Modelación y balance de masa
Felipe Sepúlveda Unda
Ingeniero Civil en Metalurgia
Doctor en Ingeniería de
Procesos de Minerales
Modelación
• Ajuste de datos.
• Hidrodinámica.
• Cinética.
• Electroquímica.
Cinética de Flotación (C. F.)
CINÉTICA FLOTACIÓN TIEMPO (τ)
EL TIEMPO DE FLOTACIÓN ES LA VARIABLE FUNDAMENTAL EN
EL DISEÑO.
CORRESPONDE AL TIEMPO MÁXIMO QUE HAY QUE DARLEA LAS
PARTÍCULAS MÁS LENTAS PARA QUE PUEDAN SER EXTRAÍDAS
DE LA PULPA.
EL TIEMPO DE RESIDENCIA ESTÁ VINCULADO AL FLUJO DE
AIRE, DE MODO TAL QUE SI ESTE FUESE PEQUEÑO, τ DEBERÍA
SER ALTO PARA COLECTAR TODAS LAS PARTÍCULAS
C. F. ES EL ESTUDIO DE LA VARIACIÓN DE CONCENTRADO QUE
REBALSA
EN LA ESPUMA CON EL TIEMPO Y LA IDENTIFICACIÓN
CUANTITATIVA DE
TODAS LAS VARIABLES QUE CONTROLAN LA VELOCIDAD.
Tipos de celdas de flotación
Discontinuas o semi continuas (Batch)
• laboratorios
Continuas
• Platas piloto e industriales
La flotación es un proceso que se caracteriza
por el cambio.
En las celdas, a medida que transcurre el
tiempo va bajando la cantidad de material
flotable y se modifican las variables físico
químicas que caracterizan la pulpa.
De allí que la forma más correcta para
modelar este proceso sea mediante una
aproximación de tipo cinética.
En las operaciones tipo Batch o semi batch
que se llevan a cabo en laboratorios, al
estudiar como van evolucionando al
transcurrir el tiempo, tanto la recuperación
como la ley de concentrado se pueden
obtener las curvas típicas.
Parámetros de control
• Recuperación: eficiencia de separación
en un proceso de concentración
(flotación)
Parámetros de control
Proceso de flotación
Cola o relave
- Flujo
Flotación - ley
Alimentación
- Flujo
- ley
Concentrado
- Flujo
- ley
Parámetros de control
Recovery = Flowrate of metal in concentrate
Flowrate of metal in feed
Parámetros de control
• Ley (%), grado de pureza de un producto
de la flotación, en relación a una especie
o elemento de interés
Grade = Mass metal in stream
Total mass of solids in stream
Mendez et al., 2009, Modeling of grinding and
classification ircuits as applied to the design of flotation
processes, computers and Chemical Engineering, 33, 97-
111.
Ley de concentrado v/s Tiempo en una prueba
discontinua de flotación
Recuperación versus tiempo en una prueba discontinua
de flotación
Modelo Fenomenológicos de Flotación Discontinua.
La forma de estas curvas sugiere que el cambio de la variable respuesta con
respecto al tiempo siguen una ley de tipo exponencial y para describirla se
han utilizados modelos fenomenologicos, empíricos y probabilísticos.
El modelo más usado es el que establece analogías con la ecuación
Fundamental de la cinética química, éste, en su forma más general, Puede
escribirse mediante la siguiente ecuación:
K: Constante cinética
n
dCi
dt
= −k C
i =1
i
ni Ec 1
Ni : representa el orden de reacción respectivo a cada una de las
concentraciones consideradas
Ci : representa una serie de concentraciones que pueden tener importancia
en el sistema, p. e. concentración de la especie flotable, concentración de
burbuja, reactivos, etc.
dC/dt : Variación de concentración de la especie flotable en la celda
La ecuación anterior se simplifica si se considera que la
velocidad de la
flotación dependiese solo de la concentración de una
especie flotable en la celda
dC Ec 2
=−k C n
C : Concentración de especie flotable
dt n : orden de la reacción
K : Constante de flotación
Suponiendo una cinética de primer orden:
− kt Ec 3
C = c0 e
C0 : Concentración de material flotable al tiempo cero
C : Concentración material flotable al tiempo t
La ecuación anterior se puede graficar de la
siguiente manera:
c0 Ec 4
ln = kt
C
Ln c0 / C
t
En la práctica no es posible recuperar el 100% por lo que se
introduce el concepto de recuperación máxima, RMÁX, ya que
no todo el material es flotable.
El término del logaritmo natural debe ser corregido, pues c0 es la
concentración de todo el material flotable menos la concentración
del material que no flotó y quedó.
c0 − c Ec 5
ln = kt
C − C
Es posible escribir :
c0 − c c0 − c
R = Ec 6 Rt = Ec 7
c0 c0
Luego la ecuación queda:
= − k (c − c )
dC n Ec 8
dt
= − k c0 (R − R )
dR n −1 n
Ec 9
dt
Denominando a la expresión k Con-1 como kn se
obtiene la sgte ecuación:
= − k n (R − R )
dR n
Ec 10
dt
Para n = 1: Rt = R 1 − e ( − kt
) Ec 11
Rt = Recuperación parcial a cualquier tiempo t
R = Recuperación máxima posible
MODELO DE 1ER ORDEN GARCÍA ZÚÑIGA
Este modelo corresponde a un caso particular de la ecuación cuando
n=1 y recibe dicho nombre en honor al Héctor García Zúñiga, ingeniero
chileno quien la propueso en el año 1935
Él realiza una descripción matemática y reconoce la analogía
algebraica con las ecuaciones de reacción de cinética de química de
primer orden.
En forma posterior, Schumann por el año 1942, puntualiza que la
constante específica de velocidad de flotación tienen el mismo
significado dimensional que la constante de reacción específica de una
reacción cinética de primer orden.
La ecuación de García zúñiga permite ajustar la gran
mayoría de los
resultados experimentales y calcular la constante de
velocidad de flotación.
(
Rt = R 1 − e − kt
) Ec 12
O bien en términos de la concentración de la
especie útil
(c0 − c )
ln = k t Ec 13
(c − c )
Trabajando la ecuación de las recuperaciones
obtenemos:
Rt − kt Rt − kt
=1− e −1= − e
R R
Rt − kt R − Rt − kt
1− =e =e
R R
(R − Rt )
ln = − kt
R
Si se grafica la última ecuación:
(R − Rt )
ln
R
-k
Tiempo de flotación
Las Rt son recuperaciones acumulativas; la R y t (datos
experimentales), sirven para obtener por medio de graficas el
valor de k.
[k] = t-1 ejemplo [k] = seg-1
Actualmente por medio de la minimización del error, uno logra
obtener el mismo resultado.
MODELO DE AGAR
EN LA PRÁCTICA LA ECUACIÓN DE GARCÍA
ZÚÑIGA ES MUY DIFÍCIL QUE SE CUMPLA
DEBIDO A LA DIFICULTAD DE AJUSTAR EL
TIEMPO CERO DE MANERA PRECISA, POR LO
CUAL A FIN DE AJUSTAR LOS DATOS
EXPERIMENTALES SE HA PROPUESTO LA
CORRECCIÓN SGTE:
(
Rt = R 1 − e − k (t + )
) Ec 14
DONDE θ ES UN PARÁMETRO QUE PERMITE QUE LA
CURVA PASE
POR EL ORIGEN
MODELO DE SEGUNDO ORDEN
En el caso que se considere cinética de segundo
orden, con n=2, el resultado que se obtiene al
integrar la ecuación es
2
R kt Ec 15
Rt =
1 + R k t
En este caso k es una constante de segundo orden
Esta ecuación es de un uso más restringido
con respecto a la de García Zúñiga pero en
algunos casos puede relacionar
adecuadamente los datos experimentales.
MODELO DE KLIMPELL
Este modelo compara los resultados de
laboratorio para los
diferentes colectores.
1
R = R 1 − 1 − e − kt
Ec 16
kt
MODELOS PARÁMETROS DISTRIBUIDO
Consideraciones
García Zúñiga propone que todas las partículas
presentan idénticas
Propiedades de flotación. Pero se ha demostrado que
esto no siempre
es así ya que la constante de velocidad depende de:
Tamaño partícula
Condiciones hidrodinámicas de la
celda.
Mineralogía
Oxidación superficial.
Por lo cual una aproximación cinética más exacta
debe de considerar a la constante de velocidad como
un parámetro distribuido
Modelo de Kelsall (1962)
Algunos científicos considerar que la especie flotable se pueden
dividir en dos grandes grupos:
1.- Baja velocidad de flotación
2.- Alta velocidad de flotación
O lo que también se puede dividir en:
a.- Mineral útil que no flota.
b.- Mineral útil de débil flotabilidad (ks)
c.- Mineral útil de alta flotabilidad. (kf)
De acuerdo a lo anterior la expresión cinética de
primer orden puede escribirse como
C = C e − ks t
+ (1 − ) e
−k f t
Ec 17
Donde Φ es la fracción de material de baja velocidad
existente en la pulpa. En términos de la ecuperación
este modelo queda descrito por la sgte ecuación:
(
R = R 1 − e
−k f t
)(1 − ) + (1 − e )
−ks t Ec 18
(
R = R 1 − e
−k f t
)(1 − ) + (1 − e )
−ks t
Otra forma de expresar la constante de velocidad
como un parámetro distribuido es mediante la
ecuación general
C = C e − kt
f (k ,0 ) dk Ec 19
0
Donde f(k,0) representa una distribución continua de
velocidad.
Inovation in flotation
technology
Paul marcos, Kostas A.
matis
Kluwer Academic
Publishers
1992.
CONDICIONES DE OPERACIÓN
MOLIENDA
muestra : 1 kg 100% bajo 10 mallas Tyler
tiempo : 11 minutos 58.8% -200 #Tyler
%sólido : 66.7
FLOTACIÓN
tiempos : 1, 3, 6, 10, 16, 20, 25, 30 minutos
acondición : 2 minutos
colectores : 6.8 [g/ton] SF-323
6.8 [g/ton] A-238
espumante : 9.07 [g/ton] M.I.B.C.
depresor : CaO libre
pH : 10
paleteo : C/10 segundos
PESOS (grs.) LEYES (% ) CuT FINOS (% ) RECUPERACIÓN
TIEMPO (min)
parcial acumulado parcial acumulado parcial acumulado parcial acumulado
1 49.5 49.5 18.36 18.36 9.088 9.088 67.49 67.49
3 29.5 79 9.71 15.13 2.864 11.953 21.27 88.76
6 23.1 102.1 2.47 12.27 0.571 12.523 4.24 93.00
10 27.8 129.9 0.85 9.82 0.236 12.760 1.75 94.75
16 22.9 152.8 0.38 8.41 0.087 12.847 0.65 95.40
20 22.5 175.3 0.15 7.35 0.034 12.880 0.25 95.65
25 32.3 207.6 0.1 6.22 0.032 12.913 0.24 95.89
30 27.5 235.1 0.09 5.50 0.025 12.937 0.18 96.08
Relave 765.8 1000.9 0.07 1.35 0.528 13.465
Cabeza Calc. 1000.9 1.345
PESOS (grs.) LEYES (% ) CuT FINOS (% ) RECUPERACIÓN
TIEMPO (seg.)
parcial acumulado parcial acumulado parcial acumulado parcial acumulado
10 21.5 21.5 26.57 26.57 5.713 5.713 22.45 22.45
20 23.8 45.3 24.56 25.51 5.845 11.558 22.97 45.42
50 24.5 69.8 19.93 23.55 4.883 16.441 19.19 64.60
90 27.5 97.3 16.21 21.48 4.458 20.898 17.52 82.12
120 18.5 115.8 11.12 19.82 2.057 22.956 8.08 90.21
240 15.7 131.5 8.42 18.46 1.322 24.278 5.19 95.40
360 13.9 145.4 4.13 17.09 0.574 24.852 2.26 97.66
Relave Cl. 73.6 219 0.81 11.62 0.596 25.448
Relave Ro. 1780.5 1999.5 0.087 1.350 1.549 26.997
Cabeza Calc. 1999.5 1.350
CINÉTICA CLEANER
30 100
(%) RECUPERACIÓN
25
(%) LEYES (CuT)
80
ley parcial
20
60 ley acumulada
15
40 % recuperación
10
5 20
0 0
10 20 50 90 120 240 360
TIEMPO (seg.)
RESULTADOS. -
Tiempo Cleaner : 101.57 segundos (1.70 min.)
Concentrado Cleaner : 20.84 (%) CuT
(%) Recuperación : 85.24 (%)
Tipos de celdas de flotación
Discontinuas o semi continuas (Batch)
• laboratorios
Continuas
• Platas piloto e industriales
Bank of Flotation Cells
TAILS 1 TAILS 2 TAILS n-1
= FEED 2 = FEED 3 = FEED n
FEED 1
FINAL TAILS
CON 1 CON 2 CON 3
CON n
CUMULATIVE CONCENTRATE AFTER 3 STAGES
CUMULATIVE CONCENTRATE AFTER n STAGES
Concentrate grade decreases down the bank
• Mineral recovery is a rate process – C=C0e-kt
• Gangue recovery is a physical process – depends on froth structure
Flotation rates – previously
Kinetics and flotation rates – without the froth
dC
= − kC First order kinetics...
dt
k
R= Mineral recovery from a single tank
k + 1 τ is the residence time (tank volume/flowrate)
This does not give a full mass balance:
Water? Hydrophilic gangue?
k1
Modelo reversible A⇒B
− kt
Rt = 1 − e
𝑗 1
𝑇𝑖 =1−
𝑗 𝑗 𝑁𝑗
1+ 𝑘𝑖 𝜏
k1
Modelo reversible A⇄B
k2
1+𝐾𝑒𝑞
−𝑘 1 𝑖 𝑖 𝑡
𝐾𝑒𝑞
𝑇𝑖 = 𝑇𝑒𝑞 𝑖 1 − 𝑒𝑥𝑝 𝑖
𝑗 𝑗 𝑁𝑗
𝑗 1 + 𝑘2 𝑖 𝜏
𝑇𝑖 = 1− 𝑗 𝑗 𝑗
1 + 𝑘1 𝑖 + 𝑘2 𝑖 𝜏
k1
Modelo Hibrido A⇒B, 𝑹∞ coef de ajuste
(
Rt = R 1 − e − kt
)
𝑗 1
𝑇𝑖 = 𝑅∞ 1 −
𝑗 𝑗 𝑁𝑗
1+ 𝑘𝑖 𝜏
𝑗 1
𝑇𝑖 = 𝑅∞ 1 −
𝑗 𝑗 𝑁𝑗
1+ 𝑘𝑖 𝜏
1 2 3 4
II Semestre 2014
Numero de celdas
Reactor de mezcla
perfecta
Yianatos J.B., Bucarey R. Lareana J. Henríquez F. Torres L., 2005, Collection
zone kinetic model for industrial flotation columns, Mineral Engineering, 18,
Wash water
1373 – 1377.
Concentrate τs
Feed
τl
Air
τs
L
Tail
1 1 𝜏𝑙 𝑘𝑚 𝜏𝑙 + 1
𝑅 = 𝑅∞ 1 − + 𝑙𝑛 −1
𝑘𝑚 𝜏𝑙 − 𝜏𝑠 𝑘𝑚 𝜏𝑠 + 1 𝜏𝑙 − 𝜏𝑠 𝑘𝑚 𝜏𝑠 + 1
II Semestre 2014
¿Que se puede comentar acá?
¿Que se puede comentar acá?
¿Que se puede comentar acá?
Circuitos de flotación
ANÁLISIS DE GRÁFICOS
Todas las variables involucradas en el proceso de flotación podemos
representarlas gráficamente asociándolas de dos o más.
Estas gráficas nos darán una pauta de cómo se está comportando
metalúrgicamente el sistema, las diferencias que puedan existir entre
distintos materiales, reactivos que se usan, parámetros involucrados, etc.
Fuertemente Hidrofóbico
% Recuperación
Gráfico Recuperación vs
tamaño de partícula para
materiales con distinto
Débilmente carácter hidrofóbico.
hidrofóbico
No hidrofóbico
Tamaño de
partícula
% Recuperación
Rougher
Limpieza
Tamaño medio de partícula (um)
Gráfico recuperación vs tamaño de partícula para distintos tipos de circuitos.
Análisis de la Superioridad del Circuito
Se trata de verificar si los circuitos son siempre
superiores a las etapas y comparar circuitos entre si
LEY
óptimo
X
X Multietapa
Una sola etapa
Recuperación
OPERACIÒN DE CIRCUITOS DE FLOTACIÒN
Circuito Rougher ( primario)
Recuperadores
Circuito Scavenger (barrido)
anterior al rougher
Circuitos alimentación por cola o relave anterior
Cleaner ( limpieza)
Limpiadores
Recleaner ( limpieza)
Operaciones complementarias:
- Molienda -
Clasificación
- Espesamiento - Filtración
Alimentación
Relave
ROUGHER SCAVENGER
RECUPERACIÓN
CLEANER LEY
Concentrado
PARÁMETROS BÁSICOS
BALANCE DE MASA :
F=C+T
BALANCE DE FINOS :
F*f=C*c + T*t
BALANCE DE AGUA :
Wf = WC + WT
PARÁMETROS BÁSICOS
RECUPERACIÓN (R)
Es el porcentaje de materia útil de mineral que es transferido por las
operaciones de concentración desde la alimentación del proceso al
concentrado
Donde :
C *c C : Masa de Concentrado
%R = * 100 c : Ley de Concentrado
F *f
F : Masa de Mineral Alimentado al Proceso
f: : Ley de Cabeza
También se puede definir por las leyes
c * (f − t )
Donde :
t : Relave
%R = * 100
f * (c − t )
LÍMITE DE CONCENTRACIÓN DE MINERALES
En la etapa de concentración de minerales, la máxima ley de concentración
que se puede lograr es la correspondiente a la de la especie mineralógica
PURA
En la realidad es casi imposible alcanzar dichas leyes ya que se estaría
realizando una concentración PERFECTA lo que implica un elevado costo
de la operación
Alcanzar una elevada ley no es muy recomendable debido a que influye
en la RECUPERACIÓN.
MIENTRAS MAYOR LEY MENOR
RECUPERACIÓN
Análisis de circuitos de Separación de Minerales
En la industria rara vez se usa una etapa individual de separación, sino que se
utilizan arreglos de dos o más etapas. En los procesos de flotación se usan
arreglos R-C-S, es decir, primarios (rougher) – limpieza (cleaner) – barrido
(scavenger).
Consideremos la siguiente configuración de circuito.
R S
F RF
CF
Esta figura se transforma en una configuración de nodos, donde cada nodo
representa una etapa de separación por flotación y las flechas con sentido a la
izquierda representan el flujo de concentrado de cada etapa de separación y con
sentido a la derecha los flujos de relave.
F
1 2R (1-R)
R
(1-R)
C S
2
(1-R)
2
R
R (1-R)
R (1-R)
CF RF
Cada nodo representa una unidad, equipo o banco de celdas.
Se considera, el caso simple de tener la misma recuperación R en cada una
de las etapas (nodos) y suponemos que la corriente de entrada a la etapa
primaria es unitaria. La recuperación de una etapa individual es R
2
La recuperación global del circuito R' es: R
R' =
1 - 2R(1 - R)
Ejemplo
Separación Cromita / Sílice, según los siguientes datos:
Ley de cabeza: 68,9 % Cr2 O3 , 31,1% Si O2
R Cr2O3 = 0.90 % (recuperación de cromita en cada etapa del circuito R-C-S).
R SiO2 = 0.15 % (recuperación de sílice en cada etapa del circuito R-C-S).
2
(0,90)
R'Cr2O 3 = * 100 = 98,8 %
1 - (2 * 0,15 * (1 - 0,90))
2
(0,15)
R'SiO2 = = 3,02 %
1 - (2 * 0,15* (1 - 0,15)
R Cr2O 3 * Ley f Cr
Ley Cr2O 3 conc. =
2O 3
(R'Cr2O 3 * Ley f Cr
2O 3
(
) + R'SiO 2 * Ley f SiO 2 )
Reemplazando en la ecuación anterior tenemos lo siguiente:
0,988 * 68,9
Ley Cr2O3 conc. = = 98,7%
( 0,988 * 68,9 ) + ( 0,0302 * 0,311)
Por lo tanto tenemos las siguientes ventajas del circuito respecto a una
etapa unitaria:
R' = 98.8 % v / s R = 90 %
Ley Cr2O3 conc. = 98.7 % v/s Ley Cr2O3 = 93 %
0,9 * 68,9
Ley conc. rougher = = 93%
( 0,9 * 68,9 ) + ( 0,15 * 31,1)
BALANCE DE MASAS
EN PROCESAMIENTO DE MINERALES.
Un aspecto de una elevada relevancia en el control, evaluación y análisis
de plantas de procesamiento de minerales es el BALANCE DE
MATERIALES.
En faena existen una variedad de tareas operacionales, para las cuales es
recomendable determinar los flujos másicos totales de las corrientes que
participan en el circuito y así compatibilizar los balances para cada uno de
los componentes químicos o de las propiedades físicas que caracterizan a
dichas corrientes.
Para lo anterior, y como una tarea rutinaria, se toman muestras de los diferentes
flujos involucrados y con dicha información se establecen los balances de masa
correspondientes. Sin embargo, las mediciones que se obtienen desde un muestreo
contienen errores aleatorios y, al resolver los balances, se generan,
inevitablemente, algunas inconsistencias que pueden distorsionar la evaluación que
se haga de la marcha operacional de la empresa o de alguna de sus secciones y
equipos en particular.
ECUACIÓN DE BALANCE DE MASA
ENTRADAS - SALIDAS = ACUMULACIÓN ó DESAPARICIÓN
SI ACUMULACIÓN ó DESAPARICIÓN ≠ 0 ESTADO NO ESTACIONARIO
COMPORTAMIENTO DEPENDE DEL TIEMPO
SI ACUMULACIÓN ó DESAPARICIÓN = 0 ESTADO ESTACIONARIO
La masa que entra es igual a la que sale
Para este curso se tomara sólo sistemas de nodos en estado estacionario
Con lo cual
ENTRADAS - SALIDAS = 0
CONCEPTOS DE NODOS Y CORRIENTES.
Para una mejor representación del circuito es conveniente representar los diversos
equipos involucrados por medio de una RED DE NODOS, los cuales estarán
comunicados por medio de las corrientes del material
Un nodo corresponde a una ubicación específica dentro del circuito, por ejemplo un
cajón de distribución, un hidrociclón o una celda de flotación
Nodo Simple: en ellos sólo participan 3 corrientes una de entrada y dos de salidas
o viceversa
F2
F1
F1
Separador UNIÓN
F3
F3
F2
Gráfico de Nodos Simples, separador y unión.
Nodos Complejos: participan en ellos más de dos corrientes en la salida o en la entrada
BALANCE DE MASA EN UN NODO
Los balances que se pueden platear suelen distinguirse dos categorías
1º Balances macroscópicos: Se refiere a alcanzar la cuadratura de masa
total, sea ésta en forma líquida, sólida o de pulpa
2º Balances Microscópicos: Se refiere a alcanzar la cuadratura con
respecto a un componente específico como compuesto químico,
granulometría, especies mineralógicas, etc. (Balances Microscópico)
Es importante mencionar que existirán oportunidades que coincidirán el número de
variables con el número de ecuaciones que se podrán plantear.
Pero en otras ocasiones habrán más ecuaciones planteadas que número de variables
por lo que el sistema tendrá más de una solución por lo que se deberá usar métodos
que permitan obtener la solución óptima.
A modo de ejemplo si se quisiera realizar un simple balance en el nodo que a
continuación se muestra
F2
F1
Separador
F3
Conociendo los valores de F1 y F2 el valor de F3 puede calcularse por simple
diferencia y así lograr la cuadratura. Pero si se conociesen los 3 flujos lo más
probable que suceda es que exista una pequeña diferencia la que no permite lograr la
cuadratura del nodo. Esta diferencia de flujos la podemos representar de la siguiente
manera:
F1 - F2 - F3 = Δ
En este caso para encontrar la cuadratura se deberá aplicar técnicas de optimización.
BALANCE DE MASA EN CIRCUITOS DE PROCESAMIENTOS
Matriz de Proceso
Todo proceso por complejo que sea puede representarse
como una red de nodos unidos por corrientes de mineral.
A modo de ejemplo, el circuito de la figura es un esquema
típico de una remolienda de concentrado, el concentrado
proveniente de la flotación primaria llega a un pozo de
bombas y es impulsado hacia el hidrociclón, en este equipo el
material se separa en dos fracciones, la gruesa que retorna al
pozo de bombas, y la fina que alimenta una celda de limpieza
de columna.
F4
1
F1 F2 F7
5
F3 W1 4 F8
W2
F5
2 3
F6
F9
F7
F2
F4
5
W1 4 F8
F1
F5
1
F6
F3 3 W2
2
F9
Red de nodos del circuito ejemplo