Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

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.

Effect of process parameters on transport phenomena during solidification in the presence of electromagnetic stirring

In the present work, a numerical study is performed to predict the effect of process parameters on transport phenomena during solidification of aluminium alloy A356 in the presence of electromagnetic stirring. A set of single-phase governing equations of mass, momentum, energy and species conservation is used to represent the solidification process and the associated fluid flow, heat and mass transfer. In the model, the electromagnetic forces are incorporated using an analytical solution of Maxwell equation in the momentum conservation equations and the slurry rheology during solidification is represented using an experimentally determined variable viscosity function. Finally, the set of governing equations is solved for various process conditions using a pressure based finite volume technique, along with an enthalpy based phase change algorithm. In present work, the effect of stirring intensity and cooling rate are considered. It is found that increasing stirring intensity results in increase of slurry velocity and corresponding increase in the fraction of solid in the slurry. In addition, the increasing stirring intensity results uniform distribution of species and fraction of solid in the slurry. It is also found from the simulation that the distribution of solid fraction and species is dependent on cooling rate conditions. At low cooling rate, the fragmentation of dendrites from the solid/liquid interface is more.

Reservoir performance under uncertainty in hydrologic impacts of climate change

Relatively few studies have addressed water management and adaptation measures in the face of changing water balances due to climate change. The current work studies climate change impact on a multipurpose reservoir performance and derives adaptive policies for possible futurescenarios. The method developed in this work is illustrated with a case study of Hirakud reservoir on the Mahanadi river in Orissa, India,which is a multipurpose reservoir serving flood control, irrigation and power generation. Climate change effects on annual hydropower generation and four performance indices (reliability with respect to three reservoir functions, viz. hydropower, irrigation and flood control, resiliency, vulnerability and deficit ratio with respect to hydropower) are studied. Outputs from three general circulation models (GCMs) for three scenarios each are downscaled to monsoon streamflow in the Mahanadi river for two future time slices, 2045-65 and 2075-95. Increased irrigation demands, rule curves dictated by increased need for flood storage and downscaled projections of streamflow from the ensemble of GCMs and scenarios are used for projecting future hydrologic scenarios. It is seen that hydropower generation and reliability with respect to hydropower and irrigation are likely to show a decrease in future in most scenarios, whereas the deficit ratio and vulnerability are likely to increase as a result of climate change if the standard operating policy (SOP) using current rule curves for flood protection is employed. An optimal monthly operating policy is then derived using stochastic dynamic programming (SDP) as an adaptive policy for mitigating impacts of climate change on reservoir operation. The objective of this policy is to maximize reliabilities with respect to multiple reservoir functions of hydropower, irrigation and flood control. In variations to this adaptive policy, increasingly more weightage is given to the purpose of maximizing reliability with respect to hydropower for two extreme scenarios. It is seen that by marginally sacrificing reliability with respect to irrigation and flood control, hydropower reliability and generation can be increased for future scenarios. This suggests that reservoir rules for flood control may have to be revised in basins where climate change projects an increasing probability of droughts. However, it is also seen that power generation is unable to be restored to current levels, due in part to the large projected increases in irrigation demand. This suggests that future water balance deficits may limit the success of adaptive policy options. (C) 2010 Elsevier Ltd. All rights reserved.

Unsteady mixed convection on the stagnation-point flow adjacent to a vertical plate with a magnetic field

An analysis is performed to study the unsteady combined forced and free convection flow (mixed convection flow) of a viscous incompressible electrically conducting fluid in the vicinity of an axisymmetric stagnation point adjacent to a heated vertical surface. The unsteadiness in the flow and temperature fields is due to the free stream velocity, which varies arbitrarily with time. Both constant wall temperature and constant heat flux conditions are considered in this analysis. By using suitable transformations, the Navier-Stokes and energy equations with four independent variables (x, y, z, t) are reduced to a system of partial differential equations with two independent variables (eta, tau). These transformations also uncouple the momentum and energy equations resulting in a primary axisymmetric flow, in an energy equation dependent on the primary flow and in a buoyancy-induced secondary flow dependent on both primary flow and energy. The resulting system of partial differential equations has been solved numerically by using both implicit finite-difference scheme and differential-difference method. An interesting result is that for a decelerating free stream velocity, flow reversal occurs in the primary flow after certain instant of time and the magnetic field delays or prevents the flow reversal. The surface heat transfer and the surface shear stress in the primary flow increase with the magnetic field, but the surface shear stress in the buoyancy-induced secondary flow decreases. Further the heat transfer increases with the Prandtl number, but the surface shear stress in the secondary flow decreases.

