Modeling re-ignition and chuffing in solid rocket motors

A nonlinear model is developed to numerically simulate dynamic combustion inside a solid rocket motor chamber. Using this model, the phenomena of re-ignition and chuffing are investigated under low-L* conditions. The model consists of two separate submodels (coupled to each other), one for unsteady burning of propellant and the other for unsteady conservation of mass and energy within the chamber. The latter yields instantaneous pressure and temperature within the chamber. The instantaneous burning rate is calculated using a one-dimensional, nonlinear, transient gas-phase model previously developed by the authors. The results presented in this paper show that the model predicts not only the critical L*, but also the various regimes of L*-instabihty. Specifically, the results exhibit (1) amplifying pressure oscillations leading to extinction, and (2) re-ignition after a dormant period following extinction. The re-ignition could be observed only when a radiation heat flux (from the combustion chamber to the propellant surface) was included. Certain high-frequency oscillations, possibly due to intrinsic instability, are observed when the pressure overshoots during re-ignition. At very low values of initial L*, successive cycles of extinction/reignition displaying typical characteristics of chuffing are predicted. Variations of the chuffing frequency and the thickness of propellant burned off during a chuff with L* are found to be qualitatively the same as that reported from experimental observations.

Spring Warming of the Eastern Arabian Sea and Bay of Bengal from Buoy Data

Observations from moored buoys during spring of 1998-2000 suggest that the warming of the mixed layer (similar to20 m deep) of the north Indian Ocean warm pool is a response to net surface heat flux Q(net) (similar to100 W m(-2)) minus penetrative solar radiation Q(pen) (similar to45 W m(-2)). A residual cooling due to vertical mixing and advection is indirectly estimated to be about 25 W m(-2). The rate of warming due to typical values of Q(net) minus Q(pen) is not very sensitive to the depth of the mixed layer if it lies between 10 m and 30 m.

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.

Land Surface Temperature with Land Cover Dynamics:Multi-Resolution, Spatio-Temporal Data Analysis of Greater Bangalore

Bangalore is experiencing unprecedented urbanisation in recent times due to concentrated developmental activities with impetus on IT (Information Technology) and BT (Biotechnology) sectors. The concentrated developmental activities has resulted in the increase in population and consequent pressure on infrastructure, natural resources, ultimately giving rise to a plethora of serious challenges such as urban flooding, climate change, etc. One of the perceived impact at local levels is the increase in sensible heat flux from the land surface to the atmosphere, which is also referred as heat island effect. In this communication, we report the changes in land surface temperature (LST) with respect to land cover changes during 1973 to 2007. A novel technique combining the information from sub-pixel class proportions with information from classified image (using signatures of the respective classes collected from the ground) has been used to achieve more reliable classification. The analysis showed positive correlation with the increase in paved surfaces and LST. 466% increase in paved surfaces (buildings, roads, etc.) has lead to the increase in LST by about 2 ºC during the last 2 decades, confirming urban heat island phenomenon. LSTs’ were relatively lower (~ 4 to 7 ºC) at land uses such as vegetation (parks/forests) and water bodies which act as heat sinks.

Sensitivity of terrestrial water and energy budgets to CO(2)-physiological forcing: an investigation using an offline land model

Increasing concentrations of atmospheric carbon dioxide (CO(2)) influence climate by suppressing canopy transpiration in addition to its well- known greenhouse gas effect. The decrease in plant transpiration is due to changes in plant physiology (reduced opening of plant stomata). Here, we quantify such changes in water flux for various levels of CO(2) concentrations using the National Center for Atmospheric Research's (NCAR) Community Land Model. We find that photosynthesis saturates after 800 ppmv (parts per million, by volume) in this model. However, unlike photosynthesis, canopy transpiration continues to decline at about 5.1% per 100 ppmv increase in CO(2) levels. We also find that the associated reduction in latent heat flux is primarily compensated by increased sensible heat flux. The continued decline in canopy transpiration and subsequent increase in sensible heat flux at elevated CO(2) levels implies that incremental warming associated with the physiological effect of CO(2) will not abate at higher CO(2) concentrations, indicating important consequences for the global water and carbon cycles from anthropogenic CO(2) emissions.

Flow visualization in turbulent free convection over horizontal smooth and grooved surfaces

