Cormix2 PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 65

Ficha 3.2.  
 
SUBSISTEMA CORMIX 2 

V. 1.1  

98

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.
La Guía Metodológica desarrollada en MEDVSA incluye los documentos:

- GUÍA METODOLÓGICA

- Anexo 1. Fichas descriptivas de los modelos comerciales

- Anexo 2. Modelos de simulación brIHne.

- 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 es parte del Anexo 1: Fichas descriptivas de los modelos


comerciales

Autores:

INSTITUTO DE HIDRÁULICA AMBIENTAL (IH CANTABRIA)

CENTRO DE ESTUDIOS DE PUERTOS Y COSTAS (CEDEX)


PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

Subsistema CORMIX2 de CORMIX:

VERTIDO SUMERGIDO MEDIANTE CHORROS MÚLTIPLES


1. INTRODUCCIÓN. AUTORES. DESARROLLO.
El subsistema CORMIX2: Submerged multiport diffusers (1) simula el vertido sumergido
mediante chorros múltiples, con tramo difusor. Permite simular diferentes configuraciones de tramo
difusor en el emisario (2) y es aplicable a chorros de flotabilidad nula, positiva y negativa. El modelo
permite tener en cuenta las condiciones en el medio receptor.
Los autores del programa estiman un error potencial (desviación típica) de los resultados
obtenidos respecto a mediciones reales “in situ”, de aproximadamente un +/-50%, para cualquier tipo
de vertido y efluente simulado (3).
La presente Ficha Descriptiva se centra en el modelado de vertidos de efluentes de
flotabilidad negativa, particularizando para el efluente salmuera. En las primeras versiones de
CORMIX, el modelado de vertidos de salmuera con CORMIX2 requería entrar en el subsistema y
considerar la concentración salina como un contaminante conservativo genérico, pero en las últimas
versiones (v.5 y v.6) el subsistema incluye explícitamente el efluente salmuera entre las opciones a
elegir y al seleccionarlo la interfaz del CORMIX2 se adapta a este tipo de efluente.

2. DESCRIPCIÓN DEL MODELO. BASE TEÓRICA.


Se trata de un modelo basado en el análisis dimensional del fenómeno (ver ficha general del
programa CORMIX), que incluye más de 31 clases de flujo diferentes y más de 25 módulos de cálculo
para la caracterización de estos flujos.
En los párrafos siguientes se definen las variables consideradas en el análisis dimensional de
CORMIX2 y los flujos (o variables integrales) que reúnen a los parámetros y fuerzas de mayor
influencia y que gobiernan el fenómeno, despreciando aquellas que resultan menos significativas. A
continuación se exponen las escalas de longitud que derivan de relacionar los flujos entre sí y con
otros parámetros determinantes en el fenómeno. Finalmente, se definen las categorías de CORMIX2
y las clases de flujo que incluyen cada una de ellas, que llevarán asociadas fórmulas de cálculo semi-
empíricas basadas en coeficientes experimentales, para caracterizar el comportamiento del efluente
en cada una de las regiones.
En los vertidos mediante chorros múltiples, CORMIX2 distingue dos casos: chorros que no
interaccionan entre sí antes de alcanzar el fondo (se comportan como chorros individuales), y chorros
que sí interaccionan. En este último supuesto, CORMIX2 asume la hipótesis de vertido mediante
ranura equivalente, que requiere la definición de nuevas escalas de longitud.

2.1. Variables axiales en CORMIX2.


La Figura 1 muestra las principales variables en la simulación de vertido con chorros múltiples
mediante el subsistema CORMIX2:

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

- Glosario de variables en el subsistema CORMIX2:

- H A : profundidad media del medio receptor en la zona de vertido (calado).


- U A : velocidad de la corriente en el medio receptor. Se considera siempre paralela a la costa.
-  A : densidad del fluido receptor
- DISTB : distancia, perpendicular a la costa, desde la orilla más cercana al emisario.
- ho : altura de las boquillas con respecto al fondo
- D : diámetro de las boquillas.
- Qo : caudal de vertido.
- U o : velocidad de salida del chorro.

-  o : densidad del efluente.

- C o : Diferencia entre la concentración inicial del efluente y la concentración en el medio receptor. En


nuestro caso: concentración salina del efluente menos salinidad del medio receptor.
-  : ángulo de inclinación del chorro respecto al fondo horizontal.
-  : ángulo horizontal entre el chorro y la corriente en el medio receptor
 
g  : gradiente de densidad relativa en el medio receptor, cuya dimensión es T 2 .  
 r z
 o : Gravedad reducida; g=valor de la gravedad (9.81m/s²).
go '  g
a
- LD : longitud del tramo difusor.
- S: separación entre boquillas.
-  : ángulo horizontal de los chorros respecto al tramo difusor.
-  : ángulo horizontal entre el tramo difusor y la corriente en el medio receptor.

2.2. Variables integrales o flujos en CORMIX2.


Las variables axiales se agrupan en variables integrales o flujos, distinguiéndose el
caso de interacción o no interacción entre chorros.

o 2.2.1 CASO DE NO INTERACCIÓN ENTRE CHORROS.

- Flujo de caudal: Q0   D 2U o , cuya dimensión es L / T


4
 3
. Representa el caudal de

efluente vertido al medio receptor.

- Flujo de cantidad de movimiento: M o  U o Qo , cuya dimensión es L / T  4 2


. Representa la
energía transmitida durante la descarga del efluente.

- Flujo de flotabilidad: J o  g o ' Qo , cuya dimensión es L / T  4 3


. Representa el efecto de la
gravedad sobre el efluente vertido.

o 2.2.1 CASO DE SI INTERACCIÓN ENTRE CHORROS.

- Flujo de caudal de vertido en ranura: q 0  Q o / L D , cuya dimensión es L / T  3


.

101
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

- Flujo de cantidad de movimiento de vertido en ranura: m o  M o / L D , cuya dimensión es


L 3

/T 2 .

- Flujo de flotabilidad de vertido en ranura: j o  J o / LD , cuya dimensión es L / T
3 3
.
2.3. Escalas de longitud en CORMIX2.
En un flujo tipo “chorro” dominan las fuerzas advectivas, de inercia, mientras que en un flujo
tipo “pluma” el movimiento viene determinado principalmente por las fuerzas de flotabilidad.
Se definen a continuación las escalas de longitud que utiliza el subsistema CORMIX2,
distinguiendo entre el caso de SI y NO interacción (solape) entre los chorros, y si se trata de un medio
receptor de densidad uniforme o estratificado:

o 2.3.1 CASO DE NO INTERACCIÓN ENTRE CHORROS.


En el caso de que no se produzca interacción entre chorros contiguos, las escalas de longitud
' '
son las mismas que para el chorro aislado tridimensional de CORMIX1: LQ , Lb , Lm , LM , Lm , Lb , ya
que el análisis dimensional conduce a los mismos monomios adimensionales, pero con una
interpretación completamente diferente, ya que los parámetros para el CORMIX2 definen el
comportamiento 3D global para el vertido total a través del difusor. Si la profundidad fuera muy
H s  ho l
grande  1 , las características del chorro global para  1 (los chorros individuales
LD LD
se habrían intermezclado formando un único chorro global) tenderán a ser iguales a las de un chorro
aislado con los mismos valores de los parámetros Lb , Lm , etc. Pero a distancias más pequeñas
l
(  1 ) estos parámetros, que seguirían marcando la diferencia entre los comportamientos 2D y 3D
LD
del difusor de longitud finita, son poco útiles y de hecho CORMIX2 no los utiliza ni para la clasificación
ni para la caracterización del flujo.

MEDIO RECEPTOR DE DENSIDAD UNIFORME.

LQ  LQ ,M : Escala flujo de caudal-cantidad de movimiento: Es la escala de longitud


representativa de la distancia a partir de la cual el comportamiento del flujo pasa de estar dominado
por el flujo de caudal a estar dominado por la cantidad de movimiento transmitida en la descarga.
Para distancias desde el punto de vertido mucho mayores que LQ, el comportamiento depende
fundamentalmente de la cantidad de movimiento, siendo despreciable el flujo de caudal.

LM  LM , B : Escala cantidad de movimiento-flotabilidad. Es la escala de longitud


representativa de la distancia a partir de la cual el comportamiento del flujo pasa de estar dominado
por la cantidad de movimiento transmitida en la descarga a estar dominado por la flotabilidad.
M o3 / 4
LM  Para distancias desde el punto de vertido mucho mayores que LM , el comportamiento
J o1/ 2
del efluente depende fundamentalmente de la flotabilidad, siendo despreciable el efecto del impulso
inicial.

Lm : Escala cantidad de movimiento-corriente del medio receptor (ambiental). Escala de


longitud representativa de la distancia a partir de la cual el chorro pasa de estar dominado por la
cantidad de movimiento transmitida en la descarga a estar dominado por la corriente en el medio

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.

Lb : Escala flotabilidad-corriente en el medio receptor. Escala de longitud representativa de la


distancia a la cual la pluma pasa de estar dominada por la flotabilidad a estar dominada por la
Jo
corriente en el medio receptor. Lb  . Para distancias desde el punto de vertido mucho mayores
U A3
que Lb, el comportamiento del flujo depende fundamentalmente de la corriente en el medio receptor,
siendo despreciable el flujo de flotabilidad.

La figura 2A muestra el significado físico de estas escalas de longitud, correspondientes a un


medio receptor homogéneo (de densidad uniforme):

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)

MEDIO RECEPTOR ESTRATIFICADO.

Lm ´ : escala cantidad de movimiento-estratificación del medio receptor. Escala de longitud


representativa de la distancia a la cual el chorro pasa de estar dominado por la cantidad de
M o1 / 4
movimiento de la descarga a estar dominado por la estratificación del medio receptor. Lm ´ .
 1/ 4
Para distancias desde el punto de vertido mucho mayores que Lm ´ , el comportamiento del flujo
depende fundamentalmente de la estratificación de densidad en el medio receptor, siendo
despreciable la influencia de la cantidad de movimiento transmitido en la descarga.

Lb ´ : Escala flotabilidad-estratificación en el medio receptor. Escala de longitud representativa


de la distancia a partir de la cual la pluma pasa de estar dominada por la flotabilidad a estar dominada
J 1/ 4
por la estratificación en el medio receptor. Lb ´ . Para distancias desde el punto de vertido
 3/8

103
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

mucho mayores que Lb ´ , el comportamiento del flujo depende fundamentalmente de la


estratificación, siendo despreciable la flotabilidad.
La figura 2B muestra el significado físico de las escalas de longitud de CORMIX2,
correspondientes a un medio receptor estratificado:

Figura 2B: Escalas de longitud de CORMIX2, para el caso de NO INTERACCIÓN entre chorros. Medio
receptor estratificado. Fuente: Manual de CORMIX (4)

o 2.3.2. CASO DE SI INTERACCIÓN ENTRE CHORROS.


