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 el enfoque está 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{aligned}{\frac {\partial \mathbf {u} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {u} &=-\nabla p+\mathbf {B} \cdot \nabla \mathbf {B} +\nu \nabla ^{2}\mathbf {u} \\[5pt]{\frac {\partial \mathbf {B} }{\partial 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{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd0712057d9430e9904ae4f147a77be4cc0776f9)
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).
El campo magnético total se puede dividir en dos partes:
(media + fluctuaciones).
Las ecuaciones anteriores en términos de variables de Elsässer (
) están
![{\displaystyle {\frac {\partial {\mathbf {z} ^{\pm }}}{\partial 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 }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/039c7a7e785a06d8416e46cec6e1e216e40139bf)
dónde
. Se producen interacciones no lineales entre las fluctuaciones alfvénicas.
.
Los parámetros adimensionales importantes para MHD son
![{\displaystyle {\begin{array}{lcl}{\text{Reynolds number }}Re&=&UL/\nu \\{\text{Magnetic Reynolds number }}Re_{M}&=&UL/\eta \\{\text{Magnetic Prandtl number }}P_{M}&=&\nu /\eta .\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d7bb6596c194fe3b3df0fa61ab37caf3f959f63)
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, sodio líquido
esta alrededor
. Pero los plasmas tienen grandes
.
El número de Reynolds es la razón del término no lineal
de la ecuación de Navier-Stokes para 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.
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).
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 disminuyendo la cascada de energía, etc. Los primeros modelos de turbulencia MHD asumían la isotropía de turbulencia, mientras que los modelos posteriores han estudiado aspectos anisotrópicos. En las siguientes secciones se resumirán estos modelos. Se pueden encontrar más discusiones sobre la turbulencia MHD en Biskamp,[1] Verma.[2] y Galtier.
Modelos isotrópicos
Iroshnikov[3] y Kraichnan[4] formularon la primera teoría fenomenológica de la turbulencia MHD. Argumentaron que en presencia de un campo magnético medio fuerte,
y
Los paquetes de ondas viajan en direcciones opuestas con la velocidad de fase de
e interactuar débilmente. La escala de tiempo relevante es el tiempo de Alfven
. Como resultado, los espectros de energía son
![{\displaystyle E^{u}(k)\approx E^{b}(k)\approx A(\Pi V_{A})^{1/2}k^{-3/2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a3ad439925e8b1d5720250528853e9995f7f422)
dónde
es la tasa en cascada de energía.
Más tarde Dobrowolny y otros derivaron las siguientes fórmulas generalizadas para las tasas en cascada.[5]
![{\displaystyle \Pi ^{+}\approx \Pi ^{-}\approx \tau _{k}^{\pm }E^{+}(k)E^{-}(k)k^{4}\approx E^{+}(k)E^{-}(k)k^{3}/B_{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6455c7b14b1f7e74d9d0ce94e79996d25c29a25d)
dónde
son las escalas de tiempo de interacción de
variables.
La fenomenología de Iroshnikov y Kraichnan sigue una vez que elegimos
.
Marsch[6] 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 tipo Kolmogorov para las variables de Elsasser:
![{\displaystyle E^{\pm }(k)=K^{\pm }(\Pi ^{\pm })^{4/3}(\Pi ^{\mp })^{-2/3}k^{-5/3}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d9bf5834212ffa95d858f2f015ae2c065a342c5)
dónde
y
son las tasas de cascada de energía de
y
respectivamente, y
son constantes.
Matthaeus y Zhou[7] 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[8] utilizando un análisis de grupo de renormalización mostrando que las fluctuaciones alfvénicas se ven afectadas por el "campo magnético medio local" dependiente de la escala. El campo magnético medio local se escala como
, cuya sustitución en la ecuación de Dobrowolny produce el espectro de energía de Kolmogorov para la turbulencia MHD.
También se han realizado análisis del grupo de renormalización para calcular la viscosidad y resistividad renormalizadas. Se demostró que estas magnitudes difusivas se escalan como
que vuelve a ceder
espectros de energía consistentes con el modelo de Kolmogorov para turbulencia MHD. El cálculo del grupo de renormalización anterior se ha realizado para helicidad cruzada cero y distinta de cero.
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.[9]
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 y col.[10] mostró usando ecuaciones cinéticas que
![{\displaystyle E(k)\sim (\Pi B_{0})^{1/2}k_{\parallel }^{1/2}k_{\perp }^{-2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64f59cc6855b48e088f3c6f9d2418cc716b22001)
dónde
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 .
Bajo el fuerte límite de turbulencia,
, Goldereich y Sridhar[11] argumentan que
("estado crítico equilibrado") lo que implica que
![{\displaystyle {\begin{aligned}E(k)&\propto k_{\perp }^{-5/3};\\[5pt]k_{\parallel }&\propto k_{\perp }^{2/3}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9e05b5549d647c912750ee4624747fd988740134)
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 de
en comparación con
, favoreciendo así la fenomenología similar a Kolmogorov para la turbulencia MHD.[12][13] Las fluctuaciones de densidad de electrones interplanetarios e interestelares también proporcionan una ventana para investigar la turbulencia MHD.
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.[14] 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.
Flujos de energía
pueden ser cantidades más confiables para validar los modelos de turbulencia MHD. Cuándo
(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 están más de acuerdo con el modelo de Kolmogorov en comparación con el modelo de Kraichnan e Iroshnikov.[15]
Los aspectos anisotrópicos de la turbulencia MHD también se han estudiado mediante simulaciones numéricas. Las predicciones de Goldreich y Sridhar[11] (
) se han verificado en muchas simulaciones.
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.[2] 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.
Véase también
Referencias
- ↑ D. Biskamp (2003), Magnetohydrodynamical Turbulence, (Cambridge University Press, Cambridge.)
- ↑ a b Verma, Mahendra K. (2004). «Statistical theory of magnetohydrodynamic turbulence: recent results». Physics Reports 401 (5–6): 229-380. ISSN 0370-1573. arXiv:nlin/0404043. doi:10.1016/j.physrep.2004.07.007.
- ↑ P. S. Iroshnikov (1964), Turbulence of a Conducting Fluid in a Strong Magnetic Field, Soviet Astronomy, 7, 566.
- ↑ Kraichnan, Robert H. (1965). «Inertial-Range Spectrum of Hydromagnetic Turbulence». Physics of Fluids (AIP Publishing) 8 (7): 1385. ISSN 0031-9171. doi:10.1063/1.1761412.
- ↑ Dobrowolny, M.; Mangeney, A.; Veltri, P. (14 de julio de 1980). «Fully Developed Anisotropic Hydromagnetic Turbulence in Interplanetary Space». Physical Review Letters (American Physical Society (APS)) 45 (2): 144-147. ISSN 0031-9007. doi:10.1103/physrevlett.45.144.
- ↑ E. Marsch (1990), Turbulence in the solar wind, in: G. Klare (Ed.), Reviews in Modern Astronomy, Springer, Berlin, p. 43.
- ↑ Matthaeus, William H.; Zhou, Ye (1989). «Extended inertial range phenomenology of magnetohydrodynamic turbulence». Physics of Fluids B: Plasma Physics (AIP Publishing) 1 (9): 1929-1931. ISSN 0899-8221. doi:10.1063/1.859110.
- ↑ Verma, Mahendra K. (1999). «Mean magnetic field renormalization and Kolmogorov's energy spectrum in magnetohydrodynamic turbulence». Physics of Plasmas (AIP Publishing) 6 (5): 1455-1460. ISSN 1070-664X. doi:10.1063/1.873397.
- ↑ Shebalin, John V.; Matthaeus, William H.; Montgomery, David (1983). «Anisotropy in MHD turbulence due to a mean magnetic field». Journal of Plasma Physics (Cambridge University Press (CUP)) 29 (3): 525-547. ISSN 0022-3778. doi:10.1017/s0022377800000933.
- ↑ Galtier, S.; Nazarenko, S. V.; Newell, A. C.; Pouquet, A. (2000). «A weak turbulence theory for incompressible magnetohydrodynamics». Journal of Plasma Physics (Cambridge University Press (CUP)) 63 (5): 447-488. ISSN 0022-3778. doi:10.1017/s0022377899008284.
- ↑ a b Goldreich, P.; Sridhar, S. (1995). «Toward a theory of interstellar turbulence. 2: Strong alfvenic turbulence». The Astrophysical Journal (IOP Publishing) 438: 763. ISSN 0004-637X. doi:10.1086/175121.
- ↑ Matthaeus, William H.; Goldstein, Melvyn L. (1982). «Measurement of the rugged invariants of magnetohydrodynamic turbulence in the solar wind». Journal of Geophysical Research (American Geophysical Union (AGU)) 87 (A8): 6011. ISSN 0148-0227. doi:10.1029/ja087ia08p06011.
- ↑ D. A. Roberts, M. L. Goldstein (1991), Turbulence and waves in the solar wind, Rev. Geophys., 29, 932.
- ↑ Müller, Wolf-Christian; Biskamp, Dieter (17 de enero de 2000). «Scaling Properties of Three-Dimensional Magnetohydrodynamic Turbulence». Physical Review Letters (American Physical Society (APS)) 84 (3): 475-478. ISSN 0031-9007. PMID 11015942. arXiv:physics/9906003. doi:10.1103/physrevlett.84.475.
- ↑ Verma, M. K.; Roberts, D. A.; Goldstein, M. L.; Ghosh, S.; Stribling, W. T. (1 de octubre de 1996). «A numerical study of the nonlinear cascade of energy in magnetohydrodynamic turbulence». Journal of Geophysical Research: Space Physics (American Geophysical Union (AGU)) 101 (A10): 21619-21625. ISSN 0148-0227. doi:10.1029/96ja01773.
Enlaces externos