0% encontró este documento útil (0 votos)
72 vistas26 páginas

Soluciones Exactas a Ecuaciones de Navier-Stokes

Este documento presenta cuatro ejemplos de soluciones exactas a las ecuaciones de Navier-Stokes para flujos viscosos unidimensionales estacionarios: (1) Flujo de corte simple entre dos placas paralelas, con una velocidad lineal; (2) Flujo en una tubería de Poiseuille, con una velocidad parabólica; (3) Película de fluido que escurre sobre una superficie inclinada, con una velocidad lineal; (4) Otros problemas unidireccionales como el flujo cortante de dos lí
Derechos de autor
© Attribution Non-Commercial (BY-NC)
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
72 vistas26 páginas

Soluciones Exactas a Ecuaciones de Navier-Stokes

Este documento presenta cuatro ejemplos de soluciones exactas a las ecuaciones de Navier-Stokes para flujos viscosos unidimensionales estacionarios: (1) Flujo de corte simple entre dos placas paralelas, con una velocidad lineal; (2) Flujo en una tubería de Poiseuille, con una velocidad parabólica; (3) Película de fluido que escurre sobre una superficie inclinada, con una velocidad lineal; (4) Otros problemas unidireccionales como el flujo cortante de dos lí
Derechos de autor
© Attribution Non-Commercial (BY-NC)
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Cap tulo 3 Flujo Viscoso: Soluciones Exactas a NS

3.1. Soluciones exactas de las ecuaciones de Navier-Stokes

Debido a la complejidad de las ecuaciones de N-S, stas solo se pueden e resolver en casos especiales. En general, buscamos que la geometr del ujo a se tal que algunos de las partes de la ecuacin se cancelen. o Lo que se busca en cada problema es encontrar el campo de velocidades: v = f (x, y, z, t). Para toda esta seccin consideraremos unicamente ujos newtonianos, o incompresibles, isotrmicos y de propiedades constantes. e Para esto consideraremos las siguientes ecuaciones: Dv = f Dt v = 0 P +
2

(3.1) (3.2)

Tambin consideraremos, para la mayor de los casos, problemas en coe a ordenadas cartesianas tal que v = (u, v, z).

3.1.1.

Flujo de corte simple o de Couette

Flujo entra dos placas paralelas innitas (2-D) Consideremos el ujo en la siguiente gura: 67

68

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS


U
Pared movil

H
y

Pared fija

Para este caso v = (u, v). Consideraremos las ecuaciones 7.1 y 7.1 para el caso 2-D: Ecuacin de continuidad: o u v + =0 x y Ecuaciones de Navier-Stokes: u u u +u +v t x y v v v +u +v t x y P + x P = + y = 2u 2u + x2 y 2 2v 2v + x2 y 2 + gx + gy

Caracter sticas del ujo: 1. Flujo estacionario (no cambia como funcin del tiempo): o /t = 0 2. Flujo desarrollado (no cambia con la posicin x, placas innitas): o /x = 0 3. La gravedad esta alineada con la direccin y y : o g = (0, gy ). 4. No hay gradiente de presin en la direccin x x : o o P/x = 0 Condiciones de contorno:

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES 1. v(x, 0) = (0, 0) 2. v(x, H) = (U, 0)

69

Considerando la caracter stica de ujo (2) en la ecuacin de conservacin o o de masa, tenemos: v =0 0+ y Esta expresin implica que o v = constante De la primera condicin de contorno sabemos que v = 0 en y = 0. Puesto o que v es constante, podemos inferir que v=0 en todas partes. Esto es un resultado muy importante. Signica que el ujo es unidireccional : v = (u, 0). Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin y y . Si v = 0 en todas partes, entonces los unicos trminos que o e sobreviven son: (0 + u(0) + 0(0)) = Reescribiendo, P + (0 + 0) + gy y

P = gy y que es la ecuacin de la hidrosttica. o a Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin x x : o u 0 + u(0) + 0( ) y entonces: 2u = (0) + 0 + 2 y + (0)

2u 0= 2 y Esta expresin se puede integrar directamente: o u = Ay + B

70

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

Sabemos de las dos condiciones de contorno que u = 0 en y = 0 y u = U en y = H: B=0 y A= Por lo tanto U H

