Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

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 devel- oped 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.

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

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.

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

采用计算流体动力学（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.

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

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