Transition from steady to oscillatory convection with chaotic feature in thermocapillary convection

The transition process from steady to turbulent convection via subharmonic bifurcation in thermocapillary convection of half floating zone was studied by numerical simulation and experimental test. Both approaches gave structure of period doubling bifurcations in the present paper, and the Feigenbaum universal law was checked for the system of thermocapillary convection.

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.

Symmetry of chaotic attractors and kink-antikink structures in CMLs

In this paper the symmetries of coupled map lattices (CMLs) and their attractors are investigated by group and dynamical system theory, as well as numerical simulation, by means of which the kink-antikink patterns of CMLs in space-amplitude plots are discussed.

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.

Excitation mechanism of thermocapillary oscillatory convection

It is suggested that the oscillation of thermocapillary convection may be excited by the buoyancy instability. By means of numerical simulation of the finite-element method, the temperature distributions in the liquid bridge are qualitatively analyzed. The temperature gradient in a certain flow region of liquid bridge may turn to be parallel to the direction of gravity when the temperature difference △T between two boundary rods of liquid bridge is larger than the critical value. The buoyancy instability may be excited, and then the thermocapillary oscillatory convection appears, as the temperature difference increases further. The distribution of the critical Marangoni number in the micro-gravity environment is derived from the data on the ground experiments. The results show that the onset of thermocapillary oscillatory convection is delayed in the case of smaller typical scale of liquid bridge and lower gravity environment.

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.

The induction plasma chemical reactor: Part I. Equilibrium model

A mathematical model is presented for the numerical simulation of the flow, temperature, and concentration fields in an rf plasma chemical reactor. The simulation is performed assuming chemical equilibrium. The extent of validity of this assumption is discussed. The system considered is the reaction of SiCl4 and NH3 for the production of Si3N4.

T形微通道中互不相溶两相流数值模拟

采用摄动有限体积(PFV)算法和水平集(level set)技术对T形微通道内互不相溶两相流动进行了数值模拟研究。考察了两相界面张力和微通道壁面润湿性对流动的影响,精确地捕捉到了油水两相流动的界面。对一些典型的T形微通道油水两相流动进行了数值计算,模拟结果和实验结果吻合较好。分析总结出了微通道内两相流动过程中的一些基本规律,为微通道内的液液两相流动实验设计和工业应用提供了新的数值预测手段.

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

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

超声速燃烧的湍流流场数值模拟

选择了6个有代表性的湍流模型，对超声速燃烧湍流流场开展了数值模拟。采用的模型包括5个线性涡粘性模型和1个非线性模型，在非线性模型中同时引入了可压缩修正。数值结果表明，本文发展的数值方法可以成功地应用于求解可压缩、含化学反应的多组分N—S方程。同时引入可压缩性修正的二阶涡粘性模型给出的结果较好。

带有相变的三维温度场数值模拟

以激光熔凝表面强韧化处理为背景，应用空间的弹塑性有限单元和高精度的数值算法、同时考虑材料组织性能的变化来模拟材料的温度场。主要研究激光熔凝加工中瞬时温度场数值模拟，同时考虑相变潜热的影响，为第二步热应力场及残余应力的数值模拟做准备。最后用算例验证了模型的正确性，并给出了不同时刻温度场的分布。

高寒草甸生态系统陆地生物圈模式研究及应用

改进了适合于高寒草甸生态系统的陆地生态圈模式，分析了模式中温度变化与水分运动分层的物理原因，说明了气候状况对地表面能量交换的影响，给出了净辐射和蒸散量新的计算方法，提出了有很差分计算中具有二阶精度的Euler隐式格式，介绍了中国科学院海北高寒草甸生态站的气候概况和野外观测情况。最后利用本模式对高寒草甸生态站地区的土壤-植被-大气间水热交换过程进行了数值模拟，模拟值与实测值吻合较好。

航天器贮箱气液自由界面追踪数值模拟

主要讨论了航天器贮箱在轨道航行时的微重力状态下其液体推进剂在贮箱内的形态分布及控制．文中采用VOF方法，加入了表面张力的效应，追踪气液两相流的自由界面，对液面在微重力条件下的位形变化进行了数值模拟．通过比较不同重力加速度及接触角下的两相流的相图，分析了影响贮箱中液体推进剂位形变化的主要因素及对其有效的控制方法．

航天运载器尾舱高空热环境数值分析

航天运载器在70km以上高空高速飞行，火箭发动机的喷流是影响其尾舱热环境的主要因素之一。采用代数方程和偏微分方程相结合的多块网格生成方法，达到热流计算中对近壁面网格的正交性和分布上的要求；通过改进的NND有限体积法高超程序，数值耦合求解多组分NS方程、湍流方程、热辐射方程和非平衡化学反应方程，提高了壁面热流计算的精度。计算结果与飞行遥测及地面实验结果对比分析表明，所采用的计算方法及得到的结果是正确有效的。