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.

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

Flow and heat transfer in a stagnation-point flow over a stretching sheet with a magnetic field

The flow and heat transfer problem in the boundary layer induced by a continuous moving surface is important in many manufacturing processes in industry such as the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheet, the cooling of an infinite metalic plate in a cooling bath (which may also be electrolyte). Glass blowing, continuous casting and spinning of fibres also involve the flow due to a stretching surface. Sakiadis [1] was the first to study the flow induced by a semi-infinite moving wall in an ambient fluid. On the other hand, Crane [2] first studied the flow over a linearly stretching sheet in an ambient fluid. Subsequently, Crane [3] also investigated the corresponding heat transfer problem. Since then several authors [4-8] have studied various aspects of this problem such as the effects of mass transfer, variable wall temperature, constant heat flux, magnetic field etc. Recently, Andersson [9] has obtained an exact solution of the Navier-Stokes equations for the MHD flow over a linearly stretching sheet in an ambient fluid. Also Chiam [10] has studied the heat transfer with variable thermal conductivity on a stretching sheet when the velocities of the sheet and the free stream are equal.

Time dependent rotational flow of a viscous fluid over an infinite porous disk with a magnetic field

Both the semi-similar and self-similar flows due to a viscous fluid rotating with time dependent angular velocity over a porous disk of large radius at rest with or without a magnetic field are investigated. For the self-similar case the resulting equations for the suction and no mass transfer cases are solved numerically by quasilinearization method whereas for the semi-similar case and injection in the self-similar case an implicit finite difference method with Newton's linearization is employed. For rapid deceleration of fluid and for moderate suction in the case of self-similar flow there exists a layer of fluid, close to the disk surface where the sense of rotation is opposite to that of the fluid rotating far away. The velocity profiles in the absence of magnetic field are found to be oscillatory except for suction. For the accelerating freestream, (semi-similar flow) the effect of time is to reduce the amplitude of the oscillations of the velocity components. On the other hand the effect of time for the oscillating case is just the opposite.

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.

Investigations of separated flow over backward facing steps in IISc hypersonic shock tunnel

Two backward facing step (2 mm and 3 mm step height) models are selected for surface heat transfer measurements. The platinum thin film gauges are deposited on the Macor inserts using both hand paint and vacuum sputtering technique. Using the Eckert reference temperature method the heating rates has been theoretically calculated along the flat plate portion of the model and the theoretical estimates are compared with experimentally determined surface heat transfer rate. Theoretical analysis of heat flux distribution down stream of the backward facing step model has been carried out using Gai’s non-dimensional analysis. Based on the measured surface heating rates on the backward facing step, the reattachment distance is estimated for 2 and 3 mm step height at nominal Mach number of 7.6. It has been found from the present study that for 2 and 3 mm step height, it approximately takes about 10 and 8 step heights downstream of the model respectively for the flow to re-attach.

Vesicle formation from dimeric ion-paired amphiphiles. Control over vesicular thermotropic and ion-transport properties as a function of intra-amphiphilic headgroup separation

A novel series of vesicle-forming ion-paired amphiphiles, bis(hexadecyldimethylammonium)alkane dipalmitate (1a-1h), containing four chains were synthesized with two isolated headgroups. In each of these amphiphiles, the two headgroup charges are separated by a flexible polymethylene spacer chain -[(CH2)(m)]- of varying lengths (m) such that the length and the conformation of the spacer chain determine the intra-"monomer" headgroup separation. Transmission electron microscopy indicated that each of these forms bilayer membranes upon dispersion in aqueous media. The vesicular properties of these aggregates have been examined by differential scanning calorimetry and temperature-dependent fluorescence anisotropy measurements. Interestingly, their T-m values decreased with the increase in the m value. Thus while the apparent T-m of the lipid with m = 2 (1a) is 74.1 degrees C, the corresponding value observed for the lipid with m = 12 (1h) is 38.9 degrees C. The fluorescence anisotropy values (r) for 1b-1g were quite high (r similar to 0.3) compared to that of 1h (r similar to 0.23) at 20-30 degrees C in their gel states. On the other hand, the r value for vesicular 1b beyond melting was higher (0.1) compared to any of those for 1c-1h (similar to 0.04-0.06). X-ray diffraction of the cast films was performed to understand the nature and the thickness of these membrane organizations. The membrane widths ranged from 30 to 51 A as the m values varied. The entrapment of a small water-soluble solute, riboflavin, by the individual vesicular aggregates, and their sustenance: under an imposed transmembrane pH gradient have also been examined. These results show that all lipid vesicles entrap riboflavin and that generally the resistance to OH- permeation decreases with the increase in m value. Finally,all the above observations were comparatively analyzed, and on the basis of the calculated structures of these lipids, it was possible to conclude that membrane propel-ties can be modulated by spacer chain length variation of the ion-paired amphiphiles.

