A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Macro-economic analysis of forestry options on carbon sequestration in India

The forestry sector provides a number of climate change mitigation options. Apart from this ecological benefit, it has significant social and economic relevance. Implementation of forestry options requires large investments and sustained long-term planning. Thus there is a need for a detailed analysis of forestry options to understand their implications on stock and flow of carbon, required investments, value of forest wealth, contribution to GNP and livelihood, demand management, employment and foreign trade. There is a need to evaluate the additional spending on forestry by analysing the environmental (particularly carbon abatement), social and economic benefits. The biomass needs for India are expected to increase by two to three times by 2020. Depending upon the forest types, ownership patterns and land use patterns, feasible forestry options are identified. It is found among many supply options to be feasible to meet the 'demand based needs' with a mix of management options, species choices and organisational set up. A comparative static framework is used to analyze the macro-economic impacts. Forestry accounts for 1.84% of GNP in India. It is characterized by significant forward industrial linkages and least backward linkage. Forestry generates about 36 million person years of employment annually. India imports Rs. 15 billion worth of forest based materials annually. Implementation of the demand based forestry options can lead to a number of ecological, economic and institutional changes. The notable ones are: enhancement of C stock from 9578 to 17 094 Mt and a net annual C-sequestration from 73 to 149 Mt after accounting for all emissions; a trebling of the output of forestry sector from Rs. 49 billion to Rs. 146 billion annually; an increase in GDP contribution of forestry from Rs. 32 billion to Rs. 105 billion over a period of 35 years; an increase in annual employment level by 23 million person years, emergence of forestry as a net contributor of foreign exchange through trading of forestry products; and an increase in economic value of forest capital stock by Rs. 7260 billion with a cost benefit analysis showing forestry as a profitable option. Implementation of forestry options calls for an understanding of current forest policies and barriers which are analyzed and a number of policy options are suggested. (C) 1997 Elsevier Science B.V.

Tracking and erosion resistance of silicone rubber nanocomposites under positive and negative dc voltages

ASTM D2303 standard provides a method for evaluating the tracking and erosion resistance of polymeric insulators under ac voltages. In this paper, the above method has been extended for evaluating the performance of the insulators under dc stresses. Tests were conducted on polymeric silicone rubber (SR) insulators under positive and negative dc stresses. Micron sized Alumina trihydrate (uATH) and nano sized Alumina (nALU) were used as fillers in SR matrix to improve the resistance to tracking and erosion. Results suggest that SR composites perform better under negative dc than under positive dc voltages. Eroded mass and leakage current data support the above result. Samples with low concentration of nano alumina fillers performed on par with the samples with large loadings of uATH.

A method of tracking the peak power points for a variable speed wind energy conversion system

In this paper, a method of tracking the peak power in a wind energy conversion system (WECS) is proposed, which is independent of the turbine parameters and air density. The algorithm searches for the peak power by varying the speed in the desired direction. The generator is operated in the speed control mode with the speed reference being dynamically modified in accordance with the magnitude and direction of change of active power. The peak power points in the P-omega curve correspond to dP/domega = 0. This fact is made use of in the optimum point search algorithm. The generator considered is a wound rotor induction machine whose stator is connected directly to the grid and the rotor is fed through back-to-back pulse-width-modulation (PWM) converters. Stator flux-oriented vector control is applied to control the active and reactive current loops independently. The turbine characteristics are generated by a dc motor fed from a commercial dc drive. All of the control loops are executed by a single-chip digital signal processor (DSP) controller TMS320F240. Experimental results show that the performance of the control algorithm compares well with the conventional torque control method.

Tracking analysis of an LMS decision feedback equalizer for a wireless channel

We consider a time varying wireless fading channel, equalized by an LMS Decision Feedback equalizer (DFE). We study how well this equalizer tracks the optimal MMSEDFE (Wiener) equalizer. We model the channel by an Autoregressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, we show via some examples that the LMS equalizer moves close to the instantaneous Wiener filter after initial transience. We also compare the LMS equalizer with the instantaneous optimal DFE (the commonly used Wiener filter) designed assuming perfect previous decisions and computed using perfect channel estimate (we will call it as IDFE). We show that the LMS equalizer outperforms the IDFE almost all the time after initial transience.

Command Tracking in High Performance Aircrafts: A New Dynamic Inversion Design

This paper proposes a new straight forward technique based on dynamic inversion, which is applied for tracking the pilot commands in high performance aircrafts.Pilot commands assumed in longitudinal mode are normal acceleration and total velocity(while roll angle and lateral acceleration are maintained at zero). In lateral mode, roll rate and total velocity are used as pilot commands (while climb rate and lateral acceleration are maintained at zero). Ensuring zero lateral acceleration leads to a better turn co-ordination. A six degree-of-freedom model of F-16 aircraft is used for both control design as well as simulation studies. Promising results are obtained which are found to be superior as compared to an existing approach (which is also based on dynamic inversion). The new approach has two potential benefits, namely reduced oscillatory response and reduced control magnitude. Another advantage of this approach is that it leads to a significant reduction of tuning parameters in the control design process.