Mixed Convection Induced in a Gas Flowing Past a Hot Horizontal Flat Plate

A boundary layer analysis of mixed convective motion over a hot horizontal flat plate is performed under the conditions of steady flow and low speed. Use of the Howarth-Dorodnytsyn transformation makes it possible to dispense with the usual Boussinesq approximation, and variable gas properties are accounted for via the assumption that dynamic viscosity and thermal conductivity are proportional to the absolute temperature. The formulation presented enables the entire mixed convection regime to be described by a single set of equations. Finite difference solutions when the Prandtl number is 0.72 are obtained over the entire range of the mixed convection parameter ξ from 0 (free convection) to 1 (forced convection) and heating parameter ▵ values from 2 to 12. The effects of both ξ and ▵on the velocity profiles, the temperature profiles, and the variation of skin friction and heat transfer functions are clearly illustrated in tables and graphs.

Radiative interactions in boundary layers

A detailed description of radiative interactions in laminar compressible boundary layers for moderate Mach numbers is presented by way of asymptotic analysis and supporting solutions. The radiation field is described by the differential approximation. While the asymptotic analysis is valid for large N (the ratio of photon mean free path to molecular mean free path) and arbitrary Boltzmann number, Bo (the ratio of convective heat flux to radiation heat flux), the solutions are obtained for Bo [double less-than sign] 1, the case of strong radiative interactions. The asymptotic analysis shows the existence of an optically thin boundary layer for large N and all Bo. For Bo [double less-than sign] 1, two outer regions are observed — one optically thin (at short distances from the leading edge) and the other optically thick (at large distances from the leading edge). An interesting feature not pointed out in the previous literature is the existence of a wall layer at large distances from the leading edge where convective heat flux can be ignored to the leading order of approximation. The radiation field in all cases can be very well approximated by a one-dimensional description. The solutions have been constructed using the ideas of matched asymptotic expansions by approximate analytical procedures and numerical methods. It is shown that, to the leading order of approximation, the radiation slip method yields exactly the same result as the more complicated matching procedure. Both the cases of linear and nonlinear radiation have been considered, the former being of interest in developing approximate methods which are subsequently generalized to handle the nonlinear problem. Detailed results are presented for both cases.

Performance of a small solar pond in the tropics

The performance of the 240 m2 solar pond in Bangalore is discussed. The problems of erosion of gradient zone and formation of internal convective zones is highlighted. The technique of passive salt addition is shown to be a viable alternative for salt recycling. Different techniques of heat extraction are discussed and the use of an immersed copper heat exchanger is shown to be most convenient. A two-zone model for prediction of the seasonal structure of the solar pond performance is proposed. The model is shown to simulate the seasonal structure of the observed variation of the temperature in the storage zone.

Indian Ocean Dipole: Processes and impacts

Equatorial Indian Ocean is warmer in the east, has a deeper thermocline and mixed layer, and supports a more convective atmosphere than in the west. During certain years, the eastern Indian Ocean becomes unusually cold, anomalous winds blow from east to west along the equator and southeastward off the coast of Sumatra, thermocline and mixed layer lift up and the atmospheric convection gets suppressed. At the same time, western Indian Ocean becomes warmer and enhances atmospheric convection. This coupled ocean-atmospheric phenomenon in which convection, winds, sea surface temperature (SST) and thermocline take part actively is known as the Indian Ocean Dipole (IOD). Propagation of baroclinic Kelvin and Rossby waves excited by anomalous winds, play an important role in the development of SST anomalies associated with the IOD. Since mean thermocline in the Indian Ocean is deep compared to the Pacific, it was believed for a long time that the Indian Ocean is passive and merely responds to the atmospheric forcing. Discovery of the IOD and studies that followed demonstrate that the Indian Ocean can sustain its own intrinsic coupled ocean-atmosphere processes. About 50% percent of the IOD events in the past 100 years have co-occurred with El Nino Southern Oscillation (ENSO) and the other half independently. Coupled models have been able to reproduce IOD events and process experiments by such models – switching ENSO on and off – support the hypothesis based on observations that IOD events develop either in the presence or absence of ENSO. There is a general consensus among different coupled models as well as analysis of data that IOD events co-occurring during the ENSO are forced by a zonal shift in the descending branch of Walker cell over to the eastern Indian Ocean. Processes that initiate the IOD in the absence of ENSO are not clear, although several studies suggest that anomalies of Hadley circulation are the most probable forcing function. Impact of the IOD is felt in the vicinity of Indian Ocean as well as in remote regions. During IOD events, biological productivity of the eastern Indian Ocean increases and this in turn leads to death of corals over a large area.Moreover, the IOD affects rainfall over the maritime continent, Indian subcontinent, Australia and eastern Africa. The maritime continent and Australia suffer from deficit rainfall whereas India and east Africa receive excess. Despite the successful hindcast of the 2006 IOD by a coupled model, forecasting IOD events and their implications to rainfall variability remains a major challenge as understanding reasons behind an increase in frequency of IOD events in recent decades.

