The computation of thermalt properties of dusty plasmas is substantial task in the area of science and technology. The thermal conductivity (λ) has been computed by applying polarization effect through molecular dynamics (MD) simulations of two dimensional (2D) strongly coupled complex dusty plasmas (SCCDPs). The effects of polarization on thermal conductivity have been measured for a wide range of Coulomb coupling (Γ) and Debye screening (κ) parameters using homogeneous non-equilibrium molecular dynamics (HNEMD) method for suitable system sizes. The HNEMD simulation method is employed at constant external force field strength (F*) and varying polarization effects. The algorithm provides precise results with rapid convergence and minute dimension effects. The outcomes have been compared with earlier available simulation results of molecular dynamics, theoretical predictions and experimental results of complex dusty plasma liquids. The calculations show that the kinetic energy of SCCDPS depends upon the system temperature (≡ 1/Г) and it is independent of higher screening parameter. Furthermore, it has shown that the presented HNEMD method has more reliable results than those obtained through earlier known numerical methods.
- Plasma thermal conductivity
- complex dusty plasma
- Homogenous non-equilibrium molecular dynamics
- force field strength
- system size
- plasma parameters etc.
Recently, thermophysical properties of complex materials are a major concern in the field of science and engineering. The term thermophysical properties used to pass on both thermodynamic and transport properties. Experimental or theoretical methods to study properties of fluids depend on microscopic and macroscopic categories [1, 2, 3, 4]. The conventional macroscopic measurements depend on the state of stress, temperature, and density. Thermodynamic properties are defined by the equilibrium conditions of the system which consist of temperature, heat capacity, entropy, pressure, internal energy, enthalpy, and density, whereas the transport properties comprise thermal conductivity, diffusion viscosity, and waves with their instabilities. For further explanation of the process in detail for these systems, data that is applicable to thermodynamics, transport, optics, transmission, light, and other features are required for non-ideal plasma [5, 6, 7, 8, 9]. In this regard, various opinions regarding computer research methods including theoretical and numerical performance have greatly improved for non-ideal Plasma . Determination for some reason, thermal conductivity is also a big problem for thermophysical researchers. Developmental aspects of heat transport in micron and nanoscale materials have shifted to the domain of technical issues as there are other areas, such as phonon heat transfer in semiconductor superlattices, which have received widespread attention from researchers. To study the internal energy of particles, their momentum, and heat transfer thus remains a crucial task. Therefore, thermal management, strategies sustainable high performance, reliability, and service life are main purposes. One such strategy is to develop new therapeutic materials based on dusty plasma that are more effective. Regulation with approval became a significant issue in modern technology . Yet similar interests are present in plasma fusion, and it can be productive radiotoxic dust in plasma-wall reactions. In many ways, this chapter provides an update literature survey on thermal transport as well as heat flow strategies to determine thermal behaviors in two-dimensional (2D) complex liquids. The coefficients were computed through the Green Kubo (GK) equilibrium molecular dynamics (EMD) simulations by Salin and Caillol  and variance procedure (VP) estimation used by the Faussurier and Murillo . Donkó and and Hartmann employed the inhomogenous non equilibrium MD (InHNEMD) method to investigate the transport and thermal conductivity . Very recently, a homogeneous NEMD (HNEMD) and homogenous perturbed MD (HPMD) schemes are introduced by Shahzad and He (current authors) for strongly coupled complex dusty plasmas (SCCDPs) to compute the thermal transport and behaviors of SCCDPs [15, 16, 17]. For the computation of transport properties, in particular, numerical models are proposed in interest to investigate thermal behavior over a suitable range of system temperature and density values (Γ, κ). Complex fluids (dusty plasma fluids) have been used for many purposes, like power generation, semiconductors industry, cosmetics, paper industry, etc. .
As we all know that 99% of matter exist in space is plasma and it is called forth state of matter. Basically plasma occurs in electrified gas form, where atoms dissociated into electrons and positive ions. It is form of matter in different areas of physics such as technical plasma, terrestrial plasma and in astrophysics. Plasma is produced artificially in laboratory used in many technical purposes likely in fluorescent lights, display, fusion energy research and other more. Term “Plasma” first time used by Irving Langmuir , who is an American physicist and defined plasma as “plasma is quasi-neutral gas of charged particles which exhibit collective behavior”. Quasi-neutral means that gas becomes electrically neutral when number of ions equal to number of electrons (
In 1922, American scientist Irving Langmuir was the only one person who defined plasma for the first time. In 1930, the study of plasma physics was started by some scholars; they are inspired by some particle problems. In 1940, hydromagnetic waves were advanced by Hanes Alfven  and these waves are called Alfven waves. Furthermore, he described that these waves would be used for the study of astrophysical plasma. At the start of 1950, the research on magnetic fusion energy was started at the same time in Soviet, Britain and USA. In 1958, the research on magnetic fusion energy was considered the branch of thermonuclear power. Primarily, this research was carried out as confidential but after the realization that controlled fusion research was not liked by military and therefore this research was publicized by above said three countries. Due to the reason, other countries may participate in fusion research based on plasma physics. At the end of 1960, plasma is created with different plasma parameters by Russian Tokomak configuration. In 1970 and 1980, various advanced tokomaks were built and approved the performance of tokomak. Moreover, fusion break almost achieved in tokomak and in 1990, the research on dusty plasma physics had begun. The dusty plasma is defined as “when charged particles absorbed in plasma, becomes four components plasma containing electrons, ions, neutral and dust particles” and dust particles alter the properties of plasma which is called “Dusty Plasma” .
1.2 Types of plasma
Plasma has complex characteristics and properties, characterized through temperature of electron and ion, density and degree of ionization. (i)
1.3 Classification of dusty plasmas
Dusty plasma is characterized by an important parameter, coulomb coupling parameter Г. The Coulomb coupling parameter is explained as, consider there are two dust particles, having same charge and separated by distance ‘
1.4 Complex (dusty) plasma and applications
Dusty plasma is generally electron ion plasma containing additional charged particulates. This charged component is sometimes termed as dust particle with size of micron. The properties of dusty plasma become more complex when charged particle immersed in plasma. Due to this reason such plasma is called dusty plasma and dusty plasma is also called complex plasma. Dust particles may be made of ice particles or it may be metallic. Dust particles are heavier than ions and their size ranging from few millimeters to nanometer. When dust particle coexists with plasma (electron, ions, neutral and dust particle) it becomes dusty plasma. Dust particle exists in different shapes and size and it presents in entire universe and also in atmosphere. Usually it is solid form but also exists in liquid and gaseous form. Dust particle can be charge by the flow of electrons and ions. Charged dust particle is affected by electric and magnetic field and their electric potential varies from 1 to 10 V. Dust particle can be grown in laboratory. Dusty plasma has attracted attention of many researchers Transport properties of dusty plasma has played a very important role in the field of science and technology. Mostly the plasma exists in universe is dusty plasma. Dusty plasmas exist in atmosphere of stars, solar wind, sun, galaxies, planetary rings, cosmic radiation, magneto and ionosphere of earth.
Human life is influenced by plasma science. It plays a very significant role in laser developments of fusion energy, sterilizing of medical instruments, plasma processing, intense particle beam, high power energy sources, lightening, high power radiation sources and development of fusion energy controlling. Plasma governs diverse important devices and technological applications. Plasma processing technologies are one of the most important technologies. Plasma processing technologies are playing important role in advance modern technologies of superconductor film growth and diamond film. In addition, the practical application of plasma physics involves the treatment of materials by means plasma technologies. The ionization of system are used to produced particular physical characteristics of plasma, which involve three types of processes, Creation of new materials, Destruction of toxic materials and Superficial modification of existing materials. For industrial process, plasma technology uses two different types of plasma, the cold plasma and thermal plasma. The first type of plasma is cold plasma. Properties of cold plasma are described by electron temperature because electron temperature is greater than ion temperature. The surfaces modification is produced due to plasma particles interact with material, as a result different functional properties of materials are achieved. Cold plasma is produced in vacuum with microwave, dc source or low power rf. The second type of plasma is thermal plasma which is produced at high pressure by radio frequency, microwave source or direct or alternating current. Mostly, thermal plasma is used to devastate toxic materials. Furthermore, plasma has become one of the fast growing research fields which have attracted many researchers. Plasma has advanced applications in the field of industry, textile, plasma chemistry, fusion devices, environmental safety and printing technology and as well as in medical field. In the past decade, plasma physics has become fast developing research in medical field due to increase the atmospheric pressure of plasma sources. Plasma used in medicine has considered the latest developing novel research field with the connection of life science and plasma physics. Moreover, in past ten to fifteen years, World wild research group has set their attention on biological materials with cold atmospheric plasma interactions. Plasma used the field of life sciences in decontamination, in therapeutic medicine and in medical implant technology. The atmospheric pressure of plasma has used to reduce the efficiency of contaminations of food containers and food products. Feed gas humidity is used to adjust the level of contamination. Furthermore, plasma created in polymer tube, used in endoscopes [15, 16, 17, 18, 19]. Operating tools such as bone saw blade, neurosurgical and endoscope are sterilized before starting the surgery or dental treatments. Plasma plays a very dominant role in diagnostic system, treatments and in medical instruments. For decontamination of germs and sterilization operation tools, the non-equilibrium discharge plasma is used, is not dangerous for environment and patient as well .
2. Molecular dynamics simulations
MD simulation is a powerful technique that can be used to solve many physical problems in atomic material research. MD simulation is handled normally all microscopic information and molecular methods have proven to be the product of applied research. Plasmas and complex liquids have various uses, ranging from semiconductor chips, colloids, thin films, and electrochemistry to biochemical films and other important areas where structures play an important role. MD simulation plays an important role in all the advanced sectors, such as textile science, engineering, physics, plasma physics, astronomy, life sciences or organic sciences, and the chemical industry. Computer simulation has become increasingly important in detecting complex motion systems. Using faster and more sophisticated computer systems, it can be studied the habitat, composition, and behavior of large complex systems. In the 1950s, Alder , Wainwright and Rahman  used the first MD simulations for liquid argon [1, 2, 3, 4, 5, 6, 7, 8, 9, 10], their references herein]. MD simulation has two basic kinds which rely on the properties so far which we can be going to calculate: one is EMD simulations (EMDS) and the other one is NEMD simulations (NEMDS). In the present work, NEMDS is applied to investigate the thermal conductivity of complex plasma at different dusty plasma parameters.
2.1 Numerical model and algorithm
NEMD simulation is used to detect the dust trajectory of an interacting system using Yukawa forces between dust particles. In present case, the HNEMD simulation (HNEMDS) method is used to calculate the thermal flow of complex (dust) plasma formed using Yukawa interaction taking in to account the charged particles with polarization effects and it is given in the form :
The first term in Eq. (1) provides the screened charge–charge interaction (form of Yukawa interaction) and the second term gives the screened dipole–dipole interaction. Yukawa potential model of dust particles interaction can be established to take into consideration the polarization effect, the temperature and the screening effects. Here
In Eq. (4), is the Yukawa interaction force acting on particle
When an external force field is selected parallel to the z-axis
3. Simulation results and discussion
In this section, we have discussed the preliminary results obtained through HNMED simulation for Coulomb coupling parameters of Г (= 10, 100), polarization values Гd = (0, 1, 10, 20, 50 and 100) and Debye screening (
There are different conditions to improve the efficient results of thermal conductivity under polarization effects. These conditions include the system size (
The polarized thermal conductivity of SCCDPs stated here may be scaled as λ0 = λ/
Figures 1–3 illustrate the normalized polarized thermal conductivity (plasma frequency,
Figure 4 shows comparisons with earlier available 2D and 3D numerical data of thermal conductivity with setting
The HNEMD scheme is used for the investigation of thermal conductivity under the influence of varying polarization values for various screening lengths
We are grateful to the National High Performance Computing Center of Xian Jiaotong University and National Advanced Computing Center, National Center of Physics (NCP), Pakistan for allocating computer time to test and run our MD code.
|Strongly coupled complex dusty plasmas||(SCCDP)|
|Homogeneous non-equilibrium molecular dynamics||(HNEMD)|
|Debye screening length||(κ)|
|External force field strength||(F*)|
|Homogenous non-equilibrium molecular dynamics||(HNEMD)|
|Non-equilibrium molecular dynamics||NEMD|
|Inhomogenous non-equilibrium molecular dynamics||(In HNEMD)|
|Strongly coupled plasma||(SCP)|
|Equilibrium molecular dynamics||(EMD)|
|Normalized thermal conductivity||(λ0)|
|Number of Particles||(N)|