Numerical Simulation Multicomponent Ion Beam Transport form ECR Ion Source

In order to simulate the transport of multi-components ion beam extracted from an ECR ion source, we have developed a multi-charged ion beam transport program named MCIBS 1.0. The program is dedicated to numerical simulation of the behavior of highly-charged ion beam and optimization of beam optics in transport lines and is realized on a PC with Windows user interface of Microsoft Visual Basic. Among all the ions with different charge states in the beam, the exchanges of electrons between highly charged ions and low charged ions or neutral,atoms of residual gas are taken into account by using classical Molecular Over-barrier Model and Monte Carlo method. An advanced Windows graphical interface makes it; comfortable and friendly for the user to operate in an interactive mode. The present program is used for the numerical calculation and optimization of beam optics in a transport line consisting of various magnetic elements, such as dipole magnet, quadrupole and so on. It is possible to simultaneously simulate 200,000 particles, in a transport line of 340 m at most, and show every particle orbit. Beam cross section graphics and emittance phase pictures can be also shown at any position in the transport line.     

Numerical Simulation on Expansion Process of Ablation Plasma Induced by Intense Pulsed Ion Beam

We present a one-dimensional time-dependent numerical model for the expansion process of ablation plasma induced by intense pulsed ion beam （IPIB）. The evolutions of density, velocity, temperature, and pressure of the ablation plasma of the aluminium target are obtained. The numerical results are well in agreement with the relative experimental data. It is shown that the expansion process of ablation plasma induced by IPIB includes strongly nonlinear effects and that shock waves appear during the propagation of the ablation plasma.     

Numerical Simulation of MHD Jet Flows

Liquid metal free surface flows （films, jets and droplets） are considered as diverter/ limiter system and first wall in fusion reactor, but the knowledge Of liquid metal free surface under a non-uniform magnetic field is very limited. In this article, the stability of a jet flow under a gradient magnetic field is investigated, and its MHD effects are the top concern. Based on numerical simulation and experimental results, a simplified model is developed to analyze the MHD effects of the jet flow and to explain the reason why it can keep stable under a strong non-uniform magnetic field.     

Interference Patterns of three-fiber Mach-Zehnder Interferometer and its Numerical Simulation

１ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｉｂｅｒｏｐｔｉｃＭａｃｈＺｅｈｎｄｅｒｉｎｔｅｒｆｅｒｏｍｅｔｅｒｉｓａｗｅｌｅｓｔａｂｌｉｓｈｅｄａｎｄｗｉｄｅｌｙｕｓｅｄｄｅｖｉｃｅ．Ｔｈｅｅｘｐｅｒｉｍｅｎｔａｌｂａｃｋｇｒｏｕｄｏｆｆｉｂｅｒｏｐｔｉｃ...     

Direct Numerical Simulation of Acoustic Turbulence: Zakharov–Sagdeev Spectrum

JETP Letters - We present the results of direct numerical simulation of three-dimensional acoustic turbulence in medium with weak positive dispersion. It is shown that at the beginning of the...     

Numerical Molecular Dynamics Simulation of the Fracture of a Ti–Al Intermetallic Nanocrystal

Doklady Physics - The results of numerical molecular dynamics simulation of the synthesis and fracture of a Ti–Al intermetallic nanocrystal under uniaxial tension are reported. It is shown...     

Lattice Boltzmann Numerical Simulation of a Circular Cylinder

The lattice Boltzmann equation (LBE) model based on the Boltzmann equation is suitable for the numerical simulation of various flow fields.The fluid dynamics equation can be recovered from the LEB model.Howeverl,compared to the Navier-Stokes transport equation,the fluid dynamics equation derived from the LBE model is somewhat different in the viscosity transport term,which contains not only the Navier-Stokes transport equation but also nonsteady pressure and momentum flux terms.The two nonsteady terms can produce the same function as the random stirring force term introduced in the direct numerical or large-eddy vortex simulation of turbulence.Through computation of a circular cylinder,it is verified that the influence of the two nonsteady terms on flow field stability cannot be ignored,which is helpful for the study of turbulence.     

Magnetic and mechanical properties of Ni–Mn–Ga/Fe–Ga ferromagnetic shape memory composite

A ferromagnetic shape memory composite of Ni–Mn–Ga and Fe–Ga was fabricated by using spark plasma sintering method. The magnetic and mechanical properties of the composite were investigated. Compared to the Ni–Mn–Ga alloy,the threshold field for magnetic-field-induced strain in the composite is clearly reduced owing to the assistance of internal stress generated from Fe–Ga. Meanwhile, the ductility has been significantly improved in the composite. A fracture strain of 26% and a compressive strength of 1600 MPa were achieved.     

Laser and Thermal Modification of Silver–Ion Implanted Glasses

We investigate sodium–calcium silicate glasses implanted with 60 keV Ag⁺ ions with a dose of 3·10¹⁶ cm^–2 at an ion current density of 10 A/cm². As a result of the ion implantation, a composite layer with silver nanoparticles is synthesized in the region near the surface. However, this layer is characterized by high nonuniformity in the size distribution of these particles over the depth of the layer. Subsequent pulsed laser irradiation in combination with equilibrium heat treatment makes it possible to modify this composite layer, improving the uniformity in the size distribution of the nanoparticles. This is particularly promising for the improvement of the technology of obtaining nonlinear optical materials. To control the parameters of the layers obtained we suggest a method based on an analysis of the optical reflection and transmission spectra measured on the side of the implanted and opposite surfaces of glass samples.     

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M∞ = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tollmien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.     

Numerical Simulation of Solitary Kinetic Alfvén Waves

