07 Ocma 07 de 13

3.
1
Captulo 3. Transporte de Sedimentos
3.- Transporte de Sedimentos

En este captulo se realiza una revisin y descripcin de algunos de los procesos fsicos
hidromorfodinmicos ms relevantes en la evolucin de las playas arenosas.
3.1.- Introduccin
Las playas arenosas han sido tradicionalmente un lugar de gran inters recreativo y
socio-econmico para el hombre. Prueba de ello es que cerca de las dos terceras partes
de la poblacin mundial vive en o cerca de la costa. Quizs la caracterstica ms
importante de las playas arenosas, desde el punto de vista fsico, es la sensibilidad y
capacidad de respuesta morfolgica a las condiciones hidrodinmicas; la arena es
transportada constantemente por la accin del oleaje, de las corrientes y por el viento.
Sin embargo, la accin del hombre sobre su entorno ha dado lugar a grandes cambios en
la zona costera, e.g. la construccin de presas, puertos, espigones y otros tipos de
estructuras de estabilizacin (la mayora, muy costosas y de escaso valor esttico). Estos
factores, junto con el hecho de que grandes cantidades de aguas residuales sean vertidas
al mar, pueden provocar un impacto econmico importante en las poblaciones costeras.
Otra funcin directa igualmente importante de una playa, es la proteccin a la zona
costera, ya que aqu se concentra la energa del oleaje, donde se disipa en una zona
relativamente estrecha cerca de la superficie del agua. En una tormenta, si la playa tiene
la cantidad suficiente de arena, puede modificar su forma, generando una serie de barras
sumergidas que pueden provocar que las olas grandes rompan y disipen su energa antes
de llegar a la costa. Las actuaciones antes mencionadas del hombre en la zona costera,
pueden minar y/o interrumpir el transporte de sedimentos, dando lugar a una
intensificacin de la accin del oleaje y a una mayor erosin. Si una playa se llega a
destruir, repararla o reconstruirla es difcil y el proceso es muy costoso.
El estudio de los procesos fsicos en la zona costera ha sido abordado tradicionalmente
por dos escuelas con puntos de vista muy diferentes. Los estudios sedimentolgicos y
de morfologa (bsicamente cualitativos, aunque no por ello menos interesantes o
importantes) han tenido un desarrollo importante, desde comienzos de la segunda mitad
del siglo XX, llevados a cabo por gelogos y geomorflogos costeros. Por otro lado, el
estudio cuantitativo de los procesos fsicos hidromorfodinmicos en la zona cercana a
la costa, ha sido realizado en profundidad slo durante los ltimos 30 aos
(aproximadamente) por investigadores dentro de los campos de la ingeniera y la
oceanografa costeras. Los fenmenos observados en el rea cercana a la costa son
extremadamente complejos y no es fcil ningn tipo de tratamiento analtico. A finales
de los aos 40 comenzaron los estudios de estos problemas a travs de experimentos de
laboratorio y de campo, y se han encontrado numerosas relaciones o correlaciones semiempricas entre los fenmenos fsicos que rigen la dinmica costera.
3.2
Desde la perspectiva de la ingeniera, los fenmenos costeros sedimentarios estn
estrechamente relacionados con varios problemas prcticos importantes, tales como la
sedimentacin de las drsenas de los puertos o la erosin de playas. Hasta no hace
mucho tiempo, la mayora de los proyectos en el litoral eran llevados a cabo mediante el
procedimiento de prueba y error, debido a la falta de conocimiento de los mecanismos
que rigen los procesos costeros. Muchas veces, el clculo del transporte de sedimentos
se realiza mediante el uso de una gran cantidad de ecuaciones que estn calibradas para
condiciones hidrodinmicas muy especficas o que no tienen un fundamento fsico real
(ecuaciones empricas). Por lo anterior y debido a que se ha extendido el desarrollo de
modelos numricos para la evaluacin del transporte de sedimentos y de evolucin del
fondo, es necesaria la obtencin y anlisis de una gran cantidad de datos
experimentales. Durante las ltimas dcadas, se han tomado grandes cantidades de datos
en estudios del transporte de sedimentos a travs de investigaciones de campo y de
laboratorio. Aunque todos estos datos son de gran ayuda, los procesos fsicos en la zona
cercana a la costa son de una complejidad extraordinaria, y la mayora siguen siendo tan
solo comprendidos en un sentido cualitativo. En este capitulo se lleva a cabo una
revisin de algunos de los estudios relacionados con la presente tesis desde el punto de
vista morfodinmico.
3.2.- Factores que regulan la variacin de la morfologa de playas

A pesar de que las playas varan ampliamente en apariencia y composicin, los
principios que gobiernan sus comportamientos son los mismos. El mecanismo bsico
del proceso de cambio es el levantamiento de millones de granos de arena del fondo
debido a la friccin (tensiones de corte) generada principalmente por el flujo oscilatorio
asociado al oleaje, la friccin de las corrientes y/o por la turbulencia generada al paso y
rotura de la ola.
Mientras el sedimento se encuentra en suspensin, es susceptible de ser transportado a
lo largo de la costa por las corrientes inducidas por el oleaje o por el viento. Por esto, las
playas constantemente estn cambiando de forma. Este transporte litoral a lo largo de la
costa juega un papel importante en la evolucin de su perfil transversal y es tambin el
responsable de la mayora de los cambios en la lnea de orilla. Puesto que el uso de las
zonas costeras por parte del hombre (usos tursticos, comerciales, etc.) se ha
incrementado de una manera importante durante el ltimo siglo, la alteracin de la
dinmica costera puede representar un grave problema para el hombre.
Se dice que una playa est en equilibrio cuando el volumen total de sedimentos movido
en la misma, en un determinado intervalo de tiempo no presenta prdidas ni ganancias,
siendo precisamente este equilibrio en playas lo que el hombre ha querido mantener por
diferentes medios para el mejor aprovechamiento de ellas. La erosin, acrecin o
estabilidad de las costas depende del balance sedimentario (las tasas de aporte o
remocin de sedimento de la costa). El concepto de balance sedimentario consiste
bsicamente en la aplicacin del principio de continuidad o conservacin de masa al
sedimento litoral de un rea conocida, y slo puede ser aplicado a aquellas situaciones
donde estn involucrados el material sedimentario y la transferencia de energa del
oleaje a distintos fenmenos hidrodinmicos. Este es un fenmeno genrico para las
playas de todo el mundo, que pueden ser consideradas como zonas de alta energa y de
3.3
migracin activa de material sedimentario mucho mas intensa que en la mayora de los
lugares en tierra o a mayor profundidad en el mar.
Los principales factores involucrados en la morfodinmica son, el rgimen
hidrodinmico al que esta sometido la zona costera, las condiciones morfolgicas y
batimtricas, y las caractersticas del material sedimentario del que se compone la playa.
Como se ha mencionado antes, las caractersticas del oleaje son bsicamente las que
determinan el rgimen hidrodinmico que se retroalimentar con las variaciones
morfodinmicas. Por ello es comnmente aceptado que hay perfiles caractersticos para
cada poca del ao, asociados a las variaciones estacionales de la energa del oleaje,
debido a que hay un balance sutil y delicado entre las fuerzas que tienden a traer arena
del mar y aquellas que se lo devuelven.
Uno de los primeros experimentos realizados para estudiar la evolucin de un perfil de
playa en funcin de las caractersticas del oleaje fue el realizado por Johnson (1956),
que encontr que la forma de equilibrio del perfil de playa depende del peralte en aguas
profundas (H0 / L0 ). Bajo trenes de oleaje con altura pequea y periodo largo (H0 / L0 ) <
0.03, la arena procedente del pie del perfil se acumula por encima del nivel medio del
mar en forma de berma, hasta alcanzar una forma de equilibrio con el perfil muy
peraltado que favorece la reflexin del oleaje incidente, cesando la tendencia
acumulativa. Por el contrario, para olas peraltadas, es decir, con mayor altura en
relacin al periodo (H0 / L0 ) > 0.03, el perfil inicial se erosiona transportando mar
adentro la arena de la berma, acumulndola por debajo del nivel medio del agua en
forma de barra sumergida. Esta barra puede provocar que el oleaje rompa y disipe su
energa a mayor distancia de la lnea de orilla, con lo que se reduce en gran medida la
accin erosiva.
3.3.- Escalas de los procesos costeros

Los procesos costeros, hidrodinmicos y sedimentarios, son extremadamente complejos
e incluyen fenmenos a escalas de espacio y tiempo completamente diferentes.
Horikawa (1970) propone una forma prctica de clasificar los fenmenos costeros en
tres categoras: macroescala, mesoescala y microescala (tabla 3.1).
Escala
Tiempo
Espacio
Macroescala
Ao
Kilmetro
Mesoescala Microescala
Da/hora
segundo
Metro
Milmetro
Tabla 3.1.- Clasificacin de las escalas de los procesos costeros.
Tericamente hablando, los fenmenos a mesoescala deben componerse de la

superposicin completa de los fenmenos a microescala, y los fenmenos a macroescala
de la de los fenmenos a mesoescala. Los sistemas morfodinmicos costeros que
componen una playa tambin se pueden clasificar jerrquicamente (Kroon, 1994), cada
uno con su propia escala espacial y temporal (tabla 3.2).
3.4
Sistema
Celda Costera Celdas Costeras Tramo Costero Lnea de Costa
Jerrquico
(unitaria)
(mltiples)
Escala
Pequea Escala Meso-Escala
Gran Escala
Meta-Escala
Tabla 3.2- Clasificacin jerrquica de sistemas morfodinmicos en la costa.
Cada uno de estos sistemas se comporta como una unidad, particular y compleja, donde
intervienen la hidrodinmica (los mecanismos impulsores), los procesos de transporte
de sedimentos y los cambios morfolgicos. Estos factores que controlan la
morfodinmica, relacionan a cada una de estas escalas espacio-temporales y se
manifiestan en forma de flujos, ciclos y transformaciones de materia y energa. Las
adaptaciones mutuas entre estos factores (hidrodinmica, procesos de transporte de
sedimentos, cambios morfolgicos y las condiciones ambientales) es lo que se conoce
como morfodinmica.
En un sistema morfodinmico costero, la hidrodinmica puede inducir un transporte de
sedimentos, lo que, a lo largo del tiempo, dar lugar a un cambio morfolgico en la
costa. Este cambio tendr una fuerte interaccin e influencia sobre la hidrodinmica, lo
que dar como resultado otra respuesta de los procesos de transporte de sedimentos, y
una nueva adaptacin morfolgica a las nuevas condiciones hidrodinmicas (figura 3.1).
Hidrodinmica
Olas y corrientes
Tte. Sedimentos
Morfologa
Volmen de Sedimento
Forma, Erosin, Acrecin
Figura 3.1.- Factores que intervienen en la morfodinmica.
Las relaciones entre la hidrodinmica, el transporte de sedimentos y la morfologa son

no- lineales. Por ejemplo, el transporte de sedimentos se lleva a cabo siempre que se
supere un valor umbral de alguna de las variables hidrodinmicas (como se ver mas
adelante) y adems, el cambio morfolgico no es inmediato, es decir, hay un desfase
temporal entre las condiciones hidrodinmicas y la respuesta morfolgica. Este periodo
de tiempo incluye un tiempo de reaccin y un tiempo de relajacin (De Boer, 1992). La
reaccin es el intervalo de tiempo entre el cambio en el mecanismo impulsor (la
hidrodinmica) y el comienzo del cambio morfolgico, mientras que la relajacin es el
intervalo de tiempo desde el comienzo en el cambio morfolgico hasta el ajuste y
3.5
acoplamiento total (forma final) de la morfologa a la hidrodinmica (en particular, para
una costa arenosa altamente energtica, la reaccin puede ser inmediata, es decir, tan
rpida que puede considerarse despreciable). Los ambientes sedimentarios costeros son
tan dinmicos que las adaptaciones morfodinmicas pueden llevarse a cabo de una
forma inmediata a un cambio en las fuerzas hidrodinmicas. El periodo de relajacin
puede considerarse como una medida de la inercia morfolgica dentro de cada uno de
los sistemas costeros y es funcin de la intensidad de la hidrodinmica y del volumen de
sedimentos involucrado en el cambio morfolgico. De todo esto se deduce que la
morfodinmica est regida por mecanismos de acoplamiento y retroalimentacin de la
hidrodinmica y los cambios morfolgicos. Los mecanismos de retroalimentacin son
determinantes en la evolucin de los sistemas costeros, pudiendo ser positivos (cuando
las inestabilidades del sistema costero tienden a amplificarse) o negativos (cuando el
sistema tiende a oponerse a las inestabilidades, permaneciendo en un estado de
equilibrio dinmico estable). Adems, por lo general, las variaciones hidrodinmicas
ocurren en periodos de tiempo menores que la duracin de los tiempos de reaccin y
relajacin morfolgica, por lo que las condiciones hidrodinmicas habrn cambiado
antes de que la morfologa se haya adaptado del todo a las condiciones iniciales. Los
ambientes hidro-morfodinmicos costeros, tienden a buscar un estado de equilibrio con
un mximo de estabilidad (Wright & Thom, 1977). No obstante, no parece que el
concepto de equilibrio dinmico pueda aplicarse sobre un periodo largo de tiempo,
debido a que los factores que intervienen en la hidro-morfodinmica son
extremadamente variables y complejos, con interacciones no- lineales.
Los sistemas costeros se encuentran cambiando la mayor parte del tiempo, tendiendo a
llegar a un estado quasi-equilibrado entre la hidrodinmica y la morfologa (debido a
que es difcil que se alcance un equilibrio dinmico en todos los componentes hidromorfodinmicos). En el caso de los procesos a escalas pequeas, las relaciones entre los
procesos impulsores y la respuesta morfolgica son muy rpidas, y actan muy
rpidamente sobre periodos de tiempo muy cortos, por lo que la morfodinmica a
escalas pequeas puede alcanzar rpidamente un estado de equilibrio dinmico. En el
caso de los procesos a escalas mayores, las relaciones entre los procesos impulsores y la
respuesta morfolgica ocurren a tasas mas bajas, en ocasiones muy bajas (casi
imperceptibles) y actan sobre periodos de tiempo muy largos.
Los ambientes costeros son muy variados, y cada uno de ellos a su vez es nico en
cuanto a sus caractersticas dinmicas, por lo que las relaciones entre los procesos a
distintas escalas espaciales y temporales son distintas tambin. Por ejemplo, una
tormenta muy severa o una marea meteorolgica de unas cuantas horas puede impactar
de tal manera en la costa que el periodo de relajacin sea de varios aos. Parece claro
que no hay una relacin nica y universal entre los procesos y los estados de equilibrio a
las distintas escalas espaciales (longitudinal y transversal a la costa) y temporales para
todos los ambie ntes costeros, aunque su importancia es evidente (figura 3.2). El estado
actual del conocimiento de estas relaciones est aun lejos de estar completo, ya que
como se ha mencionado antes, los procesos involucrados son numerosos y complejos y
hay aun muchas incertidumbres en torno a ellos.
3.6
Figura 3.2.- Escala temporal de variaciones de la lnea de costa y procesos costeros

(Adaptado de Kroon, 1994).
Los procesos morfodinmicos son especficos de cada lugar (es decir, no hay una
relacin nica entre los mecanismos y sus escalas espaciales y temporales involucradas)
por lo que las caractersticas fsicas de cada lugar son especficas. De aqu que una de
las mayores incertidumbres que conllevan los estudios numricos y experimentales de la
costa para determinar algn patrn de comportamiento o evolucin, es la seleccin de
las escalas espaciales y temporales de los procesos fsicos que ah ocurren.
Generalmente, la seleccin de estas escalas (en particular la escala espacial) se lleva a
cabo de una manera mas bien intuitiva (y por lo tanto arbitraria), pero una vez que se
han seleccionado, estas escalas establecen a su vez la escala temporal de las
predicciones del desarrollo de la costa. Las escalas longitudinales a la costa
generalmente sern mucho mayores que las escalas transversales, del orden de uno a
varias decenas de kilmetros en el primer caso, mientras que las transversales sern de
tan solo unos cuantos centenares de metros, generalmente desde la profundidad activa
hasta la lnea de dunas. Cabe sealar que las variaciones transversales en la zona
cercana a la costa suelen ser mucho m s rpidas que las variaciones longitudinales, por
lo que estas a su vez se suelen subdividir en tramos:
Desde la profundidad activa hasta la zona de rompientes.

La zona de rompientes.
Zona de swash.
Zona de playa seca.
Las posiciones y lmites entre dichas zonas no son fijas, sino mas bien graduales y
pueden cambiar con el tiempo. Esto sugiere que los distintos grupos de investigacin
(especialmente en Geologa e Ingeniera Martima) dediquen sus esfuerzos a clarificar
en paralelo, los mecanismos fsicos que intervienen en los fen menos hidromorfodinmicos a escala distintas, y que se traten de llenar los grandes huecos
existentes entre ellos.
Las caractersticas morfolgicas de una playa pueden clasificarse de acuerdo a las
distintas escalas espaciales y temporales. En la escala espacial m s pequea, se
encuentra la micromorfologa del lecho marino: los ripples generados por olas y
corrientes. La dinmica de los ripples es muy alta, el tiempo de reaccin es despreciable
3.7
y el de relajacin puede ser muy corto, del orden de segundos (asociado al flujo
oscilatorio) o de unas cuantas horas en el caso de corrientes. A mesoescala, algunos de
los rasgos morfolgicos m s caractersticos son, por ejemplo, la topografa rtmica a lo
largo de una playa (beach cusps) o los megaripples (e.g. generados por las corrientes de
marea). De cualquier manera, la dinmica de estos rasgos morfolgicos sigue siendo
muy alta y los tiempos de reaccin despreciables y de relajacin muy cortos, del orden
de minutos a horas (por ejemplo un ciclo de marea). A macroescala, los rasgos
morfolgicos son en su mayora zonas de gran erosin o deposito relacionadas con
grandes eventos de tormenta o terrazas marinas relacionadas con variaciones peridicas
en el nivel del mar. Aqu tambin existe una dinmica muy importante donde los
tiempos de reaccin son despreciables, pero los tiempos de relajacin son ms altos que
en los casos de micro y meso escala, del orden de varias horas hasta algunas semanas.
Las interacciones entre las diferentes escalas de la morfodinmica de la zona costera
pueden ser consideradas como eventos que se retro-alimentan, igual que entre la
morfodinmica y la hidrodinmica. A menudo los rasgos macro-morfolgicos del
lecho marino pueden ser interpretados como la respuesta de los cambios morfo lgicos
sucesivos a micro y meso-escala integrada en el tiempo a la vez que, los rasgos macromorfolgicos pueden determinar la posicin y ocurrencia de los rasgos micro y mesomorfolgicos. De todo esto se intuye que los estados (erosin o acrecin) as como los
tipos de playa dependern de la composicin global de los rasgos morfolgicos en todas
estas escalas. As, cuando las condiciones energticas (hidrodinmicas) sean bajas, es de
esperar que la playa sea reflejante, con una pendiente alta (y posiblemente montona,
sin barras, al menos en la parte activa del perfil). El caso opuesto, donde las condiciones
energticas sean altas, la playa ser disipativa por lo que el perfil tendr una pendiente
baja, y habr una o mas barras sumergidas en la parte activa del perfil. En definitiva,
para poder correlacionar la tasa del transporte de sedimentos con los mecanismos
hidrodinmicos impulsores, es preciso conocer los mecanismos detallados del
movimiento de fluidos y la respuesta del sedimento. Sin embargo, el conocimiento en
los mecanismos de transporte de sedimento es todava muy limitado, ya que an hay
problemas importantes que deben ser resueltos, la mayora de los cuales podran
considerarse como fenmenos a microescala. Algunos de estos problemas (Horikawa,
1981) son:
i.
ii.
iii.
iv.
v.
vi.
vii.
viii.
Caractersticas del flujo oscilatorio en la capa lmite, debido a la coexistencia de

olas y corrientes, tales como la distribucin de velocidades y la tensin de corte
de fondo bajo distintas condiciones hidrodinmicas.
Distribucin vertical del coeficiente de intercambio de momentum (o viscosidad
turbulenta) y del coeficiente de difusin en el campo de las olas.
Distribucin vertical de la velocidad del transporte de masa.
La deformacin del oleaje y el campo de velocidades en un fondo con pendiente.
Interaccin entre olas y corrientes inducidas por el oleaje.
Formacin, tamao y desaparicin de ripples, as como la formacin del llamado
sheet flow.
Distribucin vertical de la concentracin del sedimento en suspensin, dentro y
fuera de la zona de rompientes.
Los mecanismos detallados del transporte de sedimentos debido a las olas y a las
corrientes.
3.8
An cuando se han llevado a cabo numerosas investigaciones en estos temas, muchos
de los detalles se comprenden nicamente en sentido cualitativo. El conocimiento y
descripcin cuantitativa de todos estos procesos es bsico para entender y determinar los
cambios a mesoescala (que son los que mas impacto tienen para el hombre), como las
variaciones en la posicin de la lnea de costa, de la topografa del fondo marino,
formacin de barras, formas rtmicas y corrientes cercanas a la costa.
Debido a que el conocimiento de la dinmica de la zona cercana a la costa ha mejorado
mucho desde la ltima dcada, los movimientos del fluido pueden ser predichos con
bastante precisin bajo las suposiciones adecuadas. Se sabe que las partculas de
sedimento son transportadas por el fluido en movimiento, el flujo oscilatorio asociado al
oleaje y por las corrientes impulsadas por la rotura del oleaje, viento o mareas (los
efectos de estas dos ltimas quedan fuera de la presente discusin). Por lo tanto, la
transformacin del oleaje durante su propagacin hacia aguas someras (refraccin,
shoaling, friccin con el fondo, rotura, etc.) deber determinarse de antemano. El
clculo de dicha transformacin del oleaje, bajo las simplificaciones adecuadas, puede
realizarse exitosamente con los conocimientos hoy disponibles. El siguiente paso es una
prediccin precisa de las corrientes cercanas a la costa, lo cual tambin es actualmente
posible. Sin embargo, debido a la complejidad de los mecanismos de transporte, sigue
habiendo muchas incertidumbres en el campo de investigacin de los procesos costeros
sedimentarios y el tema seguir abierto durante bastante tiempo an. De aqu que, por
ejemplo, an no haya una expresin precisa para evaluar la tasa del transporte de
sedimentos en la costa. Por esta razn, aun cuando se pueden combinar las ecuaciones
fundamentales del movimiento de los fluidos y la conservacin de la masa sedimentaria
para calcular los cambios del fondo en el rea cercana a la costa, hasta ahora los
resultados tienen una precisin limitada.
Para poder superar este problema en aplicaciones prcticas, se han propuesto
numerosas frmulas de transporte longitudinal, basadas en gran medida en resultados
empricos. Combinando una de estas frmulas y la ecuacin de continuidad (o
conservacin del sedimento) obtenida bajo condiciones simplificadas, puede efectuarse
una prediccin de los cambios en la posicin de la lnea de costa. Esta clase de
aproximacin es muy til para fines prcticos, pero no siempre las suposiciones hechas
en los anlisis antes mencionados son correctas, por lo que en ocasiones las
predicciones estn rodeadas de incertidumbres no despreciables, mismas que sern
analizadas y discutidas mas adelante.
Para poder simular los cambios en la configuracin del fondo marino, se parte de la
siguiente ecuacin genrica:
(1 p )
h q x q y
t x y = 0

[3.1]
donde x e y denotan los ejes del plano horizontal, p es la porosidad del sedimento, h es
la profundidad en un punto determinado (x, y) en un tiempo determinado (t) y qx y qy las
tasas de transporte de sedimentos a travs de una longitud unitaria en cualquiera de las
direcciones (x o y respectivamente por unidad de tiempo). Para poder aplicar esta
3.9
ecuacin, qx y qy deben evaluarse considerando las caractersticas del oleaje (altura,
perodo y ngulo), profundidad del agua h y las caractersticas del sedimento (tamao
d50 , densidad, s y porosidad, p). Sin embargo, las ecuaciones disponibles para estimar
las tasas de transporte de sedimentos no son suficientemente precisas para evaluar qx y
qy.
En lo referente al estudio de los fenmenos a macroescala, el punto de vista geolgico
es de mucha ayuda para entender la tendencia general de los procesos costeros. Sin
embargo, la escala de tiempo de inters geolgico es por lo regular muy larga para los
propsitos ingenieriles; de aqu que los grupos de investigacin en Ingeniera Costera
hayan desarrollado sus propias mediciones y dispositivos para obtener datos de los
procesos costeros a un plazo relativamente corto. Este corto plazo, a grandes rasgos,
el plazo de la vida humana, es an bastante largo cuando se considera la prctica normal
de la ingeniera. No obstante, para realizar una correcta explotacin de los recursos que
representa una playa, deberan ser tomados en cuenta los cambios histricos del nivel
del mar, las variaciones en el nivel del suelo producto de los basculamientos del
continente, as como los cambios debidos a perturbaciones ocasionadas por la actividad
humana (e.g. construccin de estructuras costeras que son una de las causas principales
de la erosin de playas en costas de todo el mundo).
Por lo tanto, la combinacin de causas artificiales y naturales en una zona determinada
pueden provocar cambios, en ocasiones irreversibles, sobre la lnea de la costa. Pero en
el futuro, los efectos antropognicos pueden incrementar constantemente su
contribucin al problema de erosin de playas. Por lo tanto, los cambios costeros a corto
y largo plazo deben ser predecibles con la precisin suficiente para cualquier actuacin
en la costa.
En la figura 3.3 se muestran las relaciones entre varios de los elementos discutidos
anteriormente, tales como olas, corrientes cercanas a la costa, transporte de sedimentos,
topografa, etc. con el fin de clarificar estas complejas relaciones. Dependiendo del
problema especfico que se quiera estudiar, es necesario establecer adecuadamente las
escalas discutidas anteriormente.
Figura 3.3.- Relacin entre olas, corrientes cercanas a la costa, transporte de sedimentos y topografa
(Adaptado de Horikawa, Harikai & Kraus 1979).
3.10
3.4.- Mecanismos de transporte de sedimentos

