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

    Convection effect on an arc plasma evolution process in a two parallel contact system

    2020-05-06 05:59:14JianningYIN尹健寧ShunguiLIU劉順貴FengTANG唐峰QishenLV呂啟深andXingwenLI李興文
    Plasma Science and Technology 2020年4期
    關(guān)鍵詞:興文

    Jianning YIN (尹健寧),Shungui LIU (劉順貴),Feng TANG (唐峰),Qishen LV (呂啟深) and Xingwen LI (李興文)

    1 State Key Laboratory of Electrical Insulation and Power Equipment,Xi’an Jiaotong University,Xi’an 710049,People’s Republic of China

    2 Shenzhen Power Supply Co.Ltd,Shenzhen 518020,People’s Republic of China

    Abstract

    Keywords: air arc,circuit breaker,MHD

    1.Introduction

    A low voltage circuit breaker(LVCB)is a crucial component for controlling and protecting a power distribution network,and its movable contact consists of two parallel contacts for products with a higher rated current [1].The arc will be ignited between movable contacts and a static contact when the fault current appears,and the arc evolution process will directly affect the breaking capacity of the LVCB.Therefore,the study of the arc evolution process of a two parallel contact structure is crucial.

    Many studies have been performed on air arc behaviors.For a magneto–hydro–dynamic (MHD) simulation of an air arc,the air plasma physical properties must be obtained first.Murphy et al[2,3]calculated the transport coefficients of air,nitrogen–air and oxygen–air plasma.They also calculated the air and different metal vapor mixture plasma properties.Ander et al [4,5] and Zhang et al [6] calculated the composition and thermodynamic properties of ablated vapors of Cu,PMMA,PA6-6,PETP,POM and PE.In addition,the modeling of the plasma–sheath boundary region and the influence of the surface structure on arc voltage were studied[7,8].A sheath with a voltage–current density characteristic was proposed to represent the near-anode and near-cathode voltage,and also,the formation process of a new arc root on a splitter plate is described by it [9].Based on the above parameters,the effects of several factors,such as different mixture plasma properties [10],the ferromagnetic material[11,12],the arc ignition position [13],the erosion of the splitter plate [14] and the moving contact motion [15–17] on the arc behavior,were studied by the MHD model.Also,the behavior of the arc root was studied in[18–20].However,the quasi-steady solution for a magnetic field was assumed in the calculation.Thus,the temporal variation of the magnetic field was neglected.The studies [21,22] calculated the eddy current in the splitter plate and electrodes.

    However,all of the above studies considered single contact models and few papers have been published studying the air arc evolution process in a parallel contacts model.

    In this paper,a 3D MHD simulation model of two parallel contacts is built to investigate the arc evolution process.To represent the voltage drop of the electrodes,a 0.1 mm sheath with a voltage–current density characteristic is applied.For the numerical method,in order to take the temporal variation of the magnetic field in the whole calculation domain into account,the ANSYS Emag solver is applied to calculate the time-varying electromagnetic field.The coupling between the flow field and electromagnetic field uses the coupling server MpCCI.Based on the above methods,the evolution process of the anode arc root is analyzed by the current density distribution.The convection and conduction flux densities of air plasma are calculated and their effects on arc motion are quantitatively analyzed.The shape of the arc column is analyzed.The asymmetry of the anode and cathode arc root is also studied.The arc motion process of the two parallel contact system is preliminarily revealed.The study in this paper is beneficial for understanding the arc evolution mechanism of a two parallel contact system,and provides theoretical support to optimize high performance low voltage circuit breakers.

    2.Numerical model and governing equations

    The arc evolution process is not only determined by the electromagnetic field,but also by the gas dynamic.The arc model presented here is based on the MHD approach,where the arc plasma is described as one fluid.The plasma characteristic justifies some assumptions for the simplification.

    2.1.Hypotheses

    (1) The conditions of the local thermodynamic equilibrium and the laminar flow are satisfied.

    (2) Evaporation from the electrodes and wall materials is not considered.The physical properties of air plasma are the functions of temperature and pressure [2].

    (3) No ferromagnetic materials in the domain are presented.

    2.2.Geometry model

    Figure 1(a) shows the simulation model for the fluid.This model is a simplification of the miniature circuit breaker,which includes the electrode domain and arc column region.The dimensions of the model are 8 mm × 8 mm × 30 mm in the x–y–z direction.The origin of the coordinate is on the left wall (z = 0).In order to solve the electro-magnetic field,a far-away domain is built in the model.The thickness of the electrodes is 1.5 mm,and 0.1 mm near electrode layers surround the electrodes.Figure 1(b) shows the anode geometry of two parallel contacts.The anode consists of two parallel contacts,and 12 solid contacts connect them.The dimensions of contact 1 (contact 2) are 3 mm × 1.5 mm × 23 mm,and the dimensions of the 12 solids are 8 mm × 1.5 mm × 2 mm.The dimensions of the cathode are 8 mm × 1.5 mm ×25 mm.The material of the electrodes is copper.The distance between the two electrodes is 5 mm.

    Figure 1.The geometry model.

    Table 1.Equations solved in the fluid model.

    2.3.Governing equations

    Basically,the plasma process is modeled in detail combining the finite-element and finite-volume methods based on an MHD approach [15].A thin layer of elements with a nonlinear voltage–current density characteristic is employed to represent the voltage drop of the electrodes [13].The MHD equations are a combination of Navier–Stokes flow equations extended to include arc specific source terms with additional Maxwell equations for the electromagnetic field [13].

    2.3.1.Fluid model.The Navier–Stokes flow equations are written in the form of the generalized conservation equation (1):

    where Φ is the variable solved,ρ is the mass density,t is the time,is the velocity vector,ΓΦis the diffusion coefficient andSΦis the source term.All the equations solved for the fluid model are summarized in table 1 [23,24].

    In table 1,

    whereσE2is the Joule heat,which represents the injected energy into the air plasma.The radiation is calculated by equation (6):

    whereεnis the net emission coefficient [25].

    In the above equations,vx,vy,vzare the velocities of the x,y,z direction,jx,jy,jzare the current densities of the x,y,z direction,Bx,By,Bzare the magnetic flux densities of the x,y,z direction,is the Lorentz force produced by the current loop,η is the viscosity,λ is the thermal conductivity,cpis the specific heat,qηis the viscous dissipation,p is pressure,and h is enthalpy.

    The Navier–Stokes flow equations are solved using the ANSYS Fluent solver [15].

    2.3.2.Electromagnetic model.The motion of an arc is not only affected by the gas dynamic,but also by the magnetic field (induced by the arc itself and the current loop of the anode–arc–cathode).Therefore,the Lorentz forceis added in the source item of the momentum balance equation as shown in equations (2)–(4).

    The magnetic vector potential method is adopted to solve the nonlinear transient electro-magnetic field.The equations are shown as follows [26]

    where μ0is the permeability of vacuum,and σ is the electrical conductivity.The electrical conductivity of plasma is the function of the plasma temperature and pressure,φ is the electrical potential,andis the magnetic vector potential.

    The electromagnetic equations (7) and (8) are solved using the ANSYS Emag solver.

    During the simulation,the data transfer between the two solvers is necessary due to the interaction between the gas flow field and the electro-magnetic field.This is done by using the MpCCI [15].The electrical conductivity σ is transferred from the ANSYS Fluent to the ANSYS Emag.The Lorentz force and Joule heat are transferred from the ANSYS Emag to ANSYS Fluent by user defined codes.

    The time-varying electromagnetic field induced by the current loop,including the electrodes and arc column,is calculated in this paper.Therefore,compared with most of the previous studies [5–13],the itemis included in equations (7) and (8).

    2.3.3.Near-electrode model.The near-electrode voltage is an important part of the arc voltage and must be considered in the calculation.The 0.1 mm layers with a nonlinear resistance surrounding the electrodes are applied and its effective electrical conductivity is calculated by equation (9) [27]:

    where J is the current density,Δyis the thickness of the near electrode layer,and Ulis the voltage drop in the near electrode layer.The value of the effective electrical conductivity is calculated by the voltage–current density characteristic [13].

    A higher voltage UIhas to be exceeded at a low current density to describe the formation process of the arc root.After the formation of the new arc root,the near-electrode voltage U0is a constant.In this paper,the higher voltage UIis set to 22.3 V,and the near-electrode voltage U0is set to 10.3 V [13].

    The near-electrode voltage (near-anode and near-cathode) is calculated by a nonlinear voltage–current density characteristic [13].The difference in voltage between the near-anode and near-cathode is not considered in this model.Therefore,the arc voltage may be reduced.In addition,if the difference in voltage between the near-anode and nearcathode is considered in the model,the temperature of the anode arc root will be higher.

    2.4.Boundary conditions

    (1) Velocity: for all walls,the velocity is set to no-slip conditions.

    (2) Temperature: for the side wall (including the front surface,back surface and left surface of the arc column region,as shown in figure 1(a)) and the electrode outside planes,the environment temperature of 300 K is applied to the thermal conditions.For the pressure outlet,as shown in figure 1(a),the temperature is set to 300 K[27].At the interface between the air plasma and electrode surfaces,the temperature is determined as [13]:

    (1) Pressure:for the pressure outlet,as shown in figure 1(a),the total pressure is set to 101325 Pa.

    (2) Electrical and magnetic field: the electrical and magnetic field is solved in ANSYS Emag.In order to ensure the current is freely distributed between the two contacts according to the position of the arc root,the arc current is applied on the surface ‘inlet’.The electrical potential is set to zero on the surface ‘outlet’,as shown in figure 1(a) [14].For the interface between the arc region and the electrode,the boundary condition of the electrical field is described by equation (11):

    For the side wall of the arc region and electrode outside planes,the boundary condition is described by equation(12):

    For the boundary of the magnetic field,the zero magnetic potential is set to a far-away domain boundary to calculate the magnetic field.At the interface between the arc region and electrode,equations (13) and (14) are used [26]:

    The calculating current is set tosin (ωt+π/12) A,whereω=2πf.The equations are solved with a time step of 2.5 μs,and the coupled solver of Fluent is used.To improve the initial convergence,the simulation begins with a steady state result.The position of the initial arc at z = 5 mm in this paper and this arc ignition point are determined by the distance between the arcing position and the outlet of the model according to the structure of the miniature circuit breaker.

    The difference between the model established in this paper and the previous model is mainly that this model considers the two parallel contacts and their interactions.Also,the eddy current in the whole domain (arc domain and electrodes) is considered.In addition,the convection flux is calculated and its effect on the arc plasma evolution process in a two parallel contact system is analyzed in this paper.

    Figure 2.Current density distribution sequences in the arc column domain.

    3.Numerical results and analysis

    Figure 3.Curves of current through two contacts and arc voltage.The whole process is divided into processes(1)–(6).S1–S4 represent stage (1)–stage (4).

    The arc evolution process in a two parallel contact system is simulated based on a 3D MHD approach.The current density distribution,temperature distribution and the arc voltage are obtained.The arc shape from the anode to the cathode is analyzed by the current density distribution.The displacements of the arc roots are obtained.The convection flux density and heat conduction flux density are calculated at various instants.The arc evolution process and the influence of convection on it in a two parallel contacts model are analyzed.

    3.1.Evolution process of an arc

    The current density distribution sequences in the arc column domain are shown in figure 2.The red lines represent the position of the two parallel contacts.The current flows from the anode to the cathode.

    The curves of the current through the two contacts and arc voltage are shown in figure 3.In figure 3,i is the arc current,i1and i2are the currents through arc branch 1 and arc branch 2,respectively,u is the arc voltage.The curves of the arc resistance (including total arc resistance,resistance of arc branch 1 and arc branch 2) are shown in figure 4.The resistance trend of the two arc branches is similar to the trend of the current flowing through the two arc branches.The arc motion process can be analyzed in the following stages according to the current curves,as shown in figure 3.

    3.1.1.The anode arc root only appears on contact 1(t=0–0.38 ms).The starting position of the arc is at z = 5 mm.After that,the arc column moves along the electrodes due to the gas flow force and Lorentz force,so the arc voltage increases,as shown in figure 3.Also,the arc column expands.The heated air flows from arc branch 1 to arc branch 2 due to the convection effect,which results in the air under contact 2 being heated; therefore,the electrical conductivity of the air plasma under contact 2 increases.At this stage,the arc root of the anode only appears on contact 1,as shown in figure 2(a),but the convective effect provides favorable conditions for the formation of a new arc root on contact 2.After t = 0.3 ms,the current through the arc branch 2 increases.

    Figure 4.Curves of arc resistance. Rarc is the total arc resistance, R1 and R2 represent the resistance of arc branch 1 and arc branch 2,respectively.

    3.1.2.A new arc root is formed on contact 2(t=0.38–0.495 ms).Stage 1 (t = 0.38–0.445 ms): the electrical conductivity of the air under contact 2 is constantly increased.At t = 0.38 ms,a new arc root is formed on contact 2,as shown in figure 2(b),and the arc voltage reaches its peak.

    After t = 0.38 ms,the arc column under the anode consists of two branches.The arc column under contact 1 is defined as arc branch 1 and the arc column under contact 2 is defined as arc branch 2 for the convenience of the following analysis.At t = 0.38 ms,the current density and arc column radius of arc branch 1 are larger than those of arc branch 2.The temperature gradient between the two arc branches is the largest.The air plasma of arc branch 2 is heated and the electrical conductivity is increased due to convection and conduction.Therefore,current i2increases.In contrast,current i1decreases.The resistance of the two arc branches is equal until t = 0.445 ms.

    Stage 2 (t = 0.445–0.495 ms): the current flowing through arc branch 2 continues to increase due to thermal inertia after t = 0.445 ms.The arc resistance of branch 2 is less than that of arc branch 1.Therefore,the current flowing through branch 2 will continue to increase due to the decrease of arc resistance R2.When t = 0.495 ms,the current flowing through arc branch 2 reaches its local peak.

    3.1.3.Current i1increases(t=0.495–0.655 ms).Stage 3(t = 0.495–0.565 ms): at t = 0.495 ms,as shown in figure 2(c),the current density and radius of arc branch 2 are larger than those of arc branch 1.After that,the electrical conductivity of arc branch 1 increases,which leads to an increase in the current through arc branch 1.The arc currents of the two branches are equal again at t = 0.565 ms.

    Figure 5.Schematic diagram of the dynamic u–i curves.

    Stage 4 (t = 0.565–0.655 ms): the current flowing through arc branch 1 continues to increase at this stage due to the decrease of arc resistance R1,which leads to an increase of the temperature gradient between the two arc branches.

    The arc evolution processes (4) and (6) are similar to process (2),and the evolution process (5) is similar to process (3).

    An alternative variation of R1and R2appears in the whole process.However,the total arc resistance always reduces in the arc evolution process according to figure 4;that is to say,the arc in a two parallel contact system always develops toward the direction in which the total arc resistance decreases.

    At process (2) and process (3) of the arc evolution process,the arc voltage decreases because the cross section of the arc column becomes larger.The arc voltage increases slightly from t = 0.655 ms to t = 0.8 ms,while it decreases again at t = 0.8 ms.After t = 1.0 ms,the arc voltage increases because of the elongation of the arc column.

    Figure 5 shows the dynamic voltage–current curves of the two arcs,where the curve of the black solid line represents the dynamic u–i1characteristic of arc branch 1 and the red dash line represents the u–i2characteristic of arc branch 2 [28].

    In order to analyze all possible intersections of the two curves,curves (1)–(3) show the possible situations of arc branch 1.The points a–f are possible points of intersection.

    Through the analysis of each intersection in figure 5,the two parallel arc branches hardly co-exist at points from‘a(chǎn)’–‘c’ (one of two arc branches will be extinguished).At points from‘d’–‘g’,the two parallel arc branches can co-exist.

    In summary,if conditions (1) du/di1> 0 and du/di2> 0 or (2) du/di1< 0 and du/di2< 0 are satisfied,that is to say,when the arc resistance characteristics of two arc branches are consistent,two parallel arc branches can coexist.

    It can be seen that the above-mentioned analysis based on the dynamic current–voltage curves could not explain the calculation results as shown in figure 3.It is predicted that the interaction between the two parallel arc branches should be considered from the view of a flow field.

    3.2.Calculation results of convection and conduction

    Figure 6.Temperatures on the cross section of y = 2.4 mm along the x-axis at four instants.

    In order to further analyze the arc evolution mechanism in a two parallel contact system,the temperatures at the cross section of the anode arc root(y = 2.4 mm)along the x-axis at four instants (figures 2(c)–(f)) are obtained and shown in figure 6.

    It can be noted that the position of arc root 1 and arc root 2 is at about x = 2.2 mm and x = 5.8 mm respectively.At t = 0.655 ms and t = 1.1 ms,the central temperature of arc root 1 is significantly higher than that of arc root 2,which leads to the direction of the temperature gradient along the positive direction of the x-axis.At t = 0.495 ms and t = 0.915 ms,the central temperature of arc root 2 is higher than that of arc root 1,which leads to the direction of the temperature gradient along the –x-axis direction.The convection and conduction are directly related to the temperature gradient.The convection and conduction will affect the movement of air plasma and the arc evolution process.

    In this paper,the convective flux density is defined as[29]:

    The conduction flux density is calculated by the following formula:

    The distribution of the x component of the conduction flux density on cross section y = 2.4 mm at four instants(figures 2(c)–(f)) is shown in figure 8.

    The arc motion process is mainly affected by convection and conduction.It can be seen that the convection flux density is much larger than the conduction flux density by comparing figures 7 and 8.That is to say,convection plays a dominant role in the arc evolution process.

    Figure 7.The x component of convective flux density on cross section y = 2.4 mm ((c) = 0.495 ms,(d) = 0.655 ms,(e) = 0.915 ms,(f) = 1.1 ms).

    Figure 8.The x component of conduction flux density on cross section y = 2.4 mm ((c) = 0.495 ms,(d) = 0.655 ms,(e) = 0.915 ms,(f) = 1.1 ms).

    In figure 8,the curves of (c)–(f) are evaluated for the cross section at z = 9 mm,z = 11.5 mm,z = 14.5 mm and z = 17 mm,respectively (the position of the black dash line in figure 9).Figure 8 shows that the conduction flux density is zero at x = 2.2 mm and x = 5.8 mm,respectively,and the heated air plasma diffuses to both sides of the arc root.The heated air plasma diffuses to the middle of the two arc roots due to conduction,therefore,the positions of the two anode arc roots are maintained in a fixed position.

    3.3.The convection effect on the arc evolution process

    The flow field distribution (on the cross section of the anode arc root)at the four instants is obtained and shown in figure 9.Based on this,the effect of convection and conduction on arc evolution behavior is analyzed.

    The boxes of the black dash line represent the positions of contact 1 and contact 2,as shown in figure 9(a).

    Figure 9.Velocity distribution at different instants.

    The vertical black dash line shows the position of the anode arc root with a larger current density.It is noted that the convective flux density in the position of the arc root equals zero at about at t = 0.495 ms,t = 0.655 ms and t = 1.1 ms,as shown in figure 7,because the velocity of the airflow is along the z-axis,while the x component of the velocity is very small,as shown in figures 9(a),(b),(d).At t = 0.915 ms,the convection flux density equals 2 × 108W m-2due to the large x component of velocity in the position of the arc root,as shown in figure 9(c).

    3.3.1.The evolution behavior at stage(2).The new arc root forms on contact 2 at t = 0.38 ms.The current density of arc root 1 is higher with the arc motion,resulting in the direction of convection being from arc branch 1 to branch 2.The electrical conductivity of arc branch 2 is increased and the arc resistance of arc branch 2 begins to decrease due to convection.Therefore,current i2increases.At point ‘a(chǎn)’,as shown in fgiure 4,i1= i2.The i2current continues to increase and i1decreases after point‘a(chǎn)’ due to the thermal inertia of the arc column.Until t = 0.495 ms (curve (c) as shown in fgiure 7),the convection fulx density is less than zero when z > 11.5 mm and the convection fulx density is greater than zero when z < 9 mm,which shows that the cool air behind the arc column folws to arc branch 2 and the heated air before the arc column flows to arc branch 1,as shown in fgiure 9(a).Therefore,the electrical conductivity of arc branch 1 is increased,so that current i1increases after t = 0.495 ms,as shown in figure 3.

    3.3.2.The evolution process at stage(3).From 0.495 ms to 0.565 ms,the arc resistance R1decreases and R2increases,as shown in fgiure 4,due to the effect of convection.Therefore,current i1increases and i2decreases.At t = 0.565 ms,current i1= i2at point ‘b’,as shown in fgiure 3.From 0.565 ms to 0.655 ms,the arc resistance R1continues to decrease due to thermal inertia,resulting in current i1continuing to increase and i2continuing to decrease.Therefore,the temperature gradient between the two arc roots increases with the increase of the difference value of current density between the two arc branches,which leads to an enhancement of the convection effect and its direction is from arc root 1 to arc root 2.At t = 0.655 ms,the difference value of the current density between the two arc branches is the highest.The heated air before the arc column folws towards branch 2,as shown in fgiure 9(b).The electrical conductivity of branch 2 increases due to convection(curve(d)as shown in fgiure 7).Therefore,current i2increases and i1decreases after t = 0.655 ms.This explains the current change process at stage (3) during the arc evolution process.

    The convection effects on the arc evolution process at stages (4) and (6) are similar to stage (2),and the convection effect on the arc evolution process at stage (5) is similar to stage (3).The arc evolution process is repetitive.

    As per the above mentioned analysis,the current densities of the two arc branches under the anode do not reach the steady state in the whole arc motion process because the arc motion process is a dynamic process.Therefore,the current change process of i1and i2is a dynamic equilibrium process.

    It can be seen that two arc roots appear in the anode with an increase in the current through the above analysis.The change in the current flowing through the two branches is mainly affected by the convection and arc resistance in the whole arc evolution process in the two parallel contact system according to the above analyses.However,the electrical conductivity of the arc branches will be changed due to convection.Therefore,the arc resistance of the two arc branches will be changed.In other words,the main physical mechanism of the arc evolution process is the influence of convection.

    Figure 10.The displacement of the anode and cathode arc roots.

    3.4.The shape of the arc column and displacements of arc roots

    3.4.1.Shape of arc column.It can be noted that the size of the two arc columns varies alternately from figure 2.The new arc root forms on contact 2 after 0.38 ms.The anode has two arc roots and the cathode has one after 0.38 ms.From the anode to the cathode,the two arc branches merge into one.Therefore,the shape of the arc column is ‘Y’ type.One reason is the electromagnetic attraction of the two arc branches.The other reason is the conduction effect as shown in figure 8,which leads to the heated air plasma diffusing to the middle of the two arc roots.These two reasons also explain the phenomena of the two anode arc roots maintained in a fixed position.

    3.4.2.The displacements of arc roots.The structure of the anode and cathode is different.The anode consists of two parallel contacts and the cathode consists of one contact.It will affect the displacements of the arc roots.The displacements of the arc roots will be analyzed in this section.The displacements of the anode and cathode arc roots are shown in figure 10.

    The arc is ignited at z = 5 mm.Before t = 0.38 ms,the velocity of arc root 1 is faster than the cathode arc root because of the higher current density.Therefore,the displacement of anode arc root 1 is larger.At 0.38 ms,anode arc root 2 forms at z = 7.5 mm and it is behind anode arc root 1 because of the tailing of the arc.The equivalent cross section of the anode arc root increases after the formation of anode arc root 2,so that the current density of the cathode arc root is larger.Therefore,there is a slower velocity of the anode arc root than the cathode arc root.Thus,the cathode arc root will catch up with the anode arc root and the displacements of the anode and cathode arc roots are equal until 0.62 ms.The current density of the anode arc root increases more than the cathode arc root with the increase of current after 0.62 ms.Therefore,the displacement of the anode arc root is greater than the cathode arc root at 1.1 ms again,as shown in figure 10.This phenomenon can also be seen in figure 2(f).

    4.Conclusions

    A 3D simulation model of a two parallel contact system has been built.The anode consists of two parallel contacts in the model.A 0.1 mm near-electrode layer with nonlinear resistance is applied to simulate the formation of the arc root.The convection flux density and conduction flux density are calculated.Based on this,the arc evolution behavior of a two parallel contact system is analyzed.It can be concluded that:

    (1) The arc root will extend from one contact to two contacts on parallel contacts with the arc motion.The size of two arc branches will change alternately during the arc motion process.

    (2) Two anode arc roots can co-exist during the arc evolution behavior.The arc column of the two parallel contact system always develops toward the direction in which the total arc resistance decreases.The arc current change process through each arc branch is a dynamic process,and they try to reach the steady point in the arc motion process.The arc evolution behavior is affected by convection.

    (3) The arc column is‘Y’type after the formation of a new arc root due to the attraction of two arc branches near the anode and heat conduction.The arc root motion is asymmetrical because of the structural asymmetry.

    In future work,more attention should be paid to the plasma–sheath boundary region modelling,and the influence of the electrode motion process and splitter plates.

    Acknowledgments

    This work is supported by the Natural Science Foundation of Shaanxi Province of China (No.2017ZDJC-16) and Shenzhen Power Supply Co.Ltd (No.SZKJXM20170480).

    猜你喜歡
    興文
    我校詹興文副校長(zhǎng)申報(bào)項(xiàng)目獲第三批海南省高校思想政治工作精品項(xiàng)目立項(xiàng)
    一個(gè)放羊娃成就的“書香之家”
    火紅的日子
    宜賓試驗(yàn)站打造蠶桑產(chǎn)業(yè)強(qiáng)鎮(zhèn)——興文周家鎮(zhèn)
    興文蠶桑產(chǎn)業(yè)創(chuàng)新驅(qū)動(dòng)發(fā)展的做法與成效
    這 里 是 北 京
    青年歌聲(2019年6期)2019-06-26 09:14:42
    大澤山:從賣葡萄到興文旅
    商周刊(2018年19期)2018-10-26 03:31:20
    興文塔的建筑特色與藝術(shù)價(jià)值
    興文農(nóng)商銀行中干競(jìng)聘
    四川興文世界地質(zhì)公園景觀特征與價(jià)值評(píng)價(jià)
    国产色爽女视频免费观看| 日韩国内少妇激情av| 少妇猛男粗大的猛烈进出视频 | 国产午夜精品久久久久久一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 嫩草影院精品99| 亚洲第一区二区三区不卡| 国产精品一二三区在线看| 亚洲国产精品专区欧美| 久久久精品免费免费高清| 97热精品久久久久久| 亚洲自拍偷在线| 久久精品熟女亚洲av麻豆精品| av在线播放精品| 自拍偷自拍亚洲精品老妇| 国产成人a区在线观看| 少妇裸体淫交视频免费看高清| 麻豆久久精品国产亚洲av| 蜜桃亚洲精品一区二区三区| 一本色道久久久久久精品综合| 99re6热这里在线精品视频| 下体分泌物呈黄色| 日本wwww免费看| 不卡视频在线观看欧美| 熟女人妻精品中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人妻一区二区av| 男人爽女人下面视频在线观看| 亚洲熟女精品中文字幕| 国产伦精品一区二区三区四那| 亚洲综合精品二区| 九九久久精品国产亚洲av麻豆| 亚洲在久久综合| 一级黄片播放器| 最新中文字幕久久久久| 一个人看视频在线观看www免费| 51国产日韩欧美| 777米奇影视久久| 精品少妇黑人巨大在线播放| 婷婷色综合www| 亚洲成人一二三区av| 国产男女超爽视频在线观看| 久久精品国产a三级三级三级| 成年女人在线观看亚洲视频 | 亚洲国产av新网站| 一级毛片久久久久久久久女| 欧美极品一区二区三区四区| 欧美最新免费一区二区三区| 国产高清不卡午夜福利| 最近手机中文字幕大全| 一区二区三区精品91| 精品人妻一区二区三区麻豆| 三级男女做爰猛烈吃奶摸视频| 亚洲av免费高清在线观看| 只有这里有精品99| av线在线观看网站| 精品久久久久久电影网| 最近最新中文字幕大全电影3| 国产成人精品婷婷| 成人欧美大片| 免费观看性生交大片5| 国内精品美女久久久久久| 欧美人与善性xxx| 精品久久久精品久久久| 亚洲av中文av极速乱| 极品少妇高潮喷水抽搐| 国产亚洲最大av| 国产男女超爽视频在线观看| 免费在线观看成人毛片| 免费av毛片视频| 色5月婷婷丁香| 国产精品久久久久久av不卡| 中文天堂在线官网| 国产91av在线免费观看| 欧美另类一区| 国产熟女欧美一区二区| 超碰97精品在线观看| 国产精品精品国产色婷婷| 免费大片黄手机在线观看| 人妻少妇偷人精品九色| 日韩伦理黄色片| 特大巨黑吊av在线直播| 国产成人福利小说| 插逼视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲av.av天堂| 99视频精品全部免费 在线| 午夜精品国产一区二区电影 | 久久精品国产自在天天线| 天堂俺去俺来也www色官网| 身体一侧抽搐| 欧美日韩亚洲高清精品| 黑人高潮一二区| 亚洲无线观看免费| 日韩 亚洲 欧美在线| 久久久久久伊人网av| 国产探花极品一区二区| 亚洲人与动物交配视频| 97超视频在线观看视频| 一区二区三区乱码不卡18| 日本-黄色视频高清免费观看| 夜夜爽夜夜爽视频| 99热全是精品| 天天躁日日操中文字幕| 99热网站在线观看| 99视频精品全部免费 在线| 男女无遮挡免费网站观看| 精品人妻一区二区三区麻豆| 亚洲国产欧美在线一区| 一个人看视频在线观看www免费| 免费观看性生交大片5| 亚洲成人中文字幕在线播放| 人妻一区二区av| a级毛色黄片| 国产精品精品国产色婷婷| 久久久亚洲精品成人影院| 尤物成人国产欧美一区二区三区| 国产高清有码在线观看视频| 国产真实伦视频高清在线观看| 日韩在线高清观看一区二区三区| 丝袜脚勾引网站| 一级a做视频免费观看| 在线天堂最新版资源| 新久久久久国产一级毛片| 欧美成人一区二区免费高清观看| 免费观看av网站的网址| 欧美 日韩 精品 国产| 亚洲内射少妇av| 18禁裸乳无遮挡动漫免费视频 | 欧美最新免费一区二区三区| 精品99又大又爽又粗少妇毛片| 嫩草影院精品99| 中文欧美无线码| 亚洲欧美成人精品一区二区| 亚洲人与动物交配视频| 我要看日韩黄色一级片| 午夜日本视频在线| 欧美高清成人免费视频www| 少妇人妻精品综合一区二区| 男女边吃奶边做爰视频| 亚洲精品自拍成人| 91精品一卡2卡3卡4卡| 亚洲精品第二区| 亚洲av电影在线观看一区二区三区 | 麻豆成人午夜福利视频| 九色成人免费人妻av| 亚洲精品视频女| 国产精品不卡视频一区二区| 美女cb高潮喷水在线观看| 亚洲av男天堂| 国内精品美女久久久久久| 美女被艹到高潮喷水动态| 日韩一区二区三区影片| 狂野欧美激情性xxxx在线观看| 一区二区三区精品91| 一个人看视频在线观看www免费| 在现免费观看毛片| 欧美精品一区二区大全| 免费看日本二区| 国产黄a三级三级三级人| 成年女人看的毛片在线观看| 男女无遮挡免费网站观看| 黄色视频在线播放观看不卡| 一级毛片电影观看| 午夜日本视频在线| 搞女人的毛片| 午夜精品国产一区二区电影 | 国产综合精华液| 26uuu在线亚洲综合色| av福利片在线观看| 亚洲天堂国产精品一区在线| 久久久久九九精品影院| 久久人人爽人人爽人人片va| 成人毛片60女人毛片免费| 99久久精品一区二区三区| 成人鲁丝片一二三区免费| 最近中文字幕高清免费大全6| 美女主播在线视频| 老司机影院毛片| 午夜精品一区二区三区免费看| 搡女人真爽免费视频火全软件| 熟女电影av网| 观看免费一级毛片| 久久精品夜色国产| 嫩草影院新地址| 国内揄拍国产精品人妻在线| 欧美日韩综合久久久久久| 成人综合一区亚洲| 久久精品夜色国产| 日本av手机在线免费观看| 黄片wwwwww| 欧美三级亚洲精品| 欧美潮喷喷水| 99热这里只有是精品在线观看| 一级二级三级毛片免费看| 婷婷色综合大香蕉| 日本-黄色视频高清免费观看| 国产在线男女| 如何舔出高潮| 国产免费一区二区三区四区乱码| 精品午夜福利在线看| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜添av毛片| 欧美3d第一页| 久久6这里有精品| xxx大片免费视频| 一区二区三区精品91| 欧美最新免费一区二区三区| 免费人成在线观看视频色| 国产精品国产三级国产专区5o| 啦啦啦中文免费视频观看日本| 亚州av有码| 婷婷色麻豆天堂久久| 久久鲁丝午夜福利片| 成人亚洲精品av一区二区| 大香蕉97超碰在线| 一级毛片久久久久久久久女| 亚洲av.av天堂| 日韩欧美 国产精品| 国产欧美亚洲国产| 精品一区二区三卡| 嘟嘟电影网在线观看| 卡戴珊不雅视频在线播放| 亚洲内射少妇av| 人妻夜夜爽99麻豆av| 99re6热这里在线精品视频| 老女人水多毛片| 麻豆精品久久久久久蜜桃| 97热精品久久久久久| 国产亚洲午夜精品一区二区久久 | 欧美性猛交╳xxx乱大交人| 色哟哟·www| 91aial.com中文字幕在线观看| 久久99热这里只频精品6学生| 久久99精品国语久久久| 黄色日韩在线| kizo精华| 午夜精品一区二区三区免费看| 香蕉精品网在线| 大香蕉97超碰在线| 69av精品久久久久久| 欧美另类一区| 视频区图区小说| 久久99精品国语久久久| 精品国产露脸久久av麻豆| 欧美老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 大又大粗又爽又黄少妇毛片口| 香蕉精品网在线| 国产精品熟女久久久久浪| 国产高清有码在线观看视频| 日本猛色少妇xxxxx猛交久久| 婷婷色综合大香蕉| 国产老妇女一区| 欧美老熟妇乱子伦牲交| 联通29元200g的流量卡| 欧美成人午夜免费资源| 人人妻人人爽人人添夜夜欢视频 | 激情五月婷婷亚洲| 久久精品人妻少妇| 午夜福利在线观看免费完整高清在| 欧美激情国产日韩精品一区| 亚洲丝袜综合中文字幕| 亚洲美女搞黄在线观看| 国产人妻一区二区三区在| 日日摸夜夜添夜夜爱| 午夜激情久久久久久久| 91午夜精品亚洲一区二区三区| 国产一区二区三区综合在线观看 | 婷婷色av中文字幕| 一区二区三区乱码不卡18| 国产精品久久久久久精品电影| 少妇的逼水好多| 99热全是精品| 日韩亚洲欧美综合| 激情五月婷婷亚洲| 大码成人一级视频| 婷婷色综合大香蕉| 日韩av在线免费看完整版不卡| 高清午夜精品一区二区三区| 久久久久久久久久久免费av| 免费av观看视频| 日韩一区二区视频免费看| 极品少妇高潮喷水抽搐| 久久久久久久国产电影| 欧美xxxx黑人xx丫x性爽| 久久综合国产亚洲精品| 日本猛色少妇xxxxx猛交久久| 少妇人妻一区二区三区视频| 欧美xxⅹ黑人| 伦精品一区二区三区| 可以在线观看毛片的网站| 在线观看人妻少妇| 51国产日韩欧美| 国产淫片久久久久久久久| freevideosex欧美| 免费高清在线观看视频在线观看| 午夜福利在线在线| 欧美高清成人免费视频www| 国产 一区 欧美 日韩| 亚洲av一区综合| 精品人妻偷拍中文字幕| 亚洲国产日韩一区二区| 国产毛片在线视频| 日韩电影二区| 在线观看人妻少妇| 18禁动态无遮挡网站| 嫩草影院新地址| 九九久久精品国产亚洲av麻豆| 亚洲精品久久久久久婷婷小说| 成人亚洲欧美一区二区av| 国产亚洲av嫩草精品影院| 毛片一级片免费看久久久久| 在线观看三级黄色| 三级国产精品欧美在线观看| 亚洲欧美成人精品一区二区| 日本三级黄在线观看| 欧美 日韩 精品 国产| 欧美 日韩 精品 国产| 国产精品久久久久久精品电影小说 | 高清毛片免费看| 亚洲美女搞黄在线观看| 人妻系列 视频| 搡老乐熟女国产| 日韩大片免费观看网站| freevideosex欧美| 国产大屁股一区二区在线视频| 精品久久久久久电影网| av免费观看日本| av国产免费在线观看| 寂寞人妻少妇视频99o| 夫妻性生交免费视频一级片| 青春草视频在线免费观看| 高清在线视频一区二区三区| 女人被狂操c到高潮| 欧美成人一区二区免费高清观看| 国产精品人妻久久久影院| 色视频www国产| 久久久色成人| 亚洲av一区综合| 舔av片在线| 久久久国产一区二区| 欧美日韩视频高清一区二区三区二| 久久久久久伊人网av| 国产精品久久久久久精品电影| 国产亚洲最大av| 久久人人爽人人爽人人片va| 亚洲色图综合在线观看| 青青草视频在线视频观看| 亚洲最大成人av| 一区二区三区四区激情视频| 久久久色成人| 日日摸夜夜添夜夜添av毛片| 国产永久视频网站| 久久韩国三级中文字幕| 成人美女网站在线观看视频| 免费少妇av软件| 色5月婷婷丁香| 欧美高清成人免费视频www| 欧美一级a爱片免费观看看| 午夜精品一区二区三区免费看| 亚洲精品乱码久久久v下载方式| 午夜精品国产一区二区电影 | 国产爽快片一区二区三区| 成人亚洲精品av一区二区| 日日啪夜夜爽| 联通29元200g的流量卡| 小蜜桃在线观看免费完整版高清| 欧美一级a爱片免费观看看| 男女无遮挡免费网站观看| av专区在线播放| 国产老妇伦熟女老妇高清| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩综合久久久久久| 欧美日韩综合久久久久久| 日本免费在线观看一区| 亚洲精品aⅴ在线观看| 亚洲美女搞黄在线观看| 美女内射精品一级片tv| 丰满人妻一区二区三区视频av| 成人毛片a级毛片在线播放| 天天躁夜夜躁狠狠久久av| 午夜福利在线在线| 两个人的视频大全免费| 只有这里有精品99| 久久99蜜桃精品久久| 亚洲熟女精品中文字幕| 草草在线视频免费看| 国产一区亚洲一区在线观看| 国产视频首页在线观看| 激情 狠狠 欧美| 超碰av人人做人人爽久久| 久久97久久精品| 亚洲精品久久午夜乱码| 欧美日韩国产mv在线观看视频 | 亚洲av不卡在线观看| 干丝袜人妻中文字幕| 国产精品久久久久久精品古装| 晚上一个人看的免费电影| 国产永久视频网站| 亚洲第一区二区三区不卡| 熟妇人妻不卡中文字幕| 涩涩av久久男人的天堂| 欧美成人午夜免费资源| 免费黄频网站在线观看国产| 国产色爽女视频免费观看| 久久这里有精品视频免费| 久久久久国产网址| 亚洲国产日韩一区二区| 日韩国内少妇激情av| 五月开心婷婷网| 插逼视频在线观看| 国产欧美亚洲国产| 国产欧美亚洲国产| 国产女主播在线喷水免费视频网站| 九九在线视频观看精品| 亚洲精品乱码久久久久久按摩| 深夜a级毛片| 亚洲欧美精品专区久久| 亚洲av成人精品一二三区| 亚洲精品色激情综合| 狂野欧美激情性xxxx在线观看| 超碰av人人做人人爽久久| 男女下面进入的视频免费午夜| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜老司机福利剧场| 亚洲av电影在线观看一区二区三区 | 亚洲久久久久久中文字幕| 久热久热在线精品观看| 一级毛片久久久久久久久女| 久久午夜福利片| 精品人妻熟女av久视频| 亚洲av不卡在线观看| 极品少妇高潮喷水抽搐| 日韩欧美一区视频在线观看 | 国产女主播在线喷水免费视频网站| 男人舔奶头视频| 亚洲自拍偷在线| 精品午夜福利在线看| 少妇裸体淫交视频免费看高清| 亚洲精品国产成人久久av| 久久99热这里只频精品6学生| 欧美国产精品一级二级三级 | 简卡轻食公司| 最近手机中文字幕大全| av线在线观看网站| kizo精华| 在线a可以看的网站| 黄色配什么色好看| 亚洲最大成人手机在线| 色播亚洲综合网| 亚洲人与动物交配视频| 欧美成人午夜免费资源| 涩涩av久久男人的天堂| 在线观看美女被高潮喷水网站| 亚洲综合精品二区| 最近手机中文字幕大全| 国产淫片久久久久久久久| 18禁裸乳无遮挡动漫免费视频 | 国产黄片美女视频| 久久99精品国语久久久| 亚洲国产欧美在线一区| 日本av手机在线免费观看| 久久久久精品性色| 国产伦精品一区二区三区视频9| 亚洲国产欧美在线一区| 午夜福利视频精品| 久久久久网色| 久久久久久久大尺度免费视频| 九九爱精品视频在线观看| 久久久精品免费免费高清| 久久女婷五月综合色啪小说 | 日韩在线高清观看一区二区三区| 成人美女网站在线观看视频| 精品久久久噜噜| h日本视频在线播放| 深夜a级毛片| 国产成人免费观看mmmm| 国产亚洲最大av| 五月天丁香电影| 五月天丁香电影| 高清在线视频一区二区三区| 一个人看视频在线观看www免费| 亚洲欧美成人精品一区二区| 一级毛片黄色毛片免费观看视频| 国内少妇人妻偷人精品xxx网站| 亚洲av中文av极速乱| 亚洲真实伦在线观看| 日韩欧美精品v在线| 免费观看a级毛片全部| 国产亚洲最大av| 男人和女人高潮做爰伦理| 国产视频内射| 国产淫语在线视频| 成人亚洲欧美一区二区av| 国产一区二区三区av在线| 精品久久国产蜜桃| 嫩草影院新地址| 一本一本综合久久| av线在线观看网站| 国产欧美亚洲国产| 欧美日韩视频精品一区| 欧美精品国产亚洲| 国产精品久久久久久av不卡| 在线免费观看不下载黄p国产| 亚洲av成人精品一区久久| 人人妻人人看人人澡| 99视频精品全部免费 在线| 亚洲欧洲国产日韩| 黄色怎么调成土黄色| 99热6这里只有精品| 久久久久久久久久久丰满| 内射极品少妇av片p| 我的老师免费观看完整版| 男女下面进入的视频免费午夜| 狂野欧美激情性xxxx在线观看| 欧美少妇被猛烈插入视频| 在线观看三级黄色| 2021少妇久久久久久久久久久| 日本黄色片子视频| 联通29元200g的流量卡| 亚洲欧美精品自产自拍| 久久久久久伊人网av| 亚洲欧洲日产国产| 青青草视频在线视频观看| 嘟嘟电影网在线观看| 人人妻人人澡人人爽人人夜夜| 日日啪夜夜撸| 男的添女的下面高潮视频| 精品久久久久久久久亚洲| 精华霜和精华液先用哪个| 麻豆国产97在线/欧美| 亚洲欧美清纯卡通| 亚洲色图综合在线观看| 18禁在线无遮挡免费观看视频| 国产乱来视频区| 日韩视频在线欧美| 人人妻人人澡人人爽人人夜夜| 又大又黄又爽视频免费| 精品少妇久久久久久888优播| 热re99久久精品国产66热6| 久久人人爽人人片av| 国产一区有黄有色的免费视频| 国产精品国产av在线观看| 亚洲精品亚洲一区二区| 国产精品99久久久久久久久| 日韩,欧美,国产一区二区三区| 免费av不卡在线播放| 亚洲电影在线观看av| 国产亚洲5aaaaa淫片| 国产爱豆传媒在线观看| 国产精品无大码| 少妇人妻久久综合中文| 爱豆传媒免费全集在线观看| 狂野欧美激情性xxxx在线观看| av卡一久久| 亚洲精品成人av观看孕妇| 国产精品久久久久久久久免| 国产探花在线观看一区二区| 日本熟妇午夜| 夜夜看夜夜爽夜夜摸| 亚洲自偷自拍三级| 99久久精品国产国产毛片| 精品一区二区三卡| 交换朋友夫妻互换小说| www.色视频.com| 能在线免费看毛片的网站| 精品人妻熟女av久视频| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 男人和女人高潮做爰伦理| 午夜精品一区二区三区免费看| 两个人的视频大全免费| 嫩草影院精品99| 成人二区视频| 亚洲无线观看免费| 看十八女毛片水多多多| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 国产色婷婷99| 国内精品美女久久久久久| 亚洲在久久综合| 亚洲丝袜综合中文字幕| 国产精品一区二区在线观看99| 中文天堂在线官网| 美女xxoo啪啪120秒动态图| 丝袜脚勾引网站| 最后的刺客免费高清国语| 国产成人aa在线观看| a级一级毛片免费在线观看| 国产伦理片在线播放av一区| 亚洲av成人精品一二三区| 老司机影院成人| 日韩欧美一区视频在线观看 | 91精品国产九色| 波多野结衣巨乳人妻| 亚洲精品久久久久久婷婷小说| av一本久久久久| 国产色婷婷99| 久久久久精品久久久久真实原创| 久久精品国产亚洲网站| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 美女高潮的动态| 男女啪啪激烈高潮av片| 亚洲综合精品二区| 国产在线男女| 国产男人的电影天堂91| 观看免费一级毛片| 欧美三级亚洲精品| 欧美日韩国产mv在线观看视频 | 亚洲成人中文字幕在线播放| 午夜激情福利司机影院|