Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier--Stokes method

In this study, the Reynolds-averaged Navier--Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of $k$--$\varepsilon$ turbulence closure schemes, namely, standard model, renormalization group (RNG) model and realizable $k$--$\varepsilon$model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizable $k$--$\varepsilon$ model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable $k$--$\varepsilon$ model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.     

激光二极管端面抽运的高效率1.34μm Nd:YVO4 平凹腔型激光器

理论方面分析了端面抽运的1.34mm Nd:YVO4平凹腔型激光器的各项参数；文中的分析可为同类激光器的设计提供参考。实验方面实现了LD端面抽运的1.34mm Nd:YVO4平凹腔型激光器的运转；在输入泵浦功率为9.0W时，获得了1.34mm 激光输出功率3.51W，光－光转换效率为39.0%。     

非定常粘性空化流动模型及其数值计算

基于气液两相当地均相介质模型,本文给出了一种模拟非定常粘性空化流动的计算模型。认为空化绕流流场中流动介质是一种当地均匀的气液混相物。控制方程采用了应用两相流模型的N-S方程。基于液相和气相的状态方程推导了混合介质密度的表达式。为了保证数值计算的稳定性,控制方程的数值求解采用了TVD-MacCormack格式。为了评价计算模型的可靠性,分别计算了绕台阶和管道水锤的空化流动,所得结果是合理的,说明该方法可以用于空化流动的数值计算。     

Hopf Bifurcations,Drops in the Lid-Driven Square Cavity Flow

The lid-driven square cavity flow is investigated by numerical experiments. It is found that from $ \mathrm{Re} $$=$ $5,000 $ to $ \mathrm{Re} $$=$$ 7,307.75 $ the solution is stationary, but at $ \mathrm{Re}$$=$$7,308 $ the solution is time periodic. So the critical Reynolds number for the first Hopf bifurcation localizes between $ \mathrm{Re} $$=$$ 7,307.75 $ and $ \mathrm{Re} $$=$$ 7,308 $. Time periodical behavior begins smoothly, imperceptibly at the bottom left corner at a tiny tertiary vortex; all other vortices stay still, and then it spreads to the three relevant corners of the square cavity so that all small vortices at all levels move periodically. The primary vortex stays still. At $ \mathrm{Re} $$=$$ 13,393.5 $ the solution is time periodic; the long-term integration carried out past $ t_{\infty} $$=$$ 126,562.5 $ and the fluctuations of the kinetic energy look periodic except slight defects. However, at $ \mathrm{Re} $$=$$ 13,393.75 $ the solution is not time periodic anymore: losing unambiguously, abruptly time periodicity, it becomes chaotic. So the critical Reynolds number for the second Hopf bifurcation localizes between $ \mathrm{Re} $$=$$ 13,393.5 $ and $ \mathrm{Re} $$=$$ 13,393.75 $. At high Reynolds numbers $ \mathrm{Re} $$=$$ 20,000 $ until $ \mathrm{Re} $$=$$ 30,000 $ the solution becomes chaotic. The long-term integration is carried out past the long time $ t_{\infty} $$=$$ 150,000 $, expecting the time asymptotic regime of the flow has been reached. The distinctive feature of the flow is then the appearance of drops: tiny portions of fluid produced by splitting of a secondary vortex, becoming loose and then fading away or being absorbed by another secondary vortex promptly. At $ \mathrm{Re} $$=$$ 30,000 $ another phenomenon arises—the abrupt appearance at the bottom left corner of a tiny secondary vortex, not produced by splitting of a secondary vortex.     

Characteristics of temperature fluctuation in two-dimensional turbulent Rayleigh-Bénard convection

We study the characteristics of temperature fluctuation in two-dimensional turbulent Rayleigh–Benard convection in′a square cavity by direct numerical simulations.The Rayleigh number range is 1×10⁸≤Ra≤1×10¹³,and the Prandtl number is selected as Pr=0.7 and Pr=4.3.It is found that the temperature fluctuation profiles with respect to Ra exhibit two different distribution patterns.In the thermal boundary layer,the normalized fluctuationθrms/θrms,max is independent of Ra and a power law relation is identified,i.e.,θrms/θrms,max～(z/δ)0.99±0.01,where z/δis a dimensionless distance to the boundary(δis the thickness of thermal boundary layer).Out of the boundary layer,when Ra≤5×10⁹,the profiles ofθrms/θrms,max descend,then ascend,and finally drop dramatically as z/δincreases.While for Ra≥1×10¹⁰,the profiles continuously decrease and finally overlap with each other.The two different characteristics of temperature fluctuations are closely related to the formation of stable large-scale circulations and corner rolls.Besides,there is a critical value of Ra indicating the transition,beyond which the fluctuation hθrmsiV has a power law dependence on Ra,given by hθrmsiV～Ra?0.14±0.01.     

