MHD duct flows under hydrodynamic “slip” condition

Three classic MHD problems are revisited assuming hydrodynamic slip condition at the interface between the electrically conducting fluid and the insulating wall: (1) Hartmann flow; (2) fully developed flow in a rectangular duct; and (3) quasi-two-dimensional (Q2D) turbulent flow. The first two problems have been solved analytically. Additionally to the Hartmann number (Ha), a new dimensionless parameter S, the ratio of the slip length to the thickness of the Hartmann layer, has been identified. One of the most important conclusions of the paper is that the duct flows with the slip still exhibit Hartmann layers, whose thickness scales as 1/Ha, while the thickness of the side layers is a function of both Ha and S. In the case of Q2D flows, a new expression for the Hartmann braking time has been derived showing its increase at Ha >> 1 by factor (1+ S). Numerical simulations performed for a flow with the “M-shaped” velocity profile show that in the presence of the slip, a Q2D flow becomes more irregular as vortical structures experience less Joule and viscous dissipation in the Hartmann layers.     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.     

Empirical slip and viscosity model performance for microscale gas flow

For the simple geometries of Couette and Poiseuille flows, the velocity profile maintains a similar shape from continuum to free molecular flow. Therefore, modifications to the fluid viscosity and slip boundary conditions can improve the continuum based Navier–Stokes solution in the non‐continuum non‐equilibrium regime. In this investigation, the optimal modifications are found by a linear least‐squares fit of the Navier–Stokes solution to the non‐equilibrium solution obtained using the direct simulation Monte Carlo (DSMC) method. Models are then constructed for the Knudsen number dependence of the viscosity correction and the slip model from a database of DSMC solutions for Couette and Poiseuille flows of argon and nitrogen gas, with Knudsen numbers ranging from 0.01 to 10. Finally, the accuracy of the models is measured for non‐equilibrium cases both in and outside the DSMC database. Flows outside the database include: combined Couette and Poiseuille flow, partial wall accommodation, helium gas, and non‐zero convective acceleration. The models reproduce the velocity profiles in the DSMC database within an L₂ error norm of 3% for Couette flows and 7% for Poiseuille flows. However, the errors in the model predictions outside the database are up to five times larger. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling spurt and stress oscillations in flows of molten polymers

The paper presents an approach for modeling polymer flows with non-slip, slip and changing non-slip — slip boundary conditions at the wall. The model consists of a viscoelastic constitutive equation for polymer flows in the bulk, prediction of the transition from non-slip to sliding boundary conditions, a wall slip model, and a model for the compressibility effects in capillary polymer flows. The bulk viscoelastic constitutive equation contains a hardening parameter which is solely determined by the polymer molecular characteristics. It delimits the conditions for the onset of solid, rubber-like behavior. The non-monotone wall slip model introduced for polymer melts, modifies a slip model derived from a simple stochastic model of interface molecular dynamics for cross-linked elastomers. The predictions for the onset of spurt, as well as the numerical simulations of hysteresis, spurt, and stress oscillations are demonstrated. They are also compared with available data for a high molecular weight, narrow distributed polyisoprene. By using this model beyond the critical conditions, many of the qualitative features of the spurt and oscillations observed in capillary and Couette flows of molten polymers, are described.Notations upper convected derivative of elastic strain tensor - f, f_m, f_min dimensionless (sliding) shear friction characteristics, and its maximum and minimum - G Hookean elastic modulus - G_p plateau modulus - G, G storage and loss moduli - I₁, I₂ first and second invariant of strain tensor - I₁, I₀ capillary and barrel lengths - M non-dimensional mass flow rate - M_C critical molecular weight - M^*, M_e molecular weights of a statistical segment, and of polymer chain between entanglements - M_n, M_W number average and weight average molecular weights - m, k two fitting parameters of slip model - _s , _s ^o nominal and characteristic sliding velocities - u non-dimensional sliding velocity - u _sc initial (infinitesimal) slip velocity - u ₁ upper limit of u on the lower branch - u ₂ lower limit of u on the upper branch - u _max value of u corresponding to f_min - u _min value of u corresponding to f_max - U piston speed - Q nominal volumetric flow rate - q non-dimensional volumetric flow rate - R, R_o capillary and barrel radii - M non-dimensional mass flow rate     

Viscous dissipation effect on heat transfer characteristics of rectangular microchannels under slip flow regime and H1 boundary conditions

