Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control chemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

Voluntary Control of Multisaccade Gaze Shifts During Movement Preparation and Execution

Ramakrishnan A, Chokhandre S, Murthy A. Voluntary control of multisaccade gaze shifts during movement preparation and execution. J Neurophysiol 103: 2400-2416, 2010. First published February 17, 2010; doi: 10.1152/jn.00843.2009. Although the nature of gaze control regulating single saccades is relatively well documented, how such control is implemented to regulate multisaccade gaze shifts is not known. We used highly eccentric targets to elicit multisaccade gaze shifts and tested the ability of subjects to control the saccade sequence by presenting a second target on random trials. Their response allowed us to test the nature of control at many levels: before, during, and between saccades. Although the saccade sequence could be inhibited before it began, we observed clear signs of truncation of the first saccade, which confirmed that it could be inhibited in midflight as well. Using a race model that explains the control of single saccades, we estimated that it took about 100 ms to inhibit a planned saccade but took about 150 ms to inhibit a saccade during its execution. Although the time taken to inhibit was different, the high subject-wise correlation suggests a unitary inhibitory control acting at different levels in the oculomotor system. We also frequently observed responses that consisted of hypometric initial saccades, followed by secondary saccades to the initial target. Given the estimates of the inhibitory process provided by the model that also took into account the variances of the processes as well, the secondary saccades (average latency similar to 215 ms) should have been inhibited. Failure to inhibit the secondary saccade suggests that the intersaccadic interval in a multisaccade response is a ballistic stage. Collectively, these data indicate that the oculomotor system can control a response until a very late stage in its execution. However, if the response consists of multiple movements then the preparation of the second movement becomes refractory to new visual input, either because it is part of a preprogrammed sequence or as a consequence of being a corrective response to a motor error.

Thyroidal control of hepatic release and metabolism of vitamin A

The inverse relationship that exists between thyroxine and the vitamin A level of plasma has been examined in chicken. Thyroxine treatment leads to a decrease in the level of vitamin A carrier proteins, retinol-binding protein and prealbumin-2 in plasma and liver. There is an accumulation of vitamin A in the liver, with a greater proportion of vitamin A alcohol being present compared to that of control birds. In thyroxine treatment there is enhanced plasma turnover of retinol-binding protein and prealbumin-2, while their rates of synthesis are marginally increased. Amino acid supplementation partially counteracts effects of thyroxine treatment. Amino acid supplementation of thyroxine-treated birds does not alter the plasma turnover rates of retinol-binding protein and prealbumin-2 but increases substentially their rates of synthesis. The release of vitamin A into circulation is interfered with in hyperthyroidism due to inadequate availability of retinol-binding protein being caused by enhanced plasma turnover rate not compensated for by synthesis.

Sequential design of decentralized control systems

Our main result is a new sequential method for the design of decentralized control systems. Controller synthesis is conducted on a loop-by-loop basis, and at each step the designer obtains an explicit characterization of the class C of all compensators for the loop being closed that results in closed-loop system poles being in a specified closed region D of the s-plane, instead of merely stabilizing the closed-loop system. Since one of the primary goals of control system design is to satisfy basic performance requirements that are often directly related to closed-loop pole location (bandwidth, percentage overshoot, rise time, settling time), this approach immediately allows the designer to focus on other concerns such as robustness and sensitivity. By considering only compensators from class C and seeking the optimum member of that set with respect to sensitivity or robustness, the designer has a clearly-defined limited optimization problem to solve without concern for loss of performance. A solution to the decentralized tracking problem is also provided. This design approach has the attractive features of expandability, the use of only 'local models' for controller synthesis, and fault tolerance with respect to certain types of failure.

Control of switching surges due to energization of 765kV UHV transmission line

Reduction of switching surge over voltages allows an economic design of UHV transmission system with reduced insulation. The various means of switching surge over voltage control with pre-insertion resistors/closing resistors, shunt re-actors and controlled switching are illustrated. The switching surge over voltages during the energization of series compensated line are compared with uncompensated line. An Electromagnetic transients program has been developed for studying the effect of various means of control of switching transients during 765kV UHV transmission line energization. This paper presents the studies carried out on switching surges control in 765kV UHV transmission line energization.