Judd-Oflet analysis of spectrum and laser performance of Ho：YAP crystal end-pumped by 1.91μm Tm：YLF laser

The Ho:YAP crystal is grown by the Czochralski technique.The room temperature polarized absorption spectra of Ho:YAP crystal was measured on a c cut sample with 1 at% holmium.According to the obtained Judd-Ofelt intensity parameters Ω2 = 1.42 × 10-20 cm2,Ω4 = 2.92 × 10-20 cm2,and Ω6 = 1.71 × 10-20 cm2,this paper calculated the fluorescence lifetime to be 6 ms for 5I7 →5 I8 transition,and the integrated emission cross section to be 2.24×10-18 cm2.It investigates the room temperature Ho:YAP laser end pumped by a 1.91 μm Tm:YLF laser.The maximum output power was 4.1 W when the incident 1.91 μm pump power was 14.4 W.The slope efficiency is 40.8%,corresponding to an optical to optical conversion efficiency of 28.4%.The Ho:YAP output wavelength was centred at 2118 nm with full width at half maximum of about 0.8 nm.     

在HT-7托卡马克上密度扫描实验中的边界湍流和输运

在HT-7托卡马克上密度扫描实验中使用一个快速扫描气动朗缪尔探针来研究边界湍流和输运。随着中心弦平均等离子体密度的增加，在等离子体边界观察到径向电场剪切增强，径向电场的增强能够解释边界输运的减小和全局粒子约束的改善。粒子约束时间的增加和垂直扩散系数的减小证明粒子输运降低了。在等离子体边界区域发现电离和辐射驱动湍流的特征，表明电离和辐射在湍流驱动过程中具有重要性。     

Co-Al-W基高温合金的团簇成分式

Co-Al-W基高温合金具有类似于Ni基高温合金的γ+γ'相组织结构.根据面心立方固溶体的团簇加连接原子结构模型,Ni基高温合金的成分式即最稳定的化学近程序结构单元可以描述为第一近邻配位多面体团簇加上次近邻的三个连接原子.本文应用类似方法,首次给出了Co-Al-W基高温合金的团簇成分式.利用原子半径和团簇共振模型,可计算出Co-Al-W三元合金的团簇成分通式,为[Al-Co_(12)](Co,Al,W)_3,即以Al为中心原子、Co为壳层原子的[Al-Co_(12)]团簇加上三个连接原子.对于多元合金,需要先将元素进行分类:溶剂元素——类Co元素Co (Co, Cr, Fe, Re, Ni,Ir,Ru)和溶质元素——类Al元素Al (Al,W,Mo, Ta,Ti,Nb,V等);进而根据合金元素的配分行为,将类Co元素分为Co~γ(Cr, Fe, Re)和Co~(γ')(Ni, Ir, Ru);根据混合焓,将类Al元素分为Al, W (W, Mo)和Ta (Ta, Ti, Nb, V等).由此,任何多元Co-Al-W基高温合金均可简化为Co-Al伪二元体系或者Co-Al-(W,Ta)伪三元体系,其团簇加连接原子成分式为[Al-Co_(12)](Co_(1.0)Al_(2.0))(或[Al-Co_(12)] Co_(1.0)Al_(0.5)(W,Ta)_(1.5)=Co_(81.250)Al_(9.375)(W,Ta)_(9.375) at.%).其中,γ与γ'相的团簇成分式分别为[Al-Co_(12)](Co_(1.5)Al_(1.5))(或[Al-Co_(12)] Co_(1.5)Al_(0.5)(W,Ta)_(1.0)=Co_(84.375)Al_(9.375)(W,Ta)_(6.250) at.%)和[Al-Co_(12)](Co_(0.5)Al_(2.5))(或[Al-Co_(12)] Co_(0.5)Al_(0.5)(W, Ta)_(2.0)=Co_(78.125)Al_(9.375)(W,Ta)_(12.500)at.%).例如,Co_(82)Al_9W_9合金的团簇成分式为[Al-Co_(12)]Co_(1.1)Al_(0.4)W_(1.4)(～[Al-Co_(12)]Co_(1.0)Al_(0.5)W_(1.5)),其中γ相的团簇成分式为[Al-Co_(12)]Co_(1.6)Al_(0.4)W_(1.0)(～[Al-Co_(12)]Co_(1.5)Al_(0.5)W_(1.0)),γ'相的团簇成分式为[Al-Co_(12)]Co_(0.3)Al_(0.5)W_(2.2)(～[AlCo_(12)]Co_(0.5)Al_(0.5)W_(2.0)).     