The present paper discusses the flow visualization for turbulent free convection in a tank of water with the bottom surface being a smooth or a grooved surface and the top of the water surface exposed to ambient. The grooved surface is of parallel 90 degrees V-grooves with groove height of 10 mm and groove width of 20 mm. The experiment is carried out with aspect ratio (AR) of 2.9 and Rayleigh number (Ra) in the range, 1.3 x 10(7) - 4 x 10(7). Here AR is the aspect ratio (= width of fluid layer/height of fluid layer). Heat flux at the bottom surface is from electrical heating. From the pH-dye visualization, interesting flow structures are observed and these structures are analyzed with the help of plumes dynamics and temperature variations with time. (C) 2011 Elsevier Ltd. All rights reserved.

Processes of 30-90 days sea surface temperature variability in the northern Indian Ocean during boreal summer

During summer, the northern Indian Ocean exhibits significant atmospheric intraseasonal variability associated with active and break phases of the monsoon in the 30-90 days band. In this paper, we investigate mechanisms of the Sea Surface Temperature (SST) signature of this atmospheric variability, using a combination of observational datasets and Ocean General Circulation Model sensitivity experiments. In addition to the previously-reported intraseasonal SST signature in the Bay of Bengal, observations show clear SST signals in the Arabian Sea related to the active/break cycle of the monsoon. As the atmospheric intraseasonal oscillation moves northward, SST variations appear first at the southern tip of India (day 0), then in the Somali upwelling region (day 10), northern Bay of Bengal (day 19) and finally in the Oman upwelling region (day 23). The Bay of Bengal and Oman signals are most clearly associated with the monsoon active/break index, whereas the relationship with signals near Somali upwelling and the southern tip of India is weaker. In agreement with previous studies, we find that heat flux variations drive most of the intraseasonal SST variability in the Bay of Bengal, both in our model (regression coefficient, 0.9, against similar to 0.25 for wind stress) and in observations (0.8 regression coefficient); similar to 60% of the heat flux variation is due do shortwave radiation and similar to 40% due to latent heat flux. On the other hand, both observations and model results indicate a prominent role of dynamical oceanic processes in the Arabian Sea. Wind-stress variations force about 70-100% of SST intraseasonal variations in the Arabian Sea, through modulation of oceanic processes (entrainment, mixing, Ekman pumping, lateral advection). Our similar to 100 km resolution model suggests that internal oceanic variability (i.e. eddies) contributes substantially to intraseasonal variability at small-scale in the Somali upwelling region, but does not contribute to large-scale intraseasonal SST variability due to its small spatial scale and random phase relation to the active-break monsoon cycle. The effect of oceanic eddies; however, remains to be explored at a higher spatial resolution.

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

Nucleation dynamics and pool boiling characteristics of high pressure refrigerant using thermochromic liquid crystals

This paper reports an experimental investigation of oscillating temperature field beneath a single isolated nucleation site using a non-invasive TLC (thermochromic liquid crystal) based thermography technique. Empirical correlations are presented to demonstrate the influence of system pressure and wall heat flux on different ebullition characteristics in the nucleate pool boiling regime of refrigerant R-134a. TLC transient response and two-phase flow structure are captured using synchronized, high resolution imaging. It is observed that the area of influence of nucleation site exhibits a two-part distinct transient behavior during the bubble growth period and broadens to a maximum of 1.57 times the bubble diameter at the instant of bubble departure. This is accompanied by a sharp fall of 2.5 degrees C in the local excess temperature at the nucleation site, which results in momentary augmentation (similar to 40%) in the local heat transfer coefficient at the nucleation origin. The enhanced heat transfer rate observed during the bubble peel-off event is primarily due to transient micro-convection in the wake of the retreating bubble. Further, the results indicate that a slight increase in system pressure from 813.6 to 882.5 kPa has no considerable effect on either the wall superheat or the overall heat transfer coefficient and ebullition frequency. In addition, correlations have been obtained for bubble Reynolds number, Jackob number and the dimensionless bubble generation frequency in terms of modified boiling number.

Study of buoyant jets in natural ventilation of a model room

Buoyant jets in natural ventilation of a model room with water as the fluid medium have been studied. A constant heat flux has been maintained on the bottom surface of the room. The buoyancy causes flow to enter through the bottom opening and leave through the top opening. The shadowgraph technique is used for visualization. At the inlet, a negatively buoyant jet is observed, whereas a positively buoyant jet is observed at the outlet. The theoretical results for the centerline trajectories of these buoyant jets using both Gaussian and top-hat profiles are discussed considering the variation of the entrainment coefficient with the local Froude number and the variation of the spreading ratio of buoyancy to velocity profile with the distance from the source. The shape of the profiles is found to evolve from top-hat to Gaussian geometry.

