stringtranslate.com

Inestabilidad de Rayleigh-Taylor

Simulación hidrodinámica de un único "dedo" de la inestabilidad de Rayleigh-Taylor. [1] Nótese la formación de inestabilidades de Kelvin-Helmholtz , en la segunda instantánea y posteriores mostradas (comenzando inicialmente alrededor del nivel ), así como la formación de una "capa de hongo" en una etapa posterior en el tercer y cuarto fotograma de la secuencia.
Dedos de inestabilidad RT evidentes en la Nebulosa del Cangrejo

La inestabilidad de Rayleigh-Taylor , o inestabilidad RT (en honor a Lord Rayleigh y GI Taylor ), es una inestabilidad de una interfaz entre dos fluidos de diferentes densidades que ocurre cuando el fluido más ligero empuja al fluido más pesado. [2] [3] [4] Los ejemplos incluyen el comportamiento del agua suspendida sobre el petróleo en la gravedad de la Tierra , [3] nubes en forma de hongo como las de las erupciones volcánicas y las explosiones nucleares atmosféricas , [5] explosiones de supernovas en las que el gas del núcleo en expansión se acelera hasta convertirse en un gas de capa más denso, [6] [7] inestabilidades en reactores de fusión de plasma y [8] fusión por confinamiento inercial. [9]

El agua suspendida sobre petróleo es un ejemplo cotidiano de inestabilidad de Rayleigh-Taylor, y puede ser modelada por dos capas completamente paralelas de fluido inmiscible , el fluido más denso sobre el menos denso y ambos sujetos a la gravedad de la Tierra. El equilibrio aquí es inestable a cualquier perturbación o alteración de la interfaz: si una parcela de fluido más pesado se desplaza hacia abajo con un volumen igual de fluido más ligero desplazado hacia arriba, la energía potencial de la configuración es menor que el estado inicial. Por lo tanto, la perturbación crecerá y conducirá a una mayor liberación de energía potencial , ya que el material más denso se mueve hacia abajo bajo el campo gravitacional (efectivo), y el material menos denso se desplaza aún más hacia arriba. Esta fue la configuración estudiada por Lord Rayleigh. [3] La idea importante de GI Taylor fue su comprensión de que esta situación es equivalente a la situación en la que los fluidos se aceleran , con el fluido menos denso acelerando hacia el fluido más denso. [3] Esto ocurre en las profundidades submarinas en la superficie de una burbuja en expansión y en una explosión nuclear. [10]

A medida que se desarrolla la inestabilidad RT, las perturbaciones iniciales progresan desde una fase de crecimiento lineal a una fase de crecimiento no lineal, desarrollando eventualmente "columnas" que fluyen hacia arriba (en el sentido de flotabilidad gravitacional) y "picos" que caen hacia abajo. En la fase lineal, el movimiento del fluido puede aproximarse estrechamente mediante ecuaciones lineales , y la amplitud de las perturbaciones crece exponencialmente con el tiempo. En la fase no lineal, la amplitud de la perturbación es demasiado grande para una aproximación lineal, y se requieren ecuaciones no lineales para describir los movimientos del fluido. En general, la disparidad de densidad entre los fluidos determina la estructura de los flujos de inestabilidad RT no lineales posteriores (asumiendo que otras variables como la tensión superficial y la viscosidad son insignificantes aquí). La diferencia en las densidades del fluido dividida por su suma se define como el número de Atwood , A. Para A cercano a 0, los flujos de inestabilidad RT toman la forma de "dedos" simétricos de fluido; Para A cercano a 1, el fluido mucho más liviano "debajo" del fluido más pesado toma la forma de columnas más grandes similares a burbujas. [2]

Este proceso es evidente no sólo en muchos ejemplos terrestres, desde domos salinos hasta inversiones climáticas , sino también en astrofísica y electrohidrodinámica . Por ejemplo, la estructura de inestabilidad RT es evidente en la Nebulosa del Cangrejo , en la que la nebulosa de viento de púlsar en expansión impulsada por el púlsar del Cangrejo está barriendo material expulsado de la explosión de supernova hace 1000 años. [11] La inestabilidad RT también se ha descubierto recientemente en la atmósfera exterior del Sol, o corona solar , cuando una prominencia solar relativamente densa se superpone a una burbuja de plasma menos densa. [12] Este último caso se asemeja a las inestabilidades RT moduladas magnéticamente. [13] [14] [15]

