Stagnation flow towards a shrinking sheet

The stagnation flow towards a shrinking sheet is studied. A similarity transform reduces the Navier-Stokes equations to a set of non-linear ordinary differential equations which are then integrated numerically. Both two-dimensional and axisymmetric stagnation flows are considered. It is found that solutions do not exist for larger shrinking rates and may be non-unique in the two-dimensional case. The non-alignment of the stagnation flow and the shrinking sheet complicates the flow structure. Convective heat transfer decreases with the shrinking rate due to an increase in boundary layer thickness.     

Direct numerical simulation of a liquid sheet in a compressible gas stream in axisymmetric and planar configurations

A thin liquid sheet present in the shear layer of a compressible gas jet is investigated using an Eulerian approach with mixed-fluid treatment for the governing equations describing the gas–liquid two-phase flow system, where the gas is treated as fully compressible and the liquid as incompressible. The effects of different topological configurations, surface tension, gas pressure and liquid sheet thickness on the flow development of the gas–liquid two-phase flow system have been examined by direct solution of the compressible Navier–Stokes equations using highly accurate numerical schemes. The interface dynamics are captured using volume of fluid and continuum surface force models. The simulations show that the dispersion of the liquid sheet is dominated by vortical structures formed at the jet shear layer due to the Kelvin–Helmholtz instability. The axisymmetric case is less vortical than its planar counterpart that exhibits formation of larger vortical structures and larger liquid dispersion. It has been identified that the vorticity development and the liquid dispersion in a planar configuration are increased at the absence of surface tension, which when present, tends to oppose the development of the Kelvin–Helmholtz instability. An opposite trend was observed for an axisymmetric configuration where surface tension tends to promote the development of vorticity. An increase in vorticity development and liquid dispersion was observed for increased liquid sheet thickness, while a decreasing trend was observed for higher gas pressure. Therefore surface tension, liquid sheet thickness and gas pressure factors all affect the flow vorticity which consequently affects the dispersion of the liquid.      

Fluid flow in an arbitrary cascade on an axisymmetric stream surface in a variable thickness layer

A solution is given for the problem of flow past a cascade on an axisymmetric stream surface in a layer of variable thickness, which is a component part of the approximate solution of the three-dimensional problem for a three-dimensional cascade. Generalized analytic functions are used to obtain the integral equation for the potential function, which is solved via iteration method by reduction to a system of linear algebraic equations. An algorithm and a program for the Minsk-2 computer are formulated. The precision of the algorithm is evaluated and results are presented of the calculation of an example cascade.In the formulation of [1, 3] the problem of flow past a three-dimensional turbomachine cascade is reduced approximately to the joint solution of two-dimensional problems of the averaged axisymmetric flow and the flow on an axisymmetric stream surface in an elementary layer of variable thickness.In the following we solve the second problem for an arbitrary cascade with finite thickness rotating with constant angular velocity in ideal fluid flow: the solution is carried out on a Minsk-2 computer.Many studies have been devoted to this problem. A method for solving the direct problem for a cascade of flat plates in a hyperbolic layer was presented in [2]. Methods were developed in [1, 3] for constructing the flow for the case of a channel with variable thickness; these methods are approximately applicable for dense cascades but yield considerable error for small-load turbomachine cascades. The solution developed in [4], somewhat reminiscent of that of [2], is applicable for thin, slightly curved profiles in a layer with monotonically varying thickness. A solution has been given for a circular cascade for layers varying logarithmically [5] and linearly [6]. Approximate methods for slightly curved profiles in a monotonically varying layer with account for layer variability only in the discharge component were examined in [7–9]. A solution is given in [10] for an arbitrary layer by means of the relaxation method, which yields a roughly approximate flow pattern. The general solution of the problem by means of potential theory and the method of singularities presented in [11] is in error because of neglect of the crossflow through the skeletal line. The computer solution of [12] contains an unassessed error for the calculations in an arbitrary layer. The finite difference method is used in [13] to solve the differential equation of flow, which is illustrated by numerical examples for monotonie layers of axial turbomachines. The numerical solution of [13] is very complex.The solution presented below is found in the general formulation with respect to the geometric parameters of the cascade and the axisymmetric surface and also in terms of the layer thickness variation law.The numerical solution requires about 15 minutes of machine time on the Minsk-2 computer.     

On axisymmetric flow of Sisko fluid over a radially stretching sheet

This study presents an analysis of the axisymmetric flow of a non-Newtonian fluid over a radially stretching sheet. The momentum equations for two-dimensional flow are first modeled for Sisko fluid constitutive model, which is a combination of power-law and Newtonian fluids. The general momentum equations are then simplified by invoking the boundary layer analysis. Then a non-linear ordinary differential equation governing the axisymmetric boundary layer flow of Sisko fluid over a radially stretching sheet is obtained by introducing new suitable similarity transformations. The resulting non-linear ordinary differential equation is solved analytically via the homotopy analysis method (HAM). Closed form exact solution is then also obtained for the cases n=0 and 1. Analytical results are presented for the velocity profiles for some values of governing parameters such as power-law index, material parameter and stretching parameter. In addition, the local skin friction coefficient for several sets of the values of physical parameter is tabulated and analyzed. It is shown that the results presented in this study for the axisymmetric flow over a radially non-linear stretching sheet of Sisko fluid are quite general so that the corresponding results for the Newtonian fluid and the power-law fluid can be obtained as two limiting cases.     

