MODELAJE Y SIMULACION DE DISPERSIONES DE HUMO GENERADO POR

UN INCENDIO EN UN ESPACIO CONFINADO – TUNEL; A PARTIR DE UNA

HERRAMIENTA CFD.

JULIANA EGEA MEJIA

NICOLÁS GARZÓN GARCÍA

UNIVERSIDAD DE LOS ANDES

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
BOGOTÁ D.C.
2009

1
MODELAJE Y SIMULACION DE DISPERSIONES DE HUMO GENERADO POR
UN INCENDIO EN UN ESPACIO CONFINADO – TUNEL; A PARTIR DE UNA
HERRAMIENTA CFD.

JULIANA EGEA MEJIA

NICOLÁS GARZÓN GARCÍA

Proyecto de Grado para

Optar por El Titulo de Ingenieros Químicos

Asesor
Ingeniero Felipe Muñoz Giraldo, Ph.D.

UNIVERSIDAD DE LOS ANDES

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
BOGOTÁ D.C.
2009

2
 CONTENIDOS

Pág.

LISTA DE TABLAS
LISTA DE FIGURAS
RESUMEN
INTRODUCCIÓN …………………………………………………………………. 10
1. OBJETIVOS …………………………………………………………………… 12
1.1. OBJTIVO GENERAL …………………………………………………... 12
1.2. OBJETIVOS ESPECIFICOS ...………………………………………… 12
2. MARCO TEÓRICO …………………………………………………………… 13
2.1. FENÓMENO ……………………………………………………………... 13
2.1.1. Dispersión …………………………………………………………. 13
2.1.2. Incendios …………………………………………………………... 14
2.1.2.1. Combustión …………………………………………………... 14
2.1.2.2. Tasa de liberación de calor (HRR) ………………………… 16
2.1.3. Incendios en Espacios Confinados ….………………………….. 18
2.1.3.1. Back Layering ………………………………………………... 19
2.1.3.2. Velocidad Crítica……………………………………………… 20
2.1.3.3. Sistemas de ventilación……………………………………… 20
2.1.3.3.1 Ventilación Natural…………………………………...... 21
2.1.3.3.2 Ventilación Inducida por el Tráfico…………………... 21
2.1.3.3.3 Ventilación Mecánica………………………………….. 21
2.1.3.3.3.1 Sistemas Longitudinales…………………….…… 22
2.1.3.3.3.2 Otros Sistemas………………………………….... 22
2.1.4. Criterios de Peligrosidad………………………………………….. 23
2.1.4.1. Temperatura Alcanzada.………………………………….... 24
2.1.4.2. Calor por Radiación……………………………………........ 25
2.1.4.3. Visibilidad………………………………………………......... 26
2.1.4.4. Toxicidad de Contaminantes………………………............. 27
2.1.4.5. Tiempo de Escape…………………………………….......... 27
2.1.5. Retorno de Experiencias ………………………………………… 28
2.1.5.1. Túnel de Mont Blanc, Francia-Italia (1999)………….......... 28

3
 2.1.5.2. Túnel de Tauren, Austria (1999) …………………………… 29
2.2. MODELOS Y MÉTODOS ………………………………………………. 29
2.2.1. MODELOS QUE DESCRIBEN EL FENÓMENO ……………… 31
2.2.1.1. Turbulencia ..………………………………………………….. 31
2.2.1.2. Combustión …………………………………………………… 35
2.2.1.3. Radiación ……………………………………………………… 36
2.2.2. DINÁMICA DE FLUIDOS COMPUTACIONAL (CFD) ………… 37
2.2.2.1. Ventajas de utilizar un CFD…………………………............ 38
2.2.2.2. CFD Comerciales ……………………………………………. 38
2.2.2.2.1. COMSOL Multiphisics …………………………………. 39
2.2.2.2.2. ANSYS CFX ……………………………………………. 39
2.2.2.2.3. PHOENICS…………………..………………………….. 40
2.2.2.2.4. FDS……….…………………..………………………….. 40
2.2.3. MÉTODO DE ELEMENTOS FINITOS …………………............. 40
2.2.3.1. Elementos ……………………………………………..…….. 41
2.2.3.2. Nodos ………………………………………………………… 41
2.2.3.3. Transformación Matemática ………………………............. 41
2.2.3.4. Ensamble ……………………………………………............. 43
3. METODOLOGÍA ………………………………………………………………. 44
3.1. GEOMETRIA – Solid Edge ………………………………………......... 45
3.2. ENMALLADO – Ansys CFX – Mesh ……………………………......... 46
3.3. ESPECIFICACIÓN – Ansys CFX – Pre ………………………………. 47
3.3.1. Parámetros Modelos …………………………………………...... 47
3.3.1.1. Modelo de Turbulencia …………………………………….. 47
3.3.1.2. Modelo Transferencia de Calor …………………………… 48
3.3.2. Variable Adicional …………………………………………………. 49
3.3.3. Estado Transiente …………………........................................... 50
3.3.4. Especificación de escenarios..……………………...................... 54
3.4. RESOLVEDOR – Ansys CFX – Solver ……………………………….. 55
3.5. RESULTADOS – Ansys CFX – Post …………………….................... 55
4. RESULTADOS Y DISCUCIÓN ……………………………………………… 56
4.1. CONCENTRACIÓN …………………………………………………….. 56
4.1.1. Tamaño de la fuente ……………………………………………… 56
4.1.2. Velocidad del viento ………………………………....................... 58
4.1.3. Variación del flujo másico emitido ………………………………. 60
4.2. TEMPERATURA ………………………………………………………… 61
4.2.1. Tamaño de la fuente ……………………………………………… 61

4
 4.2.2. Velocidad del viento ………………………………………….….. 62
4.2.3. Variación de la temperatura ……..………………………….….. 63
4.3. VELOCIDAD ……………………………………………………….…… 64
4.3.1. Tamaño de la fuente ……………………………………….……. 64
4.3.2. Velocidad del viento ………………………………………….….. 65
4.4. BACK LAYERING …………………………………………………….… 65
5. CONCLUSIONES ……………………………………………………….…… 67
6. RECOMENDACIONES ………………………………………………….….. 69
BIBLIOGRAFÍA 70

5
 LISTA DE TABLAS

Pág.

Tabla 1. Concentraciones necesarias para producir inconsciencia y

muerte (Migoya Valor, 2002) …………………………………………………... 27
Tabla 2. Dimensiones fuentes de emisión………………………………… 46
Tabla 3. Composición del humo……………………………………………. 49
Tabla 4. Valores encontrados de n, r (Ingason, 2009)………………………. 53
Tabla 5. Parámetros de las curvas de diseño……………………………. 54
Tabla 6. Velocidades de viento utilizadas………………………………… 55

6
 LISTA DE FIGURAS

Pág.

Figura 1. Triángulo de Fuego …………………………………………………. 15

Figura 2. Velocidad de quemado para CH4 en diferentes medios (Migoya
Valor, 2002) ………………………………………………………………………. 16
Figura 3. Esquema de Crecimiento del HRR de un incendio…………… 17
Figura 4: Tasa de liberación de calor en el tiempo …………...………… 18
Figura 5: Esquema del fenómeno de retorno de gases. (Hua, Huo a, &
Chow b, 2008)……………………………………………………………………. 20
Figura 6. Efecto de Ventilación Longitudinal ( PIARC Committee in Road
Tunnels, 1999)………………………………………………………...………… 22
Figura 7. Tiempo de exposición para producir quemaduras en la piel
(Migoya Valor, 2002) ……………………………............................................ 24
Figura 8. Flujo de calor por radiación (Migoya Valor, 2002) ……………….. 25
Figura 9. Niveles de tolerancia a diferentes flujos de calor (Migoya Valor,
2002)……………………………………………………………………………..... 26
Figura 10. Coeficiente de Extinción (Migoya Valor, 2002)…………………... 26
Figura 11. Tipos de elementos ……………………………………………... 41
Figura 12. División de un dominio en subdominios y nodos comunes……. 41
Figura 13. Esquema de un Elemento……………………………………….. 42
Figura 14. Funcionamiento de ANSYS Workbench………………………. 45
Figura 15. (a) Imagen de la entrada del túnel Número 1 de YuanJiang,
(b)Diagrama esquemático utilizado de la sección transversal (Hua, Huo a,
& Chow b, 2008)…………………………………………………….........……… 45
Figura 16: Efectos del refinamiento de la malla………………………….. 47
Figura 17: Ejemplos de los tres (3) tipos de curvas de incendios
diseñados………………………………………………………………………. 51
Figura 18: Perfil de concentración de humo, Automóvil 0,9m/s, 3 y 5min... 57
Figura 19: Perfil de concentración de humo, Van 0,9m/s, 3 y 5min………. 57
Figura 20: Perfil de concentración de humo, HGV 0,9m/s, 3 y 7min……… 57
Figura 21: Perfil de concentración de humo, HGV 0,5m/s, 3 y 25min…… 58
Figura 22: Perfil de concentración de humo, HGV 1,85m/s, 3 y 25min… 59
Figura 23: Perfil de concentración de humo, Van 0,2m/s, 3 y 15min……... 60

7
Figura 24: Perfil de concentración de humo, Van 3m/s, 3 y 15min…………. 60
Figura 25: Perfil de concentración de humo, Van 0,9m/s; -15%w y +15%w. 61
Figura 26: Perfil de Temperatura 0,9m/s: (a) Auto 10 minutos, (b) Van 15
minutos, (c) HGV 25minutos. Incendios totalmente desarrollados…………. 62
Figura 27: Perfil de Temperatura, Van 15min; (a) 0,2m/s, (b) 3m/s.
Incendios totalmente desarrollados…………………………………………….. 63
Figura 28: Perfil de concentración de humo, Van 0,9m/s; -15%T y +15%T.. 63
Figura 29: Velocidad - 0,9m/s, (a) Auto 10 minutos, (b) Van 15 minutos,
(c) HGV 25 minutos. Incendios totalmente desarrollados…………………… 64
Figura 30: Velocidad Van 15min: (a) 0,2m/s (b) 3m/s. Incendios totalmente
desarrollados……………………………………………………………………. 65

8
 RESUMEN

Incendios vehiculares en espacios confinados son uno de los principales riesgos

que enfrentan las infraestructuras viales subterraneas. Accidentes como los
registrados en Tauern (Austria) en 1999, Mont Blanc (Francia-Italia) en 1999 y San
Gotardo (Suiza) en 2001, han dejado un balance dramático de pérdida de vidas
(85 muertes en el periodo comprendido entre 1991 y 2001, en Europa), daños a
infraestructura e impactos económicos por parálisis de actividades regionales. La
ocurrencia de este tipo de eventos acarrea la múltiple aparición de fenómenos
peligrosos asociados a la emisión de gases tóxicos y al transporte de calor por
radiación, y en algunos casos extremos la propagación de ondas de sobrepresión.
El primero de estos efectos tiene una relación directa con las afectaciones en
personas. Cerca del 67% del total de victimas registrado en este tipo de sucesos
se pueden asociar directa o indirectamente a la presencia de humo o gases
tóxicos. (Pitts, 1995) Tradicionalmente el modelaje y la simulación de este tipo de
flujos de peligro se ha realizado utilizando aproximaciones semiempìricas, las
cuales hacen una simplificación de la ocurrencia del fenómeno y solo brindan un
panorama del evento. Hoy en día, la incursión de herramientas computacionales
CFD, como Ansys CFX, brindan la oportunidad de abordar fenomenológicamente
el estudio de eventos, lo que permite la realización de múltiples escenarios. Este
trabajo contempla la utilización de Ansys CFX para la evaluación de dispersión de
gases tóxicos en un túnel vehicular, como producto de un incendio automotriz. Se
hizo uso del modelo de turbulencia κ − ε , propuesto por Launder y Spalding en el
año 1974. Este utiliza la hipótesis del gradiente de difusión para relacionar los
esfuerzos de Reynolds con los gradientes de la velocidad y la viscosidad
turbulenta. Como resultado de este estudio se han simulado los perfiles de
concentración, temperatura y velocidad, dentro del túnel, en estado estacionario y
transiente y sus efectos asociados a fenomenos como el “Back Layerning”, todo
esto haciendo uso de una herramienta CFD no especializada, lo que permite el
modelaje y simulación de este tipo de eventos, con miras a la realización de
análisis de riesgos y consecuencias con una aproximación multifísica.

9
 INTRODUCCIÓN

En las actividades industriales y en general en cualquier actividad antrópica, los

incendios constituyen cerca del 41% de los accidentes reportados, aunque no
representan la mayor causa de pérdida de vidas. La combustión de materiales
inflamables, además de producir efectos por radiación, acarrea efectos de
naturaleza toxica que se proyectan sobre el ambiente y la población expuesta.
(Lees F. P., 1996)

Los incendios en espacios confinados representan uno de los mayores peligros

conocidos, dada su capacidad destructiva, su dificultad de extinción y la gran
cantidad de humo generado por la ausencia de un suministro de aire suficiente.
(Migoya Valor, 2002)

Los túneles presentan dos características principales: 1. Espacios Confinados 2.

Alta densidad de población expuesta. La ocurrencia de un incendio dado su
condición de espacio confinado puede acarrear: 1.La dispersión hacia todo el
espacio disponible por ausencia de aire de ventilación suficiente. 2. La generación
de una capa en la parte alta del cilindro que puede fluir en contra de la dirección
de la ventilación si el flujo de aire es bajo. 3. La extracción total de la capa de
gases si la ventilación es adecuada. (Modic, 2003)

En los últimos años, se han presentado numerosas pérdidas de vidas humanas y

