Turbulencia magnetohidrodinámica - Magnetohydrodynamic turbulence

La turbulencia magnetohidrodinámica se refiere a los regímenes caóticos de flujo de magnetofluidos a un alto número de Reynolds . La magnetohidrodinámica (MHD) se ocupa de lo que es un fluido cuasi neutro con una conductividad muy alta . La aproximación fluida implica que la atención se centra en escalas macro de longitud y tiempo que son mucho más grandes que la duración de la colisión y el tiempo de colisión, respectivamente.

Ecuaciones MHD incompresibles

Las ecuaciones MHD incompresibles son

{\ Displaystyle {\ begin {alineado} {\ frac {\ parcial \ mathbf {u}} {\ parcial t}} + \ mathbf {u} \ cdot \ nabla \ mathbf {u} & = - \ nabla p + \ mathbf {B} \ cdot \ nabla \ mathbf {B} + \ nu \ nabla ^ {2} \ mathbf {u} \\ [5pt] {\ frac {\ parcial \ mathbf {B}} {\ parcial t}} + \ mathbf {u} \ cdot \ nabla \ mathbf {B} & = \ mathbf {B} \ cdot \ nabla \ mathbf {u} + \ eta \ nabla ^ {2} \ mathbf {B} \\ [5pt] \ nabla \ cdot \ mathbf {u} & = 0 \\ [5pt] \ nabla \ cdot \ mathbf {B} & = 0. \ end {alineado}}}

donde u , B , p representan los campos de velocidad, magnético y de presión total (térmico + magnético), y representan la viscosidad cinemática y la difusividad magnética . La tercera ecuación es la condición de incompresibilidad . En la ecuación anterior, el campo magnético está en unidades Alfvén (lo mismo que las unidades de velocidad). ${\ Displaystyle \ nu}$ ${\ Displaystyle \ eta}$

El campo magnético total se puede dividir en dos partes: (media + fluctuaciones). ${\ Displaystyle \ mathbf {B} = \ mathbf {B_ {0}} + \ mathbf {b}}$

Las ecuaciones anteriores en términos de variables de Elsässer ( ) son ${\ Displaystyle \ mathbf {z} ^ {\ pm} = \ mathbf {u} \ pm \ mathbf {b}}$

{\ displaystyle {\ frac {\ parcial {\ mathbf {z} ^ {\ pm}}} {\ parcial t}} \ mp \ left (\ mathbf {B} _ {0} \ cdot {\ mathbf {\ nabla }} \ right) {\ mathbf {z} ^ {\ pm}} + \ left ({\ mathbf {z} ^ {\ mp}} \ cdot {\ mathbf {\ nabla}} \ right) {\ mathbf { z} ^ {\ pm}} = - {\ mathbf {\ nabla}} p + \ nu _ {+} \ nabla ^ {2} \ mathbf {z} ^ {\ pm} + \ nu _ {-} \ nabla ^ {2} \ mathbf {z} ^ {\ mp}}

donde . Se producen interacciones no lineales entre las fluctuaciones alfvénicas . ${\ Displaystyle \ nu _ {\ pm} = {\ frac {1} {2}} (\ nu \ pm \ eta)}$ ${\ Displaystyle z ^ {\ mp}}$

Los parámetros adimensionales importantes para MHD son

{\ displaystyle {\ begin {array} {lcl} {\ text {número de Reynolds}} Re & = & UL / \ nu \\ {\ text {número de Reynolds magnético}} Re_ {M} & = & UL / \ eta \\ { \ text {Número de Prandtl magnético}} P_ {M} & = & \ nu / \ eta. \ end {array}}}

El número de Prandtl magnético es una propiedad importante del fluido. Los metales líquidos tienen pequeños números de Prandtl magnéticos, por ejemplo, el sodio líquido está alrededor . Pero los plasmas tienen grandes . ${\ Displaystyle P_ {M}}$ ${\ displaystyle 10 ^ {- 5}}$ ${\ Displaystyle P_ {M}}$

El número de Reynolds es la relación entre el término no lineal de la ecuación de Navier-Stokes y el término viscoso. Mientras que el número de Reynolds magnético es la relación entre el término no lineal y el término difusivo de la ecuación de inducción. ${\ Displaystyle \ mathbf {u} \ cdot \ nabla \ mathbf {u}}$

En muchas situaciones prácticas, el número de Reynolds del flujo es bastante grande. Para tales flujos, típicamente, la velocidad y los campos magnéticos son aleatorios. Tales flujos están llamados a exhibir turbulencias MHD. Tenga en cuenta que no es necesario que sea grande para la turbulencia del MHD. juega un papel importante en el problema de la dínamo (generación de campo magnético). ${\ displaystyle Re}$ ${\ displaystyle Re_ {M}}$ ${\ displaystyle Re_ {M}}$

