5.2.

1 LEY LOGARÍTMICA

Los perfiles de velocidad turbulentos difieren mucho del parabólico. El perfil para
la placa plana para flujo turbulento es aproximadamente logarítmico, con una
estela débil en la región exterior y una subcapa viscosa delgada. Esto fue
demostrado por C.B. Millikan en 1937. Por tanto, se supondrá que la ley
logarítmica (5.25) es válida en todo el espesor de la capa límite.
1 /2
u 1 y u¿ τ
u k
¿ = log
v
+ B u¿= 0
ρ( ) ( ) (5.25)

Donde k=0.41 y B=5.

La ecuación de Von Karman escrita en función del espesor de cantidad de
movimiento queda de la siguiente forma:
dθ
τ 0=ρ U 2 (5.26)
dx
Donde el espesor de cantidad de movimiento θes:
δ
u u
θ=∫
0 U (
1− dy (5.27)
U )
La resistencia por fricción de la placa plana D ( x ) es:
D ( x ) =ρb U 2 θ(5.28)
O igualmente:
x
D ( x ) =b∫ τ 0 ( x ) dx (5.29)
0

Siendo b la anchura de la placa plana.

Se inician los cálculos con θ, teniendo en cuenta el perfil logarítmico 5.25.
δ
u u
θ=∫
0 U
1− dy
U ( )
1 y u¿ ¿ 1 y u¿

θ=∫
δ
u
¿
( ( ) )( ( ( ) ))
k
log
v
+B
1−
u
k
log
v
+B
dy
0 U U
2
¿ 1 y u¿

θ=∫
δ

0
( 1 ¿
u log
k
y u¿

U
v ( )
+ Bu ¿

−
u
k
((
log

U
v
+B ( ) )
)) dy
Se separan las integrales.
2
1 y u¿ ¿1 y u¿
θ=∫
u¿ log
k δ
v
+ B u¿ ( )δ u
dy−∫
k
log
v
(
+ B u¿
dy
( ) )
0 U 0 U
2
¿1 y u¿ ¿1 y u¿
δ
θ=∫
u
k
log
v δ
dy+∫
( )
B u¿
δ
dy−∫
u
k
log
v
( ( )) dy +¿ ¿
0 U 0 U 0 U

¿1 y u¿
δ u
k
log
v ( )
Bu ¿ δ
Bu ¿ 2
−∫ 2
0 U2
dy −∫
0 U
dy ( )
Se integra por separado.
La integral del primer sumando se realiza por partes.
1 y u¿ ¿1 ¿1
δ u¿ log
k v
u ( )
k
δ
yu ¿ u
k δ u¿
∫
0 U
dy =
U 0
∫ log
v
dy=¿
U ( )
δlog
v
−δ ¿ ( ( ) )
Ahora el segundo sumando.
δ
Bu ¿ B u¿
∫ U dy = δ
0 U
Y el tercer sumando.
2 2
1 y u¿ 1
δ

( u¿ log
k v ( ) ) ( ) k
u¿ δ
y u¿
2
−∫
0 U
dy=−
U
∫
0
( ( ))
log
v
dy=¿ ¿

2
1
¿−
U( )[ (
k
u¿
δ log
v
2

( )) ( ( ) )]
δ u¿
−2 δlog
δ u¿
v
−δ

El cuarto sumando.
1 y u¿
δ
−∫ 2
u log
k
¿
v
Bu ¿ ( )
dy =¿ ¿
2
0 U
¿−2¿ ¿
El quinto sumando.
 δ
B u¿ 2 B u¿ 2
−∫
0
( )
U
dy=¿
U ( )
δ¿

Por tanto, el resultado final de la integral es:

1
u¿ ¿ ¿
k δu Bu
θ=
U
δlog
v [ ( ) ]
−δ +
U
δ +¿

2
1 y u¿

(
u ¿ log
k v ( ) )[ ( 2
−
U
δ log ( )) ( ( ) )]
y u¿
v
−2 δlog
δ u¿
v
−δ +¿

−2 ¿ ¿
dθ
Ahora se evalúa para poder sustituir en la ecuación 5.26.
dx
1 ¿ δ u¿ 1 ¿ Bu¿ 1 2
θ=
kU
u δlog
v ( )
−
kU
u δ+
U
δ−
kU
¿ ( )
1 2
+2 ( )
kU
¿

+2B
¿
k U2
Primer sumando.
' '
1 ¿ δ u¿ 1 ¿' δ u¿ 1 ¿ δ u¿
[ kU
u δlog
v
= ( )]
kU
u δlog
v
+
kU
u δlog ( )
v
=¿ ( ( ))
1 ¿' δ u¿ 1 ¿ ' δ u¿ 1 ¿' 1 ¿ '
kU
u δlog
v
+
kU ( )
u δ log
v
+
kU
u δ+ ( )
kU
u δ

