Flow and heat transfer in a second grade fluid over a stretching sheet

An analysis is carried out to study the flow and heat transfer characteristics in a second grade fluid over a stretching sheet with prescribed surface temperature including the effects of frictional heating, internal heat generation or absorption, and work due to deformation. In order to solve the fourth-order non-linear differential equation, associated with the flow problem, a fourth boundary condition is augmented and a proper sign for the normal stress modulus is used. It is observed that for a physical flow problem the solution is unique. The solutions for the temperature and the heat transfer characteristics are obtained numerically and presented by a table and graphs. Furthermore, it is shown that the heat flow is always from the stretching sheet to the fluid.     

Modeling of transverse self-oscillations of a circular cylinder in an incompressible fluid flow in a plane channel with circulation

Experimentally observed self-oscillations of a cylinder in a plane channel whose width is slightly greater than the cylinder diameter under the impact of the incoming fluid flow are modeled. Within the model of a nonseparated potential flow around the cylinder, the coefficients of added mass of the cylinder are calculated with the help of the generalized method of images. When the cylinder touches the channel wall, the circulation sign changes, and its value is determined by the boundary element method and the no-slip condition for the fluid at the contact point. The Joukowski force and the drag force proportional to the square of velocity are taken into account in the equations of motion of cylinder.     

A cantilever conveying fluid: coherent modes versus beam modes

A semi-analytical approach to obtain the proper orthogonal modes is described for the non-linear oscillation of a cantilevered pipe conveying fluid. Theoretically, while the spatial coherent structures are the eigenfunctions of the time-averaged spatial autocorrelation functions, it emerges that once the Galerkin projection of the proper orthogonal modes is carried out using the uniform cantilever-beam modes, the spatial dependency of the integral eigenvalue problem can be eliminated by analytical manipulation which avoids any spatial discretization error. As the solution of the integral equation is obtained semi-analytically by linearly projecting the proper orthogonal modes on the cantilever-beam modes, any linear or non-linear operation can be carried out analytically on the proper orthogonal modes. Furthermore, the reduced-order eigenvalue problem minimizes the numerical pollution which leads to spurious eigenvectors, as may arise in the case of a large-scale eigenvalue problem (without the Galerkin projection of the eigenvectors on the cantilever-beam modes). This methodology can conveniently be used to study the convergence of the numerically calculated proper orthogonal modes obtained from the full-scale eigenvalue problem.     

Field-aligned flow of a conducting fluid past a source of magnetism

Summary A prescribed source of magnetism moves at constant speed through a viscous conducting incompressible fluid with an aligned uniform magnetic field. The velocity and magnetic fields induced at a distance from the source are calculated. The induced fields are also calculated for the case in which the applied field is absent. Although no special symmetry or alignment is assumed, the source is ideal in the sense that enclosures (wires or magnets) are infinitesimal in at least two dimensions. Dynamical interactions will occur in a viscous fluid and their effect in the far field is estimated.As a consequence of finite conductivity and viscosity, the usual wakes are present which trail or lead the source depending upon the sign of (1–A ²), where A is the ratio of the source speed to the Alfvén speed in the undisturbed fluid. Outside the wake the total perturbation magnetic field due to the source is the static field plus a monopole field, divided by (1–A ²).An estimate is also made of the rate at which energy is dissipated as a consequence of viscous interactions and ohmic heating throughout the fluid, outside the immediate vicinity of the source.Geo-Astrophysics Laboratory.Plasma Physics Laboratory.     

A note on the drag for fluids of grade three

An expression is obtained for the traction vector t on a solid surface which is adjacent to an incompressible fluid of grade three which is compatible with thermodynamics. It is found that unlike fluids of grade two wherein there is no additional drag due to the non-Newtonian nature of the fluid for bodies with certain geometric symmetries (e.g. sphere), fluids of grade three provide an additional drag which is of the same sign as that provided by the viscous terms, provided certain symmetry conditions are met by the velocity field.     

Spatial local solutions of the Navier-Stokes equations

