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.

Fuzzy waste load allocation model: a multiobjective approach

Fuzzy Waste Load Allocation Model (FWLAM), developed in an earlier study, derives the optimal fractional levels, for the base flow conditions, considering the goals of the Pollution Control Agency (PCA) and dischargers. The Modified Fuzzy Waste Load Allocation Model (MFWLAM) developed subsequently is a stochastic model and considers the moments (mean, variance and skewness) of water quality indicators, incorporating uncertainty due to randomness of input variables along with uncertainty due to imprecision. The risk of low water quality is reduced significantly by using this modified model, but inclusion of new constraints leads to a low value of acceptability level, A, interpreted as the maximized minimum satisfaction in the system. To improve this value, a new model, which is a combination Of FWLAM and MFWLAM, is presented, allowing for some violations in the constraints of MFWLAM. This combined model is a multiobjective optimization model having the objectives, maximization of acceptability level and minimization of violation of constraints. Fuzzy multiobjective programming, goal programming and fuzzy goal programming are used to find the solutions. For the optimization model, Probabilistic Global Search Lausanne (PGSL) is used as a nonlinear optimization tool. The methodology is applied to a case study of the Tunga-Bhadra river system in south India. The model results in a compromised solution of a higher value of acceptability level as compared to MFWLAM, with a satisfactory value of risk. Thus the goal of risk minimization is achieved with a comparatively better value of acceptability level.

Improved Inverted V-Notch or Chimney Weir

This paper is devoted to the improvement of the measuring range of inverted V-notch (IVN) weir, a practical linear sharp-crested weir, designed earlier by the writers. The range of linearity of IVN can be considerably enhanced (by more than 200%) by the addition of a retangular weir of width 0.265W (W = half crest width) at a depth of 0.735d (d = altitude of IVN), above the crest of the weir, which is equivalent to providing at this depth two vertical straight lines to the IVN, resulting in a chimney-shaped profile; hence, the modified weir is named chimney weir. The design parameters of the weir, that is, the linearity range, base flow depth, and datum constant, which fixes the reference plane of the weir, are estimated by solving the nonlinear programming problem using a numerical optimization procedure. For flows through this weir above a depth of 0.22d, the discharges are proportional to the depth of flow measured above a reference plane situated at 0.08d above the weir crest for all heads in the range 0.22d <= h <= 2.43d, within a maximum percentage deviation of ±1.5 from the theoretical discharge. A significant result of the analysis is that the same linear head-discharge relationship governing the flow through the IVN is also valid for the extended chimney weir. Experiments with three different chimney weirs show excellent agreement with the theory by giving a constant average coefficient of discharge for each weir.

Optimal crop planning and conjunctive use of water resources in a coastal river basin

Due to increasing trend of intensive rice cultivation in a coastal river basin, crop planning and groundwater management are imperative for the sustainable agriculture. For effective management, two models have been developed viz. groundwater balance model and optimum cropping and groundwater management model to determine optimum cropping pattern and groundwater allocation from private and government tubewells according to different soil types (saline and non-saline), type of agriculture (rainfed and irrigated) and seasons (monsoon and winter). A groundwater balance model has been developed considering mass balance approach. The components of the groundwater balance considered are recharge from rainfall, irrigated rice and non-rice fields, base flow from rivers and seepage flow from surface drains. In the second phase, a linear programming optimization model is developed for optimal cropping and groundwater management for maximizing the economic returns. The models developed were applied to a portion of coastal river basin in Orissa State, India and optimal cropping pattern for various scenarios of river flow and groundwater availability was obtained.

Relevance of local parallel theory to the linear stability of laminar separation bubbles

Laminar separation bubbles are thought to be highly non-parallel, and hence global stability studies start from this premise. However, experimentalists have always realized that the flow is more parallel than is commonly believed, for pressure-gradient-induced bubbles, and this is why linear parallel stability theory has been successful in describing their early stages of transition. The present experimental/numerical study re-examines this important issue and finds that the base flow in such a separation bubble becomes nearly parallel due to a strong-interaction process between the separated boundary layer and the outer potential flow. The so-called dead-air region or the region of constant pressure is a simple consequence of this strong interaction. We use triple-deck theory to qualitatively explain these features. Next, the implications of global analysis for the linear stability of separation bubbles are considered. In particular we show that in the initial portion of the bubble, where the flow is nearly parallel, local stability analysis is sufficient to capture the essential physics. It appears that the real utility of the global analysis is perhaps in the rear portion of the bubble, where the flow is highly non-parallel, and where the secondary/nonlinear instability stages are likely to dominate the dynamics.

Analysis of Annular Plug Nozzle Flowfield

The flowfields associated with truncated annular plug nozzles of varying lengths are studied both experimentally and using computational tools. The nozzles are designed to observe wake structure transition for the range of pressure ratios considered. A classification of the open wake regime is proposed for comparing and analyzing the plug flowfields. The three-dimensional relief experienced by the annular plug flow leads to greater wave interactions on the plug surface as compared with linear plug flow, resulting in a delayed transition of the base wake. The Reynolds averaged Navier-Stokes based solvers employed in the studies could predict the plug surface flow accurately, whereas they exhibited limitations with regard to plug base flow predictions. Based on the experimental data generated, an empirical model for predicting closed wake base pressure is proposed and compared with other models available in literature.

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.

