CFD Analysis of conjugate natural convection in a vertical annulus

Conjugate natural convection in a vertical annulus with a centrally located vertical heat generating rod is studied numerically. The governing equations are discretized on a staggered mesh and are solved using a pressure-correction algorithm. A parametric study is performed by varying the Grashof number, aspect ratio, and the solid-to-fluid thermal conductivity ratio over wide ranges with the Prandtl number fixed at 0.7. Results are presented for the variation of several quantities of interest such as the local Nusselt numbers on the inner and outer boundaries, the axial variation of the centerline and interface temperatures, maximum solid, average solid and average interface temperature variations with Grashof number, and the average Nusselt number variation for the inner and outer boundaries with Grashof number. The average Nusselt number from the conjugate analysis is found to be between the Nusselt numbers of the isothermal and the isoflux cases. The average Nusselt numbers on the inner and outer boundaries show an increasing trend with the Grashof number. Correlations are presented for the Nusselt number and the dimensionless temperatures of interest in terms of the parameters of the problem.

Effect of surface radiation on conjugate natural convection in a horizontal annulus driven by inner heat generating solid cylinder

This paper reports a numerical study of the laminar conjugate natural convection heat transfer with and without the interaction of the surface radiation in a horizontal cylindrical annulus formed between an inner heat generating solid circular cylinder and an outer isothermal circular boundary. Numerical solutions are obtained by solving the governing equations with a pressure correction method on a collocated (non-staggered) mesh. Steady-state results are presented for the flow and temperature distributions and Nusselt numbers for the heat generation based Grashof number ranging from 10(7) to 10(10), solid-to-fluid thermal conductivity ratios of 1, 5, 10, 50 and 100, radius ratios of 0.226 and 0.452 and surface emissivities of 0-0.8 with air as the working medium. It is observed that surface radiation reduces the convective heat transfer in the annulus compared to the pure natural convection case and enhances the overall Nusselt number.

Unsteady natural convection from a horizontal annulus filled with a porous medium

The unsteady natural convection flow from a horizontal cylindrical annulus filled with a non-Darcy porous medium has been studied. The unsteadiness in the problem arises due to the impulsive change in the wall temperature of the outer cylinder. The Navier–Stokes equations along with the energy equation governing the unsteady natural convection flow have been solved by the finite-volume method. The effect of time variation on the heat transfer is more pronounced only in a small time interval immediately after the start of the impulsive motion and the steady state is reached after certain time. The results show that the annulus completely filled with a porous medium has the best insulating effectiveness. Convection in the horizontal annulus is confined mostly at top and bottom regions. Hence, only these regions should be insulated. In case of annulus partially filled with a porous material, insulating the region near the outer cylinder is more effective than insulating the region near the inner cylinder. The effect of Darcy number on the heat transfer is more pronounced than that of the Grashof number.

Steady flow of a maxwell fluid in a porous cylindrical annulus with circular cross-section

When a fluid with memory is injected into any flow region some assumptions regarding the initial state of stress have to be made in order to determine the state of stress at any subsequent instant. For a Maxwell fluid, it is assumed that the fluid near the surface of injection is suddenly stressed and responds by starting flow in accordance with the mechanical model chosen. The flow of a Maxwell fluid with a single relaxation time has been determined under the above assumption in the following two cases: (i) annulus between two porous concentric circular cylinders, and (ii) space between two porous and infinitely extending parallel plates. The nature of flow in the present case is similar to that of the Reiner-Rivlin fluids obtained by Narasimhan2).

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.

Exact solution for the unsteady motion of a viscous-fluid in a porous annulus

An exact solution to the problem of time-dependent motion of a viscous fluid in an annulus with porous walls is obtained under the assumption that the rate of suction at one wall is equal to the rate of injection at the other. Finite Hankel transform is used to obtain a closed-form solution for the axial velocity. The average axial velocity profiles are depicted graphically.

Non-Boussinesq conjugate natural convection in a vertical annulus

Non-Boussinesq conjugate natural convection in a vertical annulus with a centrally located vertical heat generating rod is studied numerically, taking into account variable transport properties. Results are presented for maximum solid temperatures, average Nusselt numbers and average pressure. In general, the Boussinesq model predicts higher temperatures in the solid and lower average Nusselt numbers on the inner and outer boundaries. (C) 2010 Elsevier Ltd. All rights reserved.

Studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring

Numerical and experimental studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring are performed. The alloy is electromagnetically stirred to produce semisolid slurry in a cylindrical graphite mould placed in the annulus of a linear electromagnetic stirrer. The mould is cooled at the bottom, such that solidification progresses from the bottom to the top of the cylindrical mould. A numerical model is developed for simulating the transport phenomena associated with the solidification process using a set of single-phase governing equations of mass. momentum, energy. and species conservation. The viscosity variation of the slurry, used in the model, is determined experimentally using a rotary viscometer. The set of governing equations is solved using a pressure-based finite volume technique, along with an enthalpy based phase change algorithm. The numerical study involves prediction of temperature, velocity, species and solid fraction distribution in the mould. Corresponding solidification experiments are performed, with time-temperature history recorded at key locations. The microstructures at various temperature measurement locations in the solidified billet are analyzed. The numerical predictions of temperature variations are in good agreement with experiments, and the predicted flow field evolution correlates well with the microstructures observed at various locations.

Structure and assembly of Sesbania mosaic virus

