Adaption and use of a compressible flow code for turbomachinery design

The present paper presents and discusses the use of dierent codes regarding the numerical simulation of a radial-in ow turbine. A radial-in ow turbine test case was selected from published literature [1] and commercial codes (Fluent and CFX) were used to perform the steady-state numerical simulations. An in-house compressible- ow simulation code, Eilmer3 [2] was also adapted in order to make it suitable to perform turbomachinery simulations and preliminary results are presented and discussed. The code itself as well as its adaptation, comprising the addition of terms for the rotating frame of reference, programmable boundary conditions for periodic boundaries and a mixing plane interface between the rotating and non-rotating blocks are also discussed. Several cases with dierent orders of complexity in terms of geometry were considered and the results were compared across the dierent codes. The agreement between these results and published data is also discussed.

Steady natural convection of non-Newtonian power-law fluid in a trapezoidal enclosure

Numerically investigation of free convection heat transfer in a differentially heated trapezoidal cavity filled with non-Newtonian Power-law fluid has been performed in this study. The left inclined surface is uniformly heated whereas the right inclined surface is maintained as uniformly cooled. The top and bottom surfaces are kept adiabatic with initially quiescent fluid inside the enclosure. Finite volume based commercial software FLUENT 14.5 is used to solve the governing equations. Dependency of various flow parameters of fluid flow and heat transfer is analyzed including Rayleigh number, Ra ranging from 10^5 to 10^7, Prandtl number, Pr of 100 to 10,000 and power index, n of 0.6 to 1.4. Outcomes have been reported in terms of isotherms, streamline, and local Nusselt number for various Ra, Pr, n and inclined angles. Grid sensitivity analysis is performed and numerically obtained results have been compared with those results available in the literature and found good agreement.

Effect of hydrodynamics and surface roughness on the electrochemical behaviour of carbon steel in CSG produced water

The influence of fluid flow, surface roughness and immersion time on the electrochemical behaviour of carbon steel in coal seam gas produced water under static and hydrodynamic conditions has been studied. The disc electrode surface morphology before and after the corrosion test was characterized using scanning electron microscopy (SEM). The corrosion product was examined using X-ray photoelectron spectroscopy (XPS) and X-ray diffractometry (XRD).The results show that the anodic current density increased with increasing surface roughness and consequently a decrease in corrosion surface resistance. Under dynamic flow conditions, the corrosion rate increased with increasing rotating speed due to the high mass transfer coefficient and formation of non-protective akaganeite β- FeO(OH) and goethite α- FeO(OH) corrosion scale at the electrode surface.The corrosion rate was lowest at 0 rpm.The corrosion rate decreased in both static and dynamic conditions with increasing immersion time. The decrease in corrosion rate is attributed to the deposition of corrosion products on the electrode surface. SEM results revealed that the rougher surface exhibited a great tendency toward pitting corrosion.

Flow evaluation and hemolysis analysis of BVAD centrifugal blood pump by computational fluids dynamics

Computational fluid dynamics (CFD) and particle image velocimetry (PIV) are commonly used techniques to evaluate the flow characteristics in the development stage of blood pumps. CFD technique allows rapid change to pump parameters to optimize the pump performance without having to construct a costly prototype model. These techniques are used in the construction of a bi-ventricular assist device (BVAD) which combines the functions of LVAD and RVAD in a compact unit. The BVAD construction consists of two separate chambers with similar impellers, volutes, inlet and output sections. To achieve the required flow characteristics of an average flow rate of 5 l/min and different pressure heads (left – 100mmHg and right – 20mmHg), the impellers were set at different rotating speeds. From the CFD results, a six-blade impeller design was adopted for the development of the BVAD. It was also observed that the fluid can flow smoothly through the pump with minimum shear stress and area of stagnation which are related to haemolysis and thrombosis. Based on the compatible Reynolds number the flow through the model was calculated for the left and the right pumps. As it was not possible to have both the left and right chambers in the experimental model, the left and right pumps were tested separately.

Puzzles of excited charm meson masses

We attempt a comprehensive analysis of the low lying charm meson states which present several puzzles, including the poor determination of masses of several non-strange excited mesons. We use the well-determined masses of the ground states and the strange first excited states to 'predict' the mass of the non-strange first excited state in the framework of heavy hadron chiral perturbation theory, an approach that is complementary to the well-known analysis of Mehen and Springer. This approach points to values for the masses of these states that are smaller than the experimental determinations. We provide a critical assessment of these mass measurements and point out the need for new experimental information. (c) 2007 Elsevier B.V. All rights reserved.

Terahertz wave characteristics of a single-walled carbon nanotube containing a fluid flow using the nonlocal Timoshenko beam model

This paper presents the effect of nonlocal scaling parameter on the terahertz wave propagation in fluid filled single walled carbon nanotubes (SWCNTs). The SWCNT is modeled as a Timoshenko beam,including rotary inertia and transverse shear deformation by considering the nonlocal scale effects. A uniform fluid velocity of 1000 m/s is assumed. The analysis shows that, for a fluid filled SWCNT, the wavenumbers of flexural and shear waves will increase and the corresponding wave speeds will decrease as compared to an empty SWCNT. The nonlocal scale parameter introduces certain band gap region in both flexural and shear wave mode where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite (or wave speed tends to zero). The frequency at which this phenomenon occurs is called the ``escape frequency''. The effect of fluid density on the terahertz wave propagation in SWCNT is also studied and the analysis shows that as the fluid becomes denser, the wave speeds will decrease. The escape frequency decreases with increase in nonlocal scaling parameter, for both wave modes. We also show that the effect of fluid density and velocity are negligible on the escape frequencies of flexural and shear wave modes. (C) 2010 Elsevier B.V. All rights reserved.