Los chorros contiguos interaccionan entre sí antes de alcanzar el fondo; para este caso
CORMIX2 define nuevas escalas de longitud.

MEDIO RECEPTOR DE DENSIDAD UNIFORME.

l M : Escala cantidad de movimiento-flotabilidad de flujo bidimensional, en ranura.


Es la escala de longitud representativa de la distancia a partir de la cual el comportamiento
del flujo bidimensional, en ranura, pasa de estar dominado por la cantidad de movimiento transmitida
mo
en la descarga a estar dominado por la flotabilidad. l M  .Para distancias desde el punto de
jo2 / 3
vertido mucho mayores que l M , el comportamiento del efluente depende fundamentalmente de la
flotabilidad, siendo despreciable el efecto del impulso inicial.

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

MEDIO RECEPTOR ESTRATIFICADO.

l m ´ : Escala cantidad de movimiento-estratificación del medio receptor para flujo en ranura.


Escala de longitud representativa de la distancia a la cual el chorro pasa de estar dominado por la
cantidad de movimiento de la descarga a estar dominado por la estratificación del medio receptor.
m1o / 3
lm´ . Para distancias desde el punto de vertido mucho mayores que l m ´ , el comportamiento del
 1/ 3
flujo depende fundamentalmente de la estratificación de densidad en el medio receptor, siendo
despreciable la influencia de la cantidad de movimiento transmitido en la descarga.

lb ´ : Escala flotabilidad-estratificación en el medio receptor. Escala de longitud representativa


de la distancia a partir de la cual la pluma pasa de estar dominada por la flotabilidad a estar dominada
jo1 / 3
por la estratificación en el medio receptor. lb ´ . Para distancias desde el punto de vertido
 1/ 2
mucho mayores que l b ´ , el comportamiento del flujo depende fundamentalmente de la estratificación,
siendo despreciable la flotabilidad.

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.

Las figuras 3 y 4 muestran el significado físico de algunas de las anteriores escalas de


longitud:

Figura 3: Escalas de longitud de transición de CORMIX2, chorro-pluma y chorro-ranura, en un medio receptor


de densidad uniforme. Fuente: Manual CORMIX2 (1)

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

2.4. Categorías de clasificación para efluentes de flotabilidad negativa.


Estableciendo relaciones entre las escalas de longitud anteriormente definidas y otros
parámetros fundamentales en el comportamiento del vertido de efluente hiperdenso mediante chorros
múltiples (profundidad del fondo, ángulo de inclinación del chorro, ángulo horizontal entre el chorro y
la corriente en el medio receptor, etc.), CORMIX2 establece las siguientes categorías de flujo (que
incluyen vertidos sobre el fondo, cercano a la superficie, en medio homogéneo o estratificado, en
reposo o con presencia de corrientes):
o Categoría MS: vertido cercano al fondo, flujo atrapado en una picnoclina en un medio
receptor estratificado.
o Categoría IMS: vertidos mediante chorros múltiples, con emisario cercano a la superficie,
descargados en un medio receptor estratificado. El flujo queda atrapado en el nivel de la
picnoclina.
o Categoría MNU: vertido mediante chorros múltiples, con emisario cercano al fondo,
descargados en un medio receptor de densidad uniforme.
o Categoría IMNU: vertido mediante chorros múltiples, con emisario cercano a la superficie,
descargados en un medio receptor de densidad uniforme.
o Categoría A: vertidos mediante chorros múltiples, con emisario cercano al fondo, estando
el flujo afectado por atrapamiento o adhesión (“attachment”) dinámica con el fondo.

Estas categorías de CORMIX2 a su vez se subdividen en distintas clases de flujo (MNU1,


MNU2, etc.) que llevarán asociadas distintas formulaciones de cálculo sencillas y dependientes de
coeficientes experimentales. Los vertidos de efluentes hiperdensos desde la superficie o cercanos la
misma (IMU) serán descritos en mayor detalle en la ficha correspondiente al subsistema D-CORMIX:
Las figuras 5, 6, 7, 8 y 9 muestran los esquemas de clasificación en flujos, para las distintas
categorías de CORMIX2 aplicables a efluentes de flotabilidad negativa, incluyendo vertidos
sumergidos cercanos al fondo o a la superficie, en medio receptor homogéneo o estratificado, con
medio receptor en calma o presencia de corrientes en el medio receptor.

106
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

A) VERTIDOS MEDIANTE CHORROS MÚLTIPLES SUMERGIDOS. EFLUENTE DE


FLOTABILIDAD NEGATIVA

Figura 5: MEDIO RECEPTOR HOMOGÉNEO. EMISARIO CERCANO AL FONDO. AGUAS


PROFUNDAS.

Figura 6: MEDIO RECEPTOR HOMOGÉNEO. EMISARIO CERCANO AL FONDO. AGUAS


REDUCIDAS.

107
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

Figura 7: MEDIO RECEPTOR ESTRATIFICADO. EMISARIO CERCANO AL FONDO.

Figura 8: MEDIO RECEPTOR HOMOGÉNEO. EMISARIO CERCANO A LA SUPERFICIE.

108
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

Figura 9: MEDIO RECEPTOR ESTRATIFICADO. EMISARIO CERCANO A LA SUPERFICIE.

2.5. Clases de flujo en CORMIX2 para la clasificación para efluentes de


flotabilidad negativa.
En los esquemas de las páginas anteriores se han mostrado las categorías de CORMIX2
para vertidos de flotabilidad negativa, así como las clases de flujo en que se dividen estas categorías
según los valores de los ratios entre las escalas de longitud y de éstas con otros parámetros
fundamentales.
Para el rango de valores realistas de los flujos en vertidos al Mar mediterráneo de salmuera
(procedentes de osmosis con tasas de conversión en torno al 50%), con descargas cercanas al fondo
y números de Froude entre 10 y 25, en un medio receptor en reposo y densidad homogénea, difusor
lineal y chorros unilaterales, con una sola boquilla por elevador y chorros perpendiculares al tramo
difusor, las clases más frecuentes en que CORMIX2 clasifica el flujo son:
- MNU1: vertido en una zona de profundidad suficiente, de modo que los chorros no impactan
con la superficie. En el comportamiento del flujo domina la flotabilidad frente al impulso inicial. El
medio receptor es homogéneo y se encuentra en reposo o con presencia de una corriente ambiental
débil de la misma dirección y sentido que el chorro. El impacto con el fondo se produce en un ángulo
prácticamente vertical, tras lo cual el efluente se expande radialmente por el fondo.
-MNU2: vertido en una zona de calado suficiente, de modo que los chorros no impactan con
la superficie. En el comportamiento del flujo domina el impulso inicial frente a la flotabilidad. El medio
receptor es homogéneo y la corriente en el medio receptor es significativa y de la misma dirección y
sentido que el chorro.
- MNU7. Vertido en una zona de calado limitado, de modo que los chorros impactan con la
superficie. El medio receptor es homogéneo y la corriente ambiental en el medio receptor es paralela
a los chorros.

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):

3.1. Fenómeno de interacción entre chorros.


Como se ha comentado anteriormente, cuando CORMIX2 obtiene que la separación entre
boquillas es menor que el ancho de los chorros en toda la trayectoria, asume que no hay interacción
entre chorros y los trata como chorros individuales tridimensionales.
Sin embargo, cuando el ancho del chorro iguala a la separación entre boquillas, CORMIX2
considera que se produce la interacción entre chorros (“merging”), y aplica la hipótesis de vertido
mediante ranura equivalente (“equivalente slot diffuser”). Esta hipótesis asume que el flujo de
descarga a través de un difusor con boquillas igualmente espaciadas, cuyos chorros 3D interaccionan
entre sí, es equivalente al flujo de descarga a través de una ranura con la misma longitud que el
difusor, el mismo flujo de caudal ( Q0   D 2U o ) y el mismo flujo de cantidad de movimiento
4
( M o  U o Qo ), flotabilidad y concentración de trazador. Este flujo en ranura equivale a un chorro
plano de sección rectangular, de modo que el problema 3D de chorros se transforma en un problema
bidimensional, lo que facilita su resolución numérica.
La interacción entre chorros supone una transición desde los perfiles de distribución tipo
Gauss hacia un perfil horizontal de concentraciones homogéneas a lo ancho de la pluma que se
forma al considerar la mezcla completa entre los chorros. Sin embargo, CORMIX2 no simula este
cambio suave sino que automáticamente homogeneiza la concentración y lo considera un flujo
bidimensional.
Las figuras 10, 11 y 12 muestran el paso de chorro circular tridimensional a pluma
bidimensional, que CORMIX2 realiza cuando detecta interacción entre chorros y asume la hipótesis
de ranura. La figura 10 señala el caso de difusor con una boquilla por elevador y chorros unilaterales
vertiendo en la misma dirección.

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)

La figura 12 muestra la sección transversal de un chorro individual tridimensional cuyas


variables en la sección transversal sigue una distribución de tipo Gauss (izquierda), así como la
sección transversal correspondiente a la pluma bidimensional que se asume tras la interacción entre
los chorros (derecha), donde la distribución de variables en homogénea en dirección horizontal y es
de tipo Gauss en dirección vertical.

Figura 12: Secciones transversales del chorro circular 3D y de la pluma bidimensional


(asumida en la hipótesis de vertido mediante ranura equivalente). Fuente: Manual CORMIX
(3)
Izquierda: Sección del chorro. R=Radio, S=dilución en el eje (centro de sección).
Derecha: Sección de pluma. BH=ancho lateral; BV=altura; S=dilución en eje-

112
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

3.2. Configuración del dispositivo de vertido (tramo difusor).


También es importante destacar las hipótesis de CORMIX2 en relación con las
configuraciones del dispositivo de vertido (tramo difusor). A pesar de que teóricamente el
subsistema puede simular una amplia variedad de diseños (ver la sección 4 de esta ficha), los
módulos de cálculo aplican simplificaciones que reducen muy significativamente las opciones
disponibles.
Recordamos a continuación la definición de los distintos ángulos que determinan la
configuración del dispositivo de vertido, que se ilustran además mediante la figura 13:
-  : ángulo de inclinación de salida del chorro.
-  : ángulo horizontal entre el chorro y la corriente en el medio receptor

-  : ángulo horizontal de los chorros respecto al tramo difusor.


-  : ángulo horizontal del tramo difusor con la corriente en el medio receptor.

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:

Figura 14: Simplificación de CORMIX2 para dispositivos de vertido de dos


boquillas por elevador y chorros bilaterales formando 180º. Izquierda: dispositivo
real. Derecha: simplificación de CORMIX2. Vistas en alzado

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

Simplificaciones en los dispositivos de vertido con varias boquillas por elevador


