Cormix2 PDF
Cormix2 PDF
Cormix2 PDF
Ficha 3.2.
SUBSISTEMA CORMIX 2
V. 1.1
98
”
- GUÍA METODOLÓGICA
Autores:
Figura 1: Planta (izq) y perfil (dcha) del vertido mediante chorros múltiples con CORMIX2.
100
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
101
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
102
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
M 01 / 2
receptor. Lm . Para distancias desde el punto de vertido mucho mayores que Lm, el flujo
UA
depende de la corriente ambiental y el efecto de la cantidad de movimiento es despreciable.
Figura 2A: Escalas de longitud de CORMIX2, para el caso de NO INTERACCIN de chorros. Medio receptor de
densidad uniforme. Fuente: Manual de CORMIX2 (1)
103
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Figura 2B: Escalas de longitud de CORMIX2, para el caso de NO INTERACCIÓN entre chorros. Medio
receptor estratificado. Fuente: Manual de CORMIX (4)
l m : Escala cantidad de movimiento-corriente del medio receptor para flujo en ranura. Escala
de longitud representativa de la distancia a partir de la cual el chorro bidimensional en ranura pasa de
estar dominado por la cantidad de movimiento transmitida en la descarga a estar dominado por la
mo
corriente en el medio receptor. l m Para distancias desde el punto de vertido mucho mayores
U A2
que Lm, el flujo depende de la corriente en el medio receptor y el efecto de la cantidad de movimiento
es despreciable.
104
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
UA
la : Escala corriente en el medio receptor-estratificación. Escala de longitud
1/ 2
representativa de la distancia a partir de la cual la pluma pasa de estar dominada por la corriente en
el medio receptor, a estar dominada por la estratificación del mismo.
Figura 4: Escalas de longitud para medio receptor estratificado, de CORMIX2. Fuente: Manual CORMIX2 (1).
105
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
106
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
107
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
108
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
109
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
- MNU8: Vertido en una zona de calado limitado, de modo que los chorros interaccionan con
la superficie. El medio receptor es homogéneo y la corriente ambiental es débil y perpendicular a los
chorros. El flujo es inestable.
Para cada una de las clases de flujo de CORMIX2, el subsistema lleva asociados módulos de
cálculo para caracterizar el comportamiento del vertido, teniendo en cuentas las características del
efluente, los parámetros de descarga y las condiciones en el medio receptor. En total, CORMIX2
incluye módulos de cálculo, clasificados en los siguientes siete grupos:
o Módulos para el cálculo del campo cercano en flujos de chorros múltiples convectivos.
o Módulos para el cálculo de los procesos de interacción con los contornos.
o Módulos para el cálculo de tramos difusores inestables: flujos de mezcla en campo
cercano.
o Módulos para el cálculo de los procesos de interacción con los contornos en difusores
estables.
o Módulos para el cálculo del flujo de campo intermedio en tramos difusores inestables.
o Módulos de cálculo para los procesos de esparcimiento por flotabilidad residual.
o Módulos para el cálculo de los procesos de difusión pasiva.
La descripción detallada sobre la base teórica del modelo, la formulación de cada módulo de
cálculo, las opciones disponibles y el manejo del programa se pueden encontrar en el documento
base original de CORMIX2, así como en el manual de usuario, disponibles en la página web de la
EPA (“Environmental Protection Agency”) y en la página web de CORMIX respectivamente. Al final de
la ficha se indica la bibliografía de mayor interés en relación con el subsistema CORMIX2.
110
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
3. HIPÓTESIS SIMPLIFICATIVAS.
A continuación se describen las principales hipótesis simplificativas del subsistema CORMIX2
en la simulación del fenómeno de interacción entre chorros, y el modelado del vertido bajo distintas
configuraciones del dispositivo de vertido (tramo difusor):
Figura 10: Interacción entre chorros y paso de chorro 3D a pluma 2D, vista lateral
(izquierda) y en planta (derecha). Difusor de chorros unilaterales. Fuente: Manual
CORMIX (3)
111
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
La figura 11 muestra el proceso de interacción entre chorros en el caso de un difusor con dos
boquillas por elevador y chorros bilaterales o alternos, que descargan en direcciones opuestas (en
este caso formando 180º entre sí):
Figura 11: Interacción de chorros y paso de chorro 3D a pluma 2D, vista lateral (izquierda) y
en planta (derecha). Difusor de chorros bilaterales o alternos. Fuente: Manual CORMIX2 (1)
112
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Figura 13: Angulos de configuración del tramo difusor en CORMIX2. Vista lateral
(arriba) y vista en planta (abajo)
Simplificaciones en los dispositivos de vertido con una boquilla por elevador (“one
nozzle per riser”):
Para el caso de una boquilla por elevador, CORMIX2, distingue entre chorros unilaterales
(“unidirectional”), situado en el mismo lado del difusor, y chorros bilaterales o alternos (“alternating”)
situados a ambos lados del difusor.
Para el caso de chorros unilaterales, CORMIX2 requiere que el usuario introduzca los cuatro
ángulos de configuración del vertido: , , , . Aunque además distingue entre chorros de la misma
dirección (“same direction”) y en abanico (“fanned out”), aplica exactamente las mismas fórmulas en
ambos casos, de modo que realmente no existe esta distinción en el código del modelo.
Para el caso de chorros bilaterales, CORMIX2 no ofrece la opción de introducir ningún ángulo
porque asume directamente que el vertido se realiza mediante una boquilla con ángulo de inclinación
del chorro vertical 90º , perpendicular al difusor: 90º , aplicando la hipótesis de vertido
mediante ranura equivalente. Al igual que en el caso anterior, aunque el modelo permita seleccionar
entre misma dirección y en abanico, el código no lo distingue (y de hecho es una opción absurda
porque en realidad el vertido lo considera vertical).
113
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Simplificaciones en los dispositivos de vertido con dos boquillas por elevador (“two
nozzles per riser”):
En el caso de tramo difusor con dos boquillas por elevador, chorros bilaterales (alternos) y
separados un ángulo horizontal de 180º (“Two nozzles per riser: about 180 degrees”), CORMIX2 no
requiere definir los ángulos: , , , , porque asume directamente un dispositivo con una única
boquilla por elevador, vertiendo el caudal equivalente a las dos boquillas mediante un chorro de
inclinación vertical ( 90º ) y perpendicular al tramo difusor ( 90º ). Bajo esta drástica
simplificación, aplica directamente la hipótesis de vertidos mediante ranura equivalente. En el caso de
efluentes de flotabilidad positiva, la aproximación podría ser válida para ciertos tipos de flujos, pero
nunca en el caso de efluentes hiperdensos, para los que CORMIX2 llega a una clasificación de
recirculación con inestabilidad, aunque sea erróneo.
Esta simplificación se observa directamente en la interfaz y en la hoja de resultados del
modelo. Si se introduce un dispositivo de vertido de dos boquillas por elevador formando 180º entre sí
(“Two nozzles per riser: about 180 degrees”), CORMIX2 no permite introducir ningún valor del ángulo
de inclinación del chorro: . Tras ejecutar el modelo, puede observarse en la hoja de resultados que
el modelo ha tomado por defecto: 90º y 90º y que ha asumido la hipótesis de vertido
mediante ranura equivalente. De hecho, si ejecutamos el mismo caso, para un vertido con una única
boquilla por elevador y la misma dirección, y una inclinación de 90º , los resultados serán
exactamente iguales. La figura 14 muestra la simplificación de CORMIX2 en este caso:
Lo mismo sucede al simular un tramo difusor con dos boquillas por elevador, de chorros
unilaterales y separados entre sí por un ángulo horizontal mayor a 60º ,en cuyo caso CORMIX2
ofrece únicamente la opción de considerar un ángulo de separación de 180º, aplicando entonces la
hipótesis de vertido mediante ranura equivalente.
Cuando se trata de un tramo difusor con dos boquillas por elevador, de chorros unilaterales y
separados entre sí por un ángulo horizontal inferior a 60º (“Two nozzles per riser: less than 60
degrees”), el programa solicita los cuatro ángulos de configuración: , , , ; como si fuera una
única boquilla por elevador. Esto es así porque en realidad CORMIX2 asume que es una única
boquilla, vertiendo el caudal equivalente a las dos boquillas, con el ángulo de chorro respecto al
difusor que el usuario haya introducido: ( ). La figura 15 muestra un esquema en planta de esta
simplificación:
114
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Figura 15: Simplificación de CORMIX2 para dispositivos de vertido de dos boquillas por
elevador y chorros unilaterales formando un ángulo menor a 60º entre sí. Izquierda:
dispositivo real. Derecha: simplificación de CORMIX2. Vistas en planta
La Tabla 1 agrupa las distintas configuraciones que teóricamente ofrece CORMIX2, en las
dos configuraciones que realmente, a efectos de cálculo, es posible simular:
115
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Tramo difusor de varias boquillas por elevador, flujo neto de cantidad de movimiento
horizontal distinto de cero, chorros la misma dirección o en abanico y 75º
(“several nozzles: net horizontal momentum flux Non-zero: same direction / fanned out” y
75º ). (G.8)
Tramo difusor de una boquilla por elevador, chorros unilaterales, vertiendo en la misma
dirección o en abanico y 75º
(“single nozzle: unidirectional: same direction / fanned out ”). (Casos G.1 y G.2)
Tramo difusor de una boquilla por elevador, chorros bilaterales (alternos) CONFIG 2
(“single nozzle: alternating”). (Caso G.3)
Asume chorro
Tramo difusor de dos boquillas por elevador, separadas entre sí por un ángulo horizontal único:
inferior a 60º , chorros vertiendo en la misma dirección o en abanico y 75º 90º
(“two nozzles: less than 60 degrees: same direction / fanned out” y 75º ). (Casos G.4 y
G.5)
Tramo difusor de dos boquillas por elevador, formando un ángulo de unos 180º entre sí. y aplica
(“two nozzles: about 180 degrees”). (Caso G.6)
Tramo difusor de varias boquillas por elevador, flujo neto de cantidad de movimiento Hipótesis de
horizontal cero, chorros vertiendo en la misma dirección o en abanico. vertido
(“several nozzles: net horizontal momentum flux zero: same direction / fanned out”). (Caso mediante
G.7) ranura
equivalente
Tramo difusor de varias boquillas por elevador, con flujo neto de cantidad de movimiento
horizontal distinta de cero, chorros vertiendo en la misma dirección o en abanico y 75º
(“several nozzles: net horizontal momentum flux Non-zero: same direction/fanned out” y
75º ). (G.8)
Tabla 1. Casos a los que se simplifican las configuraciones de vertido de CORMIX2.
NOTA: Los casos a los que se refiere el cuadro anterior son los modelados en el análisis de sensibilidad, ver
Anexo 3 de la presente ficha.
116
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
117
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
118
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
La tabla 3 indica el rango de valores admisible por CORMIX2 para los datos de entrada,
según los autores del modelo. Cuando el valor introducido por el usuario se encuentra fuera de este
rango, el programa avisa y no ejecuta:
Profundidad media en en la
0 – 1000m Densidad del efluente 800 – 1400 Kg/m³
columna de agua
119
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
- La clasificación de CORMIX2 incluye múltiples clases de flujo distintos, en los que aplica módulos de cálculo
diferentes, cubriendo, en teoría, un amplio abanico de casos.
- Teóricamente modela el comportamiento de un vertido mediante chorro individual tanto en un medio receptor
ilimitado como confinado, simulando en este caso la interacción del flujo con los contornos (laterales, fondo,
superficie). Sin embargo es necesario tener en cuenta las limitaciones, hipótesis y errores del CORMIX2 en la
simulación de este tipo de fenómenos.
- Considera efectos COANDA de adhesión del efluente al fondo.
- En el caso en que CORMIX2 detecta que el flujo no interacciona con los contornos, aplica el módulo de cálculo
CORJET, ofreciendo en los resultados una descripción numérica detallada de la evolución del flujo.
- Distintas posibilidades respecto al diseño del dispositivo de vertido: configuración del tramo difusor, número de
boquillas por elevador, separación, disposición de los chorros, altura, diámetro, ángulo de inclinación del chorro, etc.
- Simula el fenómeno de interacción entre chorros (“merging”) durante su trayectoria y tras el impacto con el fondo.
- El modelo aproxima el comportamiento del efluente tanto en la región de campo cercano como lejano, mediante el
acoplamiento de distintos módulos de cálculo, incluida la capa turbulenta convectiva tras el impacto con el fondo.
120
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
La Figura 17 muestra las configuraciones teóricas del subsistema CORMIX2, aunque según
lo explicado en la sección 3, las simplifica distinguiendo únicamente entre dos tipos de configuración:
UNIDIRECTIONAL DIFFUSSER STAGGED DIFFUSER ALTERNATING DIFFUSER
121
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Densidad
efluente y del A mayor diferencia de densidad entre el efluente y el medio receptor, ligeramente
Medio receptor Pequeña menor es la altura máxima y el alcance del chorro en el punto de impacto, así
como menor es la dilución.
A , o
122
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Altura boquilla
sobre el fondo A mayor altura de la boquilla, ligeramente mayor es también la altura máxima
Muy pequeña
alcanzada por el chorro, su alcance hasta el punto de impacto y la dilución.
ho
123
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Para el caso de difusores con una boquilla por elevador, para las mismas
características de chorro, si se modifican la longitud del difusor y el número de
boquillas, manteniendo la misma separación, el número de boquillas no influye en
Número de
el comportamiento del chorro.
boquillas nula
N
En el caso de elevadores con varias boquillas, los resultados son muy similares,
disminuyendo muy ligeramente la dilución con el aumento del número de
boquillas.
Tabla 5: Conclusiones del análisis de sensibilidad del subsistema CORMIX2 frente a variaciones en los datos de
entrada
124
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Figura 18: Resultados de las validación del subsistema CORMIX2 para vertidos de efluentes de
flotabilidad negativa. Chorros vertidos verticalmente, en un medio receptor homogéneo.
125
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
disponible el caso descrito en esta sección 6, correspondiente a un flujo de la clase MNU2, e incluye
únicamente geometría del chorro pero no valores de dilución.
Respecto a la validación de las otras clases de flujo de CORMIX2 para efluentes de
flotabilidad negativa, los autores no presentan ninguna validación de resultados: ni para las variables
ambientales (medio receptor en reposo, corriente con otro ángulo respecto al chorro, estratificación
en la columna de agua, etc.), ni para los distintos parámetros de descarga (inclinaciones del chorro,
altura de la boquilla, configuraciones de vertido, separación entre boquillas, etc.) Por otra parte, a
pesar de que CORMIX2 plantea como uno de sus valores añadidos la opción de simular la interacción
del flujo con los contornos, los autores no han presentado ningún estudio de validación de este
supuesto, así como tampoco para el caso de chorros que interaccionan entre sí, para los que el
subsistema aplica la hipótesis de vertido mediante ranura equivalente.
En vista de lo anterior, se puede concluir que los estudios de validación del subsistema
CORMIX2, para vertidos de efluentes de flotabilidad negativa, son muy escasos y pobres (en realidad
es solamente uno y no incluye valores de dilución), lo que resta fiabilidad al modelado de vertidos
hiperdensos con este subsistema. Esto es especialmente relevante en los casos complejos
(interacción con los contornos, flujos inestables ocupando toda la columna de agua, etc.), que
teóricamente CORMIX2 es capaz de simular, pero que en realidad no se encuentran validados.
Con el fin de remediar este vacío y de valorar la fiabilidad de las simulaciones de CORMIX2 y
el error estimado de sus resultados, en la simulación de vertidos de chorros hiperdensos, se ha
realizado un nuevo estudio de validación del modelo aparte y más completo que el presentado por los
autores del modelo. Para ello se han utilizado los resultados experimentales obtenidos por distintos
autores en ensayos en modelo físico, anteponiendo aquellos autores con mayor experiencia en este
tipo de experimentos, y aquellos que han utilizado técnicas láser de medición (LIF: Fluorescencia
Láser Inducida), dado que permite obtener resultados más precisos del comportamiento del efluente.
En la sección 7 de esta ficha se muestran los resultados y las principales conclusiones
derivadas de este estudio de validación.
126
MODELOS COMERCIALES: CORMIX2 de CORMIX
No existen datos publicados en relación con vertidos de efluentes hipersalinos mediante tramo
difusor de chorros múltiples.
127
MODELOS COMERCIALES: CORMIX2 de CORMIX
La profundidad que el subsistema toma para realizar los cálculos es “discharge detpth”, y
no “average depth”, por lo que debe introducirse el calado real en dicha casilla.
S
S
128
MODELOS COMERCIALES: CORMIX2 de CORMIX
En el Anexo 2 de esta ficha se presenta una tabla de valores recomendados para los
distintos parámetros de entrada al modelo, para el diseño del vertido al mar de salmuera,
particularizando para el Mar Mediterráneo.
Teniendo en cuenta el análisis de sensibilidad del modelo ante los distintos
parámetros de entrada (cuyas conclusiones se exponen en la sección 6 y cuyos resultados
numéricos se incluyen en el Anexo 3 de esta ficha) se recomienda:
- Evitar modelar con CORMIX2 el caso de vertidos en chorro con inclinaciones superiores a
75º, dado que iguala los resultados a una inclinación vertical.
- El modelo CORMIX2 es poco sensible a la densidad del efluente y del medio receptor (en el
rango de valores realistas para salmuera de osmosis inversa vertida en el mar
Mediterráneo), Sin embargo, la concentración salina (y con ello la densidad) del efluente y
del medio receptor sí son determinantes en el cálculo del grado de dilución necesario para
garantizar que no se superan los límites críticos de salinidad establecidos para el medio
receptor. Esta sensibilidad se refleja en la fórmula para calcular la dilución que es
necesario alcanzar, según el valor límite establecido para la concentración salina:
V V A C o C A C o g´o o A .
S o
Vo C CA C g´ lim ite A
Tabla 6: Recomendaciones respecto al modelado de vertidos de salmuera mediante chorros múltiples en tramo
difusor con el subsistema CORMIX2
Recordamos que el eje “x” en CORMIX coincide siempre con la dirección de la corriente en el
medio receptor y con la dirección de la línea de costa.
- Si el tramo difusor de chorros es perpendicular a la costa (y a la corriente) y contiguo al
emisario, los valores de distancia a introducir serán:
“Distance to one/other endpoint”: YB2=YB1+L, siendo “L” longitud del difusor: “length diffuser” y
0º
En el caso de chorros perpendiculares al tramo difusor: 90º , de modo que 90º
____________________________________________________________________
129
MODELOS COMERCIALES: CORMIX2 de CORMIX
Tabla 7: Recomendaciones respecto al uso y manejo del subsistema CORMIX2 para el modelado de vertidos de
salmuera mediante chorros múltiples.
130
MODELOS COMERCIALES: CORMIX2 de CORMIX
Uso académico - - - - -
- (2) Jirka, G.H. and P.J. Akar (1991). "Hydrodynamic Classification of Submerged Multiport
Diffusers Discharges". Journal of Hydraulic Engineering, ASCE, 117, No.9, 1113-1128.
- (3) Doneker, R.L; Jirka, G.H. (2001) “CORMIX-GI systems for mixing zone analysis of brine
wastewater disposal”. Desalination, volume 139, pp. 263–274.
- (4) Roberts, P.J.W; Tian, X. (2004). “New experimental techniques for validation of marine
discharges models”. Environmental Modelling & Software 19.
- (5)Tong, S.S., Stolzenbach, K.D., (1979). ”Submerged discharges of dense effluent”. Report
No. 243. (R79-35). Ralph M. Parsons Laboratory for Water Resources and Hydrodynamics.
Dept. of Civil Engineering. Mass. Inst. of Tech.
- (6) Zeitoun, M.A et al. (1970). “Conceptual designs of outfall systems for desalination plants”
.Research and Development Progress Rept. No 550. Office of Saline Water, U.S. Dept, of
Interior.
- (7) Roberts, P.J.W; Fellow; Ferrier, A; Daviero, G. (1997). “Mixing in inclined dense jets”.
Journal of Hydraulic Engineering, vol 123, No 8.
- (8) Autores varios, (2003): CEDEX, Aguas de la Cuenca del Segura. S.A., Universidades de
Alicante y Murcia, Instituto Oceanográfico de Murcia, y Centro de Estudios Avanzados de
Blanes-CSIC,). “Estudio de los efectos de incrementos de salinidad sobre la fanerógama
131
MODELOS COMERCIALES: CORMIX2 de CORMIX
marina Posidonia oceanica y su ecosistema, con el fin de prever y minimizar los impactos
que pudieran causar los vertidos de aguas de rechazo de plantas desaladoras”. Documento
de síntesis.
- Iso, S; Suizu, S; Maejima, A. (1994). “The Lethal Effect of Hypertonic Solutions and
Avoidance of MarineOrganisms in relation to discharged brine from a Desalination Plant”.
DESALINATION 97, pp389-399. ELSEVIER.
- Lloret, M.P; Sánchez, J.L. (2001). “Bioensayo para la evaluación de los efectos producidos
por un incremento de salinidad sobre la supervivencia del erizo de mar Paracentrotus
lividus”. II Congreso Nacional de Desalación. Alicante. Asociación Española de Desalación y
Reutilización Alicante.
- Terrados, J (1991). “Crecimiento y producción de las praderas de macrófitos del Mar Menor,
Murcia”. Tesis Doctoral. Universidad de Murcia.
- M.M. Mendez-Diaz and G.H. Jirka. (1996). "Buoyant Plumes from Multiport Diffusers in
Deep Water", J. Hydraulic Engineering, ASCE, 122, No.8.
- Jirka, G.H.; Doneker, R.L.; Hinton, S.W. (1996). “MANUAL FOR CORMIX: A Hydrodynamic
Mixing Zone Model and Decision Support System for Pollutant Discharges into Surface
Waters”. U.S.EPA, Office os Science and Technology, Washington.
- Doneker, R.L.; Jirka, G.H. (2007). “CORMIX USER MANUAL: A Hydrodynamic Mixing Zone
Model and Decision Support System for Pollutant Discharges into Surface Waters”.
Disponible en http://www.mixzon.com/downloads/
132
MODELOS COMERCIALES: CORMIX2 de CORMIX
ANEXOS
ANEXO 1.
NOTA: En CORMIX la densidad debe introducirla directamente el usuario, que habrá de calcularla previamente a
partir de los datos de temperatura y salinidad (fluido incompresible). En CORMIX el eje “x” coincide siempre con
la dirección de la corriente en el medio receptor y con la dirección de la costa.
- Parámetros de descarga de la salmuera: ver “Recomendaciones de diseño del vertido de
salmuera” y características del efluente (considerando agua de alimentación y tecnología de
desalinización). En este caso:
133
MODELOS COMERCIALES: CORMIX2 de CORMIX
134
MODELOS COMERCIALES: CORMIX2 de CORMIX
135
MODELOS COMERCIALES: CORMIX2 de CORMIX
136
MODELOS COMERCIALES: CORMIX2 de CORMIX
137
MODELOS COMERCIALES: CORMIX2 de CORMIX
CASILLA “CORSPY”
Exporta los datos introducidos en el proyecto y permite visualizar gráficamente la
configuración del dispositivo vertido elegida. La configuración diseñada tiene el siguiente aspecto
representada mediante CORSPY, se representa solamente un tramo:
138
MODELOS COMERCIALES: CORMIX2 de CORMIX
CASILLA “3DCORVUE”.
Exporta los resultados de la simulación del comportamiento del vertido y muestra
gráficamente en 3D, la geometría del chorro y pluma, coloreando en distintos tonos según la
concentración de salmuera. El efluente tiene el siguiente aspecto al ser representado por
3DCORVUE:
139
MODELOS COMERCIALES: CORMIX2 de CORMIX
ANEXO 2.
“EFFLUENT”
(CARACTERISTICAS DE LA SALMUERA)
Concentración del efluente Velocidad de salida del
Tipo contaminante Densidad del efluente
respecto M. Receptor chorro
(Pollutant type) (Effluent density)
(Discharge concentration) (Effluent velocity)
Co 30.5 psu 3-5m/s-5m/s 1052-1047Kg/m³
conservative
(*1) (*2) (*3)
“AMBIENT”
(CONDICIONES EN EL MEDIO RECEPTOR)
Profundidad media
de columna de agua Densidad del medio
Velocidad viento Velocidad corriente Coef. rozamiento
en vertido (calado) receptor
(Wind speed) (Ambient velocity) (Manning/Darcy)
(Average water (Ambient density)
depth/Depth dischar)
“DISCHARGE”
(PARAMETROS DE LA DESCARGA). CORMIX2
Longitud Distancia a Diámetro/Area Coef. Contrac Nº de Angulo entre
Altura de la
difusor costa de la boquilla diámetro boquillas tramo difusor y
boquilla
corriente
(Diffuser (Distance to (Port Diamet, (Contraction (Total number
(Port height)
length) ends points) area) ratio) of openings) (Gamma)
1-1.5m >0.15m
L YB1/YB2 1 N 0-180º
(*8) (*9)
140
MODELOS COMERCIALES: CORMIX2 de CORMIX
(Two nozzles per riser) (Same direction/ Fanned out) (Same direction/ Fanned out)
NOTAS:
(*1) La concentración salina del efluente salmuera respecto a la del medio receptor (Mar
Mediterráneo) depende de la tasa de conversión: “R” en el proceso de desalinización en planta. Para
plantas de osmosis inversa, que son las más frecuentes en nuestras costas mediterráneas, la tasa es
del orden del R=40%-50%. Para calcular la concentración salina del efluente se aplica la fórmula:
C
Co ( A ) ; siendo R: tasa de conversión de la planta; y : salinidad del agua de alimentación de la
1 R
planta (en este caso, agua del Mar Mediterráneo: 37-38psu).
(*2) Cuanto mayor es la velocidad de salida del chorro, mayor será el número de Froude
densimétrico y la dilución de la salmuera en campo cercano. Se recomienda diseñar con valores de
3.5-5m/s y números de Froude densimétrico entre 20-25.
141
MODELOS COMERCIALES: CORMIX2 de CORMIX
(*4). El calado en la zona de descarga será el que corresponda. El valor que toma el
subsistema para los cálculos es el correspondiente a “Discharge depth”. Siempre que sea posible, se
recomienda localizar el vertido a una profundidad suficiente para evitar el impacto del chorro con la
superficie del mar, dado que este hecho reduce la dilución y es difícil de simular con los modelos
comerciales disponibles. Hay que tener en cuenta que en el caso de flujo que interacciona con la
superficie, la fiabilidad de los resultados del subsistema CORMIX2 no está demostrada.
(*5) Para los valores característicos de viento en el Mar Mediterráneo, considerando las
profundidades de vertido, y el hecho de que la salmuera se hunde hacia el fondo y se propaga como
una corriente de gravedad sobre el lecho, la velocidad del viento no es un factor determinante. Desde
un punto de vista conservativo, se recomienda dejar el valor “0”.
(*6) Las corrientes en el Mar Mediterráneo no superan valores de 1m/s, siendo generalmente
mucho menores e incluso prácticamente nulas durante algunos meses del año (calma chicha). La
simulación debe incluir los distintos escenarios mensuales o estivales y considerar el valor y dirección
de la corriente en el medio receptor en las profundidades en las que se desplaza la salmuera.
(*7) Con respecto al coeficiente de rozamiento con el fondo, dependerá del tipo de fondo:
arenoso, rocoso, etc.
(*8) Se recomienda una altura de la boquilla respecto al fondo de 1m, para evitar la
realimentación del chorro con salmuera. En el caso de vertidos con limitación de calado, se
recomienda diseñar al menos con una altura de 0,5 m.
(*9) En función del caudal de salmuera a verter se jugará con el número de boquillas y el
diámetro de las mismas, de modo que la velocidad de salida de cada chorro de 3.5-5m/s (para
142
MODELOS COMERCIALES: CORMIX2 de CORMIX
garantizar una dilución significativa). Se recomiendan diámetros de boquilla superiores a 15cm para
minimizar el riesgo de obstrucción ante el fouling.
(*11) Desde un punto de vista constructivo es más recomendable, además de ser el caso de
difusor perpendicular al tramo difusor, el único validado con datos experimentales.
(*12) Los estudios experimentales llevados a cabo demuestran que si el chorro se vierte en la
misma dirección y sentido que la corriente predominante en el medio receptor, o bien, perpendicular a
la misma, se favorece la dilución. Si el vertido se realiza en la misma dirección pero sentido contrario
a la corriente, se producen fenómenos de intrusión y la dilución se reduce (7).
(*13) Depende del exceso de salinidad de los ecosistemas a proteger en cada caso (8).
143
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
ANEXO 3.
- (sigma): ángulo horizontal entre el chorro y la corriente ambiental (en el medio receptor)
- (gamma): ángulo horizontal del tramo difusor con la corriente en el medio receptor.
144
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
- N: número de boquillas.
- YB1: distancia desde el extremo más cercano del tramo difusor a la costa.
- YB2: distancia desde el extremo mas alejado del tramo difusor hasta la costa.
- Z: coordenada vertical en el sistema de referencia cartesiano.
- X: coordenada horizontal, que corresponde a la dirección y sentido de la corriente en el medio receptor, en el sistema de referencia cartesiano.
- Y: coordenada horizontal, que corresponde a la dirección y sentido perpendicular a la corriente en el medio receptor, en el sistema de referencia cartesiano.
- Zmax: Valor de la coordenada vertical (z) del eje del chorro, correspondiente al punto de máxima altura alcanzada por el eje.
- S (Zmax): Valor de la dilución en el eje del chorro, para el punto de máxima altura del eje del chorro.
- Xi: Valor de la coordenada horizontal del eje del chorro, en el punto de impacto del eje del chorro con el fondo (alcance horizontal del chorro)
- Seje (Xi): Valor de la dilución en el eje del chorro, para el punto de impacto del eje del chorro con el fondo.
- XL: Valor de la coordenada horizontal del eje del chorro en el final de la zona de campo cercano
- Seje (XL): Valor de la dilución en el eje del flujo, para el final de la zona de campo cercano.
- B: semiancho del chorro, distancia radial donde la velocidad es igual al 37% de la velocidad en el eje del chorro.
- BH: semiancho horizontal de la pluma.
- BV: semi espesor de la pluma
Caso del análisis de sensibilidad:
Para el análisis de sensibilidad se considera la salmuera derivada de un proceso de desalinización mediante osmosis inversa con tasa de conversión del
45%. El efluente se vierte al Mar Mediterráneo mediante chorros múltiples, con tramo difusor lineal.
El sistema de referencia es cartesiano, con origen la boquilla de salida del chorro. La corriente en el medio receptor (y la línea de costa) coincide en
dirección y sentido con el eje “x”, de modo que si el chorro se vierte en la misma dirección que la corriente ( 0º ) seguirá la dirección y sentido del eje “x”,
mientras que si se vierte perpendicularmente a la corriente ( 90º ) seguirá inicialmente la dirección del eje “y”, pero luego será desviado por la corriente en
dirección “x”. Para un chorro vertido en la misma dirección pero sentido opuesto a la corriente ambiental ( 180º ), se desplazará inicialmente en sentido opuesto
a la corriente (“x” negativo), hasta ser desviada en la dirección de la misma (“x” positivo).
El rango de valores de las variables corresponde a parámetros de diseño característicos de vertido de salmuera y de escenarios ambientales en el Mar
Mediterráneo.
145
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
75º 23 / 215 MNU1 X=1.5; BH=75; S=9.7 Asume vertido vertical e hipótesis ranura
90º 23 / 215 MNU1 X=1.5; BH=75; S=9.7 Asume vertido vertical e hipótesis ranura
146
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
B) INFLUENCIA DE LA DENSIDAD DEL EFLUENTE VERTIDO (ρo) Y DE LA DENSIDAD DEL FLUIDO EN EL MEDIO RECEPTOR (ρa).
U A 0 ; H A =20m ; 0º (co-flowing), ho 1 ; 60º ; U o 5m/s ; D 0.2 m; U viento 0 ; RMZ=100m; L=150m;YB1=1000; YB2=1150; N=11;s=15;
Parámetros fijos:
90º ; 90º
B.1) Para A = 1026Kg/m³, variamos la densidad del efluente: o , desde 1047Kg/m³ a 1053 Kg/m³ (con incrementos de 2kg/m³).
1049 23.8 / 222 MNU1 7.6 6.6 10.5 15.9 NO No ejecuta (Ua=0)
B.2) Para A = 1027.5Kg/m³, variamos la densidad del efluente: o , desde 1047Kg/m³ a 1053 Kg/m³ (con incrementos de 2kg/m³).
1047 25.9 / 241.5 MNU1 8.1 7.1 11.3 17.1 NO No ejecuta (Ua=0)
1049 24.7 / 230 MNU1 7.8 6.8 10.8 16.3 NO No ejecuta (Ua=0)
1051 23.6 / 220 MNU1 7.5 6.6. 10.4 15.7 NO No ejecuta (Ua=0)
147
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
B.3.) Para A = 1028Kg/m³, variamos la densidad del efluente: o , desde 1047Kg/m³ a 1053 Kg/m³ (con incrementos de 2kg/m³).
148
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
15 1051.2 1027.8 24 /220 MNU1 7.5 6.6 10.4 15.7 NO No ejecuta (Ua=0)
18 1050.2 1027 23.8 / 222 MNU1 7.6 6.6 10.5 15.9 NO No ejecuta (Ua=0)
21 1049.3 1026 23.7 / 220 MNU1 7.5 6.6 10.4 15.8 NO No ejecuta (Ua=0)
C.2.) Representando una situación primaveral/estival en el Mar Mediterráneo, con presencia de termoclina a profundidades entre 15 y 20m. Consideramos
que la picnoclina está situada a una profundidad de 12m, en un calado total de 20m, existiendo un saldo término entre la superficie y el fondo. Consideramos que el
agua de alimentación se extrae de profundidades inferiores (no contaminada por salmuera), de una zona no afectada por la picnoclina, de modo que To=24ºC y ρo =
1048.3Kg/m³;
1027.8 X=75;
1024.2 Hipótesis vertido mediante ranura equivalente. Ocupa
(T=15ºC y 25 / 235 MU2 BV=8, BH=42
(T=27ºC y S=37.5psu) toda la columna de agua, los 8metros, a 40m de origen
S=37.5psu) S=29.7
1028.15 X=75;
1024.2 Hipótesis de vertido mediante ranura equivalente. Ocupa
(T=15ºC y 25.5 / 237 MU2 BV=8, BH=42
(T=27ºC y S=37.5psu) toda la columna de agua, los 8metros, a 40m de origen
S=38psu) S=29.7
149
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
C.3.) Representando una situación primaveral/estival en el Mar Mediterráneo, con presencia de termoclina a profundidades entre 15 y 20m. Consideramos
que la picnoclina está situada a una profundidad de 15m desde superficie, en un calado total de 25m, existiendo un saldo término entre la superficie y el fondo.
Consideramos que el agua de alimentación se extrae de profundidades inferiores (no contaminada por salmuera), de una zona no afectada por la picnoclina, de
modo que To=24ºC y ρo = 1048.3Kg/m³;
1024.2 1027.8
No ejecuta
(T=27ºC y (T=15ºC y 25.3 / 235.5 MU1V 8 7 11 16.7 NO
(Ua=0)
S=37.5psu) S=37.5psu)
1024.2 1028.15
No ejecuta
(T=27ºC y (T=15ºC y 25.3 / 235.5 MU1V 8 7 11.2 16.8 NO
(Ua=0)
S=37.5psu) S=38psu)
150
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
D.1.)Variamos U A desde 0 hasta 0.9m/s (incrementos de 0.05 y0.15m/s), para 0º , 90º (corriente de la misma dirección y sentido que el chorro).
D.2.)Variamos U A , desde 0 hasta 0.9m/s (incrementos de 0.05 y 0.15m), para 90º , 0º (corriente perpendicular al chorro).
151
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
D.3.) Variamos U A , desde 0 hasta 1m/s (incrementos de 0.1m), con 180º , 270º (corriente de la misma dirección y sentido opuesto al chorro).
0.05
NO EJECUTA EL CASO DE CHORRO EN LA MISMA DIRECCION PERO SENTIDO
0.1
OPUESTO A LA CORRIENTE
0.15
152
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
0.30
NOTA: Respecto al sistema de coordenadas, el eje “x” indica la dirección y sentido de la corriente en el medio receptor. De acuerdo con esto, para
0º , la dirección y sentido del chorro y la corriente coinciden y es “x”; para 90º , ambos son perpendiculares, y el eje “y” será la dirección y sentido del
chorro, mientras que el “x” sigue correspondiente a la corriente; para 180º ¸ misma dirección pero sentido opuesto, el chorro en el vertido tendrá valores de “x”
negativos, que van haciéndose positivos, cuando el chorro invierte su dirección por efecto de la corriente, que tiene la dirección y sentido del eje “x”.
U A 0 ; A 1027.5 Kg/m³; 0º (co-flowing), ho 1 ; 60º ; o 1052 Kg/m³; U o 5m/s ; D 0.2 m; U viento 0 ; RMZ=100m;
Parámetros fijos: :
L=150m;YB1=1000; YB2=1150; N=11;s=15; 90º ; 90º
Z(*) X(*) S
153
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
3 El programa no se ejecuta porque la altura del difusor es superior a 1/3 del calado disponible
NOTA: (*). En este caso el efluente alcanza la superficie, hay un cambio en el tipo de flujo: MNU7. De acuerdo con CORMIX2, el efluente no realiza una
trayectoria de chorro sino que ocupa toda la columna de agua. Se indican los valores a la misma distancia del punto de impacto del anterior flujo MNU1
H A =20m; U A 0 ; A 1027.5 Kg/m³; 0º (co-flowing), 60º ; o 1052 Kg/m³; U o 5m/s ; D 0.2 m; U viento 0 ; RMZ=100m;
Parámetros fijos: :
L=150m;YB1=1000; YB2=1150; N=11;s=15; 90º ; 90º
F.1.) Variamos la altura de la boquilla respecto al fondo: ho entre 0 y 1.5m (incrementos de 0.5m).
154
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
U A 0 ; H A =20m ; A 1027.5 Kg/m³; 0º (co-flowing), ho 1 ; 60º ; o 1052 Kg/m³; U o 5m/s ; D 0.2 m; U viento 0 ; RMZ=100m;
Parámetros fijos: :
L=150m;YB1=1000; YB2=1150; N=11; 90º ; 90º
G.1.) Variamos la separación (s) entre boquillas entre 2m y 20m (incrementos de 3m). Variamos con ello la longitud del difusor
155
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
U A 0 ; H A =20m ; A 1027.5 Kg/m³; 0º (co-flowing), ho 1 ; 60º ; o 1052 Kg/m³; U o 5m/s ; D 0.2 m; U viento 0 ; RMZ=100m;
Parámetros fijos: :
L=150m;YB1=1000; YB2=1150; s=15; 90º ; 90º
H.1.) Variamos el numero de elevadores (boquillas), (N) entre 5 y 25 (incrementos de 5). Variamos la separación entre boquillas, sin interacción.
156
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
I.1) Configuración: tramo difusor con elevadores de una boquilla, chorros unilaterales (al mismo lado), vertiendo en la misma dirección y sentido. (“single
nozzle: unidirectional: same direction”), 90º ; 90º , N=11.
Same direction , 30º 23 /215 MNU1 3.6 6 11.7 12.2 NO Mod. CORJET
Same direction , 60º 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO Mod. CORJET
X(m) S BH BV
Hip. Vertido en
Same direction , 75º 23 /215 MNU1 1.5 9.7 75.5 0.45
ranura equiv.
Hip. Vertido en
Same direction , 90º 23 /215 MNU1 1.5 9.7 75.5 0.45
ranura equiv.
157
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
Jet maximum
Impact point
height
Multiport Diffuser design Discharge angles Calculus
Zmax Xi
Seje (Zmax) Seje (Xi)
(m) (m)
X(m) S BH BV
158
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
I.2) Configuración: tramo difusor con elevadores de una boquilla, chorros unilaterales vertiendo en abanico (“single nozzle: unidirectional: fanned
out”), 90º 90º , N=11
Fanned out, 60º 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO CORJET
X(m) S BH BV
Hip. Vertido en
Fanned out, 75º 23 / 215 MNU1 1.5 9.7 75.5 0.45
ranura equiv.
Hip. Vertido en
Fanned out, 90º 23 /215 MNU1 1.5 9.7 75.5 0.45
ranura equiv.
I.3) Configuración: tramo difusor con elevadores de una boquilla, con chorros bilaterales o alternos (“single nozzle: alternating”).
VARIABLE
RESULTADOS CON CORMIX2
MODIFICADA
Nº de Froude CLASE DE
densimétrico FLUJO
Configuración X(m) S BH(m) BV(m) Cálculos
NOTA:
(*). En ambos casos, el modelo considera, sin que el usuario lo haya introducido en los datos de entrada, los valores: 90º y 0º . Como se
observa en los resultados, los cálculos realizados por CORMIX2 son idénticos para el caso de chorros en la misma dirección (“same direction”, como en abanico:
“fanned out”. En ambos casos asume la hipótesis de vertido mediante ranura equivalente.
159
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
I.4) Configuración: tramo difusor con elevadores de dos boquillas, separadas entre sí por un ángulo horizontal inferior a 60º, y chorros vertiendo en la
misma dirección. (“two nozzles: less than 60 degrees: same direction”). 0º ; 90º , 90º . N=22
I.5) Configuración: tramo difusor con elevadores de dos boquillas, separadas entre sí por un ángulo horizontal inferior a 60º y chorros vertiendo en abanico
(“two nozzles: less than 60 degrees : fanned out”). 0º ; 90º , 90º .
160
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
I.6) Configuración: tramo difusor con elevadores de dos boquillas, separadas entre sí por un ángulo horizontal de unos 180º (“two nozzles: about 180
degrees”).
hipótesis de vertido en
About 180º: Same direction 23 / 152 (**)MNU1 1.87 7.7 75.5 0.6
ranura equivalente
NOTA:
(**). En el caso de difusor con elevadores de dos boquillas y chorros bilaterales formando unos 180º entre sí (“about 180º”), asume que se trata de un
difusor con boquilla única vertiendo verticalmente, con ángulos: 90º y asume directamente la hipótesis de vertido en ranura equivalente.
I.7) Configuración: tramo difusor con elevadores de varias boquillas, con flujo neto de cantidad de movimiento horizontal igual a cero y misma dirección o en
abanico (“several nozzles: net horizontal momentum flux zero: same direction// fanned out”) (*)
NOTA:
(*) El modelo aplica las mismas fórmulas de cálculo y obtiene los mismos resultados en los casos: misma dirección (“same direction”) y en abanico
(“Fanned out”).
(**). En este caso de varias boquillas por elevador y flujo neto de cantidad de movimiento horizontal nulo, considera, para todos los casos, que se trata de
una única boquilla por elevador vertiendo verticalmente, asumiendo: 90º . Aplica la hipótesis de vertido mediante ranura equivalente.
161
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX
I.8) Configuración de varias boquillas por elevador, con flujo neto de cantidad de movimiento horizontal distinto de cero, misma dirección y en abanico
(“several nozzles: net horizontal momentum flux non-zero: same direction/ fanned out”) (*)
NOTA
(*) El modelo aplica las mismas fórmulas de cálculo y obtiene los mismos resultados en los casos: misma dirección (“same direction”) y en abanico
(“Fanned out”).
162