IV.- FLUJO INCOMPRESIBLE NO VISCOSO pfernandezdiez.es

€ IV.- FLUJO INCOMPRESIBLE NO VISCOSO pfernandezdiez.es IV.1.- CINEMÁTICA DE FLUIDOS La Cinemática de Fluidos tiene una correspondencia biunívoca co

1 downloads 57 Views 2MB Size

Recommend Stories


VENTILACIÓN MECANICA NO INVASIVA Y OXIGENOTERAPIA DE ALTO FLUJO
VENTILACIÓN MECANICA NO INVASIVA Y OXIGENOTERAPIA DE ALTO FLUJO Modalidades especiales de tratamiento en las infecciones respiratorias agudas bajas co

CAPITULO IV CAPITULO IV
CAPITULO IV ________________________________________________________________________ CAPITULO IV 37 CAPITULO IV __________________________________

Story Transcript



IV.- FLUJO INCOMPRESIBLE NO VISCOSO pfernandezdiez.es

IV.1.- CINEMÁTICA DE FLUIDOS La Cinemática de Fluidos tiene una correspondencia biunívoca con el Primer Principio de la Termodinámica aplicado a sistemas abiertos. En un fluido en movimiento, cada partícula posee una velocidad  V que depende de la posición (x,y,z) de dicha partícula y del tiempo t, es decir: V = f(x, y, z, t) y sus proyecciones sobre los tres ejes son función también de dichas variables, viniendo representadas ⎧u = u ( x , y, z, t ) por: ⎨v = v ( x, y, z, t ) ⎩w = w ( x , y, z, t ) Se llama movimiento permanente o estacionario a aquel en que sus características, como la presión, velocidad, etc, son independientes del tiempo, es decir, son sólo función de la posición (x,y,z) p = f1 ( x,y,z) V = f 2 ( x,y,z ) ⇒

⎧u = u ( x, y , z) ⎨v = v ( x , y, z) ⎩w = w( x , y, z )

La trayectoria es el lugar geométrico de las posiciones ocupadas por una misma partícula, cuando varía el tiempo t. Si en un instante dado se asigna a cada punto un vector representando, la velocidad en dicho punto, se obtiene un conjunto de vectores llamado campo de velocidades. La línea de corriente ψ es una línea tangente, en cada uno de sus puntos, a la velocidad en ese punto y en el instante considerado; la línea de corriente satisface la condición: dx = dy = dz u v w En general, las líneas de corriente varian de un instante a otro; en un movimiento permanente son fijas y coinciden con las trayectorias. Se dice que en un punto, línea, superficie o volumen, existe un manantial, cuando en dicho punto, lípfernandezdiez.es

Flujo incompresible no viscoso.IV.-39

nea, superficie o volumen aparecen ciertas cantidades de fluido que a partir del momento en que aparecen participan en la circulación. Un sumidero en un punto, línea, superficie o volumen es aquel en que desaparecen ciertas cantidades de fluido que antes habían participado en la circulación. Un flujo se representa gráficamente mediante las líneas de corriente, que son las envolventes de los vectores velocidad de las partículas fluidas del flujo; cuando el flujo es permanente, las partículas fluidas se mueven a lo largo de trayectorias coincidentes con las líneas de corriente ψ. Si el flujo no es permanente, (régimen transitorio), una configuración de líneas de corriente indica únicamente la representación instantánea del flujo, y en este caso no existe, en general, una correspondencia tan sencilla entre las trayectorias y las líneas de corriente ϕ. El conjunto de las líneas de corriente que pasan por el contorno de un área infinitesimal, en un instante determinado, forman un tubo de fluido que se conoce como tubo de corriente o filete fluido, Fig IV.1, y es de gran utilidad en el estudio de los fenómenos fluidos. De la definición de línea de corriente es evidente que no existe paso de flujo a tra-

Fig IV.1.- Tubo de flujo

vés de la superficie lateral del tubo de corriente; un tubo de corriente se comporta como un conducto de paredes impermeables y espesor nulo, de sección recta infinitesimal.

Un número infinito de tubos de corriente adyacentes, da lugar a un tubo de sección recta finita, que se conoce frecuentemente como vena fluida . El método de estudio puede realizarse a partir del concepto de campo de velocidades V(x, y, z, t), haciendo dos tipos de consideraciones: a) Se pueden fijar las coordenadas (x1, y1, z1) de un punto en las funciones que dan el campo de velocidades, expresándose la velocidad de las partículas móviles al pasar por dicho punto en el transcurso del tiempo; matemáticamente viene expresado por (x1, y1, z 1, t). Mediante esta técnica, conocido un punto fijo del espacio, las velocidades de las diversas partículas que pasan por ese punto, forman un continuo; este punto de vista se conoce como método de Euler. b) Se puede estudiar una partícula genérica del flujo, siguiendo a dicha partícula, método de Lagrange, lo cual significa que (x, y, z) no permanecen constantes en la expresión V(x, y, z, t), sino que varían de forma continua, dando en cada instante la posición de la partícula genérica. Por lo tanto, en este caso, las coordenadas espaciales serán función del tiempo; ambas consideraciones no dependen de si el campo es permanente o no. Para el flujo bidimensional en coordenadas rectangulares: ∂ϕ ∂ϕ ; v= - El potencial de velocidades ϕ se define como: u = ∂x ∂y ∂ψ ∂ψ ; v=- La función corriente ψ se define como: u = ∂y ∂x La derivada total de ϕ se puede poner en la forma: dϕ =

∂ϕ ∂ϕ dx + dy = u dx + v dy ∂x ∂y

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-40

Las líneas equipotenciales son aquellas a lo largo de las cuales la función ϕ es constante, es decir: dy =- u ⇒ dx v

d ϕ = u dx + v dy = 0 ;

dx = - dy v u

que proporciona el gradiente de la línea de potencial. La derivada total de ψ es: dψ =

∂ψ ∂ψ dx + dy = - v dx + u dy ∂x ∂y

Para una línea de corriente, función de corriente constante ψ, resulta: d ψ = - v dx + u dy = 0

;

dy = v dx u



dy = dx v u

que es la ecuación diferencial de las líneas de corriente para el flujo bidimensional. Construcción gráfica de las líneas de corriente y de movimiento.- En la red ortogonal constituida por las líneas, ϕ = Cte, ψ = Cte, vamos a considerar una línea de corriente ψ1 y otra equipotencial ϕ1 que pasan por el punto O de la Fig IV.2, y por otra, la línea de corriente, ψ1 +dψ1, y la equipotencial,   ϕ1 +dϕ1. Si consideramos como origen de coordenadas el punto O y por él trazamos las tangentes s y n , a ψ1 y ϕ1, respectivamente, puesto que la velocidad se proyecta en verdadera magnitud sobre el eje Os por cuanto éste es tangente a la línea de corriente, las relaciones obte€ €  nidas para la velocidad u , que son: ∂ϕ ∂ψ ∂ϕ ∂ψ = , se transforman en : u = = ∂x ∂y ∂s ∂n € y escogiendo para dψ y dϕ valores iguales, los valores de ds y dn tamu=

bién serán iguales, con lo que la malla elemental será un cuadrado. Esto permite la construcción gráfica de la red de líneas de corriente y equipotenciales, a partir de una línea de corriente de la distribución del potencial a lo largo de esta línea de corriente y de las líneas equipotenciales extremas. Este método se conoce como método de Präsil, Fig IV.3. Fig IV.2

Si suponemos una línea de corriente, ψ = ψ1, y sobre la misma dos puntos A1 y A2 lo suficientemente próximos como para que el arco (A1A2) se pueda confundir con la recta, ( A1 A 2 ) = a.

Fig IV.3

Si ϕ1 es el potencial en el punto A1 y, ϕ1 +Δϕ, el potencial en el punto A2, el potencial en el punto B, punto medio del Δϕ segmento ( A1 A 2 ) es: ϕ 1 + 2

La normal al segmento ( A1 A 2 ) en el punto M1 es la tangente a la línea equipotencial que pasa por B.

Δϕ En el punto C sobre la normal, tal que BC = a , el potencial será ϕ 1 + , por lo que de acuerdo con 2 2 la ecuación: pfernandezdiez.es

Flujo incompresible no viscoso.IV.-41



u=

∂ϕ ∂ψ = ∂s ∂n

la línea de corriente que pasa por C tiene una cota de valor ϕ 1 + gráfica de la red, Fig IV.4.

Δϕ ; de aquí se deduce la construcción 2

Fig IV.4.- Trazado de las líneas de corriente

Sobre la línea de corriente ψ1 se toman los puntos A1, A2..., tales que correspondan a potenciales en progresión aritmética de razón Δϕ. Por cada uno de los puntos A1, A2..., se trazan a un mismo lado de la línea de corriente, dos semirrectas que forman 45° con la línea de corriente y que se cortan en los puntos C1, C2, ..., siendo el lugar geométrico de los puntos C la línea de corriente de cota, ϕ 1 +

Δϕ . 2

