Stability analysis of shear-thinning flow between rotating cylinders

Stability of the shear thinning Taylor–Couette flow is carried out and complete bifurcation diagram is drawn. The fluid is assumed to follow the Carreau–Bird model and mixed boundary conditions are imposed. The low-order dynamical system, resulted from Galerkin projection of the conservation of mass and momentum equations, includes additional nonlinear terms in the velocity components originated from the shear-dependent viscosity. It is observed, that the base flow loses its radial flow stability to the vortex structure at a lower critical Taylor number, as the shear thinning effects increases. The emergence of the vortices corresponds to the onset of a supercritical bifurcation which is also seen in the flow of a linear fluid. However, unlike the Newtonian case, shear-thinning Taylor vortices lose their stability as the Taylor number reaches a second critical number corresponding to the onset of a Hopf bifurcation. Complete flow field together with viscosity maps are given for different scenarios in the bifurcation diagram.     

Mathematical and numerical treatment of instabilities of non-axisymmetric confined vortices under the Dirichlet boundary conditions

This paper is concerned with an efficient analytical and numerical method for computing instability bounds of non-axisymmetric confined vortices. We propose a new approach based on a recently mathematical method of meshless collocation, when the shifted orthogonal Chebyshev base is adapted to satisfy the boundary conditions. Appropriately modified Tollmien–Schlichting modes used for the linear stability analysis offer the main advantage of dealing with a steady perturbation model of a confined vortex that leads to a variety of information describing the behavior of the full three-dimensional unsteady flow.     

Unsteady thermal boundary-layer on an infinite yawed wedge whose temperature gradient is prescribed

In this paper, the temperature distribution at the surface of an infinite yawed wedge, when temperature of the main-stream is constant, is studied under two cases; one when temperature gradient of the wall is steady and the velocity of the main-stream is unsteady, and the other when temperature gradient of the wall is unsteady and the velocity of the main-stream is steady. It is found that the heat transfer depends on the wedge angle, the angle parameter and the yaw of the wedge. The behaviour with these parameters are studied and are illustrated graphically.     

Exact solutions of problems of the three-dimensional flow of ideal viscous incompressible fluids near cylindrical surfaces

Three-dimensional flows of an incompressible fluid, the parameters of which depend on two coordinates and time, are considered. The stream surfaces of such flows are cylindrical. The equations of continuity and the Navier-Stokes equations can be transformed to relations, one of which is the equation for the stream function the other is the integral of the equations relating the pressure and the stream function, and the third is a linear equation for the projection of the velocity vector onto the axis parallel to the generatrix of the cylindrical surfaces. The problems of modelling the flows are considered on the basis of the exact solutions of the Navier-Stokes equations and Euler's equations using examples. Relations for the distribution of the flow parameters in the channel created by hyperbolical cylinders are derived for the case of unsteady inviscid flow. The streamlines of these flows are situated on the side surfaces of the hyperbolical cylinders and intercept the generatrices of the cylinders at certain indirect angles. The flow around a circular cylinder and the flow of fluid inside an elliptic cylinder are considered in the case of steady inviscid flow. The streamlines on the circular cylinder are arranged transverse to the cylinder (the projection of the velocity vector onto the coordinate axis, parallel to the generatrix of the cylinder, is equal to zero). Far from the cylinder the streamlines are also situated on a cylindrical surfaces, but not transverse to the cylinder, making certain indirect angles with the generatrix. Viscous three-dimensional flows, possessing a certain symmetry, are considered. In the case of radial symmetry the streamlines are helical lines. The non-planar Couette flow between parallel moving planes is characterized by the fact that the velocity vectors, being situated in the same plane, are collinear, while the velocity vectors in parallel planes are not collinear. Relations for viscous steady three-dimensional flows, using well-known relations, obtained for the stream function of two-dimensional flows, are given.     

Measurements of transient flow fields driven by a discontinuously applied rotating magnetic field

