微尺度串列双圆柱绕流的分子动力学并行模拟

本文采用基于空间分解算法的分子动力学并行模拟方法,研究了微尺度、低雷诺数（Re=40）下串列等大的双圆柱绕流现象。结果表明：随着间距比L＊／D＊的增加,流动存在3种特征状态：当L＊／D＊〈1．1时,同单一物体的绕钝体流动相似;当1．1〈L＊／0＊〈3．5时,涡脱落现象只在下游圆柱出现,在两圆柱之间有交替附着于下游圆柱的...     

Flow patterns of cross-flow around a varicose cylinder

The near wake of a varicose cylinder has been experimentally investigated using Laser Induced Fluorescence (LIF) and Digital Particle Imaging Velocimetry (DPIV). The work aims to provide understanding to the mechanism of the cross flow around varicose cylinder as well as to comprehend why the introduction of relatively small degrees of spanwise waviness can have a significant effect on drag reduction and suppression of the cylinder vibration. The evolution of the flow patterns and the corresponding vortex interactions are obtained. The experimental results indicated that the wake width and the formation length vary along the span of the varicose cylinder. A wider wake and a longer formation length were observed in the saddle plane. In addition, an interpretation of the three-dimensional wake structures is postulated and conceptually shown. The numerical simulation by 3-D finite volume method is successful in predicting the flow features found by the experiments.     

Lattice Boltzmann Simulation of the Cross Flow Over a Cantilevered and Longitudinally Vibrating Circular Cylinder

The multiple-relaxation-time lattice Boltzmann method （MRT-LBM） is implemented to numerically simulate the cross flow over a longitudinal vibrating circular cylinder. This research is carried out on a three-dimensional （3D） finite cantilevered cylinder to investigate the effect of forced vibration on the wake characteristics and the 3D effect of a cantilevered cylinder. To meet the accuracy of this method, the present calculation is carried out at a low Reynolds number Re =100, as well as to make the vibration obvious, we make the vibration strong enough. The calculation results indicate that the vibration has significant influence on the wake characteristics. When the vibrating is big enough, our early works show that the 2D vortex shedding would be locked up by vibration. Contrarily, this phenomenon would not appear in the present 3D case because of the end effect of the cantilevered cylinder.     

气固两相三维圆柱绕流的直接数值模拟

本文采用直接数值模拟(DNS)方法研究了气周两相三维圆柱绕流的涡量场和颗粒扩散,并着重讨论了圆柱绕流中卡门涡街的形成和涡结构的转捩过程。同时分析了圆柱绕流中不同Stokes数的颗粒在涡结构作用下的横向扩散。结果显示Stokes数为1的颗粒主要分布在流场大尺度涡结构的外边界,而Stokes数为0．01的颗粒在涡结构的作用下,在流场中充分混合。     

A numerical study of transient flow around a cylinder and aerodynamic sound radiation

The fully 3D turbulent incompressible flow around a cylinder and in its wake at a Reynolds number Re = = 9×10⁴ based on the cylinder diameter and Mach number M = 0.1 is calculated using Large Eddy Simulations (LES). Encouraging results are found in comparison to experimental data for the fluctuating lift and drag forces. The acoustic pressure in far-field is commutated through the surface integral formulation of the Ffowcs Williams and Hawkings (FWH) equation in acoustic analogy. Five different sound sources, the cylinder wall and four permeable surfaces in the flow fields, are employed. The spectra of the sound pressure are generally in quantitative agreement with the measured one though the acoustic sources are pseudo-sound regarding the incompressible flow simulation. The acoustic component at the Strouhal number related to vortex shedding has been predicted accurately. For the broad band sound, the permeable surfaces in the near wake region give qualitative enough accuracy level of predictions, while the cylinder wall surface shows a noticeable under-prediction. The sound radiation of the volumetric sources based on Lighthill tensors at vortex shedding is also studied. Its far-field directivity is of lateral quadrupoles with the weak radiations in the flow and cross-flow directions.     

