3D fluid model analysis on the generation of negative hydrogen ions for negative ion source of NBI

A radio-frequency(RF) inductively coupled negative hydrogen ion source(NHIS) has been adopted in the China Fusion Engineering Test Reactor(CFETR) to generate negative hydrogen ions.By incorporating the level-lumping method into a three-dimensional fluid model,the volume production and transportation of H^(-) in the NHIS,which consists of a cylindrical driver region and a rectangular expansion chamber,are investigated self-consistently at a large input power(40 k W) and different pressures(0.3–2.0 Pa).The results indicate that with the increase of pressure,the H^(-) density at the bottom of the expansion region first increases and then decreases.In addition,the effect of the magnetic filter is examined.It is noteworthy that a significant increase in the H^(-) density is observed when the magnetic filter is introduced.As the permanent magnets move towards the driver region,the H^(-) density decreases monotonically and the asymmetry is enhanced.This study contributes to the understanding of H-distribution under various conditions and facilitates the optimization of volume production of negative hydrogen ions in the NHIS.

A THREE-FLUID MODEL OF THE SAND-DRIVEN FLOW

The sand-driven flow is studied from the continuum viewpoint in this paper. The crux of this work is how to model the stresses of the particle phase properly. By analysing the two-fluid model which usually, works in solving gas-particle two-phase .flow,. we find that this model has many. deficiencies for studying the sand-driven flow,even for the simplest case- the steady, two-dimensional fully-developed flow.Considering this, we have proposed the three-fluid model in which the upward particles and the downward-particles ore regarded as two kinds of fluids respectively.It is shown that the three-fluid model is better than the two-fluid model in reflecting the internal structure of the flow, region and the influence of the boundary situations on the flow. and it is advantageous to find an approximate solution in that the main components of the particle-phase stresses can be explicitly expressed by those variables in the three-fluid model.In the end, the governing equations as well as the boundary. conditions for the three-fluid model are provided with a discussion.

Simulation of Variable Viscosity and Jeffrey Fluid Model for Blood Flow Through a Tapered Artery with a Stenosis

Non-Newtonian fluid model for blood flow through a tapered artery with a stenosis and variable viscosity by modeling blood as Jeffrey fluid has been studied in this paper. The Jeffrey fluid has two parameters, the relaxation time A1 and retardation time A2. The governing equations are simplified using the case of mild stenosis. Perturbation method is used to solve the resulting equations. The effects of non-Newtonian nature of blood on velocity profile, temperature profile, wall shear stress, shearing stress at the stenotsis throat and impedance of the artery are discussed. The results for Newtonian fluid are obtained as special case from this model.

Validity of the two-term Boltzmann approximation employed in the fluid model for high-power microwave breakdown in gas

The electron energy distribution function （EEDF）, predicted by the Boltzmann equation solver BOLSIG＋ based on the two-term approximation, is introduced into the fluid model for simulating the high-power microwave （HPM） breakdown in argon, nitrogen, and air, and its validity is examined by comparing with the results of particle-in-cell Monte Carlo collision （PIC/MCC） simulations as well as the experimental data. Numerical results show that, the breakdown time of the fluid model with the Maxwellian EEDF matches that of the PIC/MCC simulations in nitrogen; however, in argon under high pressures, the results from the Maxwellian EEDF were poor. This is due to an overestimation of the energy tail of the Maxwellian EEDF in argon breakdown. The prediction of the fluid model with the BOLSIG＋ EEDF, however, agrees very well with the PIC/MCC prediction in nitrogen and argon over a wide range of pressures. The accuracy of the fluid model with the BOLSIG＋ EEDF is also verified by the experimental results of the air breakdown.

The Efficient Finite Element Methods for Time-Fractional Oldroyd-B Fluid Model Involving Two Caputo Derivatives

In this paper,we consider the numerical schemes for a timefractionalOldroyd-B fluidmodel involving the Caputo derivative.We propose two efficient finite element methods by applying the convolution quadrature in time generated by the backward Euler and the second-order backward difference methods.Error estimates in terms of data regularity are established for both the semidiscrete and fully discrete schemes.Numerical examples for two-dimensional problems further confirmthe robustness of the schemes with first-and second-order accurate in time.

Stability analysis of wireless network with improved fluid model

Wireless communication is easily disturbed by unfortunate factors which drive the wireless environment unstable and complicated. Therefore, it is essential to consider these factors in stability analysis of the wireless network. However, wireless channel characteristics and packets collisions are neglected in the classical fluid model. A wireless TCP fluid model （WTFM） for stability analysis of wireless network is proposed based on cross layers, which not only makes the congestion control based on random early detection （RED） available for wireless network, but also provides a more accurate model to analyze the stability of wireless system theoretically. In the proposed model, active queue management, abnormality of wireless channels and packets collisions are taken into consideration. The comparisons between evaluating results from the WTFM and the practical performance from NS2 simulations validate the accuracy of the proposed WTFM in the perspectives of delay, dropping probability, throughput, sliding window size and queue length. A set of comparisons among the proposed WTFM, the classical fluid model and the convex optimization model are conducted. The results demonstrate that the proposed WTFM model performs better than other schemes in comprehensive aspects on capturing the characteristic of the wireless network and computing complexity.

LARGE-TIME BEHAVIOR OF SOLUTIONS TO THE INFLOW PROBLEM OF THE NON-ISENTROPIC MICROPOLAR FLUID MODEL

We investigate the asymptotic behavior of solutions to the initial boundary value problem for the micropolar fluid model in a half line R+:=(0,∞).Inspired by the relationship between a micropolar fluid model and Navier-Stokes equations,we prove that the composite wave consisting of the transonic boundary layer solution,the 1-rarefaction wave,the viscous 2-contact wave and the 3-rarefaction wave for the inflow problem on the micropolar fluid model is time-asymptotic ally stable under some smallness conditions.Meanwhile,we obtain the global existence of solutions based on the basic energy method.

Numerical investigation of radio-frequency negative hydrogen ion sources by a three-dimensional fluid model

A three-dimensional fluid model is developed to investigate the radio-frequency inductively coupled H2 plasma in a reactor with a rectangular expansion chamber and a cylindrical driver chamber,for neutral beam injection system in CFETR.In this model,the electron effective collision frequency and the ion mobility at high E-fields are employed,for accurate simulation of discharges at low pressures(0.3 Pa-2 Pa)and high powers(40 kW-100 kW).The results indicate that when the high E-field ion mobility is taken into account,the electron density is about four times higher than the value in the low E-field case.In addition,the influences of the magnetic field,pressure and power on the electron density and electron temperature are demonstrated.It is found that the electron density and electron temperature in the xz-plane along permanent magnet side become much more asymmetric when magnetic field enhances.However,the plasma parameters in the yz-plane without permanent magnet side are symmetric no matter the magnetic field is applied or not.Besides,the maximum of the electron density first increases and then decreases with magnetic field,while the electron temperature at the bottom of the expansion region first decreases and then almost keeps constant.As the pressure increases from 0.3 Pa to 2 Pa,the electron density becomes higher,with the maximum moving upwards to the driver region,and the symmetry of the electron temperature in the xz-plane becomes much better.As power increases,the electron density rises,whereas the spatial distribution is similar.It can be summarized that the magnetic field and gas pressure have great influence on the symmetry of the plasma parameters,while the power only has little effect.

Numerical study on discharge characteristics influenced by secondary electron emission in capacitive RF argon glow discharges by fluid modeling

A one-dimensional（1D） fluid model of capacitive RF argon glow discharges between two parallel-plate electrodes at low pressure is employed. The influence of the secondary electron emission on the plasma characteristics in the discharges is investigated numerically by the model. The results show that as the secondary electron emission coefficient increases,the cycle-averaged electric field has almost no change; the cycle-averaged electron temperature in the bulk plasma almost does not change, but it increases in the two sheath regions; the cycle-averaged ionization rate, electron density, electron current density, ion current density, and total current density all increase. Also, the cycle-averaged secondary electron fluxes on the surfaces of the electrodes increase as the secondary electron emission coefficient increases. The evolutions of the electron flux, the secondary electron flux and the ion flux on the powered electrode increase as the secondary electron emission coefficient increases. The cycle-averaged electron pressure heating, electron Ohmic heating, electron heating, and ion heating in the two sheath regions increase as the secondary electron emission coefficient increases. The cycle-averaged electron energy loss increases with increasing secondary electron emission coefficient.