In-cylinder investigations and analysis of a SI gas engine fuelled with H-2 and CO rich syngas fuel: Sensitivity analysis of combustion descriptors for engine diagnostics and control

The sensitivity of combustion phasing and combustion descriptors to ignition timing, load and mixture quality on fuelling a multi-cylinder natural gas engine with bio-derived H-2 and CO rich syngas is addressed. While the descriptors for conventional fuels are well established and are in use for closed loop engine control, presence of H-2 in syngas potentially alters the mixture properties and hence combustion phasing, necessitating the current study. The ability of the descriptors to predict abnormal combustion, hitherto missing in the literature, is also addressed. Results from experiments using multi-cylinder engines and numerical studies using zero dimensional Wiebe function based simulation models are reported. For syngas with 20% H-2 and CO and 2% CH4 (producer gas), an ignition retard of 5 +/- 1 degrees was required compared to natural gas ignition timing to achieve peak load of 72.8 kWe. It is found that, for syngas, whose flammability limits are 0.42-1.93, the optimal engine operation was at an equivalence ratio of 1.12. The same methodology is extended to a two cylinder engine towards addressing the influence of syngas composition, especially H-2 fraction (varying from 13% to 37%), on the combustion phasing. The study confirms the utility of pressure trace derived combustion descriptors, except for the pressure trace first derivative, in describing the MBT operating condition of the engine when fuelled with an alternative fuel. Both experiments and analysis suggest most of the combustion descriptors to be independent of the engine load and mixture quality. A near linear relationship with ignition angle is observed. The general trend(s) of the combustion descriptors for syngas fuelled operation are similar to those of conventional fuels; the differences in sensitivity of the descriptors for syngas fuelled engine operation requires re-calibration of control logic for MBT conditions. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Optimal bounded control of first-passage failure of quasi-integrable Hamiltonian systems with wide-band random excitation

The optimal bounded control of quasi-integrable Hamiltonian systems with wide-band random excitation for minimizing their first-passage failure is investigated. First, a stochastic averaging method for multi-degrees-of-freedom (MDOF) strongly nonlinear quasi-integrable Hamiltonian systems with wide-band stationary random excitations using generalized harmonic functions is proposed. Then, the dynamical programming equations and their associated boundary and final time conditions for the control problems of maximizinig reliability and maximizing mean first-passage time are formulated based on the averaged It$\ddot{\rm o}$ equations by applying the dynamical programming principle. The optimal control law is derived from the dynamical programming equations and control constraints. The relationship between the dynamical programming equations and the backward Kolmogorov equation for the conditional reliability function and the Pontryagin equation for the conditional mean first-passage time of optimally controlled system is discussed. Finally, the conditional reliability function, the conditional probability density and mean of first-passage time of an optimally controlled system are obtained by solving the backward Kolmogorov equation and Pontryagin equation. The application of the proposed procedure and effectiveness of control strategy are illustrated with an example.

Control dynamics of severe acute respiratory syndrome transmission

Severe acute respiratory syndrome (SARS) is a serious disease with many puzzling features. We present a simple, dynamic model to assess the epidemic potential of SARS and the effectiveness of control measures. With this model, we analysed the SARS epidemic data in Beijing. The data fitting gives the basic case reproduction number of 2.16 leading to the outbreak, and the variation of the effective reproduction number reflecting the control effect. Noticeably, our study shows that the response time and the strength of control measures have significant effects on the scale of the outbreak and the lasting time of the epidemic.

Control of dynamical systems with temporal logic specifications

This thesis is motivated by safety-critical applications involving autonomous air, ground, and space vehicles carrying out complex tasks in uncertain and adversarial environments. We use temporal logic as a language to formally specify complex tasks and system properties. Temporal logic specifications generalize the classical notions of stability and reachability that are studied in the control and hybrid systems communities. Given a system model and a formal task specification, the goal is to automatically synthesize a control policy for the system that ensures that the system satisfies the specification. This thesis presents novel control policy synthesis algorithms for optimal and robust control of dynamical systems with temporal logic specifications. Furthermore, it introduces algorithms that are efficient and extend to high-dimensional dynamical systems.