Semi-similar solutions of the unsteady compressible second-order boundary layer flow at the stagnation point

Semi-similar solutions of the unsteady compressible laminar boundary layer flow over two-dimensional and axisymmetric bodies at the stagnation point with mass transfer are studied for all the second-order boundary layer effects when the free stream velocity varies arbitrarily with time. The set of partial differential equations governing the unsteady compressible second-order boundary layers representing all the effects are derived for the first time. These partial differential equations are solved numerically using an implicit finite-difference scheme. The results are obtained for two particular unsteady free stream velocity distributions: (a) an accelerating stream and (b) a fluctuating stream. It is observed that the total skin friction and heat transfer are strongly affected by the surface mass transfer and wall temperature. However, their variation with time is significant only for large times. The second-order boundary layer effects are found to be more pronounced in the case of no mass transfer or injection as compared to that for suction. Résumé Des solutions semi-similaires d'écoulement variable compressible de couche limite sur des corps bi-dimensionnels thermique, sont étudiées pour tous les effets de couche limite du second ordre, lorsque la vitesse de l'écoulement libre varie arbitrairement avec le temps. Le systéme d'équations aux dérivées partielles représentant tous les effets est écrit pour la premiére fois. On le résout numériquement á l'aide d'un schéma implicite aux différences finies. Les résultats sont obtenus pour deux cas de vitesse variable d'écoulement libre: (a) un écoulement accéléré et (b) un écoulement fluctuant. On observe que le frottement pariétal total et le transfert de chaleur sont fortement affectés par le transfert de masse et la température pariétaux. Néanmoins, leur variation avec le temps est sensible seulement pour des grandes durées. Les effets sont trouvés plus prononcés dans le cas de l'absence du transfert de masse ou de l'injection par rapport au cas de l'aspiration.

Mixed convection along a vertical cone

An analysis is performed to study the flow and heat transfer characteristics for the case of laminar mixed convection along a vertical circular cone. A mixed-convection parameter is introduced in the formulation of the problem such that smooth transition from one convective limit to the other is possible. The transformed conservation equations of the nonsimilar boundary layers are solved by an efficient finite-difference method.

Greening India Mission

The Government of India has announced the Greening India Mission (GIM) under the National Climate Change Action Plan. The Mission aims to restore and afforest about 10 mha over the period 2010-2020 under different sub-missions covering moderately dense and open forests, scrub/grasslands, mangroves, wetlands, croplands and urban areas. Even though the main focus of the Mission is to address mitigation and adaptation aspects in the context of climate change, the adaptation component is inadequately addressed. There is a need for increased scientific input in the preparation of the Mission. The mitigation potential is estimated by simply multiplying global default biomass growth rate values and area. It is incomplete as it does not include all the carbon pools, phasing, differing growth rates, etc. The mitigation potential estimated using the Comprehensive Mitigation Analysis Process model for the GIM for the year 2020 has the potential to offset 6.4% of the projected national greenhouse gas emissions, compared to the GIM estimate of only 1.5%, excluding any emissions due to harvesting or disturbances. The selection of potential locations for different interventions and species choice under the GIM must be based on the use of modelling, remote sensing and field studies. The forest sector provides an opportunity to promote mitigation and adaptation synergy, which is not adequately addressed in the GIM. Since many of the interventions proposed are innovative and limited scientific knowledge exists, there is need for an unprecedented level of collaboration between the research institutions and the implementing agencies such as the Forest Departments, which is currently non-existent. The GIM could propel systematic research into forestry and climate change issues and thereby provide global leadership in this new and emerging science.