Flow instability of nanofuilds in jet

The flow instability of nanofluids in a jet is studied numerically under various shape factors of the velocity profile, Reynolds numbers, nanoparticle mass loadings,Knudsen numbers, and Stokes numbers. The numerical results are compared with the available theoretical results for validation. The results show that the presence of nanoparticles enhances the flow stability, and there exists a critical particle mass loading beyond which the flow is stable. As the shape factor of the velocity profile and the Reynolds number increase, the flow becomes more unstable. However, the flow becomes more stable with the increase of the particle mass loading. The wavenumber corresponding to the maximum of wave amplification becomes large with the increase of the shape factor of the velocity profile, and with the decrease of the particle mass loading and the Reynolds number. The variations of wave amplification with the Stokes number and the Knudsen number are not monotonic increasing or decreasing, and there exists a critical Stokes number and a Knudsen number with which the flow is relatively stable and most unstable,respectively, when other parameters remain unchanged. The perturbation with the first azimuthal mode makes the flow unstable more easily than that with the axisymmetric azimuthal mode. The wavenumbers corresponding to the maximum of wave amplification are more concentrated for the perturbation with the axisymmetric azimuthal mode.     

Stability of viscoelastic shear flows in the limit of high Weissenberg and Reynolds numbers

We consider a limit of the upper convected Maxwell model where both the Weissenberg and Reynolds numbers are large. The limiting equations have a status analogous to that of the Euler equations for the high Reynolds number limit. These equations admit parallel shear flows with an arbitrary profile of velocity and normal stress. We consider the stability of these flows. An extension of Howard’s semicircle theorem can be used to show that the flow is stabilized if elastic effects are sufficiently strong. We also show how to analyze the long wave limit in a fashion similar to the inviscid case.     

Large eddy simulation using a dynamic mixing length subgrid‐scale model

A novel dynamic mixing length (DML) subgrid‐scale model for large eddy simulations is proposed in this work to improve the cutoff length of the Smagorinsky model. The characteristic mixing length (or the characteristic wave number) is dynamically estimated for the subgrid‐scale fluctuation of turbulence by the cutoff wave‐number, k_c, and the dissipation wave‐number, k_d. The dissipation wave number is derived from the kinetic energy spectrum equation and the dissipation spectrum equation. To prove the promise of the DML model, this model is used to simulate the lid‐driven cubical cavity with max‐velocity‐based Reynolds numbers 8850 and 12,000, the channel flows with friction‐velocity‐based Reynolds numbers 180, 395, 590, and 950, and the turbulent flow past a square cylinder at the higher Reynolds number 21,400, respectively, compared with the Smagorinsky model and Germano et al.'s dynamic Smagorinsky model. Different numerical experiments with different Reynolds numbers show that the DML model can be used in simulations of flows with a wide range of Reynolds numbers without the occurrence of singular values. The DML model can alleviate the dissipation of the Smagorinsky model without the loss of its robustness. The DML model shows some advantages over Germano et al.'s dynamic Smagorinsky model in its high stability and simplicity of calculation because the coefficient of the DML model always stays positive. The characteristic mixing length in the DML model reflects the subgrid‐scale fluctuation of turbulence in nature and thus the characteristic mixing length has a spatial and temporal distribution in turbulent flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation and visualization of fiber suspension in a turbulent round jet

The turbulent properties of the fiber suspension in a turbulent round jet are numerically simulated and visualized, and some of the results are compared with the experimental data. The effects of the Reynolds number, fiber volume fraction, and aspect ratio are analyzed. The results show that the fiber injection in the flow has a delay effect on the streamwise velocity decay along the jet axis, and such an effect becomes more obvious with the increases in the fiber volume fraction and aspect ratio and the decrease in the Reynolds number. The flow with fibers shows an increase in the streamwise velocity along the radial direction, and the increase magnitude is directly proportional to the fiber volume fraction and aspect ratio and inversely proportional to the Reynolds number. The presence of fibers makes the turbulent kinetic energy and Reynolds stress increase, and the extent increases with the fiber volume fraction, Reynolds number, and fiber aspect ratio.     

Instability of plane Couette flow of viscoelastic liquids

A review of our work on the stability of plane Couette flow of a viscoelastic liquid is given. The first part of the review is based on the assumption of a “short memory” of the fluid. The Reynolds-Orr energy criterion intimates the possibility of instability at very low Reynolds numbers. A linear stability analysis for disturbances in the flow plane shows that beyond the stability limit given by the energy criterion there are always disturbances which grow with time. A critical assessment of the short memory theory shows the severe limitations of its applicability.In the second part of the paper, the assumption of short memory is dropped. The stability of plane Couette flow with respect to special disturbances perpendicular to the flow plane is investigated for a Maxwell fluid. The flow is unstable if the product of Reynolds number and Weissenberg number is higher than a certain limit, which has the value one for a simple Maxwell fluid. This result can also be interpreted as follows: The flow becomes unstable if the velocity at the boundary walls is higher than the shear wave velocity of the fluid.     

