Pauta E1

PAUTA EJERCICIO 1
CI71Q-CI6111 HIDRODINAMICA AMBIENTAL

Prof.: Y. Niño
Prof. Auxiliar: O. Sepúlveda
Sem. Otoño 2013
ENERGIA CINETICA TURBULENTA, VISCOSIDAD DE REMOLINOS Y TASA DE

DISIPACION EN EL FLUJO EN UN CANAL ABIERTO.
Partiendo desde las ecuaciones RANS y simplificando para un flujo en canales, permanente, uni-
forme, se obtienen las siguientes ecuaciones:
Continuidad
∂u ∂v
+ =0 (1)
∂x ∂z
Como el flujo es uniforme, se tiene que:
∂u
=0
∂x
por lo tanto, se obtiene que v = 0.
Movimiento
x
b
µt (∂u/∂z)
  
1 ∂ pb 1  ∂  ∂u z }| {
0=− +  µ − ρu0 v 0  (2)

ρ ∂x ρ ∂z ∂z
zb

1 ∂ pb 1 ∂ 02
0=− + −ρv ) (3)
ρ ∂z ρ ∂z
Si se considera la teorı́a de longitud de mezcla de Prandtl para el cierre dela turbulencia, nos queda
que (2) se puede expresar como:
  
esfuerzos viscosos
z}|{ 2 
1 ∂ pb 1 ∂ 
 ∂u ∂u
0=− +   µ + ρl2 (4)

ρ ∂x ρ  ∂z  ∂z ∂z


| {z }
esfuerzos turbulentos
1
CI71Q-CI6111 Hidrodinámica Ambiental
1. ANALISIS DE LA REGION CERCANA A LA PARED
Suponiendo que dentro de la subcapa viscosa los esfuerzos de Reynolds son despreciables
comparados con los esfuerzos viscosos y que en esta región el esfuerzo total puede ser
considerado constante e igual a τb , demuestre que la velocidad media en el sentido del
flujo, ū, está dada por:
ū z u∗
u+ = = z+ =
u∗ ν
p
Tomando sólo los esfuerzos viscosos (τ = µ∂u/∂y) y la expresión u∗ = τb /ρ en el fondo
del canal, nos queda:
∂u
u2+ = ν
∂z
Integrando y considerando que u(0) = 0, tenemos:
u2∗ z = νu
Lo que es equivalente a:
zu∗ u
=
ν u∗
Usando la teorı́a de longitud de mezcla y suponiendo que en una región más arriba los
esfuerzos de corte viscosos son despreciables comparados con los esfuerzos de corte de
Reynolds y que en esa región el esfuerzo de corte total todavı́a es constante e igual a τb ,
demuestre la ley logarı́tmica de velocidades.
Tomando de los esfuerzos turbulentos de la ecuación (4) y utilizando l = κz (Von Kar-

man), nos queda que:
∂u 2

τt = ρ κ z
∂z
p
Reemplazando u∗ = τ /ρ se obtiene que:
∂u u∗
=
∂z κz
du 1 dz
=
u∗ κ z
Se tiene que z+ = zu∗ /ν al mismo tiempo se tiene que dz/z = dz+ /z+ , por lo que queda:
du 1 dz+
=
u∗ κ z+
Integrando se obtiene la ley logarı́tmica de velocidades.
Departamento de Ingenierı́a Civil 2 Universidad de Chile

Suponga que en la región de equilibrio de la turbulencia, donde la ley logarı́tmica es

válida y la producción y disipación de la energı́a cinética turbulenta están en balance,
tal que: P ≈ . Usando el modelo de una ecuación y la teorı́a de longitud de mezcla,
demuestre que en esta región:
K 1 ν 1 νt
2
= 2 (a) ; 4
= (b) ; = κ z+ (c)
u∗ α u∗ κ z+ ν
asuma que las constantes CD y α del modelo de una ecuación siguen la relación: CD = α3 .
√
a) Tomando el modelo de una ecuación: νt = α Kl y elevando al cuadrado, nos queda:
νt2 = α2 Kl2
Reemplazando que l2 = νt (∂u/∂z)−1 se llega a:
∂u
νt = α2 K
∂z
Por otra parte, la derivación de la ley logarı́tmica respeto a z nos estrega la siguiente
ecuación:
∂u
u∗ = l
∂z
2
Lo que junto a νt = l (∂u/∂z), entrega que:
∂u
νt = u2∗
∂z
Ası́,
K 1
u2∗ = α2 Kl2 ⇔ 2
= 2
u∗ α
b) Según el modelo de una ecuación ε, está dado po:
K 3/2 K 3/2
ε = CD ⇔ ε = CD
l κz
Tomando la ecuación (a) y elevándola a ( )3/2 , se obtiene:
K 3/2 1 K 3/2 1
= ⇔ =
u3∗ α3 u3∗ CD
Reemplazando se obtiene:
u3∗ u3∗ ν 1
ε= ⇔ ε= u∗ ⇔ 4
=
κz κz+ ν u∗ κ z+
c) Derivando en z la ley logarı́tmica se llega a que:
∂u
u∗ = l
∂z
reemplazando en el esquma de longitud de mezla, y agregando la hipòtesis de Von
Karman queda que:
νt = κzu∗
Luego, dividiendo esto último por ν se obtiene el resultado.

