Plasma modeling

Last updated

Plasma modeling refers to solving equations of motion that describe the state of a plasma. It is generally coupled with Maxwell's equations for electromagnetic fields or Poisson's equation for electrostatic fields. There are several main types of plasma models: single particle, kinetic, fluid, hybrid kinetic/fluid, gyrokinetic and as system of many particles.

Contents

Chart for modeling plasma Plasma Modeling Table.png
Chart for modeling plasma

Single particle description

The single particle model describes the plasma as individual electrons and ions moving in imposed (rather than self-consistent) electric and magnetic fields. The motion of each particle is thus described by the Lorentz Force Law. In many cases of practical interest, this motion can be treated as the superposition of a relatively fast circular motion around a point called the guiding center and a relatively slow drift of this point.

Kinetic description

The kinetic model is the most fundamental way to describe a plasma, resultantly producing a distribution function

where the independent variables and are position and velocity, respectively. A kinetic description is achieved by solving the Boltzmann equation or, when the correct description of long-range Coulomb interaction is necessary, by the Vlasov equation which contains self-consistent collective electromagnetic field, or by the Fokker–Planck equation, in which approximations have been used to derive manageable collision terms. The charges and currents produced by the distribution functions self-consistently determine the electromagnetic fields via Maxwell's equations.

Fluid description

To reduce the complexities in the kinetic description, the fluid model describes the plasma based on macroscopic quantities (velocity moments of the distribution such as density, mean velocity, and mean energy). The equations for macroscopic quantities, called fluid equations, are obtained by taking velocity moments of the Boltzmann equation or the Vlasov equation. The fluid equations are not closed without the determination of transport coefficients such as mobility, diffusion coefficient, averaged collision frequencies, and so on. To determine the transport coefficients, the velocity distribution function must be assumed/chosen. But this assumption can lead to a failure of capturing some physics.

Hybrid kinetic/fluid description

Although the kinetic model describes the physics accurately, it is more complex (and in the case of numerical simulations, more computationally intensive) than the fluid model. The hybrid model is a combination of fluid and kinetic models, treating some components of the system as a fluid, and others kinetically. The hybrid model is sometimes applied in space physics, when the simulation domain exceeds thousands of ion gyroradius scales, making it impractical to solve kinetic equations for electrons. In this approach, magnetohydrodynamic fluid equations describe electrons, while the kinetic Vlasov equation describes ions. [1] [2]

Gyrokinetic description

In the gyrokinetic model, which is appropriate to systems with a strong background magnetic field, the kinetic equations are averaged over the fast circular motion of the gyroradius. This model has been used extensively for simulation of tokamak plasma instabilities (for example, the GYRO and Gyrokinetic ElectroMagnetic codes), and more recently in astrophysical applications.

Quantum mechanical methods

Quantum methods are not yet very common in plasma modeling. They can be used to solve unique modeling problems; like situations where other methods do not apply. [3] They involve the application of quantum field theory to plasma. In these cases, the electric and magnetic fields made by particles are modeled like a field; A web of forces. Particles that move, or are removed from the population push and pull on this web of forces, this field. The mathematical treatment for this involves Lagrangian mathematics.

Commercial plasma physics modeling codes

See also

Related Research Articles

<span class="mw-page-title-main">Force</span> Influence that can change motion of an object

In physics, a force is an influence that can cause an object to change its velocity, i.e., to accelerate, meaning a change in speed or direction, unless counterbalanced by other forces. The concept of force makes the everyday notion of pushing or pulling mathematically precise. Because the magnitude and direction of a force are both important, force is a vector quantity. The SI unit of force is the newton (N), and force is often represented by the symbol F.

<span class="mw-page-title-main">Magnetic sail</span> Proposed spacecraft propulsion method

A magnetic sail is a proposed method of spacecraft propulsion where a spacecraft shown as a purple dot in the center of the illustration has a source that generates a magnetic field, which under certain conditions, summarized in the overview section, creates a magnetopause and a bow shock that deflect a plasma wind of charged particles. A low density region forms creating an artificial magnetospheric bubble around the spacecraft. The term magnetospheric object refers to the magnetosphere, magnetopause and the bow shock. Charged particles in the plasma deflected by the bow shock and magnetopause together create an effective sail blocking area that exerts a wind force that acts on the magnetic field, which in turn exerts a magnetic force on the field source that creates a spacecraft force that accelerates the spacecraft in the same direction as the plasma wind.