Optimal transient growth in turbulent pipe flow

The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the governing equations onto the subspace spanned by the radial vorticity and radial velocity. The method is validated by comparing with the previous studies. Two peaks of the maximum transient growth amplification curve are found at different Reynolds numbers ranging from 20 000 to 250 000. The optimal flow structures are obtained and compared with the experiments and DNS results. The location of the outer peak is at the azimuthal wave number n=1, while the location of the inner peak is varying with the Reynolds number. It is observed that the velocity streaks in the buffer layer with a spacing of 100δ_v are the most amplified flow structures. Finally, we consider the optimal transient growth time and its dependence on the azimuthal wave length. It shows a self-similar behavior for perturbations of different scales in the optimal transient growth process.     

On Rayleight instability in decaying plane Couette flow

The inflexion point criterion of Rayleigh is one of the most well-known results in hydrodynamic stability theory but cannot easily be demonstrated experimentally in wall bounded flows. For plane Couette flow, where both walls move with equal speed in opposite directions, it is possible to establish a (time-dependent) inflectional velocity profile if both walls are brought momentarily to rest. If the Reynolds number is high enough a growing stationary instability develops. This situation is ideally suited for flow visualization of the instability. In this paper we show flow visualization experiments and stability calculations of the developing transverse roll cell instability in such a flow at low Reynolds numbers. Although the stability calculations are based on a quasi-stationary velocity profile, the measured and most amplified wave length obtained from the calculations are in excellent agreement.     

Effect of pulsations on the stability of a gas column

The effect of radial pulsations on the stability of a compressible cylindrical gas column surrounded by an ambient liquid is discussed. In the absence of pulsations, the stationary interface is susceptible to the Rayleigh capillary instability, which promotes the growth of longitudinal waves whose wave length is larger than 2 times the column radius, irrespective of the Reynolds number. A Floquet stability analysis for potential flow shows that the pulsations further destabilize the interface by extending the range of unstable wave numbers to a sequence of islands. A similar stability analysis for Stokes flow shows that the pulsations also have a destabilizing influence, though the presence of an insoluble surfactant has a competing stabilizing influence that may cause an overall reduction in the range of unstable wave numbers.     

Invariant solutions in a channel flow using a minimal restricted nonlinear model

Simulations using a Restricted Nonlinear (RNL) system, where mean flow distortion resulting from Reynolds stress feedback regenerates rolls, is applied in a channel flow under subcritical conditions. This quasi-linear restriction of the dynamics is used to study invariant solutions located in the bulk of the flow found recently by Rawat et al. (2016) [14]. It is shown that the RNL system truncated to a single streamwise mode for the perturbation supports invariant solutions that are found to bifurcate from a relative periodic orbit into a travelling wave solution when the spanwise size is increasing. In particular, the travelling wave solution exhibits a spanwise localized structure that remains unchanged for large values of the spanwise extent as the invariant solution lying on the lower branch found by Rawat et al. (2016) [14]. In addition, travelling wave solutions provided by this minimal RNL system are self-similar with respect to the Reynolds number based on the centreline velocity, and the half-channel height varying from 2000 to 5000.     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Characteristics of fiber suspension flow in a rectangular channel

Velocity profile of fiber suspension flow in a rectangular channel is measured by pulsed ultrasonic Doppler velocimetry (PUDV), and the effect of fiber concentration and Reynolds number on the shape of the velocity profile is investigated. Five types of flow behavior are observed when fiber concentration increases or flow rate decreases progressively. The turbulent velocity profiles of fiber suspension can be described by a correlation with fiber concentration, nl³, and Reynolds number, Re as the main parameters. The presence of fiber in the suspension will reduce the turbulence intensity and thus reduce the turbulent momentum transfer. On the other hand, fibers in the suspension have the tendency to form fiber networks, which will increase the momentum transfer. The relative contribution of these two types of momentum flux will determine the final shape of the velocity profile.     

后台阶分离流动中大涡结构演变的数值模拟

本文对后台阶分离流动中涡结构的演变进行了大涡模拟，研究了流场结构的变化规律。详细讨论了随着雷诺数的增加流场结构的典型特征的变化规律，指出流场中的涡结构随着雷诺数的增大变得十分复杂和丰富，回流区的数目、大小及其出现的位置也显著地不同。这些结果与已有的一些实验值和流场显示结果是吻合的。在此基础上，进一步研究了高雷诺数时流场中大尺度涡结构的瞬时发展和演化过程，展示了其中大涡的产生、追随、吸引、合并和破碎等过程。对于高雷诺数情况，对大涡模拟得到的数值结果进行了统计，得到的时均速度分布以及台阶后方的回流区长度与现有的其他实验结果符合得很好。本研究是针对后台阶分离流动深入开展湍流控制以及两相流动研究的基础。     