Flow of a power law fluid in an annulus with porous walls

The steady flow of a power law fluid in annuli with porous walls is investigated. The solution for the axial velocity component is obtained as a power series in terms of the cross flow Reynolds number, the first term of the series giving the solution for the case of the solid wall annulus. The cross flow is restricted to be such that the rate of injection of fluid at one wall of the annulus is equal to the rate of suction at the other wall and also we have considered only very small values of the cross flow velocity. The velocity profiles are drawn for different values of n and for different gaps and the results are discussed in detail. The behaviour of the average flux, in different eases is also discussed.

Breakage of a drop of inviscid fluid due to a pressure fluctuation at its surface

In this work, an attempt is made to gain a better understanding of the breakage of low-viscosity drops in turbulent flows by determining the dynamics of deformation of an inviscid drop in response to a pressure variation acting on the drop surface. Known scaling relationships between wavenumbers and frequencies, and between pressure fluctuations and velocity fluctuations in the inertial subrange are used in characterizing the pressure fluctuation. The existence of a maximum stable drop diameter d(max) follows once scaling laws of turbulent flow are used to correlate the magnitude of the disruptive forces with the duration for which they act. Two undetermined dimensionless quantities, both of order unity, appear in the equations of continuity, motion, and the boundary conditions in terms of pressure fluctuations applied on the surface. One is a constant of proportionality relating root-mean-square values of pressure and velocity differences between two points separated by a distance l. The other is a Weber number based on turbulent stresses acting on the drop and the resisting stresses in the drop due to interfacial tension. The former is set equal to 1, and the latter is determined by studying the interaction of a drop of diameter equal to d(max) with a pressure fluctuation of length scale equal to the drop diameter. The model is then used to study the breakage of drops of diameter greater than d(max) and those with densities different from that of the suspending fluid. It is found that, at least during breakage of a drop of diameter greater than d(max) by interaction with a fluctuation of equal length scale, a satellite drop is always formed between two larger drops. When very large drops are broken by smaller-length-scale fluctuations, highly deformed shapes are produced suggesting the possibility of further fragmentation due to instabilities. The model predicts that as the dispersed-phase density increases, d(max) decreases.

The sputtering temperature of a cooling cylindrical rod without and with an insulated core in a two-fluid medium

Utilising Jones' method associated with the Wiener-Hopf technique, explicit solutions are obtained for the temperature distributions on the surface of a cylindrical rod without an insulated core as well as that inside a cylindrical rod with an insulated inner core when the rod, in either of the two cases, is allowed to enter, with a uniform speed, into two different layers of fluid with different cooling abilities. Simple expressions are derived for the values of the sputtering temperatures of the rod at the points of entry into the respective layers, assuming the upper layer of the fluid to be of finite depth and the lower of infinite extent. Both the problems are solved through a three-part Wiener-Hopf problem of special type and the numerical results under certain special circumstances are obtained and presented in tabular forms.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Flow of a stratified fluid in a wavy channel

The flow of a stratified fluid in a channel with small and large deformations is investigated. The analogy of this flow with swirling flow in tubes with non-uniform cross-sections is studied. The flow near the wall is blocked when the Froude number takes certain critical values. The possibility of preventing the stagnation zones in the flow field is also discussed

A note on the damping and oscillations of a fluid drop moving in another fluid

The oscillations of a drop moving in another fluid medium have been studied at low values of Reynolds number and Weber number by taking into consideration the shape of the drop and the viscosities of the two phases in addition to the interfacial tension. The deformation of the drop modifies the Lamb's expression for frequency by including a correction term while the viscous effects split the frequency into a pair of frequencies—one lower and the other higher than Lamb's. The lower frequency mode has ample experimental support while the higher frequency mode has also been observed. The two modes almost merge with Lamb's frequency for the asymptotic cases of a drop in free space or a bubble in a dense viscous fluid but the splitting becomes large when the two fluids have similar properties. Instead of oscillations, aperiodic damping modes are found to occur in drops with sizes smaller than a critical size ($\sim\hat{\rho}\hat{\nu}^2/T $). With the help of these calculations, many of the available experimental results are analyzed and discussed.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Unsteady MUD rotating flow over a rotating sphere near the equator

The unsteady rotating flow of a laminar incompressible viscous electrically conducting fluid over a rotating sphere in the vicinity of the equator has been studied. The fluid and the body rotate either in the same direction or in opposite directions. The effects of surface suction and magnetic field have been included in the analysis. There is an initial steady state that is perturbed by a sudden change in the rotational velocity of the sphere, and this causes unsteadiness in the flow field. The nonlinear coupled parabolic partial differential equations governing the boundary-layer flow have been solved numerically by using an implicit finite-difference scheme. For large suction or magnetic field, analytical solutions have also been obtained. The magnitude of the radial, meridional and rotational velocity components is found to be higher when the fluid and the body rotate in opposite directions than when they rotate in the same direction. The surface shear stresses in the meridional and rotational directions change sign when the ratio of the angular velocities of the sphere and the fluid lambda greater than or equal to lambda(0). The final (new) steady state is reached rather quickly which implies that the spin-up time is small. The magnetic field and surface suction reduce the meridional shear stress, but increase the surface shear stress in the rotational direction.