Calculations of the flow resistance and heat emission of a sphere in the laminar and high-turbulent gas flows

An early drag crisis can occur at high turbulence of incoming gas flow to a sphere. To study the influence of a crisis on heat transfer from a sphere to gas, a numerical experiment was carried out in which the free gas flow around a sphere with a temperature lower than the sphere temperature was simulated for two cases. The flow was laminar in the first case and highly turbulent in the second case. To take into account turbulence, the kinematic coefficient of turbulent viscosity with a value, which is much higher (up to 2000 times) than that for physical viscosity, was introduced. The results of calculations show that the early drag crisis occurs at Reynolds numbers of about 100 and results in considerable (by four to seven times) decrease in the hydrodynamic force and sphere drag coefficient C _d . The early drag crisis is also accompanied by the crisis of heat transfer from a sphere to gas with a decrease in Nusselt numbers Nu by three to six times.     

Calculation of the flow about a sphere and the drag of the sphere under laminar and strongly turbulent conditions

To investigate the influence of a strongly turbulent incoming flow on the hydrodynamic drag of a body and occurrence of the early crisis of drag, a numerical experiment is conducted in which a free gas flow about a sphere is simulated for two cases, namely, for a laminar flow and for a strongly turbulent flow. Turbulence is simulated by assuming a high kinematic coefficient of turbulent viscosity. Calculation data lead us to conclude that the early crisis of drag at Reynolds numbers near 100, which shows up as a considerable (four-to sevenfold) decrease in the hydrodynamic force and the drag coefficient of the body, can be explained by the strong turbulence of the incoming flow.     

Translational oscillations of a microsphere in superfluid helium

The translational motion of a microsphere (radius 100 μm) in liquid helium is investigated. The sphere is levitating inside a superconducting capacitor and oscillates about its equilibrium position. The velocity amplitude and the resonance frequency are measured as a function of driving force and temperature (0.35 K up to 2.2 K). By increasing the driving force we first find a linear regime (laminar flow) which changes abruptly into a nonlinear one (turbulent flow). For temperatures below 0.7 K the linear drag is given by ballistic roton and phonon scattering whereas for temperatures above 1.1 K the hydrodynamic force on the sphere is described by Stoke's solution. In the turbulent regime, above a temperature independent threshold velocity, we find the drag force to be given by turbulence in the superfluid component plus an essentially laminar drag by the normal component.     

Calculation of the drag and heat transfer from a sphere in the gas flow in a cylindrical channel

A numerical experiment on the simulation of heat transfer from a sphere to a gas flow in a cylindrical channel in the Stokes and transient flow regimes has been described. Radial and axial profiles of the gas temperature and the dependences of drag coefficient C_d of the body and Nusselt number Nu on Reynolds number Re have been calculated and analyzed. The problem of the influence of the early drag crisis for a sphere on its heat transfer to the gas flow has been considered. The estimation of this phenomenon has shown that the early drag crisis of the sphere in a strongly turbulent flow causes a reduction in heat transfer from the sphere to the gas by three to six times (in approximately the same proportion as for its drag coefficient).     

Crisis of hydrodynamic drag of drops in the two-phase turbulent flow of a spray produced by a mechanical nozzle at transition reynolds Numbers

When processing experimental data for the hydrodynamics of a two-phase flow in a spray produced by a mechanical nozzle, we revealed an anomaly in the behavior of the hydrodynamic drag of drops: the drag coefficient turns out to be four to seven times lower than the previously known values. Several hypotheses are put forward to explain the anomaly. It is found that, when the gas flows around drops under highly turbulent conditions, an “early” (i.e., observed even at transition Reynolds numbers, Re>50) crisis of drag resistance of drops takes place. This new physical phenomenon allows us to account for a number of features of the two-phase flow that are observed in the experiment. Among these features is, in particular, the fact that the momentum transferred to the gas is roughly half the initial momentum of the liquid jet.     

The effect of riblets in rectangular duct flow

