OPTIMAL CONTROL OF A CLASS OF PARABOLIC
|
|
- Anthony Hicks
- 5 years ago
- Views:
Transcription
1 OPTIMAL CONTROL OF A CLASS OF PARABOLIC TIME-VARYING PDES James Ng a, Ilyasse Aksikas a,b, Stevan Dubljevic a a Department of Chemical & Materials Engineering, University of Alberta, Canada b Department of Mathematics, King Abdelaziz University, Jeddah, KSA 7 th Workshop on control of Distributed Parameter Systems Wuppertal, Germany July 18-22, 2011
2 INTRODUCTION MOTIVATING EXAMPLES Catalytic tubular reactor Tc T D2 ξ = v (T in T(0, t)) ξ=0 CA D1 ξ = v (CA in CA(0, t)) ξ=0 r ξ T ξ = 0 ξ=l CA ξ = 0 ξ=l Figure: Catalytic reactor Many industrial processes, e.g methanol, ammonia and other petrochemicals Diffusion-Convection-Reaction Process Tubular reactor systems with catalyst deactivation, Loss of catalyst activity Time-varying rates of reaction = Parabolic time-varying PDEs
3 INTRODUCTION MOTIVATING EXAMPLES Crystal Growth Process Important industrial process utilized for the production of semi-conductor material in the electronics and microprocessor industry. Materials produced: Silicon (Si), Germanium (Ge). Temperature dynamics: Parabolic PDE with time-varying coefficients Convective transport term is time-varying due to the motion of the domain boundary. Figure: Crystal Process diagram
4 INTRODUCTION MOTIVATING EXAMPLES Parabolic partial differential equations (PDEs) with time-varying features represent an important class of models for reaction-diffusion-convection processes. e.g. Tubular and packed bed reactor systems with catalyst deactivation, Crystal growth and annealing type processes with time-varying spatial domains These time-dependent features play an important role in the system dynamics, and therefore must be incorporated into the model based controller design. * Approach: Evolution systems representation Operator differential Riccati equation
5 INTRODUCTION RELATED WORKS Nonautonomous PDEs I. Aksikas, J.F. Forbes, and Y. Belhamadia, Optimal control design for time-varying catalytic reactors: a Riccati equation based approach, Int. J. Control., 2009 I. Aksikas and J.F. Forbes, Linear quadratic regulator for time-varying hyperbolic distributed parameter systems, IMA. J. Mathematical Control and Information, P. Acquistapace, F. Flandoli, and B. Terreni, Initial boundary value problems and optimal control for nonautonomous parabolic systems, SIAM J. Cont. & Optim., A. Smyshlyaev and M. Krstic, On control design for PDEs with space-dependent diffusivity and time-dependent reactivity, Automatica, PDEs with time-varying spatial domains A. Armaou and P. D. Christofides, Robust control of parabolic PDE systems with time-dependent spatial domains, Automatica, P.K.C.Wang, Stabilization and control of distributed systems with time-dependent spatial domains, J. Optim. Theor. & Appl., 1990.
6 GENERAL MODEL PDE DESCRIPTION Let Ω be a bounded open set of R m with smooth boundary Ω. Consider the initial and boundary value problem of the form: z(ξ, t) + A(t)z(ξ, t) = f (ξ, t) t in Ω [0, T] z(ξ, 0) = z 0 (ξ) in Ω z(ξ, t) n = 0, on Ω [0, T] The family of operators A(t) is defined as: m A(t)z := a ij (ξ, t) 2 z m + b i (t) z + c(ξ, t)z ξ i,j=1 i ξ j ξ i=1 i (i) z(ξ, t) represents, for example, temperature or concentration, with initial distribution z 0 (ξ). (ii) a ij (ξ, t) describes the heterogeneous thermal conductivity or diffusivity. (iii) b i (t) is a convective transport coefficient (e.g. time-dependent fluid superficial velocity). (iv) c(ξ, t) is a linearized reaction term (e.g. due to catalyst deactivation).
7 GENERAL MODEL PDE OPERATOR PROPERTIES Assumptions: P1. For each t [0, T], the operator A(t) is strongly elliptic, i.e. m a ij (ξ, t)η i η j ε η 2, for η R m i,j=1 P2. The coefficients c(ξ, t) L 2 ([0, T], L 2 (Ω)), b i (t) C 1 ([0, T]) and a ij (ξ, t) are sufficiently Hölder continuous, i.e. a ij (ξ, t) a ij (ξ, s) L t s β for s, t [0, T], ξ Ω and constant L > 0 and β (0, 1]. P3. The function f (ξ, t) L 2 (Ω) satisfies: ( 1 f (ξ, t) f (ξ, s) dξ) 2 2 L t s β, 0 s < t T Ω {A(t)} t [0,T] forms a family of strongly elliptic operators which admit a family of eigenfunctions {φ n(t)} t [0,T] with corresponding family of eigenvalues {λ n(t)} t [0,T].
8 INFINITE-DIMENSIONAL SYSTEM REPRESENTATION NONAUTONOMOUS PARABOLIC EVOLUTION SYSTEM Under the properties of: strong ellipticity and continuity of a ij (ξ, t), b i (t) and c(ξ, t), the operator A(t) satisfies: 1. For every t [0, T], the resolvent R(λ; A(t)) exists for all λ with Reλ 0 and there exists a constant L 1 such that R(λ; A(t)) L 1 / λ ; 2. There exists constants L 2 and β (0, 1] such that (A(t) A(s))A(τ) 1 L 2 t s β for s, t, τ [0, T]. The initial and boundary value problem is represented as a non-autonomous evolution system on L 2 (Ω): dz(t) dt for 0 s < t T and z s L 2 (Ω). = A(t)z(t) + f (t), z(s) = z s
9 INFINITE-DIMENSIONAL SYSTEM REPRESENTATION TWO-PARAMETER SEMIGROUP The solution of the nonautonomous evolution system is expressed in the form of: t z(t) = U(t, s)z s + U(t, τ)f (τ)dτ s where U(t, s) is the two-parameter semigroup generated by the operator A(t). The operator A(t) is defined as: A(t) := λ n(t), ψ n(t) φ n(t) n=1 The operator A(t) : D(A(t)) L 2 (Ω) L 2 (Ω) is the infinitesimal generator of the two-parameter semigroup: U(t, s)z(s) := e µn(t) µn(s) z(s), ψ n(s) φ n(t), for 0 s t T n=1
10 INFINITE-DIMENSIONAL SYSTEM REPRESENTATION TWO-PARAMETER SEMIGROUP The operator U(t, s) satisfies the following identities: A1. U(t, t) = I, A2. U(t, s) = U(t, r)u(r, s) for 0 s r t T A3. U(t, s) is continuous on 0 s < t T, and: U(t, s) t A4. U(t, s) L 1, = A(t)U(t, s), and A5. A(t)U(t, s) L 2 (t s) 1, and U(t, s) s A6. A(t)U(t, s)a(s) 1 L 3 for constants L i > 0. = U(t, s)a(s)
11 INFINITE-DIMENSIONAL SYSTEM REPRESENTATION OPTIMAL CONTROL PROBLEM The finite-time horizon optimal control problem is given as the following: T ( min J(u) = min C(τ)z(τ) 2 + u(τ) 2) dτ + Qz(T), z(t) u u 0 where the functional J(u) is minimized over all trajectories of { ẋ(t) = A(t)x(t) + B(t)u(t) x(0) = x 0 (1) If Q 0, B( ) C([0, T; L(U; L 2 (Ω))) and C( ) C([0, T; L(Y, L 2 (Ω))), then the optimization problem has the unique minimizing solution u min (t) such that the optimal pair u min (t) C([0, T]; U) and z min (t) C 1 ([0, T]; X) C([0, T]; D(A(t)) are related by the feedback formula, u min t [0,T] (t) = B (t)π(t)z min (t) The operator Π(t) L(X) is the strongly continuous, self adjoint, nonnegative solution of the differential Riccati equation, Π(t) + A (t)π(t) + Π(t)A(t) Π(t)B(t)B (t)π(t) + C (t)c(t) = 0, Π(T) = Q
12 EXAMPLE: CRYSTAL GROWTH PROCESS PROCESS MODEL Model of temperature dynamics: z(ξ, t) t = κ 2 z z v(t) ξ2 ξ + b(ξ)u(t) ξ = 0 z z (0, t) = 0, ξ z(ξ, 0) = z 0 (ξ) (l(t), t) = 0 ξ f(ξ, t) w(ξ, t) Q(ξ, t) κ thermal conductivity constant, v(t) is the boundary velocity, and b(ξ)u(t) is the distributed heat input along the length of the rod. The material domain is time-varying due to the motion of the boundary at the material-fluid interface, ξ = l(t). ξ = l(t) Crystal Melt Control objective: Optimally stabilize the temperature around some desired nominal distribution. Crucible
13 EXAMPLE: CRYSTAL GROWTH PROCESS TIME-VARYING SPATIAL DOMAIN AND FUNCTION SPACE DESCRIPTION The time-varying spatial domain at some time t [0, T] is denoted Ω = (0, l(t)) with 0 < ξ < l(t) l max. Spatial domain evolution is considered as a sequence of subdomains: Ω j Ω j+1 Ω Let φ(ξ, t) denote a family of functions defined on the subdomains Ω j with: φ(ξ, t) = { φ(ξ) for ξ Ωj 0 for ξ Ω c j R L 2 (Ω j ) forms a family of function spaces which are precompact in L 2 (Ω) for all t [0, T], (see, Adams, 1975): L 2 (Ω j ) L 2 (Ω j+1 ) L 2 (Ω) Enables the use of single inner product, on L 2 (Ω) for functions defined on an arbitrary subdomain Ω: φ, φ L 2 (Ω) = Ω φ(ξ, t)φ(ξ, t)dξ = φ(ξ)φ(ξ)dξ + 0 dξ = φ, φ L 2 Ω Ω c (Ω)
14 EXAMPLE: CRYSTAL GROWTH PROCESS PDE OPERATOR PROPERTIES The PDE operator A(t) is defined as: A(t) := κ 2 z z v(t) ξ2 ξ The operator A(t) is strongly elliptic for each t [0, T]; For each t [0, T] the family of eigenfunctions {φ n(t)} [0,T] are: ( ( ) φ n(ξ, t) = B n(t)e 2 1 κ 1 v(t)ξ nπ cos l(t) ξ 1 ( )) v(t) nπ 2 κ 1 (nπ/l(t)) sin l(t) ξ with coefficients: ( ( ) ) 2 v(t) B n(t) = 1 + l(t) 2κ (nπ/l(t)) φ n(ξ, t) are orthonormal to the eigenfunctions: ψ n(ξ, t) = e κ 1 v(t)ξ φ n(ξ, t) of the adjoint operator A (ξ, t) for each t [0, T]. The corresponding family of eigenvalues are: ( ) nπ 2 λ n(t) = κ 1 v(t)2 κ 1 l(t) 2 2
15 EXAMPLE: CRYSTAL GROWTH PROCESS INFINITE-DIMENSIONAL SYSTEM REPRESENTATION The associated family of linear operators A(t) in L 2 (Ω) is defined as: with domain D(A(t)): D(A(t)) := { φ L 2 (Ω) : φ, A(t)z = A(ξ, t)z for z D(A(t)) φ ξ are a.c., A(t)φ L2 (Ω), and φ } ξ (0, t) = 0, φ (l(t), t) = 0 ξ Initial and boundary value control problem is represented as the nonautonomous parabolic initial value problem: dz dt = A(t)z(t) + B(t)u(t), z(0) = z 0, 0 s t < T The solution of the initial value problem is expressed in terms of the two parameter semigroup U(t, s), t z(t) = U(t, 0)z 0 + U(t, τ)b(τ)u(τ)dτ 0 with z(s) L 2 (Ω).
16 EXAMPLE: CRYSTAL GROWTH PROCESS INFINITE-DIMENSIONAL SYSTEM REPRESENTATION The operator A(t) : D(A(t)) L 2 (Ω) L 2 (Ω) generates the two-parameter semigroup U(t, s), 0 s t T expressed as: U(t, s)z(s) = n=1 { ( exp κ 0 (nπ) 2 t l(t) 2 s ) l(s) 2 1 4κ 0 ( tv(t) 2 sv(s) 2)} z(s), ψ n(s) φ n(t) One can note that in the case where the boundary motion ceases and Ω is constant, i.e. l(t) = l, v(t) = 0, the expression for U(t, s) becomes: { ( ) nπ 2 U(t, s)z(s) = exp κ 0 (t s)} z(s), φ φ = T(t s)z(s) l n=1 where T(t), t 0 is the C 0 -semigroup of operators on L 2 (Ω) which is generated by the standard heat equation.
17 EXAMPLE: CRYSTAL GROWTH PROCESS CONTROLLER FORMULATION Let consider the weighted inner product on L 2 (Ω), for z 1, z 2 D(A(t)), z 1, z 2 r = r(ξ, t)z 1 (ξ, t)z 2 (ξ, t)dξ Ω with weight function r(ξ, t) := exp( (v(t)/κ 0 )ξ). Note that φ n, φ m r = δ nm where δ nm = 1 if n = m and 0 otherwise. The differential Riccati equation can be written under the form: φ n(t), Π(t)φ m(t) r + A(t)φ n(t), Π(t)φ m(t) r + φ n(t), Π(t)A(t)φ m(t) r Π(t)B(t)B (t)π(t)φ n(t), φ m(t) r + C(t)φ n(t), C(t)φ m(t) r = 0 with φ n(t), Π(T)φ m(t) r = φ n(t), Qφ m(t) r.
18 EXAMPLE: CRYSTAL GROWTH PROCESS CONTROLLER FORMULATION In the case where B(t) = I and C(t) = I the differential Riccati equation becomes the system of infinitely many ordinary differential equations: Π nn(t) + 2λ n(t)π nn(t) + 1 Π 2 nn (t) = 0, Πnn(T) = Q The input u min is determined from the solution of the system of ODEs, and the t [0,T] optimal state trajectory is the mild solution of the state feedback system and is expressed as: ż = (A(t) B(t)B (t)π(t))z(t), z(0) = z 0 t z(t) = e µn(t) z 0, ψ n(0) φ n(t) U(t, τ) n=1 0 n Π nn(t) z 0 (τ), ψ n(τ) φ n(τ)dτ
19 NUMERICAL RESULTS HEAT INPUT AND TEMPERATURE DISTRIBUTION l(t) v(t) t Figure: Slab domain length l(t) and boundary velocity v(t). t u(t) z(ξ, t) t t closed-loop open-loop Figure: (Top) Optimal input profile u min (t) applied to the slab at input location ξ c = (Bottom) Total open and closed loop system energy.
20 NUMERICAL RESULTS CLOSED LOOP SYSTEM z(ξ, t) t ξ 3.0 Figure: Slab temperature evolution in the time-dependent spatial domain with diffusivity κ 0 = Input applied at ξ c =
21 SUMMARY CONCLUDING REMARKS The general two-parameter semigroup representation of a class of nonautonomous parabolic PDE has been presented. The optimal control formulation for the infinite-dimensional system representation is considered. A practical example of an crystal growth process is considered with PDE model defined on time-varying spatial domain. The infinite-dimensional system representation of the PDE is determined, and the explicit two-parameter semigroup expression is provided. The corresponding optimal control problem is considered and numerical results demonstrate the stabilization of the temperature distribution in the time-dependent region.
Control Design of a Distributed Parameter. Fixed-Bed Reactor
Workshop on Control of Distributed Parameter Systems p. 1/1 Control Design of a Distributed Parameter Fixed-Bed Reactor Ilyasse Aksikas and J. Fraser Forbes Department of Chemical and Material Engineering
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationControl of Time-Varying Distributed Parameter Plug Flow Reactor by LQR
Proceedings of the 17th World Congress The International Federation of Automatic Control Control of Time-Varying Distributed Parameter Plug Flow Reactor by LQR Ilyasse Aksikas, Adrian M. Fuxman, J. Fraser
More informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
More informationStabilization of Distributed Parameter Systems by State Feedback with Positivity Constraints
Stabilization of Distributed Parameter Systems by State Feedback with Positivity Constraints Joseph Winkin Namur Center of Complex Systems (naxys) and Dept. of Mathematics, University of Namur, Belgium
More informationPositive Stabilization of Infinite-Dimensional Linear Systems
Positive Stabilization of Infinite-Dimensional Linear Systems Joseph Winkin Namur Center of Complex Systems (NaXys) and Department of Mathematics, University of Namur, Belgium Joint work with Bouchra Abouzaid
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationThe Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems
The Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems Mehdi Badra Laboratoire LMA, université de Pau et des Pays
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationAn Output Stabilization of Bilinear Distributed Systems
Int. Journal of Math. Analysis, Vol. 7, 213, no. 4, 195-211 An Output Stabilization of Bilinear Distributed Systems E. Zerrik, Y. Benslimane MACS Team; Sciences Faculty, Moulay Ismail University zerrik3@yahoo.fr,
More informationOn feedback stabilizability of time-delay systems in Banach spaces
On feedback stabilizability of time-delay systems in Banach spaces S. Hadd and Q.-C. Zhong q.zhong@liv.ac.uk Dept. of Electrical Eng. & Electronics The University of Liverpool United Kingdom Outline Background
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationDichotomy, the Closed Range Theorem and Optimal Control
Dichotomy, the Closed Range Theorem and Optimal Control Pavel Brunovský (joint work with Mária Holecyová) Comenius University Bratislava, Slovakia Praha 13. 5. 2016 Brunovsky Praha 13. 5. 2016 Closed Range
More informationOptimal Control of a Tubular Reactor via Cayley-Tustin Time Discretization. Peyman Tajik
Optimal Control of a Tubular Reactor via Cayley-Tustin Time Discretization by Peyman Tajik A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Chemical
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationChapter 2 Distributed Parameter Systems: Controllability, Observability, and Identification
Chapter 2 Distributed Parameter Systems: Controllability, Observability, and Identification 2.1 Mathematical Description We introduce the class of systems to be considered in the framework of this monograph
More informationStability of Linear Distributed Parameter Systems with Time-Delays
Stability of Linear Distributed Parameter Systems with Time-Delays Emilia FRIDMAN* *Electrical Engineering, Tel Aviv University, Israel joint with Yury Orlov (CICESE Research Center, Ensenada, Mexico)
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More informationOn Controllability of Linear Systems 1
On Controllability of Linear Systems 1 M.T.Nair Department of Mathematics, IIT Madras Abstract In this article we discuss some issues related to the observability and controllability of linear systems.
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationOptimal control of a linearized continuum model for re-entrant manufacturing production systems
Optimal control of a linearized continuum model for re-entrant manufacturing production systems Xiaodong Xu 1 and Stevan Dubljevic 1 1 Department of Chemical and Materials Engineering, University of Alberta,
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationInput-output finite-time stabilization for a class of hybrid systems
Input-output finite-time stabilization for a class of hybrid systems Francesco Amato 1 Gianmaria De Tommasi 2 1 Università degli Studi Magna Græcia di Catanzaro, Catanzaro, Italy, 2 Università degli Studi
More informationIEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 58, NO. 5, MAY invertible, that is (1) In this way, on, and on, system (3) becomes
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 58, NO. 5, MAY 2013 1269 Sliding Mode and Active Disturbance Rejection Control to Stabilization of One-Dimensional Anti-Stable Wave Equations Subject to Disturbance
More informationChapter 4 Optimal Control Problems in Infinite Dimensional Function Space
Chapter 4 Optimal Control Problems in Infinite Dimensional Function Space 4.1 Introduction In this chapter, we will consider optimal control problems in function space where we will restrict ourselves
More informationA Distributed Parameter Control Approach to Optimal Filtering and Smoothing with Mobile Sensor Networks
A Distributed Parameter Control Approach to Optimal Filtering and Smoothing with Mobile Sensor Networks J A Burns, Fellow, IEEE, E M Cliff, Member, IEEE, and C Rautenberg Center for Optimal Design and
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationEXISTENCE OF ASYMPTOTICALLY ALMOST AUTOMORPHIC MILD SOLUTIONS FOR NONAUTONOMOUS SEMILINEAR EVOLUTION EQUATIONS
Electronic Journal of Differential Equations, Vol. 208 208, No. 37, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF ASYMPTOTICALLY ALMOST AUTOMORPHIC MILD
More informationAdvanced Mechatronics Engineering
Advanced Mechatronics Engineering German University in Cairo 21 December, 2013 Outline Necessary conditions for optimal input Example Linear regulator problem Example Necessary conditions for optimal input
More informationTemporal-Spatial Optimal Tracking Control of Microwave Heating Process Based on the Spectral Galerkin Method
Temporal-Spatial Optimal Tracking Control of Microwave Heating Process Based on the Spectral Galerkin Method Jiaqi Zhong, Shan Liang, Yupeng Yuan Chongqing University College of Automation No 174, Shazhengjie,
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationVibrating Strings and Heat Flow
Vibrating Strings and Heat Flow Consider an infinite vibrating string Assume that the -ais is the equilibrium position of the string and that the tension in the string at rest in equilibrium is τ Let u(,
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationPredictive control of parabolic PDEs with boundary control actuation
Chemical Engineering Science 61 (26) 6239 6248 www.elsevier.com/locate/ces Shorter Communication Predictive control of parabolic PDEs with boundary control actuation Stevan Dubljevic, Panagiotis D. Christofides
More informationThe Role of Exosystems in Output Regulation
1 The Role of Exosystems in Output Regulation Lassi Paunonen In this paper we study the role of the exosystem in the theory of output regulation for linear infinite-dimensional systems. The main result
More informationNumerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps
Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps Yuri A. Kuznetsov joint work with W. Govaerts, A. Dhooge(Gent), and E. Doedel (Montreal) LCBIF
More informationSub-Optimal Boundary Control of Semilinear PDEs Using a Dyadic Perturbation Observer
Sub-Optimal Boundary Control of Semilinear PDEs Using a Dyadic Perturbation Observer Aditya A. Paranjape and Soon-Jo Chung Abstract In this paper, we present a sub-optimal controller for semilinear partial
More informationAn Integrodifferential Approach to Adaptive Control Design for Heat Transfer Systems with Uncertainties
An Integrodifferential Approach to Adaptive Control Design for Heat Transfer Systems with Uncertainties Vasily Saurin Georgy Kostin Andreas Rauh Luise Senkel Harald Aschemann Institute for Problems in
More informationMA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE
MA 201: Method of Separation of Variables Finite Vibrating String Problem ecture - 11 IBVP for Vibrating string with no external forces We consider the problem in a computational domain (x,t) [0,] [0,
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationGame Theory Extra Lecture 1 (BoB)
Game Theory 2014 Extra Lecture 1 (BoB) Differential games Tools from optimal control Dynamic programming Hamilton-Jacobi-Bellman-Isaacs equation Zerosum linear quadratic games and H control Baser/Olsder,
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationIntroduction to Semigroup Theory
Introduction to Semigroup Theory Franz X. Gmeineder LMU München, U Firenze Bruck am Ziller / Dec 15th 2012 Franz X. Gmeineder Introduction to Semigroup Theory 1/25 The Way Up: Opening The prototype of
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationLecture 4 Continuous time linear quadratic regulator
EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES
13 Kragujevac J. Math. 3 27) 13 26. ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES Boško S. Jovanović University of Belgrade, Faculty of Mathematics, Studentski trg 16, 11 Belgrade, Serbia
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationNotes for ENEE 664: Optimal Control. André L. Tits
Notes for ENEE 664: Optimal Control André L. Tits DRAFT August 2008 Contents 1 Motivation and Scope 3 1.1 Some Examples.................................. 3 1.2 Scope of the Course................................
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationExponential stabilization of a Rayleigh beam - actuator and feedback design
Exponential stabilization of a Rayleigh beam - actuator and feedback design George WEISS Department of Electrical and Electronic Engineering Imperial College London London SW7 AZ, UK G.Weiss@imperial.ac.uk
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationSolution of Stochastic Optimal Control Problems and Financial Applications
Journal of Mathematical Extension Vol. 11, No. 4, (2017), 27-44 ISSN: 1735-8299 URL: http://www.ijmex.com Solution of Stochastic Optimal Control Problems and Financial Applications 2 Mat B. Kafash 1 Faculty
More informationRisk-Sensitive and Robust Mean Field Games
Risk-Sensitive and Robust Mean Field Games Tamer Başar Coordinated Science Laboratory Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign Urbana, IL - 6181 IPAM
More informationThe Direct Transcription Method For Optimal Control. Part 2: Optimal Control
The Direct Transcription Method For Optimal Control Part 2: Optimal Control John T Betts Partner, Applied Mathematical Analysis, LLC 1 Fundamental Principle of Transcription Methods Transcription Method
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationHyperbolic PDEs. Chapter 6
Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationSyed Abbas. June 19th, 2009
Almost Department of Mathematics and Statistics Indian Institute of Technology Kanpur, Kanpur, 208016, India June 19th, 2009 Harald Bohr Function Almost Bohr s early research was mainly concerned with
More informationThe Fattorini-Hautus test
The Fattorini-Hautus test Guillaume Olive Seminar, Shandong University Jinan, March 31 217 Plan Part 1: Background on controllability Part 2: Presentation of the Fattorini-Hautus test Part 3: Controllability
More informationExplicit State and Output Feedback Boundary Controllers for Partial Differential Equations
JOURNAL OF AUTOMATIC CONTROL, UNIVERSITY OF BELGRADE; VOL. 13():1-9, 3 Explicit State and Output Feedback Boundary Controllers for Partial Differential Equations Andrey Smyshlyaev and Miroslav Krstic Abstract
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationStability, Pole Placement, Observers and Stabilization
Stability, Pole Placement, Observers and Stabilization 1 1, The Netherlands DISC Course Mathematical Models of Systems Outline 1 Stability of autonomous systems 2 The pole placement problem 3 Stabilization
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More information3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationA Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1
A Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1 Ali Jadbabaie, Claudio De Persis, and Tae-Woong Yoon 2 Department of Electrical Engineering
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationRecent result on porous medium equations with nonlocal pressure
Recent result on porous medium equations with nonlocal pressure Diana Stan Basque Center of Applied Mathematics joint work with Félix del Teso and Juan Luis Vázquez November 2016 4 th workshop on Fractional
More informationConvexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls
1 1 Convexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls B.T.Polyak Institute for Control Science, Moscow, Russia e-mail boris@ipu.rssi.ru Abstract Recently [1, 2] the new convexity
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Using Quadrature-Based Projection Methods Julian Köllermeier, Manuel Torrilhon 29th International Symposium on Rarefied Gas Dynamics, Xi an July 14th, 2014 Julian Köllermeier, Manuel Torrilhon 1 / 12 Outline
More informationTrajectory Smoothing as a Linear Optimal Control Problem
Trajectory Smoothing as a Linear Optimal Control Problem Biswadip Dey & P. S. Krishnaprasad Allerton Conference - 212 Monticello, IL, USA. October 4, 212 Background and Motivation 1 Background and Motivation
More informationNew phenomena for the null controllability of parabolic systems: Minim
New phenomena for the null controllability of parabolic systems F.Ammar Khodja, M. González-Burgos & L. de Teresa Aix-Marseille Université, CNRS, Centrale Marseille, l2m, UMR 7373, Marseille, France assia.benabdallah@univ-amu.fr
More informationSPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 43, Number 3, 2013 SPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS MILUTIN R. OSTANIĆ ABSTRACT. We establish the exact asymptotical behavior
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationMPC/LQG for infinite-dimensional systems using time-invariant linearizations
MPC/LQG for infinite-dimensional systems using time-invariant linearizations Peter Benner 1 and Sabine Hein 2 1 Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstr. 1, 39106 Magdeburg,
More informationHomework 7 Solutions
Homework 7 Solutions # (Section.4: The following functions are defined on an interval of length. Sketch the even and odd etensions of each function over the interval [, ]. (a f( =, f ( Even etension of
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationOPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN
Dynamic Systems and Applications 16 (2007) 393-406 OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN College of Mathematics and Computer
More informationDecentralized Disturbance Attenuation for Large-Scale Nonlinear Systems with Delayed State Interconnections
Decentralized Disturbance Attenuation for Large-Scale Nonlinear Systems with Delayed State Interconnections Yi Guo Abstract The problem of decentralized disturbance attenuation is considered for a new
More informationHIGHER ORDER MULTI-TERM TIME-FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS INVOLVING CAPUTO-FABRIZIO DERIVATIVE
Electronic Journal of Differential Equations, Vol. 27 27), No. 243, pp.. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu HIGHER ORDER MULTI-TERM TIME-FRACTIONAL PARTIAL DIFFERENTIAL
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationExponential stability of families of linear delay systems
Exponential stability of families of linear delay systems F. Wirth Zentrum für Technomathematik Universität Bremen 28334 Bremen, Germany fabian@math.uni-bremen.de Keywords: Abstract Stability, delay systems,
More information4. Special Solutions and Stability
4. Special Solutions and Stability Most PDEs in science and engineering are nonlinear. They cannot be solved either explicitly or semiexplicitly (separation of variables, transform methods, Green s functions,
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationLINEAR SYSTEMS OF GENERALIZED ORDINARY DIFFERENTIAL EQUATIONS
Georgian Mathematical Journal Volume 8 21), Number 4, 645 664 ON THE ξ-exponentially ASYMPTOTIC STABILITY OF LINEAR SYSTEMS OF GENERALIZED ORDINARY DIFFERENTIAL EQUATIONS M. ASHORDIA AND N. KEKELIA Abstract.
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationNumerics and Control of PDEs Lecture 7. IFCAM IISc Bangalore. Feedback stabilization of a 1D nonlinear model
1/3 Numerics and Control of PDEs Lecture 7 IFCAM IISc Bangalore July August, 13 Feedback stabilization of a 1D nonlinear model Mythily R., Praveen C., Jean-Pierre R. /3 Plan of lecture 7 1. The nonlinear
More information