Equilibrio Dinámico en Redes de Transporte Con Múltiples Orígenes y Destinos

UNIVERSIDAD DE CHILE
FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS

DEPARTAMENTO DE INGENIERÍA MATEMÁTICA
EQUILIBRIO DINÁMICO EN REDES DE TRANSPORTE CON MÚLTIPLES

ORÍGENES Y DESTINOS
TESIS PARA OPTAR AL GRADO DE MAGÍSTER EN CIENCIAS DE LA

INGENIERÍA, MENCIÓN MATEMÁTICAS APLICADAS
MEMORIA PARA OPTAR AL TÍTULO DE INGENIERO CIVIL MATEMÁTICO
FELIPE IGNACIO GARRIDO LUCERO
PROFESOR GUÍA:
ROBERTO M. COMINETTI COTTI-COMETTI
MIEMBROS DE LA COMISIÓN:
JOSÉ R. CORREA HAEUSSLER
JORGE G. AMAYA ARRIAGADA
Este trabajo ha sido parcialmente financiado por Centro de Modelamiento Matemático y

Fondecyt1130564
SANTIAGO DE CHILE
2017
ii
RESUMEN DE LA MEMORIA PARA OPTAR AL GRADO DE
MAGÍSTER EN CIENCIAS DE LA INGENIERÍA,
MENCIÓN MATEMÁTICAS APLICADAS,
AL TÍTULO DE INGENIERO CIVIL MATEMÁTICO
POR: FELIPE IGNACIO GARRIDO LUCERO
FECHA: 2017
PROF. GUÍA: ROBERTO M. COMINETTI COTTI-COMETTI
EQUILIBRIO DINÁMICO EN REDES DE TRANSPORTE CON MÚLTIPLES

ORÍGENES Y DESTINOS
En el presente trabajo se estudia un modelo teórico para el equilibrio dinámico de los flujos
en una red de transporte. Específicamente, se estudian las redes de colas fluidas, modelo que
permite describir la evolución temporal de los flujos que circulan dentro de esta. Este tipo de
redes han sido utilizadas como un modelo simplificado para describir situaciones tales como
el trafico en una red urbana, el movimiento de pasajeros en un sistema de transporte publico,
el flujo vehicular saliendo de Santiago el primer día de vacaciones a las distintas regiones del
país, el flujo de datos en una red de comunicaciones, o la evacuación de un estadio al término
de un espectáculo deportivo. En dichas redes, cada arco se modela como una cola fluida de
capacidad finita seguida de un segmento con tiempo de viaje constante.
El equilibrio dinámico describe la evolución temporal de los flujos bajo el supuesto de que
cada partícula se mueve desde su origen hacia su destino utilizando un camino de tiempo
mínimo, considerando tanto el tiempo fijo requerido para atravesar cada arco como el tiempo
de espera en la cola presente al momento de ingresar a este.
Dentro del trabajo se estudiaron dos modelos generales para una red de transporte. En el
Capitulo 2 se extienden los resultados recientes para un único par origen-destino, al caso en
que existen flujos concurrentes desde un origen común hacia múltiples destinos. Para ello, se
describe el modelo de la red de transporte, se discuten las distintas formas de caracterizar
un equilibrio dinámico, y se extiende la técnica de reconstrucción del equilibrio a partir de
sus derivadas, las cuales se caracterizan mediante la noción de flujos finos normalizados. Se
presentan dos algoritmos, el primero que extiende el flujo de modo que siga siendo un flujo
en equilibrio, y otro que muestra como obtener los flujos finos normalizados. Estos resultados
extienden y completan los trabajos recientes de Koch y Skutella [13] así como Cominetti,
Correa y Larre [3], los cuales consideran el caso de un único par origen-destino.
En el Capítulo 3 se aborda el caso más complejo en que coexisten varios orígenes y uno
o más destinos. Tras extender el modelo teórico se discuten las complicaciones que surgen
al intentar extender el concepto de flujo fino normalizado. Para ello se presentan algunos
ejemplos que ponen de manifiesto la perdida de linealidad en los sistemas de ecuaciones
necesarios para el desarrollo de un algoritmo, y las dificultades para analizar las interferencias
temporales entre flujos provenientes desde distintos orígenes.
iii
iv
Le dedico esto a toda la gente linda que conocí en el viaje, quienes me ayudaron a
reencontrar mi alegría.
v
vi
Agradecimientos
Agradezco a Rebeca, quien me ha ayudado y entregado su apoyo desde el primer día, a

quien le deseo lo mejor en sus proyectos, que me ha demostrado que las cosas llegan, y que
siempre se debe querer lo que uno hace, estar orgulloso del trabajo de uno y entender que es
importante. A Hernán, un increíble ejemplo a seguir, desde su humilde origen, siempre con
una sonrisa en la cara y una broma precisa, me ha enseñado que la vida se lleva mejor si
uno está contento, que sin importar lo agotadores que sean los días, siempre quedan energías
para seguir si uno está feliz haciendo lo que hace. A Francisca, una luchadora de la vida,
de quien aprendí la importancia de los amigos y que la elección correcta siempre será pasar
tiempo con ellos. Orgulloso queda pequeño al lado de lo que pienso de ella, quien desde
siempre ha querido cambiar el mundo y luchar contra las injusticias que nuestra sociedad
presenta. Se que lo podrás lograr, se que aportarás mucho más que un simple grano de arena.
A José, quien me ha demostrado que es importante siempre ser uno, no dejar de ser original
a pesar de las circunstancias en que se encuentre y que se debe vivir relajado, que no se le
debe dar importancia a las cosas que no lo requieren. A Seba, Diego, Rorro, Cata, Caro,
Phillippe, con quienes he crecido todos estos años desde cuarto medio y hemos compartido
las distintas transiciones que uno tiene en la U. Muchísimo ha pasado desde que compartimos
esas pizzas en el forestal y la Carito se comía su súper ensalada de lechuga. Hemos cambiado
de peinados, ropa, ideologías, la forma de alimentarnos y la forma de expresarnos, pero cada
vez que me junto con cada uno de ustedes, es como volver a ese parque y poder conversar
las mismas cosas de siempre. A mis amigos de primer año, con quienes vivimos todas las
primeras cosas que uno se encuentra en Beauchef, con quienes dimos los controles mechones,
jugamos taca taca y pool, comimos chaparras, fuimos a Cartagua y Coñaripe. A Dim in da
house, con quienes hemos vivido los últimos años y hemos compartido el final del proceso.
Espero que aunque nuestros caminos se separen ahora, algunos con postgrados, otros con
trabajo, siempre hayan oportunidades de compartir un almuerzo, una cerveza helada, de
jugar unos cachos y de reírnos con las historias del Pancho o las irreverencias de Pablo. Les
deseo lo mejor cabros. A la Javi y Abner, los mejores compañeros de oficina, con quienes
luchamos hasta el fin por tener el pouf, usamos pantuflas, escuchamos música, capturamos
pokemones y hablamos horas. Al profesor Cominetti, gracias por todo el apoyo entregado
este último año, por confiar y aceptar ser mi profe guía, además de darse el tiempo cada
semana para trabajar y darme siempre buenas ideas en la matraca, además de la ayuda en
las postulaciones. Finalmente, a quien fue la persona más importante para mi durante estos
años universitarios, Andrea, mi gran y mejor amiga, quien me enseñó lo linda que puede ser
la vida, lo feliz que uno puede llegar a ser con otra persona, que no se deben organizar los días
cuadradamente, que las mejores risas son aquellas que aparecen de la nada, que las mejores
conversaciones son a las tres de la mañana, y que me entregó su más honesto y puro amor.
vii
viii
Tabla de Contenido
Introducción 1
1. Preliminares. 6
2. Modelo de cola fluida, caso de un origen y múltiples destinos. 9

2.1. Modelo de red de flujo dinámico. . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1.1. Flujo dinámico. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1.2. Cola fluida. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1.3. Tiempo de viaje por un arco. . . . . . . . . . . . . . . . . . . . . . . 11
2.1.4. Existencia y unicidad de descomposición de un flujo de salida. . . . . 12
2.1.5. Caminos dinámicos mínimos. . . . . . . . . . . . . . . . . . . . . . . 14
2.1.6. Equilibrio Dinámico. . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.7. Flujo acumulativo y Flujo fino normalizado . . . . . . . . . . . . . . 17
2.2. Algoritmo de extensión de flujos en equilibrio dinámico. . . . . . . . . . . . . 22
2.2.1. Algoritmo de extensión. . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.2.2. Correctitud del algoritmo de extensión de flujo en equilibrio dinámico. 24
2.3. Cálculo de un flujo fino normalizado. . . . . . . . . . . . . . . . . . . . . . . 30
3. Modelo de cola fluida, caso múltiples orígenes y destinos. 35

3.1. Modelo de red de flujo dinámico. . . . . . . . . . . . . . . . . . . . . . . . . 35
3.1.1. Redes de flujo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.1.2. Equilibrio dinámico. . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.2. Complicaciones en el modelo teórico. . . . . . . . . . . . . . . . . . . . . . . 39
3.3. Ejemplos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.3.1. Ejemplo 1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.3.2. Ejemplo 2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Conclusión 45
Bibliografía 48
ix
x
Introducción
Una red de flujos dinámicos es un modelo matemático abstracto que permite describir
la evolución temporal de los flujos que circulan dentro de una red de transporte. Como
motivación se puede mencionar la circulación de vehículos en la red de calles y avenidas de
una ciudad, el movimiento de pasajeros en un sistema de transporte público, los flujos de
evacuación de un estadio al término de un evento deportivo, el flujo vehicular saliendo de
Santiago el primer día de vacaciones a las distintas regiones del país, la transmisión de datos
en una red de telecomunicaciones, o el movimiento de partes y piezas dentro de una cadena
de montaje.
A diferencia de los últimos dos ejemplos en que la dinámica queda determinada básica-
mente por el diseño del sistema, en las cuatro primeras situaciones las decisiones de ruteo son
tomadas de forma descentralizada por los distintos individuos, cada uno buscando minimizar
su tiempo de traslado desde su origen al destino tomando en cuenta la congestión y tiempos
de viaje que se observarán al momento de atravesar cada arco de la ruta escogida. Tales
situaciones de interacción estratégica entre múltiples agentes pueden ser modeladas usando
los conceptos y modelos de la teoría de juegos. En el contexto de sistemas de transporte,
caracterizados por la presencia de una gran cantidad de agentes y donde cada individuo tiene
un efecto marginal sobre la congestión, el enfoque clásico considera los flujos como variables
continuas y utiliza el concepto de Wardrop para describir la situación de equilibrio de los
flujos.
En términos matemáticos, los flujos en una red se describen mediante funciones

f : [0, ∞) → [0, ∞)
en que f (θ) representa la tasa de llegada en un cierto punto de la red en función del tiempo
θ ≥ 0. Así, una red de flujos dinámicos se modela como un grafo dirigido G = (V, E) en el
cual cada arco e = vw ∈ E viene caracterizado por una aplicación que transforma el flujo
fe+ (·) que ingresa al arco e a través de su nodo inicial v, en una tasa de salida fe− (·) que deja
el arco a través de su nodo terminal w. En particular, en esta tesis se estudian las redes de
colas fluida en las que cada arco se modela como una cola fluida de capacidad νe > 0 seguida
de un segmento de ruta que tiene un tiempo de viaje constante τe .
νe
τe
Cola Link
Figura 1: Arco
1
Esto quiere decir que, si una partícula infinitesimal de flujo llega a un arco e, debe esperar
en una cola que se vacía a tasa νe y, tras haber dejado esta, debe viajar τe unidades de tiempo
para alcanzar la cabeza del arco.
Ejemplo Considere el siguiente grafo con nodos orígenes s1 , s2 y nodo destino t:
ν1 = 1 r ν4 = 2
u1 = 1 s1 t
τ1 = 2 τ4 = 1
τ2 = 2
ν2 = 1 τ3 = 5
s2
ν3 = 1
u2 = 2
Figura 2: Red ejemplo
Las tasas de flujo entrante son constantes con valores u1 (θ) ≡ 1, u2 (θ) ≡ 2 respectivamente
para s1 , s2 . El origen s1 posee solo un camino para enviar su flujo hasta el destino t, de modo
que siempre utilizará los arcos 1 y 4. Por otro lado, s2 posee dos opciones, enviar su flujo
directamente por el arco 3, o bien pasar por r y compartir el arco 4 con el flujo proveniente
de s1 . En un comienzo, para s2 el camino compartido tiene menor costo, pues no existen
colas y todo se reduce a estudiar las latencias, así que escoge enviar todo su flujo por ahí.
Es importante notar que si un arco e posee capacidad νe , desde este nunca saldrá más de νe
unidades de flujo. Como s1 envía una unidad de flujo, nunca se formará cola en el primer arco,
además, en el arco 4, no se verá flujo hasta los 2 segundos que es cuando llega el proveniente
de cada fuente. Sin embargo, como llega una unidad desde cada origen, tampoco se formará
cola. De este modo, inicialmente el único arco que formará cola será el segundo. Como a este
entra flujo a tasa 2, la cola crecerá linealmente según z2 (θ) = θ. Debido al crecimiento de
esta, llegará un instante en que el costo de enviar flujo por el camino compartido se igualará
con el costo de utilizar el arco 3. Efectivamente, una partícula que ingresa a la red desde s2
en el instante θ debe esperar en la cola existente y luego viajar por ambos arcos, de modo
que su tiempo total de viaje T viene dado por
z2 (θ)
T (θ) = + τ2 + τ4 = θ + 2 + 1 = 3 + θ
ν2
En θ = 2, el tiempo de viaje por ambos caminos se iguala, y ya no resulta óptimo para s2

seguir enviando todo su flujo por el camino compartido. Para que el costo de viajar por ambas
rutas sea el mismo a partir de θ = 2, y las partículas infinitesimales no tengan incentivo a
cambiar de arcos, el nodo origen divide su flujo y comienza a enviar una unidad por cada
camino, con lo cual las colas se paralizan y el tiempo de viaje se congela en T (θ) = 5. Con
esto se alcanza un estado estacionario en el cual todas las colas y los tiempos de viaje se
congelan, de modo que no ocurren más cambios en la red.
2
Resultados previos.
Se comienza esta sección describiendo los equilibrios de Wardrop estáticos, enfoque clásico
utilizado en estos modelos, para luego pasar a ver sus limitaciones.
Considere un grafo dirigido G = (V, E), y un conjunto finito de jugadores J, donde cada
i ∈ J busca viajar desde un origen si ∈ V a un destino ti ∈ V . Cada jugador i selecciona
un subconjunto Ei de arcos de E para armar su ruta a través del grafo. Asociado a la ruta
elegida, el jugador i debe incurrir en un costo de viaje, dado por el tiempo total que necesita
para cruzar la red X
Ti = ce (ue )
e∈Ei
en donde cada ce (ue ) es el tiempo del arco e ∈ E como función creciente del flujo ue en el
arco. Es importante notar que un aumento de la congestión ue tiene externalidades negativas
pues aumenta los costos de todos los jugadores que utilizan el arco e. Finalmente, considere
que el número de jugadores |J| tiende a infinito, lo que da origen a un juego de ruteo no
atómico.
Por simplicidad, suponga que existe solo un origen s y un destino t, con demanda total
u > 0. Sea R el conjunto de todos los caminos que conectan (s, t) y xr ≥ 0 la fracción de
flujo enviado por el camino r ∈ R. El conjunto F de flujos factibles viene dado por los pares
de vectores (x, y) con x = (xr )r∈R correspondiente a los flujos en rutas, y = (ye )e∈E vector
de flujo en cada arco, que satisfacen
X X
u= xr , ye = xr .
r∈R r3e
El costo de cada arco es ce (ye ) con ce : R+ → R+ continua no decreciente, de modo que los
costos de las rutas y el costo mínimo son:
X
cr (y) = ce (ye ), γ(y) = mı́n cr (y). (1)
r∈R
e∈r
Con esto, un flujo (x, y) ∈ F es llamado equilibrio de Wardrop si solo utiliza rutas mínimas,
es decir, si para cada r ∈ R se tiene xr > 0 ⇒ cr (y) = γ(y).
Para mayores detalles sobre equilibrios de Wardrop, revisar [5]. Si bien este modelo de
equilibrios es útil para el caso de tener un flujo continuo de partículas (jugadores) ingresando
a la red, considerando la congestión que generan, sufre de limitaciones que no permiten
trabajar a lo largo de la tesis con él. En efecto, suponga que comienza a ingresar flujo a
tasa constante a la red. En el instante en que la primera partícula infinitesimal sale de s,
se decide instantáneamente cuál o cuáles son los caminos óptimos por los que debe viajar el
torrente de flujo, correspondientes a aquellos que satisfacen el mínimo en 1. Una vez tomada
la decisión de qué caminos utilizar, el flujo sigue viajando por ellos hasta alcanzar el instante
de tiempo T en que dejan de ingresar jugadores a la red. Es por esta razón que no se puede
estudiar cómo la congestión se propaga a través de la red a lo largo del tiempo, la cual puede
provocar cambios en 1 y en los caminos que satisfacen el mínimo. La razón de estos posibles
cambios es la presencia de capacidades en el modelo, generándose colas en los arcos del grafo
que modifican los tiempos de viaje.
3
Como se vio en el ejemplo de la sección anterior, la saturación de las colas provoca un
cambio en los caminos óptimos de viaje, fenómeno que se refleja en la realidad donde el uso
creciente de sistemas de ruteo guiado que utilizan información en linea del estado de la red,
tales como Tom-Tom o Waze, plantean la necesidad de considerar los fenómenos temporales
y espaciales en la propagación de los flujos y la congestión. La existencia de atochamientos
en las calles requiere un estudio del flujo que considere la variación de los caminos óptimos
a lo largo del tiempo.
El estudio de flujos dinámicos fue inicialmente hecho en el área de optimización. En los