El campo magnético medio juega un papel importante en la turbulencia MHD, por ejemplo, puede hacer que la turbulencia sea anisotrópica; suprimir la turbulencia al disminuir la cascada de energía, etc. Los primeros modelos de turbulencia MHD asumían isotropía de turbulencia, mientras que los modelos posteriores han estudiado aspectos anisotrópicos. En las siguientes discusiones se resumirán estos modelos. Se pueden encontrar más discusiones sobre la turbulencia MHD en Biskamp, Verma. y Galtier.

Modelos isotrópicos

Iroshnikov y Kraichnan formularon la primera teoría fenomenológica de la turbulencia MHD. Argumentaron que, en presencia de un fuerte campo magnético media, y paquetes de ondas viajan en direcciones opuestas con la velocidad de fase de e interactúan débilmente. La escala de tiempo relevante es la hora de Alfven . Como resultado, los espectros de energía son ${\ Displaystyle z ^ {+}}$ ${\ Displaystyle z ^ {-}}$ ${\ Displaystyle B_ {0}}$ ${\ Displaystyle (B_ {0} k) ^ {- 1}}$

{\ Displaystyle E ^ {u} (k) \ approx E ^ {b} (k) \ approx A (\ Pi V_ {A}) ^ {1/2} k ^ {- 3/2}.}

donde es la tasa de cascada de energía. ${\ Displaystyle \ Pi}$

Posteriormente Dobrowolny et al. derivó las siguientes fórmulas generalizadas para las tasas en cascada de las variables: ${\ Displaystyle z ^ {\ pm}}$

{\ Displaystyle \ Pi ^ {+} \ approx \ Pi ^ {-} \ approx \ tau _ {k} ^ {\ pm} E ^ {+} (k) E ^ {-} (k) k ^ {4 } \ approx E ^ {+} (k) E ^ {-} (k) k ^ {3} / B_ {0}}

donde están las escalas de tiempo de interacción de las variables. ${\ Displaystyle \ tau ^ {\ pm}}$ ${\ Displaystyle z ^ {\ pm}}$

La fenomenología de Iroshnikov y Kraichnan sigue una vez que elegimos . ${\ Displaystyle \ tau ^ {\ pm} \ approx 1 / (kV_ {A})}$

Marsch eligió la escala de tiempo no lineal como la escala de tiempo de interacción para los remolinos y el espectro energético derivado de Kolmogorov para las variables de Elsasser: ${\ Displaystyle T_ {NL} ^ {\ pm} \ approx (kz_ {k} ^ {\ mp}) ^ {- 1}}$

{\ Displaystyle E ^ {\ pm} (k) = K ^ {\ pm} (\ Pi ^ {\ pm}) ^ {4/3} (\ Pi ^ {\ mp}) ^ {- 2/3} k ^ {- 5/3}}

donde y son las tasas de cascada de energía de y respectivamente, y son constantes. ${\ Displaystyle \ Pi ^ {+}}$ ${\ Displaystyle \ Pi ^ {-}}$ ${\ Displaystyle z ^ {+}}$ ${\ Displaystyle z ^ {-}}$ ${\ Displaystyle K ^ {\ pm}}$

Matthaeus y Zhou intentaron combinar las dos escalas de tiempo anteriores postulando que el tiempo de interacción es la media armónica del tiempo de Alfven y el tiempo no lineal.

La principal diferencia entre las dos fenomenologías en competencia (−3/2 y −5/3) son las escalas de tiempo elegidas para el tiempo de interacción. La principal suposición subyacente es que la fenomenología de Iroshnikov y Kraichnan debería funcionar para un campo magnético medio fuerte, mientras que la fenomenología de Marsh debería funcionar cuando las fluctuaciones dominan el campo magnético medio (turbulencia fuerte).

Sin embargo, como discutiremos a continuación, las observaciones del viento solar y las simulaciones numéricas tienden a favorecer el espectro de energía de -5/3 incluso cuando el campo magnético medio es más fuerte en comparación con las fluctuaciones. Este problema fue resuelto por Verma utilizando un análisis de grupo de renormalización al mostrar que las fluctuaciones de Alfvénic se ven afectadas por el "campo magnético medio local" dependiente de la escala. El campo magnético medio local se escala como , la sustitución de la cual en la ecuación de Dobrowolny produce el espectro de energía de Kolmogorov para la turbulencia MHD. ${\ displaystyle k ^ {- 1/3}}$

También se han realizado análisis del grupo de renormalización para calcular la viscosidad y resistividad renormalizadas. Se demostró que estas cantidades difusivas se escalan a medida que eso produce nuevamente espectros de energía consistentes con el modelo similar a Kolmogorov para turbulencias MHD. El cálculo del grupo de renormalización anterior se ha realizado para helicidad cruzada cero y distinta de cero. ${\ displaystyle k ^ {- 4/3}}$ ${\ displaystyle k ^ {- 5/3}}$

