Numerical simulation for space charge dominated beam transport

To simulate the intense bunched beam transport, a computer program LEADS-3D has been developed. The particle trajectories are analysed with the Lie algebraic method. The third order approximation of the trajectory solutions is made with space charge forces off, and the second order approximation is made with space charge forces on. The particle distribution in the 3D ellipsoid is uniform or Gaussian. Most of the conventional beam optical elements are incorporated in the code. The optimization procedures are provided to fit the beam lines to satisfy the given optical conditions.     

Lie algebraic analysis for the nonlinear transportof intense pulsed beams in electrostatics lenses

The Lie algebraic method is applied to the analysis of the nonlinear transport of an intense pulsed beam in cylindrically symmetrical electrostatic lenses, and particle orbits in a six-dimensional phase space (x, p_x, y, p_y, τ, p_τ) are obtained in the second order approximation. They can also be acquired in the third or higher order approximation if needed. In the analysis, we divide the electrostatic lenses into several segments. Each segment is considered as a uniform accelerating field, and each dividing point is treated as a thin lens. The particle distribution in a three-dimensional ellipsoid is of Gaussian type.     

Thermodynamic Properties of Random Transverse Field Mixed Spin System in the Presence of Single-Ion Anisotropy

We study the thermodynamic properties of random transverse field mixed spin system in the presence of single-ion anisotropy on a square lattice. By making use of the effective field theory and a cutting approximation, the detailed phase diagrams are described and some interesting results are found under trimodal random transverse field distribution. A small single-ion anisotropy can magnify magnetic ordering region at low temperatures and existence of a large transverse field can assist the occurrence of reentrant phenomena. With increasing disorder, second-order phase transitions are shown to change into first-order phase transitions. The trajectory of the tricritical point in the phase space as a function of disorder is presented. These indicate a strong correlation with the corresponding to trimodal transverse field distribution.     

Analytical formula of free electron laser exponential gain for a non-resonant electron beam

The free electron laser(FEL) gain formulas for a non-resonant case are studied, and some new rigorous analytical formulas are given explicitly. For the mono-energetic and non-resonant electron beam, the exact expression of the solution of the FEL characteristic cubic equation is obtained with a form much more simple than that in the literatures, and the gain length as the function of the detuning parameter is explicitly given. Then the gain for different detuning parameters and from low to high can be easily calculated. A simplified approximation formula is also given for the exponential gain calculation in the non-resonant case. For the case of the electron beam with an energy spread, the solution of the characteristic cubic equation is given explicitly for rectangular energy distribution and Lorentz distribution, respectively. Moreover the explicit expression also can be used for the solution of the characteristic cubic equation including the impact of the space charge. The transition from the low gain to the high gain is analyzed. The variations of the gain bandwidth and of the detuning parameter for the maximum gain are demonstrated. The applicable ranges of the small signal gain formula and the exponential gain formula are analyzed.     

A modal analysis for the acoustic radiation problems, I. Theory

For the acoustic radiation problems from a complex vibrating body surface, a modal analysis approach is put forward. All the normal vibration velocities on a vibrating surface form the Hilbert space. In the Hilbert space, an operator is defined, which includes the radiation property of the vibrating surface and is linear, self-adjoint and positive. Using the operator, a set of basis functions in the Hilbert space are obtained, which describe the radiation patterns and are called the radiation modes. Based on the radiation modes, a set of basis functions of the radiation field are obtained by the Helmholtz simple layer potentials, which describe the distribution patterns of the radiation field and are called the field distribution modes. The radiation behavior can be expressed by expansions of the radiation modes and the field distribution modes. The modal analysis approach is introduced into the acoustic radiation problems.     

A Simple Calculation Approach for the Second－Harmonic Sound Beam Generated by an Arbitrary Distribution Source

We present a simple calculation approach for the fundamental and second-harmonic sound beams with an arbitrary distribution source in the quasilinear approximation. The analysis is based on the assumption that the source function with an arbitrary geometry and distribution is expanded into the sum of a set of two-dimensional Gaussian functions. The two- and five-dimensional integral solutions for the fundamental and second-harmonic fields are, respectively, reduced in terms of Gaussian functions and simple one-dimensional integrals. The numerical evaluation of field distributions is then greatly simplified.     