2. ANALISIS DE LA REGIÓN EXTERIOR

En esta región las escalas relevantes son la altura del flujo, H, y la velocidad de corte u∗
(outer scales). Un perfil vertical de la viscosidad de remolinos puede ser determinado por la
teorı́a de longitud de mezcla suponiendo que ella es válida a través de toda la columna de
agua.
Demuestre que en tal caso:

νt z z
= κ (1 − )
u∗ H H H
(Perfil parabólico de viscosidad de remolinos)
Tomando un perfil de corte lineal, se tiene que:

Z
τ = τb 1 −
H
Luego reemplazando que taub = ρu2∗ y que τt = νt ∂u/∂z Nos queda:

2 Z ∂u
ρu∗ 1 − = νt
H ∂z
Se demostró anteriormente que ∂u = ∂z = u∗ /l y con l = κz (Von Karman), se llega a

que:
2 Z u∗
u∗ 1 − = νt
H κz
Multiplicando lo anterior por z/H, se llega a que:
νt z z
=κ 1−
u∗ H H H
item Usando el modelo de una ecuación y el perfil parabólico de νt , determine relaciones

similares para la energı́a cinética turbulenta adimensional: K/u2∗ y la tasa de disipación
adimensional: H/u3∗ , como funciones de z/H.
El modelo de una ecuación no estrega las siguientes expresiones para nut y ε:

√ K 3/2
νt = α Kl ; ε = CD
l
Reemplazando el perfil parabólico en la primera ecuación, se llega a que:
K 1 z 2
= 1 −
u2∗ α H
Reemplazando el perfil parabólico en la segunda ecuación, se llega a que:

εH 1H z 3
= 1 −
u3∗ κ z H

3. ESCALA DE TIEMPO DE MEZCLA EN LA VERTICAL Obtenga un valor promediado en

la altura del flujo de la viscosidad de remolinos adimensional: νt /(u∗ H). Usando este valor
estime una escala de tiempo para el transporte vertical de masa en la columna de agua.
Promediando el perfil parabólico

Z H Z H
1 νt 1 z z
dz = κ 1− dz
H 0 u∗ H H 0 H H
Se llega, fácilmente a que:

νt κ
=
u∗ H 6
Del análisis dimensional del coeficiente de difusión, se sabe que la escala de tiempo es
del orden:
L2
t∼
D
Como se desea analizar el transporte vertical, para este caso la escala de longitud es H.
Por otra parte, el número de Schmidt turbulento se define como:
νt νt
σt = ⇔ Dt =
Dt σt
Con esto la escala de tiempo para el transporte vertical, del problema estará dado por:
σt
t ∼ H2
νt

Pauta E1

Cargado por

Copyright:

Formatos disponibles

Pauta E1

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

Pauta E1

Cargado por

Copyright:

Formatos disponibles

PAUTA EJERCICIO 1

CI71Q-CI6111 HIDRODINAMICA AMBIENTAL

ENERGIA CINETICA TURBULENTA, VISCOSIDAD DE REMOLINOS Y TASA DE

1. ANALISIS DE LA REGION CERCANA A LA PARED

Tomando de los esfuerzos turbulentos de la ecuación (4) y utilizando l = κz (Von Kar-

Integrando se obtiene la ley logarı́tmica de velocidades.

Departamento de Ingenierı́a Civil 2 Universidad de Chile

Suponga que en la región de equilibrio de la turbulencia, donde la ley logarı́tmica es

Departamento de Ingenierı́a Civil 3 Universidad de Chile

2. ANALISIS DE LA REGIÓN EXTERIOR

Demuestre que en tal caso:

Luego reemplazando que taub = ρu2∗ y que τt = νt ∂u/∂z Nos queda:

Se demostró anteriormente que ∂u = ∂z = u∗ /l y con l = κz (Von Karman), se llega a

item Usando el modelo de una ecuación y el perfil parabólico de νt , determine relaciones

El modelo de una ecuación no estrega las siguientes expresiones para nut y ε:

Reemplazando el perfil parabólico en la segunda ecuación, se llega a que:

Departamento de Ingenierı́a Civil 4 Universidad de Chile

3. ESCALA DE TIEMPO DE MEZCLA EN LA VERTICAL Obtenga un valor promediado en

Promediando el perfil parabólico

Se llega, fácilmente a que:

Departamento de Ingenierı́a Civil 5 Universidad de Chile

También podría gustarte