TFM - Daniel Valero PDF

Modelacin hidrulica de flujos
multifase en grandes presas.

Aplicacin al diseo de cuencos
de amortiguacin.
Tutor:
Autor:
Rafael Garca Bartual
Daniel Valero Huerta
Tesina Final de Mster, tipo B

Mster Universitario en Ingeniera Hidrulica y Medio Ambiente
Valencia, Abril de 2015
No tena miedo a las dicultades: lo que la asustaba era la obligacin de tener
que escoger un camino. Escoger un camino signicaba abandonar otros."
Paulo Coelho
vii
Agradecimientos
Son muchas las personas a las que me gustara agradecer su apoyo durante estos
ltimos aos, no solo durante la realizacin del Mster si no tambin la carrera. A
ellos les podra escribir un documento tan largo como el presente y an me faltara
espacio. Sin embargo, intentar ser mucho ms concreto.
Me gustara agradecer muy sinceramente a mi tutor, el Prof. Rafael Garca Bartual, su apoyo y su conanza que en todo momento me ha brindado. Ha sido para m
siempre un ejemplo de motivacin continua y de brillantez. Su gran conocimientos y
perspectiva de todo en general y en particular de la hidrulica y de la mecnica de
uidos va ms all de donde puede llegar mi todava joven entendimiento. Trabajar
con l ha sido sin ninguna duda una de las ms graticantes experiencias que he
vivido. Quiero agradecer tambin al resto de profesores de la unidad docente Hidrulica e Hidrologa del Departamento de Ingeniera Hidrulica y Medio Ambiente
(DIHMA) por prestarme su tiempo cuando lo he necesitado, que no han sido pocas
las veces. Empezando por el Prof. Juan Marco con el que he tenido la oportunidad
de compartir charlas interesantsimas sobre ujos multifase y mecnica de uidos en
general. No quiero olvidarme tampoco de Francisco Valls, Ignacio Andrs, Miguel
ngel Eguibar y Eduardo Albentosa con quien he podido trabajar bajo el paraguas
del proyecto EMULSIONA. Igualmente, agradecer a mparo Lpez su buen juicio
y el trato cercano con el que siempre me ha tratado. Todos ellos siempre han sabido
ofrecerme un buen consejo y me han ayudado mucho a crecer en lo profesional y en
lo personal.
No menos importante ha sido el apoyo de mis compaeros de laboratorio, Arnau
y Beatriz, con quien he compartido ideas fantsticas en todos en esos innumerables
cafs. Con ellos he aprendido mucha ms hidrulica de lo que hubiera imaginado. A
Ximo y Carlos, tcnicos capaces y resolutivos. Omar y Pedro, quienes me sufrieron
durante muchas horas. A Mario, quien siempre tiene un punto de vista muy prctico
e interesante; Ismael, Jess, Alicia y Chiara, quien me contagiaron su entusiasmo
desde mi primera beca. A mis amigos de ICCP, especialmente a Pepe y Agustn, a
quien les espera un futuro brillante y lleno de xitos.
Tambin quiero agradecer a Scherezade su innito apoyo, comprensin y dedicacin. Ha sacado de m siempre lo mejor, con su gracia natural y su buen hacer.
Como en el cuento, ha estado a mi lado ms de mil y una noches. Su apoyo ha sido
de gran ayuda. Nunca le estar lo sucientemente agradecido. A mi familia quien
son los que ms tiempo me han sufrido sin ninguna duda. Especialmente a mi padre
Jos Mara de quien he aprendido el gusto por el trabajo bien hecho, mi madre
Mari Carmen quien me ha enseado a ser valiente y no temer nunca al fracaso; y
mi hermano Eduardo que es mi hermano favorito y no deja de asombrarme cada da.
viii
ndice General
ndice General
ix
ndice de Figuras
xiii
ndice de Tablas
xv
Notacin
xvii
1. Introduccin
1.1.
Motivacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.2.
Objetivos
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2. Estado del Arte

2.1.
2.2.
Resea histrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.1.
Las presas en Espaa y en el mundo
2.1.2.
Dinmica de Fluidos Computacional (CFD)
. . . . . . . . . . . . . .
. . . . . . . . . .
11
Fundamentos fsicos . . . . . . . . . . . . . . . . . . . . . . . . . . . .
20
2.2.1.
20
Ecuaciones de gobierno . . . . . . . . . . . . . . . . . . . . . .
2.2.1.1.
Conservacin de la masa . . . . . . . . . . . . . . . .
21
2.2.1.2.
Conservacin de la cantidad de movimiento
. . . . .
22
2.2.1.3.
Sistema de ecuaciones de Navier-Stokes . . . . . . . .
22
2.2.2.
Naturaleza del trmino de presiones . . . . . . . . . . . . . . .
24
2.2.3.
La ecuacin de la vorticidad . . . . . . . . . . . . . . . . . . .
25
2.2.4.
Ecuaciones de gobierno para ujos multifase
26
2.2.5.
. . . . . . . . . .
. . . . . . . . .
29
2.2.5.1.
Naturaleza de las ecuaciones de Navier-Stokes
Comportamiento elptico . . . . . . . . . . . . . . . .
32
2.2.5.2.
Comportamiento parablico . . . . . . . . . . . . . .
33
2.2.5.3.
Comportamiento hiperblico . . . . . . . . . . . . . .
35
2.2.5.4.
Clasicacin de las ecuaciones del ujo . . . . . . . .
36
ix
2.3.
2.4.
2.2.6.
Ecuaciones de Reynolds o del ujo medio . . . . . . . . . . . .
37
2.2.7.
Problema del cierre . . . . . . . . . . . . . . . . . . . . . . . .
40
2.2.8.
Turbulencia
. . . . . . . . . . . . . . . . . . . . . . . . . . . .
41
Hidrulica de grandes presas . . . . . . . . . . . . . . . . . . . . . . .
43
2.3.1.
Aliviadero . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
45
2.3.1.1.
Supercie libre
45
2.3.1.2.
Perl de velocidades
. . . . . . . . . . . . . . . . . .
46
2.3.1.3.
Turbulencia . . . . . . . . . . . . . . . . . . . . . . .
50
2.3.1.4.
Determinacin del punto de inicio de la aireacin
. .
50
2.3.1.5.
Concentracin de equilibrio
. . . . . . . . . . . . . .
54
2.3.1.6.
Reduccin de friccin por el efecto de la aireacin . .
58
Fundamentos numricos
. . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . .
61
2.4.1.
El Mtodo de los Volmenes Finitos . . . . . . . . . . . . . . .
61
2.4.2.
Grados de aproximacin a la resolucin 3D de las ecuaciones

del ujo
2.4.3.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Modelado de la turbulencia
63
. . . . . . . . . . . . . . . . . . .
64
2.4.3.1.
Hiptesis Boussinesq o de viscosidad turbulenta . . .
64
2.4.3.2.
La ecuacin de la energa cintica turbulenta . . . . .
66
2.4.3.3.
Modelos de 0 ecuaciones . . . . . . . . . . . . . . . .
70
2.4.3.3.1.
Modelo de van Driest . . . . . . . . . . . . .
70
2.4.3.3.2.
Modelo Cebeci/Smith
71
2.4.3.3.3.
Modelo de Baldwin y Lomax
2.4.3.4.
. . . . . . . . . . . .
. . . . . . . .
72
Modelos de 1 ecuacin . . . . . . . . . . . . . . . . .
74
2.4.3.4.1.
Modelo de Baldwin-Barth
. . . . . . . . . .
75
2.4.3.4.2.
Modelo de Spalart-Allmaras . . . . . . . . .
76
Modelos de 2 ecuaciones . . . . . . . . . . . . . . . .
78
2.4.3.5.
2.4.3.5.1.
Modelo
2.4.3.5.2.
Modelo

. . . . . . . . . . . . . . . . .
80
RNG . . . . . . . . . . . . . .
81
. . . . . . . . . . . . . . . . . . .
82
2.4.4.
Modelado de la lmina libre
2.4.5.
Modelado de la reaireacin turbulenta
. . . . . . . . . . . . .
3. Calibracin del modelo de reaireacin turbulenta
86
91
3.1.
Introduccin y objetivos
. . . . . . . . . . . . . . . . . . . . . . . . .
91
3.2.
Descripcin general . . . . . . . . . . . . . . . . . . . . . . . . . . . .
92
3.3.
Opciones numricas . . . . . . . . . . . . . . . . . . . . . . . . . . . .
93
3.3.1.
Generalidades . . . . . . . . . . . . . . . . . . . . . . . . . . .
93
3.3.2.
Mallado
94
3.3.3.
Modelo de turbulencia
. . . . . . . . . . . . . . . . . . . . . .
96
3.3.4.
Modelo de lmina libre . . . . . . . . . . . . . . . . . . . . . .
97
3.4.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Resultados de la calibracin
. . . . . . . . . . . . . . . . . . . . . . .
97
xi
4. Validacin del modelo de reaireacin turbulenta
103
4.1.
Introduccin y objetivos
4.2.
Descripcin general . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
4.3.
. . . . . . . . . . . . . . . . . . . . . . . . . 103
Resultados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
4.3.1.
Resultados obtenidos directamente en el proceso de calibracin 104
4.3.2.
Opciones numricas . . . . . . . . . . . . . . . . . . . . . . . . 106
4.3.3.
Resultados obtenidos en el proceso de validacin . . . . . . . . 106
5. Aplicacin al diseo de disipadores de energa
109
5.1.
Objetivos
5.2.
Descripcin general . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109

5.2.1.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
Geometra, condiciones de contorno y opciones numricas . . . 111
5.3.
Estructura del ujo en el cuenco de amortiguacin . . . . . . . . . . . 113
5.4.
Resultados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.4.1.
Condiciones hidrulicas de diseo . . . . . . . . . . . . . . . . 115
5.4.2.
Condiciones hidrulicas desfavorables . . . . . . . . . . . . . . 118
6. Conclusiones
121
Bibliografa
124
xii
ndice de Figuras
1.1.
Flujo supercrtico reproducido en laboratorio (Pster, 2008). . . . . .
1.2.
Flujo supercrtico emulsionado sobre aliviadero y aireadores (presa de

Alicura, Argentina). Foto tomada por el profesor Dr. H.-E. Minor. . .
2.1.
2
3
Simulacin 3D de la regin no aireada sobre un perl WES (Ho et al.,

2001).
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
19
2.2.
Diferentes secciones cnicas en funcin del valor del discriminante. . .
32
2.3.
Dominio de dependencia de las ecuaciones elpticas. . . . . . . . . . .
33
2.4.
Dominio de dependencia de las ecuaciones parablicas.
. . . . . . . .
34
2.5.
Dominio de dependencia de las ecuaciones hiperblicas. . . . . . . . .
35
2.6.
Velocidad en un ujo turbulento y turbulencia instantnea en el punto
tras el cilindro. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
38
2.7.
Esquema aliviadero. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
45
2.8.
Clculo de la velocidad en un aliviadero segn USBR (Peterka, 1958).
47
2.9.
Campo de velocidades en el entorno de la cresta para
2.10. Campo de velocidades en el entorno de la cresta para

= 0,5.
= 1,0.
= 1,5.
= 2,0.
. . . .
48
. . . .
48
. . . .
49
. . . .
49
2.13. Determinacin del punto crtico de aireacin: interseccin de la capa

lmite con la lmina libre.
. . . . . . . . . . . . . . . . . . . . . . . .
52
2.14. Visualizacin de una capa lmite para ujo con Mach = 2.82. El ujo
se produce de derecha a izquierda. Adaptado de (Smith and Smits,
1995).
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.15. Estructura del ujo supercrtico aireado.
. . . . . . . . . . . . . . . .
54
55
2.16. a) Concentracin media b) Nomenclatura grcos c) Concentracin

media y de fondo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
55
2.17. Perl de concentraciones terico y medido en prototipo. . . . . . . . .
56
2.18. Perl de velocidades terico y medido en prototipo.
57
xiii
. . . . . . . . . .
NDICE DE FIGURAS
xiv
2.19. Evolucin de las concentraciones en un aliviadero. . . . . . . . . . . .
58
2.20. Reduccin de friccin por el efecto del aire en el ujo. . . . . . . . . .
59
2.21. Valores calculados con la ecuacin 2.89 frente a valores observados.
61
. . . . . . . . . . . . . .
83
2.23. Rotura de presa. Resultados del modelo fsico medidos con un PIV.
0,5
Velocidad adimensionalizada U = u/(4gd0 ) , con d0 = 0,40m y T =
0,5
t(gd1
= a) 0.50, b) 0.99, c) 1.98, e) 2.48 y f ) 2.97. Adaptado de
0 )
(Oertel and Bung, 2012). . . . . . . . . . . . . . . . . . . . . . . . . .
84
2.22. Representacin tpica de una lmina libre mediante el VOF. En el

centro de las celdas se muestra el valor de
F.
2.24. Rotura de presa. Resultados del modelo numrico empleando un m0,5

todo VOF. Velocidad adimensionalizada U = u/(4gd0 ) , con d0 =
1 0,5
0,40m y T = t(gd0 ) = a) 0.50, b) 0.99, c) 1.98, e) 2.48 y f ) 2.97.
Adaptado de (Oertel and Bung, 2012).
3.1.
. . . . . . . . . . . . . . . . .
85
Crecimiento de la capa lmite e interseccin con la lmina libre segn

el mtodo de (Wood et al., 1983; Wood, 1991). . . . . . . . . . . . . .
92
3.2.
Perles de velocidad sobre la cresta del aliviadero y en la seccin nal. 95
3.3.
Valores caractersticos de la energa cintica turbulenta sobre la cresta

del aliviadero y en la seccin nal. . . . . . . . . . . . . . . . . . . . .
3.4.
el caso
3.5.
4.1.
x = 8
cm.
. . . . . . . . . . . . . . . . . . . . . . . . . . . .
Error obtenido en la concentracin media para el caso
x = 8
cm. . .
98
99
Perles de concentracin de aire obtenidos con los valores ptimos de

los parmetros en funcin del tamao de celda.
4.2.
96
Error obtenido en la determinacin del punto crtico de aireacin para
. . . . . . . . . . . . 104
Perles de concentracin de aire descompuesto en la fraccin emulsionada y atrapada. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
4.3.
Resultados del proceso de validacin.
5.1.
Dimensiones del cuenco de amortiguacin y caractersticas del ujo

de entrada.
. . . . . . . . . . . . . . . . . . 107
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
hcb = 0. .
5.2.
Geometra en 3D para el caso
5.3.
Tamaos de diente de pie de presa ensayados.
5.4.
Estructura del ujo dentro del cuenco de amortiguacin.
5.5.
Cada de la velocidad mxima dentro del cuenco.
5.6.
Regin de cortante o shear region. . . . . . . . . . . . . . . . . . . . . 118
5.7.
. . . . . . . . . . . . . . . . . 112
. . . . . . . . . . . . . 112
. . . . . . . 114
. . . . . . . . . . . 117
Anlsis transitorio del funcionamiento hidrulico del cuenco de amortiguacin.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
ndice de Tablas
2.1.
Presas en Espaa hasta el s. XIX, adaptado de (Vallarino, 1991). . . .
2.2.
Evolucin histrica de las presas en Espaa.
2.3.
Estadstica presas en Espaa segn su altura.
. . . . . . . . . . . . . .
. . . . . . . . . . . . .
11
2.4.
Clasicacin de los principales tipos de ujo. . . . . . . . . . . . . . .
36
3.1.
Concentraciones de equilibrio en el ujo completamente desarrollado

segn los datos de (Straub and Anderson, 1958). . . . . . . . . . . . .
3.2.
3.3.
93
Valores ptimos de los parmetros del modelo de reaireacin turbulenta.100

Valores ptimos de
Cair
forzando la tensin supercial a la fsicamente
basada. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
4.1.
Concentraciones de equilibrio obtenidas con FLOW-3D, propuestas

por Chanson, 1996 (*) y propuestas por Hager, 1991 (**). Se muestran
tambin los errores relativos a cada uno de los estudios experimentales.106
5.1.
Resumen de las caractersticas de las mallas empleadas en cada una

de las simulaciones de este captulo. . . . . . . . . . . . . . . . . . . . 113
5.2.
Resultados obtenidos para el parmetro de sumergencia.
xv
. . . . . . . 116
xvi
NDICE DE TABLAS
Notacin
Notacin
Alfabeto latino
AS
rea de la lmina libre en una celda.
A+
Constante del modelo de van Driest y del modelo de Baldwin y Lomax.
C0,5
Concentracin para la cual la friccin efectiva se reduce a la mitad.
Cair
Constante de calibracin del modelo de reaireacin turbulenta.
Ce
Concentracin de equilibrio.
Cent
Concentracin de aire emulsionado.
Ce
Concentracin de aire emulsionado para el ujo completamente desarrollado.
Cmean
Co
Crough
Ctotal
Concentracin de media en el lete uido.

Constante del modelo de van Driest y del modelo de Baldwin y Lomax.
Concentracin de aire atrapado.
Concentracin de aire atrapado y emulsionado.
C1
Constante de la ecuacin de transporte para
.
C2
Constante de la ecuacin de transporte para
.
Constante de proporcionalidad para la viscosidad turbulenta.
CFD Computational Fluid Dynamics.

CFL Nmero de Courant.
CHJ Classic Hydraulic Jump.
CWJ Classic Wall Jet.
DNS Direct Numerical Simulation.
EDO Ecuacin Diferencial Ordinaria.
EDP Ecuacin en Derivadas Parciales.
EVM Eddy Viscosity Models.
F~c
Vector de ujos convectivos.
xvii
Notacin
Fr
Nmero de Froude.
F~v
Vector de ujos difusivos.
Fluid-in-Cell.
Factor de friccin de Darcy.
fe
Factor de friccin efectivo.
f~e
Vector de fuerzas externas.
gi
Proyeccin de la gravedad sobre el eje
gn
Proyeccin de la gravedad sobre la normal a la lmina libre.
HD
Hs
I
i.
Altura del embalse para el caudal de diseo.

Diferencia de alturas entre un punto a una distancia
x y el nivel del aliviadero.
Matriz identidad.
ICOLD International Commission on Large Dams.
J
K
k
ks
L
Tensor de inercia.
Constante de calibracin del modelo de reaireacin turbulenta.
Longitud de onda turbulenta.
Rugosidad del aliviadero.
Escala longitudinal.
LES Large-Eddy Simulation.
LT
L
l
Escala longitudinal turbulenta.

Subespacio de restricciones.
Longitud de onda caracterstica del torbellino de valor
2/k .
li
Longitud de mezcla turbulenta en la regin interna de la capa lmite.
lo
Longitud de mezcla turbulenta en la regin externa de la capa lmite.
Flujo msico a travs de una supercie.
Ma
Nmero de Mach.
xviii
Notacin
MAC Marker-and-Cell.
n
N ()
Rugosidad de Manning.
Operador de Navier-Stokes.
Pd
Presin de las fuerzas estabilizadoras.
Pt
Presin de las fuerzas desestabilizadoras.
PIC Particle-in-Cell.
Presin.
~
Q
Vector de trminos fuente.
Radio de curvatura de las ondulaciones de la lmina libre.
RANS Reynolds Averaged Navier-Stokes.
Re
Nmero de Reynolds.
RSM Reynolds Stress Models
Supercie de separacin entre dos fases.
Altura de lmina libre adimensionalizada
Segmento de lnea.
Altura de lmina libre sobre aliviadero.
S = s/HD ..
SPANCOLD Spanish Committe on Large Dams.
Tensin supercial.
Ti
Intensidad turbulenta.
TT
t
t0
U
u
uf s
Escala temporal turbulenta.

Variable temporal.
Instante temporal conocido.
Velocidad de corte.
Componente de la velocidad en el eje
Velocidad del ujo irrotacional.
xix
x.
Variable escalar.
Notacin
ur
urms
Velocidad de ascensin de las burbujas.

Fluctuacin de las velocidades.
ui
Valor medio de la componente
Volumen de un slido rgido.
de la velocidad.
VOF Volume of Fluid.

v Componente de la velocidad en el eje
~v
Vector velocidad con componentes
~v
Vector velocidad media.
v1
vmax
We
~
W
u, v , w .
Velocidad media mxima para una seccin situada a una distancia adimensional
v~i
y.
X.
Vector velocidad de la interfase entre dos uidos.

Velocidad media mxima a la entrada del cuenco de amortiguacin.
Nmero de Weber.
Vector de variables conservativas.
Componente de la velocidad en el eje
z.
Distancia longitudinal adimensionalizada con el calado total a la entrada del

cuenco.
~
X
~x
xlong
Vector posicin para un cuerpo.

Vector posicin con componentes
x, y , z .
Coordinada longitudinal medida desde el inicio del resalto hidrulico.
xs
Coordinada longitudinal medida desde la cresta del aliviadero.
y+
Distancia adimensional a la pared.
y90
Calado para el cual la concentracin
ylong
C(y90 ) = 0,90.
Coordinada perpendicular a la solera medida desde el inicio del resalto hidrulico.
xx
Notacin
Alfabeto griego
Tensin supercial.
Diferencia nita de una variable.
Espesor de la capa lmite.
Altura sobre la solera correspondiente al punto donde se encuentra la velocidad

