Source: http://docplayer.es/718796-Anexo-5-reglas-basicas-para-la-implementacion-de-modelos-cfd-en-el-modelo-del-vertido-en-campo-cercano-v-1-1.html
Timestamp: 2018-10-18 01:44:21+00:00

Document:
ANEXO 5 REGLAS BÁSICAS PARA LA IMPLEMENTACIÓN DE MODELOS CFD EN EL MODELO DEL VERTIDO EN CAMPO CERCANO V PDF
Download "ANEXO 5 REGLAS BÁSICAS PARA LA IMPLEMENTACIÓN DE MODELOS CFD EN EL MODELO DEL VERTIDO EN CAMPO CERCANO V. 1.1"
Ángel Moreno Quintana
1 ANEXO 5 REGLAS BÁSICAS PARA LA IMPLEMENTACIÓN DE MODELOS CFD EN EL MODELO DEL VERTIDO EN CAMPO CERCANO V. 1.1
2 Este documento se ha obtenido como resultado del proyecto de I+D+i: MEDVSA. Desarrollo e implementación de una metodología para la reducción del impacto ambiental de los vertidos de salmuera procedentes de las plantas desaladoras 045/RN08/03.3 MEDVSA es un proyecto de I+D+i ( ) financiado por el Ministerio de Medio Ambiente, Medio Rural y Marino. Los trabajos han sido coordinados por el Instituto de Hidráulica Ambiental de la UC y desarrollados en colaboración con el Centro de Estudios de Puertos y Costas (CEDEX). La Guía Metodológica desarrollada en el proyecto MEDVSA incluye la documentación: - GUÍA METODOLÓGICA - Anexo 1. Fichas descriptivas de los modelos comerciales - Anexo 2. Modelos de simulación brihne (MEDVSA). - Anexo 3. Formulación basada en los resultados experimentales - Anexo 4. Caracterización de clima marino - Anexo 5. Reglas básicas para la implementación de modelos CFDs en la simulación de campo cercano - Anexo 6. Reglas básicas para la implementación de modelos hidrodinámicos en la simulación de campo lejano - Anexo 7. Medidas preventivas y Programa de Vigilancia Ambiental El presente documento constituye el Anexo 5: Reglas básicas para la implementación de modelos CFDs (Computational Fluid Mechanics) en la simulación de campo cercano. Autor: INSTITUTO DE HIDRÁULICA AMBIENTAL Personas de contacto: Javier L. Lara Iñigo J. Losada
3 ANEXO 5 1. INTRODUCCIÓN 1 2. COMPORTAMIENTO DEL VERTIDO HIPERSALINO EN EL CAMPO CERCANO Y SU MODELADO NUMÉRICO ESTADO DEL ARTE MODELADO CON OpenFOAM DESCRIPCIÓN DEL MODELO TWOLIQUIDMIXINGFOAM CONFIGURACIÓN NUMÉRICA RESULTADOS METODOLOGÍA A SEGUIR FUTURAS LÍNEAS DE INVESTIGACIÓN REFERENCIAS... 17
5 1. Introducción El trabajo incluido en este Anexo ha sido desarrollado por el Instituto de Hidráulica Ambiental de Cantabria. La necesidad de entender en detalle los procesos físicos asociados a los fluidos ha demandando en los últimos años el uso de los modelos basados en las ecuaciones completas de Navier-Stokes. La rama de la mecánica que resuelve este tipo de problemas se conoce como Mecánica de Fluidos Computacional (CFD, Computational Fluid Mechanics) y utiliza una serie de algoritmos y técnicas numéricas que permiten analizar problemas asociados a los fluidos, con alto grado de detalle y minimizando el numero de hipótesis de partida. Para ello, la técnica CFD resuelve directamente la mayoría de los términos de las ecuaciones de Navier-Stokes y las ecuaciones de transporte asociadas a éstas. Los términos no calculados directamente están principalmente asociados al flujo turbulento y se resuelven mediante el modelado y la parametrización de los procesos físicos. Su desarrollo comenzó al final de la década de los años sesenta para resolver problemas muy sencillos, pero hoy en día, debido al perfeccionamiento del las técnicas computacionales y al fácil acceso a sistemas de cómputo masivo, el uso de estas herramientas para resolver problemas complejos es bastante común para estudiar no sólo problemas científicos, sino para resolver problemas asociados a la ingeniería. Dentro del campo de los chorros hipersalinos el uso de las técnicas CFD permite el estudio en detalle de los procesos turbulentos asociados a la descarga, la incorporación de fluido del entorno ( entrainment ), disipación de la cantidad de movimiento por fricción con el fluido del entorno, los procesos desencadenados por las fuerzas gravitacionales en la caída de la pluma hipersalina, y el impacto de ésta con el fondo marino. Este tipo de técnicas permiten por un lado, definir y caracterizar mejor el campo de velocidades y de concentraciones asociados al vertido, de una manera más exacta que los tradicionales modelos integrados que usan unas hipótesis de partida muy rígidas, y por otro complementar las medidas realizadas en los modelos físicos. El proceso desarrollado en el campo cercano de los vertidos hipersalinos es muy complejo dado que en él intervienen escalas espacio temporales diferentes entre el comportamiento inicial de chorro y en el final de pluma. Dicha variedad de escalas y su doble comportamiento chorro-pluma hace que el modelado numérico sea muy difícil de realizar, incluso con las técnicas numéricas más avanzadas en el campo de las herramientas CFD. Este hecho justifica los pocos trabajos existentes en la literatura en este tipo de aplicaciones. En el presente anexo se realiza un estudio preliminar del comportamiento del vertido de salmuera en el campo cercano y se presenta una metodología de trabajo para su estudio mediante técnicas CFD. El trabajo se centra en el vertido mediante un dispositivo en chorro sumergido. Se ha elegido esta configuración en el estudio porque se trata de uno de los sistemas de descarga más eficaces desde el punto de vista de la dilución, que es uno de los sistemas más empleados en la plantas desaladoras recientes, y dado que es la configuración modelada físicamente mediante técnicas ópticas no intrusivas, en el marco del proyecto MEDVSA, por lo que se dispone de una gran base de datos para calibración y validación de los resultados numéricos. El anexo se estructura de la siguiente manera. Inicialmente se describe brevemente el comportamiento del vertido hipersalino en el campo cercano. Posteriormente, se detallan los trabajos existentes en el estado del arte relativos al uso de técnicas CFD en chorros con flotabilidad negativa. A continuación, se describen las simulaciones numéricas realizadas mediante el código OpenFOAM usado en este proyecto. Posteriormente se detallan las ecuaciones usadas por los modelos CFD y las ecuaciones de transporte que resuelven asociadas a la salmuera. Finalmente, se propone una metodología a seguir en el uso del modelo para el estudio de este tipo de procesos. 1
6 2. Comportamiento del vertido hipersalino en el campo cercano y su modelado numérico El comportamiento de un efluente vertido al mar depende fundamentalmente de la cantidad de movimiento transmitida en la descarga y de la diferencia de densidad entre el efluente y el medio receptor. Los efluentes de flotabilidad negativa son aquellos con mayor densidad que el fluido del medio receptor donde se vierten. En el caso de la salmuera procedente de las plantas desalinizadoras, su exceso de salinidad le confiere una mayor densidad que el agua marina. Atendiendo a los procesos físicos que el efluente desencadena y los procesos hidrodinámicos que aparecen, se puede zonificar el vertido hipersalino en dos regiones: el campo cercano y el campo lejano. La región situada en el entorno del punto de vertido, también llamada zona de mezcla inicial, es la que se conoce con el nombre de campo cercano. Este Anexo se centra en el modelado numérico de la hidrodinámica y los procesos de mezcla que se desencadenan en esta región, en particular, en el caso de una descarga mediante chorro. En el campo cercano se producen diferentes procesos físicos asociados con la mezcla entre el efluente y el medio receptor, que se resumen en un aumento de la dilución del efluente y una perdida de cantidad de movimiento del flujo del efluente, debido a su mezcla y fricción con el fluido receptor, que posee una cantidad de movimiento inferior. Los parámetros geométricos y de dilución del efluente responden a escalas espaciales y temporales pequeñas, del orden de metros y de minutos respectivamente. El vertido hipersalino mediante chorro tiene dos comportamientos hidrodinámicos diferenciados. En su estado inicial, presenta un claro comportamiento de chorro, donde la descarga dota al efluente hipersalino de un cierto impulso hidrodinámico inicial y genera un flujo turbulento que facilita el mezclado. La cantidad de movimiento inicial disminuye a medida que el fluido se aleja del chorro por efecto de la fricción con el fluido del medio receptor y la incorporación de la masa de éste en el cuerpo del chorro. En esta zona se distinguen dos zonas, la zona de desarrollo de los perfiles autosemejantes (ZOFE) y la zona de perfiles autosemejantes (ZFE). La ZOFE abarca una distancia desde el punto de vertido entre 5 y 6 veces el diámetro del chorro, y supone una zona de transición entre los perfiles cuasi-horizontales, típicos de un perfil turbulento en un tubo, a un perfil gausiano, característico de un chorro. Esta zona está dominada por un flujo irrotacional en el centro del chorro, y una alta fuerza de cizalla en el exterior, que realiza el trabajo de modificar la forma del perfil y crea las inestabilidades que desencadenan los procesos de mezcla ( entrainment ) y la generación de vórtices. Las escalas temporales y especiales en esta zona son muy pequeñas, inferiores al segundo y por debajo del diámetro del chorro. El modelado numérico de esta zona requiere una discretización muy elevada y modelos de turbulencia que resuelvan las inestabilidades con esquemas numéricos de alto orden para asegurar una disipación numérica reducida. El uso de modelos basados en la técnica de simulación de grandes remolinos (LES, Large Eddy Simulation) junto con el uso de una discretización de detalle, es la más adecuada para estudiar la capa exterior del chorro y su deformación. Resolver esta zona del chorro conlleva muy altos elevados costes de simulación. La ZFE se caracteriza por presentar unos perfiles estadísticamente semejantes de tipo gausiano, que en el caso de chorros hipersalinos no son simétricos debidos al efecto de las fuerzas gravitacionales que tienden a desplazar el chorro hacia el fondo. La hidrodinámica de esta zona se caracteriza por ser la región donde se producen los mayores procesos de mezcla y estar el flujo turbulento completamente desarrollado. El chorro pierde cantidad de movimiento, y el fluido del medio receptor se mezcla completamente con el chorro. Las escalas espacio temporales de esta zona son mayores que en la ZOFE, y se observan vórtices más grandes incluso que el tamaño inicial del diámetro del tubo. Este hecho hace que las técnicas de cómputo permitan realizar una discretización mayor a la hora de resolver esta zona, aunque los procesos turbulentos reducen el esfuerzo computacional. El uso de modelos tipo RANS (Reynolds Averaged Navier- Stokes) o LES es óptimo en esta región, siendo los primeros los que menos coste computacional requieren debido al tipo modelo de cierre que usan para la turbulencia y por requerir un mallado de menor tamaño. Al final de la zona ascensional del chorro, la fuerza dominante es la gravitacional, debido a que el chorro ha perdido prácticamente la totalidad de la cantidad de movimiento transmitida en la descarga. Alcanzado el punto de máxima altura del chorro, comienza su trayectoria descendente con hidrodinámico de pluma. En su descenso el efluente continúa mezclándose con el fluido del medio receptor debido a las fuerzas de cizalla que se generan por la diferencia de 2
7 velocidad entre el fluido receptor y la masa del efluente sobre la que actúa la gravedad. La anchura de los perfiles aumenta a medida que desciende por efecto del entrainment. El proceso de descenso es más lento que el de ascenso y como consecuencia las escalas temporales son inferiores, aunque las escalas espaciales son similares a las observadas en la zona final del ascenso. Al igual que en el caso anterior, tanto los modelos tipo RANS o LES se podrían aplicar a esta región, con un mayor coste computacional de estos últimos. Finalmente, el chorro impacta con el fondo creando una zona de mayor mezcla, por efecto del gradiente de presiones que impone el fondo, que origina un cambio en el flujo de cantidad de movimiento vertical a horizontal. El impacto con el fondo de la pluma de salmuera produce la formación de un nuevo vórtice que genera una advección del efluente de manera radial. Las escalas espaciales y temporales observadas en esta región son menores que las observadas en la rama descendente, y requieren un tratamiento más refinado. Al igual que en el caso anterior, tanto los modelos tipo RANS o LES se podrían aplicar a esta región, con un mayor coste computacional de estos últimos. La figura A5-1 muestra dos imágenes de ensayos mediante modelo físico de un vertido en chorro, donde se aprecian las diferentes zonas del chorro. Se puede observar cómo en la zona cercana la boquilla las escalas espaciales de los remolinos son inferiores que en las zonas alejadas de éste, tanto en la rama ascensional del chorro como en la de descenso. Figura A5-1. Ensayos de vertido mediante chorro individual. Arriba: vista lateral. Abajo: Vista frontal 3
8 3. Estado del arte El modelado CDF de un proceso de descarga de un vertido hipersalino mediante un chorro implica la consideración de muchos tipos de escalas espacio-temporales que interaccionan entre si, confiriéndole al problema una gran complejidad por el tipo de mallas, modelos de turbulencia y métodos de resolución a utilizar. Aunque el modelado numérico de los procesos individuales de chorros, plumas o impactos de estos con contornos sólidos está presente en la literatura, el número de publicaciones en las que se modele por completo un chorro hipersalino y sus diferentes estadios es muy limitado y está ligado principalmente al uso de modelos comerciales. Uno de los primeros trabajos presentes en la literatura es el presentado por Vafeiadou et al. (2005) [1], que estudia chorros de flotabilidad negativa inclinados 45º y 90º, utilizando para el modelado del flujo ecuaciones RANS y un modelo de cierre para la turbulencia tipo SST (Shear Stress Transport). Utilizan el código CFX-5, que resuelven sobre una malla tridimensional de elementos. La validación del modelo se realiza comparando magnitudes medias obtenidas experimentalmente por Bloomfield et al. (2002) [2] y por Roberts et al. (1997) [3]. Los resultados muestran un comportamiento bastante aceptable de los resultados aunque no se facilitan comparaciones de los flujos turbulentos de velocidades o de concentraciones, ni se muestra la interacción de la pluma descendente con el fondo ni se realiza un análisis en detalle del flujo calculado. Posteriormente, Olivier et al. (2008) [4] presentan un estudio realizado con el mismo paquete comercial CFX, sobre una malla de millones de elementos sobre las que resuelven las ecuaciones del fluido y del transporte del efluente hipersalino. El modelo usado para el cierre de la turbulencia es el modelo k- estándar, sobre el que se realiza un calibrado del numero de Schmidt turbulento. Este calibrado realizado a partir de ensayos de laboratorio muestra una mejora en la predicción de las magnitudes medias y de las variables asociadas a éstas, como es el caso de la trayectoria seguida por el chorro. Tampoco en este trabajo se muestra la interacción de la pluma descendente con el fondo ni se analiza en detalle el flujo calculado. Finalmente, Seil et al. (2010) [5] muestran un trabajo en el que se presentan simulaciones numéricas de vertido con un único chorro o con múltiples usando el código ANSYS FLUENT. El trabajo, más tecnológico que científico muestra el potencial del software comercial para ser usado como herramienta de diseño. En este caso se presenta el uso del modelo SST y RNG (Renormalization Group) para el modelado del flujo turbulento. La validación de los resultados numéricos es muy somera, y compara sus resultados con las fórmulas semi-empíricas de Roberts et al. (1997) [3]. Como se puede apreciar, el número de publicaciones científicas es muy pequeño, por lo que este tipo de modelado está es sus estadios iniciales. La validación, basada en los flujos medios en la mayoría de los casos no son lo exhaustiva que se precisa y no existe una uniformidad en el modelo utilizado para el cierre de los términos turbulento del flujo. Así mismo, aunque el modelado LES parece mas preciso y presente en otras disciplinas en las que se estudia la mecánica de fluidos computacional, no se han encontrado relativos al uso de dicha aproximación. 4
9 4. Modelado con OpenFOAM En este apartado se describe el trabajo realizado dentro del proyecto MEDVSA en el campo de la mecánica de fluidos computacional. Se ha elegido utilizar el modelo OpenFOAM, de código abierto, ampliamente aceptado como herramienta de diseño en el campo de la hidromecánica y que además cubre los dos tipos de modelos, comerciales y de investigación que se incluyen en la memoria. El código es carácter gratuito, lo que hace que su uso no esté condicionado al pago de una licencia. Dado que no se ha encontrado en la literatura ninguna aplicación de este código a la resolución de vertidos hipersalinos, el estudio presentado en este apartado constituye la primera aproximación al problema mediante este código En los apartados 4.1. y 4.2. se describen brevemente el modelo y su método de resolución, mientras que los apartados 4.3 y 4.4. muestran la aplicación del código a la reproducción de un ensayos realizado en modelo físico realizado en el marco del proyecto MEDVSA, describiendo el diseño del mallado, el modelo de turbulencia utilizado y los primeros resultados obtenidos. 4.1 Descripción del modelo OpenFOAM (OpenField Operation and Manipulation) es un paquete de herramientas y modelos para la resolución de problemas de mecánica de fluidos computacional. Es de código abierto (open source), totalmente gratuito y accesible a todo el mundo a través de su portal en internet: Funciona nativamente en entornos Linux, aunque hay proyectos paralelos para portarlo tanto a Mac como a Windows. Básicamente está formado por un gran número de librerías y códigos en C++ que permiten la resolución de problemas complejos como son la turbulencia, el flujo en una o dos fases, electromagnetismo, reacciones químicas, combustión, etc., todos ellos mediante una discretización en volúmenes finitos. Asimismo incluye aplicaciones para llevar a cabo el preproceso y postproceso de los casos simulados, incluyendo malladores muy automatizados y sofisticados (blockmesh y snappyhexmesh), herramientas para la inicialización de los campos, la descomposición de la malla en subdominios y el muestreo de variables dentro del dominio. OpenFOAM está preparado para ejecutar los casos en paralelo, permitiendo un flujo de trabajo idéntico al de los casos en serie, ya que él mismo se encarga de la descomposición del dominio y de la recomposición final (opcional). Su código en C++ y su filosofía orientada a objetos permite que tenga una estructura muy modular, lo que es claramente una ventaja a la hora de programar nuevos módulos de resolución, condiciones de contorno o modelos de turbulencia, ya que permite añadir o modificar funcionalidades sin entrar en la base del código. Esto es debido a tipos de datos y clases que son específicos en OpenFOAM, y que permiten trabajar con campos de una forma más compacta y sencilla, evitando muchas veces los bucles para realizar ciertas operaciones. En contraposición al software comercial propietario, OpenFOAM no es una caja negra, el usuario tiene control sobre todos y cada uno de los pasos del proceso de resolución modificando el código fuente. Esta es una gran ventaja ya que entre los módulos incluidos existen conversores de datos provenientes de modelos comerciales (Ansys, Fluent, CFX) hacia OpenFOAM y viceversa, lo que permite validaciones cruzadas. Aparte hay programas externos que añaden funcionalidades adicionales a OpenFOAM. El caso más importante es sin duda Paraview, programa principal utilizado para el postproceso de los casos, que tiene dos intérpretes específicos para OpenFOAM, uno de ellos para manejar los casos en serie y el otro para los casos en paralelo sin necesidad de recomponerlos. 4.2 Módulo de resolución: twoliquidmixingfoam El módulo de resolución utilizado en el modelado de los vertidos hipersalinos es twoliquidmixingfoam, incluido por defecto en OpenFOAM. Este módulo resuelve las ecuaciones de Navier-Stokes promediadas por Reynolds (RANS) en tres dimensiones mediante el esquema de volúmenes finitos para dos fases inmiscibles, asumiendo un flujo incompresible. 5
10 Cada fase tiene sus propiedades (densidad y viscosidad), y adicionalmente son necesarios la difusividad molecular y el número de Schmidt (viscosidad cinemática / difusividad molecular). Asimismo se ha de especificar el modelo de turbulencia: simulación numérica directa (DNS), k-epsilon, k-omega SST, largeeddysimulation (LES) En un primer estadio de resolución se advectan las fases mediante la resolución de la ecuación 1. d 1 U D Dt 0 dt Sc (1) Donde es la fracción de fase por celda (1 implica que la celda está llena del fluido 1, mientras que 0 indica que está llena del fluido 2, siendo una mezcla de ambos para valores intermedios); es la densidad, teniendo en cuenta ambas fases ( );U es la velocidad del fluido; D es la difusividad molecular; Sc es el número de Schmidt, anteriormente definido; y Dt es la difusividad debida a la turbulencia. La dinámica general del caso se calcula mediante las anteriormente citadas ecuaciones RANS, que incluyen conservación de la masa (ecuación 2) y conservación de la cantidad de movimiento (ecuación 3). U 0 (2) d U dt * UU U U g X p (3) eff eff Donde eff es la viscosidad dinámica compuesta por la propia viscosidad del fluido más la añadida debido al modelo de turbulencia usado; g es el vector de la fuerza de la gravedad; X es el vector posición de la celda; y gradiente de presiones dinámicas, esto es, la componente de la presión en exceso sobre la hidrostática. * p es el El paso de tiempo se calcula automáticamente para cumplir una restricción impuesta al valor máximo del número de Courant. Éste se define en la ecuación (4), y se calcula celda a celda. C r U t (4) l l es la longitud media de la celda y t es el paso de tiempo En la ecuación anterior U es la velocidad de cada celda, que empleará el programa. Puesto que el algoritmo de resolución es el PISO (Pressure Implicit with Splitting of Operators) y debido a los esquemas numéricos de discretización utilizados el número de Courant no puede pasar de 1. Los modelos de turbulencia disponibles para este módulo son numerosos: k-epsilon, k-omega, k-omega SST, k-epsilon realizable o Spalart-Allmaras, entre muchos otros. 6
11 4.3 Configuración numérica El caso que se presenta a continuación es el modelado de la descarga de salmuera a través de un tubo inclinado 60º respecto al fondo e inmerso en un tanque con agua, como se aprecia en la Figura A5-2. El caso representa uno de los ensayos en modelo físico realizado mediante técnicas ópticas no intrusivas PIV (Particle Image Velocimetry) y LIF (Laser Induced Fluorescence), realizado en el marco del proyecto MEDVSA (Anexo 3). Figura A5-2. Geometría del tanque de vertido. Las constantes físicas utilizadas son las siguientes: la viscosidad cinemática de ambos fluidos es la misma e igual a 10-6 m 2 /s. La densidad del agua es de kg/m 3, mientras que la del efluente salino es mayor, kg/m 3. La difusividad molecular se ha estimado en 10-9 m 2 /s. El número de Schmidt se ha tomado como 1, ya que al no considerar un modelo explícito de turbulencia es un parámetro irrelevante. El número de Courant, anteriormente definido se ha limitado a 0.9, de forma que el paso de tiempo sea lo más largo posible sin comprometer los resultados. La zona que limita la simulación en este aspecto es la salida del tubo de vertido, ya que es el lugar de mayor velocidad (1 m/s) y de mayor resolución. El paso de tiempo en la simulación se mantiene constante e igual a 0,0002 segundos, ya que la velocidad se mantiene constante durante toda la simulación. Las ejecuciones de OpenFOAM se realizan sobre una malla que no varía a lo largo de la ejecución de forma automática. La malla es el elemento clave a la hora de obtener unos buenos resultados, de modo que ha de tener la resolución suficiente para ser capaz de representar los procesos más relevantes, y sin embargo tener un tamaño moderado, acorde con las posibilidades computacionales del equipo informático que ejecute el modelo. En este aspecto se recomienda un número mínimo de celdas de entre 5 y 10 para modelar correctamente la salida del flujo en la boquilla. Dicha resolución es demasiado fina como para mantenerla en el resto del dominio, haciendo inviable la simulación, ya que el número total de celdas sería demasiado grande. Adicionalmente, la disposición de las celdas en relación con la dirección del chorro es muy relevante, ya que la advección de la concentración de salmuera se produce óptimamente de forma perpendicular a las caras y en otro caso aparece una disipación numérica fruto de los esquemas de discretización usados, que puede llegar a enmascarar los resultados. Puesto que el efluente se descarga a través de un tubo con un ángulo de 60º con respecto al suelo, para permitir una correcta advección de la salmuera, se ha decidido realizar una malla orientada en la dirección del tubo, y modificar la dirección de la gravedad de forma que se mantenga su dirección relativa. Además, dado que el vertido es simétrico, solo se calcula la mitad del dominio, aplicando una condición de contorno específica en el plano de simetría. Existe la posibilidad de considerar dos mallas, una gruesa en las zonas lejanas, para los primeros momentos en los que se inicia la descarga del chorro, y una refinada en el comienzo del movimiento de descenso. De esta forma en los instantes iniciales no se resuelve para gran cantidad de celdas en la zona más alejada, donde no ha llegado la 7
12 concentración. El cambio entre una y otra se ha llevado a cabo mediante una técnica de mapeo por interpolación. Cada celda de la nueva malla localiza los valores de variables más próximos a ella e interpola esos valores, obteniendo un valor medio. La reanudación del caso es automática, y ni los contornos ni las condiciones de contorno varían. Las mallas en OpenFOAM tienen siempre como base un bloque. Este bloque es un paralelepípedo que está discretizado en sus tres dimensiones. En este caso se ha propuesto un bloque de dimensiones 1,20 m de largo, 0,35 m de ancho y 0,80 m de alto, con una discretización constante de cubos de 5 mm de lado, lo que resulta en 240 x 70 x 160 celdas (2,7 millones). Esta geometría inicial es entonces cortada con cualquier superficie que se defina, y se refinan las celdas intersectadas tantas veces como se haya definido. Adicionalmente se pueden incluir áreas específicas en las que también se lleve a cabo ese refinamiento. Dicho proceso se basa en la subdivisión de las celdas para obtener el doble de resolución cada vez, por lo que el número de celdas final a partir de una sola celda para cierto nivel de refinamiento n es de 2 3n. En la Tabla A5-1 se muestra la resolución de las celdas (cubos) en función de su nivel. Nivel Resolución 0 5 mm 1 2,5 mm 2 1,25 mm 3 0,625 mm Tabla A5-1: resolución del lado de las celdas en función del nivel de refinado. La malla original se ha cortado con 2 geometrías distintas, la primera es el suelo del tanque de vertido, que no se ha refinado (mantiene la discretización de 5 mm) dado que el análisis del campo lejano no es el objeto del presente estudio y dado que el chorro no se ve influido por el fondo hasta la zona del punto de impacto. La segunda geometría es el tubo de descarga. Éste tiene 5 mm de diámetro por lo que la resolución del nivel 0 sería insuficiente para representarlo. Es por ello que se opta por refinar hasta el nivel 3, de forma que la boquilla queda definida por 7 celdas, que es un número adecuado para permitir una correcta representación de los efectos de salida del chorro. Además del refinamiento de las celdas intersectadas por las superficies, se pueden definir zonas específicas en las que también se refinen las celdas. En el caso que nos ocupa se ha definido una zona próxima a la salida de la salmuera con un nivel de refinado 3. Esta zona es de nuevo un ortoedro de longitud (dimensión en la dirección de salida del chorro) 8 cm, anchura (en la perpendicular de la longitud, y en el mismo plano de la fuerza de la gravedad) de 4,4 cm y profundidad de 2 cm. Para que la transición del tamaño de celda no sea abrupta, se especifican otros 2 volúmenes de refinamiento, uno de ellos para el nivel 2 y otro para el nivel 1, de forma que el paso del nivel 3 cercano al chorro al nivel 0 en la zona lejana sea gradual. El refinamiento de nivel 2 incluye las celdas dentro de un ortoedro de longitud 18 cm, anchura 8,4 cm y profundidad 4 cm. El refinamiento de nivel 1 incluye las celdas dentro de un ortoedro de longitud 27 cm, anchura 12,4 cm y profundidad 6 cm. En zonas menos relevantes la gradación de niveles es de 1 celda de la malla base. Todas estas zonas se pueden ver en la siguiente figura, en la que se muestran las celdas del plano de simetría y las de un plano perpendicular a él. Las mallas resultantes son estructuradas y no conformes, lo que confiere un comportamiento más flexible y una mejor adaptabilidad a geometrías complejas. La malla gruesa tiene un total de 2,2 millones de celdas, mientras que la fina ronda los 5,75 millones. Esta última se ha ejecutado en un cluster, descomponiendo el caso en 40 procesadores (2,6 GHz), de forma que se simulaban unos 2 segundos cada día. 8
13 Las Figuras A5-3, A5-4, A5-5 y A5-6 muestran la malla inicial, correspondiente a los primeros instantes de vertido antes de que el chorro empiece a caer, desde diferentes perspectivas. Figura A5-3: celdas del plano de simetría y en perpendicular a él. Figura A5-4: celdas del plano de simetría en la orientación original de la malla base. Se distinguen los 4 niveles de refinamiento. La gravedad actúa formando 60º (en el sentido de las agujas del reloj) con respecto al borde superior, y por lo tanto de forma perpendicular al borde inferior, que corresponde a la intersección de la malla con la superficie del suelo del tanque. 9
14 Figura A5-5: detalle de la malla en la cercanía del tubo de vertido: La Figura A5-6 muestra la malla correspondiente al segundo estadio, donde la malla pasa de una zona de refinamiento 1 a refinamiento de nivel 2 y su anchura se extiende a la zona en la que se espera que descienda el chorro. Figura A5-6: malla fina en el plano de simetría. Las condiciones de contorno que se aplican dependen de la naturaleza de la superficie. En ese aspecto hay que distinguir la entrada, paredes, condiciones de contorno abiertas y planos de simetría. El plano de simetría es una condición específica que se aplica a todas las variables en la cara de simetría de la malla, y aplica que el gradiente en la dirección normal a dicha cara es cero. Para el campo de concentraciones las paredes ( ) tienen una condición de contorno de gradiente en la normal a la cara igual a cero. Matemáticamente se presenta en la ecuación (5). 10
15 n i 0 x i (5) donde n es el vector unitario normal al contorno. La cara de entrada mantiene un valor constante igual a 1 por lo que siempre introduce una concentración unitaria; y los contornos abiertos tienen una condición especial denominada inletoutlet que permite la salida de concentración y en caso de necesitar la entrada de fluido lo hace con valor 0 de concentración. Para las presiones, tanto la entrada como las paredes tienen la condición de contorno buoyantpressure, que calcula la presión de forma que la segunda derivada en la dirección perpendicular al contorno sea 0. En la condición de contorno abierta se especifica la presión total, igual a 0. En cuanto a las velocidades en la entrada se especifica la velocidad del chorro (1 m/s). Las paredes consideran una condición de contorno deslizante, ya que no hay resolución suficiente para calcular la capa límite. Finalmente los contornos abiertos calculan la velocidad en función de la presión ya calculada ( pressureinletoutletvelocity ). La síntesis se presenta en la tabla A5-2. Entrada Pared Abierto P buoyantpressure buoyantpressure Valor impuesto (0) U Valor impuesto (1 m/s) Deslizante pressureinletoutletvelocity α Valor impuesto (1) Gradiente = 0 inletoutlet Tabla A5-2: condiciones de contorno aplicadas. 11
16 4.4 Resultados Los resultados obtenidos se presentan en las Figuras A5-7, A5-8 y A5-9. Figura A5-7: campo medio de velocidades en el plano de simetría. Se aprecia que las velocidades de 1 m/s aparecen muy cerca del tubo de vertido y en unos primeros centímetros se desplazan sin apenas esparcimiento lateral. A partir de unas 6 anchuras de boquilla el chorro se empieza a difundir lateralmente, reduciendo la velocidad. Puesto que el efluente salino es más denso que el agua, va cayendo, por lo que aparecen velocidades a lo largo de toda la trayectoria. Figura A5-8: campo medio de concentraciones unitarias en el plano de simetría. 12
17 El anterior campo de velocidades es el responsable del movimiento de la concentración del efluente. La concentración parte de un valor unidad, que disminuye rápidamente por la dispersión en el chorro, como puede observarse claramente en la Figura A5-9, donde se representan concentraciones entre 0 y 0,20. Este campo es acorde al de velocidades y muestra más claramente la tendencia del chorro a caer, así como su punto de impacto con el suelo del tanque. Figura A5-9: evolución tridimensional del chorro. La dispersión del efluente es un proceso altamente tridimensional por lo que el caso se ha ejecutado en tres dimensiones. La Fig. A5-9 representa un instante del chorro tridimensional, que se ha obtenido como la isosuperficie de valor de concentración 0,04 y se ha coloreado en función de la velocidad en cada punto de la superficie. Además, se han representado concentraciones tanto en el suelo del tanque como en el plano de simetría, según la escala anteriormente presentada. 13
18 5. Metodología a seguir La metodología de trabajo se proporciona en el siguiente diagrama de la Figura A5-10. Figura A5-10: metodología de trabajo. Las principales etapas se señalan en rojo, mientras que los pasos que aparecen en amarillo son opcionales. Por ejemplo no es necesario crear una malla gruesa y una fina, sino que dependerá de las necesidades del problema a resolver. 1) Mallado. Esta es la etapa inicial y de su calidad depende el resultado final, por lo que es necesario acatar unas normas de calidad específicas para el método de resolución, que en este caso será con volúmenes finitos. Por ejemplo, la relación de aspecto de las celdas se debe cuidar de forma que la dimensión mayor entre la más pequeña nunca sea mayor de 5. Asimismo la intersección entre las aristas debe ser lo más cercana a 90º posible, evitando grandes coeficientes de asimetría en las caras derivados de ángulos mayores a 135º o menores a 45º, que pueden causar inestabilidades numéricas. Otra regla de buena praxis es evitar las celdas tetraédricas en las zonas de capa límite, ya que la advección 14
19 no es óptima, sustituyéndolas por celdas poliédricas cuando sea posible. Los malladores de OpenFOAM tienden a minimizar estos defectos, y disponen de una herramienta específica para evaluar la calidad de las mallas, denominada checkmesh. Las mallas en OpenFOAM parten de un paralelepípedo mallado en sus tres dimensiones que se crea mediante la aplicación blockmesh. Para optimizar la advección del efluente es conveniente que dicha malla base esté orientada en la dirección de salida del chorro. Para este caso se ha utilizado un ortoedro de dimensiones 1,20 x 0,35 x 0,80 m,con una resolución de celda de 5 mm. Esta base es después cortada con las superficies que van a definir los contornos del problema mediante el módulo snappyhexmesh ; estas superficies son en este caso el suelo y el tubo de entrada. Los niveles de refinamiento de cada superficie son: nivel 0 para el suelo, ya que no afecta la simulación hasta que el chorro impacta en él y nivel 3 en el tubo, que permite una adecuada resolución de 7 celdas en la boquilla. Además se especifican las zonas adicionales en las que se va a refinar la malla. Una primera cubre el mayor nivel de refinamiento en la zona inicial de salida del chorro, donde el chorro viaja en línea recta sin apenas variaciones. En la parte exterior, la gran diferencia de momento con respecto a la masa de agua receptora hace que los esfuerzos tangenciales sean grandes. Los vórtices que aparecen son de tamaño menor al de la boquilla de salida, por lo que la resolución necesaria es del orden de 5 a 10 celdas por el ancho de la boquilla (0,625 mm en este caso). Esta zona se recomienda que se extienda alrededor de 15 diámetros en la dimensión longitudinal para cubrir en gran parte el inicio de los vórtices, y alrededor de 8 diámetros en las direcciones perpendiculares (4 a cada lado). Es en esta zona en donde se producen las mayores velocidades (siempre hay cesión de cantidad de movimiento) y, por ser las celdas de menor tamaño, son las más restrictivas en el cálculo del número de Courant. Una vez que el chorro pierde cantidad de movimiento, la concentración comienza a expandirse lateralmente de forma muy rápida, formando estructuras de vórtices de mayor tamaño, que hacen necesaria una resolución tan fina como en la zona inicial. El paso natural en OpenFOAM es el doble de tamaño de celda. La nueva zona debe partir aproximadamente de la mitad de la anterior (en la dirección longitudinal) y cubrir unas distancias iguales al doble, es decir, unos 30 diámetros longitudinalmente y 8 diámetros a cada lado en las demás. Con esta resolución se consigue una correcta definición global. Siempre que sea posible, sería deseable mantenerla en todo el desarrollo del chorro, pero dado que las dimensiones de la celda son todavía pequeñas, llevaría a un número total de celdas muy alto, lo que haría necesario un nuevo nivel. Se puede de nuevo pasar a una zona de menor resolución según lo explicado anteriormente: cubriendo longitudinalmente la mitad de la anterior, y de dimensiones el doble de nuevo. En el momento en que se espere que el chorro pueda empezar a caer, es conveniente variar la geometría y empezar a mallar la zona inmediatamente aledaña hasta que se intercepte el suelo. Las mallas resultantes serán estructuradas y no conformes. Posteriormente se les aplican las condiciones de contorno como se recoge en la tabla 2. 2) Ejecuciones En este segundo paso se procede a la resolución temporal del problema mediante el programa de resolución twoliquidmixingfoam, siendo posible descomponer el caso en diferentes dominios, de forma que se paralelice la resolución y se obtengan resultados de una forma más rápida. Se ha fijado en 10-7 la tolerancia de los solvers para que los resultados se calculen en un tiempo razonable, y a la vez muestren precisión. El cálculo de la presión es especialmente importante, ya que es la variable de la que deriva el campo de velocidades a través de su gradiente. 15
20 Si se ha ejecutado una simulación inicial en malla gruesa, pudiéndose interpolar la solución a una malla más fina para una mayor resolución en las zonas de interés. 3) Análisis de datos. Finalmente se procede a analizar los resultados mediante el programa de postprocesado por defecto paraview u otro cualquiera, comercial o propio. El análisis permite acceder a los campos instantáneos de presión, velocidad y concentración unitaria y obtener mediante el uso de filtros todos los resultados tridimensionales que se deseen. Se pueden promediar en el tiempo los campos para obtener campos medios que representan el comportamiento global del chorro. 6. Futuras líneas de investigación. El presente documento supone un primer avance en la implementación de modelos CFDs para la simulación del campo cercano de vertidos hipersalinos, estableciendo criterios para el diseño del mallado y para la elección del modelo de turbulencia, incluyendo los primeros resultados obtenidos para el modelado de un ensayo experimental. En el futuro próximo se seguirá avanzando en esta línea de investigación, llevando a cabo las siguientes tareas: - Profundizar en la optimización de los parámetros de modelado y en el diseño de la malla para reducir el tiempo computacional necesario para el modelado. - Estudiar la idoneidad de los diferentes modelos de turbulencia existentes en la literatura para el estudio de chorros hipersalinos para una correcta caracterización de los procesos de mezcla de masa y cantidad de movimiento. - Validar los resultados numéricos obtenidos con OpenFoam con datos experimentales obtenidos, en el marco del proyecto MEDVSA, con técnicas experimentales PIF y LIF, comparando las magnitudes medias y turbulentas y el valor de cada término en las ecuaciones de gobierno. - Ampliar el estudio de modelado con OpenFoam al caso de tramo difusor con chorros múltiples que interaccionan entre sí, y a otras configuraciones de descarga de interés en los vertidos de salmuera. - Ampliar el dominio del modelado del vertido en chorro a la región completa de campo cercano, incluida la capa de esparcimiento lateral ( spreading layer ), e incluso a la región de campo lejano. - Adaptar el código OpenFoam para simular medio receptores dinámicos, con presencia de corrientes y oleaje. 16
ANEXO 4 CARACTERIZACIÓN DEL CLIMA MARINO PARA EL VERTIDO DE SALMUERA V. 1.1
ANEXO 4 CARACTERIZACIÓN DEL CLIMA MARINO PARA EL VERTIDO DE SALMUERA V. 1.1 Este documento se ha obtenido como resultado del proyecto de I+D+i: MEDVSA. Desarrollo e implementación de una metodología para

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución