Boundary layer flow of Maxwell fluid due to torsional motion of cylinder: modeling and simulation

This paper investigates the boundary layer flow of the Maxwell fluid around a stretchable horizontal rotating cylinder under the influence of a transverse magnetic field. The constitutive flow equations for the current physical problem are modeled and analyzed for the first time in the literature. The torsional motion of the cylinder is considered with the constant azimuthal velocity E. The partial differential equations (PDEs) governing the torsional motion of the Maxwell fluid together with energy transport are simplified with the boundary layer concept. The current analysis is valid only for a certain range of the positive Reynolds numbers. However, for very large Reynolds numbers, the flow becomes turbulent. Thus, the governing similarity equations are simplified through suitable transformations for the analysis of the large Reynolds numbers. The numerical simulations for the flow, heat, and mass transport phenomena are carried out in view of the bvp4c scheme in MATLAB. The outcomes reveal that the velocity decays exponentially faster and reduces for higher values of the Reynolds numbers and the flow penetrates shallower into the free stream fluid. It is also noted that the phenomenon of stress relaxation, described by the Deborah number, causes to decline the flow fields and enhance the thermal and solutal energy transport during the fluid motion. The penetration depth decreases for the transport of heat and mass in the fluid with the higher Reynolds numbers. An excellent validation of the numerical results is assured through tabular data with the existing literature.     

On the orbital motion of a rotating inner cylinder in annular flow

In this paper, numerical calculations have been performed to analyse the influence of the orbital motion of an inner cylinder on annular flow and the forces exerted by the fluid on the inner cylinder when it is rotating eccentrically. The flow considered is fully developed laminar flow driven by axial pressure gradient. It is shown that the drag of the annular flow decreases initially and then increases with the enhancement of orbital motion, when it has the same direction as the inner cylinder rotation. If the eccentricity and rotation speed of the inner cylinder keep unchanged (with respect to the absolute frame of reference), and the orbital motion is strong enough that the azimuthal component (with respect to the orbit of the orbital motion) of the flow‐induced force on the inner cylinder goes to zero, the flow drag nearly reaches its minimum value. When only an external torque is imposed to drive the eccentric rotation of the inner cylinder, orbital motion may occur and, in general, has the same direction as the inner cylinder rotation. Under this condition, whether the inner cylinder can have a steady motion state with force equilibrium, and even what type of motion state it can have, is related to the linear density of the inner cylinder. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study on the vortex motion patterns around a rotating circular cylinder and their critical characters

A hybrid finite difference and vortex method (HFDV), based on the domain decomposition method (DDM), is used for calculating the flow around a rotating circular cylinder at Reynolds number Re=1000, 200 and the angular‐to‐rectilinear speed ratio α∈(0.5, 3.25) respectively. A fully implicit third‐order eccentric finite difference scheme is adopted in the finite difference method, and the deduced large broad band sparse matrix equations are solved by a highly efficient modified incomplete LU decomposition conjugate gradient method (MILU‐CG). The long‐time, fully developed features about the variations of the vortex patterns in the wake, as well as the drag and lift forces on the cylinder, are given. The calculated streamline contours are in good agreement with the experimentally visualized flow pictures. The existence of the critical state is confirmed again, and the single side shed vortex pattern at the critical state is shown for the first time. Also, the optimized lift‐to‐drag force ratio is obtained near the critical state. Copyright © 1999 John Wiley & Sons, Ltd.     

Thermal response of a periodic boundary layer near an axisymmetric stagnation point on a circular cylinder

Boundary layer solutions are provided to study the time-mean heat transfer characteristics in a laminar flow in the vicinity of an axisymmetric stagnation point. The velocity of the oncoming flow is assumed to oscillate relative to the body. Different solutions are constructed for the small and high values of the reduced frequency parameter. Numerical solutions for the temperature functions are presented, and the wall values of the thermal gradients are tabulated.     