Si trazamos desde los puntos A1, A2..., las perpendiculares (A1 A1'), (A2 A2'), ..., a la línea de corriente Δϕ en potenciales; a partir de esta nueva línea se construirá otra, ψ 1 se podrá graduar a la línea ϕ 1 + 2 ϕ1 +Δϕ, así como los puntos equipotenciales A1”, A2”..., cubriéndose de esta forma el dominio a estudiar, con líneas de corriente y equipotenciales. El trazado de las líneas de corriente debe respetar la condición de ortogonalidad a las equipotenciales que limitan el campo y la conservación del potencial sobre las equipotenciales. El método es válido si se invierte el papel de las equipotenciales y líneas de corriente. Funciones armónicas, son aquellas que tienen sus derivadas primeras regulares en un dominio y que satisfacen la ecuación de Laplace Δϕ = 0. La suma de dos funciones armónicas es otra función armónica lo que implica que se pueden superponer dos movimientos con potencial de velocidades, sumando sus potenciales. IV.2.- CAUDAL A TRAVÉS DE UNA SUPERFICIE ELEMENTAL En general, el caudal Q de una corriente para una sección determinada, es el volumen de fluido que la atraviesa en la unidad de tiempo, m3/seg. Si se considera un tubo de corriente de sección S 1 normal en cada uno de sus puntos a la línea de corriente correspondiente, y un elemento infinitesimal de sección dS, el volumen de fluido que pasa por dS en el tiempo dt es (V dS dt) ya que (V dt) es la longitud de este tubo de corriente infinitesimal, Fig IV.1. Por la sección S1 pasa un volumen de fluido: dW =

∫S

1

V dS dt = dt

pfernandezdiez.es

∫ S V dS 1

Flujo incompresible no viscoso.IV.-42







y como se ha definido el caudal Q, como el volumen de fluido que pasa por la sección S 1 en la unidad de tiempo, resulta: Q =

dW = dt

∫S

1

V dS

Cuando S1 no sea perpendicular a la línea de corriente en cada punto, el caudal es: Q =

∫S

1

∫S V

Vn dS =

1

cos θ dS

La velocidad media correspondiente a la sección S1 es:  V =

∫S

1

Vn dS S1

=

Q S1

Un elemento de volumen no sigue las líneas de corriente, sino las de movimiento; sólo coincidirán cuando el campo de velocidades sea estacionario. IV.3.- ECUACIÓN DE CONTINUIDAD Si se considera en un fluido en movimiento un paralelepípedo fijo e invariable, Fig IV.5, de volumen, (dx dy dz), en el tiempo t contiene una masa M de fluido de la forma, ρ dx dy dz; en el transcurso de un ∂ρ tiempo dt, el valor de la masa es M + ∂M dt , por lo que el valor de ρ se modifica a ρ + , permanecien∂t ∂t do invariable el volumen (dx dy dz); la variación de masa en el volumen citado en el tiempo dt, es: ∂ρ dt dx dy dz ∂t Este incremento procede de ⎧La masa entrante por las seis caras del paralelepípedo ⎨ ⎩La masa que entra por las fuentes contenidas en el volumen infinitesimal Por la cara (dy dz) entra en el tiempo dt una masa de fluido:

ρ u dt dy dz



Por la cara opuesta sale: {ρ u +

Fig IV.5.- Paralelepípedo elemental de fluido

-

∂( ρ u ) dx } dt dy dz ∂x

En consecuencia, por el conjunto de las dos caras perpendiculares al eje Ox se produce una variación de la masa de fluido:

∂( ρ u) dt dx dy dz ∂x

siendo la variación de las masas de fluido entrantes y salientes por el conjunto de las seis caras: -(

∂( ρ u ) ∂( ρ v ) ∂( ρ w ) + + ) dt dx dy dz ∂x ∂y ∂z A su vez, la variación de masa por unidad de volumen debida a los manantiales y sumideros, en el

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-43

tiempo dt, la representamos por el caudal q, de la forma: n

q=



i=1

⎧ m ⎪ q1i - ∑ q2j , siendo: ⎨ j=1 ⎪ ⎩

n

∑ q1i , la suma de los caudales debidos a los sumideros i=1 m

j

∑ q2 , la suma de los caudales debidos a los manantiales j=1

La variación de masa debida a los manantiales y sumideros es:

ρ q dt dx dy dz por lo que: ∂ρ ∂( ρ u ) ∂( ρ v ) ∂( ρ w ) dt dx dy dz = - ( + + ) dt dx dy dz + ρ q dt dx dy dz ∂t ∂x ∂y ∂z  ∂ρ ∂( ρ u ) ∂( ρ v) ∂( ρ w ) ∂ρ + ( + + ) = ρq ⇒ + div ( ρ V ) = ρ q ∂t ∂x ∂y ∂z ∂t que es la expresión general de la ecuación de continuidad. Para el caso en que no existan manantiales ni sumideros, q = 0



∂ρ ∂( ρ u ) ∂( ρ v ) ∂( ρ w ) +( + + )= 0 ∂t ∂x ∂y ∂z Para fluidos incompresibles, ρ = Cte ; ∂ρ = 0 ∂t

;

∂u ∂v ∂w + + = 0 ∂x ∂y ∂z

 div V = 0

;

CASOS PARTICULARES €

a) Si el fluido se mueve paralelamente al eje Ox, las componentes de la velocidad son: u=u ; v=0 ; w=0 y la ecuación de continuidad: ∂( ρ u ) = 0 ∂x

;

ρ u = Cte

;

ρ g u = C´

;

γ u = C´

Como: Espacio Espacio Peso = Peso = G ⇒ G= γSu Volumen Tiempo Tiempo Volumen S   b) Si el movimiento del fluido es irrotacional, rot V = 0 , y como V = - grad ϕ = - ∇ϕ , por derivar

γu =



de un potencial, se tiene que:  div V = div (- ∇ϕ ) = Δϕ ⇒

 div V + Δϕ =€ 0 ( Ec. de Poisson € )

 ⎧div V = 0 Si ⎨ (ecuación de Laplace), el movimiento es conservativo, permanente e incompresible. ⎩Δϕ = 0



Para un movimiento plano irrotacional, las líneas de corriente tienen como ecuación, ψ (x, y)= Cte, y €

forman una red ortogonal con las líneas equipotenciales dadas por, ϕ (x, y)= Cte:

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-44

Δϕ = 0 ; Δψ = 0 ;



∂ϕ ∂ψ = ∂x ∂y

;

v =

∂ϕ ∂ψ =∂y ∂x

El caudal q que pasa entre dos líneas de corriente ψ1 y ψ2 es: q = Ψ2 -Ψ1 Fig IV.6.- Rotor

La longitud de un arco de línea equipotencial es: d ϕ = V dS

La longitud de un arco de€línea de corriente es: d ψ = V dS  El vector torbellino b es igual a la velocidad angular instantánea:  i ∂ ∂x u

 j   1 1 ∂ w= rot V € = 2 2 ∂y v  i    Como: V = w ∧ R = wx x





u =

 k ∂ ∂z w  j wy y

 k wz z

   = i (wy z - wz y) + j (wz x - wx z ) + k (wx y - wy x)

resulta: €

 rot V =

   i j k ∂ ∂ ∂ ∂x ∂y ∂x (w y z - wz y ) (wz x - wx z ) (w x y - wy x )

    = 2 wx i + 2 w y j + 2 wz k = 2 b

El vector torbellino representa la velocidad angular instantánea:    1 w = b = rot V 2 IV.4.- TORBELLINOS CILÍNDRICOS  Interesa calcular la componente z del vector b ; dicho vector se poner en la forma:  ⎛ i ⎜  1 ⎜ ∂ b = 2 ⎜ ∂x ⎜ u ⎝

 j ∂ ∂y v

 k ⎞ ⎟ ∂ ⎟ 1 ∂w ∂v € ∂u ∂w  ∂v ∂u  = {( )i + ( )j + ( )k} ∂x ⎟ 2 ∂y ∂z ∂z ∂x ∂x ∂y w ⎟⎠ por lo que la componente bz = 1 ( ∂v - ∂u ) 2 ∂x ∂y



Fig IV.7.- Torbellino cilíndrico€ pfernandezdiez.es

y ⎧⎪ u = - V sen θ = - V R De acuerdo con la Fig IV.7: ⎨ ⎪ v = V cos θ = V x ⎩ R ∂V ∂R ∂R R - V ∂u V ∂R ∂y ∂y =- y = 2 ∂y R R R2 = x 2 + y 2 V y ∂V ∂R V ∂R =+ 2 y = ∂R y R R ∂R ∂y ∂y = R ∂y R

=

Flujo incompresible no viscoso.IV.-45





y2 V y2 = - V - 2 ∂V + R R ∂R R3 ∂V ∂R R - V ∂R ∂v = V + x ∂R ∂x ∂x = V + x ∂V ∂R - V x ∂R = ∂x R R R ∂R ∂x R 2 ∂x R2

R 2 = x 2 + y 2 ; ∂R = x = ∂x R 2 2 = - V - x 2 ∂V - V x3 R R ∂R R

2 2 y2 Vy2 1 ( 2 V + x 2 + y 2 ∂V - V x 2 + y 2 ) = bz = 1 ( V + x 2 ∂V - V x3 + V + 2 ∂V ) = 2 R R 2 R ∂R R ∂R R R ∂R R3 R2 R3

