Analytical series solution for the fully developed forced convection duct flow with frictional heating and variable viscosity

Laminar forced convection flow of a liquid in the fully developed region of a circular duct with isothermal wall is analyzed. The effects of viscous dissipation as well as of temperature dependent viscosity are taken into account. The coupled momentum and energy equations are solved analytically by means of a power series method. Then, reference is made to the Poiseuille model for the temperature change of viscosity. For a fixed value of the axial pressure gradient along the duct, dual solutions are found for the velocity and temperature fields. Although dual solutions correspond to the same value of the axial pressure gradient, they lead in general to different values of the average fluid velocity, of the average fluid temperature and of the wall heat flux. It is shown that, for a given fluid and for a fixed duct radius, the absolute value of the axial pressure gradient has an upper bound above which no steady laminar solution can exist.     

Thermal radiation effect on a mixed convection flow and heat transfer of the Williamson fluid past an exponentially shrinking permeable sheet with a convective boundary condition

The thermal radiation effect on a steady mixed convective flow with heat transfer of a nonlinear (non-Newtonian) Williamson fluid past an exponentially shrinking porous sheet with a convective boundary condition is investigated numerically. In this study, both an assisting flow and an opposing flow are considered. The governing equations are converted into nonlinear ordinary differential equations by using a suitable transformation. A numerical solution of the problem is obtained by using the Matlab software package for different values of the governing parameters. The results show that dual nonsimilar solutions exist for the opposing flow, whereas the solution for the assisting flow is unique. It is also observed that the dual nonsimilar solutions exist only if a certain amount of mass suction is applied through the porous sheet, which depends on the Williamson parameter, convective parameter, and radiation parameter.     

New exact solutions of Stokes' second problem for an MHD second grade fluid in a porous space

We investigate a problem describing the oscillating flow of an incompressible magnetohydrodynamic (MHD) second grade fluid in a porous half space. Exact solutions for sine and cosine oscillations are developed by applying the Laplace transform method. The total obtained solution is a sum of steady and transient solutions. Particular attention is given to the effects of magnetic and porous medium parameters on the velocity. It is shown that previous results for a non-porous medium and hydrodynamic fluid are the limiting cases of the present problem. The results for velocity are plotted and discussed carefully.     

Dual steady solutions in natural convection between horizontal concentric cylinders

Dual steady solutions in natural convection in an annulus between two horizontal concentric cylinders are numerically investigated for a fluid of Prandtl number 0.7. It is found that, when the Rayleigh number based on the gap width exceeds a certain critical value, dual steady two-dimensional (2-D) flows can be realized: one being the crescent-shaped eddy flow commonly observed and the other the flow consisting of two counter-rotating eddies and their mirror images. The critical Rayleigh number decreases as the inverse relative gap width increases.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

Boundary-layer flow of Reiner–Philippoff fluids past a stretching wedge

This paper presents a numerical analysis of the steady boundary-layer flow of a Reiner–Philippoff fluid induced by a 90° stretching wedge in a variable free stream. The governing partial differential equations are converted into a set of two ordinary differential equations by the use of a similarity transformation. The flow is therefore governed by a stretching velocity parameter λ and two non-Newtonian fluid parameters γ and μ₀. The variation of the skin friction, as well as other flow characteristics, as a function of the governing parameters is presented graphically and tabulated. A stability analysis has also been performed for this self-similar flow based on linear disturbances to the steady similarity solutions. The results presented in this paper reveal that there are no multiple (dual) solutions for the present problem and the unique solution is stable.     

Effects of suction or blowing on the velocity and temperature distribution in the flow past a porous flat plate of a power-law fluid

An analysis is made of the steady flow of a non-Newtonian fluid past an infinite porous flat plate subject to suction or blowing. The incompressible fluid obeys Ostwald-de Waele power-law model. It is shown that steady solutions for velocity distribution exist only for a pseudoplastic (shear-thinning) fluid for which the power-law index n satisfies 0<n<1 provided that there is suction at the plate. Velocity at a point is found to increase with increase in n. No steady solution for velocity distribution exists when there is blowing at the plate. The solution of the energy equation governing temperature distribution in the flow of a pseudoplastic fluid past an infinite porous plate subject to uniform suction reveals that temperature at a given point near the plate increases with n but further away, temperature decreases with increase in n. A novel result of the analysis is that both the skin-friction and the heat flux at the plate are independent of n.     