cuantiosas pérdidas en infraestructuras por la ocurrencia de este tipo de
incidentes. Recientes desastres como los ocurridos en Baku (Azerbayan) en
1995, Mont Blanc (Francia-Italia) en 1999 y Kitzsteinhorn (Austria) en 2000, han
disparado las alarmas mundiales obligando a los ingenieros de diseño a simular
este tipo de eventos con miras a la prevención, control o mitigación de sus efectos.
Estos procesos de simulación permiten describir la dinámica de flujo del incendio y
del movimiento del humo para saber qué cantidad de oxígeno puede acceder y
alimentar el fuego, y a qué concentración de humo la gente estará expuesta.
(Neophytou & Britter, 2005)

Los software de simulación de dinámica de fluidos computacional (CFD)

constituyen una herramienta fundamental para el análisis en tres (3) dimensiones

10
de las diferentes variables involucradas en la dispersión de humo. Estos
simuladores utilizan soluciones numéricas a las ecuaciones diferenciales que
gobiernan el fenómeno, dichas soluciones generalmente requieren la manipulación
repetitiva de varios miles, incluso millones de números por lo que sus aplicaciones
se extienden a problemas de mayor detalle y sofisticación, convirtiéndose en una
herramienta práctica y útil para la ingeniería. (Anderson, 1995).

El uso de un software CFD no especializado para el modelaje y simulación de

eventos de incendio en espacios confinados, constituye una herramienta de gran
potencial dada la capacidad de diseño y generación de diferentes escenarios que
aporten información importante para la prevención y mitigación de efectos sobre
personas y estructuras.

11
 1. OBJETIVOS

1.1. Objetivo General

Modelar y simular la dispersión del humo generado por un incendio en un espacio

confinado, túnel; haciendo uso de una herramienta CFD.

1.2. Objetivos Específicos

1.2.1. Seleccionar el modelo de flujo turbulento, más adecuado, para representar

el fenómeno.

1.2.2. Seleccionar la herramienta CFD más adecuada según sus aplicaciones,

flexibilidad y disponibilidad.

1.2.3. Realizar la simulación de un caso simple para comprender el

funcionamiento de la herramienta y del modelo seleccionado.

1.2.4. Por medio del planteamiento de diferentes escenarios, estudiar la incidencia

de las variables más relevantes en la dispersión del humo dentro del túnel.

12
 2. MARCO TEÓRICO

2.1. FENÓMENO

2.1.1. Dispersión

La dispersión de una nube de gas sigue una difusión, o transporte de su masa,

producida por el movimiento molecular del aire. La nueve se diluye, a medida que
va disminuyendo su concentración por la entrada de nuevas moléculas de aire
dentro de la misma. (Casal & Montiel, 2001)

Las principales características que condicionan una situación de dispersión, están

ligadas al material emitido, a las condiciones de operación de la fuente,
condiciones meteorológicas y a la topografía del medio.

Para el material emitido se tiene:

1. La flotación: Los efectos de la flotación están relacionados con su densidad

específica, y por lo tanto, con la temperatura y el peso molecular del mismo; y de
esta forma, los fluidos pueden ser caracterizados: De flotación neutra, cuando la
densidad de la mezcla tiene un valor cercano al del aire. De flotación positiva,
generalmente gases calientes o de bajo peso molecular. Y de flotación negativa,
generalmente son gases de alto peso molecular.

Por otra parte, para la fuente, se debe considerar:

1. La cantidad de movimiento: El momentum está relacionada con la energía

cinética de la fuente, y tiene influencia sobre la densidad de la nube y su
dispersión. De acuerdo a esto, una liberación continua de material, con energía
cinética baja forma un penacho (“plume”), que tiende a comportarse como una ola;
Ahora, si la energía cinética es alta, se arrastran grandes cantidades de aire y por
lo tanto se forma un chorro (“jet”), con una forma bien definida. (Lees F. P., 1996)

2. La geometría: Los principales tipos de fuentes consideran: Puntual, lineal o de

área.

13
3. La duración: La duración de la liberación de gas de la fuente puede ser:
• Instantáneas
• Continuas
• Intermitentes

4. La elevación: Es posible clasificar la elevación de la fuente en: elevada o al nivel

del suelo, y en general, son estas últimas las relacionadas con los escapes más
peligrosos. (Lees F. P., 1996)

Y finalmente, se debe tener en cuenta la topografía (p.e. presencia de obstáculos);

y las características meteorológicas: la dirección, la velocidad y la turbulencia del
viento, además de la estabilidad atmosférica, que determina el grado de mezcla.
(La dispersión es mayor en condiciones inestables y menor en condiciones
estables).

2.1.2. Incendios

Después de presentar y comprender la importancia que tienen las características

de la fuente sobre el comportamiento del fenómeno de dispersión de gases, es
importante caracterizar la fuente de la liberación en el escenario de interés:
Incendio en espacio confinado - túnel.

2.1.2.1. Combustión (Lees F. P., 1996)

Los incendios se presentan por un proceso de combustión, que consiste

básicamente en una reacción de oxidación en la que se libera calor. Y en general,
ocurre cuando una fuente de calor entra en contacto con un material combustible,
que puede ser un líquido, un sólido o un gas. Existen tres (3) condiciones
esenciales para que se produzca un incendio, un combustible, oxígeno y calor.
Las características del medio pueden intervenir en que la combustión sea
completa o incompleta produciendo CO2 y CO respectivamente, este último es
altamente tóxico para los humanos y es una fuente considerable de riesgo.

14
 Figura 1. Triángulo de Fuego

El calor necesario para iniciar una reacción de combustión normalmente es

suministrado por una fuente externa y a su vez el calor necesario para mantenerla
lo suministra el proceso en si. La cantidad inicial de calor depende tanto de la
sustancia como de la fase en la que se encuentre, un gas o un vapor necesita
apenas de una chispa o una llama pequeña mientras que un sólido necesitará una
fuente de calor más intensa.

Las tres (3) condiciones del triangulo de combustión indican cómo se puede
atacar. El primer método es suprimir la fuente de combustible, este es usado en
casos de fugas en plantas de proceso. El segundo es remover calor, para lo cual
se usa agua sobre el fuego. El tercer método consiste en suprimir la fuente de
oxígeno lo cual se logra con el uso de espumas y gases inertes entre otros.

Como en cualquier proceso térmico en una combustión existen 3 clases de

transferencia de calor: Conducción, Convección y Radiación. Sin embargo en la
mayoría de los incendios un 75% del calor es transportado por Convección, en un
ambiente abierto este calor es transmitido a la atmósfera mientras que en un
ambiente confinado se transmite al techo del recinto. La radiación es el segundo
modo de transporte en importancia, éste es transmitido directamente a objetos
cercanos cruzando espacios abiertos sin preferencia direccional. Los aportes del
transporte conductivo no tienen tanto peso como sus contrapartes, sin embargo en
los escenarios de incendio puede existir dadas ciertas características de forma y
composición.

Los incendios han sido clasificados de acuerdo a la naturaleza de sus

combustibles y su tamaño, dicha clasificación es de especial relevancia para la
escogencia apropiada de un agente para extinguirlos. A su vez dentro de la

15
industria existe una clasificación particular en la que se tienen en cuenta otras
características.

La tasa de combustión de una mezcla inflamable depende de la velocidad de

quemado o de reacción que es una propiedad de la mezcla. Generalmente se
determina por medio de un quemador bunsen en la región laminar de la flama y se
obtiene dividiendo el flujo volumétrico del gas entre el área del cono frontal de la
llama. El valor de esta velocidad depende del tipo de material que se esté
quemando, de la fase en que se encuentre y de la cantidad de oxígeno presente
en el ambiente.

Figura 2. Velocidad de quemado para CH4 en diferentes medios (Migoya Valor, 2002)

En la figura se puede apreciar el comportamiento de la velocidad de quemado

para el Metano en Oxígeno y en Aire para diferentes relaciones estequiométricas,
allí se aprecia una composición óptima de la mezcla Oxígeno-Combustible en la
que la velocidad alcanza su máximo valor. El estudio de esta variable es
importante para análisis de la propagación de un incendio y también de la
dispersión de los productos de combustión en el ambiente.

2.1.2.2. Tasa de Liberación de Calor (HRR)

La Tasa de Calor Liberado (HRR) en un incendio, se define como la cantidad de

calor que se libera a los alrededores por unidad de tiempo y es una de las
características más importantes en el estudio de este fenómeno, ésta determina
parámetros de gran importancia como el tamaño del incendio, la temperatura y la
producción de gases tóxicos. Hoy en día el estudio y la medición del HRR se ha
16
convertido en una etapa vital de la ingeniería de seguridad contra incendios y el
desarrollo de técnicas precisas para su medición constituye una fuente de
información muy importante para la experimentación y el modelaje de estos
escenarios. (Lees F. P., 1996)

El HRR se define como la cantidad de calor que un incendio libera a sus

alrededores por unidad de tiempo. Se le considera como la fuerza motriz del
incendio, ya que de no existir tampoco existiría energía para llevar la mezcla
oxígeno-combustible a su temperatura de ignición. A su vez el HRR está
involucrado en un proceso de retroalimentación tal como se indica en el siguiente
diagrama, el calor producido en primera instancia Q1 aporta energía para crear y
calentar una mayor cantidad de mezcla inflamable que una vez llega a su
temperatura de ignición, libera un calor Q2 que cumplirá el mismo papel en un
nivel más alto (Babrauskas, 2009).

Q3
Q2
Q1

Figura 3. Esquema de Crecimiento del HRR de un incendio

Sin embargo es prudente anotar que el diagrama propone un esquema de

disponibilidad infinita de combustible lo cual no es aplicable a un caso real, los
incendios tienden a alcanzar un máximo HRR a determinado tiempo para luego
disminuir dicha tasa de liberación de calor hasta que finalmente se convierte en
cero, es decir, la combustión termina. Lo anterior indica que se pueden identificar
diferentes zonas de crecimiento, de estabilidad y de decrecimiento en la dinámica
de un incendio.

17
 Figura 4: Tasa de liberación de calor en el tiempo

• Ignición: Puede ser causada por un medio externo (chispa, llama) o de

forma espontánea debido al aumento de la temperatura de la mezcla
combustible oxígeno.
• Crecimiento: El fuego crece a un ritmo dependiente del tipo de combustible,
disponibilidad de oxígeno, configuración del compartimiento donde se
desarrolle y otros factores. El calentamiento de superficies sólidas
contiguas promueve la combustión de las mismas.
• Flashover: Es una transición rápida que existe entre el crecimiento y el
fuego totalmente desarrollado, en ella la totalidad del área superficial del
material combustible esté en llamas.
• Fuego completamente desarrollado: Durante esta etapa el ritmo de calor
liberado es máximo y el desarrollo del fuego está limitado por la
disponibilidad de oxígeno.
• Decaimiento: El combustible disminuye, por lo que la tasa de calor liberada
disminuye hasta que el incendio se apaga.

2.1.3. Incendios en Túneles

La dispersión del humo de un incendio dentro de un túnel depende de varios

factores que cambian de un incendio a otro:

• Liberación de calor
• Convección del calor
• Inclinación longitudinal
• Tipo de ventilación

18
 • Dimensiones y posibles obstrucciones
• Factores meteorológicos (dirección y magnitud del viento)

Se puede decir que el humo producido en un incendio es levantado dadas las

altas temperaturas a las que es liberado, este se mueve en la dirección en la que
la velocidad de ventilación sea alta o en ambas direcciones si la velocidad de
ventilación es baja. Esto implica que hay un pequeño espacio de algunos metros
entre la capa de humo y el suelo al menos por algunos minutos.

Conforme se libera más humo, este intercambia calor con las paredes del túnel y
con el aire de ventilación enfriándose y fluyendo hacia abajo contribuyendo a la
formación de vórtices. Esto quiere decir que en algún punto en las inmediaciones
de la fuente de emisión la capa de humo que se había elevado hacia el techo en
un principio cae hacia el suelo acabando con la zona libre de humo mencionada;
esto constituye un llenado progresivo del túnel conforme es liberado más y más
material.

2.1.3.1. Back Layering

En un incendio de un túnel de una única dirección frecuentemente existe un

sistema de ventilación longitudinal que entra en acción para crear una zona libre
de humo aguas arriba del fuego para facilitar la evacuación y extinción del mismo.
Si la velocidad de esta ventilación es lo suficientemente baja, el humo producido
puede propagarse en dirección contraria del viento, es decir, hacia la zona libre de
humo expuesta anteriormente. La velocidad crítica, es la velocidad mínima
requerida en la ventilación para eliminar el retroceso del humo hacia la zona
segura.

19
 Figura 5: Esquema del fenómeno de retorno de gases. (Hua, Huo a, & Chow b, 2008)

La figura 3 ilustra tres (3) escenarios donde la velocidad de la ventilación aumenta

progresivamente, el fenómeno es conocido como Back Layering y se ha
convertido en uno de los criterios más importantes en el diseño de sistemas de
ventilación en túneles. Su determinación tiene que ver principalmente con las
características del incendio y con la geometría del túnel.

2.1.3.2. Velocidad Crítica

En un incendio de un túnel de una única dirección frecuentemente existe un

sistema de ventilación longitudinal que entra en acción para crear una zona libre
de humo aguas arriba del fuego para facilitar la evacuación y extinción del mismo.
Si la velocidad la velocidad de esta ventilación es lo suficientemente baja, el humo
producido puede propagarse en dirección contraria del viento, es decir, hacia la
zona libre de humo expuesta anteriormente. La velocidad crítica, es la velocidad
mínima requerida en la ventilación para eliminar el retroceso del humo hacia la
zona segura.

Una velocidad de ventilación de 3 m/s se recomienda para todos los incendios