纳米流体在伸/缩楔体上的MHD流动数值研究

对纳米流体在伸/缩楔体上的磁流体(MHD)流动进行了数值研究。首先,通过相似变换将控制偏微分方程转化为非线性常微分方程组;然后,利用Matlab软件,借助打靶法,结合四阶五常龙格库塔迭代方案进行数值求解;最后,详细讨论了各控制参数对无量纲速度、温度、浓度、表面摩擦系数、局部Nusselt数和局部Sherwood数的影响。结果表明,楔体在拉伸情况下只有唯一解,理论上不会出现边界层分离;而在一定收缩强度范围内存在双解,边界层流动在壁面处可能会出现边界层分离,壁面抽吸会使边界层分离推迟;楔体在拉伸情况下,磁场参数对表面摩擦系数的影响较大,对局部Nusselt数和局部Sherwood数的影响较小。     

Hydrodynamics of flow over a transversely oscillating circular cylinder beneath a free surface

In this paper, hydrodynamic force coefficients and wake vortex structures of uniform flow over a transversely oscillating circular cylinder beneath a free surface were numerically investigated by an adaptive Cartesian cut-cell/level-set method. At a fixed Reynolds number, 100, a series of simulations covering three Froude numbers, two submergence depths, and three oscillation amplitudes were performed over a wide range of oscillation frequency. Results show that, for a deeply submerged cylinder with sufficiently large oscillation amplitudes, both the lift amplitude jump and the lift phase sharp drop exist, not accompanied by significant changes of vortex shedding timing. The near-cylinder vortex structure changes when the lift amplitude jump occurs. For a cylinder oscillating beneath a free surface, larger oscillation amplitude or submergence depth causes higher time-averaged drag for frequency ratio (=oscillation frequency/natural vortex shedding frequency) greater than 1.25. All near-free-surface cases exhibit negative time-averaged lift the magnitude of which increases with decreasing submergence depth. In contrast to a deeply submerged cylinder, occurrences of beating in the temporal variation of lift are fewer for a cylinder oscillating beneath a free surface, especially for small submergence depth. For the highest Froude number investigated, the lift frequency is locked to the cylinder oscillation frequency for frequency ratios higher than one. The vortex shedding mode tends to be double-row for deep and single-row for shallow submergence. Proximity to the free surface would change or destroy the near-cylinder vortex structure characteristic of deep-submergence cases. The lift amplitude jump is smoother for smaller submergence depth. Similar to deep-submergence cases, the vortex shedding frequency is not necessarily the same as the primary-mode frequency of the lift coefficient. The frequency of the induced free surface wave is exactly the cylinder oscillation frequency. The trends of wave length variation with the Froude number and frequency ratio agree with those predicted by the linear theory of small-amplitude free surface waves.     

Numerical study of adverse pressure gradient generation over a flat plate using a rotating cylinder

Generating an adverse pressure gradient (APG), using a rotating cylinder in the proximity of a plane wall under a laminar freestream flow, is studied numerically in this work. The magnitude of the generated APG is a function of the gap, G, between the cylinder and the wall, and the rotational speed of the cylinder, Ω. The flow in such a configuration is characterized by periodic transient vortex shedding at high Reynolds number. A numerical model for the computation of the transient flow for this configuration is developed using the ANSYS CFD simulation tool. The model is validated against published experimental and numerical data for similar flow configurations and excellent agreement is observed. A parametric study is carried out for different combinations of G and Ω for two different Reynolds numbers of 200 and 1000 to examine the development of the resulting separation bubble due to the generated APG. The mechanism of the boundary layer separation over the plane wall and the corresponding wake dynamics is investigated. Results are presented in terms of the distribution of the pressure coefficient as well as skin friction coefficient along the wall and flow patterns around and downstream of the cylinder in the proximity of the wall. The results of these computations confirm that using a rotating cylinder over a plane wall in a freestream flow is an effective technique to generate a controlled range of adverse pressure gradients.     

