Difusão ambipolar em plasmas

Equação da difusão ambipolar

Diferentemente de um gás de partículas neutras, um plasma (elétrons e íons), são menos livres ao se moverem por causa da atração coulombiana. Em um caso em que um plasma se movimenta elvolto em um gás neutro, os coeficientes de difusão dos elétrons e dos íons são tipicamente dados por

$D_{e}={\frac {k_{b}T_{e}}{m_{e}\nu _{e}}}$ e $D_{i}={\frac {h_{b}T_{i}}{m_{i}\nu _{i}}}$

onde $T_{e}$, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle T_i} , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_e} , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_i} , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nu_e} e $\nu _{i}$, são as temperaturas, massas e frequências de colisão dos elétrons e íons com os isótopos dos átomos neutros. Devido à massa do elétron ser muito menor que a massa de um íon, $D_{e}$ é maior que $D_{i}$, então quando um plasma começa a se espalhar, incialmente os elétrons se espalham mais rapidamente que os íons e isso gera um campo elétrico que freia os elétron e acelera os íons.

(botar uma figura aqui)

O fluxo de elétron e íons pode ser escrito como

${\vec {J}}_{e}=-\mu _{e}n_{e}({\vec {r}},t){\vec {E}}({\vec {r}},t)-D_{e}{\vec {\nabla }}n_{e}({\vec {r}},t)$

${\vec {J}}_{i}=\mu _{i}n_{i}({\vec {r}},t){\vec {E}}({\vec {r}},t)-D_{i}{\vec {\nabla }}n_{i}({\vec {r}},t)$

onde $\mu _{e}$ e $\mu _{i}$ são as razões entre as velocidades médias dos elétron e dos íons e o campo elétrico.

O Método

Para resolver a Difusão ambipolar em plasmas em 1D utilizamos o método FTCS(Forward Time Central Space).  

${\frac {\partial n}{\partial t}}=-D_{a}{\frac {\partial ^{2}n}{\partial x^{2}}}$

${\frac {n(x,t+dt)-n(x,t)}{dt}}=-D_{a}{\frac {n(x+dx,t)-2n(x,t)+n(x-dx,t)}{dx^{2}}}$

${\frac {\partial ^{2}f}{\partial s^{2}}}\rightarrow {\frac {f_{i-1}-2f_{i}+f_{i+1}}{(\Delta s)^{2}}}$