This experimental study considers the transient liquid metal flow which is generated inside a cylindrical container by discontinuously applying a rotating magnetic field (RMF). The focus is on the fluid motion arising from the impulsive spin–up from the resting state. The ultrasonic Doppler velocimetry (UDV) has been used to determine profiles of the fluid velocity in the ternary alloy GaInSn. The azimuthal and vertical velocity components have been measured allowing for an analysis of both the primary, swirling flow and the secondary flow in the radial–meridional plane. The experimental results show an excellent agreement with recently published numerical results. The investigations reveal that the recirculating flow in the radial–meridional plane undergoes characteristic oscillations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical computation of end plate effect on Taylor vortices between rotating conical cylinders

End plate effect on Taylor vortices between rotating conical cylinders is studied by numerical method in this paper. We suppose that the inner cone rotates together with the end plate at the top and the outer one as well as the end plate at the bottom remains at rest. It is found that the instability sets in at a critical Reynolds number about Re = 80. Increase Re to about Re = 200 the first single clockwise vortex is formed near the top of the flow system. Further increase Re to about Re = 440 another clockwise vortex is formed under the first one. At about Re = 448 the third vortex is formed which rotates in counterclockwise direction between the first two vortices. With increasing of Re the process continues. Finally, a configuration is obtained with an odd number of vortices in the annulus at about Re = 700, which confirms the experimental observation. Moreover, the local extreme values of pressure and velocity are achieved at the adjacent lines between neighboring vortices or at the medium lines of vortices. The effect of gap size on vortices is also considered. It is shown that the number of vortices increases with decreasing of the gap size and the end plates play an important role in the parity of the number of the vortices.     

Effects of chemical reaction on MHD free convection and mass transfer flow of a micropolar fluid with oscillatory plate velocity and constant heat source in a rotating frame of reference

This paper concerns with studying the steady and unsteady MHD micropolar flow and mass transfers flow with constant heat source in a rotating frame of reference in the presence chemical reaction of the first-order, taking an oscillatory plate velocity and a constant suction velocity at the plate. The plate velocity is assumed to oscillate in time with a constant frequency; it is thus assumed that the solutions of the boundary layer are the same oscillatory type. The governing dimensionless equations are solved analytically after using small perturbation approximation. The effects of the various flow parameters and thermophysical properties on the velocity and temperature fields across the boundary layer are investigated. Numerical results of velocity profiles of micropolar fluids are compared with the corresponding flow problems for a Newtonian fluid. The results show that there exists completely oscillating behavior in the velocity distribution.     

On a Control Problem for a Linear System with Measurements of a Part of Phase Coordinates

We consider a three-dimensional unsteady flow with a rotating detonation wave arising in an annular gap of an axially symmetric engine between two parallel planes perpendicular to its symmetry axis. The corresponding problem is formulated and studied. It is assumed that there is a reservoir with quiescent homogeneous propane–air combustible mixture with given stagnation parameters; the mixture flows from the reservoir into the annular gap through its external cylindrical surface toward the symmetry axis, and the parameters of the mixture are determined by the pressure in the reservoir and the static pressure in the gap. The detonation products flow out from the gap into a space bounded on one side by an impermeable wall that is an extension of a side of the gap. Through a hole on the other side of the gap and through a conical output section with a half-opening angle of 45°, the gas flows out from the engine into the external space. We formulate a model of detonation initiation by energy supply in which the direction of rotation of the detonation wave is defined by the position of the energy-release zone of the initiator with respect to the solid wall situated in a plane passing through the symmetry axis. After a while, this solid wall disappears (burns out). We obtain and analyze unsteady shock-wave structures that arise during the formation of a steady rotating detonation. The analysis is carried out within single-stage combustion kinetics by the numerical method based on the Godunov scheme with the use of an original software system developed for multiparameter calculations and visualization of flows. The calculations were carried out on the Lomonosov supercomputer at Moscow State University.     

A third-order Taylor–Galerkin model for the simulation of two-dimensional unsteady free surface flows

A Taylor–Galerkin third-order method is presented to integrate the equations describing two-dimensional, unsteady flows having a free surface. The discretization in time, with Taylor series expansions, is based on a fractional step, while the spatial approximation is obtained by the conventional Galerkin finite element method. Results are presented, in terms of the water profile history and flow velocity field, for two simulations: a partial and sudden breakage of a dam in a horizontal and frictionless river-bed, and the flow through a channel contraction.     

On the development of a wake vortex in inviscid flow