U y H Es una distribucin de velocidades lineal: o u(y) =


U
Pared movil

H
y

u(y)=U/H y

Pared fija

Otra resultado que podemos obtener de este resultado es el esfuerzo cortante en la pared. La relacin constitutiva newtoniana es o ij = P ij + ij Dkk + 2Dij El esfuerzo cortante en la pared es xy |y=o : xy = 2Dxy = Para el ujo cortante simple: w = u U = y H u v + y x

Ntese que el esfuerzo es constante a travs de todo el campo de ujo. o e

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES


P1
y
Pared fija

71
P2

H/2

H/2

Pared fija

3.1.2.

Flujo en una tuber o de Poiseuille a

Consideremos el ujo en la siguiente gura: Para este caso v = (u, v). Consideraremos las ecuaciones 7.1 y 7.1 para el caso 2-D: Ecuacin de continuidad: o u v + =0 x y Ecuaciones de Navier-Stokes: u u u +u +v t x y v v v +u +v t x y Caracter sticas del ujo: 1. Flujo estacionario (no cambia como funcin del tiempo): o /t = 0 2. Flujo desarrollado (no cambia con la posicin x, placas innitas): o /x = 0 3. La gravedad esta alineada con la direccin y y : o g = (0, gy ). 4. El gradiente de presin en la direccin x x es: o o P/x (P1 P2 )/L = G Esta cantidad es constante. Esta cantidad se dene negativa porque para que exista un ujo de izquierda a derecha la presin P2 debe se mas baja que P1 . o P + = x P = + y 2u 2u + x2 y 2 2v 2v + x2 y 2 + gx + gy

72

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS Condiciones de contorno:

1. v(x, H/2) = (0, 0) 2. v(x, H/2) = (0, 0) Paredes jas. Considerando la caracter stica de ujo (2) en la ecuacin de conservacin o o de masa, tenemos: v 0+ =0 y Esta expresin implica que o v = constante De la primera condicin de contorno sabemos que v = 0 en y = 0. Puesto o que v es constante, podemos inferir que v=0 en todas partes. Esto es un resultado muy importante. Signica que el ujo es unidireccional : v = (u, 0). Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin y y . Si v = 0 en todas partes, entonces los unicos trminos que o e sobreviven son: (0 + u(0) + 0(0)) = Reescribiendo, P + (0 + 0) + gy y

P = gy y que es la ecuacin de la hidrosttica. o a Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin x x : o 0 + u(0) + 0( entonces: G = u ) y = P 2u + 0+ 2 x y 2u y 2 + (0)

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES Esta expresin se puede integrar directamente: o u= G 2 y + Ay + B 2

73

Sabemos de las dos condiciones de contorno que u = 0 en y = H/2 y u = 0 en y = H/2: A=0 y B= Por lo tanto u(y) = G 2 H 2
2

G 2

H 2

y2

Es una distribucin de velocidades parablica: o o


P1
y Velocidad
Pared fija

Esfuerzo

P2

H/2

H/2

Pared fija

Podemos calcular la velocidad del uido en el centro de la tuber a: umax = u(y = 0) = G 2 H 2


2

Podemos tambin calcular el campo de esfuerzos cortantes: e xy = G u = y y

74

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

Ntese que el esfuerzo es var linealmente a travs de y. El esfuerzo en la o a e pared es u GH w = = . y 2 Podemos calcular el gasto Q que pasa a travs de la tuber Sabemos e a. que Q=
a

udA

Puesto que es un problema en 2-D, podemos simplicar al elemento diferencial de area como dA = dydz = dy(1), por unidad de z. As tenemos ,
H/2

Q =
H/2 H/2

u(y)dy G 2 H 2
2

=
H/2

y2

dy

G 3 H 12

3.1.3.

Pel cula de uido que escurre sobre una pared inclinada

Consideremos el ujo en la siguiente gura: Para este caso v = (u, v). Consideraremos las ecuaciones 7.1 y 7.1 para el caso 2-D. Ver en los ejemplos anteriores. Caracter sticas del ujo: 1. Flujo estacionario (no cambia como funcin del tiempo): o /t = 0 2. Flujo desarrollado (no cambia con la posicin x, placa y pel o cula innitas): /x = 0 3. La gravedad acta en las dos direcciones x x y y y : u g = (s sen , g cos ). 4. No hay gradiente de presin en la direccin x x : o o P/x = 0

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES

75

gc os

gs

en

superficie libre

Condiciones de contorno: 1. No deslizamiento en la pared: v(x, 0) = (0, 0) 2. Supercie libre de esfuerzo: xy (x, H) = 0 Considerando la caracter stica de ujo (2) en la ecuacin de conservacin o o de masa, tenemos: v = constante De la primera condicin de contorno sabemos que v = 0 en y = 0. Puesto o que v es constante, podemos inferir que v=0 en todas partes. Esto es un resultado muy importante. Signica que el ujo es unidireccional : v = (u, 0). Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin y y . Si v = 0 en todas partes, entonces los unicos trminos que o e sobreviven son: P (0 + u(0) + 0(0)) = + (0 + 0) + g cos y Reescribiendo, P = g cos y

76

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

que es la ecuacin de la hidrosttica. o a Consideremos ahora la ecuacin de conservacin de momentum en la dio o reccin x x : o u 0 + u(0) + 0( ) y entonces: 0= 2u = (0) + 0 + 2 y 2u + g sin y 2 + g sin

Esta expresin se puede integrar directamente: o u= g y2 sin + Ay + B 2

Sabemos de las dos condiciones de contorno que u = 0 en y = 0 y xy = 0 en y = H: B=0 y yx = entonces y A= Por lo tanto u(y) = u (y = H) = 0 y

g u = 0 = sin H + A y g sin H

g y2 sin Hy 2

Podemos calcular la velocidad del uido en la supercie libre: umax = u(y = H) = g sin H2 2

Podemos tambin calcular el esfuerzo en la pared: e w = u g (y = 0) = sin H. y

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES


U

77

H x

3.1.4.

Otros problemas unidireccionales estacionarios

Flujo cortante simple de dos l quidos inmiscibles Resolver el ujo para cada lado. En las paredes considerar condicin de o no-deslizamiento. En la interfaz, el esfuerzo cortante debe ser el mismo. As para el uido 1 (superior), tenemos u1 (y) = C1 y + C2 en y = H, u = U y en y = 0, xy = int . Entonces, U = C1 (H) + C2 y int = 1 u1 = 1 C1 y

Para el uido 2 (inferior), tenemos u2 (y) = C3 y + C4 en y = H, u = 0 y en y = 0, xy = int . Entonces, 0 = C3 (H) + C4 y int = 2 u2 = 2 C3 y

78

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS En la interfaz, int1 = int2 : 1 C 1 = 2 C 3

y u1 (y = 0) = u2 (u = 0), lo que implica que, C2 = C 4 Entonces C1 = C2 C3 C4 Por lo tanto u1 = u2 = y 1 U + + 1 H 2 y U +1 2 +1 H 1


1 2

1 +1 1 = U 2 +1 1 U 1 = H 2 + 1 1 1 = U 2 +1 1 U H
1 2

Podemos calcular el esfuerzo cortante: U 1 = 1 C 1 = 1 H y 2 = 2 C 3 = 2 Notemos que U H

1 2

1 +1 1 +1

2 1

U 1 2 . H 1 + 2 As podr amos calcular una viscosidad efectiva del medio: U = ef xy 2H entonces 1 2 ef = 2 1 + 2 1 = 2 =

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES Flujo de Poiseuille en una tuber circular a

79

Este problema es muy parecido al problema de ujos de Poiseuille entre dos placas (ver seccin 3.1.2). o La unica diferencia es que debemos considerar un sistema de coordenadas cil ndricas: v = (ur , u , uz ) En este caso la ecuacin de conservacin de masa, v = 0, se escribe o o en forma expl cita como: 1 1 (rur ) + u + uz = 0 r r r z Por ser un ujo desarrollado (/z) y (/z) axisimtico (/), tenemos: e 1 (rur ) = 0 r r Por lo tanto: ur = 0 Flujo unidireccional. El vector de velocidad se reduce a v = (0, 0, uz ). Resolviendo la ecuacin de conservacin de momentum unicamente en la o o direccin donde el componente de velocidad no es cero, direccin z z . Para o o coordenadas cil ndricas tenemos: uz uz u uz uz + ur + + uz t r r z = P + z 1 r r r uz r + 1 2 uz 2 uz +gz + r2 2 z 2