Waveguide laser modelling of Erbium/Ytterbium activated monoclinic double tungstates

Using the overlapping integral method, a rib waveguide laser of monoclinic potassium double tungstate, KRE $(\text{ WO}_{4})_{2}$ , co-doped with Erbium and Ytterbium has been modelled. The laser operation at 1.5 $\upmu $ m is based on an efficient pump scheme via the energy transfer from Yb $^{3+}$ to Er $^{3+}$ ions. The numerical simulation requires spectroscopic parameters of the ions involved and the waveguide geometry and index profiles. This model allows determining the laser power as a function of controllable parameters such as ions doping level, pump power, cavity length or reflectance of the input/output mirrors. It has been found that, for the standard doping level used in this matrix, the optimum cavity length is only few millimeters. Overall, using simulation tools is possible to optimize fabrication parameters, and thus saving effort in the development of experimental prototypes.     

气体中心式离心喷嘴喷雾实验与三维仿真

在实验的基础上, 基于RNG k-ε模型对常压下气体中心式同轴离心(gas-centered swirl coaxial,GCSC)喷嘴喷雾形态和破碎模式进行了三维仿真研究。采用网格自适应加密(adaptive mesh refinement,AMR)技术、耦合水平集和流体体积(coupled level-set and volume of fluid, CLSVOF)方法对气液界面进行捕捉。结果表明, 液体质量流率($\dot{m}_{\mathrm{l}}$)不变, 随着气体质量流率($\dot{m}_{\mathrm{g}}$)的增加, 中心气流的引射作用增强, 液膜内外压差增大, 雾化锥角减小, 并对其流动特性进行了分析; 而$\dot{m}_{\mathrm{g}}$不变时, 液膜在喷嘴出口的径向速度与切向速度随$\dot{m}_{\mathrm{l}}$的增大而增大, 导致雾化锥角增大。同时根据气液质量流率比(gas-liquid mass flow rate,GLR), 将喷雾的破碎模式分为穿孔破碎、气泡破碎和气动破碎。      

Steady-state and transient electronic transport properties of β-(AlxGa1-x)2O3/Ga2O3 heterostructures: An ensemble Monte Carlo simulation

The steady-state and transient electron transport properties of $\beta $-(Al$_{x}$Ga$_{1-x}$)$_{2}$O$_{3}$/Ga$_{2}$O$_{3}$ heterostructures were investigated by Monte Carlo simulation with the classic three-valley model. In particular, the electronic band structures were acquired by first-principles calculations, which could provide precise parameters for calculating the transport properties of the two-dimensional electron gas (2DEG), and the quantization effect was considered in the $\varGamma $ valley with the five lowest subbands. Wave functions and energy eigenvalues were obtained by iteration of the Schrödinger-Poisson equations to calculate the 2DEG scattering rates with five main scattering mechanisms considered. The simulated low-field electron mobilities agree well with the experimental results, thus confirming the effectiveness of our models. The results show that the room temperature electron mobility of the $\beta $-(Al$_{0.188}$Ga$_{0.812}$)$_{2}$O$_{3}$/Ga$_{2}$O$_{3}$ heterostructure at 10 kV$ \cdot$cm$^{-1}$ is approximately 153.669 cm$^{2}\cdot$V$^{-1}\cdot$s$^{-1}$, and polar optical phonon scattering would have a significant impact on the mobility properties at this time. The region of negative differential mobility, overshoot of the transient electron velocity and negative diffusion coefficients are also observed when the electric field increases to the corresponding threshold value or even exceeds it. This work offers significant parameters for the $\beta$-(Al$_{x}$Ga$_{1-x}$)$_{2}$O$_{3}$/Ga$_{2}$O$_{3}$ heterostructure that may benefit the design of high-performance $\beta$-(Al$_{x}$Ga$_{1-x}$)$_{2}$O$_{3}$/Ga$_{2}$O$_{3}$ heterostructure-based devices.     