Mixed convection heat transfer from the bottom tip of a cylinder spinning about a vertical axis in a saturated porous medium

The present work is a numerical study of heat transfer characteristics from the bottom tip of a cylinder spinning about a vertical axis in an infinitely saturated porous medium. The problem is axisymmetric. The non-dimensionalized governing equations are solved using the SIMPLER algorithm on a staggered grid. The influence of rotational Reynolds numbers and Darcy numbers on the heat transfer for a Grashof number of 104 and Prandtl number of 7.0 is studied. It is found that for very high Darcy numbers, over a wide range of rotational Reynolds numbers, the heat transfer takes place mainly due to conduction. The convective heat transfer takes place for lower Darcy numbers and for higher rotational Reynolds numbers. Moreover, there is a rapid increase in the overall Nusselt number below a certain Darcy number with increase in the rotational Reynolds numbers. The effect of the Darcy number and the rotational Reynolds number on the heat transfer and fluid flow in the porous medium is depicted in the form of streamline and isotherm plots. The variation of the overall Nusselt number with respect to the Darcy number for various rotational Reynolds numbers is plotted. The variation of the local Nusselt number with respect to the radial coordinate at the heated tip of the vertical cylinder is plotted for various Darcy and rotational Reynolds numbers.

Effect of substrate on particle arrangement in arrays formed by self-assembly of polymer grafted nanoparticles

We show that the substrate affects the interparticle spacing in monolayer arrays with hexagonal order formed by self-assembly of polymer grafted nanoparticles. Remarkably, arrays with square packing were formed due to convective shearing at a liquid surface induced by miscibility of colloidal solution with the substrate.

Critical stress intensity factors for cracked hollow pipes under transient thermal loads

Finite element analyses of a long hollow cylinder having an axisymmetric circumferential internal edge crack, subjected to convective cooling on the inner surface are performed. The transient thermal stress intensity factor is estimated using a domain version of the J-integral method. The effect of the thickness of the cylinder, crack length, and heat transfer coefficient on the stress intensity factor history are studied. The variations of critical normalized stress intensity factor with crack length-to-thickness ratio for different parameters are presented. The results show that if a small inner surface crack begins to grow, its stress intensity factor will increase with increase in crack length, reach a maximum, and then begin to drop. Based on the results, a fracture-based design methodology for cracked hollow pipes under transient thermal loads is discussed.

Effect of convection on domain growth during demixing transitions in fluids

Current analytical work on the effect of convection on the late stages of spinodal decomposition in liquids is briefly described. The morphology formed during the spinodal decomposition process depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport has a significant effect on the scaling laws in the late-stage coarsening of droplets in translational or shear flows. In addition, convective transport could result in an attractive interaction between non-Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near-symmetric quench was analysed using an area distribution function, which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to time t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Rayleigh-Benard convection during solidification of an eutectic solution cooled from the top

