The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Flow over a cylinder with a hinged-splitter plate

Previous work on rigid splitter plates in the wake of a bluff body has shown that the primary vortex shedding can be suppressed for sufficiently long splitter plates. In the present work, we study the problem of a hinged-splitter plate in the wake of a circular cylinder. The splitter plate can rotate about the hinge at the base of the cylinder due to the unsteady fluid forces acting on it, and hence the communication between the two sides of the wake is not totally disrupted as in the rigid splitter plate case. In our study, we investigate this problem in the limit where the stiffness and internal damping associated with the hinge are negligible, and the mass ratio of the splitter plate is small. The experiments show that the splitter plate oscillations increase with Reynolds numbers at low values of Re, and are found to reach a saturation amplitude level at higher Re, Re>4000. This type of saturation amplitude level that appears to continue indefinitely with Re, appears to be related to the fact that there is no structural restoring force, and has been observed previously for transversely oscillating cylinders with no restorin force. In the present case, the saturation tip amplitude level can be tip to 0.45D, where D is the cylinder diameter. For this hinged-rigid splitter plate case, it is found that the splitter plate length to cylinder diameter ratio (L/D) is crucial in determining the character and magnitude of the oscillations. For small splitter plate length (L/D <= 3.0), the oscillations appear to be nearly periodic with tip amplitudes of about 0.45D nearly independent of L/D. The nondiinensional oscillation frequencies (fD/U) on the other hand are found to continuously vary with L/D from fD/U approximate to 0.2 at L/D = 1 to fD/U approximate to 0.1 at L/D = 3. As the splitter plate length is further increased beyond L/D >= 4.0, the character of the splitter plate oscillations suddenly changes. The oscillations become aperiodic with much smaller amplitudes. In this long splitter plate regime, the spectra of the oscillations become broadband, and are reminiscent of the change in character of the wake oscillations seen in the earlier fixed-rigid splitter plate case for L/D >= 5.0. In the present case of the hinged-splitter plate, the sudden transition seen as the splitter plate length (L/D) is increased from 3 to 4 may be attributed to the fact that the wake vortices are no longer able to synchronize with the plate motions for larger splitter plate lengths. Hence, as observed in other vortex-induced vibration problems, the oscillations becomeaperiodic and the amplitude reduces dramatically.

Thermal stresses in a finite solid cylinder due to an axisymmetric temperature field at the end surface.

A three-dimensional rigorous solution for determining thermal stresses in a finite solid cylinder due to a steady state axisymmetric temperature field over one of its end surfaces is given. Numerical results for a solid cylinder having a length to diameter ratio equal to one and subjected to a symmetric temperature variation over half the radius of the cylinder at the end surfaces are included. These results have been compared with the results of the approximate solution given by W. Nowacki.

On the Role of Higher-Order Evanescent Modes in End-Offset Inlet and End-Centered Outlet Elliptical Flow-Reversal Chamber Mufflers

This paper deals with the role of the higher-order evanescent modes generated at the area discontinuities in the acoustic attenuation characteristics of an elliptical end-chamber muffler with an end-offset inlet and end-centered outlet. It has been observed that with an increase in length, the muffler undergoes a transition from being acoustically short to acoustically long. Short end chambers and long end chambers are characterized by transverse plane waves and axial plane waves, respectively, in the low-frequency range. The nondimensional frequency limit k(0)(D-1/2) or k(0)R(0) as well as the chamber length to inlet/outlet pipe diameter ratio, i.e., L/d(0), up to which the muffler behaves like a short chamber and the corresponding limit beyond which the muffler is acoustically long are determined. The limits between which neither the transverse plane-wave model nor the conventional axial plane-wave model gives a satisfactory prediction have also been determined, the region being called the intermediate range. The end-correction expression for this muffler configuration in the acoustically long limit has been obtained using 3-D FEA carried on commercial software, covering most of the dimension range used in the design exercise. Development of a method of combining the transverse plane wave model with the axial plane wave model using the impedance Z] matrix is another noteworthy contribution of this work.

