diff options
Diffstat (limited to 'content/know/concept/two-fluid-equations')
| -rw-r--r-- | content/know/concept/two-fluid-equations/index.pdc | 178 |
1 files changed, 178 insertions, 0 deletions
diff --git a/content/know/concept/two-fluid-equations/index.pdc b/content/know/concept/two-fluid-equations/index.pdc new file mode 100644 index 0000000..cd77f5e --- /dev/null +++ b/content/know/concept/two-fluid-equations/index.pdc @@ -0,0 +1,178 @@ +--- +title: "Two-fluid equations" +firstLetter: "T" +publishDate: 2021-10-19 +categories: +- Physics +- Plasma physics + +date: 2021-10-18T10:12:20+02:00 +draft: false +markup: pandoc +--- + +# Two-fluid equations + +The **two-fluid model** describes a plasma as two separate but overlapping fluids, +one for ions and one for electrons. +Instead of tracking individual particles, +it gives the dynamics of fluid elements $\dd{V}$ (i.e. small "blobs"). +These blobs are assumed to be much larger than +the [Debye length](/know/concept/debye-length/), +such that electromagnetic interactions between nearby blobs can be ignored. + +From Newton's second law, we know that the velocity $\vb{v}$ +of a particle with mass $m$ and charge $q$ is as follows, +when subjected only to the [Lorentz force](/know/concept/lorentz-force/): + +$$\begin{aligned} + m \dv{\vb{v}}{t} + = q (\vb{E} + \vb{v} \cross \vb{B}) +\end{aligned}$$ + +From here, the derivation is similar to that of the +[Navier-Stokes equations](/know/concept/navier-stokes-equations/). +We replace the time derivative with a +[material derivative](/know/concept/material-derivative/) $\mathrm{D}/\mathrm{D}t$, +and define a blob's velocity $\vb{u}$ +as the average velocity of the particles inside it, leading to: + +$$\begin{aligned} + m n \frac{\mathrm{D} \vb{u}}{\mathrm{D} t} + = q n (\vb{E} + \vb{u} \cross \vb{B}) +\end{aligned}$$ + +Where we have multiplied by the number density $n$ of the particles. +Due to particle collisions in the fluid, +stresses become important. Therefore, we include +the [Cauchy stress tensor](/know/concept/cauchy-stress-tensor/) $\hat{P}$, +leading to the following two equations: + +$$\begin{aligned} + m_i n_i \frac{\mathrm{D} \vb{u}_i}{\mathrm{D} t} + &= q_i n_i (\vb{E} + \vb{u}_i \cross \vb{B}) + \nabla \cdot \hat{P}_i{}^\top + \\ + m_e n_e \frac{\mathrm{D} \vb{u}_e}{\mathrm{D} t} + &= q_e n_e (\vb{E} + \vb{u}_e \cross \vb{B}) + \nabla \cdot \hat{P}_e{}^\top +\end{aligned}$$ + +Where the subscripts $i$ and $e$ refer to ions and electrons, respectively. +Finally, we also account for momentum transfer between ions and electrons +due to [Rutherford scattering](/know/concept/rutherford-scattering/), +leading to these **two-fluid momentum equations**: + +$$\begin{aligned} + \boxed{ + \begin{aligned} + m_i n_i \frac{\mathrm{D} \vb{u}_i}{\mathrm{D} t} + &= q_i n_i (\vb{E} + \vb{u}_i \cross \vb{B}) + \nabla \cdot \hat{P}_i{}^\top - f_{ie} m_i n_i (\vb{u}_i - \vb{u}_e) + \\ + m_e n_e \frac{\mathrm{D} \vb{u}_e}{\mathrm{D} t} + &= q_e n_e (\vb{E} + \vb{u}_e \cross \vb{B}) + \nabla \cdot \hat{P}_e{}^\top - f_{ei} m_e n_e (\vb{u}_e - \vb{u}_i) + \end{aligned} + } +\end{aligned}$$ + +Where $f_{ie}$ is the mean frequency at which an ion collides with electrons, +and vice versa for $f_{ei}$. +For simplicity, we assume that the plasma is isotropic +and that shear stresses are negligible, +in which case the stress term can be replaced +by the gradient $- \nabla p$ of a scalar pressure $p$: + +$$\begin{aligned} + m_i n_i \frac{\mathrm{D} \vb{u}_i}{\mathrm{D} t} + &= q_i n_i (\vb{E} + \vb{u}_i \cross \vb{B}) - \nabla p_i - f_{ie} m_i n_i (\vb{u}_i - \vb{u}_e) + \\ + m_e n_e \frac{\mathrm{D} \vb{u}_e}{\mathrm{D} t} + &= q_e n_e (\vb{E} + \vb{u}_e \cross \vb{B}) - \nabla p_e - f_{ei} m_e n_e (\vb{u}_e - \vb{u}_i) +\end{aligned}$$ + +Next, we demand that matter is conserved. +In other words, the rate at which particles enter/leave a volume $V$ +must be equal to the flux through the enclosing surface $S$: + +$$\begin{aligned} + 0 + &= \pdv{t} \int_V n \dd{V} + \oint_S n \vb{u} \cdot \dd{\vb{S}} + = \int_V \Big( \pdv{n}{t} + \nabla \cdot (n \vb{u}) \Big) \dd{V} +\end{aligned}$$ + +Where we have used the divergence theorem. +Since $V$ is arbitrary, we can remove the integrals, +leading to the following **continuity equations**: + +$$\begin{aligned} + \boxed{ + \pdv{n_i}{t} + \nabla \cdot (n_i \vb{u}_i) + = 0 + \qquad \quad + \pdv{n_e}{t} + \nabla \cdot (n_e \vb{u}_e) + = 0 + } +\end{aligned}$$ + +Currently, we have 8 equations (2 scalar continuity, 2 vector momentum), +but 16 unknowns $\vb{u}_i$, $\vb{u}_e$, $\vb{E}$, $\vb{B}$, $n_i$, $n_e$, $p_i$ and $p_e$. +We would like to close this system, so we need 8 more. +An obvious choice is [Maxwell's equations](/know/concept/maxwells-equations/), +in particular Faraday's and Ampère's law +(since Gauss' laws are in fact redundant; see the article on Maxwell's equations): + +$$\begin{aligned} + \boxed{ + \nabla \cross \vb{E} = - \pdv{\vb{B}}{t} + \qquad \quad + \nabla \cross \vb{B} = \mu_0 \Big( n_i q_i \vb{u}_i + n_e q_e \vb{u}_e + \varepsilon_0 \pdv{\vb{E}}{t} \Big) + } +\end{aligned}$$ + +Now we have 14 equations, so we need 2 more, for the pressures $p_i$ and $p_e$. +This turns out to be the thermodynamic **equation of state**: +for quasistatic, reversible, adiabatic compression +of a gas with constant heat capacities (i.e. a *calorically perfect* gas), +it turns out that: + +$$\begin{aligned} + \dv{t} \big( p V^\gamma \big) = 0 + \qquad \quad + \gamma + \equiv \frac{C_P}{C_V} + = \frac{N + 2}{N} +\end{aligned}$$ + +Where $\gamma$ is the *heat capacity ratio*, +and can be calculated from the number of degrees of freedom $N$ +of each particle in the gas. +In a fully ionized plasma, $N = 3$. + +The density $n \propto 1/V$, +so since $p V^\gamma$ is constant in time, +for some (known) constant $C$: + +$$\begin{aligned} + \dv{t} \Big( \frac{p}{n^\gamma} \Big) = 0 + \quad \implies \quad + p = C n^\gamma +\end{aligned}$$ + +In the two-fluid model, we thus have the following two equations of state, +giving us a set of 16 equations for 16 unknowns: + +$$\begin{aligned} + \boxed{ + p_i = C_i n_i^\gamma + \qquad \quad + p_e = C_e n_e^\gamma + } +\end{aligned}$$ + + + +## References +1. F.F. Chen, + *Introduction to plasma physics and controlled fusion*, + 3rd edition, Springer. +2. M. Salewski, A.H. Nielsen, + *Plasma physics: lecture notes*, + 2021, unpublished. |