(roseta):
En el caso de un tramo difusor con varias boquillas por elevador (“several nozzles per
riser”), se asumen las mismas simplificaciones que para los casos anteriores. Así, si el flujo de
momento neto horizontal es cero, CORMIX2 considera que el ángulo de inclinación del chorro es de
90º y aplica directamente la hipótesis de vertido mediante ranura equivalente.
Si por el contrario, el flujo de momento neto horizontal no es cero, se considera que todos los
chorros se agrupan en una dirección formando un ángulo igual o inferior a 60º y se considera una
única boquilla, de modo que sólo es necesario introducir los valores de  ,  ,  (Ver Sección 5 y
Anexo 3 de esta ficha).
Algunos autores (4) han detectado sobrestimación del orden de 45% en los resultados de
CORMIX2 para la dilución en campo cercano.

En conclusión, a pesar de las múltiples opciones de configuración que ofrece CORMIX2 en


los datos de entrada, en sus cálculos sólo distingue entre los siguientes casos:
o CONFIGURACIÓN 1: difusor con elevadores de una boquilla, chorros unilaterales, con
inclinación vertical y perpendiculares al difusor. Utiliza el módulo CORJET.
o CONFIGURACIÓN 2: difusor con elevadores de una boquilla, descargando con inclinación
vertical del chorro y perpendicular al tramo difusor (
    90º ). En este caso asume la
hipótesis de vertido mediante ranura equivalente.

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:

CONFIGURACIONES PARA LAS QUE CORMIX2 UTILIZA LAS MISMAS FÓRMULAS DE


CÁLCULO
CÁLCULO (IDÉNTICOS RESULTADOS EN CADA CASO)
Tramo difusor de una boquilla por elevador, con chorros unilaterales, vertiendo en la misma
dirección o en abanico y   75º
(“single nozzle: unidireccional: same direction / fanned out” y   75º ). (Casos G.1 y G.2) CONFIGUR 1
(*1)
Aplica
Tramo difusor de dos boquillas por elevador, separadas por un ángulo horizontal inferior a
CORJET
60º, vertiendo en la misma dirección o en abanico y   75º
(“two nozzles: less than 60 degrees: same direction / fanned out” y   75º ). (Caso G.4)

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

4. COMPONENTES Y OPCIONES DEL MODELO.


La tabla 2 indica los datos de entrada necesarios para el modelo CORMIX2, así como los
resultados que ofrece el subsistema. Entre paréntesis el nombre con que aparecen en CORMIX2:

DATOS DE ENTRADA AL SUBSISTEMA CORMIX2.


Situación de la orilla más cercana (Nearest bank is on)
Distancia desde la orilla hasta el extremo más cercano / lejano del difusor (Distance to one endpoint/ other).
Longitud del tramo difusor (diffuser length),
Altura de la boquilla con respecto al fondo (diffuser height).
Diámetro de la boquilla (Port diameter).
Coeficiente de contracción del diámetro (Contraction ratio), Número de boquillas (Total number openings).
Angulo de alineación entre el tramo difusor y la corriente (Alignment Angle GAMMA).
CONFIGURACIÓN DEL DISPOSITIVO DE VERTIDO (TRAMO DIFUSOR):
Boquillas por elevador: uno, dos o varios (Nozzles per riser: single, Two, Several).
Tipo de difusor:
1). Difusor de una boquilla por elevador (single ):
1.1. Difusor con chorros unilaterales (unidireccional):
θ Angulo de inclinación del chorro (THETA).
β: ángulo horizontal de los chorros respecto al tramo difusor (BETA).
Condiciones de
σ: ángulo horizontal entre los chorros y la corriente en el medio receptor (SIGMA).
descarga
Dirección de los chorros: misma dirección (same direction) y en abanico (fanned out)
(DISCHARGE:
CORMIX2) 1.2. Difusor con chorros bilaterales, alternos (Alternating):
Difusores con la misma dirección (Same direction) o dispuestos en abanico (Fanned out).
2). Difusores con dos boquillas por elevador (two nozzles per riser).
- Boquillas separadas entre sí un ángulo horizontal inferior a 60º (Less than 60 Degrees):
θ Angulo de inclinación del chorro (Vertical- THETA”).
β: ángulo horizontal de los chorros respecto al tramo difusor (Relative orientation- BETA).
σ: ángulo horizontal entre los chorros y la corriente en el m. receptor (Horizontal-SIGMA).
- Boquillas formando entre sí un ángulo de unos 180º (About 180 Degrees):
Chorros en la misma dirección (same direction).
Chorros vertiendo en abanico (fanned out).
3). Difusores con varias boquillas por elevador (Several Nozzles per riser).
- Flujo neto de momento nulo (“Net horizontal…zero”): chorros en la misma dirección (same direction) o en
abanico (fanned out).
- Flujo neto de momento distinto de cero (Net horizontal…Nonzero).
-Valores de los ángulos THETA, BETA, SIGMA.
Caudal/ velocidad de salida del chorro (Flow Rate/ velocity).
Efluente vertido
Concentración en la descarga (concentration). (*)
(EFFLUENT)
Densidad (density), Temperatura (Temperature).

Calado medio en la zona de vertido (average depth).


Condiciones
Profundidad del dispositivo de vertido (Dicharge depth).
Condiciones geométricas
ambientales, en Medio receptor ilimitado (Unbounded), o confinado (bounded).
el medio
Medio receptor homogéneo (Uniform) o estratificado (stratified). Temperatura
receptor Parámetros físicos
(“Temperatura”), Densidad (density).
(AMBIENT)
Condiciones Caudal / Velocidad del medio receptor (Flowrate / velocity).
hidrodinámicas Coeficiente de rozamiento con el fondo (Darcy, Manning”). Velocidad del viento (wind).
RESULTADOS EN CORMIX2.
Informe con los datos del caso simulado (Session Report).
Resultados numéricos de la simulación (Prediction file).
Clasificación y descripción de la clase de flujo (Flow class description).
Datos de salida Recomendaciones para el diseño (Design Recommendations).
(OUTPUT) Pasos intermedios, flujos, reglas de calculo, etc. (Ambietn Rulebase Browser, etc.).
-Valores y gráficos de la trayectoria del chorro: vista en planta y laterales (2D y 3D) del efluente en zona de
mezcla inicial, mediante la herramienta CORVUE.
Valores y gráficos de la concentración salina/ dilución en toda la trayectoria, y en puntos más característicos.

117
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

Tabla 2: Datos de entrada de alimentación del subsistema CORMIX2, y resultados.


NOTA (*): el programa solicita definir las unidades de la concentración de contaminante, pero
no incluye la opción “psu”. Los “psu” equivalente a ppm, sin embargo este dato de unidad no es
relevante puesto que no interviene en los cálculos, sólo se utiliza para expresar los resultados

Como sistema de referencia CORMIX2 utiliza un sistema cartesiano. En la hoja de resultados,


para la coordenada vertical (z) se considera como origen el fondo, mientras que para las coordenadas
horizontales (x,y) se toma como origen la boquilla de salida de cada chorro.

La figura 16 (coincidente con la Figura 1) muestra el esquema de los principales variables de


entrada (parámetros del efluente y del vertido y condiciones en el medio receptor) para el modelo
CORMIX2, en la simulación de vertidos con tramo difuso de chorros múltiples:

Figura 16: Vista lateral y en planta del esquema de vertido de CORMIX2.

Glosario de parámetros de entrada en CORMIX2:


- H A : profundidad media del medio receptor en la zona de vertido (calado).
- U A : velocidad de la corriente en el medio receptor. Se considera siempre paralela a la costa.
-  A : densidad del fluido receptor
- ho : altura de las boquillas con respecto al fondo
- D : diámetro de las boquillas.

- U o : velocidad de salida del chorro.

-  o : densidad del efluente.

- C o : Diferencia entre la concentración inicial del efluente y la concentración en el medio receptor.

-  : ángulo de inclinación del chorro respecto al fondo horizontal.


-  : ángulo horizontal entre el chorro y la corriente en el medio receptor

- LD : longitud del tramo difusor.


- S: separación entre boquillas.

-  : ángulo horizontal de los chorros respecto al tramo difusor.

118
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

-  : ángulo horizontal entre el tramo difusor y la corriente en el medio receptor.

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

Velocidad de corriente 0 – 10m/s Temperatura del efluente 4 – 100ºC

900 – 1200 Angulo de inclinación del chorro


Densidad del medio receptor 0 – 90º
Kg/m³ :
0 – 360º
Ángulo horizontal entre los
Temperatura del medio receptor 4 – 100ºC (excluyendo el
chorros y la corriente: 
rango 135º - 225º)
Distancia de la boquilla a la orilla
0 – 106 m Longitud del difusor 0 – 5000m
más cercana

Diámetro de boquilla 0.0001 – 20m Ratio de contracción del chorro 0.6 – 1

Altura de la boquilla respecto al


0 – 1000m Número de boquillas 3 - 10000
fondo

Ángulo horizontal entre el tramo


Velocidad de salida del chorro 0 – 50m/s difusor y la corriente:  0 – 180º

Ángulo entre los chorros y el


Caudal de efluente vertido 0 – 1000m³/s 0 – 90º
tramo difusor

Concentración de contaminante 0 - 109

Tabla 3. Rango de valores de variables

La tabla 4 indica las principales posibilidades y limitaciones del subsistema CORMIX2,


entendidas posibilidades como opciones o ventajas del programa, y limitaciones como restricciones o
carencias, para la simulación de vertidos de efluentes de flotabilidad negativa; a estas posibilidades y
limitaciones se suman las expuestas en la ficha descriptiva general del software CORMIX:

119
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

POSIBILIDADES DE CORMIX2 (aplicación SALMUERA)

- 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.

LIMITACIONES DE CORMIX2 (aplicación SALMUERA).


