Effect of desorbed gas on microwave breakdown on vacuum side of dielectric window

The gas desorbed from the dielectric surface has a great influence on the characteristics of microwave breakdown on the vacuum side of the dielectric window. In this paper, the dielectric surface breakdown is described by using the electromagnetic particle-in-cell-Monte Carlo collision(PIC-MCC) model. The process of desorption of gas and its influence on the breakdown characteristics are studied. The simulation results show that, due to the accumulation of desorbed gas, the pressure near the dielectric surface increases in time, and the breakdown mechanism transitions from secondary electron multipactor to collision ionization. More and more electrons generated by collision ionization drift to the dielectric surface, so that the amplitude of self-organized normal electric field increases in time and sometimes points to the dielectric surface. Nevertheless, the number of secondary electrons emitted in each microwave cycle is approximately equal to the number of primary electrons. In the early and middle stages of breakdown, the attenuation of the microwave electric field near the dielectric surface is very small. However, the collision ionization causes a sharp increase in the number density of electrons,and the microwave electric field decays rapidly in the later stage of breakdown. Compared with the electromagnetic PIC-MCC simulation results, the mean energy and number of electrons obtained by the electrostatic PIC-MCC model are overestimated in the later stage of breakdown because it does not take into account the attenuation of microwave electric field. The pressure of the desorbed gas predicted by the electromagnetic PIC-MCC model is close to the measured value,when the number of gas atoms desorbed by an incident electron is taken as 0.4.

Hollow cathode effect in radio frequency hollow electrode discharge in argon

Radio frequency capacitively coupled plasma source(RF-CCP)with a hollow electrode can increase the electron density through the hollow cathode effect(HCE),which offers a method to modify the spatial profiles of the plasma density.In this work,the variations of the HCE in one RF period are investigated by using a two-dimensional particle-in-cell/Monte-Carlo collision(PIC/MCC)model.The results show that the sheath electric field,the sheath potential drop,the sheath thickness,the radial plasma bulk width,the electron energy distribution function(EEDF),and the average electron energy in the cavity vary in one RF period.During the hollow electrode sheath's expansion phase,the secondary electron heating and sheath oscillation heating in the cavity are gradually enhanced,and the frequency of the electron pendular motion in the cavity gradually increases,hence the HCE is gradually enhanced.However,during the hollow electrode sheath's collapse phase,the secondary electron heating is gradually attenuated.In addition,when interacting with the gradually collapsed hollow electrode sheaths,high-energy plasma bulk electrons in the cavity will lose some energy.Furthermore,the frequency of the electron pendular motion in the cavity gradually decreases.Therefore,during the hollow electrode sheath's collapse phase,the HCE is gradually attenuated.

Direct modeling for computational fluid dynamics

All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equations. The current computational fluid dynamics （CFD） focuses on the numerical solution of partial differential equations （PDEs）, and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numer- ical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require fur- ther expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional dis- tinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of con- structing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm develop- ment. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be mod- eled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct construction of dis- crete numerical evolution equations, where the mesh size and time step will play dynamic roles in the modeling process. With the variation of the ratio between mesh size and local particle mean free path, the scheme will capture flow physics from the kinetic particle transport and collision to the hydro- dynamic wave propagation. Based on the direct modeling, a continuous dynamics of flow motion will be captured in the unified gas-kinetic scheme. This scheme can be faithfully used to study the unexplored non-equilibrium flow physics in the transition regime.

COLLISION AVOIDANCE DECISION-MAKING MODEL OF MULTI-AGENTS IN VIRTUAL DRIVING ENVIRONMENT WITH ANALYTIC HIERARCHY PROCESS

