Robust THP Transceiver Designs for Multiuser MIMO Downlink with Imperfect CSIT

We present robust joint nonlinear transceiver designs for multiuser multiple-input multiple-output (MIMO) downlink in the presence of imperfections in the channel state information at the transmitter (CSIT). The base station (BS) is equipped with multiple transmit antennas, and each user terminal is equipped with one or more receive antennas. The BS employs Tomlinson-Harashima precoding (THP) for interuser interference precancellation at the transmitter. We consider robust transceiver designs that jointly optimize the transmit THP filters and receive filter for two models of CSIT errors. The first model is a stochastic error (SE) model, where the CSIT error is Gaussian-distributed. This model is applicable when the CSIT error is dominated by channel estimation error. In this case, the proposed robust transceiver design seeks to minimize a stochastic function of the sum mean square error (SMSE) under a constraint on the total BS transmit power. We propose an iterative algorithm to solve this problem. The other model we consider is a norm-bounded error (NBE) model, where the CSIT error can be specified by an uncertainty set. This model is applicable when the CSIT error is dominated by quantization errors. In this case, we consider a worst-case design. For this model, we consider robust (i) minimum SMSE, (ii) MSE-constrained, and (iii) MSE-balancing transceiver designs. We propose iterative algorithms to solve these problems, wherein each iteration involves a pair of semidefinite programs (SDPs). Further, we consider an extension of the proposed algorithm to the case with per-antenna power constraints. We evaluate the robustness of the proposed algorithms to imperfections in CSIT through simulation, and show that the proposed robust designs outperform nonrobust designs as well as robust linear transceiver designs reported in the recent literature.

Sensitivity of the simulated precipitation to changes in convective relaxation time scale

The paper describes the sensitivity of the simulated precipitation to changes in convective relaxation time scale (TAU) of Zhang and McFarlane (ZM) cumulus parameterization, in NCAR-Community Atmosphere Model version 3 (CAM3). In the default configuration of the model, the prescribed value of TAU, a characteristic time scale with which convective available potential energy (CAPE) is removed at an exponential rate by convection, is assumed to be 1 h. However, some recent observational findings suggest that, it is larger by around one order of magnitude. In order to explore the sensitivity of the model simulation to TAU, two model frameworks have been used, namely, aqua-planet and actual-planet configurations. Numerical integrations have been carried out by using different values of TAU, and its effect on simulated precipitation has been analyzed. The aqua-planet simulations reveal that when TAU increases, rate of deep convective precipitation (DCP) decreases and this leads to an accumulation of convective instability in the atmosphere. Consequently, the moisture content in the lower-and mid-troposphere increases. On the other hand, the shallow convective precipitation (SCP) and large-scale precipitation (LSP) intensify, predominantly the SCP, and thus capping the accumulation of convective instability in the atmosphere. The total precipitation (TP) remains approximately constant, but the proportion of the three components changes significantly, which in turn alters the vertical distribution of total precipitation production. The vertical structure of moist heating changes from a vertically extended profile to a bottom heavy profile, with the increase of TAU. Altitude of the maximum vertical velocity shifts from upper troposphere to lower troposphere. Similar response was seen in the actual-planet simulations. With an increase in TAU from 1 h to 8 h, there was a significant improvement in the simulation of the seasonal mean precipitation. The fraction of deep convective precipitation was in much better agreement with satellite observations.

Vertical thermal structure of the atmosphere during active and weak phases of convection over the north Bay of Bengal: Observation and model results

Owing to the lack of atmospheric vertical profile data with sufficient accuracy and vertical resolution, the response of the deep atmosphere to passage of monsoon systems over the Bay of Bengal. had not been satisfactorily elucidated. Under the Indian Climate Research Programme, a special observational programme called 'Bay of Bengal Monsoon Experiment' (BOBMEX), was conducted during July-August 1999. The present study is based on the high-resolution radiosondes launched during BOBMEX in the north Bay. Clear changes in the vertical thermal structure of the atmosphere between active and weak phases of convection have been observed. The atmosphere cooled below 6 km height and became warmer between 6 and 13 km height. The warmest layer was located between 8 and 10 km height, and the coldest layer was found just below 5 km height. The largest fluctuations in the humidity field occurred in the mid-troposphere. The observed changes between active and weak phases of convection are compared with the results from an atmospheric general circulation model, which is similar to that used at the National Centre for Medium Range Weather Forecasting, New Delhi. The model is not able to capture realistically some important features of the temperature and humidity profiles in the lower troposphere and in the boundary layer during the active and weak spells.