Cabe señalar que la inestabilidad RT no debe confundirse con la inestabilidad Plateau-Rayleigh (también conocida como inestabilidad de Rayleigh) de un chorro de líquido. Esta inestabilidad, a veces llamada inestabilidad de manguera (o de manguera contra incendios), se produce debido a la tensión superficial, que actúa para romper un chorro cilíndrico en una corriente de gotitas que tienen el mismo volumen total pero una mayor área de superficie.

Muchas personas han presenciado la inestabilidad RT al observar una lámpara de lava , aunque algunos podrían afirmar que esto se describe con mayor precisión como un ejemplo de convección de Rayleigh-Bénard debido al calentamiento activo de la capa de fluido en la parte inferior de la lámpara.

Etapas de desarrollo y evolución final hacia la mezcla turbulenta

Esta figura representa la evolución de la inestabilidad de Rayleigh-Taylor a partir de pequeñas perturbaciones de longitud de onda en la interfaz (a) que crecen hasta convertirse en picos omnipresentes en forma de hongo (estructuras fluidas de fluido pesado en fluido ligero) y burbujas (estructuras fluidas de fluido ligero en fluido pesado) (b) y estas estructuras fluidas interactúan debido a la fusión y competencia de las burbujas (c) y finalmente se desarrollan en una región de mezcla (d). Aquí ρ2 representa el fluido pesado y ρ1 representa el fluido ligero. La gravedad actúa hacia abajo y el sistema es inestable en sentido horario.

La evolución de la RTI sigue cuatro etapas principales. [2] En la primera etapa, las amplitudes de perturbación son pequeñas en comparación con sus longitudes de onda, las ecuaciones de movimiento se pueden linealizar, lo que resulta en un crecimiento exponencial de la inestabilidad. En la primera parte de esta etapa, una perturbación inicial sinusoidal conserva su forma sinusoidal. Sin embargo, después del final de esta primera etapa, cuando comienzan a aparecer efectos no lineales, se observa el comienzo de la formación de los omnipresentes picos en forma de hongo (estructuras fluidas de fluido pesado que se transforman en fluido ligero) y burbujas (estructuras fluidas de fluido ligero que se transforman en fluido pesado). El crecimiento de las estructuras de hongo continúa en la segunda etapa y se puede modelar utilizando modelos de arrastre por flotabilidad, lo que resulta en una tasa de crecimiento que es aproximadamente constante en el tiempo. En este punto, los términos no lineales en las ecuaciones de movimiento ya no se pueden ignorar. Los picos y las burbujas comienzan a interactuar entre sí en la tercera etapa. Se produce una fusión de burbujas, en la que la interacción no lineal del acoplamiento de modos actúa para combinar picos y burbujas más pequeños para producir otros más grandes. También se produce una competencia de burbujas, en la que los picos y burbujas de longitud de onda más pequeña que se han saturado quedan envueltos por otros más grandes que aún no se han saturado. Esto finalmente se convierte en una región de mezcla turbulenta, que es la cuarta y última etapa de la evolución. En general, se supone que la región de mezcla que finalmente se desarrolla es autosimilar y turbulenta, siempre que el número de Reynolds sea lo suficientemente grande. [16]

Análisis de estabilidad lineal

Estado base de la inestabilidad de Rayleigh-Taylor. La gravedad apunta hacia abajo.

La inestabilidad bidimensional no viscosa de Rayleigh-Taylor (RT) proporciona un excelente trampolín hacia el estudio matemático de la estabilidad debido a la naturaleza simple del estado base. [17] Considere un estado base en el que hay una interfaz, ubicada en que separa medios fluidos con diferentes densidades, para y para . La aceleración gravitacional se describe mediante el vector . El campo de velocidad y el campo de presión en este estado de equilibrio, denotados con una barra superior, están dados por

donde la ubicación de referencia para la presión se toma en . Dejemos que esta interfaz se perturbe ligeramente, de modo que asuma la posición . En consecuencia, el estado base también se perturba ligeramente. En la teoría lineal, podemos escribir