Collision avoidance decision-making models of multiple agents in virtual driving environment are studied. Based on the behavioral characteristics and hierarchical structure of the collision avoidance decision-making in real life driving, delphi approach and mathematical statistics method are introduced to construct pair-wise comparison judgment matrix of collision avoidance decision choices to each collision situation. Analytic hierarchy process （AHP） is adopted to establish the agents＇ collision avoidance decision-making model. To simulate drivers＇ characteristics, driver factors are added to categorize driving modes into impatient mode, normal mode, and the cautious mode. The results show that this model can simulate human＇s thinking process, and the agents in the virtual environment can deal with collision situations and make decisions to avoid collisions without intervention. The model can also reflect diversity and uncertainly of real life driving behaviors, and solves the multi-objective, multi-choice ranking priority problem in multi-vehicle collision scenarios. This collision avoidance model of multi-agents model is feasible and effective, and can provide richer and closer-to-life virtual scene for driving simulator, reflecting real-life traffic environment more truly, this model can also promote the practicality of driving simulator.

A theoretical model of collision between soft-spheres with Hertz elastic loading and nonlinear plastic unloading

This paper presents a theoretical model on the normal(head-on) collision between soft-spheres on the basis of elastic loading of the Hertz contact for compression process and a nonlinear plastic unloading for restitution one,in which the parameters all are determined in terms of the material and geometric ones of the spheres,and the behaviors of perfect elastic,inelastic,and perfect plastic collisions appeared in the classical mechanics are fully described once a value of coefficient of restitution is specified in the region of 0 ≤ ε ≤ 1.After an empirical formula of the coefficient of restitution dependent on the impact velocity is suggested to fit the existing experimental measurements by means of the least square method,the predictions of the dependency and the collision duration are in well quantitative agreement with their experimental measurements.It is found that the measurable quantities are dependent on both the impact velocity and the parameters of spheres.Following this model,finally,an approach to determine the spring coefficient in the linear viscoelastic model of the collision is also displayed.These results obtained here will be significantly beneficial for the applications where a collision model is requested in the simulations of relevant grain flows and impact dynamics etc..

A new collision avoidance model for pedestrian dynamics

The pedestrians can only avoid collisions passively under the action of forces during simulations using the social force model, which may lead to unnatural behaviors. This paper proposes an optimization-based model for the avoidance of collisions, where the social repulsive force is removed in favor of a search for the quickest path to destination in the pedestrian＇s vision field. In this way, the behaviors of pedestrians are governed by changing their desired walking direction and desired speed. By combining the critical factors of pedestrian movement, such as positions of the exit and obstacles and velocities of the neighbors, the choice of desired velocity has been rendered to a discrete optimization problem. Therefore,it is the self-driven force that leads pedestrians to a free path rather than the repulsive force, which means the pedestrians can actively avoid collisions. The new model is verified by comparing with the fundamental diagram and actual data. The simulation results of individual avoidance trajectories and crowd avoidance behaviors demonstrate the reasonability of the proposed model.

AN IMPROVED NONADIABATIC COLLISION MODEL FOR ION-PIR FORMATION PROCESS:A+BC→A^++BC^- Zheng Ting CAI and Yu Guang MU

Hickman's fast nonadiabatic collision model for the ion-pair formation reaction A+BC→A^++BC^- was improved,where the classical trajectory has been represented by solution of motion equation UR=-dV(R)/dR, here V(R)is Morse potential.Employing this model to the CS+O_2→CS^++O_2^-reaction,a satisfactory agreement with experimental data has been obtained.

Quantum entanglement dynamics based on composite quantum collision model

By means of composite quantum collision models,we study the entanglement dynamics of a bipartite system,i.e.,two qubits S1 and S2 interacting directly with an intermediate auxiliary qubit SA,while SA is in turn coupled to a thermal reservoir.We are concerned with how the intracollisions of the reservoir qubits influence the entanglement dynamics.We show that even if the system is initially in the separated state,their entanglement can be generated due to the interaction between the qubits.In the long-time limit,the steady-state entanglement can be generated depending on the initial state of S1 and S2 and the environment temperature.We also study the dynamics of tripartite entanglement of the three qubits S1,S2,and SA when they are initially prepared in the GHZ state and separated state,respectively.For the GHZ initial state,the tripartite entanglement can be maintained for a long time when the collision strength between the environment qubits is sufficiently large.

X-ray spectra of high temperature tungsten plasma calculated with collisional radiative model