Much is known about the benefits of surface structures which mimic the riblets found on the skin of fast swimming sharks. Structures have been fabricated for study and application which replicate and improve upon the natural shape of the shark skin riblets, providing a maximum drag reduction of nearly 10% in external turbulent fluid flow. Mechanisms of fluid drag in turbulent flow and riblet-drag reduction theories from experiment and simulation are discussed. A brief review of riblet performance studies is given, and optimal riblet geometries are defined for external flow. A discussion of the structure of internal turbulent fluid flow is provided, and existing data for riblet drag reduction benefit is presented. A flow cell for studying riblet effects in internal rectangular duct flow is discussed, and data collected using several riblet surfaces fabricated for the flow cell is presented and analyzed. A discussion of the effects of the riblets on fluid flow is given, and conclusions are drawn about the possible benefits of riblets in internal fluid flow.     

Effect of groove shape on flow characteristics around a circular cylinder with grooves

In the flow around a circular cylinder, a sudden decrease in the drag force occurs at a high Reynolds number, but the same phenomenon occurs at a lower Reynolds number in the case where there exist grooves or roughness on the circular cylinder surface. In this paper, in order to make clear the flow characteristics around a circular cylinder in the case of changing the shapes of grooves, the drag coefficient, pressure distribution, velocity distribution and turbulent distribution were measured. Moreover the flow around the cylinder was analyzed by applying the RNGk · ∈ turbulent model, and the surface flow pattern was investigated using the oil-film technique. From these results, it is clear that the drag coefficient of a circular cylinder with triangular grooves decreases by about 15% compared with that of a circular cylinder with arc grooves.     

Viscous Interaction Between a Vortex Tube and a Rotating Spherical Particle

The transient advection of a cylindrical vortex tube in a viscous incompressible flow field and its interaction with a rotating/spinning spherical particle has been investigated numerically at Reynolds numbers in the range of 20≤ Re≤200 for angular velocities of 0≤Ω≤0.5. The effects of vortex parameters such as size, circulation strength and initial position relative to the particle, on the temporal behavior of the lift and drag forces are studied. Vortex‐sphere interactions bring about major changes in the flow field particularly when coupled with particle rotation. It is observed that the forces acting on the particle are significantly influenced during the time that the vortex core is in the vicinity of the particle. The extent of these local changes are about ±30% in the drag coefficient and about ±200% in the lift coefficient as compared to flow over a rotating solid sphere with no vortex interaction. It is also found that a vortex with core radius between one and two particle diameters creates the strongest temporal variations in the lift and drag coefficients. Furthermore, maximum lift variations occur for the vortex‐particle head on collision, while a vortex with an offset distance of about one diameter from the principal flow axis generates the maximum drag variations.     

The flow of a model suspension fluid past a sphere

The paper is concerned with the flow of a dilute suspension of monosized spheroids past a sphere placed at the centerline of a cylindrical tube. The suspension is modeled by the transversely isotropic fluid model and the numerical solution is obtained by a time-dependent boundary element method. No steady-state solution to the problem was found. However, the amplitude of the long-time oscillations in the drag force on the sphere is only about 1% of its mean value at a solid volume concentration of 0.01 and an aspect ratio of 10 for the spheroids. The initial orientation of the microstructure has a small effect on the drag force: it can give rise to a ±1 % variation in the drag force at large time.     

Variation of the temperature factor of models on a ballistic range

A technique for variation of the temperature factor of free-flight models by varying their initial temperature is described. An experiment on a ballistic range is carried out with a free-flying supersonic blunt cone with a half-angle of 15° at a Mach number of 2.3. The flow at the cone base is studied in the transition range (from the laminar to turbulent flow) of Reynolds numbers. The base flow pattern is determined from the shadowgraphs of the flow about the models. The drag coefficient of the blunt cone at a zero angle of attack is found by processing trajectory data. It is found that the near wake geometries and the drag coefficients of the models tested at the laboratory temperature and a temperature of 120 K differ. Explanations of this effect are given.     

介观尺度流体绕流球体的耗散粒子动力学模拟