The interaction between laminar Rayleigh-Benard convection and directional solidification is studied for the case of an eutectic solution kept in a rectangular cavity cooled from the top. Experiments and numerical simulations are carried out using an NH4Cl-H2O solution as the model fluid. The flow is visualized using a sheet of laser light scattered by neutrally buoyant, hollow-glass spheres seeded in the fluid. The numerical modeling is performed using a pressure-based finite-volume method according to the SIMPLER algorithm. The present configuration enables us to visualize flow vortices in the presence of a continuously evolving solid/liquid interface. Clear visualization of the Rayleigh-Benard convective cells and their interaction with the solidification front are obtained. It is observed that the convective cells are characterized by zones of up-flow and down-flow, resulting in the development of a nonplanar interface. Because of the continuous advancement of the solid/liquid interface, the effective liquid height of the cavity keeps decreasing. Once the height of the fluid layer falls below a critical value, the convective cells become weaker and eventually die out, leading to the growth of a planar solidification front. Results of flow visualization and temperature measurement are compared with those from the numerical simulation, and a good agreement is found.

Studies on the effects of test gas on the flow field around large angle blunt cone flying at hypersonic Mach number

The effect of the test gas on the flow field around a 120degrees apex angle blunt cone has been investigated in a shock tunnel at a nominal Mach number of 5.75. The shock standoff distance around the blunt cone was measured by an electrical discharge technique using both carbon dioxide and air as test gases. The forebody laminar convective heat transfer to the blunt cone was measured with platinum thin-film sensors in both air and carbon dioxide environments. An increase of 10 to 15% in the measured heat transfer values was observed with carbon dioxide as the test gas in comparison to air. The measured thickness of the shock layer along the stagnation streamline was 3.57 +/- 0.17 mm in air and 3.29 +/- 0.26 mm in carbon dioxide. The computed thickness of the shock layer for air and carbon dioxide were 3.98 mm and 3.02 mm, respectively. The observed increase in the measured heat transfer rates in carbon dioxide compared to air was due to the higher density ratio across the bow shock wave and the reduced shock layer thickness.

A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

The use of a fibre anemometer in turbulent flows

Quartz fibre anemometers have been used (as described in subsequent papers) to survey the velocity field of turbulent free convective air flows. This paper discusses the reasons for the choice of this instrument and provides the background information for its use in this way. Some practical points concerning fibre anemometers are mentioned. The rest of the paper is a theoretical study of the response of a fibre to a turbulent flow. An approximate representation of the force on the fibre due to the velocity field and the equation for a bending beam, representing the response to this force, form the basis of a consideration of the mean and fluctuating displacement of the fibre. Emphasis is placed on the behaviour when the spectrum of the turbulence is largely in frequencies low enough for the fibre to respond effectively instantaneously (as this corresponds to the practical situation). Incomplete correlation of the turbulence along the length of the fibre is taken into account. Brief mention is made to the theory of the higher-frequency (resonant) response in the context of an experimental check on the applicability of the low-frequency theory.

Transition to turbulence in the free convection boundary layers on an inclined heated plate

An experimental study has been made of transition to turbulence in the free convective flows on a heated plate. Observations have been made with the plate vertical and inclined at angles up to about 50° to the vertical, both above and below the plate. A fibre anemometer was used to survey the region of intermittent turbulence. Information has thus been obtained about the range of Grashof numbers over which transition takes place. Even when the plate is vertical the region of intermittent turbulence is long. When it is inclined, this region becomes still longer in the flow below the plate as a result of the stabilizing stratification, a Richardson number effect. It is possible to have a whole flow such that it should be described as transitional, not laminar or turbulent. It was noticed that in this flow and the vertical plate one, the velocity during the laminar periods could be either of two characteristic values, one of them close to zero. The behaviour above an inclined plate could be interpreted largely as a trend towards the behaviour described in a preceding paper.

Measurement of Heat Transfer Rate on Backward-Facing Steps at Hypersonic Mach Number

Two backward-facing models with step heights of 2 and 3 mm are used to measure the convective surface heat transfer rates by using platinum thin-film gauges, deposited on Macor inserts. Heat transfer rates have been theoretically calculated along the flat plate portion of a model using the Eckert reference temperature method. The experimentally determined surface heat transfer rate distributions are compared with theoretical and numerical estimations. Experimental heat flux distribution over a flat plate model showed good agreement with the reference temperature method at stagnation enthalpy range of 0.8-2 MJ/kg. Theoretical analysis has been used for downstream of a backward-facing step using Gai's nondimensional analysis. It has been found from the present study that approximately 10 and 8 step heights are required for the flow to reattach for 2 and 3 mm step height backward-facing step models, respectively, at a nominal Mach number of 7.6.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.