过渡流下近壁插入圆柱对壁面强化传热的影响

本文利用基于复合网格系统的计算方法,对Re=50～1200的近壁插入圆柱流场进行数值模拟,研究过渡流状态下在壁面附近插入圆柱对下游壁面传热强化的影响。并基于低速循环水槽流动实验台,采用粒子成像测试法(PIV)对Re=100～500的近壁插入圆柱流场进行可视化实验研究,验证了数值模拟方法的可靠性。研究结果表明:近壁插入圆柱流场在Re=100时进入过渡流状态;Re直接影响圆柱尾流中周期性涡脱和壁面涡岛的发生位置及其洗刷效应的大小,随着Re的增大,洗刷效应明显增强,因而,过渡流范围内Re越大,圆柱下游壁面传热强化越大。     

具有多孔介质覆面层的柱体绕流流场分析

建立多孔结构覆面柱体绕流模型, 采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对覆盖多孔介质层的方柱绕流进行数值模拟, 研究多孔介质对钝体绕流流场特性的影响。结果表明: 相比于不可渗透壁的柱体, 引入合适参数的多孔介质覆面层后可以有效降低其升力脉动幅值, 但阻力有所增加。同时, 较高雷诺数下多孔方柱的数值模拟表明: 多孔介质壁面使得尾迹区域的剪切层相距更远, 降低了尾流处湍动能, 并将雷诺应力的峰值移动到尾迹区域, 抑制了方柱两侧的动量交换, 使动量交换的位置发生在尾迹区域, 继而使得尾迹的涡街更加规则化。     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.     

Modification of the wake behind a circular cylinder by using synthetic jets

Experimental tests were conducted to control the flow around a cylinder by means of unsteady blowing (synthetic jet) through a single slot disposed on the wall of the model. The flow Reynolds number (based on the diameter of the model) wasR _D=10₅. The efficiency of the synthetic jet is quantified in terms of delaying separation and modifying the drag coefficient. The investigations were of three types: measurements of the mean pressure distribution, wall visualizations of the separation line position and measurements of the mean flow-field in the wake.     

Modeling the flow of liquid-metal coolant in the T-shaped mixer

The results of experimental studies on the structure of the temperature field in the tube cross section at the flow of liquid-metal coolant in a T-shaped mixer are presented. Experiments were carried out using the Rose alloy as the working fluid. To determine temperature distribution on the test section wall, infrared thermography was used; to determine temperature distribution in the channel cross section, a mobile thermocouple was used. Considerable temperature maldistribution in the mixing zone of liquid flows with different temperatures on the tube wall and in the coolant melt is shown.     

Reynolds数22 000的孤立方柱绕流的大涡模拟

方柱绕流是典型的钝体绕流问题,蕴含了丰富的流体力学现象,对这类流动的准确预测面临着诸多挑战.采用自主发展的大涡模拟程序,对来流Mach数M=0.3,Reynolds数Re_D=22 000的绕孤立方柱流动进行了细致模拟,亚格子模型使用动力涡黏模型.对计算结果的分析表明,大涡模拟所得的平均流场及Reynolds应力分布与已有实验数据和直接数值模拟结果均吻合较好,验证了预测结果的可靠性;在此基础上对瞬态流场进行了研究,展示了计算条件下方柱绕流分离转捩及尾迹区旋涡交替脱落形成Karman涡街的全过程,为更细致的流动机理探索奠定了基础.      

Effect of elasticity of a polymer solution on the Hele-Shaw flow

