Operational Optimization of GaN Thin Film Growth Employing Numerical Simulation in a Showerhead MOCVD Reactor

A detailed reaction-tran sport model was studied in a showerhead reactor for metal organic chemical vapor deposition of GaN film by using computational fluid dynamics simulation. It was found that flat flow lines without swirl are crucial to improve the uniformity of the film growth, and thin temperature gradient above the suscptor can increase the film deposition rate. By above-mentioned research, we can employ higher h (the distance from the susceptor to the inlet), P (operational pressure) and the rate of susceptor rotation to improve the film growth.

集热栅对太阳能烟囱发电系统效率的影响

采用计算流体动力学（CFD）方法对太阳能烟囱发电装置进行数值模拟，得到装置内部的温度场、速度场、压力场等分布情况。对集热棚的各种几何和物理参数进行研究和分析。结果表明，集热棚直径、太阳辐照强度、覆盖材料的透明度等诸多参数对系统效率有直接而重要的影响。

Space grown semi-insulating gallium arsenide single crystal and its application

Low noise field effect transistors and analogue switch integrated circuits (ICs) have been fabricated in semi-insulating gallium arsenide (SI-GaAs) wafers grown in space by direct ion-implantation. The electrical behaviors of the devices and the ICs have surpassed those fabricated in the terrestrially grown SI-GaAs wafers. The highest gain and the lowest noise of the transistors made from space-grown SI-GaAs wafers are 22.8 dB and 0.78 dB, respectively. The threshold back-gating voltage of the ICs made from space-grown SI-GaAs wafers is better than 8.5 V The con-elation between the characterizations of materials and devices is studied systematically. (C) 2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

CFD analysis of pneumatic conveying in a double-tube-socket (DTS (R)) pipe

A pipeline with a bypass is widely used for the pneumatic conveying of material. The double-tube-socket (DTS (R)) technology, which uses a special inner bypass, represents the current state of the art. Here, a new methodology is proposed based on the use of computational fluid dynamics (CFD) to predict the energy consumption of DTS conveying. The predicted results are in good agreement with those from pilot-scale experiments. (C) 2010 Elsevier Inc. All rights reserved.

三维风沙运动的CFD-DEM数值模拟

采用计算流体力学(Computational Fluid Dynamics,CFD)和颗粒离散元(Discrete Element Method,DEM)耦合的方法模拟三维风沙运动,并且将三维模拟结果和二维模拟结果以及实验结果进行了对比.计算结果表明:沙粒水平速度随着高度按幂函数规律增加,沙床表面附近沙粒撞击和起跳速度的概率分布均可用对数正态函数描述,沙粒撞击和起跳角度的概率分布均可用指数函数描述,沙粒水平速度,展向速度和垂直速度在不同高度处的概率分布可分别用对数正态分布,正态分布和正态分布表示.与二维计算结果的分析对比表明:二维计算得到的颗粒速度的分布规律和三维计算结果类似,但二维计算的颗粒表观密度明显偏大,由此导致输沙量计算偏大.和实验结果的对比表明:三维计算得到的颗粒速度概率分布与实验基本保持一致

Speed of the internal pellet target in CSRm

Pellet target is one of the main candidate targets in CSRm (cooler storage ring’s main ring) for hadron physics studies. Pellet speed is an important physical parameter for the target. Larger pellet speed could shorten the interacting time interval between the pellet and the cyclotron beam, and thus results in a small temperature variation for the pellet. This could make the pellet facility work in a stable condition. A fluid dynamic simulation was carried out for the pellet speed, and it was found that the maximum speed for the target pellet may be restricted to about 100 m/s even if all working parameters were set to their optimal values.

Speed of the internal pellet target in CSRm

Pellet target is one of the main candidate targets in CSRm (cooler storage ring’s main ring) for hadron physics studies. Pellet speed is an important physical parameter for the target. Larger pellet speed could shorten the interacting time interval between the pellet and the cyclotron beam, and thus results in a small temperature variation for the pellet. This could make the pellet facility work in a stable con-dition. A fluid dynamic simulation was carried out for the pellet speed, and it was found that the maxi-mum speed for the target pellet may be restricted to about 100 m/s even if all working parameters were set to their optimal values.

一种基于粒子的牛顿流体与粘弹性流体统一模拟方法

在越来越受到人们关注的基于物理流体动画领域,目前分别模拟牛顿流体或粘弹性流体的方法很多,但很少有统一模拟两者的方法.文中基于光滑粒子流体动力学方法,通过对传统纳维-斯托克斯方程添加弹性应力项,提出了一种新的统一模拟牛顿流体和粘弹性流体的方法.通过实验说明该方法不仅有效,易于实现,而且具有良好的可控性,仅仅通过调节参数就可以模拟不同粘弹性、不同类型的流体现象.

Self-organized Criticality Model for Ocean Internal Waves

