Structure-Preserving Model Reduction of Physical Network Systems

This paper considers physical network systems where the energy storage is naturally associated to the nodes of the graph, while the edges of the graph correspond to static couplings. The first sections deal with the linear case, covering examples such as mass-damper and hydraulic systems, which have a structure that is similar to symmetric consensus dynamics. The last section is concerned with a specific class of nonlinear physical network systems; namely detailed-balanced chemical reaction networks governed by mass action kinetics. In both cases, linear and nonlinear, the structure of the dynamics is similar, and is based on a weighted Laplacian matrix, together with an energy function capturing the energy storage at the nodes. We discuss two methods for structure-preserving model reduction. The first one is clustering; aggregating the nodes of the underlying graph to obtain a reduced graph. The second approach is based on neglecting the energy storage at some of the nodes, and subsequently eliminating those nodes (called Kron reduction).</p

Semi-blind robust indentification and robust control approach to personalized anemia management.

The homeostatic blood hemoglobin (Hb) content of a healthy individual varies between the range of 14-18 g/dL for a male and 12-16 g/dL for a female. This quantity provides an estimate of red blood cell (RBC) count in circulation at any given moment. RBC is a protein carrying substance that transports oxygen from the lungs to other tissues in the body and is synthesized by the kidney through a process known as erythropoiesis where erythropoietin is secreted in response to hypoxia. In this regard, the kidneys act not only as a controller but also as a sensor in regulating RBC levels. Patients with chronic kidney diseases (CKD) have dysfunctional kidneys that compromise these fundamental kidney functions. Consequently, anemia is developed. Anemics of CKD have low levels of Hb that must be controlled and properly regulated to the appropriate therapeutic range. Until the discovery of recombinant human erythropoietin (EPO) over three decades ago, treatment procedure of anemia conditions primarily involved repeated blood transfusions–a process known to be associated with several other health related complications. This discovery resulted in a paradigm shift in anemia management from blood transfusions to dosage therapies. The main objective of anemia management with EPO is to increase patients’ hemoglobin level from low to a suitable therapeutic range as defined by the National Kidney Foundation-Kidney Disease Outcomes Quality Initiative (NKF-KDOI) to be in the range of 10 - 12 g/dL while avoiding response values beyond 14 g/dL to prevent other complications associated with EPO medication. It is therefore imperative that clinicians balance dosage efficacy and toxicity in anemia management therapies. At most treatment facilities, protocols are developed to conform to NKF-KDOI recommendations. These protocols are generally based on EPO packet inserts and the expected Hb responses from the average patient. The inevitable variability within the patient group makes this “one-size-fits-all” dosing scheme non-optimal, at best, and potentially dangerous for certain group of patients that do not adhere to the notion of expected “average” response. A dosing strategy that is tailored to the individual patients’ response to EPO medication could provide a better alternative to the current treatment methods. An objective of this work is to develop EPO dosing strategies tailored to the individual patients using robust identification techniques and modern feedback control methods. First, a unique model is developed based on Hb responses and dosage EPO of the individual patients using semi-blind robust identification techniques. This provides a nominal model and a quantitative information on model uncertainty that accounts for other possible patient’s dynamics not considered in the modeling process. This is in the framework of generalized interpolation theory. Then, from the derived nominal model and the associated uncertainty information, robust controller is designed via the =H1-synthesis methods to provide a new dosing strategies for the individual patients. The H1 control theory has a feature of minimizing the influence of some unknown worst case gain disturbance on a system. Finally, a framework is provided to strategize dosing protocols for newly admitted patients

A Bayesian approach to robust identification: application to fault detection

In the Control Engineering field, the so-called Robust Identification techniques deal with the problem of obtaining not only a nominal model of the plant, but also an estimate of the uncertainty associated to the nominal model. Such model of uncertainty is typically characterized as a region in the parameter space or as an uncertainty band around the frequency response of the nominal model. Uncertainty models have been widely used in the design of robust controllers and, recently, their use in model-based fault detection procedures is increasing. In this later case, consistency between new measurements and the uncertainty region is checked. When an inconsistency is found, the existence of a fault is decided. There exist two main approaches to the modeling of model uncertainty: the deterministic/worst case methods and the stochastic/probabilistic methods. At present, there are a number of different methods, e.g., model error modeling, set-membership identification and non-stationary stochastic embedding. In this dissertation we summarize the main procedures and illustrate their results by means of several examples of the literature. As contribution we propose a Bayesian methodology to solve the robust identification problem. The approach is highly unifying since many robust identification techniques can be interpreted as particular cases of the Bayesian framework. Also, the methodology can deal with non-linear structures such as the ones derived from the use of observers. The obtained Bayesian uncertainty models are used to detect faults in a quadruple-tank process and in a three-bladed wind turbine

Isogeometric Analysis for Electromagnetism

The combination of numerical analysis with the scanning technology has been seeing increased use in many research areas. There is an emerging need for high-fidelity geometric modeling and meshing for practical applications. The Isogeometric Analysis (IGA) is a comprehensive computational framework, which integrates geometric modeling and meshing with analysis. Different from other existing numerical methods, the IGA can generate analysis ready models without loss of geometrical accuracy. In IGA, the continuity and the quality of a solution can be conveniently controlled and refined. These features enable IGA to integrate modeling, analysis, and design in a unified framework, the root idea of IGA. The IGA for electromagmetics is studied here for steady and transient electromagnetics as well as electromagnetic scattering. The solution procedure and the associated Matlab codes are developed to simulate the electromagnetic radiation on a biological tissues. The scattered and the total electrical fields are computed over the complex geometry of a brain section with realistic material properties. A perfectly matched layer (PML) is developed to model the far field boundary condition. The IGA platform developed here offers a reliable simulation within an accurate representation of the geometry. The results of this research can be used both in evaluating the potential health and safety risks of electromagnetic radiations and in optimizing the design of radiating devices used in non-invasive diagnostics and therapies

HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems

System- and Data-Driven Methods and Algorithms

An increasing complexity of models used to predict real-world systems leads to the need for algorithms to replace complex models with far simpler ones, while preserving the accuracy of the predictions. This two-volume handbook covers methods as well as applications. This first volume focuses on real-time control theory, data assimilation, real-time visualization, high-dimensional state spaces and interaction of different reduction techniques