In this study, the Hele-Shaw cell is used to examine the effect offluid elasticity on the flow patterns of two-dimensional potential flow. Flows around a circular cylinder, a square cylinder and flows through abruptly converging-diverging channels (slits) with different throat lengths are tested for water and 0.2 wt % polyacrylamide aqueous solution (PAA-solution). The viscosity of the latter is well modeled by the power law, and the first normal stress difference in the steady shear flow is around ten times higher than the shear stress. Although the PAA-solution is highly shear-thinning, the flows of PAA-solution well reproduce the two-dimensional potential flow patterns that correspond to the respective flow configurations when the flow rate is very low. The potential flow patterns ofPAA solution are disturbed in the opposite way of inertia effect observed for water. The streamlines near the upstream stagnation point of cylinders are shifted upstream separating from the cylinder surface when the flow rate is higher, while streamlines in the wake approach closer to the downstream stagnation point. Streamlines offlow through the slit at flow rates higher than the potential flow region show that a pair ofvortices is formed upstream the slit entrance, while the streamlines remain attached to the downstream wall after passing the slit.     

Numerical Simulation of Two-Dimensional Flow over Three Cylinders by LatticeBoltzmann Method

The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the movingcylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D= 0.5 and f_ν=0.4 (where A is the forced amplitude, D is the cylinder diameter, and f_ν stands for the vibrating frequency, respectively). The calculated results not onlyindicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).     

Detonation in a Three-Dimensional Elbowed Channel

The results of simulation of detonation in a curved three-dimensional channel with a circular cross section of constant width blown through by a supersonic flow of a stoichiometric propane?air mixture are presented. In the bending zone, the channel wall was toroidal. The study was carried out within the framework of the one-stage combustion kinetics by the numerical method based on Godunov’s scheme in the original software package developed for multiparameter calculations and visualization of flows. The initiation of detonation occurs as a result of the formation of the shockwave configurations associated with the flow turn in the channel. Unsteady flow patterns are obtained, and their dependence on the parameters of the problem is investigated. The flow regime without detonation, the mode with the detonation wave emerging from the channel through the input cross section, and the mode with steady detonation are obtained.     

Direct numerical simulation of the flow around a wall-mounted square cylinder under various inflow conditions

The flow around a wall-mounted square cylinder of side d is investigated by means of direct numerical simulation (DNS). The effect of inflow conditions is assessed by considering two different cases with matching momentum-thickness Reynolds numbers Re_θ ? 1000 at the obstacle: the first case is a fullyturbulent zero pressure gradient boundary layer, and the second one is a laminar boundary layer with prescribed Blasius inflow profile further upstream. An auxiliary simulation carried out with the pseudo-spectral Fourier–Chebyshev code SIMSON is used to obtain the turbulent time-dependent inflow conditions which are then fed into the main simulation where the actual flow around the cylinder is computed. This main simulation is performed, for both laminar and turbulent-inflows, with the spectral-element method code Nek5000. In both cases the wake is completely turbulent, and we find the same Strouhal number St ? 0.1, although the two wakes exhibit structural differences for x > 3d downstream of the cylinder. Transition to turbulence is observed in the laminar-inflow case, induced by the recirculation bubble produced upstream of the obstacle, and in the turbulent-inflow simulation the streamwise fluctuations modulate the horseshoe vortex. The wake obtained in our laminar-inflow case is in closer agreement with reference particle image velocimetry measurements of the same geometry, revealing that the experimental boundary layer was not fully turbulent in that dataset, and highlighting the usefulness of DNS to assess the quality of experimental inflow conditions.     

Identification of Regions with Inhomogeneous Particle Behavior in Dilute Gas‐solid Flows

The behavior of the particulate phase in a highly turbulent gas flow has been investigated in a vertical channel. Variations of the flow configuration (1. Flow past a cylinder, 2. flow past a wall‐mounted obstacle and 3. flow around a horizontally injected jet) have been subject to both experiments and numerical simulations. The velocity vector field of the solid phase has been measured by digital particle image velocimetry (DPIV). The measurements have been focused on particle‐obstacle collisions and crossflow in the vicinity of the jet nozzle using the lately developed twinpeak detection method. By application of this method regions of highly inhomogeneous particle behavior could be detected mainly upstream of the flow perturbation. Numerical results have been obtained by an Eulerian‐Lagrangian method on boundary‐fitted grids. Particle‐particle interactions as well as interphase exchange of momentum have been taken into account. The simulation results showed to be well in accordance with the experimental results.     

