Design, specification, and synthesis of aircraft electric power systems control logic

Cyber-physical systems integrate computation, networking, and physical processes. Substantial research challenges exist in the design and verification of such large-scale, distributed sensing, ac- tuation, and control systems. Rapidly improving technology and recent advances in control theory, networked systems, and computer science give us the opportunity to drastically improve our approach to integrated flow of information and cooperative behavior. Current systems rely on text-based spec- ifications and manual design. Using new technology advances, we can create easier, more efficient, and cheaper ways of developing these control systems. This thesis will focus on design considera- tions for system topologies, ways to formally and automatically specify requirements, and methods to synthesize reactive control protocols, all within the context of an aircraft electric power system as a representative application area.

This thesis consists of three complementary parts: synthesis, specification, and design. The first section focuses on the synthesis of central and distributed reactive controllers for an aircraft elec- tric power system. This approach incorporates methodologies from computer science and control. The resulting controllers are correct by construction with respect to system requirements, which are formulated using the specification language of linear temporal logic (LTL). The second section addresses how to formally specify requirements and introduces a domain-specific language for electric power systems. A software tool automatically converts high-level requirements into LTL and synthesizes a controller.

The final sections focus on design space exploration. A design methodology is proposed that uses mixed-integer linear programming to obtain candidate topologies, which are then used to synthesize controllers. The discrete-time control logic is then verified in real-time by two methods: hardware and simulation. Finally, the problem of partial observability and dynamic state estimation is ex- plored. Given a set placement of sensors on an electric power system, measurements from these sensors can be used in conjunction with control logic to infer the state of the system.

Synthesis and potency of long end-associative polymers for mist control

Long linear polymers that are end-functionalized with associative groups were studied as additives to hydrocarbon fluids to mitigate the fire hazard associated with the presence of mist in a crash scenario. These polymers were molecularly designed to overcome both the shear-degradation of long polymer chains in turbulent flows, and the chain collapse induced by the random placement of associative groups along polymer backbones. Architectures of associative groups on the polymer chain ends that were tested included clusters of self-associative carboxyl groups and pairs of hetero-complementary associative units.

Linear polymers with clusters of discrete numbers of carboxyl groups on their chain ends were investigated first: an innovative synthetic strategy was devised to achieve unprecedented backbone lengths and precise control of the number of carboxyl groups on chain ends (N). We found that a very narrow range of N allows the co-existence of sufficient end-association strength and polymer solubility in apolar media. Subsequent steady-flow rheological study on solution behavior of such soluble polymers in apolar media revealed that the end-association of very long chains in apolar media leads to the formation of flower-like micelles interconnected by bridging chains, which trap significant fraction of polymer chains into looped structures with low contribution to mist-control. The efficacy of very long 1,4-polybutadiene chains end-functionalized with clusters of four carboxyl groups as mist-control additives for jet fuel was further tested. In addition to being shear-resistant, the polymer was found capable of providing fire-protection to jet fuel at concentrations as low as 0.3wt%. We also found that this polymer has excellent solubility in jet fuel over a wide range of temperature (-30 to +70°C) and negligible interference with dewatering of jet fuel. It does not cause an adverse increase in viscosity at concentrations where mist-control efficacy exists.

Four pairs of hetero-complementary associative end-groups of varying strengths were subsequently investigated, in the hopes of achieving supramolecular aggregates with both mist-control ability and better utilization of polymer building blocks. Rheological study of solutions of the corresponding complementary associative polymer pairs in apolar media revealed the strength of complementary end-association required to achieve supramolecular aggregates capable of modulating rheological properties of the solution.

Both self-associating and complementary associating polymers have therefore been found to resist shear degradation. The successful strategy of building soluble, end-associative polymers with either self-associative or complementary associative groups will guide the next generation of mist-control technology.

Boundary-layer transition on a slender cone in hypervelocity flow with real gas effects