Considerando las mismas caracter sticas del ujo que en la seccin 3.1.2 o y adems que el ujo es axisimtrico e unidireccional, la ecuacin anterior se a e o reduce a : G uz 1 = r r r r donde G = P/z = constante Integrado una vez tenemos: uz G C1 = r+ r 2 r Integrado una vez mas: uz = G 2 C1 r + r + C2 4 ln

80

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

Sabemos que la velocidad en r = 0 debe ser nita, por lo tanto C1 = 0. Tambin sabemos que la velocidad del uido en la pared debe de ser cero e (condicin de no deslizamiento): uz (r = R) = 0. o C2 = G 2 R 4

Por lo tanto, el campo de velocidades para una tuber circular bajo un a gradiente de presin constante es: o uz = G R2 r 2 4
R

Podemos calcular el ujo volumtrico como: e Q=


A

uz dA =
0

uz (2rdr)

As : Q= La velocidad media, U = Q/A, es

G 4 R 8

G 2 U= R 8 Podemos calcular el esfuerzo en la pared es: pared = rz |r=R = Entonces, el esfuerzo en la pared es: pared = GR 2 uz r

Podemos calcular el coeciente de friccin sobre la tuber o a: Cf = Ff 1 2 U A 2

El area de contacto A es 2RL. La fuerza de friccin Ff sobre la tuber o a se puede calcular directamente del esfuerzo en la pared como Ff = pared A: Ff = R2 LG

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES Por lo tanto Cf =

81

16 16 = Re (2R)U

donde Re = D U / es el nmero de Reynolds. u Factor de friccin: o Podemos reeinterpretar el resultado del gasto volumtrico como: e G= 8Q R4

como la ca de presin que ocurre el una tuber de radio R bajo el ujo da o a Q de un uido de viscosidad . Si escribimos este resultado como prdida de e carga (altura) tenemos: hf = 8Q 1 P L = R4

El factor de rozamiento f para una tuber es una cantidad usada muy a frecuentemente por INGENIEROS para calcular prdidas en tuber e as. El factor de rozamiento se dene como: hf = f L U2 D 2

donde D = 2R. Entonces, reescribiendo la expresin anterior tenemos: o 2 U D2 1 L D4 /16 U L = 32 2 D 64 L U 2 = U D D 2 64 64 D = Re U

hf =

Por lo tanto f= Diagrama de Moody:

82

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

Figura 3.1: Diagrama de Moody. Tomado de Fox, Macdonald Flujo con lineas de corriente circulares Consideremos el siguiente ujo: Una vez mas consideramos ujo estacionario, bidimensional, axisimtrico: e v = (ur , ut heta, 0) De la ecuacin de conservacin de masa podemos demostrar que la velocio o dad ur = 0, por lo que tenemos un ujo unidireccional (l neas de corriente circulares). Si resolvemos la ecuacin de conservacin de momentum para la direccin o o o r r tenemos: ur ur u ur u2 ur + ur + + uz t r r r z = P + r
2

ur + gr

La mayor de los trminos son cero ya que ur = 0. Sin embargo sobrea e viven: P u2 = r r

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES


Espacio entre los dos cilindros lleno de un fluido viscoso W1 R1 W2

83

R2

Existe un gradiente de presin en la direccin radial como resultado de la o o !fuerza centr fuga! Ahora si, resolviendo para la direccin acimutal : o u u u u u ur u + ur + + + uz t r r r z 1 P r 1 (ru ) + r r r =

1 2 u 2 u + r2 2 z 2

+ g

Eliminando trminos tenemos: e r 1 (ru ) r r =0

Integrando dos veces tenemos: (ru ) = C1 r r C1 C2 u = r+ 2 r Las condiciones de frontera son en r = R1 , u = 1 R1 . en r = R2 , u = 2 R2 . Resolviendo para C1 y C2 tenemos: C1 = C2
2 2 2(2 R2 1 R1 ) 2 2 R2 R 1 2 2 2 1 = R1 R2 2 2 R2 R 1

84

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS As , u =

