PIV analysis of near-wake behind a sphere at a subcritical Reynolds number

The vortical structure of near-wake behind a sphere is investigated using a PIV technique in a circulating water channel at Re = 11,000. The measured velocity fields show a detailed vortical structure in the recirculation region such as recirculation vortices, reversed velocity zone, and out-of-plane vorticity distribution. The vorticity distribution of the sphere wake shows waviness in cross-sectional planes. The time-averaged turbulent structures are consistent with the visualized flow showing the onset of shear layer instability. The spatial distributions of turbulent intensities provide turbulent statistics for validating numerical predictions.     

Numerical study of turbulent plane Couette flow at low Reynolds number

A direct numerical simulation of fully developed turbulent plane Couette flow has been performed. Unsteady large-scale structures, which contributed to the instantaneous energy level, were observed. These evolving and drifting vortical structures vanished after time-averaging, and the resulting mean velocity and streamwise turbulence intensity compared favourably with recent laboratory data.     

Parallel DNS with Local Grid Refinement

A parallel finite volume method for unstructured grids is used for a direct numerical simulation of the flow around a sphere at Re = 5000 (based on the sphere diameter and undisturbed velocity). The observed flow structures are confirmed by visualization experiments. A quantitative analysis of the Reynolds averaged flow provides a data base for future model evaluations.     

高超声速尾迹流场稳定性数值研究

通过数值模拟, 对高超声速尾迹流场进行了研究, 对其尾迹流动的失稳过程进行了分析.选取计算模型为圆球,Ma= 6.0, Re = 1.71\times 10^6(Re以球头半径为参考长度). 通过数值模拟,首先得到的流动是稳定解,在底部发展出一个主分离区和一个二次分离区,流动是轴对称状态. 不添加任何扰动继续进行计算,发现底部流场缓慢发展出微弱的非定常流动. 随后,该现象继续发展,出现明显的结构失稳,得到了无量纲周期为12.0的周期解. 给出了高超声速圆球绕流尾迹结构的周期性演化过程,对其涡系结构的演化及奇点特征进行了分析. 研究表明该数值模拟方法可用于底部流动稳定性问题的研究,同时证实了高超声速底部流动也存在流动不稳定性.      

Large-eddy simulation of counter-rotating Taylor–Couette flow: The effects of angular velocity and eccentricity

In the present work, turbulent flow in the annulus of a counter-rotating Taylor-Couette (CRTC) system is studied using large-eddy simulation. The numerical methodology employed is validated, for both the mean and second-order statistics, with the direct numerical simulation (DNS) data available in the literature, for a range of Reynolds numbers from 500 to 4000. Thereafter, turbulent flow occurring in this system at Reynolds numbers of 8000 and 16000 are studied, and the results obtained are analyzed using mean and second-order statistics, vortical structures, velocity vector plots and power energy spectra. Further, the spatio-temporal variation of azimuthal velocity, extracted near the inner cylinder, shows the existence of herringbone like patterns similar to that observed in the previous studies. The effect of eccentricity of the inner cylinder with respect to the outer cylinder is studied, on the turbulent flow in the CRTC system, for two different eccentricity ratios of 0.2 and 0.5 and for two different Reynolds numbers of 1500 and 4000. The results of the eccentric CRTC are analyzed using contours of pressure, mean and second-order statistics, velocity vectors, vortical structures, and turbulence anisotropy maps. It is observed from the eccentric CRTC simulations that the smaller-gap region seems to contain higher amplitude fluctuations and more vortical structures when compared with the larger-gap region. The mean turbulent kinetic energy contours do not change qualitatively with the Reynolds number, however, quantitatively a higher turbulent kinetic energy is observed in the higher Reynolds number case of 4000.     

Vortex dynamics of a sphere wake in proximity to a wall

