• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effects of inter-electrode insertion on the performance and thermal flow fields of a hollow-electrode plasma torch

    2020-01-10 07:40:40InMokYANGJunSeokNAMMinKyuCHOIJunHoSEO3andShiYoungYANG
    Plasma Science and Technology 2020年1期

    In-Mok YANG, Jun-Seok NAM, Min-Kyu CHOI, Jun-Ho SEO,3and Shi-Young YANG

    1 Department of Applied Plasma Engineering, Chonbuk National University, Jeonju 54896, Republic of Korea

    2 Graduate School of Flexible & Printable Electronics, Chonbuk National University, Jeonju 54896,

    Republic of Korea

    Abstract The effects of inter-electrode insertion on the performance of a hollow-electrode plasma torch have been investigated by numerical analysis. Simulation results revealed that when inter-electrodes are inserted,the arc voltages and plasma powers increase due to the increase in the arc length.In addition,it was predicted that thermal efficiency can be improved with the increase in plasma power by injecting plasma gases through the gaps between inter-electrodes. These unique effects of inter-electrode insertion are a result of the plasma temperatures adjusting themselves to increase arc voltages when the arc column is contracted radially by increasing gas-flow rate or decreasing inter-electrode diameter.

    Keywords: hollow-electrode plasma torch, inter-electrodes, numerical analysis, torch performance(Some figures may appear in colour only in the online journal)

    1. Introduction

    In recent decades, hollow-electrode plasma torches have been widely used as a high powered heat source in the thermal plasma treatment of various wastes, such as PCBs (polychlorinated biphenyls)containing oil[1,2],municipal solid waste[3,4],and low-level radioactive wastes [5, 6]. By producing ionized hot flames with high temperatures of >5000 K even in oxidizing gases[7],this type of plasma torch is advantageous for not only gasification of combustible wastes by pyrolysis or oxidation,but also vitrification of noncombustibles by plasma melting. These features of a hollow-electrode plasma torch that can work with oxidizing gases come from the cold-well type of cathode primarily made of copper or copper alloys [7, 8]. However, these relatively low melting-point metals require not only current limitation(<1-2 kA)but also arc-root movement using a vortex gas injection [7, 8] in order to secure an electrode lifetime for practical applications. As a result, hollow-electrode plasma torches show the unique dynamic behaviors of arc voltages caused by the arc-root movement on the surface of the cold-well type cathode, resulting in the continuous fluctuation of the plasma flame. In addition, high gas flow rates for vortex gas injection often make it difficult to produce plasma flames with a specific enthalpy higher than 10 MJ kg?1[9].

    In order to obtain stable plasma flames with high specific enthalpy, inter-electrodes (hereinafter referred to as IEs) can be inserted between the cold-well type cathode and the cylindrical anode of a conventional hollow-electrode plasma torch,as introduced in the patent of SKF Co.[10].First of all,the insertion of electrically insulated IEs can allow highvoltage operation of the torch by forcing the arc column to extend as long as the IE length.The high-voltage operation is advantageous for scaling up the operation power level under the current limitation (≤1-2 kA) of copper or copper-alloy electrodes, permitting the extension of electrode service lifetime. In addition, the elongated arc column can alsocontribute to increasing the specific enthalpy of the plasma flame by heating up the injected gases for a longer time before they exit the torch[10].Accordingly,the insertion of IEs can bring a high-voltage/high-power operation mode together with high thermal efficiency to the conventional hollowelectrode plasma torch, allowing it to have a long electrode lifetime as well as a high specific enthalpy plasma jet.Although these advantages were revealed in the patent of SKF Co. [10] in the 1980s, there have been few papers that show the effects of IE insertion through numerical simulations.

    Table 1.Design parameters for numerical analysis of hollowelectrode plasma torches with IEs.

    In this work, we present numerical results about these effects of IE insertion on the torch performance and thermal flow fields of a hollow-electrode plasma torch. By predicting the torch performance according to the IE insertion, our numerical work can be expected to help in designing and improving this type of plasma torch for actual applications,such as the thermal plasma treatment of hazardous wastes.

    2. Numerical modeling

    2.1. Torch geometry, numerical methods, and basic assumptions

    Figures 1(a)and(b)show the electrode configuration and the computational domain of a hollow-electrode plasma torch with IEs employed in the numerical work.In these figures,the presented plasma torch has three IEs, designated as first IE,second IE, and third IE, between the hollow cathode and cylindrical anode. In order to investigate the effects of IE insertion, however, numerical analyses were conducted by changing the number of IEs from one to three.

    In figure 1, the plasma-forming gases are injected into the plasma torch through the gaps between electrodes.Considering that the gaps are distributed with the addition of IEs, plasma-forming gases were labeled after the left electrode name of the gap. For example, in figure 1(a), cathode gas with a flow rate of Gc is introduced into the torch through the gap between the cathode and the first IE. In addition, the cathode, first IE, second IE, and third IE andanode are assumed to be electrically insulated from each other. By inserting these electrically insulated IEs, the arc column can be elongated forcibly as shown in figure 1(a),illustrating the electrical connections together with the arc column between the front anode and rear cathode. In the actual system, as indicated in this electrical connection, the first IE can be used as an intermediate electrode for easy ignition [10]. For example, the arc can ignite between the cathode and the first IE when the switch is closed in figure 1(a),and the arc jet can be formed.Depending on the arc current and gas injection conditions,the generated arc jet can touch the anode surface,allowing a small amount of arc current to flow through the anode. Then, for the two arc currents flowing through the first IE and anode, a power control system (not seen in figure 1(a)) can be designed to raise up the ratio of arc current through the anode by gradually opening the switch connected to the first IE. Finally,if the power supply opens the switch completely, the long arc column can be obtained between the rear cathode and front anode.

    Table 2.Operation conditions for numerical analysis of hollowelectrode plasma torches with IEs.

    The details of torch dimensions and operating conditions illustrated in figures 1(a) and (b) are summarized in tables 1 and 2, respectively. The computational domain was divided by rectangular grids. All governing equations are discretized by the finite-volume method (FVM) and calculated using the SIMPLE algorithm [11]. Finally, the following assumptions were used for the present numerical modeling work:

    (1) Steady, 2D axisymmetric thermal plasma.

    (2) Negligible displacement current in the plasma.

    (3) Optically thin plasma.

    (4) Local thermodynamic equilibrium (LTE).

    (5) Incompressible nitrogen plasma at atmospheric pressure.

    2.2. Governing equations

    As a thermal fluid flow,arc plasma inside the hollow-electrode plasma torch can be described by the MHD (magnetohydrodynamic) equations combining the well-known conservation equations for mass, momentum, and energy with the electric potential equation inside the torch. In addition, the turbulent effects in the arc plasma can be included by selecting the appropriate turbulence model and adding the related equations to the MHD equation. In this work, we used a standard K-ε turbulence model, which was reported to produce reasonable results when simulating the arc plasma generated in the cylindrical electrodes [12]. These governing equations are expressed in the steady state, axisymmetric cylindrical coordinates as follows:

    (1) Conservation of mass

    (2) Conservation of momentum

    · Axial component

    · Radial component

    · Azimuthal component

    (1) Conservation of energy

    (2) Electrical potential equation

    (3) Standard K-ε turbulence model

    In the above equations, u, v, and w are the axial, radial, and tangential components of the velocity vector of the arc plasma.The terms ofρ,p,and μ represent the mass density,pressure and viscosity, respectively. The body forces Fzand Frin the momentum conservation equations (equations (2) and (3))indicate the axial and radial components of the Lorentz force,respectively, which are exerted on the arc currents by the selfinduced magnetic field. Gravity is neglected in this work. In equation (5), k, Cp, and h indicate the thermal conductivity,specific heat at constant pressure,and specific enthalpy of the arc plasma,respectively.The Joule heating and volumetric radiation heat loss in equation (5) are represented in terms of P and R0,respectively, where the radiation loss term R0is taken into account by using the net emission coefficient for optically thin nitrogen plasma. Additionally, Poisson's equation for the electrostatic potential φ is used to find the arc current density vector, j, in the arc column. In equation (6), σ is the electrical conductivity, and the axial and radial components of j can be expressed as the derivative forms of the electrostatic potential φ as:

    The self-induced magnetic field B by arc current is also obtained from Ampere's law, ?× B =μ0j.From the calculated results of the arc current and the magnetic field, the source terms associated with the Lorentz force and Joule heating can be obtained according to the following equations:

    In equations (7) and (8) for the description of a standard K-ε turbulence model, the viscosity term is modified by adding the turbulent viscosity, μt, to the molecular one, μ, defined as follows:

    The typical values are allocated to the turbulence constants as Cμ= 0.09, Ck= 1.00, Cε= 1.30, C1= 1.44, and C2= 1.92[13], and the turbulent generation term of G is expressed as:

    The thermodynamic and transport properties of nitrogen used for the numerical calculations of these equations can be found in [14].

    2.3. Boundary conditions

    Table 3 summarizes the boundary conditions used in this work. In a suggested type of plasma torch, plasma-forming gases are normally injected with a swirl to stabilize the arc column along the centerline [10, 15].

    In order to describe this swirling injection of plasmaforming gases, we used the swirl number, Sw, defined as the ratio of azimuthal velocity to the radial one, i.e., Sw= w v,/for the boundary conditions at gas inlets as shown in table 3.In this expression of swirl number,the radial velocity v can be calculated from the gas flow rates distributed to each gas inlet as listed in table 2. In addition, the turbulent variables at the gas inlets were calculated by using the following equations,which are widely employed for the numerical modelling of arc plasmas [12, 13].

    Table 3 also shows that there is no radial change of velocity, temperature, and turbulent variables caused by the axial symmetry along the central axis(r = 0),whereas,at the torch exit, no axial gradients for those parameters exist,assuming that the arc plasma flow is fully developed. Along the inner surface of the electrodes including the IEs,a no-slip condition and the wall function were used for boundary conditions of flow parameters(u,v,w)and turbulent variables,respectively,together with the constant temperatures as listed in table 3.

    For the boundary conditions of the electric potential equation, we used the usual assumptions widely employed in DC arc modeling [12, 15, 16]. For example, the arc current density jsis assumed to have the following exponential distribution at the cathode spot [15, 16].

    Here, the exponential coefficient b is chosen to satisfy the

    Table 3.Boundary conditions for temperatures, velocities, and turbulent variables.

    following equation, indicating that

    the arc current I0supplied from the power supply should be equal to the value obtained by integrating the current density jsover the cathode spot area, S. The location of the cathode spot was identified on the surface of the rear electrode as being 25 mm away from the closed end of the rear electrode by cold-gas analysis [16]. At the anode surface of the front electrode in figure 1,the electrostatic potential φ was assumed to be zero,taking into account that the electrical conductivity of the anode (normally Cu) is very high. In the modelling of the conventional DC plasma torch, the position of the anode spot has usually been determined by Steenbeck's minimum principle [16, 17]. As is well known, Steenbeck's minimum principle postulates that the arc length is determined naturally so as to minimize the arc voltage for a given current and boundary conditions. In this work dealing with the IE insertion,however,it is inappropriate to use Steenbeck's minimum principle to determine the anode spot position because the arc length is forced to be elongated by the electrically insulated IEs.Taking into account this difficulty in determining the arc lengths,we did not consider the formation of the electric spot on the anode surface in this work. In other words, Poisson's equation for φ(equation(6))can be solved with the boundary conditions at the anode for φ = 0 and the cathode for equation (15). According to equations (6) and (9), the calculated results for φ produce not only arc current density vector, j, but also the Joule heating, P, in the computational domain. From the computed Joule heating, plasma input power, P0(kW), can be calculated at a given arc current and other simulation conditions, such as gas flow rates and torch dimensions. In addition to plasma input power, design criteria, such as average specific enthalpy at torch exit, H0(kJ kg?1)and torch efficiency,η(%),can be also obtained.Here,H0and η are defined as the ratio of total enthalpy at torch exit to the mass flow rates of nitrogen and P0, respectively.Finally, no arc current flows through the gas inlets and the IEs, since they are electrically insulated.

    3. Results and discussion

    3.1. Effects of IE insertion

    Figures 2(a)-(c)show the calculated temperature fields of N2arc plasma simulated in the plasma torches with three different lengths of IEs, L0= 250 mm, 370 mm and 490 mm,corresponding to the number of the inserted IEs, 1, 2 and 3,respectively. These numerical results were obtained for the operating conditions of I0= 600 A and the gas-flow rates listed in table 2 with Sw= 3.0. In this figure, first, hightemperature regions of ≥10 000 K were observed to extend in the axial direction proportional to the number of the cylindrical IE. Considering that the nitrogen is ionized in earnest higher than 10 000 K, the axial extension of the high-temperature regions can be regarded as an extension of the arc column, resulting in the elevation of arc voltages by the insertion of the IE. This increase of arc voltages can be confirmed in table 4, which presents the main numerical results for three hollow-electrode plasma torches used to obtain figure 4.From table 4,one can see that the arc voltages are increasing at the rate of ~1.4 V mm?1with the increase of the total IE length L0.Because of this increase of arc voltages,it can also be observed that the plasma powers are increasing proportionally to L0. With the increase of plasma powers, H0and η were calculated to be elevated, as presented in table 4.These improvements in thermal efficiency and average specific enthalpy at torch exit can be attributed to the additional N2injected through the gaps between the existing IEs and the newly inserted ones.

    As presented in table 2 and figure 1,the total flow rate of plasma-forming gas was designed to be increased with the increase in the number of IEs by injecting additional N2through the gaps between the existing IEs and the newly inserted ones.Consequently,the raising of arc voltages by the insertion of IEs can be used to ionize this additionally injected N2, which increases not only the plasma power but also the efficient heating of plasma-forming gases shared by the gaps between IEs.

    3.2. Effects of gas-flow rates

    Figures 3(a)-(c)show the calculated temperature fields of N2arc plasmas inside the plasma torch with three IEs(L0= 490 mm)for three different gas flow rates of Gin= 750 slpm,900 slpm,and 1050 slpm with Sw= 3.0.Each value of Ginrepresents the summation of the flow rates of gas injected through the gaps between electrodes, as listed in table 2. In addition, table 5 shows the main numerical results for three cases of figures 3(a)-(c). From the comparison between figures 3(a)-(c), when Gin= 750 slpm, the temperature region of ≥10 000 K corresponding to arc column is diminished in the section of the third IE(0.1 m ≤ z ≤ 0.2 m).With the increase of total gas-flow rate Gin, however, the arc column was heating up to ≥10 000 K in the section of third IE,as shown in figures 3(b) and (c). This heating up and of arc column with the increase of Ginis to balance the generation and removal of not only heat but also charged particles in the arc column [18].

    Table 4.Calculation results for three different IE lengths (L0) and gas-flow rates (Gin) at a fixed arc current I0 = 600 A and the swirl number of Sw = 3.0.The gas flow rate Gin depends on L0 because of the torch structure, as listed in table 2.

    As explained by Pfender et al for example, if the gasflow rates of G1,G2,and G3are increased in figure 1,the arc column surface is cooled even more in the section of the third IE, increasing the loss of charged particles together with the heat loss toward the wall of the third IE. In order to compensate for these increasing losses of charged particles and energy, the ionization rate and heat generation in the arc column should be increased by raising the arc voltages [18, 19].

    Figure 4 shows how the plasma temperatures adjust themselves to meet these requirements for increasing not only the ionization rate but also heat generation by raising the arc voltage. In this figure, which shows the radial temperature profiles of arc plasma simulated for Gin= 750 slpm and Gin= 1050 slpm, one can see that the high-temperature region of ≥8000 K is expanded widely for Gin= 750 slpm,producing a relatively gradual temperature gradient. However, for Gin= 1050 slpm, the same region of ≥8000 K is contracted toward the center line, where the central temperatures are increased, making the temperature gradient steeper. In other words, if the arc column is cooled down because of the increasing gas-flow rate,the arc column can be contracted, as shown in figure 4, and then the electrical resistance can be increased. Consequently, the arc voltages will be elevated as represented in table 5, increasing the ionization rate as well as the Joule heat generation in the contracted arc column. Thanks to the increase of the ionization rate, first, the loss of charged particles can be compensated for, keeping the self-sustaining discharge condition.Next, with the increase of Joule heat generation, the core of the arc column can be heated up to temperatures higher than 10 000 K, as shown in figures 3 and 4.

    Table 5.Calculation results for three different gas-flow rates (Gin)with Sw = 3.0 at a fixed arc current of I0 = 600 A and an IE length of L0 = 490 mm.

    Table 6.Calculation results for plasma torches with different IE diameters (D3 = 37 mm and 25 mm) at the fixed gas-flow rate and arc current of Gin = 1050 slpm with Sw = 3.0 and I0 = 600 A,respectively.

    Combined with the contraction of the arc column, this increase of central temperatures brings a steep gradient to the temperature fields, allowing for the increase of radial heat transfer in accordance with the greater heat loss.Since the arc was optically thin,radiation heat loss hardly contributes to the radial heat transfer even if, above 8000 K, radiation may be strong. Thus, the radial heat transfer is dominated by the radial diffusion, i.e., the radial temperature gradient, which decreases due to the association of N atoms in the region below 8000 K [20, 21]. As a result, despite a relatively large increase of plasma power, the thermal efficiencyη is slightly changed within 1%for the gas flow rates of Ginbetween 750slpm and 1050 slpm,as listed in table 5,and then the average specific enthalpy at torch exit H0decreases proportionally to the increase in gas flow rate Gin.

    Table 7.Calculation results for plasma torch with L0 = 490 mm obtained by varying the arc currents from 400 A to 600 A at a fixed gas-flow rate, Gin = 1050 slpm with Sw = 3.0.

    3.3. Effects of IE diameters

    Figure 5 compares the temperature fields of N2arc plasmas inside the plasma torch with two different diameters of the second and third IEs of D3= 25 mm and 37 mm. In these numerical results, the arc currents and total gas-flow rates were fixed at I0= 600 A and Gin= 1050 slpm with Sw= 3.0, respectively.

    In figure 5, the temperature region of ≥8000 K is expanded more for D3= 37 mm than for D3= 25 mm,although the high-temperature region of ≥10 000 K is reduced. This change of temperature fields depending on IE diameter can also be explained by the behavior of plasma temperatures to balance generation and removal of not only heat but also charged particles in the arc column,as shown in figure 6.In this figure,comparing the radial profiles of plasma temperatures for D3= 37 mm and D3= 25 mm,for example,the central temperatures are decreasing but the temperature region of ≥7000 K is widely expanded with the increase of IE diameter from 25 mm to 37 mm.

    By changing the radial profiles of plasma temperatures as shown in figure 6,the arc column can minimize losses of both heat and charged particles, then balance their generation and removal, corresponding to the increase of IE diameter. In addition, although the temperature is as low as 7000 K, the radial expansion of this high-temperature region certainly contributes to the decrease of arc resistance. As discussed in figure 4, accordingly, arc voltage is decreasing in the arc column for D3= 37 mm,and consequently the plasma power can be lowered compared with the one for D3= 25 mm as listed in table 6, in which it should be checked that thermal efficiency is decreasing for D3= 37 mm. As shown in figure 6, the plasma temperatures higher than 6000 K are rapidly falling to the coolant temperature of 300 K near the wall regardless of IE diameter.However,the internal surfaces in the second and third IE sections are increasing with the increase in IE diameter, D3. Accordingly, a relatively large portion of heat is transferred to the wall for D3= 37 mm,despite the decrease in plasma power,leading to the decrease in thermal efficiency.

    3.4. Effects of arc current

    Table 7 presents the numerical results obtained for a plasma torch with three IEs (L0= 490 mm) by varying the arc currents from 400 A to 600 A at a fixed gas-flow rate of Gin= 1050 slpm with Sw= 3.0. In addition, figures 7 and 8 show the temperature contours in the plasma torch and their radial profiles at z = 0.015 m,respectively,for arc currents of I0= 400 A and 600 A. In table 7, first, plasma power P0is observed to increase proportionally to arc current at the rate ofabout 0.655 kW A?1. Since plasma power P0can be expressed as P0=for plasma resistance Rp,this increase in P0proportionally to I0indicates that Rpis decreasing with the increase of I0. Normally, the electrical conductivity of N2plasma is improved as the plasma temperature is increased[14]. Accordingly, the decrease in Rpseems to be primarily caused by the increase in electrical conductivity because of the increase of plasma temperatures in the arc column, as shown in figures 7 and 8. With the increase of plasma temperature, however, thermal efficiency decreases because of radiation heat loss, as listed in table 7. Despite the reduction of thermal efficiency, however, specific enthalpy is increased because of the elevation of plasma power by the increasing arc current.

    Table 8.Calculation results for plasma torch with first IE(L0 = 250)obtained by varying the swirl number from 0.0 to 3.0 at a fixed gasflow rate, Gin = 750 slpm.

    3.5. Effects of swirling

    Table 8 presents the main numerical results obtained for three different swirl numbers of Sw= 0.0,1.5 and 3.0,respectively,in the plasma torch only with the first IE. The simulated plasma torch corresponds to the one of figure 2(a) with the lengths of IEs, L0= 250 mm and the operating conditions,such as arc current I0and gas flow rates are the same with the values used for figure 2(a). Normally, the elevation of swirl number increases the radius and the electrical conductivity of the arc [22]. These effects can be confirmed in table 8, presenting the arc voltages and the plasma power decreasing with the increase in swirl number at the same operating conditions.However, the exit enthalpy is calculated to be increased with the increase in Sw,which leads to the improvement of thermal efficiency at the same operating condition. For example, the data in table 8 indicate that the swirling injection with Sw= 1.5 increases thermal efficiency by 23.5% despite the decrease of input power by about 2% in the designed plasma torch, as listed in table 5.

    The effects of swirl on the radius and the electrical conductivity can also be found in figures 9(a)-(c), showing the temperature contours obtained for three different swirl numbers,Sw= 0.0, 1.5 and 3.0, respectively. From figure 9(a) for Sw= 0.0, for example, it is observed that a cold zone with temperatures lower than 1000 K is formed in the entrance region of the first IE by the radial injection of cold nitrogen gases without a swirl. With the decrease of temperatures in the arc column developed along the center line, however, the plasma gas injection with a swirl reduces this cold zone significantly by expanding the hot region of plasma radially, as shown in figures 9(b) and (c) representing the temperature fields for Sw= 1.5 and 3.0, respectively. As discussed previously, this radial expansion of the hot region brings the increase of electrical conductivity to the arc together with the decrease of maximum temperatures in the arc column,which contributes to lowering the arc voltage along the centerline.

    Figures 10(a)-(c) show the contours of azimuthal component,w of arc plasma with different number of IEs and swirl number. For example, figure 10(a) is obtained from the calculation results at Sw= 1.5 for the plasma torch with the first IE while figures 10(b) and (c) illustrate the results computed at Sw= 1.5 and 3.0, respectively, for the plasma torch with first, second and third IE.

    From the comparison between figures 10(a) and (b), the swirl velocity field spreads from the cathode gas inlet to the next electrodes along the surface of the first IE. However, it decays rapidly in the first IE,then disappears in the section of the second IE for the plasma torch with three IEs,as shown in figure 10(b). Since the swirl velocity field provides the stabilizing radial force against the kink instability in the arc[22],it is important to keep the swirl flow field along all sections of IEs when operating this type of plasma torch. Accordingly,the effect of swirl number needs to be investigated together with optimum gas flow rates injected through the gaps between electrodes. For example, figure 10(c) shows that the swirl velocity field can be kept along the surfaces of three IEs by increasing the swirl number up to 3.0.

    4. Conclusions

    In this work,we investigated the effects of IE insertion on the performance characteristics of a hollow-electrode plasma torch, such as arc voltage, plasma power, and thermal efficiency, by using numerical analysis. First, the numerical results demonstrated that arc voltage is elevated proportionally to the number of inserted IEs,resulting in the increase in plasma power at fixed arc current.In addition,when the gases are injected additionally through the gaps between IEs,it was found that thermal efficiency can be improved simultaneously with the increases in arc voltage and plasma power.

    The numerical results for temperature fields inside the torch revealed that these unique performance characteristics primarily come from the structural advantages of IE insertion,which allows for the injection of additional gases through gaps between the IEs, together with the increase in arc voltages. For example, the additionally injected gases can cool down the arc column and contract it radially, increasing arc resistance.As a result,the arc voltage can be increased a little more to balance between the generation and the removal of heat and charged particles in the arc column, leading to the increase in not only plasma power but also thermal efficiency.However, the increases in only arc currents without the addition of plasma gases or IEs were calculated to cause no increase in arc voltages, and accordingly, the thermal efficiency decreased, as expected in the conventional hollowelectrode plasma torches. In addition, numerical results revealed that the increase in swirl number can also improve thermal efficiency by increasing exit enthalpy despite decreasing arc voltages at the fixed arc current and gas flow rate.

    It should be noted that in long arcs with a length larger than the radius,such as the arc column formed in figure 1(a),deviations from LTE or LCE (locally chemical equilibrium)can be common rather than exceptional, meaning that the thermal flow fields of arcs can be different from those computed using the LTE model.As investigated by Wu et al[21],however,the non-LCE effects are negligible in the arc column with a temperature of ≥8000 K,where most of the arc voltage is calculated. Accordingly, the numerical results obtained from the LTE model can be used to estimate the effects of IE insertion on the torch performance, such as arc voltage and plasma power, despite the inability to describe the non-LCE effects.

    By optimizing these performance characteristics according to the IEs'insertion,we expect that a hollow-electrode plasma torch with IEs can be applied to many practical applications requiring high power and high efficiency operation,such as the thermal plasma treatment of hazardous wastes.

    Acknowledgments

    This research was supported by the Technology Development Program to Solve Climate Changes of the National Research Foundation(NRF)funded by the Ministry of Science and ICT(NRF-2016M1A2A2940152).

    ORCID iDs

    午夜福利高清视频| 婷婷精品国产亚洲av| 国产熟女xx| 国产亚洲精品一区二区www| 手机成人av网站| 成人永久免费在线观看视频| 黑人操中国人逼视频| 少妇熟女aⅴ在线视频| 色婷婷久久久亚洲欧美| 少妇熟女aⅴ在线视频| 一二三四在线观看免费中文在| 又黄又爽又免费观看的视频| 亚洲色图av天堂| 久久精品91蜜桃| 久久久久久人人人人人| 91成人精品电影| 在线免费观看的www视频| 一二三四社区在线视频社区8| 久久久水蜜桃国产精品网| 国产私拍福利视频在线观看| 国产精品免费一区二区三区在线| 88av欧美| 女生性感内裤真人,穿戴方法视频| 黑丝袜美女国产一区| 性色av乱码一区二区三区2| 久久亚洲真实| 久久久久九九精品影院| 国产真实乱freesex| 97超级碰碰碰精品色视频在线观看| 欧美色欧美亚洲另类二区| 亚洲国产精品久久男人天堂| 老熟妇乱子伦视频在线观看| 在线观看一区二区三区| 久久久久国内视频| 国内精品久久久久精免费| 国产亚洲精品久久久久久毛片| 国产高清激情床上av| 少妇裸体淫交视频免费看高清 | tocl精华| 亚洲人成网站高清观看| 国产高清视频在线播放一区| 琪琪午夜伦伦电影理论片6080| 欧美日韩福利视频一区二区| 丁香六月欧美| 成人国产一区最新在线观看| 欧美激情极品国产一区二区三区| 欧美日韩福利视频一区二区| 亚洲中文日韩欧美视频| 久久人人精品亚洲av| АⅤ资源中文在线天堂| 午夜老司机福利片| 亚洲自偷自拍图片 自拍| 亚洲美女黄片视频| 久久久久久久精品吃奶| 最新美女视频免费是黄的| 侵犯人妻中文字幕一二三四区| 少妇被粗大的猛进出69影院| 又大又爽又粗| 99精品在免费线老司机午夜| 韩国精品一区二区三区| 男女之事视频高清在线观看| 国产亚洲欧美精品永久| 婷婷精品国产亚洲av| 欧美黑人精品巨大| 亚洲五月天丁香| 日本三级黄在线观看| 欧美精品亚洲一区二区| 亚洲成av片中文字幕在线观看| 两性夫妻黄色片| 国产在线精品亚洲第一网站| 亚洲成人久久性| АⅤ资源中文在线天堂| 欧美丝袜亚洲另类 | 亚洲 欧美一区二区三区| 悠悠久久av| 欧美激情极品国产一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品综合一区在线观看 | 亚洲性夜色夜夜综合| 国产亚洲欧美98| 久久精品亚洲精品国产色婷小说| 中文字幕人成人乱码亚洲影| 午夜福利在线观看吧| 午夜久久久久精精品| 夜夜看夜夜爽夜夜摸| 99精品欧美一区二区三区四区| 午夜精品在线福利| 成人精品一区二区免费| 啦啦啦免费观看视频1| 国产成人精品无人区| 一卡2卡三卡四卡精品乱码亚洲| 婷婷亚洲欧美| 日本精品一区二区三区蜜桃| 亚洲成人久久爱视频| 99久久99久久久精品蜜桃| 久久人人精品亚洲av| 国产人伦9x9x在线观看| 黑人巨大精品欧美一区二区mp4| 一进一出好大好爽视频| 久久九九热精品免费| 淫妇啪啪啪对白视频| 日韩欧美 国产精品| 女生性感内裤真人,穿戴方法视频| 国产成人av教育| 欧美日本视频| 午夜福利欧美成人| 欧美乱妇无乱码| 久久人妻av系列| 欧美性长视频在线观看| 亚洲三区欧美一区| 久久久久久九九精品二区国产 | 最好的美女福利视频网| 久久亚洲精品不卡| 免费在线观看完整版高清| 国产精品自产拍在线观看55亚洲| 日日干狠狠操夜夜爽| 青草久久国产| 成年女人毛片免费观看观看9| 免费在线观看日本一区| 无限看片的www在线观看| 日本黄色视频三级网站网址| 99re在线观看精品视频| 久久久水蜜桃国产精品网| 亚洲av第一区精品v没综合| 国产成人影院久久av| 色综合欧美亚洲国产小说| 脱女人内裤的视频| 一本大道久久a久久精品| 亚洲国产高清在线一区二区三 | 欧美av亚洲av综合av国产av| 久久婷婷人人爽人人干人人爱| 国产av一区二区精品久久| 黄频高清免费视频| 午夜激情福利司机影院| 美女扒开内裤让男人捅视频| 国产又黄又爽又无遮挡在线| 国产aⅴ精品一区二区三区波| 色播在线永久视频| 手机成人av网站| 精品一区二区三区视频在线观看免费| 国产男靠女视频免费网站| 欧美黄色片欧美黄色片| 黑丝袜美女国产一区| 国产精品久久视频播放| 久久久久久亚洲精品国产蜜桃av| 久久久国产成人精品二区| 真人一进一出gif抽搐免费| 男人的好看免费观看在线视频 | 日本免费a在线| 制服诱惑二区| 好男人在线观看高清免费视频 | 久久精品人妻少妇| 97碰自拍视频| 美女高潮到喷水免费观看| 久99久视频精品免费| 黄片播放在线免费| 精品少妇一区二区三区视频日本电影| av天堂在线播放| 国产高清激情床上av| 国产黄a三级三级三级人| 精品国产乱码久久久久久男人| 国产色视频综合| 久久伊人香网站| 欧美日韩瑟瑟在线播放| 久久欧美精品欧美久久欧美| 欧美黑人精品巨大| 亚洲九九香蕉| 日韩大码丰满熟妇| 婷婷丁香在线五月| 夜夜躁狠狠躁天天躁| 老司机在亚洲福利影院| 午夜免费成人在线视频| 欧美黑人欧美精品刺激| 大型av网站在线播放| 精品无人区乱码1区二区| 一本精品99久久精品77| 亚洲va日本ⅴa欧美va伊人久久| 精品一区二区三区视频在线观看免费| 久久久久久九九精品二区国产 | 亚洲自拍偷在线| 黄色视频,在线免费观看| 精品一区二区三区四区五区乱码| 99精品在免费线老司机午夜| 亚洲精品中文字幕在线视频| 美女高潮喷水抽搐中文字幕| 在线观看午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 757午夜福利合集在线观看| 中文字幕人成人乱码亚洲影| 久久人妻福利社区极品人妻图片| 国产国语露脸激情在线看| 精品国产乱子伦一区二区三区| 在线观看www视频免费| 美女高潮喷水抽搐中文字幕| 午夜福利在线观看吧| 丰满人妻熟妇乱又伦精品不卡| 国产黄片美女视频| 午夜a级毛片| 精品国产乱码久久久久久男人| 国产国语露脸激情在线看| 观看免费一级毛片| 欧美 亚洲 国产 日韩一| 大型黄色视频在线免费观看| 在线永久观看黄色视频| 黄色视频,在线免费观看| 久久九九热精品免费| 精品午夜福利视频在线观看一区| 久久人妻av系列| 国产熟女xx| 久久久久久人人人人人| 国产亚洲精品久久久久5区| 黄色视频不卡| 亚洲中文av在线| 一级片免费观看大全| 18美女黄网站色大片免费观看| 男人的好看免费观看在线视频 | 禁无遮挡网站| 在线av久久热| 色尼玛亚洲综合影院| 久久国产亚洲av麻豆专区| 老熟妇乱子伦视频在线观看| 18禁观看日本| 日本 av在线| 两个人视频免费观看高清| 黄色片一级片一级黄色片| 久久国产精品人妻蜜桃| 色婷婷久久久亚洲欧美| 99国产极品粉嫩在线观看| 久久久久久久午夜电影| 日本一区二区免费在线视频| 丁香六月欧美| 日韩欧美国产一区二区入口| 色综合站精品国产| 国产欧美日韩一区二区精品| 一本精品99久久精品77| 欧美成人午夜精品| 在线永久观看黄色视频| 夜夜爽天天搞| 久久午夜综合久久蜜桃| 欧美成人午夜精品| 在线天堂中文资源库| 一进一出抽搐动态| 香蕉久久夜色| 黄色女人牲交| 午夜福利视频1000在线观看| 99精品在免费线老司机午夜| 成人特级黄色片久久久久久久| 91麻豆精品激情在线观看国产| 51午夜福利影视在线观看| 亚洲免费av在线视频| 国产精品亚洲美女久久久| 国产男靠女视频免费网站| 搡老妇女老女人老熟妇| 国产97色在线日韩免费| 香蕉久久夜色| 精华霜和精华液先用哪个| 日本免费a在线| 色av中文字幕| 曰老女人黄片| av欧美777| 中亚洲国语对白在线视频| 级片在线观看| 国产99白浆流出| 90打野战视频偷拍视频| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 极品教师在线免费播放| 夜夜躁狠狠躁天天躁| 亚洲国产精品sss在线观看| 久9热在线精品视频| 香蕉久久夜色| 精品久久久久久久毛片微露脸| 欧美日韩亚洲国产一区二区在线观看| 亚洲电影在线观看av| 国产成人精品无人区| 色综合婷婷激情| 男女床上黄色一级片免费看| 在线播放国产精品三级| 草草在线视频免费看| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 国产色视频综合| 国产精品一区二区精品视频观看| 狠狠狠狠99中文字幕| 国产成人精品久久二区二区免费| 久久香蕉国产精品| 成人一区二区视频在线观看| 天堂动漫精品| 91老司机精品| 国产精品免费一区二区三区在线| 在线播放国产精品三级| 久久亚洲精品不卡| 美女 人体艺术 gogo| 中文资源天堂在线| 日韩精品免费视频一区二区三区| av视频在线观看入口| 精品欧美一区二区三区在线| 国产成+人综合+亚洲专区| 村上凉子中文字幕在线| cao死你这个sao货| √禁漫天堂资源中文www| 每晚都被弄得嗷嗷叫到高潮| av有码第一页| 欧美成人免费av一区二区三区| 色老头精品视频在线观看| 1024香蕉在线观看| a在线观看视频网站| 国内精品久久久久精免费| 亚洲一区二区三区色噜噜| 国产黄a三级三级三级人| 可以在线观看的亚洲视频| 女人爽到高潮嗷嗷叫在线视频| 村上凉子中文字幕在线| 一夜夜www| 亚洲成a人片在线一区二区| 国产伦一二天堂av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久大精品| 国产激情久久老熟女| 少妇熟女aⅴ在线视频| 男女午夜视频在线观看| 国产欧美日韩精品亚洲av| 精品国产乱码久久久久久男人| aaaaa片日本免费| 一级毛片女人18水好多| 日本免费一区二区三区高清不卡| 在线永久观看黄色视频| 麻豆av在线久日| 成人三级做爰电影| 中文字幕久久专区| 免费在线观看影片大全网站| 欧美激情极品国产一区二区三区| 老司机午夜十八禁免费视频| 婷婷精品国产亚洲av| 亚洲中文av在线| 国产亚洲欧美精品永久| 精品一区二区三区视频在线观看免费| 免费观看人在逋| 丁香六月欧美| 男人舔女人下体高潮全视频| 男人操女人黄网站| 中文字幕精品免费在线观看视频| 国产一卡二卡三卡精品| 久久这里只有精品19| av免费在线观看网站| 成人亚洲精品一区在线观看| 午夜成年电影在线免费观看| 国产高清视频在线播放一区| 亚洲七黄色美女视频| 嫩草影院精品99| 日韩欧美国产在线观看| 亚洲自拍偷在线| 欧美黑人巨大hd| 啦啦啦韩国在线观看视频| 99re在线观看精品视频| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品av在线| 国产91精品成人一区二区三区| 亚洲欧美激情综合另类| 此物有八面人人有两片| av片东京热男人的天堂| 国产野战对白在线观看| 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 可以免费在线观看a视频的电影网站| 人人妻人人澡欧美一区二区| 国产区一区二久久| 丁香欧美五月| 99久久综合精品五月天人人| 日本 av在线| √禁漫天堂资源中文www| avwww免费| 亚洲精品国产一区二区精华液| av在线播放免费不卡| 校园春色视频在线观看| 黄色女人牲交| 日本一本二区三区精品| 波多野结衣高清作品| 国内精品久久久久精免费| www.999成人在线观看| 50天的宝宝边吃奶边哭怎么回事| 一二三四社区在线视频社区8| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品国产精品久久久不卡| 99精品在免费线老司机午夜| a级毛片a级免费在线| 亚洲国产欧美网| 天天添夜夜摸| 欧美国产精品va在线观看不卡| 国产伦人伦偷精品视频| 手机成人av网站| 一级片免费观看大全| 国产精品98久久久久久宅男小说| 黄片大片在线免费观看| 12—13女人毛片做爰片一| 在线观看免费午夜福利视频| 此物有八面人人有两片| 久久久国产成人精品二区| 一级a爱片免费观看的视频| 久久草成人影院| 精品国产一区二区三区四区第35| 国产日本99.免费观看| 久久精品91蜜桃| 欧美日韩乱码在线| 国产片内射在线| 精品电影一区二区在线| 欧美成人一区二区免费高清观看 | 88av欧美| 成熟少妇高潮喷水视频| 91字幕亚洲| 男人的好看免费观看在线视频 | 久久国产精品影院| 中文字幕另类日韩欧美亚洲嫩草| 99精品在免费线老司机午夜| 免费av毛片视频| 久久草成人影院| 大型黄色视频在线免费观看| 日韩大码丰满熟妇| 女性被躁到高潮视频| 国产精品av久久久久免费| 可以在线观看毛片的网站| 中文字幕人成人乱码亚洲影| 日本一区二区免费在线视频| 黄网站色视频无遮挡免费观看| 国内毛片毛片毛片毛片毛片| 婷婷亚洲欧美| 男女视频在线观看网站免费 | 久久香蕉精品热| 久久九九热精品免费| 成人av一区二区三区在线看| 欧美色欧美亚洲另类二区| www.精华液| 巨乳人妻的诱惑在线观看| √禁漫天堂资源中文www| 久久中文看片网| 免费在线观看日本一区| 俄罗斯特黄特色一大片| 色老头精品视频在线观看| 麻豆av在线久日| 亚洲中文字幕一区二区三区有码在线看 | 熟女电影av网| 日韩有码中文字幕| 亚洲国产精品999在线| 一卡2卡三卡四卡精品乱码亚洲| 在线永久观看黄色视频| 窝窝影院91人妻| 99国产综合亚洲精品| 天天一区二区日本电影三级| 日韩高清综合在线| 少妇粗大呻吟视频| 性欧美人与动物交配| 日本 av在线| 日本一本二区三区精品| 亚洲av电影在线进入| 欧美国产精品va在线观看不卡| 久久精品91无色码中文字幕| 国产精品综合久久久久久久免费| 性欧美人与动物交配| 久热这里只有精品99| 两个人看的免费小视频| or卡值多少钱| 久久久久国产精品人妻aⅴ院| 最近最新免费中文字幕在线| 12—13女人毛片做爰片一| xxxwww97欧美| 午夜福利一区二区在线看| 午夜免费鲁丝| 成人免费观看视频高清| 日韩精品中文字幕看吧| 给我免费播放毛片高清在线观看| 亚洲精品在线观看二区| 久久伊人香网站| 亚洲av成人av| 香蕉久久夜色| 丝袜人妻中文字幕| 亚洲全国av大片| 国产v大片淫在线免费观看| 国产精品久久久久久精品电影 | 免费看美女性在线毛片视频| 又黄又爽又免费观看的视频| cao死你这个sao货| 亚洲最大成人中文| 日本免费a在线| 免费在线观看视频国产中文字幕亚洲| 香蕉丝袜av| 日日爽夜夜爽网站| 国产精品日韩av在线免费观看| 欧美日本视频| 中国美女看黄片| 亚洲专区中文字幕在线| 欧美另类亚洲清纯唯美| 长腿黑丝高跟| 欧美国产精品va在线观看不卡| 桃红色精品国产亚洲av| 免费看a级黄色片| 成人18禁在线播放| 久久久国产精品麻豆| 婷婷亚洲欧美| 欧美日韩一级在线毛片| 国产真人三级小视频在线观看| 男女那种视频在线观看| 国内毛片毛片毛片毛片毛片| 日韩中文字幕欧美一区二区| 老司机午夜十八禁免费视频| 在线永久观看黄色视频| 色播在线永久视频| 18美女黄网站色大片免费观看| 午夜福利高清视频| 亚洲国产中文字幕在线视频| 日韩三级视频一区二区三区| 国产精品国产高清国产av| 男人操女人黄网站| 久久国产精品男人的天堂亚洲| 欧美久久黑人一区二区| 91av网站免费观看| 亚洲中文字幕日韩| 欧美黄色淫秽网站| 深夜精品福利| 亚洲专区字幕在线| 午夜精品在线福利| 国产又爽黄色视频| 精品国内亚洲2022精品成人| 51午夜福利影视在线观看| 精品一区二区三区视频在线观看免费| 亚洲一区高清亚洲精品| 亚洲精品色激情综合| 亚洲国产高清在线一区二区三 | 麻豆久久精品国产亚洲av| 免费看日本二区| 99热只有精品国产| 亚洲自偷自拍图片 自拍| 国产色视频综合| 亚洲国产精品sss在线观看| 精品久久久久久久久久免费视频| 国产黄色小视频在线观看| 精品一区二区三区视频在线观看免费| 国产单亲对白刺激| 国产精品一区二区三区四区久久 | 90打野战视频偷拍视频| 一级黄色大片毛片| av片东京热男人的天堂| 久久久久久久久久黄片| 日本熟妇午夜| 国产真实乱freesex| 日韩精品免费视频一区二区三区| 久久国产亚洲av麻豆专区| av超薄肉色丝袜交足视频| 日韩精品青青久久久久久| 国产av一区在线观看免费| 嫩草影院精品99| 亚洲熟女毛片儿| 性欧美人与动物交配| 天堂动漫精品| 又黄又爽又免费观看的视频| 亚洲成人久久爱视频| 久久婷婷人人爽人人干人人爱| 看黄色毛片网站| 日韩高清综合在线| 亚洲自拍偷在线| 亚洲欧美日韩高清在线视频| 久久久久久久精品吃奶| 欧美+亚洲+日韩+国产| 国产片内射在线| 亚洲九九香蕉| 黄色 视频免费看| 91麻豆av在线| 久久久久久久午夜电影| 国产视频一区二区在线看| 性欧美人与动物交配| 久久精品影院6| 国语自产精品视频在线第100页| 国产成人一区二区三区免费视频网站| 91麻豆精品激情在线观看国产| 久久精品亚洲精品国产色婷小说| 精品欧美国产一区二区三| 午夜a级毛片| 精品卡一卡二卡四卡免费| 午夜免费激情av| 美女免费视频网站| 真人做人爱边吃奶动态| 中国美女看黄片| 天天躁狠狠躁夜夜躁狠狠躁| 成人一区二区视频在线观看| 午夜福利在线观看吧| 亚洲美女黄片视频| 国产精品电影一区二区三区| 亚洲五月天丁香| 午夜影院日韩av| 久久青草综合色| 国产精品久久视频播放| 亚洲成人精品中文字幕电影| 99精品在免费线老司机午夜| 国产精品久久视频播放| 午夜福利18| 日日爽夜夜爽网站| 99riav亚洲国产免费| 精华霜和精华液先用哪个| 精品久久久久久久人妻蜜臀av| 高清毛片免费观看视频网站| 国产精品久久久av美女十八| 精品久久久久久久人妻蜜臀av| 日日干狠狠操夜夜爽| 精品国产亚洲在线| x7x7x7水蜜桃| 啦啦啦免费观看视频1| 日韩欧美国产在线观看| 在线天堂中文资源库| 色综合婷婷激情| 狠狠狠狠99中文字幕| 黑人巨大精品欧美一区二区mp4| 久久香蕉激情| 精品久久久久久,| 正在播放国产对白刺激| 久久久水蜜桃国产精品网| 无限看片的www在线观看|