采用耗散粒子动力学(dissipative particle dynamics, DPD)方法, 对两平行平板间流体绕流三维球体进行了计算. 球体和平行平板由达到平衡状态的冻结DPD粒子组成, 流体在不同无量纲外力驱动下流动, 球体受力由组成球体的所有冻结DPD粒子求和得到. 流动达到充分发展后, 输出球体在流动方向的受力, 并计算球体的阻力系数, 与文献中的关联式进行了对比. 结果表明, 在Re≤qslant 100的范围内, DPD方法能较准确地计算出阻力系数, 在较大雷诺数时, 由于流     

Drag reduction by herringbone riblet texture in direct numerical simulations of turbulent channel flow

A bird-feather-inspired herringbone riblet texture was investigated for turbulent drag reduction. The texture consists of blade riblets in a converging/diverging or herringbone pattern with spanwise wavelength Λ_f. The aim is to quantify the drag change for this texture as compared to a smooth wall and to study the underlying mechanisms. To that purpose, direct numerical simulations of turbulent flow in a channel with height L_z were performed. The Fukagata-Iwamoto-Kasagi identity for drag decomposition was extended to textured walls and was used to study the drag change mechanisms. For Λ_f/L_z ? O(10), the herringbone texture behaves similarly to a conventional parallel-riblet texture in yaw: the suppression of turbulent advective transport results in a slight drag reduction of 2%. For Λ_f/L_z ? O(1), the drag increases strongly with a maximum of 73%. This is attributed to enhanced mean and turbulent advection, which results from the strong secondary flow that forms over regions of riblet convergence/divergence. Hence, the employment of convergent/divergent riblets in the texture seems to be detrimental to turbulent drag reduction.     

Impact of trailing wake drag on the statistical properties and dynamics of finite-sized particle in turbulence

We study by means of an Eulerian-Lagrangian model the statistical properties of velocity and acceleration of a neutrally-buoyant finite-sized particle in a turbulent flow statistically homogeneous and isotropic. The particle equation of motion, besides added mass and steady Stokes drag, keeps into account the unsteady Stokes drag force-known as Basset-Boussinesq history force-and the non-Stokesian drag based on Schiller-Naumann parametrization, together with the finite-size Faxén corrections. We focus on the case of flow at low Taylor-Reynolds number, Re_λ?31, for which fully resolved numerical data which can be taken as a reference are available [Homann H., Bec J. Finite-size effects in the dynamics of neutrally buoyant particles in turbulent flow. J Fluid Mech 651 (2010) 81-91]. Remarkably, we show that while drag forces have always minor effects on the acceleration statistics, their role is important on the velocity behavior. We propose also that the scaling relations for the particle velocity variance as a function of its size, which have been first detected in fully resolved simulations, does not originate from inertial-scale properties of the background turbulent flow but it is likely to arise from the non-Stokesian component of the drag produced by the wake behind the particle. Furthermore, by means of comparison with fully resolved simulations, we show that the Faxén correction to the added mass has a dominant role in the particle acceleration statistics even for particles whose size attains the integral scale.     

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

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

Bubbly turbulent drag reduction is a boundary layer effect

