Distributed Control and Optimization for Communication and Power Systems

We are at the cusp of a historic transformation of both communication system and electricity system. This creates challenges as well as opportunities for the study of networked systems. Problems of these systems typically involve a huge number of end points that require intelligent coordination in a distributed manner. In this thesis, we develop models, theories, and scalable distributed optimization and control algorithms to overcome these challenges.

This thesis focuses on two specific areas: multi-path TCP (Transmission Control Protocol) and electricity distribution system operation and control. Multi-path TCP (MP-TCP) is a TCP extension that allows a single data stream to be split across multiple paths. MP-TCP has the potential to greatly improve reliability as well as efficiency of communication devices. We propose a fluid model for a large class of MP-TCP algorithms and identify design criteria that guarantee the existence, uniqueness, and stability of system equilibrium. We clarify how algorithm parameters impact TCP-friendliness, responsiveness, and window oscillation and demonstrate an inevitable tradeoff among these properties. We discuss the implications of these properties on the behavior of existing algorithms and motivate a new algorithm Balia (balanced linked adaptation) which generalizes existing algorithms and strikes a good balance among TCP-friendliness, responsiveness, and window oscillation. We have implemented Balia in the Linux kernel. We use our prototype to compare the new proposed algorithm Balia with existing MP-TCP algorithms.

Our second focus is on designing computationally efficient algorithms for electricity distribution system operation and control. First, we develop efficient algorithms for feeder reconfiguration in distribution networks. The feeder reconfiguration problem chooses the on/off status of the switches in a distribution network in order to minimize a certain cost such as power loss. It is a mixed integer nonlinear program and hence hard to solve. We propose a heuristic algorithm that is based on the recently developed convex relaxation of the optimal power flow problem. The algorithm is efficient and can successfully computes an optimal configuration on all networks that we have tested. Moreover we prove that the algorithm solves the feeder reconfiguration problem optimally under certain conditions. We also propose a more efficient algorithm and it incurs a loss in optimality of less than 3% on the test networks.

Second, we develop efficient distributed algorithms that solve the optimal power flow (OPF) problem on distribution networks. The OPF problem determines a network operating point that minimizes a certain objective such as generation cost or power loss. Traditionally OPF is solved in a centralized manner. With increasing penetration of volatile renewable energy resources in distribution systems, we need faster and distributed solutions for real-time feedback control. This is difficult because power flow equations are nonlinear and kirchhoff's law is global. We propose solutions for both balanced and unbalanced radial distribution networks. They exploit recent results that suggest solving for a globally optimal solution of OPF over a radial network through a second-order cone program (SOCP) or semi-definite program (SDP) relaxation. Our distributed algorithms are based on the alternating direction method of multiplier (ADMM), but unlike standard ADMM-based distributed OPF algorithms that require solving optimization subproblems using iterative methods, the proposed solutions exploit the problem structure that greatly reduce the computation time. Specifically, for balanced networks, our decomposition allows us to derive closed form solutions for these subproblems and it speeds up the convergence by 1000x times in simulations. For unbalanced networks, the subproblems reduce to either closed form solutions or eigenvalue problems whose size remains constant as the network scales up and computation time is reduced by 100x compared with iterative methods.

Control mechanisms in the human binocular oculomotor system

A study of human eye movements was made in order to elucidate the nature of the control mechanism in the binocular oculomotor system.

We first examined spontaneous eye movements during monocular and binocular fixation in order to determine the corrective roles of flicks and drifts. It was found that both types of motion correct fixational errors, although flicks are somewhat more active in this respect. Vergence error is a stimulus for correction by drifts but not by flicks, while binocular vertical discrepancy of the visual axes does not trigger corrective movements.