= 1 ( 2 V + ∂V - V ) = 1 ( V + ∂V 2 R ∂R R 2 R ∂R que es la expresión que permite calcular la velocidad angular instantánea en cualquier tipo de torbellino cilíndrico. IV.5.- ANÁLISIS DEL MOVIMIENTO DE UN ELEMENTO DE VOLUMEN Durante el movimiento, cada elemento de volumen de fluido experimenta cambios de posición y de orientación, que analizaremos a continuación; para ello vamos a considerar un elemento de volumen cualquiera, en el instante t, que contiene a los puntos M (x, y, z) M´ (x + h, y + k, z + l) muy próximos entre sí.   Sea V (u, v, w) la velocidad en el punto M y V ʹ′(uʹ′, vʹ′, wʹ′) la velocidad en Fig IV.8

M’, por lo que se puede poner:

⎧ ∂u ∂u ∂u + k +l ⎪ u´= u (€x + h, y + k, z + l ) = u + h € ∂x ∂y ∂z ⎪  ⎪ ∂v ∂v ∂v V´ (u´, v´, w´) ⇒ ⎨ v´= v ( x + h, y + k, z + l ) = v + h +k +l ∂x ∂y ∂z ⎪ ⎪ ∂w ∂w ∂w ⎪⎩ w´= w ( x + h, y + k, z + l ) = w + h ∂x + k ∂y + l ∂z en donde se han despreciado infinitésimos de orden superior. Si consideramos la primera de estas ecuaciones, que se puede poner en la forma: u' = u + 1 { l ( ∂u - ∂w ) - k ( ∂v - ∂u )} + h ∂u + 1 { k ( ∂v - ∂u ) + l ( ∂u - ∂w )} 2 ∂z ∂x ∂x ∂y ∂x 2 ∂x ∂y ∂z ∂x y haciendo lo propio con las otras dos y llamando:

ξ = 1 ( ∂w - ∂v ) ; 2 ∂y ∂z g 1 = 1 ( ∂w + ∂v ) ; 2 ∂y ∂z

pfernandezdiez.es

η = 1 ( ∂u - ∂w ) ; 2 ∂z ∂x

ζ = 1 ( ∂v - ∂u ) 2 ∂x ∂y

g 2 = 1 ( ∂u + ∂w ) ; 2 ∂z ∂x

g 3 = 1 ( ∂v + ∂u ) 2 ∂x ∂y

Flujo incompresible no viscoso.IV.-46







⎧ u' = u + ( η l - ξ k ) + h ∂u + k g + l g 3 2 ∂x ⎪ se obtiene: ⎨ v' = v + ( ζ h - ξ l ) + k ∂v + h g 3 + l g1 ∂y ⎪ w' = w + ( ξ k - η h ) + l ∂w + h g 2 + k g1 ⎩ ∂z en las que el vector de proyecciones, ξ, η, ζ, es el vector torbellino del campo de velocidades, que sabemos vale:    1 Γ = b = rot V 2 La velocidad del punto M’ se puede considerar como el resultado de la composición geométrica de tres velocidades: a) Una velocidad en la que sus proyecciones u, v, w, se corresponden con una traslación en bloque de la  partícula, a la velocidad V . ⎧ r1 = η l - ζ k b) Una velocidad en la que sus proyecciones son: ⎨ r2 = ζ h - ξ l , y que se corresponde con una rota€ ⎩ r3 = ξ k - η h  ción en bloque de velocidad angular Γ , siendo estas proyecciones las componentes del producto vectorial:   1 rot V ∧ MMʹ′ = Γ ∧ MM ʹ′ 2 € c) Una velocidad en la que sus proyecciones son: ⎧ d = h ∂u + k g + l g 1 3 2 ∂x ⎪ ∂v ⎨ d 2 = h g 3 + k ∂y + l g1 ⎪ d = h g 2 + k g 1 + l ∂w ⎩ 3 ∂z

  que se corresponden con una deformación D . En consecuencia, V ʹ′ se puede poner en la forma:     V ʹ′ = V + Γ ∧ MM' + D € € demostrándose que ∂u , ∂v , ∂w son las velocidades de deformación lineal. ∂x ∂y ∂z Así, un elemento MM ' de anchura h paralelo a Ox, se transforma en un tiempo dt en otro elemento M1 M 1' de anchura: h ( l + ∂u dt ) ∂x En la misma forma, 2 g1, 2 g2, 2 g3, son velocidades de deformación angular. Así, dos elementos MM ' y MN ' inicialmente paralelos a los ejes Ox y Oy, formando, por lo tanto, entre ellos un ángulo de 90°, en el tiempo dt formarán otro ángulo de ( 90º- ε ) tal que:

ε = ( ∂u + ∂v ) dt ∂y ∂x

 Durante la transformación D se conserva el paralelismo entre las rectas, y un pequeño paralelepí-

pedo rectángulo se transformará en otro paralelepípedo de ángulos diferentes. Estos tres movimientos elementales se esquematizan en la Fig IV.9. € pfernandezdiez.es

Flujo incompresible no viscoso.IV.-47



Fig IV.9.- Movimientos elementales

Si se introduce una función de deformación φ, de la forma:

φ ( h, k, l ) = 1 ( h 2 ∂u + l 2 ∂v + k2 ∂w ) + k l g 1 + l h g 2 + h k g 3 2 ∂x ∂y ∂z se puede comprobar que las componentes, d1, d2 y d3 de la velocidad de deformación, son las derivadas parciales de φ, es decir: d1 =

∂φ ∂h

;

d2 =

∂φ ∂k

;

d3 =

∂φ ∂l



 D = grad φ

La descomposición que se acaba de exponer, debida a Helmholtz, es aquella para la cual la velocidad de deformación deriva de un potencial, - φ. IV.6.- POTENCIAL DE VELOCIDADES PARA FUENTES Y SUMIDEROS EN FLUJO BIDIMENSIONAL En los campos fluidodinámicos pueden existir singularidades de distintos géneros; la más sencilla es la que se presenta cuando se introduce otro fluido en el campo o se extrae del mismo. Esta singularidad toma el nombre de manantial en el primer caso o de sumidero en el segundo. En realidad se trata de una abstracción matemática que podría compararse, toscamente, con una toma o una derivación de una conducción por la que discurriese el fluido. Las singularidades modifican el campo de velocidades de la corriente fluida no perturbada. Si se suman algebraicamente tanto el potencial de la corriente no perturbada como los potenciales de las singularidades, se define el potencial en un punto cualquiera del campo fluidodinámico total. Se demuestra que la combinación de los potenciales de la corriente no perturbada y de las diversas singularidades posibles da lugar a corrientes con características iguales a las que se tendrían si se sumergiesen en el fluido cuerpos de forma determinada. En efecto, si suponemos una fuente colocada en el origen de coordenadas, y un punto A cualquiera de coordenadas (x, y), las líneas de corriente y en este tipo de flujo, Fig IV.10, son rectas radiales. pfernandezdiez.es

Flujo incompresible no viscoso.IV.-48





Fig IV.10.- Potencial de velocidades para una fuente en flujo bidimensional

 La velocidad resultante en el punto A debida al manantial es VR, radial. ⎧u = V cos θ Las componentes (u, v) de la velocidad resultante son: ⎨v = V R sen θ ⎩ R € ⎧ x = R cos θ Las coordenadas (x,y) del punto A son: ⎨ ⎩ y = R sen θ Se define el flujo Q por unidad de tiempo, (caudal), de la forma: Q = 2 π R VR

VR =



Q 2 πR

en la que, la intensidad del manantial se define, para dos y tres dimensiones, en la forma: m( 2

dim ) =

Q ; m( 3 2π

dim ) =

Q 2πb

siendo b la longitud del manantial, perpendicular al plano del dibujo. A su vez, el potencial de velocidades ϕ para la fuente se obtiene a partir de: u=

Q Q ∂ϕ x = x cos θ = = 2 2 2 2πR 2π R ∂x x +y

v=

Q Q y y ∂ϕ sen θ = = 2 = 2 2 2πR 2π R ∂y x +y

ϕ =

Q 2π



x dx 2

x +y

2

=

Q Q Q ln ( x 2 + y 2 ) = ln R 2 = ln R 4π 4π 2π ∂ϕ ∂ψ Q x = = 2 π x2+ y2 ∂x ∂y ∂ϕ ∂ψ Q y == ∂y ∂x 2 π x2+ y2

⎧ ⎪ Para las líneas de corriente se cumple: ⎨ ⎪ ⎩ y la función corriente toma la forma: ∂ψ Q x = 2 ∂y 2 π x + y2



ψ =

Q 2π



x dy 2

x +y

2

=

Q y Q arc tg = θ 2π x 2π

IV.7.- COMBINACIÓN DE UN FLUJO RECTILÍNEO Y UN MANANTIAL El potencial de velocidades para un flujo rectilíneo, como el indicado en la Fig IV.11, viene dado por, pfernandezdiez.es

Flujo incompresible no viscoso.IV.-49

ϕR = u x, y la función corriente por, ψR = u y. El potencial de velocidades es:

ϕ = u0 x +

Q ln ( x 2 + y 2 ) 4π

El potencial de la función corriente es:

ψ = u0 y +

Q y arc tg 2π x

Qx ⎧ u = ∂ϕ = u + 0 ⎪ ∂x 2 π ( x 2+ y2 ) Las componentes de la velocidad son: ⎨ ∂ϕ Qy ⎪ v = ∂y = 2 π ( x 2 + y 2 ) ⎩

Fig IV.11.- Flujo rectilíneo y fuente

Fig IV.12.- Líneas equipotenciales y de movimiento

Las líneas de corriente, ψ = Cte, son de la forma: u0 y +

Q y arc tg = Cte 2π x

y se obtienen dando a la constante los valores: 0 , ± π , ± π , ... 4 2 Se puede escoger como cuerpo sólido el dado por, ψ = 0, semicuerpo de Rankine, que es el indicado en la Fig IV.13, de la forma:

ψ=0 ;

y 2 π u0 y = tg () x Q

También se podría haber tomado como cuerpo sólido cualquiera de los, ψ = Cte, y ser las

Fig IV.13.- Semicuerpo de Rankine

u = u0 +

líneas de corriente los restantes. Un cuerpo semiinfinito separa a la corriente uniforme de la fuente; la parte superior y la parte inferior de dicho semicuerpo, coinciden en un punto de remanso, V= 0, y en él se tiene:

u=0 ; x=-a Qx = Q = u0 - m = 0 =m ; y=0 a 2 π ( x 2+ y2 ) 2π

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-50

x=-a= - m u0

Q ⇒ a= m = u0 2 π u0

⎧⎪ u = u + m cos θ 0 R Las componentes de la velocidad se pueden poner en la forma: ⎨ m sen θ ⎪⎩ v = R Para u = v = 0, se obtiene el punto de remanso A:

θ = 0º ; R = - m u0

 La velocidad V en cualquier punto viene dada por:

V 2 = u 2 + v 2 = u 02 +

2 2 2 2 m u0 2 m u0 cos θ + m2 cos 2θ + m 2 sen 2θ = u02 + cos θ + m 2 = R R R R R

€ = u02 +

2 a u02 a 2 u02 2 cos θ + = u02 (1 + 2 a cos θ + a 2 ) 2 R R R R

Un flujo rectilíneo y una fuente simulan muy bien la parte frontal de un cuerpo cilíndrico inmerso en una corriente fluida. IV.8.- OVALO DE RANKINE Cuando una fuente y un sumidero de igual intensidad se colocan equidistantes del origen de coordenadas, inmersos en una corriente uniforme (u0 x) y todo el fluido de la fuente es absorbido por el sumidero, aparece una línea de corriente divisoria, definida entre el fluido de la corriente uniforme y el fluido transferido de la fuente al sumidero, línea que puede considerarse como la intersección con el plano (x,y) de la superficie de un cilindro de forma ovoidal, conocido como ovalo de Rankine. La superposición de estos flujos da lugar a un flujo externo alrededor de un cilindro ovoidal; combinando muchas fuentes y sumideros se obtiene el flujo aproximado alrededor de un cilindro de forma arbitraria, simétrico respecto al eje Ox. El ovalo de Rankine tiene por ecuaciones, para las líneas equipotenciales y de corriente, de acuerdo con la Fig IV.14, las siguientes:

ϕ = u0x ψ = u0 y -

Q Q Q ( x + a) 2 + y 2 ln r1 + ln r2 = u0 x + ln 2π 2π 4π ( x - a )2+ y2 Q Q Q y y θ + θ = u0 y ( arc tg - arc tg )= 2π 1 2π 2 2π x-a x+a = u0 r sen θ - m ( θ 1 - θ 2 ) = u0 y - m arc tg

2ay x 2 + y2- a2

Fig IV.14.- Ovalo de Rankine pfernandezdiez.es

Flujo incompresible no viscoso.IV.-51

Los semiejes del óvalo, L y h, dependen de la intensidad relativa de la fuente y de la corriente uniforme, es decir, de la relación m ; la línea oval es, ψ = 0. u0 a m , desde 0 a valores elevados, la forma del óvalo aumenta de tau0 a maño y espesor, desde una placa plana de longitud (2 a), hasta un cilindro casi circular. Cuando se aumenta la relación

En el límite, cuando circular.

m → ∞ ⇒ ⎧⎨ L/h → 1 , correspondiente al flujo en torno a un cilindro u0 a ⎩ umáx/u0 → 2

IV.9.- DOBLETE Vamos a suponer una fuente situada en el punto A, Fig IV.15, que consideramos como origen de coordenadas, y un sumidero en B, de igual intensidad, lo cual supone el mismo valor de Q, separados una distancia ds infinitesimal. La fuente y el sumidero pueden estar tan cerca como se quiera, siempre que se mantenga constante el producto de su intensidad por la distancia que los separa (m ds = 0). La función potencial ϕD en cualquier punto P es la suma de las funciones potenciales de la fuente y del sumidero; en consecuencia se puede poner:

ϕD=

Q Q Q Q R + dR dR ln R ln ( R + dR ) = ln =ln ( 1 + ) 2π 2π 2π R 2π R

Desarrollándola en serie de potencias se obtiene:

ϕD= -

Q dR 1 dR 2 Q dR Q dS cos θ ( - ( ) + ...) = =2π R 2 R 2 πR 2πR

en la que: F =

Q dS Q = Cte, y la intensidad de la fuente o del sumidero =m 2π 2π

La función potencial para el doblete es:

ϕ D = - F cos θ = - F 2x = - 2F x 2 R R x +y La función de corriente ψD para el doblete es: u=

∂ϕ D ∂ψ D F ( x 2+ y2 ) - 2 F x 2 F ( x 2+ y2 ) = == ∂x ∂y ( x2 + y2 )2 ( x 2+ y2 )2

Fig IV.15.- Doblete (visto desde muy cerca) pfernandezdiez.es

Flujo incompresible no viscoso.IV.-52





Fig IV.16.- Líneas equipotenciales y de corriente, para un doblete en dos dimensiones (visto desde muy lejos)

v=

∂ϕ D ∂ψ D 2F xy ==∂y ∂x ( x2 + y2 )2

ψD = -



2Fxy 2

2 2

(x + y )

dx =

Fy 2

x + y2

que dice que, las líneas de corriente son círculos con centro situado sobre el eje perpendicular al x, por el punto medio de la distancia entre el manantial y el sumidero. Las líneas equipotenciales son, como hemos visto, círculos. Un doblete visto desde muy lejos, Fig IV.16, supone que las líneas de corriente son círculos tangentes al eje x en el origen, mientras que visto desde muy cerca se corresponde con la representación anterior. Las líneas equipotenciales se pueden poner en la forma: x2 + y2 + F x = 0 ϕD

⎧ a = - F 2 ϕD cuyo centro, de coordenadas (a,b), viene definido por: ⎨ ⎩ b = 0 Las líneas de corriente se pueden poner en la forma: x2 + y2 - F y = 0 ϕD

⎧ a = 0 cuyo centro, de coordenadas (a,b), viene definido por: ⎨ b = F 2 ϕD ⎩ IV.10.- COMBINACIÓN DE UN FLUJO RECTILÍNEO Y UN DOBLETE Supongamos un doblete en el origen O al cual se superpone un flujo uniforme rectilíneo, con velocidad  - u 0 a lo largo del eje Ox, Fig IV.17. El potencial de velocidades para el flujo rectilíneo es, como sabemos, de la forma (- u0 x) El potencial de velocidades total para el flujo combinado es:

ϕ = - u0 x -

Fx x2+ y2

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-53



La función de corriente total para el flujo combinado es:

ψ = - u0 y +

Fy y2

x2+

Las líneas de ψ = Cte, son las líneas de corriente; para el caso particular de, ψ = 0, Fig IV.18, se obtiene el llamado cuerpo de Rankine, de ecuación: Fy = u0 y y2

x2 +

Fig IV.17 Combinación de flujo rectilíneo y doblete

y soluciones

⎧⎪ y = 0 € ⎨ 2 F F 2 2 ⎪⎩ x + y = u0 = a , circunferencia de radio, a = u0 En consecuencia, la combinación de un flujo rectilíneo y un doblete proporciona el flujo en torno a un cilindro circular de radio a; ésto será sólo posible si, como hemos indicado anteriormente, el caudal suministrado por la fuente es igual al recogido por el sumidero. Como: F = a2 u0, el potencial de velocidad total para el flujo combinado y la función de corriente total para el mismo, se pueden poner en la forma:

Fig IV.18.- Cuerpo de Rankine; líneas de corriente y equipotenciales

ϕ = - u0 x ψ = - u0 y +

2 x u0 a2 = - u0 x ( 1 - 2a 2 ) 2 2 x +y x +y 2 y u0 a 2 = u0 y (- 1 + 2a 2 ) 2 2 x +y x +y

Distribución de velocidades en torno a un cilindro circular.- Supongamos un cilindro circular sometido a una corriente fluida uniforme, Fig IV.19; puede ser el caso de una chimenea cilíndrica sometida a la acción del viento, o el de un tubo inmerso en una corriente fluida, etc. x = a cos θ Las coordenadas de un punto cualquiera vienen dadas por: ⎧⎨ ⎩ y = a sen θ La velocidad en dicho punto es: VS =

uS2 + vS2

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-54

Fig IV.19.- Distribución de velocidades en un cilindro circular

A su vez, como: uS = uS = - u0 vS = VS =

∂ϕ ∂ϕ , vS = , se obtiene: ∂x ∂y

u0 a 2 ( y 2 - x 2 ) u a 2 ( a 2 sen 2θ - a 2 cos 2θ ) = - u0 - 0 2 = - 2 u0 sen 2θ 2 2 2 (x + y ) ( a sen 2θ + a 2 cos 2 θ ) 2

2 a 2 x u0 y 2 a 2 u0 a 2 cos θ sen θ ) = = 2 u0 sen θ cos θ ( x 2+ y2 )2 a 4 ( sen 2θ + cos 2θ ) 2 4 u02 sen 4θ + 4 u02 sen 2θ cos 2θ = 2 u0 sen θ