Hypersonic interference heating on flat plate with short three-dimensional protuberances

Experiments were conducted to measure the heat flux in the vicinity of a three-dimensional protuberance placed on a flat plate facing a hypersonic flow at zero angle of attack. The effects of flow enthalpy and height of the protuberance on the interference heating in its vicinity were studied. Evidence of disturbed flow with highly three-dimensional characteristics and heightened vorticity was observed near the protrusion. A parametric study by changing the deflection angle of the protuberance was also made. Correlations exist in the open literature for enthalpy values lower than 2  MJ/kg. This effort has yielded a new correlation that is valid for enthalpies up to 6  MJ/kg. The Z-type schlieren technique was used to visualize the flow features qualitatively for one of the flow conditions studied.

A simple model for spatial disaggregation of evaporative fraction: Comparative study with thermal sharpened land surface temperature data over India

[1] Evaporative fraction (EF) is a measure of the amount of available energy at the earth surface that is partitioned into latent heat flux. The currently operational thermal sensors like the Moderate Resolution Imaging Spectroradiometer (MODIS) on satellite platforms provide data only at 1000 m, which constraints the spatial resolution of EF estimates. A simple model (disaggregation of evaporative fraction (DEFrac)) based on the observed relationship between EF and the normalized difference vegetation index is proposed to spatially disaggregate EF. The DEFrac model was tested with EF estimated from the triangle method using 113 clear sky data sets from the MODIS sensor aboard Terra and Aqua satellites. Validation was done using the data at four micrometeorological tower sites across varied agro-climatic zones possessing different land cover conditions in India using Bowen ratio energy balance method. The root-mean-square error (RMSE) of EF estimated at 1000 m resolution using the triangle method was 0.09 for all the four sites put together. The RMSE of DEFrac disaggregated EF was 0.09 for 250 m resolution. Two models of input disaggregation were also tried with thermal data sharpened using two thermal sharpening models DisTrad and TsHARP. The RMSE of disaggregated EF was 0.14 for both the input disaggregation models for 250 m resolution. Moreover, spatial analysis of disaggregation was performed using Landsat-7 (Enhanced Thematic Mapper) ETM+ data over four grids in India for contrasted seasons. It was observed that the DEFrac model performed better than the input disaggregation models under cropped conditions while they were marginally similar under non-cropped conditions.

Hypersonic interference heating on cones with short three-dimensional protuberances

Heat fluxes around short, three-dimensional protuberances on sharp and blunt cones in hypersonic flow were experimentally measured using platinum thin-film sensors deposited on macor inserts. A parametric study of different protrusion geometries and flow conditions were conducted. Excessive heating was observed at locations near the protrusion where increased vorticity is expected, with the hottest spot being presented at the foot of the protuberance immediately upstream of it. If left unchecked, these hot spots could prove detrimental to hypersonic flight vehicles. Z-type schlieren technique was used to visualize the flow features qualitatively. New correlations to predict the heat flux at the hot spot have been proposed. (C) 2014 Elsevier Inc. All rights reserved.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Investigation of the separated region ahead of three-dimensional protuberances on plates and cones in hypersonic flows with laminar boundary layers

Heat transfer rate and pressure measurements were made upstream of surface pro-tuberances on a flat plate and a sharp cone subjected to hypersonic flow in a conventional shock tunnel. Heat flux was measured using platinum thin-film sensors deposited on macor substrate and the pressure measurements were made using fast acting piezoelectric sensors. A distinctive hot spot with highest heat flux was obtained near the foot of the protuberance due to heavy vortex activity in the recirculating region. Schlieren flow visualization was used to capture the shock structures and the separation distance ahead of the protrusions was quantitatively measured for varying protuberance heights. A computational analysis was conducted on the flat plate model using commercial computational fluid dynamics software and the obtained trends of heat flux and pressure were compared with the experimental observation. Experiments were also conducted by physically disturbing the laminar boundary layer to check its effect on the magnitude of the hot spot heat flux. In addition to air, argon was also used as test gas so that the Reynolds number can be varied. (C) 2014 AIP Publishing LLC.