小波浪高雷诺数时串列圆柱的波动压力

在２４０×１０５≤Ｒｅ≤５４０×１０５，２９４≤Ｋｃ×７８５，２５４×Ｌ／Ｄ×５０２的范围内实验研究了串列圆柱在规则波中迎浪前进时的波动压力、波动升力与波动阻力．实验模型由ＰＶＣ管加工制成，表面光滑．实验结果表明：在波流中，当Ｋｃ数较小而Ｒｅ数较高时，绕串列圆柱的流动与均匀流中的绕流有相似之处；随着Ｋｃ数和Ｒｅ数的增加，串列圆柱的柱间干扰效应会发生改变     

计算机磁头／磁盘超薄气膜润滑压强的算子分裂算法

以任意拉森数的超薄气体润滑方程为基础，采用算子分裂法和非结构三角网格的有限元法计算Ω型磁头的空气轴承气垫面(ABS)气膜压强分布、气浮力和纵倾力矩；在分析比较流量系数的各种算法的基础上，确定采用多项式拟合数据库计算流量系数．计算结果表明：采用算子分裂法可以有效克服在高轴承数时的数值不稳定性，消除数值振荡；在小轴承数时，气浮力随轴承数增加而增大，当轴承数增大到某一数值后，气浮力趋近某一稳定值，此时气膜压强分布与磁头造型基本一致；气膜的纵倾力矩在轴承数的某一临界值附近出现最大值．     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

绕旋转圆柱流动涡尾流结构和临界状态特性

采用作者提出的基于区域分解，有限差分法与涡法杂交的数值方法，结合高阶隐式差分格式，和以修正的不完全ＬＵ分解为预处理器的共轭梯度法作求解器，系统地研究了雷诺数Ｒｅ＝１０００，旋转速度比α∈（０．５，３．２５）范围内，绕旋转圆柱从突然起到充分发展，长时间内尾流旋涡结构和阻力，升力系数的变化规律，计算所得流动图案与实验流场显示符合很好。数值试验证帝了临界状态的存在，并首次给出了临界状态时的旋涡结构特性。     

应用脉冲热线技术测量含回流的高湍流流动

阐述脉冲热线测量原理、探针结构与数据处理方法；在实验雷诺数Ｒｅｃ＝６．６７×１０＾５时，用脉冲热线技术定量测量了 由楔形体模型产生的含回流的高湍流流动，得到了壁面切应力、流向时均速度、回流因子、脉动速度、偏斜因子和平坦因子等重要参数，提示了分离泡流动结构与特征。     

ON THE SECOND ORDER WAVE DIFFRACTION IN TWO LAYER FLUIDS

ＯＮ　ＴＨＥ　ＳＥＣＯＮＤ　ＯＲＤＥＲ　ＷＡＶＥ　ＤＩＦＦＲＡＣＴＩＯＮ　ＩＮ　ＴＷＯ　ＬＡＹＥＲ　ＦＬＵＩＤＳＷｕＪｉａｎｈｕａ（吴建华）；ＦａｎｇＹｉｎｇ（方颖）（ＲｅｃｅｉｖｅｄＭａｙ４，１９９５；ＲｅｓｉｖｅｄＪｕｎ．２１，１９９６；Ｃｏｍｍ...     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

平面Poisenille流动的二维Hopf分叉

本文用直接数值模拟的方法计算了二维Ｐｏｉｓｅｕｉｌｌｅ流动中扰动波的演化问题。得到了二维平衡态，在一定的波数下，Ｒｅ３９５０时，这种平衡态，将变得不稳定，模拟发现出现第二周期解，即二次分叉。     

Numerical-Analytical Investigation of Multimodal Solutions of the Jeffery-Hamel Problem for a Convergent Channel

In the Jeffery-Hamel problem for a convergent channel, for small, moderately large and asymptotically large Reynolds numbers Re, multimodal, symmetrical and asymmetrical solutions not available in the literature are obtained and studied for the entire range of admissible convergence angles. For asymptotically small Re, these solutions depend singularly on Re. A number of novel hydrodynamic effects is established and commented on.     

Second-order finite difference solutions for the flow between rotating concentric spheres

This paper describes a second-order method to calculate approximate solutions to flow of viscous incompressible fluid between rotating concentric spheres. The governing partial differential equations are presented in the stream–vorticity formulation and are written as a series of second-order equations. The technique employed makes use of second-order approximations for all terms in the governing equations and is dependent upon the direction of flow at a given point. This upwind technique has allowed us to generate approximate solutions with larger Reynolds numbers than has generally been possible for second and higher-order techniques. Solutions have been obtained with Reynolds numbers as large as 3000 and with grids as fine as a 40 × 40 mesh. Results are displayed in the form of level curves for both the stream and vorticity functions. A dimensionless quantity related to the torque acting on both spheres has been calculated from the approximate solution and compared with other results. Results with smaller Reynolds numbers such as 100 and 1000 are in excellent agreement with other published results.     

Constructing solutions of the navier-stokes equations for a fluid layer between moving parallel plates at low and moderate reynolds numbers