Aplicacin: viscosimetr o a. Consideremos el caso en que cilindro exterior se mantiene jo y el interior gira a una velocidad angular constante. Para dicho caso el perl de velocidades es: R 2 R 2 1 1 R 2 u = 2 1 2 r + 2 1 2 2 R2 R 1 R2 R 1 r para el cual el campo de esfuerzos cortantes es: r = r u r r R 2 R 2 1 = 2 1 2 2 R2 R 1 r

2 2 2 R2 1 R1 R 2 R 2 2 1 r 21 22 2 2 R2 R 1 R2 R 1 r

Si calculamos el torque sobre el cilindro interno tenemos: T = r AR1 R 2 R 2 1 (2R1 L) R1 = 2 1 2 2 R2 R 1 R1 R2 R2 = 2 1 2 2 21 R1 L R2 R 1 Si medimos el torque en el cilindro interno podr amos utilizar este sistema para medir la viscosidad del uido:
2 2 T R2 R 1 = 2 2 21 R1 L R1 R2

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES

85

3.1.5.

Flujos no-estacionarios

Hasta ahora solo hemos visto soluciones exactas para el caso de ujos estacionarios (/t = 0). Existe tambin una clases de ujos que tienen e soluciones exactas pero que no son estacionarios: ujo unidireccional v = (u, 0, 0) pero u/t = 0. Para ujos newtonianos incompresibles consideramos las ecuaciones: v = 0 Dv = f Dt (3.3) P +
2

(3.4)

Para el caso unidireccional, desarrollado y no estacionario, estas ecuaciones se reducen a: v = 0 y 2u u = G+ t y 2 donde G = P = constante. x Primer problema de Stokes Consideremos el problema de una placa innita, inicialmente en reposo, sobre la cual hay un uido viscoso. Las condiciones de frontera e iniciales para u=u(x,y,t) son: u(x, y, 0 ) = 0 u(x, 0, 0+ ) = U u(x, , t) = 0 Para este problema en particular, consideramos que G = 0, no existe gradiente de presin en la direccin x x . o o De la misma manera que para el caso de ujo uni-direccionales estacionarios, utilizando la ecuacin de continuidad podemos deducir que u = u(y, t). o La ecuacin de conservacin de momentum se reduce a: o o u 2u = 2 t y (3.7) (3.5) (3.6)

86

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS

U
a o o donde = es la viscosidad cinemtica. Esta ecuacin es una ecuacin diferencial partial de segundo orden. Esta ecuacin es la ecuacin de difusin. o o o Existen varios mtodos para resolver este tipo de ecuaciones. e 1. Mtodo de Similaridad e Consideremos la variable adimensional u = g() U donde = y(t)1/2 es la variable de similaridad que tambin es adie mensional. Debemos proceder a hacer el cambio de variable en la ecuacin de o conservacin de momentum: o u u = t t Podemos calcular u = (U g()) = U g () Entonces 1 u = (U g ())( y 1/2 t3/2 ) 2

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES simplicando tenemos u 1 = U t1 g () 2 De la misma manera podemos calcular: u u = = U (t)1/2 g () y y 2u u = = U (t)1 g () 2 y y y

87