Las fenomenologías anteriores suponen una turbulencia isotrópica que no es el caso en presencia de un campo magnético medio. El campo magnético medio normalmente suprime la cascada de energía a lo largo de la dirección del campo magnético medio.

Modelos anisotrópicos

El campo magnético medio hace que la turbulencia sea anisotrópica. Este aspecto ha sido estudiado en las últimas dos décadas. En el límite , Galtier et al. mostró usando ecuaciones cinéticas que ${\ Displaystyle \ delta z ^ {\ pm} \ ll B_ {0}}$

{\ Displaystyle E (k) \ sim (\ Pi B_ {0}) ^ {1/2} k _ {\ paralelo} ^ {1/2} k _ {\ perp} ^ {- 2}}

donde y son componentes del número de onda paralelos y perpendiculares al campo magnético medio. El límite anterior se denomina límite de turbulencia débil . ${\ Displaystyle k _ {\ paralelo}}$ ${\ Displaystyle k _ {\ perp}}$

Bajo el fuerte límite de turbulencia , Goldereich y Sridhar argumentan que ("estado crítico equilibrado"), lo que implica que ${\ Displaystyle \ delta z ^ {\ pm} \ sim B_ {0}}$ ${\ Displaystyle k _ {\ perp} z_ {k _ {\ perp}} \ sim k _ {\ paralelo} B_ {0}}$

{\ displaystyle {\ begin {alineado} E (k) & \ propto k _ {\ perp} ^ {- 5/3}; \\ [5pt] k _ {\ paralelo} & \ propto k _ {\ perp} ^ {2 / 3} \ end {alineado}}}

La fenomenología de turbulencia anisotrópica anterior se ha extendido para MHD de helicidad cruzada grande.

Observaciones del viento solar

El plasma del viento solar está en estado turbulento. Los investigadores han calculado los espectros de energía del plasma del viento solar a partir de los datos recopilados de la nave espacial. Los espectros de energía cinética y magnética, así como están más cerca en comparación , favoreciendo así la fenomenología similar a Kolmogorov para la turbulencia MHD. Las fluctuaciones de densidad de electrones interplanetarios e interestelares también proporcionan una ventana para investigar la turbulencia MHD. ${\ Displaystyle E ^ {\ pm}}$ ${\ displaystyle k ^ {- 5/3}}$ ${\ Displaystyle k ^ {- 3/2}}$

Simulaciones numéricas

Los modelos teóricos discutidos anteriormente se prueban utilizando la simulación numérica directa (DNS) de alta resolución. El número de simulaciones recientes informa que los índices espectrales están más cerca de 5/3. Hay otros que reportan índices espectrales cercanos a 3/2. El régimen de la ley del poder es típicamente de menos de una década. Dado que 5/3 y 3/2 son bastante cercanos numéricamente, es bastante difícil determinar la validez de los modelos de turbulencia MHD a partir de los espectros de energía.

Los flujos de energía pueden ser cantidades más confiables para validar los modelos de turbulencia MHD. Cuando (fluido de alta helicidad cruzada o MHD desequilibrado) las predicciones de flujo de energía del modelo de Kraichnan e Iroshnikov son muy diferentes de las del modelo similar a Kolmogorov. Se ha demostrado usando DNS que los flujos calculados a partir de las simulaciones numéricas concuerdan mejor con el modelo similar a Kolmogorov en comparación con el modelo de Kraichnan e Iroshnikov. ${\ Displaystyle \ Pi ^ {\ pm}}$ ${\ Displaystyle E ^ {+} (k) \ gg E ^ {-} (k)}$ ${\ Displaystyle \ Pi ^ {\ pm}}$

Los aspectos anisotrópicos de la turbulencia MHD también se han estudiado mediante simulaciones numéricas. Las predicciones de Goldreich y Sridhar ( ) se han verificado en muchas simulaciones. ${\ Displaystyle k_ {||} \ sim k _ {\ perp} ^ {2/3}}$

Transferencia de energía

La transferencia de energía entre varias escalas entre la velocidad y el campo magnético es un problema importante en la turbulencia MHD. Estas cantidades se han calculado tanto teórica como numéricamente. Estos cálculos muestran una transferencia de energía significativa desde el campo de velocidad a gran escala al campo magnético a gran escala. Además, la cascada de energía magnética es típicamente hacia adelante. Estos resultados tienen una relación crítica con el problema de la dínamo.

Hay muchos desafíos abiertos en este campo que, con suerte, se resolverán en un futuro próximo con la ayuda de simulaciones numéricas, modelos teóricos, experimentos y observaciones (por ejemplo, viento solar).

Languages

In other projects