Ecuaciones de la capa límite

y capa límite sobre una

superficie plana

•Jose Julian Marín Garcia

•Sebastián Morales Espina
•Anet Xcaret Campos Colima
•Rubén Rangel Figueroa
•Luis Alan Leiva Garcia
Capa Límite

• La capa límite es la capa de

fluido que está directamente en
contacto con una superficie; justo
en ese límite hay una interacción
donde la velocidad del flujo
disminuye debido a la fricción
con la superficie.
Representación gráfica
de las diferentes partes
de la capa límite
Placa plana
paralela al
flujo para
distintos
números de
Reynolds.
El número de Reynolds se
define:
Importancia en el Transporte de Fluidos

• En el transporte de fluidos, ya sea en tuberías, aeronaves o vehículos, la capa límite tiene un

impacto significativo en la resistencia al flujo (también conocida como drag) y en la eficiencia
del sistema. Dependiendo del tipo de flujo (laminar o turbulento), la capa límite puede influir en
los siguientes aspectos:

Resistencia al flujo Desprendimiento de la capa límite

La capa límite determina cómo el fluido

interactúa con las superficies, afectando la Si el flujo de un fluido es demasiado rápido o
cantidad de fuerza de fricción que se ejerce cambia de dirección, puede ocurrir un
sobre un objeto en movimiento. desprendimiento de la capa límite
Aplicaciones
prácticas
Factores
que
Afectan la
Capa
Límite
Velocidad
del flujo
•A mayores velocidades, la
capa límite tiende a
volverse más gruesa y
puede transitar de laminar
a turbulenta, lo que
incrementa la resistencia
 •La viscosidad del fluido
juega un papel importante
Propiedad en el desarrollo de la capa
límite. Fluidos más viscosos
es del tienden a mantener un flujo
laminar, mientras que los
fluido fluidos menos viscosos
pueden desarrollar un flujo
turbulento más fácilmente.
Forma de la
superficie
•La geometría de la superficie afecta el
comportamiento de la capa límite.
Superficies rugosas o irregulares pueden
causar desprendimiento más temprano de
la capa límite, mientras que superficies
lisas favorecen el flujo laminar.
Ecuaciones de la capa
límite y capa límite
sobre una superficie
plana
 Actualmente se utilizan tres técnicas para
tratar los flujos externos:

(1) soluciones numéricas en

(2) experimentación (3) teoría de la capa límite
ordenador
la teoría de la capa límite,
formulada por Ludwig Prandtl
en 1904.

• Físico alemán que se especializó en el estudio

de la mecánica de fluidos.
Ecuaciones De
La Capa Limite.
• Las ecuaciones que
gobiernan el flujo
viscoso en la capa
límite son las
ecuaciones de Navier
Stokes
En el caso de una placa plana se utiliza un flujo bidimensional por lo que podemos
simplificarlo a que sea coordenadas cartesianas.
Con esto podemos descartar la componente en z
Entonces queda de la siguiente manera
Además se cuenta con la ecuación de
continuidad
Dado que trabajamos con un
flujo bidimensional solo la
usaremos así:
Como el efecto de la gravedad es
muy pequeño lo podemos
eliminar, quedando la ecuación de
esta forma
En este caso también estaremos eliminando
los últimos términos de Z queda:
Tampoco vamos a trabajar con términos que les afecte el tiempo es
decir nuestro flujo será estacionario por lo que podemos realizar
otra simplificación
Basándonos en esta imagen diremos lo siguiente:

Que u es muy próxima a la

velocidad U de nuestra placa.
u≈U

Que ∆x es muy próxima a la

longitud L de nuestra placa.
∆x≈L

Que ∆y es muy próximo a nuestro

espesor δ de capa limite.
∆y≈δ
La velocidad ubicada en la
coordenada de y, la cual es producida
por un desviamiento de las corrientes
de nuestro fluido como se puede
observar en la siguiente imagen.
 Haciendo uso de la
ecuación de
continuidad antes • Teniendo así la ecuación de • Sustituyendo nuestras
vista vamos a darle continuidad después de aproximaciones nos quedaría
un valor a nuestra eliminar la variable Z
velocidad, el cual
será necesario para
analizar nuestros
términos.
u es muy próxima a la velocidad U de nuestra
placa.
u≈U

∆x es muy próxima a la longitud L de nuestra

placa.
∆x≈L

∆y es muy próximo a nuestro espesor δ de

capa limite.
∆y≈δ
 Recordemos que no tiene como
tal un valor dado cada termino
Dejamos nuestro termino de solo estamos trabajando con
ordenes de magnitud es decir le
velocidad solo de un lado de damos un valor relativo por lo
la ecuación para así que los signos no los vamos a
conseguir su equivalencia, tener en consideración.
dejando la ecuación de la
siguiente forma
Regresando a la ecuación de Navier-Stocks también
realizaremos un análisis con ordenes magnitud

• Empezando por nuestra direccion “x" • Sustituyendo queda:

Aquí nos podemos dar cuenta que del lado izquierdo de nuestra ecuación podemos
cancelar a nuestro espesor δ, dejándonos apreciar que estos dos términos del lado
izquierdo son iguales.

