Distributed nonlinear optimization algorithms for general network problems

There are a number of large networks which occur in many problems dealing with the flow of power, communication signals, water, gas, transportable goods, etc. Both design and planning of these networks involve optimization problems. The first part of this paper introduces the common characteristics of a nonlinear network (the network may be linear, the objective function may be non linear, or both may be nonlinear). The second part develops a mathematical model trying to put together some important constraints based on the abstraction for a general network. The third part deals with solution procedures; it converts the network to a matrix based system of equations, gives the characteristics of the matrix and suggests two solution procedures, one of them being a new one. The fourth part handles spatially distributed networks and evolves a number of decomposition techniques so that we can solve the problem with the help of a distributed computer system. Algorithms for parallel processors and spatially distributed systems have been described.There are a number of common features that pertain to networks. A network consists of a set of nodes and arcs. In addition at every node, there is a possibility of an input (like power, water, message, goods etc) or an output or none. Normally, the network equations describe the flows amoungst nodes through the arcs. These network equations couple variables associated with nodes. Invariably, variables pertaining to arcs are constants; the result required will be flows through the arcs. To solve the normal base problem, we are given input flows at nodes, output flows at nodes and certain physical constraints on other variables at nodes and we should find out the flows through the network (variables at nodes will be referred to as across variables).The optimization problem involves in selecting inputs at nodes so as to optimise an objective function; the objective may be a cost function based on the inputs to be minimised or a loss function or an efficiency function. The above mathematical model can be solved using Lagrange Multiplier technique since the equalities are strong compared to inequalities. The Lagrange multiplier technique divides the solution procedure into two stages per iteration. Stage one calculates the problem variables % and stage two the multipliers lambda. It is shown that the Jacobian matrix used in stage one (for solving a nonlinear system of necessary conditions) occurs in the stage two also.A second solution procedure has also been imbedded into the first one. This is called total residue approach. It changes the equality constraints so that we can get faster convergence of the iterations.Both solution procedures are found to coverge in 3 to 7 iterations for a sample network.The availability of distributed computer systems — both LAN and WAN — suggest the need for algorithms to solve the optimization problems. Two types of algorithms have been proposed — one based on the physics of the network and the other on the property of the Jacobian matrix. Three algorithms have been deviced, one of them for the local area case. These algorithms are called as regional distributed algorithm, hierarchical regional distributed algorithm (both using the physics properties of the network), and locally distributed algorithm (a multiprocessor based approach with a local area network configuration). The approach used was to define an algorithm that is faster and uses minimum communications. These algorithms are found to converge at the same rate as the non distributed (unitary) case.

Modified governing equations for gas-particle nozzle flows

A modified set of governing equations for gas-particle flows in nozzles is suggested to include the inertial forces acting on the particle phase. The problem of gas-particle flow through a nozzle is solved using a first order finite difference scheme. A suitable stability condition for the numerical scheme for gas-particle flows is defined. Results obtained from the present set of equations are compared with those of the previous set of equations. It is also found that present set of equations give results which are in good agreement with the experimental observation.

From Bathymetry to Bioshields: A Review of Post-Tsunami Ecological Research in India and its Implications for Policy

More than half a decade has passed since the December 26th 2004 tsunami hit the Indian coast leaving a trail of ecological, economic and human destruction in its wake. We reviewed the coastal ecological research carried out in India in the light of the tsunami. In addition, we also briefly reviewed the ecological research in other tsunami affected countries in Asia namely Sri Lanka, Indonesia, Thailand and Maldives in order to provide a broader perspective of ecological research after tsunami. A basic search in ISI Web of Knowledge using keywords ``tsunami'' and ``India'' resulted in 127 peer reviewed journal articles, of which 39 articles were pertaining to ecological sciences. In comparison, Sri Lanka, Indonesia, Thailand and Maldives had, respectively, eight, four, 21 and two articles pertaining to ecology. In India, bioshields received the major share of scientific interest (14 out of 39) while only one study (each) was dedicated to corals, seagrasses, seaweeds and meiofauna, pointing to the paucity of research attention dedicated to these critical ecosystems. We noted that very few interdisciplinary studies looked at linkages between pure/applied sciences and the social sciences in India. In addition, there appears to be little correlation between the limited research that was done and its influence on policy in India. This review points to gap areas in ecological research in India and highlights the lessons learnt from research in other tsunami-affected countries. It also provides guidance on the links between science and policy that are required for effective coastal zone management.