mitad de la velocidad mxima en la seccin.
ij
Delta de Kronecker.
Espesor de desplazamiento.
Vector derivadas con componentes
, y
, z
).
( x
Tasa de disipacin de energa turbulenta por unidad de masa.
Coordenada curvilnea.
Coordenada curvilnea.
Orientacin de un cuerpo.
Espesor de cantidad de movimiento.
Curvatura de una supercie entre dos uidos.
Energa cintica turbulenta por unidad de masa.
Constante de von Karman.
Viscosidad dinmica.
Viscosidad dinmica laminar.
Viscosidad dinmica turbulenta.
Viscosidad cinemtica.
Viscosidad cinemtica laminar.
Viscosidad cinemtica turbulenta.
Coordenada curvilnea. Componente del vector aproximacin a la solucin
Densidad.
Densidad del ujo medio.
xxi
x.
Notacin
Coeciente de tensin supercial.
Constante del modelo
.
Constante del modelo
.
Tensor de tensiones.
Tensor de esfuerzos viscosos.
Tensor de esfuerzos de Reynolds.
Velocidad de propagacin de una perturbacin.
Altura de embalse adimensionalizada
Volumen.
Tasa de disipacin de energa turbulenta por unidad de energa cintica tur-
= H/HD .
bulenta o tasa especca de disipacin turbulenta
Vector vorticidad.
xxii
/.
Introduccin
A lo largo de la historia, el agua ha tenido un protagonismo fundamental en

la vida del ser humano. Ha modelado y desplazado al hombre desde los primeros
asentamientos. La necesidad de almacenar y transportar el agua ha obligado a los
seres humanos a esforzarse para alcanzar sus objetivos. Todo esto ha derivado una
tecnologa que ha permitido el avance y progreso de las distintas civilizaciones. Precisamente del ingenio evolucionado para suplir las necesidades fundamentales del ser
humano, nacen las presas, obras faranicas con el objetivo de derivar y almacenar
el agua para su posterior uso.
Dentro de la ingeniera de presas, los desages presentan una especial complejidad tcnica. Por un lado, las tomas de explotacin estn estrechamente relacionadas
con el objetivo de la obra; ste no tendra utilidad alguna sin ellas, pues forman
un todo nico funcional. Pero junto con este objetivo de explotacin, la presa se
encuentra con otro hecho de vital importancia: la necesidad de evacuar el agua sobrante de las avenidas.
Dado que estos caudales extremos no se presentan muy localizados en el tiempo,

concentrados en avenidas de duracin relativamente corta con caudales muy grandes, su evacuacin plantea problemas de gran envergadura, no solo por los caudales
en s, sino por la elevacin del nivel producida por la presa en el cauce. Esta energa
suplementaria ha de amortiguarse de alguna forma para evitar la erosin aguas abajo de la obra que puede comprometer la estabilidad de la estructura. Cuando la obra
falla, los daos humanos y medioambientales son muy elevados. Afortunadamente,
la probabilidad de ocurrencia de estos fenmenos es muy reducida.
Los rganos destinados a la evacuacin de caudales sobrantes se llaman aliviaderos y las estructuras diseadas para garantizar el reintegro del caudal aliviado en
las condiciones energticas necesarias son los cuencos de amortiguacin.
En la ingeniera de presas, la tcnica ha avanzado considerablemente aprendiendo
Captulo 1. Introduccin
muchas veces de nuestros propios errores. Como es natural, los mtodos de clculo
empleados son cada vez ms avanzados y las posibilidades que ofrece la potencia
computacional y el desarrollo de nuevos algoritmos permite hoy en da calcular aspectos hidrodinmicos que tan solo hace una dcada eran impensables. En el terreno
de lo experimental, los mtodos tambin han avanzado mucho, tanto los intrusivos
como para los no intrusivos siendo posible medir velocidades o concentraciones sin
alterar el ujo. Esto ha permitido profundizar en estos ujos tan complejos bien en
prototipo 1.2 como en modelo reducido 1.1.
Figura 1.1: Flujo supercrtico reproducido en laboratorio (Pster, 2008).
Este avance tambin ha tenido lugar en muchos otros campos de la ingeniera, sin
embargo en el caso de la hidrulica de lmina libre el salto metodolgico es mucho
mayor. En este campo, es frecuente encontrarse con relaciones empricas obtenidas
en laboratorio ajustadas mediante coecientes experimentales. Otras veces, se hace uso de modelos numricos uni- o bidimensionales, frecuentemente acoplados con
relaciones empricas que permiten calcular fenmenos no contemplados en las hiptesis de partida de las ecuaciones empleadas. En el caso de los ujos supercrticos
de los aliviaderos; las grandes velocidades del ujo, la componente tridimensional
y el emulsionamiento del lete uido, hacen casi impensable el uso de un modelo
numrico que no recoja esta gran complejidad del ujo turbulento emulsionado.
Mientras el ingeniero hidrulico clsico manejaba variables como calados y velocidades medias, el estado de la tcnica actual obliga a hablar de reajustes en la
turbulencia, interaccin entre capa lmite y lmina libre o funciones de densidad
de probabilidad de tamaos de burbuja; cuando a la hidrulica de aliviaderos nos
referimos.
1.1. Motivacin
Es evidente, como ya se ha comentado en el captulo 2.1.1 de este trabajo, que
la evolucin de las presas tanto en nmero como en altura ha ido siempre de la
Figura 1.2: Flujo supercrtico emulsionado sobre aliviadero y aireadores (presa de

Alicura, Argentina). Foto tomada por el profesor Dr. H.-E. Minor.
mano de la evolucin en la ciencia. Las bases cientcas empleadas en la actualidad

podran quedar perfectamente recogidas en los tratados ya publicados a nales de
los 70 (USBR, 1977; Peterka, 1958), con avances casi cada 20 aos que han marcado
importantes hitos cientcos (Hager, 1992; Vischer and Hager, 1995). Por otro lado,
los pilares del anlisis de los ujos agua-aire y los problemas de cavitacin seran
estudiados poco ms tarde (Falvey, 1980, 1990; Wood, 1991; Chanson, 1996), habindose convertido la aireacin del ujo en un asunto central, como as lo revela
la ferviente actividad investigadora de los ltimos aos (Pster and Hager, 2010a;
Pster et al., 2006; Ma et al., 2011b; Meireles et al., 2014; Kramer et al., 2006;
Chanson, 2009; Bung, 2013; Oertel and Bung, 2012; Bombardelli et al., 2011; Bung,
2011), pudiendo destacar los dos Vision papers publicados muy recientemente por el
Journal of Hydraulic Research sobre ujos agua-aire (Chanson, 2013) y estructuas

hidrulicas (Hager, 2014) y que ponen todava ms de maniesto el inters actual
de la comunidad cientca por las grandes obras hidrulicas.
El diseo de los rganos de desage ha ido siempre asociado a modelos fsicos

reducidos. Sin embargo, los efectos de escala asociados al modelo reducido dejan
siempre varios aspectos sin resolver (Heller, 2011; Chanson, 2009). Este es el caso
de la aireacin del ujo. Para los nmeros de Weber y de Reynolds asociados a los
modelos reducidos, los efectos de escala no permiten una correcta reproduccin de la
aireacin del prototipo (Kramer et al., 2006; Pster and Hager, 2010a; Pster et al.,
2006). La importancia de la aireacin del ujo, tanto por la proteccin que ofrece
al aliviadero (Falvey, 1990; Kramer, 2004) como por la reduccin en la friccin de
fondo (Madavan et al., 1983; Wood, 1991; Gutirrez Serret, 1995; Chanson, 1996),
convierte a la concentracin de aire en una variable de suma importancia, al mismo
nivel que la velocidad o la presin. De hecho, esta concentracin de aire modica las
propiedades macroscpicas del ujo y da lugar a distintas distribuciones de cantidad
de movimiento. Adems se han encontrado importantes diferencias en el comportamiento de la turbulencia debido al efecto de la fase dispersa (aire) sobre la fase
portante (agua) (Balachandar and Eaton, 2010; Crowe, 2000; Kussin and Sommerfeld, 2002; Zhou, 2010; Yeo et al., 2010; Crowe, 2010). Esto, sin ninguna duda, afecta
las condiciones de entrada al cuenco y puede invalidar un diseo hidrulico que por
otro lado podra ser perfectamente aceptable para un ujo no aireado. Por ello, es
interesante emplear en el diseo otro tipo de tcnicas que permitan una aproximacin al ujo producido en el prototipo, a escala 1:1.
En los ltimos quince aos los estudios hidrulicos han ido, cada vez ms, acompaados de tcnicas CFD que permiten un anlisis detallado del ujo. Estas tcnicas
han avanzado mucho desde las primeras aplicaciones hidrulicas (Caisley et al., 1999;
Ho et al., 2001; Savage and Johnson, 2001), permitiendo actualmente simulaciones
de grandes obras hidrulicas con modelos para ujos bifase. El empleo de las tcnicas CFD puede permitir el anlisis rpido y able de mltiples alternativas, con el
aliciente del nuevo enfoque Computational Multiphase Fluid Dynamics que permite
acoplar el efecto del aire sobre el agua del mismo modo que el modelo de turbulencia
acopla el efecto de las pequeas escalas sobre el ujo medio.

Dada la situacin actual en los pases desarrollados con el recurso en gran medida
regulado, no as todava en los pases en vas de desarrollo, surgen nuevas posibilidades dentro de la ingeniera de presas como son la adaptacin de presas ya existentes
a las imposiciones de las nuevas normativas tcnicas, adems de las obras de proteccin contra avenidas an por construir. La adaptacin de obras de disipacin de
energa a nuevos caudales de diseo derivados de los nuevos estudios hidrolgicos
es una tarea que tendr que acometerse a medio plazo. La remodelacin de estas
obras hoy por hoy obliga a estudios tcnicos muy avanzados y grandes desembolsos
por parte de la administracin, siendo lo normal la construccin de una obra nueva
auxiliar que amortige estos caudales excedentes. Adems, muchos de los aliviaderos
construidos pueden perfectamente no estar protegidos contra cavitacin.
1.2. Objetivos
En este contexto y con el trabajo de investigacin que aqu se presenta, se pretende dar respuesta a:
La modelacin hidrulica tridimensional de ujos multifase que tienen lugar
sobre las grandes estructuras hidrulicas. Siendo estos ujos una mezcla con
alta concentracin de aire. El rango de aplicacin del trabajo realizado puede
ser muy amplio y con un potencial futuro todava por descubrir. Se pretende
por lo tanto desarrollar una metodologa clara para el futuro estudio de ujos
altamente aireados en grandes estructuras hidrulicas.
El estudio y calibracin de un modelo de reaireacin turbulenta ya implementado en el software comercial FLOW-3D. Se han analizado los parmetros
involucrados y se ha realizado una calibracin con datos obtenidos por otros
autores a escala prototipo (y por lo tanto, libres de efectos de escala). Esta
calibracin se ha planteado en los ltimos aos como una de las principales
tareas por acometer en el campo de la ingeniera hidrulica (Chanson, 2013;
Bombardelli and Chanson, 2009).
La validacin del modelo calibrado con el objetivo de identicar posibles puntos
dbiles y dar respuesta a estos problemas.
Finalmente, en este trabajo se propone una aplicacin prctica que se acaba
convirtiendo en una segunda validacin: la optimizacin del tamao de los
elementos de pie de presa contenidos en los cuencos disipadores de energa.
Estado del Arte

2.1. Resea histrica
2.1.1. Las presas en Espaa y en el mundo
Las presas son unas de las obras ms antiguas que existen. Ya en el antiguo
Egipto se construyeron varias, entre ellas, la hecha por el primer faran: Narmer
(hacia 3050 a.C.). sta serva para desviar el Nilo en Mens y tena una altura
de 15 metros y era de sillera. Dur unos 45 siglos. Existen evidencias tambin en
Mesopotamia, China y Mxico. En el periodo romano se producen gran cantidad
de obras hidrulicas, algunas notables por su concepto y dimensiones, y muchas de
ellas perdurables. Sus presas sobrepasaron la decena de metros.
En el medievo las persas fueron ms bien de pequea altura y destinadas a mover molinos. A partir del Renacimiento las presas van cobrando ms importancia
en nmero y altura, pero siempre sobre bases empricas y de experiencia, lo que lleva a algunos fracasos al hacer ciertas innovaciones sin el apoyo de una base cientca.
En realidad, no es hasta los siglos XVII y XIX que no se dispone de los fundamentos cientcos sucientes para proyectar las presas con un rigor mnimo. Con la
tecnologa desarrollada, la ingeniera de presas experimenta un cambio sustancial y
cualitativo y a partir de la mitad del siglo XIX se comienza a dar forma y dimensiones a las presas sobre bases fsicas y mecnicas y clculos racionales, inicindose
el periodo que genricamente puede llamarse la de las presas modernas (Vallarino,
1991). As, en un siglo las presas han avanzado incomparablemente ms que en los
40 anteriores.
En la tabla 2.1 se recoge la historia estadstica de las presas conocidas en Espaa
hasta el siglo XIX. Algunas de las presas no estn recogidas en el inventario del
MOP segn el criterio selectivo empleado en 1991 .

Lo primero que llama la atencin es el gran nmero de presas conocidas y el gran
1 Actualmente
el inventario del MAGRAMA del 2006 solo reconoce 55 presas anteriores a 1900.
Captulo 2. Estado del Arte
Siglo
Total
En servicio
Altura mx. (m)
Nm. segn MAGRAMA
II
10
24
II-V
11
10
XI
11
XIII
XIV
25
XVI
46
XVII
12
24
XVIII
25
19
49
12
XIX
11
50
33
91
63
57
Tabla 2.1: Presas en Espaa hasta el s. XIX, adaptado de (Vallarino, 1991).
porcentaje de ellas que todava estn en servicio (casi el 70 %). Estas presas todava
en servicio, como es evidente, han sido objeto de reformas y refuerzos a lo largo de
los aos. De otro modo, no podran haber subsistido.
Tambin es digno de mencin el nmero de presas romanas, que podra ser mayor
dado que algunas presas atribuidas al siglo X podran tener un origen romano.
Del periodo medieval hasta antes del siglo XIII slo se conservan presas rabes
de baja altura, de derivacin para riegos. Del XIV es probablemente la presa ms
antigua en forma de arco que se conserva.
A partir del XVI comienza una interesante actividad hidrulica, con una cierta
inexin en el XVII, y un gran impulso en el XVIII, sobre todo con Carlos III, tanto
en Espaa como en Amrica. Aumenta la altura de las presas, llegando a 49 m, y
varias presentan un diseo original. Del XVI son algunos azudes de derivacin de
corta altura, pero vertiendo caudales importantes. Son tambin de este siglo algunas
presas de arco y bveda, bveda y arco. Descatable la de Elche (construida entre
1632-1640) que fue la primera presa bveda del mundo.
El gran desarrollo y el afn innovador de la poca carolina tuvo que pagar su
tributo al pretender extrapolaciones sin la imprescindible base cientca, que permite interpretar los resultados de la experiencia y ampliar lo conocido. Este es el
caso de la presa de Puentes prxima a Lorca que se colaps en 1802 y produjo ms
de 600 vctimas; tena 50 m y se haba cimentado sobre pilotes de madera en un
lecho de gravas que fall. Y la de El Gasco (Torrelodones, Madrid) comenzada en
1787 y proyectada para 93 m de altura, casi el doble que la ms alta hasta la fecha,
colaps cuando lleg a 57 m. Estaba formada por una estructura interna de cajones
cerrados por tabiques relativamente delgados rellenos de tierra que con las lluvias
se entumecieron y la presa rompi desde dentro. De todos modos, tampoco hubiera
resistido el empuje del agua.

El gran bagaje cientco del que ya se dispona a nales del XVIII llev a la
creacin de la Escuela de Caminos en 1804, a partir de este momento las obras en
general, y las presas en particular, fueron proyectadas sobre bases ms cientcas y
no solo por meras reglas del arte.
La poltica hidrulica, siempre presente en romanos y rabes, e intensicada con
la inquietud del Renacimiento y de la Ilustracin, tiene en el trnsito de los siglos
XIX al XX sus nuevos apstoles en Joaqun Costa y Rafael Gaset, y ya en los aos
veinte, en Rafael Benjumea (I conde de Guadalhorce) y Lorenzo Pardo, que crean
en 1926 las Confederaciones Hidrogrcas, con la misin de coordinar, fomentar y
desarrollar las actividades hidrulicas de cada cuenca: estudios, proyectos, contruccin, explotacin y relaciones con los usuarios. Estas entidades subsisten con ciertas
adaptaciones para acoplarlas a las competencias de las Comunidades Autnomas.
Desde 1950 se produce un impresionante crecimiento absoluto y relativo, con un
total de 581 presas construidas en treinta aos (casi 20 por ao), y un mximo de
203 presas en la dcada 60-70. En la dcada 81-90 se han terminado 157 presas (unas
16 anuales), que todava es un buen ritmo, aunque algo menor que el anterior, lo
que es lgico despus de un largo de perodo de intesidad constructora, una vez alcanzada una mejora sustancial de la regulacin. Esta causa, unida a cierta oposicin
ecologista, se acenta en la dcada de 1990-2000, en la que se construyen 142 presas.
Decenio
Presas
Porcentaje
Antes de 1900
55
3,58 %
1900 - 1909
14
0,91 %
1910 - 1919
52
3,38 %
1920 - 1929
41
2,67 %
1930 - 1939
48
3,12 %
1940 - 1949
60
3,90 %
1950 - 1959
161
10,47 %
1960 - 1969
212
13,78 %
1970 - 1979
187
12,16 %
1980 - 1989
164
10,66 %
1990 - 1999
186
12,09 %
Posteriores al 2000
45
2,93 %
Tabla 2.2: Evolucin histrica de las presas en Espaa.

Entre 1950 y 1970 se hicieron casi todos los embalses de mayor capacidad, mul3
tiplicndose por seis la capacidad total. La existente en el 2000 es de 55.208 hm ,
un 50 % mayor que en 1970. Una buena parte de las presas de aquel perodo tenan
como objetivo la produccin de energa elctrica, que se multiplic por 10 de 1945
a 1970.
La capacidad total de los embalses en el 2000 representa el 50 % de la aportacin

anual media de las cuencas peninsulares. Si esta regulacin estuviera equitativamente
repartida, supondra la regulacin de los caudales en un ao medio, pero el desigual
reparto regional y local hace que algunas de las cuencas queden escasas de regulacin,
sobretodo en aos secos, mientras otras disfrutan incluso de cierta interanualidad.
Aproximadamente el 40 % de la capacidad actual embalsable corresponde a embalses hidroelctricos, que es una de las proporciones ms altas de Europa y del mundo.
En cuanto a actividades futuras, todava quedan presas por hacer basndonos

en dos aspectos fundamentales: la proteccin contra avenidas y la regulacin general; es decir, aumentar los mnimos y disminuir los mximos. El primer objetivo,
considerado hasta ahora complementario en algunos embalses, ha adquirido una importancia creciente. Por una parte asignando un mayor papel a los embalses actuales
en la laminacin de avenidas, incluso a costa de otros benecios directos tangibles,
o recrecindolos para este y otros objetivos. Y tambin mediante nuevas presas dedicadas exclusiva o preferentemente a la laminacin. Para ese uso podrn ser tiles
muchas cerradas de permeabilidad dudosa, no aptas para una retencin prolongada,
pero s para este tipo de regulacin espordica y breve en la que tanto da que el
caudal salga por los desages o ltrado por el terreno, pero al n regulado. Para el
segundo objetivo, podrn construirse embalses de capacidad media repartidos por
el territorio y quizs recrecer algunos existentes al objeto de mejorar los caudales
mnimos y su garanta interanual, y al propio tiempo corregir las desigualdades actuales; presas de menor entidad pueden ser tiles para regular cuencas parciales que
queden fuera de la inuencia de los grandes embalses, bien para abastecimiento de
agua o para laminacin de avenidas.
Otra actividad importante y creciente con el tiempo ser la conservacin y eventual reparacin o refuerzo del gran nmero de presas existente, que cada vez sern
ms antiguas. El envejecimiento de las presas plantea varios problemas, adems de
su deterioro, como por ejemplo que cada presa ha sido proyectada y construida segn los conocimientos de la poca. La experiencia ha enseado que algunos de estos
eran imperfectos, lo que se ha tenido en cuenta en las construidas despus. Algunas
presas consideradas sucientemente seguras hace cincuenta aos pueden plantear
dudas segn los conocimientos actuales. Adems, el cauce aguas abajo puede haber
cambiado sustancialmente, con nuevas viviendas, fbricas, cultivos, etc. Esto puede ser tan importante que obligue a una reconsideracin de la seguridad existente.
Adems, como al pasar el tiempo van envejeciendo simultneamente cada vez mayor
nmero de presas, el problema conjunto de la seguridad presenta un aspecto muy diferente del que tena hace unos aos. Es obvio que no hay una edad nica que dena
el comienzo del envejecimiento de una presa, pues este vara con el tipo, la calidad
de ejecucin y sobretodo con el mantenimiento. Sin embargo, se puede considerar
que a partir de 50 aos, como media, una presa requerir atenciones especiales (Va-
10
llarino, 1991). En el ao 2000 casi una de cada cuatro presas tena al menos 50 aos
de edad, porcentaje que aumenta considerablemente para 2020 cuando ser del 50 %.
En Espaa existen actualmente construidas (o en fase de construccin) 1.230
presas segn el ltimo inventario de presas realizado por el MAGRAMA (2006).
La propia web de la SPANCOLD mantiene un apartado especco dedicado a esa
actividad titulado presas en Espaa. Otro titulado Actividad en la construccin,
en el que se incluye la ubicacin y caractersticas ms bsicas de las presas que se
encuentran en ejecucin (altura, tipologa, etc.) y que SPANCOLD mantiene actualizado. En este sentido, cabe sealar que actualmente se encuentran en fase de
construccin un total de 18 presas, 2 de las cuales son recrecimientos de otras tantas
presas o embalses existentes.
Altura
Presas
Porcentaje
Menores de 15 m
550
35,76 %
Entre 15 y 30 m
465
30,23 %
Entre 30 y 60 m
332
21,59 %
Entre 60 y 100 m
148
9,62 %
Mayores de 100 m
43
2,8 %
Tabla 2.3: Estadstica presas en Espaa segn su altura.
Como es evidente, una presa es una obra que requiere una alta garanta de seguridad, dadas las graves consecuencias que se derivaran de un eventual fallo. La
importancia econmica y social del problema ha llevado a varios pases con gran
nmero de presas a promulgar unas normas de seguridad que se han de cumplir en
las difersas fases de la vida de una presa.
En Espaa se public en 1967 la Instruccin para el Proyecto, Construccin y
Explotacin de Grandes Presas, que ha estado vigente hasta el 30 de marzo de
1996, fecha en la que se promulg por la Direccin General de Obras Hidrulicas el
Reglamento Tcnico de Seguridad de Presas y Emblases. Actualmente, las nuevas
Normativas Tcnicas de Seguridad de Presas se encuentran en fase de discusin.
El Reglamento es de aplicacin a todas las grandes presas que se han proyectado o
construido a partir de su publicacin, fueran de propiedad pblica o privada. Para las
presas de propiedad privada existentes previamente al Reglamento, sigue vigente la
Instruccin de 1967, porque no puede exigrseles sin compensacin unas condiciones
diferentes. En cambio las presas del Estado han de adaptarse a las nuevas exigencias,
y en un futuro cercano, a las de las nuevas Normativas Tcnicas.
2.1.2. Dinmica de Fluidos Computacional (CFD)

La Dinmica de Fluidos Computacional, del ingls Computational Fluid Dyna-
mics (CFD), hace referencia a la rama de la Mecnica de Fluidos que consiste en el

11
empleo de ordenadores y tcnicas numricas para resolver todos aquellos problemas

fsicos que estn relacionados con el movimiento de los uidos y, en ocasiones, de
otros fenmenos asociados como la transferencia de calor, las reacciones qumicas,
el arrastre de slidos, etc.
En general, el CFD comprende un amplio abanico de disciplinas cientcas, entre las que cabe destacar las matemticas, la programacin, las ciencias fsicas y la
ingeniera, que deben aunarse para dar lugar al desarrollo del cdigo que sea capaz
de resolver las ecuaciones del ujo de forma satisfactoria.
Lgicamente, la consolidacin de estas tcnicas ha sido una consecuencia ms del

progresivo desarrollo de los ordenadores desde los aos 1950-1960. Hasta entonces,
los mtodos para resolver numricamente las ecuaciones del ujo no eran viables, al
no disponerse de mquinas capaces de ejecutar un gran nmero de operaciones de
clculo por unidad de tiempo.
Hace ya casi dos siglos desde que las ecuaciones de gobierno de la mecnica de
uidos quedaron denitivamente formuladas por Claude Navier (1785-1836) y George Stokes (1819-1903) cuando introdujeron los trminos de transporte viscoso a las
ecuaciones de Euler (1707-1783), dando lugar a las famosas ecuaciones de NavierStokes.
Estas ecuaciones incluyen las leyes de conservacin para la masa, la cantidad

de movimiento y la energa de un ujo. Desgraciadamente, se constituyen en un
sistema acoplado de ecuaciones, del que es prcticamente imposible obtener una solucin analtica. Por esta razn, la experimentacin y el anlisis dimensional siempre
acompaaron histricamente la va analtica, como dos herramientas esenciales en el
estudio de la Mecnica de Fluidos, para validar y contrastar los limitados estudios
tericos.
Afortunadamente, gracias a la aparicin de los modernos ordenadores, sobre todo a partir de la dcada de 1970, se dispuso de una nueva herramienta: el CFD, que
complementa tanto a los estudios tericos como a los experimentales. Basadas en
estas mismas ecuaciones, las modernas tcnicas CFD permiten su resolucin numrica en un tiempo razonable, permitiendo simular tanto ujos reales, de inters en
la industria, como ujos cannicos, para la investigacin o la docencia.
A principios del siglo XX, ya sentadas las bases tericas de la Mecnica de Fluidos, las investigaciones en este campo de la ciencia se centraron en el estudio de
capas lmite, de capas de cortadura y, en general, de la turbulencia en los ujos. De
hecho, el anlisis de la turbulencia sigue siendo uno de los mayores retos a los que
se enfrentan los investigadores en este campo cientco; su modelizacin y su tratamiento son ahora mismo la principal lnea de investigacin dentro del CFD. Osborne
Reynolds, en 1894, sienta las bases del Reynolds Averaged Navier-Stokes equations
12
(RANS), el promediado de Reynolds y la deduccin de la ecuacin turbulenta.
Los primeros trabajos en esta materia son obra de Ludwig Prandtl (1875 - 1953),
quien desarroll una teora sobre la capa lmite que ya introduca el concepto de
longitud de mezcla para el modelado de los trminos turbulentos en las ecuaciones.
Posteriormente, Theodore Von Krman (1881 - 1963) estudi el desprendimiento de
vrtices en las estelas de cuerpos en ujo externo, fenmeno que se conoce precisamente como calles de vrtices de Von Krman. Georey Ingram Taylor (1886 1975) avanz en la utilizacin de modelos estadsticos para el estudio de la turbulencia, deniendo microescalas para la isotropa de la turbulencia. Sus ideas seran
la base para la introduccin de la teora de las escalas de la turbulencia, formulada
por Andri Nikolyevich Kolmogrov (1903 - 1987) en 1941, y en la que se dena
el concepto del espectro universal de energa de la turbulencia y de la cascada de
energa.
Aunque no es fcil determinar la fecha exacta de los primeros clculos realizados

utilizando tcnicas CFD, si se puede citar al ingls Lewis Fry Richardson (1881 1953) como el precursor de la utilizacin de dichas tcnicas. De hecho, fue el primero
en utilizar un modelo numrico para la prediccin meteorolgica, dividiendo el espacio en celdas y empleando un primitivo mtodo por diferencias nitas. Su propio
intento de predecir la evolucin de las condiciones meteorolgicas durante un periodo
de 8 horas necesit de 6 semanas de clculos intensivos que terminaron en fracaso.
En posteriores estudios comenz a emplear ordenadores humanos, habitaciones repletas de personas provistas con rudimentarias calculadoras, que ejecutaban hasta
2000 procesos iterativos por semana. Aproximadamente, cada persona habra ejecutado 20000 operaciones por semana, lo que en una jornada de 40 horas semanales,
equivaldra a unos 0,140 ops. Los rudimentarios mtodos de Richardson, realizados
en las dcadas de 1910 y 1920, suponen el punto de partida del empleo del CFD,
e incluso profetizan el empleo del clculo computacional en paralelo como va para
disponer de mayores capacidades de clculo.
En la dcada de 1930, las irresolubles limitaciones de los estudios analticos fueron

una motivacin muy importante para continuar con el lento desarrollo de metodologas computacionales. As, la primera simulacin numrica (an sin computadores)
del ujo alrededor de un cilindro fue realizada en 1933 en Inglaterra por A. Thom
(1933), y comunicada por el propio G.I. Taylor. Similares resultados fueron obtenidos por M. Kawaguti en Japn, resolviendo las ecuaciones de Navier-Stokes para el
ujo alrededor de un cilindro para un nmero de Reynolds de 40. Segn palabras del
propio Kawaguti (1953), "la integracin numrica de este estudio precis de un ao
y medio, con 20 horas de clculo por semana y una considerable cantidad de trabajo
y resistencia". No en vano, l mismo realiz todos los clculos, usando simplemente
una calculadora de escritorio.
Otra gura imprescindible por su aportacin a los modelos de turbulencia es el
13
ya citado Andrei Nikolaevich Kolmogorov, cuyo artculo en 1941 sobre la teora del
equilibrio universal de la turbulencia sent una de las bases de la comprensin moderna de este fenmeno. Sin embargo, esta publicacin careci del reconocimiento
merecido debido a la explosin de la segunda guerra mundial, lo que retras varios
aos el conocimiento de su trabajo a escala mundial.
A partir de nales de la dcada de 1950 y en toda la dcada de los 60, el laboratorio nacional de Los lamos (LANL), auspiciado por la NASA, se constituy
en el verdadero impulsor de las tcnicas CFD, desarrollando los primeros cdigos y
dando los primeros pasos en el empleo de ordenadores. Hoy da, muchos mtodos
que se emplean en programas comerciales provienen de aquellos pioneros trabajos.
Dicho laboratorio contaba con los ordenadores ms potentes de la poca, como el
MANIAC (1952) que realiz los clculos para el desarrollo de la primera bomba de
hidrgeno, o el IBM 704, que fue el primer ordenador fabricado en serie y que ya
dispona de arquitectura para calcular en coma otante. Su memoria era algo limitada, pues dispona solo de 4 Kb, pero ya permita 40000 operaciones por segundo
(aprox. 40 Kiloops).
Con el centro de computacin ms grande del mundo a su servicio, el LANL cre

un gran nmero de mtodos, como son el Particle-in-Cell (PIC), el Fluid-in-Cell
(FIC), los vorticity-stream function methods o el Marker-and-Cell (MAC). Retomando los antiguos trabajos para estudiar el desprendimiento de vrtices alrededor
de cilindros, Fromm y Harlow realizan por primera vez el clculo por computador
del desprendimiento no estacionario de vrtices en 1963 para nmeros de Reynolds
por debajo de 1000 (Fromm y Harlow, 1966). Basndose en la formulacin de vorticidad y funcin de corriente, desarrollaron un modelo explcito de diferencias nitas
para ujo incompresible. Tambin se inician los primeros intentos de simular ujos
compresibles, empleando tcnicas PIC para ujos viscosos, aunque no se consiguen
reproducir los efectos de choque snico con gran precisin. Eran tiempos de aprendizaje en los que prcticamente todo estaba por hacer y en los que los resultados
obtenidos se mostraban en grcos que se hacan a mano.
Aunque los primeros mtodos para simular ujos incompresibles empleaban formulaciones basadas en la vorticidad y la funcin de corriente, a nales de los aos
60 se consiguen desarrollar las primeras simulaciones en trminos de las variables
primitivas: velocidad y presin. Esto se consigui gracias al desarrollo del mtodo
MAC por (Harlow and Welch, 1970; Amsden and Harlow, 1970), precursor del actual mallado decalado (staggered mesh ) que se implementa hoy en da en la mayora
de cdigos por volmenes nitos. Otro avance importante fue la inclusin de los
primeros modelos de turbulencia en las simulaciones a partir del desarrollo de las
bases del modelo de turbulencia
en 1967, que fue denominado entonces como
q-d por las limitaciones de los editores de texto de la poca. A nales de los aos 60
y principios de los 70, el CFD experimenta un auge muy notable, extendindose por
otros lugares del planeta abandonando el lugar que le vio prcticamente nacer como
14
lo conocemos hoy. As lo armaba Harlow en 1968, director del departamento de

CFD en el LANL durante esos aos, reconociendo que una nueva era comenzaba".
l mismo deca que ese haba sido el ltimo ao en que an estaba al corriente de
todos los desarrollos de CFD que se hacan en todo el mundo.
A partir de la dcada de 1970, un grupo liderado por Brian Spalding en el Impe-
rial College de Londres toma el relevo en la frontera del conocimiento de las tcnicas
numricas. Inspirndose en los trabajos del LANL, (Patankar and Spalding, 1972)
desarrollan una formulacin implcita (basada en el upwind dierencing ) en trminos
de velocidad y presin, en la que introducen por primera vez el mtodo de acoplamiento SIMPLE. Estos nuevos mtodos implcitos tienen una clara ventaja frente
a los explcitos, pues no presentan restricciones en la discretizacin temporal, garantizando siempre la estabilidad numrica. De este modo, el algoritmo SIMPLE se
impuso fcilmente y domin por completo durante una dcada (Van Doormaal and
Raithby, 1984). Posteriormente se desarrollaron toda una serie de modicaciones y
mejoras al modelo de acoplamiento, dando lugar a los mtodos SIMPLER (1980),
SIMPLEC (1984) o PISO (1986); una comparacin de los mtodos puede encontrarse en (Jang et al., 1986). Adems, en 1972, el propio Spalding en colaboracin con
Launder (Launder and Spalding, 1972) universaliza el modelo de turbulencia
tal y como se emplea hoy da, sin duda el modelo ms robusto y de mayor utilizacin
de entre todos los existentes hasta la fecha de hoy. Pronto queda patente la apuesta
por la economa en trminos de coste computacional, el rango de aplicabilidad y el
realismo de los modelos de turbulencia de dos ecuaciones (Launder and Spalding,
1974). En 1980, Suhas V. Patankar publica Numerical Heat Transfer and Fluid Flow
(Patankar, 1980), probablemente el primer gran libro que trata en profundidad las
metodologas de CFD y que ha servido de inspiracin para la creacin de innidad
de cdigos numricos. Patankar y el Imperial College sientan las bases denitivamente del Mtodo de Volmenes Finitos con esta notable contribucin.
El auge de estas tcnicas, junto con todos los avances conseguidos, llev a la industria aeronutica y aeroespacial a lanzarse a una carrera desbocada por el empleo
del CFD en sus fases de diseo y vericacin de modelos. A principios de los 70,
empresas punteras como el consorcio anglo-francs para la fabricacin del Concorde,
o la empresa americana Boeing (Seattle, USA) comienzan a emplear cdigos 3-D,
incompresibles de ujo potencial, obteniendo unos esperanzadores resultados con
mquinas CDC6600, que ya superaban el Megaop de capacidad. Estos xitos les
permitieron una notable reduccin del nmero de ensayos experimentales en tnel de
viento y una drstica reduccin de costes. Por aquel entonces, apenas se realizaban
100 a 200 simulaciones al ao, nada que ver con las actuales cifras. La progresiva
implantacin de estas tcnicas numricas ha reducido las pruebas experimentales
desde los 77 ensayos efectuados a nales de los 70 para el Boeing 767 a las escasas 5
utilizadas para el nuevo modelo 7E7 (base del Boeing 787, que entr en servicio en
2009). Hasta mediada la dcada de los 80, las empresas punteras an desarrollaban
sus propios cdigos, muchas veces apoyadas por centros tecnolgicos que les daban
15
soporte y nuevas soluciones.
Sin embargo, todo cambia con la llegada de la dcada de los 80. A la sombra de
la progresiva publicacin de nuevos textos especializados, comienzan a aparecer en el
mercado los primeros paquetes de CFD comerciales. Sirva de apunte que hasta 1975,
por ejemplo, ninguno de los cdigos desarrollados en el LANL haba sido distribuido
fuera de su mbito acadmico. Sin embargo, muchas empresas que se haban interesado por el uso de CFD en sus investigaciones, dejaron de desarrollar sus propios
cdigos para comenzar a utilizar los nuevos paquetes comerciales. El propio Spalding
funda Concentration Heat and Momentum LTD (CHAM) en Wimbledon, Inglaterra, en 1974 (empresa que an existe en la actualidad), comenzando a distribuir
cdigos especcos segn las necesidades de cada empresa. Pronto se dieron cuenta
de que mantener soporte para un gran nmero de cdigos distintos no era eciente,
as que en 1980 modican su losofa y presentan el primer cdigo de propsito
general, denominado Parabolic Hyperbolic Or Elliptic Numerical Integration Code
Series (PHOENICS). Fue lanzado al mercado en 1981, considerado un hito por ser
el primer paquete comercial de CFD, pensado para modelizar cualquier tipo de ujo
(dentro de las limitaciones de los modelos existentes hasta la fecha). Originalmente, este cdigo permita modelizar dominios discretizados con celdas estructuradas
(hexaedros), analizar ujos estacionarios y no estacionarios, compresibles e incompresibles, laminares y turbulentos e incluso podan modelarse reacciones qumicas y
ujos de una o dos fases.
Al mismo tiempo que Spalding y su grupo del Imperial Collage lanzaban PHOENICS al mercado, un grupo de investigadores de la Universidad de Sheeld, en Inglaterra, que estaban trabajando con cdigo CFD para simular problemas de combustin, desarrollaron, impulsados por Creare - una consultora tecnolgica ubicada
en New Hampshire, USA -, un nuevo paquete comercial de propsito general llamado
FLUENT. Su lanzamiento al mercado, en 1983, contaba tambin con discretizaciones
estructuradas (hexagonales), nicamente para ujos estacionarios, con posibilidad
de modelizar ujos laminares y turbulentos, pero sobretodo potenciaba la transferencia de calor, fenmenos de combustin y fases dispersas, as como la conveccin
natural. Uno de los mayores inconvenientes segua siendo el uso de mallas estructuradas, que limitaban notablemente la discretizacin de geometras complejas. En
1987, otro empleado de Creare demostr las grandes posibilidades de utilizar mallas
no estructuradas, mejor adaptadas a las geometras, en varios proyectos impulsados desde la NASA. Sus avances se incorporaron despus a las funcionalidades de
FLUENT, dando lugar a la creacin en 1991 del cdigo solver Rampant y del programa para crear mallas tetradricas Tgrid. Naca el primer cdigo de propsito general
que resolva mallas no estructuradas, pensado para aplicaciones aeroespaciales de
ujo compresible (con altos nmeros de Mach).
A nales de los aos 80 y principios de los 90 se produce un boom.en la creacin
de cdigos comerciales, emulando a los pioneros PHOENICS y FLUENT, entre ellos
16
cabe destacar a FIDAP (Fluid Dynamics Analysis Package ), desarrollado a partir

de los trabajos del Illinois Institute of Technology al mismo tiempo que FLUENT;
y, posteriormente adquirido por FLUENT, STAR-CD desarrollado en 1987 como
spin-o de los trabajos del Imperial Collage y centrado en mallas mviles, no estructuradas; FLOW3D (sin guin) desarrollado por la Autoridad de la Energa Atmica
(AEA) en Inglaterra a nales de los aos 80 para aplicaciones nucleares; FLOW-3D
surgido de los esfuerzos del LANL en 1985 para comercializar su potente modelo
de ujos con lmina libre (SOLA-VOF) y su modelo de interpretacin de geometras en mallas estructuradas (FAVOR); TASCow, desarrollado en Canad tambin
en 1985 y posteriormente fusionado con el FLOW3D de la AEA para dar lugar a
CFD-4 a mediados de los 90 o los cdigos FINE (Flow Integrated Environments ),
desarrollados desde 1992 por NUMECA, empresa impulsada por Charles Hirsch de
la Universidad Libre de Bruselas y autor de otro de los libros de referencia sobre
CFD (Hirsch, 2007).
Las empresas aeronuticas y automovilsticas no permanecen ajenas a estos cambios y comienzan a utilizar cdigos comerciales en sus fases de diseo. La aparicin
de mallas no estructuradas permite simulaciones con un grado de detalle que hasta
entonces no era viable. As, Boeing comienza a utilizar cdigos comerciales desde
1996, primero con CFD++ y despus con FLUENT en 1998. No tardan en seguir sus
pasos otras empresas como Airbus, utilizando Flowmaster, as como otras compaas
del sector dedicadas a la fabricacin de motores de avin, como General Electric,
Pratt & Whittney o Rolls-Royce. A partir de 1995, las tcnicas CFD comienzan a
ser aplicadas en el sector del automvil. Multitud de casos, tanto de ujos externos
como de ujos internos, son analizados desde esta nueva ptica numrica. General
Motors y Ford son pioneros en el uso de estas tcnicas. En general, la aerodinmica
del vehculo en trminos de fuerza de arrastre y penetracin al aire es el tema estrella
de estas aplicaciones. Tambin hay otros muchos estudios, a nivel de componentes,
con un gran inters como pueden ser el enfriamiento de los discos de freno, el ujo
interno y la combustin en los cilindros del motor, la inyeccin de metal fundido en
los moldes, la oscilacin del ujo y generacin de vibraciones o incluso el ujo en
las conducciones de salida de gases o en el tubo de escape. Actualmente, dentro del
sector del automvil, es el mundo de la Frmula 1 el que dispone de los medios tcnicos ms avanzados, con potentes superordenadores para analizar las caractersticas
aerodinmicas de sus monoplazas. En concreto, en 2007, el equipo BMW-Sauber F1
pas a disponer del superordenador ms potente de la F1 hasta el momento, Albert
2, con 5160 procesadores de doble ncleo y una capacidad de 12,2 Teraops, en el
que ejecutaban el programa ANSYS FLUENT, con modelos que podan llegar a
alcanzar 1000 millones de celdas.
Durante los ltimos 15 aos, el auge de los modelos numricos tridimensionales

ha permitido la simulacin de grandes obras hidrulicas. Estos modelos, comnmente conocidos por las siglas CFD (Computational Fluid Dynamics ), modelan el
17
movimiento del agua
sobre las estructuras consiguiendo por lo general errores por
debajo del 10 % tanto en calados como en velocidades, no siendo raros errores por
debajo del 5 %. Los ujos con gradientes de presiones adversos muy marcados y despegues importantes del ujo, pueden dar lugar a errores mayores. An as, su uso
permite detectar zonas con funcionamiento hidrulico deciente en fase de diseo,
con lo que se han constituido como una herramienta til en fase de diseo.
Estas tcnicas con gran tradicin aeronutica, tenan una limitada aplicacin hidrulica hasta el desarrollo de algoritmos especcos para lmina libre como son los
desarrollados en (Hirt and Nichols, 1981; Nichols et al., 1980; Amsden and Harlow,
1970; McKee et al., 2008; Gingold and Monaghan, 1977; Monaghan, 1985).
Estudios pioneros con herramientas CFD en ingeniera hidrulica de lmina libre

pueden encontrarse en (Caisley et al., 1999; Ho et al., 2001; Savage and Johnson,
2001; Oertel and Bung, 2012; Wu et al., 2000; Chanel and Doering, 2007; Bombardelli et al., 2011; Meireles et al., 2014). En la Fig. 2.1 se muestra una de las primeras
simulacioens 3D de un aliviadero en Australia. En estas simulaciones, el nmero de
celdas empleadas se encontrada en torno a 10.000. Actualmente un ordenador de
sobremesa puede simular con varios millones de celdas activas.
Sin embargo, la mayora de estudios realizados se limitan al estudio de las regiones no aireadas, donde la turbulencia producida no es capaz de atrapar y arrastrar
grandes volmenes de aire. Estos modelos no han tenido conseguido grandes xitos
en su aplicacin a las regiones marcadamente bifsicas como son la zona de ujo
completamente desarrollado de un aliviadero, donde se pueden alcanzar concentraciones de aire del 65 % (Chanson, 1996; Wood, 1991). Es por ello que la calibracin
y validacin de modelos que aproximen los fenmenos de subescala de estos procesos
es una de las claves actuales para la correcta modelacin hidrulica de las grandes
infraestructuras hidrulicas (Chanson, 2013, 2009).
Debido a la complejidad de las ecuaciones y al grado de detalle que normalmente

se desea obtener de las simulaciones, los ingenieros encargados en el desarrollo de
estas tcnicas CFD siempre han necesitado las mquinas ms potentes del mundo. Por tanto, aunque se utilice un software comercial a nivel de usuario, siempre
hay que recordar que estos cdigos estn basados en un conjunto de ecuaciones no
lineales, muy complejas y acopladas entre s, que se resuelven de forma iterativa
mediante algoritmos muy especcos incluidos en el mismo paquete. El objetivo que
se persigue es que el usuario sea capaz de resolver cualquier ujo dentro de una geometra prejada, que limitamos con unas condiciones iniciales y de contorno, y para
el que los fenmenos fsicos implicados estn identicados a priori. Los resultados
del cdigo CFD comercial pueden normalmente representarse grcamente o como
un mapa de distribuciones, tanto de variables escalares como de variables vectoriales.
2 Cuando
el ujo es monofsico: solo agua

18
Figura 2.1: Simulacin 3D de la regin no aireada sobre un perl WES (Ho et al.,
2001).
La posterior evolucin de todos estos cdigos ha estado marcada por una intensa
competitividad y una febril actividad de fusiones y cambios de manos. Tras el dominio inicial de PHOENICS, FLUENT se convirti en el referente de esta industria,
sobre todo a partir de la segunda mitad de la dcada de los 90. Sin embargo, desde
2003, se ha producido una importante concentracin de estos paquetes, iniciada por
ANSYS con la compra de CFX-4, al que pasara a llamar simplemente CFX. Finalmente, tras unos aos de intensa rivalidad comercial, en 2006 ANSYS compraba a
FLUENT por un total de 630 millones de dlares, creando as al nuevo lder mundial en el desarrollo de tcnicas numricas. Actualmente, CD-ADAPCO, heredera
del cdigo STAR-CD y participada por la importante empresa de CAD Dassault
Systems (fabricante de CATIA), se puede comparar en cierto modo al nuevo gigante
de ANSYS. Por otro lado, FLOW-3D ha conseguido alcanzar un pblico con mayor
vocacin hidrulica, sobre todo en las aplicaciones con lmina libre.
En los ltimos aos, la computacin empleando tarjetas grcas como unidades

de procesado (GPU), ha permitido que algoritmos ya conocidos desde las dcadas
19
de los 70-80 como son el mtodo SPH (Smoothed Particle Hydrodynamics) que son
ms costosos (Monaghan and Gingold, 1983; Monaghan, 1992), pero con una marcada estructura explcita, hayan recobrado una importancia casi inesperada, siendo
su punto fuerte los ujos con lmina libre.
La concentracin de los cdigos ha dado lugar a un progresivo incremento en el

precio de las licencias. Sin embargo, esta tendencia puede verse contrarrestada con
la aparicin de nuevos cdigos numricos de distribucin libre, siguiendo el espritu de sistemas operativos libres, como Linux. El pionero en la creacin de cdigos
libres de CFD es OpenCFD, que desarrolla actualmente un avanzado solver libre
denominado OpenFOAM. El principal inconveniente que presentan estos cdigos es
la ausencia de una interfaz grca y la dicultad del entorno para un usuario no
familiarizado. Sin embargo, y desde la compra de OpenFOAM por parte de ESI,
han aparecido varias inciativas para hacer este software ms amigable. As, se han
desarrollado varias interfaces grcas como HELYX-OS de engys o simFlow.
Probablemente, en un futuro no muy lejano, el papel del software libre tenga un

claro impacto en un campo tan especializado como el del CFD, obligando a reducir
el precio de las licencias de los actuales distribuidores y volviendo en cierto modo a
aquellos romnticos tiempos a mediados de los 60 en los que todo era un territorio
nuevo por explorar.
En resumen, las tcnicas CFD son ya, a todos los efectos, una herramienta ms
dentro de la ingeniera asistida por computador (CAE) utilizada universalmente en
la industria. Sus posibilidades para simular todo tipo de fenmenos y ujos permiten a los diseadores y a los analistas disponer de un tnel de viento virtual en sus
centros de computacin. As mismo, el software para el CFD ha evolucionado espectacularmente, mucho ms de lo que Claude Navier o George Stokes pudieran haber
imaginado. El CFD se ha convertido en una parte indispensable en el proceso de
diseo aerodinmico e hidrodinmico en la ingeniera naval, mecnica o aeroespacial.
2.2. Fundamentos fsicos

2.2.1. Ecuaciones de gobierno
El movimiento de un uido queda completamente descrito por las leyes de conservacin para las tres propiedades bsicas de un ujo: masa, momento y energa.
Estas tres condiciones determinan por completo el comportamiento del sistema sin
necesidad de aadir nuevas ecuaciones. La nica informacin adicional es la concerniente a las propiedades del uido: incompresible, gas perfecto, uido condensable,
material viscoelstico, etc.
20
De este modo, si adems el uido de estudio es incompresible, la ecuacin de

la energa no es necesaria pues, imponiendo la conservacin de masa y momento se
tienen cuatro ecuaciones que, junto con las tres componentes de la velocidad y la
presin, da lugar a un sistema compatible determinado. nicamente es necesaria la
ecuacin de la energa cuando se precisa una ecuacin para la densidad por ser esta
variable.
Asimismo, estas ecuaciones pueden presentarse en forma diferencial o en forma

integral. Por su aplicacin al mtodo de los volmenes nitos, se presentan primero
en forma integral, aunque tambin se ha decidido expresarlas en forma diferencial
dada la interpretacin fsica de los gradientes y del resto de componentes. De este
modo, la expresin diferencial permite tambin clasicar la naturaleza matemtica
de las ecuaciones de Navier-Stokes.
2.2.1.1. Conservacin de la masa

La ecuacin de conservacin de la masa es conocida tambin como ecuacin de
continuidad o de divergencia para ujos estacionarios. Esta ecuacin es independiente de la naturaleza del uido y supone una restriccin emprica por la cual la masa
ni se crea ni se destruye.
Por lo tanto, no hay trmino difusivo para el transporte de masa, lo que quiere
decir que la masa solo puede ser transportada por conveccin. As, en ausencia de
trminos fuente o sumidero, se obtiene:
t
Donde
~
S
es la densidad y
~v
~=0
~v dS
d +
(2.1)
la velocidad,
es el volumen del volumen de control y
la supercie con direccin y sentido normal y saliente a las caras del volumen de
control.
Para un uido incompresible, la densidad es constante y la ecuacin anterior

queda de la siguiente forma:
~ ~v = 0
21
(2.2)
2.2.1.2. Conservacin de la cantidad de movimiento

La cantidad de movimiento es una cantidad vectorial denida como el producto
de la masa por la velocidad, que expresado por unidad de volumen, es el producto
de la densidad por la velocidad.
Integrando a lo largo de la supercie en un volumen de control se puede obtener:
I
~v d +
I
I
Z

~
~
~
~v ~v dS = p dS + dS + f~e d
(2.3)
donde:
2 ~
v )ij
ij = (vj /xi + vi /xj ) (~
3

(2.4)
la viscosidad dinmica, ij la delta de Kronecker

~e el vector de fuerzas de cuerpo.
que vale 1 cuando i = j y 0 para el resto de casos, y f
es el tensor de esfuerzos viscosos,
A menudo se emplea el tensor de tensiones
= pI
sin descomponer. Donde
es la matriz identidad.
La conservacin del momento supone una ecuacin en cada eje
x, y , z .
2.2.1.3. Sistema de ecuaciones de Navier-Stokes

En los apartados anteriores se ha deducido por separado las leyes de conservacin
de masa y momento. Llegados a este punto se pueden expresar estas leyes conjuntamente en un sistema de ecuaciones con el objetivo de ver con ms claridad los
trminos involucrados. De este modo, introducimos dos vectores de ujo
El primero,
F~c ,
F~c
F~v .
contiene los trminos convectivos y se suele llamar vector de ujos
convectivos aunque tambin contiene a los trminos de presin. El segundo vector,
F~v ,
es conocido como vector de ujos difusivos y contiene los esfuerzos viscosos.
Adems, se puede denir el trmino fuente
~ , que est compuesto por las fuerzas de

Q
volumen. Con todo esto en mente y haciendo uso del vector normal a la supercie
del volumen de control
~n = (nx , ny , nz ),
Z
~ d +
W
podemos escribir:
Z

~=
~
F~c F~v dS
Qd
El vector de las variables conservativas
(2.5)
~
W
para un uido incompresible est denido
22
como:

~ = u
W
v
w
(2.6)
Para el vector de ujos convectivos se tiene:
(nx u + ny v + nz w)
u(nx u + ny v + nz w) + nx p
F~c =
v(nx u + ny v + nz w) + ny p
w(nx u + ny v + nz w) + nz p
(2.7)
Para el vector de ujos difusivos se tiene:
0
nx xx + ny xy + nz xz
F~v =
nx yx + ny yy + nz yz
nx zx + ny zy + nz zz
(2.8)
Finalmente, el trmino fuente:
~ = fe,x
Q
fe,y
fe,z
(2.9)
Este sistema de ecuaciones es conocido como sistema de ecuaciones de Navier-Stokes

para un uido incompresible. Las ecuaciones describen los ujos de masa y momento
a lo largo del contorno
de un volumen de control
Esta forma de las ecuacio-
nes se conoce como forma conservativa. Su discretizacin directa, convirtiendo las

integrales en sumatorios, asegura la conservacin de la masa entre los diferentes volmenes de control.
Se ha deducido en su forma integral por su empleo en el Mtodo de los Volmenes

Finitos (MVF), pero dada la interpretacin fsica y su transcendencia matemtica,
es importante tener en cuenta tambin la forma diferencial, y no conservativa, que
servir de base para deducciones futuras y que se obtiene en aplicacin del teorema
de Gauss al anterior sistema, teniendo en cuenta adems la incompresibilidad del
uido y denotando como
gi
a la proyeccin de la gravedad sobre el eje
23
i:
vi
= 0
xi
vi
vi
1 p
+ vj
= gi
+ v2 vi
t
xj
xi
(2.10)
La forma no conservativa del sistema de ecuaciones es adecuada para ujos suaves con bajas velocidades e incompresibles en los que las supercies de la solucin
son relativamente suaves y no discontinuas. La forma conservativa del sistema de
ecuaciones es conveniente para ujos con discontinuidades como ondas de choque en
ujos compresibles y, por lo tanto, ujos con altas velocidades.
Si en las anteriores ecuaciones se anularan los trminos de segundo orden, es decir, los viscosos, se tendran las ecuaciones de Euler o ecuaciones para un uido ideal
que seran vlidas para ujos altamente convectivos o para uidos reales lejos de los
contornos. Si se realizaran ms simplicaciones, se podran obtener otras formas como las ecuaciones de Saint Venant, de ujo potencial o de funcin de corriente, muy
empleadas todas ellas tradicionalmente para ujos marcadamente bidimensionales o
suaves con un comportamiento marcadamente laminar. Tambin podran obtenerse
las ecuaciones para capa lmite anulando los gradientes de presiones en la direccin
vertical y haciendo una serie de simplicaciones de las derivadas de la velocidad en
la direccin normal a la supercie.
Cada una de estas simplicaciones restringe el campo de aplicacin de las ecuaciones resultantes pero simplica notablemente su resolucin. Ha de tenerse en cuenta
que la resolucin del sistema de ecuaciones de Navier-Stokes completo, incluso para
un uido incompresible, es una tarea casi imposible. Las soluciones analticas solo se
pueden obtener para una serie de casos cannicos muy sencillos e incluso su resolucin numrica presenta grandes problemas dada la diferente naturaleza matemtica
que puede adquirir la ecuacin dentro del mismo dominio de resolucin, su no linealidad y la cantidad de trminos acoplados. Muestra de la gran complejidad que
supone el problema, es el enunciado de uno de los problemas del milenio: demostrar
si a partir de unas condiciones iniciales de ujo laminar, la solucin del ujo para
todos los instantes de tiempo es tambin un ujo laminar.
2.2.2. Naturaleza del trmino de presiones

El papel de las presiones en las ecuaciones de Navier-Stokes para densidad constante merece tambin un comentario especial.
Primero se observa que las tensiones normales tienen un efecto similar al de la
gravedad. Normalmente estamos acostumbrados a tratar la presin como una variable termodinmica relacionada con la densidad mediante una ecuacin de estado.
24
Sin embargo, para un uido de densidad constante no hay conexin alguna entre
presin y densidad y debe darse por tanto un enfoque distinto al campo de presiones.
De este modo, si se aplica el operador divergencia a la ecuacin de momento del
sistema de Navier-Stokes, sin asumir que el campo de velocidades es solenoidal, pero
asumiendo que
es el ratio de dilatacin (
~ ~v ),
=
el resultado es:

D
2
= R
Dt
(2.11)
donde:
1
vi vj
R = 2 p
xj xi
Considerando la solucin inicial y de contorno
si,
= 0.
p
es cero en todo el dominio, lo que implica que
(2.12)
La solucin es
=0
si, y solo
satisface la siguiente ecuacin
de Poisson:
2 p =
vi vj
xj xi
(2.13)
De lo que se concluye que para que un campo de velocidades permanezca solenoidal,

el campo escalar de presiones debe satisfacer una ecuacin de Poisson.
2.2.3. La ecuacin de la vorticidad

Una propiedad esencial de los ujos turbulentos es que son rotacionales. La vorticidad
se dene como el rotacional del campo vectorial de velocidades:
~ ~v
~ =
y es igual al doble de la velocidad de rotacin del uido en un punto
un instante
t.
(2.14)
x, y , z
para
Lo que puede ser muy til a la hora de analizar zonas problemticas
dentro del ujo. Adems, algunos mtodos numricos discretizan la ecuacin de la

vorticidad en vez de las de Navier-Stokes.
La ecuacin de la vorticidad se obtiene aplicando el operador rotacional a las
ecuaciones de Navier-Stokes:
25
Di
vi
= 2 i + j
Dt
xj
(2.15)
Donde el trmino de presin se anula para uidos de densidad constante.

Teniendo en cuenta que la ecuacin que expresa la evolucin de un elemento de
lnea material innitesimal
~s
viene dada por:
Dsi
vi
= sj
Dt
xj
(2.16)
Que, salvo por el trmino viscoso, es idntica a la ecuacin de la vorticidad.

Llegados a este punto, se puede concluir que la vorticidad en un ujo no viscoso
se comporta como un elemento lineal innitesimal, sin disiparse (teorema de Helmholtz). Se observa claramente, aunque tambin sucede para ujos viscosos, que si la
vorticidad inicial es nula, no se produce nueva vorticidad.
Si el ratio de tensiones producido por los gradientes de velocidad acta estirando el elemento de lnea alineado con el vector vorticidad, entonces la magnitud de
este vector aumenta correspondientemente. Este fenmeno es conocido como estiramiento de vrtices (vortex stretching ), que es un proceso importante en los ujos
turbulentos y por el cual al estirarse el vector vorticidad este se alarga en la direccin
del ujo. Al estirar el elemento de lnea, la velocidad en la direccin del elemento de
lnea est aumentando.
Si se particulariza, se tiene para ujos bidimensionales que el trmino de estiramiento de vrtices desaparece:
vi
=0
xj
(2.17)
y cada componente de la vorticidad pasa a comportarse como un escalar conservativo. Dada la ausencia de estiramiento de vrtices, la turbulencia bidimensional es
cualitativamente diferente de la tridimensional y, por lo tanto, no reeja elmente
lo que suele suceder en la realidad.
2.2.4. Ecuaciones de gobierno para ujos multifase

Los ujos multifase normalmente se caracterizan por la presencia de interfases
que separan y acoplan al mismo tiempo el movimiento de los dos uidos. Cuando
26
existe un ujo msico a travs de una supercie
que separa la fase 1 de la fase 2,
como por ejemplo cuando hay cambio de fase lquido-vapor, la conservacin de la

masa requiere:
m
2 (v~2 v~i ) ~n = 1 (v~1 v~i ) ~n
Donde
~n
(2.18)
es el vector unitario normal a la supercie entre los dos uidos y
v~i ~n
es
la velocidad ortogonal a esta supercie. Si representamos esta supercie como:
S (~x, t) = 0
Para un tiempo
t + dt,
tendremos
(2.19)
S (~x + v~i dt, t + dt),
S
~ =0
+ v~i S
t
~n =
S=0
en
Y el vector unitario normal, dirigido desde
S<0
con lo que:
hacia
S > 0,
~
S
~
|S|
(2.20)
viene dado por:
(2.21)
De modo que:
~n v~i =
1 S
~
t
|S|
(2.22)
S = 0 hace referencia a una supercie impermeable, como en el caso de una pared,

m
= 0 dando lugar a ~n ~v = ~n v~i . Entonces, uniendo las ecuaciones 2.18 y 2.22, se
Si
obtiene la conocida condicin de contorno cinemtica:
S
~ =0
+ ~v S
t
en
S=0
(2.23)
En las supercies slidas, para un ujo viscoso, normalmente se impone la condicin de no deslizamiento, lo que requiere que la velocidad tangencial del uido se
corresponda con la del contorno:
~n (~v v~i )
27
(2.24)
Combinando las ecuaciones 2.20 y 2.24 se puede demostrar fcilmente que para una
supercie impermeable:
~v = v~i
en
S=0
(2.25)
En el caso de dos uidos separados por una supercie, debe cumplirse la ecuacin
2.23, sin embargo en vez de denir la ecuacin 2.24, es ms indicado imponer continuidad en las velocidades tangenciales:
~n (v~1 v~2 ) = 0
Cuando
m
= 0,
(2.26)
combinando las ecuaciones 2.23 y 2.26 se impone la continuidad de
las velocidades a lo largo de toda la interfase:
v~1 = v~2
en
S=0
(2.27)
Cuando la interfase separa un lquido de un gas, los efectos dinmicos de este ltimo sobre el primero pueden tenerse en cuenta simplemente en trminos de presin,
despreciando la viscosidad. En ese caso, solo se impone la condicin 2.23, no as 2.26.
Para contornos slidos con una velocidad prescrita, la condicin 2.25 da lugar a
una solucin bien denida del sistema de ecuaciones de Navier-Stokes 2.5. Sin embargo, para una supercie libre, es necesaria una condicin ms que tenga en cuenta
el balance de cantidad de movimiento, lo que implica que el salto en las tensiones
sobre la supercie debido a
~t = ~n,
combinado con los ujos de cantidad de mo-
vimiento, deben de ser compensados por la tensin supercial.
h
i

~ I ~n~nt = I ~n~nt
~ + ~n
2 1 ~n m
(v~2 v~1 ) =
donde
(2.28)
es la tensin supercial entre dos uidos y:
~ ~n
=
Es la curvatura de la supercie. Cabe sealar tambin que
(2.29)
I ~n~nt
es el proyector
sobre el plano tangente a la interfase. En la prctica, es ms conveniente descomponer la ecuacin 2.28 en sus componentes normal y tangencial. Para la componente
normal:

p2 + p1 + ~n 2 1 ~n m
(~v2 ~v1 ~n = )
Mientras que para la componente tangencial:
28
(2.30)

~
~n 2 1 ~n = ~n ~n
(2.31)
Consideremos ahora un cuerpo rgido inmerso ene l uido con masa

inercia
J b,
volumen
Vb
y supercie
Sb .
mb ,
tensor de
Acorde a la segunda ley de Newton, el movi-
miento de este cuerpo viene dado por:
d
(mb~vb ) = F~ h + F~ e + mb~g
dt
(2.32)
Y para el momento angular:
d ~ ~h ~e
J b b = L + L
dt
Siendo en este caso
~vb
(2.33)
la velocidad del centro de gravedad del cuerpo,
cidad angular respecto al centro de gravedad tambin, y

los momentos respectivamente; los superndices
F~
~
L
~b
su velo-
son las fuerzas y
hacen referencia a las fuerzas
hidrodinmicas y a otras, de origen externo. Las hidrodinmicas vienen denidas

como:
~h
Z
~ndSb ,
F =
~h
L =
Sb
Donde
~x

~x ~n dSb
(2.34)
Sb
~n
se
el centro de masa, y la orientacin del cuerpo
se mide desde el centro de gravedad y el punto en el que el vector
dirige hacia el exterior de la supercie.

La evolucin temporal de
~,
X
vienen denidas por las relaciones cinemticas:
~
d
~b
=
dt
~
dX
= v~b ,
dt
(2.35)
2.2.5. Naturaleza de las ecuaciones de Navier-Stokes

Las ecuaciones de ujo representan un balance entre trminos convectivos y trminos difusivos, adems de la inclusin de trminos fuente. En los ujos difusivos
aparecen trminos con derivadas de segundo orden como consecuencia de la viscosidad. En los ujos convectivos aparecen derivadas de primer orden que expresan
propiedades de transporte o arrastre. Cada una de estas contribuciones tiene inuencia en la naturaleza matemtica de las ecuaciones, dando lugar a tres tipos de
sistemas de ecuaciones en derivadas parciales: elpticos, parablicos e hiperblicos.
Desde el punto de vista fsico, podemos distinguirlos a su vez segn dos comportamientos clsicos: problemas de contorno o problemas de condicin inicial, segn si el
problema es un problema de equilibrio o estacionario en el primer caso o si se trata
29
de un problema transitorio, en el segundo caso.

Desde un punto de vista estrictamente matemtico, las ecuaciones diferenciales en
derivadas parciales se pueden clasicar, a partir de una denicin genrica, en coordenadas
de las mismas como:
2u
u
2u
2u
u
+
c
+e
+ fu + g = 0
+
b
+d
2
2
x
xy
y
x
y
(2.36)
La clasicacin de estas ecuaciones est determinada por el comportamiento de sus

derivadas de orden superior, por lo que nicamente se considerarn estos trminos
y sus coecientes. Adems se supone que los coecientes
constantes o a lo sumo funciones de las coordenadas
de la variable
son
x e y , independientes por tanto
u.
Para garantizar la continuidad de la primera derivada de
u/y ,
a, b, c, d, e, f
u, ux u/x
uy
escribimos:
ux
2u
2u
ux
dx +
dy =
dx
+
dy
x
y
x2
xy
uy
uy
2u
2u
duy =
dx +
dy =
dx + 2 dy
x
y
xy
y
dux =
Donde
xy
(2.37)
conforma una supercie de soluciones por encima o por debajo del plano
y la pendiente
dy/dx
se conoce como curva caracterstica.
Todo lo anterior se puede combinar para escribir el siguiente sistema matricial:
a b c
uxx
h
dx dy 0 uxy = dux
0 dx dy
uyy
duy
(2.38)

u
u
+e
+ fu + g
h= d
x
y
(2.39)
donde:
Dado que es posible que existan discontinuidades en la segunda derivada de la variable dependiente a lo largo de las caractersticas, estas derivadas son indeterminadas.
Esto sucede cuando el determinante de la matriz de coecientes es igual a cero.

a b c

dx dy 0 = 0

0 dx dy
30
(2.40)
Lo que da lugar a:

a
dy
dx
2

b
dy
dx

+c=0
(2.41)
Resolviendo la ecuacin cuadrtica se obtiene la ecuacin de las caractersticas en

el espacio fsico:
dy
b
=
dx
Dependiendo del valor del discriminante
b2 4ac
2a
b2 4ac,
(2.42)
las curvas caractersticas pueden
ser reales o imaginarias. En los problemas en los que las curvas caractersticas son
reales, una perturbacin se propaga en una regin nita. La zona afectada se conoce
como regin de dependencia o zona de inuencia.
La ecuacin en derivadas parciales de segundo orden se clasica en funcin del
valor del discriminante:
Elptica para
b2 4ac < 0.
Parablica para
b2 4ac = 0.
Hiperblica para
b2 4ac > 0.
Es interesante hacer referencia a la homologa entre la ecuacin en derivadas

parciales y la ecuacin general de una seccin cnica:
a x2 + b xy + c y 2 + d x + e y + f = 0
(2.43)
donde se pueden identicar las siguientes propiedades geomtricas (ver Figura 2.2)
b2 4ac < 0
elipse.
b2 4ac = 0
parbola.
b2 4ac > 0
hiprbola.
En CFD, los esquemas numricos elegidos y las especicaciones de las condiciones de contorno dependen fuertemente del tipo de ecuacin en derivadas parciales
que tratemos. En muchos casos, las ecuaciones de gobierno son una mezcla de todos
los tipos.
31
Figura 2.2: Diferentes secciones cnicas en funcin del valor del discriminante.
2.2.5.1. Comportamiento elptico

La ecuacin elptica ms bsica es la ecuacin de Laplace que, expresada en un
dominio bidimensional, se obtiene de la ecuacin general cuando se ja que
a=c=1
y se igualan a cero el resto de coecientes. Por tanto, establece que:
2u 2u
+
=0
x2 y 2
(2.44)
La caracterstica ms importante de esta ecuacin es que cualquier perturbacin en

el interior del dominio conlleva un cambio en todos los puntos restantes del dominio, tal y como se observa en la Figura 2.3. Esto signica que las perturbaciones se
pueden transmitir en todas las direcciones del espacio. Si existiese simetra completa, la perturbacin se propagara idnticamente en todas las direcciones, siendo el
comportamiento de la ecuacin de tipo esfrico.
Una consecuencia directa de lo anterior es que las condiciones en las fronteras
afectan a la solucin en todos los puntos del dominio. Por tanto, en el caso de que no
existiesen trminos fuente, el valor de las variables tiene que estar necesariamente
acotado por los valores mximos en las condiciones de contorno.
Adems, puesto que cualquier perturbacin se propaga en todas las direcciones,
la solucin de todo problema fsico descrito por una ecuacin elptica es continua y
32
Figura 2.3: Dominio de dependencia de las ecuaciones elpticas.
suave, incluso si existen discontinuidades en las condiciones de contorno.

Estas dos ltimas propiedades son un buen criterio para comprobar que el mtodo numrico empleado para resolver una ecuacin elptica es apropiado. Por tanto,
es evidente que para que la informacin se propague en todas las direcciones, las
tcnicas numricas utilizadas en problemas elpticos deben permitir que cualquier
variacin en un punto del dominio pueda inuir sobre todos sus nodos vecinos.
2.2.5.2. Comportamiento parablico

La ecuacin parablica ms bsica es la ecuacin de difusin en modo no estacionario, que expresada en un dominio bidimensional, se obtiene de la ecuacin general
cuando se ja que
a =
e=1
y se igualan a cero el resto de coecientes. Por
tanto, establece que:
2u
u
= 2
t
x
(2.45)
La caracterstica fundamental de esta ecuacin es que cualquier perturbacin

en el interior del dominio slo puede tener inuencia en instantes posteriores a la
33
aparicin de la perturbacin, tal y como se muestra en la Figura 2.4, propagndose

nuevamente en todas las direcciones del espacio. Adems, la aparicin de la derivada temporal exige la denicin de unas condiciones iniciales que tendrn inuencia
en los valores futuros de la variable. Lgicamente, las condiciones iniciales afectan
tambin a la solucin en todos los puntos del dominio, aunque su impacto ir decreciendo conforme el tiempo avanza.
Figura 2.4: Dominio de dependencia de las ecuaciones parablicas.
La presencia del trmino difusivo conlleva que las condiciones de frontera afecten
a la solucin en todos los puntos del dominio, al igual que en la ecuacin elptica
(ntese que la ecuacin parablica se reduce a la elptica cuando la derivada temporal se anula). Por tanto, en ausencia de trminos fuente, el valor de las variables
debe estar acotado por los valores mximos de las condiciones de contorno y de las
condiciones iniciales.
Las ecuaciones parablicas pueden tender a un estado estacionario cuando
en el que la solucin sea independiente de los valores iniciales y por tanto se
recupere el carcter elptico. De todas formas, queda claro a la vista de la ecuacin

fundamental que en este tipo de problemas la variable temporal se comporta de manera muy diferente respecto del resto de las variables espaciales. La variacin en
slo admite inuencias en un sentido, mientras que la variable
34
x permite variaciones
en los dos sentidos, esto es, propagaciones en todas las direcciones del espacio en 3D. Por esta razn, la variable
t se conoce habitualmente como dimensin transitoria
o parablica.
2.2.5.3. Comportamiento hiperblico

La ecuacin hiperblica ms bsica es la ecuacin de onda, caracterstica en problemas vibratorios u ondulatorios y que, expresada en un dominio bidimensional
2
(t, x), se obtiene de la ecuacin general cuando se ja que a = y c = 1 y se
igualan a cero el resto de los coecientes. Por tanto, establece que:
2
2u
2 u
=
t2
x2
donde
(2.46)
representa la velocidad de propagacin de una perturbacin en el seno
de un uido. Este tipo de ecuaciones aparece en problemas no estacionarios con

fenmenos de disipacin prcticamente despreciables.
Figura 2.5: Dominio de dependencia de las ecuaciones hiperblicas.
Adems de la ecuacin de la onda y de fenmenos vibratorios anes, como el

caso de los modos de vibracin de una cuerda, el comportamiento hiperblico est
35
ntimamente ligado al caso de ujo compresible en rgimen transnico y supersnico.

Para este tipo de ujos, la velocidad del sonido (es decir, la velocidad de propagacin de las perturbaciones en el medio) puede dar lugar a la aparicin de ondas de
choque, que son claras manifestaciones de la naturaleza hiperblica de sus ecuaciones de gobierno. Las ecuaciones no viscosas de ujo compresible son hiperblicas y
ponen de relieve que las perturbaciones en un punto nicamente pueden alcanzar
zonas muy limitadas en el espacio.
2.2.5.4. Clasicacin de las ecuaciones del ujo

La ecuacin de conservacin de momento tiene mucho en comn con los distintos
tipos de comportamiento presentados anteriormente para las ecuaciones diferenciales
en derivadas parciales. As, por ejemplo, bastara con imponer condiciones estacionarias en ausencia de ujo para plantear una difusin pura con comportamiento
elptico. El mismo problema, pero en condiciones no estacionarias, mostrara un
comportamiento parablico. Por otro lado, el trmino convectivo de la ecuacin tiene un claro carcter hiperblico.
En la mayora de las aplicaciones para la ingeniera, la ecuacin presenta un comportamiento mixto, con trminos difusivos dndole un carcter elptico y trminos
temporales y convectivos otorgndole caractersticas parablicas e hiperblicas. La
correcta clasicacin de la ecuacin puede ser de inters puesto que en funcin del
comportamiento caracterstico de la ecuacin existen distintos mtodos numricos
de discretizacin y resolucin de las ecuaciones.
Las ecuaciones de Navier-Stokes, as como sus expresiones simplicadas, se pueden
clasicar segn los comportamientos vistos con anterioridad, tal y como se muestra
en la tabla 2.4.
Tipo de ujo
Flujo estacionario
Flujo no estacionario
Flujo viscoso
Elptico
Parablico
Flujo ideal
(no viscoso)
Flujo en capa lmite
M a < 1: Elptico
M a > 1: Hiperblico
Si
Si
Parablico
Hiperblico
Hiperblico
Parablico
Tabla 2.4: Clasicacin de los principales tipos de ujo.
La clasicacin de la tabla es una clasicacin formal de las ecuaciones de ujo

puesto que, en la prctica, la mayora de los ujos se comporta de forma compleja. De todas formas, queda claro que las ecuaciones de Navier-Stokes estacionarias
son formalmente elpticas, mientras que su formulacin no estacionaria es parablica.
La clasicacin matemtica de las ecuaciones para ujo ideal es notablemente
diferente debido a la falta de los trminos viscosos o disipativos de orden superior.
36
En este caso, la clasicacin depende en realidad de la importancia de los fenmenos de compresibilidad y por extensin del valor del nmero de Mach (M a). As,
la naturaleza elptica de los ujos ideales para nmeros de Mach menores que la
unidad se debe a la accin de los gradientes de presin: la presin puede propagar
perturbaciones a velocidades snicas, que son mayores que la velocidad del ujo. Por
el contrario, si la velocidad del ujo es supersnica (M a
> 1)
la presin es incapaz
de propagarse a las regiones aguas arriba, limitando las zonas de inuencia de las
perturbaciones y dando el carcter hiperblico a las ecuaciones de gobierno.
Finalmente, en el caso de ujos con elevados nmeros de Reynolds, en los que
las zonas de inuencia viscosa ocupan una extensin muy reducida dentro del problema estudiado y el resto de zonas pueden considerarse como uido ideal, todas
las derivadas de la velocidad en la direccin del ujo son mucho ms pequeas que
las derivadas en la direccin transversal. Capas lmite, chorros, capas de cortadura,
estelas y ujos en conducto totalmente desarrollados son ejemplos que responden a
este patrn. Es estas condiciones, las ecuaciones de gobierno contienen nicamente
un trmino de difusin, lo que resulta en un comportamiento claramente parablico.
Se podra concluir por tanto, que los fenmenos de difusin actan en todo el dominio, independientemente de la direccin predominante del ujo (comportamiento
elptico). Los fenmenos de conveccin actan en la direccin de propagacin, en regiones concretas del espacio (comportamiento hiperblico). Entre ambas situaciones,
una ecuacin con comportamiento parablico representa una situacin intermedia
que se puede interpretar como un proceso de difusin en todas las direcciones, pero
amortiguado en el tiempo.
2.2.6. Ecuaciones de Reynolds o del ujo medio

Como es bien sabido, los ujos reales pueden presentarse en dos regmenes completamente diferentes: bien en rgimen laminar, bien en rgimen turbulento. Los
ujos laminares se caracterizan porque las partculas de uido se mueven de manera
ordenada, en lminas, en la direccin de la corriente y sin generar ujos ni movimientos transversales. Este tipo de ujos aparece cuando la viscosidad es el fenmeno de
transporte dominante, lo cual permite atenuar cualquier tipo de perturbacin que
pudiera afectar al ujo. Se caracterizan por presentar un nmero de Reynolds bajo.
Por el contrario, los ujos turbulentos se caracterizan por presentar uctuaciones aleatorias de los campos de presiones y velocidad tanto en el tiempo como en
el espacio tal y como se observa en la Figura 2.6. Estas uctuaciones provienen de
unas inestabilidades del ujo que en forma de cascada de energa generan vrtices
cada vez ms pequeos hasta que se disipan en forma de calor por la accin de la
viscosidad. Esto puede encontrarse ms desarrollado en el Cap. 2.2.8.
37
Tpicamente, la evolucin de la variable uidodinmica en un punto del espacio

est compuesta por una serie de uctuaciones de distinta escala en el caso del ujo
turbulento, que es normalmente el tipo de ujo habitual. Es curioso observar cmo
las ecuaciones de gobierno son ecuaciones exactas, estrictamente formuladas en su
sentido matemtico y, sin embargo, sus soluciones reales dan lugar a fenmenos turbulentos, aleatorios, que slo tienen sentido prctico si se analizan desde un punto
de vista estadstico como valor medio de la velocidad e intensidad promediada de
las uctuaciones.
En concreto, para poder cuanticar el nivel de oscilacin de la velocidad es necesario realizar un promedio de la misma, con el objeto de determinar su valor medio
y la intensidad y frecuencia de las oscilaciones. De esta forma, es habitual denir el
denominado promediado de Reynolds como:
1
u = lm
T + T
t0 +T
u(t)dt
(2.47)
t0
donde la desviacin instantnea de la variable en el tiempo respecto de su valor

medio, es decir la uctuacin turbulenta se dene como:
u0 = u u
(2.48)
Figura 2.6: Velocidad en un ujo turbulento y turbulencia instantnea en el punto
tras el cilindro.
38
Partiendo del sistema de ecuaciones 2.10, y tomando el valor medio, para la

ecuacin de continuidad se obtiene:

~ ~v =
~ ~v + v~0 = 0
(2.49)
Mantenindose para el valor medio:
~ ~v = 0
(2.50)
~ v~0 = 0
(2.51)
Por diferencia se obtiene que:
Para las ecuaciones de cantidad de movimiento, las operaciones no son tan sencillas.
Para el trmino convectivo (no lineal), lo primero es escribir la derivada material en
su forma conservativa:
vj
Dvj
=
+
(vi vj )
Dt
t
xi
(2.52)
De modo que tomando el valor medio:
vj
Dvj
=
+
(vi vj )
Dt
t
xi
(2.53)
As, y sustituyendo segn la descomposicin realizada en 2.48:
vi vj = (vi + vi0 ) vj + vj0

= vi vj +
vi0 vj
= vi vj + vi0 vj0
vj0 vi
(2.54)
vi0 vj0
(2.55)
(2.56)
La convarianza de la desviacin de velocidades se les conoce como tensiones de Reynolds:
39
R = vi0 vj0
(2.57)
Sustituyendo 2.56 en 2.53:
Dvj
Dt

vj
+
vi vj + vi0 vj0
t
xi
0 0
vj
vj
+ vi
+
=
vv
t
xi xi i j
(2.58)
(2.59)
El segundo paso resulta de aplicar 2.50. Se puede ahora reescribir la media de la

derivada material como:
Dvj
0 0
Dvj
=
+
vv
Dt
Dt
xi i j
(2.60)
Que no es idntica a la derivada material media, que se podra denir como:
+ vi
t
Dt
(2.61)
Ahora es sencillo tomar la media de la ecuacin de cantidad de movimiento, dado

que los dems trminos son lineales en
vi
p.
El resultado son las ecuaciones de
Reynolds o de cantidad de movimiento media:
i
1 p
0 0
Dv
2
=
g
vv
j
j
xj
xi i j
Dt
(2.62)
En apariencia, las ecuaciones de Navier-Stokes 2.10 y las promediadas de Reynolds

2.62 son muy similares, excepto por las tensiones de Reynolds.
Las tensiones de Reynolds
vi0 vj0
juegan un papel crucial en las ecuaciones para la
v . Si vi0 vj0 se anularan, las soluciones para

velocidad media ~
~v
~v
seran idnticas.
2.2.7. Problema del cierre

Para un ujo tridimensional, existen cuatro ecuaciones independientes para la
velocidad media: tres componentes de las ecuaciones de Reynolds y una ecuacin
40
para continuidad media o una ecuacin de Poisson para la presin media

embargo, estas ecuaciones contienen ms de 4 incgnitas. Adems de
~v
p,
p.
Sin
estn
las tensiones de Reynolds.

Este es el enunciado del problema del cierre. Consecuentemente, a falta de informacin sobre las tensiones de Reynolds, el sistema no puede ser resuelto.
2.2.8. Turbulencia
La turbulencia es un estado catico e irregular del movimiento de un uido que
se establece a partir de la aparicin de irregularidades en las condiciones iniciales o
de contorno de la corriente uida. Estas inestabilidades se amplican y se retroalimentan de forma cclica, creando vrtices turbulentos.
No es fcil denir exactamente qu es la turbulencia, sin embargo todo el mundo
tiene una nocin intuitiva de lo que es: un movimiento uctuante y desordenado. Por
lo tanto, parece acertado armar que la turbulencia es un estado catico y aleatorio
del ujo, en el que la velocidad y la presin oscilan instantneamente a lo largo del
espacio y el tiempo. Este tipo de comportamiento ha propiciado tradicionalmente
un enfoque estadstico del fenmeno.
La turbulencia es una caracterstica de los ujos, no de los uidos como tales. Su
aparicin precisa la existencia de un uido en movimiento, en el que los fenmenos
de conveccin (inerciales) asociados a la velocidad sean varios rdenes de magnitud superiores a los efectos difusivos (disipativos) relacionados con la viscosidad del
uido. Recordemos que son estos los trminos no lineales (Cap. 2.5). Esta relacin
entre ujos convectivos y difusivos es el conocido nmero de Reynolds que permite
establecer la frontera aproximada entre las condiciones de ujo laminar y ujo turbulento.
La dinmica de la turbulencia es muy compleja y, dada su naturaleza, requiere
varios niveles de aproximacin para describirla. Por esta razn, no existe una teora
de la turbulencia de aplicacin general, si no mltiples modelos creados especcamente para diferentes problemas.
Sobre el tema de la turbulencia, se pueden encontrar grandes obras maestras
como pueden ser (Pope, 2000; Schlichting, 1979; Hinze, 1975).
Las propiedades ms destacables de los ujos turbulentos son:
Aleatoriedad. Es la caracterstica ms destacada de los ujos turbulentos. A

veces descrita como ruido, se maniesta por la aparicin de uctuaciones de
41
las variables uidodinmicas con tamaos y tiempos muy dispares (diferentes

escalas de la turbulencia). Estas uctuaciones instantneas no estacionarias
se desarrollan incluso en ujos estacionarios, lo cual es muestra de que solo
las propiedades estadsticas de los ujos son invariantes. Esta es la principal
razn por la que se utilizan mtodos estadsticos para su estudio y prediccin.
Adems, la existencia de escalas muy dispares hace que, aunque los ujos turbulentos parezcan caticos e impredecibles, tambin existan estructuras coherentes dentro de ese mar de irregularidades aleatorias.
Vorticidad. Es imprescindible que exista vorticidad para que un ujo pueda

ser turbulento. De hecho, todo ujo turbulento es rotacional ( ~
v
6= 0),
con
importantes niveles de vorticidad que uctan en el tiempo y en el espacio de

forma coherente, y en los que la deformacin de los vrtices supone la esencia
de la dinmica de la turbulencia.
Difusividad. Los fenmenos turbulentos intensican el transporte de masa,

cantidad de movimiento y energa. En particular, las uctuaciones a escalas
macroscpicas producen efectos de mezcla similares a los de carcter molecular
(puramente difusivos), si bien con longitudes de mezcla que dependen de los
ujos convectivos.
Transmisin de energa en forma de cascada. La turbulencia est formada por

un espectro continuo con gran rango de escalas. Con objeto de ilustrarlo, se
puede pensar en torbellinos: los ms grandes transmiten energa a los ms pequeos como si de una cascada se tratara y en estos se disipa a travs de la
disipacin molecular convirtindose en calor. Dado que la dinmica de las pequeas escalas es muy rpida en comparacin con el movimiento de las grandes
escalas, se puede asumir razonablemente que el movimiento en las pequeas
escalas es independiente de las grandes. Asimismo, la disipacin de energa
debe ser a un ratio cercano al aporte de energa, sobre todo en la turbulencia
istropa y estacionaria. Esta es la conocida teora de Kolmogorov del equilibrio
universal.
Existen tambin otras caractersticas del ujo turbulento que conviene destacar,
aunque no sean tan fundamentales como las anteriores. Por ejemplo:
Tridimensionalidad. Las escalas ms pequeas de la turbulencia tienen un

carcter muy isotrpico, lo cual implica la necesidad de modelar el ujo tridimensional, nunca bidimensional, tal y como ha quedado de maniesto en el
Cap. 2.2.3. Las escalas ms grandes, asociadas a las longitudes caractersticas del ujo analizado, pueden presentar un comportamiento bidimensional o
42
plano, pero ste se va generalizando a tridimensional segn se avanza en la

cascada de energa.
Disipacin. Los ujos turbulentos son siempre disipativos. Necesariamente han

de disipar energa en las escalas ms pequeas, energa que se obtiene del ujo
principal y que se va redistribuyendo en forma de cascada. Una vez desarrollado el ujo turbulento, la turbulencia tiende a mantenerse mediante un aporte
continuo de energa desde las escalas ms grandes. Si no existe ese suministro
de energa, la turbulencia decae.
Se dice que un fenmeno fsico es estable cuando, ante pequeas perturbaciones,

este permanece invariante, imperturbable en el tiempo. Una modicacin razonablemente pequea de las condiciones iniciales no altera signicativamente la situacin
de equilibrio.
Un ujo laminar es estable ante pequeas perturbaciones. Cuando esto no sucede,
perturbaciones innitesimales crecen espontneamente y se amplican modicando
completamente el estado inicial. Estas perturbaciones generan nuevas perturbaciones, que a su vez producen nuevas perturbaciones en un proceso degenerativo que
destruye la estructura original del ujo hasta establecer un ujo turbulento con
multitud de escalas diferentes. Esta dinmica de la inestabilidad es el origen de la
turbulencia y ocurre cuando el nmero de Reynolds es sucientemente grande.
Por tanto, la turbulencia es un fenmeno muy complejo y vortical, con escalas
muy dispares que van desde el tamao caracterstico del ujo hasta escalas disipativas muy pequeas. An siendo pequeas, estas escalas disipativas estn lejos de las
escalas de longitud molecular por lo que an se pueden emplear las ecuaciones de la
mecnica de uidos para medio continuo.
2.3. Hidrulica de grandes presas

Las bases cientcas empleadas en la actualidad podran quedar perfectamente
recogidas en los tratados ya publicados a nales de los 70 (USBR, 1977; Peterka,
1958), con avances casi cada 20 aos que han marcado importantes hitos cientcos
(Hager, 1992; Vischer and Hager, 1995). Por otro lado, los pilares del anlisis de
los ujos agua-aire y los problemas de cavitacin propios de las grandes estructuras
hidrulicas seran estudiados un poco ms tarde (Falvey, 1980, 1990; Wood, 1991;
Chanson, 1996), habindose convertido la aireacin del ujo en un asunto central,
como as lo revela la ferviente actividad investigadora de los ltimos aos (Pster
and Hager, 2010a; Pster et al., 2006; Ma et al., 2011b; Meireles et al., 2014; Kramer et al., 2006; Chanson, 2009; Bung, 2013; Oertel and Bung, 2012; Bombardelli
43
et al., 2011; Bung, 2011), pudiendo destacar los dos Vision papers recientemente
publicados por el Journal of Hydraulic Research sobre ujos agua-aire (Chanson,
2013) y estructuas hidrulicas (Hager, 2014) y que ponen todava ms de maniesto
el inters actual de la comunidad cientca por las grandes obras hidrulicas.
El diseo de los rganos de desage ha ido siempre asociado a modelos fsicos

reducidos. Sin embargo, los efectos de escala asociados al modelo reducido dejan
siempre varios aspectos sin resolver (Heller, 2011; Chanson, 2009). Este es el caso
de la aireacin del ujo. Para los nmeros de Weber y de Reynolds asociados a los
modelos reducidos, los efectos de escala no permiten una correcta reproduccin de la
aireacin del prototipo (Kramer et al., 2006; Pster and Hager, 2010a; Pster et al.,
2006). La importancia de la aireacin del ujo, tanto por la proteccin que ofrece
al aliviadero (Falvey, 1990; Kramer, 2004) como por la reduccin en la friccin de
fondo (Madavan et al., 1983; Wood, 1991; Gutirrez Serret, 1995; Chanson, 1996),
convierte a la concentracin de aire en una variable de suma importancia, al mismo
nivel que la velocidad o la presin. De hecho, esta concentracin de aire modica las
propiedades macroscpicas del ujo y da lugar a distintas distribuciones de cantidad
de movimiento. Adems se han encontrado importantes diferencias en el comportamiento de la turbulencia debido al efecto de la fase dispersa (aire) sobre la fase
portante (agua) (Balachandar and Eaton, 2010; Crowe, 2000; Kussin and Sommerfeld, 2002; Zhou, 2010; Yeo et al., 2010; Crowe, 2010). Esto, sin ninguna duda, afecta
las condiciones de entrada al cuenco y puede invalidar un diseo hidrulico que por
otro lado podra ser perfectamente aceptable para un ujo no aireado. Por ello, es
interesante emplear en el diseo otro tipo de tcnicas que permitan una aproximacin al ujo producido en el prototipo, a escala 1:1.
En los ltimos quince aos los estudios hidrulicos han ido, cada vez ms, acompaados de tcnicas CFD que permiten un anlisis detallado del ujo. Estas tcnicas
han avanzado mucho desde las primeras aplicaciones hidrulicas (Caisley et al., 1999;
Ho et al., 2001; Savage and Johnson, 2001), permitiendo actualmente simulaciones
de grandes obras hidrulicas con modelos para ujos bifase. El empleo de las tcnicas CFD puede permitir el anlisis rpido y able de mltiples alternativas, con el
aliciente del nuevo enfoque Computational Multiphase Fluid Dynamics que permite
acoplar el efecto del aire sobre el agua del mismo modo que el modelo de turbulencia
acopla el efecto de las pequeas escalas sobre el ujo medio.
En la hidrulica avanzada de grandes presas, pueden distinguirse dos regiones

muy diferenciadas de ujo: la zona de aliviadero y el cuenco de amortiguacin. Dentro del ujo sobre aliviadero, se suelen distinguir varias regiones: regin no aireada,
con aireacin gradual y completamente aireada. Si adems se usan aireadores, el
ujo puede llegar a ser ms complejo.
44
2.3.1. Aliviadero
El ujo sobre un aliviadero comprende dos regiones bien diferenciadas:
La regin no aireada, cercana a la cresta, donde el ujo es suave. Acaba con

el punto crtico de inicio de aireacin.
La regin aireada. Bien naturalmente o articialmente usando aireadores. El
ujo se airea rpidamente aumentando en calado. Esta regin se puede distinguir a simple vista pues el ujo se vuelve de un blanco opaco.
Al existir dos tipos de ujo completamente distintos, son necesarios modelos distintos para cada una de las partes (Ma et al., 2010; Bombardelli et al., 2011; Meireles
et al., 2014). Es por ello de vital importancia delimitar con precisin la zona en la
que tiene lugar cada una de las dos regiones. Dada la gran turbulencia que tiene
lugar cerca del punto crtico de aireacin, esta no es una tarea sencilla.
Desde el punto de vista de la mecnica de uidos, es el ujo en la regin aireada
el que representa un mayor reto.
2.3.1.1. Supercie libre

Conocer el perl de la lmina libre en una estructura de desage es importante
para asegurar los resguardos. Empleando las curvas facilitadas por el USBR (Peterka, 1958) se puede obtener la velocidad en el aliviadero; y por continuidad el calado.
Figura 2.7: Esquema aliviadero.
45
Otro mtodo ms general es el presentado por (Hager, 1991a). Siguiendo la notacin

de la gura 2.7, donde
y
z0
ja,
es el caudal desaguado,
hacen referencia a la posicin del labio,
el talud,
son los ejes coordenados,
x0
es el punto de transicin a pendiente
la altura de energa del embalse sobre labio jo y
la altura
de la lmina libre sobre el aliviadero; se puede denir el perl de la lmina libre

adimensionalizado como
S(X),
con
X = x/HD
siendo:
S = s/HD
Se dene
(2.63)
como la altura de energa adimensionalizada:
= H/HD
As, para el rango
2 < X/1,1 < 2
(2.64)
se puede aproximar el perl de la lmina en el
entorno de la cresta como :

S = 0,75 1,1 (1/6) X
(2.65)
2.3.1.2. Perl de velocidades

Para cualquier punto del aliviadero en el que el ujo no se encuentre aireado,
se puede emplear el baco de la gura 2.8 y aproximar el perl por uno potencial,
segn la ecuacin 2.75. Tal y como se comprueba en la fugura 2.18, la aproximacin
es razonablemente buena.
46
Figura 2.8: Clculo de la velocidad en un aliviadero segn USBR (Peterka, 1958).
La velocidad sobre la cresta
Vc (z)
es de especial inters. Se puede comprobar como
el perl se aproxima razonablemente bien por:

Vc
1 r + 0,75
=
2 r + z/H
2gH
Siendo
r ' 0,584
la curvatura adimensional de la cresta.
47
(2.66)
Figura 2.9: Campo de velocidades en el entorno de la cresta para
48
= 0,5.
= 1,0.
= 1,5.
= 2,0.
La velocidad sobre la cresta
Vc (z)
es de especial inters. Se puede comprobar como
el perl se aproxima razonablemente bien por:

Vc
1 r + 0,75
=
2 r + z/H
2gH
Siendo
r ' 0,584
la curvatura adimensional de la cresta.
49
(2.67)
2.3.1.3. Turbulencia
No existen muchos datos sobre turbulencia en aliviaderos. Aunque existen varias evidencias de que la turbulencia cambia sustancialmente cuando existe una fase
dispersa (Chanson, 2013, 2009; Prosperetti and Tryggvason, 2007; Crowe, 2010; Balachandar and Eaton, 2010; Yeo et al., 2010; Zhou, 2010; Kussin and Sommerfeld,
2002; Crowe, 2000), suele tomarse como referencia los valores propuestos por (Falvey
and Ervine, 1988):
v0 = U
Siendo
(2.68)
la velocidad de corte denida como:
U =
p
gRh S0
(2.69)
Por otro lado, (Nezu et al., 1993) propone un valor similar:
v 0 = 1,27U
(2.70)
De modo que la intensidad turbulenta quedara:
Ti =
1,27U
v0
=
v
v
(2.71)
2.3.1.4. Determinacin del punto de inicio de la aireacin

Los ujos sobre aliviadero pueden acabar dando lugar a los que comnmente se
conoce como aguas blancas", para ello tienen que concurrir dos circunstancias:
El ujo tiene que ser completamente turbulento, por ejemplo: la capa lmite
debe ocupar todo el espesor uido.
La energa cintica de los remolinos ms superciales debe ser superior a la
energa de la tensin supercial.
La ltima condicin se satisface siempre en estructuras a escala de prototipo.
Dependiendo del nivel de turbulencia, el ujo sobre aliviadero se airea naturalmente
50
en prototipo cuando las velocidades exceden de 10 - 15 m/s.
Algunos autores previos a los aos 40 sealaron las altas velocidades como la
causa del inicio de la aireacin. Estas teoras clsicas han demostrado estar obsoletas y no recoger bien la fsica del problema. Otros estudios han investigado con
cmaras de alta velocidad la formacin de burbujas (Volkart, 1980) concluyendo que
son las gotas de agua que escapan de la supercie libre las que atrapan las burbujas
al caer. Esta teora queda contradicha posteriormente (Ervine, 1987). Otros estudios
sugieren que la reaireacin en un aliviadero ocurre cuando la capa lmite intersecta
la lmina libre (Bauer, 1954). Este ltimo concepto es todava muy ampliamente
aceptado. Es necesario sin embargo que esta turbulencia que llega a la lmina libre
sea capaz de deformarla.
Otro criterio para el inicio de la aireacin basado en el valor de la turbulencia

viene expresado por (Ervine, 1987; Chanson, 1993):
r
urms =
Donde
curvatura
4
w R
(2.72)
es la tensin supercial agua-aire, w la densidad del agua y R el radio de

de las ondulaciones de la supercie libre. R puede ser sustituido por db /2
(Chanson, 1993), asumiendo implcitamente que la deformacin de la lmina libre

originara burbujas de un dimetro similar. Tambin aade una segunda condicin,
que hace relacin a la accin de la gravedad:
urms > ur cos

Siendo
ur
(2.73)
la velocidad de ascensin de las burbujas. La turbulencia debe de ser
por lo tanto mayor que la tensin supercial y los efectos de la gravedad.
Si se pretende emplear el criterio de la interseccin de la capa lmite con la lmina libre, entonces el crecimiento de la capa lmite en un aliviadero puede calcularse
segn (Wood et al., 1983):
= 0,0212
x
Donde
xs
Hs
0,11
es el espesor de la capa lmite,
desde la cresta del aliviadero,
Hs
xs
xs
ks
0,10
(2.74)
la coordenada longitudinal medida
la diferencia de nivel entre el punto donde se eva-
la el espesor de la capa lmite y la altura de agua en el embalse, y
51
ks es la rugosidad.
Figura 2.13: Determinacin del punto crtico de aireacin: interseccin de la capa

lmite con la lmina libre.
Este espesor de capa lmite
se dene como la distancia perpendicular desde la
supercie del aliviadero hasta el punto en el que la velocidad

irrotacional
uf s
es el 99 % del ujo
(Wood, 1991). El perl de velocidades puede ser aproximado como:
y 1/n
u
=
uf s
Siendo
(2.75)
la constante del perl potencial. Se recomienda que
tome el valor de 6.3
(Cain and Wood, 1981; Cain, 1978; Wood et al., 1983); otras consideraciones sobre
el valor de este parmetro pueden encontrarse en (Castro-Orgaz, 2009; Cheng, 2007;
Pope, 2000). Esta estimacin del perl de velocidades es comnmente aceptada para
ujos con
Cmean < 0,5
(Chanson, 1993).
xs = xI tiene lugar cuando

lmina libre h = h(x) (g. 2.74).
El punto crtico de inicio de aireacin

la capa lmite
= (x)
intersecta la
La ecuacin 2.74 es aplicable para taludes entre 5
y 70
el espesor de
y es especialmente til
porque muestra la inuencia de cada uno de los factores sobre el crecimiento de la

capa lmite. De este modo, la capa lmite crece con la distancia
y es reducido con la altura de energa
Hs .
xs
y la rugosidad
ks
Esta ecuacin se mantiene vlida para un
aliviadero estndar, sin embargo para un aliviadero de ancho constante y pendiente

gradualmente variable la capa lmite se puede aproximar diferenciando la ecuacin
2.74:

d
d
=
= 0,9 0,11xs cot
dxs
xs
dxs
52
(2.76)
Donde
es la pendiente del aliviadero. Esta expresin puede ser modicada tambin
para permitir cambios graduales en anchura (Wood, 1991).
La descarga (caudal especco) dentro de la capa lmite se puede obtener como:
q =
n p
2gHs
n+1
(2.77)
El calado total del ujo puede por lo tanto calcularse usando la ecuacin de Bernoulli para la velocidad fuera de la capa lmite
d=+
Siendo
uf s
como:
q q
uf s
(2.78)
el caudal especco (caudal por unidad de ancho).
El punto crtico de inicio de la aireacin tiene lugar cuando el caudal especco

completamente desaguado en la regin de capa lmite
es
q . Sin embargo, el eje externo
de la capa lmite es muy irregular (observar g. 2.14) y puede oscilar entre 0.4 y 1.2
veces el valor calculado, con lo que pueden existir regiones aireadas aguas arriba del
punto calculado y regiones sin airear aguas abajo.
53
Figura 2.14: Visualizacin de una capa lmite para ujo con Mach = 2.82. El ujo
se produce de derecha a izquierda. Adaptado de (Smith and Smits, 1995).
2.3.1.5. Concentracin de equilibrio

Deniendo la concentracin media en un perl
Cmean
Siendo
y90
equilibrio
1
=
y90
como:
y=y90
C(y)dy
(2.79)
y=0
el calado para el cual la concentracin
Ce
Cmean
C(y90 ) = 0,90. La concentracin de
se dene como la concentracin media del ujo completamente desa-
rrollado en el aliviadero. Algunos autores realizan denen el lmite superior de la

integral de la ecuacin 2.79 como
y95
y99 .
portancia pues el aumento de calado entre
54
A efectos prcticos, no tiene mayor im-
C = 0,90
C = 0,99
es mnimo.
Figura 2.15: Estructura del ujo supercrtico aireado.
Los primeros estudios experimentales sobre ujos completamente desarrollados

en aliviaderos se pueden encontrar en (Jevdjevich and Levin, 1953; Straub and Anderson, 1958). Ms adelante, (Chanson, 1993; Wood, 1991) los recogen y hacen un
importante trabajo de sntesis de los resultados.
Figura 2.16: a) Concentracin media b) Nomenclatura grcos c) Concentracin