que es la distribución de velocidades en torno a un cilindro circular, función del ángulo θ que define sobre la circunferencia de radio a, la posición del punto genérico A(x,y). IV.11.- VÓRTICE Si se analiza el caso en que se tome la función de corriente de un manantial como función potencial, se obtendrá la formulación correspondiente a un vórtice bidimensional, Fig IV.20, que satisface la ecuación de Laplace Δϕ = 0 y cuyas líneas equipotenciales son rectas radiales, mientras que sus líneas de corriente son círculos concéntricos, de centro el origen de coordenadas.  Las componentes de la velocidad total V son u y v, y para su cálculo se define el concepto de circulación Γ, como la intensidad del vórtice, que es constante, de la forma:

Γ =



∫ V dl

=

∫ u dx

+ v dy =

∫ V€R dθ

= 2πRV



V =

Γ 2πR

Si ϕ es el potencial de velocidades y ψ la función de corriente, se puede poner:



y ∂ϕ ∂ψ Γ sen θ = - Γ = = 2 2 2πR 2π x +y ∂x ∂y ∂ ϕ ∂ ψ x v = V cos θ = - Γ cos θ = Γ = =2πR 2 π x2+ y2 ∂y ∂x u = - V sen θ = -

La integración de estas ecuaciones conduce a:

Fig IV.20.- Vórtice líneas de corriente y equipotenciales pfernandezdiez.es

y ⎧⎪ ϕ = Γ arc tg = 2π x ⎨ ⎪ ψ = - Γ ln R 2 = ⎩ 4π

Γ θ 2π Γ ln ( x 2 + y 2 ) = - Γ ln R 4π 2π Flujo incompresible no viscoso.IV.-55



que son, respectivamente, las ecuaciones de las líneas equipotenciales (rectas radiales) y las de corriente; de ellas se deduce: x2 + y2 = e

-

4 πψ Γ

que son circunferencias concéntricas. IV.12.- COMBINACIÓN DE UNA FUENTE Y UN VÓRTICE En algunos estudios de maquinaria para fluidos, tales como el flujo a través de hélices o rotores de bombas centrífugas, se considera la combinación de una fuente y un vórtice; si se supone a la fuente situada en el origen de coordenadas y se la añade un vórtice, en un punto A(x,y) cualquiera del plano, las   ⎪⎧VF debida al manantial componentes de la velocidad resultante de la combinación V son: ⎨  . ⎪⎩VV debida al vórtice Si se supone que el caudal Q representa la fuerza de la fuente, y la circulación Γ representa la fuerza € del vórtice, se puede poner, Fig IV.21: € Q VF = ; VV = Γ 2 πR 2 πR

Fig IV.21.- Combinación de una fuente y un vórtice

  El ángulo α formado por VF y VV viene dado por: V Q α = arc tg F = arc tg VV Γ € € por lo que para un valor constante de Γ y de Q, el ángulo α será constante para cualquier valor de R, que tg α =

VF VV



es la distancia entre el origen de coordenadas y el punto A(x,y).La función de corriente ψ para la combinación de un vórtice y una fuente es:

ψ=

Q y Q arc tg - Γ ln ( x 2 + y 2 ) = θ - Γ ln R 2π x 4π 2π 2π

ln R =

Qθ-k Γ



R=e

Qθ -k Γ

Para, y = k = 0 R = exp (

Qθ ) Γ

€ que es la ecuación de una espiral logarítmica, Fig IV.21. pfernandezdiez.es

Flujo incompresible no viscoso.IV.-56

El valor de la circulación en torno al centro de un vórtice libre o irrotacional, es constante, e independiente del contorno que se elija, circular o no. En general, la circulación Γ es igual a la suma algebraica de las intensidades de todos los torbellinos que haya en la región interior a la curva cerrada. A continuación veremos que una región de circulación finita en una corriente está sometida a una fuerza de sustentación proporcional a u0 y Γ. IV.13.- FLUJO EN TORNO A UN CILINDRO CON VÓRTICE LIBRE La palabra elevación, o impulso ascendente, significa una fuerza en ángulo recto con la línea de flujo no perturbado. El flujo en torno a un cilindro, combinación de flujo rectilíneo y doblete, es simétrico respecto a la línea de corriente no perturbada, eje Ox, y por lo tanto no hay impulso ascendente sobre el cilindro con este tipo de flujo. Solamente se podrá desarrollar un impulso ascendente dinámico si, el conjunto de las líneas de corriente no es simétrico con respecto a la línea de corriente no perturbada; el flujo no simétrico se puede conseguir sumando un vórtice libre al doblete y al flujo rectilíneo.

Fig IV.22.- Flujo en torno a un cilindro con vórtice libre

 Sea VV la velocidad en la superficie del cilindro debida exclusivamente a un vórtice libre en sentido contrario a las agujas del reloj; la circulación en torno al cilindro es, como sabemos:

Γ€= 2 π a VV siendo a el radio del cilindro. La velocidad resultante es: V = VC + VV = 2 u0 sen θ +

Γ 2πa

Aplicando la ecuación de Bernoulli entre un punto de la corriente no perturbada, y un punto del cilindro, se tiene: p0 + 1 ρ u 02 = p + ρ V 2 2 siendo el valor de la presión:

ρ p = p0 + 1 ρ u02 ( 2 u0 sen θ + Γ ) 2 2 2 2πa

Fig IV.23.- Fuerza ascensional

pfernandezdiez.es

Para calcular la fuerza neta de impulso ascendente, en ángulo recto con el movimiento no perturbado, se hace la integración de las fuerzas de presión sobre el cilindro en la dirección Oy, Fig IV.23. Flujo incompresible no viscoso.IV.-57

Su componente vertical es: dF pvertical = p a dθ sen θ El impulso ascendente, de signo contrario, es: dFascendente = - p a d θ sen θ por lo que el impulso ascendente total sobre el cilindro Fasc es: Fasc = -



∫0

p a sen θ dθ = -

ρ u 02 ) 2

= - a ( p0 +







∫0

{ p0 + 2π

∫0

1 ρ Γ 2 ρ u 02 ( 2 u 0 sen θ + ) } a sen θ dθ = 2 2 2πa

sen θ dθ +

aρ 2



∫0

(4 u 02 sen 3θ +

2 u0 Γ sen 2θ Γ 2 sen θ + ) dθ 2 2 aπ 4π a

y como: 2π€

∫0

sen θ dθ = 0 ;



∫0

sen3θ dθ = 0

resulta: Fasc =

2 aρ 2 u0 Γ θ sen 2θ θ ] 2π = ρ u Γ [ 4 u02 {- cos θ ( sen 2θ + 2 )} + ( ) - Γ cos 0 2 2 aπ 2 4 4 π 2a2 0

que es la fuerza de impulso ascendente por unidad de longitud del cilindro; a esta ecuación se la conoce como de Kutta-Joukowski. IV.14.- EJEMPLOS RELATIVOS A LA RESOLUCIÓN ANALÍTICA DE LA ECUACIÓN DE LAPLACE EN EL PLANO Sea una función de la variable compleja (x + i y) de la forma: f(z) = ϕ (x, y) + i ψ (x, y) A todo numero complejo (x +i y) corresponde en el plano (x,y) un punto bien definido. Se dice que la función f(z) es analítica, cuando la relación: lím Δf ( z ) Δ z→0 Δz es independiente de la forma conque Δz→ 0. Esto implica que: ∂ϕ ∂ψ = ; ∂x ∂y

∂ϕ ∂ψ =∂y ∂x

pudiéndose comprobar que las funciones ϕ y ψ satisfacen la ecuación de Laplace. En consecuencia, la parte real y la parte imaginaria de una función analítica cualquiera de la variable compleja (x + i y) son dos funciones armónicas, de forma que para un flujo en el plano (x,y), la parte real se puede considerar como el potencial ϕ de velocidades y la parte imaginaria como la función ψ de corriente. pfernandezdiez.es

Flujo incompresible no viscoso.IV.-58

Las componentes u y v de la velocidad, vienen dadas, como sabemos, por las relaciones: u=

∂ϕ ∂ψ = ∂x ∂y

;

v =

∂ϕ ∂ψ =∂y ∂x

resultando: df ∂f ∂ϕ ∂ψ ∂f ∂ϕ ∂ψ = = +i = 1 =1( +i ) =u -iv dz ∂x ∂x ∂x i ∂y i ∂y ∂y a) Movimiento uniforme de velocidad u0.- Consideremos la función lineal f ( z ) = u0 z = u0 x + i u 0 y Teniendo en cuenta lo anterior, se deduce que (ϕ = u0 x), (ψ = u0 y), que representan una corriente uniforme en la dirección del eje x; la velocidad se puede calcular a partir de ϕ o de ψ, hallando la derivada de f(z) respecto de z, en la forma: df ∂ϕ ∂ψ ∂ϕ ∂ψ = +i =-i + =u - iv dz ∂x ∂x ∂y ∂y ⎧la parte real de la derivada es igual a la componente u de la velocidad en la que ⎨  ⎩la parte imaginaria es igual a - v Para conseguir un resultado práctico, la derivada