Toward getting the vortex dynamics characteristics and wake structure of a sphere in proximity to a wall, the effect of a proximal flat plate on the wake of a stationary sphere is investigated via direct numerical simulation. The vortex shedding process and the significant variation of the wake structure are described in detail. The drag coefficient reduces and the wake structure of the sphere becomes complex due to the combined effect of the wake flow and the wall. A jet flow forms between the sphere and the flat plate, which suppresses the vortex separation on the bottom of the sphere. The asymmetric distributions of the coherent structures and the recirculation region behind the sphere are discussed. Besides vortex shedding patterns, the time-averaged velocity distribution, vortex dynamics, distribution regularities of turbulent kinetic energy and enstrophy are investigated.     

Application of the ultrasonic velocity profile method to the mapping of liquid metal flows under the influence of a non-uniform magnetic field

An ultrasonic velocity profile (UVP) method has been successfully applied to the investigation of a liquid metal channel flow under the influence of an inhomogeneous magnetic field. Using velocity profiles obtained by the ultrasonic velocimeter and their numerical post processing, two-dimensional time-averaged flow maps were efficiently produced. A single transducer immersed directly into the working fluid was used in order to simplify alignment of measurement lines and avoid the undesirable refraction of the acoustic beam on the walls. An M-shaped flow and wake behind a magnetic obstacle were reconstructed as the patterns of shear and large-scale vortical flows.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.     

Shear‐improved Smagorinsky modeling of turbulent channel flow using generalized Lattice Boltzmann equation

Generalized Lattice Boltzmann equation (GLBE) was used for computation of turbulent channel flow for which large eddy simulation (LES) was employed as a turbulence model. The subgrid‐scale turbulence effects were simulated through a shear‐improved Smagorinsky model (SISM), which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. The results of mean velocity distribution across the channel showed good correspondence with direct numerical simulation (DNS) data. Negligible discrepancies were observed between the present computations and those reported from DNS for the computed turbulent statistics. Three‐dimensional instantaneous vorticity contours showed complex vortical structures that appeared in such flow geometries. It was concluded that such a framework is capable of predicting accurate results for turbulent channel flow without adding significant complications and the computational cost to the standard Smagorinsky model. As this modeling was entirely local in space it was therefore adapted for parallelization. Copyright © 2010 John Wiley & Sons, Ltd.     

Large eddy simulation of separation control over a backward-facing step flow by suction

Separation control over a backward-facing step (BFS) flow by continuous suction was numerically investigated using the turbulence model of large eddy simulation (LES). The effect of suction control on the flow fields was scrutinised by altering the suction flow coefficient, and the results indicate that suction is not only very effective in shortening the reattachment length but also very influential in reducing the tangential velocity gradient and turbulence fluctuations of the reattached flows. With increasing increments of the absolute suction flow coefficient, the effect of suction control is more significant. Furthermore, the detailed flow fields (including the time-averaged stream and velocity fields) and turbulence characteristics (including the time-averaged resolved kinetic energy and RMS velocity) for the BFS models with or without suction are presented to discuss the mechanism of suction control. Comparisons of the time-averaged statistics between the numerical simulations and corresponding experiments are conducted, and it shows that the LES based on the dynamic kinetic energy subgrid-scale model (DKEM) can acquire exact results. Therefore, feasibility of the numerical methods to simulate suction-controlled models is validated.     

Three-dimensional numerical study of jet-in-crossflow characteristics at low Reynolds number

A numerical study of a square jet in a cross flow is carried out at a Reynolds number of 100. The flow field and heat transfer characteristic downstream of the jet have been explored by solving three-dimensional unsteady Navier–Stokes equations and energy equation using higher order spatial and temporal discretization. The projection of vortical structure on a plane is seen to give the component of vortex normal to the plane. Four combinations of velocity profile namely (1) uniform crossflow and uniform jet, (2) laminar boundary layer crossflow and uniform jet, (3) uniform crossflow and parabolic jet profile, and (4) laminar boundary layer crossflow and parabolic jet are compared at same phase to see their effect on the flow field and heat transfer characteristic. All the four cases are seen to exhibit unsteadiness but the jet with parabolic profile is seen to show stronger unsteadiness. The instantaneous vortical structures of all the cases at the same phase show that the structures are more complex for the jet with parabolic velocity profile. The temperature field is seen to be correlated with the vortical structures. Comparison of the time averaged flow field reveals that the jet penetration is the highest for the jet having parabolic profile and boundary layer crossflow. The adiabatic effectiveness is observed to be more for the jet with uniform velocity profile and uniform crossflow and was least for the jet with parabolic velocity profile and boundary layer crossflow.     