Tungsten is regarded as an important candidate of plasma facing material in international thermonuclear experimental reactor （ITER）, so the determination and modeling of spectra of tungsten plasma, especially the spectra at high temperature were intensely focused on recently. In this work, using the atomic structure code of Cowan, a collisional radiative model （CRM） based on the spin-orbit-split-arrays is developed. Based on this model, the charge state distribution of tungsten ions is determined and the soft X-ray spectra from high charged ions of tungsten at different temperatures are calculated. The results show that both the average ionization charge and line positions are well agreed with others calculations and measurements with discrepancies of less than 0.63% and 1.26%, respectively. The spectra at higher temperatures are also reported and the relationship between ion abundance and temperature is predicted in this work.

The Calculation of Collision Risk on Air-Routes Based on Variable Nominal Separation

In this paper,a new method to calculate collision risk of air-routes,based on variable nominal separation,is proposed. The collision risk model of air-routes,based on the time variable and initial time interval variable,is given. Because the distance and the collision probability vary with time when the nominal relative speed between aircraft is not zero for a fixed initial time interval,the distance,the variable nominal separation,and the collision probability at any time can be expressed as functions of time and initial time interval. By the probabilistic theory,a model for calculating collision risk is acquired based on initial time interval distribution,flow rates,and the proportion of aircraft type. From the results of calculations,the collision risk can be characterized by the model when the nominal separation changes with time. As well the roles of parameters can be shown more readily.

A Theoretical Approach to Estimating Bird Risk of Collision with Wind Turbines Where Empirical Flight Activity Data Are Lacking

There are standard procedures for collecting data on numbers of birds at sites being proposed for wind farm development and evaluating collision risk for each key species. However, methods do not work well for all species. Where a local bird population is depleted, empirical data cannot provide estimates of likely collision mortality numbers if that population returns to satisfactory conservation status. Field survey methods are also inadequate for cryptic bird species. Both these problems can be important for evaluation of impacts of proposed wind farms on bird populations protected by the EU Birds Directive. We present an alternative method, based on energy constrained activity budgets and natural history, which permits assessment of likely collision numbers where empirical data are inadequate. Two case studies are presented where this approach has been successfully used to resolve disputed planning applications, one for a hen harrier population where numbers present are much below the population size at designation, and one for a cryptic species (greenshank). Our novel method helps reduce uncertainty in assessments constrained by difficulties in collecting suitable empirical data.

Modeling of Inner Surface Modification of a Cylindrical Tube by Plasma-Based Low-Energy Ion Implantation

The inner surface modification process by plasma-based low-energy ion implantation（PBLEII）with an electron cyclotron resonance（ECR）microwave plasma source located at the central axis of a cylindrical tube is modeled to optimize the low-energy ion implantation parameters for industrial applications.In this paper,a magnetized plasma diffusion fluid model has been established to describe the plasma nonuniformity caused by plasma diffusion under an axial magnetic field during the pulse-off time of low pulsed negative bias.Using this plasma density distribution as the initial condition,a sheath collisional fluid model is built up to describe the sheath evolution and ion implantation during the pulse-on time.The plasma nonuniformity at the end of the pulse-off time is more apparent along the radial direction compared with that in the axial direction due to the geometry of the linear plasma source in the center and the difference between perpendicular and parallel plasma diffusion coefficients with respect to the magnetic field.The normalized nitrogen plasma densities on the inner and outer surfaces of the tube are observed to be about 0.39 and 0.24,respectively,of which the value is 1 at the central plasma source.After a 5μs pulse-on time,in the area less than 2 cm from the end of the tube,the nitrogen ion implantation energy decreases from 1.5 keV to 1.3 keV and the ion implantation angle increases from several degrees to more than 40°;both variations reduce the nitrogen ion implantation depth.However,the nitrogen ion implantation dose peaks of about 2×10^（10）-7×10^（10）ions/cm^2 in this area are 2-4 times higher than that of 1.18×10^（10）ions/cm^2 and 1.63×10^（10）ions/cm^2 on the inner and outer surfaces of the tube.The sufficient ion implantation dose ensures an acceptable modification effect near the end of the tube under the low energy and large angle conditions for nitrogen ion implantation,because the modification effect is mainly determined by the ion implantation dose,just as the mass transfer process in PBLEII is dominated by low-energy ion implantation and thermal diffusion.Therefore,a comparatively uniform surface modification by the low-energy nitrogen ion implantation is achieved along the cylindrical tube on both the inner and outer surfaces.