años 1958 [9] y 1962 [8], LR. Ford y DR. Fulkerson estudiaron el caso de tiempo discreto
logrando encontrar su conocido algoritmo de flujo máximo en una red. En 1959 D. Gale [11]
demostró la existencia de distribuciones de un flujo de modo que en tiempo mínimo alcance
cada vértice, incluyendo al nodo destino. A finales del siglo XX, EJ. Anderson, AB. Philpott
[1] en 1994 y L. Fleischer, É. Tardos en 1998 [7] lograron extender estos resultados al caso
continuo.
Por otro lado, el estudio de flujos dinámicos considerando que las partículas se comportan
egoístamente ha sido principalmente considerado en modelos de transporte. Uno de los pri-
meros en trabajar con modelos de colas fluidas vistos desde la teoría de juegos fue Vickrey
en 1969 [19], quien buscó formas de reducir la congestión en sistemas de transporte. En 1993,
Friesz, Bernstein, Smith, Tobin y Wie [10] propusieron una nueva clase de modelos, los cuales
corresponden a una generalización dinámica de los equilibrios de Wardrop estáticos, utilizan-
do formulaciones mediante desigualdades variacionales. En 1996, Ran y Boyce [16] publicaron
un libro resumiendo sus más recientes resultados en modelos de redes de transporte para el
caso dinámico, trabajando con herramientas similares a [10]. Si bien estos últimos trabajados
logran dar modelos con reglas bien generales de propagación del flujo o dinámicas de los
arcos, incluyendo como caso particular a los modelos de colas fluidas, pocos resultados de
existencia, unicidad y caracterización de soluciones se conocen.
En el año 2000, Zhu y Marcotte [20] lograron obtener resultados de existencia de equi-
librios dinámicos para redes de transporte modelando el problema como una desigualdad
variacional en dimensión infinita. Sin embargo, los supuestos necesarios para esto dejan fuera
los modelos de colas fluidas. En 2010, Meunier y Wagner [15] también obtuvieron resulta-
dos de existencia de equilibrios dinámicos utilizando propiedades de teoría de juegos con un
continuo de jugadores. Sin embargo, este último resultado, por ser demasiado abstracto, no
entrega herramientas para calcular el equilibrio.
En 2011, Koch y Skutella [13] lograron modelar una red de cola fluida y caracterizar
las derivadas de los equilibrios dinámicos relacionándolos con flujos finos normalizados, de-
mostrando así existencia y unicidad de equilibrios. En 2015, Cominetti, Correa y Larré [3]
continuaron el trabajo hecho por Koch y Skutella, resolviendo problemas en el modelo y
obteniendo resultados de existencia y unicidad correctos. Estos últimos dos trabajos dieron
paso además a algoritmos que calculan los equilibrios dinámicos de una red de transporte.
4
Objetivos y resultados.
Los trabajos [13] y [3] se limitan al caso de redes de colas fluida con solo un origen y
un destino. Debido a esto último, el objetivo principal de esta tesis consiste en generalizar
los resultados de existencia, unicidad y caracterización de estos trabajos al caso multi-od,
entregando un algoritmo que pueda encontrar los equilibrios dinámicos de la red.
Para el caso de poseer solo un origen y múltiples destinos, se entrega el modelo de colas
fluidas de la red, se definen y caracterizan sus equilibrios dinámicos, además de demostrar
existencia de estos y unicidad en sus etiquetas, como se verá en el capítulo 2. Además,
utilizando una estrategia de flujos finos normalizados, se caracterizan las derivadas de los
equilibrios dinámicos, lo que da origen a un algoritmo de extensión de flujos. Basado en el
trabajo de R. Koch en [12], un segundo algoritmo es propuesto, el cual determina los flujos
finos normalizados de una red en un instante de tiempo θ utilizando herramientas de cortes
con congestión máxima en grafos.
En el caso de múltiples orígenes y destinos, tras entregar el modelo de colas fluidas de

la red, se describen las principales complicaciones para extender el concepto de flujo fino
normalizado al poseer múltiples fuentes. Además, se muestran ejemplos sencillos en los cuales
la red multi-od pierde linealidad en ecuaciones fundamentales utilizadas para el cálculo de
equilibrios dinámicos, dificultando el uso de un algoritmo como en el caso de un solo origen.
5
Capítulo 1
Preliminares.
En este primer capítulo se presentan algunas definiciones y resultados útiles para el resto
de la tesis. Para mayores detalles en teoría de grafos se recomienda revisar [6], mientras que
para las herramientas de teoría de la medida se recomienda revisar [17].
Definición 1.1 Una función f : I → R se dice absolutamente continua si y solo si,

∀ε > 0, ∃δ > 0 tal que para toda familia finita ((ai , bi ))i=1,...,k de intervalos de I disjuntos a
pares, se satisface
Xk k
X
|bi − ai | < δ ⇒ |f (bi ) − f (ai )| < ε
i=1 i=1
Al conjunto de todas las funciones absolutamente continuas se denota AC(R). Se define

además ACloc (R) como el conjunto de todas las funciones absolutamente continuas sobre
cada intervalo acotado.
Definición 1.2 Se define el espacio de funciones integrables sobre R como

 
 Z 
L1 (R) := f : R → R : f es medible y |f | < ∞ .
 
R
Se define además el espacio L1loc (R), como aquel que posee a todas las funciones integrables
sobre cada conjunto compacto de R, es decir,
 
 Z 
1
Lloc (R) := f : R → R : f es medible y |f | < ∞, ∀K ⊆ R, K compacto .
 
K
Este último espacio recibe el nombre de espacio de funciones localmente integrables.
Proposición 1.1 Sea f : R+ → R+ una función localmente integrable. La función F : R+ →

R+ definida por
Zθ
θ → f (η)dη
0
6
es monótona creciente y localmente absolutamente continua, con F (0) = 0. Además, F es
diferenciable c.t.p. en el eje no negativo R+ y cumple con dF
dθ
(θ) = f (θ).
Definición 1.3 Un grafo dirigido consiste en un par ordenado G := (V, E) donde V es

un conjunto de vértices y E ⊆ V × V es un conjunto ordenado de aristas. Los elementos de
E son pares de vértices de la forma e = (v, w) o simplemente e = vw. Para cada elemento
e = vw ∈ E, v corresponderá a su cola y w a su cabeza.
Definición 1.4 Dado un vértice v ∈ V en un grafo dirigido, el conjunto de arcos salientes

desde v se define por
δ + (v) := {e = vw ∈ E}
Por otro lado, el conjunto de arcos entrantes a v se define mediante
δ − (v) := {e = uv ∈ E}
Dado S ⊆ V un subconjunto de vértices, los conjuntos de arcos salientes y entrantes a S se

definen respectivamente por
δ + (S) := {e = vw ∈ E : v ∈ S, w ∈
/ S}
δ − (S) := {e = vw ∈ E : v ∈
/ S, w ∈ S}
Definición 1.5 Sea φ : X → P(Y ) una multiaplicación con X, Y conjuntos no vacíos. Se

dice que φ es cerrada en x ∈ X si ∀(xn )n ⊆ X, ∀y ∈ Y y ∀yn ∈ φ(xn )
lı́m xn = x, lı́m yn = y ⇒ y ∈ φ(x)

n→∞ n→∞
Teorema 1.1 Teorema de punto fijo de Kakutani: Sea S ⊆ Rn un conjunto no vacío,

compacto y convexo. Sea φ : S → P(S) una multiaplicación cerrada, tal que ∀x ∈ S, φ(x) es
no vacío, convexo y cerrado. Entonces existe x ∈ S tal que x ∈ φ(x).
Lema 1.1 Sea g : R → R+ una función no negativa en L1loc (R) y {(ai , bi )}i∈I una familia no
S de intervalos. Entonces, g se desvanece c.t.p en cada (ai , bi ) si y solo si se desvanece
numerable
c.t.p en i∈I (ai , bi ). Más aún, el lema es válido incluso para intervalos semi-abiertos de la
forma [ai , bi ), (ai , bi ].
Demostración. RSupongamos sin pérdida de generalidad que todos los intervalos son no vacíos.
ComoSµ(A) = A g(ξ)dξ define una medida regular en los conjuntos borelianos A ⊆ R, para
Θ = i∈I (ai , bi ) tenemos
µ(Θ) = sup{µ(K) : K compacto, K ⊆ Θ}.
Ahora, cada compacto K ⊆ Θ posee un recubrimiento finito K ⊆ nk=1 (ai , bi ) por lo tanto
S
bik
n
X n Z
X
µ(K) ≤ µ((aik , bik )) = g(ξ)dξ = 0
k=1 k=1 a
ik
7
Por lo tanto µ(Θ) = 0, lo que implica que g(ξ) = 0 para casi todo ξ ∈ Θ, lo que prueba una
de las implicancias. La otra es directa a partir de que todo intervalo (ai , bi ) es subconjunto
de Θ.
La segunda parte del lema se obtiene del hecho que las tres uniones
S difieren aSlo más en
un conjunto numerable de puntos. En efecto, consideremos N = i∈I [ai , bi ) \ i∈I (ai , bi ).
Cada punto z ∈ N debe ser un punto final z = ai y su correspondiente intervalo (ai , bi )
es disjunto de N . Dado aj ∈ N otro punto, los intervalos correspondientes deben ser de
intersección vacía, y por lo tanto hay a lo más numerable de ellos, pues toda colección
infinita
S de intervalosS abiertos en R es numerable. Un argumento similar muestra lo mismo
para i∈I (ai , bi ] \ i∈I (ai , bi ).
Proposición 1.2 Sean h ∈ ACloc (R) y g ∈ L1loc (R). Si h es monótona, se cumple la fórmula
de cambio de variable:
Zh(b) Zb
g(ξ)dξ = g(h(y))h0 (y)dy.
h(a) a
8
Capítulo 2
Modelo de cola fluida, caso de un origen

y múltiples destinos.
En este capítulo se estudia el modelo matemático para una red de flujo dinámico cuyos
arcos presentan latencias y tasas de salida o capacidad. Considere un grafo dirigido G = (V, E)
con V el conjunto de vértices y E el conjunto de arcos. Sea s un vértice origen desde el cual
se inyecta flujo hacia el grafo y {ti }i∈I ⊆ V un conjunto de nodos destino. Se considerará el
siguiente conjunto de funciones
F0 (R) := g ∈ L1loc (R) : g ≥ 0, g(θ) = 0 para casi todo θ ∈ (−∞, 0) .

(2.1)
es decir, el conjunto de funciones localmente integrables, no negativas y que desvanecen en
el eje negativo.
2.1. Modelo de red de flujo dinámico.
Dado i ∈ I, se considera ui ∈ F0 (R) la función de flujo entrante al grafo desde el nodo s

que debe llegar a ti . De este modo, el flujo total entrante al grafo está dado por
X
u(θ) := ui (θ). (2.2)
i∈I
2.1.1. Flujo dinámico.
A lo largo del modelo se trabaja con tasas de flujo en cada arista que se dirigen a algún
destino, dando paso a la siguiente definición.
Definición 2.1 Un flujo dinámico o simplemente flujo, será un par (f + , f − ) de arreglos

+ − + − +
de funciones (fi,e ), (fi,e ) con fi,e , fi,e ∈ F0 (R) para cada e ∈ E, i ∈ I, donde fi,e corresponde
al flujo entrante a la cola del arco e proveniente desde el nodo origen s que se dirige al destino
−
ti y fi,e al flujo saliente de la cabeza del arco e que debe llegar al nodo ti .
9
Se define el flujo total entrante a cada arco y el flujo total saliente de cada arco
respectivamente por
X X
+ −
fe+ := fi,e , fe− := fi,e . (2.3)
i∈I i∈I
Como |I| < ∞ claramente fe+ y fe− son funciones de F0 (R) para cada e ∈ E.
Definición 2.2 Las funciones flujo entrante acumulativo y flujo saliente acumulativo
para cada i ∈ I, e ∈ E se definen respectivamente por
Zθ Zθ
+ + − −
Fi,e (θ) := fi,e (ξ)dξ, Fi,e (θ) := fi,e (ξ)dξ. (2.4)
0 0
Estas corresponden a todo el flujo que ha entrado a la cola del arco e y a todo el flujo que
ha salido de la cabeza de e, que se dirigen a ti , en el intervalo de tiempo [0, θ]. Notar que
+ − + −
como fi,e y fi,e son localmente integrable se cumple que Fi,e y Fi,e son funciones de ACloc (R).
Dado e ∈ E, se definen los flujos acumulativos totales como

X X
+ −
Fe+ (θ) = Fi,e (θ), Fe− (θ) = Fi,e (θ). (2.5)
i∈I i∈I
las cuales, nuevamente, resultan ser elementos de ACloc (R).
2.1.2. Cola fluida.
Un arco e ∈ E, como se explicó en la introducción, es modelado por una cola fluida con
tasa de salida νe seguida por un link con tiempo de viaje o latencia constante τe .
νe
τe
Cola Link
Figura 2.1: Arco
Definición 2.3 Se define el largo de cola como la función ze : R → R dada por

ze (θ) := Fe+ (θ) − Fe− (θ + τe ). (2.6)
es decir, todo el flujo que ha entrado al arco pero aún no ha salido ni se encuentra actualmente
viajando por el link.
Definición 2.4 Se dice que las colas operan bajo capacidad si para casi todo θ ∈ R se
cumple

 νe si ze (θ) > 0,
−
fe (θ + τe ) = (2.7)
mı́n{fe+ (θ), νe } si no.

10
o equivalentemente en términos de la derivada de la función largo de cola,

 fe+ (θ) − νe si ze (θ) > 0,
ze0 (θ) = (2.8)
 +
[fe (θ) − νe ]+ si no.
La única solución de esta ecuación diferencial viene dada por (ver [2], páginas 3-11.)
Zθ
ze (θ) = máx [fe+ (ξ) − νe ] dξ. (2.9)
η∈[0,θ]
η
Esta fórmula muestra que si las colas operan bajo capacidad, el flujo entrante fe+ al arco
e, determina completamente el largo de cola ze y más aún, en virtud de 2.7, también el flujo
saliente fe− .
2.1.3. Tiempo de viaje por un arco.
Definición 2.5 Cuando una partícula de flujo llega a un arco e en un instante θ, el Tiempo
de espera que experimenta en la cola hasta poder comenzar a recorrer el arco está definido
por
Zθ+q
 
 
−
qe (θ) := mı́n q ≥ 0 : fe (ξ + τe )dξ = ze (θ) .
 
θ
es decir, el primer instante q en que el flujo saliente de la cabeza de e a partir de θ + τe iguala

el tamaño de la cola. Notar que si una partícula de flujo sale del arco en ξ + τe , entonces en
ξ se encuentra dejando la cola para comenzar a viajar por el link.
Observación 2.1 Dados fe+ , fe− los flujos de entrada y salida de un arco e fijo, la ope-
ración bajo capacidad resulta ser equivalente a que se cumplen las siguientes condiciones
simultáneamente:
(a) Restricción de capacidad: fe− (θ) ≤ νe para casi todo θ ∈ R.

(b) Restricción de no déficit: ze (θ) ≥ 0 ∀θ ∈ R.
ze (θ)
(c) Tiempo de espera: qe (θ) = νe
, ∀θ ∈ R.
Definición 2.6 Se define la función tiempo de viaje por un arco como Te : R → R dada
por la suma del tiempo de entrada θ, el tiempo de espera en la cola y la latencia, es decir,
ze (θ)
Te (θ) := θ + + τe . (2.10)
νe
La función ze (·) pertenece a ACloc (R) pues corresponde a la resta de dos funciones local-
mente absolutamente continuas, de lo cual se sigue que Te ∈ ACloc (R), ∀e ∈ E. En particular,
11
Te es derivable c.t.p y su derivada está dada para casi todo θ ∈ R por
1 +


 f (θ)
νe e
si ze (θ) > 0,
0
Te (θ) = n o (2.11)
 máx 1, 1 f + (θ)

si no.
νe e
De esto resulta que Te0 (θ) ≥ 0 por lo cual Te es no decreciente y por lo tanto las partículas
atravesando e respetan FIFO. Más aún, al igual que en [3] se tiene que todo el flujo entrante
a e en tiempo θ sale en Te (θ), es decir
Fe− (Te (θ)) = Fe+ (θ). (2.12)
Observación 2.2 Como Te (θ) → −∞ cuando θ → −∞ pues en ese caso las filas están
todas vacías, es decir, ze (θ) = 0, sumado al hecho de que Te (θ) ≥ θ por lo cual Te (θ) → ∞
cuando θ → ∞, se concluye que Te es sobreyectiva para todo e ∈ E en virtud de que son
funciones continuas.
El modelo recién descrito se utilizará para trabajar una red de flujo que posee un nodo
fuente enviando flujo a varios destinos, en la cual se busca que cada partícula infinitesimal
de este viaje siempre por el camino más corto posible considerando latencias y colas.
2.1.4. Existencia y unicidad de descomposición de un flujo de salida.
En esta sección se verá que un flujo saliente de un arco se puede descomponer de manera
única en flujos con destinos fijos respetando la propiedad FIFO para cada uno.
Lema 2.1 Sea un flujo de entrada descompuesto en nodos destinos

X
+
fe+ (θ) = fi,e (θ)
i∈I
− −
definidos para todo θ ∈ R. Entonces existen flujos únicos fi,e ∈ L∞ (R) con 0 ≤ fi,e (ξ) ≤
νe , ∀ξ ∈ R que satisfacen
TZe (θ) Zθ
− +
fi,e (ξ)dξ = fi,e (ξ)dξ (2.13)
0 0
para todo θ ∈ R, i ∈ I, e ∈ E.
Demostración. Sea A ⊆ R el conjunto de todos los tiempos θ ∈ R tal que Te0 (θ) existe y es
positiva. Recordar que como Te es absolutamente continua su derivada existe c.t.p. Se definen
las funciones fi,e , ∀e ∈ E, i ∈ I mediante
 1 +
 Te0 (θ) fi,e (θ) si θ ∈ A
