Plasma Simulator Back
Electromagnetics & Plasma

Plasma Simulator

Simulate charged particle motion under F = q(E + v×B) in real time. Observe cyclotron orbits, E×B drift, and plasma dynamics with interactive field controls.

Field Settings

Statistics

Results
1.00
ωc (electron)
1.00
Larmor radius
0.50
E×B speed
Avg KE
Visualization
Theory & Key Formulas

$$\omega_p = \sqrt{\frac{n_e e^2}{\varepsilon_0 m_e}}$$

電子プラズマ振動数:\(n_e\) 電子数密度 [m⁻³]、\(e\) 電荷、\(m_e\) 電子質量

$$\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T_e}{n_e e^2}}$$

デバイ長(遮蔽長):電場が遮蔽される特性長さ [m]

$$\mathbf{F} = q(\mathbf{E} + \mathbf{v}\times\mathbf{B})$$

ローレンツ力:荷電粒子の運動方程式の右辺

What is Plasma Motion?

🙋
What exactly is the "Lorentz force" that makes the particles in this simulator move in circles?
🎓
Basically, it's the total electromagnetic force on a charged particle. The formula is \(\vec{F}= q(\vec{E}+ \vec{v}\times \vec{B})\). The magnetic part (\(q \vec{v}\times \vec{B\)) is key for the circular motion—it's always perpendicular to the particle's velocity, causing it to constantly turn. Try setting the Electric Field slider to zero and moving the Magnetic Field slider to see a perfect circular orbit.
🙋
Wait, really? So the circle size depends on the magnetic field strength? What if I make the particle heavier?
🎓
Exactly! A stronger magnetic field (higher B) pulls the particle into a tighter turn, making a smaller circle. A heavier particle (higher mass) is harder to deflect, so it makes a wider circle. This radius is called the Larmor or gyroradius. In the simulator, try increasing the Mass slider while keeping B constant—you'll see the orbit grow larger instantly.
🙋
That's cool. But the particles sometimes drift sideways instead of going in a simple circle. What's happening there?
🎓
Great observation! That's the famous E×B drift. When you have both an electric field (E) and a magnetic field (B) at the same time, the particle's circular motion gets a net sideways push. It's like the circle itself is sliding. For instance, in the simulator, set a moderate Electric Field in the X-direction and a Magnetic Field in the Z-direction. Watch the particle's guiding center drift in the Y-direction—that's the E×B drift in action!

Physical Model & Key Equations

The motion of each charged particle is governed by Newton's second law, with the Lorentz force as the sole electromagnetic influence. This is the fundamental equation of motion solved in real-time by the simulator.

$$ m \frac{d\vec{v}}{dt}= q (\vec{E}+ \vec{v}\times \vec{B}) $$

Where \(m\) is particle mass, \(\vec{v}\) is velocity, \(q\) is charge, \(\vec{E}\) is the electric field vector, and \(\vec{B}\) is the magnetic field vector. The cross product \(\vec{v}\times \vec{B}\) ensures the magnetic force is perpendicular to both velocity and field.

From the Lorentz force, we derive two key parameters for pure magnetic field motion (when \(\vec{E}=0\)): the cyclotron frequency and the Larmor radius. These define the orbit you see.

$$ \omega_c = \frac{|q| B}{m}, \quad r_L = \frac{m v_\perp}{|q| B} $$

Here, \(\omega_c\) is the cyclotron (or gyro) frequency—how fast the particle goes around the circle. \(r_L\) is the Larmor radius—the size of the circle. \(v_\perp\) is the velocity component perpendicular to the magnetic field. Notice that \(\omega_c\) does not depend on speed, but \(r_L\) does.

Frequently Asked Questions

The initial velocity of the particles may be too high, or the settings for the magnetic and electric fields may be too strong. Please adjust the initial velocity of the particles or the strength of the electric and magnetic fields to lower values using the control panel at the bottom of the screen. You can also move the viewpoint to track the particles.
The direction of E×B drift is determined by the cross product of the electric field E and the magnetic field B (E×B). It does not depend on the sign of the charge and drifts in a direction perpendicular to both E and B. In the simulator, you can intuitively confirm the drift direction by changing the vectors of the electric and magnetic fields.
According to the Larmor radius formula r_L = mv⊥/(|q|B), to increase the radius, either increase the perpendicular velocity v⊥ of the particle or weaken the magnetic field B. In the simulator, you can adjust this by raising the slider for the perpendicular component of the initial velocity or by lowering the magnetic field strength.
Yes. You can visualize the confinement of charged particles by magnetic fields (the basics of cyclotron motion and mirror effects). However, please note that this simulator calculates single particle trajectories and does not simulate multi-particle collisions or the overall behavior of the plasma.