- Simplificación de la geometría del medio receptor a un canal recto y uniforme, de sección rectangular y fondo
plano. Las versiones más recientes del CORMIX tienen en cuenta la pendiente en el fondo, pero de un modo
simplificado y no en la simulación del comportamiento de la pluma hipersalina en campo lejano.
- Respecto a la geometría y dimensiones de cada variable, el modelo presenta las siguientes limitaciones:
* Limitación a una altura de difusor máxima igual a 1/3 del calado total.
* En vertidos prácticamente verticales (alta inclinación del chorro), el diámetro de la boquilla debe ser menor que el
calado total (D<H), mientras que en los prácticamente horizontales, debe ser menor que 1/3 del mismo (D<1/3D)
* La picnoclina ha se situarse en una posición que esté entre el 40%-90% del calado total.
* La longitud del tramo difusor debe de ser mayor a la altura de la columna de agua.
- El programa no simula todos los casos de chorro que se vierten en sentido contrario a la corriente en el medio
receptor. En particular, no simula el rango de valores  (entre chorro y corriente) de 135º-225º, por ser un vertido de
flujo complejo con recirculación.
- No permite discretizar la columna vertical en niveles. La corriente en el medio receptor, por tanto, se asume
constante en intensidad y dirección en toda la columna de agua. Respecto a la distribución de la densidad en la
columna de agua, CORMIX2 permite únicamente elegir entre tres distribuciones fijas.
- Aunque teóricamente el programa ofrece múltiples diseños en cuanto a configuración del dispositivo de vertido, a la
hora de ejecutar los cálculos simplifica estas configuraciones a un número muy reducido de casos, asumiendo
hipótesis que invalidan la fiabilidad de los resultados para vertidos de flotabilidad negativa.
- Limitaciones en el diseño de los chorros del tramo difusor: boquillas equiespaciadas, del mismo diámetro y
altura con respecto al fondo y vertiendo con la misma inclinación de chorro.
- Trabaja con valores constantes. No permite introducir regímenes temporales de ninguna de las variables.
- El algunas clases de flujo de CORMIX2, especialmente en aquellos casos en que el calado es limitado y el flujo
ocupa toda la columna de agua, los módulos de cálculo son de tipo “caja negra”, basados en fórmulas experimentales
en general poco contrastadas y que no describen la evolución del flujo.
- La hipótesis de vertido mediante ranura equivalente (“slot equivalent diffuser”) implica simplificaciones que en
muchas ocasiones llevan a resultados poco realistas.
- La mayor parte de las fórmulas de cálculo, coeficientes experimentales e hipótesis simplificativas, se han deducido
en base a ensayos experimentales con vertidos de flotabilidad positiva. Como se verá en las próximas secciones los
estudios de calibración y validación para efluentes de flotabilidad negativa con muy escasos e incluso inexistentes,
según la clase de flujo. Para el caso de interacción del efluente con los contornos, los autores no han presentado
ningún estudio de validación.
- No se han encontrado en la bibliografía disponible, estudios de validación para el caso de interacción entre chorros
de flotabilidad negativa, que demuestren la validez y fiabilidad de la hipótesis de vertido mediante ranura equivalente.

Tabla 4: Posibilidades y limitaciones en la simulación de vertidos de efluentes hiperdensos con el subsistema


CORMIX2.

120
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

La figura 16 muestra gráficamente las principales limitaciones en la geometría del cuerpo


receptor de los modelados con el subsistema CORMIX2:

Figura 16: Principales limitaciones geométricas y límites de aplicabilidad en


CORMIX2.

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

UNIDIRECTIONAL DIFFUSSER: DISPOSITIVO DE VERTIDO DE UNA BOQUILLA POR ELEVADOR.


CHORROS UNILATERALES
Dispositivo con los chorros vertiendo hacia el mismo lado del tramo difusor.
De acuerdo con las recomendaciones de los autores, esta configuración genera un momento neto directo y forma
fuertes corrientes en el medio receptor. Se recomienda este dispositivo, colocando los chorros perpendiculares a la
corriente predominante, cuando ésta se presenta con una dirección constante. Se recomienda con alineación
paralela a la corriente si la corriente es débil.

STAGGED DIFFUSER: VERTIDOS CON DIFUSOR DE CHORROS EN ABANICO.


Tramo difusor con chorros cuya dirección de descarga es aproximadamente la del tramo difusor, con pequeñas
desviaciones hacia cualquiera de los lados.
De acuerdo con los autores, este dispositivo genera también un momento directo al fluido receptor, lo que puede
generar fuertes corriente en el mismo. Para efluentes de flotabilidad negativa se desaconseja este dispositivo

ALTERNATING DIFFUSER: VERTIDOS CON DIFUSOR DE CHORROS BILATERALES, ALTERNOS


Dispositivo con chorros vertiendo hacia ambos lados del tramo difusor.
De acuerdo con los autores esta configuración no genera un momento neto directo en el fluido del medio receptor
por lo que no genera fuertes corriente en el mismo, lo que supone un menor impacto sobre el medio bentónico.
El vertido es recomendable en descarga en aguas profundas, con presencia de corrientes, o también en aguas
reducidas cuando se pretende que las corrientes tengan una influencia mínima en el comportamiento del chorro.

Figura 17: Configuraciones de vertido. Tipos de difusores en CORMIX2.

121
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

5. ANÁLISIS DE SENSIBILIDAD DEL SUBSISTEMA FRENTE A VARIACIONES


EN LOS PARÁMETROS DE ENTRADA.
Se ha realizado un análisis de sensibilidad del modelo para conocer cómo la variación del
valor de lo distintos parámetros de entrada (características del efluente, condiciones de descarga y
condiciones en el medio receptor) influyen en los resultados obtenidos con CORMIX2. La horquilla de
variación de estos parámetros responde al rango de valores realistas para un vertido de salmuera en
el Mar Mediterráneo (ver Anexo 2).
Para realizar el análisis se ha considerado un vertido sumergido mediante chorros múltiples
individual. Para todos los casos (salvo el caso de sensibilidad frente a la configuración del dispositivo
de vertido) se ha considera un tramo difusor con elevadores de una boquilla, chorros unilaterales
formando 90º con el tramo difusor.
Partiendo de este esquema, se han ido variando algunos parámetros de entrada,
manteniendo los demás parámetros constantes, para ver la influencia en la geometría y en la dilución
del efluente.
En particular se han considerando las variaciones de: inclinación del chorro, densidad del
efluente y del medio receptor, estratificación en la columna de agua, velocidad y dirección de la
corriente, número de boquillas, separación entre boquillas (considerando el caso de no interacción e
interacción entre chorros), configuración del dispositivo de vertido, calado en el medio receptor y
altura de la boquilla con respecto al fondo. Se han comparado los resultados de geometría y dilución
en el punto de máxima altura del chorro y en el punto de impacto con el fondo.
En el Anexo 3 de esta ficha, se muestran los resultados numéricos obtenidos del análisis de
sensibilidad. En todos los casos se considera que la dirección de la corriente coincide con la del eje x.
La Tabla 5 resume las principales conclusiones del análisis de sensibilidad de CORMIX2, en
el rango de valores de las variables característico de vertidos de salmuera en el Mar Mediterráneo:

VARIABLE SENSIBILIDAD COMENTARIOS RESPECTO A LOS RESULTADOS DE CORMIX2

A mayor ángulo, mayor es la cota alcanzada por el chorro y en general la


dilución.
En el caso de chorros que no interaccionan entre si durante su trayectoria
parabólica, la dilución presenta una sensibilidad muy alta, para medio receptor en
calma, con valores máximos de dilución para ángulos entre 45º y 60º. Si el medio
Ángulo receptor presenta corrientes en el medio receptor significativas, el grado de
inclinación Muy alta dilución es muy poco sensible a variaciones en el ángulo. Para ángulos
superiores a 75º, asume una inclinación vertical:   90º , aplicando la hipótesis
chorro
de vertido en ranura equivalente (“equivalent slot diffuser”).
En el caso de que exista interacción entre chorros, la dilución calculada por
CORMIX2 es muy poco sensible al ángulo de inclinación del chorro, dado que
todos se aproximan asumiendo las hipótesis: si   45º  0º y si   45º  90º

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

Si el medio es homogéneo (típica situación de invierno/otoño en el Mediterráneo,


a las profundidades de vertido) la influencia de la temperatura es pequeña.
Estratificación Si existe estratificación en la columna de agua (situación estival en el
Variable M.Mediterráneo, con un picnoclina entre 15 y 20m bajo la superficie), la influencia
y Temperatura sobre el efluente depende de si el chorro interacciona o no con la picnoclina.
- Si el chorro no alcanza la picnoclina, la influencia de la estratificación en la
columna de agua sobre el chorro es nula. Trayectoria de chorro

A mayor valor de la corriente en el medio receptor, mayor es la dilución del


efluente.
Ejemplo: para la configuración de vertido considerada en el análisis y una
Velocidad de la corriente en la dirección y sentido del chorro, CORMIX2 estima diluciones en
corriente torno a 15 en el punto de impacto en un medio receptor en reposo, mientras que
Muy alta
a velocidades superiores (0.3m/s), la dilución estimada es de 45.
UA
Superada una cierta velocidad en el medio receptor (en nuestro caso. 0.3-
0.45m/s), en que el flujo se clasifica tipo NV1, la dilución prácticamente ya no se
incrementa con la velocidad de la corriente, incluso en ciertos casos se
desfavorece.

Se ha analizado la sensibilidad para los ángulos


  0º ,90º ,180º (chorro de
la misma dirección y sentido, perpendicular y de sentido opuesto a la corriente)

Para   0 , a mayor intensidad de la corriente, menor es la altura máxima del


chorro y mayor el alcance hasta el impacto, y la dilución.
Angulo
horizontal entre Para   90 , el chorro es desviado hacia la dirección de la corriente. A mayor
la corriente y Alta intensidad, menor altura máxima y menor alcance del chorro (Y), pero mayor
chorro alcance en dirección de la corriente (X), y mayor dilución.

Para   180 º , CORMIX2 no se ejecuta, argumentando que existe fenómenos
de recirculación, que el modelo no es capaz de simular
CORMIX2 estima diluciones significativamente mayores en el caso de corrientes
perpendiculares al chorro, respecto al caso de corrientes en la misma dirección y
sentido del chorro.

Calado en el Si el chorro no interacciona con la superficie, la altura de la columna de agua por


medio receptor encima del efluente NO influye en el comportamiento del mismo.
Variable
Si el chorro SI interacciona con la superficie, los resultados indican que la altura
HA de la columna de agua apenas influye sobre la geometría y dilución del chorro.

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

Separación Si los chorros que NO interaccionan entre sí antes de alcanzar el fondo, la


entre boquillas separación entre boquillas no influye en el comportamiento del chorro.
Variable
s Si los chorros SÍ interaccionan entre sí, según CORMIX2, la separación entre
boquillas tampoco influye en el comportamiento del chorro.

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.

Como se ha descrito anteriormente, las opciones de cálculo disponibles en


CORMIX2 son en realidad muchas menos que las que ofrece en la interfaz.
CORMIX2 iguala muchas de las configuraciones disponibles, de modo que se
obtienen cálculos idénticos (Ver Sección 3).

Si en las configuraciones que lo permiten (“one nozzle per riser: undirectional


“;“two nozzles per riser: less than 60 degrees”).se introduce un ángulo de
inclinación del chorro   75º , CORMIX2 asume   90º .
El comportamiento del chorro NO es sensible a que se opte por la opción de
chorros en la misma dirección (“same direction”) o en abanico (“fanned out”), en
cualquiera de las configuraciones posibles.

El comportamiento del chorro NO es sensible a que se opte por la opción: de


Configuración tramo difusor de una boquilla por elevador y chorros unilaterales (“one nozzle per
Variable
del difusor riser: undirectional”) o la opción de Tramo difusor con dos boquillas por elevador,
separadas entre sí por un ángulo horizontal inferior a 60º. (“two nozzles per riser:
less than 60 degrees”), dado que en este último caso se simplifica unificando
ambas boquillas en una sola perpendicular al chorro.

