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.

Influence of crossed electric and quantizing magnetic fields on the Einstein relation in nonlinear optical, optoelectronic and related materials: Simplified theory, relative comparison and suggestion for experimental determination

An attempt is made to study the Einstein relation for the diffusivity-to-mobility ratio (DMR) under crossed fields' configuration in nonlinear optical materials on the basis of a newly formulated electron dispersion law by incorporating the crystal field in the Hamiltonian and including the anisotropies of the effective electron mass and the spin-orbit splitting constants within the framework of kp formalisms. The corresponding results for III-V, ternary and quaternary compounds form a special case of our generalized analysis. The DMR has also been investigated for II-VI and stressed materials on the basis of various appropriate dispersion relations. We have considered n-CdGeAs2, n-Hg1-xCdxTe, n-In1-xGaxAsyP1-y lattice matched to InP, p-CdS and stressed n-InSb materials as examples. The DMR also increases with increasing electric field and the natures of oscillations are totally band structure dependent with different numerical values. It has been observed that the DMR exhibits oscillatory dependences with inverse quantizing magnetic field and carrier degeneracy due to the Subhnikov-de Haas effect. An experimental method of determining the DMR for degenerate materials in the present case has been suggested. (C) 2010 Elsevier B.V. All rights reserved.

Sliding Modes for the Simulation of Mechanical and Electrical Systems Defined by Differential-Algebraic Equations

We offer a technique, motivated by feedback control and specifically sliding mode control, for the simulation of differential-algebraic equations (DAEs) that describe common engineering systems such as constrained multibody mechanical structures and electric networks. Our algorithm exploits the basic results from sliding mode control theory to establish a simulation environment that then requires only the most primitive of numerical solvers. We circumvent the most important requisite for the conventionalsimulation of DAEs: the calculation of a set of consistent initial conditions. Our algorithm, which relies on the enforcement and occurrence of sliding mode, will ensure that the algebraic equation is satisfied by the dynamic system even for inconsistent initial conditions and for all time thereafter. [DOI:10.1115/1.4001904]

Optimal infinite-horizon feedback laws for a general class of constrained discrete-time systems: stability and moving-horizon approximations

Stability results are given for a class of feedback systems arising from the regulation of time-varying discrete-time systems using optimal infinite-horizon and moving-horizon feedback laws. The class is characterized by joint constraints on the state and the control, a general nonlinear cost function and nonlinear equations of motion possessing two special properties. It is shown that weak conditions on the cost function and the constraints are sufficient to guarantee uniform asymptotic stability of both the optimal infinite-horizon and movinghorizon feedback systems. The infinite-horizon cost associated with the moving-horizon feedback law approaches the optimal infinite-horizon cost as the moving horizon is extended.

Exact static polarizabilities of correlated finite model systems

A finite-field method for calculating exact polarizabilities of correlated conjugated model systems within the valence bond (VB) framework is presented. The correlations reduce the polarizabilities from their noninteracting values and extend the range of linearity to higher external fields. The large nonlinear polarizabilities observed in strongly correlated conjugated organic molecules cannot be directly attributed to electron correlations. The method described can be employed to calculate static polarizabilities for any desired state of a correlated system.

Unified Galerkin- and DAE-Based Approximation of Fractional Order Systems

We consider numerical solutions of nonlinear multiterm fractional integrodifferential equations, where the order of the highest derivative is fractional and positive but is otherwise arbitrary. Here, we extend and unify our previous work, where a Galerkin method was developed for efficiently approximating fractional order operators and where elements of the present differential algebraic equation (DAE) formulation were introduced. The DAE system developed here for arbitrary orders of the fractional derivative includes an added block of equations for each fractional order operator, as well as forcing terms arising from nonzero initial conditions. We motivate and explain the structure of the DAE in detail. We explain how nonzero initial conditions should be incorporated within the approximation. We point out that our approach approximates the system and not a specific solution. Consequently, some questions not easily accessible to solvers of initial value problems, such as stability analyses, can be tackled using our approach. Numerical examples show excellent accuracy. DOI: 10.1115/1.4002516]