Viscous dissipation effect on heat transfer characteristics of a rectangular microchannel is studied. Flow is governed by the Navier–Stokes equations with the slip flow and temperature jump boundary conditions. Integral transform technique is applied to derive the temperature distribution and Nusselt number. The velocity distribution is taken from literature. The solution method is verified for the case where viscous dissipation is neglected. It is found that, the viscous dissipation is negligible for gas flows in microchannels, since the contribution of this effect on Nu number is about 1%. However, this effect should be taken into account for much more viscous flows, such as liquid flows. Neglecting this effect for a flat microchannel with an aspect ratio of 0.1 for Br=0.04 underestimates the Nu number about 5%.     

Forced convection of gaseous slip-flow in porous micro-channels under Local Thermal Non-Equilibrium conditions

Steady laminar forced convection gaseous slip-flow through parallel-plates micro-channel filled with porous medium under Local Thermal Non-Equilibrium (LTNE) condition is studied numerically. We consider incompressible Newtonian gas flow, which is hydrodynamically fully developed while thermally is developing. The Darcy–Brinkman–Forchheimer model embedded in the Navier–Stokes equations is used to model the flow within the porous domain. The present study reports the effect of several operating parameters on velocity slip and temperature jump at the wall. Mainly, the current study demonstrates the effects of: Knudsen number (Kn), Darcy number (Da), Forchheimer number (Γ), Peclet number (Pe), Biot number (Bi), and effective thermal conductivity ratio (K _R) on velocity slip and temperature jump at the wall. Results are given in terms of skin friction (C _f Re ^*) and Nusselt number (Nu). It is found that the skin friction: (1) increases as Darcy number increases; (2) decreases as Forchheimer number or Knudsen number increases. Heat transfer is found to (1) decreases as the Knudsen number, Forchheimer number, or K _R increases; (2) increases as the Peclet number, Darcy number, or Biot number increases.     

Kinetic assessment of measured mass flow rates and streamwise pressure distributions in microchannel gas flows

Measured mass flow rates and streamwise pressure distributions of gas flowing through microchannels were reported by many researchers. Assessment of these data is crucial before they are used in the examination of slip models and numerical schemes, and in the design of microchannel elements in various MEMS devices. On the basis of kinetic solutions of the mass flow rates and pressure distributions in microchannel gas flows, the measured data available are properly normalized and then are compared with each other. The 69 normalized data of measured pressure distributions are in excellent agreement, and 67 of them are within 1 ± 0.05. The normalized data of mass flow-rates ranging between 0.95 and 1 agree well with each other as the inlet Knudsen number Kn _i < 0.02, but they scatter between 0.85 and 1.15 as Kn _i > 0.02 with, to some extent, a very interesting bifurcation trend. The project supported by the National Natural Science Foundation of China (90205024, 10621202 and 10425211).     

A numerical study of pulsatile laminar flows in a pipe with a ring-type constriction

Numerical simulations have been carried out to study pulsatile laminar flows in a pipe with an axisymmetric ringtype constriction. Three types of pulsatile flows were investigated, namely a physiological flow, a pure sinusoidal flow and a non-zero mean velocity sinusoidal flow. The laminar flow governing equations were solved by the SIMPLE algorithm on a non-staggered grid and a modified Crank-Nicolson approximation was used to discretrize the momentum equations with respect to time. The maximum flow Reynolds numer (Re) is 100. The Womersley number (N_w) ranges from 0 to 50, with the corresponding Strouhal number (St) ranging from 0 to 3·98. The constriction opening ratio (d/D) and thickness ratio (h/D) are fixed at 0·5 and 0·1 respectively. Within the time period investigated, all these pulsatile flows include both forward and backward flows. The unsteady recirculation region and the recirculation points change in size and location with time. For N_w ≤ 1 and St≤ 1·56 x 10^?3 the three pulsatile flows have the same simple relation between the instantaneous flow rate and pressure loss (Δp) across the constriction and the pressure gradient in the axial direction (dp/dz) in the fully developed flow region. The phase angles between the flow rate and pressure loss and the pressure gradient are equal to zero. With increasing N_w and St, the phase angle between the flow rate and the dp/dz becomes larger and has its maximum value of 90° at N_w = 50 and St = 3·98. The three pulsatile flows also show different relations between the flow rate and the pressure gradient. The pure sinusoidal flow has the largest maximum pressure gradient and the non-zero mean velocity sinusoidal flow has the smallest. For larger N_w and St the fully developed velocity profiles in the fully developed flow region have a smaller velocity gradient along the radial direction in the central region. The maximum recirculation length increases for N_w ranging from 0 to 4·2, while this length becomes very small at N_w = 50 and St = 3·98. The deceleration tends to enlarge the recirculation region and this effect appears for N_w ≥ 3 and St ≥ 1·43×10^?2. Linear relations exist between the flow rate and the instantaneous maximum values of velocity, vorticity and shear stress.     