media y de fondo.
Otro estudio interesante sobre ujos naturalmente aireados se puede encontrar en

(Hager, 1991b), donde se establece que la concentracin de equilibrio puede ser expresada como:
55
Ce = 0,75 (sin )0,75
(2.80)
La concentracin de fondo, la que evita los daos por cavitacin, se puede expresar
como:
C0 = 1,25
3
180
0 40o
40o 80o
C0 = 0,65 sin
(2.81)
(2.82)
Los resultados de estas expresiones se muestran en la gura 2.16.

En cuanto a los perles que se desarrollan, se ha escrito mucho en la literatura
(Hager, 1991b; Chanson, 1993; Kramer et al., 2006; Pster et al., 2006). En la gura
2.17.
Figura 2.17: Perl de concentraciones terico y medido en prototipo.
56
Figura 2.18: Perl de velocidades terico y medido en prototipo.
Desde otro enfoque, (Wilhelms and Gulliver, 2005) analiza los datos histricos
y propone una descomposicin de la concentracin en una seccin en dos partes, la
atrapada y la emulsionada. Algunos autores ubican este cambio en el tipo de ujo
en
yt C(yt ) = 0,50
(Pster, 2008). La parte emulsionada es aquella que trans-
porta burbujas dispersas dentro del lete uido. Por otro lado, la parte atrapada es
arrastrada por la rugosidad de la lmina y no est compuesta de burbujas sino de
grandes masas de aire (gura 5.4).
Para la parte arrastrada por el aire, haciendo una similitud entre la rugosidad
de la capa lmite y la de la lmina libre, propone una concentracin constante de
0.23. Eliminando este valor de los datos propone unas nuevas expresiones para la
concentracin de aire emulsionado en el ujo completamente desarrollado
Ce = 0,656 1 e(0,0356(10,9))
Ce
(2.83)
Para la evolucin espacial de la concentracin, y partiendo del punto crtico de aireacin, propone:
Ce = Ce 1 e(0,010X
57
/Y )
i
(2.84)
Siendo
Yi
el calado crtico de inicio de aireacin. Estas expresiones permiten calcular
la evolucin del calado a partir de la posicin del punto crtico (gura 2.19).
Figura 2.19: Evolucin de las concentraciones en un aliviadero.
2.3.1.6. Reduccin de friccin por el efecto de la aireacin

