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.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

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.

IP方法模拟硬盘读写头气膜润滑问题

本论文中采用IP方法模拟硬盘读写头的气膜润滑问题，并选择耗时但精准的DSMC方法作为该问题中IP方法的检验标准，IP方法计算得到的读写头表面压力分布及读写头所受净升力为读写头的设计提供可靠的参考数据。 首先，将读写头复杂表面简化为平面并假定尺寸为微米（远小于真实的毫米），微米尺寸时可以得到DSMC的模拟结果用于比较。IP法对该问题的模拟结果与DSMC的模拟结果完全吻合。 其次，选定一个真实的读写头表面几何，但仍然假定尺寸为微米。IP法在该问题的应用过程中遇到四点困难，我们提出或采取了相应的解决办法。在复杂表面几何的读写头气膜润滑问题中，IP法的模拟结果与DSMC的模拟结果也相符一致。 最后完成真实几何、真实尺寸（毫米）的读写头气膜润滑问题的模拟，这也是本论文的最终研究目标。IP法模拟该问题时需要采用大网格，为此我们构造二维检验模型验证了大网格的合理性。 真实几何读写头气膜润滑问题的模拟结果表明，在流场的局部区域压力沿着高度方向变化显著，这意味着读写头问题的传统计算方法（即概括性雷诺方程）在该局部区域不成立，故用它计算得到的结果值得认真检验，但文献中尚未给出过这类检验。这正是该论文的研究背景和采用新方法（IP方法）的原因。 本论文研究的问题属于过渡领域、低速流动问题，故本论文还介绍了该领域另外一种粒子模拟方法（LVDSMC方法）以及我们对该方法提出的一些改进。 关键词：硬盘读写头，气膜润滑，概括性雷诺方程，IP方法，DSMC方法，过渡领域，LVDSMC方法

Operation of the qubit based on photon-assisted tunneling in a coupled quantum-dot system and the influence of dephasing

We derive the generalized rate equation for the coupled quantum-dot (QD) system irradiated by a microwave field in the presence of a quantum point contact. It is shown that when a microwave field is tuned in resonance with the energy difference between the ground states of two QD's, the photon-assisted tunneling occurs and, as a result, the coupled QD system may be used as the single qubit. Furthermore, we show that the oscillating current through the detector decays drastically as the dephasing rate increases, indicating clearly the influence of the dephasing effect induced by the quantum point contact used as a detecting device.

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.

Cosmological equations and thermodynamics on apparent horizon in thick braneworld

We derive the generalized Friedmann equation governing the cosmological evolution inside the thick brane model in the presence of two curvature correction terms: a four-dimensional scalar curvature from induced gravity on the brane, and a five-dimensional Gauss-Bonnet curvature term. We find two effective four-dimensional reductions of the generalized Friedmann equation in some limits and demonstrate that the reductions but not the generalized Friedmann equation can be rewritten as the first law of equilibrium thermodynamics on the apparent horizon of thick braneworld.

利用碎裂反应研究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的产生截面与实验值相近。还计算了丰质子核反应的轻带电粒子发射的同位旋效应，发现质子晕核反应中轻带电粒子的产额相对于稳定核会突然的增大。这也为寻找和判断晕核结构提供了新的手段

Thermokinetic investigation of effect of temperature on optimum NaCl concentration for petroleum bacterial growth

The power-time curves of growth of three strains of petroleum bacteria at different NaCl concentrations at 40.0 and 50.0 degreesC have been determined by using a 2277 Thermometric Thermal Activity Analyser. An equation of a power-time curve, ln[alphaP(K)/P(t) - 1] = ln[(alphaK - N-0)/N-0] - alphakt, was established based on the generalized logistic equation, where P(t) is the thermal power at time t, K the carrying capacity, P-K = P0K, P-0 the thermal power of one cell, N-0 the bacterial population at time zero, alpha = (k - D)/k. The method of four observed points with the same time interval was used to calculate the value of P-K. The growth rate constant k and the death rate constant D were calculated. The NaCl concentration of optimum growth rate of petroleum bacteria at 40.0 and 50.0 degreesC, respectively, have been obtained according to the curves k - D versus NaCl concentration, which are 0.26, 0.54 and 0.57 mol l(-1) for B-1, B-2 and B-3, respectively, at 50.0 degreesC, 0.26, 0.55 and 0.56 mol l(-1) for B-1, B-2 and B-3, respectively, at 40.0 degreesC. The results indicated that the effect of temperature on NaCl concentration of optimum growth rate was small. (C) 2002 Elsevier Science B.V. All rights reserved.

A generalized wave action conservative equation for the dissipative dynamical system in the nearshore region

To describe the various complex mechanisms of the dissipative dynamical system between waves, currents, and bottoms in the nearshore region that induce typically the wave motion on large-scale variation of ambient currents, a generalized wave action equation for the dissipative dynamical system in the nearshore region is developed by using the mean-flow equations based on the Navier-Stokes equations of viscous fluid, thus raising two new concepts: the vertical velocity wave action and the dissipative wave action, extending the classical concept, wave action, from the ideal averaged flow conservative system into the real averaged flow dissipative system (that is, the generalized conservative system). It will have more applications.

Work and energy equations and the principle of generalized effective stress for unsaturated soils

The effective stress principle has been efficiently applied to saturated soils in the soil mechanics and geotechnical engineering practice; however, its applicability to unsaturated soils is still under debate. The appropriate selection of stress state variables is essential for the construction of constitutive models for unsaturated soils. Owing to the complexity of unsaturated soils, it is difficult to determine the deformation and strength behaviors of unsaturated soils uniquely with the previous single-effective-stress variable theory and two-effective-stress-variable theory in all the situations. In this paper, based on the porous media theory, the specific expression of work is proposed, and the effective stress of unsaturated soils conjugated with the displacement of the soil skeleton is further derived. In the derived work and energy balance equations, the energy dissipation in unsaturated soils is taken into account. According to the derived work and energy balance equations, all of the three generalized stresses and the conjugated strains have effects on the deformation of unsaturated soils. For considering these effects, a principle of generalized effective stress to describe the behaviors of unsaturated soils is proposed. The proposed principle of generalized effective stress may reduce to the previous effective stress theory of single-stress variable or the two-stress variables under certain conditions. This principle provides a helpful reference for the development of constitutive models for unsaturated soils.