中等雷诺数圆柱尾流旋涡脱落的控制

介绍通过实验对圆柱尾流旋涡脱落进行抑制的方法及其结果。实验模型的展径比为38，实验的雷诺数范围为3 * 10~2～1.6 * 10~3。抑制方法是在圆柱（直径为D）表面沿展向每隔一定间距伸出一直径0.18D、长度为1.5D的小棒，实验结果表明，当棒间距小于3D，棒与来流夹角在30 °～ 90 °范围内，可有效抑制旋涡脱落。

剪切水-气界面下湍流猝发特征的实验研究

利用氢气泡时间线-脉线组合示踪技术定量地考察剪切水-气界面下的湍流猝发现象,分析猝发事件的信号特征,重点探讨猝发与湍能产生之间的联系.在猝发过程中,水面近区的瞬时流速和Reynolds切应力出现较大幅度的脉动,它们在时间和空间垂直方向上表现出高度的相干性,这是猝发事件的一个显著特征.在猝发期,猝发事件涉及的空间区域内Reynolds切应力和湍流脉动强度明显比平均值和非猝发期的情况大.其结果表明:在所考察的实验条件下,猝发是剪切水-气界面附近湍流产生的主要过程.

超声速含灰气体钝体绕流的传热增强及颗粒惯性沉积效应

研究了高、中等Reynolds数条件下，稀相含灰气体绕球体的定常超声速流动问题，对于物体迎风面前表面上发生颗粒沉积的情况，研究了由于颗粒引起的驻点热流增加，给出了传热增强的最大值，指明了驻点热流的增加取决于流动Re数、颗粒惯性参数、气固比热之比以及物面温度等控制参数。

圆柱振荡绕流的三维不稳定性研究

通过数值求解三维不可压缩Navier-Stokes方程，研究了振荡圆柱绕流的旋涡不稳定性。研究表明，在一定的参数范围内，由于旋涡不稳定性，振荡流出由二维演化成三维流动，并沿圆柱轴向形成交错排列的三维涡结构。数值计算合理地预测了三维涡结构的空间失稳波长，并与实验测试值相符很好。文中还进一步研究了圆柱的受力特性，通过求解Morison方程，计算了圆柱的阻力和惯性力特性，其计算结果与已有的实验数据相吻合。

基于槽道流LES数据库的湍流猝发检测

采用大涡模拟方法,在与两种Reynolds数情形的DNS结果进行充分验证的基础上,获得了不同Reynolds数情形槽道湍流的可靠LES数据库,由此可进一步得到任意Reynolds数时速度剖面、湍流强度、剪应力等统计量的时空分布以及猝发结构的时空特征.基于这些可靠的LES数据库,利用条件采样方法检测猝发事件的时空尺度,并提出由喷发事件时间间隔概率分布曲线确定组合参数对传统的条件采样方法进行改进,以避免检测结果的误差.检测结果表明,引入组合参数后,湍流猝发周期对门限的依赖性得到显著改善.同时,对猝发事件的空间分布进行检测,得到了平均猝发面积比.通过比较不同Reynolds数的结果发现,Reynolds数对平均猝发周期和平均猝发面积比的影响不大.

可压衰减湍流中被动标量场的直接数值模拟及谱分析

运用高精度迎风差分方法及8阶精度群速度控制型差分格式(GVC8)对初始 Reynolds数为72，初始湍流Mach数0．2-43．9的可压衰减湍流的流场及被动标量场进行了直接数值模拟，被动标量场的Schmidt数为2~10．通过不同计算网格及不同方法的数值计算对本文结果进行了验证．指出了可压衰减湍流中被动标量能谱存在的一1律，压缩性效应使得高波数成分衰减加快．对被动标量场进行了标度律分析，发现可压湍流中被动标量场具有扩展自相似性，而压缩性效应对被动标量标度指数的影响较小．发现随着被动标量Schmidt数增加，其标度指数的奇异性增强．

Statistical Simulation of Thin-Film Bearings

The information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method are used to simulate the gas flows between the write/read head and the platter of the disk drive (the slider bearing problem). The results of both methods are in good agreement with numerical solution of the Reynolds equation in the cases studied. However, the DSMC method owing to the problem of large sample size demand and the difficulty in regulating boundary conditions at the inlet and outlet was able to simulate only short bearings, while IP simulates the bearing of authentic length ~1000 m ? and can provide more detailed flow information.