df tiene que existir, y ser única; en este caso se dz

€ df tiene que: = u0 = u , que es real y, por lo tanto, como era de esperar, v = 0. dz

Las líneas equipotenciales, ϕ = Cte, son paralelas al eje Oy, y las líneas de corriente, ψ = Cte, paralelas al Ox, Fig IV.24. b) Flujo radial para una fuente (o un sumidero).- Algunas veces conviene utilizar la variable compleja en polares, en la forma: z=x +iy =

r ei θ =

⎪⎧ r = x 2 + y 2 r cos θ + i r sen θ , siendo: ⎨ y ⎪⎩θ = arc tg x

Fig IV.24.- Movimiento uniforme de velocidad u0

Fig IV.25.- Flujo radial para una fuente

Si se considera: f ( z ) = ϕ + i ψ = m ln z = m ln r e i θ = m (lnr + i θ ) pfernandezdiez.es

Flujo incompresible no viscoso.IV.-59

con la intensidad de la fuente m =

q un número real, se tiene: 2π

⎧ ϕ = m ln r = m ln x 2 + y 2 ⎨ y ψ = m θ = m arc tg ⎩ x La velocidad viene dada por: VR = (

∂ϕ ) = m dr 1 = m ∂r θ = Cte r dr r

correspondiéndose con un manantial o con un sumidero, según que el signo de m sea positivo o negativo, Fig IV.25. c) Flujo en las proximidades de un punto fijo.- Este movimiento viene definido por la función compleja: ⎧ ϕ = a ( x 2 - y 2 ) f ( z ) = a z 2 = a ( x + i y ) 2 = a ( x 2 - y 2 + 2 i x y ) ⇒ ⎨ 2 2 2 2 ⎩ ψ = a x y siendo las líneas de corriente hipérbolas equiláteras de ecuación (a x y = Cte) y las equipotenciales también hipérbolas equiláteras ortogonales a las líneas de corriente.

Fig IV.26.- Flujo en las proximidades de un punto fijo; líneas de corriente y equipotenciales

Las componentes de la velocidad en un punto vienen dadas por: ∂ϕ ⎧ ⎪ u = ∂x = a x ⎨ ∂ϕ ⎪⎩ v = ∂y = - a y siendo, por lo tanto, un movimiento plano alrededor de un punto de estancamiento A. d) Flujo entre dos paredes que forman un ángulo α.- Si una de las paredes se hace coincidir con el eje Ox, la función analítica f(z) se puede poner en la forma: f ( z ) = a z n = a ( x + i y )n , siendo: n = π n n α En polares se tiene que: z = x + i y = r ( cos θ + i sen θ ) ⇒ z n = r n ( cos θ n + i sen θ n )

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-60

⎧⎪ Líneas equipotenciales: ϕ = a r n cos θ n = Cte n ⎨ a ⎪⎩ Líneas de corriente: ψ = r n sen θ n = Cte n En las Fig IV.27-28-29-30-31 se representan algunos casos para diversos valores de a.

Fig IV.27.- a = π/4; f(z)= Az4

Fig IV.28.- a = π/2; f(z) = A z2

Fig IV.30.- a = 3π/2; f(z) = A z2/3

Fig IV.29.- a = π ; f(z) = A z

Fig IV.31.- a = 2π ; f(z) = A z1/2

IV.15.- LINEALIDAD DE LA ECUACIÓN DE LAPLACE La linealidad de la ecuación Δϕ = 0 permite hacer combinaciones lineales de soluciones conocidas de la misma. En efecto, sean f1 y f2 dos funciones, soluciones de la ecuación de Laplace. La función: f1 + f 2 , o cualquier otra combinación lineal de la forma ( λ f1 + µ f 2 ) será también solución de esta ecuación; en tal superposición, los valores de las líneas de corriente, o de las equipotenciales, se ajustan algebraicamente, mientras que las velocidades se componen geométricamente. Por lo tanto, a partir de dos redes de flujo simples, se pueden obtener nuevas redes en las que: Las nuevas líneas equipotenciales vienen dadas por: λ ϕ 1 + µ ϕ 2 Las nuevas líneas de corriente por: λ ψ 1 + µ ψ 2 A título de ejemplo, vamos a tratar el flujo alrededor de un cilindro circular con circulación. El flujo representado por la función que se propone: 2 f ( z ) = u0 ( z + a ) - i Γ ln z z 2π a

tiene una importancia aerodinámica considerable; está constituido por la combinación de dos flujos conocidos: a) El flujo alrededor de un círculo (cilindro) con velocidad uniforme u0 en el infinito b) El flujo alrededor de un vórtice (torbellino) puntual pfernandezdiez.es

Flujo incompresible no viscoso.IV.-61

por lo que representa un flujo con circulación alrededor de un círculo de radio a. La función potencial es: 2 ϕ = Γ θ + u0 ( r + a ) cos θ 2π r

La función de corriente es: 2 ψ = - Γ ln r + u0 ( r - a ) sen θ 2π a r

La línea de corriente correspondiente a, ψ = 0, se compone de un círculo de radio (r = a) y de una curva de ecuación: ln 1 Γ sen θ = Γ ln r = 2 2π a 2 π u0 r u0 ( r - a ) 1r

r ln r a = lím a → 1 2 r→ a 2 a2 a 1- 2 2 r r

=

Γ 4 π u0 r

que junto con (r = a) dan lugar a una serie de casos particulares, según que menor que la unidad, Fig IV.32, 33 y 34.

Γ sea mayor, igual o 4 π u0 r

- Si

Γ < 1, las soluciones son dos puntos reales de estancamiento A y A’ 4 π u0 r

- Si

Γ = 1, los dos puntos de estancamiento se confunden en B 4 π u0 r

- Si

Γ > 1, no hay puntos de estancamiento 4 π u0 r

Fig IV.32.- Solución con dos puntos reales de estancamiento A y A'

Fig IV.33.- Solución con dos puntos de estancamiento, que se confunden en uno solo B

Fig IV.34.- Solución en la que no existen puntos de estancamiento. Existen dos puntos conjugados de velocidad nula pfernandezdiez.es

Flujo incompresible no viscoso.IV.-62

Γ > 1 se encuentran dos puntos de velocidad nula, los puntos C y C’ 4 π u0 r conjugados con relación al círculo, Fig IV.34. Sin embargo, para el caso

Sobre el círculo la velocidad es tangente; el resultado se podía preveer de antemano: ∂ϕ V= (1 ) = Γ - 2 u0 sen θ r ∂θ r=a 2 π a IV.16.- REPRESENTACIÓN CONFORME ⎧ ζ = ξ + i η La relación compleja, ζ = F(z), en la que: ⎨ , expresa una correspondencia tal que, a cada ⎩ z = x + i y par de valores (x, y) corresponde otro par de valores (ξ, η), Fig IV.35. En consecuencia, a todo punto M del plano (x, y) corresponderá un punto P de coordenadas (ξ, η) del plano conjugado; asimismo, a una línea descrita por M corresponderá una línea descrita por P. A la intersección de dos curvas corresponde la intersección de dos curvas transformadas. Se dice que la transformación es conforme cuando conserva los ángulos, es decir, si dos curvas del plano z se cortan en M bajo un cierto ángulo α, en el plano ζ se corresponderá con dos curvas que se cortan en P, homólogo de M, bajo el mismo ángulo, (salvo para ciertos puntos singulares).

Fig IV.35.- Representación conforme

Admitiremos, sin demostración, los siguientes teoremas: - Si ζ = f(z), es una función analítica de la variable z, la transformación que permite pasar de z a ζ es conforme - Si se efectúa una transformación conforme de un flujo plano definido por su red de líneas equipotenciales y de corriente, se tiene que: * Las transformadas de las líneas equipotenciales y de corriente formarán una nueva red de líneas equipotenciales y de corriente en donde cada transformada conserva la magnitud de la línea primitiva * La circulación de las velocidades a lo largo de una línea cualquiera del plano primitivo es igual a la circulación a lo largo de la transformada. Prácticamente la resolución de un problema de flujo se reduce a buscar una transformación conforme que permita una correspondencia entre el flujo desconocido y otro más simple, o a otro ya conocido. pfernandezdiez.es