Las burbujas atrapadas a travs de la lmina libre se difunden dentro del lete
uido alcanzando la solera y reduciendo signicativamente la friccin (Wood, 1991;
Chanson, 1994, 1993; Madavan et al., 1983) tal y como se muestra en la gura 2.20.
58
Figura 2.20: Reduccin de friccin por el efecto del aire en el ujo.
Varias correlaciones empricas se han establecido para la reduccin de friccin. As,

analizando los datos de (Jevdjevich and Levin, 1953; Straub and Anderson, 1958;
Aivazyan, 1986), se puede llegar a (Chanson, 1993):

0,490 Ce
fe
= 0,5 1 + tanh 0,70
f
Ce (1 Ce )
Siendo
Ce
fe
el factor de friccin efectivo,
(2.85)
el factor de friccin para ujo monofsico y
es la concentracin de equilibrio. La correlacin para 104 puntos es de
R = 0,898.
con un enfoque ms ingenieril, existe la correlacin:
fe
= 0,307 + 0,1446 log10 (Re) 1,40Cmean
f
(2.86)
Re. En esta formulacin se emplea la

concentracin media en la seccin Cmean . Cuando se emplea en el rango Cmean > 0,25
5
7
y 2 10 < Re < 4 10 la correlacin obtenida es de R = 0,913. Una formulacin
Que incorpora tambin el nmero de Reynolds
ms genrica que funciona mejor para los valores lmite es:

fe
C0,5 Cmean
= 0,5 1 + tanh 0,70
f
Cmean (1 Cmean )
Donde
C0,5
(2.87)
representa la concentracin para la cual el factor de friccin efectivo
59
fe
toma la mitad del valor de
y viene expresado por:
C0,5 = 0,1032857 log10 (Re) 0,137857
(2.88)
La ecuacin 2.87 es aplicable para los mismos nmeros de Reynolds y da como resultado una correlacin de
R = 0,914.
Una correlacin ms sosticada que recoge el
efecto de la rugosidad puede ser expresada como:

C0,5 Cmean
fe
= 0,5 1 + tanh
f
Cmean (1 Cmean )
(2.89)
C0,5 = 0,0593 + 0,07494 log10 (Re)
(2.90)
Donde:

ks
= 0,4726 1 + (3,6644 0,4729 log10 (Re)) 2,5915 + log10
DH
Siendo
DH
(2.91)
el dimetro hidrulico. Esta ltima correlacin da lugar a un coe-
ciente de correlacin de
R = 0,920.
En la gura 2.21 se muestra una comparacin
de los datos observados con la prediccin de la ecuacin 2.89.
60
Figura 2.21: Valores calculados con la ecuacin 2.89 frente a valores observados.
2.4. Fundamentos numricos

2.4.1. El Mtodo de los Volmenes Finitos
El Mtodo de los Volmenes Finitos (MVF) es la tcnica por la cual la formulacin integral de las leyes de conservacin se discretizan directamente en el espacio
fsico. La idea de partida fundamental es el concepto de discretizacin: reemplazar
una solucin analtica en derivadas parciales que proporciona el valor de forma continua en todos los puntos del espacio por una solucin numrica aproximada que da
el valor de la solucin nicamente en una serie de puntos discretos denidos por la
malla del dominio. Una extensa revisin del mtodo y su aplicacin al CFD puede
encontrarse en (Versteeg and Malalasekera, 2007).
El mtodo de volmenes nitos propone una forma de llevar a cabo esa discretizacin. En particular, establece que los valores discretos de la solucin quedarn
descritos por un conjunto de ecuaciones algebraicas que relacionan el valor de la
variable en un punto con el valor en los puntos vecinos. La forma en la que se transmite la informacin entre esos nodos requiere de algn tipo de aproximacin, que
en el caso de volmenes nitos se realiza mediante esquemas conservativos que evalan los ujos a travs de supercies de control (Fernndez Oro, 2012; Blazek, 2005).
61
La transformacin de las ecuaciones diferenciales en un conjunto de ecuaciones

algebraicas precisa de una discretizacin espacial. Esto se consigue generando una
malla que permite dividir el dominio de inters en una serie de celdas a las cuales
se les asocia el valor de la variable discreta.
La razn por la cual el Mtodo de Volmenes Finitos es el ms empleado para

desarrollar cdigos CFD reside en su generalidad, su simplicidad conceptual y su
facilidad para ser implementado en cualquier tipo de mallado, ya sea estructurado
o no estructurado.
Existen una serie de propiedades fundamentales en el Mtodo de Volmenes Finitos que merece la pena destacar. En primer lugar, hay que recordar que el mtodo
discretiza el dominio de un nmero nito de celdas o volmenes de control. Por
tanto, el mtodo se basa en valores discretos que estn promediados en la celda,
la variable numrica fundamental en toda aplicacin de CFD. Esto distingue claramente al Mtodo de Volmenes Finitos de otros mtodos como diferencias nitas o
elementos nitos, en los que la variable numrica fundamental es el valor local de la
funcin en los nodos de la malla.
Tras la generacin del mallado, el Mtodo de Volmenes Finitos asocia un volumen nito local, tambin denominado volumen de control, a cada punto de la malla,
y a continuacin aplica las leyes integrales de conservacin a cada volumen local.
Esta es otra gran diferencia con el mtodo de diferencias nitas, donde el espacio
discretizado se considera como un conjunto de puntos, mientras que en el MVF el
espacio discretizado est formado por un conjunto de pequeas celdas y cada una
de ellas est asociada a un nodo de la malla.
El Mtodo de Volmenes Finitos tiene ms ventajas con las mallas arbitrarias, en

las cuales un gran nmero de opciones quedan abiertas para poder denir distintos
volmenes de control en los que imponer las leyes de conservacin. La posibilidad de
modicar la forma y la localizacin de los volmenes nitos asociados a los nodos, as
como las leyes y la precisin en la evaluacin de los ujos a travs de las supercies
de control, proporciona una enorme exibilidad al MVF. Esto explica la generalidad
del mtodo, an cuando existan una serie de reglas que deben respetarse durante
estas operaciones.
Una de las ventajas del Mtodo de Volmenes Finitos es que garantiza una discretizacin conservativa. En las ecuaciones discretizadas es muy importante mantener
la conservacin global de las variables bsicas del ujo: masa, cantidad de movimiento y energa; por lo que ser preciso respetar una serie de condiciones en la forma
en que se lleva a cabo la discretizacin.
El Mtodo de Volmenes Finitos tiene, por tanto, una ventaja decisiva: la discretizacin conservativa se satisface automticamente con la discretizacin directa
62
de la forma integral de las ecuaciones de gobierno.
2.4.2. Grados de aproximacin a la resolucin 3D de las ecuaciones del ujo

Las aproximaciones para resolver las ecuaciones de gobierno en ujo turbulento
se pueden clasicar en dos grandes grupos: la resolucin directa de las ecuaciones
(DNS) o bien el empleo de algn modelo de turbulencia una vez realizado un promedio temporal (RANS) o un ltrado espacial de las ecuaciones (LES).
La simulacin numrica directa resuelve todas las uctuaciones temporales y espaciales. Esto implica que para nmeros de Reynolds elevados se exigen mallados
extraordinariamente densos y pasos temporales muy bajos para poder resolver todos los tamaos de vrtice del ujo. En realidad, el procedimiento de resolucin es
anlogo al que se emplea para ujo laminar, si bien debe afrontar la resolucin de
las uctuaciones a todas las escalas. Desafortunadamente, con los medios computacionales disponibles hoy da, solamente es posible resolver con xito este problema,
obteniendo resultados tiles, para geometras muy sencillas con unos costes computacionales generalmente inasumibles. Es preciso hacer hincapi en que la forma de
trabajar de la viscosidad en las pequeas escalas es a travs de largas longitudes
de gran desarrollo, con mucha curvatura, que se entrelazan a escalas espaciales muy
pequeas. Por lo tanto, para poder resolver este mecanismo numricamente, es necesaria una gran denicin espacial, lo que implica grandes mallas con nmeros de
celdas demasiado elevados. Si no se resolviera a esta escala, solo se estara evaluando
el efecto de la viscosidad en las grandes estructuras del ujo, dando lugar a zonas de
disipacin pequeas cercanas a los grandes ujos donde los gradientes de velocidad
sean importantes desde un punto de vista macroscpico. Esto nada tiene que ver con
la estructura fractal de la turbulencia, donde una lnea de corriente aparentemente
suave y tendida se curva dando lugar a una lnea de corriente mucho ms larga de lo
que pueda parecer desde un punto de vista macroscpico y por lo tanto produciendo
mayores supercies de disipacin y mayores prdidas de energa.
La alternativa ms habitual al DNS es la resolucin de las ecuaciones promediadas, las RANS. Por tanto, en esta aproximacin se resuelven nicamente los valores
medios de las variables uidodinmicas, pero a costa de introducir una serie de trminos adicionales, las tensiones de Reynolds, que habr que modelizar para poder
garantizar el cierre del sistema de ecuaciones. Lgicamente, estos modelos son generales y pueden provocar importantes errores en los clculos numricos, pues se
construyen sobre hiptesis que muchas veces se alejan de la realidad.
63
2.4.3. Modelado de la turbulencia

2.4.3.1. Hiptesis Boussinesq o de viscosidad turbulenta
En el desarrollo histrico de cualquier campo cientco en la investigacin, es
usual que haya una sucesin de modelos propuestos para describir los fenmenos
que se estudian. A menudo, tal y como sucede en el estudio de los ujos turbulentos,
los primeros modelos son simples, pero resultan ser decientes tanto en sentido fsico como en la precisin de las predicciones. Los modelos ms avanzados pueden ser
superiores en sentido fsico y en precisin, pero carecen de simplicidad. A pesar de
sus defectos, es importante contar con una estimacin de los primeros modelos, los
ms sencillos. Los modelos ms simples pueden proporcionar una referencia sobre
cual de los fenmenos de estudios es predominante, que puede ser estudiado ms
adelante por un modelo ms complejo.
En el caso de uidos newtonianos, la Mecnica de Fluidos clsica establece que
existe una relacin lineal entre el tensor de tensiones viscosas y el tensor de deformaciones en el seno del uido. Esta evidencia, encontradada por Navier y Poisson para
casi todos los lquidos y gases, permite generalizar la ley de Newton que relaciona los
esfuerzos cortantes con los gradientes de velocidad mediante la viscosidad molecular
del uido segn la ecuacin 2.4.
Con el objetivo de tratar de resolver el problema del cierre, tal y como ha sido
enunciado en el Cap. 2.2.7, Boussinesque propuso (Boussinesq, 1877) que deba existir alguna analoga entre la interaccin de las tensiones viscosas y de las tensiones de
Reynolds. Constatando que las tensiones turbulentas aumentan cuando se incrementan las componentes del tensor promedio de deformaciones, ide que las tensiones
de Reynolds deban estar ligadas con el tensor
Sij
por medio de un coeciente de
viscosidad articial o turbulenta:
ijR
vi0 vj0

t
vi
vj
+
xj xi

= 2t Sij
(2.92)
Por analoga, el transporte turbulento de calor, masa u otra variable escalar se dene
de manera similar, suponiendo que el transporte turbulento del escalar es proporcional a los gradientes del valor medio de la variable transportada a travs de una
difusividad turbulenta.
La implicacin fsica de esta hiptesis es que el efecto de las tensiones adicionales
de Reynolds se sustituye por la contribucin de una viscosidad turbulenta adicional,
t , que se suma a la viscosidad propia del uido. De esta forma, es como si el uido,
al estar en rgimen turbulento, tuviera de forma efectiva una mayor viscosidad. Este
articio se corresponde con la idea perfectamente contrastada de que la turbulencia
favorece la mezcla y la difusin del ujo.
64
En su concepcin original, Boussinesq adopt para la viscosidad turbulenta un

valor constante, eleccin que slo proporciona buenos resultados en casos muy concretos, relacionados con ujos libres (sin efectos de pared) como chorros axisimtricos
y chorros 2-D o capas de mezcla, tratando de adaptar la idea de la viscosidad articial a esos efectos de la capa lmite en los contornos slidos. Este modelo de longitud
de mezcla (y otros derivados de l) describen las tensiones por medio de una sencilla
relacin algebraica para
en funcin de la posicin, razn por la cual se denominan
modelos de cero ecuaciones.
En los modelos algebraicos, el siguiente paso natural es introducir algn tipo de

ecuacin diferencial adicional, por ejemplo para el transporte de la energa cintica
turbulenta. Es el caso del modelo de Spalart-Allmaras, principal exponente de los
modelos de una ecuacin, en el que se resuelve una viscosidad turbulenta modicada. Sin embargo, la opcin ms habitual es la de utilizar modelos de dos ecuaciones,
como el clsico modelo
,
un modelo ms sosticado que permite la descripcin
de las propiedades turbulentas bsicas: la energa cintica turbulenta

disipacin de energa turbulenta
.
y la tasa de
Se introducen dos ecuaciones de transporte para
cada una de estas variables, que se encuentran relacionadas con la viscosidad turbulenta por medio de su cociente adimensional. Este modelo y sus variantes siguen
siendo hoy da los modelos de partida bsicos y los ms empleados en simulaciones
de tipo ingenieril.
Todos los modelos anteriores se basan en la hiptesis de Boussinesq, la cual supone implcitamente que la viscosidad turbulenta es isotrpica, es decir, la relacin
entre las tensiones de Reynolds y el tensor promedio de deformacin es la misma
en todas las direcciones. Por esta razn, se engloban todos ellos bajo la denominacin de modelos lineales de viscosidad turbulenta (Eddy Viscosity Models, EVM).
Aunque vlida en un buen nmero de situaciones, esta hiptesis pierde validez en
situaciones de ujos tridimensionales, con mucha separacin y gradientes de presin
extremadamente adversos.
Cuando la hiptesis de Boussinesq no sea vlida, ser preciso denir ecuaciones de transporte para cada una de las seis tensiones de Reynolds. La cantidad de
momento uctuante de las partculas de uido tambin es transportado por el ujo
promedio, por lo que es perfectamente razonable plantear en este caso una ecuacin de transporte para cada una de ellas. Lgicamente, al necesitar una ecuacin
de transporte para cada tensin, el coste computacional de estos modelos de tensiones de Reynolds (RSM), o modelos de segundo orden, es notablemente superior
a los EVM. Esto implica que su utilizacin frente a los anteriores debe estar bien
justicada en trminos de precisin y tiempo de clculo. Los trminos de difusin,
disipacin y deformacin por la presin que aparecen en las ecuaciones de transporte
son desconocidos a priori, por lo que ser necesario hacer algn tipo de suposicin
sobre su comportamiento en el sistema.
65
2.4.3.2. La ecuacin de la energa cintica turbulenta

La energa cintica turbulenta supone una variable de referencia en la modelacin de la turbulencia. Su implicacin directa en las RANS la hace insustituible y
los buenos resultados de estos modelos la avalan. Es por ello importante comprender
como surge la ecuacin de transporte para esta variable y las simplicaciones que
sufre en los modelos de turbulencia, pues en tanto en cuanto seamos capaces de predecir el valor real de esta variable, mejor estaremos aproximando nuestra solucin.
Los modelos de turbulencia basados en la ecuacin de la energa se desarrollaron
con la idea de incorporar los efectos no locales y los efectos de la historia del ujo
en la hiptesis de viscosidad turbulenta. Prandtl (1945) insista en la importancia
de calcular una escala de velocidades caracterstica. De hecho, tom como base de
su escala de velocidades la energa cintica turbulenta por unidad de masa, que normalmente se denota como
1
= u0i u0i
2
As, en trminos de
(2.93)
y la disipacin de energa turbulenta
,
mediante anlisis di-
mensional se puede obtener:
2

(2.94)
Lo que supone una relacin entre la viscosidad y la energa cintica turbulenta. Sin
embargo, es necesario otra variable adicional, que en este caso ha sido
relacionando
.
Asimismo,
con el tensor de esfuerzos de Reynolds se obtiene:
iiR = vi0 vi0 = 2
(2.95)
De este modo, la traza del tensor de esfuerzos de Reynolds es proporcional a la

energa cintica turbulenta de las uctuaciones de velocidad por unidad de masa.
Esta relacin ser til ms adelante. La cantidad
debera ser referida estrictamente
como energa cintica especca turbulenta al ser una cantidad por unidad de masa.
Sin embargo, y como ya se habr notado, se suele llamar energa cintica turbulenta
a secas.
Si se dene el operador de Navier-Stokes
66
N (ui )
como:
N (ui ) =
ui
ui
1 p
+ uj
gi +
2 ui
t
xj
xi
(2.96)
La ecuacin de Navier-Stokes puede ser escrita entonces como:
N (ui ) = 0
(2.97)
Multiplicando el operador por la desviacin de las velocidades de la siguiente forma

y promediando se puede obtener un resultado muy interesante:
u0i N (uj ) + u0j N (ui ) = 0
(2.98)
Operando y teniendo en cuenta la simetra del tensor de tensiones de Reynolds,

la ecuacin resultante es tambin simtrica en
y en
j.
Por claridad, procedemos
trmino a trmino de la anterior ecuacin, teniendo en cuenta adems la ecuacin