Vortex shedding and evolution induced by a solitary wave propagating over a submerged cylindrical structure

The shedding and evolution of the vortical structures generated by a solitary wave propagating over a submerged cylindrical structure are investigated experimentally and numerically. The cylindrical structure consists of two concentric cylinders and represents a simplified model for an offshore submerged intake structure typically used in coastal power plants. Flow visualization by dye injection is used to identify the dominant vortical structures near the structure. The flow visualization results show an unexpected flow reversal that causes shedding of secondary vortical structures. The experimental results are used to check a three-dimensional volume of fluid-large eddy simulation (VOF-LES) numerical model. The VOF-LES model is then used to further study the flow structure. A total of six dominant vortical structures generated by the wave motion are identified, followed by two more generated by the flow reversal. In summary, this paper provides the vorticity evolution for a complex fluid–structure interaction problem and a three-dimensional numerical simulation tool has also been validated, which can be extended to study more complex geometries and wave conditions.     

On the Validation of Large Eddy Simulation Applied to Internal Combustion Engine Flows Part II: Numerical Analysis

In internal combustion engines, the characteristic in-cylinder flow field is essential and significantly contributes to engine efficiency and performance. This paper describes the numerical investigation of the flow field in a motored 4-stroke, single-cylinder research engine. Quantitative and qualitative comparisons between experimental and numerical data have been performed at selected crank angle and results obtained in this work are discussed. Statistical flow properties are examined to analyze the averaged and instantaneous flow field. In order to investigate higher order statistical velocity moments and gain insight in the physical processes describing the engine flow structure, multi-cycle Large Eddy Simulation (LES) was carried out on two meshes with different spatial resolution. The three-dimensional structure of the flow has been also visualized by means of iso-surfaces of vortical structures, based on the Q criterion for individual cycles during intake. In order to assess the analysis and to verify that the computational mesh is applicable for the performance of LES simulations, the turbulence resolution M and the ratio of sgs-viscosity to the laminar viscosity were evaluated along the planes of interest. A direct comparison of the statistics of the flow field extracted from the numerical predictions shows a very good agreement with measurements conducted in the same configuration. Discrepancies have been however observed, in particular in the higher moments of the velocity components. Whilst this can be attributed mostly to the limited number of statistical sample (50 LES cycles) collected during the simulation, further investigation is certainly necessary to assess the relevance of modeling and spatial resolution issues.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Simulation of the Action of a Permeable Barrier on an Ideal Incompressible Channel Flow

Plane ideal incompressible flow in a rectangular channel partitioned by a thin permeable barrier (lattice) is considered. In flowing through the lattice the stream suddenly (jumpwise) changes direction and loses energy. The flow is assumed to be vortical; the vorticity is discontinuous on the lattice. A mathematical formulation of the problem for the stream function is proposed in the form of a nonlinear elliptic equation with coefficients discontinuous on the lattice line. A numerical solution is constructed using the finite-element iteration method. The results of the numerical simulation show how the flow velocity profile in the channel can be controlled by means of permeable barriers.     

VISCOUS FLOWFIELD BASED ON DISCONTINUOUS BOUNDARY ELEMENT METHOD AND VORTEX METHOD