Sustituyendo en la ecuacin de conservacin de momentum y simplio o cando tenemos 1 g () + g () = 0 (3.8) 2 la cual es una ecuacin diferencial ordinaria (a diferencia de la ecuacin o o original que era parcial. Para resolverla solo falta traducir las condiciones de frontera: Para y = 0 y t 0, u = U , entonces para = 0, g(0) = 1

Para y en cualquier t, u = 0, entonces para , g() = 0. Podemos integrar la ecuacin 3.8 o 1 g + = 0 g 2 1 ln g + 2 = C1 4 g = C 1 e 4


0
2 2 4

g = C1

d + C2

Con las condiciones de frontera podemos calcular C1 y C2 C2 = 1 C1 =


0

1 e

2 4

1 = d

88

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS Entonces u(y, t) = U Que se puede reescribir como: u(y, t) = U 1 erf y 2(t)1/2

1 1

e
0

2 4

donde erf(x) es la funcin error denida como: o 1 erf(x) =


x

e
0

x2 2

dx

2. Mtodo de Transformada de Laplace Aplicando la transformada de e Laplace L[f (t)] = f (s) = est f (t)dt
0

a la ecuacin 3.7 obtenemos: o u L t 2u = L y 2 2u s = 2 u y

Esta ecuacin es diferencial ordinaria con solucin: o o u(y, s) = C1 e s/y + C2 e s/y Para encontrar C1 y C2 debemos transformar las condiciones de frontera u(0, t) = U u(0, s) = U/s u(y, 0) = 0 u(y, 0) = 0 Considerando que u 0 en y , entonces C1 = 0. Considerando la primera condicin de frontera C2 = U/s. Por lo tanto: o u= U s/y e s

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES Aplicando la transformada inversa de Laplace L1 [f (y, s)] = f (y, t) = 1 2i
a+i

89

est f (s)ds
ai

(donde s = a + i) encontramos la velocidad como: u(y, t) = U erfc = U 3. Separacin de Variables. o Este mtodo se ver en la resolucin del siguiente problema. e a o Segundo problema de Stokes Consideremos ahora el siguiente problema: y 2 t y 2 t

1 erf

U cos nt
La ecuacin a resolver es la misma que para el caso anterior o u 2u = 2 t y Las condiciones de frontera e iniciales para u=u(x,y,t) son:

90

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS u(x, y, 0 ) = 0 u(x, 0, 0+ ) = U cos nt u(x, , t) = 0

Se puede resolver utilizando cualquiera de los dos mtodos anteriores. e Aqu usaremos el mtodo de separacin de variables para demostrar su apli e o cacin. o supondremos que la solucin, u(y, y) es el producto de dos funciones, cada o una unicamente funcin de una de las variables independientes: o u(y, t) = W (y) T (t) Puesto que la perturbacin que se est aplicado a la pared es de tipo o a cos nt podemos esperar que la funcin T (t) sea del mismo tipo. As o u(y, t) = W (y) cos nt Ms aun podemos considerar a u(y, t) = W (y) {eint } Esta expresin la podemos sustituir dentro de la ecuacin de conservacin o o o considerando u = W (y)(in)eint t u = W (y)eint y 2u = W (y)eint y 2 entonces que simplicando resulta en {W (y)(in)eint } = {W (y)eint } W (y) in W (y) = 0

que es una ecuacin diferencial ordinaria. Resolviendo esta ecuacin tenemos o o n n (1+i) (1+i) y y 2 W (y) = C1 e + C2 e 2

3.1. SOLUCIONES EXACTAS A NAVIER-STOKES

91

Considerando las condiciones de frontera, tenemos que para y , u 0, lo que implica que C2 = 0 Entonces, la solucin es o n (1+i) y int 2 e u(y, t) = C1 e n n = C1 e( 2 y) ei(nt 2 y) n n y = C1 e( 2 y) cos nt 2 Tambin sabemos que u(y = 0) = U cos nt, por lo que C1 = U . Entonces e la solucin al problema es o n 1 n y 2 u(y, t) = U e cos nt y 2 Flujo pulstil entre dos placas a Consideremos de nuevo el ujo entre dos placas paralelas, generado por un gradiente de presin: o
P1
y
Pared fija

P2

H/2

H/2

Pared fija

Pero ahora, en vez de considerar el que G = constante consideremos un gradiente de presin que sea funcin del tiempo: o o G = G0 cos nt = [Go eint ] La ecuacin a resolver sera o u 1 2u = G+ 2 t y

92

CAP ITULO 3. FLUJO VISCOSO: SOLUCIONES EXACTAS A NS Las condiciones de frontera e iniciales para u=u(x,y,t) son: u(x, H/2, t) = 0 u(x, H/2, t) = 0 este problema se puede resolver fcilmente usando separacin de variables: a o u(y, t) = [W (y)Go eint ] Que, al aplicarla a la ecuacin a resolver, nos d: o a 2W in Go W = 2 y

que es una ecuacin diferencia ordinaria no-homognea. Esta se puede reo e solver fcilmente: a W (y) = i Go + C1 cosh (1 + i) n 2 + C2 sinh (1 + i) n 2

Para las condiciones de frontera (no deslizamiento en las paredes) tenemos: C2 = 0 C1 = iGo n cosh (1 + i)
n H 2 2

Por lo tanto, el perl de velocidades esta dado por u(y, t) = i Go n 1 cosh (1 + i) cosh (1 + i)
n y 2 n H 2 2

eint

También podría gustarte