Using the two-fluid model in the case of α〉〉1 （α= β/2Q, β is the ratio of thermal pressure to magnetic pressure, and Q = me/mi ）, we numerically investigate the interactions between two solitary kinetic Alfvén waves （SKAWs） and between an SKAW and a density discontinuity. The results show that the two SKAWs would remain in their original shapes and propagate at their initiating speeds, which indicates that SKAWs behave just like standard solitons. The simulation also shows that SKAWs will reflect and refract when crossing a discontinuity and propagating into a higher density region. The transmission wave is an SKAW with increasing density, and the reverberation is a disturbance with lower amplitude.     

Numerical Simulation and Analysis on Deflection of Silicon Mirrors Irradiated by High Power Laser

ＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎａｎｄＡｎａｌｙｓｉｓｏｎＤｅｆｌｅｃｔｉｏｎｏｆＳｉｌｉｃｏｎＭｉｒｒｏｒｓＩｒｒａｄｉａｔｅｄｂｙＨｉｇｈＰｏｗｅｒＬａｓｅｒＸＩＡＪｉｎａｎ；ＣＨＥＮＧＺｕｈａｉ；ＱＩＵＪｕｎｌｉｎ（ＮａｔｉｏｎａｌＬａｂ...     

Numerical Simulation of the Movement of Frenkel–Kontorova Dislocations in Aluminum Single Crystals at Low Temperatures

The motion of Frenkel–Kontorova dislocations in the single crystals of aluminum at low temperatures has been studied, by means of the computer simulation. It is shown that the dislocation movement is realized by the quantum tunneling of the kinks of dislocations through the Peierls barriers. It is shown that the action of the Peierls high barrier is analogous to the action of low temperatures, and if the Peierls barrier overcome, the dislocation moves unevenly, accelerating under the action of the Peierls barrier and slowing down after overcoming the Peierls barrier. Based on the numerical experiment, the mean free path of dislocation, the distance between the Peierls potential barriers and the width of the Peierls barrier are calculated. The computed values correspond to the real values.     

Ion Cleaning of Facets for Improving the Reliability of High Power 980 nm Semiconductor Lasers

We report a simple and available way of improving the reliability of high power InGaAs 980 nm lasers by cleaning the facets using Ar ion before the protecting films have been coated. The Ar cleaning can remove the impurity and the oxide on the air-cleaved facets of laser diodes. It is proven that the way has marked effect on reducing the gradual degradation rate of laser diodes and improving the catastrophic-optical-damage threshold.     

Thermodynamics and surface properties of liquid Al–Ga and Al–Ge alloys

The surface properties of Al–Ga and Al–Ge liquid alloys have been theoretically investigated at a temperature of 1100 K and 1220 K respectively. For the Al–Ga system, the quasi chemical model for regular alloy and a model for phase segregating alloy systems were applied, while for the Al–Ge system the quasi chemical model for regular and compound forming binary alloys were applied. In the case of Al–Ga, the models for the regular alloys and that for the phase segregating alloys produced the same value of order energy and same values of thermodynamic and surface properties, while for the Al–Ge system, the model for the regular alloy reproduced better the thermodynamic properties of the alloy. The model for the compound forming systems showed a qualitative trend with the measured values of the thermodynamic properties of the Al–Ge alloy and suggests the presence of a weak complex of the form Al₂Ge₃. The surface concentrations for the alloys show that Ga manifests some level of surface segregation in Al–Ga liquid alloy while the surface concentration of Ge in Al–Ge liquid alloy showed a near Roultian behavior below 0.8 atomic fraction of Ge.     

Numerical Simulation of Ultrafast Laser Pulse Propagation in Tenuous Plasmas： Envelope Evolving and Modulation

We propose an effective and useful numerical simulation scheme for the investigation of the ultra-fast laser pulses in tenuous plasmas. The accuracy of the method is tested by numerical examples. We check some special examples to investigate the laser envelope evolving and modulation in plasmas. Asymmetric two-peak modulation structure is found and its underlying physics is analyzed. The advantages and shortages of the method are also discussed.     

Molecular Dynamics Simulation of Ion Implantation into HfO2 and HfO2／Si Multi-Layer Structure

We simulate the B, As and P implantations into HfO2 from 3keY to 40keY by a simulator LEACS developed based on molecular dynamics method and by the traditional Monte Carlo simulator TSUPREM4 respectively. The LEACS results accurately fit with the SIMS (secondary ion mass spectroscopy) data, while the TSUPREM4 results deviate from the SIMS data obviously except B implantation. Based on the verification of the simulator, influence of the oxide thickness on the retained range profiles in the Si layer has been quantitatively investigated in the case of HfO2/Si and SiO2/Si structures. The range profiles in the Si layer through HfO2 shift to the surface obviously for about 0.68 times of the oxide layer thickness on the average in comparison to those through SiO2. It can be predicted that this effect will have a significant impact on MOSFET (metallic oxide semiconductor field effect transistor) device performance in the integrated circuit process of the next decade if HfO2 is used to replace SiO2 as the gate dielectric.     

Spectroscopy of He–Ne afterglow plasma

Implementation of promising control schemes for the intensity and position of X-ray-laser beams with a photon energy up to several tens of kiloelectronvolts requires knowledge of the angular dependence of cross sections for photoionization of noble gas atoms by hard photons. Estimates of quadrupole corrections to the cross section for photoionization of a Kr atom by X-ray photons with an energy of about 25 keV are reported in this paper. An analytic expression for the cross section of the process is parameterized in a compact form convenient for analyzing angular distributions with an arbitrary polarization of a photon beam.     

Numerical Simulation of MHD Equilibrium Configuration for HL-2M Tokamak

In this work, the numerical simulation is employed for HL-2A equipment modification project HL-2M. The design of zero field and evolution of plasma configuration in start-up phase and single null divertor equilibrium configurations in steady-state phase has been simulated. Results show that new project can satisfy request of design.