más no para un vehículo de carga pesada que transporte sustancias muy
inflamables.

2.1.3.3. Sistemas de Ventilación ( PIARC Committee in Road Tunnels, 1999)

Los objetivos de un sistema de ventilación en un túnel en caso de incendio deben

ser:
20
 • Salvar vidas haciendo posible la evacuación de los usuarios.
• Permitir el rescate y las operaciones de lucha contra el incendio.
• Mantener zonas de escape libres de humo.
• Evitar explosiones
• Limitar los daños a la estructura del túnel y a las construcciones
circundantes.

Los posibles sistemas de ventilación son:

2.1.3.3.1 Ventilación Natural

Esta se define como la ventilación dentro de un túnel que no está equipado con
ventiladores para el control mecánico del flujo de aire en la zona de tráfico. La
ventilación natural puede ser inducida por la temperatura del aire y las condiciones
meteorológicas. El túnel puede estar ventilado con aire proveniente del exterior si
el gradiente de presión entre las entradas lo permite, y a su vez por un efecto
convectivo llamado efecto chimenea.

2.1.3.3.2 Ventilación Inducida por el Tráfico

Un flujo de aire puede ser generado por el efecto pistón de los automóviles que se
mueven dentro del túnel, una gran parte del los túneles construidos a lo largo del
globo cuentan con tan solo la ventilación exterior y la producida por el flujo pistón;
no obstante estos túneles no presentan longitudes mayores a algunos cientos de
metros. Las concentraciones promedio de emisiones de los vehículos pueden
llegar a ser diluidas en estos casos, sin embargo cuando se habla de control de
humos es necesaria la implementación de sistemas mecánicos adicionales.

2.1.3.3.3 Ventilación Mecánica

Los sistemas de ventilación instalados deben proveer aire de calidad adecuada

durante la operación normal del túnel. En un escenario de incendio, estos
sistemas deben garantizar la seguridad del público allí presente así como facilitar
las labores de emergencia y extinción del mismo; así mismo un sistema de
ventilación debe prevenir la formación de una mezcla explosiva en cercanías al
fuego.

21
2.1.3.3.3.1 Sistemas Longitudinales

La ventilación longitudinal es una forma fácil y de bajo costo de ventilar túneles de

carretera, esta implica que el sistema crea un flujo de aire uniforme a lo largo del
túnel. En este tipo de ventilación mecánica el aire entra al túnel prácticamente
limpio y gradualmente arrastra las partículas de polución generadas por los
vehículos, es relativamente fácil de instalar y en túneles de un solo sentido se ve
apoyado por el efecto pistón.

Figura 6. Efecto de Ventilación Longitudinal ( PIARC Committee in Road Tunnels, 1999)

En túneles con sistemas de ventilación longitudinales la concentración de

sustancias tóxicas aumenta en la dirección del flujo de aire y disminuye con el flujo
y por lo tanto con la velocidad del mismo. Así mismo la concentración máxima
depende de la longitud del túnel por lo que el caudal de aire necesario depende
del largo del túnel.

La velocidad máxima de flujo de aire que puede fluir por el túnel sin afectar las
actividades de los vehículos y del personal de operación se encuentra entre 8 y 9
m/s, no obstante mantener este rango de velocidad de ventilación mecánica puede
llegar a ser considerablemente costoso para túneles de mayor longitud y para
túneles de doble sentido. Generalmente la ventilación longitudinal es producida
por ventiladores de flujo axial distribuidos a lo largo del túnel, el tamaño y
distribución de estos no tiene una incidencia considerable en las velocidades
promedio por lo que solo se presentan algunas reglas técnicas al respecto.

2.1.3.3.3.2 Otros Sistemas

En sistemas diferentes a los longitudinales el aire de ventilación es proveído y/o

extraido a través de ductos. Estos sistemas son clasificados según el porcentaje
de aire fresco que es proveído y el porcentaje de aire contaminado evacuado.

22
 • Sistema Totalmente Transversal: En este tipo de sistemas hay uno o más
ductos de aire fresco paralelos al tráfico, el aire exhausto es removido de la
misma forma en el lado opuesto del túnel usando uno o más ductos de
extracción. La cantidad de aire extraído por metro de la longitud del túnel,
es la misma cantidad de de aire limpio suministrado lo cual genera una
distribución constante de la concentración de contaminantes.

• Sistema Semi-Transversal: En este sistema de ventilación existen ductos

de aire fresco a través de los cuales fluye aire proveniente del exterior, no
obstante no existen ductos de extracción por lo que el aire contaminado
fluye hacia las salidas del túnel.

• Sistema Semi-Transversal Reversible: El funcionamiento de este sistema

es muy similar al Semi-Transversal, la diferencia radica en una inversión de
los papeles ya que el aire que fluye desde el exterior lo hace de forma
longitudinal mientras que el aire contaminado fluye a través de ductos de
extracción. Generalmente este modo de operación se aplica en condiciones
de incendio dentro del túnel.

• Sistema Transversal Parcial: También llamado Pseudo-Transversal es un

sistema intermedio entre los sistemas Transversal y Semi-Transversal, por
lo que adquieren características propias de cada uno de los anteriores
dependiendo la cantidad de aire proveído o extraido.

En cualquier caso los diseños aplicados por los ingenieros generalmente son
similares a los sistemas de ventilación Transversal, lo que hace necesario tener en
cuenta tanto el aire limpio como el aire contaminado. A su vez es posible la
implementación de ventiladores que garanticen el flujo longitudinal que es
necesario sobre todo en caso de incendio.

2.1.4 Criterios de Peligrosidad (Migoya Valor, 2002)

Las normas reguladoras de seguridad en túneles se basan en diferentes criterios

de peligrosidad y tienen como objetivo asegurar que los usuarios del túnel y los
rescatistas que intervienen en la extinción de un accidente que involucre fuego
puedan escapar con un margen de seguridad. Los riesgos asociados a un
escenario de incendio en un túnel son:

23
 • Temperatura alcanzada
• Calor transmitido por radiación
• Visibilidad
• Concentración de especies tóxicas
• Tiempos de escape
• Velocidad crítica

2.1.4.1. Temperatura Alcanzada

El tiempo que un hombre puede resistir bajo cierta temperatura sin que se
produzcan daños en la piel, quemaduras graves o hipertermia, fue estudiado por
Purser (1993). La siguiente gráfica muestra su estudio experimental donde se
puede resaltar la influencia de la humedad del ambiente sobre los resultados.

Figura 7. Tiempo de exposición para producir quemaduras en la piel (Migoya Valor,

2002)

Ya que la temperatura no suele ser constante en un escenario de incendio se

acude al concepto de dosis de calor por convección. La fracción de dosis de calor
incapacitadora durante 1 minuto representa un porcentaje de la cantidad de calor
que incapacitaría a esa persona durante ese lapso de tiempo.

F th =
1 (1)
exp( 5 . 1849 − 0 . 0273 ( T − 273 . 15 ))

24
Cuando la fracción acumulada excede la unidad se alcanzan condiciones
peligrosas en las que se produce incapacidad por daños y quemaduras en la piel y
en el tracto respiratorio.

2.1.4.2. Calor por Radiación

La radiación en este tipo de escenarios dependerá de la emisividad de la llama y

de los humos calientes. Para efectos de este estudio sólo se tendrá en cuenta el
aporte de los gases de combustión.

Figura 8. Flujo de calor por radiación (Migoya Valor, 2002)

La figura anterior ilustra el flujo de calor por radiación de los gases calientes de un
incendio que recibirían las personas que huyen del mismo. La línea continua
representa un ambiente abierto mientras que la línea discontinua representa un
ambiente confinado (túnel de 8 metros de ancho, 6 metros de alto y 3 de altura),
dicho flujo de calor está dado a la altura del suelo. La siguiente ecuación describe
el comportamiento del calor por radiación dada una temperatura de gases
calientes:

E f = 5.67 ⋅ 10−8 ⋅ ε t ⋅ T 4 (2)

Donde ε୲ representa la emisividad del humo.

25
 Figura 9. Niveles de tolerancia a diferentes flujos de calor (Migoya Valor, 2002)

En la figura se ilustra el nivel de tolerancia y los tiempos de exposición de una

persona ante los efectos de un flujo de calor por unidad de área. Se aprecia que el
tiempo de exposición altera de forma dramática la consecuencia sobre la salud de
la persona.

2.1.4.3. Visibilidad

Las características ópticas del humo suelen expresarse a menudo en términos de

la transmitancia, la densidad óptica y el coeficiente de atenuación o extinción.
Experimentos realizados relacionan dicho coeficiente con la visibilidad para humo
irritante y no irritante, a su vez se aprecian relaciones empíricas asociadas a cada
experimento.

Figura 10. Coeficiente de Extinción (Migoya Valor, 2002)

26
Con base en esto se han propuesto diferentes rangos y límites para que una
persona escape con éxito de un ambiente con altas concentraciones de humo, la
visibilidad se debe encontrar entre 7 y 15 m para lograr la evacuación y llevar a
término la extinción

2.1.4.4. Toxicidad de Contaminantes

Los principales productos de la combustión pueden ser clasificados en dos (2)

principales grupos: asfixiantes e irritantes. Dentro del primer grupo solo el CO y el
HCN han sido detectados en la proporción suficiente como para causar efectos
tóxicos agudos. En el caso del CO, su combinación con la hemoglobina (Hb) de la
sangre produce carboxihemoglobina (CoHb), la cual es responsable de la
disminución de oxígeno suministrado a los diferentes órganos y tejidos del cuerpo.
Este efecto es de relevante importancia en el cerebro, donde se produce
inconsciencia. En estudios recientes se ha demostrado que los efectos del HCN
son similares a los del CO, con una diferencia importante, y es el efecto conocido
como ‘Knock down’ que tiene el HCN, es decir, el nivel de incapacitación del CO
se alcanza de manera más gradual. A su vez el nivel del CoHb (cuyo efecto es
diferente en personas activas y personas sedentarias), aumenta conforme el
organismo incrementa la actividad respiratoria, por lo que cuando una persona
corre para alejarse de un incendio está contribuyendo al aumento de los niveles de
CoHb hasta el peligroso punto de la inconsciencia.

Tabla 1: Concentraciones necesarias para producir inconsciencia y muerte

(Migoya Valor, 2002)
Tiempo (minutos)
Producto 5 30
Inconsciencia Muerte Inconsciencia Muerte
CO 6000-8000 ppm 12000-16000 ppm 1400 - 1700 ppm 2500-4000 ppm
HCN 150-200 ppm 250-400 ppm 90-120 ppm 170-230 ppm
O2 10 - 13 % < 5% < 12 % 6-7%
CO2 7 -8 % < 10% 6-7% >9%

2.1.4.5. Tiempo de Escape

Diferentes estudios sobre la velocidad de escape en varios escenarios como

edificios y estaciones de metro subterráneas muestran que la velocidad promedio

27
en condiciones libres de humo varía entre 1 y 2 m/s, sin embargo bajo condiciones
de alta irritabilidad y baja visibilidad esta velocidad puede disminuir a 0.3 m/s. No
existen datos puntuales para túneles de carretera, sin embargo una buena
suposición es se encuentre entre 0.5 y 1.5 m/s dependiendo de las condiciones
del confinamiento, la iluminación y el diseño de las señales de salida.

De acuerdo a lo anterior el tiempo total que transcurre desde que se origina el

incendio hasta que aquellas personas que se encuentran en sus inmediaciones
lleguen a un lugar seguro puede variar entre 30 s y 5 min.

2.1.5. Retorno de Experiencias

A continuación se presenta un resumen de los incendios ocurridos en el túnel de

Mont Blanc, y en el túnel de Tauren, accidentes que hacen parte de los incendios
en túneles de carretera más estudiados en los últimos años.

2.1.5.1. Túnel de Mont Blanc, Francia-Italia (1999) (Rapporteur Alfred Haack,

2004)

Antes del gran incendio ocurrido el 24 de Marzo de 1999, ya habían ocurrido un

total de 17 incendios de camión desde su inauguración en 1965, la mayoría de
estos fueron controlados con los extinguidores propios de cada camión, mientras
que en al menos 5 fue necesaria la intervención de la brigada contra incendios.

El incendio ocurrido en 1999 empezó en un camión articulado que transportaba 9

toneladas de margarina y 12 toneladas de harina, aparte de la carga otros
materiales inflamables en el vehículo eran aproximadamente 550 Lt de
combustible y la espuma altamente inflamable usada en el aislamiento térmico del
vagón de refrigeración. Este vehículo se detuvo a 6.5 km de la entrada Francesa,
un total de 14 vehículos de carga pesada se detuvieron detrás del vehículo
incendiado.

La velocidad de la ventilación en el sitio donde el vehículo se detuvo se

encontraba en el rango de 1 a 1.5 m/s, el sistema consistía en una ventilación
transversa con un nivel de extracción de aire mayor al de suministro. El HRR
máximo en los alrededores del primer camión se estimó en el rango de 75 a 110
MW. Durante el incendio la carga de margarina se derritió convirtiéndose en un
aceite peligrosamente inflamable, contribuyendo así a la propagación del incendio
hacia los alrededores del vehículo; a su vez la radiación producida por las llamas
28
que se desviaban con el techo en la dirección del viento, contribuyó a que los
vehículos contiguos llegaran a su ignición.

Como resultado de las 53 horas que duró el incendio, 39 personas fallecieron; 29

fueron encontradas dentro de los vehículos, 9 personas por fuera y un bombero
murió por heridas que adquirió dentro del túnel.

2.1.5.2. Túnel de Tauren, Austria (1999) (Rapporteur Alfred Haack, 2004)