• Ecuación • Quedándonos de esta forma:

 • Nuestra velocidad U esta
dividida entre dos
diferentes magnitudes y
Ahora enfocándonos por el análisis que
hicimos en la imagen
en el lado derecho antes vista sabemos que
L es mucho mayor a δ.
de la ecuacíon Entonces podemos
despreciar nuestro
termino dividido por L
Hecho esto podemos regresarnos a la
ecuación original
• Vamos a despejar del lado izquierdo a 𝜌 la cual está
multiplicando para que nos quede de la siguiente forma.
Sabemos que la viscosidad cinemática equivale a µ / 𝜌 por lo
que podemos sustituirla en nuestra ecuación.
Ecuación de
Bernoulli
•Para la presión nos
ayudaremos de esta
ecuación
 Tomando en consideración que la presión afuera de la capa
limite no tiene variaciones y adentro de la capa limite si puede
tener variaciones respecto de x, nos quedaría lo siguiente
Siendo esta la ecuación de Bernoulli la cual ya fue
deducida para las condiciones que se tiene en la
capa limite nos deja con la siguiente igualdad
Sustituiremos en la ecuación que hemos venido trabajando dejándonos con la
siguiente expresión.
 • Para hacer más fácil este análisis solo
Para la sustituiremos la orden de magnitud que
está dividiendo a nuestros términos, debido
• Haciendo la sustitución en la
ecuación nos queda lo siguiente:

dirección en a que ya sabemos que nuestra υ es

pequeña.
"y"
• Como podemos observar al
ser nuestra υ pequeña y
dividirla por alguna otra
cantidad, nos resulta algo
aún más pequeño por lo
que lo podemos despreciar.

• Quedándonos así solamente el

término de la presión la cual
también la vamos a despreciar
debido a que no tiene
variaciones en esa dirección
ósea que es constante y suele
ser muy pequeña por lo que la
aproximamos a cero.
Dejándonos sin ningún termino
para la dirección “y”. con esto
finalizamos los análisis en cada una
de las direcciones propuestas
dejándonos finalmente con las
ecuaciones de Prandtl que serían
las siguientes:
Ecuaciones de
Blasius
Físico alemán especializado en mecánica de fluidos. Fue
uno de los primeros estudiantes de Prandtl.
Blasius comenzó con el caso mas simple que seria una placa plana de espesor
cero alineada con la corriente junto con una velocidad constante, sin
perturbaciones es decir laminar y sin gradiente de presión
• Blasius toma directamente estas
ecuaciones y se observa que existe un
sistema de ecuaciones diferenciales
de primer orden en x y de segundo
orden en y, por definición
matemática sabemos que se puede
llevar dos ecuaciones de segundo
orden a una de tercer orden y
viceversa.
Esta hipótesis es que cualquier perfil
de velocidades que tengamos en la
placa, podemos hacerlo coincidir
con otro existente a lo largo de la
misma usando una escala adecuada.
 Para u es sencillo pues esta va a ser la velocidad de corriente libre representada por U
Pero para Y, debería ser uno que me de él orden de magnitud del espesor δ de la capa
limite. Pero como no lo conocemos debemos calcular una longitud que nos de ese
orden de magnitud de delta δ y para esto vamos a realizar lo siguiente.

Sabemos que δ/L puede ser escalado

al orden de magnitud de reinols es
decir lo siguiente:

Donde:
ReL= al Reinols de la longitud
característica.
 Si nosotros queremos calcular en un punto cualquiera X de la capa limite el espesor δ el cual
va aumentando a lo largo de X entonces vamos a utilizar el Reinols local, es decir

• Pero el Reinols local es igual a lo • Por conveniencia diremos que nuestro espesor
siguiente: delta δ se aproxima a lo que equivale el
reinols local, es decir
 Esto siendo una orden de
magnitud no una igualdad y •u es la velocidad en la capa
vamos a tomar esa dimensión δ limite y U∞ es la velocidad de la
que es un factor de escala para corriente libre siendo esta
la coordenada Y para realizar expresión adimensional.
una relación adimensional la
cual es la siguiente
•g siendo una función la cual
tiene que ser universal que
describa la distribución de
velocidad en la capa limite
laminar de una placa plana.

•η (eta) la distancia a la pared Y,

Adimencionalizada con el factor
antes visto
Acomodando la ecuación que describe la distancia a la pared, para
poderla dejar de una forma conveniente para más adelante.

• Lo único que se hizo es invertir el denominador para usarlo de la manera que necesitamos.
• Primer paso lo que hizo blasius es que para • Donde ψ (epsi) es la función de corriente del
resolver las ecuaciones de prandtl definió flujo la cual a su vez seria equivalente a lo
una función de corriente adimensional f(η) siguiente
 Tenemos que empezar a calcular la velocidad de u y la velocidad de v y sus de sus derivadas
a partir de la función de corriente ψ para después remplazarlas en la ecuación de cantidad de
movimiento

• El cual va a ser igual a lo