基于转捩SST模型凸起圆柱绕流数值研究

为了研究局部凸起对边界层转捩的影响,采用转捩SST模型分别对亚临界、临界和超临界状态下带突起的圆柱绕流问题进行了数值模拟,分析了不同Reynolds数下带突起的圆柱绕流问题的近壁面流动特征以及表面时均压力与摩擦力系数的分布和凸起对圆柱表面流动分离以及转捩的影响,对比了有无凸起两侧圆柱表面时均压力、摩擦力系数的不同. 结果表明:当来流Reynolds数处于临界区时,气流在圆柱上表面凸起处形成了3个反向旋转的漩涡,之后随着θ的增大,发生了流动分离和流动转捩现象;对于不同Reynolds数下的来流,圆柱上表面的凸起可以使气流发生转捩的位置提前;圆柱上表面的凸起使流速增大、压强降低,从而导致圆柱产生升力,随着来流Reynolds数的增大,其升力逐渐变大.      

基于转捩模型及声比拟方法的高精度圆柱分离涡/涡致噪声模拟

基于七阶加权紧致非线性格式(WCNS-E8T7),结合延迟分离涡模拟(DDES)和Ffowcs WilliamsHawkings声比拟方法,对亚临界雷诺数下单圆柱、圆柱-翼型的分离涡/涡致噪声问题进行了数值模拟.针对亚临界雷诺数下圆柱尾迹中的转捩问题,发展了基于γ-Reθ模型高精度转捩-延迟分离涡模拟(Tran-DDES)方法,并与传统基于全湍流剪切应力输运(SST)模型的高精度DDES方法进行了对比.单圆柱模拟结果表明:传统SST-DDES方法会造成平均流场的回流区增大,压差阻力偏小等问题;而添加转捩模型的Tran-DDES方法与实验符合得很好.圆柱尾迹中添加翼型后,翼型对圆柱附近流场产生影响,使SST-DDES方法造成的圆柱后回流区偏大的问题减弱,并与Tran-DDES模拟结果差异变小.但在脉动量预测以及脉动产生的噪声预测方面, Tran-DDES方法仍与实验符合得更好.     

再入尾迹湍流对雷达散射截面影响分析

通过探讨高超声速再入湍流尾迹等离子体场中电磁波的散射机制,推导出在工程上描述湍流亚密等离子体雷达散射的一阶畸变波Born近似模型,分析了该模型在充分发展湍流尾迹等离子体场中的适用性,完成了适用于三维尾迹等离子体场的程序设计.以已有的湍流尾迹等离子体流场数据为基础,分析了再入尾迹湍流等离子体流动对雷达散射截面的影响.选取考察的几个有代表性的因素为:湍流模型、转捩过程、湍流尺度、电子组分脉动初始条件等.由结果可以看到,湍流转捩过程和湍流尺度对雷达散射截面值影响不大;电子组份脉动强度初始值影响较明显;湍流模型在特定条件下影响亦不大.     

电磁力控制湍流边界层分离圆柱绕流场特性数值分析

在亚临界区高雷诺数Re=1.4×105下,采用脱体涡模拟结合湍流分离的方法对弱电解质中电磁力作用下湍流边界层分离圆柱绕流场及其升(阻)力特性进行了数值模拟和分析.结果表明,电磁力可以提高圆柱体湍流边界层内的流体动能,延缓圆柱体湍流边界层的流动分离,减弱圆柱体湍流绕流场中在流向和展向上大尺度漩涡的强度,减小圆柱体阻力时均值及其升力脉动幅值.当电磁力作用参数大于某个临界值后,湍流边界层流动分离消失,在圆柱体尾部产生射流现象,从而电磁力对圆柱体产生净推力作用,出现负阻力现象,而且升力脉动幅值接近于零,出现圆柱体升力消失现象.