Performance evaluation of chemistry transport models over India

Using continuous and near-real time measurements of the mass concentrations of black carbon (BC) aerosols near the surface, for a period of 1 year (from January to December 2006) from a network of eight observatories spread over different environments of India, a space-time synthesis is generated. The strong seasonal variations observed, with a winter high and summer low, are attributed to the combined effects of changes in synoptic air mass types, modulated strongly by the atmospheric boundary layer dynamics. Spatial distribution shows much higher BC concentration over the Indo-Gangetic Plain (IGP) than the peninsular Indian stations. These were examined against the simulations using two chemical transport models, GOCART (Goddard Global Ozone Chemistry Aerosol Radiation and Transport) and CHIMERE for the first time over Indian region. Both the model simulations significantly deviated from the measurements at all the stations; more so during the winter and pre-monsoon seasons and over mega cities. However, the CHIMERE model simulations show better agreement compared with the measurements. Notwithstanding this, both the models captured the temporal variations; at seasonal and subseasonal timescales and the natural variabilities (intra-seasonal oscillations) fairly well, especially at the off-equatorial stations. It is hypothesized that an improvement in the atmospheric boundary layer (ABL) parameterization scheme for tropical environment might lead to better results with GOCART.

Unsteady separation and vortex shedding from a laminar separation bubble over a bluff body

Boundary layers are subject to favorable and adverse pressure gradients because of both the temporal and spatial components of the pressure gradient. The adverse pressure gradient may cause the flow to separate. In a closed loop unsteady tunnel we have studied the initiation of separation in unsteady flow past a constriction (bluff body) in a channel. We have proposed two important scalings for the time when boundary layer separates. One is based on the local pressure gradient and the other is a convective time scale based on boundary layer parameters. The flow visualization using a dye injection technique shows the flow structure past the body. Nondimensional shedding frequency (Strouhal number) is calculated based on boundary layer and momentum thicknesses. Strouhal number based on the momentum thickness shows a close agreement with that for flat plate and circular cylinder.

Experimental Investigation of Heat Fluxes in the Vicinity of Protuberances on a Flat Plate at Hypersonic Speeds

Heat transfer rates measured in front and to the side of a protrusion on an aluminum flat plate subjected to hypersonic flow at zero angle of attack are presented for two flow enthalpies of approximately 2 MJ/kg and 4.5 MJ/kg. Experiments were conducted in the hypersonic shock tunnel (HST2) and free piston driven HST3 at a freestream Mach number of 8. Heat transfer data was obtained for different geometries of the protrusion of a height of 4 mm, which is approximately the local boundary layer thickness. Comparatively high rates of heat transfer were obtained at regions of flow circulation in the separated region, with the hottest spot generally appearing in front of the protuberance. Experimental values showed moderate agreement with existing empirical correlations at higher enthalpy but not at all for the lower enthalpy condition, although the correlations were coined at enthalpy values nearer to the lower value. Schlieren visualization was also done to investigate the flow structures qualitatively.

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

Transonic flutter, shocks, and parameter effects

There is a drop in the flutter boundary of an aeroelastic system placed in a transonic flow due to compressibility effects and is known as the transonic dip. Viscous effects can shift the lo-cation of the shock and depending on the shock strength the boundary layer may separate leading to changes in the flutter speed. An unsteady Euler flow solver coupled with the structural dynamic equations is used to understand the effect of shock on the transonic dip. The effect of various system parameters such as mass ratio, location of the center of mass, position of the elastic axis, ratio of uncoupled natural frequencies in heave and pitch are also studied. Steady turbulent flow results are presented to demonstrate the effect of viscosity on the location and strength of the shock.

Discrete time second order sliding mode observer for uncertain linear multi-output system

