stringtranslate.com

Espesor de la capa límite

Esta página describe algunos de los parámetros utilizados para caracterizar el espesor y la forma de las capas límite formadas por un fluido que fluye a lo largo de una superficie sólida. La característica definitoria del flujo en la capa límite es que en las paredes sólidas, la velocidad del fluido se reduce a cero. La capa límite se refiere a la delgada capa de transición entre la pared y el flujo de fluido a granel. El concepto de capa límite fue desarrollado originalmente por Ludwig Prandtl [1] y se clasifica en términos generales en dos tipos, acotado e ilimitado. [2] La propiedad diferenciadora entre capas límite limitadas e ilimitadas es si la capa límite está siendo influenciada sustancialmente por más de una pared. Cada uno de los tipos principales tiene un subtipo laminar , transicional y turbulento . Los dos tipos de capas límite utilizan métodos similares para describir el espesor y la forma de la región de transición con un par de excepciones que se detallan en la Sección Capa límite ilimitada. Las caracterizaciones que se detallan a continuación consideran el flujo estacionario, pero se extienden fácilmente al flujo inestable.

La descripción de la capa límite acotada

Capas límite delimitadas es un nombre utilizado para designar el flujo de fluido a lo largo de una pared interior de manera que las otras paredes interiores inducen un efecto de presión sobre el flujo de fluido a lo largo de la pared bajo consideración. La característica definitoria de este tipo de capa límite es que el perfil de velocidad normal a la pared a menudo tiene una asíntota suave a un valor de velocidad constante denominado u e ( x ). El concepto de capa límite limitada se representa para un flujo constante que ingresa a la mitad inferior de un canal 2-D de placa plana delgada de altura H en la Figura 1 (el flujo y la placa se extienden en la dirección positiva/negativa perpendicular al plano xy ). Ejemplos de este tipo de flujo de capa límite ocurren en el flujo de fluidos a través de la mayoría de tuberías, canales y túneles de viento. El canal 2-D representado en la Figura 1 está estacionario y el fluido fluye a lo largo de la pared interior con una velocidad promedio en el tiempo u ( x , y ), donde x es la dirección del flujo e y es la normal a la pared. La línea discontinua H /2 se agrega para reconocer que se trata de una situación de flujo de canal o tubería interior y que hay una pared superior ubicada sobre la pared inferior que se muestra en la imagen. La Figura 1 muestra el comportamiento del flujo para valores de H que son mayores que el espesor máximo de la capa límite pero menores que el espesor en el que el flujo comienza a comportarse como un flujo exterior. Si la distancia de pared a pared, H , es menor que el espesor de la capa límite viscosa, entonces el perfil de velocidad, definido como u ( x , y ) en x para todo y , adquiere un perfil parabólico en la dirección y y la El espesor de la capa límite es de sólo H /2.

En las paredes sólidas de la placa, el fluido tiene velocidad cero ( condición de frontera sin deslizamiento ), pero a medida que se aleja de la pared, la velocidad del flujo aumenta sin alcanzar un máximo y luego se acerca a una velocidad media constante u e ( x ). . Esta velocidad asintótica puede cambiar o no a lo largo de la pared dependiendo de la geometría de la pared. El punto donde el perfil de velocidad alcanza esencialmente la velocidad asintótica es el espesor de la capa límite. El espesor de la capa límite se representa como la línea discontinua curva que se origina en la entrada del canal en la Figura 1. Es imposible definir una ubicación exacta en la que el perfil de velocidad alcanza la velocidad asintótica. Como resultado, se utilizan varios parámetros de espesor de la capa límite, generalmente indicados como , para describir escalas de espesor características en la región de la capa límite. También es de interés la forma del perfil de velocidad, que es útil para diferenciar los flujos de capa límite laminares de los turbulentos. La forma del perfil se refiere al comportamiento y del perfil de velocidad a medida que pasa a u e ( x ).

Figura 1: Dibujo esquemático que representa el flujo de fluido que ingresa a la mitad inferior de un canal 2D con una separación entre placas de H. El flujo y el canal se extienden perpendicularmente al plano xy .

El 99% del espesor de la capa límite

El espesor de la capa límite, es la distancia normal a la pared hasta un punto donde la velocidad del flujo ha alcanzado esencialmente la velocidad "asintótica" . Antes del desarrollo del Método del Momento, la falta de un método obvio para definir el espesor de la capa límite llevó a gran parte de la comunidad de flujo en la segunda mitad del siglo XX a adoptar la ubicación , denotada y dada por

como el espesor de la capa límite.