La razn principal de que la fsica del transporte de sedimentos no est del todo
aclarada se debe a la casi insuperable dificultad de medir la tasa y la direccin del
transporte de sedimentos en la zona costera, conjuntamente con la medicin de olas,
corrientes, batimetra y la interaccin entre ellos. Por lo tanto, en esta seccin se
discutirn brevemente algunos temas relacionados con los mecanismos del transporte de
sedimentos en la costa desde la perspectiva de la mecnica de fluidos.
En el anlisis del transporte de sedimentos en zonas costeras, por simplicidad, es muy
comn considerar por separado el transporte de perpendicular y el paralelo a la lnea de
la costa. Es obvio que ambos estn muy relacionados entre si. No obstante, el transporte
transversal es quizs ms importante para los procesos costeros a corto plazo, mientras
que el transporte longitudinal es ms importante para la variacin a medio y largo plazo
de la topografa costera.
En el mbito de la Ingeniera de Costas, el problema del transporte de sedimentos y
evolucin costera tradicionalmente se ha abordado desde una perspectiva individual
(por separado de la hidrodinmica) y en muchos casos, emprica. Sin embargo, est
claro que no se debe seguir tratando ambos aspectos por separado. La hidrodinmica es
el mecanismo impulsor del transporte de sedimentos. Los granos del sedimento en el
lecho se movern cuando las tensiones de corte en el fondo ejercidas por el flujo
excedan un cierto valor crtico. Estas tensiones de corte son provocadas por la friccin
entre el flujo (no idealizado) y el fondo, creando una capa lmite donde habr una
disipacin de energa del oleaje y donde se producir la interaccin con el fondo y el
inicio del movimiento del sedimento. Esta capa lmite se puede definir intuitivamente
como la capa dentro de la cual el flujo est totalmente influenciado por la presencia del
fondo. Debido a la naturaleza oscilatoria del flujo asociado al oleaje (hay un cambio de
direccin cada periodo), la capa lmite del oleaje es una capa muy delgada, confinada a
las cercanas del fondo, mientras que la capa lmite de una corriente puede crecer hasta
alcanzar la regin donde el fluido se mueve, libre de la presencia del fondo. La altura de
la capa lmite del oleaje es por lo tanto mucho menor, lo que provoca que las tensiones
de corte debido a la accin del oleaje sean mucho mayores que las debidas a una
corriente, aun cuando las velocidades sean de magnitud comparable. Por lo tanto, el
inicio del movimiento del sedimento en los ambientes costeros est mucho mas
influenciado por el flujo oscilatorio asociado al oleaje.
La propia movilidad del sedimento puede afectar al desarrollo de la misma capa lmite.
Partiendo de un hipottico fondo plano, la rugosidad efectiva ser funcin nicamente
del tamao de grano (d50 ), sin embargo, una vez que se ha superado un valor crtico de
las tensiones de corte y se inicia el movimiento, se crean ripples en el lecho, que a su
vez aumentarn la rugosidad efectiva del fondo afectando al flujo. En este caso, la
magnitud de la rugosidad efectiva depender del tamao (alt ura) del ripple. De hecho,
mediciones en campo de la atenuacin de la energa del oleaje han arrojado resultados
de los valores del factor de friccin del fondo de hasta un orden de magnitud mayor de
lo esperado (Iwagaki & Kakinuma, 1963 y 1967; Treolar & Abernathy, 1978). Esto es
debido a la presencia de ripples y otras formas en el lecho, que traen como resultado
3.11
una resistencia mucho mayor al movimiento del oleaje que en lechos planos y son
consecuencia del mismo oleaje. Este incremento en la rugosidad efectiva provocar
mayores tensiones de fondo y ms disipacin de energa dentro de la capa lmite. Si
aumenta la intensidad del oleaje, entonces se forman ripples bidimensionales muy
regulares con alturas, longitudes de onda y peraltes muy bien definidos. Conforme se
incrementa la amplitud orbital del flujo oscilatorio en el fondo (que es la escala de
longitud del oleaje), tambin lo hace la altura y longitud del ripple, mientras que el
peralte permanece aproximadamente constante. Un incremento importante en la accin
del oleaje podra traer como resultado una disminucin en el peralte de los ripples, que
bajo estas condiciones dejan de ser regulares y comienzan a mostrar caractersticas
tridimensionales y una cresta ms redondeada (con respecto a su estado bidimensional),
debido a que el flujo oscilatorio del oleaje la sobrepasa y la suaviza. Por lo tanto, se
puede intuir que la geometra exacta de los ripples depender de las propiedades del
oleaje, del fluido y de los sedimentos. Si la actividad del oleaje se incrementa, los
ripples tienden a desaparecer, formndose una mezcla de agua y sedimento que se
mueve en un vaivn constante con el paso del oleaje. Este fenmeno se conoce como
sheet flow. La disipacin de la energa de la ola traer como consecuenc ia una
disminucin en su altura, por lo que para predecir la altura del oleaje en la zona costera,
entre otras cosas, es importante cuantificar esta prdida de energa, que est relacionada
con las tensiones de corte en el fondo. Otro punto importante relacionado con la capa
lmite del oleaje es que provoca que las corrientes tengan mas resistencia con el fondo
(Madsen 1994). En lo referente a los mecanismos de transporte de sedimento, las
tensiones de corte en el fondo provocadas por el flujo oscilatorio del oleaje son mayores
que las de las corrientes. Por lo tanto, el flujo oscilatorio del oleaje es capaz de
movilizar el sedimento por si solo, de tal manera que en presencia de una corriente
costera, se puede generar un transporte de sedimento importante (que puede ser por
arrastre y/o en suspensin). El transporte por arrastre se lleva a cabo en la regin
cercana al fondo, donde existe la mayor concentracin de sedimento, y es debido
principalmente a la dislocacin de los granos producida por las tensiones de corte, que
ruedan y saltan sobre el fondo. Por su parte, el transporte de sedimento en suspensin se
lleva a cabo en la columna de agua, donde la cantidad de sedimento que entra en la
columna depende de las tensiones de corte y la turbulencia del flujo (que es el principal
mecanismo de suspensin), mientras que la velocidad del flujo regular la tasa de
transporte. Al mismo tiempo, la velocidad del fluido y las tensiones de corte dependen a
su vez de las formas de lecho, lo que demuestra la complicada naturaleza de la
interaccin entre el flujo y un fondo mvil. El levantamiento del sedimento del fondo
puede provocar la generacin de una capa estratificada de agua y sedimento, un
gradiente importante de concentracin de sedimentos y el amortiguamiento de la
turbulencia dentro de la capa lmite (Soulsby, 1997). Cabe mencionar que si el
movimiento del sedimento se debiera a un oleaje lineal, no debera resultar un transporte
neto de sedimento, aunque en la realidad, cerca de la costa el oleaje es no lineal y una
parte importante del transporte transversal a la costa se debe precisamente a la no
linealidad (asimetra) del oleaje.
De aqu se desprende la importancia de las tensiones de corte en el fondo para la
cuantificacin de la mayor parte de los procesos fsicos que ocurren en la zona cercana a
la costa. Una forma de hacerlo es relacionar las tensiones de corte en el fondo con las
velocidades del flujo oscilatorio en las cercanas del lecho usando un factor de friccin
similar al usado en capas lmite de flujos estacionarios (Kajiura, 1964; Jonsson, 1966).
Este factor de friccin puede ser determinado con la rugosidad relativa del fondo, que es
3.12
la relacin entre la escala de rugosidad k n y la amplitud A de la semi-excursin de la
partcula del agua (que es la escala de longitud caracterstica del oleaje). k n se puede
determinar experimentalmente a partir de un ajuste logartmico de las velocidades
medidas en la columna de agua y cerca del fondo como k n = 30z0 (para flujos
turbulentos), donde z0 es la ordenada en el origen del ajuste logartmico. Si el flujo es
totalmente turbulento, como es en la mayora de los casos en la naturaleza, el factor de
friccin en el fondo depender nicamente de la rugosidad relativa del fondo.
Una suposicin fundamental para esta aproximacin (igual que como se hace para
capas lmite estacionarias turbulentas), es que la geometra del fondo puede ser
representada por una nica escala de longitud (Madsen, 1993). Esto, aunque es una
simplificacin muy grosera para un fondo con una geometra muy compleja, es una
aproximacin suficiente, ya que se ha visto que la dispersin de los datos
experimentales es tal que no se justifica la utilizacin de esquemas ms complejos para
describir la geometra del fondo. Por lo tanto, cuando se conoce con antelacin la
geometra del lecho, se puede estimar la longitud de la rugosidad del fondo y para un
oleaje dado, se puede calcular las tensiones de corte del oleaje.
Una vez se han generado, la rugosidad del fondo pasar a depender de las formas
geomtricas en el fondo (ya no del tamao de grano). Bajo las condiciones de sheet
flow, la rugosidad equivalente en el fondo ser funcin (estar escalada) del grosor de la
capa mezcla de agua y sedimentos cerca del fondo, por lo que habr un incremento en la
friccin entre el flujo y el fondo. Por lo tanto, la rugosidad del lecho para fondos
arenosos no es constante y depende del oleaje, el fluido y los parmetros del sedimento,
adems de que la presencia de ripples puede cambiar la resistencia al flujo, tanto del
oleaje, las corrientes como la accin combinada de ambos.
En lo referente al transporte de sedimento, este se puede llevar a cabo por la accin de
las corrientes costeras (de marea, viento o inducidas por la rotura del oleaje) o por la
accin del oleaje, aunque lo mas comn es que se vea transportado por la accin
conjunta de las corrientes y del oleaje. El primer mecanismo de transporte es la entrada
del sedimento del fondo en la columna de agua, como resultado de la friccin entre el
fondo marino y el flujo (de las corrientes y/o el oleaje) junto con la difusin turbulenta
que tiende a mantener los granos de arena en suspensin. Una vez movilizados (como
respuesta a la friccin con el flujo y la turbulencia), los sedimentos se pueden
transportar, dependiendo de las caractersticas del flujo y de las suyas propias, por
arrastre, saltacin o suspensin. El transporte por arrastre es el modo de transporte
dominante para flujos lentos y/o granos de arena gruesos. En caso de que el flujo sea
mas rpido (oleaje alto o corrientes fuertes) el sedimento ser puesto en suspensin en la
columna de agua y ser susceptible de ser transportado grandes distancias. De aqu se
define la tasa de transporte de sedimentos como la cantidad de sedimento (en masa o
volumen) que se mueve a travs de un plano vertical con una anchura unitaria que es
perpendicular a la direccin del flujo. As, la tasa de erosin / deposito de sedimento en
la zona costera depender de la diferencia entre las tasas de transporte de sedimento que
entra y que sale de una zona en particular.
3.13
3.5.- Interaccin Fluido-Sedimento

Los efectos de los mecanismos hidrodinmicos impulsores la sobre la dinmica
sedimentaria tienen lugar principalmente en las proximidades del fondo debido a la
friccin que se genera entre el lecho marino y el flujo. Esto es lo que se conoce como
tensi n de corte en el fondo, que es la fuerza que ejerce el flujo por unidad de rea del
fondo. De aqu la importancia de la revisin de los procesos que nacen de la interaccin
del fondo con el flujo, que son determinantes en la dinmica costera. Hay varias
maneras de definir la tensin de corte en el fondo b , que dimensionalmente es una
fuerza por unidad de rea [N/m2 ] y que tambin puede expresarse en trminos de
velocidad [m/s], ya que la velocidad de corte u* , se define como b = u*2 , donde es
la densidad del fluido. De aqu se define la velocidad de corte como u* = b . Esta
definicin debe tomarse con cierta cautela, ya que, aunque puede relacionarse con las
fluctuaciones turbulentas de las componentes de la velocidad real del flujo, la
velocidad de corte no representa una velocidad real del flujo.
En ocasiones es necesario utilizar las tensiones de corte de una forma adimensional
para lo cual se las relaciona con el sedimento, dando lugar al parmetro adimensional de
Shields, definido en trminos de las tensiones de corte en el fondo o de la velocidad de
corte respectivamente (ecuacin [3.142]), involucrando la aceleracin de la gravedad g,
la densidad del agua (), la densidad del sedimento (s) (o la relacin entre ambas,
definida como s = s / ) y el dimetro del sedimento d.
Las tensiones de corte generadas no solo dependen de la velocidad del flujo, tambin
dependen de la rugosidad del fondo. Esta rugosidad se puede medir, ya sea por la
rugosidad equivalente de Nikuradse, k n (que esta relacionada con el tamao del
sedimento), o por la escala de rugosidad z0 (que se puede derivar del perfil de
velocidades) y que se pueden relacionar a travs de la expresin z 0 = kn /30 , por lo que
k n = 30 z 0 . Esta notacin generalmente se emplea para definir valores no estacionarios
(variables en el tiempo), por lo que es muy til distinguir los diferentes mecanismos por
los cuales se pueden generar las tensiones de corte: olas, corrientes y/o la accin
combinada de ambas. Los valores del oleaje estn relacionados con las valores de los
parmetros del flujo oscilatorio, las corrientes involucran valores medios, y en el caso
del flujo combinado de olas y corrientes, se pueden obtener valores promediados y/o los
valores mximos en el periodo del oleaje. Aunque las fluctuaciones turbulentas pueden
provocar grandes variaciones de los valores obtenidos en cualquiera de los flujos
mencionados, para aplicaciones ingenieriles, generalmente es suficiente con los clculos
promediados.
Puesto que la naturaleza de los flujos debidos a corrientes p uras y oleaje puro son
totalmente distintos, los valores de z0 y k n pueden ser totalmente distintos entre ambos
casos, sobre todo si la orientacin del oleaje, las corrientes y la de los ripples del fondo
son distintos (Mathisen & Madsen, 1996). As, el valor de las tensiones de corte en el
fondo 0 es producto de 3 contribuciones bsicas: la friccin generada por y sobre el
sedimento 0s, el arrastre producido por los campos de presin asociados al flujo sobre
los ripples 0f y el transporte de sedimento 0t provocado por la transferencia de
3.14
momentum al movilizar al sedimento. Esto es igualmente aplicable al concepto de la
velocidad de corte u* , el parmetro adimensional de Shields y la escala de rugosidad
del fondo z0 .
Es importante recalcar que la friccin es la contribucin relacionada con el sedimento
(Van Rijn, 1993) y la manera usual de calcular esta contribucin es por medio del perfil
logartmico de velocidades, donde el valor de z0 (o de k n ) esta relacionado slo con el
tamao de grano. El valor resultante de las tensiones de corte en el fondo no es
realmente el factor de friccin que estara medido en la superficie de un fondo con
ripples, aunque es similar, y es una medida convencional (k n = 2.5 d50 ) y prctica contra
la cual la respuesta del sedimento puede correlacionarse. En el caso de que se trate de
un fondo razonablemente plano (el caso mas comn para fondos con sedimentos
gruesos, 0.8 mm > d50 ), las tensiones de corte debidas a la friccin sern las tensiones de
corte totales (es decir, solo tendrn como componente la friccin con el sedimento). En
caso contrario, cuando el sedimento es mas fino, generalmente se encuentra un fondo
con ripples y/o puede presentar un flujo en capa mas intenso, con lo cual las tensiones
de corte sern las debidas a la friccin con el sedimento, mas las debidas a las formas
del fondo y/o las debidas a la transferencia de momentum del flujo al sedimento
movilizado.
Ya que solo la friccin acta directamente sobre el grano, es importante establecer su
magnitud para evaluar el inicio del movimiento, transporte y tasa de captacin de
sedimento que se pone en suspensin. Por otra parte, la determinacin de las tensiones
de corte en el fondo sirve para calcular la resistencia al flujo con el fondo (que
determina las intensidades turbulentas que influyen directamente en la difusin del
sedimento en suspensin en la columna de agua).
En definitiva, la dinmica de la zona cercana a la costa involucra la complicada
interaccin de un gran nmero de fenmenos fsicos, muchos de ellos comprendidos
slo desde una perspectiva cualitativa, por lo cual su estudio presenta muchas
incertidumbres e imprecisiones. Las razones de dichas imprecisiones son muy variadas
y se pueden agrupar genricamente en las siguientes categoras:
Los mtodos predictivos no incluyen todos los procesos o mecanismos

implicados.
Algunos de los mecanismos no estn bien comprendidos y/o las formulaciones
de dichos mecanismos son imprecisas, muy simplificadas o incompletas.
Uso de formulaciones (particularmente las empricas) fuera de su propio rango
de validez.
Uso de datos de campo de una localidad distinta.
As, cuando se hacen predicciones de la dinmica de la zona cercana a la costa, es

importante llevar a cabo un anlisis de sensibilidad para evaluar los errores en las
estimaciones, que nacen a partir de las incertidumbres de las simplificaciones asumidas
y de los parmetros y datos de entrada, as como un anlisis comparativo de las distintas
aproximaciones utilizadas (comparaciones entre distintos modelos predictivos).
Algunas de las incertidumbres tpicas en los parmetros de entrada se deben a la
combinacin de los errores en la instrumentacin de medida y en la representatividad
3.15
de un solo valor (e.g. en la altura del oleaje) sobre una zona muy grande o sobre
periodos de tiempo largos (tabla 3.3).
Parmetro
Densidad del agua
Viscosidad cinemtica
Densidad del sedimento
Tamao del sedimento
Profundidad (batimetra)
Velocidad de la corriente
Direccin de la corriente
Altura del oleaje
Periodo del oleaje
Direccin de aproximacin del oleaje
Unidades % Error
[kg/m3 ]
0.2
2
[m /s]
10
3
[kg/m ]
2
[m]
20
[m]
5
[m/s]
10
grados
10
[m]
10
[s]
10
grados
15
Tabla 3.3.- Incertidumbres asociadas a algunos de los parmetros fsicos involucrados en la dinmica
costera (Soulsby, 1993).
Las incertidumbres relacionadas con la densidad y viscosidad del agua provienen de las
variaciones espaciales y estacinales de la temperatura y salinidad del agua. Las del
sedimento provienen de las variaciones espaciales y temporales de la composicin y
distribucin de los granos de arena en una playa. Las debidas a la profundidad se deben
a las variaciones en el nivel medio del mar, ya sea por la marea astronmica o
meteorolgica, o bien por la propia imprecisin del aparato de medida. En lo que se
refiere al oleaje, las incertidumbres que hay provienen de las propias medidas (precisin
del aparato), de las tcnicas de anlisis y de las variaciones interanuales del clima de
oleaje. De todo esto se intuye que las diferencias entre distintos mtodos o modelos
predictivos son relativamente pequeas cuando se trata de predecir mecanismos cuya
fsica se conoce o se entiende bien, y por otra parte, las diferencias son relativamente
grandes cuando se trata de mecanismos cuya fsica no se conoce del todo. As por
ejemplo, para la densidad y la viscosidad del agua, el error es del orden del 1%; para la
velocidad orbital y velocidad de cada del sedimento es del orden del 10%; para el
umbral de las tensiones de corte, velocidad de la corriente o altura de ola es del orden
del 20%; para las tensiones de corte media y mxima bajo el flujo combinado de olas y
corrientes es del orden del 50% y para las tasas de transporte de sedimento en
suspensin para corrientes y para corrientes y oleaje combinados es de un factor de 3 a 5
respectivamente.
Otra fuente de incertidumbres asociada a la prediccin de la dinmica de la zona
cercana a la costa proviene de la misma seleccin de los parmetros representativos
utilizados para describir el oleaje y los procesos de rotura. Las diferencias que se
pueden encontrar al utilizar distintos parmetros pueden dar lugar a variaciones de ms
del 50% en las predicciones, aun cuando el mtodo utilizado sea el mismo. Por otra
parte, la variacin de algunas de las propiedades del sedimento y del agua del mar puede
presentar otra fuente de incertidumbre. Se sabe, por ejemp lo, que la densidad del agua
de mar generalmente decrece con la temperatura y aumenta con la salinidad. As la
mxima densidad del agua dulce se alcanza a los 4 C mientras que la del agua de mar
3.16
es a -1.9 C. En el proceso de movilizacin de sedimentos, cerca del fondo se puede
formar una capa estratificada de aguasedimento con una densidad diferente, lo que
puede provocar un amortiguamiento de la turbulencia dentro de la capa lmite (Soulsby,
1997) e incluso provocar una corriente de turbidez si hay una pendiente en el fondo. La
viscosidad cinemtica es una propiedad molecular del agua, cuyo valor decrece con la
temperatura y se incrementa con la salinidad. Esta propiedad del agua de mar est
fuertemente relacionada con las tensiones de corte.
3.5.1.- Relacin entre la hidrodinmica y los mecanismos de transporte

Teniendo en cuenta que los procesos de transporte de sedimento se han estudiado
durante muchos aos y que an no es posible predecir las tasas de transporte sin un
grado muy alto de incertidumbre, se puede intuir la complejidad de la naturaleza de
estos procesos. La gran mayora de los estudios de laboratorio de transporte de
sedimento se han llevado a cabo bajo condiciones del flujo y lecho totalmente
idealizadas. En la realidad, la interaccin del flujo con el sedimento (y su transporte) es
un problema no lineal, complicado por las formas geomtricas en el lecho, las
caractersticas del material del fondo, la variabilidad de las corrientes, la superposicin
del oleaje, el desarrollo de flujos turbulentos en la capa lmite del fondo y por la rotura
del oleaje. El problema por lo tanto, es m s complejo que un caso idealizado, y se
complica an mas teniendo en cuenta que el transporte de sedimentos se lleva a cabo
tanto por arrastre de fondo como en suspensin, y la separacin cuantitativa de estos dos
mecanismos de transporte es tambin, altamente compleja.
La forma mas apropiada de abordar el problema de la prediccin del transporte de
sedimentos y la evolucin del fondo, es relacionando la friccin ejercida por el fluido
(debida a la viscosidad del fluido) y la generacin de una capa lmite (con la
consecuente tensin de corte en el fondo b), con la tasa de transporte de sedimento q.
Para conseguir esto, hay bsicamente dos caminos. Por una parte, las aproximaciones
empricas que usan medidas experimentales de los perfiles de velocidades del fluido
(sus variaciones en la profundidad), as como la rugosidad del fondo y el tamao de
grano para determinar las relaciones de proporcionalidad que pueda haber entre los
distintos parmetros sobre una gran variedad de distintas condiciones ambientales. Por
otra parte, estn las aproximaciones tericas, donde normalmente se intenta usar la
dinmica turbulenta del fluido para determinar los valores de proporcionalidad
directamente. No obstante y a pesar de su gran utilidad, el anlisis estocstico por si
solo no proporciona una teora predictiva de flujo turbulento, y por lo tanto, tampoco de
la mecnica de transporte de sedimento en flujos turbulentos. De aqu la necesidad de
llevar a cabo experimentos de validacin para determinar los coeficientes (por ejemplo,
de difusin) que nacen de las distintas suposiciones incorporadas en cualquier
formulacin terica predictiva sobre el transporte de sedimentos. Una de las mayores
dificultades que trae consigo el estudio del transporte de sedimento (tanto emprica
como analticamente) es la dificultad de medir adecuadamente las tasas de transporte y
las condiciones hidrodinmicas cerca del fondo. Consecuentemente las formulaciones
existentes para predecir las tasas de transporte de sedimento muestran grandes
discrepancias, an cuando se apliquen a los mismos datos.
3.17
Debido a la naturaleza oscilatoria de las velocidades orbitales del oleaje, la capa lmite
tiene poco tiempo para desarrollarse (concretamente, la mitad del periodo del oleaje).
Esto trae como resultado el desarrollo de una capa muy delgada justo encima del fondo,
la capa lmite del oleaje, dentro de la cual, la velocidad del fluido cambia desde su valor
libre (a una cierta distancia z sobre el fondo) a cero (en el fondo mismo), donde hay
condiciones de friccin. La gran velocidad de corte dentro de la capa lmite del oleaje,
que es del orden de unos cuantos centmetros de espesor, produce grandes niveles de
turbulencia y grandes tensiones de corte en el fondo. En contraste con el movimiento
inducido por el oleaje, la variacin temporal de una corriente ocurre sobre una escala de
tiempo mucho mayor, hasta del orden de unas cuantas horas. Por lo tanto, an cuando la
corriente no sea del todo estacionaria, la capa lmite generada por dicha corriente tendr
una mayor escala vertical, del orden de varios metros. En consecuencia, la velocidad de
corte, las intensidades turbulentas y las tensiones de corte generadas en el fondo sern
mucho menores para una corriente que para el flujo asociado al oleaje, aun cuando
tengan velocidades similares.
3.5.2.- Naturaleza de los flujos en la zona cercana a la costa