Geodesic acoustic mode in a reduced two-fluid model

A reduced two-fluid model is constructed to investigate the geodesic acoustic mode（GAM）. The ion dynamics is sufficiently considered by including an anisotropic pressure tensor and inhibited heat flux vector, whose evolutions are determined by equations derived from the 16-momentum model. Electrons are supposed to obey the Boltzmann distribution responding to the electrostatic oscillation with near ion acoustic velocity. In the large safety factor limit, the GAM frequency is identical with the kinetic one to the order of 1 q2 when zeroing the anisotropy. For general anisotropy, the reduced two-fluid model generates the frequency agreeing well with the kinetic result with arbitrary electron temperature. The present simplified fluid model will be of great use and interest for young researchers and students devoted to plasma physics.

Investigation of the diocotron instability of an infinitely wide sheet electron beam by using the macroscopic cold-fluid model theory

This paper investigates the diocotron instability of an infinitely wide relativistic sheet electron beam in conducting wMls propagating through a uniform magnetic field by using the macroscopic cold-fluid model theory. Assuming low- frequency perturbations with long axial wavelengths, the eigenvalue equation and the dispersion relation are acquired for a sheet electron beam with sharp boundary profile and uniform density. The results presented in this paper has developed the use of the macroscopic cold-fluid model theory by extending the parameter of the electron cyclotron frequency ωc to a wider usage range, which is restricted to be much larger than the plasma frequency ωp in the previous research work. Theoretical analyses and numerical calculations indicate that the transport of the sheet electron beam will be completely stabilized by augmenting the normalized beam thickness to a conductor gap larger than a threshold λb, which is greatly dependent on the parameter ωc/ωp. The larger ωc/ωp is, the smaller λb will be needed. Moreover, the system parameters, including the wave number kx of the perturbations and the relativistic mass factor γb, will also influence the growth rate of diocotron instability obviously.

Investigation of intense sheet electron beam transport using the macroscopic cold-fluid model and the single-particle orbit theory

The focusing and the stable transport of an intense elliptic sheet electron beam in a uniform magnetic field are investigated thoroughly by using the macroscopic cold-fluid model and the single-particle orbit theory.The results indicate that the envelopes and the tilted angles of the sheet electron beam obtained by the two theories are consistent.The single-particle orbit theory is more accurate due to its treatment of the space-charge fields in a rectangular drift tube.The macroscopic cold-fluid model describes the collective transport process in order to provide detailed information about the beam dynamics,such as beam shape,density,and velocity profile.The tilt of the elliptic sheet beam in a uniform magnetic field is carefully studied and demonstrated.The results presented in this paper provide two complete theories for systemically discussing the transport of the sheet beam and are useful for understanding and guiding the practical engineering design of electron optics systems in high power vacuum electronic devices.

Numerical Simulation of Innovative Operation of Blast Furnace Based on Multi-Fluid Model

A multi-fluid blast furnace model was simply introduced and was used to simulate several innovative ironmaking operations. The simulation results show that injecting hydrogen bearing materials, especially injecting natural gas and plastics, the hydrogen reduction is enhanced, and the furnace performance is improved simultaneously. Total heat input shows obvious decrease due to the decrease of heat consumption in direct reduction, solution loss and silicon transfer reactions. If carbon composite agglomerates are charged into the furnace, the temperature of thermal reserve zone will obviously decrease, and the reduction of iron-bearing burden materials will be retarded. However, the efficiency of blast furnace is improved just due to the decrease in heat requirements for solution loss, sinter reduction, and silicon transfer reactions, and less heat loss through top gas and furnace wall. Finally, the model is used to investigate the performance of blast furnace under the condition of top gas recycling together with plastics injection, cold oxygen blasting and carbon composite agglomerate charging. The lower furnace temperature, extremely accelerated reduction rate, drastically decreased CO2 emission and remarkably enhanced heat efficiency were obtained by using the innovative operations, and the blast furnace operation with superhigh efficiency can be realized.

Exact Solutions of a Power Law Fluid Model in Posttreatment Analysis of Wire Coating with Linearly Varying Boundary Temperature

In this paper, analysis of post-treatment of wire coating is presented. Coating material satisfies power law fluid model. Exact solutions for the velocity field, volume flow rate and average velocity are obtained. Moreover, the heat transfer results are presented for different cases of linearly varying on the boundaries. The variations of velocity, volume flow rate, radius of coated wire, shear rate and the force on the total wire are presented graphically and discussed.

