DOCSLIB.ORG

Control of Error and Convergence in ODE Solvers

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Control of Error and Convergence in ODE Solvers Control of Error and Convergence in ODE Solvers Kjell Gustafsson Integration method + Differential equation Control of Error and Convergence in ODE Solvers Kjell Gustafsson Tekn. Lic, Civ. Ing. Kristianstads Nation Akademisk avhandling som för avläggande av teknisk doktors- examen vid Tekniska fakulteten vid Universitetet i Lund kommer att offentligen försvaras i sal M:A, Maskinhuset, Lunds Tekniska Högskola, fredagen den 22 maj 1992, kl 10.15. A PhD thesis to be defended in public in room M:A, the M-building, Lund Institute of Technology, on Friday May 221992, at 10.15. Jr Control of Error and Convergence in ODE Solvers Kjell Gustafsson Integration method + Differential equation Lund 1992 Cover figure: A block diagram demonstrating the feedback control viewpoint of an implicit in- tegration method, cf. Chapter 5. Department of Automatic Control Lund Institute of Technology Box 118 S-221 00 LUND Sweden © Kjell Gustafsson, 1992 Printed in Sweden by Studentlitteratur's printing-office Lund 1992 . Document nmme Department of Automatic Control DOCTORAL DISSERTATION Lund Institute of Technology : Date of iuuc P.O. Box 118 i March 1992 S-221 00 Lund Sweden , Document Number ISRN LUTFD2/TFRT--1036--SE Author(s) Supervisor Kjell Gustafsson Karl Johan Åström, Gustaf Söderlind, Per Hagander Sponsoring organisation Swedish Research Council for Engineering Science» (TFR), contract 222/91-405. Nordiska Forskarut- bildningsakademien, contract 31065.34.024/91). Title and subtitle Control of Error and Convergence in ODE Solvers Abstract Feedback is a general principle that can be used in ma:iy different contexts. In this thesis it is applied to numerical integration of ordinary differential equations. Ar. ..ivanced integration method includes parameters and variables that should be adjusted during the execu'or In addition, the integration method should be able to automatically handle situations such as: initia!:z -. ;on, restart after failures, etc. In this thesis we regard the algorithms for parameter adjustment and s.. division as a controller. The controller measures different variables that tell the current status of the inte;- ..tion, and based on this information it decides how to continue. The design of the controller is vital in or •: lo accurately and efficiently solve i- 'arge clas» of ordinary differential equations. The application of feedback control may appear farf tched, but numerical integration methods are in fact dynamical systems. This is often overlooked in traditional numerical analysis. We derive ciynamic models that describe the behavior of the integration metho as well as the standard control algori'.ims in use today. Using these models it is possible to analyze proper';.es of current algorithms, and also expl iin some generally observed misbehaviors. Further, we use the acquireri insight to derive new and improved control algorithms, both for explicit and implicit Runge-Kutta methods. In the explicit case, the new controller gives good overall performance. In particular it • vercomes the problem with an oscillating ctepsize sequence that is often experienced when the stepsize is r-stricted by numerical stability. The controller for implicit methods is designed so that it tracks changes in Ke differential equation better than current algorithms. In addition, it includes a new strategy for the equaion solver, which allows the stepsize to vary more freely. This leads to smoother error control without excessive operations on the iteration matrix. Key words numerical integration, ordinary differential equations, it :al value problems, simulation, Runge-Kutta meth- ods, stepsize selection, error control, convergence cont.-i•', numerical analysis Classification system and/or index terms (if any) Supplementary bibliographical information ISSN and key title •l-N I 0280-5316 j Language Number of pages Heapient '* notes English 184 Security classification The report may be ordered from the Department of Automatic Control or borrowed through f S.- JJniver§ity Library 2, Box 1010, S-221 03 Lund, Sweden. Telex 33248 lubbia fond t Contents Preface 7 Acknowledgements 7 1. Introduction 9 1.1 Is This Thesis Really Needed? 10 1.2 Scope of the Thesis 12 1.3 Thesis Outline 14 2. Runge-Kutta Methods 16 2.1 Integration Methods in General 16 2.2 Runge-Kutta Methods 23 2.3 The Standard Stepsize Selection Rule 37 2.4 Equation Solvers 38 2.5 Standard Error and Convergence Control 44 2.6 Control Issues in Runge-Kutta Methods 47 3. The Stepsize-Error Relation 50 3.1 The Limitations of the Asymptotic Model 50 3.2 Different Operating Regions 58 3.3 An Improved Asymptotic Model 61 3.4 Outside the Asymptotic Region 64 3.5 Stepsize Limited by Stability 67 3.6 Experimental Verification of the Models 74 3.7 Conclusions 82 4. Control of Explicit Runge-Kutta Methods 84 4.1 Problems with the Standard Controller 84 4.2 Control Objectives 89 4.3 Changing the Controller Structure 90 4.4 Tuning the Controller Parameters 93 4.5 Restart Strategies 97 4.6 Numerical Examples 99 4.7 The Complete Controller 110 5. Control of Implicit Runge-Kutta Methods 113 5.1 A Predictive Error Controller 115 5.2 The Iteration Error in the Equation Solver 129 5.3 Convergence Control 134 5.4 The Complete Controller 153 6. Conclusions 158 7. References 163 A. Some Common Runge-Kutta Methods 171 A.1 Explicit Methods 172 A.2 Implicit Methods 177 A.3 Method Properties 179 B. List of Notations 183 Preface "It's one small step for man, one giant leap for mankind." Neil Armstrong The size of a step should be related to its environment. This is partic- ularly true in numerical integration of ordinary differential equations, where the stepsize is used to control the quality of the produced solu- tion. A large stepsize leads to a large error, while a too small stepsize is inefficient due to the many steps needed to carry out the integration. What is regarded as large or small, depends on the differential equation and the required accuracy. "The choice of stepsize in an integration method is a problem you control guys should have a look at," said Professor Gustaf Söderlind, when handing out project suggestions at the end of a course on numeri- cal integration methods in 1986/87. Michael Lundh and I got interested in the problem and started working. We soon realized that the standard rule for selecting stepsize could be interpreted as a pure integrating con- troller. A generalization to PI or PID control is then rather reasonable. We experimented with PI control and achieved good results [Gustafsson et a]., 1988]. After having finished the project both Michael and I continued work- ing on other research topics. I was, however, rather annoyed that we did not have a good explanation of why a PI controller made such a differ- ence in performance in the case where numerical stability restricts the stepsize. What was it in the integration method that changed so dra- matically? After a couple of months I returned to the problem, strongly encouraged by Professors Karl Johan Åström, Gustaf Söderlind, and Do- cent Per Hagander. Soon I was hooked. Pursuing the feedback control viewpoint turned out to be very rewarding. Not only was it possible to ex- plain the superiority of the PI controller, but the approach also provided means to ..., well, continue on reading, and you will learn all about it. Acknowledgements This will probably be the most read section of this thesis. I have looked forward to writing it, since it gives me an opportunity to express my gratitude to all dear friends. They have been of invaluable assistance during the years I have been working with this material. I have benefited immensely from their feedback on my ideas, from their help to get the right perspectives, and from their encouragement and support in times of intellectual standstill. Thank you! There are a few friends I would like to mention in particular: Working in an interdisciplinary area I have enjoyed the pleasure of several supervisors. I want to express my gratitude to Karl Johan Åström for teaching me about control, and enthusiastically backing me up when I have tried to apply feedback ideas in a somewhat different area. He has taught me to pose the right questions, and I have benefited immensely from his vast research experience. With admirable teaching skill Gustaf Söderlind introduced me to the subject of numerical integra- tion. He has been a tremendous source ofinformation and has teasingly presented one challenging research problem after another. I am also in- debted to Per Hagander. He has provided many constructive pieces of advice. I appreciate him for taking the time to penetrate even minute details, and for allowing ample time for all our discussions. Among my fellow graduate students I would like to mention Michael Lundh, Bo Bernhardsson, and Ola Dahl. The early part of the work was done together with Michael and I am grateful for his contributions. Bo has been a dear and close friend for more than a decade. We have a common fascination for life in general and engineering in particular. Thanks for sharing your insight and ideas with me. Thanks Ola for surviving in the same office with me, patiently ove^ooking me blocking the computer and spreading my things everywhere. My colleagues at the Department of Automatic Control in Lund form a friendly and inspiring community. The creative and open-minded at- mosphere has influenced me a lot. In addition, the extra-ordinary access to computing power and high quality software, excellently maintained by Leif Andersson and Anders Blomdell, has made work very stimulating. I have plagued Gustaf Söderlind, Karl Johan Åström, Per Hagan- der, and Björn Wittenmark with several versions of the manuscript. Their constructive criticism have greatly improved the quality of the final product.
Load more