Segundo sumando.

−1 ¿ ' −1 ¿' 1 ¿ '

( kU
uδ =
kU )
u δ−
kU
u δ

Tercer sumando.

Bu¿ ' B ¿' B

( U U )
δ = u δ + u¿ δ'
U

Cuarto sumando.
¿¿
 2
1 2 ¿ ¿' u¿ δ 1 2
¿−
kU( )
2 u u δ log
v [ ( )] − ( )
kU
¿

2
1 2 ¿ ¿' u¿ δ 1 2
¿−
kU( )
2u u δ log
v [ ( )] − ( )
kU
¿

1 2 ¿ u¿ δ ¿' 1 2 ¿ u¿ δ ¿ '
−( ) kU
u 2 log
v ( )
u δ−
kU ( )
u 2 log
v
uδ ( )
Quinto sumando.
¿¿

1 2 ¿ ¿' u¿ δ 1 2
¿4 ( )
kU
u u δlog ( ) ( )
v
+2
kU
¿

2 ¿ 2
1 u δ 1
¿4( ) u u δlog ( )+ 2(
kU )
¿ ¿'
¿
kU v
2 2
1 1
+2 ( ) u u δ+2 (
kU )
¿ ¿'
¿
kU
Sexto sumando.
¿¿
Séptimo sumando.
¿¿

2B u¿ δ 2B
¿−
kU 2
2u ¿ ¿'
u δlog ( )
v
−
k U2
¿

Octavo sumando.
¿¿
Noveno sumando.
¿¿
Reagrupando todo, queda:

dθ 1 ¿' δ u¿ 1 ¿ ' δ u¿ 1 ¿' 1 ¿ ' 1 ¿'

=
dx kU
u δlog
v
+
kU ( )
u δ log
v
+
kU
u δ+ ( )
kU
u δ−
kU
u δ+¿

2
−1 ¿ ' B ¿' B ¿ ' 1 2 ¿ ¿' u¿ δ
kU
u δ + u δ+ u δ −
U U kU ( )
2u u δ log
v [ ( )] +¿
 1 2
−( ) kU
¿

2 ¿ 2 ¿
1 uδ 1 u δ
u 2 log (
v )
u u δlog (
v )
−(
kU )
u δ +4 (
kU )
¿ ¿ ' ¿ ¿'
+¿

2
1
+2 (
kU )
¿

