Mathematical Modelling of the Dynamics of Granular Materials in a Rotating Cylinder

The current study is aimed at the development of a theoretical simulation tool based on Discrete Element Method (DEM) to 'interpret granular dynamics of solid bed in the cross section of the horizontal rotating cylinder at the microscopic level and subsequently apply this model to establish the transition behaviour, mixing and segregation.The simulation of the granular motion developed in this work is based on solving Newton's equation of motion for each particle in the granular bed subjected to the collisional forces, external forces and boundary forces. At every instant of time, the forces are tracked and the positions velocities and accelarations of each partcle is The software code for this simulation is written in VISUAL FORTRAN 90 After checking the validity of the code with special tests, it is used to investigate the transition behaviour of granular solids motion in the cross section of a rotating cylinder for various rotational speeds and fill fraction.This work is hence directed towards a theoretical investigation based on Discrete Element Method (DEM) of the motion of granular solids in the radial direction of the horizontal cylinder to elucidate the relationship between the operating parameters of the rotating cylinder geometry and physical properties ofthe granular solid.The operating parameters of the rotating cylinder include the various rotational velocities of the cylinder and volumetric fill. The physical properties of the granular solids include particle sizes, densities, stiffness coefficients, and coefficient of friction Further the work highlights the fundamental basis for the important phenomena of the system namely; (i) the different modes of solids motion observed in a transverse crosssection of the rotating cylinder for various rotational speeds, (ii) the radial mixing of the granular solid in terms of active layer depth (iii) rate coefficient of mixing as well as the transition behaviour in terms of the bed turnover time and rotational speed and (iv) the segregation mechanisms resulting from differences in the size and density of particles.The transition behaviour involving its six different modes of motion of the granular solid bed is quantified in terms of Froude number and the results obtained are validated with experimental and theoretical results reported in the literature The transition from slumping to rolling mode is quantified using the bed turnover time and a linear relationship is established between the bed turn over time and the inverse of the rotational speed of the cylinder as predicted by Davidson et al. [2000]. The effect of the rotational speed, fill fraction and coefficient of friction on the dynamic angle of repose are presented and discussed. The variation of active layer depth with respect to fill fraction and rotational speed have been investigated. The results obtained through simulation are compared with the experimental results reported by Van Puyvelde et. at. [2000] and Ding et at. [2002].The theoretical model has been further extended, to study the rmxmg and segregation in the transverse direction for different particle sizes and their size ratios. The effect of fill fraction and rotational speed on the transverse mixing behaviour is presented in the form of a mixing index and mixing kinetics curve. The segregation pattern obtained by the simulation of the granular solid bed with respect to the rotational speed of the cylinder is presented both in graphical and numerical forms. The segregation behaviour of the granular solid bed with respect to particle size, density and volume fraction of particle size has been investigated. Several important macro parameters characterising segregation such as mixing index, percolation index and segregation index have been derived from the simulation tool based on first principles developed in this work.

Analytical Studies on Diffusion and lntermittency in Chaotic Maps

The study of simple chaotic maps for non-equilibrium processes in statistical physics has been one of the central themes in the theory of chaotic dynamical systems. Recently, many works have been carried out on deterministic diffusion in spatially extended one-dimensional maps This can be related to real physical systems such as Josephson junctions in the presence of microwave radiation and parametrically driven oscillators. Transport due to chaos is an important problem in Hamiltonian dynamics also. A recent approach is to evaluate the exact diffusion coefficient in terms of the periodic orbits of the system in the form of cycle expansions. But the fact is that the chaotic motion in such spatially extended maps has two complementary aspects- - diffusion and interrnittency. These are related to the time evolution of the probability density function which is approximately Gaussian by central limit theorem. It is noticed that the characteristic function method introduced by Fujisaka and his co-workers is a very powerful tool for analysing both these aspects of chaotic motion. The theory based on characteristic function actually provides a thermodynamic formalism for chaotic systems It can be applied to other types of chaos-induced diffusion also, such as the one arising in statistics of trajectory separation. It was noted that there is a close connection between cycle expansion technique and characteristic function method. It was found that this connection can be exploited to enhance the applicability of the cycle expansion technique. In this way, we found that cycle expansion can be used to analyse the probability density function in chaotic maps. In our research studies we have successfully applied the characteristic function method and cycle expansion technique for analysing some chaotic maps. We introduced in this connection, two classes of chaotic maps with variable shape by generalizing two types of maps well known in literature.

Studies on the Characteristics of Vertical Circulation and Equatorial Waves using Indian MST Radar

The aim of the present study is to understand the characteristics and properties of different wave modes and the vertical circulation pattern in the troposphere and lower stratosphere over Indian region using data obtained from the Indian Mesosphere-Stratosphere Troposphere (MST) radar, National Center for Environmental Prediction/National Centres of Atmospheric Research (NCEP/NCAR) reanalysed data and radiosonde observations.Studies on the vertical motion in monsoon Hadley circulation are carried out and the results are discussed . From the analysis of MST radar data, an overall picture of vertical motion of air over Indian region is explained and noted that there exists sinking motion both during winter and summer. Besides, the study shows that there is an anomalous northerly wind in the troposphere over the southern peninsular region during southwest monsoon season.The outcome of the study on intrusion of mid-latitude upper tropospheric trough and associated synoptic-scale vertical velocity over the tropical Indian latitudes are reported and discussed . It shows that there is interaction between north Indian latitudes and tropical easterly region, when there is an eastward movement of Western Disturbance across the country. It explains the strengthening of westerlies and a change of winter westerlies into easterlies in the tropical troposphere and lower stratosphere. The divergence field computed over the MST radar station shows intensification in the downward motion in association with the synoptic systems of the northwest Indian region.

SST–convection relation over tropical oceans

According to current knowledge, convection over the tropical oceans increases with sea surface temperature (SST) from 26 to 29 °C, and at SSTs above 29 °C, it sharply decreases. Our research shows that it is only over the summer warm pool areas of Indian and west Pacific Oceans (monsoon areas) where the zone of maximum SST is away from the equator that this kind of SST-convection relationship exists. In these areas (1) convection is related to the SST gradient that generates low-level moisture convergence and upward vertical motion in the atmosphere. This has modelling support. Regions of SST maxima have low SST gradients and therefore feeble convection. (2) Convection initiated by SST gradient produces strong wind fields particularly cross-equatorial low-level jetstreams (LLJs) on the equator-ward side of the warm pool and both the convection and LLJ grow through a positive feedback process. Thus, large values of convection are associated with the cyclonic vorticity of the LLJ in the atmospheric boundary layer. In the inter-tropical convergence zone (ITCZ) over the east Pacific Ocean and the south Pacific convergence zone (SPCZ) over the west Pacific Ocean, low-level winds from north and south hemisphere converge in the zone of maximum SST, which lies close to the equator producing there elongated bands of deep convection, where we find that convection increases with SST for the full range of SSTs unlike in the warm pool regions. The low-level wind divergence computed using QuikSCAT winds has large and significant linear correlation with convection in both the warm pool and ITCZ/SPCZ areas. But the linear correlation between SST and convection is large only for the ITCZ/SPCZ. These findings have important implications for the modelling of largescale atmospheric circulations and the associated convective rainfall over the tropical oceans