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.

Steady mixed convection stagnation-point flow of upper convected Maxwell fluids with magnetic field

The steady MHD mixed convection flow of a viscoelastic fluid in the vicinity of two-dimensional stagnation point with magnetic field has been investigated under the assumption that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary layer theory is used to simplify the equations of motion. induced magnetic field and energy which results in three coupled non-linear ordinary differential equations which are well-posed. These equations have been solved by using finite difference method. The results indicate the reduction in the surface velocity gradient, surface heat transfer and displacement thickness with the increase in the elasticity number. These trends are opposite to those reported in the literature for a second-grade fluid. The surface velocity gradient and heat transfer are enhanced by the magnetic and buoyancy parameters. The surface heat transfer increases with the Prandtl number, but the surface velocity gradient decreases.

Vertical distribution of aerosols over the east coast of India inferred from airborne LIDAR measurements

The information on altitude distribution of aerosols in the atmosphere is essential in assessing the impact of aerosol warming on thermal structure and stability of the atmosphere.In addition, aerosol altitude distribution is needed to address complex problems such as the radiative interaction of aerosols in the presence of clouds. With this objective,an extensive, multi-institutional and multi-platform field experiment (ICARB-Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out under the Geosphere Biosphere Programme of the Indian Space Research Organization (ISRO-GBP) over continental India and adjoining oceans during March to May 2006. Here, we present airborne LIDAR measurements carried out over the east Coast of the India during the ICARB field campaign. An increase in aerosol extinction (scattering + absorption) was observed from the surface upwards with a maximum around 2 to 4 km. Aerosol extinction at higher atmospheric layers (>2 km) was two to three times larger compared to that of the surface. A large fraction (75-85%) of aerosol column optical depth was contributed by aerosols located above 1 km. The aerosol layer heights (defined in this paper as the height at which the gradient in extinction coefficient changes sign) showed a gradual decrease with an increase in the offshore distance. A large fraction (60-75%) of aerosol was found located above clouds indicating enhanced aerosol absorption above clouds. Our study implies that a detailed statistical evaluation of the temporal frequency and spatial extent of elevated aerosol layers is necessary to assess their significance to the climate. This is feasible using data from space-borne lidars such as CALIPSO,which fly in formation with other satellites like MODIS AQUA and MISR, as part of the A-Train constellation.

Transition of turbulent pipe flow

The velocity profile in turbulent pipe flow is usually divided into two regions, a wall or inner region and a core or outer region. For the inner region, the viscosity and wall shear stress are the important parameters governing the velocity distribution whereas for the outer region, the wall reduces the velocity below the maximum velocity independent of viscosity. In the present work, a velocity model is proposed for turbulent flow in the wall region of a pipe covering the entire transition from smooth to rough flows. Coupling this model for the wall region with the power law velocity model for the core region, an equation for the friction factor is obtained. The model constants are evaluated by using Nikuradse's experiments in the fully smooth and rough turbulent flows. The model shows good agreement with the friction factor and the velocity profiles obtained by Nikuradse for the transition region of turbulent flow.

Unsteady Natural Convection Flow over a Heated Cylinder Buried in a Fluid Saturated Porous Medium

Transient natural convection flow on a heated cylinder buried in a semi-infinite liquid-saturated porous medium has been studied. The unsteadiness in the problem arises due to the cylinder which is heated (cooled) suddenly and then maintained at that temperature. The coupled partial differential equations governing the flow and heat transfer are cast into stream function-temperature formulation, and the solutions are obtained from the initial time to the time when steady state is reached. The heat transfer is found to change significantly with increasing time in a small time interval immediately after the start of the impulsive change, and steady state is reached after some time. The average Nusselt number is found to increase with Rayleigh number When the surface of the cylinder is suddenly cooled, there is a change in the direction of the heat transfer in a small time interval immediately after the start of the impulsive change in the surface temperature;however when the surface temperature is suddenly increased, no such phenomenon is observed.

Spatial and vertical heterogeneities in aerosol properties over oceanic regions around India: Implications for radiative forcing

The influence of atmospheric aerosols on Earth's radiation budget and hence climate, though well recognized and extensively investigated in recent years, remains largely uncertain mainly because of the large spatio-temporal heterogeneity and the lack of data with adequate resolution. To characterize this diversity, a major multi-platform field campaign ICARB (Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out during the pre-monsoon period of 2006 over the Indian landmass and surrounding oceans, which was the biggest such campaign ever conducted over this region. Based on the extensive and concurrent measurements of the optical and physical properties of atmospheric aerosols during ICARB, the spatial distribution of aerosol radiative forcing was estimated over the entire Bay of Bengal (BoB), northern Indian Ocean and Arabian Sea (AS) as well as large spatial variations within these regions. Besides being considerably lower than the mean values reported earlier for this region, our studies have revealed large differences in the forcing components between the BoB and the AS. While the regionally averaged aerosol-induced atmospheric forcing efficiency was 31 +/- 6 W m(-2) tau(-1) for the BoB, it was only similar to 18 +/- 7 W m(-2) tau(-1) for the AS. Airborne measurements revealed the presence of strong, elevated aerosol layers even over the oceans, leading to vertical structures in the atmospheric forcing, resulting in significant warming in the lower troposphere. These observations suggest serious climate implications and raise issues ranging from the impact of aerosols on vertical thermal structure of the atmospheric and hence cloud formation processes to monsoon circulation.

The dynamics of magnetically trapped fluids. I - Implications for umbral dots and penumbral grains

A study of the magnetohydrodynamic system in which a nonmagnetized fluid in a gravitational field is surrounded by a fluid carrying a vertical magnetic field is presented. It is pointed out that this study can throw some light on the fine-structural features of a sunspot. The equilibrium configuration of the field-free fluid is a tapering column ending at an apex. The regions away form the apex can be studied by the slender flux tube approximation. A scheme developed to treat the apex indicates that, just below the apex, the radius of the tapering column opens up with a 3/2 power dependence on the depth below the apex. If the internal pressure of the field-free fluid is increased, the apex rises, and a static equilibrium may not be possible beyond a limit if the magnetic pressure drops quickly above a certain height. The nature of steady-flow solutions beyond this limit is investigated. Under conditions inside a sunspot, a column of field-free gas is found to rise with a velocity of about 100 km/hr. If umbral dots and penumbral grains are interpreted as regions where the field-free gas ultimately emerges, a very natural explanation of most of their observed properties is obtained.

Two-temperature accretion around rotating black holes: a description of the general advective flow paradigm in the presence of various cooling processes to explain low to high luminous sources

We investigate viscous two-temperature accretion disc flows around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter of 0.998. It is found that the flow, during its infall from the Keplerian to sub-Kepleria transition region to the black hole event horizon, passes through various phases of advection: the general advective paradigm to the radiatively inefficient phase, and vice versa. Hence, the flow governs a much lower electron temperature similar to 10(8)-10(9.5) K, in the range of accretion rate in Eddington units 0.01 less than or similar to (M) over dot less than or similar to 100, compared to the hot protons of temperature similar to 10(10.2)-10(11.8) K. Therefore, the solution may potentially explain the hard X-rays and gamma-rays emitted from active galactic nuclei (AGNs) and X-ray binaries. We then compare the solutions for two different regimes of viscosity. We conclude that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart, which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the underfed AGNs and quasars (e. g. Sgr A*) to ultraluminous X-ray sources (e. g. SS433), at different combinations of input parameters, such as the mass accretion rate and the ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in highly luminous AGNs and ultraluminous quasars (e. g. PKS 0743-67).

Analyzing information flow in brain networks with nonparametric Granger causality

Multielectrode neurophysiological recording and high-resolution neuroimaging generate multivariate data that are the basis for understanding the patterns of neural interactions. How to extract directions of information flow in brain networks from these data remains a key challenge. Research over the last few years has identified Granger causality as a statistically principled technique to furnish this capability. The estimation of Granger causality currently requires autoregressive modeling of neural data. Here, we propose a nonparametric approach based on widely used Fourier and wavelet transforms to estimate both pairwise and conditional measures of Granger causality, eliminating the need of explicit autoregressive data modeling. We demonstrate the effectiveness of this approach by applying it to synthetic data generated by network models with known connectivity and to local field potentials recorded from monkeys performing a sensorimotor task.

Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

Influence of spatial variability of permeability property on steady state seepage flow and slope stability analysis

In recent years, spatial variability modeling of soil parameters using random field theory has gained distinct importance in geotechnical analysis. In the present Study, commercially available finite difference numerical code FLAC 5.0 is used for modeling the permeability parameter as spatially correlated log-normally distributed random variable and its influence on the steady state seepage flow and on the slope stability analysis are studied. Considering the case of a 5.0 m high cohesive-frictional soil slope of 30 degrees, a range of coefficients of variation (CoV%) from 60 to 90% in the permeability Values, and taking different values of correlation distance in the range of 0.5-15 m, parametric studies, using Monte Carlo simulations, are performed to study the following three aspects, i.e., (i) effect ostochastic soil permeability on the statistics of seepage flow in comparison to the analytic (Dupuit's) solution available for the uniformly constant permeability property; (ii) strain and deformation pattern, and (iii) stability of the given slope assessed in terms of factor of safety (FS). The results obtained in this study are useful to understand the role of permeability variations in slope stability analysis under different slope conditions and material properties. (C) 2009 Elsevier B.V. All rights reserved.

Boundary-layer effects on internal flow choking in dual-thrust solid rocket motors

Theoretical studies have been carried out to examine internal flow choking in the inert simulators of a dual-thrust motor. Using a two-dimensional k-omega turbulence model, detailed parametric studies have been carried out to examine aerodynamic choking and the existence of a fluid throat at the transition region during the startup transient of dual-thrust motors. This code solves standard k-omega turbulence equations with shear flow corrections using a coupled second-order-implicit unsteady formulation. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-averaged, Navier-Stokes equations is employed. It was observed that, at the subsonic inflow conditions, there is a possibility of the occurrence of internal flow choking in dual-thrust motors due to the formation of a fluid throat at the beginning of the transition region induced by area blockage caused by boundary-layer-displacement thickness. It has been observed that a 55% increase in the upstream port area of the dual-thrust motor contributes to a 25% reduction in blockage factor at the transition region, which could negate the internal How choking and supplement with an early choking of the dual-thrust motor nozzle. If the height of the upstream port relative to the motor length is too small, the developing boundary layers from either side of the port can interact, leading to a choked,flow. On the other hand, if the developing boundary layers are far enough apart, then choking does not occur. The blockage factor is greater in magnitude for the choked case than for the unchoked case. More tangible explanations are presented in this paper for the boundary-layer blockage and the internal flow choking in dual-thrust motors, which hitherto has been unexplored.

Constitutive equation for elevated temperature flow behavior of plain carbon steels using dimensional analysis

Dimensional analysis using π-theorem is applied to the variables associated with plastic deformation. The dimensionless groups thus obtained are then related and rewritten to obtain the constitutive equation. The constants in the constitutive equation are obtained using published flow stress data for carbon steels. The validity of the constitutive equation is tested for steels with up to 1.54 wt%C at temperatures: 850–1200 °C and strain rates: 6 × 10−6–2 × 10−2 s−1. The calculated flow stress agrees favorably with experimental data.

Boson-fermion duality in SU(m vertical bar n) supersymmetric Haldane-Shastry spin chain

By using the Y(gl(m|n)) super Yangian symmetry of the SU(m|n) supersymmetric Haldane-Shastry spin chain, we show that the partition function of this model satisfies a duality relation under the exchange of bosonic and fermionic spin degrees of freedom. As a byproduct of this study of the duality relation, we find a novel combinatorial formula for the super Schur polynomials associated with some irreducible representations of the Y(gl(m|n)) Yangian algebra. Finally, we reveal an intimate connection between the global SU(m|n) symmetry of a spin chain and the boson-fermion duality relation. (C) 2007 Elsevier B.V. All rights reserved.

Continuous-Flow Surface Aeration Systems

An aeration process in ail activated sludge plant is a continuous-flow system. In this system, there is a steady input flow (flow from the primary clarifier or settling tank with some part from the secondary clarifier or secondary settling tank) and output flow connection to the secondary clarifier or settling tank. The experimental and numerical results obtained through batch systems can not be relied on and applied for the designing of a continuous aeration tank. In order to scale up laboratory results for field application, it is imperative to know the geometric parameters of a continuous system. Geometric parameters have a greater influence on the mass transfer process of surface aeration systems. The present work establishes the optimal geometric configuration of a continuous-flow surface aeration system. It is found that the maintenance of these optimal geometric parameters systems result in maximum aeration efficiency. By maintaining the obtained optimal geometric parameters, further experiments are conducted in continuous-flow surface aerators with three different sizes in order to develop design curves correlating the oxygen transfer coefficient and power number with the rotor speed. The design methodology to implement the presently developed optimal geometric parameters and correlation equations for field application is discussed.