Next Article in Journal
A Few Points of the Engineering Logic Discussed in ITER EDA on Evaluation of Halo- and AVDE-Induced Loads in Tokamaks
Next Article in Special Issue
Structural Characteristics of Ion Holes in Plasma
Previous Article in Journal
An Ion Source’s View of Its Plasma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formation of Kinetics Coherent Structures in Weakly Collisional Media

1
Department of Electrophysical Facilities, National Research Nuclear University MEPhI, Kashirskoye hwy, 31, Moscow 115409, Russia
2
Institute for High Temperatures, Russian Academy of Sciences, Izhorskaya 13/19, Moscow 127412, Russia
3
Research Institute of Quality, Safety and Technologies of Specialized Food Products of the PRUE G.V. Plekhanov, Stremyanny Lane, 36, Building 6, Moscow 117997, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 27 May 2021 / Revised: 5 June 2021 / Accepted: 9 June 2021 / Published: 18 June 2021

Abstract

:
The formation of nonlinear, nonstationary structures in weakly collisional media with collective interactions are investigated analytically within the framework of the kinetic description. This issue is considered in one-dimensional geometry using collision integral in the Bhatnagar-Gross-Krook form and some model forms of the interparticle interaction potentials that ensure the finiteness of the energy and momentum of the systems under consideration. As such potentials, we select the Yukawa potential, the δ -potential, which describes coherent structures in a plasma. For such potentials we obtained a dispersion relation which makes it possible to estimate the size and type of the forming structures.

1. Introduction

There are a lot of natural systems which are capable of exchanging energy and negentropy with the immediate environment [1,2]. In this case, the steady time-dependent states of the systems appear to be possible as a result of two effects: the formation of coherent states due to the self-consistent potential field and the formation of dissipative equilibria due to collisions. A well-known example of dissipative processes is the Beloysov-Ghabotinskyl reaction [2,3], whereas the formation of nonlinear wave-like vortex-like structures in the collisionless Coulomb plasmas [4,5] is the most vivid and vital sample of the system with the purely potential interaction between particles, which have been described by soliton-like distributions. In this case, one should keep in mind that such wave-like solutions possess some properties of particles [6,7]. For example, the polar molecules clouds can combine into hyperparticle structures which one may consider to be macroscopic atoms.
From this perspective, we could suggest that the internal structure as well as the physical properties of such hyperparticle samples is similar to real atoms. In particular, we may expect that they are capable of accumulating some power inside themselves and to produce (gate out) it during the fusion. It is clear that such macroscopic entities can be formed in various systems with collective interacting particles. (We shall call such media as plasma-like media).
Therefore, it would be interesting to study the possibility of such macroscopic, time-dependent entities formation in weakly dissipative media with different collective interaction between particles. Therefore, the special question within this study is dependence of coherent states on the type of potential – Lennard-Jones, Yukawa, Coulomb potential, etc. To understand this, we are going to consider the influence of the potential type on the formation of coherent states in ideal, weakly collision plasma-like media.
This article discusses an approach to describing the behavior of particles in the approximation of physical kinetics. We examined the formation of coherent structures in the environment, consisting of a Maxwellian ensemble of particles. At the same time, various potentials of interaction between particles are considered.

2. Basic Equations