donde es el número de onda real en la dirección y es la tasa de crecimiento de la perturbación. Entonces, el análisis de estabilidad lineal basado en las ecuaciones de gobierno de la no viscosa muestra que

Por lo tanto, si , el estado base es estable y mientras que si , es inestable para todos los números de onda. Si la interfaz tiene una tensión superficial , entonces la relación de dispersión se convierte en

lo que indica que la inestabilidad ocurre solo para un rango de números de onda donde ; es decir, la tensión superficial estabiliza números de onda grandes o escalas de longitud pequeñas. Entonces la tasa de crecimiento máxima ocurre en el número de onda y su valor es

Detalles del análisis de estabilidad lineal. [17] Una derivación similar aparece en Chandrasekhar 1981, §92, pp. 433–435.

La perturbación introducida en el sistema se describe mediante un campo de velocidad de amplitud infinitesimalmente pequeña. Debido a que se supone que el fluido es incompresible, este campo de velocidad tiene la representación de función de corriente

donde los subíndices indican derivadas parciales . Además, en un fluido incompresible inicialmente estacionario, no hay vorticidad y el fluido permanece irrotacional , por lo tanto . En la representación de la función de corriente, A continuación, debido a la invariancia traslacional del sistema en la dirección x , es posible hacer el ansatz

donde es un número de onda espacial. Por lo tanto, el problema se reduce a resolver la ecuación

El dominio del problema es el siguiente: el fluido con la etiqueta 'L' vive en la región , mientras que el fluido con la etiqueta 'G' vive en el semiplano superior . Para especificar la solución completamente, es necesario fijar las condiciones en los límites y la interfaz. Esto determina la velocidad de onda c , que a su vez determina las propiedades de estabilidad del sistema.

La primera de estas condiciones se proporciona mediante detalles en el límite. Las velocidades de perturbación deben satisfacer una condición de ausencia de flujo, de modo que el fluido no se escape en los límites . Por lo tanto, en , y en . En términos de la función de corriente, esto es

Las otras tres condiciones se proporcionan mediante detalles en la interfaz .

Continuidad de la velocidad vertical: En , las velocidades verticales coinciden, . Utilizando la representación de la función de corriente, esto da

Ampliando sobre da

donde HOT significa "términos de orden superior". Esta ecuación es la condición de interfase requerida.

Condición de superficie libre: En la superficie libre , se cumple la condición cinemática:

Linealizando, esto es simplemente

donde la velocidad se linealiza sobre la superficie . Utilizando las representaciones de modo normal y función de corriente, esta condición es , la segunda condición interfacial.

Relación de presión a través de la interfaz: Para el caso de tensión superficial , la diferencia de presión sobre la interfaz en está dada por la ecuación de Young-Laplace :

donde σ es la tensión superficial y κ es la curvatura de la interfaz, que en una aproximación lineal es

De este modo,

Sin embargo, esta condición se refiere a la presión total (base+perturbada), por lo tanto

(Como es habitual, las cantidades perturbadas se pueden linealizar sobre la superficie z=0 ). Utilizando el equilibrio hidrostático , en la forma

Esto se convierte en

Las presiones perturbadas se evalúan en términos de funciones de corriente, utilizando la ecuación de momento horizontal de las ecuaciones de Euler linealizadas para las perturbaciones, con el resultado

Poniendo esta última ecuación y la condición de salto juntas,

Sustituyendo la segunda condición interfacial y utilizando la representación en modo normal, esta relación se convierte en

donde no hay necesidad de etiquetar (solo sus derivados) porque en

Solución

Ahora que se ha establecido el modelo de flujo estratificado, la solución está al alcance de la mano. La ecuación de la función de corriente con las condiciones de contorno tiene la solución

La primera condición interfacial establece que en , lo que obliga a La tercera condición interfacial establece que

Introduciendo la solución en esta ecuación se obtiene la relación

La A se cancela de ambos lados y nos quedamos con

Para comprender en su totalidad las implicaciones de este resultado, es útil considerar el caso de tensión superficial cero.

y claramente

  • Si , y c es real. Esto sucede cuando el líquido para encendedores se encuentra en la parte superior;
  • Si , y c es puramente imaginario. Esto sucede cuando el fluido más pesado se encuentra en la parte superior.