The first contribution of this thesis is the generalization of a classical linear temporal logic (LTL) control synthesis approach to optimal and robust control. We show how we can extend automata-based synthesis techniques for discrete abstractions of dynamical systems to create optimal and robust controllers that are guaranteed to satisfy an LTL specification. Such optimal and robust controllers can be computed at little extra computational cost compared to computing a feasible controller.

The second contribution of this thesis addresses the scalability of control synthesis with LTL specifications. A major limitation of the standard automaton-based approach for control with LTL specifications is that the automaton might be doubly-exponential in the size of the LTL specification. We introduce a fragment of LTL for which one can compute feasible control policies in time polynomial in the size of the system and specification. Additionally, we show how to compute optimal control policies for a variety of cost functions, and identify interesting cases when this can be done in polynomial time. These techniques are particularly relevant for online control, as one can guarantee that a feasible solution can be found quickly, and then iteratively improve on the quality as time permits.

The final contribution of this thesis is a set of algorithms for computing feasible trajectories for high-dimensional, nonlinear systems with LTL specifications. These algorithms avoid a potentially computationally-expensive process of computing a discrete abstraction, and instead compute directly on the system's continuous state space. The first method uses an automaton representing the specification to directly encode a series of constrained-reachability subproblems, which can be solved in a modular fashion by using standard techniques. The second method encodes an LTL formula as mixed-integer linear programming constraints on the dynamical system. We demonstrate these approaches with numerical experiments on temporal logic motion planning problems with high-dimensional (10+ states) continuous systems.

Temporal and spacial control of RNA stability in the early embryo of Drosophila melanogaster

Early embryogenesis in metazoa is controlled by maternally synthesized products. Among these products, the mature egg is loaded with transcripts representing approximately two thirds of the genome. A subset of this maternal RNA pool is degraded prior to the transition to zygotic control of development. This transfer of control of development from maternal to zygotic products is referred to as the midblastula transition (or MBT). It is believed that the degradation of maternal transcripts is required to terminate maternal control of development and to allow zygotic control of development to begin. Until now this process of maternal transcript degradation and the subsequent timing of the MBT has been poorly understood. I have demonstrated that in the early embryo there are two independent RNA degradation pathways, either of which is sufficient for transcript elimination. However, only the concerted action of both pathways leads to elimination of transcripts with the correct timing, at the MBT. The first pathway is maternally encoded, is triggered by egg activation, and is targeted to specific classes of mRNAs through cis-acting elements in the 3' untranslated region (UTR}. The second pathway is activated 2 hr after fertilization and functions together with the maternal pathway to ensure that transcripts are degraded by the MBT. In addition, some transcripts fail to degrade at select subcellular locations adding an element of spatial control to RNA degradation. The spatial control of RNA degradation is achieved by protecting, or masking, transcripts from the degradation machinery. The RNA degradation and protection events are regulated by distinct cis-elements in the 3' untranslated region (UTR). These results provide the first systematic dissection of this highly conserved process in development and demonstrate that RNA degradation is a novel mechanism used for both temporal and spatial control of development.

Distributed Optimal Control of Cyber-Physical Systems: Controller Synthesis, Architecture Design and System Identification

The centralized paradigm of a single controller and a single plant upon which modern control theory is built is no longer applicable to modern cyber-physical systems of interest, such as the power-grid, software defined networks or automated highways systems, as these are all large-scale and spatially distributed. Both the scale and the distributed nature of these systems has motivated the decentralization of control schemes into local sub-controllers that measure, exchange and act on locally available subsets of the globally available system information. This decentralization of control logic leads to different decision makers acting on asymmetric information sets, introduces the need for coordination between them, and perhaps not surprisingly makes the resulting optimal control problem much harder to solve. In fact, shortly after such questions were posed, it was realized that seemingly simple decentralized optimal control problems are computationally intractable to solve, with the Wistenhausen counterexample being a famous instance of this phenomenon. Spurred on by this perhaps discouraging result, a concerted 40 year effort to identify tractable classes of distributed optimal control problems culminated in the notion of quadratic invariance, which loosely states that if sub-controllers can exchange information with each other at least as quickly as the effect of their control actions propagates through the plant, then the resulting distributed optimal control problem admits a convex formulation.