A numerical study of flow through wavy-walled channels

A numerical procedure is developed for the analysis of flow in a channel whose walls describe a travelling wave motion. Following a perturbation method, the primitive variables are expanded in a series with the wall amplitude as the perturbation parameter. The boundary conditions are applied at the mean surface of the channel and the first-order perturbation quantities are calculated using the pseudospectral collocation method. Although limited by the linear analysis, the present approach is not restricted by the Reynolds number of the flow and the wave number and frequency of the wavy-walled channel. Using the computed wall shear stresses, the positions of flow separation and reattachment are determined. The variations in velocity and pressure with frequency of excitation are also presented. © 1998 John Wiley & Sons, Ltd.     

Magnetohydrodynamic peristaltic flow of a hyperbolic tangent fluid in a vertical asymmetric channel with heat transfer

In the present paper we discuss the magnetohydrodynamic (MHD) peristaltic flow of a hyperbolic tangent fluid model in a vertical asymmetric channel under a zero Reynolds number and long wavelength approximation. Exact solution of the temperature equation in the absence of dissipation term has been computed and the analytical ex- pression for stream function and axial pressure gradient are established. The flow is analyzed in a wave frame of reference moving with the velocity of wave. The expression for pressure rise has been computed numerically. The physical features of pertinent parameters are analyzed by plotting graphs and discussed in detail.     

Stability of wave forms of motion of a viscous fluid layer under tangential stress

We study the stability of wave flow of a viscous incompressible fluid layer subjected to tangential stress and an inclined gravity force with respect to long-wave disturbances.An asymptotic solution is constructed for the equations of the disturbed motion and the problem is reduced to the study of a second-order ordinary differential equation. It is shown that after loss of stability by a Poiseuille flow the laminar nature of the flow is not destroyed, but the form of the free surface acquires a wave-like profile. The Poiseuille regime is stable for low Reynolds numbers. The critical Reynolds number for wave flow is found, and the stability and instability regions are determined.     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

Experiments on the linear instability of flow in a wavy channel

The effects of wall corrugation on the stability of wall-bounded shear flows have been examined experimentally in plane channel flows. One of the channel walls has been modified by introduction of the wavy wall model with the amplitude of 4% of the channel half height and the wave number of 1.02. The experiment is focused on the two-dimensional travelling wave instability and the results are compared with the theory [J.M. Floryan, Two-dimensional instability of flow in a rough channel, Phys. Fluids 17 (2005) 044101 (also: Rept. ESFD-1/2003, Dept. of Mechanical and Materials Engineering, The University of Western Ontario, London, Ontario, Canada, 2003)]. It is shown that the flow is destabilized by the wall corrugation at subcritical Reynolds numbers below 5772, as predicted by the theory. For the present corrugation geometry, the critical Reynolds number is decreased down to about 4000. The spatial growth rates, the disturbance wave numbers and the distribution of disturbance amplitude measured over such wavy wall also agree well with the theoretical results.     

On the Lower Critical Reynolds Number for Flow in a Circular Pipe

Flow in a circular pipe is investigated experimentally at Reynolds numbers higher than that at which the resistance coefficients calculated from the Blasius formula for laminar flow and from the Prandtl formula for turbulent flow are equal. The corresponding Reynolds number based on the mean-flow velocity and the pipe diameter is about 1000. The experiments were performed at a high level of inlet pulsations produced by feeding gas into the pipe through a hole with a diameter several times smaller than the pipe diameter. In our experiments the critical Reynolds number was determined as the value, independent of the distance from the inlet, at which the ratio of the axial to the mean-flow velocity as a function of the Reynolds number deviated from 2. At the maximum ratio of the pipe cross-sectional area to the area of the hole through which the gas entered the pipe, equal to 26, the critical Reynolds number was about 2300. After a fivefold increase in the hole area the critical Reynolds number increased by approximately 4%.At Reynolds numbers below 2000, after at a high level of the inlet pulsations an almost laminar flow had developed in the pipe, a perturbation was introduced by inserting a diametrically oriented cylindrical rod with a diameter 10–20 times smaller than the pipe diameter. In these experiments, at Reynolds numbers higher than 1000, at a distance from the rod equal to 50 pipe diameters the axial to mean-flow velocity ratio was less than 2, approaching this value again at large distances from the rod. The insertion of the rod led to a decrease in the critical Reynolds number by approximately 12%.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.