Flujo incompresible no viscoso.IV.-63

IV.17.- TRANSFORMACIÓN DE JOUKOWSKI Esta transformación es de la forma: 2 ζ = ξ +i η = 1 (z+ b ) 2 z

Si se aplica a las líneas de corriente que derivan del potencial complejo: 2 f ( z ) = u0 ( z + a ) - i Γ ln z z 2π a

se puede encontrar la forma de las líneas de corriente alrededor de obstáculos cuyo perfil sea el de ala de avión, y conocer así el reparto de velocidades alrededor de los mismos. Los resultados teóricos así obtenidos concuerdan perfectamente con la experiencia. Se demuestra que esta transformación equivale a las relaciones siguientes: ⎧⎪ξ = ⎨ ⎪⎩η =

1 (r + 2 1 (r 2

b2 ) cos θ r b2 ) sen θ r

Fig IV.36.- Transformación de Joukowski

En la Fig IV.36 se tiene que: Los puntos A y B se transforman en los puntos A’ y B’ Al círculo C1 de diámetro (AB) le corresponde el segmento (A’B’) A los círculos de ecuación, x2 + y2 = R2, de centro O, corresponden las elipses homofocales de focos A y B Al círculo C2 de centro I, que pasa por los puntos A y B, corresponde el arco de círculo (A’I’B’) Al círculo Q tangente en B al círculo C2, corresponde el perfil Q’ tangente en B’ al arco (A’I’B’) Si se traza el flujo, con o sin circulación, alrededor del círculo, la transformación de Joukowski permitirá obtener el flujo alrededor del perfil de ala de avión, con lo que se entraría en el estudio de perfiles, mas propio de la aerodinámica teórica. Si mediante un procedimiento de transformación conforme pasamos de perfiles circulares a perfiles alares, las presiones teóricas sobre cada punto del perfil se corresponderán con las de los correspondientes puntos del círculo y, por lo tanto, al integrarse las mismas se puede llegar a la sustentación de un ala que se extiende indefinidamente (ala infinita).

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-64

IV.18.- PARADOJA DE D’ALEMBERT Cuando un cuerpo sólido se mueve en el interior de un fluido real, se originan unas fuerzas debidas a la viscosidad. Por el Principio de acción y reacción, el cuerpo ejerce sobre el fluido una fuerza igual y de sentido contrario a la que el fluido ejerce sobre el solido, es decir, el fenómeno de la resistencia que un solido experimenta al moverse en el seno de un fluido, es análogo al de la resistencia que un fluido experimentaría al moverse en el interior de un solido, como una tubería, por ejemplo.

Fig IV.37.- Cilindro y placa en el seno de un fluido perfecto en movimiento

Sea un cilindro, una esfera, etc, Fig IV.37, colocados en el seno de un flujo bidimensional, carente de viscosidad; si el fluido se mueve en régimen laminar, su velocidad en los puntos A y A’ será nula, (puntos de estancamiento), siendo, por lo tanto, su presión máxima; en los puntos B y B' la velocidad alcanza su valor máximo, y la presión su valor mínimo; en consecuencia, la distribución de presiones en torno al cilindro, esfera, etc., es simétrica, llegándose a la conclusión de que estos sólidos introducidos en un fluido ideal, moviéndose en régimen laminar, no están sometidos a ninguna fuerza, permaneciendo quietos, resultado que está en contradicción con la realidad, por cuanto la experiencia dice que este hecho de la distribución simétrica de presiones no se cumple, hecho afortunado atribuible a la viscosidad del fluido, por cuanto los sólidos introducidos en el seno de un fluido en movimiento presentan una resistencia al avance, tanto mayor, cuanto mayor sea la diferencia de presiones existente entre los puntos A y A'.  Si el cilindro, esfera, etc., se mueve con velocidad constante u0 en el seno del fluido ideal en reposo, será equivalente a considerar el cilindro circular fijo, y el fluido ideal moviéndose en sentido contrario con  velocidad - u 0 . € Al considerar el fluido ideal se cumple que: €

- La energía será constante en todos los puntos de una misma línea de corriente. - Es irrotacional, lo cual implica el que la energía es constante en todos los puntos, aunque no estén en la misma línea de corriente. La velocidad del fluido en un punto de la superficie del cilindro es: Vs = 2 u0 sen θ, tal como se muestra en la Fig IV.38. Si la gravedad no interviene, bien porque el plano del dibujo sea horizontal, o porque el fluido sea un gas, aplicando Bernoulli entre la sección O, corriente sin perturbar, y un punto cualquiera S del cilindro, se tiene: Fig IV.38

Circulación de un fluido perfecto en torno a un cilindro pfernandezdiez.es

Flujo incompresible no viscoso.IV.-65

u2 p0 p V2 + 0 = s + s γ 2g γ 2g Despejando ps y teniendo en cuenta que (γ = ρ g), resulta: ps = p0 +

ρ u02 ρ ρ ( u02 - Vs2 ) = p0 + ( u02 - 4 u02 sen 2θ ) = p0 + ( 1 - 4 sen 2θ ) 2 2 2

que se puede poner también en la forma: p s - p0 Δp = = 1 - 4 sen 2θ 2 ρ u0 ρ u02 2 2 Las fuerzas debidas a la presión son normales en cada punto del cilindro; dando valores a θ se obtienen los correspondiente valores de Δp que se manifiestan en la Fig IV.39, en la que el diagrama polar de presiones nos permite visualizar la paradoja de D'Alembert, y en donde por la simetría de las presiones, la resultante según el eje horizontal, fuerza de arrastre, es nula, es decir, un cilindro se mueve en el seno de un fluido ideal sin experimentar resistencia alguna. En el fenómeno, la presión puede ser mayor, igual o menor que cero, obteniéndose los siguientes resultados:

Fig IV.39 Diagrama de presiones en la circulación de un fluido perfecto en torno a un cilindro

⎧⎪ Δp = 0 ; 1 - 4 sen 2θ = 0 ; θ = ± 30º Para: ⎨ Δp > 0 ; 1 - 4 sen 2θ > 0 ; θ < ± 30º ⎪ Δp < 0 ; 1 - 4 sen 2θ < 0 ; θ > ± 30º ⎩

Esto es lo mismo que si el potencial de una singularidad formada por un par de manantiales de signo contrario y muy cercanos se combina con el potencial ϕ = u0 de la corriente unidimensional; se tiene un flujo igual al que se tendría si en la corriente se sumergiese una esfera. En este flujo el movimiento relativo de la esfera estaría privado de resistencia (paradoja de D'Alembert). En cambio, si el potencial del par indicado anteriormente se combina con el de un movimiento horizontal plano (bidimensional), se obtiene un flujo igual al que se tendría si se sumergiese en la corriente un cilindro circular de eje normal al plano y de altura infinitesimal como el espesor de la corriente. Incluso en este flujo el movimiento relativo del cilindro resultaría privado de resistencia porque la distribución de las presiones en la parte anterior sería igual a la distribución de las mismas en la parte posterior. Además no se manifestarían fuerza ni tan siquiera en sentido transversal al movimiento, porque la distribución de las presiones es simétrica incluso en aquella dirección. Si además se añade a esta configuración el potencial de un torbellino se tiene un flujo igual al que se tendría si el cilindro descrito anteriormente empezase a girar sobre su propio eje, Fig IV.40, en la que se observa el comportamiento de los filetes de fluido en un cilindro que gira en sentido contrario a las agujas del reloj alrededor de su eje. Fig IV.40

En estas condiciones la resistencia todavía es nula, pero en sentido transversal la distribución de las presiones no es ya simétrica sobre ambos lados; aparece entonces una resultante normal, perpendicupfernandezdiez.es

Flujo incompresible no viscoso.IV.-66

lar al movimiento, que se denomina sustentación, que tiene un valor que viene dado por:    Vectorialmente: Fasc = ρ ( u0 ∧ Γ ) Escalarmente: Fasc = ρ u0 Γ siendo € Γ la circulación (teorema de Kutta-Joukowski). El cilindro representa el núcleo del torbellino. Puesto que, en circunstancias reales, el fluido no es perfecto, detrás del cilindro se da una distribución de presiones superficiales distinta de la que se da en la parte anterior, porque la vena fluida se separa de la superficie nace de esta forma una resistencia. Conviene entonces achatar la sección circular hasta reducirla a un perfil no simétrico, de manera que en la parte posteFig IV.41 Diagrama de presiones en la circulación de un fluido perfecto en torno a un cilindro en rotación

rior los hilos de fluido se separen del perfil lo más tarde posible.

La asimetría de la distribución de presiones sobre dicho perfil recuerda la del cilindro que gira alrededor de su eje; por eso, el perfil continúa siendo la representación física del torbellino del mismo modo que el cilindro. En la Fig IV.41 se presenta la distribución de presiones alrededor de un cilindro en movimiento de rotación, inmerso en un fluido en movimiento relativo respecto al baricentro del cilindro. IV.19.- FLUJO PARA UN FLUIDO REAL EN TORNO A UN CILINDRO CIRCULAR El flujo alrededor de un cilindro circular es un excelente ejemplo de los efectos causados por el desprendimiento de la capa límite sobre la fuerza de arrastre. Si se considera al cilindro inmerso en la corriente, desde un punto de vista macroscópico la configuración de la misma en torno al cilindro sería la indicada en la Fig IV.40, pero si nos detenemos en un punto cualquiera de la superficie del cilindro, y se pasa a la observación microscópica, se encuentra que la distribución de velocidades, como se sabe, tiene otra forma, según la cual, la capa de fluido contigua al cilindro se adhiere al mismo, por su viscosidad, por lo que su velocidad en el punto de contacto con el cilindro se reduce a cero, Fig IV.42. Esta velocidad aumenta rapidísimamente hasta que, pasada una cierta película de fluido de espesor  δ, la velocidad que adquiere el fluido es V0 correspondiente a las líneas de corriente. En la ecuación de Newton τ = η du se tiene que, aunque η sea muy pequeña, caso del agua y del aidy € re, resulta que el gradiente de velocidades du es muy grande, por cuanto la variación de velocidad tiene dy lugar dentro del espesor δ, muy pequeño, y en consecuencia, el esfuerzo cortante τ es muy grande. En la práctica, y salvo aquellos casos en que la velocidad del fluido sea muy pequeña, no se suele dar la configuración de las lineas de corriente como las hemos visto hasta ahora sino que a partir de un cierto punto del cilindro, las líneas de corriente se separan, desprendiéndose la capa límite, y creándose aguas abajo del cilindro unos remolinos que configuran la estela, y que originan una depresión. pfernandezdiez.es

Flujo incompresible no viscoso.IV.-67

Fig IV.42.- Cilindro inmerso en una corriente y observación microscópica del flujo en sus proximidades

Al estudiar, como caso particular, el flujo bidimensional de un fluido incompresible en torno a un cilindro circular de diámetro D, el coeficiente de arrastre Cw es a su vez función del numero de Reynolds, Fig IV.43, según hemos visto por análisis dimensional, pudiéndose considerar 6 casos para distintos valores de Cw y Re, tal como se indican a continuación.

Fig IV.43.- Valores del coeficiente de arrastre en función del número de Re

El tipo a, Fig IV.44, tiene un numero de Reynolds muy pequeño, Re < 1, por lo que las fuerzas de inercia son pequeñas en comparación con las de viscosidad, cerrándose el flujo suavemente detrás del cilindro, aguas abajo del mismo

Fig IV.44 a) Flujo bidimensional en torno a un cilindro, Re< 1

