~(133)Sm和~(149)Yb的β缓发质子衰变

利用40Ca+96Ru融合蒸发反应产生了近质子滴线核133Sm,配合氦喷嘴带传输系统采用“质子-γ”符合方法观测了它们的β缓发质子衰变,其中包括半衰期、质子能谱、第二代子核低位态之间的γ跃迁,并估计出衰变到第二代子核不同低位态的分支比．通过统计理论拟合上述实验数据,指认了133Sm的自旋宇称的可能范围．并用Woods-Saxon Strutinsky方法计算了限制组态的133Sm的核势能面,通过对比发现133Sm的自旋宇称可能有两种成分:5／2+和1／2-．这一结果与2001年发表的133Sm(EC+β+)衰变的简单衰变纲图是相容的．此外用同一方法分析了2001年Eur．Phys．J．A12: 1-4中发表的有关149Yb的β缓发质子衰变实验数据,由此指认了149Yb的基态自旋宇称为1／2-．

beta-delayed proton decays of Sm-133 and Yb-149

Sm-133 was produced via fusion evaporation in the reaction Ca-40+Ru-96. Its P-delayed proton decay was studied by means of "p-gamma" coincidence in combination with a He-jet tape transport system, including half-lives, proton energy spectra, gamma-transitions following the proton emissions, and the branching ratios to the low-lying states in the grand-daughter nuclei. The possible spins and parities of 133Sm were extracted by fitting the experimental data with a statistical model calculation. The configuration-constrained nuclear potential energy surfaces of Sm-133 were calculated by using the Woods-Saxon Strutinsky method. Comparing the experimental and calculated results, the spins and parities Of Sm-133 were assigned to be 5/2(+) and 1/2(-), which is reconciled with our published simple (EC+beta(+)) decay scheme Of Sm-113 in 2001. In addition, our experimental data on the beta-delayed proton decay of Yb-149 reported in Eur. Phys. J., 2001, A12: 1-4 was also analyzed by using the same method. The spin and parity of Yb-149 was assigned to be 1/2-.

beta-delayed proton decays and spin assignments for Sm-133 and Yb-149

The proton-rich isotope Sm-133 was produced via the fusion evaporation reaction Ca-40 + Ru-96. Its beta-delayed proton decay was studied by p-gamma coincidence in combination with a He-jet tape transport system, and half-lives, proton energy spectra, gamma-transitions following the proton emission, as well as beta-delayed proton branching ratios to the low-lying states in the grand-daughter nucleus were determined. Comparing the observed beta-delayed proton branching ratios with statistical model calculations, the best agreement is found assuming that only one level with the spin of 3/2 in Sm-133 decays or two levels with the spins of 1/2 and 5/2 decay with similar half-lives. The configuration-constrained nuclear potential energy surfaces of Sm-133 were calculated using the Woods-Saxon-Strutinsky method, which suggests a 1/2-ground state and a 5/2(+) isomer with an excitation energy of 120 keV. Therefore, the simple(EC+beta(+)) decay scheme of Sm-133 in Eur. Phys. J.A 11,277(2001) has been revised. In addition, our previous experimental data on the beta-delayed proton decay of Yb-149 reported in Eur. Phys. J. A 12,1 ( 2 0 0 1) was also analyzed using the same method. The spin-parity of Yb-149 is suggested to be 1/2(-).

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

煤与瓦斯突出的层裂－粉碎模型

提出了描述煤与瓦斯突出的１维模型。该模型把煤体破坏分为２个阶段──层裂和层裂片的粉碎，并用粉碎率的概念来描述煤体的这种非均匀破坏。利用上述模型和从模拟实验得来的参数，对突出过程进行数值计算，得到了破裂阵面推进速度与瓦斯初始压力和煤初始粉碎率之间的拟合关系．计算结果与实验结果符合得很好。在煤激波管的实验中，揭示了一系列重要现象。值得提到的是，在严格控制实验煤样参数不变，仅改变充气气体压力的情况下，发现有３种典型现象：当充气气压ｐｇ小于某一临界值ｐｃｌ时，煤样不发生破坏；当ｐｇ＞ｐｃｌ但小于另一临界值ｐｃ２时，可观察到部分煤样破裂；而当ｐｇ＞ｐｃ２时，观察到煤样全部破裂并被抛出，破裂阵面近似恒稳地推进。在建立本文的物理模型时，将利用这些实验事实。１物理模型比照郑哲敏等的１维实验，研究图１所示装置破膜后的过程。煤与瓦斯突出研究的关键问题是对煤体破裂的研究，即对煤的破裂准则和破裂阵面运动规律的研究。对此本文提出以下模型：（１）把最大有效拉应力准则作为破裂准则，认为破裂是不连续的，破裂阵面跳跃式发展；（２）引入混相流和粉碎率的概念来描述层裂片的运动；（３）层裂片的形态在运动中不稳定，满足一定条件后将发生进一步的粉碎。