Slip effects on unsteady mixed convection of hybrid nanofluid flow near the stagnation point

The unsteady mixed convection of the Al2O3-Cu/H2O hybrid nanofluid flow near the stagnation point past a vertical plate is analyzed. The bvp4c technique is used to solve the resulting ordinary differential equations. The combined effects of the velocity and thermal slip are addressed. The effects of different relevant physical parameters are studied numerically. The results show that the heat transfer rate is reduced when the volume fraction of the nanoparticles increases, while the unsteadiness...     

On Low-Dimensional Modeling of Channel Turbulence

We investigate a sequence of low-dimensional models of turbulent channel flows. These models are based on the extraction of the Karhunen–Loève (KL) eigenfunctions from a large-scale simulation in a wide channel with R _*=180 (based on the friction velocity). KL eigenfunctions provide an optimal coordinate system in which to represent the dynamics of the turbulent flow. The hierarchy of KL modes identifies the most energetic independent phenomena in the system. A series of Galerkin projections is then used to derive a dynamical approximation to flows. In order to capture essential aspects of the flow in a low-dimensional system, a careful selection of modes is carried out. The resulting models satisfy momentum and energy conservation as well as yielding the amount of viscous dissipation found in the exact system. Their dynamics includes modes which characterize the flux, rolls, and propagating waves. Unlike previous treatments the instantaneous streamwise flow is time dependent and represented by KL flux modes. The rolls correspond to the streaks observed in experiments in the viscous sublayer. Propagating waves which first appear in the model flow at low Reynolds number (R _*∼ 60) persist through the chaotic regime that appears as the Reynolds number is increased. Statistical measures of the chaotic flows which have been generated by the models compare favorably with those found in full-scale simulations. Received 13 July 1998 and accepted 8 January 1999     

Experimental estimation of wall slip for filled rubber compounds

The phenomenon of wall slip during flow of rubber compounds through capillaries is investigated for a typical styrene-butadiene elastomer with carbon black. It was found that at low temperature (110°C) the dependencies of slip velocity V _c on shear stress are described by the power law but, additionally, V _c depends on radius of a channel. At high temperatures there is a critical shear stress below which sliding is absent. Sliding appears only at higher shear stresses where, again, V _c depends on shear stress and the radius of a channel.     

Analysis of Streamwise Velocity Fluctuations in Turbulent Pipe Flow with the Use of an Ultrasonic Doppler Flowmeter

Data collected from several studies of experimental and numerical nature in wall-bounded turbulent flows and in particular in internal flows (channel and pipe flows, Mochizuki and Nieuwstadt [1]) at different Reynolds numbers R ⁺(Ru _*/ν), indicate that: (i) the peak of the rms-value (normalized by u _*) of the streamwise velocity fluctuations (σ _u ⁺|_peak) is essentially independent of the Reynolds number, (ii) the position of the rms peak value (y ⁺|_peak) is weakly dependent of the Reynolds number, (iii) the skewness of the streamwise velocity fluctuations (S _u ) is close to zero at the position in which the variance has its peak. A series of measurements of streamwise velocity fluctuations has been performed in turbulent pipe flow with the use of an Ultrasonic Doppler Velocimeter and our results support those reported in [1]. This revised version was published online in July 2006 with corrections to the Cover Date.     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part II Buoyant flow

This paper presents a derivation of an explicit algebraic model for two-dimensional (2-D) buoyant flows. It is an extension of the work reported in Part I (So et al. [27]). The tensor representation method of Jongen and Gatski [14] is extended to derive an explicit algebraic Reynolds stress model (EASM) for 2-D buoyant flow invoking the Boussinesq approximation. The projection methodology is further extended to treat the heat flux transport equation in the derivation of an explicit algebraic heat flux model (EAHFM) for buoyant flow. Again, the weak equilibrium assumption is invoked for the scaled Reynolds stress and scaled heat flux equation. An explicit algebraic model for buoyant flows is then formed with the EASM and EAHFM. From the derived EAHFM, an expression for the thermal diffusivity tensor in buoyant shear flows is deduced. Furthermore, a turbulent Prandtl number (Pr_T) for each of the three heat flux directions is determined. These directional Pr_T are found to be a function of the gradient Richardson number. Alternatively, a scalar Pr_T can be derived; its value is compared with the directional Pr_T. The EASM and EAHFM are used to calculate 2-D homogeneous buoyant shear flows and the results are compared with direct numerical simulation data and other model predictions, where good agreement is obtained. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