Micropolar fluid flow over a shrinking sheet

An analysis is carried out for the steady two-dimensional flow of a micropolar fluid over a shrinking sheet in its own plane. The shrinking velocity is assumed to vary linearly with the distance from a fixed point on the sheet. The features of the flow and heat transfer characteristics are analyzed and discussed. It is found that the solution exists only if adequate suction through the permeable sheet is introduced. Moreover, stronger suction is necessary for the solution to exist for a micropolar fluid compared to a classical Newtonian fluid. Dual solutions are obtained for certain suction and material parameters.     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

Resistance of rough plates in a rotating fluid

Generalizing Navier’s partial slip condition, the flow due to a rough or striated plate moving in a rotating fluid is studied. It is found that the motion of the plate, the fluid surface velocity, and the shear stress are in general not in the same direction. The solution is extended to the case of finite depth, or Couette slip flow in a rotating system. In this case an optimum depth for minimum drag is found. The solutions are also closed form exact solutions of the Navier–Stokes equations. The results are fundamental to flows with Coriolis effects.     

MHD flow of a viscous fluid on a nonlinear porous shrinking sheet with homotopy analysis method

The present paper investigates the magnetohydrodynamic（MHD） flow of a viscous fluid towards a nonlinear porous shrinking sheet.The governing equations are simplified by similarity transformations.The reduced problem is then solved by the homotopy analysis method.The pertinent parameters appearing in the problem are discussed graphically and presented in tables.It is found that the shrinking solutions exist in the presence of MHD.It is also observed from the tables that the solutions for f（0） with different values of parameters are convergent.     

Hall effects on the unsteady hydromagnetic oscillatory flow of a second-grade fluid

An exact solution of an oscillatory flow is constructed in a rotating fluid under the influence of an uniform transverse magnetic field. The fluid is considered as second-grade (non-Newtonian). The influence of Hall currents and material parameters of the second-grade fluid is investigated. The hydromagnetic flow is generated in the uniformly rotating fluid bounded between two rigid non-conducting parallel plates by small amplitude oscillations of the upper plate. The exact solutions of the steady and unsteady velocity fields are constructed. It is found that the steady solution depends on the Hall parameter but is independent of the material parameter of the fluid. The unsteady part of the solution depends upon both (Hall and material) parameters. Attention is focused upon the physical nature of the solution, and the structure of the various kinds of boundary layers is examined. Several results of physical interest have been deduced in limiting cases.     

Four-Component Gas/Water/Oil Displacements in One Dimension: Part III,Development of Miscibility

The structure of the system of conservation laws in fully compositional three-phase, four-component flow is examined for the first time. Two of the eigenvalues can be found analytically for this flow regime, regardless of the equation of state used to model the phase behavior. A cubic equation of state is used to calculate the gas/oil phase behavior and Henry’s law is used to represent the partitioning of hydrocarbons between the hydrocarbon and water phases. Sample analytical solutions are found for Riemann problems modeling the injection of carbon-dioxide and water into an oil reservoir. Finally, the structure of the solutions at the minimum miscibility pressure (MMP) for the hydrocarbon system is studied. We show that when water is injected simultaneously with gas at the MMP, multicontact miscible displacement of the oil by the injected gas develops only if the fraction of water in the injected fluid is below a critical value. For water fractions above the critical value, water flowing at a high velocity forces the composition path of the solution to remain in the three-phase region in which two hydrocarbon phases are present. We study both condensing and vaporizing gas drives to demonstrate that this result is general.     

Do general viscoelastic stresses for the flow of an upper convected Maxwell fluid satisfy the momentum equation?

Given a general velocity field consistent with the stagnation point flow, can the viscoelastic stresses arising in the flow of an upper convected Maxwell fluid found by solving the constitutive equation also satisfy the momentum equation? Consideration is given to the study of the stress tensor arising in the steady flow of an upper convected Maxwell (UCM) fluid with a velocity field consistent with the stagnation point flow. By the method of characteristics, exact solutions to the partial differential equations arising in the approximating model of the viscoelastic stresses in the flow of an upper convected Maxwell (UCM) fluid are obtained for the three components of the stress tensor, for reasonably general velocity fields. We are able to account for the effects of variable boundary data at the inflow by considering the viscoelastic stresses over two spatial variables. Furthermore, we assume a relatively general velocity field. As a special case, some results present in the recent literature are obtained; it is known that these special case solutions do not satisfy the momentum equation. In the general case we consider, we find that the general solution will not satisfy the momentum equation except in a limited restricted case. We discuss how this shortcoming might be rectified by use of a more general velocity field.     