THERMAL MODELLING OF COLLISIONAL OROGENY: IMPLICATIONS FOR THE ULTRA-HIGH PRESSURE METAMORPHISM

The petrological research on the ultra high pressure metamorphism (UHP) of collisional orogen indicates that the upper crustal rocks is subducted to depths exceeding 100 km, and returned to the surface rapidly. In this study, we investigate the thermal structure of collisional orogen as a slab of continental lithosphere being subducted beneath an overriding wedge of continental lithosphere by the 2 D finite element method. The advection heat transfer due to the accretion of orogenic wedge is considered. The wedge is composed of the upper crust materials through the accretion from the down going plate to the upper plate. For identifying the significance of the geometric and/or kinetic factors on the thermal structure of continental subduction, the different combinations of parameters, including dip angle of subduction zone, accretion or erosion rates, and the convergence velocity etc., are used in modelling. The time span of continental subduction in our calculation is less than 30 Ma, according to the short duration of ultra deep subduction of continental slab suggested by the preservation of metastable pre peak low pressure mineralogy assemblage in the garnet of UHP rocks. Therefore, the steep dip angle of down going plate and/or low rate of accretion favour the ultra deep subduction of upper crust materials, especially for the slower down going slab. Meanwhile, taking the erosion rate as the level of exhumation rate of UHP rocks in some orogens (i.e., 1-2 km/Ma or more) does not result in the anatexis melting of crust of the overriding plate, due to the cooling effect of the rapid down going slab. However, the temperature structures of all models are generally cooler than those recovered by thermobarometric studies of the UHP rocks. This implies the significant increase of temperature after the rapid subduction of continental slab. Following the method of Davies and von Blackenburg (1998), we show that the slab breakoff can occur at the depth exceeding 100 km. Thermal modelling on the post subduction stage shows the heating related to the plate breakoff can cause the higher temperature recorded by the exhumed UHP rocks. The higher geotherm during post subduction stage leads to the weak strength of the orogenic wedge, and favours the faster upward movement of the UHP rock slices as ductile agents. The lower temperature gradient of the subduction slab predicted by modelling suggests the cold subducting slab could have transported significant fluids to mantle depth, not released during subduction. Accordingly, the absence of coeval calc alkalic magmatism in UHP orogens might resulted from the lower temperature as well as the fluid free circumstance, both are related to the rapid subduction of cold plate. Therefore, shear heating is not needed for explanation the thermal evolution of UHP orogen. On the other hand, the post collisional or late stage granitic plutonism is closely related to the deep seated heat producing materials of the accretion wedge.

Elliptic Flow Splitting between Particles and their Antiparticles in Au＋Au Collisions from a Multiphase Transport Model

The elliptic flow v2, for π±, K±, p and p in Au＋Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the v2 values for p and p is observed, and the values of v2 splitting are larger compared with π＋ and π-, K＋ and K-. The difference increases with decreasing the center-of-mass energy. The effect of the quark coalescence mechanism in a multi-phase transport model to the value of elliptic difference △v2 between p and p- has been discussed. The simulation of Au＋Au collisions at 14.5 GeV shows that the effect of hadron cascade to △v2 is not obvious, and a larger patton-scattering cross section can lead to a larger △v2.

Non-invasive optical characterization and estimation of Zn porosity in gas tungsten arc welding of Fe–Al joints using CR model and OES measurements