Real-World Applications

Magnetic Confinement Fusion (Tokamaks): This is the principle behind reactors like ITER. Powerful magnetic fields are used to confine a hot plasma (like the particles in the simulator) in a donut-shaped ring, preventing it from touching and melting the reactor walls. The particles' circular orbits keep them trapped within the magnetic "bottle."

Mass Spectrometers: Devices that identify chemical substances by their mass. Ions are shot into a known magnetic field. Heavier ions (larger m) bend in larger radii (\(r_L \propto m\)), while lighter ones bend more sharply. By measuring the bend, the instrument determines the mass, just like changing the Mass slider changes the orbit size here.

Spacecraft Shielding (Active Concepts): To protect astronauts from deadly cosmic radiation on long Mars missions, one proposed idea is to generate a powerful magnetic field around the spacecraft. Like the magnetic field in the simulator deflecting a charged particle, this "magnetic shield" would deflect incoming high-energy ions away from the crew cabin.

Industrial Plasma Processing: Used in manufacturing computer chips and coating materials. Plasma (a soup of charged particles) is guided and controlled inside a vacuum chamber using precisely shaped electric and magnetic fields. Engineers use E×B drifts and confinement principles to direct the plasma exactly where it's needed for etching or deposition.

Common Misconceptions and Points to Note

First, keep in mind that this simulator shows the motion of a single particle. The plasmas you deal with in practice are collections of a vast number of particles, where collisions between particles and collective behavior (like plasma oscillations) are fundamentally important. While you learn from this tool that E×B drift is "independent of particle type," in actual plasmas, other drifts due to density gradients or magnetic field curvature occur, which make confinement more difficult.

Next, the pitfall of "real-time" display. The simulation speed depends on the computational load. For instance, if you set the particle mass extremely low (like an electron) or the magnetic field very strong, the cyclotron frequency \(\omega_c\) becomes very large. The numerical calculation might not keep up, causing the display to stutter or leading to numerical instability. In practical CAE, knowing how to stably solve such "stiff problems" is key.

Finally, the importance of initial conditions. In this tool, the particle's initial position and velocity are often fixed or set randomly. However, in simulations for something like a fusion reactor, particles are typically set to be generated according to a specific distribution function (e.g., a Maxwellian distribution). Don't jump to conclusions like "this is how plasma moves" just from looking at the trajectory of a single particle.

How to Use

  1. Set magnetic field strength (B) using BvalNum or Bslider—typical range 0.1–2.0 Tesla for lab plasmas.
  2. Define electric field components via ExvalNum, EyvalNum sliders in V/m (e.g., 1000–10000 V/m for ionosphere studies).
  3. Adjust particle count (NvalNum) from 10 to 1000 ions to observe collective drift behavior and density gradients.
  4. Click simulate to observe real-time trajectories governed by F = q(E + v × B).

Worked Example

Argon ions (m=40 amu, q=+1.6e-19 C) in B=0.5 T perpendicular field, Ex=2000 V/m: cyclotron frequency ωc ≈ 4.8 MHz, Larmor radius rL ≈ 0.18 mm at 1 km/s drift velocity. E×B drift velocity vE×B = E/B ≈ 4000 m/s perpendicular to both fields. Simulating N=500 particles reveals circular orbits superimposed on bulk drift, gyroradius expansion under acceleration.

Practical Notes

  1. Crossed E and B fields (perpendicular orientation) produce steady E×B drift independent of particle mass—used in Hall thrusters and tokamak confinement diagnostics.
  2. Reduce timestep if particles escape boundaries; instability occurs when ωc·Δt > 0.1 for accurate orbit resolution.
  3. High N (>500) and strong B reveal diamagnetic effects and pressure-gradient-driven instabilities relevant to magnetospheric substorms.
  4. Zero electric field isolates pure cyclotron motion—validate against ωc = qB/m before adding drifts.