Synthesis of ZnO nanorods on a flexible Phynox alloy substrate: influence of growth temperature on their properties

A novel flexible alloy substrate (Phynox, 50 mm thick) was used for the synthesis of zinc oxide (ZnO) nanorods via a low-temperature solution growth method. The growth of ZnO nanorods was observed over a low temperature range of 60-90 degrees C for a growth duration of 4 hours. The as-synthesized nanorods were characterized using field-emission scanning electron microscopy (FE-SEM), X-ray diffraction (XRD), transmission electron microscopy (TEM), and X-ray photoelectron spectroscopy (XPS) for their morphology, crystallinity, microstructure and composition. The as-grown ZnO nanorods were observed to be relatively vertical to the substrate. However, the morphology of the ZnO nanorods in terms of their length, diameter and aspect ratio was found to vary with the growth temperature. The morphological variation was mainly due to the effects of the various relative growth rates observed at the different growth temperatures. The growth temperature influenced ZnO nanorods were also analyzed for their wetting (either hydrophobic or hydrophilic) properties. After carrying out multiple wetting behaviour analyses, it has been found that the as-synthesized ZnO nanorods are hydrophobic in nature. The ZnO nanorods have potential application possibilities in self-cleaning devices, sensors and actuators as well as energy harvesters such as nanogenerators.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

The problem of liquid droplet combustion-reexamination

The simple quasi-steady analysis of the combustion of a liquid fuel droplet in an oxidising atmosphere provides unsatisfactory explanations for several experimental observations. It's prediction of values for the burning constant (K), the flame-to-droplet diameter ratio ( ) and the flame temperature (Tf) have been found to be amgibuous if not completely inaccurate. A critical survey of the literature has led us to a detailed examination of the effects of unsteadiness and variable properties. The work published to date indicates that the gas-phase unsteadiness is relatively short and therefore quite insignificant.A new theoretical analysis based on heat transfer within the droplet is presented here. It shows that the condensed-phase unsteadiness lasts for about 20â??25% of the total burning time. It is concluded that the discrepancies between experimental observations and the predictions of the constant-property quasi-steady analysis cannot be attributed either to gas-phase or condensed-phase unsteadiness.An analytical model of quasi-steady droplet combustion with variable thermodynamic and transport properties and non-unity Lewis numbers will be examined. Further findings reveal a significant improvement in the prediction of combustion parameters, particularly of K, when consideration is given to variations of cp and λ with the temperature and concentrations of several species. Tf is accurately predicted when the required conditions of incomplete combustion or low ( ) at the flame are met. Further refinement through realistic Lewis numbers predicts ( ) meaningfully.

Thermal stresses in a finite hollow cylinder due to an axisymmetric temperature field at the end surface

A three-dimensional exact solution for determining the thermal stresses in a finite hollow cylinder subject to a steady state axisymmetric temperature field over one of its end surfaces has been given. Numerical results for a hollow cylinder, having lenght to outer diameter ratio equal to one and inner to outer diameter ratio equal to 0.75, subjected to a symmetric temperature variation over the end surfaces of the cylinder have been given.

Transient Momentum and Enthalpy Transfer in Packed Beds at High Reynolds Numbers

An experimental study for transient temperature response and pressure drop in a randomly packed bed at high Reynolds numbers is presented.The packed bed is used as a compact heat exchanger along with a solid-propellant gas generator, to generate room-temperature gases for use in control actuation, air bottle pressurization, etc. Packed beds of lengths 200 and 300 mm were characterized for packing-sphere-based Reynolds numbers ranging from 0.8 x 10(4) to 8.5 x 10(4).The solid packing used in the bed consisted of phi 9.5 mm steel spheres. The bed-to-particle diameter ratio was with the average packed-bed porosity around 0.43. The inlet flow temperature was unsteady and a mesh of spheres was used at either end to eliminate flow entrance and exit effects. Gas temperature and pressure were measured at the entry, exit,and at three axial locations along centerline in the packed beds. The solid packing temperature was measured at three axial locations in the packed bed. A correlation based on the ratio of pressure drop and inlet-flow momentum (Euler number) exhibited an asymptotically decreasing trend with increasing Reynolds number. Axial conduction across the packed bed was found to he negligible in the investigated Reynolds number range. The enthalpy absorption rate to solid packing from hot gases is plotted as a function of a nondimensional time constant for different Reynolds numbers. A longer packed bed had high enthalpy absorption rate at Reynolds number similar to 10(4), which decreased at Reynolds number similar to 10(5). The enthalpy absorption plots can be used for estimating enthalpy drop across packed bed with different material, but for a geometrically similar packing.