Nanocrystallization behaviour of a ternary amorphous alloy during isothermal annealing： a Monte Carlo simulation

The nanocrystallization behaviour of Zr70Cu20Ni10 metallic glass during isothermal annealing is studied by employing a Monte Carlo simulation incorporating with a modified Ising model and a Q-state Potts model. Based on the simulated microstructure and differential scanning calorimetry curves, we find that the low crystal-amorphous interface energy of Ni plays an important role in the nanocrystallization of primary Zr2Ni. It is found that when T〈T1max （where T1max is the temperature with maximum nucleation rate）, the increase of temperature results in a larger growth rate and a much finer mierostrueture for the primary Zr2Ni, which accords with the microstructure evolution in ＂flash annealing＂. Finally, the Zr2Ni/Zr2Cu interface energy σG contributes to the pinning effect of the primary nano-sized Zr2Ni grains in the later formed normal Zr2Cu grains.     

Solving forward and inverse problems of the nonlinear Schrödinger equation with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential via PINN deep learning

In this paper, based on physics-informed neural networks (PINNs), a good deep learning neural network framework that can be used to effectively solve the nonlinear evolution partial differential equations (PDEs) and other types of nonlinear physical models, we study the nonlinear Schrödinger equation (NLSE) with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential, which is an important physical model in many fields of nonlinear physics. Firstly, we choose three different initial values and the same Dirichlet boundary conditions to solve the NLSE with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential via the PINN deep learning method, and the obtained results are compared with those derived by the traditional numerical methods. Then, we investigate the effects of two factors (optimization steps and activation functions) on the performance of the PINN deep learning method in the NLSE with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential. Ultimately, the data-driven coefficient discovery of the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential or the dispersion and nonlinear items of the NLSE with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential can be approximately ascertained by using the PINN deep learning method. Our results may be meaningful for further investigation of the nonlinear Schrödinger equation with the generalized ${ \mathcal P }{ \mathcal T }$-symmetric Scarf-II potential in the deep learning.     

Mutual recombination in slow Si^＋ ＋ H^- collisions

This paper studies the process of mutual neutralization of Si^＋ and H^- ions in slow collisions within the multichannel Landau-Zener model. All important ionic-covalent couplings in this collision system are included in the collision dynamics. The cross sections for population of specific final states of product Si atom are calculated in the CM energy range 0.05 e∨/u-5 ke∨/u. Both singlet and triplet states are considered. At collision energies below -10 e∨/u, the most populated singlet state is Si（3p4p, ^1S0）, while for energies above -150e∨/u it is the Si（3p, 4p, ^1P1） state. In the case of triplet states, the mixed 3p4p（^3S1 ＋^3P0） states are the most populated in the entire collision energy range investigated. The total cross section exhibits a broad maximum around 200 300e∨/u and for ECM ≤ 10e∨/u it monotonically increases with decreasing the collision energy, reaching a value of 8 × 10^-13 cm^2 at ECM = 0.05 e∨/u. The ion-pair formation process in Si（3p^2 ^3PJ）＋H（1s） collisions has also been considered and its cross section in the considered energy range is very small （smaller than 10^-20 cm^2 in the energy region below 1 ke∨/u）.     

The fabrication of nickel silicide ohmic contacts to n-type 6H-silicon carbide

This paper reports that the nickel silicide ohmic contacts to n-type 6H-SiC have been fabricated. Transfer length method test patterns with NiSi/SiC and NiSi硅化镍;欧姆触点;n型碳化硅;制造;能带;带隙Project supported by the National Basic Research Program of China (Grant No~2002CB311904), the National Defense Basic Research Program of China (Grant No~51327010101) and the National Natural Science Foundation of China (Grant No~60376001).2006-09-192006-10-30This paper reports that the nickel silicide ohmic contacts to n-type 6H-SiC have been fabricated. Transfer length method test patterns with NiSi/SiC and NiSi2/SiC structure axe formed on N-wells created by N^＋ ion implantation into Si-faced p-type 6H-SiC epilayer respectively. NiSi and NiSi2 films are prepared by annealing the Ni and Si films separately deposited. A two-step annealing technology is performed for decreasing of oxidation problems occurred during high temperature processes. The specific contact resistance Pc of NiSi contact to n-type 6H-SiC as low as 1.78× 10^-6Ωcm^2 is achieved after a two-step annealing at 350 ℃for 20 min and 950℃ for 3 min in N2. And 3.84×10-6Ωcm^2 for NiSi2 contact is achieved. The result for sheet resistance Rsh of the N＋ implanted layers is about 1210Ω/□. X-ray diffraction analysis shows the formation of nickel silicide phases at the metal/n-SiC interface after thermal annealing. The surfaces of the nickel silicide after thermal annealing are analysed by scanning electron microscope.     

