Capa límite de Blasius
En física y mecánica de fluidos, una capa límite de Blasius, llamada así por Paul Richard Heinrich Blasius, describe la capa límite laminar bidimensional constante que se forma en una placa semi-infinita que se mantiene paralela a un flujo unidireccional constante. Falkner y Skan generalizaron posteriormente la solución de Blasius al flujo en cuña (capa límite de Falkner-Skan), es decir, a los flujos en los que la placa no es paralela al flujo.
Ecuaciones de la capa límite de Prandtl
Utilizando argumentos de escala, Ludwig Prandtl[1] argumentó que aproximadamente la mitad de los términos de las ecuaciones de Navier-Stokes son despreciables en los flujos de capa límite (excepto en una pequeña región cerca del borde de ataque de la placa). Esto conduce a un conjunto reducido de ecuaciones conocidas como ecuaciones de la capa límite. Para un flujo incompresible constante con viscosidad y densidad constantes, dichas ecuaciones son:
Continuidad de masas:
-Momento:
-Momento:
Aquí el sistema de coordenadas se elige con apuntando paralelo a la placa en la dirección del flujo y la coordenada apuntando normal a la placa, y son las componentes y de la velocidad, es la presión, es la densidad y es la viscosidad cinemática.
Se han encontrado varias soluciones de similitud para este conjunto de ecuaciones para diversos tipos de flujo, incluido el flujo sobre una placa plana delgada. El término similitud se refiere a la propiedad de que los perfiles de velocidad en diferentes posiciones del flujo son los mismos, aparte de los factores de escala. Los factores de escala de similitud reducen el conjunto de ecuaciones diferenciales parciales a un conjunto relativamente fácil de resolver de ecuaciones diferenciales ordinarias no lineales. Paul Richard Heinrich Blasius, uno de los alumnos de Prandtl, desarrolló el modelo de similitud correspondiente al flujo para el caso en que el gradiente de presión, /, a lo largo de una placa plana delgada es despreciable en comparación con cualquier gradiente de presión en la región de la capa límite.[2]
Ecuación de Blasius - capa límite de primer orden
Blasius demostró que para el caso en que , la ecuación de Prandtl -momentum tiene una solución autosimilar. La solución autosimilar existe porque las ecuaciones y las condiciones de contorno son invariantes bajo la transformación
donde es cualquier constante positiva. Introdujo las variables autosimilares
![](http://upload.wikimedia.org/wikipedia/commons/thumb/9/9c/Blasius.png/220px-Blasius.png)
donde
- es el espesor de la capa límite,
- es la velocidad de la corriente libre,
- y es la función de corriente.
La función de corriente es directamente proporcional a la función normalizada, , que es sólo una función de la variable de espesor de similitud. Esto nos lleva directamente a las componentes de velocidad:[3]: 136
Donde el primero denota derivación con respecto a . La sustitución en la ecuación del -momento da la ecuación de Blasius
Las condiciones de contorno son la condición de no deslizamiento, la impermeabilidad de la pared y la velocidad de la corriente libre fuera de la capa límite
Se trata de una ecuación diferencial ordinaria no lineal de tercer orden que puede resolverse numéricamente, por ejemplo, con el método de disparo.
Con la solución para y obtenidas sus derivadas, la ecuación de Prandtl -momentum puede ser adimensionalizada y reordenada para obtener el gradiente de presión , /, como[4]: 46
donde es el espesor de desplazamiento de Blasius.
La velocidad normal de Blasius y el gradiente de presión ascienden a un valor de 0.86 y 0.43, respectivamente, a valores grandes de mientras que asciende a la velocidad de la corriente libre . Como tiende a cero, la escala -gradiente de presión va a 0,16603.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/4/49/Graph_Blasius_velocities_and_y-pressure_gradient_versus_normal_height_above_wall.png/300px-Graph_Blasius_velocities_and_y-pressure_gradient_versus_normal_height_above_wall.png)
La forma límite para pequeña es
y la forma límite para grandes es
- [5]
Los parámetros característicos de las capas límite son el espesor de la capa límite viscosa de dos sigmas,,[6] , el espesor de desplazamiento , el espesor del momento , el esfuerzo cortante de la pared y la fuerza de arrastre actuando sobre una longitud de la placa. Para la solución de Blasius, vienen dadas por
El factor en la fórmula de la fuerza de arrastre es para tener en cuenta ambos lados de la placa.
La integral de momento de Von Kármán y la integral de energía para el perfil de Blasius se reducen a
donde
- es el esfuerzo cortante de pared,
- es la velocidad de inyección/succión de pared,
- es la tasa de disipación de energía,
- es el espesor de momento y
- es el espesor de energía.
Unicidad de la solución de Blasius
La solución de Blasius no es única desde una perspectiva matemática,[7]: 131 como el propio Ludwig Prandtl señaló en su transposición y analizado por una serie de investigadores como Keith Stewartson, Paul A. Libby.[8] A esta solución puede añadirse cualquiera del conjunto discreto infinito de funciones propias, cada una de las cuales satisface la ecuación linealmente perturbada con condiciones homogéneas y decaimiento exponencial en el infinito. La primera de estas funciones resulta ser la derivada de primer orden respecto a de la solución de Blasius, que representa la incertidumbre en la localización efectiva del origen.
Capa límite de segundo orden
Esta aproximación de la capa límite predice una velocidad vertical distinta de cero lejos de la pared, que debe tenerse en cuenta en la capa externa no viscosa del siguiente orden y en la solución correspondiente de la capa límite interna, que a su vez predice una nueva velocidad vertical, y así sucesivamente. La velocidad vertical en el infinito para el problema de la capa límite de primer orden a partir de la ecuación de Blasius es:
La solución para la capa límite de segundo orden es cero. La solución para la capa límite externa no viscosa e interna son[7]: 134
De nuevo como en el problema de frontera de primer orden, a esta solución se le puede añadir cualquiera del conjunto infinito de eigensoluciones. En todas las soluciones puede considerarse como un número de Reynolds.
Capa límite de tercer orden
Dado que el problema interno de segundo orden es cero, las correcciones correspondientes al problema de tercer orden son nulas, es decir, el problema externo de tercer orden es el mismo que el problema externo de segundo orden.[7]: 139 La solución para la corrección de tercer orden no tiene una expresión exacta, pero la expansión de la capa límite interior es de la forma siguiente:
donde es la primera solución propia de la solución de primer orden de la capa límite (que es derivada de la solución de primer orden de Blasius) y la solución para no es única y el problema queda con una constante indeterminada.
Capa límite de Blasius con succión
La succión es uno de los métodos habituales para posponer la separación de la capa límite.[9] Considérese una velocidad de succión uniforme en la pared . Bryan Thwaites[10] demostró que la solución para este problema es la misma que la solución de Blasius sin succión para distancias muy próximas al borde de ataque. Introduciendo la transformación
en las ecuaciones de la capa límite conduce a
con las siguientes condiciones de contorno:
Transformación de Von Mises
Iglisch obtuvo la solución numérica completa en 1944.[11] Si además se introduce la transformación de von Mises[12]
entonces las ecuaciones se convierten en:
con las siguientes condiciones de contorno:
Esta ecuación parabólica en derivadas parciales se puede plantear partiendo de numéricamente.
Perfil asintótico de succión
Dado que la convección debida a la succión y la difusión debida a la pared sólida actúan en sentido contrario, el perfil alcanzará solución estable a gran distancia, a diferencia del perfil de Blasius en el que la capa límite crece indefinidamente. La solución fue obtenida por primera vez por Griffith y F.W. Meredith.[13] Para distancias desde el borde de ataque de la placa , tanto el espesor de la capa límite como la solución son independientes de dada por
Stewartson estudió la correspondencia de la solución completa con el perfil de succión asintótica.[14]
Capa límite de Blasius compresible
Aquí se estudia la capa límite de Blasius con una entalpía específica especificada en la pared es estudiada. La densidad , viscosidad y la conductividad térmica ya no son constantes aquí. La ecuación de conservación de la masa, el momento y la energía serán las siguientes:
donde es el número de Prandtl con sufijo que representa propiedades evaluadas en el infinito. Las condiciones de contorno se convierten en
- ,
- .
A diferencia de la capa límite incompresible, la solución de similitud sólo existe si la transformación
se cumple y esto sólo es posible si .
Transformación de Howarth
![](http://upload.wikimedia.org/wikipedia/commons/thumb/7/72/CompBlasius.jpg/200px-CompBlasius.jpg)
Introduciendo las variables autosimilares mediante la transformación Howarth-Dorodnitsyn
las ecuaciones se reducen a
donde es el coeficiente de dilatación adiabática y es el número de Mach, donde es la velocidad del sonido. La ecuación puede resolverse una vez especificados . Las condiciones de contorno son
Las expresiones comúnmente utilizadas para el aire son . Si es constante, entonces . La temperatura dentro de la capa límite aumentará aunque la temperatura de la placa se mantenga a la misma temperatura que la ambiente, debido al calentamiento disipativo y, por supuesto, estos efectos de disipación sólo son pronunciados cuando el número de Mach es grande.
Capa límite de Blasius de primer orden en coordenadas parabólicas
Dado que las ecuaciones de la capa límite son ecuaciones parabólicas en derivadas parciales, las coordenadas naturales para el problema son coordenadas parabólicas.[7]: 142 La transformación desde coordenadas cartesianas a coordenadas parabólicas viene dada por
- .
Referencias
- ↑ Prandtl, L. (1904). «Über Flüssigkeitsbewegung bei sehr kleiner Reibung». Verhandlinger 3. Int. Math. Kongr. Heidelberg: 484-491.
- ↑ Blasius, H. (1908). «Grenzschichten in Flüssigkeiten mit kleiner Reibung». Z. Angew. Math. Phys. 56: 1-37.
- ↑ Schlichting, H., (1979). Boundary-Layer Theory, 7th ed., McGraw-Hill, New York.
- ↑ Weyburne, David (2022). Aspects of Boundary Layer Theory. ISBN 978-0-578-98334-9.
- ↑ Boyd, J. (2008). «The Blasius function: computations before computers, the value of tricks, undergraduate projects, and open research problems». SIAM Rev. 50: 791-804.
- ↑ Weyburne, D. (2014). «New thickness and shape parameters for the boundary layer velocity profile». Experimental Thermal and Fluid Science 54: 22-28.
- ↑ a b c d Van Dyke, Milton (1975). Perturbation methods in fluid mechanics. Parabolic Press. ISBN 9780915760015.
- ↑ Libby, Paul A., y Herbert Fox. "Algunas soluciones de perturbación en la teoría de la capa límite laminar". Journal of Fluid Mechanics 17.3 (1963): 433-449.
- ↑ Rosenhead, Louis, ed. Laminar boundary layers. Clarendon Press, 1963.
- ↑ Thwaites, Bryan. On certain types of boundary-layer flow with continuous surface suction. HM Stationery Office, 1946.
- ↑ Iglisch, Rudolf. Cálculo exacto de la capa límite laminar en la placa plana de flujo longitudinal con succión homogénea. Oldenbourg, 1944.
- ↑ Von Mises, Richard. "Observaciones sobre hidrodinámica". Z. Angew. Math. Mech 7 (1927): 425-429
- ↑ Griffith, A. A., y F. W. Meredith. "La posible mejora del rendimiento de las aeronaves debido al uso de la succión de la capa límite". Royal Aircraft Establishment Report Nº E 3501 (1936): 12.
- ↑ Stewartson, K. "On asymptotic expansions in the theory of boundary layers". Studies in Applied Mathematics 36.1-4 (1957): 173-191.
Bibliografía
- Parlange, J. Y.; Braddock, R. D.; Sander, G. (1981). «Analytical approximations to the solution of the Blasius equation». Acta Mech. 38 (1–2): 119-125. Bibcode:1981AcMec..38..119P. doi:10.1007/BF01351467.
- Pozrikidis, C. (1998). Introduction to Theoretical and Computational Fluid Dynamics. Oxford. ISBN 978-0-19-509320-9.
- Schlichting, H. (2004). Boundary-Layer Theory. Springer. ISBN 978-3-540-66270-9.
- Wilcox, David C. Basic Fluid Mechanics DCW Industries Inc. 2007
- Boyd, John P. (1999), «The Blasius function in the complex plane», Experimental Mathematics 8 (4): 381-394, ISSN 1058-6458, MR 1737233, doi:10.1080/10586458.1999.10504626 .
Enlaces externos
- [1] - Traducción al inglés del documento original de Blasius - NACA Technical Memorandum 1256.
Datos: Q4925317