A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Dynamics of pulsed flow in an elastic tube

Internal haemorrhage, often leading to cardio-vascular arrest happens to be one of the prime sources of high fatality rates in mammals. We propose a simplistic model of fluid flow in our attempt to specify the location of the haemorrhagic spot, which, if located accurately, could possibly be operated leading to an instant cure. The model we employ for the purpose is basically fluid mechanical in origin and consists of a viscous fluid, pumped by a periodic force and flowing through an elastic tube. The analogy is with that of blood, pumped from the heart and flowing through an artery or vein. Our results, aided by graphical illustrations, match reasonably well with experimental observations.

Albedo enhancement of marine clouds to counteract global warming: impacts on the hydrological cycle

Recent studies have shown that changes in solar radiation affect the hydrological cycle more strongly than equivalent CO(2) changes for the same change in global mean surface temperature. Thus, solar radiation management ``geoengineering'' proposals to completely offset global mean temperature increases by reducing the amount of absorbed sunlight might be expected to slow the global water cycle and reduce runoff over land. However, proposed countering of global warming by increasing the albedo of marine clouds would reduce surface solar radiation only over the oceans. Here, for an idealized scenario, we analyze the response of temperature and the hydrological cycle to increased reflection by clouds over the ocean using an atmospheric general circulation model coupled to a mixed layer ocean model. When cloud droplets are reduced in size over all oceans uniformly to offset the temperature increase from a doubling of atmospheric CO(2), the global-mean precipitation and evaporation decreases by about 1.3% but runoff over land increases by 7.5% primarily due to increases over tropical land. In the model, more reflective marine clouds cool the atmospheric column over ocean. The result is a sinking motion over oceans and upward motion over land. We attribute the increased runoff over land to this increased upward motion over land when marine clouds are made more reflective. Our results suggest that, in contrast to other proposals to increase planetary albedo, offsetting mean global warming by reducing marine cloud droplet size does not necessarily lead to a drying, on average, of the continents. However, we note that the changes in precipitation, evaporation and P-E are dominated by small but significant areas, and given the highly idealized nature of this study, a more thorough and broader assessment would be required for proposals of altering marine cloud properties on a large scale.

Unified Vdd Vth Optimization Based DVFM Controller for a Logic Block

In this paper analytical expressions for optimal Vdd and Vth to minimize energy for a given speed constraint are derived. These expressions are based on the EKV model for transistors and are valid in both strong inversion and sub threshold regions. The effect of gate leakage on the optimal Vdd and Vth is analyzed. A new gradient based algorithm for controlling Vdd and Vth based on delay and power monitoring results is proposed. A Vdd-Vth controller which uses the algorithm to dynamically control the supply and threshold voltage of a representative logic block (sum of absolute difference computation of an MPEG decoder) is designed. Simulation results using 65 nm predictive technology models are given.

Evolution of 1/f α noise during electromigration stressing of metal film: Spectral signature of electromigration process

In this paper we report a systematic study of low-frequency 1/fα resistance fluctuation in a metal film at different stages of electromigration. The resistance fluctuation (noise) measurement was carried out in presence of a dc electromigration stressing current. We observe that in addition to the increase in the spectral power SV(f), the frequency dependence of the spectral power changes as the electromigration process progresses and the exponent α starts to change from 1 to higher value closer to 1.5. We interpret this change in α as arising due to an additional contribution to the spectral power with a 1/f3/2 component, which starts to contribute as the electromigration process progresses. This additional component SV(f) ∼ 1/f3/2 has been suggested to originate from long range diffusion that would accompany any electromigration process. The experimental observation finds support in a model simulation, where we also find that the enhancement of noise during electromigration stressing is accompanied by a change in spectral power frequency dependence.