The medium is assumed to be infinite, i.e., the pure initial value problem is to be studied. Within the framework of a kinetic approach, the governed equations for the one-particle distribution function f = f ( t , r , v ) describing polar molecules dynamics are
f t + v · f F m · f v = I c ,
where F is the force influencing the medium particles of mass m and I c is the collision integral.
To show the basic features of formation of coherent structures in collisionless and weakly collision media we start from the simplest possible form of collision integral
I c = ν ( f f e )
where ν is the collision frequency and f e is some equilibrium distribution to be determined.
Here we suppose that the interaction between particles depends only on distance between particles r r but not their velocities v and v . So one can represent self-consistent force F in terms of the scalar potential Φ = Φ ( t , r ) ,
F = Φ ,
where potential is defined by the relation [4]:
Φ ( t , r ) = K 1 , 2 ( r r ) f ( t , r , v ) d r d v + + K 1 , 2 , 3 ( r r , r r , r r ) f ( t , r , v ) f ( t , r , v ) d r d v d r d v + + Π ( t , r ) ,
where Π ( t , r ) is some known function. In the general case, this function reflects the influence of external superposed factors.
Then, let us discuss the influence of separate terms on the nonlinear properties of the system; however, now for simplicity we restrict our study to the case when
Φ ( t , r ) = K ( r r ) f ( t , r , v ) d r d v + Π ( t , r ) .
The kernel of integral relation is determined by the nature of interaction between particles of the system, for example, that may be dipolar interaction. However, the present approach is valid only when the kernel K ( r , r ) satisfies the relation
γ = + K ( r r ) d r < ,
whereas Coulomb-type potentials cannot be defined by such way since for these fields the relation (5) is violated. In this case we have to use Poison equation.
We are going to construct some partial solution of the issue (1), (4) using [8]:
f ( t , r , v ) = k = 0 L F k ( t , v ) Φ k ( t , r ) ,
where F k ( t , v ) are unknown functions to be determined, and Φ ( t , r ) is determined by relation (4). We assume that L is a finite integer or infinite value. We are going to seek the local equilibrium distribution function in similar form
f e ( t , r , v ) = k = 0 L C k F 0 ( t , v ) Φ k ,
where C k are some constants.
Additionally, we suppose that Φ ( t , r ) satisfies
Φ t + V · Φ = 0 ,
where drift velocity V is a constant.
By applying (6)–(8) to Equation (1), we find
k = 0 L F k t + ν ( F k C k F 0 ) + ( k + 1 ) F k + 1 ( v V ) · Φ + v F k · Φ Φ k = 0 .
If there is no connection between Φ ( t , r ) and F k ( t , v ) , relation (9) has to be satisfied for any power of Φ / q j and Φ (here q 1 = x , q 2 = y , q 3 = z ). So we obtain
v F k + ( k + 1 ) F k + 1 ( v V ) = 0 ,
F k t = ν ( F k C k F 0 ) .
The solution of (11) can be written as
F k ( t , v ) = W k ( v ) exp ( ν t ) + C k F 0 ( v ) ,
where W k ( v ) is some functions. Then substituting (12) into Equation (10) we have
( k + 1 ) ( v V ) W k + 1 e ν t + C k + 1 F 0 + C k v F 0 + e ν t v W k = 0
from which we can obtain
( k + 1 ) C k + 1 F 0 ( v V ) + C k v F 0 = 0 ,
( k + 1 ) W k + 1 ( v V ) + v W k = 0 .
For
C k + 1 = C k ( k + 1 ) T ,
where T > 0 is a constant which can be considered to be temperature; from Equation (14) it follows
F 0 = C 0 exp ( v V ) 2 2 T ,
where C 0 is an arbitrary constant. From (15) we obtain
W k = 1 k ! D k W 0 ( v ) , k = 1 , 2 , , L ,
where
D = V v ( v V ) 2 v .
Finally, the substitution of (16) and (17) into (12) leads to
F ( t , v ) = e ν t k ! D k W 0 ( v ) + C 0 k ! exp ( v V ) 2 2 T
and the distribution function can then be written as
f ( t , r , v ) = k = 0 L e ν t D k W 0 ( v ) + C 0 T k e ( v V ) 2 / 2 T Φ k k ! .
The function W 0 ( v ) can be any function of velocity v such that all moments of the distribution function (19),
I n ( t , r ) = + v n f ( t , r , v ) d v ,
must be finite, i.e.,
I n ( t , r ) < , n = 0 , 1 , 2 , .
and for this kernel K ( r , r ) provide the finite Φ ( t , r ) , i.e.,
K ( r r ) f ( t , r , v ) d r d v < .
Moreover, we have to require
D k = + D k W 0 ( v ) d v < .
It is essential to stress conditions (20)–(22) determine the allowable choices of W 0 ( v ) .
To define the admissible form of Φ ( t , r ) we put (19) into (6). As a result we obtain
Φ ( t , r ) = k = 0 L + K ( r r ) k ! C + D k e ν t Φ k ( t , r ) d r .
Here we used C 0 = C / ( 2 π T ) 3 / 2 .
Relations (23) and (19) with conditions (20)–(22) form the basis for most of the following analysis.

3. Maxwellian Type Distributions