<span class="mw-page-title-main">Magnetohydrodynamics</span> Model of electrically conducting fluids

Magnetohydrodynamics is a model of electrically conducting fluids that treats all interpenetrating particle species together as a single continuous medium. It is primarily concerned with the low-frequency, large-scale, magnetic behavior in plasmas and liquid metals and has applications in numerous fields including geophysics, astrophysics, and engineering.

Plasma diagnostics are a pool of methods, instruments, and experimental techniques used to measure properties of a plasma, such as plasma components' density, distribution function over energy (temperature), their spatial profiles and dynamics, which enable to derive plasma parameters.

In plasma physics, the particle-in-cell (PIC) method refers to a technique used to solve a certain class of partial differential equations. In this method, individual particles in a Lagrangian frame are tracked in continuous phase space, whereas moments of the distribution such as densities and currents are computed simultaneously on Eulerian (stationary) mesh points.

<span class="mw-page-title-main">Magnetic reconnection</span> Process in plasma physics

Magnetic reconnection is a physical process occurring in electrically conducting plasmas, in which the magnetic topology is rearranged and magnetic energy is converted to kinetic energy, thermal energy, and particle acceleration. Magnetic reconnection involves plasma flows at a substantial fraction of the Alfvén wave speed, which is the fundamental speed for mechanical information flow in a magnetized plasma.

<span class="mw-page-title-main">Guiding center</span>

In physics, the motion of an electrically charged particle such as an electron or ion in a plasma in a magnetic field can be treated as the superposition of a relatively fast circular motion around a point called the guiding center and a relatively slow drift of this point. The drift speeds may differ for various species depending on their charge states, masses, or temperatures, possibly resulting in electric currents or chemical separation.

In plasma physics, waves in plasmas are an interconnected set of particles and fields which propagate in a periodically repeating fashion. A plasma is a quasineutral, electrically conductive fluid. In the simplest case, it is composed of electrons and a single species of positive ions, but it may also contain multiple ion species including negative ions as well as neutral particles. Due to its electrical conductivity, a plasma couples to electric and magnetic fields. This complex of particles and fields supports a wide variety of wave phenomena.

The Vlasov equation is a differential equation describing time evolution of the distribution function of plasma consisting of charged particles with long-range interaction, such as the Coulomb interaction. The equation was first suggested for the description of plasma by Anatoly Vlasov in 1938 and later discussed by him in detail in a monograph.

The gyroradius is the radius of the circular motion of a charged particle in the presence of a uniform magnetic field. In SI units, the non-relativistic gyroradius is given by

In plasma physics, the Hasegawa–Mima equation, named after Akira Hasegawa and Kunioki Mima, is an equation that describes a certain regime of plasma, where the time scales are very fast, and the distance scale in the direction of the magnetic field is long. In particular the equation is useful for describing turbulence in some tokamaks. The equation was introduced in Hasegawa and Mima's paper submitted in 1977 to Physics of Fluids, where they compared it to the results of the ATC tokamak.

Gyrokinetic ElectroMagnetic (GEM) is a gyrokinetic plasma turbulence simulation that uses the particle-in-cell method. It is used to study waves, instabilities and nonlinear behavior of tokamak fusion plasmas. Information about GEM can be found at the GEM web page. There are two versions of GEM, one is a flux-tube version and the other one is a global general geometry version. Both versions of GEM use a field-aligned coordinate system. Ions are treated kinetically, but averaged over their gyro-obits and electrons are treated as drift-kinetic.

Gyrokinetics is a theoretical framework to study plasma behavior on perpendicular spatial scales comparable to the gyroradius and frequencies much lower than the particle cyclotron frequencies. These particular scales have been experimentally shown to be appropriate for modeling plasma turbulence. The trajectory of charged particles in a magnetic field is a helix that winds around the field line. This trajectory can be decomposed into a relatively slow motion of the guiding center along the field line and a fast circular motion, called gyromotion. For most plasma behavior, this gyromotion is irrelevant. Averaging over this gyromotion reduces the equations to six dimensions rather than the seven. Because of this simplification, gyrokinetics governs the evolution of charged rings with a guiding center position, instead of gyrating charged particles.