The evolution of an initial perturbation in an axisymmetric subsonic normal inviscid gas flow through a pipe is directly simulated. The basic (unperturbed) flow has a zero radial velocity component, while its axial velocity component (along the axis of symmetry) increases or decreases linearly with the radius. The perturbation is specified as a swirl (rotation about the axis) with a positive or negative velocity vanishing on the central axis and the lateral surface. Irrespective of its direction, the swirl gives rise to a steady-state vortex carried by the flow. It shape is spherical (contiguous to the rotation axis) or circular (sliding along the impermeable lateral surface).     

On the influence of plasma DBD actuator on the flow in a rectangular channel

Previously, the vortical structures generated by plasma DBD actuator working in unsteady regime were investigated in detail. The generalized model describing the behaviour of these vortices in dependency on input power parameters was introduced. This paper should reveal how the wall-jet-like-flow generated by that actuator will affect the developed flow in a rectangular channel with cross-section dimension of 250 × 100 mm. The actuator is considered in spanwise configuration where the induced flow has the same or opposite orientation as the main flow. The flow control will be tested both for steady regime and for unsteady regime. The dynamic of that complex phenomenon will be studied and results in qualitatively and quantitatively meaning will be presented. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

One-dimensional viscoelastic fluid model where viscosity and normal stress coefficients depend on the shear rate

We study the unsteady motion of a viscoelastic fluid modeled by a second-order fluid where normal stress coefficients and viscosity depend on the shear rate by using a power-law model. To study this problem, we use the one-dimensional nine-director Cosserat theory approach which reduces the exact three-dimensional equations to a system depending only on time and on a single spatial variable. Integrating the equation of conservation of linear momentum over the tube cross-section, with the velocity field approximated by the Cosserat theory, we obtain a one-dimensional system. The velocity field approximation satisfies both the incompressibility condition and the kinematic boundary condition exactly. From this one-dimensional system we obtain the relationship between average pressure and volume flow rate over a finite section of the tube with constant and variable radius. Also, we obtain the correspondent equation for the wall shear stress which enters directly in the formulation as a dependent variable. Attention is focused on some numerical simulation of unsteady/steady flows for average pressure, wall shear stress and on the analysis of perturbed flows.     

A modified discrete-vortex method algorithm with shedding criterion for aerodynamic coefficients prediction at high angle of attack

Low-order methods require less computing power than classical computational fluid dynamics and can be implemented on a laptop computer, which is needed for engineering tasks. Discrete vortex methods are such low order methods that can describe the unsteady separated flow around an airfoil. After a presentation of the leading edge suction parameter discrete vortex method, a modified algorithm is proposed, in order to reduce the computing cost, and compared with the previous one. Several reference unsteady airfoil motions are discussed in terms of gain in the computation time with comparisons between the previous scheme and the present one. The accuracy of the new method is demonstrated through aerodynamic coefficients. The application of the present discrete vortex method to a transient pitching motion of an airfoil is also presented, in order to understand the leading edge vortex formation, and its implication in terms of lift and drag coefficients. The method is not limited to unsteady or transient motions but can also simulate the flow around a constant angle of attack airfoil. In that case, an original method of fast summation of the vortices located far away from the airfoil, allows a linear dependence of the computation time versus the number of vortices shed, which is a great improvement over the quadratic dependence observed in the classical discrete vortex methods. The development of the aerodynamic coefficients with angle of attack, from values ranging between −10° and 90°, is obtained for a purely two-dimensional flow. In particular, the shape of the lift coefficient of the airfoil in the fully detached flow region is established. Comparisons with relevant experimental or computational fluid dynamics data are discussed in order to grasp the influence of upstream turbulence level and three-dimensional effects in the measured data in the fully detached flow region.     

Numerical solution of steady and unsteady flow over a vibrating airfoil in a channel