2 2
1 1
−2 ( ) 2 u u δ−2 (
kU )
¿ ¿'
¿
kU
−2 B
¿
k U2
+2B
¿
k U2
Se simplifica y se saca factor común.
2
dθ ¿'
dx
=u δ
1
kU
log
δ u¿ B
v
+ −
U kU[
1 2 ¿
2 u log
u¿ δ
v ( ) ( ) [ ( )] +2
1 2 ¿
( )
kU
u log
u¿ δ
( ) ( )
v
−2
1 2 ¿ 2B
kU
u−
kU 2
2u ¿
log
u¿ δ
( )
v
+
2
kU
2

+u δ ¿1
kU
'
log
[
δ u¿ B
v
+ −
U
1 2 ¿
kU
u log
u¿ δ
v ( ) ( ) [ ( )] −
2B ¿
k U2
u log
u¿ δ
v
−
B 2 ¿
( )( )
U ]
u (5.31)

Ahora se utiliza la ecuación empírica de Blasius 5.32 como línea de fricción, pues
∂u
al intentar utilizar por definición τ 0= μ
∂ y y=0 ( )
evaluado con el perfil logarítmico el
resultado sería ∞ .
−1/ 4
Uδ
c f =0.046 ( ) v
(5.32)

Que relacionando τ 0 con c f puede ser escrita también como:

1/ 4
v
τ 0=0.0233 ρ U
Uδ
2
( ) (5.33)

Se despeja δ de y se sustituye en la ecuación (5.31), también se sustituye u¿ , k y B.

1/ 2
τ
u= 0
¿
ρ ( ) ; k=0.41 , B=5(5.34)

v
δ =0.0233 4 ρ 4 U 7 (5.35)
τ 40
Así:
4 4 7 '
dθ 0.0233 ρ U v τ 0
= 1/ 2 9/ 2
¿
dx 2 ρ τ0
1 2
1 2 τ0 0.02334 ρ 4 U 7
−
kU ( ) ( )[ (
2
ρ
2
log
ρ1/ 2 τ 7/0 2 )] +¿

1
1 2
τ0 0.02334 ρ4 U 7 1 2
τ 0 12
+2
kU( )( ) ( ρ
2
log 1/ 2 7/ 2
ρ τ0
−2
kU) ( )( ) ρ
+¿

1 1 1
−2 B τ 0 0.02334 ρ4 U 7 2 B τ0 2 B 2 τ0 2
k U2
2 ( ) (
ρ ) 2
log
ρ 1/2 τ 70 /2
+
k U2 ρ
−
U
2
ρ ( ) ( )( )
¿

1 '
τ τ
−4 ( ) 0.0233 ρ U
0 2 4 4 7 0
¿ 5
ρ τ 0

1 2 1
2
τ0 0.0233 4 ρ4 U 7 2 B τ0 0.02334 ρ4 U 7
−
1
kU ( )( ) [ ( ρ
2
log
ρ1 /2 τ 70 /2 )] −
k U2 ρ ( ) ( 2
log
ρ1/ 2 τ 7/0 2 )
B 2
τ 0 12
−
U ( )( ) ρ
¿ (5.36)

Se saca factor común τ '0 y se sustituye en la ecuación de Von Karman (5.26).

' τ0
τ 0= /¿
ρU 2
1 2
1 2 τ0 0.02334 ρ 4 U 7 2
τ 0 12 0.0233 4 ρ4 U 7
−
kU ( ) ( )[ (
2
ρ
2
log
)] ( ) ( ) (
ρ1/ 2 τ 7/0 2
+2
) 1
kU ρ
log
ρ1 /2 τ 70 /2
+¿

τ 0 12 2 B
2 1 1 1
1 τ0 0.0233 4 ρ4 U 7 2 B τ 0 2 B 2 τ 0
−2
kU( )( ) ρ
−
kU 2
2
ρ ( ) (
2
log
) ( ) ( )()
ρ1 /2 τ 70 /2
+
k U2 ρ
−
U
2
ρ
2
¿

1
τ0 0.02334 ρ4 U 7
−4 ( )
ρ
2
5
τ0
¿

1 2 1
2
τ0 0.0233 4 ρ4 U 7 2 B τ0 0.02334 ρ4 U 7
−
1
kU ( )( ) [ ( ρ
2
log
ρ1 /2 τ 70 /2 )] −
k U2 ρ ( ) ( 2
log
ρ1/ 2 τ 7/0 2
+¿
)
B 2
τ 0 12
−
U ( )( ) ρ
¿ ¿ (5.37)
Para resolver esta ecuación diferencial ordinaria se va a utilizar el programa
Matlab, y en concreto la rutina OD45.

5.2.2 LEY POTENCIAL

Como muestra la Figura 5.1, los perfiles turbulentos difieren mucho del perfil
laminar. En este apartado, se van a ajustar los datos del perfil de velocidades con
una ley de potencias, suponiendo que esta ley es válida en todo el espesor de la
capa límite. La forma de ley de potencias es:

7 Rn<1 07

{ }
1 /n
u y
U
= ()
δ
n= 8 1 07 < Rn<1 08 (5.38)
8
9 1 0 < Rn<1 0
9

Figura 5.1. Comparación de los perfiles adimensionales de velocidad de la capa límite.

Ahora se puede aplicar la ecuación integral de von Karman siguiendo los pasos
que se describieron para el flujo laminar, excepto en el momento de evaluar el

∂u
esfuerzo cortante en la pared. La forma de ley de potencias da τ 0=μ =∞ ; para
∂y
y=0, por tanto, el perfil no produce buenos resultados para el esfuerzo cortante

cerca de la pared. Así, en lugar de usar τ 0= μ ( ∂∂ uy ) y=0

, se utilizará la relación

empírica de Blasius que relaciona el coeficiente de fricción superficial local con el

espesor de la capa límite:
−1/ 4
Uδ
c f =0.046 ( ) v
(5.39)

O bien, relacionando τ 0 con c f :

1/ 4
v
τ 0=0.0233 ρ U 2 ( )
Uδ
(5.40)

5.2.2.1 LEY DE POTENCIA (n=7)

Se parte del perfil de velocidades para Reynolds menores a 10 7, de acuerdo a la

ley de potencias.
1 /7
u y
U
=
δ () (5.41)

Se sustituye el perfil de velocidades en la ecuación de Von Karman.

δ 1/7 1 /7
d
τ 0= ∫ ρU 2 δy
dx 0 () [ ()]
1−
y
δ
dy=
7
72
ρ U2
dδ
dx

Al combinar la expresión de Blasius con la anterior se tiene:

1 /4
7 dδ v
72
ρ U 2 =0.0233 ρ U 2
dx Uδ ( )
Eliminando los términos de velocidad del flujo exterior al cuadrado y densidad y
luego reagrupando:
1/4
v
δ 1/ 4 dδ=0.2396 ( )
U
dx

Se integran ambos términos.

1/4
5 5/ 4 v
4
δ =0.2396 x
U ( )
Se despeja δ.
 1 /5
v
δ =0.3812 x 4 /5 ( )
U

1
Como x 4 /5 =x /(x 5 ), la expresión se puede poner de la siguiente forma:
1/ 5
v
δ =0.3812 x
Ux ( )
O lo que es lo mismo:

δ =0.3812 x ℜx−1 /5 (5.42)

Ahora se evalúa el esfuerzo cortante en la pared τ 0, los coeficientes de fricción

local c f y medio C F y la fuerza de arrastre D F.
Al sustituir la expresión anterior para δ en la ecuación de Blasius queda:
c f =0.059 3 ℜ x−1/ 5 (5.43)
Integrando para el cálculo del coeficiente medio de fricción:
L
1
CF= ∫ c ( x ) dx (5.44)
L 0 f
L L −1/ 5
1 0.0593 Ux
C F = ∫ 0.0593 R n−1/
L 0 x
5
dx=
L
∫
0 v ( ) dx=¿ ¿

L L
0.0593U −1/ 5 0.0593 U −1/ 5 5 4 /5
¿
L v−1 /5
∫
0
x−1/ 5
dx=
L v−1 /5 4
x [ ] 0

5
C F =0.0741 R n−1/
L (5.45)
Para el esfuerzo cortante τ 0, al sustituir la expresión para δ, queda:

τ 0=0.0296 ρ U 2 R n−1/ 5
x (5.46)
Y para la fuerza de arrastre, integrando se obtiene:
L
D F=∫ τ 0 ( x ) Edx ( 5.47 )
0

L −1 L −1
2 5 2 5
D F=E ∫ 0.0296 ρ U R n x dx=0.0296 ρ U E ∫ R n x dx=¿ ¿
0 0

L L
0.0296 ρ v 1 /5 U 2 E 5 4 /5 0.0371 ρ v 1 /5 U 2 L4 / 5 E
−1 /5
Ux
¿ 0.0296 ρ U 2 E ∫
0
( )v
=
U 1/5 [ ]
4
x
0
=
U 1/5
 1
Para dejar la ecuación en función de Rn, L4 /5 =L/(L 5 ) y la ecuación queda:
1 /5
v
D F=0.0371 ρ U 2≤ ( )UL

D F=0.0371 ρ U 2 LER n−1/ 5

L (5.48)

5.2.2.2 LEY DE POTENCIA (n=8)

Se parte del perfil de velocidades para Reynolds menores a 10 8, de acuerdo a la

ley de potencias.
1 /8
u y
U
=
δ () (5.49)

Se sustituye el perfil de velocidades en la ecuación de Von Karman.

δ 1/ 8 1 /8
d
τ 0= ∫ ρU 2 δy
dx 0 () [ ()]
1−
y
δ
dy=
8
90
ρU 2
dδ
dx

Al combinar la expresión de Blasius con la anterior se tiene:

1/ 4
8 dδ v
90
ρ U 2 =0.0233 ρ U 2
dx Uδ ( )
Eliminando los términos de velocidad del flujo exterior al cuadrado y densidad y
luego reagrupando:
1/ 4
v
δ 1/ 4
dδ=0.2 621
U ( ) dx

Se integran ambos términos.

1/ 4
5 5/ 4 v
4
δ =0.2621 x
U ( )
Se despeja δ.

1/ 5
v
δ =0. 4096 x 4 / 5 ( )
U

1
Como x 4 /5 =x /(x 5 ), la expresión se puede poner de la siguiente forma:
1 /5
v
δ =0. 4096 x
Ux ( )
O lo que es lo mismo:

δ =0. 4096 x ℜ x−1/ 5( 5.50)

Ahora se evalúa el esfuerzo cortante en la pared τ 0, los coeficientes de fricción

local c f y medio C F y la fuerza de arrastre D F.
Al sustituir la expresión anterior para δ en la ecuación de Blasius queda:
c f =0.0583 ℜx−1/ 5 (5.51)
Integrando para el cálculo del coeficiente medio de fricción:
L
1
CF= ∫ c ( x ) dx
L 0 f
L L −1 /5
1 0.058 3 Ux
C F = ∫ 0.0583 R n−1/
L 0 x
5
dx =
L
∫
0 v ( ) dx=¿ ¿

L L
0.0583 U −1/ 5 0.058 3 U −1 /5 5 4/ 5
¿
L v−1 /5
∫
0
x−1/ 5
dx =
L v−1/ 5 4
x [ ] 0

C F =0.0728 R n−1 /5
L (5.52)
Para el esfuerzo cortante τ 0, al sustituir la expresión para δ, queda:

τ 0=0.0291 ρU 2 R n−1 /5
x (5.53)
Y para la fuerza de arrastre, integrando se obtiene:
L
D F=∫ τ 0 ( x ) Edx
0

L −1 L −1
2 5 2 5
D F=E ∫ 0.0291 ρ U R n x dx=0.0291 ρ U E∫ R n x dx =¿ ¿
0 0

L −1/ 5 L
Ux 0.0291 ρ v 1/5 U 2 E 5 4 / 5 0.0364 ρ v 1/ 5 U 2 L4 /5 E
¿ 0.0291 ρU 2 E ∫
0
( )
v
=
U 1 /5 4 [ ]
x
0
=
U 1 /5
1
Para dejar la ecuación en función de Rn, L4 /5 =L/(L 5 ) y la ecuación queda:
1/ 5
v
D F=0.0364 ρU 2≤ ( )
UL

D F=0.0364 ρU 2 LER n−1 /5

L (5.54)

5.2.2.3 LEY DE POTENCIA (n=9)

Se parte del perfil de velocidades para Reynolds menores a 10 9, de acuerdo a la
ley de potencias.
1 /9
u y
U
=
δ() (5.55)

Se sustituye el perfil de velocidades en la ecuación de Von Karman.

δ 1/ 9 1 /9
d y y 9 dδ
τ 0= ∫ ρU 2
dx 0 δ () [ 1− ()
δ ] dy=
110
ρU 2
dx

Al combinar la expresión de Blasius con la anterior se tiene:

1/ 4
9 dδ v
110
ρ U 2 =0.0233 ρU 2
dx Uδ ( )
Eliminando los términos de velocidad del flujo exterior al cuadrado y densidad y
luego reagrupando:
1 /4
v
δ 1/ 4 dδ=0.2 848 ( ) U
dx

Se integran ambos términos.

1 /4
5 5/ 4 v
4
δ =0.2848 x
U ( )
Se despeja δ.

1/ 5
v
δ =0. 4377 x 4 / 5 ( )
U

1
Como x 4 /5 =x /(x 5 ), la expresión se puede poner de la siguiente forma:
1 /5
v
δ =0. 4377 x ( )
Ux

O lo que es lo mismo:

δ =0. 4377 x ℜ x−1/ 5( 5.56)

Ahora se evalúa el esfuerzo cortante en la pared τ 0, los coeficientes de fricción

local c f y medio C F y la fuerza de arrastre D F.
Al sustituir la expresión anterior para δ en la ecuación de Blasius queda:
 c f =0.057 3 ℜx−1/ 5 (5.57)
Integrando para el cálculo del coeficiente medio de fricción:
L
1
C F = ∫ c f ( x ) dx
L 0
L L −1/ 5
1 0.057 3
CF= ∫ 0.057 3 R n−1/
L 0 x
5
dx=
L
∫ Ux
v ( ) dx =¿ ¿
0

L L
0.057 3U −1/ 5 0.057 3 U −1/ 5 5 4 /5
¿
Lv
−1 /5 ∫
0
x−1/ 5
dx=
Lv
−1 /5
4
x [ ] 0

C F =0.0716 R n−1 /5
L (5.58)
Para el esfuerzo cortante τ 0, al sustituir la expresión para δ, queda:

τ 0=0.0286 ρ U 2 R n−1/ 5
x (5.59)
Y para la fuerza de arrastre, integrando se obtiene:
L
D F=∫ τ 0 ( x ) Edx
0

L −1 L −1
D F=E ∫ 0.0286 ρ U 2 R nx5 dx=0.0286 ρ U 2 E ∫ R n x5 dx=¿ ¿
0 0

L L
0.0286 ρ v 1 /5 U 2 E 5 4 /5 0. 0358 ρ v 1 /5 U 2 L4 / 5 E
−1 /5
Ux
2
¿ 0.0286 ρ U E ∫
0
( )
v
=
U
1/5 [ ]
4
x
0
=
U
1/5

1
Para dejar la ecuación en función de Rn, L4 /5 =L/(L 5 ) y la ecuación queda:
1 /5
v
D F=0.0358 ρ U 2≤ ( )
UL

D F=0.0358 ρ U 2 LERn−1/ 5
L (5.60)