Mechanisms of summer intraseasonal sea surface temperature oscillations in the Bay of Bengal

Intraseasonal variations (ISV) of sea surface temperature (SST) in the Bay of Bengal (BoB) is highest in its northwestern part. An Indian Ocean model forced by QuikSCAT winds and climatological river discharge (QR run) reproduces ISV of SST, albeit with weaker magnitude. Air-sea fluxes, in the presence of a shallow mixed layer, efficiently effect intraseasonal SST fluctuations. Warming during intraseasonal events is smaller (<1°C) for June - July period and larger (1.5° to 2°C) during September, the latter due to a thinner mixed layer. To examine the effect of salinity on ISV, the model was run by artificially increasing the salinity (NORR run) and by decreasing it (MAHA10 run). In NORR, both rainfall and river discharge were switched off and in MAHA10 the discharge by river Mahanadi was increased tenfold. The spatial pattern of ISV as well as its periodicity was similar in QR, NORR and MAHA10. The ISV was stronger in NORR and weaker in MAHA10, compared to QR. In NORR, both intraseasonal warming and cooling were higher than in QR, the former due to reduced air-sea heat loss as the mean SST was lower, and the latter due to enhanced subsurface processes resulting from weaker stratification. In MAHA10, both warming and cooling were lower than in QR, the former due to higher air-sea heat loss owing to higher mean SST, and the latter due to weak subsurface processes resulting from stronger stratification. These model experiments suggest that salinity effects are crucial in determining amplitudes of intraseasonal SST variations in the BoB.

Dependence of climate forcing and response on the altitude of black carbon aerosols

Black carbon aerosols absorb solar radiation and decrease planetary albedo, and thus can contribute to climate warming. In this paper, the dependence of equilibrium climate response on the altitude of black carbon is explored using an atmospheric general circulation model coupled to a mixed layer ocean model. The simulations model aerosol direct and semi-direct effects, but not indirect effects. Aerosol concentrations are prescribed and not interactive. It is shown that climate response of black carbon is highly dependent on the altitude of the aerosol. As the altitude of black carbon increases, surface temperatures decrease; black carbon near the surface causes surface warming, whereas black carbon near the tropopause and in the stratosphere causes surface cooling. This cooling occurs despite increasing planetary absorption of sunlight (i.e. decreasing planetary albedo). We find that the trend in surface air temperature response versus the altitude of black carbon is consistent with our calculations of radiative forcing after the troposphere, stratosphere, and land surface have undergone rapid adjustment, calculated as ``regressed'' radiative forcing. The variation in climate response from black carbon at different altitudes occurs largely from different fast climate responses; temperature dependent feedbacks are not statistically distinguishable. Impacts of black carbon at various altitudes on the hydrological cycle are also discussed; black carbon in the lowest atmospheric layer increases precipitation despite reductions in solar radiation reaching the surface, whereas black carbon at higher altitudes decreases precipitation.

Bose-Hubbard model in a strong effective magnetic field: Emergence of a chiral Mott insulator ground state

Motivated by experiments on Josephson junction arrays, and cold atoms in an optical lattice in a synthetic magnetic field, we study the ``fully frustrated'' Bose-Hubbard model with half a magnetic flux quantum per plaquette. We obtain the phase diagram of this model on a two-leg ladder at integer filling via the density matrix renormalization group approach, complemented by Monte Carlo simulations on an effective classical XY model. The ground state at intermediate correlations is consistently shown to be a chiral Mott insulator (CMI) with a gap to all excitations and staggered loop currents which spontaneously break time-reversal symmetry. We characterize the CMI state as a vortex supersolid or an indirect exciton condensate, and discuss various experimental implications.

Monsoon circulation interaction with Western Ghats orography under changing climate - Projection by a 20-km mesh AGCM