Quantum Dynamics of Cooled Atoms in the Presence of Bose—Einstein Condensates

Under the Markov approximation,the quantum dynamics of cooled atoms in the presence of Bose-Einstein condensates is studied.A master equation governing the evolution of such a system is derved.Using this master equation,the distribution of the atoms in the excited states at finite temperature and the dynamics of the excited atom at zero temperature are given and discussed.     

Frequency/time domain modeling of a direct drive point absorber wave energy converter

Two different methods to model a point absorber wave energy converter(WEC)with direct drive linear power take-off(PTO)are proposed in the present study:the frequency domain(FD)method and the time domain(TD)method.In the FD analysis,the frequency response function(FRF)of the WEC device is obtained via the equation of motion,and the expressions of power capture width in regular and random waves are derived as well.In the TD modeling,based on a state space approximation of the convolution term in the motion equation,both regular wave and random wave simulations are carried out.The regular wave simulation results indicate that the state space approximation is sufficiently accurate and the capture width reaches the maximum in the vicinity of the natural frequency.In the random wave simulations,the effects of buoy size,the PTO damping and wave climate on the power capture width are discussed in detail,which leads to the conclusion that the capture widths are influenced by the natural frequency of the WEC device,peak frequency of the wave spectrum,the amplitude of FRF and PTO damping.Furthermore,the increase of the capture width is at the cost of a relatively large buoy size and PTO damping when control is not included.     

Why Do Deconfinement and Chiral Restoration Coincide at High Temperature？

The global colour model in free space is extended to finite temperature to study the deconfinement and chlral phase transitions at high temperature T and zero chemical potential in the mean field approximation. Both possibilities of coincidence and non-coincidence of the two distinct phase transitions are found when the model     

Quasispecies distribution of Eigen model

We have studied sharp peak landscapes of the Eigen model from a new perspective about how the quasispecies are distributed in the sequence space. To analyse the distribution more carefully, we bring in two tools. One tool is the variance of Hamming distance of the sequences at a given generation. It not only offers us a different avenue for accurately locating the error threshold and illustrates how the configuration of the distribution varies with copying fidelity $q$ in the sequence space, but also divides the copying fidelity into three distinct regimes. The other tool is the similarity network of a certain Hamming distance $d_{0}$, by which we can gain a visual and in-depth result about how the sequences are distributed. We find that there are several local similarity optima around the centre (global similarity optimum) in the distribution of the sequences reproduced near the threshold. Furthermore, it is interesting that the distribution of clustering coefficient $C(k)$ follows lognormal distribution and the curve of clustering coefficient $C$ of the network versus $d_{0}$ appears to be linear near the threshold.     

Lie algebraic analysis for the nonlinear transport of intense bunched beam in electrostatic quadrupoles

In this paper, the nonlinear transport of intense bunched beams in electrostatic quadrupoles is analyzed using the Lie algebraic method, and the results are briefly presented of the linear matrix approximation and the second order correction of particle trajectory in the state space. Beam having K-V distribution and Gaussian distribution approximation are respectively considered. A brief discussion is also given of the total effects of the quadrupole and the space charge forces on the evolution of the beam envelope.     

离散速度方向模型在微尺度气体流动中的应用

离散速度方向模型是一种简化Boltzmann方程的新方法。该方法通过减少Boltzmann方程的维数来降低数值求解的计算量。在DVD模型中,分子速度的方向是离散的,而分子的速率仍然是连续的,这样就可以用一组三维的速率分布函数来代替Boltzmann方程中六维的速度分布函数。由于减少了三个动量维,同Boltzmann方程相比,DVD模型的数值计算量可以降低几个数量级。本文用数值的方法对DVD模型进行了研究。数值结果显示,在广泛的Knudsen数下,DVD方法可给出精确的计算结果。同线性化Boltzmann方程的计算结果相比,最大的误差不超过6%,在连续介质领域中,误差甚至不超过1%。     

Non-equilibrium approach to transport phenomena in quantum crystals