• trabajaremos principalmente • Empezaremos con el termino siguiente
con la ecuación siguiente; u de nuestra ecuación
• Realizando un acomodo con • Ahora sustituiremos a lo que • Vamos a derivar respecto a lo
la regla de la cadena quedaría equivale ψ y η; que tenemos
así:

• Vamos a simplificar nuestro

resultado eliminando v y x debido a
que los tenemos en el denominador
y en el numerador
Ahora
continuaremos
con el • Quedaría de la • Sustituyendo los valores
siguiente siguiente manera
haciendo un arreglo
término de la con la regla de la
ecuación de cadena de la misma
manera que en el
cantidad de anterior caso • Ahora vamos a realizar la derivación
movimiento correspondiente para cada una de
nuestras expresiones en partes
• Quedándonos para η
para x vamos a sacar lo que es constante y nos
quedaremos simplemente con x que es nuestra
variable la cual sacamos de la raíz
Ahora que
tenemos nuestro
termino solo lo
derivaremos
Ahora vamos a acomodar los
términos para poder simplificar más
nuestro resultado quedándonos de la
siguiente manera;
Ahora vamos a acomodar los
términos para poder simplificar
más nuestro resultado
quedándonos de la siguiente
manera;

el siguiente termino es igual a η por lo

que podemos remplazarlo por la
misma para hacer una mayor
simplificación como se muestra a
continuación
Una vez terminada la simplificación de nuestros términos
podemos juntar ambos resultados quedándonos
c
 Continuaremos con el Sustituyendo los valores
siguiente termino de le
ecuación

En esta ocasión las derivadas van a

ser mas sencillas debido que en la
El cual también vamos a realizarle segunda derivada es respecto de y lo
un acomodo con la regla de la cadena cual simplifica mucho todo
de la siguiente forma
Siguiente término
• Tenemos que derivar dos veces respecto de y por lo que para hacer las cosas mas simples vamos a utilizar el
resultado anterior que ya fue derivado respecto a y, también aremos el mismo arreglo con la regla de la cadena

Derivando nos queda

Vamos a realizar el siguiente análisis respecto al producto de las raíces
que tenemos para simplificar nuestro resultado

• Aplicándolo a las raíces que tenemos

• simplificando más términos • ahora vamos a derivar para
podemos multiplicar nuestras U la componente en v la cual
dejándonos con una U2 equivale a lo siguiente
quedándonos debido a la definición de
corriente.

Sustituimos el valor de ψ en la ecuación

Como podemos observar existe un producto de funciones por lo que vamos a
realizar la derivada de acuerdo a esa fórmula de cuando tenemos un producto de dos
funciones
• Primero tomamos el término de la raíz y • Ahora le sumamos a(f(η)) y multiplicamos a la
derivamos respecto a(f(η)) y multiplicamos misma por la derivada de la raíz respecto de x
por la derivada de η respecto de x como se
muestra a continuación;
En la parte donde tenemos que derivar η lo que se hizo es que se
remplazó el valor de la misma y después se sacó la x de la raíz y se
invirtió para poderla derivar de manera más sencilla, al igual en la otra
parte se sacó x de la raíz para facilitar su derivación
Procedemos a derivar nuestros términos y nos
quedaría
Ahora con los términos arreglados de este modo podemos darnos cuenta que
podemos remplazar a η en nuestro resultado como se muestra a continuación
Tenemos raíces en ambas partes por lo que vamos a factorizarlas junto con el
1/2x el cual también lo tenemos en ambas partes

nuestro resultado
•Una vez encontrado todos los valores de
nuestras incógnitas vamos a sustituir en la
ecuación original
Así es como queda la ecuación sacando términos comunes
• Como podemos observar existen dos términos • Y a su vez eliminamos U∞ /x para quedar con la
iguales y de distinto signo por lo que vamos a expresión lo más reducido posible y Despejando y
eliminarlos agrupando nuestros términos
como es una ecuación de tercer orden vamos necesitar 3 condiciones para resolverla, teniendo en
cuenta esto tenemos que recordar que f es una función de corriente adimensional y que η es la
distancia adimensional a la pared

• Por lo tanto, las condiciones que nosotros conocemos para u y v son;

Todo esto representado para la variable Eta η queda de la siguiente manera;

tabla con los resultados numéricos de la solución de la ecuación de blasius por medio del
método de kutta, esta tabla es proveniente del libro de mecánica de fluidos del cengel;
Si nosotros adoptamos la convención de que el espesor de la capa limite es aquel donde
u/U∞=99% podemos observar en la tabla de resultados que “eta” η=5

despejamos delta δ y le remplazamos el valor de eta η

Dejándonos con la expresión para calcular el espesor de la

capa limite
 Distribución de tensiones de corte sobre la placa

vamos a derivar con cuando eta η vale 0, dejándonos

el arreglo de la regla sustituiremos y así que la solución para las
de la cadena derivaremos tensiones de corte seria
Coeficiente de fricción adimensional
PROGRAM_SOLUCION DE BLASIUS POR EL
METODO RUNGE KUTTA