Some salient features of the atmosphere observed over the north Bay of Bengal during BOBMEX

This paper describes the near surface characteristics and vertical variations based on the observations made at 17.5degreesN and 89degreesE from ORV Sagar Kanya in the north Bay of Bengal during the Bay of Bengal Monsoon Experiment (BOBMEX) carried out in July-August 1999. BOBMEX captured both the active and weak phases of convection. SST remained above the convection threshold throughout the BOBMEX. While the response of the SST to atmospheric forcing was clearly observed, the response of the atmosphere to SST changes was not clear. SST decreased during periods of large scale precipitation, and increased during a weak phase of convection. It is shown that the latent heat flux at comparable wind speeds was about 25-50% lower over the Bay during BOBMEX compared to that over the Indian Ocean during other seasons and tropical west Pacific. On the other hand, the largest variations in the surface daily net heat flux are observed over the Bay during BOBMEX. SST predicted using observed surface fluxes showed that 1-D heat balance model works sometime but not always, and horizontal advection is important. The high resolution Vaisala radiosondes launched during BOBMEX could clearly bring out the changes in the vertical structure of the atmosphere between active and weak phases of convection. Convective Available Potential Energy of the surface air decreased,by 2-3 kJ kg(-1) following convection, and recovered in a time period of one or two days. The mid tropospheric relative humidity and water vapor content, and wind direction show the major changes between the active and weak phases of convection.

Surface pressure and heat transfer measurements in the unsteady separated hypersonic flow field over double cones

Reliable bench mark experimental database in the separated hypersonic flow regime is necessary to validate high resolution CFD codes. In this paper we report the surface pressure and heat transfer measurements carried out on double cones (first cone semi-apex angle = 15, 25 deg.; second cone semi-apex angle= 35, 68 deg.) at hypersonic speeds that will be useful for CFD code validation studies. The surface pressure measurements are carried out at nominal Mach number of 8.35 in the IISc hypersonic wind tunnel. On the other hand the surface heat transfer measurements are carried out at a nominal Mach number of 5.75 in the IISc hypersonic shock tunnel. The flow separation point on the first cone, flow reattachment on the second cone and the wild fluctuation of the transmitted shock on the second cone surface (25/68 deg. double cone) in the presence of severe adverse pressure gradient are some of the flow features captured in the measurements. The results from the CFD studies indicate good agreement with experiments in the attached flow regime while considerable differences are noticeable in the separated flow regime.

Estimation of heavy rainfall during cyclonic storms from microwave observations using nonlinear approach over Indian Ocean

In this study, an effort has been made to study heavy rainfall events during cyclonic storms over Indian Ocean. This estimate is based on microwave observations from tropical rainfall measuring mission (TRMM) Microwave Imager (TMI). Regional scattering index (SI) developed for Indian region based on measurements at 19-, 21- and 85-GHz brightness temperature and polarization corrected temperature (PCT) at 85 GHz have been utilized in this study. These PCT and SI are collocated against Precipitation Radar (PR) onboard TRMM to establish a relationship between rainfall rate, PCT and SI. The retrieval technique using both linear and nonlinear regressions has been developed utilizing SI, PCT and the combination of SI and PCT. The results have been compared with the observations from PR. It was found that a nonlinear algorithm using combination of SI and PCT is more accurate than linear algorithm or nonlinear algorithm using either SI or PCT. Statistical comparison with PR exhibits the correlation coefficients (CC) of 0.68, 0.66 and 0.70, and root mean square error (RMSE) of 1.78, 1.96 and 1.68 mm/h from the observations of SI, PCT and combination of SI and PCT respectively using linear regressions. When nonlinear regression is used, the CC of 0.73, 0.71, 0.79 and RMSE of 1.64, 1.95, 1.54 mm/h are observed from the observations of SI, PCT and combination of SI and PCT, respectively. The error statistics for high rain events (above 10 mm/h) shows the CC of 0.58, 0.59, 0.60 and RMSE of 5.07, 5.47, 5.03 mm/h from the observations of SI, PCT and combination of SI and PCT, respectively, using linear regression, and on the other hand, use of nonlinear regression yields the CC of 0.66, 0.64, 0.71 and RMSE of 4.68, 5.78 and 4.02 mm/h from the observations of SI, PCT and combined SI and PCT, respectively.