The identification of quadratic invariance as an appropriate means of "convexifying" distributed optimal control problems led to a renewed enthusiasm in the controller synthesis community, resulting in a rich set of results over the past decade. The contributions of this thesis can be seen as being a part of this broader family of results, with a particular focus on closing the gap between theory and practice by relaxing or removing assumptions made in the traditional distributed optimal control framework. Our contributions are to the foundational theory of distributed optimal control, and fall under three broad categories, namely controller synthesis, architecture design and system identification.

We begin by providing two novel controller synthesis algorithms. The first is a solution to the distributed H-infinity optimal control problem subject to delay constraints, and provides the only known exact characterization of delay-constrained distributed controllers satisfying an H-infinity norm bound. The second is an explicit dynamic programming solution to a two player LQR state-feedback problem with varying delays. Accommodating varying delays represents an important first step in combining distributed optimal control theory with the area of Networked Control Systems that considers lossy channels in the feedback loop. Our next set of results are concerned with controller architecture design. When designing controllers for large-scale systems, the architectural aspects of the controller such as the placement of actuators, sensors, and the communication links between them can no longer be taken as given -- indeed the task of designing this architecture is now as important as the design of the control laws themselves. To address this task, we formulate the Regularization for Design (RFD) framework, which is a unifying computationally tractable approach, based on the model matching framework and atomic norm regularization, for the simultaneous co-design of a structured optimal controller and the architecture needed to implement it. Our final result is a contribution to distributed system identification. Traditional system identification techniques such as subspace identification are not computationally scalable, and destroy rather than leverage any a priori information about the system's interconnection structure. We argue that in the context of system identification, an essential building block of any scalable algorithm is the ability to estimate local dynamics within a large interconnected system. To that end we propose a promising heuristic for identifying the dynamics of a subsystem that is still connected to a large system. We exploit the fact that the transfer function of the local dynamics is low-order, but full-rank, while the transfer function of the global dynamics is high-order, but low-rank, to formulate this separation task as a nuclear norm minimization problem. Finally, we conclude with a brief discussion of future research directions, with a particular emphasis on how to incorporate the results of this thesis, and those of optimal control theory in general, into a broader theory of dynamics, control and optimization in layered architectures.

Management of Biesboch reservoirs for quality control with special reference to eutrophication

The three Biesbosch Reservoirs are pumped storage reservoirs, fed with rather polluted and highly eutrophic water from the River Meuse. Air injection at the bottom of the reservoirs prevents thermal stratification, which would otherwise result in serious water quality deterioration. Reservoir mixing also serves as an economic algal control measure; mixing over sufficient depth causes light to play the role of limiting factor and this, combined with zooplankton grazing, keeps the biomass of phytoplankton at acceptable levels. Special problems are caused by benthic, geosmin-producing Oscillatoria species growing on the inner embankment. Rooting up the bottom with a harrow is used as the method of control, based on underwater observations by biological staff trained as SCUBA-divers. With regard to pollutant behaviour the three reservoirs act as a series of fully mixed reactors. This enables the application of kinetic models to describe their behaviour and allows the use of a selective intake policy, e.g. for suspended solids with associated contaminants, ammonia and polynuclear aromatic hydrocarbons. A combination of selective intake and self- purification processes - enhanced by the compartmentalisation of the storage volume in three reservoirs - leads to a striking improvement for many water-quality parameters.

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.

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.

Feedback control and the arrow of time

The purpose of this paper is to highlight the central role that the time asymmetry of stability plays in feedback control. We show that this provides a new perspective on the use of doubly-infinite or semi-infinite time axes for signal spaces in control theory. We then focus on the implication of this time asymmetry in modeling uncertainty, regulation and robust control. We point out that modeling uncertainty and the ease of control depend critically on the direction of time. We finally discuss the relationship of this control-based time arrow with the well-known arrows of time in physics. © 2008 IEEE.