The laminar to turbulent transition process in boundary layer flows in thermochemical nonequilibrium at high enthalpy is measured and characterized. Experiments are performed in the T5 Hypervelocity Reflected Shock Tunnel at Caltech, using a 1 m length 5-degree half angle axisymmetric cone instrumented with 80 fast-response annular thermocouples, complemented by boundary layer stability computations using the STABL software suite. A new mixing tank is added to the shock tube fill apparatus for premixed freestream gas experiments, and a new cleaning procedure results in more consistent transition measurements. Transition location is nondimensionalized using a scaling with the boundary layer thickness, which is correlated with the acoustic properties of the boundary layer, and compared with parabolized stability equation (PSE) analysis. In these nondimensionalized terms, transition delay with increasing CO₂ concentration is observed: tests in 100% and 50% CO₂, by mass, transition up to 25% and 15% later, respectively, than air experiments. These results are consistent with previous work indicating that CO₂ molecules at elevated temperatures absorb acoustic instabilities in the MHz range, which is the expected frequency of the Mack second-mode instability at these conditions, and also consistent with predictions from PSE analysis. A strong unit Reynolds number effect is observed, which is believed to arise from tunnel noise. N_Tr for air from 5.4 to 13.2 is computed, substantially higher than previously reported for noisy facilities. Time- and spatially-resolved heat transfer traces are used to track the propagation of turbulent spots, and convection rates at 90%, 76%, and 63% of the boundary layer edge velocity, respectively, are observed for the leading edge, centroid, and trailing edge of the spots. A model constructed with these spot propagation parameters is used to infer spot generation rates from measured transition onset to completion distance. Finally, a novel method to control transition location with boundary layer gas injection is investigated. An appropriate porous-metal injector section for the cone is designed and fabricated, and the efficacy of injected CO₂ for delaying transition is gauged at various mass flow rates, and compared with both no injection and chemically inert argon injection cases. While CO₂ injection seems to delay transition, and argon injection seems to promote it, the experimental results are inconclusive and matching computations do not predict a reduction in N factor from any CO₂ injection condition computed.

Distributed Load Control in Multiphase Radial Networks

The current power grid is on the cusp of modernization due to the emergence of distributed generation and controllable loads, as well as renewable energy. On one hand, distributed and renewable generation is volatile and difficult to dispatch. On the other hand, controllable loads provide significant potential for compensating for the uncertainties. In a future grid where there are thousands or millions of controllable loads and a large portion of the generation comes from volatile sources like wind and solar, distributed control that shifts or reduces the power consumption of electric loads in a reliable and economic way would be highly valuable.

Load control needs to be conducted with network awareness. Otherwise, voltage violations and overloading of circuit devices are likely. To model these effects, network power flows and voltages have to be considered explicitly. However, the physical laws that determine power flows and voltages are nonlinear. Furthermore, while distributed generation and controllable loads are mostly located in distribution networks that are multiphase and radial, most of the power flow studies focus on single-phase networks.

This thesis focuses on distributed load control in multiphase radial distribution networks. In particular, we first study distributed load control without considering network constraints, and then consider network-aware distributed load control.

Distributed implementation of load control is the main challenge if network constraints can be ignored. In this case, we first ignore the uncertainties in renewable generation and load arrivals, and propose a distributed load control algorithm, Algorithm 1, that optimally schedules the deferrable loads to shape the net electricity demand. Deferrable loads refer to loads whose total energy consumption is fixed, but energy usage can be shifted over time in response to network conditions. Algorithm 1 is a distributed gradient decent algorithm, and empirically converges to optimal deferrable load schedules within 15 iterations.

We then extend Algorithm 1 to a real-time setup where deferrable loads arrive over time, and only imprecise predictions about future renewable generation and load are available at the time of decision making. The real-time algorithm Algorithm 2 is based on model-predictive control: Algorithm 2 uses updated predictions on renewable generation as the true values, and computes a pseudo load to simulate future deferrable load. The pseudo load consumes 0 power at the current time step, and its total energy consumption equals the expectation of future deferrable load total energy request.