In this study,we employed a non-invasive approach based on the collisional radiative(CR)model and optical emission spectroscopy(OES)measurements for the characterization of gas tungsten arc welding(GTAW)discharge and quantification of Zn-induced porosity during the GTAW process of Fe–Al joints.The OES measurements were recorded as a function of weld current,welding speed,and input waveform.The OES measurements revealed significant line emissions from Zn-I in 460–640 nm and Ar-I in 680–800 nm wavelength ranges in all experimental settings.The OES coupled CR model approach for Zn-I line emission enabled the simultaneous determination of both essential discharge parameters i.e.electron temperature and electron density.Further,these predictions were used to estimate the Zn-induced porosity using OES-actinometry on Zn-I emission lines using Ar as actinometer gas.The OES-actinometry results were in good agreement with porosity data derived from an independent approach,i.e.x-ray radiography images.The current study shows that OES-based techniques can provide an efficient route for real-time monitoring of weld quality and estimate porosity during the GTAW process of dissimilar metal joints.

Development of the PARASOL Code and Full Particle Simulation of Tokamak Plasma with an Open-Field SOL-Divertor Region Using PARASOL

The PARASOL code and the simulation by using PARASOL are introduced briefly. The PARASOL code with particle-in-cell （PIC） method and binary collision model was developed in JAERI and JAEA. Simulations using PARASOL code were carried out in order to investigate the power and particle control with diveror system in fusion reactors. The one-dimensional （1D） version of PARASOL was adopted to investigate the Bohm criterion, the supersonic flow, the SOL heat conduction, and so on. The heat propagation due to edge localized mode （ELM） was studied with the 1D-dynamic PARASOL. The two-dimensional version of PARASOL for the whole tokamak plasma including scrape-off-layer （SOL）-divertor region was useful for simulating the SOL flow pattern, the electric field formation etc. Based on PARASOL simulation results, improved physics modeling for the fluid simulation was built up.

Soft X-Ray Laser Gain from Neon like Gallium and Germanium

Gain coefficients are calculated for neon-like gallium and germanium ions. Shorter wavelengths are calculated and predicted to be emitted. The gain coefficients are calculated among 457 energy levels of the neon-like ions. Collisional excitations were calculated through the distorted wave approximations through five electron temperatures Te = 300, 500, 700, 1000 and 1500 eV.

Nonadiabatic collision model calculations of ion-pair formation cross sections of Cs+O2→Cs++O2-

This paper has improved Hickman's nonadiabatic collision model by substituting Hickman's constant velocity classical straight line trajectory approximation with the solution of motion equation mR=-dV(R)/dR, and has calculated the cross sections of ion-pair formation Cs+O2 -Cs++O2- with the improved nonadiabatic collision model (INCM). A comparison of our results with other theoretical and experimental results has been made.

Effect of entanglement embedded in environment on quantum non-Markovianity based on collision model

By means of collision models(CMs)where the environment is simulated by a collection of ancillas consisting of two entangled qubits,we investigate the effects of entanglement in the environment on the non-Markovianity of an open quantum system.Two CMs are considered in this study,in the first one the open quantum system S directly collides with the environment,while in the second one the system interacts with two intermediate qubits which,in turn,are coupled to the environment.We show that it is possible to enhance the non-Markovianity by environment entanglement in both models.In particular,in the second model,we show that the initial state of the auxiliary qubits can also affect the non-Markovianity of the system and there exists the optimal combination of the initial environmental state and the initial state of auxiliary qubits.In this case,the non-Markovianity can be greatly enhanced.

Collision mitigation and vehicle transportation safety using integrated vehicle dynamics control systems

The aim of this paper is to investigate the effect of vehicle dynamics control systems （VDCS） on both the collision of the vehicle body and the kinematic behaviour of the ve- hicle＇s occupant in case of offset frontal vehicle-to-vehicle collision. A unique 6-degree-of- freedom （6-DOF） vehicle dynamics/crash mathematical model and a simplified lumped mass occupant model are developed. The first model is used to define the vehicle body crash parameters and it integrates a vehicle dynamics model with a vehicle front-end structure model. The second model aims to predict the effect of VDCS on the kinematics of the occupant. It is shown from the numerical simulations that the vehicle dynamics/crash response and occupant behaviour can be captured and analysed quickly and accurately. Yurthermore, it is shown that the VDCS can affect the crash characteristics positively and the occupant behaviour is improved.