湍动雾化射流液雾粒径分布的数值模拟

对水在空气中湍动雾化射流的气液两相流场进行了数值模拟.其中,气体流场采用k-ε湍流模型进行模拟,给出了载气轴向速度的分布情况.对喷雾粒子的运动采用颗粒轨道法,建立粒子破碎和碰撞的数值模型,研究了液雾粒径在不同工况下沿轴向的变化趋势.数值模拟结果与实验结果在多种气液比下进行了比较,两者吻合较好.同时,分析了初始粒径和粒子总数及喷嘴尺寸对液雾粒径变化趋势的影响,讨论了有助于粒径均匀分布的条件.     

HL-1M装置弹丸加料等离子体中的蛇形振荡

在弹丸注入实验中，人们观测到穿透至q=1面以内的弹丸能在等离子体芯部激发出蛇形振荡。蛇形振荡是由一个局域的高密度区域在等离子体芯部旋转而反映在软X射线、密度信号上的一种振荡现象。在JET，ASDEX，Tore Supra等装置上都观测到了弹丸注入条件下的蛇形振荡现象。     

混合层中柱状粒子运动的研究

本文在直接数值模拟混合层流场的基础上,求解柱状粒子运动方程而得到粒子运动轨迹,说明粒子在混合层中的不同区域,其运动轨迹也不同,有的呈现波状轨迹,有的呈现环状轨迹,有的是混合轨迹。粒子的混合程度与St数有关。文中结论对工程应用有指导意义。     

Orientation distribution of fibres in a channel flow of fibre suspension

The orientation and concentration distributions of fibres in laminar and turbulent channel flows are investigated numerically. The obtained results are in good agreement with the experimental data. In the laminar flow regime, more fibres orient to the flow direction as the Reynolds number increases. The shear rate of fluid around a fibre plays an important role in determining the orientation distribution of fibres, while the fibre density and the fibre aspect-ratio have marginal influence on the orientation distribution. In the turbulent regime, the orientation distribution of fibres becomes more homogeneous with the increase of Reynolds number, and the concentration profile is flatter than that in the laminar regime. The fluctuating intensity of fibre velocity in the downstream direction is larger than that in the lateral directions.     

70MHz连续波质子束脉冲化装置

为了进行强流回旋加速器关键技术研究, 中国原子能科学研究院建立了一个强流回旋加速器综合试验装置. 中国原子能研究院将在这个回旋加速器综合试验装置上建立强流脉冲化实验装置, 目标是实现几十至百keV量级的强流束的脉冲化. 具体是将70MHz连续波负氢束脉冲化为重复频率1—8MHz, 脉冲宽度约为10ns的脉冲质子束. 脉冲化装置将主要包括束流切割器和聚束器两大系统. 聚束器采用频率为70.487MHz的双间隙单漂移栅网结构, 可以将直流束压缩到±30°的回旋加速器高频接受相宽之内. 束流切割器将采用频率为2.2MHz的正弦波, 切割后的脉冲宽度将小于8ns, 最后得到的脉冲束的重复频率为4.4MHz.     

Absence of Charge-Resonance-Enhanced Ionization in Attosecond Pulse Photoionization： Numerical Result on One-Dimensional H2^＋

We present a numerical result of photoionization rate for the one-dimensional molecular hydrogen ion model exposed to intense light of 1 × 10^16-2×10^16 W/cm^2, 55-as pulse duration, and 800nm wavelength. In contrast to the previous calculation result of charge-resonance-enhanced ionization for lower intensity and much longer pulse, our result exhibits an ionization saturation. The numerical results are interpreted in the field-dressed potential picture as over-the-barrier liberation of electrons. This extremely short pulsewidth and relatively high field phenomenon requests experimental demonstration.     

Behaviour of Self-Standing CVD Diamond Film with Different Dominant Crystalline Surfaces in Thermal-Iron Plate Polishing

Self-standing CVD diamond films with different dominant crystalline surfaces are polished by the thermal-iron plate polishing method. The influence of the dominant crystalline surfaces on polishing etfficiency is investigated by measuring the removal rate and final roughness. The smallest rms roughness of 0.14 μm is measured with smallest removal rate in the films with the initial （220） dominant crystalline surface. Activation energy for the polishing is analysed by the Arrhenius relation. It is found that the values are 170kJ/mol, 222kJ/mol and 214kJ/mol for the film with three different dominant crystalline surfaces. Based on these values, the polishing cause is regarded as the graphitization-controlling process. In the experiment, we find that transformation of the dominant crystalline surfaces from （111） to （220） always appears in the polishing process when we polish the （111） dominant surface.     

The most common friend first immunization

In this paper, a standard susceptible-infected-recovered-susceptible(SIRS) epidemic model based on the Watts–Strogatz(WS) small-world network model and the Barabsi–Albert(BA) scale-free network model is established, and a new immunization scheme — "the most common friend first immunization" is proposed, in which the most common friend's node is described as being the first immune on the second layer protection of complex networks. The propagation situations of three different immunization schemes — random immunization, high-risk immunization, and the most common friend first immunization are studied. At the same time, the dynamic behaviors are also studied on the WS small-world and the BA scale-free network. Moreover, the analytic and simulated results indicate that the immune effect of the most common friend first immunization is better than random immunization, but slightly worse than high-risk immunization. However,high-risk immunization still has some limitations. For example, it is difficult to accurately define who a direct neighbor in the life is. Compared with the traditional immunization strategies having some shortcomings, the most common friend first immunization is effective, and it is nicely consistent with the actual situation.     

Characteristics of Flow around an Impulsively Rotating Square Cylinder via LB-DF/FD Method

Flow around an impulsively rotating square cylinder in a viscous fluid in range of 1 ≤ Re ≤ 300 is numerically investigated by the previously developed LB-DF/FD method, which combines the lattilce Boltzmann method （LBM） and direct-forcing fictitious domain （DF/FD） scheme. Results show that in total three kinds of transient characteristics depending on Re are observed： 1 ≤ Re ≤ 20, four vortices arising from the corners of the square cylinder separate from the surfaces and gradually become stable without vortex integration or shedding; 20 ≤ Re 〈 100, vortices integration is observed when they grow long enough, then separated from each other; Re ≥ 100, vortex shedding takes place in this regime. The shedding vortex joins the downstream vortex to form a new one. It is also found that vortex shedding happens more than one time when Re ≥ 160. Furthermore, each vortex shedding induces a fluctuation in the torque exerted on the cylinder.     

Stripping extraction calculation and simulation for CYCIAE-100

A 100 MeV H- compact cyclotron is under construction at China Institute of Atomic Energy (CYCIAE-100). The proton beams of 75 MeV-100 MeV at an intensity of 200 μA will be extracted in dual opposite directions by charge exchange stripping devices. The crossing point at the switching magnet center is fixed inside the magnet yoke and the stripping points for various extraction energies are calculated by the code CYCTRS. With the code GOBLIN, we can calculate the transfer matrix including the dispersion effects from the stripping point to the switch magnet. The beam distribution just after stripping foil can be obtained from the multi-particle tracking code COMA and the extracted beam parameters after the switch magnet such as emittance, envelope, dispersion, energy spread, bunch length, etc. are given by the extraction orbit simulations.