Assessment of Five SGS Models for Low-Rem MHD Turbulence

Assessment of three regularization-based and two eddy-viscosity-based subgrid-scale (SGS) turbulence models for large eddy simulations (LES) are carried out in the context of magnetohydrodynamic (MHD) decaying homogeneous turbulence (DHT) with a Taylor scale Reynolds number (Re_λ) of 120 and a MHD transition-to-turbulence Taylor-Green vortex (TGV) problems with a Reynolds number of 3000, through direct comparisons to direct numerical simulations (DNS). Simulations are conducted using the low-magnetic Reynolds number approximation (Re_m<<1). LES predictions using the regularization-based Leray- α,LANS- α, and Clark- α SGS models, along with the eddy viscosity-based non-dynamic Smagorinsky and the dynamic Smagorinsky models are compared to in-house DNS for DHT and previous results for TGV. With regard to the regularization models, this work represents their first application to MHD turbulence. Analyses of turbulent kinetic energy decay rates, energy spectra, and vorticity fields made between the varying magnetic field cases demonstrated that the regularization models performed poorly compared to the eddy-viscosity models for all MHD cases, but the comparisons improved with increase in magnitude of magnetic field, due to a decrease in the population of SGS eddies within the flow field.     

Kovasznay Mode Decomposition of Velocity-Temperature Correlation in Canonical Shock-Turbulence Interaction

The correlation coefficient R_uT between the streamwise velocity and temperature is investigated for the case of canonical shock-turbulence interaction, motivated by the fact that this correlation is an important component in compressible turbulence models. The variation of R_uT with the Mach number, the turbulent Mach number, and the Reynolds number is predicted using linear inviscid theory and compared to data from DNS. The contributions from the individual Kovasznay modes are quantified. At low Mach numbers, the peak post-shock R_uT is determined by the acoustic mode, which is correctly predicted by the linear theory. At high Mach numbers, it is determined primarily by the vorticity and entropy modes, which are strongly affected by nonlinear and viscous effects, and thus less well predicted by the linear theory.     

Aerodynamic characteristics of a thin airfoil cascade in an ideal incompressible flow with separation from the leading edges

A new method of solving the problem of separation flow past a cascade of thin airfoils is proposed. The method is based on the model of separated flow past a single airfoil in which vortex wakes shed from the airfoil edges are modeled by vortex layers with time-averaged intensities. Under the assumption of a small deviation of the freestream angle α from the “impactless entry” angle α ₀ the singular integral equation governing the flow kinematics is closed by a dynamic equation. It is shown that the separation effect on the time-averaged aerodynamic characteristics of the cascade is associated with the total pressure loss due to the flow energy expenditure on the vortex wake formation. The aerodynamic characteristics of the cascade calculated with and without account for flow separation differ by the second-order quantity ? = α ? α ₀.     

Numerical study on evolution of subharmonic varicose low-speed streaks in turbulent channel flow

The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing \( + \frac{{\partial u}}{{\partial x}}\). When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.     

Simulation of Receptivity and Induced Transition From Discrete Roughness Elements

Simulations have been carried out to predict the receptivity and growth of crossflow vortices created by Discrete Roughness Elements (DREs) The final transition to turbulence has also been examined, including the effect of DRE spacing and freestream turbulence. Measurements by Hunt and Saric (2011) of perturbation mode shape at various locations were used to validate the code in particular for the receptivity region. The WALE sub-grid stress (SGS) model was adopted for application to transitional flows, since it allows the SGS viscosity to vanish in laminar regions and in the innermost region of the boundary layer when transition begins. Simulations were carried out for two spanwise wavelengths: λ= 12mm (critical) and λ= 6mm (control) and for roughness heights (k) from 12 μm to 42 μm. The base flow considered was an ASU (67)-0315 aerofoil with 45 ⁰ sweep at -2.9 ⁰ incidence and with onset flow at a chord-based Reynolds number Re _c= 2.4x10 ⁶. For λ= 12mm results showed, in accord with the experimental data, that the disturbance amplitude growth rate was linear for k = 12 μm and 24 μm, but the growth rate was decreased for k = 36 μm Receptivity to λ= 6mm roughness showed equally good agreement with experiments, indicating that this mode disappeared after a short distance to be replaced by a critical wavelength mode. Analysis of the development of modal disturbance amplitudes with downstream distance showed regions of linear, non-linear, saturation, and secondary instability behaviour. Examination of breakdown to turbulence revealed two possible routes: the first was 2D-like transition (probably Tollmien-Schlichting waves even in the presence of crossflow vortices) when transition occurred beyond the pressure minimum; the second was a classical crossflow vortex secondary instability, leading to the formation of a turbulent wedge.     