Charge orderings in the atomic limit of the extended Hubbard model

The extended Hubbard model in the atomic limit (AL-EHM) on a square lattice with periodic boundary conditions is studied with use of the Monte Carlo (MC) method. Within the grand canonical ensemble the phase and order-order boundaries for charge orderings are obtained. The phase diagrams include three types of charge ordered phases and the nonordered phase. The system exhibits very rich structure and shows unusual multicritical behavior. In the limiting case of t_ij=0, the EHM is equivalent to the pseudospin model with single-ion anisotropy , exchange interaction W in an effective magnetic field . This classical spin model is analyzed using the MC method for the canonical ensemble. The phase diagram is compared with the known results for the Blume-Capel model.     

Effective action for phase fluctuations in <Emphasis Type="Italic">d</Emphasis>-wave superconductors near a Mott transition

Phase fluctuations of a d-wave superconducting order parameter are theoretically studied in the context of high-T_c cuprates. We consider an extended t-J model describing electrons in a layer which also contains long-range Coulomb interactions. The constraint of having at most singly occupied sites is enforced by an additional Hubbard term. The Heisenberg interaction is decoupled by a d-wave order parameter in the particle-particle channel. Assuming first that the equilibrium state has long-range phase order, the effective action is derived perturbatively for small fluctuations within a path integral formalism, in the presence of the Coulomb and Hubbard interaction terms. In a second step, a more general derivation of is performed in terms of a gradient expansion which only assumes that the gradients of the order parameter are small whereas the value of the phase may be large. We show that in the phase-only approximation the resulting reduces in leading order in the field gradients to the perturbative one which thus allows to treat also the case without long-range phase order or vortices. Our result generalizes previous expressions for to the case of interacting electrons, is explicitly gauge invariant, and avoids problematic singular gauge transformations.     

Guiding of 150 keV O6＋ ions through nanocapillaries in an uncoated Al2O3 membrane： special time dependence of the transmission profile width

<正>This paper reports that the transmission of O⁶⁺ ions with energy of 150keV through capillaries in an uncoated Al₂O₃ membrane was measured,and agreements with previously reported results in general angular distribution of the transmitted ions and the transmission fractions as a function of the tilt angle well fitted to Gaussian-like functions were observed.Due to using an uncoated capillary membrane,ourψ_c is larger than that using a gold-coated one with a smaller value of（?）,which suggests a larger equilibrium charge Q_∞in our experiment.The observed special width variation with time and a larger width than that using a smaller（?） were qualitatively explained by using mean-field classical transport theory based on a classical-trajectory Monte Carlo simulation.     

Sideband manipulation of population inversion in a three-level Λ atomic system

Sideband manipulation of population inversion in a three-level A atomic configuration is investigated theoretically. Compared with the case of a nearly monochromatic field, a population inversion between an excited state and a ground state has been found in a wide sideband intensity range by increasing the difference in frequency between three components. Furthermore, the population inversion can be controlled by the sum of the relative phases of the sideband components of the trichromatic pump field with respective to the phase of the central component. Changing the sum phase from 0 to π, the population inversion between the excited state and the ground state can increase within nearly half of the sideband intensity range. At the same time, the sideband intensity range that corresponds to the system exhibiting inversion ρ00 〉 ρ11 also becomes wider evidently.     

Kinetic Ising model in a time-dependent oscillating external magnetic field: effective-field theory

Recently,Shi et al.[2008 Phys.Lett.A 372 5922] have studied the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field and presented the dynamic phase diagrams by using an effective-field theory(EFT) and a mean-field theory(MFT).The MFT results are in conflict with those of the earlier work of Tom’e and de Oliveira,[1990 Phys.Rev.A 41 4251].We calculate the dynamic phase diagrams and find that our results are similar to those of the earlier work of Tom’e and de Oliveira;hence the dynamic phase diagrams calculated by Shi et al.are incomplete within both theories,except the low values of frequencies for the MFT calculation.We also investigate the influence of external field frequency(ω) and static external field amplitude(h0) for both MFT and EFT calculations.We find that the behaviour of the system strongly depends on the values of ω and h0.