Compensator for constant relative stability in process control systems

Process control systems are designed for a closed-loop peak magnitude of 2dB, which corresponds to a damping coefficient () of 0.5 approximately. With this specified constraint, the designer should choose and/or design the loop components to maintain a constant relative stability. However, the manipulative variable in almost all chemical processes will be the flow rate of a process stream. Since the gains and the time constants of the process will be functions of the manipulative variable, a constant relative stability cannot be maintained. Up to now, this problem has been overcome either by selecting proper control valve flow characteristics or by gain scheduling of controller parameters. Nevertheless, if a wrong control valve selection is made then one has to account for huge loss in controllability or eventually it may lead to an unstable control system. To overcome these problems, a compensator device that can bring back the relative stability of the control system was proposed. This compensator is similar to a dynamic nonlinear controller that has both online and offline information on several factors related to the control system. The design and analysis of the proposed compensator is discussed in this article. Finally, the performance of the compensator is validated by applying it to a two-tank blending process. It has been observed that by using a compensator in the process control system, the relative stability could be brought back to a great extent despite the effects of changes in manipulative flow rate.

Exact traveling- wave solution for model geophysical systems

Exact traveling-wave solutions of time-dependent nonlinear inhomogeneous PDEs, describing several model systems in geophysical fluid dynamics, are found. The reduced nonlinear ODEs are treated as systems of linear algebraic equations in the derivatives. A variety of solutions are found, depending on the rank of the algebraic systems. The geophysical systems include acoustic gravity waves, inertial waves, and Rossby waves. The solutions describe waves which are, in general, either periodic or monoclinic. The present approach is compared with the earlier one due to Grundland (1974) for finding exact solutions of inhomogeneous systems of nonlinear PDEs.

Exact nonlinear travelling hydromagnetic wave solutions

Exact travelling wave solutions for hydromagnetic waves in an exponentially stratified incompressible medium are obtained. With the help of two integrals it becomes possible to reduce the system of seven nonlinear PDE's to a second order nonlinear ODE which describes an one dimensional harmonic oscillator with a nonlinear friction term. This equation is studied in detail in the phase plane. The travelling waves are periodic only when they propagate either horizontally or vertically. The reduced second order nonlinear differential equation describing the travelling waves in inhomogeneous conducting media has rather ubiquitous nature in that it also appears in other geophysical systems such as internal waves, Rossby waves and topographic Rossby waves in the ocean.

A structure-preserving energy function for stability analysis of AC/DC systems

Direct stability analysis ofAC/DC power systems using a structure-preserving energy function (SPEF) is proposed in this paper. The system model considered retains the load buses thereby enabling the representation of nonlinear voltage dependent loads. TheHVDC system is represented with the same degree of detail as is normally done in transient stability simulation. The converter controllers can be represented by simplified or detailed models. Two or multi-terminalDC systems can be considered. The stability analysis is illustrated with a 3-machine system example and encouraging results have been obtained.

Conditional linearization in nonlinear random vibration

In this paper, an improved probabilistic linearization approach is developed to study the response of nonlinear single degree of freedom (SDOF) systems under narrow-band inputs. An integral equation for the probability density function (PDF) of the envelope is derived. This equation is solved using an iterative scheme. The technique is applied to study the hardening type Duffing's oscillator under narrow-band excitation. The results compare favorably with those obtained using numerical simulation. In particular, the bimodal nature of the PDF for the response envelope for certain parameter ranges is brought out.

Dynamic analysis of voltage instability in AC-DC systems

This paper presents the analysis and study of voltage collapse at any converter bus in A C-DC systems considering the dynamics of DC system. The problem of voltage instability is acute when HVDC links are connected to weak AC systems, the strength determined by short circuit ratio (SCR) at the converter bus. The converter control strategies are important in determining voltage instability. Small signal analysis is used to identify critical modes and evaluate the effect of AC system strength and control parameters. A sample two-terminal DC system is studied and the results compared with those obtained from static analysis. Also, the results obtained from small signal analysis are validated with nonlinear simulation.