Este accidente ocurrió el 29 de Mayo de 1999 y fue originado por un choque entre
un camión, 4 vehículos y un segundo camión que se encontraba más adelante, el
incendio involucró un total de 16 camiones y 24 carros.

La ventilación del túnel de Tauren consiste en un sistema totalmente transversal

con 4 secciones. El máximo volumen de aire fresco es de aproximadamente 190
m3/s km mientras que el de aire exhausto es de 114 m3/s km. Una vez iniciado el
incendio el sistema de extracción de la sección involucrada fue operado a un
caudal de 230 m3/s km durante aproximadamente 15 minutos ya que si bien el
humo era extraído sin mayor problema, la temperatura del mismo no permitió que
el sistema siguiera operando, promoviendo un flujo de humo hacia el portal norte
del túnel. El HRR máximo del incendió se ubica en el rango de 300 a 400 MW.

Se concluyó que el sistema de ventilación allí instalado era completamente

inadecuado para lidiar con un incendio de estas características, lo cual produjo un
desafortunado balance de 12 personas muertas y 57 heridas.

2.2. MODELOS Y MÉTODOS

El comportamiento de las masas de humo emitidas por un incendio dentro de un

túnel ha sido estudiado por diferentes autores, cuyos intereses han abarcado tanto
la parte experimental como la de modelación y simulación. Sin embargo, la
dificultad asociada a la experimentación a escala real, en términos de seguridad,
esfuerzo requerido y altos costos de ejecución, ha llevado a que el modelaje vaya
tomando más fuerza con el pasar de los años.

No obstante, también se ha encontrado que el modelaje matemático del

crecimiento del fuego y la dinámica del humo no es fácil debido a los procesos
físicos involucrados como la turbulencia, la flotación, la convección, la
transferencia de calor por radiación, y la combustión; además de otras variables
29
como la ubicación del fuego, las condiciones externas de viento, y los sistemas de
ventilación disponibles dentro del túnel, que también afectan los resultados. Esto
ha llevado al desarrollo de dos tipos de modelos, conocidos como modelos
zonales y de campo.

Los modelos zonales utilizan expresiones empíricas, para describir el fenómeno;

son modelos unidimensionales que dividen cada espacio en un pequeño número
de volúmenes (zonas), como la capa superior-caliente, la capa inferior-fría, el
penacho de fuego y la frontera del compartimento. Cada una de estas zonas se
asume de propiedades uniformes. Las ecuaciones de conservación se solucionan
numéricamente junto con las ecuaciones empíricas obtenidas previamente de
experimentos. Algunos ejemplos de estos modelos son los desarrollados por el
NIST1: HARVARD6 code, FAST, CCFM, y CFAST2.

Los modelos de campo, por otra parte, utilizan la aproximación por los “principios
básicos” en la solución de las ecuaciones clásicas de conservación de masa,
momento y energía en puntos discretos en tiempo y espacio. Con esta
aproximación, se puede entonces modelar la naturaleza multidimensional y
dependiente del tiempo, de un incendio, en túneles obstruidos o libres y de
cualquier geometría; por lo cual, solo se hace posible con la llegada de
computadores modernos con gran capacidad y alta velocidad de procesamiento.
Del gran número de modelos de campo que se han desarrollados, algunos son
códigos de propósito general como CFX (Módulo CDF3 de ANSYS, Inc.),
Phoenics, Flow3D, los cuales son aplicables a una gran variedad de problemas;
mientras que otros códigos como FDS (desarrollado por el NIST), FIRE, JASMINE,
desarrollado por la FRS4, han sido creados especialmente para problemas
relacionados con incendios.

Las herramientas CFD, son una forma rigurosa de estudiar el comportamiento del
humo en estos escenarios (incendio dentro de un túnel). Actualmente, se deben
incluir ciertas suposiciones en los modelos de campo, debido a que el
entendimiento los procesos físicos involucrados, es límitado, y por esto, la
precisión de los resultados obtenidos con la herramienta CFD depende de la
pertinencia y la exactitud de los modelos físicos empleados en los códigos. Así,
1
NIST, por sus siglas en Inglés. National Institute of Standards and Technology.
2
CFAST, Por sus siglas en Inglés. The Consolidated Model of Fire and Smoke Transport.
3
CFD, por sus siglas en Inglés. Dinámica de fluidos computaional.
4
FRS, por sus siglas en Inglés. Fire Research Station
30
por ejemplo, para realizar esta simulación se requieren los submodelos de
combustión, turbulencia y transferencia de calor por radiación, para poder describir
de forma completa el fenómeno de un incendio en un espacio confinado. (Shorab,
Shashi, & Surendra, 2008).

A continuación se presenta una breve introducción de los modelos de flujo

turbulento y los métodos numéricos de solución de ecuaciones diferenciales
empleados en el desarrollo del proyecto.

2.2.1. MODELOS QUE DESCRIBEN EL FENÓMENO

El grupo de ecuaciones, que describen los procesos de transferencia de masa,

calor y momento son ecuaciones diferenciales parciales, que no tienen una
solución general de forma analítica y es por esto que deben ser solucionadas
numéricamente; por otra parte, las ecuaciones generales (conocidas como las
ecuaciones de Navier-Stokes), son modificadas/ simplificadas de acuerdo a las
condiciones del fenómeno estudiado.

Teniendo en cuenta esto, y la complejidad de los procesos involucrados en un

incendio que se presenta dentro de un túnel, además de la importancia que tienen
los modelos utilizados en la simulación del fenómeno, y por lo tanto en los
resultados obtenidos con la misma; se va a puntualizar en los modelos
considerados más relevantes: turbulencia, combustión y radiación.

2.2.1.1. Turbulencia

La turbulencia consiste en fluctuaciones en el campo de flujo en tiempo y espacio,

es un proceso complejo, principalmente porque estan involucradas las tres (3)
dimensiones, es inestable y contiene un amplio rango de escalas de tiempo y
longitud. La turbulencia ocurre cuando las fuerzas inerciales en el fluido se hacen
significativas comparadas a las fuerzas viscosas, y se caracteriza por un alto
número de Reynolds; por lo tanto puede tener un efecto significativo sobre las
características del flujo.

Para incendios en espacios confinados, las escalas dependerán de las

características del incendio y de las dimensiones del recinto, pero usualmente las
las escalas de longitud puede variar desde menos de un milímetro hasta unos

31
metros, y es por esta razón que el modelaje de flujos turbulentos se convierte en
un reto.

Para esto, se han propuesto varias aproximaciones, las cuales pueden ser
clasificadas básicamente en tres (3) categorias:

• Simulación Numérica Directa (DNS5): Todos los movimientos turbulentos se

resuleven solucionando directamente las ecuaciones de Navier-Stokes que
gobiernan los flujos de fluidos, sin incluir ecuaciones adicionales. Esta
aproximación requiere un número muy grande de elementos y por esto es
poco práctico para usos industriales. No se ha utilizado en la Ingeniería de
seguridad contra incendio.

• Simulación de Remolinos Grande (LES6): Todos, exceptuando los

movimientos turbulentos más pequeños son resueltos con las ecuaciones
de Navier-Stokes. Los remolinos más finos se pueden ignorar o también
pueden ser modelados. Sin embargo, y aunque esta aproximación permita
emplear menos elementos que la descrita anteriormente, todavía es poco
práctica para ser aplicada a incendios y al movimiento de humo en espacios
complejos. Solo hasta hace poco se viene aplicando a incendios, pero su
uso se ha limitado a escenarios simples.

• Navier-Stokes con Reynolds promediado (RANS7): Las ecuaciones de

Navier-Stokes son modificadas con la introducción de cantidades
fluctuantes promediadas en el tiempo y asi las ecuaciones obtenidas no
apuntan a resolver los movimientos turbulentos, pero si proporcionan las
características más relevantes del flujo. Teniendo en cuenta esto, solo se
requiere que la malla sea suficientemente fina para captar las
carasteristicas importantes del flujo; convirtiendose en la aproximación más
práctica de las tres, y por lo tanto la más utilizada en la ingeniería de
seguridad contra incendios, al igual que en otras aplicaciones de la
ingeniería. (Gobeau, Ledin, & Lea, 2002)

5
DNS, por sus siglas en Inglés, Direct Numerical simulation
6
LES, por sus siglas en Inglés, Large Eddy Simulation
7
RANS, por sus siglas en Inglés, Reynolds Averaged Navier-Stokes
32
Los modelos de turbulencia basados en las ecuaciones RANS son conocidos
como Modelos de Turbulencia Estadísticos. Entre estos, se ha encontrado
adecuado utilizar el modelo de turbulencia κ − ε para el modelaje y simulación de
incendios.
En este modelo, tanto la velocidad como la escala de longitud son solucionadas
con ecuaciones separadas de transporte. Estos modelos de turbulencia de dos
ecuaciones son extensamente utilizados, ya que con estos se tiene un buen
equilibrio entre el esfuerzo computacional y la precisión de los resultados.

El modelo κ − ε utiliza la hipótesis del gradiente de difusión para relacionar los

esfuerzos de Reynolds con los gradientes de la velocidad y la viscosidad
turbulenta. La viscosidad turbulenta es modelada como el producto de una
velocidad turbulenta y una escala de longitud turbulenta. La escala de velocidad
de turbulencia es calculada de la energía cinética turbulenta, que proporciona la
solución de su respectiva ecuación de transporte.

La κ representa la energía cinética turbulenta y está definida como la varianza de

las fluctuaciones en la velocidad y ε es la disipación turbulenta de los remolinos (la
tasa a la que las fluctuaciones de la velocidad se disipan).

El nuevo sistema de ecuaciones, con las dos nuevas variables que introduce el
modelo, es entonces:

Ecuación de Continuidad

∂ρ
+ ∇ • ( ρU ) = 0 (3)
∂t

Ecuación de Momento
∂ρU
+ ∇ • ( ρU ⊗ U ) − ∇ • ( µ eff ∇U ) = ∇p'+∇ • ( µ eff ∇U ) T + B (4)
∂t

Donde B es la suma de las fuerzas, µ eff es la viscosidad efectiva que aporta a la

turbulencia, y p ' es la presión modificada dada por:

33
 2
p' = p + ρκ (5)
3

El modelo κ − ε , está basado en el concepto de viscosidad de los remolinos y por

esto, la definición de µ eff , en la Ecuación 6, en la que aparece el término de la

viscosidad turbulenta ( µt ), la cual está relacionada con la energía cinética

turbulenta y la disipación por la Ecuación 7 ( C µ
es una constante).

µ eff = µ + µ t (6)

κ2
µt = Cµ ρ (7)
ε

Las ecuaciones diferenciales de transporte, para la energía cinética turbulenta y la

tasa de disipación de la turbulencia, de las cuales se obtienen los valores de κ y ε
son:

∂(ρκ)  µ  
+ ∇ • (ρUκ ) = ∇ •  µ + t ∇κ  + pk − ρε (8)
∂t  σκ  

∂(ρε )  µ   ε
+ ∇ • ( ρUε ) = ∇ •  µ + t ∇ε  + (Cε 1 pk − Cε 2 ρε ) (9)
∂t  σε   κ

Donde Cε1 , Cε 2 ,σ k ,σ ε son constantes y pk es la producción de turbulencia debido

a las fuerzas viscosas y de flotación, que se modela utilizando la siguiente
expresión. (ANSYSInc)

2
pκ = µt + ∇U • (∇U + ∇U T ) − ∇ • U (3µt ∇ • U + ρκ ) + pkb (10)
3

34
2.2.1.2. Combustión

Los productos de la combustión son muchos y varían. Generalmente, la mezcla

resultante de productos es conocida como humo; el humo producido en un
incendio, es una combinación volátil de partículas sólidas de hollín, combustible no
quemado, o parcialmente quemado, y una mezcla de aire caliente y gases tóxicos.
(Rapporteur Alfred Haack, 2004).

Como no es fácil conocer con exactitud la cinética de la reacción y la velocidad de

la combustión en incendios en espacios confinados, (Xue, Ho, & Cheng, 2001) no
es sorprendente ver que la selección de modelos de combustión entre
investigadores varíe dependiendo de su aplicación; sin embargo es importante
mencionar que hay un pequeño efecto en la predicción cualitativa de los patrones
de flujo de aire en los incendios en túneles, dependiendo del modelo utilizado.

Se tienen básicamente dos formas para representar la combustión: primero un

modelo de fuente de calor volumétrica (VHS8), que no predice la liberación de
calor y de humo en la llama, pero si su transporte a partir de ella. Las cantidades
de calor y de humo liberadas por el incendio, y el volumen de llama donde ocurren
las liberaciones tienen que ser especificadas por el usuario. La distribución de
calor y de humo se asumen uniformes en el volumen de la llama; otra
aproximación, son los modelos de combustión que se enfocan en predecir, aunque
de una manera simple, las reacciones químicas que pasan en la llama, se tiene en
cuenta la distribución no uniforme de calor en la región de la llama y la influencia
del flujo local.

El acercamiento más simple para representar un incendio generalmente es un

modelo VHS, y por lo tanto es ampliamente utilizado para investigar las
consecuencias de un incendio en grandes espacios. Esta aproximación no tiene
en cuenta la reacción química, y por lo tanto, no predice el proceso de combustión,
pero en lugar de esto representa el efecto resultante - la producción de calor y
humo, imponiendo una distribución típica uniforme de estos, sobre un volumen
preestablecido. La única condición necesaria para que el modelo VHS sea válido
consiste en que el área del incendio debe ser restringida. (Gobeau, Ledin, & Lea,
2002).

