Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part I. Mathematical modeling

The conventional Newton and fast decoupled power flow (FDPF) methods have been considered inadequate to obtain the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. It is well known that the PV and Q-theta decoupling assumptions of the fast decoupled power flow formulation no longer hold in the vicinity of the critical point. Moreover, the Jacobian matrix of the Newton method becomes singular at this point. However, the maximum loading point can be efficiently computed through parameterization techniques of continuation methods. In this paper it is shown that by using either theta or V as a parameter, the new fast decoupled power flow versions (XB and BX) become adequate for the computation of the maximum loading point only with a few small modifications. The possible use of reactive power injection in a selected PV bus (Q(PV)) as continuation parameter (mu) for the computation of the maximum loading point is also shown. A trivial secant predictor, the modified zero-order polynomial which uses the current solution and a fixed increment in the parameter (V, theta, or mu) as an estimate for the next solution, is used in predictor step. These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained with the new approach for the IEEE test systems (14, 30, 57 and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that parameters can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. (C) 2003 Elsevier B.V. All rights reserved.

Mathematical modeling of the ester oil-refrigerant R134A mixture two-phase flow with foam formation through a small diameter tube

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

Mathematical glimpse on the Y chromosome degeneration

The Y chromosomes are genetically degenerate and do not recombine with their matching partners X. Non-recombination of XY pairs has been pointed out as the key factor for the degeneration of the Y chromosome. The aim here is to show that there is a mathematical asymmetry in sex chromosomes which leads to the degeneration of Y chromosomes even in the absence of XX and XY recombination. A model for sex-chromosome evolution in a stationary regime is proposed. The consequences of their asymmetry are analyzed and lead us to a couple of conclusions. First, Y chromosome degeneration shows up v 2 more often than X chromosome degeneration. Second, if nature prohibits female mortalities from beeing exactly 50%, then Y chromosome degeneration is inevitable.

A mathematical model for wave propagation in elastic tubes with inhomogeneities: Application to blood waves propagation

We study the propagation of waves in an elastic tube filled with an inviscid fluid. We consider the case of inhomogeneity whose mechanical and geometrical properties vary in space. We deduce a system of equations of the Boussinesq type as describing the wave propagation in the tube. Numerical simulations of these equations show that inhomogeneities prevent separation of right-going from left-going waves. Then reflected and transmitted coefficients are obtained in the case of localized constriction and localized rigidity. Next we focus on wavetrains incident on various types of anomalous regions. We show that the existence of anomalous regions modifies the wavetrain patterns. (c) 2007 Elsevier B.V. All rights reserved.

Mathematical models of generalized diffusion

We discuss in this paper equations describing processes involving non-linear and higher-order diffusion. We focus on a particular case (u(t) = 2 lambda (2)(uu(x))(x) + lambda (2)u(xxxx)), which is put into analogy with the KdV equation. A balance of nonlinearity and higher-order diffusion enables the existence of self-similar solutions, describing diffusive shocks. These shocks are continuous solutions with a discontinuous higher-order derivative at the shock front. We argue that they play a role analogous to the soliton solutions in the dispersive case. We also discuss several physical instances where such equations are relevant.

On non-linear dynamics and control designs applied to the ideal and non-ideal variants of the Fitzhugh-Nagumo (FN) mathematical model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mathematical modeling and control of population systems: Applications in biological pest control

The aim of this paper is to apply methods from optimal control theory, and from the theory of dynamic systems to the mathematical modeling of biological pest control. The linear feedback control problem for nonlinear systems has been formulated in order to obtain the optimal pest control strategy only through the introduction of natural enemies. Asymptotic stability of the closed-loop nonlinear Kolmogorov system is guaranteed by means of a Lyapunov function which can clearly be seen to be the solution of the Hamilton-Jacobi-Bellman equation, thus guaranteeing both stability and optimality. Numerical simulations for three possible scenarios of biological pest control based on the Lotka-Volterra models are provided to show the effectiveness of this method. (c) 2007 Elsevier B.V. All rights reserved.

On non-linear dynamics and an optimal control synthesis of the action potential of membranes (ideal and non-ideal cases) of the Hodgkin-Huxley (HH) mathematical model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effectiveness-ntu computation with a mathematical model for cross-flow heat exchangers

Due to the wide range of design possibilities, simple manufactured, low maintenance and low cost, cross-flow heat exchangers are extensively used in the petroleum, petrochemical, air conditioning, food storage, and others industries. In this paper a mathematical model for cross-flow heat exchangers with complex flow arrangements for determining epsilon -NTU relations is presented. The model is based on the tube element approach, according to which the heat exchanger outlet temperatures are obtained by discretizing the coil along the tube fluid path. In each cross section of the element, tube-side fluid temperature is assumed to be constant because the heat capacity rate ratio C*=Cmin/Cmax tends toward zero in the element. Thus temperature is controlled by effectiveness of a local element corresponding to an evaporator or a condenser-type element. The model is validated through comparison with theoretical algebraic relations for single-pass cross-flow arrangements with one or more rows. Very small relative errors are obtained showing the accuracy of the present model. epsilon -NTU curves for several complex circuit arrangements are presented. The model developed represents a useful research tool for theoretical and experimental studies on heat exchangers performance.

Mathematical models for Isoptera (Insecta) mound growth

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Mathematical equation correction to spectral and transport interferences in high-resolution continuum source flame atomic absorption spectrometry: determination of lead in phosphoric acid

Um método de correção de interferência espectral e de transporte é proposto, e foi aplicado para minimizar interferências por moléculas de PO produzidas em chama ar-acetileno e de transporte causada pela variação da concentração de ácido fosfórico. Átomos de Pb e moléculas de PO absorvem a 217,0005 nm, então Atotal217,0005 nm = A Pb217,0005 nm + A PO217,0005 nm. Monitorando o comprimento de onda alternativo de PO em 217,0458 nm, é possível calcular a contribuição relativa de PO na absorbância total a 217,0005 nm: A Pb217,0005 nm = Atotal217,0005 nm - A PO217,0005 nm = Atotal217,0005 nm - k (A PO217,0458 nm). O fator de correção k é a razão entre os coeficientes angulares de duas curvas analíticas para P obtidas a 217,0005 e 217,0458 nm (k = b217,0005 nm/b217,0458 nm). Fixando-se a taxa de aspiração da amostra em 5,0 ml min-1, e integrando-se a absorbância no comprimento de onda a 3 pixels, curvas analíticas para Pb (0,1 - 1,0 mg L-1) foram obtidas com coeficientes de correlação típicos > 0,9990. As correlações lineares entre absorbância e concentração de P nos comprimentos de onda 217,0005 e 217,0458 foram > 0,998. O limite de detecção de Pb foi 10 µg L-1. O método de correção proposto forneceu desvios padrão relativos (n=12) de 2,0 a 6,0%, ligeiramente menores que os obtidos sem correção (1,4-4,3%). As recuperações de Pb adicionado às amostras de ácido fosfórico variaram de 97,5 a 100% (com correção pelo método proposto) e de 105 a 230% (sem correção).

Mathematical tool to size rural digesters

Os biodigestores têm sido objetos de grande destaque devido a atual crise de energia e conseqüente busca de fontes alternativas. Outro fator que coloca os biodigestores em evidência é o intenso processo de modernização da agropecuária, que além da grande demanda de energia, produz um volume de resíduos animais e de culturas, que ocasiona muitas vezes problemas de ordem sanitária. O objetivo deste trabalho é fornecer uma ferramenta matemática para determinação de parâmetros para projetos de construção de biodigestores rurais, levando-se em consideração o atendimento de necessidades energéticas, obedecendo os dimensionamentos dos sistemas, fatores de rendimento e garantindo a funcionalidade. Para isto, foram formulados modelos de otimização não lineares, de fácil resolução, para os três principais tipos de biodigestores rurais. Com a resolução destes modelos são determinados a altura e o diâmetro que levem a um volume mínimo para cada tipo, com isto reduz-se a quantidade necessária de materiais de alvenaria e consequentemente o custo do biodigestor é diminuído.