Articulo de referencia

Capa límite

La capa límite alrededor de una mano humana, fotografía schlieren . La capa límite es el borde verde brillante, más visible en el dorso de la mano (haga clic para ver la imagen ...

La capa límite alrededor de una mano humana, fotografía schlieren . La capa límite es el borde verde brillante, más visible en el dorso de la mano (haga clic para ver la imagen en alta resolución).

En física y mecánica de fluidos , una capa límite es la delgada capa de fluido que se forma en las inmediaciones de una superficie límite debido al flujo del fluido a lo largo de dicha superficie. La interacción del fluido con la pared induce una condición de contorno de no deslizamiento (velocidad cero en la pared). La velocidad del flujo aumenta monótonamente por encima de la superficie hasta alcanzar la velocidad del flujo principal. La delgada capa de fluido cuya velocidad aún no ha regresado a la velocidad del flujo principal se denomina capa límite de velocidad.

El aire que rodea a una persona se calienta, lo que genera un flujo de aire convectivo inducido por la gravedad, que a su vez crea una capa límite tanto de velocidad como térmica. Una brisa perturba esta capa límite, y el cabello y la ropa la protegen, haciendo que la persona sienta más o menos calor. En el ala de un avión , la capa límite de velocidad es la parte del flujo cercana al ala, donde las fuerzas viscosas distorsionan el flujo no viscoso circundante. En la atmósfera terrestre , la capa límite atmosférica es la capa de aire (~1 km) cercana al suelo. Se ve afectada por la superficie; los flujos de calor diurnos y nocturnos causados ​​por el calentamiento del suelo por el sol, la humedad o la transferencia de momento hacia o desde la superficie. 

Tipos de capas límite

Visualización de la capa límite, que muestra la transición de la condición laminar a la turbulenta.

Las capas límite laminares pueden clasificarse, a grandes rasgos, según su estructura y las circunstancias en las que se crean. La delgada capa de cizallamiento que se desarrolla en un cuerpo oscilante es un ejemplo de capa límite de Stokes , mientras que la capa límite de Blasius se refiere a la conocida solución de similitud cerca de una placa plana adherida a un flujo unidireccional incidente, y la capa límite de Falkner-Skan , una generalización del perfil de Blasius. Cuando un fluido rota y las fuerzas viscosas se equilibran con el efecto Coriolis (en lugar de la inercia convectiva), se forma una capa de Ekman . En la teoría de la transferencia de calor, se produce una capa límite térmica. Una superficie puede tener varios tipos de capas límite simultáneamente.

La naturaleza viscosa del flujo de aire reduce las velocidades locales en una superficie y es responsable de la fricción superficial. La capa de aire sobre la superficie del ala que se ralentiza o se detiene debido a la viscosidad es la capa límite. Existen dos tipos diferentes de flujo de capa límite: laminar y turbulento. [ 1 ]

flujo de capa límite laminar

La capa límite laminar es un flujo muy suave, mientras que la capa límite turbulenta contiene remolinos o vórtices. El flujo laminar genera menos resistencia por fricción superficial que el flujo turbulento, pero es menos estable. El flujo de la capa límite sobre la superficie de un ala comienza como un flujo laminar suave. A medida que el flujo se extiende desde el borde de ataque hacia atrás, la capa límite laminar aumenta de espesor.

Flujo de capa límite turbulenta

A cierta distancia del borde de ataque, el flujo laminar se interrumpe y se transforma en un flujo turbulento. Desde el punto de vista de la resistencia aerodinámica, es recomendable que la transición del flujo laminar al turbulento se produzca lo más atrás posible en el ala, o que una gran parte de la superficie del ala se encuentre dentro de la porción laminar de la capa límite. Sin embargo, el flujo laminar de baja energía tiende a interrumpirse de forma más abrupta que la capa turbulenta.

El concepto de capa límite de Prandtl

Ludwig Prandtl
Perfil de velocidad de la capa límite laminar

La capa límite aerodinámica fue propuesta por primera vez por Ludwig Prandtl en un artículo presentado el 12 de agosto de 1904 en el tercer Congreso Internacional de Matemáticos en Heidelberg, Alemania . Esta teoría simplifica las ecuaciones del flujo de fluidos al dividir el campo de flujo en dos áreas: una dentro de la capa límite, dominada por la viscosidad y que genera la mayor parte de la resistencia que experimenta el cuerpo en dicha capa; y otra fuera de la capa límite, donde la viscosidad puede despreciarse sin que ello afecte significativamente a la solución. Esto permite obtener una solución analítica para el flujo en ambas áreas mediante simplificaciones importantes de las ecuaciones completas de Navier-Stokes . La misma hipótesis es aplicable a otros fluidos (además del aire) con viscosidad moderada a baja, como el agua. En el caso de que exista una diferencia de temperatura entre la superficie y el fluido, se observa que la mayor parte de la transferencia de calor hacia y desde un cuerpo tiene lugar en las proximidades de la capa límite de velocidad. Esto, a su vez, permite simplificar las ecuaciones en el campo de flujo fuera de la capa límite. La distribución de presión a lo largo de la capa límite en la dirección normal a la superficie (como en un perfil aerodinámico ) permanece relativamente constante a lo largo de la capa límite, y es la misma que en la superficie misma.

El espesor de la capa límite de velocidad se define normalmente como la distancia desde el cuerpo sólido hasta el punto en el que la velocidad del flujo viscoso es el 99 % de la velocidad de la corriente libre (la velocidad superficial de un flujo no viscoso ). [ 2 ] El espesor de desplazamiento es una definición alternativa que establece que la capa límite representa un déficit en el flujo másico en comparación con el flujo no viscoso con deslizamiento en la pared. Es la distancia por la que la pared tendría que desplazarse en el caso no viscoso para dar el mismo flujo másico total que en el caso viscoso. La condición de no deslizamiento requiere que la velocidad del flujo en la superficie de un objeto sólido sea cero y que la temperatura del fluido sea igual a la temperatura de la superficie. La velocidad del flujo aumentará entonces rápidamente dentro de la capa límite, regida por las ecuaciones de la capa límite, a continuación.

El espesor de la capa límite térmica es, de forma similar, la distancia desde el cuerpo a la cual la temperatura alcanza el 99 % de la temperatura del flujo libre. La relación entre ambos espesores viene determinada por el número de Prandtl . Si el número de Prandtl es igual a 1, ambas capas límite tienen el mismo espesor. Si el número de Prandtl es mayor que 1, la capa límite térmica es más delgada que la capa límite de velocidad. Si el número de Prandtl es menor que 1, como ocurre en el aire en condiciones estándar, la capa límite térmica es más gruesa que la capa límite de velocidad.

En diseños de alto rendimiento, como planeadores y aeronaves comerciales, se presta mucha atención al control del comportamiento de la capa límite para minimizar la resistencia. Hay que considerar dos efectos. Primero, la capa límite aumenta el espesor efectivo del cuerpo, a través del espesor de desplazamiento , lo que incrementa la resistencia de presión. Segundo, las fuerzas de corte en la superficie del ala crean resistencia por fricción superficial .

A números de Reynolds elevados , típicos de aeronaves de tamaño real, es deseable tener una capa límite laminar . Esto resulta en una menor fricción superficial debido al perfil de velocidad característico del flujo laminar. Sin embargo, la capa límite inevitablemente se engrosa y se vuelve menos estable a medida que el flujo se desarrolla a lo largo del cuerpo, y finalmente se vuelve turbulenta , proceso conocido como transición de la capa límite . Una forma de abordar este problema es succionar la capa límite a través de una superficie porosa (véase Succión de la capa límite ). Esto puede reducir la resistencia, pero suele ser poco práctico debido a su complejidad mecánica y la potencia necesaria para mover el aire y eliminarlo. Las técnicas de flujo laminar natural (NLF) desplazan la transición de la capa límite hacia atrás al remodelar el perfil aerodinámico o el fuselaje de manera que su punto más grueso se encuentre más atrás y sea menos grueso. Esto reduce las velocidades en la parte delantera y se logra el mismo número de Reynolds con una mayor longitud.

A números de Reynolds bajos , como los que se observan en los aeromodelos, es relativamente fácil mantener un flujo laminar. Esto produce una baja fricción superficial, lo cual es deseable. Sin embargo, el mismo perfil de velocidad que confiere a la capa límite laminar su baja fricción superficial también provoca que se vea gravemente afectada por gradientes de presión adversos . A medida que la presión comienza a recuperarse en la parte posterior de la cuerda del ala, una capa límite laminar tiende a separarse de la superficie. Dicha separación del flujo provoca un gran aumento de la resistencia de presión , ya que incrementa considerablemente el tamaño efectivo de la sección del ala. En estos casos, puede ser ventajoso inducir deliberadamente la turbulencia en la capa límite en un punto anterior a la ubicación de la separación laminar, utilizando un turbulador . El perfil de velocidad más completo de la capa límite turbulenta le permite soportar el gradiente de presión adverso sin separarse. Por lo tanto, aunque aumenta la fricción superficial, la resistencia total disminuye. Este es el principio que explica el efecto de hoyuelos en las pelotas de golf, así como el funcionamiento de los generadores de vórtices en las aeronaves. También se han diseñado secciones de ala especiales que optimizan la recuperación de presión, reduciendo o incluso eliminando la separación laminar. Esto representa un equilibrio óptimo entre la resistencia de presión derivada de la separación del flujo y la fricción superficial causada por la turbulencia inducida.

Cuando se utilizan modelos a escala reducida en túneles de viento, a veces se emplea una peniche para reducir o eliminar el efecto de la capa límite.

Ecuaciones de la capa límite

La deducción de las ecuaciones de la capa límite fue uno de los avances más importantes en la dinámica de fluidos. Mediante un análisis de orden de magnitud , las conocidas ecuaciones de Navier-Stokes que rigen el flujo de fluidos viscosos pueden simplificarse enormemente dentro de la capa límite. En particular, la característica de las ecuaciones diferenciales parciales (EDP) se vuelve parabólica, en lugar de la forma elíptica de las ecuaciones completas de Navier-Stokes. Esto simplifica enormemente la solución de las ecuaciones. Al realizar la aproximación de la capa límite, el flujo se divide en una porción no viscosa (que es fácil de resolver mediante varios métodos) y la capa límite, que se rige por una EDP más fácil de resolver. Las ecuaciones de continuidad y de Navier-Stokes para un flujo incompresible estacionario bidimensional en coordenadas cartesianas vienen dadas por

incógnita+υy=0{\displaystyle {\partial u \over \partial x}+{\partial \upsilon \over \partial y}=0}
incógnita+υy=1ρpagincógnita+ν(2incógnita2+2y2){\displaystyle u{\partial u \over \partial x}+\upsilon {\partial u \over \partial y}=-{1 \over \rho }{\partial p \over \partial x}+{\nu }\left({\partial ^{2}u \over \partial x^{2}}+{\partial ^{2}u \over \partial y^{2}}\right)}
υincógnita+υυy=1ρpagy+ν(2υincógnita2+2υy2){\displaystyle u{\partial \upsilon \over \partial x}+\upsilon {\partial \upsilon \over \partial y}=-{1 \over \rho }{\partial p \over \partial y}+{\nu }\left({\partial ^{2}\upsilon \over \partial x^{2}}+{\partial ^{2}\upsilon \over \partial y^{2}}\right)}

dónde{\displaystyle u}yυ{\displaystyle \upsilon }son los componentes de la velocidad,ρ{\displaystyle \rho }es la densidad,pag{\displaystyle p}es la presión, yν{\displaystyle \nu }es la viscosidad cinemática del fluido en un punto.

La aproximación establece que, para un número de Reynolds suficientemente alto , el flujo sobre una superficie se puede dividir en una región externa de flujo no viscoso no afectado por la viscosidad (la mayor parte del flujo) y una región cercana a la superficie donde la viscosidad es importante (la capa límite).{\displaystyle u}yυ{\displaystyle \upsilon }sean las velocidades longitudinal y transversal (normal a la pared) respectivamente dentro de la capa límite. Utilizando el análisis de escala , se puede demostrar que las ecuaciones de movimiento anteriores se reducen dentro de la capa límite a convertirse en

incógnita+υy=1ρpagincógnita+ν2y2{\displaystyle u{\partial u \over \partial x}+\upsilon {\partial u \over \partial y}=-{1 \over \rho }{\partial p \over \partial x}+{\nu }{\partial ^{2}u \over \partial y^{2}}}
1ρpagy=0{\displaystyle {1 \over \rho }{\partial p \over \partial y}=0}

y si el fluido es incompresible (como lo son los líquidos en condiciones estándar):

incógnita+υy=0{\displaystyle {\partial u \over \partial x}+{\partial \upsilon \over \partial y}=0}

El análisis de orden de magnitud supone que la escala de longitud longitudinal es significativamente mayor que la escala de longitud transversal dentro de la capa límite. De ello se deduce que las variaciones en las propiedades en la dirección longitudinal son generalmente mucho menores que las de la dirección normal a la pared. Al aplicar esto a la ecuación de continuidad se observa queυ{\displaystyle \upsilon }, la velocidad normal de la pared, es pequeña en comparación con{\displaystyle u}la velocidad longitudinal.

Dado que la presión estáticapag{\displaystyle p}es independiente dey{\displaystyle y}Entonces, la presión en el borde de la capa límite es la presión en toda la capa límite en una posición longitudinal dada. La presión externa se puede obtener mediante la aplicación de la ecuación de Bernoulli . SeaU{\displaystyle U}sea ​​la velocidad del fluido fuera de la capa límite, donde{\displaystyle u}yU{\displaystyle U}son ambos paralelos. Esto da como resultado al sustituir porpag{\displaystyle p}el siguiente resultado

incógnita+υy=UdUdincógnita+ν2y2{\displaystyle u{\partial u \over \partial x}+\upsilon {\partial u \over \partial y}=U{\frac {dU}{dx}}+{\nu }{\partial ^{2}u \over \partial y^{2}}}

Para un flujo en el que la presión estáticapag{\displaystyle p}tampoco cambia en la dirección del flujo

dpagdincógnita=0{\displaystyle {\frac {dp}{dx}}=0}

entoncesU{\displaystyle U}permanece constante.

Por lo tanto, la ecuación de movimiento se simplifica a:

incógnita+υy=ν2y2{\displaystyle u{\partial u \over \partial x}+\upsilon {\partial u \over \partial y}={\nu }{\partial ^{2}u \over \partial y^{2}}}

Estas aproximaciones se utilizan en diversos problemas prácticos de flujo de interés científico y de ingeniería. El análisis anterior se aplica a cualquier capa límite laminar o turbulenta instantánea, pero se utiliza principalmente en estudios de flujo laminar, ya que el flujo medio también es el flujo instantáneo debido a la ausencia de fluctuaciones de velocidad. Esta ecuación simplificada es una EDP parabólica y puede resolverse mediante una solución de similitud, a menudo denominada capa límite de Blasius .

Teorema de transposición de Prandtl

Prandtl observó que de cualquier solución(incógnita,y,t), v(incógnita,y,t){\displaystyle u(x,y,t),\ v(x,y,t)}que satisface las ecuaciones de la capa límite, solución adicional(incógnita,y,t), v(incógnita,y,t){\displaystyle u^{*}(x,y,t),\ v^{*}(x,y,t)}, que también satisface las ecuaciones de la capa límite, se puede construir escribiendo [ 3 ]

(incógnita,y,t)=(incógnita,y+F(incógnita),t),v(incógnita,y,t)=v(incógnita,y+F(incógnita),t)F(incógnita)(incógnita,y+F(incógnita),t){\displaystyle u^{*}(x,y,t)=u(x,y+f(x),t),\quad v^{*}(x,y,t)=v(x,y+f(x),t)-f'(x)u(x,y+f(x),t)}

dóndeF(incógnita){\displaystyle f(x)}es arbitrario. Dado que la solución no es única desde una perspectiva matemática, [ 4 ] a la solución se le puede agregar cualquiera de un conjunto infinito de autofunciones como lo muestran Stewartson [ 5 ] y Paul A. Libby . [ 6 ] [ 7 ]

Integral de momento de von Kármán

Von Kármán derivó la ecuación integral integrando la ecuación de la capa límite a través de la capa límite en 1921. [ 8 ] La ecuación es

τwρU2=1U2t(Uδ1)+δ2incógnita+2δ2+δ1UUincógnita+vwU{\displaystyle {\frac {\tau _{w}}{\rho U^{2}}}={\frac {1}{U^{2}}}{\frac {\partial }{\partial t}}(U\delta _{1})+{\frac {\partial \delta _{2}}{\partial x}}+{\frac {2\delta _{2}+\delta _{1}}{U}}{\frac {\partial U}{\partial x}}+{\frac {v_{w}}{U}}}

dónde

τw=μ(y)y=0,vw=v(incógnita,0,t),δ1=0(1U)dy,δ2=0U(1U)dy{\displaystyle \tau _{w}=\mu \left({\frac {\partial u}{\partial y}}\right)_{y=0},\quad v_{w}=v(x,0,t),\quad \delta _{1}=\int _{0}^{\infty }\left(1-{\frac {u}{U}}\right)\,dy,\quad \delta _{2}=\int _{0}^{\infty }{\frac {u}{U}}\left(1-{\frac {u}{U}}\right)\,dy}
τw{\displaystyle \tau _{w}}es la tensión cortante de la pared,vw{\displaystyle v_{w}}es la velocidad de succión/inyección en la pared,δ1{\displaystyle \delta _{1}}es el espesor de desplazamiento yδ2{\displaystyle \delta _{2}}es el espesor de momento. La aproximación de Kármán-Pohlhausen se deriva de esta ecuación.

Integral de energía

La integral de energía fue derivada por Wieghardt . [ 9 ] [ 10 ]

2ερU3=1Ut(δ1+δ2)+2δ2U2Ut+1U3incógnita(U3δ3)+vwU{\displaystyle {\frac {2\varepsilon }{\rho U^{3}}}={\frac {1}{U}}{\frac {\partial }{\partial t}}(\delta _{1}+\delta _{2})+{\frac {2\delta _{2}}{U^{2}}}{\frac {\partial U}{\partial t}}+{\frac {1}{U^{3}}}{\frac {\partial }{\partial x}}(U^{3}\delta _{3})+{\frac {v_{w}}{U}}}

dónde

ε=0μ(y)2dy,δ3=0U(12U2)dy{\displaystyle \varepsilon =\int _{0}^{\infty }\mu \left({\frac {\partial u}{\partial y}}\right)^{2}dy,\quad \delta _{3}=\int _{0}^{\infty }{\frac {u}{U}}\left(1-{\frac {u^{2}}{U^{2}}}\right)\,dy}
ε{\displaystyle \varepsilon }es la tasa de disipación de energía debida a la viscosidad a través de la capa límite yδ3{\displaystyle \delta _{3}}es el espesor de energía. [ 11 ]

Transformación de Von Mises

Para capas límite bidimensionales estacionarias, von Mises [ 12 ] introdujo una transformación que tomaincógnita{\displaystyle x}yψ{\displaystyle \psi }( función de flujo ) como variables independientes en lugar deincógnita{\displaystyle x}yy{\displaystyle y}y utiliza una variable dependienteχ=U22{\displaystyle \chi =U^{2}-u^{2}}en lugar de{\displaystyle u}La ecuación de la capa límite se convierte entonces en:

χincógnita=νU2χ2χψ2{\displaystyle {\frac {\partial \chi }{\partial x}}=\nu {\sqrt {U^{2}-\chi }}\,{\frac {\partial ^{2}\chi }{\partial \psi ^{2}}}}

Las variables originales se recuperan de

y=U2χdψ,=U2χ,v=incógnita(1)dψ.{\displaystyle y=\int {\sqrt {U^{2}-\chi }}\,d\psi ,\quad u={\sqrt {U^{2}-\chi }},\quad v=u\int {\frac {\partial }{\partial x}}\left({\frac {1}{u}}\right)\,d\psi .}

Esta transformación es extendida posteriormente a la capa límite compresible por von Kármán y HS Tsien . [ 13 ]

La transformación de Crocco

Para una capa límite compresible bidimensional estacionaria, Luigi Crocco [ 14 ] introdujo una transformación que tomaincógnita{\displaystyle x}y{\displaystyle u}como variables independientes en lugar deincógnita{\displaystyle x}yy{\displaystyle y}y utiliza una variable dependienteτ=μ/y{\displaystyle \tau =\mu \partial u/\partial y}(esfuerzo cortante) en lugar de{\displaystyle u}La ecuación de la capa límite se convierte entonces en:

μρincógnita(1τ)+2τ2μdpagdincógnita(1τ)=0,si dpagdincógnita=0, entonces μρτ2τincógnita=12τ2.{\displaystyle {\begin{aligned}&\mu \rho u{\frac {\partial }{\partial x}}\left({\frac {1}{\tau }}\right)+{\frac {\partial ^{2}\tau }{\partial u^{2}}}-\mu {\frac {dp}{dx}}{\frac {\partial }{\partial u}}\left({\frac {1}{\tau }}\right)=0,\\[5pt]&{\text{if }}{\frac {dp}{dx}}=0,{\text{ then }}{\frac {\mu \rho }{\tau ^{2}}}{\frac {\partial \tau }{\partial x}}={\frac {1}{u}}{\frac {\partial ^{2}\tau }{\partial u^{2}}}.\end{aligned}}}

La coordenada original se recupera de

y=μdτ.{\displaystyle y=\mu \int {\frac {du}{\tau }}.}

Capas límite turbulentas

El tratamiento de las capas límite turbulentas es mucho más difícil debido a la variación temporal de las propiedades del flujo. Una de las técnicas más utilizadas para abordar los flujos turbulentos es la descomposición de Reynolds . En este método, las propiedades instantáneas del flujo se descomponen en un componente medio y uno fluctuante, suponiendo que el valor medio del componente fluctuante es siempre cero. La aplicación de esta técnica a las ecuaciones de la capa límite proporciona las ecuaciones completas de la capa límite turbulenta, que no suelen aparecer en la literatura.

¯incógnita+v¯y=0{\displaystyle {\partial {\overline {u}} \over \partial x}+{\partial {\overline {v}} \over \partial y}=0}
¯¯incógnita+v¯¯y=1ρpag¯incógnita+ν(2¯incógnita2+2¯y2)y(v¯)incógnita(2¯){\displaystyle {\overline {u}}{\partial {\overline {u}} \over \partial x}+{\overline {v}}{\partial {\overline {u}} \over \partial y}=-{1 \over \rho }{\partial {\overline {p}} \over \partial x}+\nu \left({\partial ^{2}{\overline {u}} \over \partial x^{2}}+{\partial ^{2}{\overline {u}} \over \partial y^{2}}\right)-{\frac {\partial }{\partial y}}({\overline {u'v'}})-{\frac {\partial }{\partial x}}({\overline {u'^{2}}})}
¯v¯incógnita+v¯v¯y=1ρpag¯y+ν(2v¯incógnita2+2v¯y2)incógnita(v¯)y(v2¯){\displaystyle {\overline {u}}{\partial {\overline {v}} \over \partial x}+{\overline {v}}{\partial {\overline {v}} \over \partial y}=-{1 \over \rho }{\partial {\overline {p}} \over \partial y}+\nu \left({\partial ^{2}{\overline {v}} \over \partial x^{2}}+{\partial ^{2}{\overline {v}} \over \partial y^{2}}\right)-{\frac {\partial }{\partial x}}({\overline {u'v'}})-{\frac {\partial }{\partial y}}({\overline {v'^{2}}})}

Utilizando un análisis de orden de magnitud similar, las ecuaciones anteriores se pueden reducir a términos de orden principal. Al elegir escalas de longitudδ{\displaystyle \delta }para cambios en la dirección transversal, yL{\displaystyle L}para cambios en la dirección del flujo, conδ<<L{\displaystyle \delta <<L}, la ecuación de momento en x se simplifica a:

¯¯incógnita+v¯¯y=1ρpag¯incógnita+ν2¯y2y(v¯).{\displaystyle {\overline {u}}{\partial {\overline {u}} \over \partial x}+{\overline {v}}{\partial {\overline {u}} \over \partial y}=-{1 \over \rho }{\partial {\overline {p}} \over \partial x}+{\nu }{\partial ^{2}{\overline {u}} \over \partial y^{2}}-{\frac {\partial }{\partial y}}({\overline {u'v'}}).}

Esta ecuación no satisface la condición de no deslizamiento en la pared. Al igual que Prandtl hizo con sus ecuaciones de capa límite, se debe utilizar una nueva escala de longitud más pequeña para permitir que el término viscoso se convierta en el término de orden principal en la ecuación de momento. Al elegirη<<δ{\displaystyle \eta <<\delta }En la escala y , la ecuación de momento de primer orden para esta "capa límite interna" viene dada por:

0=1ρpag¯incógnita+ν2¯y2y(v¯).{\displaystyle 0=-{1 \over \rho }{\partial {\overline {p}} \over \partial x}+{\nu }{\partial ^{2}{\overline {u}} \over \partial y^{2}}-{\frac {\partial }{\partial y}}({\overline {u'v'}}).}

En el límite de número de Reynolds infinito, se puede demostrar que el término de gradiente de presión no tiene efecto en la región interna de la capa límite turbulenta. La nueva "escala de longitud interna"η{\displaystyle \eta }es una escala de longitud viscosa y es de ordenν{\displaystyle {\frac {\nu }{u_{*}}}}, con{\displaystyle u_{*}}siendo la escala de velocidad de las fluctuaciones turbulentas, en este caso una velocidad de fricción .

A diferencia de las ecuaciones de la capa límite laminar, la presencia de dos regímenes regidos por diferentes conjuntos de escalas de flujo (es decir, la escala interna y externa) ha hecho que encontrar una solución de similitud universal para la capa límite turbulenta sea difícil y controvertido. Para encontrar una solución de similitud que abarque ambas regiones del flujo, es necesario igualar asintóticamente las soluciones de ambas regiones. Dicho análisis dará como resultado la llamada ley logarítmica o ley potencial .

Enfoques similares al análisis anterior también se han aplicado a capas límite térmicas, utilizando la ecuación de energía en flujos compresibles. [ 15 ] [ 16 ]

El término adicionalv¯{\displaystyle {\overline {u'v'}}}En las ecuaciones de la capa límite turbulenta, la tensión cortante de Reynolds es desconocida a priori . Por lo tanto, la solución de dichas ecuaciones requiere el uso de un modelo de turbulencia , cuyo objetivo es expresar la tensión cortante de Reynolds en función de variables de flujo conocidas o sus derivadas. La falta de precisión y generalidad de estos modelos constituye un obstáculo importante para la predicción exitosa de las propiedades del flujo turbulento en la dinámica de fluidos moderna.

En la región cercana a la pared existe una capa de tensión constante. Debido a la amortiguación de las fluctuaciones de velocidad vertical cerca de la pared, el término de tensión de Reynolds se vuelve despreciable y se observa un perfil de velocidad lineal. Esto solo se cumple en la región muy cercana a la pared .

Transferencia de calor y masa

En 1928, el ingeniero francés André Lévêque observó que la transferencia de calor por convección en un fluido en movimiento se ve afectada únicamente por los valores de velocidad muy cercanos a la superficie. [ 17 ] [ 18 ] Para flujos con un número de Prandtl elevado, la transición de temperatura/masa desde la superficie hasta la temperatura de la corriente libre tiene lugar a través de una región muy delgada cercana a la superficie. Por lo tanto, las velocidades del fluido más importantes son aquellas dentro de esta región muy delgada en la que el cambio de velocidad puede considerarse lineal con la distancia normal desde la superficie. De esta manera, para

(y)=U[1(yh)2h2]=Uyh[2yh],{\displaystyle u(y)=U\left[1-{\frac {(y-h)^{2}}{h^{2}}}\right]=U{\frac {y}{h}}\left[2-{\frac {y}{h}}\right]\;,}

cuandoy0{\displaystyle y\rightarrow 0}, entonces

(y)2Uyh=θy,{\displaystyle u(y)\approx 2U{\frac {y}{h}}=\theta y,}

donde θ es la tangente de la parábola de Poiseuille que interseca la pared. Aunque la solución de Lévêque era específica para la transferencia de calor en un flujo de Poiseuille, su intuición ayudó a otros científicos a encontrar una solución exacta del problema de la capa límite térmica. [ 19 ] Schuh observó que en una capa límite, u es nuevamente una función lineal de y , pero que en este caso, la tangente de la pared es una función de x . [ 20 ] Expresó esto con una versión modificada del perfil de Lévêque,

(y)=θ(incógnita)y.{\displaystyle u(y)=\theta (x)y.}

Esto da como resultado una muy buena aproximación, incluso para valores bajos.PAGr{\displaystyle Pr}números, de modo que solo los metales líquidos conPAGr{\displaystyle Pr}mucho menor que 1 no puede tratarse de esta manera. [ 19 ] En 1962, Kestin y Persen publicaron un artículo que describía soluciones para la transferencia de calor cuando la capa límite térmica está contenida completamente dentro de la capa de momento y para varias distribuciones de temperatura de la pared. [ 21 ] Para el problema de una placa plana con un salto de temperatura enincógnita=incógnita0{\displaystyle x=x_{0}}Proponen una sustitución que reduce la ecuación parabólica de la capa límite térmica a una ecuación diferencial ordinaria . La solución de esta ecuación, la temperatura en cualquier punto del fluido, puede expresarse como una función gamma incompleta . [ 18 ] Schlichting propuso una sustitución equivalente que reduce la ecuación de la capa límite térmica a una ecuación diferencial ordinaria cuya solución es la misma función gamma incompleta. [ 22 ] Se pueden derivar soluciones analíticas con el Ansatz autosimilar dependiente del tiempo para las ecuaciones de la capa límite incompresible que incluyen la conducción de calor. [ 23 ]

Como es bien sabido en varios libros de texto, la transferencia de calor tiende a disminuir con el aumento de la capa límite. Recientemente, se observó a gran escala que el viento que fluye a través de un generador fotovoltaico tiende a "atrapar" el calor en los paneles fotovoltaicos bajo un régimen turbulento debido a la disminución de la transferencia de calor. [ 24 ] A pesar de que frecuentemente se asume que es inherentemente turbulento, esta observación accidental demuestra que el viento natural se comporta en la práctica muy cerca de un fluido ideal, al menos en una observación que se asemeja al comportamiento esperado en una placa plana, lo que potencialmente reduce la dificultad de analizar este tipo de fenómeno a mayor escala.

Constantes de transferencia convectiva a partir del análisis de la capa límite

Paul Richard Heinrich Blasius derivó una solución exacta a las ecuaciones de la capa límite laminar anteriores . [ 25 ] El espesor de la capa límiteδ{\displaystyle \delta }es una función del número de Reynolds para flujo laminar.

δ5.0incógnitaRmi{\displaystyle \delta \approx 5.0{x \over {\sqrt {Re}}}}
δ{\displaystyle \delta }= el espesor de la capa límite: la región del flujo donde la velocidad es menor al 99% de la velocidad en el campo lejano.v{\displaystyle v_{\infty }};incógnita{\displaystyle x}es posición a lo largo de la placa semiinfinita, yRmi{\displaystyle Re}es el número de Reynolds dado porρvincógnita/μ{\displaystyle \rho v_{\infty }x/\mu }(ρ={\displaystyle \rho =}densidad yμ={\displaystyle \mu =}viscosidad dinámica).

La solución de Blasius utiliza condiciones de contorno en forma adimensional:

vincógnitavSvvS=vincógnitav=vyv=0{\displaystyle {v_{x}-v_{S} \over v_{\infty }-v_{S}}={v_{x} \over v_{\infty }}={v_{y} \over v_{\infty }}=0}     en     y=0{\displaystyle y=0}
vincógnitavSvvS=vincógnitav=1{\displaystyle {v_{x}-v_{S} \over v_{\infty }-v_{S}}={v_{x} \over v_{\infty }}=1}     en     y={\displaystyle y=\infty }yincógnita=0{\displaystyle x=0}
Las capas límite de velocidad y temperatura comparten una forma funcional.
La capa límite de velocidad (arriba, naranja) y la capa límite de temperatura (abajo, verde) comparten una forma funcional debido a la similitud en los balances de momento/energía y las condiciones de contorno.

Tenga en cuenta que en muchos casos se cumple la condición de contorno de no deslizamiento.vS{\displaystyle v_{S}}La velocidad del fluido en la superficie de la placa es igual a la velocidad de la placa en todos los puntos. Si la placa no se mueve, entoncesvS=0{\displaystyle v_{S}=0}Se requiere una derivación mucho más compleja si se permite el deslizamiento del fluido. [ 26 ]

De hecho, la solución de Blasius para el perfil de velocidad laminar en la capa límite sobre una placa semiinfinita puede extenderse fácilmente para describir las capas límite térmica y de concentración para la transferencia de calor y masa, respectivamente. En lugar del balance diferencial de momento en el eje x (ecuación de movimiento), se utiliza un balance de energía y masa derivado de forma similar:

Energía:        vincógnitaTincógnita+vyTy=kρdopag2Ty2{\displaystyle v_{x}{\partial T \over \partial x}+v_{y}{\partial T \over \partial y}={k \over \rho C_{p}}{\partial ^{2}T \over \partial y^{2}}}

Masa:          vincógnitadoAincógnita+vydoAy=DAB2doAy2{\displaystyle v_{x}{\partial c_{A} \over \partial x}+v_{y}{\partial c_{A} \over \partial y}=D_{AB}{\partial ^{2}c_{A} \over \partial y^{2}}}

Para el balance de momento, viscosidad cinemáticaν{\displaystyle \nu }puede considerarse como la difusividad de momento . En el balance de energía, esto se reemplaza por la difusividad térmica.α=k/ρdoPAG{\displaystyle \alpha ={k/\rho C_{P}}}y por difusividad de masaDAB{\displaystyle D_{AB}}en el balance de masas. En la difusividad térmica de una sustancia,k{\displaystyle k}es su conductividad térmica,ρ{\displaystyle \rho }es su densidad ydoPAG{\displaystyle C_{P}}es su capacidad calorífica. El subíndice AB denota la difusividad de la especie A que se difunde en la especie B.

Bajo el supuesto de queα=DAB=ν{\displaystyle \alpha =D_{AB}=\nu }, estas ecuaciones se vuelven equivalentes al balance de momento. Por lo tanto, para el número de PrandtlPAGr=ν/α=1{\displaystyle Pr=\nu /\alpha =1}y número de SchmidtSdo=ν/DAB=1{\displaystyle Sc=\nu /D_{AB}=1}La solución de Blasius se aplica directamente.

En consecuencia, esta derivación utiliza una forma relacionada de las condiciones de contorno, reemplazandov{\displaystyle v}conT{\displaystyle T}odoA{\displaystyle c_{A}}(temperatura absoluta o concentración de la especie A). El subíndice S denota una condición de superficie.

vincógnitavSvvS=TTSTTS=doAdoASdoAdoAS=0{\displaystyle {v_{x}-v_{S} \over v_{\infty }-v_{S}}={T-T_{S} \over T_{\infty }-T_{S}}={c_{A}-c_{AS} \over c_{A\infty }-c_{AS}}=0}     en     y=0{\displaystyle y=0}
vincógnitavSvvS=TTSTTS=doAdoASdoAdoAS=1{\displaystyle {v_{x}-v_{S} \over v_{\infty }-v_{S}}={T-T_{S} \over T_{\infty }-T_{S}}={c_{A}-c_{AS} \over c_{A\infty }-c_{AS}}=1}     en     y={\displaystyle y=\infty }yincógnita=0{\displaystyle x=0}

Utilizando la función de línea de corriente, Blasius obtuvo la siguiente solución para el esfuerzo cortante en la superficie de la placa.

τ0=(vincógnitay)y=0=0,332vincógnitaRmi1/2{\displaystyle \tau _{0}=\left({\partial v_{x} \over \partial y}\right)_{y=0}=0.332{v_{\infty } \over x}Re^{1/2}}

Y a través de las condiciones de contorno, se sabe que

vincógnitavSvvS=TTSTTS=doAdoASdoAdoAS{\displaystyle {v_{x}-v_{S} \over v_{\infty }-v_{S}}={T-T_{S} \over T_{\infty }-T_{S}}={c_{A}-c_{AS} \over c_{A\infty }-c_{AS}}}

Se nos proporcionan las siguientes relaciones para el flujo de calor/masa que sale de la superficie de la placa.

(Ty)y=0=0,332TTSincógnitaRmi1/2{\displaystyle \left({\partial T \over \partial y}\right)_{y=0}=0.332{T_{\infty }-T_{S} \over x}Re^{1/2}}
(doAy)y=0=0,332doAdoASincógnitaRmi1/2{\displaystyle \left({\partial c_{A} \over \partial y}\right)_{y=0}=0.332{c_{A\infty }-c_{AS} \over x}Re^{1/2}}

Entonces paraPAGr=Sdo=1{\displaystyle Pr=Sc=1}

δ=δT=δdo=5.0incógnitaRmi{\displaystyle \delta =\delta _{T}=\delta _{c}={5.0x \over {\sqrt {Re}}}}

dóndeδT,δdo{\displaystyle \delta _{T},\delta _{c}}son las regiones de flujo dondeT{\displaystyle T}ydoA{\displaystyle c_{A}}son menos del 99% de sus valores de campo lejano. [ 27 ]

Debido a que el número de Prandtl de un fluido en particular no suele ser igual a la unidad, el ingeniero alemán E. Polhausen, que trabajó con Ludwig Prandtl, intentó extender empíricamente estas ecuaciones para aplicarlas aPAGr1{\displaystyle Pr\neq 1}Sus resultados pueden aplicarse aSdo{\displaystyle Sc}Asimismo. [ 28 ] Encontró que para un número de Prandtl mayor que 0,6, el espesor de la capa límite térmica venía dado aproximadamente por:

El número de Prandtl afecta el espesor de la capa límite térmica. Cuando Prandtl es menor que 1, la capa térmica es más gruesa que la capa de velocidad. Cuando Prandtl es mayor que 1, la capa térmica es más delgada que la capa de velocidad.
Gráfico que muestra el espesor relativo en la capa límite térmica frente a la capa límite de velocidad (en rojo) para varios números de Prandtl. ParaPAGr=1{\displaystyle Pr=1}, ambos son iguales.
δδT=PAGr1/3{\displaystyle {\delta \over \delta _{T}}=Pr^{1/3}}          y por lo tanto          δδdo=Sdo1/3{\displaystyle {\delta \over \delta _{c}}=Sc^{1/3}}

A partir de esta solución, es posible caracterizar las constantes de transferencia de calor/masa por convección en función de la región del flujo de la capa límite. La ley de conducción de Fourier y la ley de enfriamiento de Newton se combinan con el término de flujo derivado anteriormente y el espesor de la capa límite.

qA=k(Ty)y=0=hincógnita(TST){\displaystyle {q \over A}=-k\left({\partial T \over \partial y}\right)_{y=0}=h_{x}(T_{S}-T_{\infty })}
hincógnita=0,332kincógnitaRmiincógnita1/2PAGr1/3{\displaystyle h_{x}=0.332{k \over x}Re_{x}^{1/2}Pr^{1/3}}

Esto proporciona la constante convectiva local.hincógnita{\displaystyle h_{x}}en un punto del plano semiinfinito. La integración a lo largo de la placa da un promedio

hL=0,664kincógnitaRmiL1/2PAGr1/3{\displaystyle h_{L}=0.664{k \over x}Re_{L}^{1/2}Pr^{1/3}}

Siguiendo la derivación con términos de transferencia de masa (k{\displaystyle k}= constante de transferencia de masa convectiva,DAB{\displaystyle D_{AB}}= difusividad de la especie A en la especie B,Sdo=ν/DAB{\displaystyle Sc=\nu /D_{AB}}), se obtienen las siguientes soluciones:

kincógnita=0,332DABincógnitaRmiincógnita1/2Sdo1/3{\displaystyle k'_{x}=0.332{D_{AB} \over x}Re_{x}^{1/2}Sc^{1/3}}
kL=0,664DABincógnitaRmiL1/2Sdo1/3{\displaystyle k'_{L}=0.664{D_{AB} \over x}Re_{L}^{1/2}Sc^{1/3}}

Estas soluciones se aplican para flujo laminar con un número de Prandtl/Schmidt mayor que 0,6. [ 27 ]

Muchos de los principios que se aplican a las aeronaves también se aplican a los barcos, submarinos y plataformas marinas, siendo el agua el fluido principal de interés en lugar del aire. Dado que el agua no es un fluido ideal, los barcos que se mueven en ella experimentan resistencia. Las partículas del fluido se adhieren al casco del barco debido a la fuerza de adhesión entre el agua y el barco, creando una capa límite donde la velocidad de flujo del fluido forma un gradiente de velocidad pequeño pero pronunciado , con el fluido en contacto con el barco idealmente con una velocidad relativa de 0, y el fluido en el borde de la capa límite con la velocidad de flujo libre , o la velocidad relativa del fluido alrededor del barco. [ 29 ]

Mientras que la parte delantera del buque se enfrenta a fuerzas de presión normales debido al fluido que la rodea, la parte trasera experimenta una componente de presión menor debido a la capa límite. Esto genera una mayor resistencia a la presión, conocida como «resistencia viscosa por presión» o « resistencia de forma ». [ 29 ]

En los buques, a diferencia de las aeronaves, se trabaja con flujos incompresibles, donde la variación de la densidad del agua es insignificante (un aumento de presión cercano a 1000 kPa produce una variación de tan solo 2-3  kg/m³ ) . Este campo de la dinámica de fluidos se denomina hidrodinámica. Un ingeniero naval diseña priorizando la hidrodinámica y, posteriormente, la resistencia. El desarrollo, la ruptura y la separación de la capa límite se vuelven cruciales debido a que la alta viscosidad del agua genera elevadas tensiones de corte.

Turbina de capa límite

Este efecto se aprovechó en la turbina Tesla , patentada por Nikola Tesla en 1913. Se la conoce como turbina sin álabes porque utiliza el efecto de la capa límite y no un fluido que incide sobre los álabes como en una turbina convencional. Las turbinas de capa límite también se conocen como turbinas de cohesión, turbinas sin álabes y turbinas de capa de Prandtl (en honor a Ludwig Prandtl ).

Predicción del espesor transitorio de la capa límite en un cilindro mediante análisis dimensional

Al utilizar las ecuaciones de fuerza transitoria y viscosa para un flujo cilíndrico, se puede predecir el espesor de la capa límite transitoria hallando el número de Womersley (nortew{\displaystyle N_{w}}).

Fuerza transitoria =ρvw{\displaystyle \rho vw}

Fuerza viscosa =μvδ12{\displaystyle {\mu v \over \delta _{1}^{2}}}

Igualándolos entre sí se obtiene:

ρvw=μvδ12{\displaystyle \rho vw={\mu v \over \delta _{1}^{2}}}

Al despejar delta se obtiene:

δ1=μρw= v w{\displaystyle \delta _{1}={\sqrt {\mu \over \rho w}}={\sqrt {\ v \over \ w}}}

En forma adimensional:

Lδ1=Lw v=nortew{\displaystyle {L \over \delta _{1}}={L{\sqrt {w \over \ v}}}=N_{w}}

dóndenortew{\displaystyle N_{w}}= Número de Womersley;ρ{\displaystyle \rho }= densidad;v{\displaystyle v}= velocidad;w={\displaystyle w=}frecuencia de oscilaciones;δ1{\displaystyle \delta _{1}}= longitud de la capa límite transitoria;μ{\displaystyle \mu }= viscosidad;L{\displaystyle L}= longitud característica.

Predicción de las condiciones de flujo convectivo en la capa límite de un cilindro mediante análisis dimensional.

Al utilizar las ecuaciones de fuerza convectiva y viscosa en la capa límite para un flujo cilíndrico, se pueden predecir las condiciones de flujo convectivo en la capa límite hallando el número de Reynolds adimensional (Rmi{\displaystyle Re}).

Fuerza convectiva:ρv2 L{\displaystyle \rho v^{2} \over \ L}

Fuerza viscosa:μvδ22{\displaystyle {\mu v \over \delta _{2}^{2}}}

Igualándolos entre sí se obtiene:

ρv2 L=μvδ22{\displaystyle {\rho v^{2} \over \ L}={\mu v \over \delta _{2}^{2}}}

Al despejar delta se obtiene:

δ2=μLρv{\displaystyle \delta _{2}={\sqrt {\mu L \over \rho v}}}

En forma adimensional:

Lδ2=ρvLμ=Rmi{\displaystyle {L \over \delta _{2}}={\sqrt {\rho vL \over \mu }}={\sqrt {Re}}}

dóndeRmi{\displaystyle Re}= Número de Reynolds;ρ{\displaystyle \rho }= densidad;v{\displaystyle v}= velocidad;δ2{\displaystyle \delta _{2}}= longitud de la capa límite convectiva;μ{\displaystyle \mu }= viscosidad;L{\displaystyle L}= longitud característica.

Ingestión de la capa límite

La ingestión de la capa límite promete un aumento en la eficiencia del combustible de las aeronaves mediante un propulsor montado en la parte trasera que ingiere la capa límite lenta del fuselaje y reenergiza la estela para reducir la resistencia y mejorar la eficiencia propulsora . Para operar en un flujo de aire distorsionado, el ventilador es más pesado y su eficiencia se reduce, y su integración es compleja. Se utiliza en conceptos como el Aurora D8 o el Nova de la agencia de investigación francesa Onera , ahorrando un 5 % en crucero al ingerir el 40 % de la capa límite del fuselaje. [ 30 ]

Airbus presentó el concepto Nautilus en el congreso ICAS en septiembre de 2018: para absorber toda la capa límite del fuselaje, minimizando la distorsión del flujo azimutal , el fuselaje se divide en dos ejes con ventiladores de relación de derivación de 13-18:1. La eficiencia propulsora alcanza hasta el 90%, como la de los rotores abiertos contrarrotatorios, con motores más pequeños, ligeros, menos complejos y ruidosos. Podría reducir el consumo de combustible en más de un 10% en comparación con un motor convencional bajo el ala con una relación de derivación de 15:1. [ 30 ]

Véase también

Referencias

  1. Young, AD (1989). Capas límite (1.ª  ed. publicada). Washington, DC: Instituto Americano de Aeronáutica y Astronáutica. ISBN 0930403576.
  2. Schlichting, Hermann; Gersten, Klaus (2017). "2.1 Concepto de capa límite" . Teoría de la capa límite (Novena ed.). Berlín Heidelberg: Springer. pag. 29.doi : 10.1007 /978-3-662-52919-5_2 . ISBN   978-3-662-52917-1Consultado el 5 de agosto de 2023. Con frecuencia, el límite se define arbitrariamente como el punto donde la velocidad alcanza un cierto porcentaje de la velocidad exterior, por ejemplo, el 99 %. Para mayor claridad, se suele utilizar un índice, por ejemplo, δ99.
  3. ^ Prandtl, L. (1938). "Zur Berechnung der Grenzschichten". Zeitschrift für Angewandte Mathematik und Mechanik . 18 (1): 77– 82. Bibcode : 1938ZaMM...18...77P . doi : 10.1002/zamm.19380180111 .
  4. Van Dyke, Milton. Métodos de perturbación en mecánica de fluidos. Parabolic Press, Incorporated, 1975.
  5. Stewartson, K. (1957). "Sobre expansiones asintóticas en la teoría de capas límite". Journal of Mathematics and Physics . 36 ( 1– 4): 173– 191. doi : 10.1002/sapm1957361173 .
  6. Libby, Paul A.; Fox, Herbert (1963). "Algunas soluciones de perturbación en la teoría de la capa límite laminar". Journal of Fluid Mechanics . 17 (3): 433. doi : 10.1017/S0022112063001439 . S2CID 123824364 . 
  7. Fox, Herbert; Libby, Paul A. (1964). "Algunas soluciones de perturbación en la teoría de la capa límite laminar Parte 2. La ecuación de energía". Journal of Fluid Mechanics . 19 (3): 433– 451. Bibcode : 1964JFM....19..433F . doi : 10.1017/S0022112064000830 . S2CID 120911442 . 
  8. von Kármán, T. (1921). "Über laminare und turbulente Reibung" . Zeitschrift für Angewandte Mathematik und Mechanik . 1 (4): 233– 252. Bibcode : 1921ZaMM....1..233K . doi : 10.1002/zamm.19210010401 .
  9. Wieghardt, K. Sobre una ecuación de energía para el cálculo de capas límite laminares. Joint Intelligence Objectives Agency, 1946.
  10. ^ Wieghardt, K. (1948). "Über einen Energiesatz zur Berechnung laminarer Grenzschichten". Ingenieur-Archiv . 16 ( 3– 4): 231– 242. Bibcode : 1948AAM....16..231W . doi : 10.1007/BF00548007 . S2CID 119750449 . 
  11. Rosenhead, Louis, ed. Capas límite laminares. Clarendon Press, 1963.
  12. ^ Tollmien, Walter; Schlichting, Hermann; Görtler, Henry; Riegels, FW (1961). "Bemerkungen zur Hydrodynamik". Ludwig Prandtl Gesammelte Abhandlungen . págs. 627– 631. doi : 10.1007/978-3-662-11836-8_49 . ISBN  978-3-662-11837-5.{{cite book}}: Incompatibilidad de ISBN/Fecha ( ayuda )
  13. von Kármán, T.; Tsien, HS (1938). "Capa límite en fluidos compresibles". Journal of the Aeronautical Sciences . 5 (6): 227– 232. doi : 10.2514/8.591 .
  14. Crocco, L. "Una transformación característica de las ecuaciones de la capa límite en gases." ARC 4582 (1939): 1940.
  15. von Karman, T. (1939). "La analogía entre la fricción de fluidos y la transferencia de calor". Transactions of the American Society of Mechanical Engineers . 61 (8): 705– 710. doi : 10.1115/1.4021298 . S2CID 256805665 . 
  16. Guo, J.; Yang, XIA; Ihme, M. (marzo de 2022). "Estructura de la capa límite térmica en flujos de canal turbulentos en condiciones transcríticas" . Journal of Fluid Mechanics . 934 A45. Bibcode : 2022JFM...934A..45G . doi : 10.1017/jfm.2021.1157 . ISSN 0022-1120 . S2CID 246066677 .  
  17. ^ Levêque, A. (1928). "Las leyes de la transmisión de calor por convección". Annales des Mines ou Recueil de Mémoires sur l'Exploitation des Mines et sur les Sciences et les Arts qui s'y Rattachent, Mémoires (en francés). XIII (13): 201-239 .
  18. 1 2 Niall McMahon. "André Lévêque p285, una revisión de su aproximación del perfil de velocidad" . Archivado del original el 4 de junio de 2012.
  19. 1 2 Martin, H. (2002). "La ecuación generalizada de Lévêque y su uso práctico para la predicción de tasas de transferencia de calor y masa a partir de la caída de presión". Chemical Engineering Science . 57 (16): 3217– 3223. Bibcode : 2002ChEnS..57.3217M . doi : 10.1016/S0009-2509(02)00194-X .
  20. Schuh, H. (1953). "Sobre soluciones asintóticas para la transferencia de calor a temperaturas variables de la pared en una capa límite laminar con perfiles de velocidad de Hartree". Journal of the Aeronautical Sciences . 20 (2): 146– 147. doi : 10.2514/8.2566 .
  21. Kestin, J. y Persen, LN (1962). "La transferencia de calor a través de una capa límite turbulenta a números de Prandtl muy altos". International Journal of Heat and Mass Transfer . 5 (5): 355– 371. Bibcode : 1962IJHMT...5..355K . doi : 10.1016/0017-9310(62)90026-1 .
  22. Schlichting, H. (1979). Teoría de la capa límite (7.ª ed.). Nueva York (EE. UU.): McGraw-Hill. 
  23. Barna, Imre Ferenc; Bognár, Gabriella; Mátyás, László; Hriczó, Krisztián (2022). "Análisis autosimilar de las capas límite compresibles e incompresibles dependientes del tiempo, incluida la conducción de calor" . Revista de Análisis Térmico y Calorimetría . 147 : 13625–13632 . arXiv : 2101.08990 . doi : 10.1007/s10973-022-11574-3 .
  24. Rossa, Carlos (2023). "Pérdidas de energía en generadores fotovoltaicos debido a patrones de viento" . Nature Communications Engineering . 2 (66) 66. Bibcode : 2023CmEng...2...66R . doi : 10.1038/s44172-023-00119-7 . PMC 10956078 . 
  25. ^ Blasio, H. (1908). "Grenzschichten in Flüssigkeiten mit kleiner Reibung" . Zeitschrift für Mathematik und Physik . 56 : 1-37 .( Traducción al inglés )
  26. Martin, Michael J. (2001). "Solución de la capa límite de Blasius con condiciones de flujo deslizante". AIP Conference Proceedings . Vol. 585. pp. 518– 523. doi : 10.1063/1.1407604 . hdl : 2027.42/87372 .  
  27. 1 2 Geankoplis, Christie J. Procesos de transporte y principios de procesos de separación: (incluye operaciones unitarias). Cuarta ed. Upper Saddle River, NJ: Prentice Hall Professional Technical Reference, 2003. Impreso.
  28. Pohlhausen, E. (1921). "Der Wärmeaustausch zwischen festen Körpern und Flüssigkeiten mit kleiner reibung und kleiner Wärmeleitung" . Zeitschrift für Angewandte Mathematik und Mechanik . 1 (2): 115– 121. Bibcode : 1921ZaMM....1..115P . doi : 10.1002/zamm.19210010205 .
  29. 1 2 "Resistencia y propulsión de buques" (PDF) . usna.edu . Consultado el 14 de febrero de 2024 .
  30. 1 2 Graham Warwick (19 de noviembre de 2018). "La semana en tecnología, del 19 al 23 de noviembre de 2018" . Aviation Week & Space Technology .
  • Chanson, H. (2009). Hidrodinámica aplicada: Introducción a los flujos de fluidos ideales y reales . CRC Press, Taylor & Francis Group, Leiden, Países Bajos, 478 páginas. ISBN 978-0-415-49271-3.
  • AD Polyanin y VF Zaitsev, Manual de ecuaciones diferenciales parciales no lineales , Chapman & Hall/CRC Press, Boca Raton – Londres, 2004. ISBN 1-58488-355-3
  • AD Polyanin, AM Kutepov, AV Vyazmin y DA Kazenin, Hidrodinámica, transferencia de masa y calor en ingeniería química , Taylor & Francis, Londres, 2002. ISBN 0-415-27237-8
  • Hermann Schlichting, Klaus Gersten, E. Krause, H. Jr. Oertel, C. Mayes "Teoría de la capa límite" octava edición Springer 2004 ISBN 3-540-66270-7
  • John D. Anderson Jr., "La capa límite de Ludwig Prandtl" , Physics Today , diciembre de 2005
  • Anderson, John (1992). Fundamentos de aerodinámica (2.ª  ed.). Toronto: SSCHAND. págs. 711–714 . ISBN  0-07-001679-8.
  • H. Tennekes y JL Lumley , "Un primer curso sobre turbulencia", The MIT Press, (1972).
  • Conferencias sobre turbulencia para el siglo XXI, por William K. George
  • Biblioteca Digital Nacional de Ciencias – Capa límite
  • Moore, Franklin K., " Efecto de desplazamiento de una capa límite tridimensional ". Informe NACA 1124, 1953.
  • Benson, Tom, " Capa límite ". Tecnologías de aprendizaje Glenn de la NASA.
  • Separación de la capa límite
  • Ecuaciones de la capa límite: Soluciones exactas – de EqWorld
  • Jones, TRANSFERENCIA DE CALOR EN LA CAPA LÍMITE DE TV
  • "El revolucionario concepto de "capa límite" y su prevalencia en la aeronáutica, por Sourabh S. Diwan" . YouTube . Centro Internacional de Ciencias Teóricas. 18 de febrero de 2022.