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.

Alveolar nitric oxide and its role in pediatric asthma control assessment

Resumen Background: Nitric oxide can be measured at multiple flow rates to determine proximal (maximum airway nitric oxide flux; Jaw(NO)) and distal inflammation (alveolar nitric oxide concentration; CA(NO)). The main aim was to study the association among symptoms, lung function, proximal (maximum airway nitric oxide flux) and distal (alveolar nitric oxide concentration) airway inflammation in asthmatic children treated and not treated with inhaled glucocorticoids. Methods: A cross-sectional study with prospective data collection was carried out in a consecutive sample of girls and boys aged between 6 and 16 years with a medical diagnosis of asthma. Maximum airway nitric oxide flux and alveolar nitric oxide concentration were calculated according to the two-compartment model. In asthmatic patients, the asthma control questionnaire (CAN) was completed and forced spirometry was performed. In controls, differences between the sexes in alveolar nitric oxide concentration and maximum airway nitric oxide flux and their correlation with height were studied. The correlation among the fraction of exhaled NO at 50 ml/s (FENO50), CA(NO), Jaw(NO), forced expiratory volume in 1 second (FEV1) and the CAN questionnaire was measured and the degree of agreement regarding asthma control assessment was studied using Cohen's kappa. Results: We studied 162 children; 49 healthy (group 1), 23 asthmatic participants without treatment (group 2) and 80 asthmatic patients treated with inhaled corticosteroids (group 3). CA(NO) (ppb) was 2.2 (0.1-4.5), 3 (0.2-9.2) and 2.45 (0.1-24), respectively. Jaw(NO) (pl/s) was 516 (98.3-1470), 2356.67 (120-6110) and 1426 (156-11805), respectively. There was a strong association (r = 0.97) between FENO50 and Jaw(NO) and the degree of agreement was very good in group 2 and was good in group 3. There was no agreement or only slight agreement between the measures used to monitor asthma control (FEV1, CAN questionnaire, CA(NO) and Jaw(NO)). Conclusions: The results for CA(NO) and Jaw(NO) in controls were similar to those found in other reports. There was no agreement or only slight agreement among the three measure instruments analyzed to assess asthma control. In our sample, no additional information was provided by CA(NO) and Jaw(NO).

Rapid quantitative detection of Yersinia pestis by lateral-flow immunoassay and up-converting phosphor technology-based biosensor

Up-converting phosphor technology (UPT)-based lateral-flow immunoassay has been developed for quantitative detection of Yersinia pestis rapidly and specifically. In this assay, 400 nm up-converting phosphor particles were used as the reporter. A sandwich immumoassay was employed by using a polyclonal antibody against F1 antigen of Y. pestis immobilized on the nitrocellulose membrane and the same antibody conjugated to the UPT particles. The signal detection of the strips was performed by the UPT-based biosensor that could provide a 980 nm IR laser to excite the phosphor particles, then collect the visible luminescence emitted by the UPT particles and finally convert it to the voltage as a signal. V-T and V-c stand for the multiplied voltage units for the test and the control line, respectively, and the ratio V-T/V-C is directly proportional to the number of Y pestis in a sample. We observed a good linearity between the ratio and log CFU/ml of Y pestis above the detection limit, which was approximately 10(4) CFU/mI. The precision of the intra- and inter-assay was below 15% (coefficient of variation, CV). Cross-reactivity with related Gram-negative enteric bacteria was not found. The UPT-LF immunoassay system presented here takes less than 30 min to perform from the sample treatment to the data analysis. The current paper includes only preliminary data concerning the biomedical aspects of the assay, but is more concentrated on the technical details of establishing a rapid manual assay using a state-of-the-art label chemistry. (c) 2006 Elsevier B.V. All rights reserved.

Fluidic flow delay by ionogel passive pumps in microfluidicpaper-based analytical devices

[EN] A new concept for fluid flow manipulation in microfluidic paper-based analytical devices ( µPADs) is presented by introducing ionogel materials as passive pumps. µPADs were fabricated using a new doubleside contact stamping process and ionogels were precisely photopolymerised at the inlet of the µPADs.The ionogels remain mainly on the surface of the paper and get absorbed in the superficial paper-fibers allowing for the liquid to flow from the ionogel into the paper easily. As a proof of concept the fluid flow and mixing behaviour of two different ionogels µPADs were compared with the non-treated µPADs.It was demonstrated that both ionogels highly affect the fluid flow by delaying the flow due to their different physical and chemical properties and water holding capacities.