CONTINUOUS THERMODYNAMICS FOR POLYMER SOLUTIONS II . LATTICE FLUID MODEL

Using lattice-fluid model, a continuous thermodynamic framework is presented for phase-equilibrium calculations for binary solutions with a polydisperse polymer solute. A two-step process is designed to form a real polymer solution containing a solvent and a polydisperse polymer solute occupying a volume at fixed temperature and pressure. In the first step, close-packed pure components including solvent and polymers with different molar masses or different chain lengths are mixed to form a closed-packed polymer solution. In the second step, the close-packed mixture, considered to be a pseudo-pure substance is mixed with holes to form a real polymer solution with a volume dependent on temperature and pressure. Revised Freed's model developed previously is adopted for both steps. Besides pure-component parameters, a binary size parameter cr and a binary energy parameter e12 are used. They are all temperature dependent. The discrete-multicomponent approach is adopted to derive expressions for chemical potentials, spinodals and critical points. The continuous distribution function is then used in calculations of moments occurring in those expressions. Computation procedures are established for cloud-point-curve, shadow-curve, spinodal and critical-point calculations using standard distribution or arbitrary distribution on molar mass or on chain length. Illustrative examples are also presented.

Asymptotic Stability of Combination of Contact Discontinuity with Rarefaction Waves for the One-Dimensional Viscous Micropolar Fluid Model

In this paper, we consider with the large time behavior of solutions of the Cauchy problem to the one-dimensional compressible micropolar fluid model, where the far field states are prescribed. When the corresponding Riemann problem for the Euler system admits the solution consisting of contact discontinuity and rarefaction waves, it is shown that the combination wave corresponding to the contact discontinuity, with rarefaction waves is asymptotically stable provided that the strength of the combination wave and the initial perturbation are suitably small. This result is proved by using elementary L2-energy methods.

STRONG SOLUTIONS FOR THE INCOMPRESSIBLE FLUID MODELS OF KORTEWEG TYPE

In this article, we are concerned with the strong solutions for the incompress- ible fluid models of Korteweg type in a bounded domain Ω СR^3. We prove the existence and uniqueness of local strong solutions to the initial boundary value problem. We point out that in this article we allow the existence of initial vacuum provided initial data satisfy a compatibility condition.

Fluid Model Modulated by an M/M/1 Working Vacation Queue with Negative Customer

This paper investigates a fluid model driven by an M/M/1 queue with working vacations and RCE （Removal of customer in the end） policy of negative customer. In the external environment, the negative customer is not served by the server and only removes the positive customer in the end one-to-one. We establish a fluid flow model based on this stochastic process, and obtain the mean buffer content and the probability of empty buffer for this fluid queue using the LT （Laplace transform） method. Moreover, several special cases of the model here are obtained. Finally, some numerical examples are presented to demonstrate the effects of parameters on the performance indices of the fluid model.

Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.

Two-fluid modeling of Geldart A particles in gas-fluidized beds

We have investigated the effect of cohesion and drag models on the bed hydrodynamics of Geldart A particles based on the two-fluid （TF） model. For a high gas velocity U0 = 0.03 m/s, we found a transition from the homogeneous fluidization to bubbling fluidization with an increase of the coefficient C1, which is used to account for the contribution of cohesion to the excess compressibility. Thus cohesion can play a role in the bed expansion of Geldart A particles. Apart from cohesion, we have also investigated the influence of the drag models. When using the Wen and Yu drag correlation with an exponent n = 4.65, we find an under-prediction of the bed expansion at low gas velocities （U0 = 0.009 m/s）. When using a larger exponent （n = 9.6）, as reported in experimental studies of gas-fluidization, a much better agreement with the experimental bed expansion is obtained. These findings suggest that at low gas velocity, a scale-down of the commonly used drag model is required. On the other hand, a scale-up of the commonly used drag model is necessary at high gas velocity （U0 = 0.2 and 0.06 m/s）. We therefore conclude that scaling the drag force represent only an ad hoc way of repairing the deficiencies of the TF model, and that a far more detailed study is required into the origin of the failure of the TF model for simulating fluidized beds of fine powders.