Sesbania mosaic virus (SeMV) is a ss-RNA (4149 nt) plant sobemovirus isolated from farmer's field around Tirupathi, Andhra Pradesh. The viral capsid (30 nm diameter) consists of 180 copies of protein subunits (MW 29 kDa) organized with icosahedral symmetry. In order to understand the mechanism of assembly of SeMV, a large number of deletion and substitution mutants of the coat protein (CP) were constructed. Recombinant SeMV CP (rCP) as well as the N-terminal rCP deletion mutant Delta N22 were found to assemble in E. coli into virus-like particles (VLPs). Delta N36 and Delta N65 mostly formed smaller particles consisting of 60 protein subunits. Although particlem assembly was not affected due to the substitution of aspartates (D14 and D149) that coordinate calcium ions by asparagines, the stability of the resulting capsids was drastically reduced. Deletion of residues forming a characteristic beta-annulus at the icosahedral 3-folds did not affect the assembly of VLPs. Mutation of a single tryptophan, which occurs near the icosahedral fivefold axis to glutamate or lysine, resulted in the disruption of the capsid leading to soluble dimers that resembled the quasi-dimer structure of the native virus. Replacement of positively charged residues in the amino terminal segment of CP resulted in the formation of empty shells. Based on these observations, a plausible mechanism of assembly is proposed.

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.

Determination of solid fraction evolution during solidification of aluminum alloy in presence of linear electromagnetic stirring

In many industrial casting processes, knowledge of the solid fraction evolution during the solidification process is a key factor in determining the process parameters such as cooling rate, stirring intensity and in estimating the total solidification time. In the present work, a new method of estimating solid fraction is presented, which is based on calorimetric principles. In this method, the cooling curve data at each point in the melt, along with the thermal boundary conditions, are used to perform energy balance in the mould, from which solid fraction generation during any time interval can be estimated. This method is applied to the case of a rheocasting process, in which Al-Si alloy (A356 alloy) is solidified by stirring in a cylindrical mould placed in the annulus of a linear electromagnetic stirrer. The metal in the mould is simultaneously cooled and stirred to produce a cylindrical billet with non-dendritic globular microstructure. Temperature is measured at key locations in the mould to assess the various heat exchange processes prevalent in the mould and to monitor the solidification rate. The results obtained by energy balance method are compared with those by the conventional procedure of calculating solid fraction using the Schiel equation.

Performance Evaluation of Distance-Hop Proportionality on Geometric Graph Models of Dense Sensor Networks

Wireless sensor networks can often be viewed in terms of a uniform deployment of a large number of nodes on a region in Euclidean space, e.g., the unit square. After deployment, the nodes self-organise into a mesh topology. In a dense, homogeneous deployment, a frequently used approximation is to take the hop distance between nodes to be proportional to the Euclidean distance between them. In this paper, we analyse the performance of this approximation. We show that nodes with a certain hop distance from a fixed anchor node lie within a certain annulus with probability approach- ing unity as the number of nodes n → ∞. We take a uniform, i.i.d. deployment of n nodes on a unit square, and consider the geometric graph on these nodes with radius r(n) = c q ln n n . We show that, for a given hop distance h of a node from a fixed anchor on the unit square,the Euclidean distance lies within [(1−ǫ)(h−1)r(n), hr(n)],for ǫ > 0, with probability approaching unity as n → ∞.This result shows that it is more likely to expect a node, with hop distance h from the anchor, to lie within this an- nulus centred at the anchor location, and of width roughly r(n), rather than close to a circle whose radius is exactly proportional to h. We show that if the radius r of the ge- ometric graph is fixed, the convergence of the probability is exponentially fast. Similar results hold for a randomised lattice deployment. We provide simulation results that il- lustrate the theory, and serve to show how large n needs to be for the asymptotics to be useful.

EXPERIMENTAL INVESTIGATION OF UNSTEADY FLOW IN A TRANSONIC UNI-STAGE AXIAL COMPRESSOR

Detailed steady and unsteady experimental measurements and analysis were performed on a Single stage Transonic Axial Compressor with asymmetric rotor tip clearance for studying the compressor stall phenomena. The installed compressor had asymmetric tip clearance around the rotor casing varying from about 0.65mm to 1.25mm. A calibrated 5-hole aerodynamic probe was traversed radially at exit of rotor and showed the characteristics of increased flow angle at lower mass flow rates for all the speeds. Mach number distribution and boundary layer effects were also clearly captured. Unsteady measurements for velocity were carried out to study the stall cell behavior using a single component calibrated hotwire probe oriented in axial and tangential directions for choke/free flow and near stall conditions. The hotwire probe was traversed radially across the annulus at inlet to the compressor and showed that the velocity fluctuations were dissimilar when probe was aligned axial and tangential to the flow. Averaged velocities across the annulus showed the reduction in velocity as stall was approached. Axial mean flow velocity decreased across the annulus for all the speeds investigated. Tangential velocity at free flow condition was higher at the tip region due to larger radius. At stall condition, the tangential velocity showed decreased velocities at the tip and slightly increased velocities at the hub section indicating that the flow has breakdown at the tip region of the blade and fluid is accelerated below the blockage zone. The averaged turbulent intensity in axial and tangential flow directions increased from free flow to stall condition for all compressor rated speeds. Fast Fourier Transform (FFT) of the raw signals at stall flow condition showed stall cell and its corresponding frequency of occurrence. The stalling frequency of about half of rotational speed of the rotor along with large tip clearance suggests that modal type stall inception was occurring.