This paper presents a second order sliding mode observer (SOSMO) design for discrete time uncertain linear multi-output system. The design procedure is effective for both matched and unmatched bounded uncertainties and/or disturbances. A second order sliding function and corresponding sliding manifold for discrete time system are defined similar to the lines of continuous time counterpart. A boundary layer concept is employed to avoid switching across the defined sliding manifold and the sliding trajectory is confined to a boundary layer once it converges to it. The condition for existence of convergent quasi-sliding mode (QSM) is derived. The observer estimation errors satisfying given stability conditions converge to an ultimate finite bound (within the specified boundary layer) with thickness O(T-2) where T is the sampling period. A relation between sliding mode gain and boundary layer is established for the existence of second order discrete sliding motion. The design strategy is very simple to apply and is demonstrated for three examples with different class of disturbances (matched and unmatched) to show the effectiveness of the design. Simulation results to show the robustness with respect to the measurement noise are given for SOSMO and the performance is compared with pseudo-linear Kalman filter (PLKF). (C) 2013 Published by Elsevier Ltd. on behalf of The Franklin Institute

Effect of angle of incidence on mixed convective wake dynamics and heat transfer past a square cylinder in cross flow at Re=100

In this paper, a numerical investigation is performed to study the mixed convective flow and heat transfer characteristics past a square cylinder in cross flow at incidence. Utilizing air (Pr = 0.71) as an operating fluid, computations are carried out at a representative Reynolds number (Re) of 100. Angles of incidences are varied as, 0 degrees <= alpha <= 45 degrees. Effect of superimposed positive and negative cross-flow buoyancy is brought about by varying the Richardson number (RI) in the range -1.0 <= Ri <= 1.0. The detail features of flow topology and heat transport are analyzed critically for different angles of incidences. The thermo fluidic forces acting on the cylinder during mixed convection are captured in terms of the drag (C-D), lift (C-L), and moment (C-M) coefficients. The results show that the lateral width of the cylinder wake reduces with increasing alpha and the isotherms spread out far wide. In the range 0 degrees < alpha < 45 degrees, C-D reduces with increasing Ri. The functional dependence of C-M with Ri reveals a linear relationship. Thermal boundary layer thickness reduces with increasing angle of incidences. The global rate of heat transfer from the cylinder increases with increasing alpha. (C) 2014 Elsevier Ltd. All rights reserved.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Performance of WRF-Chem over Indian region: Comparison with measurements