In turbulent Taylor-Couette flow, the injection of bubbles reduces the overall drag. On the other hand, rough walls enhance the overall drag. In this work, we inject bubbles into turbulent Taylor-Couette flow with rough walls (with a Reynolds number up to 4 x 10(5), finding an enhancement of the dimensionless drag as compared to the case without bubbles. The dimensional drag is unchanged. As in the rough-wall case no smooth boundary layers can develop, the results demonstrate that bubbly drag reduction is a pure boundary layer effect.     

Effectiveness of blowing for improving the high-speed trains aerodynamics

The promising method of drag reduction with the use of micro-blowing through the streamlined surface has been proposed for its use to the external surface of high-speed train. The advantages of high-speed train as an object of micro-blowing application are introduced. The corresponding RANS-based mathematical model is elaborated, and the computations of the external flow around a long train body are performed. Predictions of the turbulent boundary layer over penetrable surface with different modes of micro-blowing have been presented and analyzed. The developed modifications of mathematical model of turbulence have been used to take into account the micro-blowing influence in the inner region of turbulent boundary layer. The obtained results of parametric analysis of drag reduction depending on the area of permeable sections, intensity of micro-blowing, and high-speed train length have been analyzed. In particular, the dependence between drag reduction effect and length of train body with realized micro-blowing as well as its intensity is established. Realization of micro-blowing with blowing velocity just 0.25 % of train speed (V = 100 m/s) on the 70 % of the streamlined surface area for just one train carriage (L = 25 m) allows one to reduce the aerodynamic drag (including the most actual friction and head-tail pressure components) of the whole train (L = 200 m) by about 5.25 %, so in case of micro-blowing realization on all its 8 carriages, the train’s aerodynamic drag can be reduced approximately by 42 %.     

Use of DES in mildly separated internal flow: dimples in a turbulent channel

Detached eddy simulation (DES) is investigated as a means to study an array of shallow dimples with depth to diameter ratios of 1.5% and 5% in a turbulent channel. The DES captures large-scale flow features relatively well, but is unable to predict skin friction accurately due to flow modelling near the wall. The current work instead relies on the accuracy of DES to predict large-scale flow features, as well as its well-documented reliability in predicting flow separation regions to support the proposed mechanism that dimples reduce drag by introducing spanwise flow components near the wall through the addition of streamwise vorticity. Profiles of the turbulent energy budget show the stabilising effect of the dimples on the flow. The presence of flow separation however modulates the net drag reduction. Increasing the Reynolds number can reduce the size of the separated region and experiments show that this increases the overall drag reduction.     

确定分布的展向Lorentz力调制下的槽道湍流涡结构

采用直接数值模拟方法,对槽道湍流中确定分布的Lorentz力的流动控制与减阻问题进行研究.讨论了Lorentz力作用于槽道湍流后,流场的特性和涡结构的特性,并对此类Lorentz力对槽道湍流的控制与减阻机理进行了讨论.研究发现:1)Lorentz力诱导的层流流场壁面附近存在梯度极大的展向速度剪切层,该剪切层容易形成流向涡结构;2)在给定合适参数的确定分布的Lorentz力作用下,湍流流场仅剩周期分布的准流向涡;3)与未控制流场相比,控制后的流场中,准流向涡的抬升高度大大降低,从而减小猝发强度,使壁面阻力下降.     

Flow characteristics around a square cylinder with changing chamfer dimensions

Abstract  

It is known that for a square cylinder subjected to uniform flow, the drag force changes with the angle of attack. To clarify the flow characteristics around a square cylinder with corner cutoffs, we measured the drag coefficient and the Strouhal number for changing chamfer dimensions. We analyzed the flow around a square cylinder with corner cutoffs by applying the RNG k–ε turbulent model, and investigated the surface flow pattern using visualization by means of the oil film and mist flow method. From these results, we obtained the surface flow patterns by the oil film method and numerical analysis. The numerical results agreed well with the experimental values. The drag coefficient of the square cylinder with corner cutoffs decreased suddenly at an angle of attack of about α = 0°– 10° when compared with the drag coefficient for a square cylinder. The minimum value of the drag coefficient for the square cylinder with corner cutoffs decreased by about 30% compared with that for the square cylinder. The drag coefficient of the square cylinder with 10% corner cutoffs was found to be smallest, since the wake area of this square cylinder was smaller compared with that of the other square cylinder.     

介观尺度流体绕流球体的耗散粒子动力学模拟

采用耗散粒子动力学(dissipative particle dynamics, DPD)方法, 对两平行平板间流体绕流三维球体进行了计算. 球体和平行平板由达到平衡状态的冻结DPD粒子组成, 流体在不同无量纲外力驱动下流动, 球体受力由组成球体的所有冻结DPD粒子求和得到. 流动达到充分发展后, 输出球体在流动方向的受力, 并计算球体的阻力系数, 与文献中的关联式进行了对比. 结果表明, 在Re≤qslant 100的范围内, DPD方法能较准确地计算出阻力系数, 在较大雷诺数时, 由于流体的压缩性导致计算结果出现差异.