Para flujos de capa límite laminar a lo largo de un canal de placa plana que se comportan de acuerdo con las condiciones de la solución de Blasius , el valor se aproxima mucho por [3]

donde es constante y donde

es el número de Reynolds ,
es la velocidad de la corriente libre,
es la velocidad asintótica,
es la distancia aguas abajo desde el inicio de la capa límite, y
es la viscosidad cinemática.

Para capas límite turbulentas a lo largo de un canal de placa plana, el espesor de la capa límite, viene dado por [4]

Esta fórmula del espesor de la capa límite turbulenta supone 1) el flujo es turbulento desde el inicio de la capa límite y 2) la capa límite turbulenta se comporta de una manera geométricamente similar [5] (es decir, los perfiles de velocidad son geométricamente similares junto con el flujo en la dirección x, diferenciándose solo por los parámetros de escala en y ). Ninguno de estos supuestos es cierto para el caso general de la capa límite turbulenta, por lo que se debe tener cuidado al aplicar esta fórmula.

Espesor de desplazamiento

El espesor del desplazamiento, o , es la distancia normal a un plano de referencia que representa el borde inferior de un hipotético fluido no viscoso de velocidad uniforme que tiene el mismo caudal que ocurre en el fluido real con la capa límite. [6]

El espesor de desplazamiento modifica esencialmente la forma de un cuerpo sumergido en un fluido para permitir, en principio, una solución no viscosa si los espesores de desplazamiento fueran conocidos a priori .

La definición del espesor de desplazamiento para flujo compresible , basada en el caudal másico, es

¿Dónde está la densidad? Para flujo incompresible , la densidad es constante, por lo que la definición basada en el caudal volumétrico se convierte en

Para los cálculos de la capa límite turbulenta, se utilizan la densidad y la velocidad promediadas en el tiempo.

Para flujos de capa límite laminar a lo largo de una placa plana que se comportan de acuerdo con las condiciones de la solución de Blasius , el espesor del desplazamiento es [7]

donde es constante.

El espesor del desplazamiento no está directamente relacionado con el espesor de la capa límite, pero se da aproximadamente como . [8] Tiene un papel destacado en el cálculo del factor de forma. También aparece en varias fórmulas del Método del Momento.

Espesor del momento

El espesor del momento, o , es la distancia normal a un plano de referencia que representa el borde inferior de un hipotético fluido no viscoso de velocidad uniforme que tiene el mismo caudal de impulso que ocurre en el fluido real con la capa límite. [9]

La definición de espesor de momento para flujo compresible basada en el caudal másico es [10] [11] [12]

Para flujo incompresible , la densidad es constante, por lo que la definición basada en el caudal volumétrico se convierte en

donde están la densidad y es la velocidad "asintótica".

Para los cálculos de la capa límite turbulenta, se utilizan la densidad y la velocidad promediadas en el tiempo.

Para flujos de capa límite laminar a lo largo de una placa plana que se comportan de acuerdo con las condiciones de la solución de Blasius , el espesor del momento es [13]

donde es constante.

El espesor del momento no está directamente relacionado con el espesor de la capa límite, pero se da aproximadamente como . [14] Tiene un papel destacado en el cálculo del factor de forma.

Un parámetro relacionado llamado Espesor de Energía [15] a veces se menciona en referencia a la distribución de energía turbulenta, pero rara vez se utiliza.

factor de forma

Se utiliza un factor de forma en el flujo de la capa límite para ayudar a diferenciar el flujo laminar y el turbulento. También aparece en varios tratamientos aproximados de la capa límite, incluido el método Thwaites para flujos laminares. La definición formal viene dada por

donde es el factor de forma, es el espesor del desplazamiento y es el espesor del momento.

Convencionalmente, = 2,59 (capa límite de Blasius) es típico de flujos laminares, mientras que = 1,3 - 1,4 es típico de flujos turbulentos cerca de la transición laminar-turbulenta. [16] Para flujos turbulentos próximos a la separación, 2.7. [17] La ​​línea divisoria que define los valores laminar-transicional y transicional-turbulento depende de varios factores, por lo que no siempre es un parámetro definitivo para diferenciar las capas límite laminar, transicional o turbulenta.

Método del momento

Un método relativamente nuevo [18] [19] para describir el espesor y la forma de la capa límite utiliza la metodología matemática del momento que se usa comúnmente para caracterizar funciones de probabilidad estadística . El método del momento de la capa límite se desarrolló a partir de la observación de que la gráfica de la segunda derivada de la capa límite de Blasius para el flujo laminar sobre una placa se parece mucho a una curva de distribución gaussiana. La implicación de la forma gaussiana de la segunda derivada es que la forma del perfil de velocidad para el flujo laminar se aproxima mucho a una función gaussiana dos veces integrada. [20]