Ahora, cuando el fluido más pesado se asienta en la parte superior, y

donde es el número de Atwood . Al tomar la solución positiva, vemos que la solución tiene la forma

y esto está asociado a la posición de la interfaz η por: Ahora defina

Cuando se permite que las dos capas del fluido tengan una velocidad relativa, la inestabilidad se generaliza a la inestabilidad de Kelvin-Helmholtz-Rayleigh-Taylor, que incluye tanto la inestabilidad de Kelvin-Helmholtz como la inestabilidad de Rayleigh-Taylor como casos especiales. Recientemente se descubrió que las ecuaciones de fluidos que gobiernan la dinámica lineal del sistema admiten una simetría de paridad-tiempo , y la inestabilidad de Kelvin-Helmholtz-Rayleigh-Taylor ocurre cuando y solo cuando la simetría de paridad-tiempo se rompe espontáneamente. [18]

Explicación de la vorticidad

Visualización de una configuración inestable de Rayleigh-Taylor donde el par baroclínico en la interfaz crea vorticidad e induce un campo de velocidad que aumenta el par baroclínico. Aquí ω es la vorticidad, p es la presión, ρ es la densidad, u es la velocidad y g es la gravedad. Las flechas circulares gruesas representan el campo de velocidad creado por el vórtice.

La inestabilidad RT puede verse como el resultado del torque baroclínico creado por la desalineación de los gradientes de presión y densidad en la interfaz perturbada, como se describe en la ecuación de vorticidad no viscosa bidimensional , donde ω es la vorticidad, ρ la densidad y p es la presión. En este caso, el gradiente de presión dominante es hidrostático , resultante de la aceleración.

Cuando se encuentra en la configuración inestable, para un componente armónico particular de la perturbación inicial, el torque en la interfaz crea vorticidad que tenderá a aumentar la desalineación de los vectores de gradiente . Esto a su vez crea vorticidad adicional, lo que conduce a una mayor desalineación. Este concepto se representa en la figura, donde se observa que los dos vórtices contrarrotativos tienen campos de velocidad que se suman en el pico y el valle de la interfaz perturbada. En la configuración estable, la vorticidad, y por lo tanto el campo de velocidad inducido, estarán en una dirección que disminuye la desalineación y, por lo tanto, estabiliza el sistema. [16] [19]

Una explicación mucho más simple de la física básica de la inestabilidad de Rayleigh-Taylor se puede encontrar en la Ref.20.

Comportamiento tardío

El análisis de la sección anterior falla cuando la amplitud de la perturbación es grande. El crecimiento se vuelve entonces no lineal a medida que los picos y las burbujas de la inestabilidad se enredan y se enrollan formando vórtices. Entonces, como en la figura, se requiere una simulación numérica del problema completo para describir el sistema.

Véase también