Detached eddy simulation of flow past an isolated inclined solar panel

Detached Eddy Simulations (DES) were carried out to analyze the influence of ground clearance and to estimate the unsteady wind loading on solar panels. A generic type of solar panel was considered at three different ground clearances. For the range of ground clearance values considered in the study, the results show that an increase in the clearance produces stronger vortex shedding fluctuations, higher shedding frequencies, and larger mean and unsteady wind loading. For the smallest ground clearance, an irregular vortex shedding was observed. The flow structures formed in the wake are analyzed using various indicators of turbulence characteristics including the vorticity and the λ₂ criterion. The flow structures suggest a complex formation and subsequent shedding of vortices. The absence of regular vortex shedding at the smallest ground clearance and the presence of the same at the mid-ground clearance value suggests a reference for minimizing the unsteady wind loading of solar panels. Furthermore, the unsteady results reveal that the lower sections of the panel are more vulnerable to wind loading.     

Mach waves produced in the supersonic jet mixing layer by shock/vortex interaction

The noise emission of free jets has been extensively investigated for many decades. At subsonic jet velocities, coherent structures of the mixing layer move at subsonic speed and emit sound waves. Free jets blowing at supersonic speeds, however, can emit weak shock waves, called Mach waves. At supersonic speeds, two cases must be distinguished: the structures move either subsonically or supersonically relative to the inside and/or outside speed of sound. In the case of supersonic movement, the Mach waves exist inside as well as outside the jet. At subsonic speeds, no Mach waves appear. Although numerous theories have been established to find the origin of the Mach waves, to the authors’ best knowledge, the mechanism of the Mach wave formation has not yet been clearly explained. Recently another theory of Mach waves in supersonic jets was developed, as described herein, which outlines the causes for the Mach wave production and stability as well as their dynamics. The theory’s principle is that the Mach waves are initiated by vortices which move downstream at three speeds w, \({w}'\) and \({w}''\) inside of the mixing layer. These three types of vortices and Mach waves are described in a comprehensive manner by the theory and are called the “w-, \({w}'\)- and \({w}''\)-vortices” and “w-, \({w}'\)- and \({w}''\)-Mach waves,” respectively.     

Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Study of the Spectral Difference Numerical Dissipation for Turbulent Flows Using Unstructured Grids

In this paper, the numerical dissipation properties of the Spectral Difference (SD) method are studied in the context of vortex dominated flows and wall-bounded turbulence, using uniform and distorted grids. First, the validity of using the SD numerical dissipation as the only source of subgrid dissipation (the so-called Implicit-LES approach) is assessed on regular grids using various polynomial degrees (namely, p = 3, p = 4, p = 5) for the Taylor-Green vortex flow configuration at R e = 5 000. It is shown that the levels of numerical dissipation greatly depend on the order of accuracy chosen and, in turn, lead to an incorrect estimation of the viscous dissipation levels. The influence of grid distortion on the numerical dissipation is then assessed in the context of finite Reynolds number freely-decaying and wall-bounded turbulence. Tests involving different amplitudes of distortion show that highly skewed grids lead to the presence of small-scale, noisy structures, emphasizing the need of explicit subgrid modeling or regularization procedures when considering coarse, high-order SD computations on unstructured grids. Under-resolved, high-order computations of the turbulent channel flow at R e _τ = 1000 using highly-skewed grids are considered as well and present a qualitatively similar agreement to results obtained on a regular grid.     

Dual yielding in capillary suspensions

Rheological measurements were performed to examine the yielding behavior of capillary suspensions prepared by mixing cocoa powder as dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid. Here, we investigated the yielding behavior of solid-fluid-fluid systems with varying particle volume fraction, ?, spanning the regime from a low volume fraction (? = 0.25) to a highly filled regime (? = 0.65) using dynamic oscillatory measurements. While for ? ≤ 0.4 with a fixed water volume fraction (? _w ) of 0.06 as the secondary fluid, capillary suspensions exhibited a single yield point due to rupturing of aqueous capillary bridges between the particles, while capillary suspensions with ? ≥ 0.45 showed a two-step yielding behavior. On plotting elastic stress (G^′ γ) as a function of applied strain (γ), two distinct peaks, indicating two yield stresses, were observed. Both the yield stresses and storage modulus at low strains were found to increase with ? following a power law dependence. With increasing ? _w (0 – 0.08) at a fixed ? = 0.65, the system shifted to a frustrated, jammed state with particles strongly held together shown by rapidly increasing first and second yield stresses. In particular, the first yield stress was found to increase with ? _w following a power law dependence, while the second yield stress was found to increase exponentially with ? _w . Transient steady shear tests were also performed. The single stress overshoot for ? ≤ 0.4 with ? _w = 0.06 reflected one-step yielding behavior. In contrast, for high ? (≥ 0.45) values with ? _w = 0.06, two stress overshoots were observed in agreement with the two-step yielding behavior shown in the dynamic oscillatory measurements. Experiments on the effect of resting time on microstructure recovery demonstrated that aggregates could reform after resting under quiescent conditions.     

Experimental study of the velocity field in a vortex-flow heat exchanger with isothermal conditions

The velocity field in a vortex heat cell was investigated experimentally using laser Doppler velocimetry for a wide range of flow conditions. Experimental results point out the three dimensionality of the exchanger's flow, which is composed into a main vortex flow developing along the side walls. The strength of the flow increases up to a limiting value reached for a Reynolds number ranging between 15,000 and 30,000; a secondary flow, caused by interaction between centrifugal and inertial forces, extends perpendicularly to the main flow and remains Reynolds number dependent. It is composed of multiple counter-rotating structures occurring at the exchanger periphery with low inlet Reynolds numbers, thus reducing the rate of centripetal momentum transfer. With increasing inlet Reynolds number, the secondary flow extends across the whole exchanger radius, thus increasing the rate of mixing of the treated fluid. The appearance of so-called Taylor–Görtler vortices tends to reduce the z- and r-axis vorticity transfer.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Manipulation of Large-Scale Vortical Structures and Mixing in a Coaxial Jet with Miniature Jet Actuators Toward Active Combustion Control

Manipulation of large-scale vortical structures and associated mixing in a methane-air coaxial jet is carried out by using miniature jet actuators installed on the inner surface of the annular nozzle. The periodic radial miniature jet injections are achieved with a rapid-response servo-valve. The spatio-temporal primary jet structures are investigated through phase-locked 2C-PIV (2 Component Particle Image Velocimetry) and stereoscopic-PIV. It is found that intense ring-like vortices are produced perfectly in phase with the periodic miniature jet injections regardless of the valve-driven frequency f_v examined. When the Strouhal number St_v, which is defined with f_v, is larger than unity, the ring-like vortices are densely formed and thus methane/air mixing is prompted with low periodic fluctuation. The diameter of the vortices becomes small as St_v is increased, so that the transport range of the inner methane and outer air fluids can be controlled by changing St_v. In addition, the evolution of counter-rotating vortex pair is also observed in the cross-sectional plane. These streamwise vortices are directly formed as a result of the radial miniature jet injection, which leads to entrainment of the ambient fluid near the primary jet shear layer, and they also contribute to the mixing enhancement. Moreover, it is demonstrated that coaxial jet flame characteristics such as carbon monoxide (CO) emission and flame holding can be drastically improved under different equivalence ratios by the present jet control scheme.     

Effect of Flow Structures on Turbulence Statistics of Taylor-Couette Flow in the Torque Transition State

Direct numerical simulations of Taylor-Couette flow from Re= 8000 to 25000 have been conducted to investigate changes of turbulence statistics in the transition of the Reynolds number dependency of the mean torque near Re= 10000. The velocity fluctuations are decomposed into the contributions of the Taylor vortex and remaining turbulent fluctuations. Significant Reynolds number dependencies of these components are observed in the radial profiles of the Reynolds stress and the transmission of the mean torque. The contributions of Taylor vortex and turbulent components in the net amount of mean torque are evaluated. The Taylor vortex component is overtaken by the turbulent counterpart around Re= 15000 when they are defined as the azimuthally averaged component and the remnants. The results show that the torque transition can be explained by the competition between the contributions of azimuthally averaged Taylor vortex and the remaining turbulent components.     

Near-wake flow structure downwind of a wind turbine in a turbulent boundary layer

Wind turbines operate in the surface layer of the atmospheric boundary layer, where they are subjected to strong wind shear and relatively high turbulence levels. These incoming boundary layer flow characteristics are expected to affect the structure of wind turbine wakes. The near-wake region is characterized by a complex coupled vortex system (including helicoidal tip vortices), unsteadiness and strong turbulence heterogeneity. Limited information about the spatial distribution of turbulence in the near wake, the vortex behavior and their influence on the downwind development of the far wake hinders our capability to predict wind turbine power production and fatigue loads in wind farms. This calls for a better understanding of the spatial distribution of the 3D flow and coherent turbulence structures in the near wake. Systematic wind-tunnel experiments were designed and carried out to characterize the structure of the near-wake flow downwind of a model wind turbine placed in a neutral boundary layer flow. A horizontal-axis, three-blade wind turbine model, with a rotor diameter of 13 cm and the hub height at 10.5 cm, occupied the lowest one-third of the boundary layer. High-resolution particle image velocimetry (PIV) was used to measure velocities in multiple vertical stream-wise planes (x–z) and vertical span-wise planes (y–z). In particular, we identified localized regions of strong vorticity and swirling strength, which are the signature of helicoidal tip vortices. These vortices are most pronounced at the top-tip level and persist up to a distance of two to three rotor diameters downwind. The measurements also reveal strong flow rotation and a highly non-axisymmetric distribution of the mean flow and turbulence structure in the near wake. The results provide new insight into the physical mechanisms that govern the development of the near wake of a wind turbine immersed in a neutral boundary layer. They also serve as important data for the development and validation of numerical models.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

FEM solution to natural convection flow of a micropolar nanofluid in the presence of a magnetic field

The two-dimensional, laminar, unsteady natural convection flow in a square enclosure filled with aluminum oxide (\(\hbox {Al}_{2} \hbox {O}_{3}\))–water nanofluid under the influence of a magnetic field, is considered numerically. The nanofluid is considered as Newtonian and incompressible, the nanoparticles and water are assumed to be in thermal equilibrium. The mathematical modelling results in a coupled nonlinear system of partial differential equations. The equations are solved using finite element method (FEM) in space, whereas, the implicit backward difference scheme is used in time direction. The results are obtained for Rayleigh (Ra), Hartmann (Ha) numbers, and nanoparticles volume fractions (\(\phi\)), in the ranges of \(10^3 \le Ra \le 10^7\), \(0\le Ha \le 500\) and \(0 \le \phi \le 0.2\), respectively. The streamlines and microrotation contours are observed to show similar behaviors with altering magnitudes. For low Ra values, when \(Ha=0\), symmetric vortices near the walls and a central vortex in opposite direction are observed in vorticity. As Ra increases, the central vortex splits into two due to the circulation in the effect of the buoyant flow. Boundary layer formation is observed when Ha increases for almost all Rayleigh numbers in both streamlines and vorticity. The isotherms have horizontal profiles for high Ra values owing to convective dominance over conduction. As Ha is increased, the convection effect is reduced, and isotherms tend to have vertical profiles. This study presents the first FEM application for solving highly nonlinear PDEs defining micropolar nanofluid flow especially for large values of Rayleigh and Hartmann numbers.     

DNS of Turbulent Flows in Ducts with Complex Shape

We carry out Direct Numerical Simulation (DNS) of flows in closed straight ducts with complex peripheral shape. To perform the simulations the Navier-Stokes equations in cylindrical coordinates are discretized by a second-order finite difference scheme, and the immersed-boundary technique is used to resolve the flow close to walls of complex shape. The basic geometry is a circular pipe of radius R, with imposed sinusoidal perturbations of the type \(\eta R \sin (N_{w}\theta )\). Simulations by varying N _w at fixed η were performed to investigate the effect of the perturbation wavenumber. Additional simulations by fixing N _w and varying η also allow to investigate the influence of the amplitude of the wall corrugations. The modifications of the near-wall structures due to change in the shape of the walls are well depicted through contour plots of the radial component of the vorticity. The presence of geometrical disturbances anchors the structures at the locations where curvature changes, and the shape of the structures is strongly linked to the amplitude of the wall corrugation. Our interest is also in understanding the influence of the shape of the surface on wall friction. We were expecting some changes in the profile of the total stress with respect to that of the circular pipe, which instead were not found. This is a first indication that changes in the near-wall region do not affect the outer region, and that Townsend’s similarity hypothesis holds.     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.