El comportamiento del chorro NO es sensible a que se opte por las


configuraciones: Tramo difusor de una boquilla por elevador y chorros bilaterales
o alternos (“single nozzle: alternating”), o tramo difusor con dos boquillas por
elevador separadas 180º entre sí (“two nozzles: about 180 degrees”). Para los
tres dispositivos, CORMIX2 considera el vertido a través de una única boquilla
(con el caudal suma de todas las boquillas del elevador) descargando con
inclinación vertical y chorro perpendicular al difusor (     90º ). CORMIX2
aplica la hipótesis de vertido mediante ranura equivalente. Si el caudal total por
elevador y velocidad son los mismos, se obtendrán idénticos resultados.

Tabla 5: Conclusiones del análisis de sensibilidad del subsistema CORMIX2 frente a variaciones en los datos de
entrada

En conclusión, a pesar de las múltiples opciones de configuración que ofrece CORMIX2 en


los datos de entrada, en sus cálculos sólo distingue entre los siguientes casos:
o CONFIGURACIÓN 1: vertido en tramo difusor de una boquilla por elevador, chorros
unilaterales, perpendiculares al difusor. Para los cálculos utiliza el módulo CORJET, y en
sensible a todos los parámetros a los que CORJET es sensible (ver ficha descriptiva de
CORJET).
o CONFIGURACIÓN 2: vertido en tramo difusor de una boquilla por elevador, descargando
con inclinación vertical del chorro y perpendicular al tramo difusor (     90º ). En
este caso asume la hipótesis de vertido mediante ranura equivalente.

124
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

6. CALIBRACIÓN Y VALIDACION DEL MODELO, PRESENTADA POR LOS


AUTORES, PARA EFLUENTES DE FLOTABILIDAD NEGATIVA.
A pesar de que el subsistema CORMIX2 es capaz de simular un gran número de clases de
flujo, los datos de validación presentados por los autores son bastante limitados, principalmente en el
caso de efluentes de flotabilidad negativa.
La única validación que se presenta para este tipo de efluentes utiliza los resultados del
estudio experimental realizado por Tong et al (1979), (5), que corresponde a un vertido mediante
tramo difusor con elevadores de una boquilla descargando con inclinación vertical (   90º ),
descargando en un medio receptor homogéneo y con presencia de una corriente en la misma
dirección y sentido que los chorros (   0º ). El ratio entre la longitud del difusor y el diámetro de las
boquillas es igual a 57.7. El número de Froude densimétrico de cada chorro es Fo=153, y el número
de Froude densimétrico bajo la hipótesis de vertido mediante ranura equivalente es Fro=208.1. El
cociente entre la velocidad de salida del chorro y la de la corriente en el medio receptor es R=9.36 (lo
que equivale, para una velocidad de salida del chorro de 4.5m/s, a una corriente en el medio receptor
de 0.46m/s, en prototipo).
Para este caso de ensayo, CORMIX2 clasifica el flujo en la clase MNU2, sin interacción de
los chorros entre sí y con presencia de una corriente en el medio receptor significativa, aplicando para
los cálculos el módulo CORJET.
La figura 18 representa gráficamente los resultados de trayectoria y geometría del chorro,
obtenidos visualmente de los ensayos, y numéricamente de la simulación con CORMIX2. La
validación no incluye valores de dilución.

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.

Comparando los resultados obtenidos experimental y numéricamente, se observa cómo la


dimensión del chorro estimada por CORMIX2 es menor que la observada experimentalmente, lo que
implicaría, en principio, predicción de diluciones por debajo de las diluciones reales. El alcance del
chorro es mayor en las predicciones de CORMIX2 que en los resultados experimentales.
Es importante señalar que de los más de quince clases de flujo que incluye el CORMIX2 para
vertidos mediante chorros sumergidos desde el fondo de efluentes hiperdensos, y de los más diez
que incluye para vertidos hiperdensos desde superficie (ver Figuras 5, 6, 7, 8 y 9), los autores no
presentan ningún dato de calibración de las fórmulas de cálculo para este tipo de flujos de flotabilidad
negativa. Respecto a los estudios de validación de las fórmulas del subsistema, sólo se encuentra

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

7. VALIDACION DEL MODELO CON DATOS EXPERIMENTALES DE ENSAYOS


EN MODELO FÍSICO.

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

8. RECOMENDACIONES DE APLICACIÓN Y MANEJO DEL MODELO.


Revisado el subsistema CORMIX2, para el caso de vertido de salmuera mediante chorros
múltiples sumergidos, con emisario cercano al fondo, se establecen las recomendaciones expuestas
en las tablas 6 y 7 respecto al modelado y uso y manejo del programa. A estas recomendaciones hay
que añadir las expuestas en la ficha general del programa CORMIX:

RECOMENDACIONES RESPECTO AL MODELADO

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.

Respecto a la configuración del dispositivo de vertido (tramo difusor), para el caso


de efluentes hiperdensos, se recomienda seguir los consejos, recordando siempre que para
CORMIX el eje “x” coincide con la dirección de la corriente y de la línea de costa:
o 1.). Si se pretende construir un tramo difusor paralelo a la costa y perpendicular al
emisario se recomienda utilizar un diseño de tramo difusor de elevadores con una
boquilla, chorros unilaterales, vertiendo perpendicularmente al tramo difusor, hacia la
dirección de mayor profundidad y chorros con ángulos de inclinación   45º 60º .
En CORMIX2, el modelado se realizaría eligiendo las opciones de configuración:
“single nozzles per riser: unidireccional: same direction” y los ángulos:   90º ;
  270º ,   0º y   45º 60º . En este caso YB1=YB2. En las versiones
actuales el valor de  está limitado al rango 0-90º, por lo que no distinguen entre
vertido dirigido hacia costa o hacia mar adentro.
Es importante en esta simulación asegurarnos de que CORMIX2 clasifica el flujo como
estable (MNU1 o MNU2), de modo que aplique el modulo de calculo de CORJET. Se
recomienda diseñar con separación suficiente entre boquillas para evitar la interacción
entre chorros y que el subsistema aplique la hipótesis de vertido mediante ranura, dado
que no ha sido validada por los autores para efluentes hiperdensos.
o 2.). Si se pretende construir un emisario con tramo difusor contiguo al emisario y
perpendicular a la costa (y a la corriente), se recomienda utilizar el diseño de chorros
unilaterales ya descrito, o la configuración de doble boquilla por elevador, formando
un ángulo de 180º entre sí (configuración en T), y descargando en ambos casos con
un ángulo de inclinación entre 45º-60º. En teoría, CORMIX2 permite modelar este
último diseño mediante la opción (“Two nozzles: about 180º, same direction”), pero
en realidad, con esta opción, asume que el vertido se realiza a través de una única
boquilla, con inclinación vertical; hipótesis que invalida completamente los resultados
en el caso de efluentes de flotabilidad negativa. Por tanto, para realizar una
simulación rigurosa, se debe modelar cada lado del tramo difusor por separado,
como un difusor de una boquilla por elevador y chorros unilaterales, calculando el
caudal de modo que el caudal por boca sea el que corresponde al diseño original. En
caso de que existan corrientes que pudieran desviar el chorro y producir interacción
entre los chorros de distinto lado del difusor, esta simplificación no sería correcta.

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

- Dado el alto grado de sensibilidad al valor de la intensidad de las corrientes en el medio


receptor, es importante introducir valores realistas de corrientes en el medio receptor a la
profundidad de influencia del chorro. Además, se recomienda modelar con varias
velocidades (alrededor de la velocidad de la corriente) con el fin de asegurarnos de que el
resultado es acorde al flujo y que no nos encontramos en una de las zonas de transición
donde el modelo simula saltos bruscos de una clase de flujo a otra.
- Evitar el modelado con CORMIX2 de los casos en que los chorros impactan con la
superficie, dado que los resultados indican, en algunos casos, diluciones superiores a las
obtenidas para el mismo diseño pero un mayor calado, de modo que no exista interacción,
y dado que una vez que hay interacción los resultados son independientes del calado en el
medio receptor. Estos resultados se contradicen con lo obtenido experimentalmente.
- Evitar el modelado con CORMIX2 de chorros que interaccionan entre sí, dado que para
CORMIX2, una vez que existe interacción, los resultados son independientes de la
separación entre chorros, lo cual está en contradicción con los resultados experimentales.
- Respecto a la configuración de vertido, evitar el modelado de dispositivos de chorros
bilaterales (a distintos lados del difusor), o de dispositivo con más de una boquilla por
elevador, a la vista de las hipótesis simplificativas que el modelo asume (una boquilla y
    90º ) y que suponen un fuerte error en el modelado de chorros hiperdensos.

Tabla 6: Recomendaciones respecto al modelado de vertidos de salmuera mediante chorros múltiples en tramo
difusor con el subsistema CORMIX2

RECOMENDACIONES RESPECTO AL USO Y MANEJO DEL PROGRAMA

Para evitar incompatibilidad en los ángulos de entrada al modelo:    .

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

- Si el tramo difusor es paralelo a la costa (y a la corriente) y perpendicular al emisario, los


valores de distancia a introducir serán:
“Distance to one/other endpoint”: YB2=YB1.
 Si el vertido se realiza alejándose de la costa:   270º
 Si el vertido se realiza hacia la costa:   90º
En el caso de chorros perpendiculares al tramo difusor:   90º y   0º

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

9. ACCESO AL MODELO Y A LA INFORMACIÓN, PRECIO, INSTALACIÓN,


MANEJO, TIEMPO DE COMPUTACIÓN, ENTRADA DE DATOS Y
PRESENTACIÓN DE RESULTADOS.
Acceso al modelo y a la información: Documentos y trámite de solicitud de licencia en la
página web: www.cormix.com .
- Idioma: La información disponible (manuales, artículos, etc.), así como la interfaz del modelo
con el usuario está disponible únicamente en lengua inglesa.
- Precio: El precio de la licencia varía según el modelo se adquiera temporal o
definitivamente. El subsistema CORMIX2 se incluye en la versión más básica (v6.0 G), de
herramientas generales de CORMIX, cuya licencia tiene los siguientes precios (Mayo 2010):

ALQUILER DEL MODELO ($) COMPRA DEL MODELO


V6.0G: Herramientas generales, que
incluyen al CORMIX2 V6.0G: Herramientas generales, que
incluyen al CORMIX2
1 mes 3 meses 6meses 1 año

Uso académico - - - - -

Uso comercial 1.099 1.599 1.899 2.599 8.999

Tabla 8: Precios de la Licencia de herramientas generales, para el uso de CORMIX2

El resto de características son idénticas a las descritas en la ficha general de CORMIX. .

10. REFERENCIAS Y BIBLIOGRAFÍA RECOMENDADA.


- (1) Akar, P.J. and G.H. Jirka (1991). “CORMIX2: An Expert System for Hydrodynamic Mixing
Zone Analysis of Conventional and Toxic Submerged Multiport Diffuser Discharges", U.S.
Environmental Protection Agency, Office of Research and Development, Washington, DC.