8
VHS, por sus siglas en Inglés, Volumetric Heat Sourse
35
Como se mencionó anteriormente, en el modelo VHS, el HRR, y las características
del volumen (tamaño y forma), equivalentes al incendio, deben ser especificadas
como datos de entrada y para esto, se utilizan datos experimentales o
correlaciones previamente establecidas.

2.2.1.3. Radiación

Muchos flujos de interés práctico, ocurren en situaciones donde el fluido y/o las
fronteras de inclusión están expuestos a altas temperaturas. En tales situaciones,
el efecto de transferencia de calor por radiación es significativo. (ANSYSInc)

La radiación puede jugar un papel muy importante en la transferencia total de calor

en flujos de combustión, generalmente cuando las temperaturas están por encima
de los 600K. Las principales fuentes de radiación son el CO2 y el H2O, que emiten
la energía en bandas discretas, y hollín, que a su vez emiten radiación en todas
las longitudes de onda.

La transferencia de calor por radiación ocurre entre los emisores y receptores, p.e.
entre superficies sólidas y las mezclas de fase hollín/gas de las llamas.

Actualmente se cuenta con varias aproximaciones para modelar la radiación.

Salvo declaración contraria, estos modelos proporcionan una predicción de la
transferencia de calor por radiación, que para los modelos más sofisticados,
significará la solución de un grupo de ecuaciones diferenciales. El valor obtenido
es entonces adicionado como un término de la fuente en la ecuación de transporte
para la temperatura, puntualmente la ecuación de energía y de entalpia. (Gobeau,
Ledin, & Lea, 2002)

En el las simulaciones realizadas se utiliza el modelo P1 o la aproximación

diferencial. Este modelo asume que el medio participante es opaco y emplea la
aproximación P1 armónico esférico para la solución de la ecuación de la
transferencia de radiación (RTE9). (Shorab, Shashi, & Surendra, 2008)

El modelo P1 es también una simplificación de la ecuación de transporte de

radiación, en la que se asume que la intensidad de la radiación es isotrópica, o
independiente de la dirección en una posición dada en el espacio. En este, las

9
RTE, Por sus siglas en ingles, Radiative Transfer Equation
36
fronteras abiertas (entradas, salidas y aperturas), son tratadas como fronteras
totalmente transparentes, es decir, ellos absorben toda la energía saliente, y la
energía entrante es calculada como un cuerpo negro en la temperatura local o en
una temperatura especificada por el usuario. De manera general se tiene entonces
que el flujo de calor radioactivo espectral en el límite de difusión para un emisor,
receptor, y un medio de dispersión lineal, puede ser calculado con la Ecuación 11;
y finalmente en la Ecuación 12 se tiene la ecuación para la radiación incidente
espectral que resulta de la sustitución de los términos mencionados en la ecuación
de transporte de radiación.

1
q rv = − ∇Gv
3( K av − K sv ) − AK sv (11)
 1 
− ∇ ⋅  ∇Gv  = Kav ( Ebv − Gv ) (12)
 3( Kav − K sv ) − AKsv 

donde A es el coeficiente de anisotropía lineal. (ANSYSInc).

2.2.2. DINÁMICA DE FLUIDOS COMPUTACIONAL (CFD)

La dinámica de Fluidos computacional es una de las ramas de la mecánica de

fluidos que utiliza métodos numéricos y algoritmos para solucionar y analizar los
problemas que involucran flujos de fluidos (gases y líquidos), la transferencia de
calor, y otros procesos físicos relacionados. (ANSYSInc)

Las características y el comportamiento de cualquier fluido están gobernadas por

tres (3) principios fundamentales: 1. La Conservación de Masa; 2. La
Conservación de Energía y 3. La segunda ley de Newton (
Fuerza = Masa × Aceleració n ). (Anderson, 1995) Las leyes de conservación
generalmente se expresan en términos de los cambios de unas variables con
respecto a otras (ecuaciones diferenciales parciales - EDP). Un CFD se encarga
de encontrar una solución numérica a dichas EDP mediante métodos de algebra
discreta de manera tal que provee una predicción de la características de un
fenómeno en un campo puntual de espacio y tiempo.

37
2.2.2.1. Ventajas de utilizar CFD

Utilizar un CFD, para predecir la dinámica del humo, tiene un gran número de
ventajas comparadas con otras técnicas, más sencillas, que han sido
desarrolladas en el pasado.

Los modelos empíricos, usualmente están limitados a escenarios muy simples, y

los resultados obtenidos a partir de los mismos, no pueden ser extrapolados más
allá de la configuración en la cual los experimentos fueron realizados. Por otra
parte, un CFD permite predecir el movimiento del humo, sin las restricciones
impuestas por otras ténicas de simulación más simples, en particular puede ser
usado para representar diferentes configuraciones de incendio en geometrías
complejas.

Un CFD es una herramienta predictiva que puede ser utilizada tanto en la etapa de
diseño, como en la evaluación de una estructura ya terminada y por esto es
potencialmente útil para evaluar los efectos de cambios en el diseño estructural y
de la ventilación, y evaluar el funcionamiento de medidas de seguridad tomadas
sobre diferentes incendios; diferenciando tamaño, duración y ubicación, entre
otras caracteristicas relevantes en la dinámica del mismo. Adicionalmente, los
resultados obtenidos con un CFD, contienen gran cantidad de información sobre
las variables asociadas a los flujos estudiados; que puede incluir la velocidad, la
temperatura, el humo y concentraciones de gas en decenas o aún cientos de miles
de puntos dentro del espacio.

No obstante, los CFDs también tienen limitaciones, es una técnica de gran

potencial y complejidad; que requiere de conocimientos previos del fenómeno para
poder aplicarse y producir resultados confiables. Además es importante recordar
que cualquier ténica de modelaje, no pasa de ser una aproximación de la realidad
y por lo tanto un uso inadecuado de la misma, puede producir predicciones que se
alejen significativamente de la situación real. (Gobeau, Ledin, & Lea, 2002)

2.2.2.2. CFD Comerciales

Como se mencionó anteriormente, y debido a que cada día los software de

Dinámica de Fluidos Computacional toman más fuerza, se ha desarrollado una
gran variadad de CFD; algunos de proposito general, los cuales pueden ser

38
aplicados a diferentes situaciones relacionadas con el flujo de fluidos, la
transferencia de calor entre otros procesos físicos; y otros de propósito especifico
aplicables unicamente a escenarios de incendio.

2.2.2.2.1. COMSOL Multiphysics

Anteriormente conocido como FEMLAB este es un paquete de software de análisis

y resolución por elementos finitos aplicable a múltiples ramas de la física y la
ingeniería, especialmente fenómenos acoplados, o multifísicos. COMSOL
Multiphysics también ofrece una amplia y bien gestionada interfaz a MATLAB y
sus toolboxes que proporcionan una amplia variedad de posibilidades de
programación, preprocesado y postprocesado. Los modelos de resolución que
emplea pueden interactuar con un ambiente de programación numérico llamado
Comsol Script.

2.2.2.2.2. ANSYS CFX

Dentro del paquete tecnológico de ANSYS Inc., se encuentra, ANSYS CFX, un

software CFD, de propósito general que ha sido aplicado con un alto rendimiento
durante más de 20 años. La clave que tiene ANSYS CFX para alcanzar soluciones
confiables y exáctas, de una manera rápida y robusta, es la avanzada tecnología
que utiliza su moderno resolvedor, éste es a su vez la base para la escogencia de
gran cantidad de modelos físicos con los cuales es posible abordar prácticamente
cualquier tipo de fenómeno relacionado con el flujo de fluidos: laminar a turbulento
(incluyendo la transición), incompresible a totalmente compresible, subsónico a
trans y super sónico, isotérmico o con transferencia de calor por convección y/o
radiación, no reaccionante a combustible, un solo fluido o una mezcla de fluidos en
una o varias fases (incluyendo superficies libres), entre muchas otras
configuraciones. Todo ésto incluido en un moderno, intuitivo, y flexible ambiente
para el usuario, con amplias posibilidades de personalización, escritura, y un
poderoso lenguaje de expresión.

Pero ANSYS CFX es más que solo un código CFD: con su integración en la
plataforma de ANSYS Workbench, se puede obtener el beneficio de tener
conexiones bidireccionales a todos los módulos (geometría, Mallado,
Especificación, Resolverdor, Visualización de resultados), lo cual facilita, y hace
más práctico su uso. (ANSYSInc)

39
2.2.2.2.3. PHOENICS

Es un software CFD y el producto principal de CHAM10, es una aplicación de

propósito general pensada para la resolución de problemas complejos en el
ámbito de la dinámica de fluidos. PHOENICS proporciona un entorno preciso y
robusto para la resolución y simulación de procesos que incluyan el tratamiento de
flujos de fluidos, calor o transferencia de masa, reacciones químicas y/o
fenómenos de combustión en aplicaciones de ingeniería y para el tratamiento de
problemáticas medioambientales. (Aertia Software, 2005)

2.2.2.2.4. FDS

El Simulador de la Dinámica de Incendios (FDS11), es un software CFD,

desarrollado por el Instituto Nacional de Estadares y Tecnología (NIST),
especificamente para el flujo de fluidos conducido por incendios. El FDS soluciona
numéricamente una forma de las ecuaciones de Navier-Stokes apropiada para un
flujo lento, termicamente conducido, haciendo énfasis en el transporte de humo y
de calor a partir del incendio.
Smokeview es un programa de visualización separado que se utiliza para mostrar
los resultados de una simulación obtenida con el FDS. (NIST, 2009)

2.2.3. METODO DE ELEMENTOS FINITOS

Este método usa formulaciones integrales para crear un sistema de ecuaciones

algebraicas a partir de las ecuaciones diferenciales ordinarias (ODE’s)
provenientes del modelo κ − ε . Los pasos básicos para la implementación de este
método se plantean a continuación (Moaveni, 1999):

1. Discretizar el dominio en elementos finitos.

2. Asumir una función continua como solución.
3. Desarrollar las ecuaciones para el elemento.
4. Ensamblar los elementos y construir la “Stiffness Matrix”.
5. Aplicar condiciones de frontera y valores iniciales
6. Resolver las ecuaciones lineales o no lineales de forma simultánea.

10
CHAM, por sus siglas en Inglés, Concentration Heat and Momentum, Ltd. Es una compañía de
consultoría y desarrollo de aplicaciones verticales especializadas en la simulación de procesos de
flujo de fluidos y transferencia de calor.
11
FDS, por sus siglas en Inglés, Fire Dynamics Simulator
40
2.2.3.1. Elementos

Dependiendo de la geometría y de las características del fenómeno los elementos

pueden ser de línea, de área o de volumen, estos están identificados por una
secuencia global y específica de números de nodos.

Figura 11. Tipos de elementos

2.2.3.2. Nodos

Los elementos están conectados mediante nodos. Los nodos están definidos
como el lugar en el espacio donde los grados de libertad del problema existen,
existen a
su vez los grados de libertad de los nodos dependen de la naturaleza física del
fenómeno y del tipo de elemento; por ejemplo, para un problema de flujo el grado
de libertad de los nodos es la velocidad, mientras que para un problema de
conducción de calor
lor es la temperatura. (Madenci & Guven, 2006).
2006)

Figura 12. División de un dominio en subdominios y nodos comunes

2.2.3.3. Transformación Matemática

A continuación se hará una descripción simple y muy general del funcionamiento

funcionami
de la transformación que hace en el proceso de determinación de las ecuaciones
de cada elemento:

41
 Figura 13. Esquema de un Elemento

En la figura 13 se aprecian 2 nodos y un elemento de dimensión Y.

Suponiendo que allí ocurre un fenómeno físico regido por la siguiente ecuación:

dr
S =C (13)
dy

Donde C es una constante, la forma discretizada de esta ecuación sobre el

elemento quedaría de la forma:

∆r
S = C (14)
Y

Donde ∆r = r2 − r1 y S 1 = − S 2 , por lo tanto la Ecuación 14 se puede reescribir para

los nodos 1 y 2.

( r2 − r1 )
S1 = C , S 2 = C ( r1 − r 2 ) (15)
Y Y

De forma matricial ambas ecuaciones para este elemento se plantean de la

siguiente forma:

 − 1 1  r1   S1 
K    =   (16)
 1 − 1  r2   S 2 

La Ecuación 16 representa el sistema de ecuaciones algebraicas a resolver para

hallar los valores de r en cada nodo y por lo tanto los correspondientes valores de
s. De forma general esta ecuación se expresa como:

42
k (e) r (e) = s(e) (17)

(e )
Donde k es la matriz característica del elemento, θ (e ) es el vector de
(e)
desconocidos nodales y q es el vector de del lado derecho de la igualdad.
(Madenci & Guven, 2006)

2.2.3.4. Ensamble

Para resolver el problema sobre todo el dominio es necesario realizar un ensamble

de las ecuaciones de cada elemento para obtener un sistema global de
ecuaciones:

Ku = F (18)

Donde K es el ensamblaje de las matrices características de los elementos, u es

el vector de desconocidos nodales y F es el ensamblaje de las matrices del lado
derecho de la igualdad de cada elemento. Dicho ensamblaje consiste en la
sumatoria sobre todo el dominio de todas las matrices involucradas.

De manera que se tiene un conjunto de ecuaciones algebraicas que han de

resolverse numéricamente para encontrar la solución del problema y de esa
manera lograr una caracterización del fenómeno. (Madenci & Guven, 2006)

43
 3. METODOLOGÍA

Haciendo uso de la herramienta CFD, ANSYS CFX, se diseñaron diferentes