Second, we investigated the non-linearities of the oculomotor system by examining the eye movement responses to point targets moving in two dimensions in a subjectively unpredictable manner. Such motions consisted of hand-limited Gaussian random motion and also of the sum of several non-integrally related sinusoids. We found that there is no direct relationship between the phase and the gain of the oculomotor system. Delay of eye movements relative to target motion is determined by the necessity of generating a minimum afferent (input) signal at the retina in order to trigger corrective eye movements. The amplitude of the response is a function of the biological constraints of the efferent (output) portion of the system: for target motions of narrow bandwidth, the system responds preferentially to the highest frequency; for large bandwidth motions, the system distributes the available energy equally over all frequencies. Third, the power spectra of spontaneous eye movements were compared with the spectra of tracking eye movements for Gaussian random target motions of varying bandwidths. It was found that there is essentially no difference among the various curves. The oculomotor system tracks a target, not by increasing the mean rate of impulses along the motoneurons of the extra-ocular muscles, but rather by coordinating those spontaneous impulses which propagate along the motoneurons during stationary fixation. Thus, the system operates at full output at all times.

Fourth, we examined the relative magnitude and phase of motions of the left and the right visual axes during monocular and binocular viewing. We found that the two visual axes move vertically in perfect synchronization at all frequencies for any viewing condition. This is not true for horizontal motions: the amount of vergence noise is highest for stationary fixation and diminishes for tracking tasks as the bandwidth of the target motion increases. Furthermore, movements of the occluded eye are larger than those of the seeing eye in monocular viewing. This effect is more pronounced for horizontal motions, for stationary fixation, and for lower frequencies.

Finally, we have related our findings to previously known facts about the pertinent nerve pathways in order to postulate a model for the neurological binocular control of the visual axes.

Uptake and factors that influence the use of 'sit less, move more' occupational intervention strategies in Spanish office employees

Background: Little is known about the types of 'sit less, move more' strategies that appeal to office employees, or what factors influence their use. This study assessed the uptake of strategies in Spanish university office employees engaged in an intervention, and those factors that enabled or limited strategy uptake. Methods: The study used a mixed method design. Semi-structured interviews were conducted with academics and administrators (n = 12; 44 +/- 12 mean SD age; 6 women) at three points across the five-month intervention, and data used to identify factors that influenced the uptake of strategies. Employees who finished the intervention then completed a survey rating (n = 88; 42 +/- 8 mean SD age; 51 women) the extent to which strategies were used [never (1) to usually (4)]; additional survey items (generated from interviewee data) rated the impact of factors that enabled or limited strategy uptake [no influence (1) to very strong influence (4)]. Survey score distributions and averages were calculated and findings triangulated with interview data. Results: Relative to baseline, 67% of the sample increased step counts post intervention (n = 59); 60% decreased occupational sitting (n = 53). 'Active work tasks' and 'increases in walking intensity' were the strategies most frequently used by employees (89% and 94% sometimes or usually utilised these strategies); 'walk-talk meetings' and ` lunchtime walking groups' were the least used (80% and 96% hardly ever or never utilised these strategies). 'Sitting time and step count logging' was the most important enabler of behaviour change (mean survey score of 3.1 +/- 0.8); interviewees highlighted the motivational value of being able to view logged data through visual graphics in a dedicated website, and gain feedback on progress against set goals. 'Screen based work' (mean survey score of 3.2 +/- 0.8) was the most significant barrier limiting the uptake of strategies. Inherent time pressures and cultural norms that dictated sedentary work practices limited the adoption of 'walk-talk meetings' and ` lunch time walking groups'. Conclusions: The findings provide practical insights into which strategies and influences practitioners need to target to maximise the impact of 'sit less, move more' occupational intervention strategies.

Efetividade da aplicação semestral de verniz fluoretado na redução da incidência de cárie dentária em pré-escolares: resultados após 18 meses