平面粘性流体扰动与哈密顿体系

通过变分原理,将哈密顿体系的理论引入到平面粘性流体扰动的问题中,导出一套哈密顿算子矩阵的本征函数向量展开求解问题的方法。基于直接法求解流体力学基本方程,导出流场一般特征关系,通过本征值的求解及本征向量的叠加,得到波扰动解,继可分析流场端部效应。从而在该领域用在哈密顿体系下辛几何空间中研究问题的方法代替了传统在拉格朗日体系欧几里德空间分析问题的方法。为流体力学的研究提供一条新途径。     

On relaxation times in the Navier-Stokes-Voigt model

We study analytically and numerically the relaxation time of flow evolution governed by the Navier-Stokes-Voigt (NSV) model. We first show that for the Taylor–Green vortex decay problem, NSV admits an exact solution which evolves slower than true fluid flow. Secondly, we show numerically for a channel flow test problem using standard discretisation methods that although NSV provides more regular solutions compared to usual Navier-Stokes solutions, NSV approximations take significantly longer to reach the steady state.     

Mixed convection boundary layer flow near stagnation-point on vertical surface with slip

This paper considers the steady mixed convection boundary layer flow of a viscous and incompressible fluid near the stagnation-point on a vertical surface with the slip effect at the boundary. The temperature of the sheet and the velocity of the external flow are assumed to vary linearly with the distance from the stagnation-point. The governing partial differential equations are first transformed into a system of ordinary differential equations, which are then solved numerically by a shooting method. The features of the flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed. Both assisting and opposing flows are considered. The results indicate that for the opposing flow, the dual solutions exist in a certain range of the buoyancy parameter, while for the assisting flow, the solution is unique. In general, the velocity slip increases the heat transfer rate at the surface, while the thermal slip decreases it.     

Mixed convection boundary-layer flow over a vertical surface embedded in a porous medium

In this paper we have numerically investigated the existence and uniqueness of a vertically flowing fluid passed a model of a thin vertical fin in a saturated porous media. We have assumed the two-dimensional mixed convection from a fin, which is modelled as a fixed, semi-infinite vertical surface, embedded in a fluid-saturated porous media under the boundary-layer approximation. We have taken the temperature, in excess of the constant temperature in the ambient fluid on the fin, to vary as  , where is measured from the leading edge of the plate and λ is a fixed constant. The Rayleigh number is assumed to be large so that the boundary-layer approximation may be made and the fluid velocity at the edge of the boundary-layer is assumed to vary as . The problem then depends on two parameters, namely λ and , the ratio of the Rayleigh to Péclet numbers. It is found that when λ>0 (<0) there are (is) dual (unique) solution(s) when is grater than some negative values of (which depends on λ). When λ<0 there is a range of negative value of (which depends on λ) for which dual solutions exist and for both λ>0 and λ<0 there is a negative value of (which depends on λ) for which there is no solution. Finally, solutions for 0<1 and 1 have been obtained.     

Quantitative NMR velocity imaging of a main-chain liquid crystalline polymer flowing through an abrupt contraction

 The flow of isotropic and liquid crystalline (LC) hydroxypropylcellulose (HPC) aqueous solutions into an abrupt axisymmetric contraction has been quantitatively measured by pulsed-field-gradient NMR techniques. Steady-state axial velocity profiles, acquired upstream of the contraction, reveal a large contraction entry length for the LC solution. This entry flow field exists over an order of magnitude change in flow rate and is attributed to elasticity that is associated with polydomain liquid crystallinity. Pronounced, off-centerline velocity maxima (in an axisymmetric flow field) were present upstream of the contraction, in the entry flow region. Apparently, a more viscous and elastic core of fluid was present along the centerline; this fluid resisted elongational strain more than the fluid closer to the walls. Quantitative velocity profiles were extracted from displacement distributions and corrected for elongational dispersion. The isotropic solution velocity profiles matched those obtained from viscoelastic simulations using an approximate Doi-Edwards model, parameterized with independent rheological data. Received: 29 April 1999/Accepted: 30 August 1999