In this paper, we present a simple spring-block model for ocean internal waves based on the self-organized criticality (SOC). The oscillations of the water blocks in the model display power-law behavior with an exponent of -2 in the frequency domain, which is similar to the current and sea water temperature spectra in the actual ocean and the universal Garrett and Munk deep ocean internal wave model [Geophysical Fluid Dynamics 2(1972) 225; J. Geophys. REs. 80 (1975) 291]. The influence of the ratio of the driving force to the spring coefficient to SOC behaviors in the model is also discussed.

Numerical simulation of principal tidal constituents in the South China Sea, Gulf of Tonkin and Gulf of Thailand

The principal tidal constituents M-2, S-2, K-1 and O-1 in the South China Sea, Gulf of Tonkin and Gulf of Thailand are simulated simultaneously using the numerical scheme of Kwok et al. (1995 Proceedings of the 1st Asian Computational Fluid Dynamics Conference, pp. 16-19). The average differences between the computed and observed harmonic constants are mostly within 5 cm and 10 degrees for amplitudes and phase-lags, respectively. The simulated tidal regimes in the present model are believed to be more accurate than the previous numerical results. Our studies confirm that a clockwise rotating M-2 amphidromic system lies in the southeast of the Gulf of Thailand and an S-2 amphidromic system at the near-shore area of the northeast South China Sea. The linear tidal energy equation developed by Garrett (1975 Deep-Sea Research 22, 23-35) is generalized to the nonlinear case. Based on the numerical results, the energy budgets in the South China Sea and its subareas, namely the Taiwan Strait, the Gulf of Tonkin, the Gulf of Thailand and the remaining area are investigated. The tidal motion in the Taiwan Strait is maintained mainly by the energy fluxes from the East China Sea for both semidiurnal and diurnal species and partially from the Luzon Strait for semidiurnal species. For the other parts of the South China Sea, the tidal motion is mainly maintained by the energy fluxes through the Luzon Strait. The energy inputs from the tide-generating force are negative for semidiurnal species and positive for diurnal species. (C) 1999 Elsevier Science Ltd. All rights reserved.

Rossby waves with linear topography in barotropic fluids

Rossby waves are the most important waves in the atmosphere and ocean, and are parts of a large-scale system in fluid. The theory and observation show that, they satisfy quasi-geostrophic and quasi-static equilibrium approximations. In this paper, solitary Rossby waves induced by linear topography in barotropic fluids with a shear flow are studied. In order to simplify the problem, the topography is taken as a linear function of latitude variable y, then employing a weakly nonlinear method and a perturbation method, a KdV (Korteweg-de Vries) equation describing evolution of the amplitude of solitary Rossby waves induced by linear topography is derived. The results show that the variation of linear topography can induce the solitary Rossby waves in barotropic fluids with a shear flow, and extend the classical geophysical theory of fluid dynamics.

CFX与USAERO的水下机器人操纵性仿真计算研究

研究了借助计算流体力学软件CFX与USAERO对水下机器人的操纵性进行仿真计算的方法,并将仿真计算结果与模型试验结果对比,验证仿真计算的可行性,指出仿真计算的适用范围.该方法可以对某些现实中难以进行的操纵性试验进行仿真计算,对于研究水下机器人的操纵性尤其是在方案设计阶段的操纵性设计和降低试验成本有一定的现实意义.

APPLICATION OF THE KHOKHLOV-ZABOLOTSKAYA-KUZNETSOV EQUATION TO MODELING HIGH-INTENSITY FOCUSED ULTRASOUND BEAMS

High-intensity focused ultrasound is a form of therapeutic ultrasound which uses high amplitude acoustic waves to heat and ablate tissue. HIFU employs acoustic amplitudes that are high enough that nonlinear propagation effects are important in the evolution of the sound field. A common model for HIFU beams is the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation which accounts for nonlinearity, diffraction, and absorption. The KZK equation models diffraction using the parabolic or paraxial approximation. Many HIFU sources have an aperture diameter similar to the focal length and the paraxial approximation may not be appropriate. Here, results obtained using the “Texas code,” a time-domain numerical solution to the KZK equation, were used to assess when the KZK equation can be employed. In a linear water case comparison with the O’Neil solution, the KZK equation accurately predicts the pressure field in the focal region. The KZK equation was also compared to simulations of the exact fluid dynamics equations (no paraxial approximation). The exact equations were solved using the Fourier-Continuation (FC) method to approximate derivatives in the equations. Results have been obtained for a focused HIFU source in tissue. For a low focusing gain transducer (focal length 50λ and radius 10λ), the KZK and FC models showed excellent agreement, however, as the source radius was increased to 30λ, discrepancies started to appear. Modeling was extended to the case of tissue with the appropriate power law using a relaxation model. The relaxation model resulted in a higher peak pressure and a shift in the location of the peak pressure, highlighting the importance of employing the correct attenuation model. Simulations from the code that were compared to experimental data in water showed good agreement through the focal plane.