This paper considers solutions of the Navier-Stokes equations polynomial in the coordinates, which. are called local solutions. For an incompressible fluid, all higher-order terms (sums of higher-order. monomials) of degree 2 are found and it is proved that nontrivial axisymmetric higher-order terms. of degree higher than 2 do not exist. Nonsolenoidal axisymmetric solutions are listed, which can be. treated as steady-state barotropic gas flows in a potential external-force field. All elliptic vortices. generalizing the well-known Kirchhoff solution are calculated. All solutions of degree 3 with the. higher-order term of partial form are found. Some of these solutions break down in a finite time. regardless of the value and sign of viscosity. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 109–119, March–April, 2009.     

Studies on the law of laminar helical flow of power-law fluid in eccentric annuli

According to the principle of fluid mechanics,the law of laminar,helical flow of power-law fluid in eccentric annuli is studied extensively in this paper.The apparent viscosity,velocities distribution of laminar helical flow of power-law fluid are discussed and calculating methods of flow rate and pressure loss are presented.The factors influencing flow are also analysed.On the basis of theoretical studies some new results of the present paper are compared with the theories of the helical flow of the power-law fluid in concentric annuli.The test verified that theoretical formulas in this article are proper and general.They can be used for calculating hydraulic parameters in drilling engineering.     

Heat transfer and entropy generation analysis of non-Newtonian fluid flow through vertical microchannel with convective boundary condition

The entropy generation and heat transfer characteristics of magnetohydrodynamic(MHD) third-grade fluid flow through a vertical porous microchannel with a convective boundary condition are analyzed. Entropy generation due to flow of MHD non-Newtonian third-grade fluid within a microchannel and temperature-dependent viscosity is studied using the entropy generation rate and Vogel's model. The equations describing flow and heat transport along with boundary conditions are first made dimensionless using proper non-dimensional transformations and then solved numerically via the finite element method(FEM). An appropriate comparison is made with the previously published results in the literature as a limiting case of the considered problem.The comparison confirms excellent agreement. The effects of the Grashof number, the Hartmann number, the Biot number, the exponential space-and thermal-dependent heat source(ESHS/THS) parameters, and the viscous dissipation parameter on the temperature and velocity are studied and presented graphically. The entropy generation and the Bejan number are also calculated. From the comprehensive parametric study, it is recognized that the production of entropy can be improved with convective heating and viscous dissipation aspects. It is also found that the ESHS aspect dominates the THS aspect.     

A domain perturbation study of steady flow in a cone-and-plate rheometer of non-ideal geometry

The method of domain perturbation developed by Joseph is used to calculate velocity and stress profiles in a slightly misaligned cone-and-plate rheometer where the cone is spinning and the plate is stationary. Results for a Newtonian fluid, a Criminale-Ericksen-Filbey fluid, an upper-convected Maxwell fluid, and a White-Metzner fluid are presented and compared with earlier results in which the cone is stationary and the plate is spinning (Dudgeon and Wedgewood, 1993). Streamlines calculated for the Newtonian fluid show a very small recirculation region near the stationary plate. Velocity and stress contours are symmetric around the plane of largest gap width. For the elastic fluids studied, streamlines are asymmetric. The fluid response lags where the fluid is dominated by memory effects. Much larger recirculation regions are calculated for fluids dominated by shear thinning. These recirculation regions contain a large fraction of the fluid in the apparatus and have the effect of changing the shape of the flow domain for the remaining fluid that rotates around the cone's axis. Elasticity also has a pronounced effect on the stress profile, indicating that the accuracy of the cone and plate may be compromised even for small mis-alignments.     

Announcements

Recent numerical investigations on pressure surges during pump trip in pumping installations showed that by including an air entrainment variable wave speed model, reasonable predictions of transient pressure surges with proper phasing and attenuation of pressure peaks can be obtained. These calculated results are consistent with similar field measurements made with the pumps operating at low pump cut-out levels, when air entrainment due to an attached surface vortex was observed. However, in the numerical calculation procedures it is assumed that the inertia of the moving elements of the check valve is small and that the check valve closes at zero reverse flow velocity. In practice, check valves seldom close precisely at zero reverse flow velocity. With the check valves not closing at zero reverse velocity, the present numerical computations show that the air content in a fluid system can adversely affect the check valve performance. With the fluid system operating within a critical range of air entrainment values, the present analysis showed that there is a possibility of ‘check valve slamming’ when the check valves are selected based only on the analysis of an air-free system. This phenomenon is confirmed through field observations.     