En la zona cercana a la costa, el flujo cercano al fondo estar determinado por la
interaccin no lineal de las olas y las corrientes. La superposicin de flujos de diferentes
escalas de tiempo y por lo tanto, de diferentes escalas espaciales de la capa lmite traen
como resultado que la capa lmite del oleaje, considerablemente mas delgada, quede
anidada dentro de la capa lmite de las corrientes. El aumento de la turbulencia generada
dentro de la capa lmite debido al oleaje trae como consecuencia que las corrientes
experimenten una mayor resistencia con el fondo en presencia del oleaje que en
ausencia de ste. Por su parte, el flujo cercano al fondo, dentro de la capa lmite del
oleaje, tambin se ver afectado por la presencia de las corrientes, pero en mucho menor
medida que lo que la corriente se ve afectada por la presencia del oleaje.
Debido a la naturaleza oscilatoria del flujo asociado al oleaje, la capa lmite del fondo
tiene un periodo de tiempo limitado para desarrollarse, aproximadamente la mitad del
periodo del oleaje. Esto trae como resultado el desarrollo de una capa lmite delgada, de
unos cuantos centmetros de grosor, inmediatamente por encima del fondo (la llamada
capa lmite del oleaje), dentro de la cual las velocidades del fluido cambian desde su
valor en la corriente libre a cero en el fondo de la capa lmite. La gran velocidad de
corte dentro de la capa lmite del oleaje produce grandes niveles de turbulencia de gran
intensidad y una gran tensin de corte en el fondo.
En contraste al movimiento del oleaje, una corriente (inducida por viento, marea, etc.)
ser mucho mas estacionaria, con lo cual su periodo de variacin es mucho mayor,
pudiendo ser hasta del orden de horas. Por lo tanto, an si la corriente es fuerte y vara
poco en el tiempo, la capa lmite del fondo asociada a la corriente tend r una escala
longitudinal mucho mayor que la capa lmite del fondo asociada al oleaje, hasta del
orden de varios metros. As, como se ha mencionado antes, la velocidad de corte, las
intensidades turbulentas y las tensiones de corte en el fondo sern mucho menores para
una corriente que para el flujo asociado al oleaje, an cuando tengan una velocidad
comparable.
3.18
De acuerdo con lo anterior, parece claro que los procesos que rigen la fsica de la hidromorfodinmica cerca del fondo en la zona cercana a la costa (la interaccin fluidosedimento), estarn dominados por el oleaje. Por lo tanto, el oleaje es el responsable de
iniciar el movimiento del sedimento (tanto por las tensiones de corte en el fondo, como
por la turbulencia asociada a su rotura). Aunque las corrientes no tienen un papel tan
importante en el inicio del movimiento del sedimento, si que lo tienen como mecanismo
de transporte, ya que en principio y debido a la naturaleza oscilatoria del oleaje, el
transporte neto de sedimento por accin del oleaje es casi nulo. Adems de ser el
mecanismo responsable del inicio del movimiento del sedimento para que las corrientes
(que por si mismas no tienen gran capacidad de iniciar el movimiento) sean capaces de
transportarlo, el movimiento del oleaje provoca que el fondo cambie su forma,
generando ripples. La aparicin de los ripples en un fondo arenoso incrementa la
rugosidad fsica de la interfase fluido-sedimento, que para un fondo plano viene dada
por el tamao del grano, y pasa a depender de la geometra del ripple. Es decir, al
rugosidad aumenta en algunos rdenes de magnitud. Por lo tanto, el oleaje no slo
afecta a la resistencia de las corrientes con el fondo al aumentar las intensidades
turbulentas dentro de la capa lmite del oleaje, si no que adems, al propiciar la
generacin de ripples, incrementa la rugosidad fsica del fondo.
Una vez que se han alcanzado las condiciones para que el sedimento se ponga en
movimiento, habr un movimiento oscilatorio del sedimento acorde al oleaje con una
componente neta de transporte asociada a la presencia de una corriente o de un fondo
con pendiente. Para unas condiciones de flujo justo por encima del umbral de
movimiento, el movimiento del sedimento se ve confinado a arrastre y un poco de
saltacin. A este modo de transporte se le conoce como transporte por arrastre de fondo.
Cuando el transporte por arrastre de fondo se hace m s intenso, los granos que se
encuentran en saltacin (y cuya tendencia a volver al fondo esta asociada con su
velocidad de cada) pueden escapar lo suficientemente lejos del fondo como para ser
capturados por los vrtices turbulentos que transportan al sedimento hacia la columna
de agua. Balanceando el transporte hacia el fondo asociado con la velocidad de cada
del sedimento y el transporte hacia la columna de agua provocado por la mezcla
turbulenta, se produce un perfil de concentracin de sedimentos en la columna del agua.
Una vez el sedimento ha sido puesto en suspensin, puede ser advectado por las
corrientes en lo que se conoce como transporte de sedimentos en suspensin.
En buena parte de las aproximaciones para el estudio de estos fenmenos fsicos se
utiliza un modelo simple de viscosidad turbulenta propuesta por Grant & Madsen
(1986) con el cual se obtiene una expresin de cierre analtica para flujos en la capa
lmite del fondo bajo la accin combinada de olas y corrientes. Otras teoras alternativas
de interaccin olas-corrientes que se aplican al transporte de sedimentos se pueden
encontrar en Bijker (1967), Fredse (1984), Madsen (1993) y Davies et al. (1988).
Para describir la formacin de una capa lmite se considerarn slo sedimentos no
cohesivos, teora lineal del oleaje antes de la rotura y corrientes unidireccionales y
estacionarias. Esto implica que esta descripcin dentro de la zona de rompientes es
discutible, y no se puede aplicar si las corrientes presentan cambios apreciables en el
vector velocidad en la columna del agua.
3.19
3.5.3.- Definicin de tensiones de corte y capa lmite
La mayora de los flujos que transportan sedimentos presentan una capa lmite con flujo
de corte turbulento, donde las fuerzas que el fluido ejerce sobre el sedimento son
debidas bsicamente a la turbulencia. En una situacin simplificada bidimensional, estos
flujos se caracterizan por tener mayor variacin horizontal que vertical. Esta diferencia
en las escalas de longitud permite despreciar la aceleracin vertical dentro de la capa
lmite (Schlichting, 1960). El concepto de capa lmite, que es el vnculo entre el flujo de
un fluido ideal (no viscoso) a un fluido real (viscoso) fue ideado en 1925 por Ludwig
Prandtl, que propone que para fluidos con poca viscosidad, el efecto de la friccin
interna en el fluido es apreciable solo en una delgada capa en las vecindades de los
contornos del flujo. Sin embargo, es posible que la capa lmite pueda desarrollarse en la
totalidad de la columna de agua en algunos casos, como en corrientes estacionarias bien
desarrolladas, donde las propiedades del flujo (e.g. densidad, velocidad, presin) en
cualquier punto permanecen constantes en el tiempo. En la figura 3.4 se muestra el
desarrollo de una capa lmite, cuyo grosor () se define como la distancia desde el
contorno hasta un punto tal en la columna de agua donde u = 0.995 U (donde U es la
velocidad del flujo, libre de la influencia del fondo).
Figura 3.4.- Desarrollo de una capa lmite (Liu, 2001).
El desarrollo de la capa lmite puede expresarse como:
flujo laminar
Ux
= 5
x

0.5
dondeRe ( x) =
0.2
flujo turbulento
Ux
= 0.4
x

Ux
< 5 x105
Ux
dondeRe ( x ) =
> 5 x105
[3.2]
3.20
Prandtl (1925) propone el concepto de longitud de mezcla (teora de mezcla de Prandtl)
para poder calcular las tensiones de corte turbulentas en un flujo. En la figura 3.5 se
muestra la estructura de la capa lmite turbulenta y el perfil medio de velocidades para
un flujo horizontal bidimensional fluyendo en el plano x - z, siendo x el plano horizontal
y z el plano vertical.
Figura 3.5.- Estructura de la capa lmite turbulenta y del perfil medio de velocidades (CEM, 1998).
La capa lmite turbulenta se compone de tres subcapas; la subcapa viscosa, la capa de

generacin de turbulencia y la capa exterior. Hidrodinmicamente, la capa de
generacin de turbulencia se caracteriza por una turbulencia muy energtica y de
pequea escala y gran cizalladura. Los vrtices turbulentos que se generan en esta
regin se mueven hacia dentro y hacia fuera de la subcapa viscosa. Si los elementos
rugosos del fondo (los granos de arena o las formas del lecho) tienen una mayor altura
que la subcapa viscosa, entonces la capa de generacin de turbulencia se extiende hasta
el fondo.
La capa exterior constituye la mayor parte de la capa turbulenta y se caracteriza por
vrtices ms grandes, que son ms eficientes para el transporte de momentum. Es esta
gran eficiencia en el transporte de momentum lo que produce que el perfil medio de
velocidad sea mucho ms suave que la capa de generacin de turbulencia.
La subcapa viscosa se caracteriza por tener fluctuaciones de velocidad (turbulencia),
pero las fluctuaciones de velocidad en direccin normal al fondo (contorno) tienden a
cero en sus proximidades. Por lo tanto, cerca del fondo, el transporte de momentum del
fluido est dominado por la difusin molecular sobre la difusin turbulenta, y la tensin
de corte en el fondo se puede suponer en principio como una aproximacin a la
tensin de corte de una capa lmite de flujo laminar (aquellos que fluyen a muy poca
velocidad). El flujo ser entonces en capas que se deslizan suavemente unas sobre otras
sin mezcla macroscpica de partculas de fluido entre ellas. Entonces las tensiones de
corte en un flujo de esta naturaleza estn dadas por la segunda ley de Newton (donde
interviene la viscosidad del fluido):
3.21
u
z
[3.3]
donde es la densidad del fluido, es la viscosidad cinemtica ( 10-6 m2 /s) y u es la

velocidad horizontal. Para flujos con una capa lmite turbulenta, los modelos de las
tensiones de corte varan desde modelos muy simples (ecuacin anterior) donde el valor
de la viscosidad puede aumentarse artificialmente (con respecto a su valor laminar)
para tomar en consideracin el incremento de la mezcla asociada a los flujos turbulentos
(viscosidad turbulenta t), hasta modelos numricos de cierre para la turbulencia muy
sofisticados. Una revisin mas extensa de otros modelos alternativos que emplean
esquemas mas sofisticados de cierre turbulento puede verse en Madsen &
Wikramanayake (1991).
Los modelos de capa lmite turbulenta pueden simplificarse artificialmente escogiendo
un valor constante de la viscosidad turbulenta t, que es mucho mayor que su valor
laminar (molecular) y que refleja el tamao y estructura de los vrtices asociados al
flujo turbulento. Esta suposicin resulta en un modelo muy simple y conceptual de
viscosidad turbulenta para la tensin de corte turbulenta, que puede explicarse de la
siguiente manera.
En la figura 3.6 se puede ver al fluido movindose en una serie de capas de espesor l
con velocidad variable u = u(z) en direccin longitudinal al contorno del fondo,
concretamente se muestran dos capas vecinas, en la que la capa inferior se mueve a u1 , y
la capa superior se mueve a u2 u1 ( u z ) l ' . En este flujo turbulento, u(z) representa
la velocidad promediada en el tiempo (es decir, sin las fluctuaciones de las velocidades
turbulentas). Ya que el flujo medio es paralelo al contorno, la velocidad media en la
componente vertical necesariamente es cero. No obstante, las fluctuaciones turbulentas
en la vertical (w) darn lugar a un intercambio de fluido a travs del lmite entre ambas
capas. Debido a que no hay transporte neto de fluido entre las dos capas (si lo hubiera,
la velocidad vertical sera no nula), las tasas de transferencia vertical ascendente y
descendente de volumen por unidad de rea horizontal deben ser iguales.
Figura 3.6.- Transferencia de momentum y mezcla (tensiones de corte) en capas lmite turbulentas
(Madsen, 1993).
3.22
Suponiendo que una parcela de fluido viaja una cierta longitud l antes de que se
transfiera su momentum, es decir, para que este fluido recin llegado forme parte de la
capa superior, de acuerdo con la teora de la longitud de mezcla de Prandtl, su velocidad
debe incrementarse en:
u
u ' = u2 u1 l '
z
[3.4]
Por lo tanto, la capa superior acelerar al fluido proveniente de la capa inferior. De

acuerdo con la segunda ley de Newton, la tasa de intercambio de momentum es igual a
una fuerza, por lo que la transferencia vertical ascendente de fluido desde la capa
inferior (con menor velocidad) hacia la capa superior (con mayor velocidad) es
equivalente a una fuerza por unidad de rea (es decir, la tensin de corte) que acta
sobre la capa superior en direccin opuesta la flujo. De igual manera, considerando el
problema a la inversa, si hay una transferencia vertical descendente del fluido desde la
capa superior (con mayor velocidad) hacia la capa inferior (con menor velocidad) se
producir una tensin de corte de igual magnitud, pero actuando en la misma direccin
del flujo. Para simplificar las cosas, en el anlisis de capas lmite para flujos turbulentos,
generalmente se adopta un modelo simple de viscosidad turbulenta. Este modelo puede
ser derivado de la hiptesis de longitud de mezcla turbulenta clsica de Prandtl
(Schlichting, 1960), en donde se supone que la fluctuacin de velocidad instantnea
vertical tiene la misma magnitud que la horizontal (es decir, las fluctuaciones
turbulentas de la velocidad horizontal (u) y las fluctuaciones de la velocidad vertical
(w) son comparables). A partir de la discusin anterior de la diferencia de velocidades
horizontales asociadas con la transferencia vertical de fluido entre dos capas vecinas, se
desprende que las fluctuaciones de la velocidad estn determinadas por la diferencia de
velocidades entre las dos capas, se tiene que
u ' = ( u1 u2 ) = l '
u
w'
z
[3.5]
Por lo tanto, el producto de la magnitud de las fluctuaciones de velocidad turbulentas