Steady flow of couple stress fluid through tubes of slowly varying cross-sections--application to blood flows

Steady laminar flow of a non-Newtonian fluid based on couple stress fluid theory, through narrow tubes of varying cross-sections has been studied theoretically. Asymptotic solutions are obtained for the basic equations and the expressions for the velocity field and the wall shear stress are derived for a general cross-section. Computation and discussions are carried out for the geometries which occur in the context of physiological flows or in particular blood flows. The tapered tubes and constricted tubes are of special importance. It is observed that increase in certain parameters results in erratic flow behaviour proximal to the constricted areas which is further enhanced by the increase in the geometric parameters. This elucidates the implications of the flow in the development of vascular lesions.

Studies on Unsteady Pressure Fields in the Region of Separating and Reattaching Flows

Experimental studies on the measurement of pressure fields in the region of separating and reattaching flows behind several two-dimensional fore-bodies and one axisymmetric body are reported. In particular, extensive measurements of mean pressure, surface pressure fluctuation, and pressure fluctuation within the flow were made for a series of two-dimensional fore-body shapes consisting of triangular nose with varying included angle. The measurements from different bodies are compared and one of the important findings is that the maximum values of rms pressure fluctuation levels in the shear layer approaching reattachment are almost equal to the maximum value of the surface fluctuation levels.

Biorthogonal series method for Oseen type flows

A biorthogonal series method is developed to solve Oseen type flow problems. The theory leads to a new set of eigenfunctions for a specific class of linear non-selfadjoint operators containing the biharmonic one. These eigenfunctions differ from those given earlier in the literature for the biharmonic operator. The method is applied to the problem of thermocapillary flow in a cylindrical liquid bridge of finite length with axial through flow. Flow and temperature distributions are obtained at leading order of an expansion for small surface tension Reynolds number and Prandtl number. Another related problem considered is that of cylindrical cavity flow. Solutions for both cases are presented in terms of biorthogonal series. The effect of axial through flow on velocity and temperature fields is discussed by numerical evaluation of the truncated analytical series. The presence of axial through flow not only convectively shifts the vortices induced by surface forces in the direction of the through flow, but also moves their centers toward the outer cylindrical boundary. This process can lead to significantly asymmetric flow structures.

Impact of climate change on Indian forests: a dynamic vegetation modeling approach

We make an assessment of the impact of projected climate change on forest ecosystems in India. This assessment is based on climate projections of the Regional Climate Model of the Hadley Centre (HadRM3) and the dynamic global vegetation model IBIS for A2 and B2 scenarios. According to the model projections, 39% of forest grids are likely to undergo vegetation type change under the A2 scenario and 34% under the B2 scenario by the end of this century. However, in many forest dominant states such as Chattisgarh, Karnataka and Andhra Pradesh up to 73%, 67% and 62% of forested grids are projected to undergo change. Net Primary Productivity (NPP) is projected to increase by 68.8% and 51.2% under the A2 and B2 scenarios, respectively, and soil organic carbon (SOC) by 37.5% for A2 and 30.2% for B2 scenario. Based on the dynamic global vegetation modeling, we present a forest vulnerability index for India which is based on the observed datasets of forest density, forest biodiversity as well as model predicted vegetation type shift estimates for forested grids. The vulnerability index suggests that upper Himalayas, northern and central parts of Western Ghats and parts of central India are most vulnerable to projected impacts of climate change, while Northeastern forests are more resilient. Thus our study points to the need for developing and implementing adaptation strategies to reduce vulnerability of forests to projected climate change.

Impact of climate change at species level: a case study of teak in India

In this study, we model the long-term effect of climate change on commercially important teak (Tectona grandis) and its productivity in India. This modelling assessment is based on climate projections of the regional climate model of the Hadley Center (HadRM3) and the dynamic vegetation model, IBIS. According to the model projections, 30% of teak grids in India are vulnerable to climate change under both A2 and B2 SRES scenarios because the future climate may not be optimal for teak at these grids. However, the net primary productivity and biomass are expected to increase because of elevated levels of CO2. Given these directions of likely impacts, it is crucial to further investigate the climate change impacts on teak and incorporate such findings into long-term teak plantation programs. This study also demonstrates the feasibility and limitations of assessing the impact of projected climate change at the species level in the tropics.