El método del momento se basa en integrales simples del perfil de velocidad que utilizan todo el perfil, no solo unos pocos puntos de datos de la región de la cola como lo hace . El método del momento introduce cuatro nuevos parámetros que ayudan a describir el espesor y la forma de la capa límite. Estos cuatro parámetros son la ubicación media, el ancho de la capa límite , la asimetría del perfil de velocidad y el exceso del perfil de velocidad . La asimetría y el exceso son verdaderos parámetros de forma a diferencia de los parámetros de relación simples como el H 12 . La aplicación del método del momento a las derivadas primera y segunda del perfil de velocidad genera parámetros adicionales que, por ejemplo, determinan la ubicación, la forma y el espesor de las fuerzas viscosas en una capa límite turbulenta. Una propiedad única de los parámetros del método del momento es que es posible demostrar que muchos de estos parámetros de espesor de velocidad también son parámetros de escala de similitud. Es decir, si hay similitud en un conjunto de perfiles de velocidad, entonces estos parámetros de espesor también deben ser parámetros de escala de longitud de similitud. [21]

Es sencillo convertir el perfil de velocidad adecuadamente escalado y sus dos primeras derivadas en núcleos integrales adecuados.

Los momentos centrales basados ​​en los perfiles de velocidad escalados se definen como

donde es el espesor del desplazamiento y la ubicación media, está dada por

Hay algunas ventajas al incluir también descripciones de los momentos de las derivadas del perfil de la capa límite con respecto a la altura sobre el muro. Considere los momentos centrales del perfil de velocidad de la primera derivada dados por

donde la ubicación media de la primera derivada es el espesor del desplazamiento .

Finalmente, los momentos centrales del perfil de velocidad de la segunda derivada están dados por

donde la ubicación media de la segunda derivada, está dada por

¿Dónde está la viscosidad y dónde está el esfuerzo cortante de la pared ? La ubicación media, , para este caso se define formalmente como u e ( x ) dividido por el área bajo la curva de la segunda derivada.

Las ecuaciones anteriores funcionan tanto para capas límite laminares como turbulentas siempre que se utilice la velocidad promedio en el tiempo para el caso turbulento.

Con los momentos y las ubicaciones medias definidas, el espesor y la forma de la capa límite se pueden describir en términos de anchos de la capa límite ( varianza ), asimetrías y excesos ( curtosis excesiva ). Experimentalmente, se encuentra que el espesor definido como donde , rastrea el pozo en busca de flujos turbulentos en la capa límite. [22]

Siguiendo el ejemplo de las ecuaciones de equilibrio de momento de la capa límite , los momentos de la segunda derivada de la capa límite, rastrean el espesor y la forma de esa porción de la capa límite donde las fuerzas viscosas son significativas. Por lo tanto, el método del momento hace posible rastrear y cuantificar la capa límite laminar y la región viscosa interna de las capas límite turbulentas usando momentos, mientras que el espesor de la capa límite y la forma de la capa límite turbulenta total se rastrean usando momentos y .

El cálculo de los momentos de la segunda derivada puede ser problemático ya que bajo ciertas condiciones las segundas derivadas pueden volverse positivas en la región más cercana a la pared (en general, es negativa). Este parece ser el caso del flujo interior con un gradiente de presión adverso (APG). Los valores del integrando no cambian de signo en el marco de probabilidad estándar, por lo que la aplicación de la metodología del momento al caso de la segunda derivada dará como resultado medidas de momento sesgadas. Una solución simple [23] es excluir los valores problemáticos y definir un nuevo conjunto de momentos para un perfil de segunda derivada truncado que comienza en el mínimo de la segunda derivada. Si el ancho, , se calcula utilizando el mínimo como ubicación media, entonces el espesor de la capa límite viscosa, definido como el punto donde el perfil de la segunda derivada se vuelve insignificante sobre la pared, se puede identificar adecuadamente con este enfoque modificado.

Para momentos derivados cuyos integrandos no cambian de signo, los momentos se pueden calcular sin la necesidad de tomar derivadas utilizando la integración por partes para reducir los momentos a integrales simples basadas en el núcleo de espesor de desplazamiento dado por

Por ejemplo, el valor de la segunda derivada es y la asimetría de la primera derivada, se puede calcular como

Se demostró que este parámetro rastrea los cambios de forma de la capa límite que acompañan a la transición de la capa límite laminar a turbulenta. [24]