Entropy analysis of unsteady magneto-nanofluid flow past accelerating stretching sheet with convective boundary condition

The unsteady laminar magnetohydrodynamics(MHD) boundary layer flow and heat transfer of nanofluids over an accelerating convectively heated stretching sheet are numerically studied in the presence of a transverse magnetic field with heat source/sink. The unsteady governing equations are solved by a shooting method with the Runge-KuttaFehlberg scheme. Three different types of water based nanofluids, containing copper, aluminium oxide, and titanium dioxide, are taken into consideration. The effects of the pertinent parameters on the fluid velocity, the temperature, the entropy generation number, the Bejan number, the shear stress, and the heat transfer rate at the sheet surface are graphically and quantitatively discussed in detail. A comparison of the entropy generation due to the heat transfer and the fluid friction is made with the help of the Bejan number. It is observed that the presence of the metallic nanoparticles creates more entropy in the nanofluid flow than in the regular fluid flow.     

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder atRe=200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number. The project supported by the State Key Fundamental Research Project of “Large Scale Scientific Computation Research” (G199903281)     

A computational solution of hydromagnetic‐free convective flow past a vertical plate in a rotating heat‐generating fluid with Hall and ion‐slip currents

Taking Hall and ion‐slip current into account, the unsteady magnetohydrodynamic heat‐generating free convective flow of a partially ionized gas past an infinite vertical plate in a rotating frame of reference is investigated theoretically. A computer program using finite elements is employed to solve the coupled non‐linear differential equations for velocity and temperature fields. The effects of Hall and ion‐slip currents as well as the other parameters entering into the problem are discussed extensively and shown graphically. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of laminar flow past a circular cylinder with slip conditions

The no‐slip condition is an assumption that cannot be derived from first principles and a growing number of literatures replace the no‐slip condition with partial‐slip condition, or Navier‐slip condition. In this study, the influence of partial‐slip boundary conditions on the laminar flow properties past a circular cylinder was examined. Shallow‐water equations are solved by using the finite element method accommodating SU/PG scheme. Four Reynolds numbers (20, 40, 80, and 100) and six slip lengths were considered in the numerical simulation to investigate the effects of slip length and Reynolds number on characteristic parameters such as wall vorticity, drag coefficient, separation angle, wake length, velocity distributions on and behind the cylinder, lift coefficient, and Strouhal number. The simulation results revealed that as the slip length increases, the drag coefficient decreases since the frictional component of drag is reduced, and the shear layer developed along the cylinder surface tends to push the separation point away toward the rear stagnation point so that it has larger separation angle than that of the no‐slip condition. The length of the wake bubble zone was shortened by the combined effects of the reduced wall vorticity and wall shear stress which caused a shift of the reattachment point closer to the cylinder. The frequency of the asymmetrical vortex formation with partial slip velocity was increased due to the intrinsic inertial effect of the Navier‐slip condition. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical study of vortex shedding from a rotating cylinder immersed in a uniform flow field

A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=10³ and α≤3; the other with Re=10⁴ and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd.     

Melting heat transfer in the stagnation‐point flow of an upper‐convected Maxwell (UCM) fluid past a stretching sheet

The steady laminar boundary layer flow and heat transfer past a stretching sheet arre considered. Upper‐convected Maxwell (UCM) fluid is treated as a rheological model. The resulting nonlinear differential system is solved by homotopy analysis method (HAM). The influence of melting parameter (M), Prandtl number (Pr), Deborah number (β) and stretching ratio (A = a/c) on the velocity and temperature profiles is thoroughly examined. It is noticed that fields are effected appreciably with the variation of parameters. Furthermore, it is seen that the local Nusselt number is a decreasing function of melting parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

Transient analytical solution for the motion of a vibrating cylinder in the Stokes regime using Laplace transforms