O objetivo deste trabalho foi avaliar se, após 18 meses de acompanhamento, a aplicação semestral de verniz fluoretado em pré-escolares reduz a incidência de cárie na dentição decídua. Os possíveis efeitos adversos e a aceitabilidade da intervenção também foram investigados. Desenho do estudo: ensaio clínico controlado, randomizado, duplo-cego, com uso de placebo. Ingressaram no estudo 200 crianças na faixa etária de 12 a 48 meses, recrutadas em uma unidade da rede pública de saúde. Cem crianças foram alocadas no grupo teste (verniz fluoretado Duraphat, Colgate Oral Pharmaceuticals, New York, NY, USA) e 100 no grupo controle (verniz placebo, Colgate Oral Pharmaceuticals, New York, NY, USA). Os exames clínicos, seguidos das aplicações semestrais de verniz nos grupos teste e controle, foram realizados por duas odontopediatras treinadas (Kappa=0,85) que utilizaram o International Caries Detection and Assessment System II (ICDAS II) para registrar o índice ceo-s nos níveis c1 (lesão ativa em esmalte não cavitado, lesão em esmalte cavitada e lesão em dentina) e c3 (lesão em dentina) de detecção de doença. Verificou-se também a ocorrência de efeitos adversos e outras queixas por entrevistas telefônicas com os responsáveis 7 a 10 dias após a primeira aplicação de verniz. Como resultado aos 18 meses de acompanhamento, 85 crianças do grupo teste e 92 do grupo controle foram examinadas. As diferenças nos incrementos médios de cárie, entre os grupos teste e controle, foram -0,9 (IC95%: -2,4;0,4) no nível c1 e -0,8 (IC95%:-1,9;0,2) no nível c3. O risco de desenvolvimento de novas lesões de cárie em dentina foi 30% (RR=0,70; IC 95%:0,49;0,99) menor nas crianças do grupo teste em relação às do grupo controle. Foram feitas entrevistas com 183 responsáveis sobre os efeitos adversos e a aceitabilidade da intervenção. Um relatou que a sua criança sentiu ardência na cavidade bucal após a aplicação do placebo e outro ficou incomodado com a coloração amarelada dos dentes da criança imediatamente após a aplicação do verniz fluoretado. A aplicação semestral de verniz fluoretado na dentição decídua de pré-escolares reduz o risco de desenvolver lesão de cárie dentária em dentina. Além disso, essa intervenção é segura e bem aceita pelos pais e pelas crianças.

Influência do comportamento semi-rígido de placas de base e de ligações viga-coluna na resposta dinâmica de pórticos de aço.

Tradicionalmente, na análise e dimensionamento de estruturas de aço, assume-se que as ligações viga-coluna são rígidas ou flexíveis (rotuladas). Por outro lado, é de conhecimento geral que a grande maioria das ligações viga-coluna apresenta um comportamento intermediário, ou seja, semi-rígido. Inúmeros trabalhos de pesquisa têm sido desenvolvidos nos últimos vinte e cinco anos, de forma a estudar o comportamento desse tipo de ligação. Um dos principais objetivos desta investigação é o de propor uma metodologia de análise que represente de forma apropriada a influência do comportamento semi-rígido de placas de base e de ligações viga-coluna, sobre a resposta dinâmica (linear e não-linear) de estruturas de aço. Outra contribuição desta dissertação diz respeito à investigação do comportamento dinâmico (linear e não-linear) de pórticos de aço, a partir da consideração de ligações viga-coluna simétricas e não-simétricas e especialmente das placas de base. A análise estrutural é desenvolvida com base no emprego do programa de elementos finitos ANSYS [27]. Nos modelos em elementos finitos foram considerados os efeitos de não-linearidade geométrica (efeitos de segunda ordem), o comportamento não-linear das placas de base e das ligações viga-coluna e, bem como, o efeito de histerese que ocorre quando a estrutura é submetida a cargas cíclicas. Os resultados alcançados indicaram que o fenômeno físico da ressonância não ocorre no que se refere à resposta dinâmica dos modelos semi-rígidos não-lineares. A ressonância não ocorre na resposta dos modelos devido ao fato de que, na análise dinâmica não-linear, o efeito de histerese presente nas ligações (placas de base e viga-coluna), essencialmente com comportamento não-linear, provoca um amortecimento na resposta dinâmica da estrutura.

The cost of quinine Cinchona pubescens control on Santa Cruz Island, Galapagos