The work deals with a numerical solution of 2D steady and unsteady inviscid incompressible flow over the profile NACA 0012 in a channel. The finite volume method (FVM) in a form of cell-centered explicit schemes at quadrilateral C-mesh is used. Governing system of equations is the system of incompressible Euler equations. The method of artificial compressibility and time dependent method is applied to steady computations. The small disturbance theory (SDT) applied to a numerical solution of flow over a rotated profile by a small angle only is mentioned. Brief introduction is given to the Arbitrary (Semi) Lagrangian-Eulerian (ALE) method used for unsteady computations. Some numerical results of unsteady flow over a vibrating profile achieved by both SDT and ALE method are presented. Unsteady flow is caused by prescribed oscillations of the profile (one degree of freedom) fixed in an elastic axis. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The stress state of an elastic composite cone with centre of rotation at the apex

The torsion of a composite cone that has a centre of rotation at its apex is investigated in a spherical system of coordinates. A composite cone is a cone with one shear modulus, inserted into a conical funnel having another shear modulus and with ideal mechanical contact between its surface and the inner surface of the conical funnel. The auxiliary problem of a composite cone with its apex truncated by a spherical surface is considered first. The outer surface of such a conical body is not loaded, but a load that reduces to a torque is applied to its spherical surface. The auxiliary problem is reduced to a one-dimensional discontinuous boundary-value problem using a specially constructed integral transformation. The exact solution of this boundary-value problem is constructed. The limit is then taken in the solution obtained as the radius of the spherical surface tends to zero for the purpose of obtaining an exact solution of the problem of the torsion of a composite cone that has a centre of rotation at the apex.     

圆锥形血管中的振荡发展流动

本文在小锥度角的假设下,研究了圆锥形血管的非定常振荡的发展流动问题.导得了相应的速度分布公式.分析表明,所有收缩的圆锥形血管的流动都是发展流动,而且锥度角对发展流动的影响随着锥度角的增大而增大.     

Computations of transport of pollutant in shallow water

The focus of this paper is to simulate the transport of a passive pollutant by a flow modelled by the two-dimensional shallow water equations. Considering the friction terms, new model for simulating the steady and unsteady transport of pollutant is established. Then the adaptive semi-discrete central-upwind scheme based on central weighted essentially non-oscillatory reconstruction is utilized for simulating the two-dimensional steady and unsteady transport of pollutant. The non-oscillatory behavior and accuracy of the scheme are demonstrated by the numerical result.     

Unsteady viscous flow over a rotating stretchable disk with deceleration

In this work, the laminar unsteady flow over a stretchable rotating disk with deceleration is investigated. The three dimensional Navier–Stokes (NS) equations are reduced into a similarity ordinary differential equation group, which is solved numerically using a shooting method. Mathematically, two solution branches are found for the similarity equations. The lower solution branch may not be physically feasible due to a negative velocity in the circumferential direction. For the physically feasible solution branch, namely the upper solution branch, the fluid behavior is greatly influenced by the disk stretching parameter and the unsteadiness parameter. With disk stretching, the disk can be friction free in both the radial and the circumferential directions, depending on the values of the controlling parameters. The results provide an exact solution to the whole unsteady NS equations with new nonlinear phenomena and multiple solution branches.     

圆截面螺旋管道内非定常流动研究

以血液流动为背景,利用双参数摄动法研究了圆截面螺旋管内低频振荡流动,得到问题的二阶摄动解,分析了轴向速度、二次流、壁面剪应力在不同时刻的特点及随时间和Womersley数的变化情况。研究表明:挠率对圆截面曲线管道内低频振荡流动的影响不可忽略,尤其是轴向压力梯度绝对值很小时,挠率将对二次流动结构起主要影响作用。流函数的剧烈变化只发生在正、负数值发生转变的很小的时间段内,大部分时间段内变化平缓。壁面剪应力随θ的变化也很大。     

A numerical simulation of flow between two rotating coaxial frustum cones

The behavior of flow between two coaxial frustum cones, with the inner one rotating and the outer stationary, is studied in this paper. It is found that the fluid at the outlet does not flow out directly, but flows up till a certain height. This reflux generates a vortex area with a quite large velocity and pressure magnitude. This reflux area, between Z/H = 0.05 and 0.30, has the trend to move up with increasing Reynolds number Re. The velocity magnitude is linear in the radial direction if the Re is small. This linear relation converts to quasi-quadratic function as the Re increasing. If the frustum cone inclination is small, the flow will tend to be unstable with a quite large velocity and pressure magnitude. Finally, a comparison is made with Taylor–Couette flow.