Network constraints, e.g., transformer loading constraints and voltage regulation constraints, bring significant challenge to the load control problem since power flows and voltages are governed by nonlinear physical laws. Remarkably, distribution networks are usually multiphase and radial. Two approaches are explored to overcome this challenge: one based on convex relaxation and the other that seeks a locally optimal load schedule.

To explore the convex relaxation approach, a novel but equivalent power flow model, the branch flow model, is developed, and a semidefinite programming relaxation, called BFM-SDP, is obtained using the branch flow model. BFM-SDP is mathematically equivalent to a standard convex relaxation proposed in the literature, but numerically is much more stable. Empirical studies show that BFM-SDP is numerically exact for the IEEE 13-, 34-, 37-, 123-bus networks and a real-world 2065-bus network, while the standard convex relaxation is numerically exact for only two of these networks.

Theoretical guarantees on the exactness of convex relaxations are provided for two types of networks: single-phase radial alternative-current (AC) networks, and single-phase mesh direct-current (DC) networks. In particular, for single-phase radial AC networks, we prove that a second-order cone program (SOCP) relaxation is exact if voltage upper bounds are not binding; we also modify the optimal load control problem so that its SOCP relaxation is always exact. For single-phase mesh DC networks, we prove that an SOCP relaxation is exact if 1) voltage upper bounds are not binding, or 2) voltage upper bounds are uniform and power injection lower bounds are strictly negative; we also modify the optimal load control problem so that its SOCP relaxation is always exact.

To seek a locally optimal load schedule, a distributed gradient-decent algorithm, Algorithm 9, is proposed. The suboptimality gap of the algorithm is rigorously characterized and close to 0 for practical networks. Furthermore, unlike the convex relaxation approach, Algorithm 9 ensures a feasible solution. The gradients used in Algorithm 9 are estimated based on a linear approximation of the power flow, which is derived with the following assumptions: 1) line losses are negligible; and 2) voltages are reasonably balanced. Both assumptions are satisfied in practical distribution networks. Empirical results show that Algorithm 9 obtains 70+ times speed up over the convex relaxation approach, at the cost of a suboptimality within numerical precision.

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.

Micro-ramp control for oblique shock wave / boundary layer interactions

Supersonic engine intakes operating supercritically feature shock wave / boundary layer interactions (SBLIs), which are conventionally controlled using boundary layer bleed. The momentum loss of bleed flow causes high drag, compromising intake performance. Micro-ramp sub-boundary layer vortex generators (SBVGs) have been proposed as an alternative form of flow control for oblique SBLIs in order to reduce the bleed requirement. Experiments have been conducted at Mach 2.5 to characterise the flow details on such devices and investigate their ability to control the interaction between an oblique shock wave and the naturally grown turbulent boundary layer on the tunnel floor. Micro-ramps of four sizes with heights ranging from 25% to 75% of the uncontrolled boundary layer thickness were tested. The flow over all sizes of microramp was found to be similar, featuring streamwise counter-rotating vortices which entrain high momentum fluid, locally reducing the boundary layer displacement thickness. When installed ahead of the shock interaction it was found that the positioning of the micro-ramps is of limited importance. Micro-ramps did not eliminate flow separation. However, the previously two-dimensional separation was broken up into periodic three-dimensional separation zones. The interaction length was reduced and the pressure gradient across the interaction was increased.

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.

Studies on the chemical control of psychrophilic bacterial spoilage of fish. 4 - The effect of chemical preservatives on the growth of psychrophilic bacteria isolated from marine fish

Dehydroacetic acid and ammonia were found to be very effective in checking the growth of all the cultures at all concentrations tried. The two nitrofuran derivatives namely, semicarbazone and AF-2 were fairly effective, semicarbazone being more effective than AF-2. Sodium nitrite was found to be totally ineffective against the cultures at all concentrations tried.