Communication: Unusual dynamics of hybrid nanoparticles and their binary mixtures

We present the results on the evolution of microscopic dynamics of hybrid nanoparticles and their binary mixtures as a function of temperature and wave vector. We find unexpectedly a nonmonotonic dependence of the structural relaxation time of the nanoparticles as a function of the morphology. In binary mixtures of two of the largest nanoparticles studied, we observe re-entrant vitrification as a function of the volume fraction of the smaller nanoparticle, which is unusual for such high diameter ratio. Possible explanation for the observed behavior is provided. (C) 2010 American Institute of Physics. doi:10.1063/1.3495480]

Fatigue crack propagation model for plain concrete - An analogy with population growth

In this work, an analytical model is proposed for fatigue crack propagation in plain concrete based on population growth exponential law and in conjunction with principles of dimensional analysis and self-similarity. This model takes into account parameters such as loading history, fracture toughness, crack length, loading ratio and structural size. The predicted results are compared with experimental crack growth data for constant and variable amplitude loading and are found to capture the size effect apart from showing a good agreement. Using this model, a sensitivity analysis is carried out to study the effect of various parameters that influence fatigue failure. (C) 2010 Elsevier Ltd. All rights reserved.

Mean and fluctuating pressure in boat-tail separated flows at transonic speeds

Experiments were carried out investigating the features of mean and unsteady surface pressure fluctuations in boat-tail separated flows relevant to launch vehicle configurations at transonic speeds. The tests were performed on a generic axisymmetric body in the Mach-number range of 0.7-1.2, and the important geometrical parameters, namely, the boat-tail angle and diameter ratio, were varied systematically. The measurements made included primarily the mean and unsteady surface-pressure fluctuations on nine different model configurations. Flow-visualization studies employing a surface oil flow, and schlieren techniques were carried out to infer features like boundary-layer separation, reattachment, and shock waves in the flow. The features of mean and fluctuating surface pressures are discussed in detail including aspects of similarity. It has been observed that, on a generic configuration employed in the present study, the maximum levels of surface-pressure fluctuations in the reattachment zone are appreciably lower than those found on launch vehicle configurations having a bulbous or hammerhead nose shape. A simple correlation is suggested for the maximum value of rms pressure fluctuations in the reattachment zone at different freestream Mach numbers.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

Effective thermal conductivity of a heat generating rod bundle dissipating heat by natural convection and radiation

A numerical study of conjugate natural convection and surface radiation in a horizontal hexagonal sheath housing 19 solid heat generating rods with cladding and argon as the fill gas, is performed. The natural convection in the sheath is driven by the volumetric heat generation in the solid rods. The problem is solved using the FLUENT CFD code. A correlation is obtained to predict the maximum temperature in the rod bundle for different pitch-to-diameter ratios and heat generating rates. The effective thermal conductivity is related to the heat generation rate, maximum temperature and the sheath temperature. Results are presented for the dimensionless maximum temperature, Rayleigh number and the contribution of radiation with changing emissivity, total wattage and the pitch-to-diameter ratio. In the simulation of a larger system that contains a rod bundle, the effective thermal conductivity facilitates simplified modelling of the rod bundle by treating it as a solid of effective thermal conductivity. The parametric studies revealed that the contribution of radiation can be 38-65% of the total heat generation, for the parameter ranges chosen. Data for critical Rayleigh number above which natural convection comes into effect is also presented. (C) 2011 Elsevier B.V. All rights reserved.