Evolution Equation of the Free Expansion State

In the field of plasma physics, free expansion refers to the spontaneous outward movement of a plasma cloud into a vacuum or a surrounding medium of significantly lower density. This process occurs in the absence of external confining magnetic fields or driving forces, driven primarily by the internal thermal pressure of the plasma. Understanding the evolution of this state is critical for various applications, ranging from laser-plasma interactions and inertial confinement fusion to the development of plasma-based propulsion systems.

To mathematically describe the complex dynamics of an expanding plasma, researchers typically employ two primary frameworks: the fluid model, which treats the plasma as a continuous medium, and the kinetic model, which tracks the evolution of particle distribution functions in phase space.

Fundamental Assumptions and Physical Regime

To render the governing equations tractable for analytical and numerical studies, several simplifying assumptions are commonly adopted:

• Collisionless Approximation: The plasma is assumed to be in a regime where the mean free path between particle collisions is much larger than the characteristic system scale. This allows the use of the Vlasov equation or inviscid fluid equations.
• One-Dimensional Planar Symmetry: The expansion is modeled along a single axis (typically the $x$-axis), assuming that all physical quantities vary only with position $x$ and time $t$.
• Quasi-neutrality: On scales much larger than the Debye length ($\lambda_D$), the plasma is assumed to maintain local charge neutrality, such that $n_e \approx Z n_i$.
• Single-Temperature Approximation: For simplicity, the electrons and ions are often assumed to share the same temperature, or their respective temperatures are treated independently while ignoring complex heat flux terms between species.

The Fluid Dynamics Framework

The macroscopic evolution of the free expansion can be described by a set of coupled conservation laws: the continuity, momentum, and energy equations.

2.1 Conservation of Mass (Continuity Equation)

The evolution of the particle density $n(x,t)$ is governed by the continuity equation, which ensures the conservation of mass:
$$\frac{\partial n}{\partial t} + \frac{\partial}{\partial x}(n u) = 0$$
where $u(x,t)$ represents the local fluid velocity.

2.2 Conservation of Momentum

The momentum equation describes how the velocity field evolves under the influence of pressure gradients:
$$\frac{\partial}{\partial t}(n u) + \frac{\partial}{\partial x}(n u^2 + p) = 0$$
Here, $p$ is the thermal pressure. In an ideal plasma, this is closed by the equation of state $p = n k_{\mathrm{B}} T$.

2.3 Conservation of Energy and the Adiabatic Assumption

In a free expansion where no external work is performed and thermal conduction is negligible, the process is considered adiabatic. The energy evolution is expressed as:
$$\frac{\partial}{\partial t}\left(\frac{p}{\gamma-1}\right) + \frac{\partial}{\partial x}\left(\frac{p}{\gamma-1}u + pu\right) = 0$$
where $\gamma$ is the adiabatic index (for a monoatomic electron-ion plasma, $\gamma = 5/3$).

2.4 Normalization for Scaling Analysis

To facilitate numerical simulations and reveal universal behaviors, variables are typically normalized using characteristic plasma scales:

• Position: $\tilde{x} = x/\lambda_D$
• Time: $\tilde{t} = \omega_{pi} t$ (where $\omega_{pi}$ is the ion plasma frequency)
• Velocity: $\tilde{u} = u/c_s$ (where $c_s = \sqrt{\gamma k_{\mathrm{B}} T_0/m_i}$ is the sound speed)
• Density and Pressure: $\tilde{n} = n/n_0$ and $\tilde{p} = p/(n_0 k_{\mathrm{B}} T_0)$

Self-Similar Solutions

One of the most profound features of free expansion is its self-similarity. As the plasma expands, the spatial profiles of density, velocity, and pressure maintain their functional form, merely scaling with time. This allows us to transform the partial differential equations (PDEs) into ordinary differential equations (ODEs) by introducing the dimensionless similarity variable:
$$\xi = \frac{x}{c_s t}$$

