El presente trabajo consiste en el estudio numérico del confinamiento de un campo de alta temperatura, producido por una fuente de calor intenso, por medio de cortinas de aire tipo doble jet-doble flujo (DJ-DF), emulando las condiciones de operación de la instalación experimental construida en el marco del proyecto Fondecyt 1040498. Para este efecto se utilizó el modelo de turbulencia Reynolds Stress Model (RSM). El diseño de un dispositivo de confinamiento de escalares activos mediante cortinas de aire y la simulación numérica de sus condiciones de operación con modelos turbulentos de orden superior han motivado este estudio, cuyo objetivo ulterior es proponer soluciones para la seguridad en túneles viales. La eficacia de confinamiento por barreras gaseosas abre las puertas para diseñar su aplicación tecnológica en el campo de la seguridad dentro de túneles viales, en los cuales se busca controlar el transporte de calor y los gases tóxicos producidos en caso de incendio. La atención del estudió se centró tanto en las cortinas como en el comportamiento global del dispositivo. Simulaciones isotérmicas y no-isotérmicas fueron realizadas. Se efectuó una comparación del modelo RSM con los modelos k- standard y V2F. Se realizaron simulaciones no-isotérmicas considerando diferentes condiciones de operación del dispositivo, variando los parámetros relevantes que gobiernan los fenómenos estudiados (número de Reynolds, potencia térmica). Estas comparaciones fueron complementadas con mediciones de temperatura obtenidas de la instalación experimental. Por último, se varió el dominio de cálculo de las simulaciones: se realizó una simulación tridimensional del dispositivo y se realizó una simulación bidimensional, pero disminuyendo el largo de los circuitos de recirculación. Las simulaciones permitieron caracterizar las distintas zonas de las cortinas DJ-DF. En las zonas cerca de las boquillas de salida las cortinas se comportan como doble jet, pero aguas abajo se comportan como un jet plano simple. Los términos de transferencia de calor turbulento permitieron identificar las zonas de inhibición y transferencia de calor por estos mecanismos. Se observó que las capas de mezcla y la zona de impacto favorecen de manera importante el transporte de calor turbulento. Al aumentar el número de Reynolds de las cortinas se obtuvo una disminución de la actividad turbulenta en el eje de la cortina favoreciendo así el confinamiento. Se concluye que el modelo RSM se adapta mejor que los modelos citados, justificando el aumento de los recursos computacionales por la calidad de los resultados. Las simulaciones 3D indican cambios poco significativos en el comportamiento de las cortinas con respecto a los 2D, sin embargo, se predicen temperaturas mayores en la zona confinada debido a efectos radiativos en las paredes laterales. Las cortinas son en esencia bidimensionales y los efectos 3D se limitan a la zona próxima al penacho generado por la fuente térmica. La disminución de la longitud de los circuitos de recirculación del dispositivo no afecta la habilidad de confinamiento de las cortinas, pero produce un aumento global de la temperatura en la zona confinada. Dada la dispersión de los resultados al comparar el campo térmico experimental, etapas futuras de investigación numérica deben considerar las pérdidas de calor por las paredes (no adiabaticidad). Además, si se desea simular condiciones reales de incendio dentro del dispositivo, el modelo de radiación ocupado debe considerar el efecto radiativo que produce el humo dentro del túnel, que debería afectar la temperatura y dinámica de los flujos.