<span class="mw-page-title-main">Adolfo Figueroa Viñas</span> Puerto Rican astrophysicist

Adolfo Figueroa Viñas is the first Puerto Rican astrophysicist at the National Aeronautics and Space Administration (NASA) and is an expert in solar and space plasma physics at the Heliophysics Science Division. As a staff scientist his research interests include studying plasma kinetic physics and magnetohydrodynamics of the solar wind, heliosphere, shock waves, MHD and kinetic simulation of plasma instabilities, and turbulent processes associated with space, solar and astrophysical plasmas.

<span class="mw-page-title-main">Liu Chen (physicist)</span> American physicist

Liu Chen is an American theoretical physicist who has made original contributions to many aspects of plasma physics. He is known for the discoveries of kinetic Alfven waves, toroidal Alfven eigenmodes, and energetic particle modes; the theories of geomagnetic pulsations, Alfven wave heating, and fishbone oscillations, and the first formulation of nonlinear gyrokinetic equations. Chen retired from University of California, Irvine (UCI) in 2012, assuming the title professor emeritus of physics and astronomy.

<span class="mw-page-title-main">Plasma (physics)</span> State of matter

Plasma is one of four fundamental states of matter, characterized by the presence of a significant portion of charged particles in any combination of ions or electrons. It is the most abundant form of ordinary matter in the universe, mostly in stars, but also dominating the rarefied intracluster medium and intergalactic medium. Plasma can be artificially generated, for example, by heating a neutral gas or subjecting it to a strong electromagnetic field.

This glossary of physics is a list of definitions of terms and concepts relevant to physics, its sub-disciplines, and related fields, including mechanics, materials science, nuclear physics, particle physics, and thermodynamics. For more inclusive glossaries concerning related fields of science and technology, see Glossary of chemistry terms, Glossary of astronomy, Glossary of areas of mathematics, and Glossary of engineering.

The Farley–Buneman instability, or FB instability, is a microscopic plasma instability named after Donald T. Farley and Oscar Buneman. It is similar to the ionospheric Rayleigh-Taylor instability.

Bernstein–Greene–Kruskal modes are nonlinear electrostatic waves that propagate in an unmagnetized, collisionless plasma. They are nonlinear solutions to the Vlasov–Poisson equation in plasma physics, and are named after physicists Ira B. Bernstein, John M. Greene, and Martin D. Kruskal, who solved and published the exact solution for the one-dimensional case in 1957.

Harold Weitzner is an American applied mathematician and physicist whose primary research is plasma physics. He is Professor Emeritus of Mathematics at the Courant Institute of Mathematical Sciences and has served as Director of the Magneto-Fluid Dynamics Division at Courant since 1981, succeeding Harold Grad. He has published over 120 research articles on the topics of plasma physics, magnetohydrodynamics, fluid mechanics, fractional equations and kinetics, and chaos.

References

  1. Pokhotelov, D.; von Alfthan, S.; Kempf, Y.; Vainio, R.; Koskinen, H.E.J.; Palmroth, M. (2013). "Ion distributions upstream and downstream of the Earth's bow shock: First results from Vlasiator". Ann. Geophys. 31 (12): 2207–2212. doi: 10.5194/angeo-31-2207-2013 . hdl: 10138/161497 .
  2. von Alfthan, S.; Pokhotelov, D.; Kempf, Y.; Hoilijoki, S.; Honkonen, I.; Sandroos, A.; Palmroth, M. (2014). "Vlasiator: First global hybrid-Vlasov simulations of Earth's foreshock and magnetosheath". J. Atmos. Sol. Terr. Phys. 120: 24–35. doi:10.1016/j.jastp.2014.08.012.
  3. Hedditch, John (2018). "A different approach to the MHD equilibrium". arXiv: 1808.00622 [physics.plasm-ph].