DRIFTS studies on CO and NO adsorption and NO plus CO reaction over Pd2+-substituted CeO2 and Ce0.75Sn0.25O2 catalysts

Sequential adsorption of CO and NO as well as equimolar NO + CO reaction with variation of temperature over Pd2+ ion-substituted CeO2 and Ce0.75Sn0.25O2 supports has been studied by DRIFTS technique. The results are compared with 2 at.% Pd/Al2O3 containing Pd-0. Both linear and bridging Pd-0-CO bands are observed over 2 at.% Pd/Al2O3. But, band positions are shifted to higher frequencies in Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts that could be associated with Pd delta+-CO species. In contrast, a Pd2+-CO band at 2160 cm(-1) is observed upon CO adsorption over Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts pre-adsorbed with NO and a Pd+-CO band at 2120 cm(-1) is slowly developed on Ce(0.73)Srl(0.25)Pd(0.02)O(2-delta) over time. An intense linear Pd-0-NO band at 1750 cm(-1) found upon NO exposure to CO pre-adsorbed 2 at.% Pd/Al2O3 indicates molecular adsorption of NO. On the other hand, a weak Pd2+-NO band at 1850 cm(-1) is noticed after NO exposure to Ce0.98Pd0.02O2-delta catalyst pre-adsorbed with CO indicating dissociative adsorption of NO which is crucial for NO reduction. Pd-0-NO band is initially formed over CO pre-adsorbed Ce0.73Sn0.25Pd0.02O2-delta which is red-shifted over time along with formation of Pd2+-NO band. Several intense bands related to nitrates and nitrites are observed after exposure of NO to fresh as well as CO pre-adsorbed Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts. Ramping the temperature in a DRIFTS cell upon NO and CO adsorption shows the formation of N2O and NCO surface species, and N2O-formation temperature is comparable with the reaction done in a reactor.

Reply to Discussion on ``Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis'' Comput. Geotech. 70 (2015) 138-149]

A discussion has been provided for the comments raised by the discusser (Clausen, 2015)1] on the article recently published by the authors (Chakraborty and Kumar, 2015). The effect of exponent alpha for values of GSI approximately smaller than 30 becomes more critical. On the other hand, for greater values of GSI, the results obtained by the authors earlier remain primarily independent of alpha and can be easily used. (C) 2015 Elsevier Ltd. All rights reserved.

Anthropogenic aerosol fraction over the Indian region: model simulations versus multi-satellite data analysis