Los errores numéricos encontrados al calcular los momentos, especialmente los de orden superior, son una preocupación seria. Pequeños errores experimentales o numéricos pueden hacer que explote la porción de corriente nominalmente libre de los integrandos. Existen ciertas recomendaciones de cálculo numérico [25] que se pueden seguir para mitigar estos errores.

La descripción de la capa límite ilimitada

Las capas límite ilimitadas, como su nombre lo indica, son típicamente flujos de capas límite exteriores a lo largo de paredes (y algunos flujos interiores con espacios muy grandes en canales y tuberías). Aunque no se aprecia ampliamente, la característica definitoria de este tipo de flujo es que el perfil de velocidad pasa por un pico cerca del borde de la capa límite viscosa y luego lentamente asíntota a la velocidad de la corriente libre u 0 . Un ejemplo de este tipo de flujo de capa límite es el flujo de aire cerca de la pared sobre un ala en vuelo. El concepto de capa límite ilimitada se representa para un flujo laminar estable a lo largo de una placa plana en la Figura 2. La curva discontinua inferior representa la ubicación de la velocidad máxima u max ( x ) y la curva discontinua superior representa la ubicación donde u ( x , y ) esencialmente se convierte en u 0 , es decir . la ubicación del espesor de la capa límite. Para el caso de placa plana muy delgada, el pico es pequeño dando como resultado que la capa límite exterior de la placa plana se parezca mucho al caso del canal plano de flujo interior. Esto ha llevado a que gran parte de la literatura sobre flujo de fluidos trate incorrectamente los casos acotados e ilimitados como equivalentes. El problema con este pensamiento de equivalencia es que el valor máximo máximo puede exceder fácilmente el 10-15% de u 0 para el flujo a lo largo de un ala en vuelo. [26] Las diferencias entre la capa límite limitada e ilimitada se exploraron en una serie de Informes de la Fuerza Aérea. [27] [28] [29]

El pico de la capa límite ilimitada significa que algunos de los parámetros de forma y espesor del perfil de velocidad que se utilizan para los flujos interiores de la capa límite limitada deben revisarse para este caso. Entre otras diferencias, el caso de la capa límite laminar ilimitada incluye regiones viscosas e inerciales dominadas similares a los flujos turbulentos de la capa límite.

Figura 2: Representación de la capa límite laminar “ilimitada” a lo largo de una placa plana 2-D con el flujo y la placa extendiéndose perpendicular al plano xy .

Método del momento

Para flujos exteriores ilimitados de la capa límite, es necesario modificar las ecuaciones de momento para lograr el objetivo deseado de estimar las diversas ubicaciones del espesor de la capa límite. El comportamiento de pico del perfil de velocidad significa que la normalización del área de los momentos se vuelve problemática. Para evitar este problema, se ha sugerido [30] que la capa límite ilimitada se divida en regiones viscosas e inerciales y que el espesor de la capa límite se pueda calcular luego utilizando integrales de momento separadas específicas para esa región. Es decir, la región viscosa interna de las regiones de la capa límite ilimitada laminar y turbulenta se puede rastrear usando momentos modificados, mientras que el espesor de la capa límite inercial se puede rastrear usando momentos y modificados . La lenta velocidad a la que las asíntotas máximas de la velocidad de la corriente libre significa que los valores calculados del espesor de la capa límite son típicamente mucho mayores que en el caso de la capa límite acotada.

Los momentos y modificados para la región de la capa límite inercial se crean: 1) reemplazando el límite integral inferior por la ubicación del pico de velocidad designado por , 2) cambiando el límite integral superior a h donde h está ubicado en lo profundo de la corriente libre, y 3) cambiar la escala de velocidad de a . El espesor del desplazamiento en los momentos modificados debe calcularse utilizando los mismos límites integrales que las integrales de momento modificadas. Al tomar como ubicación media, el espesor de la capa límite 3-sigma modificado se convierte en donde está el ancho modificado .

Los momentos de la segunda derivada modificada se pueden calcular usando las mismas integrales definidas anteriormente pero reemplazando H /2 para el límite integral superior. Para evitar errores numéricos se deben seguir ciertas recomendaciones de cálculo [31] . Las mismas preocupaciones para los momentos de la segunda derivada con respecto a las capas límite acotadas APG para el caso acotado anterior también se aplican a los momentos modificados para el caso no acotado.