Eurobalise-Train communication modelling to assess interferences in railway control signalling systems

The evolution of the railway sector depends, to a great extent, on the deployment of advanced railway signalling systems. These signalling systems are based on communication architectures that must cope with complex electromagnetical environments. This paper is outlined in the context of developing the necessary tools to allow the quick deployment of these signalling systems by contributing to an easier analysis of their behaviour under the effect of electromagnetical interferences. Specifically, this paper presents the modelling of the Eurobalise-train communication flow in a general purpose simulation tool. It is critical to guarantee this communication link since any lack of communication may lead to a stop of the train and availability problems. In order to model precisely this communication link we used real measurements done in a laboratory equipped with elements defined in the suitable subsets. Through the simulation study carried out, we obtained performance indicators of the physical layer such as the received power, SNR and BER. The modelling presented in this paper is a required step to be able to provide quality of service indicators related to perturbed scenarios.

Shock / boundary-layer interaction control using three-dimensional bumps for transonic wings

Three-dimensional bumps have been developed and investigated, aiming at the two major objectives of shock-wave / boundary-layer interaction control, i.e. drag reduction and suppression of separation, simultaneously. An experimental investigation has been conducted for a default rounded bump in channel now at University of Cambridge and a computational study has been performed for a spanwise series of rounded bumps mounted on a transonic aerofoil at University of Stuttgart. Observed in both cases are wave drag reduction owing to A-shock structures produced by three-dimensional surface bumps and mild control effects on the boundary layer. The effects of rough surface and tall extension have been investigated as well as several geometric variations and multiple bump configurations. A double configuration of narrow rounded bumps has been found to best perform amongst the tested, considerably reducing wave drag through a well-established A-shock structure with little viscous penalty and thus achieving substantial overall drag reduction. Counter-rotating streamwise vortex pairs have been produced by some configurations as a result of local flow separation, but they have been observed to be confined in relatively narrow wake regions, expected to be beneficial in suppressing large-scale separation under off-design condition despite increase of viscous drag. On the whole a large potential of three-dimensional control with discrete rounded bumps has been demonstrated both experimentally and numerically, and experimental investigation of bumps fitted on a transonic aerofoil or wing is suggested toward practical application.

Experimental investigation of 3D shock / boundary layer interaction control in transonic flows

A novel supersonic wind tunnel setup is proposed to enable the investigation of control on a normal shock wave. Previous experimental arrangements were found to suffer from shock instability. Wind tunnel tests with and without control have confirmed the capability of the new setup to stabilise a shock structure at a target position without changing the nature of the shock wave / boundary layer interaction flow at M∞ = 1.3 and M ∞ = 1.5. Flow visualisation and pressure measurements with the new setup have revealed detailed characteristics of shock wave / boundary layer interactions and a λ-shock structure as well as benefits of control in total drag reduction in the presence of 3D bump control.

The time-resolved flow within an unsteady ejector

The flow field within an unsteady ejector has been investigated using experimental and computational techniques. The experimental results show a peak thrust augmentation of 1.4; numerical simulation gives a value of 1.37. It is shown that the vortex ring dominates the flow field. At optimal thrust augmentation the vortex ring acts like a fluid piston accelerating the fluid inside the ejector. A model is proposed for the operation of unsteady ejectors, based on the vortex ring acting like a fluid piston. Control volume analysis is presented showing that mass entrainment is responsible for thrust augmentation. It is proposed that the spacing of successive vortex rings determines the mass entrainment and therefore thrust augmentation. The efficiency of unsteady ejectors was found to vary between 28% and 32% depending on the L/D ratio of the unsteady jet source. Copyright © 2008 by J H Heffer.

Flow of wet powder in a conical centrifugal filter-an analytical model

