Analytical solution of a multidimensional Langevin equation at high friction limits and probability passing over a two-dimensional saddle

The analytical solution of a multidimensional Langevin equation at the overdamping limit is obtained and the probability of particles passing over a two-dimensional saddle point is discussed. These results may break a path for studying further the fusion in superheavy elements synthesis.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Production cross section of neutron cluster from the reaction of light neutron rich nuclei

Within the Boltzmann-Langevin equation, the neutron cluster production cross sections in the reactions induced by Be-14, He-8, He-6, Li-11, B-17, Be-11, C-19 on C-12 at 35MeV/u were studied. The experimental data for (4)n production cross section from Be-14+C-12 at 35MeV/u can be reproduced. It is found that the production cross section of neutron cluster is large in the reaction that the projectile has more halo nucleons. And the projectiles with big mass number are easy to produce the neutron cluster, when they have the same number of halo nucleons. The neutron cluster is probably mainly from the halo nucleons of projectile.

利用碎裂反应研究6He和6Li的团簇结构

本文在碎裂反应实验的基础上研究了原子核的团簇结构以及晕核的特性。由兰州放射性束流线(RIBLL)提供的6He和6Li次级束与Be靶的反应，用24单元的CsI（Tl）阵列探测器测量了6He和6Li的碎裂产物的能量以及小范围内的角分布。用脉冲形状的方法鉴别出了6He与6Li的轻带电粒子碎片。通过重构不变质量谱的方法得到了RGM预言的6He和6Li的3核子共振态。拟合实验数据发现6Li在Ex=22.9MeV和Ex=30.3MeV出现t+3He团簇共振态，6He在Ex=19.2MeV和Ex=29.8MeV出现t+t团簇共振态。并且与理论以及其它的实验的结果进行了比较。这是在实验上首次通过重构不变质量谱的方法直接观测到6He和6Li的三核子共振态。比较了6He和6Li的碎片产生几率，结果表明丰中子的6He产生丰中子碎片的几率要比6Li大，而产生缺中子碎片的几率要比6Li低，结论与BLE的计算是一致的。通过比较6He和6Li反应产生3He和3H的截面角分布，发现丰中子核反应容易产生丰中子的碎片这种现象在前角区更为敏感，随着角度的增大差别逐渐变小。因此前角区的3H与3He的发射产额比是一个鉴别晕核的灵敏探针。用双碎片关联和电粒子产额估算了6He碎裂各出射道的分支比，还估算了6Li碎裂成双带电粒子出射道的相对比例。由出射道的分支比估算了自由发射的中子质子比，并且与QMD的结果进行了比较，发现晕核的自由发射中子质子比要比稳定核大，这也可作为鉴别晕核的一种探针。另外用同位旋相关的BLE(Boltzmann-Langevin Equation)模型计算了轻丰中子核反应中子团簇集团产生的截面，发现3，4n团簇的截面最大，其中4n的产生截面与实验值相近。还计算了丰质子核反应的轻带电粒子发射的同位旋效应，发现质子晕核反应中轻带电粒子的产额相对于稳定核会突然的增大。这也为寻找和判断晕核结构提供了新的手段

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.

A lattice Boltzmann equation for waves

We propose a lattice Boltzmann model for the wave equation. Using a lattice Boltzmann equation and the Chapman-Enskog expansion, we get 1D and 2D wave equations with truncation error of order two. The numerical tests show the method can be used to simulate the wave motions.