Recommended publications
  • Coping with the Bounds: a Neo-Clausewitzean Primer / Thomas J
    About the CCRP The Command and Control Research Program (CCRP) has the mission of improving DoD’s understanding of the national security implications of the Information Age. Focusing upon improving both the state of the art and the state of the practice of command and control, the CCRP helps DoD take full advantage of the opportunities afforded by emerging technologies. The CCRP pursues a broad program of research and analysis in information superiority, information operations, command and control theory, and associated operational concepts that enable us to leverage shared awareness to improve the effectiveness and efficiency of assigned missions. An important aspect of the CCRP program is its ability to serve as a bridge between the operational, technical, analytical, and educational communities. The CCRP provides leadership for the command and control research community by: • articulating critical research issues; • working to strengthen command and control research infrastructure; • sponsoring a series of workshops and symposia; • serving as a clearing house for command and control related research funding; and • disseminating outreach initiatives that include the CCRP Publication Series. This is a continuation in the series of publications produced by the Center for Advanced Concepts and Technology (ACT), which was created as a “skunk works” with funding provided by the CCRP under the auspices of the Assistant Secretary of Defense (NII). This program has demonstrated the importance of having a research program focused on the national security implications of the Information Age. It develops the theoretical foundations to provide DoD with information superiority and highlights the importance of active outreach and dissemination initiatives designed to acquaint senior military personnel and civilians with these emerging issues.
    [Show full text]
  • Numerical Solution of Ordinary Differential Equations
    NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS Kendall Atkinson, Weimin Han, David Stewart University of Iowa Iowa City, Iowa A JOHN WILEY & SONS, INC., PUBLICATION Copyright c 2009 by John Wiley & Sons, Inc. All rights reserved. Published by John Wiley & Sons, Inc., Hoboken, New Jersey. Published simultaneously in Canada. No part of this publication may be reproduced, stored in a retrieval system, or transmitted in any form or by any means, electronic, mechanical, photocopying, recording, scanning, or otherwise, except as permitted under Section 107 or 108 of the 1976 United States Copyright Act, without either the prior written permission of the Publisher, or authorization through payment of the appropriate per-copy fee to the Copyright Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, (978) 750-8400, fax (978) 646-8600, or on the web at www.copyright.com. Requests to the Publisher for permission should be addressed to the Permissions Department, John Wiley & Sons, Inc., 111 River Street, Hoboken, NJ 07030, (201) 748-6011, fax (201) 748-6008. Limit of Liability/Disclaimer of Warranty: While the publisher and author have used their best efforts in preparing this book, they make no representations or warranties with respect to the accuracy or completeness of the contents of this book and specifically disclaim any implied warranties of merchantability or fitness for a particular purpose. No warranty may be created ore extended by sales representatives or written sales materials. The advice and strategies contained herin may not be suitable for your situation. You should consult with a professional where appropriate. Neither the publisher nor author shall be liable for any loss of profit or any other commercial damages, including but not limited to special, incidental, consequential, or other damages.
    [Show full text]
  • Observation of Nonlinear Verticum-Type Systems Applied to Ecological Monitoring
    2nd Reading June 4, 2012 13:53 WSPC S1793-5245 242-IJB 1250051 International Journal of Biomathematics Vol. 5, No. 6 (November 2012) 1250051 (15 pages) c World Scientific Publishing Company DOI: 10.1142/S1793524512500519 OBSERVATION OF NONLINEAR VERTICUM-TYPE SYSTEMS APPLIED TO ECOLOGICAL MONITORING S. MOLNAR´ Institute of Mathematics and Informatics Szent Istv´an University P´ater K. u. 1., H-2103 Godollo, Hungary [email protected] M. GAMEZ´ ∗ and I. LOPEZ´ † Department of Statistics and Applied Mathematics University of Almer´ıa La Ca˜nada de San Urbano 04120 Almer´ıa, Spain ∗[email protected] †[email protected] Received 4 November 2011 Accepted 13 January 2012 Published 7 June 2012 In this paper the concept of a nonlinear verticum-type observation system is introduced. These systems are composed from several “subsystems” connected sequentially in a particular way: a part of the state variables of each “subsystem” also appears in the next “subsystem” as an “exogenous variable” which can also be interpreted as a con- trol generated by an “exosystem”. Therefore, these “subsystems” are not observation systems, but formally can be considered as control-observation systems. The problem of observability of such systems can be reduced to rank conditions on the “subsystems”. Indeed, under the condition of Lyapunov stability of an equilibrium of the “large”, verticum-type system, it is shown that the Kalman rank condition on the linearization of the “subsystems” implies the observability of the original, nonlinear verticum-type system. For an illustration of the above linearization result, a stage-structured fishery model with reserve area is considered.
    [Show full text]
  • Nonlinear System Stability and Behavioral Analysis for Effective
    S S symmetry Article Nonlinear System Stability and Behavioral Analysis for Effective Implementation of Artificial Lower Limb Susmita Das 1 , Dalia Nandi 2, Biswarup Neogi 3 and Biswajit Sarkar 4,* 1 Narula Institute of Technology, Kolkata 700109, India; [email protected] 2 Indian Institute of Information Technology Kalyani, Kalyani 741235, India; [email protected] 3 JIS College of Engineering., Kalyani 741235, India; [email protected] 4 Department of Industrial Engineering, Yonsei University, Seoul 03722, Korea * Correspondence: [email protected] or [email protected] Received: 8 September 2020; Accepted: 14 October 2020; Published: 19 October 2020 Abstract: System performance and efficiency depends on the stability criteria. The lower limb prosthetic model design requires some prerequisites such as hardware design functionality and compatibility of the building block materials. Effective implementation of mathematical model simulation symmetry towards the achievement of hardware design is the focus of the present work. Different postures of lower limb have been considered in this paper to be analyzed for artificial system design of lower limb movement. The generated polynomial equations of the sitting and standing positions of the normal limb are represented with overall system transfer function. The behavioral analysis of the lower limb model shows the nonlinear nature. The Euler-Lagrange method is utilized to describe the nonlinearity in the field of forward dynamics of the artificial system. The stability factor through phase portrait analysis is checked with respect to nonlinear system characteristics of the lower limb. The asymptotic stability has been achieved utilizing the most applicable Lyapunov method for nonlinear systems.
    [Show full text]
  • Spectral Properties of the Koopman Operator in the Analysis of Nonstationary Dynamical Systems
    UNIVERSITY of CALIFORNIA Santa Barbara Spectral Properties of the Koopman Operator in the Analysis of Nonstationary Dynamical Systems A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Mechanical Engineering by Ryan M. Mohr Committee in charge: Professor Igor Mezi´c,Chair Professor Bassam Bamieh Professor Jeffrey Moehlis Professor Mihai Putinar June 2014 The dissertation of Ryan M. Mohr is approved: Bassam Bamieh Jeffrey Moehlis Mihai Putinar Igor Mezi´c,Committee Chair May 2014 Spectral Properties of the Koopman Operator in the Analysis of Nonstationary Dynamical Systems Copyright c 2014 by Ryan M. Mohr iii To my family and dear friends iv Acknowledgements Above all, I would like to thank my family, my parents Beth and Jim, and brothers Brennan and Gralan. Without their support, encouragement, confidence and love, I would not be the person I am today. I am also grateful to my advisor, Professor Igor Mezi´c,who introduced and guided me through the field of dynamical systems and encouraged my mathematical pursuits, even when I fell down the (many) proverbial (mathematical) rabbit holes. I learned many fascinating things from these wandering forays on a wide range of topics, both contributing to my dissertation and not. Our many interesting discussions on math- ematics, research, and philosophy are among the highlights of my graduate school career. His confidence in my abilities and research has been a constant source of encouragement, especially when I felt uncertain in them myself. His enthusiasm for science and deep, meaningful work has set the tone for the rest of my career and taught me the level of research problems I should consider and the quality I should strive for.
    [Show full text]
  • Polarization Fields and Phase Space Densities in Storage Rings: Stroboscopic Averaging and the Ergodic Theorem
    Physica D 234 (2007) 131–149 www.elsevier.com/locate/physd Polarization fields and phase space densities in storage rings: Stroboscopic averaging and the ergodic theorem✩ James A. Ellison∗, Klaus Heinemann Department of Mathematics and Statistics, University of New Mexico, Albuquerque, NM 87131, United States Received 1 May 2007; received in revised form 6 July 2007; accepted 9 July 2007 Available online 14 July 2007 Communicated by C.K.R.T. Jones Abstract A class of orbital motions with volume preserving flows and with vector fields periodic in the “time” parameter θ is defined. Spin motion coupled to the orbital dynamics is then defined, resulting in a class of spin–orbit motions which are important for storage rings. Phase space densities and polarization fields are introduced. It is important, in the context of storage rings, to understand the behavior of periodic polarization fields and phase space densities. Due to the 2π time periodicity of the spin–orbit equations of motion the polarization field, taken at a sequence of increasing time values θ,θ 2π,θ 4π,... , gives a sequence of polarization fields, called the stroboscopic sequence. We show, by using the + + Birkhoff ergodic theorem, that under very general conditions the Cesaro` averages of that sequence converge almost everywhere on phase space to a polarization field which is 2π-periodic in time. This fulfills the main aim of this paper in that it demonstrates that the tracking algorithm for stroboscopic averaging, encoded in the program SPRINT and used in the study of spin motion in storage rings, is mathematically well-founded.
    [Show full text]
  • Approximations in Numerical Analysis, Cont'd
    Jim Lambers MAT 460/560 Fall Semester 2009-10 Lecture 7 Notes These notes correspond to Section 1.3 in the text. Approximations in Numerical Analysis, cont'd Convergence Many algorithms in numerical analysis are iterative methods that produce a sequence fαng of ap- proximate solutions which, ideally, converges to a limit α that is the exact solution as n approaches 1. Because we can only perform a finite number of iterations, we cannot obtain the exact solution, and we have introduced computational error. If our iterative method is properly designed, then this computational error will approach zero as n approaches 1. However, it is important that we obtain a sufficiently accurate approximate solution using as few computations as possible. Therefore, it is not practical to simply perform enough iterations so that the computational error is determined to be sufficiently small, because it is possible that another method may yield comparable accuracy with less computational effort. The total computational effort of an iterative method depends on both the effort per iteration and the number of iterations performed. Therefore, in order to determine the amount of compu- tation that is needed to attain a given accuracy, we must be able to measure the error in αn as a function of n. The more rapidly this function approaches zero as n approaches 1, the more rapidly the sequence of approximations fαng converges to the exact solution α, and as a result, fewer iterations are needed to achieve a desired accuracy. We now introduce some terminology that will aid in the discussion of the convergence behavior of iterative methods.
    [Show full text]
  • Rethinking the International System As a Complex Adaptive System
    Munich Personal RePEc Archive A New Taxonomy for International Relations: Rethinking the International System as a Complex Adaptive System Scartozzi, Cesare M. University of Tokyo, Graduate School of Public Policy 2018 Online at https://mpra.ub.uni-muenchen.de/95496/ MPRA Paper No. 95496, posted 12 Aug 2019 10:50 UTC 1SFQVCMJDBUJPOJournal on Policy and Complex Systems • Volume 4, Number 1 • Spring 2018 A New Taxonomy for International Relations: Rethinking the International System as a Complex Adaptive System Cesare Scartozzi Graduate School of Public Policy, University of Tokyo [email protected] 7-3-1, Hongo, Bunkyo, Tokyo 113-0033, Japan ORCID number: 0000-0002-4350-4386. Abstract The international system is a complex adaptive system with emer- gent properties and dynamics of self-organization and informa- tion processing. As such, it is better understood with a multidis- ciplinary approach that borrows methodologies from the field of complexity science and integrates them to the theoretical perspec- tives offered by the field of international relations (IR). This study is set to formalize a complex systems theory approach to the study of international affairs and introduce a new taxonomy for IR with the two-pronged aim of improving interoperability between differ- ent epistemological communities and outlining a formal grammar that set the basis for modeling international politics as a complex adaptive system. Keywords: international politics, international relations theory, complex systems theory, taxonomy, adaptation, fitness, self-orga- nization This is a prepublication version of: Scartozzi, Cesare M. “A New Taxonomy for International Relations: Rethinking the International System as a Complex Adaptive System.” Journal on Policy and Complex Systems 4, no.
    [Show full text]
  • Chapter 8 Nonlinear Systems
    Chapter 8 Nonlinear systems 8.1 Linearization, critical points, and equilibria Note: 1 lecture, §6.1–§6.2 in [EP], §9.2–§9.3 in [BD] Except for a few brief detours in chapter 1, we considered mostly linear equations. Linear equations suffice in many applications, but in reality most phenomena require nonlinear equations. Nonlinear equations, however, are notoriously more difficult to understand than linear ones, and many strange new phenomena appear when we allow our equations to be nonlinear. Not to worry, we did not waste all this time studying linear equations. Nonlinear equations can often be approximated by linear ones if we only need a solution “locally,” for example, only for a short period of time, or only for certain parameters. Understanding linear equations can also give us qualitative understanding about a more general nonlinear problem. The idea is similar to what you did in calculus in trying to approximate a function by a line with the right slope. In § 2.4 we looked at the pendulum of length L. The goal was to solve for the angle θ(t) as a function of the time t. The equation for the setup is the nonlinear equation L g θ�� + sinθ=0. θ L Instead of solving this equation, we solved the rather easier linear equation g θ�� + θ=0. L While the solution to the linear equation is not exactly what we were looking for, it is rather close to the original, as long as the angleθ is small and the time period involved is short. You might ask: Why don’t we just solve the nonlinear problem? Well, it might be very difficult, impractical, or impossible to solve analytically, depending on the equation in question.
    [Show full text]
  • NONLINEARITIES, FEEDBACKS and CRITICAL THRESHOLDS WITHIN the EARTH's CLIMATE SYSTEM 1. Introduction Nonlinear Phenomena Charac
    NONLINEARITIES, FEEDBACKS AND CRITICAL THRESHOLDS WITHIN THE EARTH’S CLIMATE SYSTEM JOSÉ A. RIAL 1,ROGERA.PIELKESR.2, MARTIN BENISTON 3, MARTIN CLAUSSEN 4, JOSEP CANADELL 5, PETER COX 6, HERMANN HELD 4, NATHALIE DE NOBLET-DUCOUDRÉ 7, RONALD PRINN 8, JAMES F. REYNOLDS 9 and JOSÉ D. SALAS 10 1Wave Propagation Laboratory, Department of Geological Sciences CB#3315, University of North Carolina, Chapel Hill, NC 27599-3315, U.S.A. E-mail: [email protected] 2Atmospheric Science Dept., Colorado State University, Fort Collins, CO 80523, U.S.A. 3Dept. of Geosciences, Geography, Univ. of Fribourg, Pérolles, Ch-1700 Fribourg, Switzerland 4Potsdam Institute for Climate Impact Research, Telegrafenberg C4, 14473 Potsdam, P.O. Box 601203, Potsdam, Germany 5GCP-IPO, Earth Observation Centre, CSIRO, GPO Box 3023, Canberra, ACT 2601, Australia 6Met Office Hadley Centre, London Road, Bracknell, Berkshire RG12 2SY, U.K. 7DSM/LSCE, Laboratoire des Sciences du Climat et de l’Environnement, Unité mixte de Recherche CEA-CNRS, Bat. 709 Orme des Merisiers, 91191 Gif-sur-Yvette, France 8Dept. of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139-4307, U.S.A. 9Department of Biology and Nicholas School of the Environmental and Earth Sciences, Phytotron Bldg., Science Dr., Box 90340, Duke University, Durham, NC 27708, U.S.A. 10Dept. of Civil Engineering, Colorado State University, Fort Collins, CO 80523, U.S.A. Abstract. The Earth’s climate system is highly nonlinear: inputs and outputs are not proportional, change is often episodic and abrupt, rather than slow and gradual, and multiple equilibria are the norm.
    [Show full text]
  • Modeling Vascular Homeostasis and Improving Data Filtering Methods for Model Calibration
    Modeling Vascular Homeostasis and Improving Data Filtering Methods for Model Calibration by Jiacheng Wu A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Engineering − Mechanical Engineering and the Designated Emphasis in Computational Science and Engineering in the Graduate Division of the University of California, Berkeley Committee in charge: Professor Shawn C. Shadden, Chair Professor Grace D. O'Connell Professor Peter L. Bartlett Spring 2019 Modeling Vascular Homeostasis and Improving Data Filtering Methods for Model Calibration Copyright 2019 by Jiacheng Wu 1 Abstract Modeling Vascular Homeostasis and Improving Data Filtering Methods for Model Calibration by Jiacheng Wu Doctor of Philosophy in Engineering − Mechanical Engineering and the Designated Emphasis in Computational Science and Engineering University of California, Berkeley Professor Shawn C. Shadden, Chair Vascular homeostasis is the preferred state that blood vessels try to maintain against external mechanical and chemical stimuli. The vascular adaptive behavior around the homeostatic state is closely related to cardiovascular disease progressions such as arterial aneurysms. In this work, we develop a multi-physics computational framework that couples vascular growth & remodeling (G&R), wall mechanics and hemodynamics to describe the overall vascular adaptive behavior. The coupled simulation is implemented in patient-specific geometries to predict aneurysm progression. Lyapunov stability analysis of the governing
    [Show full text]
  • Nonlinear Control Systems
    1. Introduction to Nonlinear Systems Nonlinear Control Systems Ant´onioPedro Aguiar [email protected] 1. Introduction to Nonlinear Systems IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7Epedro/NCS2012/ 2011/2012 1 1. Introduction to Nonlinear Systems Objective The main goal of this course is to provide to the students a solid background in analysis and design of nonlinear control systems Why analysis? (and not only simulation) • Every day computers are becoming more and more powerful to simulate complex systems • Simulation combined with good intuition can provide useful insight into system's behavior Nevertheless • It is not feasible to rely only on simulations when trying to obtain guarantees of stability and performance of nonlinear systems, since crucial cases may be missed • Analysis tools provide the means to obtain formal mathematical proofs (certificates) about the system's behavior • results may be surprising, i.e, something we had not thought to simulate. 2 1. Introduction to Nonlinear Systems Why study nonlinear systems? Nonlinear versus linear systems • Huge body of work in analysis and control of linear systems • most models currently available are linear (but most real systems are nonlinear...) However • dynamics of linear systems are not rich enough to describe many commonly observed phenomena 3 1. Introduction to Nonlinear Systems Examples of essentially nonlinear phenomena • Finite escape time, i.e, the state can go to infinity in finite time (while this is impossible to happen for linear systems) • Multiple isolated equilibria, while linear systems can only have one isolated equilibrium point, that is, one steady state operating point • Limit cycles (oscillation of fixed amplitude and frequency, irrespective of the initial state) • Subharmonic, harmonic or almost-periodic oscillations; • A stable linear system under a periodic input produces an output of the same frequency.
    [Show full text]