A finite element analysis of laminar unsteady flow of viscoelastic fluids through channels with non-uniform cross-sections

The effects of non-Newtonian behaviour of a fluid and unsteadiness on flow in a channel with non-uniform cross-section have been investigated. The rheological behaviour of the fluid is assumed to be described by the constitutive equation of a viscoelastic fluid obeying the Oldroyd-B model. The finite element method is used to analyse the flow. The novel features of the present method are the adoption of the velocity correction technique for the momentum equations and of the two-step explicit scheme for the extra stress equations. This approach makes the computational scheme simple in algorithmic structure, which therefore implies that the present technique is capable of handling large-scale problems. The scheme is completed by the introduction of balancing tensor diffusivity (wherever necessary) in the momentum equations. It is important to mention that the proper boundary condition for pressure (at the outlet) has been developed to solve the pressure Poisson equation, and then the results for velocity, pressure and extra stress fields have been computed for different values of the Weissenberg number, viscosity due to elasticity, etc. Finally, it is pertinent to point out that the present numerical scheme, along with the proper boundary condition for pressure developed here, demonstrates its versatility and suitability for analysing the unsteady flow of viscoelastic fluid through a channel with non-uniform cross-section.     

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

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

Where is the rudder of a fish？： the mechanism of swimming and control of self-propelled fish school

Numerical simulation and control of self- propelled swimming of two- and three-dimensional biomimetic fish school in a viscous flow are investigated. With a parallel computational fluid dynamics package for the two- and three-dimensional moving boundary problem, which combines the adaptive multi-grid finite volume method and the methods of immersed boundary and volume of fluid, it is found that due to the interactions of vortices in the wakes, without proper control, a fish school swim with a given flap- ping rule can not keep the fixed shape of a queue. In order to understand the secret of fish swimming, a new feedback con- trol strategy of fish motion is proposed for the first time, i,e., the locomotion speed is adjusted by the flapping frequency of the caudal, and the direction of swimming is controlled by the swinging of the head of a fish. Results show that with this feedback control strategy, a fish school can keep the good order of a queue in cruising, turning or swimming around circles. This new control strategy, which separates the speed control and direction control, is important in the construction of biomimetic robot fish, with which it greatly simplifies the control devices of a biomimetic robot fish.     

Numerical simulation of chemically reactive Powell-Eyring liquid flow with double diffusive Cattaneo-Christov heat and mass flux theories

numerical study is reported for two-dimensional flow of an incompressible Powell-Eyring fluid by stretching the surface with the Cattaneo-Christov model of heat diffusion. Impacts of heat generation/absorption and destructive/generative chemical reactions are considered. Use of proper variables leads to a system of non-linear dimensionless expressions. Velocity, temperature and concentration profiles are achieved through a finite difference based algorithm with a successive over-relaxation (SOR) method. Emerging dimensionless quantities are described with graphs and tables. The temperature and concentration profiles decay due to enhancement in fluid parameters and Deborah numbers.     

Experimental and numerical study on flow characteristics and heat transfer of an oscillating jet in a channel

A fluidic oscillator can produce self-induced and self-sustaining oscillating jet by fluid supply without moving parts. This device has attracted research interest in heat and mass transfer enhancement in recent years. In the current study, a double-feedback fluidic oscillator was numerically investigated based on three-dimensional unsteady Reynolds-averaged Navier-Stokes equations (3D-URANS) while the operating fluid is an incompressible flow. Then, the results were validated with experimental data by two-dimensional time-resolved particle image velocimetry (2D-TR-PIV) and thermographic phosphor thermometry (TPT) for the velocity and temperature field, respectively. A grid sensitivity study was done by comparison of instantaneous and time-averaged flow fields. Additionally, the proper orthogonal decomposition (POD) method was used to find the phase information of the oscillating jet, and fast Fourier transform (FFT) analysis was used to find the frequency of the oscillating jet to validate the numerical results. The effect of the working fluid was also studied. Finally, in order to determine the effect of the Reynolds number on heat transfer enhancement, the Q-criterion was calculated to provide detailed insight into the oscillating mechanism. The results show that the non-dimensional frequency of oscillation is independent of either the working fluid or mass flow rate. Additionally, for a given fluid, increasing Re causes strong vortices and increases the frequency of oscillation. However, the convection heat transfer did not change significantly when varying the mass flow rate because the convection velocity of vortices increases as the mass flow rate is enhanced. A comparison with a free jet reveals that the oscillating jet in a channel is useful in terms of covering a larger area.     