Effects of wall permeability on turbulence

In order to understand the effects of the wall permeability on turbulence near a porous wall, flow field measurements are carried out for turbulent flows in a channel with a porous bottom wall by a two-component particle image velocimetry (PIV) system. The porous media used are three kinds of foamed ceramics which have almost the same porosity (0.8) but different permeability. It is confirmed that the flow becomes more turbulent over the porous wall and tends to be turbulent even at the bulk Reynolds number of Re_b=1300 in the most permeable wall case tested. Corresponding to laminar to turbulent transition, the magnitude of the slip velocity on the porous wall is found to increase drastically in a narrow range of the Reynolds number. To discuss the effects of the wall roughness and the wall permeability, detailed discussions are made of zero-plane displacement and equivalent wall roughness for porous media. The results clearly indicate that the turbulence is induced by not only the wall roughness but the wall permeability. The measurements have also revealed that as Re_b or the wall permeability increases, the wall normal fluctuating velocity near the porous wall is enhanced due to the effects of the wall permeability. This leads to the increase of the turbulent shear stress resulting in higher friction factors of turbulence over porous walls.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

Near critical phenomena in laminar boundary layers

Asymptotic analysis of boundary layer separation in the limit of large Reynolds number Re→∞ has shown that in a number of cases which are of importance from a practical point of view solutions of the resulting interaction equations describing two-dimensional (2-D) steady flows exist up to a limiting value Γ_c of the relevant controlling parameter Γ only while two branches of solutions exist in a regime Γ<Γ_c. The present study aims at a better understanding of near critical flows |Γ-Γ_c|→0 and in particular the changes of the flow behaviour associated with the passage of Γ through Γ_c.     

The “Twin capillary” a simple device to separate shear- and slip-flow of fluids

When the flow behaviour of fluids is investigated with capillary-or rotational rheometers, adhesion of the fluid to the wall is normally one of the boundary conditions. For many fluids, especially for suspensions, this assumption is not valid. These fluids tend to slip at the wall. Therefore the normal evaluation of rheometer measurements leads to apparent but not compatible flow functions. The flow behaviour of these fluids can be characterized with two material functions which describe separately slipping in the boundary layer and shearing within the fluid. Only if both functions are known, correct predictions of flow processes are possible. A simple equipment to separate the shear function and the slip function is described.List of symbols Y^* apparent shear rate - Y _w ^* apparent wall shear rate - Y_w wall shear rate corrected with Rabinowitsch and Weissenberg correction - Y_s reduced shear rate (slip corrected) - Y_ws reduced wall shear rate (slip corrected) - * (r) velocity distribution in a capillary - _G slip velocity (at the wall) - * (r) velocity distribution in a capillary (without slip) - shear stress - _w wall shear stress - V_S total volume rate - V_G shear volume rate - V_G slip volume rate - p ₁ pressure in the reservoir channel of the capillary rheometer - p ₀ athmospheric pressure - L capillary length - R capillary radius     

Effects of Reynolds number on physiological‐type pulsatile flows in a pipe with ring‐type constrictions

The effects of Reynolds number on the physiological‐type of laminar pulsatile flow fields within the vicinity of mechanical ring‐type constriction in small pipes were studied numerically. The parameters considered are: the Reynolds number (Re) in the range of 50–1500; Strouhal number (St) in the range of 0.00156–3.98; Womersley number (Nw) from 0.0 to 50.0. The pulsatile flows considered were physiological‐type of simulated flows. Within a pulsating cycle, detailed flow characteristics were studied through the pulsating contours of streamline (ψ), vorticity (Ω), shear stress (τ) and isobar. The relations between the instantaneous flow rate (Q) and instantaneous pressure gradients (dp/dz) are observed to be elliptic. The relations between the instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations were observed between the instantaneous flow rate (Q) and the maximum velocity, maximum vorticity and maximum shear stress. The Reynolds number of the flow in a pulsating cycle was found to have significant effects on the recirculation length and the pressure gradient within the pulsatile flow regime. Copyright © 1999 John Wiley & Sons, Ltd.     

