Question 1

¿Qué es la inestabilidad de Kelvin-Helmholtz en un plasma?

Accepted Answer

La inestabilidad de Kelvin-Helmholtz (KH) en un plasma ocurre en la interfaz entre dos capas de plasma que se mueven a distintas velocidades. La cizalladura de velocidad en el límite amplifica pequeñas perturbaciones: la presión cae donde el flujo es más rápido (efecto Bernoulli), lo que tira de la interfaz hacia dentro, generando vórtices. En un plasma magnetizado, el campo magnético puede estabilizar o desestabilizar según su orientación respecto al flujo de cizalladura.

Question 2

¿Qué es la fórmula de la tasa de crecimiento γ = k·ΔV/2?

Accepted Answer

Para dos capas de fluido incompresible de igual densidad con una diferencia de velocidad ΔV = V₁ − V₂, la tasa de crecimiento KH es γ = k·ΔV/2, donde k = 2π/λ es el número de onda de la perturbación superficial. Esto significa que las longitudes de onda más cortas (k mayor) crecen más rápido en el límite de densidad igual. Cuando las densidades difieren, la tasa de crecimiento se convierte en γ = k·ΔV·√(ρ₁ρ₂)/(ρ₁+ρ₂), reduciendo el crecimiento para grandes contrastes de densidad.

Question 3

¿Dónde ocurre la inestabilidad KH del plasma en la naturaleza?

Accepted Answer

La inestabilidad KH del plasma es omnipresente en los plasmas espaciales. En la magnetopausa de la Tierra, el viento solar que se desliza junto a la magnetosfera impulsa vórtices KH que transportan plasma solar hacia el interior. Los cinturones de nubes de Júpiter muestran ondulaciones KH en los límites de las corrientes en chorro. En la corona solar, se han captado imágenes directas de vórtices KH en plumas y lluvia coronal. Las nebulosas de viento de púlsar y los chorros astrofísicos muestran patrones KH característicos en sus límites.

Question 4

¿Cómo afecta la relación de densidad a la inestabilidad KH?

Accepted Answer

La relación de densidad ρ₂/ρ₁ modifica la tasa de crecimiento a través del factor √(ρ₁ρ₂)/(ρ₁+ρ₂). Con densidad igual (relación = 1) este factor equivale a 1/2, dando la tasa de crecimiento máxima k·ΔV/2. A medida que aumenta el contraste de densidad (ρ₂/ρ₁ grande), la tasa de crecimiento efectiva disminuye y la capa pesada resiste la deflexión. Una relación de densidad muy grande (ρ₂/ρ₁ ≫ 1) suprime la inestabilidad: la capa más pesada apenas se mueve mientras la capa ligera se deforma a su alrededor.

Question 5

¿Qué es el enrollamiento de vórtices y cómo se desarrolla?

Accepted Answer

El enrollamiento de vórtices es la saturación no lineal de la inestabilidad KH. Inicialmente, las perturbaciones sinusoidales crecen exponencialmente (fase lineal). A medida que las amplitudes se hacen comparables a la longitud de onda, la vorticidad se concentra en núcleos discretos: la interfaz se enrolla formando vórtices en forma de ojo de gato. Los vórtices adyacentes acaban fusionándose (emparejamiento de vórtices) para formar estructuras más grandes, una cascada de energía inversa que impulsa la mezcla turbulenta a través de la capa de cizalladura.

Question 6

¿Cómo afecta el modo de perturbación k a la simulación?

Accepted Answer

El modo de perturbación entero k establece cuántas longitudes de onda sinusoidales completas caben a lo largo del dominio de la simulación. El modo k=1 muestra un único enrollamiento de vórtice grande; k=4 muestra cuatro vórtices simultáneos. Los modos más altos crecen más rápido (γ ∝ k) en la fase lineal, pero también se saturan a una amplitud menor en relación con el tamaño del dominio. El enrollamiento visualmente más espectacular suele producirse en modos intermedios k=2–4.

Question 7

¿Qué estabiliza la inestabilidad de Kelvin-Helmholtz en un plasma magnetizado?

Accepted Answer

Una componente del campo magnético paralela al flujo proporciona tensión magnética que resiste la flexión de las líneas de campo durante la formación de vórtices. La inestabilidad se estabiliza cuando la velocidad de Alfvén V_A = B/√(μ₀ρ) supera ΔV/2, es decir, cuando la energía magnética supera la energía cinética de la cizalladura. Un campo magnético perpendicular no tiene efecto estabilizador porque las líneas de campo pueden deslizarse entre sí, pero añade efectos de compresibilidad.

Question 8

¿Cuál es la tasa de crecimiento de la capa de mezcla en la fase no lineal?

Accepted Answer

Tras el crecimiento exponencial inicial, el ancho de la capa de mezcla δ crece linealmente con el tiempo como δ = α·ΔV·t, donde α ≈ 0,06–0,1 es el coeficiente de mezcla turbulenta (constante de Brown-Roshko). Esta fase de crecimiento lineal persiste hasta que la capa de cizalladura llena el dominio. El espesor de momento de la capa de mezcla θ crece como dθ/dt ≈ 0,017·ΔV, un resultado confirmado tanto por experimentos como por simulaciones numéricas directas.

Question 9

¿Cómo es relevante la inestabilidad KH del plasma para la energía de fusión?

Accepted Answer

En los plasmas de borde de los tokamaks, la inestabilidad KH impulsa el transporte turbulento a través de la separatriz, mezclando el plasma caliente del núcleo con el plasma más frío del borde. Los modos localizados de borde (ELM) implican inestabilidades tipo KH en el límite del pedestal. En la fusión por confinamiento inercial (ICF), las inestabilidades KH en la interfaz ablador-combustible degradan la simetría de la implosión y mezclan material ablador frío en el combustible caliente, apagando la ignición. Controlar la inestabilidad KH es, por tanto, central para lograr la energía de fusión.

Question 10

¿Qué método numérico utiliza esta simulación?

Accepted Answer

Esta simulación utiliza un método de lámina de vórtices 2D en un dominio periódico. La interfaz entre las dos capas de plasma se representa mediante N puntos marcadores, cada uno con una circulación proporcional al salto de velocidad local ΔV. Sus posiciones evolucionan mediante la integral de Birkhoff-Rott, discretizada usando la ley de Biot-Savart con un parámetro de desingularización (regularización de vórtice tipo blob). Un integrador temporal de Runge-Kutta de cuarto orden avanza la lámina, capturando el enrollamiento hasta la fase de saturación.

⚡ Inestabilidad de Plasma — Kelvin-Helmholtz

Acerca de la simulación

¿Qué encontraste?