- (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.

- Fernández-Torquemada, Y; Sánchez-Lisazo, J.L; (2006). “Effect of salinity on growth and


survival of Cymodocea nodosa ascherson and Zostera noltii Hornermann”. Biology Marine
Mediterranean 13, pp.46-47.

- 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.

- Palomar, P; Losada, I. (2008). ”Desalinización de agua marina en España: aspectos a


considerar en el diseño del sistema de vertido para protección del medio marino”. Revista de
Obras Públicas, nº 3486, pp. 37-52.

- Más información sobre el manejo y sistema de clasificación del subsistema CORMIX2:

- 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.

EJEMPLO DE APLICACIÓN DEL MODELO CORMIX2 A UN VERTIDO


PROTOTIPO DE SALMUERA AL MAR MEDITERRANEO. PASOS E INTERFAZ.
CASO DE APLICACIÓN:
Se desea modelar el comportamiento de un vertido de salmuera procedente de una nueva
planta desalinizadora proyectada en la costa Mediterránea. Se opta por una configuración de
descarga mediante tramo difusor de chorros múltiples, con elevadores de una boquilla y chorros
unilaterales, descargando con cierta incitación. La salmuera procede de un proceso de osmosis
inversa de tasa de conversión del 45%. Los valores de datos de entrada responden a un caso dentro
del rango realista de valores de diseño y de valores de las condiciones en el medio receptor.
METODOLOGIA PARA LA SIMULACION DEL VERTIDO MEDIANTE CORMIX2.
1. Recopilar los valores de los datos de entrada al modelo, que incluyen:
- Condiciones en el medio receptor, correspondientes a distintos escenarios, estacionales, en
el Mar Mediterráneo, y teniendo en cuenta el clima local.

CONDICIONES AMBIENTALE EN EL MEDIO RECEPTOR.


(entre paréntesis el nombre con que se definen en CORMIX2)
Profundidad de vertido, m. Medio receptor Mar
15
(Depth) (Bounded/Unbounded) (ilimitado)

Velocidad corriente en el medio


receptor, m/s. (Ambient velocity)
0.15
Densidad del medio receptor,
kg/m3 (correspondiente a 1027
Exceso de salinidad máximo en T=15ºC y S=37.5psu)
1.5
zona a proteger, psu

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:

PARÁMETROS DE DISEÑO DEL VERTIDO DE SALMUERA


(entre paréntesis el nombre con que se definen en CORMIX2)

Diámetro de boquilla, m. Velocidad de salida del chorro, m/s.


0.2 5
(Port diameter) (Effluent velocity)

Altura difusor, m. Densidad del efluente, kg/m3


1 1052
(Port height) (Effluent density)

Angulo inclinación descarga, Temperatura del efluente, ºC.


60º 15/24
(Vertical angle: THETA) (Effluent temp) (*)

Angulo de los chorros respecto a Concentración, exceso de salinidad


la corriente 0º respecto M Receptor, psu, 30.5
(Horizontal angle: SIGMA) (Concentration)

Longitud del difusor Configuración del dispositivo de vertido:tramo


150
(diffuser length) difusor lineal y paralelo a la costa (a 500m), con

133
MODELOS COMERCIALES: CORMIX2 de CORMIX

una boquilla por elevador, chorros unilaterales


Número de elevadores (boquillas) vertiendo en perpendicular al tramo difusor:
11
(Number of openings)
(GAMMA=90º, BETA=90º)
NOTA1: Se supone desalinización por Osmosis inversa con tasa de conversión, R= 45%.
Qproduccion ; C
Qsalmuera  Qo  Csalmuera  Co  A ; siendo C=concentración salina; CA=salinidad del agua
1 R 1 R
alimentación;

En el Anexo 2 de esta Ficha Descriptiva se ha incluido una tabla de valores realistas y


recomendados para los datos de entrada de CORMIX2 (ver Sección 4: “opciones del modelo”),
en el modelado de un vertido de salmuera en el Mar Mediterráneo

2. INTERFAZ DEL SUBSISTEMA CORMIX2.


Las imágenes que se muestran en las páginas siguientes corresponden a la versión más
actual (v.6) del sistema CORMIX
Se inicia en primer lugar el programa con doble click en el icono de CORMIX. La interfaz del
programa presenta el siguiente aspecto:

En la zona superior existen siete


pestañas para acceder a las distintas opciones
del modelo:
- PROJECT: proyecto.
- AMBIENT: medio receptor.
- EFFLUENT: efluente de salmuera.
- DISCHARGE: descarga o vertido.
- MIXING ZONE: Zona de mezcla.
- OUTPUT: Resultados.
- PROCESSING: Verificación y ejecución

3. PESTAÑA DE DEFINICIÓN DEL PROYECTO: “PROJECT”


Incluye los datos sobre el proyecto y la descripción del mismo.
Para comenzar a utilizar CORMIX puede cargarse un proyecto ya existente o crearse
directamente uno nuevo. Los proyectos en CORMIX2 se guardan con extensión “cmx”.
Para cargar un proyecto ya existente: Project: Load. Se rellena la información con nuestros
datos y se guarda con el nombre deseado: Project: save as.
Para un proyecto nuevo: Se rellena la casilla “Design case” y se guarda: Project: save as.

134
MODELOS COMERCIALES: CORMIX2 de CORMIX

4. PESTAÑA DE CARACTERÍSTICAS DEL EFLUENTE: “EFFLUENT”.


Incluye las características del efluente salmuera. Se rellenan las casillas con los datos de
entrada del modelo (ver la sección 4 de la ficha: “Opciones del modelo”, y el anexo 2 de valores
realistas y recomendados). Para simular un vertido de salmuera se puede seleccionar la opción de
contaminante conservativo (“Conservative pollutant”), o bien seleccionar directamente salmuera
(“brine”). Las opciones son muy similares. La interfaz de “EFFLUENT” es:

135
MODELOS COMERCIALES: CORMIX2 de CORMIX

5. PESTAÑA DE CONDICIONES EN EL MEDIO RECEPTOR: “AMBIENT”


Incluye las variables (condiciones) en el medio receptor.
Se completan las casillas con los datos de entrada (ver la sección 4 de la ficha: “Opciones del
modelo” y el Anexo 2, de valores realistas y recomendados). Siguiendo estos pasos para el caso
propuesto, la interfaz de “AMBIENT” resulta:

6. PESTAÑA DE DISEÑO DEL DISPOSITIVO DE VERTIDO: “DISCHARGE”.


Incluye los parámetros de diseño del sistema de vertido. En este caso de vertido mediante
chorro individual, lo primero es seleccionar la pestaña de “CORMIX2”, y a continuación ir rellenando
las casillas con los datos de entrada (ver la sección 4 de la ficha: “Opciones del modelo” y el anexo 2
de valores realistas y recomendados).
Siguiendo estos pasos para el caso propuesto, la interfaz de “DISCHARGE” resulta:

136
MODELOS COMERCIALES: CORMIX2 de CORMIX

7. PESTAÑA DE ZONA DE MEZCLA: “MIXING ZONE”


Incluye las características de la zona de mezcla. La salmuera no se considera un efluente
tóxico. La norma de calidad ambiental podrá establecer como el exceso de salinidad máxima
permitido en la zona a proteger. En la normativa vigente aún no se ha definido las dimensiones de la
zona de mezcla en el caso de vertidos de salmuera. Se indica la región de interés, que será el ámbito
de aplicación del CORMIX. La interfaz de “MIXING ZONE” resulta:

8. PESTAÑA DE RESULTADOS: “OUTPUT”


Incluye las opciones para generar los distintos tipos de informes de resultados (ver sección 4
de la ficha; “Opciones del modelo”). La interfaz de “OUTPUT” es la siguiente:

El bloque de opciones “output options” muestra los distintos ficheros de resultados,


destacando la opción de “prediction file” que ofrece los resultados numéricos de la simulación.
El bloque de opciones “Rulebase Browser Display Options” incluye las reglas de clasificación
que se han seguido para definir la categoría, clase de flujo, módulo de cálculo, etc.
El bloque de “coordinate system:corUCS” está disponible en las últimas versiones de
CORMIX y permite georreferenciar el vertido.

137
MODELOS COMERCIALES: CORMIX2 de CORMIX

9. PESTAÑA DE PROCESADO DE EJECUCION: “PROCESSING”.


Incluye las opciones de ejecución del modelo, que se pueden realizar de un solo paso con la
opción: “1,2,3&RUN”, o bien , paso a paso, del siguiente modo:
- Validar los datos de entrada (Validate input): Detecta posibles errores o inconsistencias.
- Cálculo de los parámetros fundamentales del modelo (Calculate parameters): Calcula los flujos
de caudal, de cantidad de movimiento y de flotabilidad, las escalas de longitud que utilizará
CORMIX2 y el número de Froude densimétrico.
- Clasificación del flujo (Classify flow)
- Ejecutas mediante “1,2,3&RUN”.
El programa generará archivos de resultados que pueden guardarse para ser legibles
mediante cualquier lector de textos. Para guardarlo se pulsa la casilla: Save, en la hoja de resultados.
Siguiendo estos pasos para el caso propuesto, la interfaz de “PROCESSING” resulta:

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.

RANGO DE VALORES REALISTAS EN LOS PARAMETROS DE ENTRADA,


PARA LA SIMULACION DEL VERTIDO DE SALMUERA AL MAR MEDITERRANEO.
DESALINIZACION POR ORMOSIS INVERSA. SUBSISTEMA CORMIX2
Es importante señalar que la horquilla de valores admitidos por CORMIX2 es más amplia,
pero en muchas ocasiones está fuera del rango de diseño de nuestras plantas.

VALORES REALISTAS Y RECOMENDADOS DE PARAMETROS DE ENTRADA EN SIMULACION DE


VERTIDOS DE SALMUERA AL MAR MEDITERRANEO (CORMIX2).

“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)

La que corresponda. 0 0.01-1m/s 1028- 1024kg/m³ No afecta al chorro


(*4) (*5) (*6) (*3) hasta impacto (*7)

“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)

Configuración del dispositivo de vertido (difusor)


(Diffuser configuration)

Chorros unilaterales: (unidireccional).


Misma dirección/ Abanico: Chorros bilaterales: (Alternating).
(Same direction/Fanned out) Misma dirección/Abanico,
Indiferente una y otra opción. (Same direction/ Fanned ou)t
Una boquilla por elevador
CORMIX2 realiza los mismos cálculos.
(Single nozzle per riser)
Theta Betta (*) Sigma (*)

SI: 0º a 90º Indiferente una y otra opción.


45º-60º (<75º) 90º CORMIX2 realiza los mismos cálculos
NO: 180º
(*10) (*11)
(*12)
Ángulo horizontal entre chorros menor 60º Ángulo entre chorros de unos 180º
Dos boquillas por
(Less than 60º). (About 180º).
elevador
Misma dirección/ abanico Misma dirección/Abanico