We analyse the cost of controlling the invasive quinine tree Cinchona pubescens Vahl in the highlands of Santa Cruz Island, Galapagos. Control costs in ten 400 m2 plots formed the basis for estimating the cost of control over the whole island. In the plots, densities were 2100–24,000 stems/ha (stems >150 cm tall) and 55,000–138,000 stems/ha (all size classes combined). Control involved uprooting small plants, and applying of a mix of metsulfuron methyl and picloram to cut stumps or to machete cuts in the bark of larger trees. These methods are presently used by Galapagos National Park field crews to control quinine. Costs (in man hours, herbicide and US$) were related to stem density; the density of stems summed across four height classes was a better predictor of costs than density of any one size class. Regressions (on all size classes combined) formed the basis for predictive models of costs. Costs ranged from $14 to $2225 per ha depending on stem density. The amount of herbicide (active ingredient/ha) that must be applied to high density stands of quinine is higher than typical rates of application in an agricultural setting. The cost of treating all existing plants once across quinine’s known range on Santa Cruz Island (c. 11,000 ha) was estimated at c. US$1.65 million. CDF Contribution Number 1013.

Model-based control of a Rijke tube combustion instability

A Rijke tube is used to demonstrate model-based control of a combustion instability, where controller design is based on measurement of the unstable system. The Rijke tube used was of length 0.75m and had a grid-stabilised laminar flame in its lower half. A microphone was used as a sensor and a loudspeaker as an actuator for active control. The open loop transfer function (OLTF) required for controller design was that from the actuator to the sensor. This was measured experimentally by sending a signal with two components to the actuator. The first was a control component from an empirically designed controller, which was used to stabilise the system, thus eliminating the non-linear limit cycle. The second was a high bandwidth signal for identification of the OLTF. This approach to measuring the OLTF is generic and can be applied to large-scale combustors. The measured OLTF showed that only the fundamental mode of the tube was unstable; this was consistent with the OLTF predicted by a mathematical model of the tube, involving 1-D linear acoustic waves and a time delay heat release model. Based on the measured OLTF, a controller to stabilise the instability was designed using Nyquist techniques. This was implemented and was seen to result in an 80dB reduction in the microphone pressure spectrum. A robustness study was performed by adding an additional length to the top of the Rijke tobe. The controller was found to achieve control up to an increase in tube length of 19%. This compared favourably with the empirical controller, which lost control for an increase in tube length of less than 3%.

Impedance control reduces instability that arises from motor noise.

There is ample evidence that humans are able to control the endpoint impedance of their arms in response to active destabilizing force fields. However, such fields are uncommon in daily life. Here, we examine whether the CNS selectively controls the endpoint impedance of the arm in the absence of active force fields but in the presence of instability arising from task geometry and signal-dependent noise (SDN) in the neuromuscular system. Subjects were required to generate forces, in two orthogonal directions, onto four differently curved rigid objects simulated by a robotic manipulandum. The endpoint stiffness of the limb was estimated for each object curvature. With increasing curvature, the endpoint stiffness increased mainly parallel to the object surface and to a lesser extent in the orthogonal direction. Therefore, the orientation of the stiffness ellipses did not orient to the direction of instability. Simulations showed that the observed stiffness geometries and their pattern of change with instability are the result of a tradeoff between maximizing the mechanical stability and minimizing the destabilizing effects of SDN. Therefore, it would have been suboptimal to align the stiffness ellipse in the direction of instability. The time course of the changes in stiffness geometry suggests that modulation takes place both within and across trials. Our results show that an increase in stiffness relative to the increase in noise can be sufficient to reduce kinematic variability, thereby allowing stiffness control to improve stability in natural tasks.

The direct field boundary impedance of two-dimensional periodic structures with application to high frequency vibration prediction.

Large sections of many types of engineering construction can be considered to constitute a two-dimensional periodic structure, with examples ranging from an orthogonally stiffened shell to a honeycomb sandwich panel. In this paper, a method is presented for computing the boundary (or edge) impedance of a semi-infinite two-dimensional periodic structure, a quantity which is referred to as the direct field boundary impedance matrix. This terminology arises from the fact that none of the waves generated at the boundary (the direct field) are reflected back to the boundary in a semi-infinite system. The direct field impedance matrix can be used to calculate elastic wave transmission coefficients, and also to calculate the coupling loss factors (CLFs), which are required by the statistical energy analysis (SEA) approach to predicting high frequency vibration levels in built-up systems. The calculation of the relevant CLFs enables a two-dimensional periodic region of a structure to be modeled very efficiently as a single subsystem within SEA, and also within related methods, such as a recently developed hybrid approach, which couples the finite element method with SEA. The analysis is illustrated by various numerical examples involving stiffened plate structures.