u0i
u
= 0). Adems se ha empleado la notacin tensorial para
de continuidad ( i =
xi
xi
las derivadas por economa de espacio. As, para el trmino temporal se tiene:
u0i (uj ),t + u0j (ui ),t = u0i (uj + u0j ),t + u0j (ui + u0i ),t
= u0i uj,t + u0i u0j,t + uj ui,t + u0j u0i,t
= u0i u0j,t + u0j u0i,t
= (u0i u0j ),t
ijR
=
t
(2.99)
Desarrollando el trmino convectivo se obtiene:
u0i uk uj,k + u0j uk ui,k = u0i (uk + u0k )(uj + u0j ),k + u0j (uk + u0k )(ui + u0i ),k
= u0i uk u0j,k + u0i u0k (uj + u0j ),k + u0j uk u0i,k + u0j u0k (ui + u0i ),k
= uk (u0i u0j ),k + u0i u0k uj,k + u0j u0k ui,k + uk (u0i u0j ),k
= uk
ijR
uj
R ui
ijR
jk
+
(u0 u0 u0 )
xk
xk
xk
xk i j k
67
(2.100)
Para el trmino del gradiente de presiones se obtiene:
u0i p,j + u0j p,i = u0i (p + p0 ),j + u0j (p + p0 ),i

= u0i p0,j + u0j p0,i
= u0i
p0
p0
+ u0j
xj
xi
(2.101)
Finalmente, de los trminos viscosos se obtiene:
(u0i uj,kk + u0j ui,kk ) = u0i (uj + u0j ),kk + u0j (ui + u0i ),kk
= u0i u0j,kk + u0j u0i,kk
= (u0i u0j,k ),k + (u0j u0i,k ),k 2u0i,k u0j,k
= (u0i u0j ),kk 2u0i,k u0j,k
2 ijR
u0i u0j
=
2
xk xk
xk xk
(2.102)
Recogiendo los trminos obtenidos en los desarrollos anteriores se obtiene la ecuacin:
"
#
0
R
0 u0
0
0
0
ijR
u
ijR
u
u
u
u
p
p
j
i
j
j
ij
+uk
= ikR
R
+2 i
+ i
+
+
+ u0i u0j u0k

t
xk
xk jk xk
xk xk xj xi xk
xk
(2.103)
Lo que supone 6 nuevas ecuaciones, una por cada componente del tensor de tensiones de Reynolds. Sin embargo, se obtienen tambin 22 nuevas incgnitas:
u0i u0j u0k 10 inc

ognitas
u0i u0j
6 inc
ognitas
xk xk
u0j p0
u0i p0
+
6 inc
ognitas
xj
xi
2
68
Agrupando trminos, se puede escribir de forma ms compacta:
"
#
R
ijR
ijR
u
u
j
i
ij
R
+ uk
= ikR
jk
+ ij ij +
+ Cijk
t
xk
xk
xk
xk
xk
(2.104)
donde:
ij
p
=
u0i u0j
+
xj
xi
u0i u0j
xk xk
= u0i u0j u0k + p0 u0i jk + p0 u0j ik
ij = 2
Cijk
(2.105)
El desarrollo anterior ilustra perfectamente el problema del cierre pues segn se toman momentos de mayor orden, se generan nuevas incgnitas y ecuaciones. Despus
de todo, no es tan sorprendente, pues las operaciones realizadas son puramente matemticas y no introducen ninguna condicin fsica nueva, lo que no podra cerrar
el problema por s mismo. En esencia, el promediado de Reynolds (que da lugar al
problema del cierre) es una simplicacin brutal que da lugar a la prdida de gran
parte de la informacin de las ecuaciones de Navier-Stokes. La funcin de la modelacin de la turbulencia es dar un valor a los trminos desconocidos que se producen al
promediar, basndose en las propiedades del ujo y estableciendo nuevas ecuaciones
para satisfacer las nuevas incgnitas introducidas. Haciendo estas aproximaciones se
cierra el sistema y se puede resolver el problema.
Tomando la traza de ecuacin anterior y teniendo en cuenta la incompresibilidad
y ijR se obtiene la siguiente ecuacin de transporte
para la energa cintica turbulenta:
del ujo y la relacin entre

ui
1 0 0 0
1 0 0
+ uj
= ij
+
ui ui uj p uj
t
xj
xj
xj
xj
2
La cantidad
(2.106)
es la disipacin por unidad de masa y se dene como:
=
u0i u0i
xk xk
69
(2.107)
Los trminos que aparecen en la anterior ecuacin de transporte para la energa

turbulenta representan procesos fsicos que ocurren segn la turbulencia se mueve
en el seno del ujo. La suma de los dos trminos de la izquierda, el trmino temporal
y el convectivo, forman la derivada material de
que viene a ser la variacin de
siguiendo una partcula uida. El primer trmino del lado derecho es conocido como
la produccin de energa turbulenta y representa el ratio de transferencia de energa
turbulenta como consecuencia del ujo medio. La disipacin es el proceso por el cual
la energa turbulenta se convierte en calor. La difusin molecular que depende de la
viscosidad representa el ratio de transferencia de energa turbulenta por difusin. El
trmino de la triple correlacin de velocidades se conoce como transporte turbulento
y maniesta el ratio por el cual la turbulencia es transportada por las desviaciones
de las velocidades. El ltimo trmino se conoce como la difusin de la presin.
De este modo, la ecuacin estrictamente exacta de la energa cintica turbulenta ha sido formulada, pero como ya se ha mencionado, presenta nuevas incgnitas.
El modo en el que los diferentes modelos de turbulencia de una o dos ecuaciones
reformulen la ecuacin de transporte para la cantidad
cerrando nuevamente el
problema, har ms vlida o menos la aproximacin.
2.4.3.3. Modelos de 0 ecuaciones

Los modelos de turbulencia de 0 ecuaciones son modelos mediante los cuales la
viscosidad turbulenta se relaciona con las variables medias del ujo mediante un
ecuacin algebraica.
La hiptesis fundamental de estos modelos es que la produccin de energa turbulenta y el ratio de disipacin son aproximadamente iguales. Esta hiptesis es evidentemente restrictiva y aade adems que no se modela el transporte de la energa
cintica turbulenta, lo que no dota de cierta historia al ujo. Todo esto es contrario
a la fsica de la mayora de ujos, sin embargo supone una forma fcil de cerrar las
RANS.
Muchos modelos suelen incluir una formulacin para la regin interna y otra para
la externa de la capa lmite para representar la longitud de mezcla.
2.4.3.3.1. Modelo de van Driest

Un modelo usado comnmente utiliza una funcin exponencial (funcin de amortiguamiento de van Driest) para la regin interna, mientras que la regin externa es
proporcional al espesor de la capa lmite. Matemticamente se extresa como:
70
li = (1 ey
+ /A+
)y
(2.108)
lo = Co
donde
(2.109)
es la constante de von Karman ( 0,41), y A+ es un parmetro que depende
del gradiente de presiones. Para un valor nulo del gradiente de presiones, tiene un
valor de 26. La constante
0,08 0,09,
mientras que
Co en la ecuacin 2.109 normalmente

es el espesor de la capa lmite.
toma un valor entre
2.4.3.3.2. Modelo Cebeci/Smith

Otra formulacin muy comn para la viscosidad turbulenta de la regin externa
es la del modelo de Cebeci/Smith expresada como:
nut = ue
(2.110)
toma un valor de 0.0168 para ujos donde el nmero de Reynolds del espesor
cantidad de movimiento Re es mayor de 5000 y el espesor de desplazamiento
donde
de la
se dene como:
=
0

u
dy
1
ue
P ara f lujo incompresible.
(2.111)
P ara f lujo compresible.
(2.112)
y:
=
0

u
1
dy
e ue
Recordar que el nmero de Reynolds basado en el momento del espesor se dene

como:
Re =
e ue
e
71
(2.113)
donde el espesor de cantidad de movimiento se dene como:
Z
=
0
u
ue
u
1
ue

dy
(2.114)
y:
Z
=
0
u
e ue

u
1
dy
ue
(2.115)
El modelo algebraico descrito requiere informacin del espesor y las propiedades

del ujo de la capa lmite. Cuando se resuelven las ecuaciones de Navier-Stokes, es
difcil determinar el espesor y las propiedades requeridas. Esto sucede especialmente
cuando existe separacin del ujo.
2.4.3.3.3. Modelo de Baldwin y Lomax

Un modelo de turbulencia que no es escrito en trminos de las propiedades de la
capa lmite es el introducido por Baldwin y Lomax. La regin interna es aproximada
mediante:
t = l2 | |
donde
es la vorticidad que en 2 dimensiones tiene una sola componente:
=
y
(2.116)
v u
x y
se calcula mediante la ecuacin 2.108. La coordenada adimensional
(2.117)
y+
se puede
escribir como:
y
y + = u
=
w
o:
72
| w | w
y
w w
(2.118)
y+ =
donde el subndice
p
y
w | w |
w
(2.119)
hace referencia a las cantidades de la pared. La regin exterior
es aproximada mediante:
t = CCP Fwake FKleb
donde
toma un valor de 0.0168 (como en el modelo de Cebeci/Smith) y:
"
Fwake
Un valor tpico de
Cwake
(V )2
= min ymax Gmax , Cwake ymax
Gmax
es 1.0. El parmetro

Gmax = max
Donde
(2.120)
es la constante de von Karman y
Gmax
#
(2.121)
se dene como:

l
||
(2.122)
la longitud de mezcla que queda deter-
minada por la funcin de amortiguamiento de la ecuacin 2.108 del modelo de Van

Driest. La diferencia entre los valores absolutos de las velocidades mxima y mnima
dentro de la regin viscosa se denota
V .
Para ujos afectados por una pared, la
velocidad mnima ocurre en la supercie y vale cero, entonces:
1/2
V = (u2 + v 2 )max
y para ujos de corte,

la velocidad en
(2.123)
se dene como la diferencia entre la velocidad mxima y
ymax :
1/2
V = (u2 + v 2 )max
(u2 + v 2 )1/2
ymax
FKleb
(2.124)
es el factor de intermitencia denido como:
"
FKleb = 1 + 5,5
73
CKleb y
ymax
6 #1
(2.125)
ymax
es el valor de
para el cual
Gmax
ocurre.
Los valores tpicos para la constante de Klebano
CKleb
CCP
son respectiva-
mente 0.3 y 1.6 para gradientes de presin suaves o nulos. Sin embargo, para incluir
la inuencia del gradiente de presin, se pueden emplear las siguientes expresiones:
CKleb =
0,01312
2
3 0,1724
(2.126)
Donde:
=
y
ymax V
u x
es la velocidad de cortadura. El gradiente de velocidad en
(2.127)
se calcula fuera de
la regin viscosa. Una vez se evala la constante de Klebano, se puede determinar
CCP
mediante:
cCP =
3 4CKleb
3
2CKleb (2 3CKleb + CKleb
)
(2.128)
Finalmente, la viscosidad turbulenta a lo largo de la capa lmite se obtiene como:
t = min(ti , to )
Lo que signica que
a
to ,
donde
se ja a
pasa a valer
ti
desde la pared hasta un lugar en el que
(2.129)
ti
supera
to .
2.4.3.4. Modelos de 1 ecuacin

Los modelos de una ecuacin utilizan una ecuacin en derivadas parciales para
deducir la escala de velocidades turbulenta, que normalmente es la energa cintica
turbulenta
desarrollada en profundidad en el apartado 2.4.3.2, y una ecuacin
algebraica para la escala espacial turbulenta.

La viscosidad turbulenta se puede escribir por lo tanto como:
t = 1/2 l
74
(2.130)
2.4.3.4.1. Modelo de Baldwin-Barth

El modelo de turbulencia de una ecuacin de Baldwin-Barth se obtiene a partir
del
estndar de dos ecuaciones. Este modelo se desarrolla ms adelante en el
apartado 2.4.3.5.1.
El procedimiento para obtener el modelo de una ecuacin de Baldwin-Barth es com-
binar ambas ecuaciones del modelo
para dar lugar a una sola ecuacin en
trminos del nmero de Reynolds turbulento
ReT =
ReT
denido como:
2

(2.131)
Las variables empleadas para denir este nmero de Reynolds son

ntica turbulenta,
la viscosidad cinemtica y
la energa ci-
la tasa de disipacin turbulenta.
Con objeto de tener en cuenta la regin cercana a la pared, el nmero de Reynolds

turbulento se divide en dos partes de modo que
funcin de amortiguamiento siendo
ReT ' ReT
ReT = ReT f (ReT ),
f es
ReT .
donde
para grandes nmeros de
La ecuacin de transporte para el nmero de Reynolds turbulento
ReT
una
tiene la
forma:

t
ReT + ~v ReT =
+
2 ReT
t

q

1
(t ) ReT + (C2 f2 C1 ) ReT P

(2.132)
La produccin de turbulencia se expresa de la forma:

P = t
ui uj
+
xj
xi
La viscosidad cinemtica turbulenta
ui
2
t
xj
3
D1
D2
2
= t S 2
(2.133)
queda determinada por:
t = c (ReT )D1 D2
donde
uk
xk
(2.134)
son funciones de amortiguamiento que extienden la validez del mo-
75
delo a las regiones de pared y vienen dados por:
D1 = 1 ey
D2 = 1 e
+ /A+
y + /A+
2
(2.135)
(2.136)
Las funciones y constantes de la ecuacin 2.132 son las siguientes:
= (C2 C1 ) c /2

t = t

hp
C1
C1
1
+
f2 (y ) =
+ 1
+ D1 D2
D1 D2 +
C2
C2
y +

y+
1 y+ /A+
1 y+ /A+2
+
e
D2 + + e
D1
A2
D1 D2 A+
= 0,41,
C1 = 1,2,
C2 = 2,0
+
C = 0,09,
A = 26,
A+
2 = 10
(2.137)
2.4.3.4.2. Modelo de Spalart-Allmaras

El modelo de Spalart-Allmaras resuelve una ecuacin de transporte para una variable
que est relacionada con la viscosidad turbulenta. La ecuacin de transporte
empleada es la siguiente:

2
Cb1
d
1
2
=
(( + )) + Cb2 () +cb1 S(1ft2 ) Cw1 fw 2 ft2
+ft1 (q)2
dt
d
(2.138)
donde la viscosidad turbulenta viene dada por:
t = f1
Y:
76
(2.139)
3
3
3 + C1
f1 =
(2.140)
(2.141)
Varias funciones y constantes de la ecuacin 2.138 se denen como:
S=S+
donde
es la distancia a la pared,
d2
2 d2
(2.142)
es la constante de von Karman,
es la magni-
tud de la vorticidad:

v u

S=
x x
(2.143)
Y:
f2 = 1
La funcin
fw
1 + f1
(2.144)
es:
fw (r) = g
6
1 + Cw3
6
g 6 + Cw3
1/6
!
(2.145)
Donde:
g = r + Cw2 (r6 r)
(2.146)
Y:
r=
Grandes valores de
S2 d2
se truncan a 10. La funcin
77
(2.147)
ft2
viene dada por:
ft2 = Ct3 eCt4

y la funcin
ft1
(2.148)
es:
2 2
2 2
ft1 = Ct1 gt e[Ct2 (t /q) (d +gt dt )]
donde
dt
es la distancia de cualquier punto a la supercie,
supercie,
(2.149)
es la vorticidad en esta
la diferencia de velocidades entre el punto y la supercie y
mnimo de 1.0 y
q/t x,
donde
gt
es el
es el espaciamiento de la malla a lo largo de
la pared.
Las constantes empleadas son:
=
Cw1 =
C1
Cb1
+ (1 + Cb2 )/
2
= 7,1
Ct1 = 1,0
2
3
Cb1 = 0,1355
Cb2 = 0,622
Cw2 = 0,3
Cw3 = 2,0
= 0,41
Ct2 = 2,0
Ct3 = 1,1
Ct4 = 2,0
(2.150)
2.4.3.5. Modelos de 2 ecuaciones

Los modelos de una ecuacin se basan en la resolucin de la ecuacin de la energa turbulenta o ms bien de su aproximacin. Los de dos ecuaciones trabajan con
la misma ecuacin de la energa turbulenta y con otra adicional. Los modelos de dos
ecuaciones han sido la base de la investigacin en turbulencia durante las ltimas 3
dcadas. Se consideran modelos completos (Wilcox et al., 1998; Huang et al., 1997;
Patankar, 1980), pues sirven para estudiar ujos turbulentos sin conocer a priori la
estructura de la turbulencia.
El punto de partida para seguramente todos los modelos de dos ecuaciones es la
hiptesis de Boussinesq (Cap. 2.4.3.1), que relaciona una variable escalar conocida
como viscosidad turbulenta
T , mediante anlisis dimensional, con y otra variable
del ujo.
Por denicin, la disipacin turbulenta es no negativa y tiene por dimensiones
[L2 /T 3 ]. Una escala longitudinal turbulenta podra deducirse a partir de y como:
78
LT
3/2

(2.151)
Del mismo modo, la escala temporal turbulenta sera:
TT
(2.152)
La viscosidad turbulenta ser proporcional al producto de la velocidad turbulenta y

la escala espacial turbulenta. As se obtiene:
3/2
(2.153)
Siendo estrictos, esa proporcionalidad puede expresarse como:
T = C
2

(2.154)
Kolmogorov (1942) ya expresaba la necesidad de formular dos ecuaciones, una para
y otra para otra variable .

Si formulramos y desarrollramos:
u0i
[N (ui )]
xj xj
Se obtendra una ecuacin matemticamente exacta para
(2.155)
procediendo del mismo
modo que con la ecuacin de la energa cintica turbulenta en el apartado anterior:

ui

2 ui
+ uj
= 2 u0i,k u0j,k + u0k,i u0k,j
2u0k u0i,j
t
xj
xj
xk xj
2u0i,k u0i,m u0k,m 2 2 ui,km u0i,km

0 0
0 0
0
+
uj ui,m ui,m 2 pm uj,m

xj
xj
3 La
(2.156)
eleccin de esta otra variable determinar el modelo de turbulencia de dos ecuaciones que
se emplea (p.e.: , ,...
79
Esta ecuacin es todava ms compleja que la de la energa turbulenta y presenta

nuevas incgnitas. Por ello, es preciso aproximar estas ecuaciones teniendo en cuenta
los procesos fsicos reales maniestos en los ujos.
De este modo, se pretende en los modelos de dos ecuaciones, modelar la fsica

de la turbulencia ms que las ecuaciones diferenciales exactas que gobiernan el comportamiento de
y de
,
dado que estas siguen siendo inabordables. Es importante
tener todo esto en cuenta para comprender las limitaciones ms que importantes de
estos modelos y lo que se pretende conseguir mediante ellos.
2.4.3.5.1. Modelo
El modelo de turbulencia
es el modelo de dos ecuaciones ms empleado
tradicionalmente. Se basa en la resolucin de una ecuacin para la energa cintica

turbulenta
y otra para el ratio de disipacin de energa turbulenta
.
Finalmente
estas cantidades se emplean para calcular la viscosidad turbulenta segn la ecuacin

2.154.
La ecuacin empleada para aproximar la energa cintica turbulenta es la siguiente:
ui
+ uj
= ijR
+
t
xj
xj
xj

T
+
xj
(2.157)
Y la ecuacin para la aproximacin de la disipacin de energa turbulenta es:

R ui
2
+ uj
= C1 ij
C2
t
xj
xj
xj

T
+
xj
(2.158)
Se observa claramente la analoga entre las ecuaciones empleadas para
y las
matemticamente exactas deducidas para estos valores en los anteriores apartados.

Aunque no se resuelven las ecuaciones exactas debido al nuevo problema del cierre
que surge en su obtencin, se modela el comportamiento fsico de
y de
y se
calibran los coecientes introducidos mediante modelacin fsica.
As se emplean los siguientes valores de los coecientes para el modelo

estndar:
80
C1
C2
C
=
=
=
=
=
1,44
1,92
0,09
1,0
1,3
(2.159)
2.4.3.5.2. Modelo RNG

La moderna versin RNG del clsico modelo
,
est basada en la teora ma-
temtica de Renormalizacin de Grupos (Yakhot et al., 1992; Yakhot and Orszag,

1986). Esta teora pretende recoger las interacciones entre distintos grados de libertad en sistemas, tanto lineales como no lineales. Este modelo de turbulencia tiene
ciertas ventajas sobre el anterior, como por ejemplo que es ms preciso para ujos
rpidamente variados, con curvaturas mayores de las lneas de corriente y para ujos
de menor nmero de Reynolds.
Las ecuaciones empleadas para
sin embargo, se modica el coeciente
son
C2 :
las mismas que en el anterior modelo ,
C 3 (1 /0 )
C2 = C2 +
1 + 3
(2.160)
Donde:
p
2Sij Sji

1 ui uj
=
+
2 xj
xi
=
Sij
(2.161)
Los valores de los coecientes toman los siguientes valores:
4 Aunque
en este caso las ecuaciones son derivadas a partir de las ecuaciones de Navier-Stokes
mediante el uso de la teora RNG(Pope, 2000).
81
C1
C2
C
=
=
=
=
=
=
=
1,42
1,68
0,085
0,72
0,72
0,012
4,38
(2.162)
2.4.4. Modelado de la lmina libre

El tratamiento de la lmina libre requiere una especial atencin, pues supone
desde el punto de vista numrico una condicin de contorno mvil. El tratamiento
de esta asegura un correcto campo de presiones y velocidades en los ujos en los
que la lmina libre juegue un papel importante. La lmina libre se considera una
supercie de discontinuidad de las variables, lo que da ciertos problemas a la hora
de evaluar los gradientes de las variables.
Son tres los problemas que surgen del tratamiento numrico de la lmina libre:
su representacin discreta, su evolucin en el tiempo y la forma en la que las condiciones de contorno se imponen.
Hay varias formas de llevar a cabo este tratamiento, tanto desde un punto de
vista euleriano como lagrangiano. La representacin lagrangiana del uido es conceptualmente ms sencilla ya que cada zona del mallado que subdivide el uido
permanece identicada con el mismo elemento uido durante todo el tiempo. Las
fuerzas de volumen y supercie son fciles de denir, por lo que es fcil de calcular
la respuesta dinmica de los elementos. En una representacin euleriana el mallado
permanece jo y la identidad de cada elemento uido no se mantiene. Sin embargo,
cada celda de la malla se puede interpretar como un elemento uido en cuyas caras
se calculan las fuerzas de supercie de una forma anloga al enfoque lagrangiano.
Los dos mtodos dieren en la manera en que los elementos uidos se mueven de una
posicin a otra despus de que se calculen sus velocidades. En el caso lagrangiano,
la malla simplemente se mueve con las velocidades calculadas, mientras que en la
euleriana es necesario calcular el ujo del uido entre las celdas de la malla. Este
ujo, o clculo de la conveccin, requiere un promediado de las propiedades del ujo
de todos los elementos uidos que se encuentran en una misma celda despus de
un periodo de tiempo. El promediado de la conveccin da lugar a un suavizado de
la solucin, en particular en las zonas de discontinuidad como sucede en la lmina
libre. El nico modo de compensar esta falta de precisin en los contornos mviles
es introducir un tratamiento especial que reconozca la discontinuidad y evite el promedio a travs de esta.
82
Una tcnica muy ecaz que incorpora las ventajas de los mtodos lagrangianos
para los esquemas eulerianos es el Volume Of Fluid (VOF) (Hirt and Nichols, 1981;
Nichols et al., 1980). Un buen contraste numrico experimental se pueden encontrar
en (Oertel and Bung, 2012), algunos resultados se muestran en las guras 2.23 y 2.24.
En este mtodo, se dene una funcin
en cada punto cuyo valor vale 0 para
puntos sin uido y 1 para puntos con uido. El promedio de estos valores para una
celda da la fraccin de uido que la llena. De este modo, las celdas con valores distintos de 0 o 1 son celdas que contienen una supercie libre.
Figura 2.22: Representacin tpica de una lmina libre mediante el VOF. En el centro
de las celdas se muestra el valor de
F.
La direccin normal al contorno se obtiene evaluando la direccin en la que vara

ms rpidamente el valor de
F . Dado que F
es una funcin escalonada, sus derivadas
deben ser calculadas de una manera especial. Una vez el valor de
y sus derivadas
son conocidas, se puede trazar una lnea en cada celda que aproxima la interfaz. La
localizacin de este contorno puede utilizarse entonces para imponer las condiciones
de contorno.
La dependencia temporal de
viene dada por la ecuacin:
F
F
F
F
+u
+v
+w
=0
t
x
y
z
83
(2.163)
Figura 2.23: Rotura de presa. Resultados del modelo fsico medidos con un PIV.
1 0,5
0,5
Velocidad adimensionalizada U = u/(4gd0 ) , con d0 = 0,40m y T = t(gd0 )
=
a) 0.50, b) 0.99, c) 1.98, e) 2.48 y f ) 2.97. Adaptado de (Oertel and Bung, 2012).
84
Figura 2.24: Rotura de presa. Resultados del modelo numrico empleando un m0,5
todo VOF. Velocidad adimensionalizada U = u/(4gd0 ) , con d0 = 0,40m y
0,5
T = t(gd1
= a) 0.50, b) 0.99, c) 1.98, e) 2.48 y f ) 2.97. Adaptado de (Oer0 )
tel and Bung, 2012).
85
Esta ecuacin de transporte comprende que
se mueve con el uido sin afectar
primeramente al campo de velocidades y de presiones. De hecho, el campo de velocidades se toma como conocido al inicio de cada iteracin en la evolucin temporal
de la lmina libre. Evidentemente, la evolucin de la lmina libre s que afectar al
campo de velocidades y de presiones, pero esto sucede en la siguiente iteracin.
En un esquema lagrangiano, todo esto se reduce a que
elemento uido. En el caso euleriano
es constante en cada
sirve nicamente de marcador. El enfoque
es por lo tanto el mismo, salvando las diferencias evidentes entre cada mtodo.
En resumen, este mtodo permite calcular la lmina libre de una forma able
con poco coste computacional. Dado que este mtodo marca la evolucin de zonas
de uido y no de supercies, cualquier problema relacionado con la interseccin
de supercies queda completamente resuelto. De este modo, se puede resolver un
problema que es evidentemente bifsico (agua-aire), resolviendo exclusivamente las
ecuaciones de la masa y del momento de un solo uido. As, este mtodo ser mucho menos caro computacionalmente que resolver las ecuaciones del ujo en ambos
dominios uidos pero presenta la desventaja de que no permite resolver el efecto del
aire sobre la lmina, ms all de especicar una condicin tipo presin en las celdas
de la interfaz. Sin embargo, y a pesar de este problema evidente, el mtodo funciona
muy bien para la mayora de ujos con lmina libre, si analizamos detenidamente el
efecto del arrastre del aire sobre el agua, nos damos cuenta de que la gran diferencia
de densidades permite obtener buenas aproximaciones despreciando este efecto del
aire sobre el agua.
2.4.5. Modelado de la reaireacin turbulenta

Dada la importancia que tiene el efecto del aire sobre el agua en los ujos extremadamente turbulentos con lmina libre, puede ser de inters evaluar las tasas de
entrada de gas en el lquido. Si se modela la entrada de aire a travs de la lmina
libre, debe tenerse en cuenta el efecto del cambio de densidad de cada celda debido a
la mezcla de dos uidos de diferente densidad. Se debe por lo tanto tener en cuenta
los ujos de densidad y el efecto de emulsionamiento de la lmina libre.
Para evaluar la tasa de entrada de gas a travs de la supercie del lquido, se
puede realizar un balance de fuerzas estabilizadoras y desestabilizadoras: por un
lado gravedad y tensin supercial y, por otro lado, la energa cintica turbulenta
en la lmina de libre. Dado que es necesario conocer el valor de la energa cintica
turbulenta, es necesario emplear un modelo de turbulencia con ecuacin de transporte para
La idea bsica es que los torbellinos turbulentos en la supercie levantan pequeas masas de agua que pueden atrapar aire y arrastrarlo hacia el interior del
86
caracteriza la turbulencia mediante

y la disipacin de energa turbulenta . Una medida
volumen uido. El modelo de turbulencia

la energa cintica turbulenta
caracterstica de los torbellinos turbulentos se puede obtener por anlisis dimensional segn la ecuacin 2.151 y se puede hacer exacta mediante:
r
LT = C
3 3/2
2
(2.164)
Con esta escala se caracterizan la altura de las perturbaciones de la lmina debido

al efecto de la turbulencia. Se puede denir por lo tanto una presin asociada a las
fuerzas estabilizadoras como:
Pd = gn LT +
Donde
es la densidad del lquido,
K
LT
(2.165)
es el coeciente de tensin supercial y
gn
es
la componente de la gravedad normal a la supercie libre.

Asimismo, se puede denir una medida de la presin de las fuerzas perturbadoras:
Pt =
(2.166)
Para que exista entrada de aire debe suceder que el valor de
Pt
sea mayor que
Pd .
De otro modo, las fuerzas estabilizadoras predominan y la lmina no se perturba lo

suciente para atrapar el gas. El volumen de aire entrado por unidad de tiempo,
debe ser proporcional al rea de la celda
As ,
V ,
y la altura de la perturbacin sobre el
nivel medio de la lmina. Se puede escribir por lo tanto:
s
V = Cair AS
Donde
Cair
Pt Pd
(2.167)
es un coeciente de proporcionalidad. El valor de este coeciente debe
ser menor que la unidad, pues solo una porcin del volumen perturbado es ocupado
por aire.
La experiencia demuestra que este modelo es apropiado para distintos tipos de aplicaciones. Sin embargo, hay una situacin donde el modelo debe ser modicado:
cuando un chorro intercepta la lmina de agua, se necesitara una resolucin exageradamente pequea tanto temporal como espacialmente para capturar la energa
turbulenta que realmente sucede. Para estos casos, se aumenta un 10 % la energa
87
cintica turbulenta asociada a la velocidad normal a la supercie. Adems, se modica la escala longitudinal turbulenta para estos casos.
Es necesario tener en cuenta que cuando se modela la entrada de aire, se est
empleando lo que se conoce como un modelo de submalla, es decir, un modelo que
recoge los efectos que tendra el clculo con una malla mejor. Ha de tenerse en cuenta que los procesos involucrados en el atrapamiento de burbujas en la lmina tienen
lugar en una escala espacial muy pequea, muy por debajo de lo que nos permite
calcular nuestra potencia computacional. Es decir, se debera resolver a una escala
casi innitamente pequea la interfaz bifsica para poder recoger las salpicaduras y
los atrapamientos sin necesidad de modelo de reaireacin.
88
Calibracin del modelo de reaireacin turbulenta

3.1. Introduccin y objetivos
En las grandes estructuras hidrulicas es comn encontrar ujos supercrticos

con grandes cantidades de aire atrapadas. Las burbujas de aire son transportadas
aguas abajo con el ujo medio, cambiando las propiedades macroscpicas del ujo:
aumentando el calado, reduciendo la densidad y cambiando las propiedades turbulentas. La aireacin del ujo induce cierto emulsionamiento (Chanson, 1996; Falvey,
1980) que cambia signicativamente la distribucin de cantidad de movimiento del
ujo de entrada a la obra de disipacin de energa. Afortunadamente, esta aireacin protege la obra de posibles daos por cavitacin lo que convierte la aireacin
en esencial para una correcta operacin de la estructura hidrulica (Wood, 1991).
La aireacin natural del ujo es la caracterstica ms importante de los ujos supercrticos (Novak et al., 2007), aunque tambin pueden emplearse aireadores para
emulsionar articialmente el ujo en los aliviaderos (Pster and Hager, 2010a,b).
Estos ujos supercrticos emulsionados son difciles de reproducir en laboratorio
y se encuentran normalmente sujetos a importantes efectos de escala (Heller, 2011;
Chanson, 2009, 2013). Se puede encontrar ejemplos de ujos preaireados articialmente en laboratorio mediante un jetbox en (Kramer et al., 2006; Schmocker et al.,
2008; Valero Huerta et al., 2013; Valero et al., 2014). Las tcnicas de medicin son
adems complejas y costosas. De este modo, la combinacin de herramientas CFD
se presentan como una herramienta muy til con posibilidad de reproducir las variables hidrodinmicas, tanto primitivas como derivadas, de los ujos multifase. Sin
embargo, la falta de validacin y comprobacin con datos libres de efectos de escala
es todava un problema por resolver (Chanson, 2013; Chanson and Lubin, 2010; Ma
et al., 2011b,a).
En esta primera parte del trabajo de investigacin, se ha llevado a cabo la calibracin del modelo de reaireacin turbulenta tal y como se ha presentado en el Cap.
1 Computational Fluid Dynamics.
Captulo 3. Calibracin del modelo de reaireacin turbulenta Daniel Valero Huerta
2.4.5. El objetivo principal es obtener unos valores de los parmetros que permitan
al modelo ser empleado en el diseo de grandes obras hidrulicas. Para ello, se han
empleado datos de presas a escala 1:1, libres de efectos de escala.
La alta sensibilidad del modelo a los parmetros de entrada aportados por otros
modelos como son la energa cintica turbulenta
(modelo de turbulencia, Cap.
2.4.3) o la supercie libre (modelo de lmina libre, Cap. 2.4.4 ) son tambin analizados.
De este modo, el principal objetivo de este captulo es obtener unos valores de
los parmetros del modelo de reaireacin turbulenta que permitan posteriormente
simular problemas relacionados con la aireacin del ujo en grandes presas.
3.2. Descripcin general

En el proceso de calibracin del modelo de reaireacin turbulenta se han llevado
a cabo ms de 200 simulaciones variando los parmetros del modelo y los tamaos
de celda del modelo. Con esto se pretende realizar un anlisis de sensibilidad no
solo al valor de los parmetros (K ,
Cair ),
si no tambin al grado de convergencia
numrica de las variables hidrodinmicas que, como se ver ms adelante, juega un
papel especialmente importante en el uso de este modelo .
Figura 3.1: Crecimiento de la capa lmite e interseccin con la lmina libre segn el
mtodo de (Wood et al., 1983; Wood, 1991).
2 Ms
informacin en el Cap. 3.3.2.

92
Para la calibracin del modelo se ha elegido un problema lo ms cannico posible

dentro de la ingeniera de presas. La geometra modelada consiste en un aliviadero
de 20 metros de altura con un perl tipo USBR (Peterka, 1958) y una pendiente
de 53.1 grados (0.75H:1V) . La altura de energa en el embalse es de 1 metro sobre

la cresta del aliviadero, lo que se corresponde con un caudal especco desaguado
2
de 2.18 m /s. Para esta geometra el punto crtico de aireacin tiene lugar a 12.1
metros (medido en la proyeccin horizontal) desde el intrados de la presa; la concentracin de equilibrio esperada es de 0.65 acorde a (Chanson, 1993). El punto crtico
de inicio de aireacin se produce cuando el ujo es sucientemente turbulento como
para vencer las fuerzas estabilizadoras (gravedad y tensin supercial) y as atrapar pequeos volmenes de aire. Esto generalmente sucede cuando la capa lmite
intersecta la lmina libre. Para obtener este valor se ha empleado la ecuacin 2.74
para la capa lmite y la ecuacin 2.78 desarrolladas en el Cap. 2.3.1.4. Este mtodo
de obtencin del punto crtico de aireacin es el conocido mtodo de Wood (Wood
et al., 1983; Wood, 1991). La rugosidad empleada es de 1 mm, tal y como se sugiere
en trabajos previos (Gutirrez Serret, 1995). Se puede encontrar un esquema del
resultado en la gura 3.1.
Para el clculo del valor objetivo de la concentracin de equilibrio se han empleado los datos experimentales recogidos por (Chanson, 1993). Aunque los datos no
estn tomados no son originales de este trabajo del Dr. Chanson, hace un esfuerzo
de recopilacin, sntesis y anlisis importante . De aqu se extrae un valor esperable

de 0.65 para la concentracin de aire en agua. En la tabla 3.1 se recogen los valores
de concentraciones de aire de equilibrio en el ujo completamente desarrollado.
Pendiente
7.5
Cmean
0.161
15.0 22.5
0.241
0.310
30
37.5
45
60
75
0.410
0.569
0.622
0.680
0.721
Tabla 3.1: Concentraciones de equilibrio en el ujo completamente desarrollado segn los datos de (Straub and Anderson, 1958).
3.3. Opciones numricas