A new analytical solution for the motion of an elastic cylinder in a viscous fluid is derived using Laplace transforms. Unlike previously available solutions, full expressions for transient terms are given. The solution is compared with conventional treatments of this problem. It is expected to have particular value for applications related to viscosity measurement using vibrating-wire viscometers applied to higher viscosity fluids.     

Slip effects and heat transfer analysis in a viscous fluid over an oscillatory stretching surface

In this study, we investigate the heat transfer problem in a viscous fluid over an oscillatory infinite sheet with slip condition. The sheet is moved back and forth in its own plane. The derived problem involves a dimensionless parameter indicating the relative magnitude of frequency to sheet stretching rate. A system of non‐linear partial differential equations is solved numerically using the finite‐difference scheme, in which a coordinate transformation is employed to transform the semi‐infinite physical space to a bounded computational domain. The physical features of interesting parameters on the velocity and temperature distributions are shown graphically and discussed. The values of the skin‐friction coefficient and the local Nusselt number are given in tabular form. Copyright © 2008 John Wiley & Sons, Ltd.     

Weakly nonlinear stability analysis of triple diffusive convection in a Maxwell fluid saturated porous layer

The weakly nonlinear stability of the triple diffusive convection in a Maxwell fluid saturated porous layer is investigated. In some cases, disconnected oscillatory neutral curves are found to exist, indicating that three critical thermal Darcy-Rayleigh numbers are required to specify the linear instability criteria. However, another distinguishing feature predicted from that of Newtonian fluids is the impossibility of quasi-periodic bifurcation from the rest state. Besides, the co-dimensional two bifurcation points are located in the Darcy-Prandtl number and the stress relaxation parameter plane. It is observed that the value of the stress relaxation parameter defining the crossover between stationary and oscillatory bifurcations decreases when the Darcy-Prandtl number increases. A cubic Landau equation is derived based on the weakly nonlinear stability analysis. It is found that the bifurcating oscillatory solution is either supercritical or subcritical, depending on the choice of the physical parameters. Heat and mass transfers are estimated in terms of time and area-averaged Nusselt numbers.     

Simulation of the flow past a circular cylinder in the supercritical regime by blending RANS and variational-multiscale LES models

A strategy which blends a variational multiscale large eddy simulation (VMS-LES) model and a RANS model in a hybrid approach is investigated. A smooth blending function, which is based on the value of a blending parameter, is used for switching from VMS-LES to RANS. Different definitions of the blending parameter are investigated. The capabilities of the novel hybrid approach are appraised in the simulation of the flow around a circular cylinder at a Reynolds number 1.4×10⁵, based on the freestream velocity and on the cylinder diameter, in the presence of turbulent boundary-layer due to turbulent inflow conditions. A second study at Reynolds numbers from Re=6.7×10⁵ to 1.25×10⁶ is also presented. The effect of using the VMS-LES approach in the hybrid model is evaluated. Results are compared to those of other RANS, LES and hybrid simulations in the literature and with experimental data     

The Eckert number phenomenon: Experimental investigations on the heat transfer from a moving wall in the case of a rotating cylinder

The Eckert number phenomenon was investigated theoretically by Geropp in 1969 and describes a reversal in heat transfer from a moving wall at an Eckert number Ec ≈ 1. In this report the Eckert number phenomenon is confirmed experimentally for the first time. For that purpose the heat transfer from a heated, vertically rotating cylinder in a crossflow was investigated. In order to perform the experiments in a range where the predicted phenomenon occurs, extreme rotational speeds were necessary. A heating concept had to be developed which allowed an input of heating power independent of the speed and which therefore had to be contact-free. The results show, among other things, that the temperature difference between the wall and the surrounding fluid has a significant effect on the predicted reversal of heat transfer at the wall. Moreover, maximum heat transfer occurs at an Eckert number Ec ≈ 0.3, which is of great importance for the cooling of hot surfaces in a gas-flow.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.