Linear and nonlinear optical properties of cumulenes and polyenynes: a model exact study

Model exact static and frequency-dependent polarizabilities, static second hyperpolarizabilities and THG coefficents of cumulenes and polyenynes, calculated within the correlated Pariser-Parr-Pople (PPP) model defined over the pi-framework are reported and compared with the results for the polyenes. It is found that for the same chain length, the polarizabilities and THG coefficients of the cumulenes are largest and those of the polyenynes smallest with the polyenes having an intermediate value. The optical gap of the infinite cumulene is lowest (0.75 eV) and is associated with a low transition dipole moment for an excitation involving transfer of an electron between the two orthogonal conjugated pi-systems. The polyenynes have the largest optical gap (4.37 eV), with the magnitude being nearly independent of the chain length. This excitation involves charge transfer between the conjugated bonds in the terminal triple bond. Chain length and frequency dependence of alpha(ij) and gamma(ijkl) of these systems are also reported. The effect of a heteroatom on the polarizability and THG coefficients of acetylenic systems is also reported. It has been found that the presence of the heteroatom reduces the polarizability and THG coefficients of these systems, an effect opposite to that found in the polyenes and cyanine dyes. This result has been associated with the different nature of the charge transfer in the acetylenic systems.

Relaxation behavior of twin nonlinear optical chromophores: Effect of the spacer length

A new series of twin nonlinear optical (NLO) molecules, having two 4-nitrophenol chromophores that are linked via a flexible polymethylene spacer of varying length [(CH2)(n), n = 1-12], were synthesized. Powder second harmonic generation measurements of these twin samples indicated a pronounced odd-even oscillation, with the odd twins exhibiting a high SHG value while the even ones gave no measurable SH signal. This behavior reflects the crystal packing preferences in such twin NLO systems that have odd and even numbers of atoms linking them - the even ones appear to prefer a centrosymmetric packing arrangement. The orientational/disordering dynamics of these twin NLO molecules, doped in a polymer (poly(methyl methacrylate)) matrix, has also been studied using SHG in electric field poled samples. Interestingly, the maximum attainable SH signal, chi((2)), in, the poled samples also showed an odd-even oscillation; the odd ones again having a higher value of chi((2)) This unprecedented odd-even oscillation in such molecularly doped systems is rationalized as being due to the intrinsically greater ease of a parallel alignment of the two chromophores in the twins with an odd spacer than in those with an even one. Further, the temporal stability of the SHG intensity at 70 degrees C, after the removal of the applied corona, was also studied. The relaxation of all the twin chromophores followed a biexponential decay; the characteristic relaxation time (tau(2)) for the slow decay component suggests that while the twin with a single methylene unit relaxes relatively slowly, the relaxation is significantly faster in cases where n = 2 and 3. In the twins with even longer spacer segments, the relaxation again becomes slower and reaches a saturation value. The observed minimum appears to reflect the interplay of two competing factors that affect the chromophore alignment in such twin systems, namely, the electrostatic repulsion between neighboring oriented dipoles and the intrinsic flexibility of the spacer.

Microwave heating in multiphase systems: evaluation of series solutions

The 1D electric field and heat-conduction equations are solved for a slab where the dielectric properties vary spatially in the sample. Series solutions to the electric field are obtained for systems where the spatial variation in the dielectric properties can be expressed as polynomials. The series solution is used to obtain electric-field distributions for a binary oil-water system where the dielectric properties are assumed to vary linearly within the sample. Using the finite-element method temperature distributions are computed in a three-phase oil, water and rock system where the dielectric properties vary due to the changing oil saturation in the rock. Temperature distributions predicted using a linear variation in the dielectric properties are compared with those obtained using the exact nonlinear variation.