3.3.1. Generalidades
El modelo numrico empleado es FLOW-3D, que resuelve numricamente las
ecuaciones de Navier-Stokes. Para ello, discretiza las ecuaciones por el Mtodo de
los Volmenes Finitos y resuelve el campo de velocidades en las caras de las celdas,
almacenando las presiones en el centro. Una explicacin ms detallada del mtodo
3 Dicha
pendiente corresponde a uno de los valores ms elevados que se suelen emplear.

valores fueron obtenidos principalmente en el estudio de (Straub and Anderson, 1958),
aunque han quedado ampliamente vericados posteriormente.
4 Los
93
se encuentra en el Cap. 2.4.1. Las presiones se acoplan con las velocidades mediante
la resolucin de una ecuacin de Poisson parecida a la desarrollada en 2.2.2. Esto
da lugar a un sistema de ecuaciones lineales que se resuelve de forma iterativa. La
discretizacin del trmino advectivo se realiza empleando un esquema de segundo orden con conservacin de los gradientes (Hirsch, 2007). Esto asegura un rpido grado
de convergencia en el clculo de las velocidades sin perder precisin. Las tensiones
tangenciales, sin embargo, se han calculado implcitamente. El mtodo empleado
para la resolucin del sistema de ecuaciones lineales es el mtodo de proyeccin sobre el subespacio de Krylov; este tipo de mtodos presenta una convergencia rpida
con un muy buen grado de precisin. Se puede encontrar una discusin sobre los
mtodos de resolucin de sistemas de ecauciones lineales (tanto iterativos como de
proyeccin) en el libro (Saad, 2003).
El ncleo central de la modelacin numrica del ujo multifase es el constituido
por el tringulo: modelo de reaireacin turbulenta, modelo de turbulencia y modelo
de lmina libre (Cap. 2.4.5), 2.4.3 y 2.4.4. El objetivo de la calibracin es obtener
con un grado de precisin aceptable las concentraciones
que tienen lugar en el
aliviadero, para ello es necesario calcular la energa cintica turbulenta
y su evo-
lucin dentro del lete uido, nalmente es necesario determinar correctamente la

localizacin de la lmina libre. Para esta investigacin, los modelos empleados han
sido para la turbulencia el RNG
(descrito con mayor detalle en el Cap. 2.4.3.5.2)
especialmente adecuado para aplicaciones hidrulicas; el VOF segn se ha descrito

en el Cap. 2.4.4; y el modelo de reaireacin turbulenta (tal y como se dene en el
Cap. 2.4.5) en el que los parmetros
Cair
han tomado valores comprendidos
en los rangos 0 - 1 y 0 - 5 N/m correspondientemente.
3.3.2. Mallado
Se ha elegido una malla ortogonal estructurada con tamaos de celda uniforme,
lo que garantiza que no habr reduccin del orden de convergencia de los modelos
numricos empleados debido a las diferencias de tamaos de celda y que el acceso
a la informacin de cada volumen de control se realizar de una forma eciente. El
tamao de las celdas en cada una de las simulaciones ha sido fruto de una amplia
discusin. Por uno lado, se pretenda garantizar la convergencia del modelo numrico
y por otro no emplear tamaos de celda prohibitivos para el nmero de simulaciones
necesario en un proceso de calibracin de estas caractersticas.
Como se ha mencionado anteriormente, se han llevado a cabo ms de 200 simulaciones con varios tamaos de celda. Con esto se ha observado la sensibilidad del
modelo al tamao de celda para condiciones en las que el modelo numrico se hubiera considerado convergido en condiciones menos exigentes. Se ha decidido emplear
un mallado 2D, pese a que la turbulencia es eminentemente tridimensional.
Para el anlisis de la convergencia numrica, se han elegido dos variables y dos
94
secciones de control: velocidad y energa cintica turbulenta. Su eleccin se ha basado

en criterios de representatividad del ujo que tiene lugar. Las secciones a estudiar
son la cresta del aliviadero, donde el fenmeno de propagacin de la capa lmite
empieza a tener lugar; y el nal del aliviadero, donde el fenmeno se encuentra
completamente desarrollado.
Figura 3.2: Perles de velocidad sobre la cresta del aliviadero y en la seccin nal.
En el caso de la velocidad se observa como los perles convergen rpidamente

para resoluciones ms nas de 8 cm. En la cresta, la convergencia es todava ms
rpida. Sin embargo, a la salida del aliviadero las diferencias se han ido acumulando
y propagando dando lugar a valores ms desiguales. Pueden observarse estos resultados en la gura 3.2.
Para el caso de la energa cintica turbulenta (), se han elegido: el valor mximo
sobre la cresta y el valor en la celda que forma parte de la lmina libre al nal del
aliviadero. Pueden observarse estos resultados en la gura 3.3. Se observa que el
valor sobre la cresta es invariable. Sin embargo, la convergencia del valor al nal
del aliviadero en la celda de la lmina libre es ms sensible. Este valor, recordemos
que es la principal variable de entrada del modelo de reaireacin turbulenta y por
lo tanto este resultado explicara la alta sensibilidad del modelo al tamao de celda.
No es hasta el tamao de celda
x = 4cm
que el valor representativo al nal del
aliviadero se puede considerar convergido. Como se describir ms adelante, esto

tendr un impacto ms que signicativo en el estudio realizado.
95
Figura 3.3: Valores caractersticos de la energa cintica turbulenta sobre la cresta

del aliviadero y en la seccin nal.
Para la calibracin, se han llevado a cabo simulaciones con tamaos de celda de
x = 8
cm, donde los perles de velocidad empiezan a converger, de
x = 6
cm
donde los perles de velocidad ya han convergido pero no todava la turbulencia; y

de
x = 4 cm donde tanto las velocidades como la turbulencia se muestran estables.
Por lo tanto se recoge una situacin con bajo grado de convergencia, una con nivel
medio y otra con una convergencia muy estable.
3.3.3. Modelo de turbulencia

El modelo de turbulencia juega un papel fundamental en el problema analizado.
La turbulencia, generada en la capa lmite, se transporta con el ujo creciendo el
espesor de esta hasta que alcanza la lmina libre. En este punto, la turbulencia tiene
energa suciente para atrapar pequeas porciones de aire que pasan a ser transportadas con el ujo aumentando la concentracin de aire gradualmente hasta que
se alcanza un equilibrio. La correcta reproduccin de estas cantidades turbulentas
asegurar el buen funcionamiento del modelo de reaireacin turbulenta, descrito en
96
el Cap. 2.4.5.
El modelo de turbulencia completo ms sencillo es el modelo de turbulencia de

dos ecuaciones (Pope, 2000). Esto implica que es el modelo con mayor aplicabilidad
para un amplio rango de problemas, aunque no tiene porque ser el modelo de mayor
precisin. El modelo ms empleado es generalmente el
que utiliza una ecuacin

5
para la energa cintica turbulenta y otra para la disipacin de esta . Como ya se

ha discutido anteriormente la versin RNG de este modelo de turbulencia presenta
ciertas ventajas, por lo que es el modelo de turbulencia empleado en este estudio.
3.3.4. Modelo de lmina libre

El modelado de la lmina libre es otro de los aspectos clave del trabajo realizado.
Sin una correcta reproduccin de la evolucin de la lmina libre, la deteccin de la
localizacin del punto crtico de aireacin sera realmente imposible y por lo tanto,
podra no llegar a existir aireacin en el modelo numrico. En el desarrollo de este
trabajo de investigacin se ha empleado el modelo VOF tal y como queda descrito
en el Cap. 2.4.4.
Todas las simulaciones realizadas son 2D, al haberse comprobado previamente la

idoneidad del mtodo para reproducir correctamente la lmina libre en condiciones
ms complejas que la de este estudio (Oertel and Bung, 2012; Meireles et al., 2014;
Bombardelli et al., 2011). El modelo de turbulencia, por tratarse de un modelo
RANS, puede captar bien en un modelo 2D la turbulencia generada .
3.4. Resultados de la calibracin

La calibracin del modelos de reaireacin turbulenta se ha realizado con distintos
valores del tamao de celda, tal y como se muestra en la grca 3.3 por ejemplo. Los
x = 8 cm, x = 6 cm y x = 4 cm. En las guras

sipercies de error para el caso x = 8 cm ajustadas a los
tamaos de celda empleados son

3.4 y 3.5 se observan las
errores obtenidos con cada simulacin. Estas supercies de error permiten obtener
los valores ptimos de los parmetros.
5 Se puede encontrar una explicacin ms detallada en la seccin 2.4.3.5, en particular en 2.4.3.5.1
y 2.4.3.5.2. Se recomienda adems la lectura del captulo 10 de (Pope, 2000) y el libro (Wilcox
et al., 1998)
6 Recordemos que una de las propiedades de la turbulencia es que es esencialmente tridimensional.
97
Figura 3.4: Error obtenido en la determinacin del punto crtico de aireacin para
el caso
x = 8
cm.
98
Figura 3.5: Error obtenido en la concentracin media para el caso
x = 8
cm.
Es importante destacar que no existe ninguna combinacin de parmetros

y
Cair
que anule el error en la determinacin del punto crtico de aireacin y en
la concentracin media de aire simultneamente. Por ello, se ha dado prioridad a

la anulacin del error de la concentracin media sobre el error de la otra variable,
que simplemente se ha minimizado. De este modo, se ha obtenido la curva de error
nulo para el error en la determinacin de la concentracin media (Figura 3.5) y se
ha buscado el valor de los parmetros pertenecientes a esta curva que minimizan el
error en la determinacin del punto de inicio de aireacion (Figura 3.4).
Lo primero que se observa es como los valores de los parmetros calibrados muestran una sensibilidad muy alta al tamao de celda, pese a que la convergencia en
trminos de velocidades queda demostrada a partir de tamaos de celda
x = 8 cm
(segn se observa en la Figura 3.2).

Los valores ptimos de los parmetros calibrados se encuentran en la tabla 3.2.
Un resultado muy interesante es que cuando el valor de la turbulencia en la lmina libre alcanza un valor estable - convergencia ms exigente que en el caso de los
99
perles de velocidad nicamente - el valor de la tensin supercial amplicada por

el coeciente
toma el valor de la tensin supercial fsicamente basada del agua.
Si eliminamos este parmetro de la calibracin y por lo tanto calibramos solo
Cair
obtendremos por lo tanto otros valores distintos. El resultado queda recogido en la

tabla 3.3.
x (cm)
Cair (-)
K (N/m)
0.525
0.619
0.631
0.073
0.570
0.978
Tabla 3.2: Valores ptimos de los parmetros del modelo de reaireacin turbulenta.
x (cm)
Cair (-)
K (N/m)
Tabla 3.3: Valores ptimos de
Cair
0.525
0.291
0.334
0.073
0.073
0.073
forzando la tensin supercial a la fsicamente
basada.
De estos resultados, puede concluirse que el parmetro
no requiere de
calibracin y que una exhaustiva convergencia numrica de las simulaciones es imprescindible para el correcto funcionamiento del modelo. Se recomienda por lo tanto
realizar un test de independencia de malla basndonos siempre en los valores de la
turbulencia, no solo velocidades y calados.
100
Validacin del modelo de reaireacin turbulenta

4.1. Introduccin y objetivos
Con la validacin realizada se pretende comprobar la bondad del ajuste del modelo a otros casos no contemplados durante el proceso de calibracin. Para ello, se
han empleado nuevamente datos de prototipo libres de efectos de escala para aliviaderos con distintas pendientes. El hecho de emplear distintas pendientes garantiza
que el balance entre la turbulencia y la tensin supercial de un resultado distinto
que para el aliviadero empleado en la calibracin.
El set-up numrico de la validacin es evidentemente el mismo que el empleado
en el proceso de calibracin. No tendra sentido emplear un modelo de turbulencia
ni un modelo de lmina libre distintos.

Las opciones numricas empleadas en el proceso de validacin son las mismas
que las descritas en el Cap. 3.3, con lo que la conguracin de las simulaciones es
idntica cambiando solo la geometra de estudio. Las geometras ensayadas en las
simulaciones de validacin se correspondan con aliviaderos con inclinaciones entre
30 y 48.5 grados. Son por lo tanto, aliviaderos con menor pendiente que el caso
de calibracin y con una concentracin de aire menor. La preparacin del modelo
geomtrico se ha realizado nuevamente en AutoCAD y se ha exportado en formato
STL (estereolitogrco) con el n de poder importarlo despus en FLOW-3D.
Las simulaciones realizadas para la validacin se han llevado a cabo todas con
mallas con
x = 4
cm, que es la resolucin que asegura la convergencia de la
turbulencia a la altura de la lmina libre hacia el nal del aliviadero. Los valores
empleados de los parmetros del modelo de reaireacin turbulenta son por lo tanto:
Cair = 0,525
K = 0,073
N/m.
Captulo 4. Validacin del modelo de reaireacin turbulenta
4.3. Resultados
4.3.1. Resultados obtenidos directamente en el proceso de
calibracin
Aunque no se trata de una validacin estricta, algunos de los resultados obtenidos
durante el proceso de calibracin son objeto de un inters especial. Este es el caso de
los perles de concentracin de aire. Durante la calibracin, lo que se ha intentado
obtener es el valor exacto de la concentracin media, pero en ningn caso se ha
forzado el perl a encajar con el observado en prototipo. Se observa en la Figura
4.1 como al mejorar la resolucin del mallado se obtienen perles que recogen mejor
la forma del perl observado en prototipo. Este perl queda muy bien recogido y
descrito en (Chanson, 1993).
Figura 4.1: Perles de concentracin de aire obtenidos con los valores ptimos de los
parmetros en funcin del tamao de celda.
Otro resultado interesante se puede observar cuando se promedia la variable fraccin de uido y concentracin. Esta concentracin a la que se hace referencia en todo
104
momento, no es ms que la fraccin emulsionada de aire
Cent
(la que forma burbu-
jas). Sin embargo, existe otro tipo de aire transportado con el ujo: el aire atrapado
Crough . El aire atrapado tiene lugar debido a la gran rugosidad de la lmina libre que
permite que grandes porciones de aire sean arrastradas sin llegar a formar burbujas
(Pster, 2008). Evidentemente, tambin existe algo de aire emulsionado dentro de
la fraccin uida que rodea la fraccin de aire atrapado pero no es la fraccin predominante en la regin superior del lete uido.
En el modelo numrico, el aire atrapado no es ms que la ausencia de fraccin

de uido - y por lo tanto, ausencia tambin de cantidad de movimiento del ujo. Si
analizamos los resultados promediados, se observa como la fraccin de aire atrapado
se encuentra localizada en la regin superior del lete uido, mientras que la mayor
parte del aire emulsionado se produce en la regin inferior del perl. Adems, se
observa una transicin brusca cuando la concentracin total supera el 50 %, tal y
como se ha observado previamente en estudios experimentales (Wilhelms and Gulliver, 2005; Pster, 2008). Este cambio brusco separa las dos regiones: aire atrapado
y aire emulsionado. El resultado se puede observar en la gura 4.2.
Figura 4.2: Perles de concentracin de aire descompuesto en la fraccin emulsionada

y atrapada.
105
4.3.2. Opciones numricas
Para el proceso de validacin se han empleado las mismas opciones numrica que
en la calibracin (3.3), con el objetivo de garantizar un resultado comparable. En
cuanto al mallado empleado, se han empleado la opcin con
x = 4
cm, que ha
demostrado unos buenos indicadores de convergencia numrica.
4.3.3. Resultados obtenidos en el proceso de validacin
En el proceso de validacin se han simulado varios aliviaderos con pendientes

menores al caso de calibracin. As, las pendientes simuladas y las concentraciones
obtenidas se muestran en la tabla 4.1. En esta tabla se presentan los valores sugeridos
por varios autores y el error relativo a cada uno de ellos. Una presentacin ms visual
de los resultados se puede encontrar en la Figura 4.3.
pendiente
30.00
34.62
39.25
43.88
48.50
Cmean (FLOW-3D)
Cmean (*)
0.38
0.47
0.52
0.56
0.63
0.41
0.51
0.58
0.61
0.63
Error % (*)
7.04
8.51
10.32
8.40
0.48
Cmean
(**)
Error % (**)
0.45
0.49
0.53
0.57
0.60
14.54
4.77
1.86
1.92
-3.84
Tabla 4.1: Concentraciones de equilibrio obtenidas con FLOW-3D, propuestas por

Chanson, 1996 (*) y propuestas por Hager, 1991 (**). Se muestran tambin los
errores relativos a cada uno de los estudios experimentales.
106
Figura 4.3: Resultados del proceso de validacin.
107
Aplicacin al diseo de disipadores de energa

5.1. Objetivos
El diseo de los elementos disipadores de una obra hidrulica es, tanto desde
el punto de vista tcnico como econmico, uno de los aspectos ms costosos. Las
grandes obras hidrulicas frecuentemente deben superar grandes desniveles en un
corto recorrido haciendo necesaria una obra de disipacin de energa donde el ujo
es altamente turbulento y multifase.
El objetivo de esta parte del trabajo de investigacin es aplicar el modelo de
reaireacin turbulenta calibrado y validado al diseo de cuencos de amortiguacin.
Para ello se ha elegido un cuenco de amortiguacin segn USBR tipo II, dado su
gran rango de aplicacin. Ms informacin sobre el diseo de este tipo de cuencos
puede encontrarse en (Peterka, 1958; Hager, 1992; Vischer and Hager, 1995).
El cuenco tipo II del USBR es uno de los cuencos tipicados ms empleados. El
diseo de este cuenco es fruto de una exhaustiva campaa experimental realizada
en torno a los aos 50. En ella se consideraron diferentes opciones geomtricas que
se modelaron a escala. Como es bien sabido, los ujos aireados son muy difciles
de reproducir en laboratorio. De hecho, es fsicamente imposible escalar dos uidos
al mismo tiempo y mantener sus propiedades intactas. Con el modelo calibrado, se
pretende estudiar la estructura del ujo en este cuenco y determinar la dimensin
ptima de los dientes de pie de presa - por ser los elementos ms en contacto con el
ujo aireado.

El cuenco tipo II del USBR es uno de los cuencos tipicados ms empleados. El
diseo de este cuenco es fruto de una exhaustiva campaa experimental realizada
Captulo 5. Aplicacin al diseo de disipadores de energa
en torno a los aos 50. En ella se consideraron diferentes opciones geomtricas que
se modelaron a escala. Como es bien sabido, los ujos aireados son muy difciles de
reproducir en laboratorio. De hecho, es fsicamente imposible escalar dos uidos al
mismo tiempo y mantener sus propiedades intactas.
Figura 5.1: Dimensiones del cuenco de amortiguacin y caractersticas del ujo de

entrada.
Con el modelo de reaireacin turbulenta calibrado, se pretende estudiar la estructura del ujo desarrollada en este cuenco y determinar la dimensin ptima de
los dientes de pie de presa - por ser los primeros elementos en contacto con el ujo
supercrtico emulsionado. Ser por ello el elemento de la estructura de disipacin
ms sensible al funcionamiento del modelo de reaireacin empleado.
La forma de proceder ha sido la siguiente:
1. Anlisis en estado estacionario para las condiciones hidrulicas de diseo. Da-
dos varios tamaos de diente de pie de presa (hcb ) , se analiza su funcionamiento en estado estacionario para las condiciones de diseo. Estas condiciones
deben dar lugar a un resalto estacionario y estable dentro de la estructura de
disipacin de energa. Para alcanzar el estacionario se han simulado en una
primera etapa 60 segundos quedando en evidencia que 20 - 25 eran sucientes
para establecerse el estacionario. Se han simulado 5 segundos adicionales para
promediar la solucin y analizar sus variables hidrodinmicas.
2. Anlisis en estado transitorio para condiciones hidrulicas desfavorables. Con
el objetivo de analizar un rango mayor de situaciones, se ha realizado un
1 El
subndice de hcb hace referencia a

pie de presa.
chute block,
110
que es el nombre en ingls de los dientes de
anlisis transitorio donde el nivel de energa de aguas abajo
va decreciendo
gradualmente. Para ello se ha introducido una condicin de contorno aguas

abajo de nivel exponencial decreciente en el tiempo. Esto tiene lugar durante
30 segundos.
5.2.1. Geometra, condiciones de contorno y opciones numricas
La geometra modelada queda descrita en la gura 5.1, que est obtenida acorde
a la gua de diseo del USBR (Peterka, 1958). El modelo 3D construdo en AutoCAD se muestra en la gura 5.2. Cabe destacar que los dientes de pie de presa no
constituyen un elemento con una dimensin ja si no que su dimensin ptima es
objeto del estudio realizado. Las dimensiones de diente de pie de presa ensayadas
(hcb ) se muestran en la gura 5.3.
Las dimensiones de los dientes de pie de presa se muestran adimensionalizadas

con el calado
hw ,
que hace referencia a la fraccin de agua del calado supercrtico a
la entrada del cuenco. El valor de
hw
para el caso de estudio es de 0.141 m, calculado
acorde con las recomendaciones anteriormente referidas. El muro de nal de cuenco

se ha diseo acorde a las recomendaciones comentadas anteriormente, sin ningn
tipo de modicacin.
2 El
nivel de energa de aguas abajo viene denido por el calado y2 .

111
Figura 5.2: Geometra en 3D para el caso
hcb = 0.
Figura 5.3: Tamaos de diente de pie de presa ensayados.
Las condiciones de entrada al cuenco son especialmente particulares, tal y como

se muestra en la gura 5.1. El ujo se airea de forma natural a lo largo del aliviadero
y llega al cuenco con una concentracin dada y por lo tanto una distribucin de cantidad de movimiento muy diferente a la del resalto hidrulico clsico. El encuentro
de este perl de cantidad de movimiento con los dientes de pie de presa dar lugar
a ligeros cambios en la estructura del resalto contenido en el cuenco. Por ello es tan
importante denir la altura ptima de estos elementos que permitirn modicar la
estructura del resalto.
112
Las opciones numricas empleadas son las mismas que las descritas en el Cap.
3.3; empleadas tanto en el proceso de calibracin del modelo de reaireacin turbulenta como en su posterior validacin.
El mallado sin embargo, ha sido necesario cambiarlo por limitaciones de tiempos

de clculo. De este modo, se han empleado varias mallas 3D con tamaos de celda
entre 4 cm cerca del aliviadero y 16 cm cerca de la condicin de contorno de aguas
abajo. La progresin del tamao de celda entre mallas es de 1:2 siempre con lo que
varios ejes coinciden y se evitan as algunos problemas de redondeo en las interpolaciones. En la tabla 5.1 se muestra un resumen de las caractersticas de las 4 mallas
empleadas en cada simulacin.
Malla no.
Tipo
2D
2
3
4
(m)
No. celdas
Largo x ancho x alto (m)
0.04
262500
21.0 x 5.0 x 20.0
3D
0.04
2340000
4.8 x 5.0 x 6.0
3D
0.08
1124240
10.2 x 5.0 x 7.36
3D
0.16
84000
12.0 x 5.0 x 5.6
Tabla 5.1: Resumen de las caractersticas de las mallas empleadas en cada una de
las simulaciones de este captulo.
La malla 1 recoge la regin del aliviadero y embalse, lo que permite reproducir

las condiciones de entrada al cuenco en trminos de concentracin de aire, velocidad
y turbulencia. La malla 2, recoge los dientes de pie de presa y la zona de impacto,
ya en 3D. La malla 3 recoge la mayor parte del cuenco de amortiguacin y la malla
4 nicamente pretende alejar la condicin de contorno de aguas abajo del cuenco
con el objetivo de no alterar articialmente el ujo dentro del cuenco.
5.3. Estructura del ujo en el cuenco de amortiguacin

El uso de modelos numricos 3D permite no solo analizar los diseos hidrulicos
libres de efectos de escala si no tambin descubrir con mayor profundidad lo que
pasa exactamente dentro del ujo. De este modo, se han detectado varios mecanismos dentro de la estructura que por su inters caben destacar. En la gura 5.4 se
encuentran esquematizados estos mecanismos.
113
Figura 5.4: Estructura del ujo dentro del cuenco de amortiguacin.
Los mecanismos detectados, segn la numeracin seguida en la gura 5.4, son

los siguientes:
1. Impacto del ujo supercrtico. El ujo impacta contra la masas de agua contenida dentro del resalto, a mucha menor velocidad. La diferencia de velocidades
da lugar a una regin de cortante importante que ayuda a frenar el ujo. La
turbulencia generada en esta regin produce una aireacin mucho mayor que
la propia del ujo supercrtico entrante.
2. Deexin del chorro. El chorro se ve obligado a cambiar su trayectoria al
acercarse al codo rpida-cuenco. Adems, en esta zona el se encuentra con
los dientes de pie de presa que, como se ver ms adelante, pueden provocar
distintos efectos en funcin de su tamao.
3. Cada de la velocidad del chorro. El chorro va perdiendo velocidad a lo largo del
cuenco hasta llegar llegar a unas velocidades de salida del cuenco aceptables.
4. Regin de ujo cortante o shear region. La lnea intermitente roja de la gura
5.4 hace referencia a la zona del ujo en la que la velocidad toma el valor mitad
de la velocidad mxima. Varios autores sealan esta zona como la regin de
cortante (Rajaratnam, 1976; Hager, 1992; Chanson and Brattberg, 2000). Esta
zona es una fuente de turbulencia y disipacin de energa.
5. Recirculacin o roller. La regin de cortante da lugar a una estructura vortical
de grandes dimensiones que se sustenta sobre el chorro. Esta gran masa de
agua queda connada idealmente dentro del cuenco aunque, dentro de sus
oscilaciones, puede llegar a salirse parcialmente fuera sin que esto comprometa
114
la estabilidad de la obra de disipacin. Por ello es especialmente importante la

localizacin de los dientes de pie de presa, que adems ayudan al connamiento
de la zona de recirculacin.
6. Proteccin aguas abajo del cuenco. El chorro, aunque ha perdido mucha velocidad, sigue presentando sus velocidades ms elevadas cerca de la zona de solera.
Los dientes de nal de cuenco ayudan a cortar el chorro y deectarlo, elevando
el centro de gravedad del perl de velocidades y estableciendo una pequea
zona de recirculacin inmendiatamente aguas abajo del cuenco, con lo que el
ro se encuentra con una zona localizada de proteccin con velocidades bajas
y poca turbulencia antes de establecerse el rgimen subcrtico completamente
desarrollado.
Todos los mecanismos detectados estn de algn modo interconectados desde
que el chorro impacta (1) con la masa de agua ms lenta dentro del cuenco, hasta
que abandona el cuenco (6) en condiciones aceptables para la estabilidad del ro y de
la obra de disipacin de energa. Se trata por lo tanto de un problema de mecnica
de uidos multifase altamente complejo.
5.4. Resultados
3
Las caractersticas del resalto hidrulico clsico parcialmente sobre pendiente
han sido estudiados en detalle por varios autores (Hager, 1992; Kawagoshi and Hager, 1990; Ohtsu and Yasuda, 1991; Hager, 1988; Peterka, 1958). Este caso coincide
en gran medida con el caso de estudio
hcb = 0,
siendo las principales diferencias que
el ujo a la entrada del cuenco est aireado y que a la salida existen unos dientes que
deectan el chorro remanente. Por lo tanto, los datos experimentales documentados
por otros autores servirn de comparacin.
De este modo, cabe esperar un funcionamiento hidrulico comprendido entre el
del resalto hidrulico clsico (CHJ - Classic Hydraulic Jump ) (Chanson and Brattberg, 2000) y el del chorro de pared (CWJ - Classic Wall Jet )(Rajaratnam, 1976).
5.4.1. Condiciones hidrulicas de diseo

Para el anlisis del funcionamiento hidrulico del cuenco en condiciones de di-
seo, se han simulado 25 segundos hasta establecerse el estacionario . Una vez obtenidos estos 25 segundos, se han simulado 5 segundos adicionales para promediar
la solucin. An as, puede observarse cierta intermitencia en algunos de los perles
3 El
resalto parcialmente sobre pendiente es comnmente conocido como B-jump.

ha comprobado previamente en una simulacin de 60 segundos que a partir de 20 - 25
segundos la solucin era estadsticamente estacionaria.
4 Se
115
analizados.
Una de las variables ms relevantes que mejor resumen el comportamiento del
cuenco de amortiguacin es el parmetro de sumergencia
E,
denida como:
y2 z1
E=
y2
Donde
z1
y2
entre z1
(5.1)
queda denida en la gura 5.1 e
es el calado de aguas abajo. De

y2 , mayor ser la eciencia del
resalto (Hager, 1992; Kawagoshi and Hager, 1990). El valor del parmetro de sueste modo, cuando mayor es la diferencia
mergencia se encuentra comprendido entre 0 y 1. As, cuando
E 1,
el cuenco se
comportar como un resalto hidrulico (CHJ), por el contrario cuando
E 0
el
comportamiento se asemejar ms al de un chorro de pared (CWJ). Para las condiciones de diseo, las sumergencias tienden a ser nulas, dado que hay suciente nivel
energtico aguas abajo como para mantener el resalto en una posicin estable. En
la tabla 5.2 se muestran las distintas sumergencias obtenidas para cada tamao de
diente analizado.
hcb
0
1hw
3,1hw
7,1hw
10,7hw
y2
2.558
2.672
0.043
2.277
2.672
0.148
2.533
2.672
0.052
2.696
2.672
-0.009
2.787
2.672
-0.043
z1
(m)
Tabla 5.2: Resultados obtenidos para el parmetro de sumergencia.
Se observa como para las condiciones de diseo, el tamao de diente ptimo es

el correspondiente a
hcb = 1hw ,
tal y como se sugiere en (Peterka, 1958). En ese
estudio, sin embargo, se analizaba un rango menos amplio de valores (entre

y
2,72hw )
0,81hw
Otro de los parmetros que mejor caracteriza el comportamiento del ujo dentro
del cuenco es la cada de velocidades. Se suele denir la cada de velocidades como
la velocidad mxima en cada seccin adimensionalizada con la velocidad mxima
a la entrada al cuenco. El resultado obtenido para los escenarios ensayados queda
recogido en la gura 5.5. En esta gura, la coordenada
longitudinal adimensionalida con el calado
y95
5 En
es la distancia en el eje
de entrada al cuenco.
el estudio de Peterka se concluye que normalmente se utilizan dientes ms grandes de lo