escenarios de incendio en un túnel con el objetivo de apreciar el comportamiento
de las variables más importantes involucradas en el desarrollo de un incendio en
un espacio confinado. Como se ha expuesto anteriormente las variables críticas
que determinan el desarrollo y los efectos de un incendio en una geometría
confinada tipo túnel son: la geometría de la fuente de emisión, la cantidad de
materia liberada, la cantidad de energía térmica liberada y las condiciones de
ventilación; por lo tanto se crearon diferentes combinaciones de dichos parámetros
para establecer la sensibilidad del fenómeno ante su variación.

Se hizo un análisis en estado transitorio con el objetivo de caracterizar el

comportamiento de variables potencialmente peligrosas para personas y
estructuras, como lo son la temperatura y la concentración de humo. El
comportamiento de dichas variables en el tiempo arroja información mucho más
relevante, dado que en un escenario de incendio en espacios confinados, la
evolución del evento es de vital importancia para la toma de decisiones en pro de
la prevención y mitigación de efectos.

ANSYS CFX, fue escogido por su gran trayectoria, versatilidad y disponibilidad;

éste CFD permite llevar a cabo dicho análisis de una forma adecuada, en la que
las diferentes etapas del proceso de simulación se presentan ordenadamente en
una interfaz amable y clara. Funciona mediante una plataforma de proyecto
llamada ANSYS Workbench en la que se divide el proceso en 5 etapas principales
que permiten que el proceso sea más funcional y asi se pueda llevar a término una
simulación de forma exitosa.

44
 1
Solid Edge

3 2

CFX ANSYS
Mesh
PRE WORKBENCH

4 5
CFX
Solver
POST

Figura 14. Funcionamiento de ANSYS Workbench

3.1. GEOMETRIA – Solid Edge

Por medio del software Solid Edge se generó la geometría en 3-D de un túnel de
carretera. Las dimensiones utilizadas corresponden al túnel número 1 de
Yuangjiang ubicado al oeste de China, en el cual se desarrollaron diferentes
experimentos de incendio a escala real, el túnel tiene un área transversal
correspondiente a medio circulo, tiene una longitud de 1032 m, 7.2 m y 10.8 m de
alto y ancho respectivamente. (Hua, Huo a, & Chow b, 2008). Es importante
mencionar que por efectos prácticos no se utilizó la longitud total del túnel, y en
lugar de esto, se trabajó sobre la tercera parte correspondiente a la zona central
del mismo.

(a) (b)

Figura 15. (a) Imagen de la entrada del túnel Número 1 de YuanJiang, (b)Diagrama
esquemático utilizado de la sección transversal (Hua, Huo a, & Chow b, 2008)

45
En el centro del túnel se ubicaron las tres (3) diferentes fuentes de emisión de
humo, las cuales corresponden a un automóvil de 4 pasajeros, una van y a un
vehículo de carga pesada (HGV12). Las dimensiones se presentan a continuación.

Tabla 2: Dimensiones fuentes de emisión

Alto (m) Ancho (m) Largo (m)
Automóvil 1,48 1,88 5,15
Van 2,08 2,02 5,69
HGV 4,25 2,50 12,60

De esta forma se especifica el área transversal a través de la cual fluye el caudal

de humo de cada caso y así mismo se especifica la altura de la fuente, variables
que influyen en el comportamiento y características de la dispersión.

3.2. ENMALLADO – Ansys CFX – Mesh

Ya que el modelo de turbulencia que se utiliza en esta simulación implica el

planteamiento de ecuaciones diferenciales que van a ser resueltas por medio de
un método numérico, el volumen generado previamente debe ser discretizado, es
decir, dividido en pequeños sub volúmenes llamados elementos. En el ámbito de
la comunidad CFD, este procedimiento recibe el nombre de enmallado.

Tal como se expuso anteriormente la cantidad de elementos es directamente

proporcional a la precisión de los resultados y a la demanda computacional, por lo
que el tamaño de los elementos y por ende la cantidad de los mismos tuvo que
balancearse con el recurso de cómputo disponible.

ANSYS CFX cuenta con un módulo de enmallado en el que se pueden ajustar

diferentes parámetros para lograr una discretización acorde con la deseada, entre
esos parámetros se encuentra el tamaño de cada elemento, la geometría, la
uniformidad y distribución. De forma análoga se puede configurar el motor de
enmallado de forma tal que se haga un refinamiento en sectores determinados
(disminución localizada del tamaño de los elementos).

12
HGV, por sus siglas en inglés, Heavy Goods Vehicle
46
Se creó un refinamiento en las zonas contiguas a la fuente y se creó una capa de
elementos semicilíndricos contra las paredes y el piso del túnel para obtener
resultados más precisos. A su vez se llevó a cabo un proceso iterativo para
encontrar la configuración más adecuada, en el que se analizaban las diferencias
de las mallas generadas sobre los resultados de la simulación, de esa manera se
encontró un balance entre los parámetros de mallado que dio como resultado una
división de la geometría de aproximadamente 118.000 elementos y 29.000 nodos.

En la figura 16 se puede apreciar el cambio de los resultados cuando se

implementa un refinamiento de la malla en la zona contigua a la fuente de emisión
de humo, y se tiene una transición básicamente de 26.600 a 118.000 elementos
en total:

Figura 16. Efectos del refinamiento de la malla

3.3. ESPECIFICACIÓN – Ansys CFX – Pre

En este módulo se ingresa la información necesaria para la simulación. Se

especifican submodelos, materiales, componentes, dominios, subdominios,
interfases y condiciones de frontera. Dentro de la información ingresada, se va a
precisar básicamente en cuatro (4) temas que se consideran los más relevantes
dentro de la caracterización de los escenarios estudiados: Parámetros de los
modelos, Variable adicional, estado transiente, especificación de escenarios.

3.3.1. Parámetros Modelos

3.3.1.1. Modelo de Turbulencia κ −ε

47
El modelo de turbulencia utilizado requiere la especificación de algunos
parámetros relacionados con las características del fluido y de la geometría en la
que se encuentre, en este caso se especificaron 2 parámetros dentro del modelo
κ − ε .

El primero es conocido como la escala de longitud de la turbulencia, que es una

cantidad física que describe el tamaño de los remolinos grandes en un flujo
turbulento. Generalmente esta escala de longitud se usa para estimar las
propiedades turbulentas de las entradas en una simulación CFD, su valor está
relacionado con la dimensión del problema y se define como una relación entre los
parámetros propios del modelo con una constante de proporcionalidad. (CFD
Online)

Con base en la escala de longitud turbulenta para el flujo totalmente desarrollado

en una tubería éste parámetro para el caso de la dispersión de humo en un túnel,
se relacionó con las dimensiones promedio de la fuente de emisión como un 10%
de su longitud.

El segundo parámetro especificado fue la intensidad de la turbulencia también

llamada nivel de turbulencia, es la relación entre la media de las fluctuaciones de
velocidad turbulenta y la velocidad media. Para ingresar el valor de este parámetro
se tomaron como referencia rangos experimentales de acuerdo al tipo de fluido y a
las características de la geometría como sigue: (CFD Online)

• Turbulencia Alta: Flujo a alta velocidad en geometrías complejas (Turbinas,

Compresores). Valores de intensidad entre 5% y 20%.
• Turbulencia Media: Flujo en geometrías no complicadas como tuberías o
ductos de ventilación. Valores de intensidad entre 1% y 5%
• Turbulencia Baja: Velocidades de flujo muy bajas. Valores de intensidad por
debajo de 1%.

3.3.1.2. Modelo de Transferencia de Calor.

Estos modelos son usados para predecir la temperatura en los diferentes nodos
de la geometría de estudio, en ellos se incluye la transferencia de calor por

48
conducción, por convección y de ser necesario se implementa el mezclado
turbulento y el trabajo viscoso.

Los dos modelos más relevantes en cuanto a las características de este proyecto
se conocen como: (1) Energía Térmica y (2) Energía Total.

Entre los anteriores el modelo más adecuado y por lo tanto el que se seleccionó
para la simulación, es el modelo de Energía Térmica. Este difiere del modelo de
Energía Total en la omisión de términos de aportes de energía cinética del flujo,
por lo que es aplicable a flujos de velocidades bajas donde los efectos cinéticos
son despreciables. Diferentes mediciones de velocidad de flujo en túneles han
sido reportadas, mostrando valores no mayores a 10 m/s por lo que la aplicación
modelo expuesto anteriormente es válida.

Si bien los efectos cinéticos no se tienen en cuenta, el modelo de Energía Térmica

es más robusto en el planteamiento y desarrollo matemático de las ecuaciones de
conservación, lo cual es una ventaja sobre el modelo de Energía Total.

3.3.2. Variable Adicional

Como el humo producido en un incendio, es una combinación volátil de partículas

sólidas de hollín, combustible parcialmente quemado, una mezcla de aire caliente
y gases tóxicos, y las proporciones en la que se encuentra cada uno de los
componentes varía significativamente de un caso a otro; para la simulación se
planteó el humo como un gas resultado de la mezcla homogénea de diferentes
gases, utilizando los datos expuestos por (Commission on Engineering and
Technical Systems, National Research Council, 1991) para la combustión de
poly(n-phnyleneisophthalamide).

La siguiente tabla ilustra los compuestos con mayores aportes en la ponderación

de propiedades físicas como la densidad.

Tabla 3: Composición del humo.

Concentración Concentración 3
Componente ρ (Kg/m )
(mg/g) % másico
CO 32 18,76% 1,145
CO2 56 32,83% 1,842

49
Agua 53 31,08% 1,607
Benceno 7,1 4,16% 2,778
Tolueno 1,9 1,11% 3,041
Benzonitrilo 13,6 7,97% 2,790

Teniendo estos datos, se creó una variable adicional volumétrica, con densidad
constante 1,691kg/m3, para representar el humo producido por el incendio,
teniendo además en cuenta que el 18,76% de la mezcla es Monóxido de Carbono
(CO).

3.3.3. Estado Transiente

Cuando se tienen situaciones en las que el comportamiento de las variables

cambia con el tiempo, se hace necesario un análisis en estado no estacionario, ya
que en muchos casos resulta más relevante conocer la transición de un estado
inicial al estado estable, como es el caso de la dispersión de un gas tóxico en un
espacio confinado.

En este caso, el tiempo es dividido en pasos discretos (de la misma manera que el
espacio es subdividido en pequeños volúmenes). El paso de tiempo, debe ser
escogido con cuidado, ya que éste debe ser suficientemente pequeño para
reproducir las características transientes importantes del flujo, y además debe ser
acorde con el refinamiento de la malla; entre más pequeño sea el tamaño de los
volúmenes, más corto debe ser el paso del tiempo.

De esta forma es posible representar el período de crecimiento de un incendio y

cualquier cambio en los patrones de flujo de aire con el tiempo. Estas
simulaciones son útiles en particular para investigar la eficiencia de los sistemas
de ventilación de emergencia, la extensión de humo y de calor, entre otras
características de interés para los investigadores de seguridad contra incendio. Sin
embargo, se debe tener en cuenta que para llevar a cabo la simulación, es
necesario conocer previamente información que caracterice de alguna manera el
incendio con el que se va a trabajar, ya que el CFD no puede predecir, entre otros,
la tasa de liberación de calor intrínseca de un incendio, y es por esto, que se debe
recurrir a utilizar los escenarios de incendio diseñados. (Gobeau, Ledin, & Lea,
2002).

50
Un escenario de incendio diseñado en un túnel proporciona los datos de entrada
para cualquier tipo de cálculo requerido para realizar predicciones y evaluación de
los peligros térmicos y tóxicos para las personas, así como de los daños térmicos
y no térmicos a estructuras y equipos relacionados. El diseño de un escenario de
incendio se puede obtener a partir de métodos determinísticos, donde la curva de
diseño representa el peor incendio concebible, o puede estar basado en el análisis
de riesgo que representa un escenario que está relacionado con datos
experimentales y estadísticos. (Ingason, 2009).

La forma habitual de representar un escenario de incendio, es escogiendo una

HRR máximo de una tabla de diseño.

En realidad los incendios diseñados, representan una idealización de un incendio