−
fi,e (Te (θ)) =
0 si no.

12
Como la función Te es sobreyectiva, la definición se tiene sobre todo R. Lo primero es ver
que estos flujos de salida descompuestos en destinos están efectivamente bien definidos. Sean
θ, θ0 tal que θ < θ0 , Te0 existe en θ y θ0 , y Te (θ) = Te (θ0 ) = α. Como la función Te es creciente,
se tiene que Te (ξ) = α, ∀ξ ∈ [θ, θ0 ], luego Te0 = 0 en [θ, θ0 ] en aquellos puntos donde existe la
derivada, con lo cual θ, θ0 ∈ Ac y luego
− −
fi,e (Te (θ)) = fi,e (Te (θ0 )) = 0
por lo tanto la definición es consistente.
Sea θ ∈ A, entonces
X
− 1 X + fe+ (θ)
fi,e (Te (θ)) = f (θ) =
i∈I
Te0 (θ) i∈I i,e Te0 (θ)
y utilizando 2.11 se obtiene que

X  νe si ze (θ) > 0,
−
fi,e (Te (θ)) =
mı́n{fe+ (θ), νe } si no.

i∈I
− −
P
De este modo todo fi,e (Te (θ)) ≤ νe , ∀θ ∈ R. Como cada fi,e es no negativo se concluye
i∈I
− −
que fi,e ≤ νe , ∀i ∈ I, luego fi,e ∈ L∞ (R), ∀i ∈ I.
−
Dado θ ∈ R con Te (θ) < 0, si θ ∈ Ac se tendrá que fi,e (Te (θ)) = 0. Si θ ∈ A entonces
+
como Te (θ) < 0 ≤ Te (0) se concluye que θ < 0 y luego como fi,e desvanece en el eje negativo
−
se obtiene que fi,e también pues Te es sobreyectiva. Con esto,
TZe (θ) TZe (θ) Zθ

− − −
fi,e (ξ)dξ = fi,e (ξ)dξ = fi,e (Te (ξ))Te0 (ξ)dξ
0 −∞ −∞
donde la última igualdad sale de utilizar un cambio de variable. Luego, si ξ ∈ Ac , Te0 (ξ) = 0
ya que no se consideran aquellos tiempos ξ donde la derivada no está definida pues es un
conjunto de medida nula, entonces, de 2.11 se concluye que fe+ (ξ) = 0 c.t.p. y en particular
+
fi,e (ξ) = 0 c.t.p. Por lo tanto para casi todo ξ ∈ R se cumple la igualdad
− +
fi,e (Te (ξ))Te0 (ξ) = fi,e (ξ)
ya que si ξ ∈ A se obtiene directo de la definición. Se concluye entonces que
TZe (θ) Zθ Zθ
− + +
fi,e (ξ)dξ = fi,e (ξ)dξ = fi,e (ξ)dξ
0 −∞ 0
+
ya que fi,e desvanece en el eje negativo, demostrando así la existencia.
La unicidad sale del hecho de que como Te (·) es sobreyectiva, si dos funciones g, h cumplen
con 2.13 entonces
Zα Zα
g(ξ)dξ = h(ξ)dξ, ∀α ≥ 0
0 0
con lo cual g = h c.t.p, es decir, son iguales vistas como funciones de L1loc (R).
13
2.1.5. Caminos dinámicos mínimos.
Si una partícula proveniente de s ingresa a un camino P = (e1 , e2 , ..., ek ) en tiempo θ,

llega al final de este en tiempo
lP (θ) = Tek ◦ · · · ◦ Te1 (θ).
Sea Pw el conjunto de todos los s-w caminos en G. El tiempo más corto en que una partícula
saliendo de s en θ logra llegar a w viene dado por
lw (θ) = mı́n lP (θ). (2.14)

P ∈Pw
Estas funciones corresponden a los caminos más cortos sobre arcos que poseen latencias y
retraso de colas en los momentos correspondientes considerando el tiempo que toma alcanzar
cada arco. Reciben el nombre de caminos dinámicos mínimos.
Gracias a que las funciones Te son absolutamente continuas y no decrecientes lo mismo se