Notas

  1. ^ Li, Shengtai y Hui Li. "Código AMR paralelo para ecuaciones MHD o HD compresibles". Laboratorio Nacional de Los Álamos . Consultado el 5 de septiembre de 2006 .
  2. ^ abc Sharp, DH (1984). "Una visión general de la inestabilidad de Rayleigh-Taylor". Physica D . 12 (1): 3–18. Bibcode :1984PhyD...12....3S. doi :10.1016/0167-2789(84)90510-4.
  3. ^ abcd Drazin (2002) págs.
  4. ^ David Youngs (ed.). "Inestabilidad y mezcla de Rayleigh-Taylor". Scholarpedia .
  5. ^ "Por qué las bombas nucleares crean nubes de hongos". 20 de noviembre de 2013.
  6. ^ Wang, C.-Y. y Chevalier RA (2000). "Inestabilidades y aglutinamiento en remanentes de supernovas de tipo Ia". The Astrophysical Journal . 549 (2): 1119–1134. arXiv : astro-ph/0005105v1 . Código Bibliográfico :2001ApJ...549.1119W. doi :10.1086/319439. S2CID  15244583.
  7. ^ Hillebrandt, W.; Höflich, P. (1992). "Supernova 1987a en la Gran Nube de Magallanes". En RJ Tayler (ed.). Astrofísica estelar . CRC Press . pp. 249–302. ISBN. 978-0-7503-0200-5.. Véase la página 274.
  8. ^ Chen, HB; Hilko, B.; Panarella, E. (1994). "La inestabilidad de Rayleigh-Taylor en el pinch esférico". Journal of Fusion Energy . 13 (4): 275–280. Bibcode :1994JFuE...13..275C. doi :10.1007/BF02215847. S2CID  122223176.
  9. ^ Betti, R.; Goncharov, VN; McCrory, RL; Verdon, CP (1998). "Tasas de crecimiento de la inestabilidad ablativa de Rayleigh-Taylor en la fusión por confinamiento inercial". Física de plasmas . 5 (5): 1446–1454. Código Bibliográfico :1998PhPl....5.1446B. doi :10.1063/1.872802.
  10. ^ John Pritchett (1971). "EVALUACIÓN DE VARIOS MODELOS TEÓRICOS DE EXPLOSIÓN SUBACUÁTICA" (PDF) . Gobierno de EE. UU., pág. 86. Archivado desde el original (PDF) el 18 de octubre de 2012. Consultado el 9 de octubre de 2012 .
  11. ^ Hester, J. Jeff (2008). "La Nebulosa del Cangrejo: una quimera astrofísica". Revista anual de astronomía y astrofísica . 46 : 127–155. Código Bibliográfico :2008ARA&A..46..127H. doi :10.1146/annurev.astro.45.051806.110608.
  12. ^ Berger, Thomas E.; Slater, Gregory; Hurlburt, Neal; Shine, Richard; et al. (2010). "Dinámica de prominencias quiescentes observadas con el telescopio óptico solar Hinode. I. Columnas de flujo ascendente turbulento". The Astrophysical Journal . 716 (2): 1288–1307. Bibcode :2010ApJ...716.1288B. doi : 10.1088/0004-637X/716/2/1288 .
  13. ^ Chandrasekhar 1981, Cap. X.
  14. ^ Hillier, A.; Berger, Thomas; Isobe, Hiroaki; Shibata, Kazunari (2012). "Simulaciones numéricas de la inestabilidad magnética de Rayleigh-Taylor en el modelo de prominencia de Kippenhahn-Schlüter. I. Formación de flujos ascendentes". The Astrophysical Journal . 716 (2): 120–133. Bibcode :2012ApJ...746..120H. doi : 10.1088/0004-637X/746/2/120 .
  15. ^ Singh, Chamkor; Das, Arup K.; Das, Prasanta K. (2016), "Inestabilidad monomodo de una interfaz ferrofluido-mercurio bajo un campo magnético no uniforme", Physical Review E , 94 (1): 012803, Bibcode :2016PhRvE..94a2803S, doi :10.1103/PhysRevE.94.012803, PMID  27575198
  16. ^ ab Roberts, MS; Jacobs, JW (2015). "Los efectos de las perturbaciones iniciales forzadas de longitud de onda pequeña y ancho de banda finito y la miscibilidad en la inestabilidad turbulenta de Rayleigh Taylor". Journal of Fluid Mechanics . 787 : 50–83. Bibcode :2016JFM...787...50R. doi :10.1017/jfm.2015.599. OSTI  1436483. S2CID  126143908.
  17. ^ desde Drazin (2002) págs. 48–52.
  18. ^ Qin, H.; et al. (2019). "La inestabilidad de Kelvin–Helmholtz es el resultado de la ruptura de la simetría de paridad-tiempo". Física de plasmas . 26 (3): 032102. arXiv : 1810.11460 . Código Bibliográfico :2019PhPl...26c2102Q. doi :10.1063/1.5088498. S2CID  53658729.
  19. ^ Roberts, MS (2012). Experimentos y simulaciones sobre la inestabilidad incompresible de Rayleigh-Taylor con perturbaciones iniciales de longitud de onda pequeña (tesis doctoral). Tesis doctorales de la Universidad de Arizona. Bibcode :2012PhDT.......222R. hdl : 10150/265355 .

20.^ AR Piriz, OD Cortazar, JJ López Cela y NA Tahir, "The Rayleigh-Taylor instability", Am. J. Física. 74 , 1095 (2006)

Referencias

Artículos de investigación originales

Otro

Enlaces externos