140
MODELOS COMERCIALES: CORMIX2 de CORMIX

(Two nozzles per riser) (Same direction/ Fanned out) (Same direction/ Fanned out)

Indiferente una u otra opción.


CORMIX2 realiza los mismos cálculos

Theta Beta Sigma Indiferente una u otra opción.


CORMIX2 realiza los mismos cálculos
SI: 0º a 90º
45º-60º (<75º) 90º
NO: 180º
(*10) (*11)
(*12)
Flujo de C. movimiento distintos de cero
Flujo de cantidad de movimiento nulo (Momentum flux non Zero).
Momentum flux zero. Misma dirección/ Abanico
Misma dirección /Abanico (Same direction/ Fanned out)
Varias boquillas por
elevador Same direction /Fanned out Indiferente una u otra opción.
(Several nozzles per CORMIX2 realiza los mismos cálculos
riser)
Theta Beta Sigma
Indiferente una u otra opción. SI: 0º a 90º
45º-60º (<75º) 90º
CORMIX2 realiza los mismos cálculos NO: 180º
(*10) (*11)
(*12)
“MIXING ZONE”
(ZONA DE MEZCLA)
Norma de calidad M. Receptor,
Zona de mezcla Intervalos hasta
exc. salinidad Zona de interés
(Mixing zone specified. representación
(Non toxic effluent. WQ (Region of Interest)
Distance/With/Area) (Intervals for display)
standard)

Exceso de salinidad permitido: Dimensiones de la zona


2.5psu, 2 psu, 1.5psu. de mezcla (donde no se La que corresponda 10-20
(*13) cumplen WQ estándar)

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.

(*3) Considerando el rango de densidades realistas en el Mar Mediterráneo en las


profundidades de vertido y tasas de conversión de la planta del 45%, y suponiendo (como
prácticamente ocurre en la realidad) que la temperatura de la salmuera es prácticamente la misma
que la del agua de alimentación, el siguiente cuadro indica el rango de densidades aproximadas.
Para una salinidad media del agua de alimentación de 37.5psu y una temperatura media de
21ºC, y R=45%, la densidad aproximada de la salmuera sería de 1049Kg/m³.

141
MODELOS COMERCIALES: CORMIX2 de CORMIX

AGUA DE ALIMENTACIÓN (OBRA


SALMUERA
DE TOMA)
Concentración
Salinidad Temperatura Densidad Temperatura Densidad
Salinidad
(psu) (ºC) (Kg/m³) (ºC) (Kg/m³)
(psu)
37 1027.4 67.3 1050.6
37.5 15 1027.8 68 15 1051.2
38 1028.1 OSMOSIS 69.1 1052
37 1026.6 INVERSA. 67.3 1049.7
37.5 18 1027 Tasa de 68 18 1050.2
38 1027.3 conversión 69.1 1051.1
37 1025.7 del 67.3 1048.8
37.5 21 1026 45% 68 21 1049.3
38 1026.5 69.1 1050.2
37 1024.8 67.3 1047.8
37.5 24 1025.2 68 24 1048.3
38 1025.6 69.1 1049.2
37 1023.8 67.3 1046.8
37.5 27 1024.2 68 27 1047.3
38 1024.6 69.1 1048.2

(*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.

(*10) Acorde a la opinión de la comunidad científica (6)

(*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.

RESULTADOS NUMÉRICOS DEL ANÁLISIS DE SENSIBILIDAD PARA EL MODELO CORMIX2 DE CORMIX.


En este anexo se presentan los resultados del análisis de sensibilidad del subsistema CORMIX2 a los valores de los datos de entrada al modelo. Las
conclusiones a este análisis se exponen en la sección 5 de esta ficha descriptiva.
Glosario de términos para los datos de entrada (características del efluente, parámetros de descarga y condiciones en el medio receptor:

-  A : densidad del medio receptor.

- H A : profundidad media de la columna de agua en la zona de vertido.


- U A ; velocidad de la corriente en el medio receptor.
- S A : salinidad del medio receptor
- TA : temperatura del medio receptor

-  (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.

-  (beta): ángulo horizontal de los chorros respecto al tramo difusor.

- o : densidad del efluente a la salida.

- ho : altura de la boquilla con respecto al fondo.

-  (theta): ángulo de inclinación del chorro en la descarga.

- U o : velocidad de salida del chorro.


- D : diámetro de la boquilla.
- To : temperatura del efluente.
- S o : concentración salina del efluente

- L: longitud del tramo difusor.


- s: separación entre boquillas.

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.

- Gravedad reducida: g ´ g  o  g  o   a . - Número de Froude densimétrico: F  Uo


a a
o
g o ´D

- 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

A) INFLUENCIA DEL ÁNGULO DE INCLINACIÓN DEL CHORRO (  ).

U A  0 ; H A =20m ;  A  1027.5 Kg/m³;   0º (co-flowing), ho  1 ;  o  1052 Kg/m³; U o  5m/s ; D  0.2 m; U viento  0 ;


Parámetros fijos: Parámetros fijos: :
RMZ=100m; L=150m;YB1=1000; YB2=1150; N=11;s=15;
  90º ;   90º
A.1) Variamos  , desde 0 a 90º (con incrementos de 15º).

RESULTADOS CON CORMIX2


VARIABLE
Nº Froude Punto de máxima altura del Punto de impacto del eje del
MODIFICADA CLASE DE
densimétrico eje del chorro chorro con el fondo
FLUJO Interacción
Boquilla/ranura Capa turbulenta
Chorros
 (º) Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

0º 23 / 215 MNU1 1 1 5.4 5 NO No ejecuta (Ua=0)

15º 23 / 215 MNU1 2 4.7 9.3 9 NO No ejecuta (Ua=0)

30º 23 / 215 MNU1 3.6 6 11.7 12.2 NO No ejecuta (Ua=0)

45º 23 / 215 MNU1 5.6 6.5 12 14.3 NO No ejecuta (Ua=0)

60º 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO No ejecuta (Ua=0)

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³).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura del Punto de impacto del eje
Nº de Froude
CLASE DE FLUJO eje del chorro del chorro con el fondo Interacción
densimétrico Capa turbulenta
de chorros
O (Kg/m )
3 Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

1047 25 / 232 MNU1 7.9 6.9 11 16.5 NO No ejecuta (Ua=0)

1049 23.8 / 222 MNU1 7.6 6.6 10.5 15.9 NO No ejecuta (Ua=0)

1051 22.9 / 213 MNU1 7.3 6.4 10 15.3 NO No ejecuta (Ua=0)

1053 22 / 205 MNU1 7.1 6.2 9.7 14.8 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³).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura Punto de impacto del eje
Nº de Froude
CLASE DE FLUJO del eje del chorro del chorro con el fondo Interacción
densimétrico Capa turbulenta
de chorros
O (Kg/m3)
Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

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)

1053 22.7 / 211 MNU1 7.3 6.3 10 15 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³).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura Punto de impacto del eje
Nº de Froude
CLASE DE FLUJO del eje del chorro del chorro con el fondo Interacción
densimétrico Capa turbulenta
de chorros
O 3
(Kg/m )
Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

1047 26 / 245 MNU1 8.2 7.2 11.5 17.3 NO No ejecuta (Ua=0)

1049 25 / 232 MNU1 7.9 6.9 11 16.5 NO No ejecuta (Ua=0)

1051 24 / 222 MNU1 7.6 6.6 10.5 15.9 NO No ejecuta (Ua=0)

1053 23 / 213 MNU1 7.3 6.4 10 15.3 NO No ejecuta (Ua=0)

148
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

C) INFLUENCIA DE LA ESTRATIFICACIÓN EN LA COLUMNA DE AGUA.


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º
C.1.) Representando una situación otoñal/invernal en el Mar Mediterráneo, sin presencia de Termoclina y una temperatura uniforme en la columna de agua
del medio receptor, con valores de 15º, 18º y 21ºC. Se considera la misma temperatura en el efluente y el medio receptor

RESULTADOS CON CORMIX2


VARIABLES MODIFICADAS
Nº de Froude CLASE Punto de máxima altura del Punto de impacto del eje del
densimétrico DE eje del chorro chorro con el fondo
Interacción
Chorro/ranura FLUJO Capa turbulenta
T (ºC) A O Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)
chorros
3 3
(Kg/m ) (Kg/m )

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³;

RESULTADOS CON CORMIX2


VARIABLES MODIFICADAS
Nº de Froude CLASE DE
densimétrico FLUJO
 A _ sup  A _ fondo Valores Simulación
3
(Kg/m ) (Kg/m3)

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³;

RESULTADOS CON CORMIX2


VARIABLES MODIFICADA
Punto de máxima altura del Punto de impacto del eje del
Nº de Froude CLASE DE eje del chorro chorro con el fondo
densimétrico FLUJO Interacción
Capa turbulenta
 A _ sup  A _ fondo
Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)
chorros
3
(Kg/m ) (Kg/m3)

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) INFLUENCIA DE LA VELOCIDAD DE LA CORRIENTE ( U A ) Y DEL ÁNGULO HORIZONTAL ENTRE LOS CHORROS Y LA


CORRIENTE EN EL MEDIO RECEPTOR (  ).
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;s=15;   90º ;   90º

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).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura del Punto de impacto del eje
Nº de Froude CLASE DE Capa turbulenta
eje del chorro del chorro con el fondo Interacción de
densimétrico FLUJO
chorros
U A (m/s) Zmax (m) Seje (Zmax) Xi (m) Seje (Xi) X S

0 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO No ejecuta (Ua=0)

0.05 23 / 215 MNU1 6.6 6.6 11.5 19 NO Intrusion

0.10 23 / 215 MNU2 6.4 8 13 24.5 NO 17 35

0.15 23 / 215 MNU2 6.15 9.7 15.1 30.9 NO 19.4 43.7

0.30 23 / 215 MNU2 5.1 11.7 21.3 44.9 NO 25 63.5

0.45 23 / 215 MNU2 4.6 14.3 27 51.3 NO 30.9 72.6

0.6 23 / 215 MNU2 4.1 15.6 33 54.3 NO 35.8 76.8

0.75 23 / 215 MNU2 3.6 13.5 37 55.5 NO 40 78.4

0.9 23 / 215 MNU2 3.3 13.9 41.4 55.6 NO 44.3 78.6

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

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura del Punto de impacto del eje del
Nº de Froude CLASE DE Capa turbulenta
eje del chorro chorro con el fondo Interacción
densimétrico FLUJO
de chorros
UA (m/s)
Zmax (m) Seje (Zmax) Xi (m) Seje (Xi) Y/X S

0 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO No ejecuta (Ua=0)

0.05 23 / 215 MNU1 7.3 6.7 10/ 2 16.4 NO Intrusión