En la Figura 3 se muestra un ejemplo de los momentos modificados para el flujo de la capa límite ilimitada a lo largo de una sección de ala. [32] Esta figura se generó a partir de una simulación 2-D [33] para flujo de aire laminar sobre una sección de ala NACA_0012. En esta figura se incluyen el 3-sigma modificado , el 3-sigma modificado y las ubicaciones. El valor de la relación modificada es 311, el valor de la relación modificada es ~2 y el valor es un 9% mayor que el valor. La gran diferencia entre el valor y el valor demuestra lo inadecuado del espesor de la capa límite. Además, el gran pico de velocidad demuestra el problema de tratar las capas límite interiores limitadas como equivalentes a las capas límite exteriores ilimitadas.

Figura 3: El perfil de velocidad de una simulación de perfil aerodinámico NACA0012 en x/c = 0,3. [34]

δmáximoespesor

La ubicación del pico de velocidad, indicada como, es una ubicación de demarcación obvia para la capa límite ilimitada. El principal atractivo de esta elección es que esta ubicación es aproximadamente la ubicación divisoria entre las regiones viscosa e inercial. Para la simulación de flujo laminar a lo largo de un ala, se encuentra que [35] u max ubicado en δ max se aproxima al espesor de la capa límite viscosa dado como +, que indica los picos de velocidad justo por encima del espesor de la capa límite viscosa δ v . Para las regiones inerciales de flujos laminares y turbulentos, es un límite inferior conveniente para las integrales de momento. Si el ancho, , se calcula utilizando como ubicación media, entonces se puede identificar adecuadamente el espesor de la capa límite, definido como el punto donde la velocidad esencialmente se vuelve u 0 por encima de la pared.

El 99% del espesor de la capa límite

Una implicación significativa del comportamiento de pico es que el espesor del 99%, NO se recomienda [36] como parámetro de espesor para la capa límite ilimitada de flujo exterior, ya que ya no corresponde a una ubicación de capa límite de importancia. Solo es útil para flujo laminar ilimitado a lo largo de una placa plana muy delgada con un ángulo de incidencia cero con respecto a la dirección del flujo, ya que el pico para este caso será muy pequeño y el perfil de velocidad se aproximará mucho al caso de la capa límite limitada. Para placas-paredes gruesas, ángulos de incidencia distintos de cero o flujo alrededor de la mayoría de las superficies sólidas, el exceso de flujo debido al arrastre de forma da como resultado un pico cerca de la pared en el perfil de velocidad, lo que no es útil.

Espesor de desplazamiento, espesor de momento y factor de forma.

El espesor del desplazamiento, el espesor del momento y el factor de forma pueden, en principio, calcularse utilizando el mismo enfoque descrito anteriormente para el caso de la capa límite acotada. Sin embargo, la naturaleza puntiaguda de la capa límite ilimitada significa que la sección inercial del espesor de desplazamiento y el espesor de momento tenderán a cancelar la porción cercana a la pared. Por lo tanto, el espesor del desplazamiento y el espesor del momento se comportarán de manera diferente para los casos acotados e ilimitados. Una opción para hacer que el espesor de desplazamiento ilimitado y el espesor de momento se comporten aproximadamente como el caso acotado es usar u max como parámetro de escala y δ max como límite integral superior.

Notas

  1. ^ L.Prandtl, 1904
  2. ^ Weyburne, 2017
  3. ^ Schlichting, p.140
  4. ^ Schlichting, pag. 638
  5. ^ Schlichting, p.152
  6. ^ Schlichting, pag. 140
  7. ^ Schlichting, pag. 141
  8. ^ Schlichting, pag. 28
  9. ^ Schlichting, pag. 141
  10. ^ Schlichting, pag. 354
  11. ^ Whitfield, pag. 13
  12. ^ Schlichting, pag. 258
  13. ^ Schlichting, pag. 141
  14. ^ Schlichting, pag. 161
  15. ^ Schlichting, pag. 354
  16. ^ Schlichting, pag. 454.
  17. ^ X. Wang, W. George, L. Castillo, 2004
  18. ^ Weyburne, 2006
  19. ^ Weyburne, 2014
  20. ^ Weyburne, 2006, pág. 1678
  21. ^ Weyburne, 2017
  22. ^ Weyburne, 2014, pág. 26
  23. ^ Weyburne, 2020a
  24. ^ Weyburne, 2014, pág. 25
  25. ^ Weyburne, 2014
  26. ^ Weyburne, 2020a
  27. ^ Weyburne, 2020a
  28. ^ Weyburne, 2020b
  29. ^ Weyburne, 2020c
  30. ^ Weyburne, 2020a
  31. ^ Weyburne, 2014
  32. ^ Weyburne, 2020a
  33. ^ R. Swanson y S. Langer, 2016
  34. ^ R. Swanson y S. Langer, 2016
  35. ^ Weyburne, 2020a
  36. ^ Weyburne, 2020a

Referencias

Lectura adicional