In this study, the authors have investigated the likely future changes in the summer monsoon over the Western Ghats (WG) orographic region of India in response to global warming, using time-slice simulations of an ultra high-resolution global climate model and climate datasets of recent past. The model with approximately 20-km mesh horizontal resolution resolves orographic features on finer spatial scales leading to a quasi-realistic simulation of the spatial distribution of the present-day summer monsoon rainfall over India and trends in monsoon rainfall over the west coast of India. As a result, a higher degree of confidence appears to emerge in many aspects of the 20-km model simulation, and therefore, we can have better confidence in the validity of the model prediction of future changes in the climate over WG mountains. Our analysis suggests that the summer mean rainfall and the vertical velocities over the orographic regions of Western Ghats have significantly weakened during the recent past and the model simulates these features realistically in the present-day climate simulation. Under future climate scenario, by the end of the twenty-first century, the model projects reduced orographic precipitation over the narrow Western Ghats south of 16A degrees N that is found to be associated with drastic reduction in the southwesterly winds and moisture transport into the region, weakening of the summer mean meridional circulation and diminished vertical velocities. We show that this is due to larger upper tropospheric warming relative to the surface and lower levels, which decreases the lapse rate causing an increase in vertical moist static stability (which in turn inhibits vertical ascent) in response to global warming. Increased stability that weakens vertical velocities leads to reduction in large-scale precipitation which is found to be the major contributor to summer mean rainfall over WG orographic region. This is further corroborated by a significant decrease in the frequency of moderate-to-heavy rainfall days over WG which is a typical manifestation of the decrease in large-scale precipitation over this region. Thus, the drastic reduction of vertical ascent and weakening of circulation due to `upper tropospheric warming effect' predominates over the `moisture build-up effect' in reducing the rainfall over this narrow orographic region. This analysis illustrates that monsoon rainfall over mountainous regions is strongly controlled by processes and parameterized physics which need to be resolved with adequately high resolution for accurate assessment of local and regional-scale climate change.

Ship-Mounted Real-Time Surface Observational System on board Indian Vessels for Validation and Refinement of Model Forcing Fields

A network of ship-mounted real-time Automatic Weather Stations integrated with Indian geosynchronous satellites Indian National Satellites (INSATs)] 3A and 3C, named Indian National Centre for Ocean Information Services Real-Time Automatic Weather Stations (I-RAWS), is established. The purpose of I-RAWS is to measure the surface meteorological-ocean parameters and transmit the data in real time in order to validate and refine the forcing parameters (obtained from different meteorological agencies) of the Indian Ocean Forecasting System (INDOFOS). Preliminary validation and intercomparison of analyzed products obtained from the National Centre for Medium Range Weather Forecasting and the European Centre for Medium-Range Weather Forecasts using the data collected from I-RAWS were carried out. This I-RAWS was mounted on board oceanographic research vessel Sagar Nidhi during a cruise across three oceanic regimes, namely, the tropical Indian Ocean, the extratropical Indian Ocean, and the Southern Ocean. The results obtained from such a validation and intercomparison, and its implications with special reference to the usage of atmospheric model data for forcing ocean model, are discussed in detail. It is noticed that the performance of analysis products from both atmospheric models is similar and good; however, European Centre for Medium-Range Weather Forecasts air temperature over the extratropical Indian Ocean and wind speed in the Southern Ocean are marginally better.

An approach for assessing the accuracy of accelerometer responses during axial impact tests

The primary objective of the present study is to show that for the most common configuration of an impactor system, the accelerometer cannot exactly reproduce the dynamic response of a specimen subjected to impact loading. An equivalent Lumped Parameter Model (LPM) of the given impactor set-up has been formulated for assessing the accuracy of an accelerometer mounted in a drop-weight impactor set-up for an axially loaded specimen. A specimen under the impact loading is represented by a non-linear spring of varying stiffness, while the accelerometer is assumed to behave in a linear manner due to its high stiffness. Specimens made of steel, aluminium and fibre-reinforced composite (FRC) are used in the present study. Assuming the force-displacement response obtained in an actual impact test to be the true behaviour of the test specimen, a suitable numerical approach has been used to solve the governing non-linear differential equations of a three degrees-of-freedom (DOF) system in a piece-wise linear manner. The numerical solution of the governing differential equations following an explicit time integration scheme yields an excellent reproduction of the mechanical behaviour of the specimen, consequently confirming the accuracy of the numerical approach. However, the spring representing the accelerometer predicts a response that qualitatively matches the assumed force-displacement response of the test specimen with a perceptibly lower magnitude of load.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.