Edge fracture in cone-plate and parallel plate flows

Edge fracture is an instability of cone-plate and parallel plate flows of viscoelastic liquids and suspensions, characterised by the formation of a `crack' or indentation at a critical shear rate on the free surface of the liquid. A study is undertaken of the theoretical, experimental and computational aspects of edge fracture. The Tanner-Keentok theory of edge fracture in second-order liquids is re-examined and is approximately extended to cover the Criminale-Ericksen-Filbey (CEF) model. The second-order theory shows that the stress distribution on the semi-circular crack is not constant, requiring an average to be taken of the stress; this affects the proportionality constant, K in the edge fracture equation −N _2c = KΓ/a, where N _2c is the critical second normal stress difference, Γ is the surface tension coefficient and a is the fracture diameter. When the minimum stress is used, K = 2/3 as found by Tanner and Keentok (1983). Consideration is given to the sources of experimental error, including secondary flow and slip (wall effect). The effect of inertia on edge fracture is derived. A video camera was used to record the inception and development of edge fracture in four viscoelastic liquids and two suspensions. The recorded image was then measured to obtain the fracture diameter. The edge fracture phenomenon was examined to find its dependence on the physical dimensions of the flow (i.e. parallel plate gap or cone angle), on the surface tension coefficient, on the critical shear rate and on the critical second normal stress difference. The critical second normal stress difference was found to depend on the surface tension coefficient and the fracture diameter, as shown by the theory of Tanner and Keentok (1983); however, the experimental data were best fitted by the equation −N _2c = 1.095Γ/a. It was found that edge fracture in viscoelastic liquids depends on the Reynolds number, which is in good agreement with the inertial theory of edge fracture. Edge fracture in lubricating grease and toothpaste is broadly consistent with the CEF model of edge fracture. A finite volume method program was used to simulate the flow of a viscoelastic liquid, obeying the modified Phan-Thien-Tanner model, to obtain the velocity and stress distribution in parallel plate flow in three dimensions. Stress concentrations of the second normal stress difference (N ₂) were found in the plane of the crack; the velocity distribution shows a secondary flow tending to aid crack formation if N ₂ is negative, and a secondary flow tending to suppress crack formation if N ₂ is positive. Received: 4 January 1999 Accepted: 19 May 1999     

Simulation of turbulent flows around a circular cylinder using nonlinear eddy-viscosity modelling: steady and oscillatory ambient flows

Steady incident flow past a circular cylinder for sub- to supercritical Reynolds number has been simulated as an unsteady Reynolds-averaged Navier–Stokes (RANS) equation problem using nonlinear eddy-viscosity modelling assuming two-dimensional flow. The model of Craft et al. (Int. J. Heat Fluid Flow 17 (1996) 108), with adjustment of the coefficients of the ‘cubic’ terms, predicts the drag crisis at a Reynolds number of about 2×10⁵ due to the onset of turbulence upstream of separation and associated changes in Strouhal number and separation positions. Slightly above this value, at critical Reynolds numbers, drag is overestimated because attached separation bubbles are not simulated. These do not occur at supercritical Reynolds numbers and drag coefficient, Strouhal number and separation positions are in approximate agreement with experimental measurements (which show considerable scatter). Fluctuating lift predictions are similar to sectional values measured experimentally for subcritical Reynolds numbers but corresponding measurements have not been made at supercritical Reynolds numbers. For oscillatory ambient flow, in-line forces, as defined by drag and inertia coefficients, have been compared with the experimental values of Sarpkaya (J. Fluid Mech. 165 (1986) 61) for values of the frequency parameter, β=D²/νT, equal to 1035 and 11240 and Keulegan–Carpenter numbers, KC=U₀T/D, between 0.2 and 15 (D is cylinder diameter, ν is kinematic viscosity, T is oscillation period, and U₀ is the amplitude of oscillating velocity). Variations with KC are qualitatively reproduced and magnitudes show best agreement when there is separation with a large-scale wake, for which the turbulence model is intended. Lift coefficients, frequency and transverse vortex shedding patterns for β=1035 are consistent with available experimental information for β≈250−500. For β=11240, it is predicted that separation is delayed due to more prominent turbulence effects, reducing drag and lift coefficients and causing the wake to be more in line with the flow direction than transverse to it. While these oscillatory flows are highly complex, attached separation bubbles are unlikely and the flows probably two dimensional.