Absolute/Convective Instability Transition in a Backward Facing Step Combustor: Fundamental Mechanism and Influence of Density Gradient

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms-(1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are colocated. Regions in the space of parameters characterizing the base flow velocity profile, i.e., shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on understanding prior experimental studies of combustion instability in backward facing step combustors and hydrodynamic instability in other flows such as heated jets and bluff body stabilized flames is discussed.

ABSOLUTE/CONVECTIVE INSTABILITY TRANSITION IN A BACKWARD FACING STEP COMBUSTOR: FUNDAMENTAL MECHANISM AND INFLUENCE OF DENSITY GRADIENT

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms - (1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are co-located. Regions in the space of parameters characterizing the base flow velocity profile, i.e. shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on prior observations of flow instability in other flows such as heated jets and bluff-body stabilized flames is discussed.

Tensile flow behaviour of 2.25Cr-1Mo ferritic steel base metal and simulated heat affected zone structures of 2.25Cr-1Mo weld joint

Tensile tests in the temperature range 298 to 873 K have been performed on 2.25Cr-1Mo base metal and simulated heat affected zone (HAZ) structures of its weld joint, namely coarse grain bainite, fine grain bainite and intercritical structure. Tensile flow behaviour of all the microstructural conditions could be adequately described by the Hollomon equation (sigma = K-1 epsilon(n1)) at higher (> 623 K) temperatures. Deviation from the Hollomon equation was observed at low strains and lower (< 623 K) temperatures. The Ludwigson modification of Hollomon's equation, sigma = K-1 epsilon(n1) + exp (K-2 + n(2) epsilon), was found to describe the flow curve. In general, the flow parameters n(1), K-1, n(2) and K-2 were found to decrease with increase in temperature except in the intermediate temperature range (423 to 623 K). Peaks/plateaus were observed in their variation with temperature in the intermediate temperature range coinciding with the occurrence of serrated flow in the load-elongation curve. The n(1) Value increased and the K-1 value decreased with the type of microstructure in the order: coarse grain bainite, fine grain bainite, base metal and intercritical structure. The variation of nl with microstructure has been rationalized on the basis of mean free path (MFP) of dislocations which is directly related to the inter-particle spacing. Larger MFP of dislocations lead to higher strain hardening exponents n(1).

Effect of base dissipation on the granular flow down an inclined plane

The effect of base dissipation on the granular flow down an inclined plane is examined by altering the coefficient of restitution between the moving and base particles in discrete element (DE) simulations. The interaction laws between two moving particles are kept fixed, and the coefficient of restitution (damping constant in the DE simulations) between the base and moving particles are altered to reduce dissipation, and inject energy from the base. The energy injection does result in an increase in the strain rate by up to an order of magnitude, and the temperature by up to two orders of magnitude at the base. However, the volume fraction, strain rate and temperature profiles in the bulk (above about 15 particle diameters from the base) are altered very little by the energy injection at the base. We also examine the variation of h(stop), the minimum height at the cessation of flow, with energy injection from the base. It is found that at a fixed angle of inclination, h(stop) decreases as the energy dissipation at the base decreases.

Transition due to base roughness in a dense granular flow down an inclined plane

Particle simulations based on the discrete element method are used to examine the effect of base roughness on the granular flow down an inclined plane. The base is composed of a random configuration of fixed particles, and the base roughness is decreased by decreasing the ratio of diameters of the base and moving particles. A discontinuous transition from a disordered to an ordered flow state is observed when the ratio of diameters of base and moving particles is decreased below a critical value. The ordered flowing state consists of hexagonally close packed layers of particles sliding over each other. The ordered state is denser (higher volume fraction) and has a lower coordination number than the disordered state, and there are discontinuous changes in both the volume fraction and the coordination number at transition. The Bagnold law, which states that the stress is proportional to the square of the strain rate, is valid in both states. However, the Bagnold coefficients in the ordered flowing state are lower, by more than two orders of magnitude, in comparison to those of the disordered state. The critical ratio of base and moving particle diameters is independent of the angle of inclination, and varies very little when the height of the flowing layer is doubled from about 35 to about 70 particle diameters. While flow in the disordered state ceases when the angle of inclination decreases below 20 degrees, there is flow in the ordered state at lower angles of inclination upto 14 degrees. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4710543]

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.

Dynamic strain ageing in Ni-base superalloy 720Li

An experimental investigation into the dynamic strain ageing (DSA) of a wrought Ni-base superalloy 720Li was conducted. Characteristics of jerky, flow have been studied at intermediate temperatures of 350, 400 and 450 degrees C at strain-rates between 10(-3) and 10(-5) s(-1). Serrations of Type C are predominant within the temperature/strain-rate range explored. The major characteristics of the serrations-i.e. (a) critical plastic strain for onset of serrations, epsilon(c); (b) average stress decrement, Delta sigma(avg); and (c) strain increment between serrations. Delta epsilon(BS)-have been examined at selected temperatures and strain-rates. Negative strain-rate sensitivity was observed in the DSA regime. However. temperature did not influence tensile properties such as yield strength, ultimate strength. elongation, reduction in area, and work hardening rate or fracture features in DSA regime. Analysis of the results Suggests that locking of the mobile dislocations by substitutional alloying elements is responsible for the DSA in alloy 720Li.

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.