0.1 23 / 215 MNU2 7.2 7.8 9.4/4.3 20.5 NO 9.4/80 29

0.15 23 / 215 MNU2 6.6 7.1 8.6/7.1 29.9 NO 8.6/82 42

0.30 23 / 215 MNU2 6 12.4 6.8/17 59 NO 6.8/92 84

0.45 23 / 215 MNU2 5.7 20.3 5.7/27 63 NO 5.7/101 103

0.6 23 / 215 MNU2 5.3 24.6 5.2/36.4 80.9 NO 5.2/111 114

0.75 23 / 215 MNU2 4.9 27.3 4.9/45.6 85.9 NO 4.9/120 121

0.9 23 / 215 MNU2 4.5 29.1 4.5/54 89 NO 4.5/ 129 126

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).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto máxima altura eje Punto de impacto del eje
Nº de Froude CLASE DE
del chorro del chorro con el fondo Interacción de
densimétrico FLUJO Capa turbulenta
Seje chorros
UA (m/s)
Zmax (m)
(Zmax)
Y/Xi (m) Seje (Xi)

0 23 / 215 MNU1 9 7.9 12.7 18.9 NO No ejecuta (Ua=0)

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”.

E) INFLUENCIA DEL CALADO DEL MEDIO RECEPTOR ( H A ).

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º

E.1) Variamos el calado en el medio receptor: H A entre 35m y 5m (incrementos de 5m).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura del Punto de impacto del eje del
Nº de Froude CLASE DE
eje del chorro chorro con el fondo
densimétrico FLUJO Interacción de chorros Capa turbulenta
H A (m) Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

35 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO

30 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO

25 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO

20 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO

15 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO

10 23 / 215 MNU1 7.4 6.4 10.2 15.4 NO No ejecuta (Ua=0)

Z(*) X(*) S

7 23 / 215 MNU7 1.33 10 11.1 NO

5 23 / 215 MNU7 1.2 10 11 NO

3.5 23 / 215 MNU7 1.1 10 10.9 NO

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

F) INFLUENCIA DE LA ALTURA DE LA BOQUILLA CON RESPECTO AL FONDO ( ho ).

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).

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Punto de máxima altura del Punto de impacto del eje del
Nº de Froude CLASE DE
eje del chorro chorro con el fondo Interacción de
densimétrico FLUJO Capa turbulenta
chorros
ho (m) Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

0 23 / 215 MNU1 6.4 6.4 9.8 13.9 NO

0.5 23 / 215 MNU1 6.7 6 10 14.6 NO


No ejecuta (Ua=0)
1 23 / 215 MNU1 7.3 6.1 10.2 15.4 NO

1.5 23 / 215 MNU1 7.9 6.2 10.3 16.2 NO

154
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

G) INFLUENCIA DE LA SEPARACIÓN ENTRE BOQUILLAS (S).

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

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Nº de Froude CLASE DE Punto de máxima altura del eje Punto de impacto del eje del
densimétrico FLUJO del chorro chorro con el fondo Interacción
Capa turbulenta
de chorros
S (m) Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

20 (L=200m) 23 / 249 MNU1 7.3 6.4 10.2 15.4 NO

15 (L=150m= 23 /215 MNU1 7.3 6.4 10.2 15.4 NO

10 (L=100m) 23 / 176 MNU1 7.3 6.4 10.2 15.4 NO

5 (L=50m) 23 /124 MNU1 7.3 6.4 10.2 15.4 NO

3 (l=30m) 23 /96 MNU1 7.3 6.4 10.2 15.4 NO No ejecuta (Ua=0)

2.5 (l=25m) 23 /88 MNU1 7.5 5.7 10.3 17.8 SI (X=8.2)

2 (L=20m) 23 / 78 MNU1 7.5 5.7 10.7 14.7 SI (X=6)

1 (L=10m= El programa no se ejecuta porque L<ha

155
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

H) INFLUENCIA DEL NÚMERO DE BOQUILLAS (ELEVADORES), (N).

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.

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Nº de Froude CLASE DE Punto de máxima altura del eje Punto de impacto del eje del
Capa turbulenta
densimétrico FLUJO del chorro chorro con el fondo Interacción de
chorros
N Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

26 (s=6m) 23 /140 MNU1 7.4 6.4 10.2 15.4 NO

21 (s=7.5m) 23 / 155 MNU1 7.4 6.4 10.2 15.4 NO


No ejecuta (Ua=0)
16 (s=10m) 23 / 178 MNU1 7.4 6.4 10.2 15.4 NO

11 (s=15m) 23 /215 MNU1 7.4 6.4 10.2 15.4 NO

6 (s=30) 23 / 291 MNU1 7.4 6.4 10.2 15.4 NO

156
PROYECTO I+D+i “MEDVSA”
MODELOS COMERCIALES: CORMIX2 de CORMIX

I) INFLUENCIA DE LA CONFIGURACION DEL DISPOSITIVO DE VERTIDO (TRAMO DIFUSOR).


Parámetros fijos: : U A  0 ; H A =20m ;  A  1027.5 Kg/m³;   0º (co-flowing), ho  1 ;  o  1052 Kg/m³; U o  5m/s ; D  0.2 m; U viento  0 ; RMZ=100m;
L=150m;YB1=1000; YB2=1150; N=11;s=15;

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.

RESULTADOS CON CORMIX2


VARIABLE MODIFICADA
Punto de máxima altura del Punto de impacto del eje del
Nº de Froude CLASE DE eje del chorro chorro con el fondo
densimétrico FLUJO Interacción
Configuración, ángulo de Cálculos
chorros
inclinación  Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

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)

Single nozzle: unidirectional:


  90º   90 º 7.4 6.4 10.2 15.4 CORJET
Same direction   60º
Single nozzle: unidirectional:
  90º   90º 7.4 6.4 10.2 15.4 CORJET
Fanned out   60º

Two nozzles:less than 60 degrees


  90º ;   90 º 7.4 6.4 10.2 15.4 CORJET
Same direction/ fanned out

X(m) S BH BV

Single nozzle: unidirectional Equivalent slot


Same direction/ fanned out ,
  90º 1.5 9.7 75.5 0.45
diffuser hyp.

Single nozzle: unidirectional Equivalent slot


Same direction/ fanned out ,
  90º 1.5 9.7 75.5 0.45
diffuser hyp.

Single nozzle: alternating Equivalent slot


- 1.5 9.7 75.5 0.45
Same direction/ fanned out diffuser hyp.

Two nozzles:less than 60 degre: Equivalent slot


Same direction/fanned out   90º - 1.87 7.7 75.5 0.6
diffuser hyp.

Two nozzles: about 180 degrees Equivalent slot


- 1.87 7.7 75.5 0.6
Same direction/ fanned out diffuser hyp.

Several nozzles: net horizontal


momentum flux zero (3ports) Equivalent slot
- 2.1 6.7 75.5 0.65
diffuser hyp.
Same direction// fanned out

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

RESULTADOS CON CORMIX2


VARIABLE MODIFICADA Punto de máxima altura del eje Punto de impacto del eje del
Nº de Froude CLASE DE
densimétrico FLUJO del chorro chorro con el fondo Interacción
Cálculos
de chorros
Configuración Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

Fanned out:,   30º 23 / 215 MNU1 3.6 6 11.7 12.2 NO CORJET

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

Hipótesis de vertido en ranura


Same direction, 23 / 215 (*) MNU1 1.5 9.7 75 0.45
equivalente

Hipótesis de vertido en ranura


Fanned out 23 /215 (*) MNU1 1.5 9.7 75.5 0.45
equivalente

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

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Nº de Froude Punto de máxima altura Punto de impacto del eje del
CLASE DE FLUJO del eje chorro chorro con el fondo Interacción de
densimétrico Cálculos
chorros
Configuración Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

Same direction Hipótesis de vertido mediante ranura.


23 / 152 MNU7 X =75; BV=20; BH=37; S=43.8 Ocupa toda columna(con 35º aplica
  30º CORJET)
Same direction
23 /152 MNU1 7.4 6.4 10.2 15.4 NO CORJET
  60º
X(m) S BH BV

Same direction Hip. vertido


23 /152 MNU1 1.87 7.7 75.5 0.6
  90º ranura equiva

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º .

RESULTADOS CON CORMIX2


VARIABLE
MODIFICADA Nº de Froude CLASE DE Punto de máxima altura del Punto de impacto del eje del chorro
densimétrico FLUJO eje del chorro con el fondo Interacción de
Cálculos
chorros
Configuración Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

Fanned out Hipótesis de vertido mediante ranura.


23 / 152 MNU7 X =75; BV=20; BH=37; S=43.8 Ocupa toda columna(con 35º aplica
  30º CORJET)
Fanned out
23/152 MNU1 7.4 6.4 10.2 15.4 NO CORJET
  60º
X(m) S BH BV

Fanned out Hipótesis de vertido mediante ranura


23/152 MNU1 1.87 7.7 75.5 0.6
  90º equivalente

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”).

VARIABLE MODIFICADA RESULTADOS CON CORMIX2


Nº de Froude CLASE DE
Configuración densimétrico FLUJO X(m) S BH BV Cálculos

hipótesis de vertido en
About 180º: Same direction 23 / 152 (**)MNU1 1.87 7.7 75.5 0.6
ranura equivalente

Hipótesis de vertido en ranura


Abou 180º: Fanned out 23 /152 (**)MNU1 1.87 7.7 75.5 0.6
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”) (*)

VARIABLE MODIFICADA RESULTADOS CON CORMIX2


Nº de Froude
CLASE DE FLUJO
densimétrico
Configuración X(m) S BH BV Cálculos

Same direction // Fanned out Hipótesis de vertido en ranura


23 /96 (**) MNU1 2.5 5.7 75.6 0.76
(5 bocas/difusor. N=55) equivalente

Same direction // Fanned out Hipótesis de vertido en ranura


23 / 124 MNU1 2.1 6.7 75.5 0.65
(3 bocas/difusor. N=33) equivalente

ºSame direction // Fanned out Hipótesis de vertido en ranura


23 / 102 MNU1 2.4 5.9 75.7 0.7
(7 bocas/difusor. N=77) equivalente

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”) (*)

RESULTADOS CON CORMIX2


VARIABLE MODIFICADA Punto de máxima altura del eje Punto de impacto del eje
Nº de Froude CLASE DE
densimétrico FLUJO del chorro del chorro con el fondo Interacción Módulo de
de chorros cálculo
Configuración Zmax (m) Seje (Zmax) Xi (m) Seje (Xi)

Same direction // fanned out


MNU1 3.6 6 11.7 12.2 NO CORJET
  30º
Same direction // fanned out
MNU1 7.4 6.4 10.2 15.4 NO CORJET
  60º
X(m) S BH BV

Same direction // fanned out Hipótesis de vertido en ranura


MNU1 1.65 8.7 75.5 0.5
  90º equivalente

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

También podría gustarte