By substituting this variable into the fluid equations, we derive a closed system of ODEs:
$$\begin{aligned}
(1-\xi) \frac{d\tilde{n}}{d\xi} + \tilde{n}\frac{d\tilde{u}}{d\xi} &= 0 \
(1-\xi) \frac{d\tilde{u}}{d\xi} + \frac{1}{\tilde{n}}\frac{d\tilde{p}}{d\xi} &= 0 \
(1-\xi) \frac{d\tilde{p}}{d\xi} + \gamma \tilde{p}\frac{d\tilde{u}}{d\xi} &= 0
\end{aligned}$$

Using the relation $\tilde{p} = \tilde{n}^\gamma$, these equations yield the following analytical solutions:

• Velocity Profile: $\tilde{u}(\xi) = \frac{2}{\gamma+1}(\xi - 1)$
• Density Profile: $\tilde{n}(\xi) = \left[ 1 - \frac{\gamma-1}{\gamma+1}\xi \right]^{\frac{2}{\gamma-1}}$
• Pressure Profile: $\tilde{p}(\xi) = \tilde{n}(\xi)^\gamma$

These solutions are valid for $\xi \le \xi_{\max} = \frac{2}{\gamma-1}$. At the boundary $\xi_{\max}$, the density drops to zero, marking the rarefaction front—the interface between the expanding plasma and the vacuum.

Example: The $\gamma = 5/3$ Case
For a standard electron-ion plasma, the solutions simplify to:
$$\tilde{u}(\xi) = \frac{3}{4}(\xi-1), \quad \tilde{n}(\xi) = \left(1 - \frac{1}{4}\xi\right)^3$$
At the expansion origin ($\xi=0$), the velocity is negative ($\tilde{u} = -3/4$), indicating a residual inward flow relative to the expanding frame, while the density reaches zero at $\xi=4$.

Kinetic Perspective: The Vlasov-Poisson System

While the fluid model captures the macroscopic "bulk" motion, it cannot describe fine-grained phenomena such as electron escape or ion acceleration at the expansion front. For these effects, a kinetic description is required using the Vlasov equation for each species $s$ (electrons and ions):
$$\frac{\partial f_s}{\partial t} + v\frac{\partial f_s}{\partial x} + \frac{q_s}{m_s}E\frac{\partial f_s}{\partial v} = 0$$
The self-consistent electric field $E(x,t)$ is determined by the Poisson equation:
$$\frac{\partial E}{\partial x} = \frac{e}{\varepsilon_0}(Zn_i - n_e)$$
In the limit of rapid expansion, the electric field is often concentrated in a very thin layer at the expansion front (the thin-sheath approximation). This field acts as a driver that accelerates ions forward while allowing electrons to expand more rapidly, creating a charge separation that sustains the expansion.

Numerical Implementation Strategies

Simulating the evolution of the free expansion state requires careful attention to numerical stability and physical accuracy. Key considerations include:

• Conservation Laws: It is essential to use conservative schemes (such as Finite Volume Methods) to ensure that mass, momentum, and energy are strictly preserved during the simulation.
• Grid Management: Because the rarefaction front moves and sharp gradients develop, adaptive mesh refinement (AMR) or Lagrangian tracking is highly effective for maintaining resolution at the expansion edge.
• Boundary Conditions: To prevent unphysical reflections of waves, the vacuum boundary should be treated with open or absorbing boundary conditions.
• Numerical Dissipation: High-order limiters (e.g., WENO or TVD schemes) should be employed to suppress oscillations near the steep density gradients of the expansion front.

Summary

The evolution of a free expansion state is a fundamental problem characterized by a transition from high-pressure equilibrium to a self-similar rarefaction wave. While the fluid model provides elegant analytical solutions through self-similar scaling, the kinetic model remains indispensable for capturing the micro-physics of particle acceleration and sheath formation. Mastery of these evolution equations provides the theoretical foundation necessary for exploring advanced plasma phenomena in high-energy-density physics.