(w) con el grosor (l) de la capa individual (que est fsicamente relacionado con la
longitud de mezcla asociada con las fluctuaciones de la velocidad vertical) es un
equivalente de la viscosidad cinemtica (Madsen, 1993). Estos son los mismos
argumentos presentados para conceptuar la naturaleza de las tensiones de corte en flujos
laminares (Bird, Stewart and Ligthfoot, 1960), en el cual, w y l estn asociados al
movimiento molecular. Sin embargo, para flujos turbulentos es necesario relacionar la
escala de velocidades turbulentas w y la longitud de mezcla l con las propiedades del
flujo turbulento:
3.23
= u ' w ' = ( l ')
u u
2 du
= ( l ')
z z
dz
[3.6]
Tanto fsicamente como dimensionalmente, la viscosidad turbulenta se puede

relacionar con las fluctuaciones de velocidad u y w como t = ( l ' ) ( u z ) .
2
Dado que los flujos de corte en una capa lmite sobre un fondo horizontal estn
caracterizados por una mayor escala de variacin en el plano horizontal con respecto a
la escala de variacin en el plano vertical, es posible despreciar las aceleraciones
verticales dentro de la capa lmite, por lo que el gradiente de presin horizontal es aquel
impuesto en el borde exterior de la capa lmite. Para flujos bidimensionales horizontales
en el plano x-z siendo x la horizontal y z la vertical, la forma linealizada de la ecuacin
de la capa lmite (figura 3.4) es:
u
p
=-
t
x z
[3.7]
donde es la densidad del fluido, u es la velocidad horizontal, p es la presin en el

borde exterior de la capa lmite, es decir en z = , es la tensin de corte, y t es el
tiempo. Integrando la ecuacin [3.7] desde el fondo, donde z = 0, hasta una elevacin z
por encima del fondo, es decir:
z u
p
dz = b =
+ dz
0
z
t x
[3.8]
Si se impone que el lmite superior de integracin est inmediatamente por encima del
fondo (es decir en z = 0+) el factor de la derecha se hace muy pequeo, con lo que
inmediatamente por encima del fondo:
b = tensin de corte en el fondo
[3.9]
que es lo que se conoce como la ley de la pared.

Combinando las ecuaciones [3.6] y [3.9], se obtiene la expresin para la escala de las
fluctuaciones de velocidades turbulentas:
w' = l '
u
=
z
b
= u*
[3.10]
3.24
donde u* es la velocidad de corte. Por lo tanto, la escala de fluctuaciones de velocidades

turbulentas est relacionada con la magnitud de la tensin de corte en el fondo, que a su
vez depende del flujo.
Para terminar de definir las tensiones de corte turbulentas, es necesario parametrizar la
longitud de mezcla l. Tomando en cuenta que la validez de la ley de la pared se limita
a las inmediaciones del fondo, donde el movimiento vertical de las partculas del fluido
est limitado por la presencia del contorno (es decir, conforme z 0), es de esperar que
l 0. Por lo tanto, la longitud de mezcla en las vecindades del fondo debe ser
proporcional a la distancia z por encima del fondo, es decir, l = k z, donde la constante
de proporcionalidad k es la constante de von Karman.
Con l = k z y w = u* el modelo para las tensiones de corte turbulentas puede
expresarse como:
= t
u
z
[3.11]
donde la viscosidad turbulenta es :
t = k u *z
u
z
[3.12]
donde k es la constante de von Karman ( 0.41), y
u* =
[3.13]
es la llamada velocidad de corte, donde b es la tensin de corte en el fondo (z = 0).

Propiamente dicho, u* no es una velocidad real ( b tiene unidades de [m/s]) pero
est relacio nada con el factor de friccin (de hecho, se considera una forma alternativa
de expresar la friccin en trminos de una velocidad).
Como se ha mencionado antes, las tensiones de corte en el fondo b son debidas a las
aportaciones de las tensiones de corte turbulentas y viscosas. No obstante, en el fondo
no hay turbulencia (dado que u = w = 0, u = w= 0). Por lo tanto, de la ecuacin [3.6] se
tiene que = ( u ' w ') = 0 . Por lo tanto, en una capa muy delgada cerca del fondo, las
tensiones de corte viscosas son dominantes sobre las tensiones de corte turbulentas y
por lo tanto, el flujo es laminar. Esta delgada capa es la llamada subcapa viscosa. Por
3.25
encima de la subcapa viscosa, es decir, en la mayor parte del flujo, las tensiones de corte
turbulentas son dominantes. Algunos resultados de medidas de las tensiones de corte en
la subcapa viscosa muestran que esta es constante e igual a las tensiones de corte en el
fondo (Liu, 2001).
Aunque el desarrollo de este modelo conceptual para las tensiones de corte fue llevado
a cabo considerando las condiciones en las proximidades del fondo, en la capa lmite (es
decir, para valores pequeos de z), los resultados de los clculos de las tensiones de
corte, con t como se expresa en la ecuacin [3.12] aplicada para todo el grosor de la
capa lmite, tienen un alto grado de precisin (Madsen, 1993).
3.6.- Corrientes
Las corrientes pueden ser provocadas por mareas, el viento, gradientes de presin
atmosfrica, la rotura del oleaje, descargas de ros en el mar, gradientes de densidad,
etc. En la zona cercana a la costa, las corrientes inducidas por la rotura del oleaje son las
dominantes y tienen la capacidad de transportar grandes cantidades de sedimento en
suspensin. Por ello es fundamental describir el perfil de velocidades de las corrientes.
El flujo (en la columna del agua) se puede clasificar desde el punto de vista fsico y
comenzando desde el fondo, (figura 3.7) de la siguiente manera:
Subcapa viscosa. Una delgada capa por encima del fondo donde casi no hay
turbulencia y las tensiones de corte viscosas son constantes. Por lo tanto, el flujo
es laminar. Por encima de esta capa el flujo es turbulento.
Capa de transicin (tambin llamada capa buffer) donde la viscosidad y la

turbulencia tienen la misma importancia relativa.
Capa turbulenta con perfil logartmico, donde las tensiones de corte viscosas son
despreciables en comparacin con las tensiones de corte turbulentas. Aqu las
tensiones turbulentas son constantes e iguales que las tensiones de corte en el
fondo. Es en esta capa donde Prandtl introdujo el concepto de longitud de
mezcla y obtuvo el perfil de velocidades logartmico.
Capa exterior turbulenta, donde las velocidades son casi constantes debido a la
presencia de grandes vrtices que producen una fuerte mezcla en el flujo.
En la capa turbulenta con perfil logartmico de velocidades, las tensiones de corte son
constantes e iguales a las tensiones de corte en el fondo. Prandtl obtuvo el perfil de
velocidades logartmico suponiendo que la longitud de mezcla es proporcional a la
distancia con el fondo (es decir, l = k z ).
Aun cuando se han propuesto varias expresiones para la distribucin (perfil) de
velocidades en la capa de transicin y en la capa turbulenta exterior, en la actualidad
ninguna de ellas est totalmente aceptada. No obstante, modificando la suposicin de la
longitud de mezcla (siguiente seccin), el perfil logartmico de velocidades se puede
3.26
aplicar tambin a la capa de transicin y la capa exterior turbulenta, lo cual es
razonablemente preciso cuando se comparan las medidas con los clculos de
velocidades en esta regin.
Figura 3.7.- Clasificacin de las regiones de flujo. Los grosores no estn a escala y la capa exterior
turbulenta es del orden del 80 al 90% del flujo (Liu, 2001).
As, desde el punto de vista ingenieril, la capa de transicin y la capa exterior turbulenta
pueden ser descritas con un perfil de velocidades logartmico (figura 3.8) (como con la
subcapa viscosa).
Nikuradse (1933) fue el primero en estudiar el efecto de la rugosidad de los contornos
en la distribucin de velocidades y fue l quien introdujo el concepto de la rugosidad
equivalente de grano k n (o Rugosidad de Nikuradse o rugosidad del fondo).
Experimentalmente verific los siguientes puntos:
u * kn
5 . La rugosidad del fondo es
mucho mas pequea que el grosor de la subcapa viscosa. Por lo tanto, la

rugosidad del fondo no afectar la distribucin de velocidades.
Flujos moderadamente turbulentos para
Flujos altamente turbulentos para
Flujos de transicin para 5
u * kn
70 . La rugosidad del fondo es tan
grande que produce vrtices turbulentos cerca del fondo, por lo que no existe
una subcapa viscosa y la velocidad del flujo no depende de la viscosidad.
u * kn
70 . La distribucin de velocidades se ve
afectada por la rugosidad del fondo y por la viscosidad.
3.27
Figura 3.8.- Clasificacin ingenierl de las distintas regiones del flujo (Liu, 2001). Nota: el grosor de las
capas no est a escala.
En la capa turbulenta, las tensiones de corte se deben exclusivamente a la aportacin de

las tensiones de corte turbulentas t (figura 3.9), donde las tensiones de corte totales se
incrementan linealmente con la profundidad, es decir:
z
t ( z) = b 1
h
[3.14]
De acuerdo con la teora de la longitud de mezcla de Prandtl:
du
t = l
dz
[3.15]
y suponiendo que la longitud de mezcla es:
z
l = kz 1
h
0.5
[3.16]
donde k es la constante de vo n Karman (k = 0.4), podemos obtener:
b
du
u
=
= *
dz
kz
kz
[3.17]
3.28
Integrando la expresin anterior se obtiene el perfil de velocidades, cuya forma es
logartmica:
u (z ) =
u* z
ln
k z0
[3.18]
donde el significado fsico de la constante de integracin z0 es la elevacin

correspondiente a un valor nulo de velocidad, es decir, donde u z= z = 0 , que segn
0
Nikuradse (1933) tiene los siguientes valores:

u k
para flujos moderadamente turbulentos * n 5

9u*
k
u k
z 0 = n para flujos altamente turbulentos * n 70

30
k
u k
+ n para flujos de transicion 5 < * n < 70

9u* 30

[3.19]
donde k n es la rugosidad de equivalente de grano de Nikuradse. Es interesante notar que

la velocidad de corte, que por definicin no tiene nada que ver con la velocidad del
flujo, es la velocidad del flujo a una elevacin de z = z0 e k , es decir:
uz =z
0e
= u*
[3.20]
Figura 3.9.- Distribucin de las componentes de las tensiones de corte en el fondo.
En los estudios relacionados con la cuantificacin del transporte de sedimentos, es

importante notar que la velocidad de corte, adems de estar relacionada con la b y la
3.29
friccin del flujo con el fondo, es la velocidad del flujo muy cerca del fondo (figura
3.10).
Por otra parte, en el caso de los flujos moderadamente turbulentos existe la subcapa
viscosa. La tensin de corte viscosa es constante en esta capa e igual a las tensiones de
corte en el fondo, es decir:
Aplicando la definicin que u
z =0
du
= b
dz
[3.21]
= 0 e integrando con respecto de z, se obtiene:
u2
u (z) = z = * z
[3.22]
lo que significa que el perfil de velocidades dentro de la subcapa viscosa es lineal. Esta
distribucin lineal de velocidades se intersecta con el perfil logartmico de velocidades a
una elevacin de z = 11.6 ( u* ) . Estos perfiles de velocidades se muestran en la figura
3.10 con una descripcin detallada de la velocidad del fluido cerca del fondo.
Figura 3.10.- Perfil de velocidades en flujos moderadamente y altamente turbulentos (Liu, 2001).
3.30
3.6.1.- Perfil de velocidades y rugosidad del fondo
Una corriente fluyendo sobre el fondo del mar experimenta una friccin con el fondo,
lo que provoca la generacin de una capa lmite turbulenta, que en este caso puede ser
del orden de varios metros de grosor. En aguas someras, la capa lmite puede ocupar la
totalidad de la profundidad, mientras que en aguas profundas, la capa lmite puede
ocupar tan solo la parte mas profunda de la columna de agua, por debajo de una capa
donde la corriente no se ve afectada por la presencia del fondo. Dentro de la capa lmite,
la velocidad de la corriente se incrementa desde cero (su valor en el fondo) a un mximo
en o cerca de la superficie del agua (perfil de velocidades), donde el mayor gradiente de
velocidades se encuentra precisamente en las cercanas del fondo.
Debido a que la forma mas comn de definir la velocidad de la corriente es su valor
verticalmente integrado U , se puede relacionar esta velocidad con el perfil de
1 h
velocidades mediante su definicin U = u( z ) dz , donde U es la velocidad de la
h 0
corriente verticalmente integrada, h es la profundidad y u(z) es la velocidad de la
corriente a una cierta elevacin sobre el fondo z. Cabe mencionar que si el perfil de
velocidades tiende a cero a una profundidad z = z0 (por ejemplo, para un perfil
logartmico de velocidades), el lmite inferior de esta integral debe ser cambiado de 0 a
z0 .
En comparacin con el flujo oscilatorio del oleaje, una corriente puede considerarse
como un flujo razonablemente estacionario. Si la tensin del fondo debida a las
corrientes es c, el perfil de la velocidad de la corriente puede expresarse como:
u (z ) =
u*c z
ln
k z0
[3.23]
donde u*c = bc es la velocidad de corte de la corriente. Por lo tanto, suponiendo un

flujo estacionario no estratificado sobre un fondo horizontal (que puede tener ripples),
la velocidad de la corriente u vara con la altura sobre el fondo (z) en forma logartmica.
Esta ecuacin es un perfil logartmico clsico de velocidades donde la constante de
integracin (que representa una escala de longitud y es totalmente arbitraria) se expresa
en trminos z0 , que es el valor de z en el cual el perfil logartmico de velocidades
predice una velocidad igual a cero. Resulta curioso que esta condicin de friccin no se
encuentre en z = 0, que es donde se localiza el contorno slido.
Para un fondo no rugoso, esta condicin tiene una posicin bien definida, que es el
propio contorno, es decir, z0 = z = 0. Para este tipo de fondo, la viscosidad turbulenta
(ecuacin [3.12]) se hace mas pequea que la viscosidad molecular para valores de
z < ( ( ku* ) ) . Por lo tanto, el flujo en las inmediaciones del contorno, concretamente
en la subcapa viscosa, es laminar y no turbulento. Para un contorno rugoso, la
localizacin del contorno no es un valor nico de z, por lo que fijar el valor z = 0, puede
ser una definicin un tanto ambigua de la posicin terica del fondo, ya que algunas
3.31
partes estarn localizadas en z > 0 y otras en z < 0. Por lo tanto, la aplicacin de la
ecuacin [3.23] en las inmediaciones de un fondo slido es nicamente una formalidad,
y por lo tanto, la prediccin de uc = 0 en z = z0 no tiene significado fsico.
El rango de posiciones (alturas) sobre el fondo para las cuales la ecuacin [3.23] es
vlida (lejos de estructuras costeras y fuera de la zona de rompientes) es desde unos
cuantos centmetros por encima del fondo hasta una altura aproximada del 20 al 30% de
la columna de agua (es decir, aproximadamente a z = 2 a 3 m en aguas someras o z = 20
a 30 m en aguas profundas). Soulsby (1990) propone una ecuacin para describir el
perfil de velocidades de una corriente de marea (que puede considerarse razonablemente
estacionaria):
z
U ln
z 0 para z < z < 0.5
u (z) =
0

ln
2 z0 2h
[3.24]a, b
z
U ln
2z0
u (z) =
para0.5 < z < h

ln
2 z0 2h
donde U es la velocidad verticalmente integrada, h es la profundidad del agua y es el

grosor de la capa lmite. La ecuacin [3.24]b es compatible con la ecuacin [3.23] cerca
del fondo.
Alternativamente a estas ecuaciones, se puede encontrar el perfil de velocidades de las
corrientes de marea mediante las siguientes frmulas empricas:
z 7
U (z) =
U para z0 < z < 0.5
0.32 h
[3.25]a, b
U ( z ) = 1.07U para0.5 < z < h
La rugosidad del fondo k n , conocida tambin como la rugosidad equivalente de arena

de Nikuradse es una forma de medir la escala o longitud de la rugosidad del fondo.
En realidad la nica situacin donde se puede obtener directamente la rugosidad del
fondo es aquella donde este consista de esferas uniformes, donde k s = dimetro de la
esfera. Pero en la naturaleza, los fondos estn compuestos de granos con distintos
tamaos (irregulares) y rara vez son planos, con formas geomtricas como ripples o
3.32
dunas que aparecern dependiendo del tamao del grano y de las caractersticas del
flujo. En este caso, la rugosidad del fondo puede obtenerse indirectamente a partir de la
medida de las velocidades, realizando un ajuste logartmico del perfil de velocidades
por mnimos cuadrados para obtener la velocidad de corte u* , utilizando k s = 30 z0
(ecuacin [3.19]) y b = u* 2 . Madsen (1993) propone una pequea simplificacin del
criterio del rgimen del flujo, considerando nicamente flujos moderadamente
turbulentos para valores de k n (u* / ) 3.3 z0 = ( / 9u* ) y flujos altamente
turbulentos para valores de k n (u* / ) 3.3 z0 = (kn / 30).
La suposicin del perfil logartmico de velocidades sugiere que la velocidad mxima
ocurre en la superficie del flujo. No obstante, es comn que las velocidades medidas
indiquen que el mximo se encuentra a una cierta distancia de la superficie, debido a la
friccin de la superficie con el aire. De hecho, el perfil de velocidades logartmico se
desarrolla bsicamente para la capa turbulenta (logartmica) cerca del fondo. Por lo
tanto, para la determinacin de las tensiones de corte en el fondo (b ) y de la rugosidad
equivalente de Nikuradse (k n ) a partir de las velocidades deben realizarse
preferentemente con medidas realizadas a una elevacin de k n < z < 0.2 h.
Por otra parte, Christoffersen & Jons son (1985) proponen una expresin de ajuste entre
la viscosidad cinemtica del agua, el nmero de Reynolds del Sedimento
Re s = ( u* k n ) y la escala (longitud) de la rugosidad del fondo z0 :
z0 =
kn
u k
1 exp * n +
30
27 9u*
[3.26]a
Anteriormente, Colebrook & White (1937) propusieron una versin simplificada de la

expresin anterior:
z0 =
kn
+
30 9u*
[3.26]b
u * kn
> 70 ), la expresin de Colebrook &
White (1937) se reduce a una forma idntica a las ecuaciones [3.19]:
que para un flujo altamente turbulento (para
z0 =
kn
30
mientras que si el flujo es moderadamente turbulento, (para

Colebrook & White (1937) se reduce a:
[3.26]c
u * kn
< 5 ), la expresin de
3.33
z0 =
9u*
[3.26]d
Si el flujo se encuentra en transicin, se debe usar la expresin completa.

An cuando las caractersticas propias de los materiales (potencialmente distintos) que
compone n el fondo pueden provocar que el rgimen del flujo sea distinto (fondos de
sedimentos finos o muy finos estn asociadas a flujos suaves o de transicin, mientras
que las arenas gruesas y las gravas estn asociados a flujos hidrodinmicamente ms
intensos), es comn considerar todos los flujos sobre arenas como hidrodinmicamente
intensos ya que es matemticamente mas simple realizar los clculos. Las
aproximaciones de la expresin de Colebrook & White (1937) conllevan un error de
menos del 10% en los clculos de u* en la regin posterior al umbral del movimiento
para granos mayores a 60 m (Soulsby, 1997).
Se han propuesto una gran cantidad de expresiones para la escala de rugosidad z0 para
un lecho de arena sin ripples en trminos de la rugosidad de Nikuradse k n y en funcin
de la distribucin granulomtrica del sedimento (Sleath, 1984, Van Rijn, 1993, etc.)
aunque quizs la ms frecuente sea en funcin del dimetro medio del sedimento:
k n = 2.5 d 50
[3.27]
Es interesante notar que, an cuando hay una discrepancia importante entre las distintas
alternativas para el clculo de k n , y por lo tanto, del clculo de las tensiones de corte, la
incertidumbre ser menor, ya que stas dependen logartmicamente de k n .
Combinando las ecuaciones [3.26]c y [3.27], se obtiene una relacin directa de z0 con el
tamao del sedimento para flujos hidrodinmicamente intensos:
z0 =
d 50
12
[3.28]
Madsen (1993) sostiene que para flujos turbulentos sobre un fondo plano consistente en
material granular, se suele tomar k n = d50 = tamao de grano del sedimento del fondo.
En cambio, para un flujo sobre un fo ndo irregular con elementos rugosos (por ejemplo,
los ripples) el valor k n = 30 z0 se refiere a la rugosidad equivalente de los granos de
arena de Nikuradse, donde z0 se obtiene extrapolando el perfil logartmico de
velocidades por encima del fondo hasta un valor de z = z0 , donde la velocidad de la
corriente u (z) se desvanece.
Mathisen (1993) realiz unos experimentos de laboratorio donde obtuvo un perfil de
velocidades de una corriente fluyendo sobre un fondo con barras triangulares de 1.5 cm
3.34
de altura y con una separacin de 10 cm (mostrado en un grfico semilogartmico en la
figura 3.11). Sus medidas mostraron que sobre las crestas y en medio de los elementos
rugosos, las medidas de velocidades se desviaban de la lnea recta prevista por debajo
de los primeros 2.5 a 3 cm sobre el fondo y tambin por encima de 10 cm. Las
desviaciones de las velocidades por encima de los 10 cm. (el perfil no- logartmico de
velocidades) se deben a que el flujo por enc ima de esta profundidad no se encuentra
afectado por la resistencia del flujo con el contorno y la turbulencia de las paredes del
canal. La desviacin de las velocidades (del perfil logartmico) cerca del fondo se debe
a la propia configuracin del fondo, es decir, a la proximidad de los elementos rugosos.
Aqu, la velocidad es funcin de la posicin, es decir, el flujo no es uniforme en las
inmediaciones de los elementos rugosos.
Figura 3.11.- Perfil de velocidades (semi logartmico) de una corriente fluyendo sobre elementos rugosos
bidimensionales (Mathisen, 1993).
El aumento logartmico de la velocidad (uc vs. log z) queda bien definido en la regin
entre 3 cm < z < 10 cm. A partir del perfil de velocidades de la corriente de esta regin,
se puede obtener, por extrapolacin a u (z) = 0, el valor de z = z0 0.7 cm o k n = 30 z0
21 cm (debido a que el flujo es altamente turbulento). Este ejemplo muestra claramente
que k n , la rugosidad equivalente de Nikuradse, es funcin de la configuracin de la
rugosidad del fondo, y no necesariamente refleja la escala fsica de las formas del fondo.
As, en los experimentos realizados por Mathisen (1993), donde la escala de rugosidad
(la altura de las barras triangulares) de 1.5 cm, result en una rugosidad equivalente de
Nikuradse de 21 cm. Este mayor valor de k n est asociado a una mayor resistencia del
flujo producido por las barras triangulares.
3.35
3.6.2.- Fuerzas del flujo sobre el sedimento
3.6.2.1.- Fuerzas de arrastre y sustentacin

La tensin de corte en el fondo (o la friccin con el fondo) es la fuerza de friccin por
unidad de rea que ejerce una corriente al fluir. Por lo tanto, para los clculos del
transporte de sedimentos es de particular importancia su evaluacin, ya que representa
la fuerza inducida por el flujo sobre los granos de arena del fondo (aunque generalmente
es ms conveniente expresar la tensin de corte en el fondo asociada con el flujo en la
capa lmite en trminos de la friccin de la corriente con el fondo a travs de un factor
de friccin f c). Suponiendo que no hay formas geomtricas del fondo y que el transporte
de sedimento no es muy intenso, es razonable suponer que las tensiones de corte totales
sern nicamente debidas a las provocadas por la friccin entre el fluido y el fondo. No
obstante, considerando un fluido real movindose alrededor de un cuerpo slido, el flujo
ejerce una fuerza de arrastre sobre dicho cuerpo (figura 3.12).
Figura 3.12.- Fuerzas de arrastre y sustentacin ejercidas por el flujo sobre un cuerpo slido (Liu, 2001).
La fuerza de arrastre se compone de la friccin de arrastre y de la friccin de la forma

geomtrica. La componente de la friccin de arrastre tiene su origen en la proyeccin de
la fuerza de friccin en la direccin del flujo, mientras que la componente de la friccin
de la forma geomtrica tiene su origen en la proyeccin de la presin ejercida por el
flujo sobre las formas (irregularidades) geomtricas del cuerpo (sobre la misma
direccin del flujo). Esta fuerza de arrastre total se puede escribir como:
3.36
1
FD = CD AU 2
2
[3.29]
Por su parte, la componente de sustentacin se puede expresar de la misma manera:
FL =
1
CL AU 2
2
[3.30]
donde U es la velocidad media del flujo, A es el rea proyectada del cuerpo al plano
perpendicular a la direccin del flujo y CD y CL son los coeficientes de arrastre y
sustentacin, que dependen de la forma y de la rugosidad de la superficie del cuerpo en
cuestin y del nmero de Reynolds, todos ellos, por lo regular, determinados
experimentalmente.
3.6.2.2.- Coeficiente de friccin

El esquema de definicin del coeficiente de friccin del flujo sobre los granos de arena
depositados en el fondo se muestra en la figura 3.13.
Figura 3.13.- Fuerzas del fluido actuando sobre los granos del fondo (Liu, 2001).
La fuerza de arrastre del flujo sobre el grano se puede expresar como:
2
1
FD = CD A ( U )
2
[3.31]
3.37
donde es un trmino que se incluye debido a que la velocidad del fluido sobre el
grano se desconoce, aunque es razonable suponer que es funcin de la velocidad media
del flujo y de otros parmetros. De la misma manera se puede suponer que los granos
ejercen una fuerza de resistencia FD sobre el flujo. Considerando que el rea del grano
proyectada en el plano horizontal es A, entonces las tensiones de corte en el fondo se
pueden expresar como:
bc =
FD 1
1
2 A
2
2
= CD
U = f c U
A' 2
A'
2
[3.32]
donde f c es el coeficiente (adimensional) de friccin del fondo. sta expresin de bc

tambin se puede expresar en trminos de una cierta velocidad conocida uc a un cierto
nivel de referencia z r, que indica que el factor de friccin de la corriente f c es funcin del
nivel de referencia, z = zr, en el cual se define un valor para la velocidad de la corriente
(Madsen, 1993).
bc =
2
1
f c ( uc ( zr ) )
2
[3.33]
Combinando las ecuaciones [3.13] y [3.32], se obtiene el valor de las tensiones de corte
debidas a las corrientes, que tambin se pueden expresar en trminos del factor de
friccin f c y la velocidad media:
bc
f
= c U
u*c =
[3.34]
o de la velocidad a un cierto nivel de referencia uc (zr):
f
u*c = c
2
uc ( z r )
[3.35]
Introduciendo esta ltima expresin en la ecuacin [3.23], con z = zr y k = 0.4, se

obtiene una ecuacin para el factor de friccin de la corriente en trminos de la
elevacin del nivel de referencia y la escala de rugosidad del fondo:
1
4 fc
log 10
zr
z0
[3.36]
3.38
En flujos altamente turbulentos (donde se ha definido que z 0 = kn 30 ), la ecuacin
[3.36] es una ecuacin explcita para f c en trminos de la rugosidad relativa z r k n . En
flujos moderadamente turbulentos (donde se ha definido que z 0 = 9u*c ), a partir de la
ecuacin [3.36] se deriva una ecuacin implcita para f c (Madsen 1993):
1
z u (z )
= log 10 r c r + 0.20
+ log 10
4 fc
4 fc
1
[3.37]
que muestra que cuando el flujo es moderadamente turbulento, el factor de friccin de la

corriente f c depende del nmero de Reynolds ( z ru c ) .
Para la determinacin del factor de friccin de una corriente definida por su velocidad
uc (zr) y una cierta rugosidad del fondo k n , se puede suponer en primera instancia, que el
flujo es altamente turbulento (como sucede en la mayor parte de los casos en la
naturaleza), por lo que z 0 = kn 30 . Con esta suposicin, se puede obtener fcilmente el
factor de friccin de las corrientes f c a partir de la ecuacin [3.36] y u*c de la ecuacin
[3.34]. Evaluando el rgimen del flujo, si se confirma que k n u* > 3.3 , entonces el
flujo es efectivamente altamente turbulento. Sin embargo, si k n u* < 3.3 , entonces el
flujo es moderadamente turbulento, y se puede proceder de dos maneras. En primer
lugar se puede resolver la ecuacin implcita para f c, (ecuacin [3.37], de la misma
manera que se determina la ecuacin de friccin del oleaje (como se ver mas adelante),
cuya forma es muy similar a la de f c.
La segunda opcin es un proceso iterativo para determinar z 0 = 9u*c , donde (en
principio y suponiendo un flujo altamente turbulento) se obtiene u*c, y con este valor de
z0 se resuelven las ecuaciones [3.36] y [3.34] para obtener un nuevo valor de u*c. Este
u*c se puede usar de nuevo para recalcular z 0 = 9u*c y as continuar con un proceso
iterativo, en el cual, es de esperar que la convergencia sea rpida, ya que z0 es parte del
argumento de un logaritmo.
3.6.2.3.- Coeficiente de Chzy

Por otra parte, las tensiones de corte en el fondo se pueden relacionar con la velocidad
de la corriente verticalmente integrada U a travs de un coeficiente de arrastre CD por
medio de una ley de friccin cuadrtica:
bc = CDU
[3.38]
3.39
donde el coeficiente de arrastre puede expresarse de distintas formas, incluyendo el
coeficiente de Darcy-Weisbach (f) el coeficiente de Chzy (C), o el coeficiente de
Manning (n), que aplicados al agua de mar, se pueden relacionar matemticamente de la
manera siguiente:
CD =
f
g
gn 2
= 2 = 1
8 C
h3
[3.39]
El coeficiente de Chzy se plante originalmente como un coeficiente emprico para la

velocidad media en un flujo uniforme en un canal:
U = C RS
[3.40]
donde R es el radio hidrulico, es decir el rea de la seccin transversal dividida por el

parmetro mojado, S es la pendiente del fondo y C es el coeficiente emprico llamado
coeficiente de Chzy (que en principio es constante) y para el cual se han propuesto
varias frmulas. Una de las formas en la que es posible obtener la expresin terica del
coeficiente de Chzy es a partir del perfil logartmico de velocidades. Considerando que
en el caso de flujos en canales, la velocidad de corte se expresa como:
u* = gRS
[3.41]
Entonces, combinando ambas ecuaciones ([3.40] y [3.41]), se obtiene una expresin

para el coeficiente de Chzy en funcin de la velocidad media y la velocidad de corte:
C=
U
u*
[3.42]
Ya que la velocidad media se puede obtener a partir de integracin del perfil

logartmico de velocidades, es decir:
U=
1 h
u h z
u ( z ) dz = * ln
h z0
kh z0 z0
dz
u h
z u h
* ln 1 + 0 * ln
k z 0
h k z 0 e
[3.43]
3.40
e insertando la ecuacin [3.43] en la ecuacin [3.42], y utilizando la expresin de z0 y
convirtiendo el logaritmo natural en logaritmo en base 10, se obtiene:
12h
18log
3.3
g h
u*
C=
ln
k
z
e
0
12h
18log
kn
flujos hidrulicamente suaves u*kn 5

[3.44]
u k
flujos hidrulicamente fuertes * n 70

m
de donde se desprende que C tiene unidades de
.
s
Igualando el coeficiente de Chzy (ecuaciones [3.42] y [3.44]) con el factor de friccin

de la corriente (ecuaci n [3.34]), ambos en trminos de la velocidad de friccin u* , se
obtiene una expresin para el factor de friccin en funcin del coeficiente de Chzy:
0.06
flujos hidraulicamente suaves
12
h
log

3.3
2 g
u*
f c = 2 =
C
0.06
flujos hidraulicamente intensos

2
12
h
log

kn
u* kn
[3.45]
u* k n
70

De esta manera, (al igual que con el coeficiente de Chzy), se puede calcular el
coeficiente de friccin (ecuacin [3.45]) a partir de conocer la velocidad media y
obtener as la tensin de corte en el fondo (ecuacin [3.32]). Cabe recalcar que, debido a
que el coeficiente de friccin f c es adimensional, por lo regular se prefiere su uso sobre
el coeficiente de Chzy.
As, conociendo la velocidad media (ecuacin [3.43]), que generalmente es el caso y las
caractersticas del fondo (por ejemplo, el tamao medio de grano y/o la altura estimada
de ripples presentes), se puede obtener el coeficiente de Chzy (ecuacin [3.44]), la
velocidad de corte en el fondo a partir de este coeficiente (ecuacin [3.42]) y la tensin
de corte en el fondo (ecuacin [3.13]).
As, el trmino de friccin CD depende entonces de la escala de rugosidad del lecho z0 y
de la profundidad del agua h. Una manera simple de establecer esta dependencia es
mediante una ley potencial (Soulsby, 1997):
3.41
z
C D = 0
h
[3.46]
donde CD es el coeficiente de friccin involucrado en las ecuaciones [3.38] y [3.39] y

y son coeficientes para los que existen diversas expresiones:
Manning Strickler (Soulsby, 1997): = 0.0474 y = 1/3
Dawson et al., (1983): = 0.0190 y = 0.208
Mediante el ajuste de datos experimentales de
u*
d
contra 50 se obtiene la siguiente
U
h
expresin:
u* 1 d50 7
=
U 7 h
[3.47]
Sustituyendo la ecuacin [3.28] (z0 = d50 /12) en la ecuacin [3.47] se obtienen valores
de = 0.0415 y = 2/7 en la ecuacin [3.46]. Estos valores son recomendados
(Soulsby, 1997) para calcular las tensiones de corte debidas a la friccin con el
sedimento en lugar de los recomendados por Manning Strickler (Soulsby, 1997) y
Dawson et al. (1983) que fueron derivados para otros propsitos. La ecuacin [3.47] es
compatible con la ecuacin para el perfil de velocidades para una corriente dada en la
ecuacin [3.25]. Igualando las ecuaciones [3.25] y [3.47] se obtiene una ecuacin para
la velocidad de corte a partir de la velocidad medida u(z) a una altura z cerca del fondo:
d 7
u* = 0.121 50 u ( z )
z
[3.48]
Alternativamente, para evaluar el factor de friccin se puede utilizar la siguiente

ecuacin:
CD =
z
B + ln 0
[3.49]
Por lo general, se acepta que el perfil logartmico de velocidades (ecuacin [3.23]) se

mantiene en la columna de agua, con lo cual los factores de la ecuacin [3.49] sern k
3.42
(constante de von Karman) = 0.4 y B = 1, dando lugar a una ecuacin ampliamente
utilizada del coeficiente de arrastre:
0.40
CD =
z
1 + ln 0
[3.50]
En aguas profundas, el perfil de velocidades en la columna de agua se obtiene con la

ecuacin [3.24] y el coeficiente de arrastre con la ecuacin [3.49] donde k es la
constante de Von Karman (k = 0.4) y B = ( 2 h ) ln ( 2h ) .
El mtodo escogido depende del tipo de aplicacin. La forma logartmica (ecuacin
[3.49]) es la que tiene la justificacin fsica ms fuerte, sin embargo la forma de le y
potencial es por lo regular ms conveniente para cualquier manipulacin matemtica, y
est mas o menos bien apoyada por datos. En caso de no tener informacin disponible, o
de ser necesaria solo una estimacin aproximada, se sugiere un valor del coeficiente de
arrastre de CD = 0.0025 (Soulsby, 1997).
Debido a que las corrientes pueden tener una variacin vertical compleja, adems de no
ser estacionarias, y a que las tensiones de corte en el fondo estn mas relacionadas con
la corriente cerca del fondo, se puede relacionar 0 (de la ecuacin [3.38]) con la
velocidad de la corriente a 1 m (U100 ) y un coeficiente de arrastre C100 > CD, cuyos
valores dependern del tipo de fondo (tabla 3.4).
b = C100U 100 2
[3.51]
3.6.2.4.- Tensiones de corte totales asociadas a una corriente

La mayor parte de las veces, el lecho sobre el cual fluye una corriente no es plano, sino
que tendr formas geomtricas, como dunas, ripples, etc., que es el caso ms comn
fuera de la zona de rompientes. Es evidente que la resistencia que pueda experimentar el
flujo con el fondo, depender de la geometra del mismo. Sobre fondos no planos, (e.g.
en presencia de ripples) donde el transporte de sedimentos sea limitado, las tensiones de
corte totales estarn compuestas por dos contribuciones: la tensin de corte efectiva bcs
debida a la friccin que acta directamente sobre los granos individuales, y la tensin de
corte de arrastre bcf debido al campo de presin que acta sobre las formas del lecho
(figura 3.14), de manera que:
bc = bcs + bcf
[3.52]
3.43
En el caso de un fondo plano, bcf = 0 y la rugosidad del fondo estara definida por el
tamao de grano, como k s = 2.5 d50 . En este caso, las tensiones de corte efectivas para
un flujo hidrodinmicamente intenso estaran dadas por:
bcs
0.06
2
1
1
2
2
= b = fc U = 12h U
2
2 log

2.5 d 50
[3.53]
donde h es la profundidad del agua, U es la velocidad media de la corriente y el

coeficiente de friccin f c est dado por la ecuacin [3.45]. En el caso de un fondo con
ripples, la componente bcs es la misma (ecuacin [3.53]), pero las tensiones de corte
totales son mayores debidas a la aportacin de la presin ejercida por el flujo sobre las
formas del fondo (donde la escala de rugosidad estar dada por la altura de las formas):
0.06
2
1
2
b =
12 h U
log
2
[3.54]
donde se asume que la rugosidad del lecho es igual a la altura de los ripples r.
Figura 3.14.- Resistencia al flujo debida a un fondo con ripples (Liu, 2001).
La distincin entre bc y bcs explica el fenmeno de que, la aparicin de los ripples en

el fondo (y por lo tanto, un incremento en b ) no necesariamente trae consigo un
incremento en el transporte de sedimento por arrastre de fondo (Liu, 2001).
3.44
Segn Soulsby (1997), la relacin entre bc / bcs es generalmente en el rango de 2 10
para un fondo con ripples. El mecanismo impulsor del movimiento de los granos de
sedimento es bs, de manera que es necesario calcular esta componente cuando se quiere
determinar el umbral del inicio del movimiento, el transporte por arrastre de fondo, etc.
Sin embargo, en el mar, la componente debida nicamente a la friccin es difcil de
calcular cuando hay formas geomtricas presentes. Por lo tanto, generalmente las
tensiones de corte y velocidad de friccin debida a los granos, se calcula como se
describi anteriormente igual que si no hubiera formas geomtricas presentes, y as se
usa en los clculos del transporte de sedimento (como si fuera el autntico valor de la
friccin con los granos).
El valor de la tensin de corte total bc puede calcularse asignando una longitud de
rugosidad total z0 que incluya tanto las componentes de la friccin con los granos, como
del arrastre debido a las formas geomtricas (y la componente debida al transporte de
sedimento, si es necesario). Estos valores se sustituyen en las ecuaciones descritas
anteriormente (ecuaciones [3.46] o [3.50]) para obtener el coeficiente de arrastre total,
las tensiones de corte en el fondo y la velocidad de corte. La velocidad de friccin total
y la longitud de la rugosidad total se pueden obtener a partir de medidas de campo de
los perfiles de velocidad cerca del fondo. Los valores de z0 y C100 (tabla 3.4) dados por
Soulsby (1983) pueden usarse para calcular las tensiones de corte totales.
Material del fondo

z0 (mm)
Lodos
0.2
Lodos / Arena
0.7
Limos / Arena
0.05
Arena (s/ripples)
0.4
Arena (c/ripples)
6
Arena / Conchas
0.3
Arena / Gravas
0.3
Lodos / Arenas / Gravas
0.3
Gravas
3
C100
0.0022
0.0030
0.0016
0.0026
0.0061
0.0024
0.0024
0.0024
0.0047
Tabla 3.4 .- Valores medios de z0 y C100 para distintos materiales del fondo.
A velocidades de flujo altas, cuando se presentan condiciones de sheet flow intenso,

aparece una tercera componente de la rugosidad, que proviene del intercambio de
momentum entre el flujo y el sedimento, que se pone en movimiento. Esta componente
de la escala de rugosidad debida al transporte de sedimento z0t, est relacionada con la
intensidad del transporte, que por su parte, est relacionada con las tensiones de corte
asociadas a los granos del sedimento bcs. Wilson (1989a) encontr una relacin
experimental entre bcs y z0t :
z 0t =
5bcs
30 g (s )
[3.55]
3.45
Las tensiones de corte totales pueden calcularse obteniendo primero la longitud de

rugosidad total z0 sumando las componentes relacionadas al sedimento, a las formas
geomtricas del fondo y al transporte de sedimento:
z 0 = z 0s + z0 f + z0t
[3.56]
(donde z0s est dado en la ecuacin [3.28], z0t por la ecuacin [3.55] y z 0 f = ar ( r2 r )
donde ar es un coeficiente del orden de 1, aunque puede estar entre 0.3 y 3; r y r son
la altura y longitud del ripple respectivamente).
Debe hacerse hincapi en que, debido a la naturaleza no lineal de las ecuaciones, el
valor de b obtenido a partir de calcular z0 (ecuacin [3.56]) y luego la ecuacin [3.50]
ser diferente del valor obtenido calculando por separado las tres componentes las
tensiones de corte, bcs, bcf y bct usando la ecuacin [3.50] para obtener, z0s, z0f y z0t y
posteriormente sumarlos. El primer procedimiento es el correcto.
3.7.- Oleaje
Los efectos de friccin del flujo asociado al oleaje con el fondo producen una capa
lmite oscilatoria dentro de la cual, la magnitud de la velocidad orbital se incrementa
rpidamente con la altura, desde cero en el lecho hasta el valor uw en el borde exterior
de la capa lmite. Para un fondo plano y velocidades orbitales relativamente pequeas,
la capa lmite puede ser laminar, aunque por lo regular, cuando hay sedimento en
movimiento involucrado, la capa lmite ser turbulenta.
En ausencia de una corriente, la turbulencia se ve confinada dentro de la capa lmite,
que en el caso del flujo asociado al oleaje es de tan solo unos cuantos milmetros o
centmetros de grosor, en contraste con la capa lmite de un flujo estacionario que puede
ser del orden de varios (hasta decenas) metros. Como consecuencia de esto, se produce
una velocidad de corte mucho mayor en la capa lmite del oleaje, lo que de hecho
provoca que las tensiones de corte en el fondo debidas al flujo oscilatorio del oleaje
cerca del fondo sean mucho mayores que las producidas por una corriente estacionaria.
Precisamente, al igual que con las corrientes, la propiedad hidrodinmica mas
importante del oleaje en el transporte de sedimentos es la tensin de corte en el fondo
que produce, que en este caso tambin es oscilatoria y cuya amplitud es bwm .
Debido a la naturaleza oscilatoria de las velocidades orbitales, la capa lmite del oleaje
ser no estacionaria. No obstante, es vlido suponer que el grosor de la capa lmite w
debida al oleaje es lo suficientemente pequeo de tal manera que puede calcularse con
teora lineal. Debido a la naturaleza no estacionaria del flujo oscilatorio, es de esperar
que las tensiones de corte en el fondo varen en el tiempo y por lo tanto, la velocidad de
corte sea funcin del tiempo:
3.46
u* =
b
= u* (t )
[3.57]
por lo que es complicado aplicar un modelo simple de viscosidad turbulenta, ya que esta
viscosidad turbulenta del oleaje es una variable dependiente del tiempo. Sin embargo, el
efecto de la dependencia en el tiempo de la viscosidad turbulenta es muy pequeo si se
comparan los resultados obtenidos con la aplicacin de un modelo simple y nodependiente en el tiempo, de la viscosidad turbulenta, en donde la velocidad de corte
est determinada por la mxima tensin de corte en el fondo bwm durante un ciclo de ola
(Towbridge & Madsen, 1984):
u* = u*wm =
bwm
[3.58]
3.7.1.- Perfil de velocidades del oleaje cerca del fondo

De acuerdo con la teora lineal (figura 3.15), el desplazamiento del flujo oscilatorio en
el fondo est dado por la ecuacin:
Abm =
H
1
2 sinh ( kh )
[3.59]
y la velocidad horizontal de la partcula de agua es:
ub = ubwm sin ( t )
[3.60]
donde ubwm es la velocidad horizontal mxima:
ubwm =
H
1
= Abm
T sinh ( kh )
[3.61]
La teora lineal asume que el fluido es ideal (es decir, no viscoso), de tal manera que
hay movimiento de las partculas de agua en el fondo (lo cual no sucede en la realidad).
3.47
Figura 3.15.- Perfil de velocidades horizontales y flujo oscilatorio predicho por la teora lineal de oleaje.
En el estudio del transporte de sedimento, el flujo cerca del fondo es de gran inters,
por lo que nace el concepto de la capa lmite (del oleaje, corrientes e interaccin olascorrientes) en el fondo. El concepto (genrico) de capa lmite fue desarrollado por
Prandtl (1925), que estableci que para los fluidos que tengan viscosidad, el efecto de la
friccin interna en el flujo es apreciable solo en una delgada capa que se genera en las
proximidades de los contornos del flujo. En el caso de los ambientes costeros, la capa
lmite desarrollada bajo la accin del oleaje es de especial inters en los procesos
litorales.
Lundgren & Sorensen (1956) comprobaron experimentalmente el desarrollo de la capa
lmite con un flujo oscilatorio (no realizado con oleaje sino medido en un tnel
oscilante, como se muestra en la figura 3.16) donde el agua era impulsada por un pistn
cuyo movimiento simula el desplazamiento de las partculas del agua en el fondo del
mar predicho por la teora lineal del oleaje. El flujo en la seccin del tnel es horizontal
y uniforme. El grosor de la capa lmite, aunque pequeo, vara con el tiempo, mientras
que fuera de la capa lmite, el flujo permanece imperturbado.
Figura 3.16.- Tnel con flujo oscilante (Lundgren & Sorensen, 1956).
La deduccin de la expresin para la capa lmite la realizan a partir de las ecuaciones

de Navier-Stokes en la direccin horizontal, suponiendo que el flujo en la seccin del
3.48
tnel es horizontal (w = 0) y uniforme ( u x = 0 ) , y que fuera de la capa lmite se
tiene que u = u0 y que = 0. Dado que el flujo es laminar, las tensiones de corte se
pueden expresar por la ley de Newton de la viscosidad (ecuacin [3.3]):
u
z
De aqu, se obtiene la ecuacin que gobierna el movimiento del fluido dentro del tnel:
( u u0 )
2u
= 2
t
z
[3.62]
Imponiendo las condiciones de contorno:

u z =0 = 0
[3.63]
u z = = ubwm sin ( t )
La solucin para la ecuacin [3.62] es:
uw ( z ) = ubwm sin ( t ) ubwm exp
z
t
2 sin

[3.64]
El primer trmino es el de la velocidad horizontal de la partcula de agua (ecuacin

[3.60]) y el segundo trmino es una onda amortiguada, que decae rpidamente lejos del
fondo, como se puede ver en la figura 3.17. En este caso, las tensiones de corte en el
fondo estn dadas por:
bw
u
=
z
U m
= 2 ( sin ( t ) + cos ( t ) )
z =0

= bwm sin ( t + 45 )
[3.65]
3.49
donde bwm
ubwm
=
.
Ya que la velocidad del flujo es u0 = ubwm sin ( t), se puede ver que existe un desfase
entre bw y u0 , es decir, bwm aparecer 45 por delante de ubwm .
Figura 3.17.- Velocidad local del flujo oscilatorio alrededor de u bwm (ecuacin [3.61]).
Por otra parte, de acuerdo con Madsen (1993), la tensin de corte del oleaje obtenida
del modelo simple de viscosidad turbulenta, es decir, sin considerar que la velocidad de
corte sea dependiente del tiempo (ecuacin [3.57]) sino funcin de la tensin de corte
mxima (ecuacin [3.58]), es:
bwm = t
uw
u
= k u*w m z w
z
z
[3.66]
donde uw es la velocidad orbital del oleaje.

Usando la ecuacin anterior, aplicando teora lineal del oleaje a las condiciones de
contorno y simplificando, se obtiene la solucin para el perfil de velocidades (Madsen,
1993):
uw ( z ) =
2
z
sin ubwm ln cos(t + )
z0
[3.67]
3.50
donde ubwm es la velocidad orbital mxima cerca del fondo y es el desfase de la
velocidad orbital del oleaje (con respecto de la superficie libre) cerca del fondo, dada
por:
tan =
2
ku
ln *w m 1.15
z0
[3.68]
Cabe mencionar que el perfil de velocidades cerca del fondo para un flujo de capa
lmite de una onda peridica es logartmico, igual que en el caso de una corriente
estacionaria. Esto sugiere que hay un cierto grado de similitud entre los flujos
estacionarios y los no estacionarios en la capa lmite turbulenta inmediatamente por
encima del fondo (por lo tanto, es vlido usar la condicin de friccin deducida
experimentalmente para corrientes estacionarias).
3.7.2.- Grosor de la capa lmite del oleaje

Con la inclusin del concepto de la capa lmite del oleaje, el perfil de velocidades
correspondiente a la cresta del oleaje se muestra en la figura 3.18
Figura 3.18.- Capa lmite del oleaje para flujo laminar sobre un fondo no rugoso (Liu, 2001).
El grosor de la capa lmite depende de como se defina el lmite superior de la capa

lmite:
3.51
i.
Jonsson (1966) defini el lmite superior de la capa lmite como la elevacin

mnima donde la magnitud de la velocidad es igual a ubwm , que de acuerdo con la
ecuacin [3.64], es de:
j =
ii.
2
2
Sleath (1987) define el lmite superior de la capa lmite como la elevacin donde
la magnitud de la velocidad es el 95% de ubwm , es decir, existe una diferencia
relativa del 5%, lo que da un grosor de la capa lmite de:
0.05 = 3
iii.
[3.69]
[3.70]
Madsen, (1993) propone que a partir de la aproximacin logartmica dada por la

ecuacin [3.67] (nicamente vlida solo cerca del fondo y para valores altos de
(k u*wm / z0 ), la posicin en la cual |ubw| = ubwm es z (k u*wm / ). A pesar
que el uso de la ecuacin [3.67] para la prediccin del grosor de la capa lmite
excede su validez, diversos resultados (Grant & Madsen, 1986) sugieren que:
w =
ku*wm
=k
fw
A
2 bm
[3.71]
En la realidad, el flujo sobre un fondo con ripples es turbulento, por lo que el grosor de
la capa lmite se ve afectado por la rugosidad del fondo (figura 3.19). Por esto, Sleath
(1987) propone la siguiente ecuacin emprica para obtener el grosor de la capa lmite
en funcin de la rugosidad del fondo:
A
0.05
= 0.26 bm
kn
kn
0.70
[3.72]
donde k n es la rugosidad del fondo y Abm es la amplitud de la excursin de la partcula

del agua en el fondo (ecuacin [3.59]).
As, conociendo las caractersticas del oleaje (altura y longitud), la profundidad del
agua y la altura de los ripples, es posible obtener el grosor de la capa lmite del oleaje,
encontrando la amplitud del desplazamiento de la partcula del agua (asociada al flujo
oscilatorio) en el fondo (ecuacin [3.59]) y aplicando las expresiones de Jonsson (1966)
3.52
o Sleath (1987), utilizando la rugosidad del fondo como la altura de los ripples. Cabe
recalcar que el grosor de la capa lmite del oleaje vara con el tiempo.
Figura 3.19.- Grosor de la capa lmite del oleaje sobre un fondo con ripples (Liu, 2001).
El mtodo descrito aqu es para determinar su anchura en el paso de la cresta de la ola,

es decir, el grosor mximo de dicha capa. Si bien la capa lmite generada por una
corriente puede cubrir el total de la profundidad, la capa lmite generada por el flujo
oscilatorio del oleaje es muy delgada debido precisamente a la oscilacin de las
partculas de agua, ya que en medio periodo del oleaje, las partculas cerca del fondo se
mueven en el sentido de propagacin del oleaje (desarrollando una capa lmite). Justo a
la mitad del periodo del oleaje, las partculas del agua en el fondo se detienen por
completo para comenzar a moverse en direccin contraria al avance del oleaje, por lo
que la capa lmite comienza a desarrollarse otra vez desde el principio.
3.7.2.1.- Factor de friccin del oleaje (f w) en fondos no rugosos

De manera similar a la definicin del coeficiente de friccin debido a las corrientes
(ecuaciones [3.32] y [3.53]), Jonsson (1966) define el coeficiente de friccin del oleaje
(tomando en cuenta que las tensiones de corte en el fondo debidas al oleaje son funcin
del tiempo) como:
bwm =
1
f w ubwm2
2
[3.73]
donde ubwm es la velocidad horizontal mxima de la partcula de agua en el fondo, dada

por la teora lineal del oleaje y f w es el coeficiente de friccin del oleaje. Es importante
recalcar que bwm y ubwm no ocurren al mismo tiempo.
De forma anloga a la velocidad de corte debida a las corrientes (ecuacin [3.34]), la
velocidad de friccin del oleaje se define como:
3.53
u* w =
bwm
=
fw
ubwm
2
[3.74]
donde f w es el coeficiente de friccin del oleaje. Esta expresin es idntica a la del

coeficiente de friccin de la corriente f c.
El coeficiente de friccin del oleaje (para una capa lmite laminar y sobre un fondos no
rugosos) se obtiene de la expresin terica para la tensin de corte mxima en el fondo
bwm para una capa lmite laminar de oleaje y un fondo no rugoso es
bwm
ubwm
[3.75]
combinando las ecuaciones [3.75] y [3.73], se obtiene una expresin del coeficiente de
friccin:
2
f w = ubwm
0.5
2

= A
= 2
2 bm
2
Abm

[3.76]
De acuerdo con resultados experimentales, se ha encontrado que la capa lmite laminar

en fondos no rugosos corresponde a:
Abm2
< 3x105
[3.77]
Por otra parte, para obtener el coeficiente de friccin del oleaje (para una capa lmite
turbulenta y sobre fondos no rugosos), Justensen (1988) sugiere que:
0.123

f w = 0.024
2
Abm
[3.78]
para
106 <
Abm2
< 108
[3.79]
3.54
3.7.2.2.- Factor de friccin del oleaje (f w) en fondos rugosos
En la naturaleza, los flujos son en su mayora, turbulentos y sobre fondos rugosos. Por
lo tanto, las expresiones descritas anteriormente no suelen ser aplicadas en problemas
costeros. Jonsson, (1966) propone una ecuacin emprica implcita para el coeficiente
de friccin para un flujo turbulento sobre un fondo rugoso, que fue posteriormente
aproximada de forma explcita por Swart (1974):
0.194
kn
f w = exp 5.213
5.977
Abm
[3.80]
Nielsen (1992) sostiene que la ecuacin de Swart (1974) tiende a sobreestimar el

coeficiente de friccin f w para valores pequeos de k n / Abm (figura 3.20), y propone una
expresin que proporciona un mejor ajuste:
k 0.2
f w = exp 5.5 n 6.3

Abm
[3.81]
Con estas expresiones es posible obtener el coeficiente de factor del oleaje a partir de la
altura y longitud del oleaje, la profundidad y la altura de los ripples, encontrando la
amplitud del desplazamiento de la partcula en el fondo (ecuacin [3.59]) y evaluar
directamente la ecuacin [3.81].
Figura 3.20.- Factores de friccin del oleaje observados para flujos turbulentos (Nielsen, 1992).
3.55
Por otra parte, Madsen (1993) sugiere que a partir del perfil de velocidades cerca del
fondo uw (z) (ecuacin [3.67]), la tensin de corte en el fondo puede obtenerse de
bw = bwm cos( t + )
[3.82]
en trminos del factor de friccin del oleaje f w, definido por Jonsson (1966), la mxima
tensin de corte en el fondo (para el caso de un fondo plano) est dado por la
ecuacin[3.73]:
bwm =
1
2
f w ubwm
2
[3.83]
que, en trminos de la velocidad de corte es:
u*wm =
bwm
=
fw
ubwm
[3.84]
Este es generalmente el caso para la zona de rompientes, ya que por lo regular, el flujo
es tan intenso que no permite la existencia de ripples. En este caso, la magnitud de la
tensin de corte en el fondo es la debida a la friccin con el sedimento.
Usando la solucin del modelo de Madsen (1993) del perfil de velocidades cerca del
fondo (ecuaciones [3.66] a [3.68]) y las tensiones de corte y factor de friccin del oleaje
derivadas de dicho perfil (ecuaciones [3.82] a [3.84]), se obtiene una ecuacin implcita
para el factor de friccin del oleaje:
f
2
2
k w ubwm

=
+
2
1.15 2
fw
ln
z0
[3.85]
donde k = 0.4 y la amplitud de la trayectoria de la partcula de agua predicha por la

teora lineal:
Abm =
ubwm
[3.86]
3.56
As, la ecuacin [3.85] puede aproximarse por una relacin implcita del factor de
friccin del fondo propuesta por Grant y Madsen (1986) para flujos altamente
turbulentos (cuando k n =30 z0 ):
1
4 fw
+ log10
1
4 fw
= log10
Abm
0.17 + 0.24 4 f w
kn
Si el flujo es suavemente turbulento, k n se sustituye por 30 z 0 =
[3.87]
3 .3
y la ecuacin
u*wm
[3.87] puede rescribirse como:
1
1
Re
+ log10
= log10
0.17 + 0.06 4 4 f w
50
4 4 fw
4 4 fw
[3.88]
donde
Re w =
ubwm Abm
[3.89]
es un nmero de Reynolds del oleaje. Finalmente, se puede usar el factor de friccin del
oleaje para una capa lmite laminar (Jonsson, 1966):
fw =
2
Re w
[3.90]
3.7.2.3.- Evaluacin del factor de friccin

Myrhaug (1989) propone una expresin implcita para f w que hace uso de la ecuacin
[3.26] y es vlida para todos los regmenes turbulentos del oleaje:
1
2
0.32
Re
f
= ln(6.36 r f w 2 ) ln 1 exp( 0.0262 w w
fw
r
+ 4.71 r
1
Re f 2
w w
+ 1.64 [3.91]
3.57
donde Rew es el nmero de Reynolds del oleaje (ubw Abm / ) y r es la rugosidad relativa
del fondo (Abm / kn ).
Se han definido una gran cantidad de expresiones para definir el factor de friccin del
oleaje f wr para flujos altamente turbulentos, entre las que se destacan las siguientes:
Swart (1974):
f wr = 0.3 para r 1.57

f wr = 0.00253exp(5.21 r 0.19 ) para r > 1.57
[3.92]a, b
Nielsen (1992), para todos los posibles valores de r:

f wr = exp( 5.5r 0. 2 6.3)
[3.93]
Soulsby (1993), para todos los posibles valores de r:

0.52
A
f wr = 1.39 bm
z0
que en trminos de z0 = ks 30
[3.94]a, b
f wr = 0.237r 0.52
El factor de friccin para fondos no-rugosos f ws se puede obtener a partir de la siguiente

ecuacin (Soulsby, 1997):
f ws = BRw
[3.95]
donde:
B = 2 y N = 0.5 para Rew 5 x 105 (flujo laminar).
B = 0.0521 y N = 0.187 para Rew > 5 x 105 (flujo suavemente turbulento) (Soulsby,
1997).
Valores alternativos para B y N:
B = 0.035 y N = 0.16 para Rew > 5 x 105 (flujo suavemente turbulento) (Fredse &
Deigaard, 1992).
B = 0.0450 y N = 0.175 para Rew > 5 x 105 (flujo suavemente turbulento) (Myrhaug,
1995).
3.58
Por su parte, Madsen (1993) propone que, conociendo la rugosidad equivalente k n , y las
condiciones del oleaje descritas por Abm y ubwm , es posible calcular el factor de friccin
del oleaje f w por medio de 3 ecuaciones distintas. El f w que se debe seleccionar es el
mayor de los 3 valores. En este contexto cabe mencionar que en el calculo de f w para
flujos altamente turbulentos o suavemente turbulentos, los resultados predichos por las
ecuaciones [3.88] y [3.90] son idnticos para Re 3 x 104 , por lo cual, segn Madsen
(1993), puede considerarse un valor razonable para la transicin del rgimen del flujo,
de laminar (Re < 3 x 104 ) a turbulento.
La evaluacin de los factores de friccin del oleaje de las ecuaciones [3.87] y [3.88]
puede realizarse mediante un proceso iterativo, que para la ecuacin [3.87] es como
sigue:
1
x
( n +1 )
A
1
= log 10 bm 0.17 log 10 ( n ) + 0.24 x ( n)
kn
x
[3.96]
donde el suprandice denota el paso iterativo x = 4 f w y la iteracin comienza desde

x (0) = 0.4.
Aqu se supone que el perfil de velocidades orbitales dentro de la capa lmite puede
expresarse por medio de la ecuacin [3.67] y que esta expresin tan solo es vlida para
valores pequeos de z 0 ( k u*wm ) = 0.12
fw
)(k
Abm ) , es decir, para valores
relativamente grandes de Abm k n . Por esta razn y debido a que la solucin iterativa de
la expresin del factor de friccin del oleaje (ecuacin [3.87]) converge lentamente para
Abm kn < 100 , es posible obtener una solucin exacta para f w y (cuando Abm kn < 100 )
a partir de la figura 3.21a,b.
Para capas limite de flujos moderadamente turbulentos, la ecuacin [3.88] puede
expresarse de una forma similar a la ecuacin [3.96] y la solucin iterativa se realiza
como se ha sealado anteriormente, con x = 4 4 f w .
3.59
Figura 3.21.- Diagrama del factor de friccin del oleaje (a) y ngulo de fase de la friccin del fondo (b)
(CEM , 1998).
3.7.2.4.- Tensiones de corte totales del oleaje

Los fondos en zonas de profundidades intermedias y someras (con excepcin de la zona
de rompientes) estn generalmente formados por ripples, que pueden haberse generado
por la accin del oleaje principalmente y mas ocasionalmente por las corrientes.
Al igual que en el caso de las corrientes sobre fondos no-planos con transporte de
sedimento limitado, las tensiones totales inducidas por el oleaje bw estn compuestas
por la debida a la friccin con la superficie del fondo bws y por la fuerza de arrastre
debida a las formas geomtricas bwf :
bw = bws + bwf
[3.97]
Se han propuestos varios mtodos para calcular las dimensiones de la rugosidad de los
ripples generados por la accin del oleaje que, por lo regular, involucran un valor de z0
que para ser calculado puede necesitar los trminos correspondientes a los efectos de un
transporte de sedimentos intenso, as como los trminos de la rugosidad por friccin de
superficie (o granular) y la rugosidad por arrastre (asociado a las formas geomtricas).
La dificultad en el clculo de estas componentes es igual que en el caso de las
corrientes, por lo que lo usual es aproximarlas a la tensin de corte de superficie o
granular bws obtenida usando el valor de z 0 = d50 12 en los clculos de las fricciones.
La razn por la que se realiza de esta manera es que bws puede ser calculada de una
forma relativamente simple y sin ambigedades y puede servir como variable
independiente a otras variables dependientes tales como las dimensiones de los ripples o
las concentraciones de referencia del sedimento, que se pueden relacionar con bws.
3.60
3.7.3.- Perfil de velocidades del oleaje fuera de la capa lmite
Teniendo en cuenta la limitacin de la validez del perfil de velocidades dada por la
ecuacin [3.67], de la ecuacin [3.85] se deduce que el perfil de velocidades dentro de
la capa lmite del oleaje puede expresarse de forma simple, con un lmite superior de
validez dado por el valor de z en el cual |uw| = ubm o aproximadamente como z = w/
(Madsen, 1993), es decir:
uw ( z ) =
fw
u
2 bwm ln z cos(t + ) = u*wm ln z cos( t + )
k
z0
k
z0
[3.98]
En la figura 3.22 se comparan las amplitudes de las velocidades orbitales dentro de la

capa lmite obtenidas experimentalmente (Jons son & Carlson, 1976) con la prediccin
obtenida a partir de la teora de Madsen (1993) presentada aqu, y por la solucin
aproximada de la ecuacin [3.98]. Claramente, el perfil aproximado de velocidades es
muy similar al perfil (cerca del fondo) exacto y al perfil experimental. An cuando hay
algunas diferencias entre los perfiles por encima del fondo, la solucin aproximada
proporciona una representacin adecuada del perfil medio, en particular cuando se
proyecta su uso para la prediccin de las tasas de transporte de sedimento.
Figura 3.22.- Comparacin entre el modelo de la prediccin de la velocidad orbital del oleaje dentro de la
capa lmite del oleaje (ec. [3.98]) con medidas realizadas por Jonsson & Carlson (1976) (CEM, 1998).
3.61
El anlisis descrito de la mecnica de la capa lmite se basa en la suposicin de ondas

simples y peridicas. En la naturaleza, el oleaje irregular (real) se puede representar por
la superposicin de varias ondas individuales, con distintas alturas, frecuencias y
direcciones de propagacin, es decir, por un espectro direccional de oleaje. Madsen,
Poon y Graber (1988) presentaron un modelo de capa lmite para olas descritas por su
espectro de velocidades orbitales cerca del fondo. En este anlisis se utiliza una
viscosidad turbulenta invariante en el tiempo, similar a la discutida anteriormente para
ondas monocromticas y las velocidades orbitales cercanas al fondo, ubwm,i, y la
frecuencia i de la i- sima componente del espectro del oleaje para calcular un valor
representativo de las caractersticas del oleaje. As por ejemplo, la velocidad orbital
mxima en el fondo es
ubwr = ubwm =
2
b w m, i
[3.99]
y la frecuencia (y por ende, el periodo) representativo es
( u
==
u
i
2
b w m, i
[3.100]
b w m, i
Aqu ubwm,i2 y i son la velocidad orbital mxima cerca del fondo y la frecuencia de la
i-sima componente del espectro.
Estas caractersticas de una onda representativa pueden ser usadas junto con las
ecuaciones derivadas anteriormente para olas peridicas, para obtener las caractersticas
ms importantes de flujos en capa lmite turbulenta para olas espectrales. Dicho de otra
forma, el uso de las condiciones representativas ubwr = ubwm y r = , junto con la
informacin de rugosidad equivalente del fondo, k n , permitir la determinacin de un
factor de friccin representativo, f wr = f w, de las ecuaciones [3.87] o [3.88]. Con estos
valores del factor de friccin del oleaje, las tensiones de corte mximas en el fondo
asociadas con la componente espectral estarn dadas por:
bwmi, =
1
f wrubwrub w m, i
2
[3.101]
y para calcular el perfil de velocidades cerca del fondo (dentro de la capa lmite del
oleaje) de las componentes individuales del oleaje ubw,i (ecuacin [3.98]), el trmino
ubwm se sustituye por ubwm,i.
Es importante destacar que en este anlisis, la caracterstica representativa del oleaje
descritas por Madsen (1993) no es la altura de ola significante, sino la altura media
3.62
cuadrtica. Por lo tanto, si las olas se describen mediante su altura y periodo significante
Hs y Ts debe obtenerse la altura de ola media cuadrtica H rms = H s 2 y T = Ts. para
as obtener la velocidad orbital representativa cercana al fondo.
3.8.- Flujos combinados de olas y corrientes

Dada la naturaleza de la interaccin del flujo con el fondo, donde las corrientes generan
una capa lmite mas grande (que puede incluso abarcar la totalidad de la profundidad)
que el oleaje, se desprende que la friccin ejercida por las corrientes (ms estacionarias)
es mucho menor que la friccin que ejerce el flujo oscilatorio del oleaje. De tal manera
que en las regiones costeras, donde coexisten olas y corrientes, an cuando la velocidad
de la corriente sea mucho mayor, que la velocidad cerca del fondo inducida por el flujo
oscilatorio, ser sta la que domine la situacin en las proximidades del fondo (figura
3.23). Esto trae como consecuencia que el oleaje pueda poner en suspensin una mayor
cantidad que las corrientes, aunque por su naturaleza oscilatoria, la capacidad de
transporte es casi nula (salvo en la zona de rompientes), siendo las corrientes las que
cumplen con esa funcin. Por lo tanto, en la gran mayora de las reas costeras, los
flujos combinados de olas y corrientes juegan un papel protagonista en la dinmica de
los sedimentos. La mayor dificultad en el estudio de estos temas es que las olas y las
corrientes interactan hidrodinmicamente de forma no lineal, de manera que los flujos
resultantes no son una suma lineal de sus comportamientos. Las principales
interacciones son:
Modificacin de la velocidad de fase y longitud de onda en su interaccin con

una corriente, dando lugar a una refraccin importante.
Interaccin de las capas limite de olas y corrientes, dando lugar a un incremento
en las tensiones de corte de ambas componentes (la estacionaria y la oscilatoria).
Generacin de corrientes costeras inducidas por la rotura del oleaje.
Figura 3.23.- Comparacin entre perfiles verticales de velocidades de oleaje y corrientes (Liu, 2001).
Cabe notar sin embargo, que en las situaciones donde se presenten flujos combinados
de olas y corrientes, la propagacin del oleaje y las transformaciones de la longitud de
3.63
onda y velocidades orbitales requieren del uso de la relacin de dispersin de olas en
presencia de corrientes, pero al hacer los clculos del comportamiento de la capa lmite,
es la frecuencia absoluta y no la relativa (en el caso de la interaccin con una corriente)
la que debe usarse (Soulsby, 1997).
Las tensiones de corte en el fondo bajo flujos combinados de olas y corrientes se ven
fortalecidas con respecto a los valores que resultaran de una simple suma de sus valores
individuales, ya que su interaccin es no lineal. Se han propuesto una gran cantidad de
teoras y modelos que describen este proceso. Soulsby (1997) realiza una comparacin
entre los valores predichos de las tensiones de corte medias bwc y mximas bwcm
durante un ciclo de ola, donde se muestra que las diferencias entre los modelos son del
orden de un 30 a 40%, y hasta de un factor de 3 para casos cuyas condiciones estaban
fuertemente dominadas por el oleaje. Soulsby (1995) realiza una serie de comparaciones
entre modelos parametrizados contra un gran nmero de datos de bwc obtenidos en
laboratorio y en campo, usando varios criterios para mejorar la bondad de los ajustes,
encontrando que ninguno de los resultados de los modelos comparados resultaba mejor,
de acuerdo con todos los criterios. En el caso de los modelos analticos, sugiere que el
de Grant & Madsen (1979) y el Fredse (1984) son de los modelos analticos que, en
trminos generales, ofrecen mejores resultados.
Por lo tanto, es posible encontrar las tensiones de corte del oleaje y las corrientes,
conociendo la altura y longitud de onda del oleaje, la profundidad y tamao de los
ripples, encontrando la amplitud del desplazamiento de la partcula en el fondo
(ecuacin [3.59]), f w (ecuacin [3.81]), f c (ecuacin [3.45]), ubwm (ecuacin [3.61]) y
evaluando bc (ecuacin [3.32]) y bwm (ecuacin [3.73]).
3.8.1.- Modelo de Madsen

El modelo conceptual de Madsen (1993) de la turbulencia ge nerada cerca del fondo,
sugiere que la viscosidad turbulenta en las inmediaciones del fondo es funcin del
tiempo siempre que el flujo sea no-estacionario. Sin embargo, en el caso de un oleaje
puro, es posible suponer una viscosidad turbulenta invariante (no dependiente) en el
tiempo, basada en la velocidad de corte mxima. Esto da como resultado una gran
simplificacin en los anlisis, sin sacrificar apreciablemente la exactitud de los clculos
(Towbridge & Madsen, 1984). En el caso de las capas lmite de flujos combinados de
olas y corrientes, la magnitud de la viscosidad turbulenta en las inmediaciones del fondo
ser funcin de la mxima velocidad de corte en dicho flujo. Ya que la extensin
vertical de la turbulencia asociada a las olas est limitada por el grosor de la capa lmite
del oleaje, la contribucin del oleaje a la mezcla turbulenta debe desvanecerse en algn
nivel en el cual solo la velocidad de corte de la corriente contribuir a la mezcla
turbulenta.
Para el caso general del flujo combinado de olas y corrientes, con la corriente fluyendo
a un cierto ngulo wc con respecto a la direccin de propagacin del oleaje, la mxima
tensin de corte combinada en el fondo bwcm puede expresarse como:
3.64
bwcm = bwm bc =
= bwm
bwm
bc cos wc
1 + 2 bc cos wc + bc
bwm
bwm
) +(
2
bc
sin wc )
[3.102]
lo que muestra que bwcm bwm ya que es de esperar que bc << bwm . Obteniendo bwcm a
partir de la ecuacin [3.102]) y las tensiones de corte en el fondo debidas a las
corrientes bc (ecuaciones [3.32] o [3.33]), el modelo de viscosidad turbulenta para
flujos combinados de olas y corrientes de Grant y Madsen (1979, 1986):
k u*wm z para z < cw

t =
k u*c z para z > cw
[3.103]
bwm
, u*c = bc y cw es el nivel al cual la turbulencia debida al oleaje
se desvanece, es decir, el grosor de la capa lmite del oleaje. Es interesante notar que
para este caso, hay una discontinuidad en z = cw. Para evitar esta inconsistencia, se
puede imponer una transicin suave entre los dos segmentos no lineales. No obstante, la
precisin lograda en los clculos no justifica lo engorroso y complicado que se vuelven
stos clculos, exactamente como ocurre con la consideracin de una viscosidad
turbulenta t variable en el tiempo. Por lo tanto, la utilizacin del modelo de t de Grant
& Madsen (1979, 1986) es perfectamente justificable.
donde u*wm =
Debido a que la capa lmite para flujos combinados de olas y corrientes interactuando
con un cierto ngulo es tridimensional, el problema debe, en principio, resolverse para
las 2 componentes del flujo horizontal. Sin embargo, cuando se escoge la direccin x
como la direccin de propagacin, slo en esta direccin habr flujo no estacionario. En
este caso, se puede demostrar que todas las formulaciones derivadas para el caso del
oleaje puro son vlidas tambin para oleaje en presencia de una corriente cuando u*wm
(obtenida de la ecuacin [3.102] con u*wm = bwm ) es reemplazada por u*wcm . Como
se ha mencionado antes, u*wm u*wcm , lo que indica que el movimiento del oleaje est
muy poco influenciado por la presencia de una corriente.
El flujo estacionario estar completamente confinado a la direccin de la corriente, es
decir, a un ngulo wc con respecto a la direccin del oleaje. Tomando en cuenta la ley
de la pared, la velocidad de la corriente se obtiene de
bc
= u*c 2
uc
ku*wcm z z para z < cw

=
ku z uc para z >
cw
*c z
[3.104]
3.65
que para z < cw, tiene como solucin:
uc ( z ) =
u*c u*c
z
ln
k u*wcm z 0
[3.105]
donde z0 est dado por la ecuacin [3.19] con u*wcm sustituyendo u*

Para z > cw, la solucin es:
uc ( z ) =
u*c
z
ln
k
z 0a
[3.106]
donde z0a es la rugosidad aparente del fondo (obtenida como la constante de integracin
arbitraria). Igualando las dos soluciones en z = cw e introduciendo la expresin para z0a
en la ecuacin [3.106], se obtiene una solucin alternativa para z > cw
uc ( z ) =
u*c
z
u

+ *c ln cw
ln
k cw u*wcm z 0
[3.107]
De aqu se desprende el importante efecto que las olas pueden tener en los perfiles de
velocidad de una corriente. En primer lugar, los gradientes de velocidades dentro de la
capa lmite para flujos combinados de olas y corrientes, se reducen en un factor de
u*c u*wcm con respecto a su valor en ausencia de olas. Esto es consecuencia del
incremento en la intensidad de la turbulencia dentro de la capa lmite del oleaje que
proviene precisamente de su accin.
En el caso (extremo) de que u*c u*wcm 0 , entonces el valor de uc ser
aproximadamente cero a lo largo de la capa lmite del oleaje. Esto significa que la
velocidad de la corriente fuera de la capa lmite del oleaje debe alcanzar un valor
cercano a cero a una profundidad z = cw, es decir, cw z0a, que es el valor de la
rugosidad aparente. Esta misma conclusin puede obtenerse de la ecuacin [3.107] en el
lmite de u*c u*wcm 0 . Por lo tanto, las corrientes en presencia del oleaje
experimentarn una mayor rugosidad de fondo (es decir, para una corriente con la
misma velocidad u (z) a un nivel dado por encima del fondo, la velocidad de corte de la
corriente y las tensiones de corte de dicha corriente se incrementarn en presencia del
oleaje).
Si se compara el perfil de velocidades de una corriente en presencia del oleaje, predicha
por esta teora de interaccin olas-corrientes con un perfil medio de velocidades de una
3.66
corriente (Bakker & Van Doorn, 1978), se concluye que la definicin de cw dada por la
ecuacin [3.71] puede adoptarse para la teora de interaccin olas-corrientes. Esta teora
tambin es aplicable para un oleaje descrito por su espectro de velocidades orbitales
cerca del fondo, usando los parmetros representativos del tren de oleaje, con la
direccin de propagacin del oleaje tomada como la direccin pico (Madsen, 1993).
3.8.1.1.- Perfil de velocidades para flujos combinados de olas y corrientes

Imponiendo como marco de referencia que la direccin de propagacin del oleaje se
defina como el eje x, el perfil de velocidades justo por encima del fondo est definido
por la ecuacin [3.67], donde el ngulo de fase definido por la ecuacin [3.68] se evala
con u*wm = u*wcm.
El vector velocidad de la corriente est dado por:
uc = u c ( z ){cos wc , sin wc }
[3.108]
donde uc (z) est dado por las ecuaciones [3.105] y [3.107].
3.8.1.2.- Tensiones de corte para flujos combinados de olas y corrientes
De forma similar a lo expuesto anteriormente, la tensin de corte asociada a un flujo no

estacionario combinado de olas y corrientes puede obtenerse de:
bwc ( t ) = {bwm cos(t + ) + bc cos wc, bc sin wc }
[3.109]
Para la aplicacin de esta teora para la interaccin de olas y corrientes para oleaje
irregular, se usa el espectro de velocidades orbitales cerca del fondo, definido por su
velocidad orbital mxima en el fondo (ecuacin [3.99]), es decir
ubwr = ubwm =
2
b w m, i
y su periodo representativo, obtenido a partir de la frecuencia
representativa (ecuacin [3.100]) = r = ( i ubwmi, 2 )
2
b w m, i
3.8.2.- Modelo de Fredse

Para el modelado del transporte de sedimentos, es necesario estimar las tensiones de
corte en el fondo y el perfil de velocidades bajo flujos combinados de olas y corrientes.
3.67
Hay varios modelos en el estado del arte que describen estos flujos, siendo uno de los
ms importantes el modelo de Fredse (1981). Este modelo sostiene que el grosor de la
capa lmite del oleaje es , y que el flujo se divide en dos zonas, fuera de la capa lmite
(z < ) donde el flujo no se ve sujeto a la friccin y dentro de la capa lmite (z > ). Con
la superposicin de una corriente dbil, se produce una cierta cantidad de turbulencia
fuera de la capa lmite del oleaje. Dentro de la capa lmite del oleaje, la turbulencia es
debida a la aportacin tanto del oleaje como de las corrientes. Sin embargo, la corriente
es tan dbil que el grosor de la capa lmite del oleaje es aproximadamente el mismo.
3.8.2.1.- Caso 1.- El oleaje y la corriente se propagan en la misma direccin

En el caso de oleaje puro, la mxima tensin de corte en el fondo est dada por la
ecuacin [3.73] y la velocidad de corte [3.74] (utilizando ubwm que es la velocidad
mxima en el fondo obtenida con teora lineal). Ya que la capa lmite del oleaje es muy
delgada, ubwm puede ser considerada como la velocidad en el borde superior de la capa
lmite. La velocidad instantnea de la capa lmite y las tensiones de corte en el fondo
estn dadas por las ecuaciones [3.60] y [3.73] (modificada, es decir, cambiando ubwm por
ubw) respectivamente.
Superponiendo una corriente, cuya magnitud en el borde superior de la capa lmite es
u, la velocidad instantnea del flujo combinado es:
uwc = ( ubw + u ) = ( ubwm sin ( t ) + u )
[3.110]
donde las tensiones de corte combinadas son:
1
bwc = f w uwc 2
2
[3.111]
y las tensiones de corte medias son:
wc =
1
T
bwc dt
2
f w ubwm u
[3.112]
As, conociendo las caractersticas del oleaje (H, T, h), y la velocidad media de la
corriente (U) se pueden calcular , f w y Um (con los datos del oleaje), pero u no es
inmediato, ya que el perfil de velocidades de la corriente se ve afectado (distorsionado)
por el oleaje.
3.68
3.8.2.1.1.- Perfil de velocidades fuera de la capa lmite del oleaje
El perfil de velocidades de la corriente (sin oleaje) se puede describir con una
combinacin de las ecuaciones [3.18] y [3.19]:
u (z ) =
u*c
z
ln
k 0.033 kn
[3.113]
donde u*,c es la velocidad de corte de la corriente, k n es la rugosidad del fondo (que para
flujos altamente turbulentos es z0 = 0.033 k n ) y k es la constante de von Karman (= 0.4).
En presencia del oleaje, Grant & Madsen (1979) sugieren que el perfil de velocidades
es:
u (z ) =
u*wc
z
ln
k
0.033 knw
[3.114]
donde u*wc es la velocidad de corte del flujo combinado de olas y corrientes, y est dado
por la ecuacin:
u*wc =
bwc
2
f w ubwm u
[3.115]
dnde k nw es la rugosidad del fondo bajo la accin combinada de olas y corrientes.

Aplicando la condicin de contorno u z= = u en la ecuacin [3.114] se obtiene una
expresin para la rugosidad del fondo k nw:
ku
k nw = 30 exp
u*wc
[3.116]
El perfil promediado de velocidades es:

z
ln
dz
k nw
1 h
u*wc 6.2 + ln
k knw
U=
1 h
u*wc
u
(
z
)
dz
=
h
kh
[3.117]
3.69
Combinando las ecuaciones [3.115] y [3.116] y [3.117] se obtiene una expresin

explcita para u :
u = C C 2 U 2
[3.118]
donde
1 h
C = U + f w ubwm 6.2 + ln
k 30
[3.119]
3.8.2.1.2.- Perfil de velocidades dentro de la capa lmite del oleaje

Dentro de la capa lmite del oleaje, la turbulencia es debida a la aportacin del oleaje y
de las corrientes. La viscosidad turbulenta combinada es:
twc = ( tc + tw ) = k ( u*wc + u*w ) z
[3.120]
Por lo tanto, (anlogamente a la ecuacin [3.11]) la tensin media de corte en el fondo

es:
du
bwc =twc wc
dz
[3.121]
Combinando las ecuaciones [3.120] y [3.121], se obtiene:
duwc 1 u*wc 2
=
dz
z u*wc + u*w
[3.122]
Integrando la ecuacin [3.122] y resolviendo para uwc :
z
1 u*wc 2
uwc ( z ) =
ln
k u*wc + u*w z0
[3.123]
3.70
donde z0 es una constante de integracin correspondiente a la altura donde la velocidad
es nula, es decir u z= z = 0 , que de acuerdo con Nikuradse, z0 = 0.033 kn (ecuacin
0
[3.19]). La figura 3.24 muestra un ejemplo del perfil de velocidades con y sin oleaje.
De esta manera, conociendo la velocidad media de una corriente y las caractersticas
del oleaje (H, T, h) y la altura de los ripples, ambos viajando en la misma direccin, es
posible calcular las tensiones de corte y el perfil de velocidades en flujos combinados de
olas y corrientes, calculando la longitud de onda con teora lineal, la amplitud del
desplazamiento de la partcula de agua en el fondo (ecuacin [3.59]) y la velocidad
mxima del flujo oscilatorio (ecuacin [3.61]). Considerando la escala de rugosidad del
fondo como la altura de los ripples, se calcula el grosor de la capa lmite (ecuacin
[3.72]), el factor de friccin del oleaje f w (ecuacin [3.81]), la velocidad de friccin u*,w
(ecuacin [3.74]) y la velocidad de la corriente en el borde superior de la capa lmite
(ecuaciones [3.118] y [3.119]). Con esto se obtiene la tensin de corte media en el fondo
wc (ecuacin [3.112]), la velocidad de corte asociada (ecuacin [3.115]), la rugosidad
bajo el flujo combinado de olas y corrientes k nw (ecuacin [3.116]) y finalmente el perfil
de la corriente dentro de la capa lmite (ecuacin [3.114]) y el perfil de velocidades
fuera de la capa lmite (ecuacin [3.123], con z0 = 0.033 k n .).
Si solo se tratara de una corriente, el procedimiento sera calcular el factor de friccin
de la corriente f c (ecuacin [3.45] ), la velocidad de corte asociada u*c (ecuacin [3.34])
y finalmente el perfil de velocidades (ecuacin [3.113]).
Figura 3.24.- Perfil de velocidades con y sin oleaje para una misma corriente (Liu, 2001).
3.8.2.2.- Caso 2.- El oleaje y la corriente se propagan con un ngulo

Cuando el oleaje y las corrientes interactan con un cierto ngulo , (figura 3.25 que
muestra el vector velocidad horizontal en z = , es decir, en el lmite superior de la capa
lmite del oleaje), la velocidad combinada de olas y corrientes u wc, es:
3.71
uwc = u2 + u w2 + 2 u uw cos
[3.124]
En este caso, las tensiones de corte (ecuacin [3.111]) actan en la misma direccin
que u (ecuacin [3.124]). Por lo tanto, las tensiones de corte que actan en la misma
direccin que la corriente son:
bwc =
1
1
u + u cos
f w uwc 2 cos = f w uwc 2 w
2
2
u wc
[3.125]
mientras que la tensin de corte media en la direccin de la corriente es:
bwc =
1 T
2
1 + cos 2
dt
=
f
u
u
bwc
w bwm
T 0
[3.126]
Figura 3.25.- Velocidad instantnea en el lmite superior de la capa lmite del oleaje (Liu, 2001).
3.8.3.- Modelo de Bijker

A diferencia de lo que asume Fredse, Bijker (1971) asume que el oleaje es tan dbil
que no afectar el grosor de la subcapa viscosa de la corriente.
Considerando nicamente la corriente, Bijker (1971) asume que la subcapa viscosa
comienza desde z = 0 a z = z0 e, donde el perfil lineal de la velocidad es tangente al
perfil logartmico (figura 3.26).
3.72
Figura 3.26.- Subcapa viscosa en flujos hidrodinmicamente intensos (Liu, 2001).
Del perfil logartmico de velocidades se obtiene que:
uc = uwc
z =z 0e
u*
k
[3.127]
y que las tensiones de corte en el fondo son:

bc = u*c 2 = k 2 uc 2
[3.128]
Ahora, superponiendo el oleaje (figura 3.26), la velocidad instantnea del flujo

combinado en la parte superior de la subcapa viscosa es.
uwc = uc 2 + ubw 2 + 2uc ubw cos
[3.129]
y las tensiones de corte combinadas son:

bwc = k 2 uwc 2
[3.130]
y la tensin de corte media (que acta en la misma direccin que la velocidad

combinada uwc) es:
3.73
bwc =
1
T
1
bwc dt bc + bwm
2
[3.131]
donde T es el periodo del oleaje.

De esta manera, conociendo la velocidad media de una corriente y las caractersticas
del oleaje (H, T, h) y la altura de los ripples, ambos viajando en la misma direccin, es
posible calcular las tensiones de corte en flujos combinados de olas y corrientes,
calculando la longitud de onda con teora lineal, la amplitud del desplazamiento de la
partcula de agua en el fondo (ecuacin [3.59]) y la velocidad mxima del flujo
oscilatorio (ecuacin [3.61]). Considerando la escala de rugosidad del fondo como la
altura de los ripples, se calcula el factor de friccin del oleaje f w (ecuacin [3.81]), las
tensiones de corte mximas debidas a la velocidad mxima del flujo oscilatorio
(ecuacin [3.73]), el factor de friccin de la corriente (ecuacin [3.45]), las tensiones
debidas a la corriente (ecuacin [3.32]) y la tensin media en el fondo wc (ecuacin
[3.131]).
Cabe mencionar que el modelo de Bijker (1971) da una mayor tensin de corte debida
al flujo combinado de olas y corrientes que el modelo de Fredse (1981).
3.8.4.- Modelo de Soulsby

Por otra parte, Soulsby (1993) propone una aproximacin algebraica para las tensiones
de corte generadas bajo la accin combinada de olas y corrientes (con una precisin de
5% en la mayora de los casos):
bwc
3.2
bw
= bc 1 + 1.2

bc bw
[3.132]
donde bc y bw son respectivamente las tensiones de corte que tendran las corrientes y
las olas por separado.
La expresin correspondiente a las tensiones de corte mximas bwcm est dada por la
suma vectorial de bw (obtenida a partir de las ecuaciones de f wr [3.94]a y bw [3.73]) y
bwc (ecuacin [3.132]):
bwcm = ( bwc + bw cos ) + ( bw sin )

2
1
2 2
[3.133]
3.74
Para determinar el umbral de inicio de movimiento del sedimento y la difusin del
sedimento, es necesario el clculo de bwcm y bwc respectivamente.
Los mtodos pueden aplicarse igualmente para calcular las tensiones de corte totales en
el fondo y la contribucin de la friccin de superficie. En el primer caso, se usa la z0
total, mientras que en el segundo caso se usa z 0 = z0 s = d50 12 .
3.9.- Umbral de inicio de movimiento

En los puntos anteriores se describi la naturaleza de los flujos en la capa lmite
turbulenta del fondo desde un punto de vista puramente hidrodinmico, de donde se
dedujo que la interaccin fluido-sedimento puede representarse por las tensiones de
corte en el fondo b , que para el caso general de flujos inducidos por la accin
combinada de olas y corrientes puede expresarse por la ecuacin [3.109] (Madsen,
1993) o la ecuacin [3.133] (Soulsby, 1997). Debido a que, tanto la componente del
oleaje bwm cos(t + ) como de la corriente bc, estn asociadas con perfiles
logartmicos de velocidad inmediatamente por encima del fondo, (ecuaciones [3.98] y
[3.105] respectivamente), se deduce que la velocidad combinada inmediatamente por
encima del fondo es logartmica.
La figura 3.27 muestra, para un fondo consistente en material granular, las condiciones
fsicas del fluido en las cercanas de los granos individuales.
Figura 3.27.- Esquema de un flujo turbulento con perfil logartmico sobre un fondo granular
rugoso (CEM, 1998).
Tomando en cuenta que el perfil logartmico en esta regin es una mera extrapolacin
de las condiciones del flujo obtenidas fuera de la frontera, es evidente que la resistencia
real experimentada por el fluido no es una fuerza por unidad de rea distribuida de
forma uniforme, sino una suma de fuerzas de arrastre sobre los granos de arena
individuales, es decir:
3.75
F D,grano
b = unidad
de area
donde FD es la fuerza de arrastre sobre los granos.

No es sino hasta un cierto nivel por encima de los elementos rugosos individuales que
es razonable asumir que los vrtices, que se forman alrededor de los granos individuales
y que son arrojados hacia el flujo, se han combinado para producir un flujo de corte
turbulento, como el considerado en el modelo conceptual de las tensiones de corte
analizado anteriormente.
Adems de la variabilidad espacial de la interaccin real entre un flujo turbulento y un
fondo granular, tambin debe tenerse en cuenta que el flujo turbulento (incluso cuando
se trate de una corriente estacionaria), es bsicamente no estacionario. Por lo tanto, las
fuerzas de arrastre sobre los granos individuales no son constantes, sino que vara n con
el tiempo debido a las fluctuaciones turbulentas. No obstante, debido a que el perfil
logartmico representa el perfil de velocidades promediado en el tiempo (y por lo tanto,
las fluctuaciones turbulentas han sido eliminadas durante el proceso de promediado), las
fuerzas de arrastre sobre los granos individuales deben evaluarse como sus valores
promediados (como lo indica la barra superior en al expresin anterior).
Parece claro que para describir la compleja naturaleza de la interaccin entre el fluido y
el sedimento en la interfase entre ambos, es necesario un riguroso tratamiento
matemtico a menos que se emplee un modelo hidrodinmico que resuelva la estructura
del flujo alrededor de los granos individuales, incluyendo la variacin temporal de los
procesos fsicos asociados a las fluctuaciones turbulentas. Debido a la obvia dificultad
que entraa el desarrollo de un modelo de estas caractersticas, es necesario adoptar
modelos conceptuales heursticos (basados en la fsica que rige el problema), que sean
consistentes con el nivel de descripcin del flujo para derivar relaciones cuantitativas
que describan la interaccin fluido-sedimento y contar con evidencias empricas para
determinar las constantes necesarias para obtener las relaciones predictivas.
Considerando un flujo estacionario sobre un fondo con sedimento no cohesivo, las
fuerzas que actan sobre el grano se ilustran en la figura 3.28.
En esta situacin, el mecanismo impulsor del movimiento del grano es la fuerza de
friccin (arrastre) que ejerce el flujo sobre ste (similar a las ecuaciones [3.29] y [3.31]):
1
d 2
2
FD = CD
( u* )
2
4
[3.134]
donde u* es la velocidad de corte (o velocidad de friccin) que es la velocidad del flujo

en las inmediaciones del fondo, es un coeficiente usado para modificar el valor de u*
de tal manera que el trmino (u* ) indica la velocidad del flujo caracterstica que pasa
por el grano. La fuerza estabilizadora puede considerarse como la fuerza de friccin que
acta sobre el grano.
3.76
Figura 3.28.- Fuerzas que actan sobre el sedimento en el fondo debidas a una corriente estacionaria
(Liu, 2001).
El umbral del movimiento del sedimento en el fondo marino es uno de los factores ms
importantes en los clculos relacionados con la respuesta del sedimento a la accin de
las olas y las corrientes. Es de especial importancia en lo s problemas que involucran
erosin alrededor de las estructuras, clculos de la movilidad del sedimento en rellenos
o extracciones de arena y transporte del sedimento en general.
Para flujos con poca velocidad sobre fondos arenosos, la arena permanece bsicamente
inmvil. Al incrementarse la velocidad del flujo de forma paulatina es evidente que
llega un punto en el cual el fluido es capaz de movilizar al sedimento. Este punto es el
llamado umbral del movimiento del sedimento. Este umbral se puede alcanzar
igualmente por el flujo oscilatorio asociado al oleaje y al flujo combinado de olas y
corrientes.
Para una corriente estacionaria, la velocidad umbral verticalmente integrada U cr
necesaria para mover un grano de dimetro d en un fondo plano, horizontal, sin ripples,
a una profundidad h, puede ser predicha por distintos mtodos. Por ejemplo, Van Rijn
(1984) propone una ecuacin vlida para agua dulce a 15 de temperatura y una
densidad del sedimento de s = 2650 kg / m3 :
4h
U cr = 0.19(d 50 )0.1 log10
para 100 d 50 500m
d90
[3.135]a, b
4h
U cr = 8.5(d 50) 0.6 log10
para 500 d 50 2000m
d
90
3.77
Soulsby (1997) propone una expresin para calcular esta velocidad crtica en funcin
de una tensin de corte umbral en el fondo y de la ley de friccin (ecuacin [3.47]),
vlida para sedimentos no cohesivos y para condiciones para en las cuales el tamao
adimensional del grano D* > 0.1:
h
U cr = 7
d50
1
7
( g ( s 1) d50 f ( D* ) ) 2 ( para D*
> 0.1)
donde:
f ( D* ) =
[3.136]a, b
1
3
0.30
g ( s 1)
+ 0.055 (1 exp(0.020 D* ) ) y D* =
2
d 50
1 + 1.2D*

donde s = relacin entre las densidades del agua y del sedimento y es la viscosidad
cinemtica del agua.
Bajo la accin del oleaje, el umbral de movimiento del sedimento depende de la
magnitud de la velocidad orbital en el fondo ubwm , el periodo del oleaje T, el dimetro
del tamao de grano d y su densidad s. La velocidad orbital umbral ubwcr puede
obtenerse a partir de las ecuaciones propuestas por Komar & Miller (1974):
2
3
1
3
1
3
ubwcr = ( 0.11g ( s 1) ) d T para d < 0.5 mm
ubwcr = (1.09 g ( s 1) ) 7 d 7T 7 para d > 0.5 mm
[3.137]a
[3.137]b
donde g es la gravedad y s la relacin entre las densidades del grano y del agua.
Estas ecuaciones son ampliamente utilizadas, pero tienen la desventaja de tener una
gran discontinuidad en d = 0.5 mm. Como alternativa, la velocidad orbital umbral
tambin puede deducirse a partir de la tensin umbral de corte en el fondo.
Soulsby (1997) propone un mtodo continuo en todo el rango de tamaos de sedimento
y bajo condiciones de flujo laminar, suave y altamente turbulento que es totalmente
compatible con el mtodo del umbral de movimiento inducido por corrientes (ecuacin
[3.136]). Este mtodo consiste en una familia de curvas de ubwcr como funcin del
tamao del sedimento para un rango de periodos de oleaje en caso especfico de granos
de cuarzo en agua de mar a 10 C y 35 ppm. Estas curvas fueron calculadas para cada
tamao de grano, incrementando paulatinamente la velocidad orbital ubw y obteniendo
de ah las tensiones de corte. Cuando las tensiones de corte calculadas son iguales al
valor umbral, entonces ubwcr = ubw. Estos clculos, sin embargo, no permiten la
obtencin de una expresin analtica debido a la complejidad de incluir las expresiones
del factor de friccin bajo flujos laminares, suave o altamente turbulentos. La figura
3.78
3.29 muestra una comparacin entre las predicciones con una serie de datos
experimentales (Van Rijn, 1989). Cabe notar que las curvas tericas sub-predicen los
datos experimentales, debido posiblemente a que es ms representativo comparar los
umbrales bajo corrientes con las tensiones de corte medias del oleaje en lugar de
compararlas con las tensiones de corte de pico. Como es de esperar, ubwcr se incrementa
con el periodo, tal como lo indican las curvas tericas (tal como sucede con las
ecuaciones de Komar y Miller, 1974).
Figura 3.29.- Velocidad oscilatoria umbral asociada al oleaje para el movimiento de sedimentos
(Van Rijn, 1989).
3.9.1.- Inicio del movimiento. Parmetro de Shields

A partir de la descripcin anterior de la naturaleza de los flujos turbulentos en las
cercanas del fondo, se puede deducir que las tensiones de corte en el fondo se pueden
considerar como una magnitud fsica relacionada con las fuerzas de arrastre, es decir, la
fuerza que al actuar sobre el sedimento en la superficie del fondo lo pone en
movimiento. Aceptando que para un sedimento con un dimetro d hay
aproximadamente 1/d 2 granos por unidad de rea, entonces la fuerza que acta sobre los
granos del fondo es:
FD,grano b d 2
[3.138]
la doble barra indica que est espacial y temporalmente promediado.

Para un sedimento no cohesivo, los granos individuales de sedimento yacen y
permanecen en el fondo debido a su propio peso, y resisten al movimiento horizontal
debido a la interaccin (friccin) con los granos adyacentes. Por lo tanto, para un
3.79
sedimento no cohesivo, la fuerza estabilizadora est asociada con el peso sumergido de
los granos individuales:
Wgrano ( s ) g d 3
[3.139]
donde s es la densidad del sedimento ( 2650 kg / m3 para cuarzo).

La relacin entre la fuerzas movilizadoras (fuerzas de arrastre) y la fuerza
estabilizadora (peso sumergido) es un concepto fsico fundamental en la interaccin
fluido-sedimento para sedimento no cohesivo. Esta relacin es lo que se conoce como el
parmetro de Shields (Shields, 1936).
Considerando que, a partir de una cierta tensin de corte (fuerza) crtica (que se puede
expresar en trminos de la velocidad de corte crtica u*,c), el sedimento puede comenzar
a moverse, el parmetro de Shields es una expresin que denota la situacin donde el
sedimento esta a punto de ponerse en movimiento, justo donde la fuerza de arrastre es
igual a la velocidad de friccin, es decir:
1
d2
d3 1
d2
2
C D
C L
( u* c ) = f ( s ) g
( u*c )
2
4
6
2
4
[3.140]
sta ecuacin puede ser re-ordenada como:
u*c 2
f
4
= 2
2
( s 1) g d CD + f CL 3 2
[3.141]
donde el parmetro de Shields se define como:

u*2
b
b
=
=
=
( s 1) g d ( s 1) g d (s ) g d
[3.142]
donde s = (s /) es la densidad relativa del sedimento con respecto al fluido.

De tal manera que, el sedimento comienza a moverse cuando:
u* > u*cr, es decir, que la velocidad de corte del flujo es mayor que la velocidad
de corte crtica.
b > bcr, es decir, que la tensin de corte crtica es mayor que la tensin de corte
umbral (b = u*,c2 ).
3.80
> cr, es decir, que el parmetro de Shields es mayo r que el parmetro de

Shields crtico (ecuacin [3.142]).
Para un flujo turbulento sobre un fondo no rugoso, consistente en sedimento no

cohesivo, con dimetro d, la rugosidad de grano equivalente de Nikuradse es k n = d. De
la descripcin de la naturaleza de flujos turbulentos en la capa lmite, se vio que la
escala o magnitud de la rugosidad experimentada por el flujo, z0 (ecuacin [3.19])
depende del rgimen (moderadamente o altamente turbulento) del flujo. El parmetro
que determina las caractersticas (rgimen) del flujo en las cercanas del fondo y por lo
tanto de la fuerza movilizadora que acta sobre los granos individuales del sedimento es
el nmero de Reynolds en el fondo:
Re s =
u*k n u*d
=
[3.143]
De aqu se deduce que las condiciones de estabilidad de un grano de sedimento en una

interfase fluido sedimento puede expresarse como un valor crtico del parmetro de
Shields, el cual a su vez, es funcin de un nmero de Reynolds crtico (Re* ) que ha de
ser determinado a partir de modelos numricos o datos experimentales:
cr =
bcr
u*cr 2
=
= f (Re s )
g ( s 1) d g ( s 1) d
[3.144]
La figura 3.30 muestra el diagrama tradicional del parmetro de Shields (ecuacin

[3.144]) contrastado con datos experimentales obtenidos de un flujo uniforme y
estacionario (Raudkivi, 1990) y que relaciona al parmetro de Shields con el nmero de
Reynolds. Esto muestra que cr 0.06 es aproximadamente constante para Re s > 100, es
decir, para flujos totalmente turbulentos, mientras que el parmetro crtico de Shields,
cr aumenta desde un valor mnimo de 0.035 para valores de Re s < 10, correspondiente a
un rgimen de flujo suavemente turbulento. En la figura 3.30 se distinguen 3 zonas bien
diferenciadas, correspondientes a 3 situaciones de flujo distintas:
Flujo hidrodinmicamente suave, para Res 2, donde dn es mucho mas pequeo

que el grosor de la subcapa viscosa, por lo que el sedimento se encuentra dentro
de la subcapa viscosa y por lo tanto, cr es independiente del dimetro del grano.
Experimentalmente se ha encontrado que cr = 0.1/ Res.
Flujos en transicin, para 2 Res 500, donde el tamao del grano es del
mismo orden que el grosor de la subcapa viscosa. Hay un valor mnimo para c
= 0.032 correspondiente a un Re s = 10.
Flujos hidrodinmicamente fuertes, para Re s 500, donde no existe la subcapa

viscosa y por lo tanto cr es independiente de la viscosidad del fluido. cr tiene
un valor constante de 0.06
3.81
Figura 3.30.- Diagrama de Shields para inicio de movimiento en un flujo turbulento estacionario
(Raudkivi, 1990) de donde se obtiene el valor crtico c en funcin del nmero de Reynolds del grano
(sedimento uniforme y no cohesivo).
Para establecer las condiciones crticas (umbrales) de inicio de movimiento, (ecuacin

[3.144]), es importante tener en cuenta la naturaleza de la interaccin entre el fluido y el
sedimento. Para valores de > cr las fuerzas de movimiento son mayores que las
fuerzas estabilizadoras, por lo que el sedimento estar en movimiento. Sin embargo,
esto no significa que para valores de ligeramente menores que cr, el sedimento no se
mueva en absoluto. Para un flujo turbulento, en las proximidades de cr, las fuerzas
de arrastre sobre algunos granos de sedimento pueden incrementarse momentneamente
por vrtices turbulentos (especialmente fuertes) provenientes de la capa de generacin
de turbulencia, por lo que si el tiempo de respuesta de los granos individuales es menor
que la duracin de la fluctuacin (pulso turbulento), puede provocarse una dislocacin
de los granos desde su posicin de reposo e iniciarse el movimiento. La curva emprica
de un diagrama de vs. Res que exprese la condicin crtica de movimiento (la curva de
Shields o el criterio de Shields), debe de ser interpretada como un rea difusa alrededor
de la propia curva. Para un flujo cuyas condiciones se encue ntren dentro de esta rea, el
movimiento del sedimento puede tener lugar, pero ser cada vez ms espordico
conforme caiga debajo de la curva que define cr.
3.9.2.- Diagrama modificado de Shields

Para poder predecir las condiciones del flujo que provocarn el movimiento del
sedimento con un cierto tamao (dimetro) conocido, la utilizacin del diagrama
(convencional) de Shields puede ser un tanto engorrosa, debido a que la caracterstica
ms importante del flujo involucrada en el movimiento del sedimento, la velocidad de
corte crtica u*cr est incluida en ambos ejes. Este problema puede solucionarse si se
3.82
toma en cuenta que la curva de Shields define una nica relacin entre cr y Re* . As, a
partir de la definicin de cr (ecuacin [3.144]) se obtiene que :
u*cr = ( s 1) gd cr
[3.145]
Aplicando esto en la definicin del nmero de Reynolds en el fondo Re* , se obtiene el

parmetro:
S* =
d ( s 1) gd
Re s
=
4
4 cr
[3.146]
Tomando los valores correspondientes de Re s y cr del diagrama tradicional de

Shields, se puede obtener el valor del parmetro sedimento - fluido S* que luego puede
usarse para reemplazar a Re s en el diagrama de Shields. De esta manera se obtiene el
diagrama modificado de Shields (Madsen y Grant, 1976) tal como se muestra en la
figura 3.31. Cabe mencionar que la razn de que aparezca el factor 4 en la definicin de
S* es para obtener valores numricos de S* comparables con los valores de Re* en el
diagrama de Shields tradicional. Esto est hecho meramente por conveniencia sin tener
ningn significado fsico (Madsen, 1993).
Figura 3.31.- Diagrama modificado de Shields que muestra c en funcin de S * (Madsen & Grant, 1976).
El lmite inferior del diagrama modificado de Shields, en S* = 1, es una regin que,

para sedimento compuesto por cuarzo en agua de mar, corresponde a un dimetro de
sedimento 0.1 mm, es decir, una arena muy fina. Aunque la cohesin puede ser
3.83
importante para dimetros menores que 0.1 mm, si estos se encuentran limpios de
materia orgnica, pueden ser considerados como no cohesivos y por lo tanto, estar
gobernados por el criterio de inicio de movimiento de Shields.
Ajustando los datos experimentales sobre condiciones del umbral de inicio de
movimiento presentados por Raudkivi (1990), obtenidos para valores bajos del nmero
de Reynolds del sedimento Re s (en el rango de 0.03 a 1), se obtiene que:
cr = 0.1Re s
1
3
para Re s < 1
[3.147]
que interpretado en el diagrama modificado de Shields (propuesto por Madsen y Grant,

1976), puede rescribirse como:
cr = 0.1 S*
2
7
para S* < 0.1
[3.148]
Para valores grandes de S* , el rgimen del flujo correspondiente a las condiciones de

inicio del movimiento son totalmente turbulentas, es decir, Re s > 100, y se puede usar
un valor de cr 0.06 para valores de S* mayores que el rango cubierto en la figura 3.32.
Figura 3.32.- Comparacin de la curva de Shields con datos experimentales de las condiciones de inicio
de movimiento en flujos turbulentos oscilatorios (Madsen & Grant, 1976).
El criterio modificado de Shields para el inicio de movimiento para sedimento no

cohesivo mostrado en la figura 3.32 se obtuvo a partir de experimentos con flujos
3.84
estacionarios turbulentos. No obstante, no debe sorprender el hecho de que este criterio
pueda ser tambin aplicable a flujos turbulentos en capa lmite no estacionarios, como
se demuestra en la misma figura 3.32 (Madsen y Grant, 1976) si se considera la
naturaleza de la puesta en movimiento del sedimento. Para las condiciones de flujo en
las proximidades de cr, el movimiento espordico de algunos granos se debe a las
fluctuaciones turbulentas de alta frecuencia de la fuerza movilizadora que acta sobre
los granos. Debido a que el perfil medio de velocidades cerca del fondo es logartmico,
tanto para el flujo asociado al oleaje como para corrientes, es fsicamente razonable
esperar que la turbulencia cercana al fondo sea similar para ambos casos, es decir, que
sea funcin de los valores instantneos de la velocidad de corte en el fondo. Por lo tanto,
dado que el tiempo de respuesta de los granos individuales del sedimento es
relativamente corto con respecto a la escala de tiempo de las fluctuaciones turbulentas
(que se espera que sea la misma para olas y corrientes, si sus velocidades de corte son
iguales), el inicio de movimiento tan slo se ver afectado por las inestabilidades si el
periodo del oleaje es del mismo orden que la escala temporal de las fluctuaciones
turbulentas. Debido a que esto en general no sucede, los efectos de las inestabilidades
del flujo medio son despreciables y los granos del sedimento reaccionan de forma casi
instantnea a las tensiones de corte, es decir, el inicio de movimiento se obtiene de:
m =
u*m 2
= cr
( s 1) gd
[3.149]
donde m es el parmetro de Shields que denota movimiento, u*m denota la mxima

velocidad de corte, es decir, u*m = u*c, u*wm y u*cwm para corrientes puras, oleaje puro y
flujos combinados de olas y corrientes respectivamente.
El efecto de la pendiente del fondo en el inicio de movimiento de sedimento en la
interfase fluido-sedime nto puede ser tomada en cuenta considerando una modificacin
(Madsen, 1991) en el balance entre la fuerzas de arrastre, gravedad y friccin
(resistencia al movimiento) a lo largo de un fondo inclinado con un cierto ngulo con
respecto a la horizontal en la direccin del flujo, donde es positiva si la pendiente del
fondo es cuesta arriba en la direccin del flujo. Este balance de fuerzas es (Madsen,
1993):
1
2 2
3
3
CD d ucr, ( s ) g d sin = ( s ) g d cos tan s [3.150]
2
4
6
6
donde CD y s son el coeficiente de arrastre y el ngulo de friccin respectivamente,

para un grano de arena esfrico en la interfase fluido-sedimento, y ucr, es una velocidad
representativa del fluido que se usa en la evaluacin de la fuerza de arrastre.
Reordenando la ecuacin [3.150] se obtiene:
3.85
2
ucr ,
( s 1) g d
4
tan
tan s cos 1 +
3CD
tan s
[3.151]
donde los trminos en el parntesis representan los efectos de la pendiente del fondo,
con positivo si la pendiente del fondo es cuesta arriba en la direccin del flujo,
mientras que los trminos restantes expresan el balance de fuerzas para un fondo
horizontal. Aceptando que este balance de fuerzas sobre un fondo horizontal expresa el
criterio de Shields, este resultado sugiere que el parmetro crtico de Shields para un
flujo sobre un fondo con pendiente puede expresarse como:
tan
cr, = cr cos 1 +

tan s
[3.152]
De hecho, este esquema simple (ecuacin [3.151]) puede ser usado para obtener una
expresin terica para el parmetro crtico de Shields, tomando ucr como la velocidad
dada por el perfil logartmico de velocidades evaluado en z = zr. Entonces, la ecuacin
[3.151] para un fondo horizontal con = 0 cambia a:
k
ucr 2
4
zr
= cr =
( s 1) gd
3CD ln
z0
tan
s
[3.153]
Aunque de manera menos intuitiva que los argumentos fsicos presentados para la
deduccin de la ecuacin [3.144] (parmetro crtico de Shields), esta expresin muestra
que cr es funcin de Re s, ya que z0 (ecuacin [3.19]) tambin es funcin de Re * . Para
una flujo totalmente turbulento (con z0 = k n / 30 = d / 30), la ecuacin [3.153] lleva a un
valor constante de cr para Re s > 100, que coincide cualitativamente con diversas
observaciones realizadas. Tomando valores razonables de zr = 0.8 y CD = 0.4, este
modelo simple sugiere un valor de cr = 0.051 tan s, el cual tambin concuerda
cuantitativamente con las observaciones para s 50 (Madsen, 1991), ya que cr
(0.051 tan 50) 0.06 para un flujo totalmente turbulento.
El diagrama modificado de Shields (figura 3.31) es mas fcil e intuitivo de usar que el
diagrama tradicional ya que, conociendo las caractersticas del fluido y del sedimento, s,
d y t, se puede obtener S* a partir de la ecuacin [3.146]. A partir de aqu, conociendo
S* , puede obtenerse cr a partir de la figura 3.31 o de la ecuacin [3.148] para S* < 0.8,
mientras que la velocidad de corte crtica puede obtenerse de la ecuacin [3.145]. De
igual manera, conociendo la densidad del sedimento, y la viscosidad del fluido, se puede
estimar la tensin de corte en el fondo crtica para el inicio de movimiento bcr para un
dimetro dado, encontrando el parmetro del fluido sedimento (S* ) (ecuacin [3.146]),
encontrando u*cr con la figura 3.31 y finalmente la bcr (ecuacin [3.13]).
3.86
3.10.- Formas geomtricas y rugosidad del fondo

Para situaciones de flujo estacionario, como una corriente, el sedimento comienza a
moverse cuando se supera un valor umbral del parmetro de Shields (parmetro crtico
de Shields), definido anteriormente en la ecuacin [3.142] y en la figura 3.31.
Para estimar el inicio del movimiento bajo la accin del oleaje, es posible utilizar
directamente el diagrama de Shields, sustituyendo las tensiones de corte inducidas por
la corriente por las tensiones de corte mximas (definidas por Jonsson, 1966) inducidas
por el flujo oscilatorio del oleaje (ecuacin [3.73]), as como el Parmetro del fluidosedimento S* (ecuacin [3.146]). Para calcular las tensiones de corte mximas inducidas
por el oleaje, se calcula la velocidad mxima de la partcula de agua cerca del fondo y la
amplitud del desplazamiento por medio de las ecuaciones [3.61] y [3.59]
respectivamente, y el coeficiente de friccin del oleaje por medio de la ecuacin [3.81].
Una vez el sedimento ha comenzado a moverse, comienzan a generarse formas
geomtricas del fondo. En estudios de laboratorio se ha visto que, en funcin del
rgimen hidrodinmico, la geometra evoluciona de ser plano a tener ripples, a dunas, a
condiciones de sheet flow y a antidunas (figura 3.33):
i.
Ripples.- Se forman en regimenes hidrodinmicos relativamente suaves, de baja

intensidad y estn asociados a sedimentos relativamente finos, con d50 < 0.7 mm.
El tamao de los ripples est controlado principalmente por el tamao del grano.
Las alturas y longitudes tpicas de los ripples estn dadas por las siguientes
expresiones empricas:
H ripples 100 d50

Lripples 1000 d50
[3.154]
En flujos de baja intensidad, los ripples tienen formas muy regulares (2D), con una
pendiente de 6 en la cara expuesta del ripple y una pendiente de 32 en la cara
protegida del ripple. Con un incremento en la intensidad del flujo, los ripples se vuelven
tridimensionales.
ii.
Dunas.- La forma de las dunas es muy similar a la de los ripples, pero son
mucho ms grandes. El tamao de las dunas est controlado por la profundidad
del fondo y se asocia a tamaos de grano gruesos, con un d50 mayor que 0.6 mm.
Con el incremento en la intensidad del flujo, las dunas crecen y la profundidad
del agua donde estn las dunas disminuye. Esto da lugar a un incremento en la
velocidad del flujo en las crestas de las dunas, lo que puede provocar que, en
caso de un aumento en el rgimen hidrodinmico, las dunas sean lavadas, el
fondo se vuelva plano y se generen condiciones de flujo en capas (sheet flow).
3.87
iii.
Antidunas.- Se forman cuando el nmero de Froude es mayor a 1. En estas

condiciones, la altura del oleaje formado en la superficie es del mismo orden de
magnitud que la altura de la anti-duna. La superficie ondulada del fondo se
vuelve inestable y puede crecer y romper en contra de la direccin del flujo, lo
que provoca que las antidunas se muevan en direccin contraria al flujo.
Figura 3.33.- Flujo sobre ripples, dunas y antidunas (Liu, 2001).
De tal manera que conociendo la velocidad media de la corriente, la profundidad del

agua y el tamao del sedimento, es posible predecir las posibles formas generadas en el
fondo por medio de diagramas empricos (e.g. Znamenskaya, 1969), donde el tamao
del sedimento se representa por medio de la velocidad de cada del sedimento (s). La
velocidad de movimiento de los ripples tambin se muestra en dicho diagrama (figura
3.34).
Figura 3.34.- Formas del fondo en funcin de la velocidad de cada, velocidad de la corriente y el nmero
de Fraude (Raudkivi, 1990).
3.88
3.10.1.- Generacin de formas del fondo por accin del oleaje

Cuando la intensidad del flujo (expresada en trminos del parmetro de Shields, )
excede una condicin crtica de inicio de movimiento, cr, los granos de sedimento
comenzarn a moverse de manera instantnea. Sin embargo, para valores de que
excedan ligeramente (del 10 al 20%) el valor de cr, comenzar la formacin de formas
geomtricas en el lecho.
La aparicin de formas en la interfase fluido sedimento en el fondo marino provoca
un cambio obvio en su geometra. Mas que provocar un incremento en la resistencia del
flujo asociada a las fuerzas de arrastre que se generan sobre los granos individuales (es
decir, la rugosidad dada por el dimetro del sedimento), el flujo se separar en la cresta
de las formas del lecho y entonces, la resistencia del flujo se deber principalmente a las
fuerzas de presin del flujo sobre las formas geomtricas del lecho. En otras palabras, la
rugosidad pasar a depender de la geometra de los ripples. La aparicin de los ripples
en el fondo cambia la rugosidad del fondo en algunos rdenes de magnitud.
A pesar del incremento en la rugosidad del fondo debida a la presencia de los ripples,
son las fuerzas de arrastre que actan sobre los granos individuales (no las fuerzas de
arrastre sobre las formas del lecho) las responsables del movimiento del sedimento
inducido por el flujo.
La friccin con el fondo, en trminos de las tensiones de corte en el fondo, es la tensin
de corte promedio que acta sobre el sedimento. La idea de dividir las tensiones de corte
en el fondo en las debidas a la friccin y las debidas a la presin sobre las formas del
fondo, es decir:
b = bs + bf
[3.155]
(donde bs denota la componente debida a la friccin y bf la componente de presin del

flujo sobre las formas del fondo, como se puede observar en la figura 3.35) ha recibido
una atencin considerable en el contexto de la mecnica del transporte de sedimentos en
flujos turbulentos estacionarios (Raudkivi, 1990).
Para un oleaje puro, las tensiones de corte debidas a la friccin, se pueden obtener de
forma simple (Madsen y Grant, 1976) evaluando las tensiones de corte en el fondo a
partir de las condiciones de oleaje dadas por ubm y , y suponiendo que la escala de
rugosidad del fondo sea el dimetro del sedimento d. As, las tensiones de corte debidas
a la friccin se obtienen como:
1
2
bwsm = f ws ubwm
2
[3.156]
3.89
donde f ws es el factor de friccin del oleaje obtenido de las ecuaciones [3.87], [3.88] o
[3.90] para una rugosidad de:
kn = kn ' = d
[3.157]
Figura 3.35.- Esquematizacin de las componentes de friccin con el fondo y presin de arrastre en un
flujo turbulento sobre un fondo arenoso con formas geomtricas (CEM, 1998).
Para flujos combinados de olas y corrientes, Glenn & Grant (1987) proponen la
evaluacin de las tensiones de corte debidas a la friccin de una forma similar, es decir,
usando el modelo de interaccin olas y corrientes descrito anteriormente, con la
rugosidad del fondo definida por la ecuacin [3.157]. Sin embargo, si un flujo de capa
lmite combinado por olas y corrientes se define por la velocidad de la corriente fuera de
la capa lmite del oleaje, es decir, por uc (zr), con zr > cw, la utilizacin directa de esta
informacin para describir la capa lmite para flujos combinados de olas y corrientes,
puede dar lugar a una determinacin un tanto ambigua de la friccin. Ello es debido a
que esta aproximacin supone que la variacin de la velocidad de la corriente fuera de
la capa lmite del oleaje estar en funcin de la rugosidad determinada por el tamao de
grano, es decir, por k n = kn = d, lo cual es slo cierto para un fondo plano (sin formas
geomtricas), es decir, en ausencia de fuerzas de presin sobre las formas geomtricas
del fondo.
Para solucionar este problema, primero se aplica la teora desarrollada para flujos
combinados de olas y corrientes utilizando k n como la rugosidad total de los granos de
arena equivalente de Nikuradse, para tomar en cuenta la presencia de formas del lecho
al predecir la velocidad de la corriente en el borde exterior de la capa lmite:
uc ( zr ) con z r = cw
[3.158]
se obtiene a partir de un modelo que considere las tensiones totales de corte en el fondo.
3.90
Con la corriente definida por la ecuacin [3.158], las tensiones de corte debidas a la
friccin del flujo con el fondo se calculan a partir de la teora descrita para flujos
combinados, usando la rugosidad del grano k n = k n = d. De esta manera se evita la
ambigedad en la prediccin de las tensiones de corte debidas a la friccin
(Wikramanayake & Madsen, 1994).
3.10.2.- Generacin de ripples

Wikramanayake & Madsen (1994) realizaron un estudio comparando datos de la
geometra de los ripples en condiciones de campo contra el parmetro de Shields basado
en la tensin de corte mxima debida a la friccin m , y encontraron que, aunque hay
una dispersin importante en los resultados, el parmetro que proporcionaba una mejor
correlacin entre la geometra de los ripples y los datos de campo era:
Z=
'm
S*
[3.159]
donde S* es el parmetro fluido sedimento (ecuacin [3.146]). Wikramanayake &

Madsen (1994) realizaron una comparacin de este parmetro contra la altura
normalizada Abmr (donde Abmr = ubwr r ), y el peralte de los ripples, adquiridos
en experimentos de campo, obteniendo las siguientes relaciones empricas
1.8x102 Z 0.5 para 0.0016 < Z < 0.012
=
Abmr
4 1.23
para 0.012 < Z < 0.18
7.0x10 Z
[3.160]
1.5x101 Z 0.009 para 0.0016 < Z < 0.012

=

2 0.65
para 0.012 < Z < 0.18
1.05x10 Z
[3.161]
donde Abm r es la amplitud de la trayectoria de la partcula de agua predicha por la teora

lineal (ecuacin [3.86] evaluada con los parmetros representativos del oleaje irregular),
y son la altura y longitud de onda del ripple respectivamente, y los rangos superior e
inferior de validez indican el rango cubierto por los datos experimentales.
3.91
Los datos de campo muestran una dispersin considerablemente mas grande con
respecto a los datos de laboratorio. Esto en realidad es una situacin esperable si se
toman en consideracin las dificultades asociadas a la obtencin de datos de campo con
respecto a los datos de laboratorio, aunque los datos de campo (altura y peralte de los
ripples) muestran la misma tendencia general que los datos de laboratorio. Para el rango
inferior de las intensidades del flujo, 0.0012 < Z < 0.016, la variacin de la geometra
normalizada de los ripples con respecto de Z es relativamente baja, mientras que la
geometra vara rpidamente para valores altos de Z. Finalmente, cuando Z excede el
valor de 0.018, el fondo se vuelve plano de nuevo, dando lugar a lo que se conoce como
sheet flow.
3.10.3.- Geometra de los ripples

A partir de la descripcin del oleaje irregular mediante sus parmetros representativos
(periodo representativo y la altura de ola media cuadrtica) y las caractersticas del
sedimento, se puede obtener las tensiones de corte bwm,s debidas a la friccin del flujo
oscilatorio (ecuacin [3.156]).
Una vez conocido el parmetro bwm,s, se puede obtener el parmetro de Shields basado
en la mxima tensin de corte debida a la friccin m,s:
m, s =
b w m, s
( s 1)gd
u*wm 2
( s 1) gd
[3.162]
y si m,s > 2cr (obtenido a partir de la curva modificada de Shields, figura 3.31), se
puede suponer que se inicia el movimiento del sedimento y se puede obtener el
parmetro Z (ecuacin [3.159]).
As, una vez conocido Z, se pueden usar las relaciones empricas dadas en las
ecuaciones [3.160] y [3.161] para obtener la geometra de los ripples. Si m,s > 0.35 o si
Z > 0.18, se supone que el fondo se vuelve plano de nuevo y que hay condiciones de
sheet flow.
3.10.4.- Rugosidad del fondo en presencia de ripples

La rugosidad del fondo k n , tambin llamada rugosidad de grano equivalente de
Nikuradse es una forma de medir la irregularidad del fondo. La nica situacin en la que
se puede obtener directamente la rugosidad del fondo es en un lecho plano consistente
en esferas uniformes, donde k n = dimetro de la esfera. Como se mencion
anteriormente, la rugosidad del fondo se puede obtener indirectamente a partir de
medidas del perfil de velocidades (realizando un ajuste logartmico). A partir de una
gran cantidad de datos sobre los valores de la escala de rugosidad del fondo, obtenidos
por medio de medidas del perfil de velocidades, en distintas situaciones de rgimen
3.92
hidrodinmico y tamaos de sedimento, se han obtenido los valores aproximados de la
escala de rugosidad del fondo:
(1 10 ) d 50
kn
100d 50 H ripple
fondo plano
fondo con ripples
[3.163]
En un flujo estacionario, la secuencia de formas del fondo que se forman una vez se
han superado las condiciones umbrales para el inicio del movimiento:
i.
ii.
iii.
iv.
Fondo plano.
Rolling grain ripples (que se forman en las proximidades del umbral de inicio de
movimiento, cuando los granos comienzan a moverse, pero no se levantan del
fondo).
Vortex ripples (se dan en condiciones de > 2cr ).
flujo en capas (en condiciones de > 0.413 d[cm]. El fondo se vuelve plano).
Los ripples inducidos por el flujo oscilatorio del oleaje se pueden distinguir de los
ripples formados por las corrientes gracias a sus caractersticas geomtricas. Los ripples
inducidos por el oleaje tienden a ser ms simtricos debido a la oscilacin de la
partcula del agua (figura 3.36).
Figura 3.36.- Geometra de los ripples inducidos por el oleaje y por las corrientes respectivamente.
Desde el punto de vista hidrodinmico, la rugosidad equivalente de Nikuradse de un

fondo con ripples es ms importante que la propia geometra de los ripples. Sin
embargo, es obvio que la geometra de los ripples debe jugar un papel importante en la
determinacin de la resistencia al flujo que experimenta el fluido en su interaccin con
el fondo. Para traducir la geometra de los ripples en trminos de la rugosidad
equivalente de Nikuradse (es decir, la determinacin de la rugosidad de un fondo
mvil), es necesario tener informacin y datos experimentales de la resistencia del flujo
sobre un fondo natural mvil con presencia de ripples. Desafortunadamente la falta de
datos experimentales bajo condiciones de oleaje natural limitan seriamente la
posibilidad de la determinacin de la rugosidad de un fondo mvil.
3.93
Para el caso de los ripples generados por el oleaje, su longitud (Lr) es aproximadamente
2 veces la amplitud del desplazamiento de la partcula del agua Abm (ecuacin [3.59]):
Lr = 2 Abm
[3.164]
y el peralte est dado por (Nielsen, 1979):
3
Hr
= 0.182 0.24 ( m, s ) 2
Lr
[3.165]
donde Hr es la altura del ripple y m,s es el parmetro de Shields efectivo debido a la

friccin, es decir, el parmetro de Shie lds calculado tomando en cuenta la rugosidad del
fondo como 2.5 d50 . La escala de rugosidad del fondo (con ripples) estar escalada por
la altura de los ripples:
k nr = (1 4) H r
[3.166]
Para determinar la rugosidad del fondo en condiciones de campo, Madsen, (1993)

sugiere que para condiciones de oleaje irregular definido en trminos del periodo
representativo, la altura de los ripples se determine usando la ecuacin [3.160] y la
rugosidad del fondo mvil k nr a partir de la ecuacin [3.166]. Si el valor de m,s no
indica inicio de movimiento, es decir, 2 m,s < cr, a menos que haya informacin del
fondo disponible, la rugosidad del fondo ser k n = kn = d.
Si m,s > 0.35 o si Z > 0.18, es de esperar que el fondo sea plano, es decir, que haya
condiciones de sheet flow. Algunos resultados de Madsen (1993) sugieren un valor
aproximado de k n = (15 d) en condiciones de sheet flow. Debido a que las condiciones a
partir de las cuales fueron obtenidos estos resultados corresponden a valores de m,s 1
y a que es fsicamente razonable suponer que la rugosidad de un fondo mvil para
condiciones de sheet flow depende de la intensidad del transporte de sedimento sobre
un fondo plano (que de hecho est relacionado con la magnitud de m,s), es vlido
suponer un valor tentativo de:
k n = 15 m, s d
[3.167]
Esta expresin tiene la misma forma de otras propuestas (e.g. Wilson, 1989), aunque
los coeficientes difieren.
As, a partir de conocer la densidad y dimetro del sedimento, la densidad y viscosidad
del agua, y las caractersticas del oleaje (H, T, h) es posible calcular la rugosidad del
3.94
fondo, calculando la longitud de onda (teora lineal), la amplitud del desplazamiento de
la partcula Abm (ecuacin [3.59]), la velocidad mxima de la partcula de agua ubwm
(ecuacin [3.61]), el factor de friccin del oleaje correspondiente a k n = 2.5 d50
(ecuacin [3.81]), las tensiones de corte efectivas en el fondo bwm (ecuacin [3.73]), el
parmetro de Shields efectivo m,s (ecuacin [3.142] para bwm ), la longitud del ripple
(ecuacin [3.164]), el peralte del ripple (ecuacin [3.165]) con lo que se encuentra la
altura r y finalmente la rugosidad con ripples k nr (ecuacin [3.166]).

07 Ocma 07 de 13

Cargado por

Copyright:

Formatos disponibles

07 Ocma 07 de 13

Cargado por

Información del documento

Descripción original:

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

07 Ocma 07 de 13

Cargado por

Copyright:

Formatos disponibles

3.

3.- Transporte de Sedimentos

3.2.- Factores que regulan la variacin de la morfologa de playas

3.3.- Escalas de los procesos costeros

Tabla 3.1.- Clasificacin de las escalas de los procesos costeros.

Tericamente hablando, los fenmenos a mesoescala deben componerse de la

Forma, Erosin, Acrecin

Figura 3.1.- Factores que intervienen en la morfodinmica.

Las relaciones entre la hidrodinmica, el transporte de sedimentos y la morfologa son

Figura 3.2.- Escala temporal de variaciones de la lnea de costa y procesos costeros

Desde la profundidad activa hasta la zona de rompientes.

Caractersticas del flujo oscilatorio en la capa lmite, debido a la coexistencia de

3.4.- Mecanismos de transporte de sedimentos

3.5.- Interaccin Fluido-Sedimento

Los mtodos predictivos no incluyen todos los procesos o mecanismos

As, cuando se hacen predicciones de la dinmica de la zona cercana a la costa, es

3.5.1.- Relacin entre la hidrodinmica y los mecanismos de transporte

3.5.2.- Naturaleza de los flujos en la zona cercana a la costa

Figura 3.4.- Desarrollo de una capa lmite (Liu, 2001).

El desarrollo de la capa lmite puede expresarse como:

La capa lmite turbulenta se compone de tres subcapas; la subcapa viscosa, la capa de

donde es la densidad del fluido, es la viscosidad cinemtica ( 10-6 m2 /s) y u es la

Por lo tanto, la capa superior acelerar al fluido proveniente de la capa inferior. De

Por lo tanto, el producto de la magnitud de las fluctuaciones de velocidad turbulentas

Tanto fsicamente como dimensionalmente, la viscosidad turbulenta se puede

donde es la densidad del fluido, u es la velocidad horizontal, p es la presin en el

que es lo que se conoce como la ley de la pared.

donde u* es la velocidad de corte. Por lo tanto, la escala de fluctuaciones de velocidades

donde la viscosidad turbulenta es :

donde k es la constante de von Karman ( 0.41), y

es la llamada velocidad de corte, donde b es la tensin de corte en el fondo (z = 0).

Capa de transicin (tambin llamada capa buffer) donde la viscosidad y la

mucho mas pequea que el grosor de la subcapa viscosa. Por lo tanto, la

Flujos moderadamente turbulentos para

Flujos altamente turbulentos para

Flujos de transicin para 5

afectada por la rugosidad del fondo y por la viscosidad.

En la capa turbulenta, las tensiones de corte se deben exclusivamente a la aportacin de

De acuerdo con la teora de la longitud de mezcla de Prandtl:

y suponiendo que la longitud de mezcla es:

donde k es la constante de vo n Karman (k = 0.4), podemos obtener:

donde el significado fsico de la constante de integracin z0 es la elevacin

Nikuradse (1933) tiene los siguientes valores:

z 0 = n para flujos altamente turbulentos * n 70

+ n para flujos de transicion 5 < * n < 70

donde k n es la rugosidad de equivalente de grano de Nikuradse. Es interesante notar que

Figura 3.9.- Distribucin de las componentes de las tensiones de corte en el fondo.

En los estudios relacionados con la cuantificacin del transporte de sedimentos, es

Aplicando la definicin que u

= 0 e integrando con respecto de z, se obtiene:

donde u*c = bc es la velocidad de corte de la corriente. Por lo tanto, suponiendo un

donde U es la velocidad verticalmente integrada, h es la profundidad del agua y es el

U ( z ) = 1.07U para0.5 < z < h

La rugosidad del fondo k n , conocida tambin como la rugosidad equivalente de arena

Anteriormente, Colebrook & White (1937) propusieron una versin simplificada de la

White (1937) se reduce a una forma idntica a las ecuaciones [3.19]:

que para un flujo altamente turbulento (para

mientras que si el flujo es moderadamente turbulento, (para