Calculation of the decay rate of oscillations of a low-viscosity fluid rotating in a cylindrical vessel

A study is made of the problem of small oscillations of a low-viscosity fluid rotating under conditions of weightlessness in a cylindrical vessel. For fixed volume of the fluid and different values of the angular velocity and the wetting angle the decay rate and frequency of the characteristic oscillations are calculated. For an ideal fluid, the shapes of the characteristic oscillations of the free surface of the fluid are also calculated.     

Bandgap calculation for mixed in-plane waves in 2D phononic crystals based on Dirichlet-to-Neumann map

In this paper, a method based on the Dirichlet-to-Neumann map is developed for bandgap calculation of mixed in-plane waves propagating in 2D phononic crystals with square and triangular lattices. The method expresses the scattered fields in a unit cell as the cylindrical wave expansions and imposes the Bloch condition on the boundary of the unit cell. The Dirichlet-to-Neumann (DtN) map is applied to obtain a linear eigenvalue equation, from which the Bloch wave vectors along the irreducible Brillouin zone are calculated for a given frequency. Compared with other methods, the present method is memory-saving and time-saving. It can yield accurate results with fast convergence for various material combinations including those with large acous- tic mismatch without extra computational cost. The method is also efficient for mixed fluid-solid systems because it considers the different wave modes in the fluid and solid as well as the proper fluid-solid interface condition.     

Three-dimensional flow structure in highly buoyant jet by scanning stereo PIV combined with POD analysis

The flow characteristics and the structure of highly buoyant jet of low density fluid issuing into a stagnant surrounding of high density fluid is studied by scanning stereo PIV combined with proper orthogonal decomposition (POD) analysis. The experiment is carried out at Froude number of 0.3 and Reynolds number of 200, which satisfies the inflow condition due to the unstable density gradient near the nozzle exit. An increase in the maximum mean velocity occurs and the vertical velocity fluctuation is highly amplified near the nozzle exit, which suggests the influence of inflow due to the unstable density gradient. The POD analysis indicates that the vertical velocity fluctuation is the major source of fluctuating energy contributing to the development of the highly buoyant jet. The examination of the POD modes show that the longitudinal structure of the vertical velocity fluctuation is generated along the jet axis having the opposite sign of velocity fluctuation on both sides of the jet axis. The vertical scale of the POD mode decreases with increasing the mode number and results in the frequent appearance of cross-flow across the buoyant jet. The reconstruction flow from the POD modes indicates that the vortex structure is caused by the highly sheared layer between the upward and downward velocity and the inflow is induced by the vortex structure. The magnitude of the vortex structure seems to be weakened with an increase in the distance from the nozzle and the buoyant jet approaches to an asymptotic state in the further downstream.     

The combined effects of inlet fluid flow and temperature nonuniformity in cross flow plate-fin compact heat exchanger using finite element method

An analysis of a crossflow plate-fin heat exchanger accouning for the combined effects of inlet fluid flow nonuniformity and temperature nonuniformity on both hot and cold fluid sides is carried out using a Finite Element Model. A mathematical equation is developed to generate different types of fluid flow/temperature maldistribution models considering the possible deviations in inlet fluid flow. Using these fluid flow maldistribution models, the exchanger effectiveness and its deteriorations due to flow/temperature nonuniformity are calculated for entire range of design and operating conditions. It was found that the performance deteriorations are quite significant in some typical applications due to inlet fluid flow/temperature nonuniformity.