Two Bubbles Thermocapillary Migration Experimental System

An experimental investigation will be performed on the thermocapillary motion of two bubbles in Chinese return-satellite. The experiment will study the migration process of bubble caused by thermocapillary effect in microgravity environment, and their interaction between two bubbles. The bubble is driven by the thermocapillary stress on the surface on account on the variation of the surface tension with temperature. The interaction between two bubbles becomes significant as the separation distance between them is reduced drastically so that the bubble interaction has to be considered. Recently, the problem has been discussed on the method of successive reflections, and accurate migration velocities of two arbitrarily oriented bubbles were derived for the limit of small Marangoni and Reynolds numbers. Numerical results for the migration of the two bubbles show that the interaction between two bubbles has significant influence on their thermocapillary migration velocities with a bubble approaching another. However, there is a lack of experimental validate for the theoretic results. Now the experimental facility is designed for experimenting time after time. A cone-shaped top cover is used to expel bubble from the cell after experiment. But, the cone-shaped top cover can cause temperature uniformity on horizontal plane in whole cell. Therefore, a metal board with multi-holes is fixed under the top cover. The board is able to let the temperature distribution on the board uniform because of their high heat conductivity, and the bubble can pass through it. In the system two bubbles are injected into the test cell respectively by two sets of cylinder. And the bubbles sizes are controlled by two sets of step-by-step motor. It is very important problem that bubble can be divorced from the injecting mouth in microgravity environment. Thus, other two sets of device for injecting mother liquid were used to push bubble. The working principle of injecting mother liquid is to utilize pressure difference directly between test cell and reservoir

Information Preservation (IP) Method in Simulation

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 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 ( ) 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 ( ) 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.

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.

Statistics of Turbulence in Benthal Boundary Layers

This paper deals with turbulence behavior inbenthalboundarylayers by means of large eddy simulation (LES). The flow is modeled by moving an infinite plate in an otherwise quiescent water with an oscillatory and a steady velocity components. The oscillatory one aims to simulate wave effect on the flow. A number of large-scale turbulence databases have been established, based on which we have obtained turbulencestatisticsof the boundarylayers, such as Reynolds stress, turbulence intensity, skewness and flatness ofturbulence, and temporal and spatial scales of turbulent bursts, etc. Particular attention is paid to the dependences of those statistics on two nondimensional parameters, namely the Reynolds number and the current-wave velocity ratio defined as the steady current velocity over the oscillatory velocity amplitude. It is found that the Reynolds stress and turbulence intensity profile differently from phase to phase, and exhibit two types of distributions in an oscillatory cycle. One is monotonic occurring during the time when current and wave-induced components are in the same direction, and the other inflectional occurring during the time when current and wave-induced components are in opposite directions. Current component makes an asymmetrical time series of Reynolds stress, as well as turbulence intensity, although the mean velocity series is symmetrical as a sine/cosine function. The skewness and flatness variations suggest that the turbulence distribution is not a normal function but approaches to a normal one with the increasing of Reynolds number and the current-wave velocity ratio as well. As for turbulent bursting, the dimensionless period and the mean area of all bursts per unit bed area tend to increase with Reynolds number and current-wave velocity ratio, rather than being constant as in steady channel flows.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Single curved fiber sedimentation under gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

DNS and LES of turbulent channel flow with hydrophobic surface

Hydrophobic surface benefits for drag reduction. Min and Kim[1] do the first Direct Numerical Simulation on drag reduction in turbulent channel flow. And Fukagata and Kasagi[2] make some theoretical analysis based on Dean[3]'s formula and some observations in the DNS results. Using their theory, they conclude that drag reduction is possible in large Reynolds number. Both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed in our research. How the LES behaving in the turbulent channel flow with hydrophobic surface is examined. Original Smagorinsky model and its Dynamical model are used in LES. The slip velocities predicted by LES using Dynamical model are in good agreement with DNS as shown in the Figure. Although the percentage of drag reduction predicted by LES shows some discrepancies, it is in the error limit for industrial flow. First order and second order moments of LES are also examined and compared with DNS's results. The first-order moments is calculated well by LES. But there are some discrepancies of second-order moments between LES and DNS. [GRAPHICS]