necesario y que un tamao de diente igual a 1hw es suciente. Sin embargo no se habla en ningn
momento de un valor ptimo.
116
Figura 5.5: Cada de la velocidad mxima dentro del cuenco.
Se puede observar como el comportamiento de la cada de velocidades es ms

parecida a la del chorro de pared como cabe de esperar para este valor del parmetro
de sumergencia (Hager, 1992; Kawagoshi and Hager, 1990). Se oberva tambin en
la gura 5.5 una zona de aceleracin del ujo en la zona de impacto, tal y como se
observa en estudios anteriores (Ohtsu and Yasuda, 1991). El ajuste de la curva para
hcb = 0 a la del CWJ es muestra del correcto funcionamiento del modelo numrico.
Es destacable que para los dientes mayores a
hcb = 1hw
la cada de velocidades es
pronunciada en el primer tramo del cuenco, donde el ujo encuentra un obstculo que
decelera su movimiento. Sin embargo, para el caso de
hcb = 1hw
el ujo se comporta
en este primer tramo de una forma similar al caso sin diente. Pero, a diferencia de
este caso ms bsico, el chorro cambia su comportamiento a mitad de cuenco y la
cada de velocidades se hace ms pronunciada. De esta gura es difcil extraer un
candidato para el valor ptimo, pero permite identicar distintos comportamientos
del chorro:
Para dientes grandes, el chorro impacta con el diente. Esto produce una cada
localizada de las velocidades que no se acompaa a lo largo del cuenco.
Para dientes pequeos, el chorro se ve corrugado por el diente, introduciendo
cierta cantidad de turbulencia que acaba frenando el chorro a lo largo del
cuenco.
Otro conclusin interesante es la que se puede concluir del estudio de la regin de
cortante o shear region, mostrada en la gura 5.6. En esta gura,
hace referencia
a la altura sobre la solera donde se encuentra la velocidad mitad de la velocidad
117
mxima
v1
para la seccin analizada.
Se observa como todos los casos presentan un comportamiento parecido, salvo

el caso
hcb = 1hw
donde el chorro queda ms pegado a la solera. Esto quiere decir
que en este caso contamos con una recirculacin o roller de mayores dimensiones.
Al mismo tiempo, al encontrarse ms bajo el centro de gravedad del perl de velocidades, la distribucin de cantidad de movimiento impacta de una forma ms directa
con los dientes del nal de cuenco.
Figura 5.6: Regin de cortante o shear region.
5.4.2. Condiciones hidrulicas desfavorables

Uno de los papeles de los dientes de pie de presa es estabilizar el resalto (Peterka,
1958), por ello es especialmente interesante analizar el funcionamiento del cuenco en
condiciones desfavorables. Para ello, y dadas las limitaciones de tiempo de clculo
que supondra analizar innitas combinaciones, se han llevado a cabo simulaciones
en estado estacionario con una condicin de contorno de aguas abajo variable en
el tiempo. Partiendo del calado de aguas abajo correspondiente a la situacin de
diseo y de la solucin obtenida para el apartado anterior, la condicin de contorno
de aguas abajo va perdiendo calado exponencialmente durante 30 segundos hasta
terminar por debajo del 40 % del valor de diseo. El resultado obtenido se muestra
en la gura 5.7.
118
Figura 5.7: Anlsis transitorio del funcionamiento hidrulico del cuenco de amortiguacin.
Se observa claramente como los valores extremos del tamao de diente dan lugar
a las soluciones menos estables. Sin embargo, el resto de valores intermedios dan unos
valores parecidos. Teniendo en cuenta que el resalto contenido tiene unas oscilaciones
con cierto periodo caracterstico asimilable a las diferencias de tiempo de despegue
de las 3 soluciones intermedias, sera imposible concluir con gran certeza cual de las
3 soluciones es la ms estable. Si bien se puede concluir que los 3 tamaos de diente
intermedios presentan un tiempo de despegue razonablemente parecido (entre 26 y
28 segundos). En la gura 5.7 tambin se muestra la serie temporal del parmetro de
sumergencia. Se observa como la lnea correspondiente a
hcb = 1hw
muestra durante
mayor tiempo una sumergencia claramente superior y por lo tanto una eciencia
mejor a la del resto de casos. Esto permitira concluir que el tamao de diente
hcb = 1hw
es el ms eciente para condiciones hidrulicas adversas.
119
Conclusiones
Los ujos turbulentos multifase son comunes en muchsimas aplicaciones ingenieriles y medioambientales. La naturaleza estocstica tanto de la fase portante como
de la fase dispersa hace del problema multifase un problema innitamente ms complejo que su anlogo de una sola fase (Balachandar and Eaton, 2010). La evidencia
tanto numrica como experimental de la existencia de un reajuste de la turbulencia
debido a la presencia de una fase portante hace que el problema adquiera un aspecto todava enigmtico y difcil de predecir. Es recomendable la lectura de trabajos
como (Zhou, 2010; Yeo et al., 2010; Crowe, 2010; Prosperetti and Tryggvason, 2007;
Chen and Faeth, 2001; Crowe, 2000; Kussin and Sommerfeld, 2002) donde se dan
evidencias sucientes de este fenmeno.
No hay que perder de vista que la modelacin numrica de ujos multifase es

sin duda un campo muy amplio, que puede abarcar desde la ingeniera qumica,
petroqumica, trmica o hidrulica. Desde el punto de vista hidrulico, son probablemente el transporte de sedimentos y el ujo aireado en aliviaderos las dos
aplicaciones de mayor inters. Las caractersticas de los ujos aireados en aliviaderos son especialmente interesantes pues nos encontramos ante ujos supercrticos
altamente aireados con una granulometra de burbujas muy amplia que a su mismo
tiempo puede evolucionar en el tiempo debido a la turbulencia y presiones locales
que sufren las burbujas. Los ujos supercrticos de por s ya presentan un reto para
la modelacin numrica. As, esta vertiente del problema multifase en la ingeniera
hidrulica hace sin ninguna duda los ujos en aliviaderos uno de los mayores retos
con los que nos encontramos en el panorama hidrulico internacional.
Es por lo tanto, el enfoque que se ha dado al presente trabajo, un enfoque ingenieril donde prima la obtencin de un resultado al conocimiento cientco completo del
fenmeno estudiado. El empleo de un modelo de subescala permite aproximar lo que
ocurre dentro de una celda sin tener completamente resuelta la hidrodinmica alrededor de las burbujas (o partculas). El acoplamiento con un modelo de turbulencia
de dos ecuaciones, el modelo completo ms sencillo (Pope, 2000), pone el broche a
Captulo 6. Conclusiones
este enfoque ingenieril del problema de una forma elegante y efectiva. Sin embargo,
la falta de validacin es sin ninguna duda uno de los problemas ms apremiantes en
el campo de la modelacin numrica de los ujos aireados (Chanson, 2013; Chanson
and Lubin, 2010). Por ello este trabajo ha partido desde un principio de la slida base que proporciona una calibracin y una posterior validacin a un modelo numrico.
El modelo de reaireacin turbulenta, descrito en el captulo 2.4.5, es calibrado
en el captulo 3 y posteriormente validado en el captulo 4. Los resultados de este
ltimo captulo muestran un buen ajuste para el rango de pendientes analizado (30
- 53 ). Lo que nos anima a testearlo en una aplicacin con cierto inters ingenieril.
Esta aplicacin, y como se ve en los resultados del captulo 5, podramos decir que
se convierte en una segunda validacin del modelo.
La aplicacin prctica planteada es la de encontrar una altura ptima de los
dientes de pie de presa, que juegan un papel estabilizador del resalto hidrulico en
el cuenco de amortiguamiento (Vallarino, 1991). El rango de tamaos ensayados
queda recogido en la gura 5.3. Las simulaciones realizadas son 3D manteniendo
una relacin entre el calado y el ancho de 1:5 con el objetivo de no perder la tridimensionalidad del ujo; tal y como sugiere el Dr. Hager tanto en (Hager, 2014)
como en correspondencia mantenida personalmente. Esto ha servido tambin para
validar la calibracin, que se haba llevado a cabo en 2D, al no observar diferencias
apreciables en las velocidades y concentraciones de aire (tal y como caba de esperar).
Sorprendentemente, de todas las alturas de diente de pie de presa analizadas, la
ptima ha resultado ser la misma que se sugera aos atrs en los primeros estudios
sobre elementos de disipacin de energa (Peterka, 1958). En ellos, se examinaba
un rango ms pequeo de dientes y no se pretenda optimizar si no dar una altura
mnima vlida. Las simulaciones conducidas tanto en condiciones de diseo (Cap.
5.4.1) como en condiciones desfavorables (Cap. 5.4.2) han mostrado la idoneidad
de este tamao de diente y su papel estabilizador del resalto. Ha servido tambin
para identicar el papel clave que juega al no funcionar como un bloque de impacto
si no como una macrorrugosidad. A modo de calibracin ltima del modelo se han
comparado los resultados con los datos experimentales para resalto hidrulico sobre
pared lisa y resalto hidrulico sobre pendiente tipo B.
122
Referencias
Aivazyan, O. (1986). Stabilized aeration on chutes. Hydrotechnical Construction,

20(12):713722.
Amsden, A. A. and Harlow, F. H. (1970). A simplied mac technique for incompressible uid ow calculations. Journal of Computational Physics, 6(2):322325.
Balachandar, S. and Eaton, J. K. (2010).
Turbulent dispersed multiphase ow.
Annual Review of Fluid Mechanics, 42:111133.

Bauer, W. J. (1954). Turbulent boundary layer on steep slopes. Transactions of the
American Society of Civil Engineers, 119(1):12121233.

Blazek, J. (2005).
Computational Fluid Dynamics: Principles and Applications:
Principles and Applications. Second Edition. Elsevier.

Bombardelli, F. A. and Chanson, H. (2009). Progress in the observation and modeling of turbulent multi-phase ows. Environmental Fluid Mechanics, 9(2):121123.
Bombardelli, F. A., Meireles, I., and Matos, J. (2011). Laboratory measurements and
multi-block numerical simulations of the mean ow and turbulence in the nonaerated skimming ow region of steep stepped spillways.
Environmental Fluid
Mechanics, 11(3):263288.
Boussinesq, J. (1877). Essai sur la thorie des eaux courantes, volume 2. Imprimerie
nationale.
Bung, D. B. (2011).
Developing ow in skimming ow regime on embankment
stepped spillways. Journal of Hydraulic Research, 49(5):639648.

Bung, D. B. (2013). Non-intrusive detection of airwater surface roughness in selfaerated chute ows. Journal of Hydraulic Research, 51(3):322329.
Cain, P. (1978). Measurements within self-aerated ow on a large spillway.
125
Cain, P. and Wood, I. R. (1981). Measurements of self-aerated ow on a spillway.
Journal of the Hydraulics Division, 107(11):14251444.

Caisley, M. E., Bombardelli, F. A., and Garcia, M. H. (1999).
Hydraulic model
study of a canoe chute for low-head dams in illinois. Civil Engineering Studies,
Hydraulic Engineering Series No-63, University of Illinois at Urbana-Champaign.

Castro-Orgaz, O. (2009). Velocity prole and ow resistance models for developing
chute ow. Journal of Hydraulic Engineering, 136(7):447452.
Chanel, P. and Doering, J. (2007). An evaluation of computational uid dynamics for
spillway modelling. In 16th Australasian Fluid Mechanics Conference (AFMC),
pages 12011206. School of Engineering, The University of Queensland.
Chanson, H. (1993). Self-aerated ows on chutes and spillways. Journal of hydraulic
engineering, 119(2):220243.
Chanson, H. (1994). Drag reduction in open channel ow by aeration and suspended
load. Journal of Hydraulic Research, 32(1):87101.
Chanson, H. (1996). Air bubble entrainment in free-surface turbulent shear ows.
Access Online via Elsevier.
Chanson, H. (2009).
Turbulent airwater ows in hydraulic structures: dynamic
similarity and scale eects. Environmental Fluid Mechanics, 9(2):125142.

Chanson, H. (2013). Hydraulics of aerated ows: qui pro quo? Journal of Hydraulic
Research, 51(3):223243.
Chanson, H. and Brattberg, T. (2000). Experimental study of the airwater shear
ow in a hydraulic jump. International Journal of Multiphase Flow, 26(4):583
607.
Chanson, H. and Lubin, P. (2010). Discussion of verication and validation of a
computational uid dynamics (cfd) model for air entrainment at spillway aerators
appears in the canadian journal of civil engineering 36 (5): 826-838.
Canadian
Journal of Civil Engineering, 37(1):135138.

Chen, J.-H. and Faeth, G. (2001).
Continuous-phase properties of homogeneous
particle-laden turbulent ows. AIAA journal, 39(1):180183.

Cheng, N.-S. (2007).
Power-law index for velocity proles in open channel ows.
Advances in water Resources, 30(8):17751784.

Crowe, C. T. (2000). On models for turbulence modulation in uidparticle ows.
International Journal of Multiphase Flow, 26(5):719727.

Crowe, C. T. (2010). Multiphase ow handbook. CRC Press.
Ervine, DA, F. H. (1987). Behaviour of turbulent water jets in the atmosphere and
in plunge pools. In ICE proceedings, volume 83, pages 295314. Thomas Telford.
Falvey, H. R. and Ervine, D. A. (1988). Aeration in jets and high velocity ows. In
Model-prototype correlation of hydraulic structures, pages 2255. ASCE.

Falvey, H. T. (1980). Air-water ow in hydraulic structures. Bureau of Reclamation,
Engineering Monograph No. 41.
Falvey, H. T. (1990). Cavitation in chutes and spillways. Bureau of Reclamation,
Engineering Monograph No. 42.
Fernndez Oro, J. M. (2012). Tcnicas numricas en Ingeniera de Fluidos.
Gingold, R. A. and Monaghan, J. J. (1977).
Smoothed particle hydrodynamics:
theory and application to non-spherical stars. Monthly notices of the royal astro-
nomical society, 181(3):375389.

Gutirrez Serret, R. (1995). Aireacin en las estructuras hidrulicas de las presas:
aliviaderos y desages profundos.
PhD thesis, Tesis doctoral, Biblioteca ETS.
Ingenieros de Caminos, Universidad Politcnica de Madrid. Documento 360718.

Hager, W. (1991a). Experiments on standard spillway ow (including appendix). In
ICE Proceedings, volume 91, pages 399416. Thomas Telford.

Hager, W. H. (1988). B-jump in sloping channel. Journal of Hydraulic Research,
26(5):539558.
Hager, W. H. (1991b). Uniform aerated chute ow. Journal of Hydraulic Enginee-
ring, 117(4):528533.
Hager, W. H. (1992). Energy dissipators and hydraulic jump, volume 8. Springer.
Hager, W. H. (2014). Hydraulic structures: a positive outlook into de future. Journal
of Hydraulic Research.
Harlow, F. H. and Welch, E. J. (1970).
Numerical calculation of time dependent
viscous incompressible ow of uid with free surface. Journal of Computational
Physics, 6(2):322325.
Heller, V. (2011). Scale eects in physical hydraulic engineering models. Journal of
Hydraulic Research, 49(3):293306.

Hinze, J. O. (1975). Turbulence. New York: McGraw-Hill.
Hirsch, C. (2007).
Numerical Computation of Internal and External Flows: The
Fundamentals of Computational Fluid Dynamics: The Fundamentals of Computational Fluid Dynamics. Butterworth-Heinemann.
Hirt, C. W. and Nichols, B. D. (1981). Volume of uid (vof ) method for the dynamics
of free boundaries. Journal of computational physics, 39(1):201225.
Ho, D., Boyes, K., and Donohoo, S. (2001). Investigation of spillway behavior under
increased maximum ood by computational uid dynamics technique. In Proc.
Conf. 14th Australasian Fluid Mechanics, pages 577580.

Huang, P., Bardina, J., and Coakley, T. (1997).
Turbulence modeling validation,
testing, and development. NASA Technical Memorandum, 110446.

Jang, D., Jetli, R., and Acharya, S. (1986). Comparison of the piso, simpler, and
simplec algorithms for the treatment of the pressure-velocity coupling in steady
ow problems. Numerical Heat Transfer, Part A: Applications, 10(3):209228.
Jevdjevich, V. and Levin, L. (1953). Entrainment of air in owing water and technical problems connected with it. In Proceedings@ sMinnesota International Hy-
draulic Convention, pages 439454. ASCE.

Kawagoshi, N. and Hager, W. (1990).
B-jump in sloping channel, ii.
Journal of
Hydraulic Research, 28(4):461480.

Kramer, K. (2004). Development of aerated Chute Flow. VAW, PhD Thesis.
Kramer, K., Hager, W. H., and Minor, H.-E. (2006). Development of air concentration on chute spillways. Journal of Hydraulic Engineering, 132(9):908915.
Kussin, J. and Sommerfeld, M. (2002). Experimental studies on particle behaviour
and turbulence modication in horizontal channel ow with dierent wall roughness. Experiments in Fluids, 33(1):143159.
Launder, B. and Spalding, D. (1972). Mathematical models of turbulence, 1972.
Launder, B. E. and Spalding, D. (1974). The numerical computation of turbulent
ows. Computer methods in applied mechanics and engineering, 3(2):269289.
Ma, J., Oberai, A., Drew, D., Lahey Jr, R., and Moraga, F. (2010). A quantitative
sub-grid air entrainment model for bubbly owsplunging jets.
Computers &
Fluids, 39(1):7786.
Ma, J., Oberai, A. A., Drew, D. A., Lahey, R. T., and Hyman, M. C. (2011a). A
comprehensive sub-grid air entrainment model for rans modeling of free-surface
bubbly ows. The Journal of Computational Multiphase Flows, 3(1):4156.
Ma, J., Oberai, A. A., Lahey Jr, R. T., and Drew, D. A. (2011b).
Modeling air
entrainment and transport in a hydraulic jump using two-uid rans and des turbulence models. Heat and mass transfer, 47(8):911919.
Madavan, N., Deutsch, S., and Merkle, C. (1983).
Reduction of turbulent skin
friction by microbubbles. Technical report, DTIC Document.
MAGRAMA (2006). Inventario de presas y embalses de espaa.

McKee, S., Tom, M., Ferreira, V., Cuminato, J., Castelo, A., Sousa, F., and Mangiavacchi, N. (2008). The mac method. Computers & Fluids, 37(8):907930.
Meireles, I. C., Bombardelli, F. A., and Matos, J. (2014).
Air entrainment onset
in skimming ows on steep stepped spillways: an analysis. Journal of Hydraulic
Research, (ahead-of-print):111.
Monaghan, J. (1985). Particle methods for hydrodynamics. Computer Physics Re-
ports, 3(2):71124.
Monaghan, J. and Gingold, R. (1983). Shock simulation by the particle method sph.
Journal of computational physics, 52(2):374389.

Monaghan, J. J. (1992). Smoothed particle hydrodynamics. Annual review of as-
tronomy and astrophysics, 30:543574.

Nezu, I., Tominaga, A., and Nakagawa, H. (1993). Field measurements of secondary
currents in straight rivers. Journal of Hydraulic Engineering, 119(5):598614.
Nichols, B., Hirt, C., and Hotchkiss, R. (1980). Sola-vof: A solution algorithm for
transient uid ow with multiple free boundaries. Technical report, Los Alamos
Scientic Lab., NM (USA).
Novak, P., Moat, A., Nalluri, C., and Narayanan, R. (2007). Hydraulic structures.
Taylor & Francis.
Oertel, M. and Bung, D. B. (2012).
Initial stage of two-dimensional dam-break
waves: laboratory versus vof. Journal of Hydraulic Research, 50(1):8997.

Ohtsu, I. and Yasuda, Y. (1991). Hydraulic jump in sloping channels. Journal of
Hydraulic Engineering, 117(7):905921.

Patankar, S. (1980). Numerical heat transfer and uid ow. CRC Press.
Patankar, S. V. and Spalding, D. B. (1972).
A calculation procedure for heat,
mass and momentum transfer in three-dimensional parabolic ows. International
Journal of Heat and Mass Transfer, 15(10):17871806.

Peterka, A. J. (1958).
Hydraulic design of stilling basins and energy dissipators.
Bureau of Reclamation, Engineering Monograph No. 25.

Pster, M. (2008).
Discussion to 'bubbles and waves description of self-aerated
spillway ow'. Journal of Hydraulic Research, 46(3):420423.

Pster, M. and Hager, W. H. (2010a). Chute aerators. i: Air transport characteristics. Journal of Hydraulic Engineering, 136(6):352359.
Pster, M. and Hager, W. H. (2010b). Chute aerators. ii: Hydraulic design. Journal
of Hydraulic Engineering, 136(6):360367.

Pster, M., Hager, W. H., and Minor, H.-E. (2006).
Bottom aeration of stepped
spillways. Journal of Hydraulic Engineering, 132(8):850853.

Pope, S. B. (2000). Turbulent ows. Cambridge university press.
Prosperetti, A. and Tryggvason, G. (2007). Computational methods for multiphase
ow. Cambridge University Press.

Rajaratnam, N. (1976). Turbulent jets. Elsevier.
Saad, Y. (2003). Iterative methods for sparse linear systems. Siam.
Savage, B. M. and Johnson, M. C. (2001). Flow over ogee spillway: Physical and
numerical model case study. Journal of Hydraulic Engineering, 127(8):640649.
Schlichting, H. (1979). Boundary-Layer Theory. New York: McGraw-Hill.
Schmocker, L., Pster, M., Hager, W. H., and Minor, H.-E. (2008). Aeration characteristics of ski jump jets. Journal of Hydraulic Engineering, 134(1):9097.
Smith, M. and Smits, A. (1995). Visualization of the structure of supersonic turbulent boundary layers. Experiments in uids, 18(4):288302.
Straub, L. G. and Anderson, A. G. (1958). Experiments on self-aerated ow in open
channels. Journal of the Hydraulics Division, 84(7):135.
USBR (1977). Design of small dams. Washington. DC: Govt. Print. O.
Valero, D., Fullana, O., Garca-Bartual, R., Andrs-Domnech, I., Valls, F., and
Marco, J. (2014). Analytical formulation for the aerated hydraulic jump and physical modeling comparison. In 3rd Europe IAHR Congress, Book of Proceedings,
2014, Porto - Portugal.

Valero Huerta, D., Fullana Montoro, O., Garca Bartual, R., Valls Morn, F. J., and
Andrs-Domnech, I. (2013). Ejecucin de un nuevo dispositivo experimental para
la caracterizacin de resaltos hidrulicos emulsionados. In VI Seminario sobre las
Lneas Prioritarias de Investigacin de la Red de Laboratorios de Hidrulica de

Espaa, RLHE.
Vallarino, E. (1991). Tratado bsico de presas. Colegio de Ingenieros de Caminos,
Canales y Puertos.
Van Doormaal, J. and Raithby, G. (1984). Enhancements of the simple method for
predicting incompressible uid ows. Numerical heat transfer, 7(2):147163.
Versteeg, H. K. and Malalasekera, W. (2007). An introduction to computational uid
dynamics: the nite volume method. Pearson Education.
Vischer, D. and Hager, W. (1995). Energy dissipators, hydraulic structures design
manual. AA Balkema.
Volkart, P. (1980). The mechanism of air bubble entrainment in self-aerated ow.
International Journal of Multiphase Flow, 6(5):411423.

Wilcox, D. C. et al. (1998). Turbulence modeling for CFD, volume 2. DCW industries
La Canada, CA.
Wilhelms, S. C. and Gulliver, J. S. (2005). Bubbles and waves description of selfaerated spillway ow. Journal of Hydraulic Research, 43(5):522531.
Wood, I. (1991).
Free surface air entrainment on spillways.
Number 4. IAHR
Hydraulic Structures Design Manual, Balkema, Rotterdam, The Netherlands.

Wood, I. R., Ackers, P., and Loveless, J. (1983). General method for critical point
on spillways. Journal of Hydraulic Engineering, 109(2):308312.
Wu, W., Rodi, W., and Wenka, T. (2000).
3d numerical modeling of ow and
sediment transport in open channels. Journal of Hydraulic Engineering, 126(1):4

15.
Yakhot, V., Orszag, S., Thangam, S., Gatski, T., and Speziale, C. (1992).
Deve-
lopment of turbulence models for shear ows by a double expansion technique.
Physics of Fluids A: Fluid Dynamics (1989-1993), 4(7):15101520.

Yakhot, V. and Orszag, S. A. (1986). Renormalization group analysis of turbulence.
i. basic theory. Journal of scientic computing, 1(1):351.
Yeo, K., Dong, S., Climent, E., and Maxey, M. R. (2010). Modulation of homogeneous turbulence seeded with nite size bubbles or particles. International Journal
of Multiphase Flow, 36(3):221233.

Zhou, L. (2010). Advances in studies on two-phase turbulence in dispersed multiphase ows. International Journal of Multiphase Flow, 36(2):100108.

TFM - Daniel Valero PDF

Cargado por

Copyright:

Formatos disponibles

TFM - Daniel Valero PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

TFM - Daniel Valero PDF

Cargado por

Copyright:

Formatos disponibles

Modelacin hidrulica de flujos

multifase en grandes presas.

Rafael Garca Bartual

Daniel Valero Huerta

Tesina Final de Mster, tipo B

Valencia, Abril de 2015

2. Estado del Arte

Las presas en Espaa y en el mundo

Dinmica de Fluidos Computacional (CFD)

Conservacin de la cantidad de movimiento

Sistema de ecuaciones de Navier-Stokes . . . . . . . .

Naturaleza del trmino de presiones . . . . . . . . . . . . . . .

Ecuaciones de gobierno para ujos multifase

Naturaleza de las ecuaciones de Navier-Stokes

Clasicacin de las ecuaciones del ujo . . . . . . . .

Ecuaciones de Reynolds o del ujo medio . . . . . . . . . . . .

Problema del cierre . . . . . . . . . . . . . . . . . . . . . . . .

Hidrulica de grandes presas . . . . . . . . . . . . . . . . . . . . . . .

Determinacin del punto de inicio de la aireacin

Reduccin de friccin por el efecto de la aireacin . .

El Mtodo de los Volmenes Finitos . . . . . . . . . . . . . . .

Grados de aproximacin a la resolucin 3D de las ecuaciones

Hiptesis Boussinesq o de viscosidad turbulenta . . .

La ecuacin de la energa cintica turbulenta . . . . .

Modelo de van Driest . . . . . . . . . . . . .

Modelo de Baldwin y Lomax

Modelado de la lmina libre

Modelado de la reaireacin turbulenta

3. Calibracin del modelo de reaireacin turbulenta

Modelo de lmina libre . . . . . . . . . . . . . . . . . . . . . .

4. Validacin del modelo de reaireacin turbulenta

Descripcin general . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103

Resultados obtenidos directamente en el proceso de calibracin 104

Opciones numricas . . . . . . . . . . . . . . . . . . . . . . . . 106

Resultados obtenidos en el proceso de validacin . . . . . . . . 106

5. Aplicacin al diseo de disipadores de energa

Descripcin general . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109

Geometra, condiciones de contorno y opciones numricas . . . 111

Estructura del ujo en el cuenco de amortiguacin . . . . . . . . . . . 113

Condiciones hidrulicas de diseo . . . . . . . . . . . . . . . . 115

Condiciones hidrulicas desfavorables . . . . . . . . . . . . . . 118

Flujo supercrtico reproducido en laboratorio (Pster, 2008). . . . . .

Flujo supercrtico emulsionado sobre aliviadero y aireadores (presa de

Simulacin 3D de la regin no aireada sobre un perl WES (Ho et al.,

Diferentes secciones cnicas en funcin del valor del discriminante. . .

Dominio de dependencia de las ecuaciones elpticas. . . . . . . . . . .

Dominio de dependencia de las ecuaciones parablicas.

Dominio de dependencia de las ecuaciones hiperblicas. . . . . . . . .

Velocidad en un ujo turbulento y turbulencia instantnea en el punto

Clculo de la velocidad en un aliviadero segn USBR (Peterka, 1958).

Campo de velocidades en el entorno de la cresta para

2.10. Campo de velocidades en el entorno de la cresta para

2.13. Determinacin del punto crtico de aireacin: interseccin de la capa

2.15. Estructura del ujo supercrtico aireado.

2.16. a) Concentracin media b) Nomenclatura grcos c) Concentracin

2.17. Perl de concentraciones terico y medido en prototipo. . . . . . . . .

2.18. Perl de velocidades terico y medido en prototipo.

2.19. Evolucin de las concentraciones en un aliviadero. . . . . . . . . . . .

2.20. Reduccin de friccin por el efecto del aire en el ujo. . . . . . . . . .

Ecuaciones de gobierno para ujos multifase

Clasicacin de las ecuaciones del ujo . . . . . . . .

Ecuaciones de Reynolds o del ujo medio . . . . . . . . . . . .

Estructura del ujo en el cuenco de amortiguacin . . . . . . . . . . . 113

Flujo supercrtico reproducido en laboratorio (Pster, 2008). . . . . .

Simulacin 3D de la regin no aireada sobre un perl WES (Ho et al.,

Velocidad en un ujo turbulento y turbulencia instantnea en el punto

2.15. Estructura del ujo supercrtico aireado.

2.16. a) Concentracin media b) Nomenclatura grcos c) Concentracin

2.17. Perl de concentraciones terico y medido en prototipo. . . . . . . . .

2.18. Perl de velocidades terico y medido en prototipo.

2.20. Reduccin de friccin por el efecto del aire en el ujo. . . . . . . . . .

Perles de velocidad sobre la cresta del aliviadero y en la seccin nal. 95

Perles de concentracin de aire obtenidos con los valores ptimos de

Perles de concentracin de aire descompuesto en la fraccin emulsionada y atrapada. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105

Dimensiones del cuenco de amortiguacin y caractersticas del ujo

Estructura del ujo dentro del cuenco de amortiguacin.

Clasicacin de los principales tipos de ujo. . . . . . . . . . . . . . .

Concentraciones de equilibrio en el ujo completamente desarrollado

forzando la tensin supercial a la fsicamente

Concentracin de aire emulsionado para el ujo completamente desarrollado.

Concentracin de media en el lete uido.

Vector de ujos convectivos.

Vector de ujos difusivos.

Flujo msico a travs de una supercie.

Supercie de separacin entre dos fases.