Difusão ambipolar em plasmas

De Física Computacional
Edição feita às 01h24min de 1 de abril de 2021 por Leirbag (Discussão | contribs)

Ir para: navegação, pesquisa

Equação da difusão ambipolar

A difusão é o modo como um fluido de dilui em um meio. Estudar as equações que governam esse fenômeno e as formas de resolvê-las é de extremo interesse para a física de fluidos e de plasmas, entre outras áreas. Aqui mostramos uma resolução numérica para o caso unidimensional da difusão ambipolar de um plasma (gás formado de elétrons e íons) envolto em um gás neutro, ou seja, o caso de um plasma se espalhando por um tubo.

Diferentemente de um gás de átomos/moléculas neutros(as), os plasmas são menos livres ao se moverem por causa das interações eletromagnéticas envolvidas no movimento das cargas, como a força de Coulomb e a força magnética. Na difusão de plasmas em um gás neutro, os coeficientes de difusão dos elétrons e dos íons são tipicamente dados por

e

onde , , , , e , são as temperaturas, massas e frequências de colisão dos elétrons e íons com os átomos neutros. Devido à massa do elétron ser muito menor que a massa de um íon, é maior que , então quando um plasma começa a se difundir, incialmente os elétrons se espalham mais rapidamente que os íons, o que gera um campo elétrico que freia os elétron e acelera os íons. Chamamos esse processo de difusão ambipolar.

Difusao ambipolar.png [1]

Como mostrado por Shimony e Cahn[2], esse problema é descrito por uma equação de onda amortecida

onde e , sendo a frequência de colisão ambipolar e o coeficiente de difusão ambipolar, que pode ser escrito como [3].

Como tratamos do caso unidimensional, a equação 1 torna-se

O Método

A resolução numérica do problema foi baseada no artigo de Najafi e Izadi [4]. Começamos com a forma mais usual de escrever a equação da onda amortecida unidimencional


No nosso caso e .

Discretizando as variáveis do problema, temos que

Substituindo as derivadas que aparecem na equação por diferenças finitas, obtemos

Substituindo essas relações na equação 3, obtemos

Omitindo todos os temos de ordem e isolando , obtemos

sendo .

Essa é a equação para resolver o problema para , mas nessecitamos ainda de uma maneira de determinar a paritr de . Para isso assumimos que a função é inicialmente estacionária e fazemos

Substituindo na equação 4 para obtemos

Com as equações 4 e 5, e tomando as devidas condições de contorno nas bordas, podemos calcular a evolução temporal da função de densidade. Esse método é estável para e seu erro é

Resultados e Discussão

Evolução temporal da densidade do plasma para diferentes valores de e
Alt text
Alt text
Alt text
Alt text

Programas Utilizados

Referências

  1. http://www.enigmatic-consulting.com/semiconductor_processing/CVD_Fundamentals/plasmas/ambipolar_diffusion.html
  2. Z. Shimony and J. H. Cahn, "Time-dependent ambipolar diffusion waves", The Physics of Fluids 8, 1704 (1965)
  3. http://uigelz.eecs.umich.edu/classes/pub/eecs517/handouts/derivation_ambipolar_diffusion_v02.pdf
  4. H. Najafi and F. Izadi, "Comparison of two finite-difference methods for solving the damped wave equation", viXra, 2016