INTRODUCCIÓN
La característica principal de ondas de Alfvén de corte en la magnetohidrodinámica ideal, cuando se considera la ley de Ohm simple, es que se propagan con una velocidad de grupo paralela al campo de inducción magnética de fondo en el sistema de referencia en el que plasma está en reposo. Además, las ondas de Alfvén son perturbaciones incompresibles, existe una relación entre las perturbaciones en velocidad y en campo de inducción magnética, y la presión total (plasma más magnética) es constante [1]. Las mismas han sido detectadas desde el comienzo de la era espacial, en una amplia variedad de plasmas como el viento solar [2] y en la magnetósfera de la Tierra [3,4]. En algunos de estos sistemas la densidad y el campo de inducción magnética son tales que los términos de gradiente de presión electrónica y de Hall en la ley de Ohm no pueden despreciarse [5,6].
Las ondas de Alfvén con término Hall han sido estudiadas por otros autores como Mattei [7], Woodward y McKenzie [8,9], Pokhotelov et al. [10], entre otros. Ellos linealizan el sistema de ecuaciones de la magnetohidrodinámica o imponen dependencias particulares para las perturbaciones.
En lugar de linealizar las ecuaciones magnetohidrodinámicas y buscar ondas monocromáticas, Sallago y Platzeck impusieron las condiciones características de ondas de Alfvén tanto para el caso con ley de Ohm simple en plasmas no uniformes[11] como en el caso con término de Hall y gradiente de presión electrónica en plasmas uniformes [12]. Este último trabajo en adelante será mencionado como (SP04). Además Sallago y Platzeck mostraron que las ondas de Alfvén pueden propagarse, cuando se tienen en cuenta la presión electrónica y el término de Hall, si se cumple una condición sobre la dependencia espacial de la densidad de corriente. Esta condición, llamada “condición de polarización”, relaciona la densidad de corriente y su rotor. En el caso linealizado, con pequeñas perturbaciones y ondas monocromáticas, esta condición significa que la perturbación en campo de inducción magnética debe estar polarizada circularmente [7,10], por este motivo no es posible imponer la condición de adiabaticidad ya que el sistema resultaría sobredeterminado.
Por otra parte, cuando se estudian plasmas parcialmente ionizados existe una cantidad relativa de neutros que se pondera de diversas maneras mediante las siguientes relaciones:
la fracción de ionización [13], fracción del plasma sin ionizar [1]
o mediante el grado de ionización [14]
donde n i y n n son la densidad número de los iones y de los neutros respectivamente. Distintos autores denominan de manera distinta al plasma. Algunos consideran dos fluidos, uno cargado y otro neutro. Estas consideraciones deberán ajustarse dependiendo del grado de acoplamiento de los iones y los neutros. Mayores detalles pueden verse el trabajo de Khomenko et al. (2012) [15].
En relación con los plasmas parcialmente ionizados, pueden encontrarse tanto en regiones del Sol quieto como del Sol activo, así como en lo ionosfera terrestre. Vinculado con las zonas solares activas, el estudio de los fenómenos que soportan los tubos de flujo magnético son de interés para el análisis del comportamiento en las estructuras solares, tanto en el análisis de las manchas solares como en los loops coronales [16] y que en su mayoría han sido realizados considerando uniformes las variables físicas no perturbadas. Cuando se admiten variaciones espaciales, el término correspondiente a la derivada de la cantidad de estado de equilibrio multiplicada por la perturbación es despreciado por considerarlo de segundo orden [17]. La importancia relativa de los distintos términos de la ley de Ohm en los ejemplos de PIP mencionados pueden verse en los trabajos de Leake et al. [18] y de Ballester et al. [14]. Estos valores son analizados a partir de modelos como VALC [19] y C7 [20], entre otros.
En el presente trabajo se tiene un gas compuesto por la suma del gas cargado (eléctricamente cuasi neutro) y el gas neutro, para el que se considera que alcanzó la temperatura de equilibrio una vez completada su capacidad de ionización. En la Sec. 2 se demuestra que las ondas de Alfvén de gran amplitud pueden propagarse, cuando se tienen en cuenta la presión electrónica, los términos de Hall y ambipolar en la ley de Ohm. En la Sec. 3 se muestra que esta solución tiende a la encontrada por Sallago y Platzeck [12] si el término ambipolar se desvanece y que la solución linealizada es similar a la solución linealizada de De Pontieu et al. (2001) [21] (en adelante, DP01), si en esta última se considera conductividad infinita.
MÉTODOS
Para poder discutir la solución del sistema de ecuaciones de los PIP en la aproximación magnetohidrodinámica, primero se recordarán algunas ideas.
Ley de Ohm generalizada
Se llega a la expresión de la ley de Ohm generalizada desde el modelo de plasma compuesto por iones, electrones y neutros y, considerando que el plasma sea eléctricamente neutro, resulta [1]:
donde e es la carga del protón, n i el número de protones por unidad de volumen, p e la presión electrónica y F es igual a una expresión en función de la frecuencia de colisión ión-neutro y de la fracción de ionización [13]
Ecuaciones de la PIP-MHD
Se suele designar como PIP-MHD a la magnetohidrodinámica cuando en la ley de Ohm se tienen en cuenta la fracción de ionización y los términos que hemos discutido anteriormente y cuando a pesar de la presencia significativa de neutros los componentes del plasma están acoplados [22].
Las ecuaciones de continuidad, movimiento y la ecuación 0, son iguales a las de la magnetohidrodinámica ideal. La ecuación que resulta modificada es la ecuación de inducción, ya que se obtiene de reemplazar la ley de Ohm en la ley Faraday.
Para un plasma en las condiciones precedentes las ecuaciones son:
Esta última ecuación puede reescribirse en función de la fracción de ionización (5),
donde ε = mi/e.
Ondas de Alfvén
En esta subsección se buscan soluciones de las ecuaciones de la PIP-MHD, no linealizadas, que puedan ser identificadas como ondas de Alfvén. Esta identificación tendrá en cuenta que dichas soluciones tengan las siguientes propiedades: las perturbaciones se propagan con una velocidad de grupo que, en el sistema de referencia en que el plasma se encuentra en reposo, sea paralela al campo magnético de fondo, las perturbaciones de los campos de velocidad y de campo magnético están relacionados, la perturbación es incompresible y que existe una magnitud que permanece constante. Además, se impondrá que las soluciones tiendan, en el límite cuando el término de ambipolar es despreciable, a las ondas de Alfvén con término de Hall del tipo correspondiente. Por otro lado, para que las ondas de Alfvén con término de Hall puedan satisfacer todas las ecuaciones, debe imponerse una condición adicional de “polarización” (ver SP04), por lo que también deberán estarlo las ondas en PIP-MHD.
Se supone que en un plasma con campos de fondo uniformes , se propagan perturbaciones incompresibles, con dependencia espacio temporal
donde
es la velocidad de grupo y a∗ es una constante que habrá que determinar para mostrar la influencia de los términos de la ley de Ohm.
Con estas suposiciones, de la ecuación de continuidad (6) resulta nula la perturbación en densidades:
Debido a ésto, la ecuación de inducción (10) se simplifica. Como ⃗∇ρ = 0, por ser ρ 0 uniforme y ρ 1 nulo, resulta que
con lo que puede verse que el término gradiente de presión electrónica no influye en las soluciones buscadas.
Ahora se deben encontrar las perturbaciones y p 1 que satisfagan a las ecuaciones.
Debido a que , y proponiendo la siguiente relación entre las perturbaciones en velocidad y campo magnético
la ecuación de inducción toma la siguiente forma:
que es
y resulta
Como puede verse, cuando no hay perturbación y , cuando la fracción de ionización tiende a la unidad (o F tiende a cero), resulta que , donde corresponde a la solución totalmente ionizada con término de Hall (ver SP04). Tomando divergencia de (15) puede verse que debido a la condición de incompresibilidad, se satisface también si , esto implica una condición entre .
Finalmente, reemplazando la relación entre las perturbaciones en velocidad y campo magnético y el valor de la velocidad de grupo en la ecuación de movimiento, se obtiene:
El miembro de la izquierda del igual puede reescribirse utilizando identidades vectoriales, con lo que la ecuación de movimiento (20) resulta así
Tomando el rotor de esta expresión (21), se tiene la siguiente ecuación
En el caso totalmente ionizado con término de Hall, el rotor de la densidad de corriente y la densidad de corriente estaban relacionadas, , resultando en una condición de polarización (SP04). Para el presente caso proponiendo que , reemplazando en la ecuación (22), se obtiene
que implica que la única solución posible es la siguiente
y, en consecuencia se tiene
que es una “condición de polarización”.
Luego de reemplazarla en la ecuación (24), se obtiene la siguiente relación entre las constantes a∗, b∗ y C
o lo que es lo mismo, se obtienen los valores de a∗
Por lo tanto, la velocidad de grupo resulta dependiente del parámetro de proporcionalidad entre la densidad de corriente y su rotor, así como del factor de ionización
Retornando a la ecuación de movimiento, reemplazando la condición de polarización y el valor de a∗, resulta la constancia en la región perturbada de una cantidad : la “presión total generalizada modificada”
La “presión total generalizada modificada” tiende a la “presión total generalizada” del caso totalmente ionizado con término de Hall cuando x → 1 ya que a∗ → a .
Nótese además que a partir de la condición de “polarización” se obtiene
Una relación entre la vorticidad y la densidad de corriente similar a ésta también existe para las ondas de Alfvén en campos uniformes cuando se tiene en cuenta el término de Hall (ver SP04).
Además, como la presión del plasma debe ser siempre positiva, la amplitud de las perturbaciones no es arbitraria. Se debe remarcar que hay un flujo de calor que proviene del gradiente de presión electrónica, como ya sucedía en el caso totalmente ionizado con término de Hall.
Comparación con el resultado de DP01
Con el objetivo de realizar la comparación con los resultados de DP01 [21], primero se repasará el significado de la condición de “polarización” (25) para la densidad de corriente en el límite linealizado cuando se consideran ondas planas. Si la perturbación del campo magnético se puede expresar de la forma
la densidad de corriente cumplirá la condición solicitada (25) si la amplitud b1 satisface la relación
Esto se puede satisfacer únicamente si el campo magnético tiene componentes en x e y desfasadas en π/2,
resultando b∗ = ±k. Por lo tanto, la condición que es necesario imponer sobre para tener una perturbación del tipo alfvénica se reduce en este caso a que la onda esté polarizada circularmente. En este caso para la frecuencia se obtiene esta expresión
Por otra parte, en el límite linealizado, cuando el término de Hall sea mucho menor que el ambipolar, considerando el caso con 0, la frecuencia ωLim se escribe así
Reemplazando el valor de C y F puede verse que existe una cantidad γ que evalúa el valor del amortiguamiento de las ondas por efecto del término ambipolar:
Puede verse que γ depende principalmente de la relación entre la “frecuencia de Alfvén” y la frecuencia de colisión entre los iones y los neutros. Debido a que para hallar el valor de ω Lim se consideró
para estos resultados se espera que las ondas se propaguen con amortiguación pequeña.
En el caso analizado en DP01, estos autores consideran conductividad finita y ondas linealmente polarizadas. Si se toma en SP01 el límite de conductividad infinita, el resultado del cálculo del valor del amortiguamiento γ de DP01 resulta idéntico a (36). Es importante remarcar que el motivo por el que se adoptó la polarización circular en este análisis reside en que las ondas de Alfvén con término de Hall se propagan si se encuentran circularmente polarizadas [7,12].
CONCLUSIONES
Se ha encontrado que para plasmas parcialmente ionizados, bajo ciertas condiciones, existen soluciones no linealizadas que se comportan con las características o propiedades de las ondas de Alfvén de corte. La perturbación en densidad de corriente y su rotor satisfacen una relación de proporcionalidad como sucedía en el caso de ondas de Alfvén en la magnetohidrodinámica con término de Hall. Cuando el plasma tiende a estar totalmente ionizado, la solución tiende a la hallada para ondas de Alfvén en la magnetohidrodinámica con término de Hall. En el límite linealizado se encuentra que, a pesar de que el término ambipolar es disipativo, bajo ciertas condiciones las ondas pueden propagarse con amortiguación pequeña.