tiene para lP y luego para lw , ∀w ∈ V . Además, por argumentos similares a los dados para
Te también se obtiene que lw (θ) → ±∞ cuando θ → ±∞ además de la sobreyectividad de
estas funciones.
Asumiendo que la suma de las latencias de cualquier ciclo es siempre positiva, además
de la monotonía de Te y la restricción de no déficit se obtiene que los caminos dinámicos
mínimos no poseen ciclos y por lo tanto pueden ser calculados resolviendo
(
θ si w = s,
lw (θ) = mı́n− Te (lv (θ)) si w =
6 s. (2.15)
e=vw∈δ (w)
para cualquier w ∈ V .
Para θ ∈ R se define el conjunto
Eθ0 := {e = vw ∈ E : Te (lv (θ)) ≤ lw (θ)} (2.16)
es decir toda arista e que pertenece a un camino dinámico mínimo en el momento θ. Un arco
e ∈ Eθ0 se dirá activo. Notar que es equivalente Te (lv (θ)) ≤ lw (θ) a Te (lv (θ)) = lw (θ).
Definición 2.7 Se define el grafo de caminos θ mínimos como el grafo acíclico Gθ =

(V, Eθ0 ).
Un arco e = vw ∈ E inactivo cumple Te (lv (θ)) > lw (θ) y gracias a la continuidad de estas
funciones la desigualdad se seguirá cumpliendo durante un lapso de tiempo más. Se definen
también los conjuntos
Θe := {θ ∈ R : e ∈ Eθ0 }.
es decir todos los tiempos θ donde e es activo.
14
2.1.6. Equilibrio Dinámico.
Se comienza esta sección dando la definición de flujo factible.
+ −
Definición 2.8 Un flujo dinámico f = (fi,e , fi,e )(i,e)∈I×E se dirá factible si para todo i ∈ I
cumple la siguiente restricción de conservación de flujo para casi todo θ ∈ R

X
+
X
− ui (θ) si v = s,
fi,e (θ) − fi,e (θ) = (2.17)
+ −
0 si v 6= {s, ti }.
e∈δ (v) e∈δ (v)
Definición 2.9 Un flujo dinámico factible f es llamado Equilibrio Dinámico si para cada
+
e = vw ∈ E, i ∈ I se cumple fi,e (ξ) = 0 para casi todo ξ ∈ lv (Θce ).
El siguiente lema da una caracterización de equilibrio dinámico.
Lema 2.2 Un flujo dinámico factible f es un equilibrio dinámico si y solo si ∀e = vw ∈ E,

+
∀i ∈ I, y casi todo ξ ∈ R se cumple fi,e (ξ) > 0 ⇒ ξ ∈ lv (Θe ).
+
Demostración. Sean i ∈ I, e ∈ E fijos. Notar que fi,e (ξ) > 0 ⇒ ξ ∈ lv (Θe ) c.t.p. es
+ c
equivalente a fi,e (ξ) = 0 para ξ ∈ lv (Θe ) c.t.p. Luego, bastará con que ver los conjuntos
lv (Θce ) y lv (Θe )c difieren en un conjunto de medida nula. Sea ξ ∈ lv (Θe )c , como la función
lv (·) es sobreyectiva, existe θ ∈ R tal que ξ = lv (θ). Se cumple que θ ∈ Θce ya que sino,
lv (θ) ∈ lv (Θe ) es decir ξ ∈ lv (Θe ) lo que es una contradicción. Por lo tanto lv (θ) ∈ lv (Θce ) con
lo cual se concluye que lv (Θe )c ⊆ lv (Θce ).
Sea ahora ξ ∈ lv (Θce ) \ lv (Θe )c = lv (Θce ) ∩ lv (Θe ), entonces existen θ, θ̄ con θ ∈ Θce , θ̄ ∈ Θe
tal que ξ = lv (θ) = lv (θ̄). Como lv (·) es no decreciente, en realidad lv (η) = ξ para todo η
entre θ y θ̄, los cuales, como viven en conjuntos disjuntos son distintos, así existe η ∈ Q tal
que ξ = lv (η), es decir ξ ∈ lv (Q) que es un conjunto numerable y por lo tanto de medida nula.
Se concluye entonces que lv (Θce ) \ lv (Θe )c está contenido en un conjunto de medida nula.
Intuitivamente es claro que en un equilibrio todo arco con largo de cola distinto de 0 debe
ser activo. Sean θ ∈ R y el conjunto
Eθ∗ := {e = vw ∈ E : ze (lv (θ)) > 0} (2.18)
La siguiente proposición demuestra que el resultado intuitivo es cierto.
Proposición 2.1 Si f es un equilibrio dinámico entonces Eθ∗ ⊆ Eθ0 y además
Eθ0 = {e = vw ∈ E : lw (θ) ≥ lv (θ) + τe } (2.19)
Eθ∗ = {e = vw ∈ E : lw (θ) > lv (θ) + τe } (2.20)
Demostración. Sean e = vw ∈ Eθ∗ y θ0 ≤ θ el último instante en que e fue activo. Como f es

+
equilibrio dinámico, fi,e (ξ) = 0 para todo i ∈ I, c.t.p. en el intervalo (lv (θ0 ), lv (θ)], luego, como
15
ze (lv (θ)) > 0, el arco debe haber tenido cola a lo largo de todo el intervalo. En particular, de
2.11 se obtiene que Te0 (ξ) = 0 c.t.p. en (lv (θ0 ), lv (θ)], o equivalentemente Te es constante en el
intervalo. Entonces
Te (lv (θ)) = Te (lv (θ0 )) = lw (θ0 ) ≤ lw (θ)
donde la primera igualdad es debido a que Te es constante, la segunda a que e es activo en θ0
y la desigualdad gracias a que lw (·) es una función creciente. En particular Te (lv (θ)) ≤ lw (θ),
concluyendo que e ∈ Eθ0 y probando así la inclusión.
Ahora para probar 2.19, notando que si e ∈ Eθ0 se tiene que

lw (θ) = Te (lv (θ)) ≥ lv (θ) + τe
donde la desigualdad sale de 2.10 y de la restricción de no déficit. Por otro lado, si lw (θ) ≥
lv (θ) + τe , suponga que ze (lv (θ)) = 0. Entonces lw (θ) ≥ Te (lv (θ)) y con esto e ∈ Eθ0 . Si
ze (lv (θ)) > 0, entonces e ∈ Eθ∗ y gracias a que Eθ∗ ⊆ Eθ0 , se obtiene que e ∈ Eθ0 .
De forma similar se puede demostrar 2.20. Sea e ∈ Eθ∗ , entonces ze (lv (θ)) > 0 y como
e ∈ Eθ0 , lw (θ) = Te (lv (θ)) > lv (θ) + τe . Recíprocamente, si lv (θ) + τe < lw (θ), como lw (θ) ≤
Te (lv (θ)) por definición, 2.10 entrega que necesariamente ze (lv (θ)) > 0, luego e ∈ Eθ∗ .
El siguiente teorema da una caracterización para el equilibrio dinámico que ayudará a

continuar el estudio de este.
Teorema 2.1 Sea f un flujo dinámico factible. Son equivalentes:
(a) f es un equilibrio dinámico.

+ −
(b) Para todo e = vw ∈ E, i ∈ I y para todo θ ∈ R se tiene Fi,e (lv (θ)) = Fi,e (lw (θ)).
+
/ Eθ ⇒ fi,e (lv (θ))lv0 (θ) = 0.
(c) Para todo e = vw ∈ E, i ∈ I y para casi todo θ ∈ R, e ∈ 0
Demostración. Sean i ∈ I, e = vw ∈ E y θ ∈ R. Sea θ0 ≤ θ el máximo instante tal que

Te (lv (θ0 )) = lw (θ). Notar que θ0 existe puesto que, si e es activo para θ entonces θ0 = θ, por
otro lado si e es inactivo, lw (θ) < Te (lv (θ)), y como Te (lv (θ)) → −∞ cuando θ → −∞ y
Te (lv (·)) es continua, existe un θ0 que cumple con lo pedido.
Sea Iθ = (θ0 , θ]. Cabe decir que θ ∈ Θe si y solo si θ0 = θ con lo cual Iθ = ∅. Lo primero
es demostrar que Θce = ∪θ Iθ . En efecto. si θ ∈ Θce , se tiene que θ0 < θ y θ ∈ Iθ lo que entrega
una primera inclusión. Recíprocamente, sea θ00 ∈ Iθ . Por la definición de θ0 y la monotonía
de Te (·), lw (·) se cumple que
lw (θ00 ) ≤ lw (θ) = Te (lv (θ0 )) < Te (lv (θ00 ))
donde la última desigualdad debe ser estricta pues si Te (lv (θ0 )) = Te (lv (θ00 )) se contradice que
θ0 < θ00 y la definición de θ0 . Luego, θ00 ∈ Θce concluyendo la otra inclusión.
Recordando el Lema 2.1, para todo θ ∈ R se cumple que

TZe (θ) Zθ
− +
fi,e (ξ)dξ = fi,e (ξ)dξ
0 0
16
por lo tanto
lZ
v (θ)
+ − + − + + +
Fi,e (lv (θ))−Fi,e (lw (θ)) = Fi,e (lv (θ))−Fi,e (Te (lv (θ0 ))) = Fi,e (lv (θ))−Fi,e (lv (θ0 )) = fi,e (ξ)dξ
lv (θ0 )
+
donde esta integral es siempre mayor o igual a 0 pues fi,e ∈ F0 (R). Luego,
+ −
Fi,e (lv (θ)) − Fi,e (lw (θ)) = 0
+
si y solo si fi,e se desvanece c.t.p en el intervalo (lv (θ0 ), lv (θ)] = lv (Iθ ). En virtud del Lema
+
1.1, esto se tiene si y solo si fi,e (ξ) = 0 para ξ ∈ ∪θ lv (Iθ ) = lv (Θce ) c.t.p, lo que prueba que
(a) y (b) son equivalentes.
Ahora, partiendo desde

lZ
v (θ)
+ − +
Fi,e (lv (θ)) − Fi,e (lw (θ)) = fi,e (ξ)dξ
lv (θ0 )
y realizando un cambio de variable el cual vale gracias a que las funciones lv son absolutamente
continuas y monótonas, se obtiene
Zθ
+ − +
Fi,e (lv (θ)) − Fi,e (lw (θ)) = fi,e (lv (ξ))lv0 (ξ)dξ
θ0
+
Al igual que antes esta integral es igual a 0 si y solo si fi,e (lv (ξ))lv0 (ξ) = 0 para todo ξ ∈ Iθ
c.t.p. Nuevamente, en virtud del Lema 1.1, esto ocurre si y solo si esta función se desvanece
en ∪θ Iθ = Θce c.t.p, probando la equivalencia entre (b) y (c).
Este teorema entrega una manera útil de caracterizar un equilibrio dinámico dentro de
una red de flujo, tal como se verá en la siguiente sección.
2.1.7. Flujo acumulativo y Flujo fino normalizado
La parte (b) del Teorema 2.1 entrega una forma de sincronizar el flujo dinámico descom-
puesto en destinos fijos usando el tiempo de salida desde s como un reloj común. Esto motiva
la siguiente definición:
Definición 2.10 Se define el flujo acumulativo inducido por un equilibrio dinámico f

+ −
como x(θ) = (xi,e (θ))i,e para cada i ∈ I, e = vw ∈ E con xi,e (θ) = Fi,e (lv (θ)) = Fi,e (lw (θ)),
∀θ ∈ R.
Estos flujos acumulativos satisfacen la restricción de conservación de flujo. En efecto,

integrando la ecuación 2.17 en el intervalo [0, lv (θ)], se obtiene que para cada θ ∈ R, i ∈ I:

X X Ui (θ) si v = s
xi,e (θ) − xi,e (θ) = (2.21)
+ −
0 si v 6= {s, ti }
17
Rθ
donde Ui (θ) = ui (ξ)dξ y utilizando el hecho que ls (θ) = θ. Es decir, cada xi (θ) corresponde
0
a un s-ti flujo estático de valor Ui (θ). Como las funciones xi,e también pertenecen a ACloc (R),
derivando a ambos lados de la igualdad anterior se obtiene que para casi todo θ ∈ R

X
0
X
0 ui (θ) si v = s
xi,e (θ) − xi,e (θ) = (2.22)
+ −
0 si v 6= {s, ti }
Es decir, los flujos x0i (θ) := (x0i,e (θ))e∈E para aquellos θ donde la derivada existe, son s-ti
flujos estáticos de valor ui (θ). Además, en virtud de la parte (c) del Teorema 2.1, para todo
e∈ / Eθ0 se tiene x0i,e (θ) = 0.
La estrategia para abordar el problema de obtener equilibrios dinámicos de una red será
relacionar estos con las derivadas de los flujos acumulativos. Para esto, primero se caracterizan
las derivadas de las funciones xi,e y lw ya que al ser estas funciones absolutamente continuas,
pueden ser recuperadas a partir de integrar x0i,e y lw0 respectivamente.
Como ya se dijo, las funciones x0i (θ) son s-ti flujos estáticos de valor ui (θ) para aquellos θ
donde la derivada existe, con x0i,e (θ) = 0 si e ∈
/ Eθ0 . Además se tiene claramente que lw0 (θ) = 1
para w = s. Para w 6= s, comenzando de 2.15 y considerando la regla de derivación de un
mínimo de funciones:
lw0 (θ) = mı́n Te0 (lv (θ))lv0 (θ) (2.23)
e=vw∈Eθ0
donde el mínimo se realiza en realidad solo sobre los arcos activos pues sobre aquellos arcos
fuera de Eθ0 nunca se alcanzará el óptimo. Utilizando 2.11 se obtiene:
1 +
f (lv (θ))lv0 (θ) si e ∈ Eθ∗


 νe e
lw0 (θ) = mı́n 0 n o
e=vw∈Eθ 
 máx l0 (θ), 1 f + (lv (θ))l0 (θ) si e ∈
v νe e v / Eθ∗ .
P
Si se considera xe (θ) := xi,e (θ), sigue que
i∈I
1 0
si e ∈ Eθ∗


 x (θ)
νe e
lw0 (θ) = mı́n
e=vw∈Eθ0
n o
 máx l0 (θ), 1 x0 (θ) si e ∈

/ Eθ∗ .
v νe e
Por comodidad se definirá ρe (lv0 , x0e ) como

1 0
si e ∈ Eθ∗


 x
νe e
ρe (lv0 , x0e ) := n o (2.24)
 máx l0 , 1 x0

/ Eθ∗ .
si e ∈
v νe e
De este modo,
lw0 (θ) = mı́n ρe (lv0 (θ), x0e (θ)). (2.25)
e=vw∈Eθ0
18
Las observaciones previas conducen a introducir las siguientes definiciones.
Dado un par de conjuntos de aristas E ∗ , E 0 contenidos en E, se dirá que estos cumplen la

propiedad (H) si:
(H) E ∗ ⊆ E 0 ⊆ E, E 0 no posee ciclos y ∀v ∈ V, hay un s − v camino en E 0 .
Se denota por Ki (E 0 , ui ) al conjunto no vacío, compacto y convexo de todos los s-ti flujos
/ E 0.
estáticos x0i = (x0i,e )e∈E ≥ 0 de valor ui con x0i,e = 0 si e ∈
Dado un conjunto I finito, un elemento (x0i )i∈I ∈ K := Ki (E 0 , ui ), se define x0 = (x0e )e∈E

Q
i∈I
mediante
X
x0e = x0i,e , ∀e ∈ E (2.26)
i∈I
El vector x0 cumple con


X X  u si v = s
x0e − x0e = 0 si v 6= {s, ti }, ∀i ∈ I (2.27)
−u si v = ti , ∀i ∈ I

e∈δ + (v) e∈δ − (v)
Para cada x0 ∈ K, se asocia una etiqueta única l0 = (lv0 )v∈V dada por ls0 = 1,
lw0 = mı́n 0 ρe (lv0 , x0e )

e=vw∈E
Definición 2.11 Dado un arreglo de flujos estáticos (x0i )i∈I ∈ K y el flujo x0 = (x0e )e∈E
que estos definen, se dirá que el vector (x0i )i∈I es un flujo fino normalizado o N T F (de
sus siglas en inglés) de valor (ui )i∈I con reseteo en E ∗ ⊆ E 0 si y solo si x0e = 0 para cada
e = vw ∈ E 0 tal que lw0 < ρe (lv0 , x0e ).
Teorema 2.2 Sea f un equilibrio dinámico y θ ∈ R tal que las derivadas por la derecha
ui = (dUi /dθ+ )(θ), lv0 = (dlv /dθ+ )(θ) y x0i,e = (dxi,e /dθ+ )(θ) existen, para todo i ∈ I, e ∈ E,
v ∈ V . Entonces (x0i )i∈I es un NTF de valor (ui )i∈I con reseteo en Eθ∗ ⊆ Eθ0 y las etiquetas
correspondientes son l0 .
Demostración. Como xi es un s-ti flujo estático de valor Ui , derivando se obtiene que x0i es
un s-ti flujo estático de valor ui . Además, si e ∈ / Eθ0 entonces e permanece inactivo en un
+
intervalo de tiempo [θ, θ + ε) para cierto ε > 0, como x0i,e (ξ) = fi,e (lv (ξ))lv0 (ξ) y f es equilibrio
dinámico, x0i,e = 0 en [θ, θ + ε) c.t.p, lo que muestra que x0i ∈ K(E 0 , ui ).
El siguiente paso es ver que l0 corresponde a las etiquetas de x0 =

P 0
xi , es decir, que ls0 = 1
i∈I
y que lw0 = mı́n 0 ρe (lv0 , x0e ). Lo primero claramente se cumple pues ls (θ) = θ. Para ver la
e=vw∈E
otra parte, considere dos nuevos subconjuntos de Eθ0 :
E+∗ := {e = vw ∈ Eθ0 : ze (lv (θ)) > 0 o ∃ε > 0, ze (ξ) > 0, ∀ξ ∈ (lv (θ), lv (θ) + ε)},
19
E+0 := {e = vw ∈ Eθ0 : ze (lv (θ)) = 0 y ∃(θn )n con θn ↓ θ tal que e ∈ Eθ0 n , ∀n}.
Se demostrarán las siguientes cuatro propiedades:
(a) lw0 ≥ lv0 para e = vw ∈ E+0 . En efecto, para cada e ∈ E+0 se tiene que lw (θn ) =
Te (lv (θn )) ≥ lv (θn ) + τe y lw (θ) = lv (θ) + τe , luego lw (θn ) − lw (θ) ≥ lv (θn ) − lv (θ).
Dividiendo por θn − θ y haciendo n → ∞, se concluye que lw0 ≥ lv0 .
(b) lw0 ≤ lv0 para e = vw ∈ Eθ0 \ E+∗ . En efecto, para e ∈ Eθ0 \ E+∗ , considere θn ↓ θ con
ze (lv (θn )) = 0 y entonces lw (θn ) ≤ Te (lv (θn )) = lv (θn ) + τe , y como lw (θ) = lv (θ) + τe
pues e es activo en θ y por continuidad de ze (lv (·)) se tiene que ze (lv (θ)) = 0, de la
misma forma que antes se concluye que lw0 ≤ lv0 .
(c) lw0 ≥ x0e /νe para e = vw ∈ Eθ0 . En efecto, para e = vw ∈ Eθ0 la restricción de operar
bajo capacidad entrega que para todo θ0 ≥ θ:
lZ 0
w (θ )
xe (θ0 ) − xe (θ) = fe− (ξ)dξ ≤ νe (lw (θ0 ) − lw (θ))

lw (θ)
dividiendo por θ0 − θ y haciendo θ0 → θ se concluye que lw0 ≥ x0e /νe .

(d) lw0 = x0e /νe para e = vw ∈ E+∗ . En efecto, para e = vw ∈ E+∗ , existe θ0 > θ lo
suficientemente cerca para que ze (lv (θ0 )) > 0. Por lo tanto, la fila en e permanecerá
no vacía durante [lv (θ0 ), lv (θ0 ) + ze (lv (θ0 ))/νe ) pues no hay que olvidar que ze (lv (θ0 ))/νe
corresponde al tiempo que debe esperar una partícula que llega a la fila en el instante
lv (θ0 ). Entonces, como la fila es no vacía se tiene que fe− (ξ) = νe para ξ ∈ [lv (θ0 ) +
τe , lv (θ0 ) + ze (lv (θ0 ))/νe + τe ) = [lv (θ0 ) + τe , lw (θ0 )) c.t.p, y por lo tanto fe− (ξ) = νe c.t.p.
en un pequeño intervalo hacia la derecha de lw (θ), pues lw (θ0 ) ≥ lw (θ). De este modo
la igualdad en
lZ 0
w (θ )
xe (θ0 ) − xe (θ) = fe− (ξ)dξ ≤ νe (lw (θ0 ) − lw (θ))

lw (θ)
se mantiene para θ0 suficientemente cercano a θ, con lo cual lw0 = x0e /νe para e ∈ E+∗ .
Utilizando 2.24, al combinar (b) y (d) se obtiene que lw0 ≤ mı́n ρe (lv0 , x0e ), con igualdad
e=vw∈Eθ0
si existe un arco e ∈ Eθ∗ incidente a w. En caso que no exista una arista con cola positiva en
tiempo θ incidente a w, se considera una sucesión de tiempos θn ↓ θ estrictamente monótona
y una sucesión de arcos activos en ∈ Eθ0 n , y de esta última se toma una subsucesión (enk )k
con enk = vw constante para todo k ∈ N. Luego, e ∈ E+0 por la definición del conjunto,
y entonces combinando (a) y (c) se obtiene que lw0 ≥ ρe (lv0 , x0e ) y como ya se tiene la otra
desigualdad, se cumple que lw0 = mı́n 0 ρe (lv0 , x0e ) para w 6= s.
e=vw∈Eθ
Finalmente, queda ver que (x0i )i∈I corresponde a un NTF. Suponga que x0i,e > 0 para algún
e = vw ∈ Eθ0 con lw0 < ρe (lv0 , x0e ). Si e ∈ Eθ∗ , entonces e ∈ E+∗ , luego se debería cumplir (d),
pero de la definición de ρe (lv0 , x0e ) y lo asumido, se llega a una contradicción. Por lo tanto
e∈/ Eθ∗ . Por otro lado, como x0i,e > 0 entonces x0e > 0 y por lo tanto xe (θ0 ) > xe (θ) para todo
20
θ0 > θ, es decir, Fe+ (lv (θ0 )) > Fe+ (lv (θ)), por lo que e debe ser activo para una sucesión θn ↓ θ
y entonces e ∈ E+0 . Pero ya se vio que esto, junto con (a) y (c) indica que lw0 ≥ ρe (lv0 , x0e ),
contradiciendo lo supuesto en un comienzo.
El teorema anterior entrega la existencia de NTF’s a partir de un equilibrio dinámico.

Ahora se procede a ver la existencia de NTF’s en general, sin que esto provenga necesaria-
mente de un equilibrio dinámico, y en la siguiente sección la forma de obtener equilibrios
dinámicos a partir de estos.
Teorema 2.3 Sean (ui )i∈I ≥ 0 y (E ∗ , E 0 ) que satisfacen (H). Entonces existe un NTF de
valor (ui )i∈I con reseteo en E ∗ ⊆ E 0 .
Ki (E 0 , ui ), que resulta ser un conjunto compacto, convexo y no

Q
Demostración. Sea K =
i∈I
vacío. Note que los NTF’s son justamente los puntos fijos de la multiaplicación Γ : K → 2K
tal que a cada x0 ∈ K le entrega el conjunto no vacío, compacto y convexo:
Γ(x0 ) = y 0 ∈ K : yi,e
0
= 0 ∀i ∈ I, e ∈ E 0 tal que lw0 < ρe (lv0 , x0e )

(2.28)
con l0 las etiquetas correspondientes de x0 y E ∗ . Como la aplicación x0 → l0 es continua, Γ

es semicontinua superior y en virtud del teorema de punto fijo de Kakutani, existe un punto
fijo x0 ∈ Γ(x0 ) el cual resulta ser, por la definición, un NTF.
En general, podrían existir varios NTF’s, cada uno con sus respectivas etiquetas. El si-
guiente teorema muestra que las etiquetas en todos estos coinciden.
Teorema 2.4 Sean (ui )i∈I ≥ 0 y (E ∗ , E 0 ) que satisfacen (H). Entonces las etiquetas l0 son
las mismas para todos los NTF’s de valor (ui )i∈I con reseteo en E ∗ ⊆ E 0 .
Demostración. Sean (x0i )i∈I , (yi0 )i∈I dos NTF’s con diferentes etiquetas l0 6= h0 y suponga sin
pérdida de generalidad que el conjunto S := {v ∈P V : lv0 > h0v } es no vacío. Como x0 e
y 0 satisfacen conservación de flujo, definiendo bs = ui , bti = −ui , bv = 0 para todo otro
i∈I
vértice v, se tiene que
X X
x0 (δ + (S)) − x0 (δ − (S)) = x0e − x0e
e=vw∈δ + (S) e=vw∈δ − (S)
X
= bv
v∈S
X X
= ye0 − ye0 = y 0 (δ + (S)) − y 0 (δ − (S))
e=vw∈δ + (S) e=vw∈δ − (S)
Sea e = vw ∈ δ + (S) y suponga que x0e > ye0 . En particular, x0e > 0, luego, por ser NTF, lw0 =
ρe (lv0 , x0e ), pero ρe (lv0 , x0e ) > ρe (h0v , ye0 ) ≥ h0w por lo que se asumió. Esto es una contradicción
con que w ∈ / S, por lo tanto x0e ≤ ye0 . Sea ahora e = vw ∈ δ − (S) y suponga que ye0 > x0e .
Entonces, en particular ye0 > 0, luego h0w = ρe (h0v , ye0 ) ≥ ρe (lv0 , x0e ) ≥ lw0 , pero como w ∈ S esto
es una contradicción, por lo tanto, ye0 ≤ x0e . En resumen
21
(a) x0 (δ + (S)) − x0 (δ − (S)) = y 0 (δ + (S)) − y 0 (δ − (S))
(b) Para todo e = vw ∈ δ + (S), x0e ≤ ye0 .
(c) Para todo e = vw ∈ δ − (S), ye0 ≤ x0e .
Uniendo las tres se obtiene que x0e = ye0 para todo e ∈ δ(S) y más aún, ye0 = 0 para todo
e ∈ δ − (S) pues si ye0 > 0 se llega a una contradicción con argumentos similares a los usados
antes. Como E 0 es acíclico, existe w ∈ S tal que todo arco e = vw ∈ E 0 pertenece a δ − (S),
luego lw0 > h0w ≥ 0 pues w ∈ S y x0e = ye0 = 0 pues e ∈ δ − (S), por lo tanto para cada uno de
estos arcos se tiene que e ∈ / E ∗.
Entonces, de 2.24 se obtiene que lv0 = ρe (lv0 , x0e ) y h0v = ρe (h0v , ye0 ), luego
h0w = mı́n 0 h0v

vw∈E
lw0 = mı́n 0 lv0

vw∈E
Como cada e = vw ∈ δ − (S) cumple con que v ∈

/ S, entonces
h0w = mı́n 0 h0v ≥ mı́n 0 lv0 = lw0

vw∈E vw∈E
Pero como w ∈ S, se llega a una contradicción.
Demostrar la existencia y unicidad de los NTF es un paso crucial para la estrategia de

flujos estáticos, pues cumplen un rol protagónico en el algoritmo de extensión que se define
en la siguiente sección.
2.2. Algoritmo de extensión de flujos en equilibrio diná-

mico.
El algoritmo entregado en esta sección describe como extender un flujo factible que es
equilibrio dinámico hasta un instante θ0 , a partir de un NTF con reseteo en los arcos activos
de θ0 . Además de extender las funciones de flujo de entrada y salida en cada arco, extiende
las funciones de caminos dinámicos mínimos utilizando las etiquetas (únicas) del NTF. La
única condición que se necesita para poder utilizar el algoritmo, es que las funciones de flujo
entrante (ui )i∈I al grafo sean constantes, simplificando la integración.
2.2.1. Algoritmo de extensión.
P Dado un arreglo de flujos de entrada (ui )i∈I tal que cada ui (θ) ≡ ui , de modo que u(θ) =
ui (θ) ≡ u0 y f un flujo factible y equilibrio dinámico sobre un intervalo [0, θk ], el equilibrio
i∈I
22
se extenderá de la siguiente manera:
Algorithm 1: Extensión de flujo en equilibrio dinámico

0 ∗ 0
1 Encontrar (xi )i∈I un NTF de valor (ui )i∈I con reseteo en Eθ ⊆ Eθ , con sus etiquetas
k k
correspondientes l0 .
2 Calcular α > 0 el máximo valor tal que
/ Eθ0 k
lw (θk ) + αlw0 − lv (θk ) − αlv0 ≤ τe , ∀e = vw ∈ (2.29)
lw (θk ) + αlw0 − lv (θk ) − αlv0 ≥ τe , ∀e = vw ∈ Eθ∗k (2.30)
Definir θk+1 = θk + α.
3 Extender las funciones de caminos dinámicos mínimos y el flujo dinámico mediante:
lv (θ) = lv (θk ) + (θ − θk )lv0 , ∀v ∈ V y θ ∈ [θk , θk+1 ]

+
fi,e (ξ) = x0i,e /lv0 , ∀e = vw ∈ E, i ∈ I y ξ ∈ [lv (θk ), lv (θk+1 ))
−
x0i,e
fi,e (ξ) = ∀e = vw ∈ E, i ∈ I y ξ ∈ [Te (lv (θk )), Te (lv (θk+1 )))
ρe (lv0 , x0e )
A una extensión de f como la anterior se le denominará una α-extensión.
Observación 2.3 En el paso (3) del algoritmo, perfectamente puede ocurrir que lv0 sea cero.
En este caso, los intervalos [lv (θk ), lv (θk+1 )) y [Te (lv (θk )), Te (lv (θk+1 ))) serán vacíos, por lo
que las funciones de flujo de entrada y salida no se extienden.
Observación 2.4 Koch y Skutella (ver [13]) entregan un método para construir un equilibrio
dinámico a partir de un NTF, en el caso de un origen y un destino y un flujo de entrada cons-
tante u(θ) ≡ u0 . Larré (ver [14]) propone un algoritmo ligeramente modificado especificando
de que forma extender el flujo para arcos inactivos, además de como extender las colas en
cada iteración. El algoritmo 1 recién propuesto se diferencia de estos en primera instancia,
debido a que considera el caso de múltiples destinos. Sin embargo, incluso en el caso reducido
a un solo destino, se tiene la diferencia de que los dos algoritmos mencionados en [13] y [14]
extienden el flujo de salida por:
x0e
fe− (ξ) = ∀e = vw ∈ E, ξ ∈ [lw (θk ), lw (θk+1 ))
lw0
mientras que en el paso 3 del algoritmo 1 se pasa a extender mediante:
x0e
fe− (ξ) = ∀e = vw ∈ E, ξ ∈ [Te (lv (θk )), Te (lv (θk+1 )))
ρe (lv0 , x0e )
Las razones para estos cambios son las siguientes. Dado ξ ∈ [lv (θk ), lv (θk+1 )), para el cálculo
de las colas y para estudiar la operación bajo capacidad, se necesita evaluar fe− (ξ + τe ), sin
embargo este instante puede perfectamente ser mayor a lw (θk+1 ), por ejemplo, para un arco
inactivo sin cola justo antes de θk+1 , por lo tanto se necesitaría saber el flujo de salida en
23
un instante que no se encuentra aún definido. Por esto se realiza el cambio del intervalo de
extensión a [Te (lv (θk )), Te (lv (θk+1 ))).
La otra razón es que en el caso de Koch y Skutella si lw0 = 0 el intervalo de extensión del
flujo de salida no existe. Distinto es si se extiende el flujo de salida mediante fe− (ξ) = x0e /lw0
y se utiliza como intervalo [Te (lv (θk )), Te (lv (θk+1 ))), pues este puede ser no degenerado pero
lw0 se puede anular. Para evitar esto se cambia la definición de fe− (·) de modo que cuando
ρe (lv0 , x0e ) se anule, no exista en realidad intervalo de extensión.
Cabe destacar que para aquellos arcos con x0e > 0 se recupera la extensión dada por Koch,
Skutella y Larré en sus respectivos trabajos.
2.2.2. Correctitud del algoritmo de extensión de flujo en equilibrio

dinámico.
En la presente sección se estudia la correctitud del algoritmo 1. Para esto, se procede de la

siguiente forma: La Proposición 2.2 muestra que las colas existentes en θk no se vacían durante
[θk , θk+1 ). La Proposición 2.3 demuestra que para aquellos arcos e con ρe (lv0 , x0e ) = 0 en θk ,
no existe extensión del flujo de salida pues su intervalo de extensión es vacío. El Lema 2.3
demuestra que las extensiones lineales de los caminos dinámicos mínimos lo siguen siendo.
La Proposición 2.4 demuestra que los arcos inactivos en θk no se vuelven activos durante
[θk , θk+1 ). El Lema 2.4 estudia que sucede con los arcos que no pertenecen a los casos de las
Proposiciones 2.2 y 2.4. El Lema 2.5 demuestra que la α-extensión de un equilibrio dinámico
sigue operando bajo capacidad. Finalmente, el Teorema 2.5 cierra esta sección demostrando
que la α-extensión de un equilibrio dinámico cumple la restricción de conservación de flujo y
es efectivamente, un equilibrio dinámico.
Proposición 2.2 Sea e = vw ∈ Eθ∗k . Entonces para todo ξ ∈ [θk , θk+1 ), ze (lv (ξ)) > 0, es
decir, si un arco posee cola en θk , esta no se vaciará durante el intervalo [θk , θk + α).
Demostración. Sea ξ ∈ [θk , θk+1 ) y e = vw ∈ Eθ∗k . La cola ze (lv (ξ)) se puede calcular según
lZ
v (ξ)
ze (lv (ξ)) = ze (lv (θk )) + [fe+ (η) − fe− (η + τe )]dη

lv (θk )
Se cumple que fe (η + τe ) ≤ νe pues si η + τe ∈ [Te (lv (θk )), Te (lv (θk+1 ))) la desigualdad se sigue
−
de la extensión de cada fi,e junto con 2.24, y si η + τe < Te (lv (θk )) entonces la desigualdad
resulta de operar bajo capacidad, luego
0
xe
ze (lv (ξ)) ≥ ze (lv (θk )) + 0 − νe (lv (ξ) − lv (θk ))
lv
Pero la cola en θk cumple con νe (lw (θk ) − lv (θk ) − τe ) = ze (lv (θk )), pues e es activo. Además,
lv (ξ) − lv (θk ) = lv0 (ξ − θk ) y x0e ≥ νe lw0 nuevamente por 2.24, luego
ze (lv (ξ)) ≥ νe (lw (θk ) − lv (θk ) − τe + lw0 (ξ − θk ) − lv0 (ξ − θk )) > 0
24
donde la desigualdad mayor a 0 se obtiene gracias a que ξ − θk < α y 2.30.
Proposición 2.3 Suponga que en [0, θk ] f = (f + , f − ) es un equilibrio dinámico y que lv0

en θk es positivo, de modo que exista la α-extensión de f . Sea e = vw con ρe (lv0 , x0e ) = 0.
Entonces el flujo de salida fe− (·) no se extiende pues [Te (lv (θk )), Te (lv (θk+1 ))) = ∅.
Demostración. Como lv0 > 0 y ρe (lv0 , x0e ) = 0, de 2.24 se obtiene que necesariamente e ∈ Eθ∗k
y x0e = 0. Luego, dado θ ∈ (θk , θk+1 ), como Te (lv (θ)) = lv (θ) + τe + ze (lv (θ))/νe , y todas las
funciones involucradas son absolutamente continuas, al derivar a ambos lados
d l0
[Te (lv (θ))] = lv0 + v fe+ (lv (θ)) − fe− (lv (θ) + τe )

dθ νe
l0
= lv0 − v fe− (lv (θ) + τe )
νe
Ahora, como e ∈ Eθ∗k , utilizando 2.30 más el hecho que θ < θk+1 , se sigue que
lv (θ) + τe < lw (θ).
Por otra parte, notando que
lw0 = mı́n ρh (lv0 , x0h ) ≤ ρe (lv0 , x0e ) = 0
h=vw∈Eθ0
k
se tiene que lw0 = 0, por lo tanto lw (·) es constante en [θk , θk+1 ), por lo tanto lw (θk ) = lw (θ),
luego
lv (θ) + τe < lw (θ) = lw (θk ) ≤ Te (lv (θk ))
Como e ∈ Eθ∗k , en virtud de la proposición anterior, e ∈ Eθ∗ y como en [0, θk ] se opera bajo
d
capacidad, fe− (lv (θ) + τe ) = νe , concluyendo que dθ [Te (lv (θ))] = 0. Por lo tanto
[Te (lv (θk )), Te (lv (θk+1 ))) = ∅
es decir, el flujo de salida en realidad no se extiende.
Lema 2.3 Sea f un flujo factible y equilibrio dinámico en el intervalo de tiempo [0, θk ]. Sea
α > 0 un real positivo que cumple 2.29 y 2.30. Entonces la extensión de las funciones de
caminos dinámicos mínimos
ˆlv (θ) = lv (θk ) + (θ − θk )l0 , ∀v ∈ V y θ ∈ [θk , θk+1 ]
v
al intervalo [0, θk+1 ] definen caminos dinámicos mínimos.
Demostración. Se quiere ver que para θ ∈ (θk , θk+1 ) se sigue cumpliendo 2.15. Es sencillo
ver que para v = s, ˆlv (θ) = θ simplemente utilizando que lv (θk ) = θk y que ls0 = 1. Para
w 6= s, dado e = vw se tiene que Te (ˆlv (θ)) = ˆlv (θ) + τe + ze (l̂νve(θ)) y como todas las funciones
involucradas son absolutamente continuas, derivando a ambos lados
d h ˆ i ˆl0 h i
Te (lv (θ)) = ˆlv0 + v fe+ (ˆlv (θ)) − fe− (ˆlv (θ) + τe )
dθ νe
25
Suponga que ˆlv (θ) + τe ∈ [Te (lv (θk )), Te (lv (θk+1 ))), entonces
d h ˆ i ˆl0 x0 x0e

ˆ
Te (lv (θ)) = lv +0 v e
−
dθ νe lv0 ρe (lv0 , x0e )
  ∗

ν e si e ∈ Eθ
ˆl0  x0

 k
= ˆlv0 + v  0e −

νe lv
n 0 o 
 mı́n x0e , νe

si e ∈ ∗
/ Eθk
l v
Notando que ˆlv0 = lv0 ,

x0e
 
 νe
si e ∈ Eθ∗k 
d h ˆ i  
0 0
Te (lv (θ)) = = ρe (lv , xe )
dθ n
0 xe
0
o
∗ 
máx lv , νe si e ∈
/ Eθk

 
Por otro lado, si ˆlv (θ) + τe < Te (lv (θk )), necesariamente se debe cumplir que e ∈ Eθ∗k pues
/ Eθ∗k entonces,
si e ∈
ˆlv (θ) + τe ≥ ˆlv (θk ) + τe = lv (θk ) + τe = Te (lv (θk ))
lo que es una contradicción. Luego, como en [0, θk ] se opera bajo capacidad, fe− (ˆlv (θ)+τe ) = νe ,
se sigue que
d h ˆ i ˆl0 x0
Te (lv (θ)) = ˆlv +
0 v e
− νe
dθ νe lv0
ˆ0 x0e ˆ0
= lv + − lv
νe
x0
= e = ρe (lv0 , x0e )
νe
donde la última igualdad sale nuevamente del hecho que e ∈ Eθ∗k . Usando la regla de derivación
de un mínimo de funciones se obtiene
" #
d
mı́n Te (ˆlv (θ)) = mı́n 0 ρe (lv0 , x0e ).
dθ e=vw∈Eθ0 k e=vw∈Eθ
k
Como en [0, θk ] las etiquetas son caminos dinámicos mínimos, se tiene que lw0 = mı́n ρe (lv0 , x0e ),
e=vw∈Eθ0
k
luego " #
d
mı́n Te (ˆlv (θ)) = lw0 = ˆlw0
dθ e=vw∈Eθ0
k
es decir, las funciones mı́n Te (ˆlv (·)) y ˆlw (·) varían de la misma forma. Como en θk las
e=vw∈Eθ0
k
funciones son iguales pues ˆl y l coinciden en ese instante, se concluye que
mı́n Te (ˆlv (θ)) = ˆlw (θ).

e=vw∈Eθ0
k
26
Finalmente, como para un arco e ∈ / Eθ0 k no se alcanzará el mínimo anterior pues de 2.30
se tiene que ˆlw (θ) < τe + ˆlv (θ) ≤ Te (ˆlv (θ)), donde la desigualdad estricta es debido a que
θ < θk+1 , se puede reescribir el mínimo anterior como
mı́n Te (ˆlv (θ)) = ˆlw (θ)

e=vw∈E
concluyendo que 2.15 aún se cumple.
/ Eθ0 k . Entonces para todo ξ ∈ [θk , θk+1 ), se tiene que e ∈

Proposición 2.4 Sea e = vw ∈ / Eξ0 .
Demostración. / Eθ0 k y ξ ∈ [θk , θk+1 ). De 2.29 se tiene

Sea e = vw ∈
lw (θk ) + (ξ − θk )lw0 − lv (θk ) − (ξ − θk )lv0 < τe
donde la desigualdad estricta sale de que ξ − θk < α, luego
lw (ξ) < lv (ξ) + τe ≤ Te (lv (ξ))
/ Eξ0 .
concluyendo así que e ∈
Según lo visto en las proposiciones 2.2 y 2.4, durante [θk , θk+1 ), ningún arco con cola verá
esta vaciada, así como ningún arco inactivo se volverá activo. Sin embargo, hay dos casos que
si se pueden dar:
(i) Existe e ∈ Eθ0 k \ Eθ∗k tal que e ∈

/ Eθ0 para θ > θk .
(ii) Existe e ∈ Eθ0 k \ Eθ∗k tal que e ∈ Eθ∗ para θ > θk .
es decir, un arco activo sin cola se vuelve inactivo o bien forma una cola durante el intervalo
[θk , θk+1 ). Afortunadamente, en ambas situaciones el NTF encontrado en el paso 1 no se ve
modificado con esta variación de los conjuntos Eθ0 k , Eθ∗k , por lo que no es necesario recalcularlo
hasta θk+1 . Esto es demostrado en el siguiente lema.
α > 0 un real positivo que cumple 2.29 y 2.30. Sea e ∈ Eθ0 k \ Eθ∗k tal que, para θ ∈ (θk , θk+1 ),
e∈/ Eθ0 o e ∈ Eθ∗ . Entonces el NTF (x0i )i∈I encontrado en 1 no se ve modificado.
Demostración. Suponga que e = vw ∈ Eθ∗ , es decir, ze (θk ) = 0 y ze (θ) > 0 para θ > θk .
Como e comienza a armar una cola justo en θ, debe cumplirse que fe+ (lv (θ)) > νe de modo
que no pueda pasar todo el flujo que está llegando. Es decir, x0e /lv0 > νe , o equivalentemente
x0e /νe > lv0 . Como e ∈
/ Eθ∗k , por 2.24 se tiene
lw0 = máx{lv0 , x0e /νe } = x0e /νe .
y como ze (θ) > 0, la ecuación correspondiente al NTF en el instante θ sería
lw0 = x0e /νe .
por lo tanto el sistema a resolver no cambia en e, por lo que tampoco lo hace el NTF.
27
Por otro lado, si e ∈ / Eθ0 para θ > θk , al calcular el NTF en θ se debe imponer x0e (θ) = 0.
Ahora bien, si xe (θk ) > 0, como x0 es NTF en θk , lw0 = ρe (lv0 , x0e ). Además, e ∈
0
/ Eθ∗k luego
/ Eθ0 , lw (θ) < lv (θ) + τe , gracias al Lema 2.3
lw (θk ) = Te (lv (θk )) = lv (θk ) + τe , pero como e ∈
que demuestra que las extensiones lineales siguen definiendo caminos dinámicos mínimos, por
lo tanto necesariamente se debe cumplir que lw0 < lv0 , luego ρe (lv0 , x0e ) < lv0 , lo que contradice
2.24 en el caso ze (θk ) = 0. Por lo tanto x0e = 0.
α > 0 un real positivo que cumple 2.29 y 2.30. Para todo e = vw ∈ E y para casi todo
ξ ∈ [θk , θk + α) se tendrá que

 νe si ze (lv (ξ)) > 0,
−
fe (lv (ξ) + τe ) =
mı́n {fe+ (lv (ξ)), νe } si no.

Demostración. Sea e = vw ∈ E, ξ ∈ [θk , θk+1 ). Si lv (ξ)+τe ≤ Te (lv (θk )) entonces la propiedad

se obtiene del hecho que hasta lv (θk ) el flujo opera bajo capacidad. Suponga que lv (ξ) + τe ∈
(Te (lv (θk )), Te (lv (θk+1 )). A partir de 2.24 se obtiene

x0i,e


 x0e e
ν si e ∈ Eθ∗k
−
fi,e (lv (ξ) + τe ) = n 0 x0i,e
o
 mı́n xi,e / Eθ∗k

 ν ,
x0 e l 0
si e ∈
e v
Luego, sumando sobre todo i ∈ I,



 νe si e ∈ Eθ∗k
fe− (lv (ξ) + τe ) = n o
 mı́n νe , x00e

/ Eθ∗k
si e ∈
l v
Suponga que e ∈ Eξ∗ . Si e ∈ Eθ∗k entonces fe− (lv (ξ) + τe ) es de la forma buscada. Por otro
lado, si e ∈ Eθ0 k \ Eθ∗k , necesariamente debe cumplirse que x0e /lv0 ≥ νe para que pueda formarse
n o
− x0e
cola, luego fe (lv (ξ) + τe ) = mı́n νe , l0 = νe .
v
Si e ∈ Eξ0 \ Eξ∗ , como las colas no se vacían durante [θk , θk+1 ), se cumple que e ∈
/ Eθ∗k , por
n 0
o
lo tanto fe− (lv (ξ) + τe ) = mı́n νe , xl0e .
v
Teorema 2.5 Sea f un flujo factible y equilibrio dinámico en el intervalo de tiempo [0, θk ].
Sea α > 0 un real positivo que cumple 2.29 y 2.30. Entonces la α-extensión de f es un flujo
factible y equilibrio dinámico en [0, θk+1 ].
Demostración. Lo primero es probar que 2.17 se sigue cumpliendo. Sea v ∈ V , y lv0 > 0, de
modo que tenga sentido hablar de extensión. Sea ξ ∈ [lv (θk ), lv (θk+1 )). Como (x0i )i∈I es un
NTF, en particular es un flujo estático de valor (ui )i∈I , por lo tanto, dado i ∈ I,

X
0
X
0 ui (ξ) si v = s
xi,e − xi,e =
+ −
0 si v 6= {s, ti }
28
Luego, dividiendo por lv0 , recordando que ls0 = 1:

X X ui (ξ) si v = s
x0i,e /lv0 − x0i,e /lv0 =
0 si v =
6 {s, ti }
e∈δ + (v) e∈δ − (v)
Sea e = uv ∈ δ − (v). Si x0i,e = 0, entonces el término realmente no aparece en la suma. Por

otro lado, si x0i,e > 0, se tiene que lv0 = ρe (lu0 , x0e ), por lo tanto

X
0 0
X
0 0 0 ui (ξ) si v = s
xi,e /lv − xi,e /ρe (lu , xe ) = (2.31)
+ −
0 si v 6 = {s, ti }
e∈δ (v) e=uv∈δ (v)
Además, si x0i,e > 0, en particular e es un arco activo en θk , por lo tanto lv (θk ) = Te (lu (θk )),
+ −
con lo cual ξ ∈ [Te (lu (θk )), Te (lu (θk+1 ))), luego por como se definen fi,e (·) y fi,e (·) en la
α-extensión, se tiene que 2.31 queda

X
+
X
− ui (ξ) si v = s
fi,e (ξ) − fi,e (ξ) = (2.32)
+ −
0 si v 6= {s, ti }
probando lo buscado.
Para ver que el flujo extendido sigue siendo un equilibrio dinámico, se verá que la parte
(c) del Teorema 2.1 se cumple. Sea i ∈ I, ξ ∈ [0, θk+1 ) y e ∈ / Eξ0 . Si ξ ≤ θk entonces
+ 0
fi,e (lv (ξ))lv (ξ) = 0 pues hasta θk el flujo es un equilibrio dinámico. Por otra parte, si ξ ∈
(θk , θk+1 ) se consideran los siguientes casos:
1. Si lv0 = 0 en [θk , θk+1 ), entonces trivialmente se cumple lo buscado.

2. Si lv0 > 0 en [θk , θk+1 ) y e ∈/ Eθ0 k entonces xi,e0 = 0, y por la forma como se extiende el
+
flujo de entrada se obtiene nuevamente que fi,e (lv (ξ))lv0 (ξ) = 0.
3. Si lv0 > 0 en [θk , θk+1 ) y e ∈ Eθ0 k , el lema 2.4 entrega que el NTF no cambia al pasar el
arco activo en θk a ser inactivo en ξ, y además x0i,e = 0 en θk , por lo tanto, nuevamente
+
fi,e (lv (ξ))lv0 (ξ) = 0.
De aquí se concluye la correctitud del algoritmo de extensión de flujos en equilibrio diná-

mico.
Observación 2.5 Es importante notar que en el paso 2 del algoritmo 1, α se puede elegir
estrictamente mayor que cero. Para ver esto, si e ∈ / Eθ0 k , entonces lw (θk ) < lv (θk ) + τe
y por continuidad existe α1 > 0 tal que 2.29 se cumple para todo α ≤ α1 . Asimismo, si
e = vw ∈ Eθ∗k , se tiene que lw (θk ) − lv (θk ) − τe = ze (lv (θk ))/νe > 0 y entonces existe α2 > 0
tal que 2.30 se cumple para todo α ≤ α2 . En virtud de lo anterior, existe α > 0 que cumple
con ambas desigualdades al mismo tiempo.
Los Teoremas 2.3 y 2.4 aseguran la existencia de un NTF para el paso 1 de cada iteración
y la unicidad de sus etiquetas.
En virtud de disminuir el trabajo computacional se recomienda siempre considerar la má-

xima constante α que cumpla con ambas desigualdades, de modo de maximizar el crecimiento
del intervalo [0, θk ] minimizando así el número de iteraciones.
29
Observación 2.6 El trabajo hecho por Koch y Skutella en [13], además de las diferencias
explicadas al comienzo con el algoritmo propuesto en la sección anterior, tiene problemas en
las demostraciones hechas para ver la correctitud de este. Primero, no se entregan demostra-
ciones formales de que las desigualdades 2.29 y 2.30 implican que los arcos no se activan ni
las colas no se vacían respectivamente, únicamente se menciona. Segundo, para estudiar la
operación bajo capacidad de la extensión, se utiliza la caracterización dada en la observación
2.1, sin embargo solo se demuestran las partes (a) y (b) dejando de lado la condición de
tiempo de espera. Un estudio formal de esta última da a conocer el problema que puede surgir
al hacer el cálculo del largo de cola como se explicó en la observación 2.4, haciendo necesario
el cambio del intervalo de extensión para el flujo de salida. En detalle, el problema surge al
estudiar aquellos arcos que poseen x0e = 0.
Larré, en [14], presenta los mismos problemas en sus demostraciones. Si bien demuestra
que aquellos arcos con x0e > 0 satisfacen la operación bajo capacidad, no discute en detalle
que sucede con los otros, limitándose al igual que Koch y Skutella a decir que estos arcos
cumplen trivialmente las propiedades, sin embargo, y como ya se mencionó anteriormente,
es necesario calcular sus colas lo que puede llevar a evaluar el flujo de salida en instantes de
tiempo donde aún no se ha extendido.
Debido a los problemas presentes en los estudios hechos anteriormente al algoritmo 1,

además de la falta de demostración de ciertas propiedades importantes, se redemostraron
todos los resultados. El orden que siguen las demostraciones hechas es de gran importancia
ya que solo se utilizan resultados previos evitando cualquier tipo de argumento circular.
2.3. Cálculo de un flujo fino normalizado.
Como se puede ver en la sección anterior, uno de los elementos necesarios en el Algoritmo
1 es la obtención de un NTF. La siguiente sección revisa un algoritmo, basado en el trabajo
de R. Koch en [12], sección 6. Se comienza dando definiciones y resultados preliminares
importantes en el cálculo, para luego entregar el algoritmo.
Desde ahora en adelante, sea la tupla (G, ν, s, {ti }i∈I ) dada por G := (V, E) un grafo
dirigido, ν las capacidades de los arcos, s ∈ V un nodo fuente y {ti }i∈I ⊆ V \ {s} un conjunto
de nodos destino.
Definición 2.12 A un conjunto X ⊆ V se le denominará corte. Si s ∈ X, se dirá s-corte,

y si además ∃t ∈ V \ X, entonces X se dirá s-t-corte.
Definición 2.13 Se entenderá por balance de nodos a cualquier vector b ∈ RV tal que
bs > 0, bti < 0, ∀i ∈ I y bv = 0 para el resto de nodos v ∈ V , que además cumpla
X
b(V ) := bv = 0.
v∈V
Definición 2.14 Dado X ⊆ V un corte y b ∈ RV un balance de nodos, la congestión de X
30
estará dada por
P
bv
b(X) v∈X
q := = P .
ν(δ + (X)) νe
e∈δ + (X)
Además, un corte de congestión máxima se denominará corte disperso o ralo.
Definición 2.15 Sea y := (yi )i∈I un arreglo tal que, para cada iP∈ I, yi ∈ RE+ . Al vector y
se le llamará flujo estático, y a partir de este se definirá ye := yi,e para cada e ∈ E. Un
i∈I
elemento yi,e de un flujo estático corresponde a la fracción de flujo que se dirige al destino
ti que está pasando por el arco e. Por otro lado, ye corresponde al flujo total que circula por
el arco. Si se consideran parejas de valores (bsi , bti ) para todo i ∈ I que simbolizan el flujo
enviado desde s hacia ti , se dice que y es factible si cumple
yi (δ + (s)) − yi (δ − (s)) = bsi , ∀i ∈ I

yi (δ + (ti )) − yi (δ − (ti )) = bti , ∀i ∈ I
yi (δ + (tj )) − yi (δ − (tj )) = 0, ∀i 6= j
yi (δ + (v)) − yi (δ − (v)) = 0, ∀i ∈ I, ∀v ∈ V : v 6= s, v 6= {tj }j∈I
P
Del mismo modo, definiendo bs := bsi se tendrá que el vector (ye )e∈E es factible o que
i∈I
respeta el balance de nodos b si
y(δ + (s)) − y(δ − (s)) = bs ,

y(δ + (ti )) − y(δ − (ti )) = bti , ∀i ∈ I
y(δ + (v)) − y(δ − (v)) = 0, ∀v ∈ V : v 6= s, v 6= {tj }j∈I
Proposición 2.5 La congestión q ∗ de un corte ralo en una red (G, ν, s, {ti }i∈I ) puede ser
calculada mediante el siguiente programa lineal:
mı́n q ∗
subject to :
y(δ + (s)) − y(δ − (s)) = bs ,
y(δ + (ti )) − y(δ − (ti )) = bti , ∀i ∈ I
y(δ + (v)) − y(δ − (v)) = 0, ∀v ∈ V : v 6= s, v 6= {tj }j∈I
y ≤ ν · q∗
y ∈ RE
+
Este problema de optimización lineal entrega un flujo factible que viaja por la red satisfa-
ciendo las demandas bs > 0 y bti < 0, ∀i ∈ I. A partir de este se puede obtener sencillamente
una descomposición por destinos (yi,e )i∈I,e∈E tal como se explica en la siguiente rutina. Sea
(ye )e∈E un flujo total obtenido según la proposición 2.5 y T := (ti )i∈I el conjunto de nodos
destinos.
31
Algorithm 2: Descomposición de un flujo por destinos.
1 Sea ti ∈ T y considere Pi el conjunto de todos los s-ti caminos con flujo positivo en
todos sus arcos.
2 Para cada P ∈ Pi , se calcula λP := mı́n{ye : e ∈ P }, se define yi,e = λP , ∀e ∈ P y se
actualiza ye := ye − λp para todo e ∈ P .
3 Se setea T := T \ {ti }.
4 Si T 6= ∅, se vuelve a 1. Sino, se define yi,e = 0 para toda pareja (i, e) que no se haya
definido aún y se retorna y := (yi,e )i∈I,e∈E .
Definición 2.16 Sea y un flujo estático. Sea además X un corte ralo y q ∗ su congestión. Se
dirá que y respeta el valor del corte ralo si se satisface
ye ≤ νe · q ∗ , ∀e ∈ E
Definición 2.17 Sea una tupla (G, ν, s, {ti }i∈I ) y E1 ⊆ E un conjunto de arcos. Se dirá red
con reseteo a la red (G, ν, s, {ti }i∈I , E1 ) y a un arco e ∈ E1 , arco con reseteo.
Definición 2.18 Sea (G, ν, s, {ti }i∈I , E1 ) una red con reseteo e y = (yi )i∈I un flujo estático
que respeta un balance de nodos b. Se dirá que y es un flujo fino con reseteo si existe un
vector h := (hv )v∈V que cumple con:
hs = 1
hw ≤ hv ∀e = vw ∈ E \ E1 , ye = 0
hw = máx{hv , ye /νe } ∀e = vw ∈ E \ E1 , ye > 0
hw = ye /νe ∀e = vw ∈ E1
A los elementos del vector h se les dirá etiquetas. En caso que E1 sea vacío, y se dirá
flujo fino sin reseteo.
Observación 2.7 Esta definición de flujo fino corresponde a la de un NTF al considerar la

red con reseteo (G, ν, s, {ti }i∈I , Eθ∗ ) y el balance de nodos bsi = ui , bti = −ui , ∀i ∈ I, donde
las etiquetas corresponden a lv0 .
Proposición 2.6 Sea b un balance de nodos, y = (yi )i∈I un flujo fino sin reseteo que respeta
a b y a la congestión q ∗ de un corte ralo X ⊆ V . Sean (hv )v∈V sus respectivas etiquetas.
Entonces, para todo nodo v ∈ V \ X, su etiqueta hv vale exactamente q ∗ .
Definición 2.19 Sea y un flujo estático que respeta un balance de nodos b. Sea además
X ⊆ V un corte. Se define la red inducida por X como la red (H, ν H , bH ) dada por el
subgrafo inducido H := G[X], ν H la restricción de ν a H y el nuevo balance de nodos bH
dado por
bH + + − −
v := bv − y(δ (v) ∩ δ (X)) + y(δ (v) ∩ δ (X)), ∀v ∈ X
En el caso que X sea un corte ralo de congestión q, en virtud de la proposición 2.6, el nuevo
balance de nodos simplemente queda
bH + +
v := bv − q · ν(δ (v) ∩ δ (X)), ∀v ∈ X
32
Proposición 2.7 Sea y un flujo estático que respeta un balance de nodos b, e y H su restric-
ción como vector a una red inducida (H, ν H , bH ) por X un corte. Entonces y H es un flujo
estático que respeta a bH .
Definición 2.20 Sea (G, ν, s, {ti }i∈I , E1 ) una red con reseteo y e = vw ∈ E1 un arco con
reseteo. Se construye una nueva red con reseteo mediante borrar el arco e = vw ∈ E1 de G y
agregar el arco es := sw con capacidad νe . Realizando este procedimiento para cada arco de
E1 se obtiene una red sin reseteo llamada red anclada.
La siguiente proposición explica como a partir de un flujo fino sin reseteo se puede obtener
uno con reseteo y viceversa.
Proposición 2.8 Sea N := (G, ν, s, {ti }i∈I , E1 ) una red con reseteo y Na := (Ga , ν a , s, {ti }i∈I )
la red anclada correspondiente. Sea además y un flujo fino con reseteo en la red N que respeta
cierto balance de nodos b. Si se define
ye0 s := ye ∀e ∈ E1 ,
ye0 := ye ∀e ∈ E \ E1 ,
b0s := bs + y(E1 ),
X
b0v := bv − ye ∀v ∈ V \ {s}.
e∈δ + (v)∩E1
se obtiene que y 0 es un flujo fino sin reseteo en la red anclada Na que respeta el balance de
nodos b0 . Viceversa, si y 0 es un flujo fino sin reseteo en Na respetando el balance de nodos b0 ,
al definir
ye := ye0 s ∀e ∈ E1 ,
ye := ye0 ∀e ∈ E \ E1 ,
X
bs := bs − ye0 s ,
e∈E1
X
bv := b0v + ye0 s ∀v ∈ V \ {s}.
e∈δ + (v)∩E1
se obtiene que y es un flujo fino con reseteo en N que respeta el balance de nodos b.
Finalmente, el siguiente algoritmo entrega una forma de encontrar un flujo fino sin reseteo,
es decir, cuando E1 = ∅, dada una red (G, ν, s, {ti }i∈I ) y b un balance de nodos.
33
Algorithm 3: Flujo fino sin reseteo
1 Se inicializa J := (G, ν, s, {ti }i∈I ), Jold := ∅.
∗ 0
2 Encontrar un corte ralo X en I de congestión q junto con un flujo estático y que lo
respete utilizando la proposición 2.5.
0
3 Dado X := V (G) \ X, definir
hv := q ∗ , ∀v ∈ X 0
+
ye := ye0 , ∀e ∈ δG (X) ∪ E(G[X 0 ])
4 Asignar Jold := J . Construir la red inducida por X y asignarla a J . Calcular el nuevo

balance de nodos en J .
5 Si J es igual que Jold , definir hv := 0, ∀v ∈ V (G) y retornar (y, h). Sino ir a 2.
Utilizando la Proposición 2.8, el Algoritmo 3 también sirve para encontrar flujos finos en
una red con reseteo. Sin embargo existe una complicación, pues se necesita saber el valor del
flujo fino en cada arco para poder pasar de una red con reseteo a una sin, y luego utilizar
el algoritmo. Esto es obviamente una contradicción pues se necesita utilizar el objeto que
se busca calcular. Para poder resolver este problema, se realiza una rutina de punto fijo.
Sea y0 un flujo fino con reseteo dado y b0 un balance de nodos en una red con reseteo
(G, ν, s, {ti }i∈I , E1 ). Sea Ga el grafo anclado dado por G y E1 y considere y := y0 , b := b0 .
Algorithm 4: Rutina de punto fijo para el cálculo de un flujo fino con reseteo.
0 0
1 Sean (y , b ) el flujo fino sin reseteo y el balance de nodos definidos por (y, b) según la
Proposición 2.8.
0 0
2 Sea z un flujo fino sin reseteo en la red (Ga , ν, s, {ti }i∈I ) con b como balance de nodos
entregado por el algoritmo 3.
0 0
3 Se calcula z el flujo fino con reseteo y b̂ el balance de nodos definidos por z y b en la
red (G, ν, s, {ti }i∈I , E1 ) según la Proposición 2.8.
4 Si b = b̂, se retorna z. Si no, se asigna y := z y se vuelve a 1.
Esta rutina logra encontrar, en general en pocas iteraciones, un flujo fino con reseteo.
Como primer flujo y0 se utiliza el valor del flujo fino que poseen los arcos en la iteración
anterior del Algoritmo 1, pues no hay que olvidar que la rutina de punto fijo se usará en cada
iteración del algoritmo de extensión de flujo en equilibrio dinámico y el flujo fino con reseteo
encontrando corresponde a un NTF cuando se utiliza la red (G, ν, s, {ti }i∈I , Eθ∗k ).
34
Capítulo 3
Modelo de cola fluida, caso múltiples

orígenes y destinos.
En este capítulo se muestra brevemente el modelo teórico para una red de flujo dinámico
con múltiples nodos origen y destino. Se explican las complicaciones presentes en el modelo
y que no permiten trabajar mediante una estrategia de flujos estáticos como se hizo en el
capítulo anterior. Finalmente, se muestran ejemplos que estudian la evolución del flujo a
través de un grafo, dejando en evidencia otras dificultades que implica el tener múltiples
fuentes.
3.1. Modelo de red de flujo dinámico.
El modelo teórico de una red de flujo con múltiples orígenes y destinos sigue la misma
línea que el de un solo origen.
3.1.1. Redes de flujo.
Sea G = (V, E) un grafo dirigido con V un conjunto de nodos y E un conjunto de arcos.

Sea (si , ti )i∈I ⊆ V una familia finita de parejas de nodos, donde si corresponde a un nodo
origen y ti a su respectivo nodo destino.
Las nociones de flujo sobre el tiempo, flujo total entrante a cada arco, flujo total saliente
de cada arco, flujo entrante acumulativo y flujo saliente acumulativo se definen de manera
análoga al caso de un solo origen, salvo que ahora corresponden al flujo que proviene del nodo
si y se dirige al nodo ti , para cada i ∈ I.
La función Largo de cola, operación bajo capacidad y la función de tiempo de viaje por un
arco se definen de la misma forma. Es importante notar que si se considera la derivada de la
35
función de tiempo de viaje por un arco, se obtiene nuevamente (2.11)
1 +


 f (θ)
νe e
si ze (θ) > 0,
0
Te (θ) = n o
 máx 1, 1 f + (θ)

si no.
νe e
de donde se puede concluir que Te0 ≥ 0 y por lo tanto que las partículas respetan FIFO.
Al igual que en el capítulo anterior, si una partícula proveniente de si ingresa en θ a un

camino P = (e1 , e2 , ..., ek ), saldrá de este en
lP (θ) := Tek ◦ · · · ◦ Te1 (θ)
Considerando Pi,w el conjunto de todos los si -w caminos en G, el mínimo tiempo en que una
partícula saliendo de si en θ puede llegar a w ∈ V está dado por
li,w (θ) := mı́n lP (θ)

P ∈Pi,w
Estas funciones corresponden a los i-caminos dinámicos mínimos.
Al igual que el caso de un solo origen, las funciones de i-caminos dinámicos mínimos y
tiempo de viaje por un arco cumplen con Te (θ) → ±∞, li,w (θ) → ±∞ cuando θ → ±∞, son
continuas y sobreyectivas.
Además se puede encontrar una recurrencia al asumir que la suma de latencias de cualquier
ciclo es siempre positiva:
(
θ si w = si ,
li,w (θ) = mı́n Te (li,v (θ)) si w 6= si .
e=vw∈E
y con esto, se definen los conjuntos de arcos i-activos en el instante θ además de los arcos
con cola mediante
0
Ei,θ := {e = vw ∈ E : Te (li,v (θ)) ≤ li,w (θ)} ,
∗
Ei,θ := {e = vw ∈ E : ze (li,v (θ)) > 0} .
3.1.2. Equilibrio dinámico.
Se define el conjunto de instantes en que un arco e es i-activo como:

0
Θi,e := {θ ∈ R : e ∈ Ei,θ }.
Además, un flujo f se dirá factible si cumple con

X
+
X
− ui (θ) si v = si
fi,e (θ) − fi,e (θ) =
0 si v 6= si
e∈δ + (v) e∈δ − (v)
con ui la tasa de flujo entrante por el nodo si al grafo.
36
Definición 3.1 Un flujo factible f se dirá Equilibrio Dinámico si para cada e = vw ∈ E
+
se cumple fi,e (ξ) = 0 para ξ ∈ li,v (Θci,e ) c.t.p.
Hay tres resultados que se recuperan del capítulo anterior, los Lemas 2.1 y 2.2, y el
Teorema 2.1. Para el primer lema, la demostración hecha en el caso de un solo origen se
ajusta perfectamente a este caso, por lo que se omite.
Lema 3.1 Sea un flujo de entrada descompuesto en los pares de nodos origen-destino
X
+
fe+ (θ) = fi,e (θ)
i∈I
− −
definidos para todo θ ∈ R. Entonces existen flujos únicos fi,e ∈ L∞ (R) con 0 ≤ fi,e (ξ) ≤
νe , ∀ξ ∈ R que satisfacen
TZe (θ) Zθ
− +
0 0
para todo θ ∈ R, i ∈ I, e ∈ E.
Lema 3.2 Un flujo sobre el tiempo factible f es un equilibrio dinámico si y solo si ∀e =

+
vw ∈ E, i ∈ I, se cumple fi,e (ξ) > 0 ⇒ ξ ∈ li,v (Θi,e ) para casi todo ξ ∈ R.
+
Demostración. Sean i ∈ I, e ∈ E fijos. Notar que fi,e (ξ) > 0 ⇒ ξ ∈ li,v (Θi,e ) c.t.p. es
+
equivalente a fi,e (ξ) = 0 para ξ ∈ li,v (Θi,e )c c.t.p. Luego, basta con ver que los conjuntos
li,v (Θci,e ) y li,v (Θi,e )c difieren en un conjunto de medida nula. Sea ξ ∈ li,v (Θi,e )c , como la
función li,v (·) es sobreyectiva, existe θ ∈ R tal que ξ = li,v (θ). Se cumple que θ ∈ Θci,e ya
que sino, li,v (θ) ∈ li,v (Θi,e ) es decir ξ ∈ li,v (Θi,e ) lo que es una contradicción. Por lo tanto
li,v (θ) ∈ li,v (Θci,e ) con lo cual se concluye que li,v (Θi,e )c ⊆ li,v (Θci,e ).
Sea ahora ξ ∈ li,v (Θci,e ) \ li,v (Θi,e )c = li,v (Θci,e ) ∩ li,v (Θi,e ), entonces existen θ, θ̄ con θ ∈ Θci,e ,
θ̄ ∈ Θi,e tal que ξ = li,v (θ) = li,v (θ̄). Como li,v (·) es no decreciente, en realidad li,v (η) = ξ para
todo η entre θ y θ̄, los cuales, como viven en conjuntos disjuntos son distintos, entonces en
particular existe η ∈ Q tal que ξ = li,v (η), es decir ξ ∈ li,v (Q) que es un conjunto numerable
y por lo tanto de medida nula. Se concluye entonces que li,v (Θci,e ) \ li,v (Θi,e )c está contenido
en un conjunto de medida nula.
Teorema 3.1 Sea f un flujo sobre el tiempo factible. Son equivalentes:
(a) f es un equilibrio dinámico.

+ −
(b) Para todo e = vw ∈ E, i ∈ I y para todo θ ∈ R se tiene Fi,e (li,v (θ)) = Fi,e (li,w (θ)).
0 + 0
(c) Para todo e = vw ∈ E, i ∈ I y para casi todo θ ∈ R, e ∈
/ Ei,θ ⇒ fi,e (li,v (θ))li,v (θ) = 0.
Demostración. Sean i ∈ I, e = vw ∈ E y θ ∈ R. Sea θ0 ≤ θ el máximo instante tal que

Te (li,v (θ0 )) = li,w (θ). Notar que θ0 existe gracias a que, si e es i-activo para θ entonces θ0 = θ,
por otro lado si e es i-inactivo, li,w (θ) < Te (li,v (θ)), y como Te (li,v (θ)) → −∞ cuando θ → −∞
y Te (li,v (·)) es continua, existe un θ0 que cumple con lo pedido.
37
Sea Iθ = (θ0 , θ] (notar que este intervalo también depende de i, sin embargo se omite el
índice). Cabe decir que θ ∈ Θi,e si y solo si θ0 = θ con lo cual Iθ = ∅. Lo primero es demostrar
que Θci,e = ∪θ Iθ . En efecto. si θ ∈ Θci,e , se tiene que θ0 < θ y entonces Iθ 6= ∅ con lo cual θ ∈ Iθ
lo que entrega una primera inclusión. Recíprocamente, sea θ00 ∈ Iθ . Por la definición de θ0 y
la monotonía de Te (·), li,w (·) se cumple que
li,w (θ00 ) ≤ li,w (θ) = Te (li,v (θ0 )) < Te (li,v (θ00 ))
donde la última desigualdad debe ser estricta pues si Te (li,v (θ0 )) = Te (li,v (θ00 )) se contradice
que θ0 < θ00 y la definición de θ0 . Luego, θ00 ∈ Θci,e concluyendo la otra inclusión.
Utilizando el Lema 3.1 se cumple para todo θ ∈ R que
TZe (θ) Zθ
− +
0 0
por lo tanto
li,v
Z (θ)
+ − + − + + +
Fi,e (li,v (θ))−Fi,e (li,w (θ)) = Fi,e (li,v (θ))−Fi,e (Te (li,v (θ0 ))) = Fi,e (li,v (θ))−Fi,e (li,v (θ0 )) = fi,e (ξ)dξ
li,v (θ0 )
+
donde esta integral es siempre mayor o igual a 0 pues fi,e ∈ F0 (R). Luego,
+ −
Fi,e (li,v (θ)) − Fi,e (li,w (θ)) = 0
+
si y solo si fi,e se desvanece c.t.p en el intervalo (li,v (θ0 ), li,v (θ)] = li,v (Iθ ), es decir, en virtud
+
del Lema 1.1, si y solo si fi,e (ξ) = 0 para ξ ∈ ∪θ li,v (Iθ ) = li,v (Θce ) c.t.p, lo que prueba que (a)
y (b) son equivalentes.
Ahora, partiendo desde
li,v
Z (θ)
+ − +
Fi,e (li,v (θ)) − Fi,e (li,w (θ)) = fi,e (ξ)dξ
li,v (θ0 )
y realizando un cambio de variable el cual vale gracias a que las funciones li,v son absoluta-
mente continuas y monótonas, se obtiene
Zθ
+ − + 0
Fi,e (li,v (θ)) − Fi,e (li,w (θ)) = fi,e (li,v (ξ))li,v (ξ)dξ
θ0
+ 0
Al igual que antes esta integral es igual a 0 si y solo si fi,e (li,v (ξ))li,v (ξ) = 0 para todo ξ ∈ Iθ
c.t.p, es decir, nuevamente en virtud del Lema 1.1, si y solo si esta función se desvanece en
∪θ Iθ = Θci,e c.t.p, probando la equivalencia entre (b) y (c).
38
3.2. Complicaciones en el modelo teórico.
Existen varias razones por las cuales el modelo recién presentado no es fácil de analizar.
0 ∗
Uno de los mayores problemas ocurre en la relación que existe entre los conjuntos Ei,θ y Ei,θ .
∗
A diferencia de la proposición 2.1, no existe una caracterización para Ei,θ y no necesariamente
se cumple la inclusión
∗ 0
Ei,θ ⊆ Ei,θ .
Intuitivamente, si un arco e cumple con ze (li,v (θ)) > 0 para cierto instante θ, no necesaria-
mente es porque esté llegando flujo desde el origen si , sino que podría perfectamente provenir
flujo desde otro origen. Por lo tanto, la inclusión pasa a ser
∗ ∗
Si e ∈ Ei,θ , entonces ∃j ∈ I tal que e ∈ Ej,θ
Otra gran complicación surge al definir un flujo fino. Recordando la definición 2.10 se
puede definir un flujo acumulativo de forma similar gracias al teorema 3.1, de este modo
+ −
xi,e (θ) := Fi,e (li,v (θ)) = Fi,e (li,w (θ))
y luego al derivar se obtiene
+
x0i,e (θ) := fi,e 0
(li,v (θ))li,v (θ)
Luego si se considera
X X
+
x0e (θ) := x0i,e (θ) = fi,e 0
(li,v (θ))li,v (θ) 6= fe+ (li,v (θ))li,v
0
(θ) (3.1)
i∈I i∈I
lo que conlleva a un problema al momento de hacer la deducción para obtener un resultado

similar a (1.24) pues se querría llegar a
0 0
li,w (θ) = mı́n ρe (li,v (θ), x0e (θ))
e=vw∈Ei,θ0
lo cual no es posible debido a 3.1. Esto último sugiere que un estudio mediante flujos estáticos
no es posible para el caso de múltiples orígenes.
3.3. Ejemplos.
En esta sección se presentan dos ejemplos. El primero muestra la forma en que se puede
estudiar la dinámica del flujo dentro de una red multi-od, mientras que el segundo corresponde
a un caso pequeño donde queda en evidencia una de las mayores dificultades que presentan
estas redes: La pérdida de linealidad en las ecuaciones para determinar un NTF.
39
3.3.1. Ejemplo 1.
Considere el siguiente grafo dirigido, un vector de latencias τ y uno de capacidades ν.

Para aquellos arcos sin capacidad se entenderá que poseen tasa infinita de salida, de modo
que nunca forman cola. Estos corresponden a los arcos sin una caja negra al comienzo.
u1 = 4 s1 νs t = 2
1 1
τs1 t1 = 17
t1
τs1 r1 = 1 τr2 t1 = 1
τs2 r2 = 2 τr1 r2 = 1 τr2 t2 = 1

u2 = 2 s2 r1 r2 t2
νr1 r2 = 6
νs2 t2 = 1 τr2 t3 = 1
τs3 r1 = 3
u3 = 3 s3 τs3 t3 = 15 t3
νs3 t3 = 2
τs2 t2 = 14
Figura 3.1: Ejemplo 1
Los nodos {s1 , s2 , s3 } son nodos fuentes, mientras que {t1 , t2 , t3 } son sus respectivos nodos
destino. Las tasas de flujo constante entrantes al grafo desde cada nodo respectivamente son
u1 = 4, u2 = 2, u3 = 3
En θ = 0, cada nodo origen decide enviar su flujo por el camino más corto posible, de
modo que los 3 escogen el camino interior. Con esto se tiene que


 0 θ ∈ [0, 1)
4 θ ∈ (1, 2)

fr+1 r2 (θ) =

 6 θ ∈ (2, 3)
9 θ ∈ (3, θ∗ )

donde θ∗ corresponde al instante en que algún nodo fuente decide cambiar de camino. Como
νr1 r2 = 6, a partir del segundo θ = 3 se comienza a armar cola según

0 θ ∈ [0, 3) 0 θ ∈ [0, 3]
żr1 r2 (θ) = ⇒ zr1 r2 (θ) =
3 θ ∈ [3, θ∗ ) 3θ − 9 θ ∈ [3, θ∗ ]
Teniendo esta función se puede ver cuanto tiempo le cuesta a cada partícula de flujo viajar
desde su origen a su destino por el camino central y compararlo con el camino externo:
(s1 ) Si una partícula ingresa en θ1 al grafo, sus tiempos de viaje (sin considerar lo tardado
en entrar al grafo) serán
1 3(θ1 + 1) − 9
lcentral (θ1 ) = 3 +
6
40
1
lexterno (θ1 ) = 17
Notar que si bien la partícula ingresa en θ1 al grafo, no es hasta θ = θ1 + 1 que se coloca
a la cola en el arco r1 r2 . El camino central seguirá siendo más conveniente mientras
1 1
lcentral ≤ lexterno , es decir,
θ1 ≤ 30
(s2 ) Si una partícula ingresa en θ2 al grafo, sus tiempos de viaje serán
2 3(θ2 + 2) − 9
6
2
2 2
Luego, el camino central seguirá siendo más conveniente mientras lcentral ≤ lexterno , es
decir,
θ2 ≤ 21
(s3 ) Si una partícula ingresa en θ3 al grafo, sus tiempos de viaje serán
3 3(θ3 + 3) − 9
6
3
3 3
Luego, el camino central seguirá siendo más conveniente mientras lcentral ≤ lexterno , es
decir,
θ3 ≤ 20
Si bien s2 y s3 deciden cambiar en instantes diferentes, estos cambios se verán reflejados

en r1 r2 al mismo tiempo. De este modo θ∗ = 23 corresponderá al instante en que el flujo
entrante a r1 r2 cambiará.
El nodo s1 seguirá comportándose igual mientras tanto ya que en θ1 = 22 aún le conviene

enviar flujo por el centro. Sea α ∈ [0, 2] el flujo de s2 enviado por el medio y β ∈ [0, 3]
análogamente el de s3 . De este modo,


 0 θ ∈ [0, 1)
4 θ ∈ (1, 2)



fr+1 r2 (θ) = 6 θ ∈ (2, 3)
9 θ ∈ (3, 23)




α + β + 4 θ ∈ (23, θ∗ )

Con esto la cola crece según


 0 θ ∈ [0, 3]
zr1 r2 (θ) = 3θ − 9 θ ∈ [3, 23]
(α + β − 2)(θ − 23) + 60 θ ∈ [23, θ∗ ]

Además, como los arcos s2 t2 y s3 t3 poseen tasas de salida, se comienzan a formar cola según
las funciones

0 θ ∈ [0, 21] 0 θ ∈ [0, 20]
zs2 t2 (θ) = ∗ zs3 t3 (θ) =
(1 − α)(θ − 21) θ ∈ [21, θ2 ] (1 − β)(θ3 − 20) θ ∈ [20, θ3∗ ]
Ahora, α y β deben ser tales que ir por el centro o por fuera tenga igual costo durante un
intervalo no vacío de tiempo.
41
2 2
(s2 ) Se hace un calculo similar al anterior y se igualan lcentral y lexterno
(α + β − 2)(θ2 + 2 − 23) + 60 (1 − α)(θ2 − 21)
4+ = 14 +
| {z 6 } | {z 1 }
2
lcentral 2
lexterno
Considerando que θ2 > 21, se llega a una primera ecuación 7α + β = 8.

(s3 ) De manera análoga se impone
(α + β − 2)(θ3 + 3 − 23) + 60 (1 − β)(θ3 − 20)
5+ = 15 +
| {z 6 } | {z 2 }
3
lcentra 3
lexterno
Considerando que θ3 > 20, se llega a una segunda ecuación α + 3β = 4.
Las dos ecuaciones encontradas forman un sistema cuyas soluciones corresponde a los flujos
enviados por cada nodo fuente por el camino central. Resolviendo el sistema se obtiene
α=β=1
Por lo tanto, en θ = 21, s2 comienza a enviar 1 unidad de flujo por fuera y 1 por dentro.
En θ = 20, s3 comienza a enviar 2 unidades de flujo por fuera y 1 por dentro. Con esto la
cola en r1 r2 queda constante con un largo zr1 r2 (θ) = 60, para θ ≥ 23. Los tiempos de viaje
también se congelan siendo estos 13, 14 y 15 para una partícula que ingresa a s1 , s2 o s3
respectivamente, de modo que el sistema llega a un estado estacionario.
3.3.2. Ejemplo 2.
Considere el siguiente grafo con dos orígenes {s1 , s2 } y solo un destino t. Al igual que en el
ejemplo anterior, aquellos arcos que no presentan la caja negra en el comienzo corresponden
a arcos con capacidad infinita, de modo que nunca forman cola. Las tasas de entrada al grafo
son iguales valiendo u1 = u2 = 8. El grafo corresponde a una red simétrica, de modo que las
fuentes se comportarán de la misma forma a medida que avanza el tiempo.
νs1 r1 = 4 νr1 v 1 = 2
u1 = 8 s1 r1 v1
τs1 r1 = 2 τr1 v 1 = 2
τs1 r2 = 3 τv1 t = 1
τs2 r1 = 3 τv2 t = 1
τs2 r2 = 2 τr2 v 2 = 2
u2 = 8 s2 r2 v2
νs2 r2 = 4 νr2 v 2 = 2
Figura 3.2: Ejemplo 2
42
Inicialmente cada origen decide enviar todo el flujo por los caminos exteriores, pues al no
haber colas, solo se deben revisar la suma de latencias. De este modo,
4θ θ ∈ [0, θ∗ ]

zs1 r1 (θ) =
4θ θ ∈ [0, θ∗ ]

zs2 r2 (θ) =
donde θ∗ corresponderá al instante de tiempo en que ocurra un cambio de caminos óptimos.
A partir de los 2 segundos, 4 unidades de flujo comienzan a llegar a los arcos r1 v1 y r2 v2 , los
cuales comienzan a armar cola según

0 θ ∈ [0, 2]
zr1 v1 (θ) =
2(θ − 2) θ ∈ [2, θ∗∗ ]

0 θ ∈ [0, 2]
zr2 v2 (θ) =
2(θ − 2) θ ∈ [2, θ∗∗ ]
con θ∗∗ el instante en que los cambios efectuados por las fuentes afecten a estos arcos. Con
esto ya se puede calcular el instante en que los orígenes decidirán cambiar sus caminos. Se
realiza el cálculo solo para s1 pues el del otro origen es análogo.
(s1 ) Primero se ve el tiempo necesario para recorrer el camino inferior. Si una partícula
ingresa a s1 en θ, llega a t en
Tv2 t ◦ Tr2 v2 ◦ Ts1 r2 (θ) = Tv2 t ◦ Tr2 v2 (θ + 3)

2(θ + 3 − 2)
= Tv2 t θ + 3 + +2
2
= Tv2 t (2θ + 6)
= 2θ + 6 + 1 = 7 + 2θ
Por otro lado, si la partícula ingresa en θ a s1 y utiliza el camino superior,

4θ
Tv1 t ◦ Tr1 v1 ◦ Ts1 r1 (θ) = Tv1 t ◦ Tr1 v1 θ + +2
4
= Tv1 t ◦ Tr1 v1 (2θ + 2)

2(2θ + 2 − 2)
= Tv1 t 2θ + 2 + +2
2
= Tv1 t (4 + 4θ) = 5 + 4θ
Con esto, el camino superior seguirá siendo óptimo mientras
5 + 4θ ≤ 7 + 2θ ⇒ θ ≤ 1
Por lo tanto, en θ = 1 las fuentes varían sus flujos. Sean α, β ∈ [0, 8] la fracción de flujo que
s1 envía por arriba y s2 por abajo, respectivamente. Notar que como s1 r1 y s2 r2 poseen colas,
el flujo saliente de estos arcos será 4 unidades, independiente del valor de α y β. Además no
43
es hasta el segundo 4 que el flujo cruzado llega pues las partículas respetan FIFO. De este
modo, las colas crecen según

4θ θ ∈ [0, 1]
zs1 r1 (θ) =
(α − 4)(θ − 1) + 4 θ ∈ [1, θ∗ ]

4θ θ ∈ [0, 1]
zs2 r2 (θ) =
(β − 4)(θ − 1) + 4 θ ∈ [1, θ∗ ]

 0 θ ∈ [0, 2]
zr1 v1 (θ) = 2(θ − 2) θ ∈ [2, 4]
(10 − β)(θ − 4) + 8 θ ∈ [4, θ∗∗ ]


 0 θ ∈ [0, 2]
zr2 v2 (θ) = 2(θ − 2) θ ∈ [2, 4]
(10 − α)(θ − 4) + 8 θ ∈ [4, θ∗∗ ]

Para determinar α y β se debe imponer que el tiempo tomado por cada camino sea el
mismo durante un intervalo no vacío de tiempo. Se hace el cálculo para s1 .
(s1 ) Si una partícula ingresa al grafo en θ y utiliza el camino inferior, alcanza t en
Tv2 t ◦ Tr2 v2 ◦ Ts1 r2 (θ) = Tv2 t ◦ Tr2 v2 (θ + 3)

(10 − β)(θ + 3 − 4) + 8
= Tv2 t θ + 3 + +2
2

(10 − β)(θ − 1)
= Tv2 t θ + 9 +
2
(10 − β)(θ − 1) (10 − β)(θ − 1)
=θ+9+ + 1 = θ + 10 +
2 2
Por otro lado, si la partícula escoge el camino superior, alcanza t en

(α − 4)(θ − 1) + 4
Tv1 t ◦ Tr1 v1 ◦ Ts1 r1 (θ) = Tv1 t ◦ Tr1 v1 θ + +2
4

(α − 4)(θ − 1)
= Tv1 t ◦ Tr1 v1 θ + 3 +
4
 
(α−4)(θ−1)
(α − 4)(θ − 1) (10 − β) · θ + 3 + 4
−4 +8
= Tv1 t θ + 3 + + + 2
4 2

(α−4)(θ−1)
(α − 4)(θ − 1) (10 − β) · θ − 1 + 4
= θ + 10 + +
4 2
44
Igualando los tiempos de cada camino hacia t, además de realizar simplificaciones con-
siderando que θ > 1, se obtiene

10 − β α − 4 (10 − β) α−4
= + · 1+
2 4 2 4
0 = 2(α − 4) + (10 − β)(α − 4)
0 = (α − 4)(12 − β)
Esta corresponde a la primera ecuación del sistema, cuya solución entrega la fracción
de flujo que s1 y s2 envían por cada camino. Es importantísimo notar que se pierde la
linealidad al obtener multiplicación de variables.
Gracias a la simetría que posee la red, la segunda ecuación sale simplemente de cambiar el
rol de α y β en la primera, de este modo el sistema será
(1) 0 = (α − 4)(12 − β)
(2) 0 = (β − 4)(12 − α)
En este caso, como α, β ∈ [0, 8], la solución del sistema se obtiene fácilmente y corresponde
simplemente a
α = β = 4.
Tras este cambio, se alcanza un estado estacionario en la red, pues las fuentes no tienen
incentivos a modificar sus caminos óptimos.
Los ejemplos recién estudiados entregan una referencia de cómo analizar la dinámica en
una red multi-od. La principal observación de este capítulo es que, incluso en un caso tan
sencillo con dos orígenes y un solo destino, se pierde la linealidad en el sistema de ecuaciones
correspondiente a encontrar un NTF, pues, si bien, no es explícitamente lo que uno realiza
en los cálculos, imponer esta igualdad de tiempos a lo largo de cada camino para encontrar
las tasas de flujo entrante no es nada más que pasar de trabajar con x0 , l0 a x, l.
45
Conclusión
A través de esta tesis se ha podido generalizar el trabajo realizado por R. Koch, M.

Skutella [13] y R. Cominetti, J. Correa, O. Larré [3], quienes lograron entregar los resultados
más completos hasta la fecha en modelos equilibrio dinámico para redes de colas fluidas.
Además, se completaron algunas demostraciones que habían sido omitidas en estos trabajos.
Primero se desarrolló el modelo teórico para el caso de una red con un solo origen y
múltiples destinos. Tras detallar la forma de ver una red de flujo dinámico, se caracterizaron
los equilibrios dinámicos de esta, dando varias formas de verlos, siendo el Teorema 2.1 la
principal. Utilizando la estrategia de flujos estáticos y flujos finos normalizados se logró
pasar el problema de encontrar equilibrios dinámicos a encontrar un NTF con reseteo en los
arcos activos y los arcos con cola positiva. Tras esto, se propuso el Algoritmo 1, uno de los
principales aportes a la teoría de colas fluidas, el cual permite encontrar efectivamente flujos
dinámicos en una red para el caso de tasas de entrada constantes.
Alternativamente a encontrar un NTF mediante fuerza bruta, lo cual es posible pero

costoso computacionalmente, se propusieron algoritmos basados en el trabajo de R. Koch
[12] que encuentran flujos finos normalizados mediante técnicas de cortes con congestión
máxima en grafos.
Como segunda parte en este trabajo, se estudió el caso de una red con múltiples orígenes
y destinos. Tras desarrollar brevemente su modelo teórico, se mostraron las dificultades que
surgen para extender las técnicas basadas en flujos finos normalizados. Mediante un ejemplo
pequeño de solo dos orígenes y un destino, se evidenció el mayor problema presente en la red
multi-od que corresponde a la pérdida de linealidad en las ecuaciones para determinar un
NTF, lo que dificulta ostensiblemente encontrar uno mediante fuerza bruta.
Si bien se explican las complicaciones existentes en el modelo con múltiples orígenes, aún
se debe encontrar una forma de caracterizar los equilibrios dinámicos. La razón de que el caso
con un solo origen tenga éxito es la capacidad de coordinar los tiempos locales de cada nodo
con respecto al del origen s y solo preocuparse del instante en que una partícula ingresa a la
red. Debido a la presencia de múltiples orígenes, claramente no se puede encontrar un reloj
común, sin embargo, modelos que no utilicen flujos estáticos deben ser desarrollados.
Como se pudo apreciar en los ejemplos estudiados, tras un par de iteraciones en el cambio
de caminos óptimos, los grafos siempre llegaron a estados estacionarios, donde los tiempos
de viaje, o bien se congelaban, o bien crecían de igual forma, de modo que las partículas
infinitesimales de flujo no tenían incentivos a cambiar sus rutas. Este fenómeno no es al
46
azar y resultados recientes obtenidos por R. Cominetti, J. Correa y N. Olver [4] encuentran
condiciones necesarias para que el estado estacionario se alcance. Sin embargo, sigue abierta
la pregunta de si este resultado se puede generalizar, al menos al caso de tener múltiples
destinos y un solo origen.
Si bien los dos problemas abiertos mencionados son interesantes, el más apremiante a ser
resuelto es encontrar condiciones que aseguren el término del Algoritmo 1. En la prática el
algoritmo de extensión de flujo en equilibrio dinámico termina sin problemas, sin embargo no
ha sido demostrado que los tiempos de extensión αk encontrados entre 0 y T no acumulen,
haciendo que las funciones de tasa de flujo comiencen a variar infinitamente dando origen
a oscilacions del tipo Zeno [18] y eventualmente nunca se pueda alcanzar el horizonte T . A
priori no es sencillo encontrar una forma de estudiar esto último, sin embargo herramientas de
sistemas lineales complementarios podrían ser de gran ayuda para demostrar la no existencia
de estados Zeno.
47
Bibliografía
[1] EJ Anderson and AB Philpott. Optimisation of flows in networks over time. Probability,
statistics and optimisation, pages 369–382, 1994.
[2] Narayanaswamy Balakrishnan. Advances on Methodological and Applied Aspects of Pro-

bability and Statistics. CRC Press, 2004.
[3] Roberto Cominetti, José Correa, and Omar Larré. Dynamic equilibria in fluid queueing
networks. Operations Research, 63(1):21–34, 2015.
[4] Roberto Cominetti, José Correa, and Neil Olver. Long term behavior of dynamic equi-
libria in fluid queuing networks. Sin publicar, 2017.
[5] José R Correa and Nicolás E Stier-Moses. Wardrop equilibria. Wiley encyclopedia of
operations research and management science, 2011.
[6] Reinhard Diestel. Graph theory {graduate texts in mathematics; 173}. Springer-Verlag
Berlin and Heidelberg GmbH & amp, 2000.
[7] Lisa Fleischer and Éva Tardos. Efficient continuous-time dynamic network flow algo-
rithms. Operations Research Letters, 23(3):71–80, 1998.
[8] LR Ford and DR Fulkerson. Flows in networks. 1962. Princeton U. Press, Princeton,
NJ, 1962.
[9] Lester R Ford Jr and Delbert Ray Fulkerson. Constructing maximal dynamic flows from
static flows. Operations research, 6(3):419–433, 1958.
[10] Terry L Friesz, David Bernstein, Tony E Smith, Roger L Tobin, and Byung-Wook Wie.
A variational inequality formulation of the dynamic network user equilibrium problem.
Operations Research, 41(1):179–191, 1993.
[11] David Gale et al. Transient flows in networks. The Michigan Mathematical Journal,
6(1):59–63, 1959.
[12] Ronald Koch. Routing games over time. PhD, TU Berlin., 2012.
[13] Ronald Koch and Martin Skutella. Nash equilibria and the price of anarchy for flows
over time. Theory of Computing Systems, 49(1):71–97, 2011.
48
[14] O Larré. Equilibrios dinámicos de flujo en redes. Memoria Ingeniería Civil Matemática,
Universidad de Chile, 2010.
[15] Frédéric Meunier and Nicolas Wagner. Equilibrium results for dynamic congestion ga-
mes. Transportation science, 44(4):524–536, 2010.
[16] B Ran. Ran, b. and de boyce (1996) modeling dynamic transportation networks, 1996.
[17] Jaime San Martín. Teoría de la Medida Apunte del Curso. Universidad de Chile, Facultad
de Ciencias Físicas y Matemáticas, 2013.
[18] Jinglai Shen and Jong-Shi Pang. Linear complementarity systems: Zeno states. SIAM
Journal on Control and Optimization, 44(3):1040–1066, 2005.
[19] William S Vickrey. Congestion theory and transport investment. The American Econo-
mic Review, 59(2):251–260, 1969.
[20] Daoli Zhu and Patrice Marcotte. On the existence of solutions to the dynamic user
equilibrium problem. Transportation Science, 34(4):402–414, 2000.
49

Equilibrio Dinámico en Redes de Transporte Con Múltiples Orígenes y Destinos

Cargado por

Copyright:

Formatos disponibles

Equilibrio Dinámico en Redes de Transporte Con Múltiples Orígenes y Destinos

Cargado por

Información del documento

Descripción original:

Título original

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

Equilibrio Dinámico en Redes de Transporte Con Múltiples Orígenes y Destinos

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS

EQUILIBRIO DINÁMICO EN REDES DE TRANSPORTE CON MÚLTIPLES

TESIS PARA OPTAR AL GRADO DE MAGÍSTER EN CIENCIAS DE LA

FELIPE IGNACIO GARRIDO LUCERO

Este trabajo ha sido parcialmente financiado por Centro de Modelamiento Matemático y

EQUILIBRIO DINÁMICO EN REDES DE TRANSPORTE CON MÚLTIPLES

Agradezco a Rebeca, quien me ha ayudado y entregado su apoyo desde el primer día, a

2. Modelo de cola fluida, caso de un origen y múltiples destinos. 9

3. Modelo de cola fluida, caso múltiples orígenes y destinos. 35

En términos matemáticos, los flujos en una red se describen mediante funciones

Ejemplo Considere el siguiente grafo con nodos orígenes s1 , s2 y nodo destino t:

Figura 2: Red ejemplo

En θ = 2, el tiempo de viaje por ambos caminos se iguala, y ya no resulta óptimo para s2

El estudio de flujos dinámicos fue inicialmente hecho en el área de optimización. En los

En el caso de múltiples orígenes y destinos, tras entregar el modelo de colas fluidas de

Definición 1.1 Una función f : I → R se dice absolutamente continua si y solo si,

Al conjunto de todas las funciones absolutamente continuas se denota AC(R). Se define

Definición 1.2 Se define el espacio de funciones integrables sobre R como

Este último espacio recibe el nombre de espacio de funciones localmente integrables.

Proposición 1.1 Sea f : R+ → R+ una función localmente integrable. La función F : R+ →

Definición 1.3 Un grafo dirigido consiste en un par ordenado G := (V, E) donde V es

Definición 1.4 Dado un vértice v ∈ V en un grafo dirigido, el conjunto de arcos salientes

Dado S ⊆ V un subconjunto de vértices, los conjuntos de arcos salientes y entrantes a S se

Definición 1.5 Sea φ : X → P(Y ) una multiaplicación con X, Y conjuntos no vacíos. Se

lı́m xn = x, lı́m yn = y ⇒ y ∈ φ(x)

Teorema 1.1 Teorema de punto fijo de Kakutani: Sea S ⊆ Rn un conjunto no vacío,

µ(Θ) = sup{µ(K) : K compacto, K ⊆ Θ}.

Modelo de cola fluida, caso de un origen

2.1. Modelo de red de flujo dinámico.

Dado i ∈ I, se considera ui ∈ F0 (R) la función de flujo entrante al grafo desde el nodo s

2.1.1. Flujo dinámico.

Definición 2.1 Un flujo dinámico o simplemente flujo, será un par (f + , f − ) de arreglos

Dado e ∈ E, se definen los flujos acumulativos totales como

las cuales, nuevamente, resultan ser elementos de ACloc (R).

2.1.2. Cola fluida.

Figura 2.1: Arco

Definición 2.3 Se define el largo de cola como la función ze : R → R dada por

2.1.3. Tiempo de viaje por un arco.

es decir, el primer instante q en que el flujo saliente de la cabeza de e a partir de θ + τe iguala

(a) Restricción de capacidad: fe− (θ) ≤ νe para casi todo θ ∈ R.

Fe− (Te (θ)) = Fe+ (θ). (2.12)

2.1.4. Existencia y unicidad de descomposición de un flujo de salida.

Lema 2.1 Sea un flujo de entrada descompuesto en nodos destinos

TZe (θ) TZe (θ) Zθ

Si una partícula proveniente de s ingresa a un camino P = (e1 , e2 , ..., ek ) en tiempo θ,

lP (θ) = Tek ◦ · · · ◦ Te1 (θ).

lw (θ) = mı́n lP (θ). (2.14)

Gracias a que las funciones Te son absolutamente continuas y no decrecientes lo mismo se

Para θ ∈ R se define el conjunto

Eθ0 := {e = vw ∈ E : Te (lv (θ)) ≤ lw (θ)} (2.16)

Definición 2.7 Se define el grafo de caminos θ mínimos como el grafo acíclico Gθ =

Se comienza esta sección dando la definición de flujo factible.

El siguiente lema da una caracterización de equilibrio dinámico.

Lema 2.2 Un flujo dinámico factible f es un equilibrio dinámico si y solo si ∀e = vw ∈ E,

Eθ∗ := {e = vw ∈ E : ze (lv (θ)) > 0} (2.18)

La siguiente proposición demuestra que el resultado intuitivo es cierto.

Proposición 2.1 Si f es un equilibrio dinámico entonces Eθ∗ ⊆ Eθ0 y además