Anthropogenic aerosols play a crucial role in our environment, climate, and health. Assessment of spatial and temporal variation in anthropogenic aerosols is essential to determine their impact. Aerosols are of natural and anthropogenic origin and together constitute a composite aerosol system. Information about either component needs elimination of the other from the composite aerosol system. In the present work we estimated the anthropogenic aerosol fraction (AF) over the Indian region following two different approaches and inter-compared the estimates. We espouse multi-satellite data analysis and model simulations (using the CHIMERE Chemical transport model) to derive natural aerosol distribution, which was subsequently used to estimate AF over the Indian subcontinent. These two approaches are significantly different from each other. Natural aerosol satellite-derived information was extracted in terms of optical depth while model simulations yielded mass concentration. Anthropogenic aerosol fraction distribution was studied over two periods in 2008: premonsoon (March-May) and winter (November-February) in regard to the known distinct seasonality in aerosol loading and type over the Indian region. Although both techniques have derived the same property, considerable differences were noted in temporal and spatial distribution. Satellite retrieval of AF showed maximum values during the pre-monsoon and summer months while lowest values were observed in winter. On the other hand, model simulations showed the highest concentration of AF in winter and the lowest during pre-monsoon and summer months. Both techniques provided an annual average AF of comparable magnitude (similar to 0.43 +/- 0.06 from the satellite and similar to 0.48 +/- 0.19 from the model). For winter months the model-estimated AF was similar to 0.62 +/- 0.09, significantly higher than that (0.39 +/- 0.05) estimated from the satellite, while during pre-monsoon months satellite-estimated AF was similar to 0.46 +/- 0.06 and the model simulation estimation similar to 0.53 +/- 0.14. Preliminary results from this work indicate that model-simulated results are nearer to the actual variation as compared to satellite estimation in view of general seasonal variation in aerosol concentrations.

On the measurement of the surface energy budget over a land surface during the summer monsoon

The measurement of surface energy balance over a land surface in an open area in Bangalore is reported. Measurements of all variables needed to calculate the surface energy balance on time scales longer than a week are made. Components of radiative fluxes are measured while sensible and latent heat fluxes are based on the bulk method using measurements made at two levels on a micrometeorological tower of 10 m height. The bulk flux formulation is verified by comparing its fluxes with direct fluxes using sonic anemometer data sampled at 10 Hz. Soil temperature is measured at 4 depths. Data have been continuously collected for over 6 months covering pre-monsoon and monsoon periods during the year 2006. The study first addresses the issue of getting the fluxes accurately. It is shown that water vapour measurements are the most crucial. A bias of 0.25% in relative humidity, which is well above the normal accuracy assumed the manufacturers but achievable in the field using a combination of laboratory calibration and field intercomparisons, results in about 20 W m(-2) change in the latent heat flux on the seasonal time scale. When seen on the seasonal time scale, the net longwave radiation is the largest energy loss term at the experimental site. The seasonal variation in the energy sink term is small compared to that in the energy source term.

Spatial distribution of aerosol black carbon over India during pre-monsoon season

Aerosol black carbon (BC) mass concentrations ([BC]), measured continuously during a multi-platform field experiment, Integrated Campaign for Aerosols gases and Radiation Budget (ICARB, March-May 2006), from a network of eight observatories spread over geographically distinct environments of India, (which included five mainland stations, one highland station, and two island stations (one each ill Arabian Sea and Bay of Bengal)) are examined for their spatio-temporal characteristics. During the period of study, [BC] showed large variations across the country, with values ranging from 27 mu g m(3) over industrial/urban locations to as low as 0.065 mu g m(-3) over the Arabian Sea. For all mainland stations, [BC] remained high compared to highland as well as island stations. Among the island stations, Port Blair (PBR) had higher concentration of BC, compared to Minicoy (MCY), implying more absorbing nature of Bay of Bengal aerosols than Arabian Sea. The highland station Nainital (NTL), in the central Himalayas, showed low values of [BC], comparable or even lower than that of the island station PBR, indicating the prevalence of cleaner environment over there. An examination of the changes in the mean temporal features, as the season advances from winter (December-February) to pre-monsoon (March-May), revealed that: (a) Diurnal variations were pronounced over all the mainland stations, with all afternoon low and a nighttime high: (b) At the islands, the diurnal variations, though resembled those over the mainlands, were less pronounced; and (c) In contrast to this, highland station showed an opposite pattern with an afternoon high and a late night or early morning low. The diurnal variations at all stations are mainly caused by the dynamics of local Atmospheric Boundary Layer (ABL), At the entire mainland as well as island stations (except HYD and DEL), [BC] showed a decreasing trend from January to May, This is attributed to the increased convective mixing and to the resulting enhanced vertical dispersal of species in the ABL. In addition, large short-period modulations were observed at DEL and HYD, which appeared to be episodic, An examination of this in the light of the MODIS-derived fire count data over India along with the back-trajectory analysis revealed that advection of BC from extensive forest fires and biomass-burning regions upwind were largely responsible for this episodic enhancement in BC at HYD and DEL.