实时TV光电子全息散斑干涉仪及其在折射率场中的应用

介绍了一种不需PC微机支持图象卡的实时TV光电子全息散斑干涉仪,成功地应用于折射率流场的显示定量测量,并给出了实验结果。

Study of the shock motion in a hypersonic shock system turbulent boundary-layer interaction

Wall pressure fluctuations and surface heat transfer signals have been measured in the hypersonic turbulent boundary layer over a number of compression-corner models. The distributions of the separation shock oscillation frequencies and periods have been calculated using a conditional sampling algorithm. In all cases the oscillation frequency distributions are of broad band, but the most probable frequencies are low. The VITA method is used for deducing large scale disturbances at the wall in the incoming boundary layer and the separated flow region. The results at present showed the existence of coherent structures in the two regions. The zero-cross frequencies of the large scale structures in the two regions are of the same order as that of the separation shock oscillation. The average amplitude of the large scale structures in the separated region is much higher than that in the incoming boundary layer. The length scale of the separation shock motion region is found to increase with the disturbance strength. The results show that the shock oscillation is of inherent nature in the shock wave/turbulent boundary layer interaction with separation. The shock oscillation is considered to be the consequence of the coherent structures in the separated region.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Crystal-growth in floating zone with phase-change and thermosolutal convections

Floating zone crystal growth in microgravity environment is investigated numerically by a finite element method for semiconductor growth processing, which involves thermocapillary convection, phase change convection, thermal diffusion and solutal diffusion. The configurations of phase change interfaces and distributions of velocity, temperature and concentration fields are analyzed for typical conditions of pulling rates and segregation coefficients. The influence of phase change convection on the distribution of concentration is studied in detail. The results show that the thermocapillary convection plays an important role in mixing up the melt with dopant. The deformations of phase change interfaces by thermal convection-diffusion and pulling rods make larger variation of concentration field in comparison with the case of plane interfaces.

激波风洞中超声速燃烧现象的初步实验研究

介绍了启用激波风洞进行超燃研究所作的改造和得到的初步实验结果。实验除测量常规流场参数和模型表面的动态参数外，主要使用二维实验模型，利用高温气体的自发光和燃烧伴随的发光现象，采用高速摄影技术来观测燃烧现象，记录到高焓流场中叠加燃料喷射和燃烧的流场，观察到自由边界条件下扩散控制的超声速燃烧现象。结果表明激波风洞用于超声速燃烧研究的一些必要的技术问题已基本解决。

试论农业工程中的力学问题及主要研究方向

综述了农业工程中所存在的主要力学问题或现象,如:应力与植物的生长、种子损伤与激活、土壤粒子(泥沙)在雨水和江河中的流动、农作物体内水分、矿质元素和光合产物的运输、农作物种植、农产品运输和土壤耕作过程中所存在的力学问题或现象,并在概述了相关研究现状的同时,着重探讨了今后农业生物力学研究的主要方向,强调了加强生物力学研究的重要性.

纤维束张紧力缠绕复合材料飞轮初应力的三维数值分析

利用平面应力型全弹性模型的思想(即将纤维束张紧力缠绕看成多层复合材料薄环连续过盈装配的过程),建立了三维纤维束张紧力缠绕复合材料飞轮初应力分析模型,并给出了基于面-面接触算法求解张紧力缠绕复合材料飞轮初应力的三维数值方法.算例分析表明,三维数值分析得到的飞轮的环向初应力及径向初始压应力(数值)均略低于平面应力模型的结果,且这种差距随着飞轮轴向长度的增加而缓慢增大;三维分析证实了平面应力模型关于张紧力缠绕复合材料飞轮的初应力分析有足够的精度.最后给出了三维模型轴向效应的表征方法.

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.