具有热源项的驻点流动和传热问题的近似解析解

研究持续拉伸变形表面上二维平面和轴对称驻点流的动量和热量传输问题。利用同伦分析方法获得速度分布和温度分布的级数解,讨论了级数的敛散性。通过图形分析主流速度与拉伸速度的比率参数,普朗特参数,热源参数和流动类型指标对速度边界层和温度边界层的影响。结果表明,这些参数对二维平面驻点流动和传热有较大的影响。     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.     

Analytic solution for axisymmetric flow over a nonlinearly stretching sheet

An analysis is performed for the boundary-layer flow of a viscous fluid over a nonlinear axisymmetric stretching sheet. By introducing new nonlinear similarity transformations, the partial differential equations governing the flow are reduced to an ordinary differential equation. The resulting ordinary differential equation is solved using the homotopy analysis method (HAM). Analytic solution is given in the form of an infinite series. Convergence of the obtained series solution is explicitly established. The solution for an axisymmetric linear stretching sheet is obtained as a special case.     

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     

Stresses distributions in a rotating functionally graded piezoelectric hollow cylinder

An analytic solution to the axisymmetric problem of a long, radially polarized, hollow cylinder composed of functionally graded piezoelectric material (FGPM) rotating about its axis at a constant angular velocity is presented. For the case that electric, thermal and mechanical properties of the material obey different power laws in the thickness direction, distributions for radial displacement, stresses and electric potential in the FGPM hollow cylinder are determined by using the theory of electrothermoelasticity. Some useful discussions and numerical examples are presented to show the significant influence of material nonhomogeneity, and adopting suitable graded indexes and applying suitable geometric size and rotating velocity ω may optimize the rotating FGPM hollow cylindrical structures. This will be of particular importance in modern engineering application.     

Axisymmetric flow and heat transfer of the Carreau fluid due to a radially stretching sheet: Numerical study

The prime objective of this article is to study the axisymmetric flow and heat transfer of the Carreau fluid over a radially stretching sheet. The Carreau constitutive model is used to discuss the characteristics of both shear-thinning and shear-thickening fluids. The momentum equations for the two-dimensional flow field are first modeled for the Carreau fluid with the aid of the boundary layer approximations. The essential equations of the problem are reduced to a set of nonlinear ordinary differential equations by using local similarity transformations. Numerical solutions of the governing differential equations are obtained for the velocity and temperature fields by using the fifth-order Runge–Kutta method along with the shooting technique. These solutions are obtained for various values of physical parameters. The results indicate substantial reduction of the flow velocity as well as the thermal boundary layer thickness for the shear-thinning fluid with an increase in the Weissenberg number, and the opposite behavior is noted for the shear-thickening fluid. Numerical results are validated by comparisons with already published results.     

On the spreading radius of surface tension driven oil on deep water

The spread of a thin oil film by surface tension gradients from an oil source of unlimited mass on deep water is considered. A similarity solution for the velocity fields of the oil and water, the oil thickness and the rate at which each grow is obtained both for axisymmetric and the previously explored planar spreading. The dimensionless size of the spread, which is oil type independent, is shown to be 1.0754 and 1.4150 for axisymmetric and planar spreading respectively. It is further shown that the oil film equation of state, which relates surface tension to oil thickness, is unique to each oil or oil-surfactant mixture.     

Thermo-magnetic analysis of thick-walled spherical pressure vessels made of functionally graded materials

This study presents an analytical solution of thermal and mechanical displacements, strains, and stresses for a thick-walled rotating spherical pressure vessel made of functionally graded materials(FGMs). The pressure vessel is subject to axisymmetric mechanical and thermal loadings within a uniform magnetic field. The material properties of the FGM are considered as the power-law distribution along the thickness. Navier's equation, which is a second-order ordinary differential equation, is derived from the mechanical equilibrium equation with the consideration of the thermal stresses and the Lorentz force resulting from the magnetic field. The distributions of the displacement,strains, and stresses are determined by the exact solution to Navier's equation. Numerical results clarify the influence of the thermal loading, magnetic field, non-homogeneity constant, internal pressure, and angular velocity on the magneto-thermo-elastic response of the functionally graded spherical vessel. It is observed that these parameters have remarkable effects on the distributions of radial displacement, radial and circumferential strains, and radial and circumferential stresses.     

Liquid film flow due to an unsteady stretching sheet