A pure dielectric quantum crystal subjected to an external mechanical force is described by non-equilibrium Green’s functions. In equilibrium the leading approximation leads to the definition of elementary excitations, the phonons in the renormalized harmonic approximation. Their temperature dependent energies are to be determined as solutions of an integral equation. For hydrodynamic disturbances a generalized transport equation for a phonon number density is derived. A similar approximation for the spectral function yields an integral equation for space and time dependent quasiparticle energies which are expressed as functionals of the displacement field and the phonon distribution. The Boltzmann equation for the latter includes the quasi-particle interaction.     

Propagation properties of vectorial elliptical Gaussian beams beyond the paraxial approximation

Starting from the vectorial Rayleigh diffraction integral formula and without using the far-field approximation, a solution of the wave equation beyond the paraxial approximation is found, which represents vectorial non-paraxial elliptical Gaussian beams in free space. The far-field expressions for non-paraxial Gaussian beams and elliptical Gaussian beams can be regarded as special cases treated in this paper. Some basic propagation properties of vectorial non-paraxial elliptical Gaussian beams, including the irradiance distribution, phase term, beam widths and divergence angles are studied. Numerical results are given and illustrated.     

Classical projected phase space density of billiards and its relation to the quantum neumann spectrum

A comparison of classical and quantum evolution usually involves a quasiprobability distribution as a quantum analogue of the classical phase space distribution. In an alternate approach that we adopt here, the classical density is projected on to the configuration space. We show that for billiards, the eigenfunctions of the coarse-grained projected classical evolution operator are identical to a first approximation to the quantum Neumann eigenfunctions. However, even though there exists a correspondence between the respective eigenvalues, their time evolutions differ. This is demonstrated numerically for the stadium and lemon-shaped billiards.     

The Δ- Method for the Boltzmann Equation

A new numerical method is proposed to solve the Boltzmann equation. A frame is set up by using a discrete velocity approximation in the infinite velocity space, but by considering only those distribution function points which are not too small. The distribution function points may occur anywhere in the infinite discrete velocity space and are not constrained to a pre-specified region. A fourth-order finite difference is used for the convection terms. A Monte Carlo-like method is applied to the discrete velocity model of the collision integral. The effort of the method is proportional to the number of discrete points. Numerical examples are given for the full Boltzmann equation and results for some benchmark problems are compared with analytical or prior solutions.     

An identity for the non-equilibrium phase space distribution function

A new identity satisfied by the phase space distribution function is derived and the Zubarev statistical operator is obtained as an approximation to this identity.     

电推力器流动模拟中的电子处理方法

 分析了电子的准中性假设、玻耳兹曼分布假设、粒子模型在电推力器流动模拟中的适用性和优劣性,提出了一种新的电子处理方法——电子漂移扩散近似,采用该方法模拟了离子发动机栅极光学系统等离子体运动过程。结果表明:该方法得出的电势分布、离子相空间分布及电子数密度分布与经典的电子玻耳兹曼分布假设处理方法计算结果一致,验证了该方法可以很好地应用于电推力器栅极光学系统模拟。     

Fundamentals of a Technique for Determining Electron Distribution Functions by Multi-Term Even-Order Expansion in Legendre Polynomials. 2. Comprehensive Investigation of a Model Plasma

Applying the new technique for finding the converged solution of the Boltzmann equation in a weakly ionized plasma, which was developed in the first part of this paper, a comprehensive study of the electron velocity distribution function for a model plasma with elastic and exciting collisions is performed by solving the Boltzmann equation with increasing order of approximation. The purpose of this investigation is that of calculating the isotropic distribution f₀, the first contribution f₁ to the anisotropy of the velocity distribution, the important macroscopic quantities and, more generally, that of studying the total anisotropy as well as the changes of all these quantities when the approximation degree is enlarged beyond the 2 terms of the conventional Lorentz approximation. By varying some parameters of the model plasma, that is the electric field strength, the magnitude of the excitation cross section and the excitation threshold, the main features of plasmas in inert as well as molecular gases are modelled and the impact of these parameters on the mentioned quantities is analysed. Some of the converged results are compared with results of corresponding Monte Carlo simulations. The approximation degree required to find the converged values of isotropic distribution, main macroscopic quantities and electron distribution in the velocity space (and thus its real anisotropy) is estimated by solving the Boltzmann equation over wide parameter ranges.