Ir al contenido

Ecuación de Rayleigh-Plesset

De Wikipedia, la enciclopedia libre
La ecuación de Rayleigh–Plesset se suele aplicar al estudio de burbujas durante procesos de cavitación. En la imagen se representa la formación de estas burbujas tras una hélice.

En mecánica de fluidos, la ecuación de Rayleigh–Plesset es una ecuación diferencial ordinaria que gobierna la mecánica de una burbuja de gas inmersa en un líquido infinito.[1][2][3][4]​ Se suele escribir en su forma general como:

Símbolo Nombre
Presión en el interior de la burbuja, asumiéndose que esta es lo bastante pequeña como para que sea uniforme
Presión externa en el fluido a una distancia infinita
Densidad del líquido que rodea a la burbuja, asumiéndose constante
Radio de la burbuja
Viscosidad cinemática del fluido que rodea a la burbuja, considerándose constante
Tensión superficial de la burbuja

Siendo y conocidas, la ecuación de Rayleigh–Plesset puede ser usada para obtener el radio de la burbuja en función del tiempo.

Historia

[editar]

La ecuación de Rayleigh–Plesset se deriva de las ecuaciones de Navier-Stokes asumiendo simetría esférica.[4]​ Fue obtenida por primera vez por John Strutt, tercer barón Rayleigh en 1917, sin efectos de viscosidad ni tensión superficial. Fue usada por primera vez al estudio de burbujas viajeras en fenómenos de cavitación por Milton S. Plesset en 1949.[5]

Derivación

[editar]

La ecuación de Rayleigh–Plesset se puede obtener de las ecuaciones de Navier-Stokes con el radio de la burbuja como un parámetro dinámico.[3]​ Considerando simetría esférica en una burbuja de radio variable en el tiempo, se puede asumir que contiene vapor homogéneamente distribuido con una temperatura uniforme. En el exterior de la burbuja existe un dominio líquido de tamaño infinito con densidad constante y viscosidad dinámica , siendo la temperatura y presión lo bastante alejados de la burbuja como para que esta no afecte al líquido y . La temperatura se asume constante. Sin embargo, en las proximidades de la burbuja el fluido se ve afectado por esta, por lo que se puede definir parámetros en función de la distancia radial al centro de la burbuja. Son , y la velocidad . Es importante recordar que estos parámetros solo están definidos en el exterior de la burbuja, .

Conservación de la masa

[editar]

Aplicando la conservación de la masa, se obtiene una ley de la inversa del cuadrado para la velocidad que debe ser inversamente proporcional a la distancia del centro de la burbuja.[5]​ Así, se deduce que debe ser variable en el tiempo

Si la transferencia de masa a través de la superficie de la burbuja es nula, la velocidad en la interfase debe ser

lo que da:

Si hay transporte de masa, el incremento de masa encerrada en la burbuja es:

con representando el volumen de la burbuja. Si es la velocidad relativa del líquido con la burbuja en , la masa entrante a esta viene dada por:

con siendo la superficie de la burbuja. Aplicando la conservación de la masa, , se obtiene . Luego:

Así:

En muchos casos la densidad del líquido en mucho mayor que la del vapor, , por lo que puede ser aproximado por el primer resultado para transferencia de masa nula , por lo que[5]

Conservación del Momentum

[editar]

Si se asume un fluido newtoniano, las ecuaciones de Navier-Stokes en coordenadas esféricas para el movimiento en dirección radial son:

Sustituyendo por la viscosidad cinemática y reordenando los términos se obtiene:

Donde al sustituir por el resultado obtenido del apartado anterior:

Se debe notar que los términos viscosos se cancelan durante la sustitución:.[5]Separando variables e integrando desde la frontera de la burbuja hasta resulta:

Condiciones de frontera

[editar]

Designando como a la tensión normal en el líquido dirigida desde el centro hacia el exterior de la burbuja, tenemos para un fluido con densidad y viscosidad constantes:

Luego en una fracción infinitesimal de la superficie de la burbuja hay una fuerza resultante neta de:

donde es la tensión superficial.[5]​ Si no hay transferencia de masa en la frontera, la fuerza por unidad de área debe ser cero, luego:

y así resulta:

donde si se reordenan los términos y se define se obtiene la ecuación de Rayleigh–Plesset[5]

Usando la notación de Newton de indicar con un punto una derivada temporal, se puede representar más sucintamente como:

Soluciones

[editar]
Integración numérica de la ecuación de Rayleigh-Plesset completa. Se considera reposo en condiciones atmosféricas con R0=50 um, y la burbuja se somete a presión oscilatoria a su frecuencia natural hasta que se produce su colapso.
Integración numérica de la ecuación de Rayleigh-Plesset completa. Se considera reposo en condiciones atmosféricas con R0=50 um, y la burbuja se somete a subpresión hasta que se produce su colapso.

No se conocen soluciones cerradas para la ecuación de Rayleigh–Plesset. Sin embargo, se pueden obtener fácilmente soluciones numéricas con la precisión que se desee. Mención expresa merece el caso de tensión superficial y viscosidad negligibles, para el que hay aproximaciones analíticas de orden elevado.[6]

Para el caso estático, en cambio, la ecuación se simplifica a la conocida como ecuación de Laplace-Young:

Cuando sólo hay variaciones infinitesimales en el radio y presión, la ecuación da comoresultado la frecuencia natural de la burbuja, un valor de interés en los flujos con cavitación.

Referencias

[editar]
  1. Rayleigh, Lord (1917). «On the pressure developed in a liquid during the collapse of a spherical cavity». Phil. Mag. 34: 94-98. 
  2. Plesset, M.S. (1949). «The dynamics of cavitation bubbles». ASME J. Appl. Mech. 16: 228-231. 
  3. a b Leighton, T. G. (17 de abril de 2007). Derivation of the Rayleigh–Plesset equation in terms of volume. Southampton, UK: Institute of Sound and Vibration Research. 
  4. a b Lin, Hao; Brian D. Storey; Andrew J. Szeri (2002). «Inertially driven inhomogeneities in violently collapsing bubbles: the validity of the Rayleigh–Plesset equation». Journal of Fluid Mechanics 452. ISSN 0022-1120. doi:10.1017/S0022112001006693. Archivado desde el original el 8 de junio de 2019. Consultado el 27 de agosto de 2015. 
  5. a b c d e f Brennen, Christopher E. (1995). Cavitation and Bubble Dynamics. Oxford University Press. ISBN 0-19-509409-3. 
  6. Obreschkow, D.; Bruderer M.; Farhat, M. (5 de junio de 2012). «Analytical approximations for the collapse of an empty spherical bubble». Physical Review E 85. arXiv:1205.4202. doi:10.1103/PhysRevE.85.066303.