The aerosol mass concentrations over several Indian regions have been simulated using the online chemistry transport model, WRF-Chem, for two distinct seasons of 2011, representing the pre-monsoon (May) and post-monsoon (October) periods during the Indo-US joint experiment `Ganges Valley Aerosol Experiment (GVAX)'. The simulated values were compared with concurrent measurements. It is found that the model systematically underestimates near-surface BC mass concentrations as well as columnar Aerosol Optical Depths (AODs) from the measurements. Examining this in the light of the model-simulated meteorological parameters, we notice the model overestimates both planetary boundary layer height (PBLH) and surface wind speeds, leading to deeper mixing and dispersion and hence lower surface concentrations of aerosols. Shortcoming in simulating rainfall pattern also has an impact through the scavenging effect. It also appears that the columnar AODs are influenced by the unrealistic emission scenarios in the model. Comparison with vertical profiles of BC obtained from aircraft-based measurements also shows a systematic underestimation by the model at all levels. It is seen that concentration of other aerosols, viz., dust and sea-salt are closely linked with meteorological conditions prevailing over the region. Dust is higher during pre-monsoon periods due to the prevalence of north-westerly winds that advect dust from deserts of west Asia into the Indo-Gangetic plain. Winds and rainfall influence sea-salt concentrations. Thus, the unrealistic simulation of wind and rainfall leads to model simulated dust and sea-salt also to deviate from the real values; which together with BC also causes underperformance of the model with regard to columnar AOD. It appears that for better simulations of aerosols over Indian region, the model needs an improvement in the simulation of the meteorology.

Influence of fuel hydrogen fraction on syngas fueled SI engine: Fuel thermo-physical property analysis and in-cylinder experimental investigations

Hydrogen, either in pure form or as a gaseous fuel mixture specie enhances the fuel conversion efficiency and reduce emissions in an internal combustion engine. This is due to the reduction in combustion duration attributed to higher laminar flame speeds. Hydrogen is also expected to increase the engine convective heat flux, attributed (directly or indirectly) to parameters like higher adiabatic flame temperature, laminar flame speed, thermal conductivity and diffusivity and lower flame quenching distance. These factors (adversely) affect the thermo-kinematic response and offset some of the benefits. The current work addresses the influence of mixture hydrogen fraction in syngas on the engine energy balance and the thermo-kinematic response for close to stoichiometric operating conditions. Four different bio-derived syngas compositions with fuel calorific value varying from 3.14 MJ/kg to 7.55 MJ/kg and air fuel mixture hydrogen fraction varying from 7.1% to 14.2% by volume are used. The analysis comprises of (a) use of chemical kinetics simulation package CHEMKIN for quantifying the thermo-physical properties (b) 0-D model for engine in-cylinder analysis and (c) in-cylinder investigations on a two-cylinder engine in open loop cooling mode for quantifying the thermo-kinematic response and engine energy balance. With lower adiabatic flame temperature for Syngas, the in-cylinder heat transfer analysis suggests that temperature has little effect in terms of increasing the heat flux. For typical engine like conditions (700 K and 25 bar at CR of 10), the laminar flame speed for syngas exceeds that of methane (55.5 cm/s) beyond mixture hydrogen fraction of 11% and is attributed to the increase in H based radicals. This leads to a reduction in the effective Lewis number and laminar flame thickness, potentially inducing flame instability and cellularity. Use of a thermodynamic model to assess the isolated influence of thermal conductivity and diffusivity on heat flux suggests an increase in the peak heat flux between 2% and 15% for the lowest (0.420 MW/m(2)) and highest (0.480 MW/m(2)) hydrogen containing syngas over methane (0.415 MW/m(2)) fueled operation. Experimental investigations indicate the engine cooling load for syngas fueled engine is higher by about 7% and 12% as compared to methane fueled operation; the losses are seen to increase with increasing mixture hydrogen fraction. Increase in the gas to electricity efficiency is observed from 18% to 24% as the mixture hydrogen fraction increases from 7.1% to 9.5%. Further increase in mixture hydrogen fraction to 14.2% results in the reduction of efficiency to 23%; argued due to the changes in the initial and terminal stages of combustion. On doubling of mixture hydrogen fraction, the flame kernel development and fast burn phase duration decrease by about 7% and 10% respectively and the terminal combustion duration, corresponding to 90%-98% mass burn, increases by about 23%. This increase in combustion duration arises from the cooling of the near wall mixture in the boundary layer attributed to the presence of hydrogen. The enhancement in engine cooling load and subsequent reduction in the brake thermal efficiency with increasing hydrogen fraction is evident from the engine energy balance along with the cumulative heat release profiles. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

EFFECT OF PLUME DYNAMICS FOR HEAT TRANSFER ENHANCEMENT AT SOLID-LIQUID INTERFACE

The present paper analyzes the effects of plumes for heat transfer enhancement at solid-liquid interface taking both smooth and grooved surfaces. The experimental setup consists of a tank of dimensions 265 x 265 x 300 (height) containing water. The bottom surface was heated and free surface of the water was left open to the ambient. In the experiments, the bottom plate had either a smooth surface or a grooved surface. We used 90 V-grooved rough surfaces with two groove heights, 10mm and 3mm. The experiment was done with water layer depths of 90mm and 140mm, corresponding to values of aspect ratio(AR) equal to 2.9 and 1.8 respectively. Thymol blue, a pH sensitive dye, was used to visualize the flow near the heated plate. The measured heat transfer coefficients over the grooved surfaces were higher compared that over the smooth surface. The enhanced heat transport in the rough cavities cannot be ascribed to the increase in the contact area, rather it must be the local dynamics of the thermal boundary layer that changes the heat transport over the rough surface.

EXPERIMENTAL INVESTIGATION OF UNSTEADY FLOW IN A TRANSONIC UNI-STAGE AXIAL COMPRESSOR

Detailed steady and unsteady experimental measurements and analysis were performed on a Single stage Transonic Axial Compressor with asymmetric rotor tip clearance for studying the compressor stall phenomena. The installed compressor had asymmetric tip clearance around the rotor casing varying from about 0.65mm to 1.25mm. A calibrated 5-hole aerodynamic probe was traversed radially at exit of rotor and showed the characteristics of increased flow angle at lower mass flow rates for all the speeds. Mach number distribution and boundary layer effects were also clearly captured. Unsteady measurements for velocity were carried out to study the stall cell behavior using a single component calibrated hotwire probe oriented in axial and tangential directions for choke/free flow and near stall conditions. The hotwire probe was traversed radially across the annulus at inlet to the compressor and showed that the velocity fluctuations were dissimilar when probe was aligned axial and tangential to the flow. Averaged velocities across the annulus showed the reduction in velocity as stall was approached. Axial mean flow velocity decreased across the annulus for all the speeds investigated. Tangential velocity at free flow condition was higher at the tip region due to larger radius. At stall condition, the tangential velocity showed decreased velocities at the tip and slightly increased velocities at the hub section indicating that the flow has breakdown at the tip region of the blade and fluid is accelerated below the blockage zone. The averaged turbulent intensity in axial and tangential flow directions increased from free flow to stall condition for all compressor rated speeds. Fast Fourier Transform (FFT) of the raw signals at stall flow condition showed stall cell and its corresponding frequency of occurrence. The stalling frequency of about half of rotational speed of the rotor along with large tip clearance suggests that modal type stall inception was occurring.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Shock tunnel measurements of surface pressures in shock induced separated flow field using MEMS sensor array

Characterized not just by high Mach numbers, but also high flow total enthalpies-often accompanied by dissociation and ionization of flowing gas itself-the experimental simulation of hypersonic flows requires impulse facilities like shock tunnels. However, shock tunnel simulation imposes challenges and restrictions on the flow diagnostics, not just because of the possible extreme flow conditions, but also the short run times-typically around 1 ms. The development, calibration and application of fast response MEMS sensors for surface pressure measurements in IISc hypersonic shock tunnel HST-2, with a typical test time of 600 mu s, for the complex flow field of strong (impinging) shock boundary layer interaction with separation close to the leading edge, is delineated in this paper. For Mach numbers 5.96 (total enthalpy 1.3 MJ kg(-1)) and 8.67 (total enthalpy 1.6 MJ kg(-1)), surface pressures ranging from around 200 Pa to 50 000 Pa, in various regions of the flow field, are measured using the MEMS sensors. The measurements are found to compare well with the measurements using commercial sensors. It was possible to resolve important regions of the flow field involving significant spatial gradients of pressure, with a resolution of 5 data points within 12 mm in each MEMS array, which cannot be achieved with the other commercial sensors. In particular, MEMS sensors enabled the measurement of separation pressure (at Mach 8.67) near the leading edge and the sharply varying pressure in the reattachment zone.

Meridional gradients in aerosol vertical distribution over Indian Mainland: Observations and model simulations

Multi-year observations from the network of ground-based observatories (ARFINET), established under the project `Aerosol Radiative Forcing over India' (ARFI) of Indian Space Research Organization and space-borne lidar `Cloud Aerosol Lidar with Orthogonal Polarization' (CALIOP) along with simulations from the chemical transport model `Goddard Chemistry Aerosol Radiation and Transport' (GOCART), are used to characterize the vertical distribution of atmospheric aerosols over the Indian landmass and its spatial structure. While the vertical distribution of aerosol extinction showed higher values close to the surface followed by a gradual decrease at increasing altitudes, a strong meridional increase is observed in the vertical spread of aerosols across the Indian region in all seasons. It emerges that the strong thermal convections cause deepening of the atmospheric boundary layer, which although reduces the aerosol concentration at lower altitudes, enhances the concentration at higher elevations by pumping up more aerosols from below and also helping the lofted particles to reach higher levels in the atmosphere. Aerosol depolarization ratios derived from CALIPSO as well as the GOCART simulations indicate the dominance of mineral dust aerosols during spring and summer and anthropogenic aerosols in winter. During summer monsoon, though heavy rainfall associated with the Indian monsoon removes large amounts of aerosols, the prevailing southwesterly winds advect more marine aerosols over to landmass (from the adjoining oceans) leading to increase in aerosol loading at lower altitudes than in spring. During spring and summer months, aerosol loading is found to be significant, even at altitudes as high as 4 km, and this is proposed to have significant impacts on the regional climate systems such as Indian monsoon. (C) 2015 Elsevier Ltd. All rights reserved.