基于非协调边界元方法和涡方法的联合应用, 模拟了二维和三维黏性不可压缩流场. 计算中利用离散涡元对漩涡的产生、凝聚和输送过程进行模拟, 并将整体计算域分解为采用涡泡模拟的内部区域和用涡列模拟的数字边界层区域. 计算域中涡量场的拉伸和对流由Lagrangian涡方法模拟, 用随机走步模拟涡量场的扩散. 内部区域涡元涡量场速度由广义Biot-Savart公式计算, 势流场速度则采用非协调边界元方法计算. 非协调边界元将所有节点均取在光滑边界处, 从而避免了法向速度的不连续现象; 而对于系数矩阵不对称的大型边界元方程组,引入了非常高效的预处理循环型广义极小残余(the generalized minimum residual, GMRES)迭代算法, 使得边界元法的优势得到了充分发挥, 同时, 在内部涡元势流场计算中对近边界点采用了正则化算法, 该算法将奇异积分转化为沿单元围道上一系列线积分, 消除了势流计算中速度及速度梯度的奇异性. 二维、三维流场算例证明了所用方法的正确性, 也验证了该算法可以大幅度提高模拟精度和效率.     

In-cylinder flows of a motored four-stroke engine with flat-crown and slightly concave-crown pistons

The temporal and spatial evolution processes of the in-cylinder flow structures and turbulence intensities in the symmetry and offset planes of a motored four-valve, four-stroke engine during the intake and compression strokes are diagnosed by using a particle image velocimeter. Two pistons of different crown shapes (flat-crown and slightly concave-crown pistons) are studied. The inception, establishment, and evolution of the tumbling vortical flow structures during the intake and compression strokes are clearly depicted. Quantitative strengths of the rotating vortical flow motions are presented by a dimensionless parameter, the tumble ratio, which can represent the mean angular velocity of the vortices in the target plane. The turbulence intensity of the in-cylinder flow is also calculated by using the measured time-varying velocity data. The results show that the flat-crown piston induces higher bulk-averaged tumble ratio and turbulence intensity than the slightly concave-crown piston does because the tumble ratio and turbulence generated by the flat-crown piston in the offset planes during the compression stroke are particularly large. The engine with the flat-crown piston also presents larger torque and power outputs and lower hydrocarbon emission than that with the slightly concave-crown piston. This might be caused by the enhanced combustion in the engine cylinder due to the stronger tumble ratio and turbulence intensity.     

三角翼大迎角不可压粘流的数值模拟

研究了人工压缩法拟压缩性系数β的选取，采用函数形式的β有效地加速了收敛过程．采用求解不可压Ｎ－Ｓ方程，对三角翼大迎角绕流进行了数值模拟，得到了与实验吻合很好的结果．分析和讨论了大迎角旋涡流动的复杂物理现象     

Unsteady separated and reattaching turbulent flow over a two-dimensional square rib

The spatio-temporal characteristics of the separated and reattaching turbulent flow over a two-dimensional square rib were studied experimentally. Synchronized measurements of wall-pressure fluctuations and velocity fluctuations were made using a microphone array and a split-fiber film, respectively. Profiles of time-averaged streamwise velocity and wall-pressure fluctuations showed that the shear layer separated from the leading edge of the rib sweeps past the rib and directly reattaches on the bottom wall (x/H=9.75) downstream of the rib. A thin region of reverse flow was formed above the rib. The shedding large-scale vortical structures (fH/U₀=0.03) and the flapping separation bubble (fH/U₀=0.0075) could be discerned in the wall-pressure spectra. A multi-resolution analysis based on the maximum overlap discrete wavelet transform (MODWT) was performed to extract the intermittent events associated with the shedding large-scale vortical structures and the flapping separation bubble. The convective dynamics of the large-scale vortical structures were analyzed in terms of the autocorrelation of the continuous wavelet-transformed wall pressure, cross-correlation of the wall-pressure fluctuations, and the cross-correlation between the wall pressure at the time-averaged reattachment point and the streamwise velocity field. The convection speeds of the large-scale vortical structures before and after the reattachment point were U_c=0.35U₀ and 0.45U₀, respectively. The flapping motion of the separation bubble was analyzed in terms of the conditionally averaged reverse-flow intermittency near the wall region. The instantaneous reattachment point in response to the flapping motion was obtained; these findings established that the reattachment zone was a 1.2H-long region centered at x/H=9.75. The reverse-flow intermittency in one period of the flapping motion demonstrated that the thin reverse flow above the rib is influenced by the flapping motion of the separation bubble behind the rib.