A one-dimensional analytical model is developed for the steady state, axisymmetric, slender flow of saturated powder in a rotating perforated cone. Both the powder and the fluid spin with the cone with negligible slip in the hoop direction. They migrate up the wall of the cone along a generator under centrifugal force, which also forces the fluid out of the cone through the powder layer and the porous wall. The flow thus evolves from an over-saturated paste at inlet into a nearly dry powder at outlet. The powder is treated as a Mohr-Coulomb granular solid of constant void fraction and permeability. The shear traction at the wall is assumed to be velocity and pressure dependent. The fluid is treated as Newtonian viscous. The model provides the position of the colour line (the transition from over- to under-saturation) and the flow velocity and thickness profiles over the cone. Surface tension effects are assumed negligible compared to the centrifugal acceleration. Two alternative conditions are considered for the flow structure at inlet: fully settled powder at inlet, and progressive settling of an initially homogeneous slurry. The position of the colour line is found to be similar for these two cases over a wide range of operating conditions. Dominant dimensionless groups are identified which control the position of the colour line in a continuous conical centrifuge. Experimental observations of centrifuges used in the sugar industry provide preliminary validation of the model. © 2011 Elsevier Ltd.

A theoretical approach to the passive control of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this.We study how this spiral mode of vortex breakdown might be suppressed or promoted. We use a Lagrangian approach to identify regions of the flow which are sensitive to small open-loop steady and unsteady (harmonic) forces. We find these regions to be upstream of the vortex breakdown bubble. We investigate passive control using a small axisymmetric control ring. In this case, the steady and unsteady control forces are caused by the drag force on the control ring. We find a narrow region upstream of the bubble where the control ring will stabilise the flow and we verify this using numerical simulations. © 2012 IEEE.

An investigation of the effects of pneumatic actuator design on slip control for heavy vehicles

Progress in reducing actuator delays in pneumatic brake systems is opening the door for advanced anti-lock braking algorithms to be used on heavy goods vehicles. However, little has been published on slip controllers for air-braked heavy vehicles, or the effects of slow pneumatic actuation on their design and performance. This paper introduces a sliding mode slip controller for air-braked heavy vehicles. The effects of pneumatic actuator delays and flow rates on stopping performance and air (energy) consumption are presented through vehicle simulations. Finally, the simulations are validated with experiments using a hardware-in-the-loop rig. It is shown that for each wheel, pneumatic valves with delays smaller than 3ms and orifice diameters around 8mm provide the best performance. © 2013 Copyright Taylor and Francis Group, LLC.

Long range exhaustion-a mathematical model for the axisymmetric air flow of a local exhaust ventilation hood assisted by a turbulent radial jet

Abstract-Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art Banged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air Bow to be reduced and provide a means of identifying and assessing the various parameters that control the air Bow. The mathematical model is formulated in terms of the Stokes steam function, ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained Bow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions. | Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art flanged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air flow to be reduced and provide a means of identifying and assessing the various parameters that control the air flow. The mathematical model is formulated in terms of the Stokes steam function, Ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained flow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Control of aircraft in the terminal manoeuvring area using parallelised sequential Monte Carlo

This paper reports on the use of a parallelised Model Predictive Control, Sequential Monte Carlo algorithm for solving the problem of conflict resolution and aircraft trajectory control in air traffic management specifically around the terminal manoeuvring area of an airport. The target problem is nonlinear, highly constrained, non-convex and uses a single decision-maker with multiple aircraft. The implementation includes a spatio-temporal wind model and rolling window simulations for realistic ongoing scenarios. The method is capable of handling arriving and departing aircraft simultaneously including some with very low fuel remaining. A novel flow field is proposed to smooth the approach trajectories for arriving aircraft and all trajectories are planned in three dimensions. Massive parallelisation of the algorithm allows solution speeds to approach those required for real-time use.

Investigation of a radial-inflow bleed as a potential for compressor clearance control

The mismatch in thermal response between a High Pressure Compressor (HPC) drum and casing is a limiting factor in the reduction of compressor clearance. An experimental test rig has been used to demonstrate the concept of radial inflow to reduce the thermal time constant of HPC discs. The testing uses a simulated idle - Maximum Take Off (MTO) - idle transient in order to measure the thermal response directly. The testing is fully scaled in the dimensionless sense to engine conditions. A simple closure model based on lumped capacitance is used to illustrate the scope of potential benefits. The proof-of-concept testing shows that HPC disc time constant reductions of the order 2 are feasible with a radial-inflow bleed of only 4% of bore flow at scaled MTO conditions. Using the experimental results, the simple closure modelling suggests that for a stage with a significant mismatch in thermal response, reductions in 2D axis-symmetric clearance of as much as 50% at MTO conditions may be possible along with significant scope for improvements at cruise conditions. Copyright © 2013 by ASME.