que podría ocurrir de verdad, y por esto, las curvas de incendio diseñadas, que se
utilizan para túneles, incluyen diferentes velocidades de crecimiento del incendio
2
(p.e. crecimiento lineal ( α t ), crecimiento cuadrático ( α t ), o exponencial
(α(1− e
−t
)
)e−t . Estas funciones de crecimiento se pueden combinar a su vez con un
valor máximo de HRR y una función de decaimiento (α −t o e −t ).

Las curvas de incendio diseñadas se pueden dividir en tres (3) tipos principales:
curvas lineales, curvas cuadráticas y curvas exponenciales. En cada caso esta
descripción puede referirse tanto al periodo de crecimiento como al de
decaimiento y puede ser combinada de diferentes formas. En la figura 17 se
muestra un ejemplo de cada uno de los tipos de curvas mencionados.

Figura 17. Ejemplos de los tres (3) tipos de curvas de incendios diseñados.
51
Una forma simple, flexible y confiable para describir una curva de incendio
diseñada, mediante una sola expresión matemática, fue presentada por (Ingason,
2009), y es la utilizada en este trabajo. El diseño consiste en una curva de
incendio exponencial que incluye una velocidad de crecimiento y un periodo de
decaimiento, también exponenciales y un periodo de tiempo con un valor
constante de HRR.

La ecuación general que describe una curva de diseño de incendio presentada es:

(
Q& ( t ) = Q& max n 1 r1 (1 − e − k1t ) n1 −1 e − k1t + n 2 r2 (1 − e − k 2 t ) n 2 −1 e − k 2 t ) (19)

Q&
Donde max es el máximo HRR en MW y para determinar las 6 constantes
desconocidas ( n1 , n 2 , r1 , r2 , k 1 , k 2 ), se integra la Ecuación 19 y se encuentra una
correlación entre Etot (Liberación total de energía), Q& max , n, r (Ecuación 20), por otro
lado, se asume que existe una constanteη que relaciona k 1 y k 2 , y con esto se
tiene para k 1 y k 2 , las Ecuaciones 21 y 22 respectivamente.
r r 
E tot = E tot ,1 E tot , 2 = Q& max  1 + 2  (20)
 k1 k 2 
Q& max
k1 = (r1 + η r2 ) (21)
E tot
Q& max  r1 
k2 =  + r2  (22)
Etot  η 

En las ecuaciones anteriores, k representa el coeficiente de ancho del tiempo, n

es el índice de retardo y r el coeficiente de amplitud de la curva, teniendo en
cuenta esto, y para hallar los valores de n , r con los que se logra obtener una
curva con un periodo máximo constante, se utiliza una hoja de cálculo clásica, y
finalmente se encuentran 10 casos (Tabla 4), de los cuales, el número 3 (
n1 = 3; n 2 = 9 , 710 ; r1 = 2 ,135 ; r2 = 1, 730 ), es el que mejor se adapta a la curva
cuadrática y por lo tanto fueron los valores utilizados durante este trabajo; el valor
de η que se encontró válido para los 10 casos, es igual a 10/7.

52
Tabla 4: Valores encontrados de n, r (Ingason, 2009)
Curva nr n1 n2 r1 r2
1 1,5 4,757 1,654 1,530
2 2 6,495 1,894 1,612
3 3 9,710 2,135 1,730
4 4 12,220 2,231 1,741
5 5 14,610 2,290 1,760
6 6 16,865 2,330 1,774
7 7 19,040 2,360 1,788
8 8 20,980 2,376 1,789
9 9 22,853 2,387 1,789
10 10 24,675 2,398 1,792

Finalmente cabe mencionar que la curva exponencial presentada es una

herramienta muy poderosa para la creación de curvas, cuando los parámetros
físicos ( Q& max , E tot ), son conocidos, y de esta forma, variando los parámetros
E tot , n , es posible ajustar el periodo de crecimiento a casi cualquier tipo de
incendio (Ingason, 2009).

Ahora, teniendo en cuenta que las funciones que se debían especificar en la

simulación eran las que relacionaban el comportamiento de la temperatura y el
flujo másico con el tiempo, se utilizó el mismo análisis planteado anteriormente
para realizar una modificación a la curva descrita en la Ecuación 19, y así poder
enocontrar las curvas buscadas (T vs t y w vs t), esto es posible debido a que el
comportamiento de las tres variables mencionadas está directamente relacionado.
(Babrauskas, 2009)

De esta forma es como la Ecuación 19 se convierte en la Ecuación 23 presentada

a continuación, en la cual se tienen dos nuevos parámetros a X y bX , los cuales
deben ser encontrados para ajustar la curva a los datos experimentales con los
que se cuente.

( (
X ( t ) = a X ⋅ Q& max n 1 r1 (1 − e − k1t ) n1 −1 e − k1t + n 2 r2 (1 − e − k 2 t ) n 2 −1 e − k 2 t )) + b X (23)

53
En esta ecuación, X puede representar, el HRR (Qmax en MW), la temperatura ( T
en K), o el Flujo másico ( w en Kg/s).

Finalmente es con este procedimiento como se especifican las ecuaciones en la

simulación, para el comportamiento de la temperatura y el flujo másico en el
tiempo. Que son las variables más determinantes en la dinámica del humo dentro
del túnel. Los valores de Q& max , Etot ,wmax , para cada uno de los casos fueron
tomados de (Ingason, 2009) y (Mégret & Vauquelin, 2000). Estos valores, al igual
que el tiempo total de duración del incendio, para cada uno de los casos, y el
momento en el que se desarrolla totalmente cada incendio se presentan en la
Tabla 5.

Tabla 5. Parámetros de las curvas de diseño

Qmax w max Tiempo Tiempo HRR
Etot (KJ) (MW) Tmax (°C) (Kg/s) Total (min) max (min)
Automóvil 17 8 210 55,264 120 10
Van 38 15 400 92,106 120 15
HGV 450 100 1282 460,53 180 25

3.3.4. Especificación de escenarios

Se diseñaron diferentes escenarios de incendio variando parámetros de entrada

como: Tamaño de la fuente, Flujo Máximo de material emitido, Temperatura
máxima de los gases de combustión y Velocidad de Ventilación. Es importante
resaltar y tener en cuenta que todos los valores con excepción de la velocidad del
viento, son valores tomados de la literatura, los cuales han sido resultados de la
experimentación, y se utilizan para poder tener un orden de magnitud para cada
uno de los casos, pero en ninguna circunstancia podran ser tomados como regla
general para cada tamaño de vehículo; como se ha mencionado a lo largo del
documento, estos datos varían de manera significativa dependiendo de la
combinación de factores presentes en cada caso.
Los datos utilizados para la velocidad del viento son tres (3) valores considerados
dentro del rango “normal”, esto referido a que son velocidades de viento
comunmente encontradas al interior de los túneles de carretera, y adicional a esto,
se utilizaron dos (2) valores correspondientes a condiciones de ventilación críticas,
una velocidad muy baja (0,2m/s) y una muy alta (3m/s), para ilustrar de una mejor
manera la importancia que tiene esta variable dentro de la dinámica del fenómeno.

54
En las Tablas 5 y 6 es se muestran los valores utilizados para los diferentes casos
estudiados.

Tabla 6. Velocidades de viento utilizadas

Velocidad Ventilación 1 2 3 4 5
(m/s) 0,2 0,5 0,9 1,85 3

3.4. RESOLVEDOR – Ansys CFX – Solver

En el modulo de Solver se especifican parámetros relacionados con grado de

precisión deseado, lo cual quiere decir que los criterios de convergencia allí
ingresados, afectan directamente el tiempo necesario para obtener la solución, por
lo cual hay que especificar valores con los cuales se logre un balance entre el
tiempo de cómputo y el nivel de confiabilidad de los resultados.

Las simulaciones se corrieron en un computador con procesador de doble núcleo

y 2 GB de memoria ram (características que definen la capacidad de cómputo),
tardando un promedio de 4 minutos para estado estacionario y 3 horas para
estado transiente, estos tiempos pueden disminuir considerablemente si la
capacidad de cómputo de la máquina donde se corran es incrementada.

3.5. RESULTADOS – Ansys CFX – Post

El módulo de visualización de ANSYS CFX se caracteriza por su versatilidad en

cuanto a la presentación de los resultados obtenidos. Es posible relacionar todas
las variables involucradas en el fenómeno de una forma práctica para darle una
buena orientación al posterior análisis.

Con base en los criterios de peligrosidad mencionados anteriormente, se

escogieron tres (3) variables principales para analizar los resultados: Temperatura,
Concentración y Velocidad. De acuerdo a la naturaleza de cada variable se
utilizaron visualizaciones tales como planos, volúmenes, vectores, y se realizaron
animaciones para observar el comportamiento en el tiempo de cada una de las
mismas.

55
 4. RESULTADOS Y DISCUCIÓN

Es importante resaltar que esta simulación se lleva a cabo bajo condiciones de

total confinamiento, ya que es un primer análisis en el cual no se tiene en cuenta
la inclusión de sistemas de ventilación y/o extracción de humo sino una sola
velocidad longitudinal del viento, lo que se ve reflejado en los resultados
mostrados y analizados a continuación.

4.1. CONCENTRACIÓN

Para analizar el perfil de concentración, se tuvo en cuenta el tiempo comprendido

entre el inicio del incendio y el momento en el que se desarrolla totalmente,
teniendo en cuenta que la información obtenida de este lapso de tiempo es de
vital importancia para el diseño de sistemas de alarma eficientes que permitan
identificar la presencia de fuego/humo, y así llegar a tiempos cortos de
evacuación. La evacuación se debe realizar inmediatamente después del inicio del
incendio ya que durante ese tiempo los niveles de concentración de humo dentro
del túnel, aguas abajo del incendio, hasta ahora están comenzando a subir, lo que
podría permitir que una persona se desplace a un lugar seguro.

Ya que los resultados se presentan en unidades de kg/m3 y los límites de

peligrosidad presentados en el numeral 2.1.4.4. están dados en partes por millón
(ppm), se estableció una relación lineal entre ambas magnitudes teniendo en
cuenta que el 18,76% del humo es CO (Commission on Engineering and Technical
Systems, National Research Council, 1991).

 kg 
ppm = 163 .912 ⋅  3  + 25 .537 (23)
m 

4.1.1. Tamaño de la Fuente de Emisión.

Se comparó la influencia del tamaño de la fuente bajo las mismas condiciones de

ventilación, es importante resaltar que el tamaño de esta, está relacionado con la
cantidad de materia que se libera y con la temperatura a la que ésta es emitida.

56
 Figura 18. Perfil de concentración de humo, Automóvil 0,9m/s, 3 y 5min.

Figura 19. Perfil de concentración de humo, Van 0,9m/s, 3 y 5min.

Figura 20. Perfil de concentración de humo, HGV 0,9m/s, 3 y 7min.

De los resultados obtenidos, y teniendo en cuenta los límites de tolerancia

mencionados anteriormente, se puede determinar un rango de peligrosidad para
cada uno de los casos: Para el automóvil, se tiene que a los 3 minutos de haber
iniciado el incendio, los niveles peligrosos de concentración ya se encuentran a

57
unos 80 metros de la fuente aguas abajo aproximadamente, lo cual quiere decir
que las personas que estén cerca a esta zona, tienen alta probabilidad de sufrir
narcolepsia y una posterior intoxicación letal, por efectos del CO. Se puede
observar como esta distancia es mucho mayor a medida que el tamaño de la
fuente se incrementa; como es en el caso del HGV, para el cual la liberación de
humo es un 1250% mayor que la del automóvil, y en este el rango peligroso ya
está más allá de 172m de la fuente.

Adicionalmente, se puede observar que antes de alcanzar el flujo máximo de

liberación de humo en el incendio, ninguna persona podría sobrevivir más de 5
minutos, en los 172 metros, aguas abajo del incendio, que son visibles en estas
imágenes.

Esta información es potencialmente útil para la ubicación de refugios y para el

diseño de los planes de evacuación.

4.1.2. Velocidad del viento.

Teniendo en cuenta que uno de los parámetros de diseño más importante en el

campo de la ingeniería de seguridad contra incendio en túneles, es la velocidad de
la ventilación, se plantearon diferentes configuraciones.

Figura 21. Perfil de concentración de humo, HGV 0,5m/s, 3 y 25min.

58
 Figura 22. Perfil de concentración de humo, HGV 1,85m/s, 3 y 25min.

En el caso del HGV, es posible observar que trascurridos 3 minutos, una velocidad
de viento de 1,85 m/s produce una distribución más uniforme, que una velocidad
menor (0,9 o 0,5m/s), lo cual indica que hay una mayor área del túnel expuesta a
una concentraciones peligrosas de humo, sin embargo a medida que avanza el
tiempo y se continua desarrollando el incendio, las concentraciones globales
correspondientes serán menores.

Es interesante observar que pasados 25 minutos (incendio totalmente

desarrollado), el perfil de concentración de humo, con una velocidad de 1,85m/s
es muy similar al comportamiento observado transcurridos solo 7 minutos de
iniciado el incendio cuando se tiene una velocidad de 0,9m/s, lo que confirma que
una velocidad mayor, contribuye a la evacuación de humo, de una manera
significativa.

Para hacer más objetivo el análisis, se tomo el caso de la van como un escenario
base, para a partir de este, agregar más niveles a las variables y de esta forma
producir resultados más relevantes. Entre los casos adicionales estudiados se
encuentran 2 velocidades extremas reportadas en la literatura (0,2m/s (Hua, Huo
a, & Chow b, 2008), 3m/s (Migoya Valor, 2002)).

59
 Figura 23. Perfil de concentración de humo, Van 0,2m/s, 3 y 15min.

Figura 24. Perfil de concentración de humo, Van 3m/s, 3 y 15min.

Analizando los perfiles de concentración obtenidos con condiciones extremas de

ventilación, se confirma lo descrito anteriormente con el caso de HGV.

Finalmente es importante destacar la importancia de la ventilación, ya que

teniendo una velocidad de viento pequeña, las consecuencias de este incendio
pueden ser equivalentes a tener un incendio de dimensión mucho mayor.

4.1.3. Incremento del Flujo másico de humo

Debido a que las características físicas de cada vehículo varían, y por lo tanto la
dinámica de su combustión, se planteó una variación de ± 15% en el flujo másico
máximo emitido por la fuente, para evaluar la sensibilidad del comportamiento del
humo ante esta variable. A partir de esto se observó que es necesario una
variación mayor ( ±30%) para poder tener cambios significativos; tanto en el
llenado progresivo del túnel, como en su concentración máxima.
60
 Figura 25. Perfil de concentración de humo, Van 0,9m/s; -15%w y +15%w.

4.2. TEMPERATURA

Para analizar el perfil de temperatura, se tuvo en cuenta todo el periodo de tiempo

que dura el incendio desde la ignición hasta su extinción, esto, con el objetivo de
obtener información que permita llevar a cabo un eventual análisis no solo a las
consecuencias para vidas humanas sino a los posibles daños sobre la
infraestructura por efectos térmicos.

Los resultados de los perfiles de concentración y temperatura tienden a ser

similares, lo cual se debe en principio a que la función que modela el estado no
estacionario de cada una, cambia en unidades y valores, mas no en tendencia (ver
numeral 4.3.2), a su vez se debe a que físicamente cuando se tiene una mayor
concentración de humo la proporción que está mezclada con el aire es baja y por
lo tanto el transporte de calor también lo es, por lo que una nube de humo muy
concentrada tendera a tener elevadas temperaturas.

Es importante tener en cuenta que la temperatura obtenida en estos perfiles,

incluye los efectos por radiación y convección.

4.2.1. Tamaño de la Fuente de Emisión.

Se comparó la influencia del tamaño de la fuente sobre el perfil de temperatura.

61
 (a) (b)

(c)
Figura 26. Perfil de Temperatura 0,9m/s: (a) Auto 10 minutos, (b) Van 15 minutos, (c)
HGV 25minutos. Incendios totalmente desarrollados

Tal como se esperaba, entre mayor sea el tamaño de la fuente, mayores van a ser
las temperaturas alcanzadas a lo largo del túnel.

Un criterio adecuado para establecer daños sobre seres humanos por

temperatura, está descrito en la figura 7, y según esto, se puede afirmar que para
el único caso en el que las personas puede resistir al menos 5 minutos sin sufrir
quemaduras en la piel bajo condiciones de incendio totalmente desarrollado,
estando a una distancia no menor de 40 metros, aguas debajo de la fuente, es
cuando se incendia el automóvil, ya que en los otros casos, las temperaturas
alcanzadas son demasiado altas y se requiere estar a una distancia que supera el
área de análisis (172 m, aguas abajode la fuente), para evitar sufrir de
quemaduras.

4.2.2. Velocidad del viento.

Se realizó la comparación de los casos extremos de de velocidad de viento, tal

como se describió en el análisis de la concentración. Y se determinó que en el
62
caso de una ventilación de 0,2m/s, la temperatura promedio a partir de 40 m
aguas debajo de la fuente, en el momento en el que el incendio está totalmente
desarrollado, aumenta en un 57% aproximadamente, mientras que en el caso de
una velocidad mayor de ventilación (3m/s), la temperatura obtenida en la misma
área disminuyó un 28%, esto, con respecto al caso base (velocidad 0,9 m/s).

(a) (b)
Figura 27. Perfil de Temperatura, Van 15min; (a) 0,2m/s, (b) 3m/s. Incendios totalmente
desarrollados

4.2.3. Incremento del la temperatura.

De manera análoga a lo que se observó con los perfiles de concentración, al variar

± 15% los valores máximos de temperatura, no se tuvieron cambios relevantes
en los perfiles obtenidos, de hecho, el cambio es aún menos notorio con la
temperatura que con la concentración. Ya teniendo en cuenta la variación de 30%,
se evidencia un cambia de un poco menos del 25 % en las temperaturas globales,
valores que ya comienzan a tener importancia en las consecuencias observadas.

Figura 28. Perfil de concentración de humo, Van 0,9m/s; -15%T y +15%T.

63
4.3. VELOCIDAD

Para analizar la velocidad (magnitud y dirección) en diferentes puntos del túnel, se

construyeron dos planos: un transversal ubicado exactamente en el centro de la
fuente y otro longitudinal, que también la atraviesa, y sobre estos se observan los
vectores.

4.3.1. Tamaño de la Fuente de Emisión.

(a)

(b)

(c)

Figura 29. Velocidad - 0,9m/s, (a) Auto 10 minutos, (b) Van 15 minutos, (c) HGV 25
minutos. Incendios totalmente desarrollados
64
La carácteristica más evidente es que la magnitud de la velocidad es proporcional
al tipo de fuente, ya que el caudal aumenta en mayor medida que el área
transversal de un tamaño de fuente a otro.

Otra característica importante tiene que ver con la altura de la fuente de emisión,
ya que como se puede observar, el HGV se encuentra a tan solo 3m del techo,
mientras que el vehículo se encuentra a 5,7, lo que se manifiesta en un aumento
en la velocidad cercana a las paredes para la fuente que está más cercana al
techo (ver figura 19).

4.3.2. Velocidad del viento.

(a) (b)
Figura 30. Velocidad Van 15min: (a) 0,2m/s (b) 3m/s. Incendios totalmente desarrollados

En la figura se aprecia claramente la influencia que tiene la velocidad del viento

sobre la dirección y la magnitud de la velocidad de las partículas en el interior del
túnel.

4.4. BACK LAYERING

En cada uno de los casos analizados anteriormente, se puede evidenciar de una

manera clara la presencia del fenómeno conocido como Back Layering (ver
sección 2.1.3.1.), y la influencia que tiene el mismo en el comportamiento del
humo aguas abajo del incendio.

En las figuras donde se presentan los perfiles de concentración, se observa un

plano ubicado a 10m aguas arriba de la emisión, en el que se puede apreciar la

65
ocurrencia de este fenómeno para los casos en los que se tiene una velocidad de
viento baja, un flujo de materia emitida alto. Este es el caso de HGV con una
velocidad de 0,5m/s presentado en la figura 21, en el cual una alta concentración
de humo atraviesa este plano, indicando que la capa superior de humo que se
devuelve va mucho más allá de los 10 metros que se están observando; sin
embargo, este fenómeno puede ocurrir cuando las velocidades de ventilación son
mayores, si el flujo de humo proveniente de la fuente es lo suficientemente alto
(ver figura 22).

La velocidad crítica, varía de acuerdo a cada configuración de incendio, es decir,

que para cada combinación de HRR y velocidad de ventilación, va a existir una
velocidad crítica diferente.

Como ya se ha mencionado, la concentración de humo y la temperatura tienen un

comportamiento análogo, por lo que los casos en los que se presenta un
incremento de la concentración aguas arriba, también se va a evidenciar un
incremento de la temperatura en esta misma zona. (figuras 23 y 27 (a)).

En las figuras de velocidad también se evidencia de manera clara este fenómeno,

y éste es representado por un remolino en la zona superior del túnel, aguas arriba
de la fuente. El tamaño del remolino formado depende del tamaño de la fuente de
emisión y de la elevación de la misma con respecto al nivel del suelo.

La representación más clara de este fenómeno, se puede observar en la figura 30,

en donde por un lado, con una velocidad de 0,2m/s el remolino formado tiene un
tamaño considerable y alcanza a cubrir un gran porcentaje de la altura total del
túnel; y por otro lado, con una velocidad de 3 m/s no se forma remolino alguno, y
por el contrario, la corriente de humo solo se mueve aguas abajo.

66
 5. CONCLUSIONES

ANSYS-CFX demuestra ser una excelente herramienta para la simulación de

fenómenos físicos tales como la mecánica de fluidos turbulentos. Para el caso
específico que se presenta en este documento resultó ser adecuado, y se
encontró que es una herramienta potencialmente confiable.

Dado que este estudio se realizó bajo condiciones de confinamiento total, los
resultados evidencian la importancia de los sistemas de ventilación y extracción en
un túnel. Un diseño adecuado de estos sistemas garantiza que las estrategias de
prevención y mitigación de consecuencias sean funcionales y efectivas.

La velocidad crítica es diferente en cada escenario dadas las diferentes

combinaciones entre flujo emitido y velocidad del viento, sin embargo se mostró
que una ventilación de 3m/s garantiza que el fenómeno de Back Layering no se
presente. Este valor si bien aplica para todos los casos puede exceder la
velocidad crítica requerida en diferentes configuraciones.

El tiempo que transcurre desde el inicio del fuego hasta el momento en que este
se encuentra totalmente desarrollado es de vital importancia para el análisis de
consecuencias por intoxicación, ya que es en este lapso de tiempo donde las
personas tendrían la oportunidad de buscar refugio. Si bien en un principio una
alta velocidad de ventilación distribuye el humo de forma homogénea, más
adelante cuando el incendio se haya desarrollado completamente, las
concentraciones globales serán menores que en el caso de una velocidad de
ventilación baja.

Para una geometría fija los cambios en los resultados con alteraciones de 15% en
variables de entrada como el flujo máximo de material emitido y la temperatura
máxima alcanzada, no son relevantes. Sin embargo para un cambio del 30% ya
se aprecian efectos sobre la dinámica del fenómeno, a su vez se observó que la
sensibilidad de la concentración a lo largo del túnel es mayor que la de la
temperatura.

67
Las simulaciones numéricas tienen un gran potencial en el área de la modelación
en ingeniería, sin embargo se debe ser muy cuidadoso en su ejecución y posterior
análisis. Es necesario comprender el fenómeno físico así cómo los métodos de
resolución para desarrollar una simulación confiable.

68
 6. RECOMENDACIONES

Los modelos y métodos presentados durante el documento no dejan de ser una

aproximación y por lo tanto, los resultados obtenidos no serán exactos bajo
ninguna circuntancia. Errores numéricos y físicos se pueden presentar en cada
una de la etapas de la simulación (Gobeau, Ledin, & Lea, 2002): datos mal
ingresados (propiedades físicas, dimensiones, geometría), selección errada de
elementos (falta de comprensión de las bases teóricas), selección de tamaño y
forma de los elementos (precisión de los resultados) y manejo de condiciones de
frontera. (Moaveni, 1999)

Por lo tanto, una forma de controlar estos errores, y obtener resultados coherentes
para realizar un análisis ajustado fenomenológicamente, es comprender de
manera clara el fenómeno físico, las suposiciones del modelo y el método de
resolución aplicados. (Gobeau, Ledin, & Lea, 2002)

Finalmente, y teniendo en cuenta que el tema trabajado es muy amplio, además

que la informacion arrojada por las simulaciones realizadas es de gran importancia
en la busqueda de la prevención y mitigación de los efectos asociados a estos
eventos, se considera interesante proponer algunos contenidos con los que se
puede continuar trabajando en el tema, a partir de los resultados obtenidos con
este trabajo. En principio incluir los sistemas de ventilación y extracción de humo,
para observar como cambian las condiciones en el interior del túnel, seguido a
esto, estudiar la influencia de obstaculos (otros vehículos) en la dinámica del
humo, contexto que se aproxima mejor a un escenario real; tercero, estudiar la
propagación del incendio, y por último, contando con toda esta información,
realizar un análisis riguroso de consecuencias sobre personas y estructuras.

69
 BIBLIOGRAFÍA

Aertia Software. (2005). Especial CFD y Phoenics. AXIS (6).

Anderson, J. (1995). Computational Fluid Dynamics: The Basics with Aplications.

Mcgraw Hill.

ANSYSInc. (s.f.). Ansys CFX 12.0 Manual. ANSYS CFX-Solver, Release 12.0:
Theory .

Babrauskas, V. (Octubre de 2009). InterFIRE Online. Obtenido de

http://www.interfire.org/features/heat_release.asp

Casal, J., & Montiel, H. (2001). Dispersión de nubes tóxicas o inflamables. En A. d.

Industriales.

CFD Online. (s.f.). Turbulence Intensity. Recuperado el Octubre de 2009, de CFD

Online: http://www.cfd-online.com/Wiki/Turbulence_intensity

CFD Online. (s.f.). Turbulence Length Scale. Recuperado el Octubre de 2009, de

CFD Online: http://www.cfd-online.com/Wiki/Turbulence_length_scale

Commission on Engineering and Technical Systems, National Research Council.

(1991). Fires in Mass Transit Vehicles: Guide for Evaluation of Toxic Hazards.
United States of America.

Gobeau, N., Ledin, H., & Lea, C. (2002). Guidance for HSE Inspectors: Smoke
Movement in Complex Enclosed Spaces - Assesment of Computational Fluid
Dynamics.

Hua, L., Huo a, R., & Chow b, W. (2008). Studies on buoyancy-driven back-
layering flow in tunnel fires. Experimental Thermal and Fluid Science , 1468 - 1483.

Ingason, H. (2009). Design fire curves for tunnels. Fire Safety Journal , 44, 259-
265.

Lees, F. P. (1996). Emission and Dispersion . En Lees' Loss Prevention in the

Process Industries. Hazard Identification, Assessment and Control (págs. 136 -
159).

70
Lees, F. P. (1996). Fire. En Lees' Loss Prevention in the Process Industries.
Hazard Identification, Assesment and Control (págs. 2-15).

Madenci, E., & Guven, I. (2006). The Finite Method and Aplications in Engineering
Using ANSYS. New York: Springer.

Mégret, O., & Vauquelin, O. (2000). A model to evaluate tunnel fire characteristics.
Fire Safety Journal , 34, 393-401.

Migoya Valor, E. (2002). Modelo Zonal para la Simulación del Movimiento de

Humos y Gases Calientes en Incendios: Aplicación a Túneles de Carretera. Tesis
Doctoral . Madrid, España.

Moaveni, S. (1999). Finite Elemen Analysis. New Yersey: Prencite Hall.

Modic, J. (2003). Fire Simulation in road tunnels. Tunnelling and Underground

Space Technology. 18, 525–530.

Neophytou, M. .., & Britter, R. E. (2005). A Simple Model for the Movement of
Smoke in a Confined Tunnel. Pure and Applied Geophysics , 162, 1941-1954.

NIST. (2009). Fire Dynamics Simulator (Version 5) User's Guide.

PIARC Committee in Road Tunnels. (1999). Fire and Smoke control in road
tunnels. PIARC.

Pitts, W. M. (1995). The Global Equivalence Ratio Concept and the Formation
Mechanisms of Carbon Monoxide in Enclosure Fires. Progress in Energy and
Combustion Science , 21, 197-237.

Rapporteur Alfred Haack, S. (2004). Design Fire Scenarios. FIT Technical Report -
Part 1 . Bruselas, Bélgica.

Shorab, J., Shashi, K., & Surendra, K. (2008). Numerical simulation of fire in a
tunnel. Comparative study of CFAST and CFX predictions. Tunnelling and
Underground Space Technology , 160-170.

Xue, H., Ho, J., & Cheng, Y. (2001). Comparison of different combustion models in
enclosure fire simulation. Fire Safety Journal , 36, 37-54.

71