We have studied two-dimensional flow of a thin liquid film over an impulsively stretching sheet under assumption of uniform initial film thickness. Using singular perturbation technique both momentum and film evolution equations are solved analytically for small Reynolds number and these solutions are verified numerically. Numerical computation for large Reynolds number shows an anomalous behaviour of film thinning rate in different time zone. These results are explained physically and the crucial role-played by viscosity in this case is highlighted. It is found that faster rate of thinning can be obtained if the sheet is stretched impulsively with continuously increasing stretching speed.     

Heat transfer in the flow of a cold,two-dimensional draining sheet over a hot,horizontal cylinder

The paper considers heat transfer characteristics of thin film flow over a hot horizontal cylinder resulting from a cold vertical sheet of liquid falling onto the surface. The underlying physical features of the developing film thickness, velocity and temperature distributions have been illustrated by numerical solutions of high accuracy for large Reynolds numbers using the modified Keller box method. The solutions for film thickness distribution are good agreement with those obtained using the Pohlhausen integral momentum technique thus providing a basic confirmation of the validity of the results presented.     

Extended homotopy perturbation method and the axisymmetric flow past a porous stretching sheet

The extended homotopy perturbation method, which is an extension of the celebrated homotopy perturbation method (HPM), is applied to obtain a solution to the problem of the steady, laminar, axisymmetric flow of a viscous, incompressible fluid past a porous stretching sheet. The solution so obtained is totally analytical and is expressible in terms of the cross‐flow velocity of the fluid past the stretching sheet. Its hallmark is that it does not depend upon computation of any auxiliary parameter for enlarging the convergence region of the solution. Rather, it calculates the solution automatically adjusting the scaling factor of the independent similarity variable normal to the sheet. The results obtained by the extended HPM are in excellent agreement with the exact numerical solution. Also, an asymptotic solution valid for large suction parameter is developed, which matches well with the exact solution even for moderate values of the suction parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

EFFECTS OF A TRANSLATING LOAD ON A FLOATING PLATE—STRUCTURAL DRAG AND PLATE DEFORMATION

The elastic deformation of a structural plate floating on water caused by a translating three-dimensional load is investigated. The problem is akin to the landing and take-off of aircraft on a structural or ice sheet. The initial-boundary-value problem is solved analytically using a free-surface condition that incorporates the flexural rigidity of the plate. The three-dimensional load is modeled as an axisymmetric, translating pressure distribution. The time-dependent analytical solution is used to obtain the unsteady drag of this moving pressure, if it exists, as well as its asymptotic behavior at large time. The behavior of the transition of the drag near a critical speed related to the minimum celerity of the free waves of the hydroelastic system is examined. Asymptotic analysis shows that the drag attains a discontinuous but finite value as the translation speed approaches the critical speed, an essential difference from some existing two-dimensional results. The growth rate of the plate slope is found to be weakly singular, like log t, for large time. Comparisons with published experimental data for plate deformation are made for the case of an ice sheet. The agreement is very favorable. Implications on the operation of floating runways are discussed.     

传热传质对有抽吸的收缩薄片上的非线性磁流体动力学边界层流动的影响

This work is concerned with Magnetohydrodynamic viscous flow due to a shrinking sheet in the presence of suction. The cases of two dimensional and axisymmetric shrinking are discussed. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are numerically solved by using an advanced numeric technique. Favorability comparisons with previously published work are presented. Numerical results for the dimensionless velocity, temperature and concentration profiles as well as for the skin friction, heat and mass transfer and deposition rate are obtained and displayed graphically for pertinent parameters to show interesting aspects of the solution.     

Non-Isotropic Length Scales During the Compression Stroke of a Motored Piston Engine

For the case of axial compression the two-point velocity correlation equations of axisymmetric homogeneous turbulence are derived. Appropriate integrations then lead to equations for the components of the Reynolds stress tensor as well as to those for the two independent integral length-scales characterizing axisymmetric homogeneous turbulence. These equations contain a certain number of empirical constants. Values for these constants are taken from the literature, or were adjusted from the present data.The resulting model is validated using data from a motored piston engine. The flow field, which has negligible swirl and tumble, has been measured using particle image velocimetry (PIV). Since turbulence is axisymmetric and homogeneous in the counter region, two-dimensional PIV provides the time history of the axial and radial length-scales. The experimental data are compared with the mathematical model.     

Some properties of multi-fluid stokes flows

The solution of Stokes' equations for a rotating axisymmetric body which possesses reflection symmetry about a planar interface between two infinite immiscible quiescent viscous fluids is shown to be independent of the viscosities of the fluids and identical with the solution when the fluids have the same viscosity. The result is generalized to a rotating axisymmetric system of bodies which possesses reflection symmetry about each interface of a plane stratified system of fluids. An analogous result for two-fluid systems with a nonplanar static interface is also derived. The effect on torque reduction produced by the presence of a second fluid layer adjacent to a rotating axisymmetric body is considered and explicit calculations are given for the case of a sphere. A proof of uniqueness for unbounded multi-fluid Stokes' flow is given and the asymptotic far field structure of the velocity field is determined for axisymmetric flow caused by the rotation of axisymmetric bodies.