For simplicity, we shall analyze Equation (23) in the case when L and when F 0 belongs the class of Maxwellian functions,
F 0 ( v ) = 1 ( 2 π θ ) 3 / 2 exp ( v V ) 2 2 θ ,
where θ > 0 is a constant which plays role of second temperature. It should be noted that in general case we can suppose T θ . Then from (19) and (24) it follows
f ( t , r , v ) = 1 ( 2 π θ ) 3 / 2 e ν t exp ( v V ) 2 + 2 θ + Φ θ + 1 ( 2 π T ) 3 / 2 exp ( v V ) 2 2 T + Φ θ
In this case, instead of (23) we obtain
Φ ( t , r ) = e ν t + K ( r r ) exp ( Φ / θ ) d r + C + K ( r r ) exp ( Φ / T ) d r
The nonlinear equation (26) admits spatially uniform solution Φ 0 = Φ 0 ( t ) which is determined from
Φ 0 = γ e ν t + Φ 0 / θ + C e Φ 0 / T ,
where γ is determined by (5). There is its nontrivial solution. In particular, if we set θ = T and Φ 0 / θ 1 , one can obtain
Φ 0 = γ e ν t + C 1 C e ν t + C / θ .
Let us consider (28) as a state of local equilibrium, near which the basic equation is linearized, for this we substitute
Φ = Φ 0 ( t ) + δ Φ 0 ( t , r )
into (26). As result we have
δ Φ e Φ 0 / θ θ e ν t + C + K ( r r ) δ Φ ( r ) d r = 0 .
Let us analyze the behavior of Equation (28) with respect to different values of handling parameters ν , θ . In particular, in the limit ν t , Equation (29) admits stationary solution of the form
δ Φ = A e i k r .
Substituting (30) into the (29), we obtain
1 λ + K ( r r ) exp ( i k ( r r ) d r = 0 .
where
λ ( θ , C ) = C exp ( Φ 0 / θ ) θ .
One can estimate the integral
+ K ( r r ) exp ( i k ( r r ) d r = 0 0 π 0 2 π K ( ρ ) exp i k ρ cos ( ϑ ) ρ 2 sin ( ϑ ) d φ d ϑ d ρ
= 4 π k 0 + K ( ρ ) ρ sin ( k ρ ) d ρ .
Thus, we obtain the following dispersion relation
1 4 π λ k 0 + K ( ρ ) ρ sin ( k ρ ) d ρ = 0 .
which defines the existence of spatial wavelet-structures for the real values k. Equation (33) has nontrivial, real root k if and only if the integral in (33) is negative. Such behavior is possible owing to the initial conditions (see Equation (32)) or the kernel of interaction K ( ρ ) .

4. Different Types of Kernels

As an illustration, we now calculate the roots of Equation (33) for the kernel of type
K ( r , r ) = μ r r exp q r r ,
where μ and q are some constants depending on the kind of the interactions which admits elementary analytical consideration and has some physical meaning. Inserting (34) into (33) we obtain
1 4 π λ μ k 2 + q 2 = 0 .
It is evident that this equation has real roots in the case
4 π λ μ q 2 = 0 .
The spatial length of periodic structures
k 1 = ( 4 π λ μ q 2 ) 1 / 2
is defined by the initial conditions over constant λ and the parameters of interaction q and μ . This example shows that for some initial conditions and kernels one can expect the formation of spatially nonuniform structures from initial uniform states [9]. The relation (33) or (35) determines only value of k but its direction is unknown. Using (37) we can determine the size of the vortex depending on the initial parameters of the system.
As another result of the application of this method, we give an example of the formation of coherent structures in the "billiard balls" model. In that case, the kernel of interaction is shown by
K ( r , r ) = μ δ ( r r a )
where μ is some constant depending of the kind of interactions and | a | is the radius of every particle of the media. Inserting (38) into dispersion relation (33) we obtain
1 4 π λ μ a sin ( a k ) k = 0 .
We can transform Equation (39) into a more convenient form
k = 4 π λ μ a sin ( a k ) .
Equation (40) has some sets of solutions. This shows us that if we imagine that particles of a medium interact only through collisions, coherent vortexes still appear. They can have different sizes that correspond to different positive roots of the Equation (40). Additionally, we add that if
q > > 1 ,
where q is the constant included in (34), the kernel (34) is greatly approximated by a delta-function such as (38). That is, the Yukawa kernel model under the condition of a low intensity of particle interaction (41) could be approximated by the model of billiard balls.

5. Discussion and Conclusions

To describe the formation of kinetic coherent structures in weakly collisional medium with collective interactions we used one-dimensional, time-dependent nonlinear Vlasov equations with simplest collision integral in the Bhatnagar-Gross-Krook form. Additionally, in our model of collisional medium we have used some physically admissible forms of the interparticle interaction potentials providing the finiteness of the energy and momentum of the systems under consideration. As such potentials, we select the Yukawa potential, the δ -potential, and the soliton potential, which may describe electron holes in a plasma [10].
We applied the method of constructing time-dependent solutions of the kinetics equations based on the expansion of the distribution function as a series in positive powers of the interparticle interaction potentials [8]. The key point of this method consists in the representation of nonstationary solution and local equilibrium distribution function in similar form (see Equations (6) and (7)). As a result, we managed to show that a Maxwellian time-dependent distribution function (19) appears as a natural solution of the considered initial-value problem (see Equations (10) and (11)) when such Maxwellian time-dependent distribution can be formed for any interparticle interaction potential.
In this class of Maxwellian distributions for the considered potentials, we obtained a dispersion relation Equation (33) which makes it possible to estimate the size and type of the forming structures in locally equilibrium media. Analyzing dispersion relation (33), we presented coherent structures as wave solutions of the equation, which describes variations with respect to the equilibrium states of system. Such an approximation is valid if the spatial scales of the resulting structures are small with respect to the characteristic scale of the system [11].
Analyzing relation (33), we presented coherent structures as wave solutions of equation, which describes variations with respect to the equilibrium states of system [12]. Taking into account such high sensitivity of our pattern to perturbations in the initial conditions, one can come to the conclusion that there exists a small interval in the system parameters in which it is possible to see the effect described in reality. Our consideration was limited to the case of Maxwellian distributions. In this regard it would be interesting to study the Lorenz-Fermi distributions which may be more relevant for some natural cases.

Author Contributions

Conceptualization, A.K., V.B.; methodology, A.K.; validation, A.K., V.B.; investigation, V.B., writing review and editing, V.B.; supervision, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research has been supported by the grant “Governing for the processes of structure formation in high-molecular compounds under non-equilibrium conditions” of Plekhanov Russian University of Economics.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Hans Schamel for his support in developing the theory. Specifically, for providing information and access to existing works on the topic, as well as for verifying the truth of the methods presented in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kadomtsev, B.B. Role of Weak Energy Fluxes in the Formation of Coherent Structures in Macromolecular Media. Uspekhi Fiz. Nauk. 1994, 164, 449. [Google Scholar] [CrossRef]
  2. Nicolis, G.; Prigogine, I. Self-Organization in Nonequilibrium Systems; Wiley: New York, NY, USA, 1977. [Google Scholar]
  3. Kuramoto, Y. Chemical Oscillations, Waves and Turbulence; Springer: Berlin, Germany, 1984. [Google Scholar]
  4. Vlasov, A.A. Nonlocal Statistical Physics [in Russian]; Nauka: Moscow, Russia, 1979. [Google Scholar]
  5. Davidson, R. Methods in Nonlinear Plasma Theory; Plasmas, Academic: New York, NY, USA, 1972. [Google Scholar]
  6. Ablowitz, M.J.; Segur, H. Solitons and the Inverse Scattering Transform; SIAM: Philadelphia, PA, USA, 1981. [Google Scholar]
  7. Born, M.; Infeld, I. Foundations of the new field theory. Proc. R. Soc. 1934, A144, 425. [Google Scholar] [CrossRef]
  8. Karimov, A.R.; Lewis, H.R. Nonlinear solutions of a Maxwellian type for the Vlasov–Poisson equations. Phys. Plasmas 2001, 6, 759. [Google Scholar] [CrossRef]
  9. Schamel, H. Hole equilibria in Vlasov–Poisson systems: A challenge to wave theories of ideal plasmas. Phys. Plasmas 2000, 7, 4831. [Google Scholar] [CrossRef]
  10. Luque, A.; Schamel, H. Electrostatic trapping as a key to the dynamics of plasmas, fluids and other collective systems. Phys. Rep. 2005, 415, 261–359. [Google Scholar] [CrossRef]
  11. Karimov, A.R.; Yu, M.Y.; Stenflo, L. Properties and evolution of anisotropic structures in collisionless plasmas. J. Plasma Phys. 2016, 82, 905820502. [Google Scholar] [CrossRef] [Green Version]
  12. Schamel, H.; Mandal, D.; Sharma, D. Diversity of solitary electron holes operating with non-perturbative trapping. Phys. Plasmas 2020, 27, 062302. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karimov, A.; Bogdanov, V. Formation of Kinetics Coherent Structures in Weakly Collisional Media. Plasma 2021, 4, 359-365. https://0-doi-org.brum.beds.ac.uk/10.3390/plasma4020024

AMA Style

Karimov A, Bogdanov V. Formation of Kinetics Coherent Structures in Weakly Collisional Media. Plasma. 2021; 4(2):359-365. https://0-doi-org.brum.beds.ac.uk/10.3390/plasma4020024

Chicago/Turabian Style

Karimov, Alexander, and Vladislav Bogdanov. 2021. "Formation of Kinetics Coherent Structures in Weakly Collisional Media" Plasma 4, no. 2: 359-365. https://0-doi-org.brum.beds.ac.uk/10.3390/plasma4020024

Article Metrics

Back to TopTop