In this paper, we study exact solutions of the Navier-Stokes equations for a layer between parallel plates, the distance between which varies according to an arbitrary power law and whose boundary has a no-slip condition, are under study. A solution in the form of a power series of the Reynolds number is obtained. Comparison with the exact solutions is performed, and high accuracy of expansions for Reynolds numbers Re = 1–10 is shown. An accurate estimate of the error of the Reynolds thin layer approximation is obtained.     

Quasiperiodic solutions of the Navier-Stokes equations induced by the quasiperiodic shape of the boundaries of a two-dimensional flow domain

Steady quasiperiodic solutions of the Navier-Stokes equations in an infinite two-dimensional layer with quasiperiodic boundaries are obtained on the Reynolds number range 0 < Re* < 200. The calculations are performed using a spectral-difference method based on the representation of the quasiperiodic solutions in the form of convergent double Fourier series. The properties of these solutions and the distinctive features of their spectra are studied and their fundamental differences from periodic solutions are shown. The possibility of applying the quasiperiodic solutions for modeling flows in fractal layers is discussed.     

Large Eddy Simulation of a plane impinging jet

Large Eddy Simulations of a plane turbulent impinging jet have been carried out using the dynamic Smagorinsky model. The statistical results are first validated with the measurements from the literature: mean and turbulent quantities along the jet axis and at different vertical locations are presented. This study is completed by the analysis of the wall shear stress at the impingement wall. The effect of the jet Reynolds number (3000Re13500) on the kinematic development of the jet is also discussed. To cite this article: F. Beaubert, S. Viazzo, C. R. Mecanique 330 (2002) 803–810.     

鱼类自主游动水动力学特性的数值模拟

通过对推力和阻力进行重新定义, 从根本上解决了鱼游研究中推力和阻力无法区分的难题.在此基础上, 利用自适应网格下的ghost-cell浸没边界方法, 模拟了鱼类以鲹科模式在黏性流体(309 \le Re \le 14\,581)和无黏流体 (相当于雷诺数无穷大情形)中的二维自主游动.结果表明: (1) Strouhal数随雷诺数增大而减小,当雷诺数趋向于无穷时, Strouhal数趋向于0.25; (2)在所有雷诺数情况下, 推力主要来源于压力分量; 当Re<3000时, 阻力的压力分量小于黏性力分量, 而当Re>3000, 二者的关系就会反过来; (3)推进效率随着雷诺数的增大而增大,当雷诺数趋向于无穷大时, 推进效率最高可以达到70%, 说明鲹科模式适用于较高雷诺数下的游动.      

On the Reynolds-number independence of mixed convection in a vertical channel subjected to asymmetric wall temperatures with and without flow reversal

The problem of mixed convection in a vertical channel with asymmetric wall temperatures including situations of flow reversal is studied numerically. The SIMPLER algorithm with a staggered grid system is employed to solve the corresponding numerical equations formulated by the finite-volume method. A second-order upwind scheme is used to model the convective term, and a suitable grid distribution is introduced. The ranges of the parameters studied are 0 r_t 1, 1 Re 1000, and 0 Gr/Re 500.

The numerical results, with the streamwise coordinate scaled by the Reynolds number (Re), show that solutions for the velocity and temperature fields are independent of the Reynolds number when Re 50, even in the presence of flow reversal. These solutions, however, are dependent on r_t and Gr/Re. Subsequently, correlations are proposed for the bulk temperature distribution and the local Nusselt numbers along the hot wall and the cold wall.     

Mixed convection about a non-isothermal vertical surface in a porous medium

The problem of mixed convection about non-isothermal vertical surfaces in a saturated porous medium is analysed using boundary layer approximations. The analysis is made assuming that the surface temperature varies as an arbitrary function of the distance from the origin. A perturbation technique has been applied to obtain the solutions. Using the differentials of the wall temperature, which are functions of distance along the surface, as perturbation elements, universal functions are derived for various values of the governing parameter Gr/Re. Both aiding and opposing flows are considered. The universal functions obtained can be used to estimate the heat transfer and fluid velocity inside the boundary layer for any type of wall temperature variation. As a demonstration of the method, heat transfer results have been presented for the case of the wall temperature varying as a power function of the distance from the origin. The results have been studied for various combinations of the parameters Gr/Re and the power index m, taking both aiding and opposing flows into consideration. On comparing these results with those obtained by a similarity analysis, the agreement is found to be good.     

Simulation of peristaltic flow of chyme in small intestine for couple stress fluid

In this article couple stress fluid have been considered for the peristaltic flow of chyme in intestine. Problem under consideration have been formulated assuming that two non-periodic sinusoidal wave of different wavelength propagate with same speed c along the outer wall of the tube. Governing equations have been simplified under the assumptions of long wavelength and low Reynolds number approximation (such assumption are consistent that Re (Reynold number) is very small and long wavelength approximation also exists in the small intestine). Exact solutions have been evaluated for velocity and pressure rise. Physical behaviour of different parameter of couple stress fluid have been presented graphically for velocity, pressure rise, pressure gradient and frictional forces. The stream lines are also made against different parameters.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.