Parkes problems revisited - Application of response number

Response number R-n(n), proposed in [3, 4], is an important independent dimensionless number for the dynamic response of structures [2]. In this paper, the response number is applied to the dynamic plastic response of the well-known Parkes' problem, i.e., beams struck by concentrated mass.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

A finite deformation analysis of polycrystal by considering the influence of grain boundary

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.

A characteristics study on the performance of a 2-stage light gas gun

In order to obtain an overall and systematic understanding of the performance of a two-stage light gas gun (TLGG), a numerical code to simulate the process occurring in a gun shot is advanced based on the quasi-one-dimensional unsteady equations of motion with the real gas effect,;friction and heat transfer taken into account in a characteristic formulation for both driver and propellant gas. Comparisons of projectile velocities and projectile pressures along the barrel with experimental results from JET (Joint European Tons) and with computational data got by the Lagrangian method indicate that this code can provide results with good accuracy over a wide range of gun geometry and loading conditions.

Behavior of ar plasma formed in a high-density plasma source - an ecr reactor

In order to develop the ultra-large scale integration(ULSI), low pressure and high density plasma apparatus are required for etching and deposit of thin films. To understand critical parameters such as the pressure, temperature, electrostatic potential and energy distribution of ions impacting on the wafer, it is necessary to understand how these parameters are influenced by the power input and neutral gas pressure. In the present work, a 2-D hybrid electron fluid-particle ion model has been developed to simulate one of the high density plasma sources-an Electron Cyclotron Resonance (ECR) plasma system with various pressures and power inputs in a non-uniform magnetic field. By means of numerical simulation, the energy distributions of argon ion impacting on the wafer are obtained and the plasma density, electron temperature and plasma electrostatic potential are plotted in 3-D. It is concluded that the plasma density depends mainly on both the power input and neutral gas pressure. However, the plasma potential and electron temperature can hardly be affected by the power input, they seem to be primarily dependent on the neutral gas pressure. The comparison shows that the simulation results are qualitatively in good agreement with the experiment measurements.

Influence of convections on dopant segregation in floating-zone crystal-growth system

The steady and axisymmetric crystal growth process of floating zone model was studied numerically to concern with the influence of convection and phase change on effective segregation. An iteration method of numerical simulation considering both thermocapillary and buoyancy effects for GaAs crystal growth gave the effective segregation coefficient, which was compared with the space experiment of GaAs on board the Chinese recoverable satellite. The calculated segregation coefficient of a two-dimensional model was found to be smaller than the one suggested by space experiment with the simplified assumption of an one-dimensional model.

A study of the mechanism of consolidating metal-powder under explosive-implosive shock-waves

A study of the two-dimensional flow pattern of particles in consolidation process under explosive-implosive shock waves has been performed to further understand the mechanism of shock-wave consolidation of metal powder, in which bunched low-carbon steel wires were used instead of powder. Pressure in the compact ranges from 6 to 30 GPa. Some wires were electroplated with brass, some pickled. By this means, the flow pattern at particle surfaces was observed. The interparticle bonding and microstructure have been investigated systematically for the consolidated specimens by means of optical and electron microscopy, as well as by microhardness. The experimental results presented here are qualitatively consistent with Williamson's numerical simulation result when particle arrangement is close packed, but yield more extensive information. The effect of surface condition of particle on consolidation quality was also studied in order to explore ways of increasing the strength of the compacts. Based on these experiments, a physical model for metal powder shock consolidation has been established.

A ground hydrologic model with inclusion of a layer of vegetation canopy that can interact with general-circulation model

In this paper, a ground hydrologic model(GHM) is presented in which the vapor, heat and momentum exchanges between ground surface covers (including vegetation canopy) and atmosphere is described more realistically. The model is used to simulate three sets of field data and results from the numerical simulation agree with the field data well. GHM has been tested using input data generated by general circulation model (GCM) runs for both the North American regions and the Chinese regions, The results from GHM are quite different from those of GHMs in GCMs. It shows that a more active concerted effort on the land surface process study to provide a physically realistic GHM for predicting the exchange between land and atmosphere is important and necessary.

变截面流化床理论模型及数值模拟

基于双流体模型,在一定假设的条件下推导得到了用于描述变截面床流化过程的局部平衡模型。对于局部平衡模型,空间导数项的离散采用五阶精度的WENO有限差分格式,时间导数项的离散采用TVD Runge-Kutta型的离散格式,对流量突变后的瞬态过程进行了模拟,得到固相体积分数在整个变化过程中沿床高的分布以及床高变化规律和床层表面颗粒速度变化曲线。对于流量突增过程,在床内不同位置形成了一系列的连续波,当所有的连续波到达床面整个变化过程结束。而对于流量突然减小过程,将会有固相体积分数间断在分布板处形成,当所有间断到达床面时,塌落过程结束。

激光微弯曲成形机理的数值研究

利用自行研制的含热传导、冲击动力学大、变形有限元程序，模拟了小尺寸梁在脉冲激光加热条件下的变形过程。在此基础上，利用商用程序模拟了冷却及残余应力的产生，研究了激光参数（强度及分布）等对于微弯曲的影响。数值模拟结果与文献中的实验观察相吻合。

气相爆轰波数值模拟中化学反应模型研究

基于改进的时-空守恒元解元算法对气相爆轰波数值模拟中3种常用化学反应模型(二步模型,基元反应模型和Sichel的二步模型)进行了考察。对平面爆轰波和具有胞格结构的爆轰波进行了数值模拟,并对数值结果进行了比较和讨论。结果表明:3种化学反应模型得到的爆轰参数准确性有所差异,但得到的胞格结构均能和实验结果较好吻合。3种化学反应模型在爆轰波数值模拟中各有优缺点,应视具体问题决定使用哪种化学反应模型。

微重力蜡烛火焰特征数值模拟

建立了微重力蜡烛火焰的数学模型。计算与分析表明，火焰的形状由空气动力学特征决定，火焰是温度取决于化学反应动力学特征和火焰的热损失。在静止微重力环境中，自然对流的消失使火焰为半球形。辐射热损失对烛火焰温度（颜色）特征的形成有重要贡献，在静止重力环境下，化学反应放热速率受氧气扩散速率控制，辐射热损失的冷却使火焰温度低于正常重力温度值。但当环境气体的流动速度加大时，辐射热损失的影响逐渐减小，蜡烛火焰的温度逐渐接近正常重力蜡烛火焰的温度。当氧浓度较小时，火焰峰值温度小于烟黑形成的阈值温度（1300 K）；当氧浓度较大时，火焰温度大于烟黑形成的阈值温度。

应力波在多层介质中传播特性数值分析

通过与相关文献结果的对比,在验证了数值模拟应力波传播可行性的基础上,比较了应力波通过3层花岗岩和夹层为泡沫铝的3层介质后,发现后者应力波幅值的衰减远大于前者,应变能增加为前者的1.6倍,证明了软夹层在研究的速度量级上对能量耗散具有显著作用;在入射波波长为3层介质总厚度1/2的条件下,当泡沫铝的厚度占总厚度约0.2时,得出了3层介质的应变能约为系统总能量的60%,此时入射波波长为泡沫铝厚度的2.5倍,组合介质获得较佳的衰减性能.

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率，对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序。该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想，验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合，说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动，为了提高计算规模及计算速度，作者对程序进行了并行化并针对微机机群系统进行了并行优化，实际计算表明本程序具有较高的并行效率。