Fig IV.45 b) Flujo bidimensional en torno a un cilindro, Re = 20

Fig IV.46.- c) Flujo bidimensional en torno a un cilindro, Re > 20 pfernandezdiez.es

Flujo incompresible no viscoso.IV.-68

Fig IV.47.- d) Flujo bidimensional en torno a un cilindro, 5.000 < Re < 15.000

Fig IV.48.- e) Flujo bidimensional en torno a un cilindro, 50.000 < Re < 200.000

Fig IV.49.- f) Flujo bidimensional en torno a un cilindro, Re > 200.000

Al aumentar el numero de Reynolds, el coeficiente de arrastre Cw disminuye; este caso b Fig IV.45, se presenta para, Re ≈ 20. El flujo es laminar a lo largo de la superficie del cilindro aguas arriba; en el punto S se produce la separación de la capa límite laminar. Aguas abajo aparece une región en la que se encuentran dos remolinos estacionarios. Para valores del número de Reynolds, Re > 20, Fig IV.46, los remolinos se hacen inestables, comienzan a vibrar irregularmente, y posteriormente se separan alternativamente de los dos lados, periódicamente, describiendo lo que se conoce como la trayectoria del vórtice de Kàrmàn, existiendo después del cilindro una disposición estable y escalonada de vórtices. La capa límite del lado de la corriente arriba del cilindro es laminar, separándose detrás del cilindro. En la estela posterior aparecen una serie de capas libres laminares. Para valores del numero de Reynolds comprendidos entre, 5.000 < Re < 15.000, Fig IV.47, el flujo se separa del cilindro para formar una estela simétrica con capas libres turbulentas caso d. El tipo e Fig IV.48, aparece para números de Reynolds comprendidos entre, 50.000 < Re < 200.000, siendo, para este caso, constante el coeficiente de arrastre Cw. El punto S es la separación de la capa límite laminar, produciéndose la transición al mismo tiempo que la separación. Para el tipo f Fig IV.49, el numero de Reynolds es, Re > 200.000, y la transición se produce en la capa límite laminar, por delante del punto de separación. En la capa límite turbulenta existirán una intensa mezcla de porciones de fluido, siendo el perfil de velocidades mas brusco, retrasándose, en consecuencia, la separación. Detrás del cilindro, la capa límite turbulenta se separara para formar una estela turbulenta. Como resumen, se puede asegurar que, para el flujo alrededor de un cilindro, y para números de Re muy pequeños, el flujo es laminar en todos sus puntos; para números grandes de Re, el flujo se puede considerar como potencial, salvo en la capa límite y en la estela. La capa límite se forma a partir del punto de estancamiento, aguas arriba del cilindro, y suele ser laminar, para la cual, un gradiente de presión adverso, precipita la separación antes que en una capa lípfernandezdiez.es

Flujo incompresible no viscoso.IV.-69

mite turbulenta, por ser la cantidad de movimiento relativamente pequeña en la capa laminar. Se encuentra que el punto de separación para la capa límite laminar se sitúa mas arriba del cilindro que el punto de separación correspondiente al caso en que la capa límite se hace turbulenta antes de que se separe. La transferencia de la cantidad de movimiento incrementada en la capa límite turbulenta retrasa el punto de separación, deduciéndose la estela en tamaño. IV.20.- PERFIL AERODINÁMICO Para el caso del perfil aerodinámico inmerso en un fluido, Fig IV.50, en el extradós se produce una fuerte depresión (aspiración) y en el intradós se produce una pequeña presión que se extiende hasta el borde de ataque (borde frontal) del ala, para pequeños valores de su ángulo α de inclinación (ángulo de ataque) en relación a la dirección de la corriente no perturbada. El valor medio de la diferencia de presiones entre el intradós y el extradós, para una superficie de ala igual a: S (Lcuerda x 1 metro lineal de ala), proporciona la fuerza de sustentación, de la forma: Fsust =

Cw ρ S u02 2

Fig IV.50.- Coeficientes aerodinámicos de un perfil de ala de avión

en la que Cw es el coeficiente de sustentación, que depende de la forma del perfil (ala) y del ángulo α. Se puede imaginar el ala infinita sustituida por un serie de torbellinos, con el eje paralelo al borde de ataque, como tramos ds de fluido adherente que entran en la longitud L del ala; cuando se examina el ala de longitud L finita, estos torbellinos, que están adheridos a la superficie, no terminan donde el ala está cortada, sino que se repliegan hacia atrás a 90º (torbellinos en forma de herradura o de estribo), formando torbellinos libres. Estos torbellinos, por efecto de la inducción, desvían hacia abajo la corriente detrás del ala, con velo  cidad v menor que la velocidad V del fluido. Por tal motivo, el ala finita puede considerarse equivalente a un ala infinita con el mismo perfil, pero con una incidencia reducida, y puesto que para el ala infinita la sustentación es función lineal de la incidencia, el ala finita tiene una sustentación menor. € € Para un ala infinita, la resultante aerodinámica es perpendicular a la velocidad, (es decir, la resistencia es nula), pero para un ala finita se origina una resistencia inducida Ri, que se añade a la resistencia de perfil Rp debida al rozamiento y a la falta de recuperación de presiones provocada por la separación de la corriente; la suma de estas dos resistencias es: pfernandezdiez.es

Flujo incompresible no viscoso.IV.-70

R=

Cw ρ S V 2 2

En la práctica, la sustentación no crece indefinidamente con α sino que el coeficiente de sustentación alcanza un valor y después disminuye bruscamente para α = 12º÷16º . Esto es debido a que la vena fluida, cuando el ángulo α llega a estos valores, no puede permanecer unida a la superficie y se produce una reducción de la depresión, que es el fenómeno de mayor influencia en la sustentación del ala. Para obtener la resistencia total, es necesario añadir a la resistencia del perfil deducida anteriormente, la resistencia de rozamiento, como la que se tendría sobre una cara de una superficie plana lisa y cuyo coeficiente de arrastre puede obedecer a: C=

0 ,455 ( lg Re ) 2 ,58

;

Re = V L ν

En la capa límite, la corriente es laminar en la primera parte del cuerpo, y después turbulenta debajo de una subcapa laminar (analogía de Prandtl). La resistencia de rozamiento depende de la capa límite, en la que se forman los torbellinos que, al pasar al restante campo del fluido, originan la estela.

pfernandezdiez.es

Flujo incompresible no viscoso.IV.-71

Get in touch

Social

© Copyright 2013 - 2024 MYDOKUMENT.COM - All rights reserved.