skip to main content


Title: A Generalized LMI Formulation for Input-Output Analysis of Linear Systems of ODEs Coupled with PDEs
In this paper, we consider input-output properties of linear systems consisting of PDEs on a finite domain coupled with ODEs through the boundary conditions of the PDE. This framework can be used to represent e.g. a lumped mass fixed to a beam or a system with delay. This work generalizes the sufficiency proof of the KYP Lemma for ODEs to coupled ODE-PDE systems using a recently developed concept of fundamental state and the associated boundary-condition-free representation. The conditions of the generalized KYP are tested using the PQRS positive matrix parameterization of operators resulting in a finite-dimensional LMI, feasibility of which implies prima facie provable passivity or L2-gain of the system. No discretization or approximation is involved at any step and we use numerical examples to demonstrate that the bounds obtained are not conservative in any significant sense and that computational complexity is lower than existing methods involving finite-dimensional projection of PDEs.  more » « less
Award ID(s):
1935453
NSF-PAR ID:
10113773
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Proceedings of the IEEE Conference on Decision & Control
ISSN:
2576-2370
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Modern control theory provides us with a spectrum of methods for studying the interconnection of dynamic systems using input-output properties of the interconnected subsystems. Perhaps the most advanced framework for such inputoutput analysis is the use of Integral Quadratic Constraints (IQCs), which considers the interconnection of a nominal linear system with an unmodelled nonlinear or uncertain subsystem with known input-output properties. Although these methods are widely used for Ordinary Differential Equations (ODEs), there have been fewer attempts to extend IQCs to infinitedimensional systems. In this paper, we present an IQC-based framework for Partial Differential Equations (PDEs) and Delay Differential Equations (DDEs). First, we introduce infinitedimensional signal spaces, operators, and feedback interconnections. Next, in the main result, we propose a formulation of hard IQC-based input-output stability conditions, allowing for infinite-dimensional multipliers. We then show how to test hard IQC conditions with infinite-dimensional multipliers on a nominal linear PDE or DDE system via the Partial Integral Equation (PIE) state-space representation using a sufficient version of the Kalman-Yakubovich-Popov lemma (KYP). The results are then illustrated using four example problems with uncertainty and nonlinearity. 
    more » « less
  2. We present a framework for stability analysis of systems of coupled linear Partial-Differential Equations (PDEs). The class of PDE systems considered in this paper includes parabolic, elliptic and hyperbolic systems with Dirichelet, Neuman and mixed boundary conditions. The results in this paper apply to systems with a single spatial variable and assume existence and continuity of solutions except in such cases when existence and continuity can be inferred from existence of a Lyapunov function. Our approach is based on a new concept of state for PDE systems which allows us to express the derivative of the Lyapunov function as a Linear Operator Inequality directly on L2 and allows for any type of suitably well-posed boundary conditions. This approach obviates the need for integration by parts, spacing functions or similar mathematical encumbrances. The resulting algorithms are implemented in Matlab, tested on several motivating examples, and the codes have been posted online. Numerical testing indicates the approach has little or no conservatism for a large class of systems and can analyze systems of up to 20 coupled PDEs. 
    more » « less
  3. It has recently been shown that the evolution of a linear Partial Differential Equation (PDE) can be more conveniently represented in terms of the evolution of a higher spatial derivative of the state. This higher spatial derivative (termed the `fundamental state') lies in $L_2$ - requiring no auxiliary boundary conditions or continuity constraints. Such a representation (termed a Partial Integral Equation or PIE) is then defined in terms of an algebra of bounded integral operators (termed Partial Integral (PI) operators) and is constructed by identifying a unitary map from the fundamental state to the state of the original PDE. Unfortunately, when the PDE is nonlinear, the dynamics of the associated fundamental state are no longer parameterized in terms of PI operators. However, in this paper we show that such dynamics can be compactly represented using a new tensor algebra of partial integral operators acting on the tensor product of the fundamental state. We further show that this tensor product of the fundamental state forms a natural distributed equivalent of the monomial basis used in representation of polynomials on a finite-dimensional space. This new representation is then used to provide a simple SDP-based Lyapunov test of stability of quadratic PDEs. The test is applied to three illustrative examples of quadratic PDEs. 
    more » « less
  4. Abstract

    In this study we consider a coupled system of partial differential equations (PDE's) which describes a certain structural acoustics interaction. One component of this PDE system is a wave equation, which serves to model the interior acoustic wave medium within a given three dimensional chamber Ω. This acoustic wave equation is coupled on a boundary interface Γ0to a two dimensional system of thermoelasticity: this thermoelastic PDE is composed in part of a structural beam or plate equation, which governs the vibrations of flexible wall portion Γ0of the chamber Ω. Moreover, this elastic dynamics is coupled to a heat equation which also evolves on Γ0, and which imparts a thermal damping onto the entire structural acoustic system. As we said, the interaction between the wave and thermoelastic PDE components takes place on the boundary interface Γ0, and involves coupling boundary terms which are above the level of finite energy. We analyze the stability properties of this coupled structural acoustics PDE model, in the absence of any additive feedback dissipation on the hard walls Γ1of the boundary. Under a certain geometric assumption on Γ1, an assumption which has appeared in the literature in connection with structural acoustic flow, and which allows for the invocation of a recently derived microlocal boundary trace estimate, we show that classical solutions of this thermally damped structural acoustics PDE decay uniformly to zero, with a rational rate of decay.

     
    more » « less
  5. In this work, we present a Linear Matrix Inequality (LMI) based method to synthesize an optimal H1 estimator for a large class of linear coupled partial differential equations (PDEs) utilizing only finite dimensional measurements. Our approach extends the newly developed framework for representing and analyzing distributed parameter systems using operators on the space of square integrable functions that are equipped with multipliers and kernels of semi-separable class. We show that by redefining the state, the PDEs can be represented using operators that embed the boundary conditions and input-output relations explicitly. The optimal estimator synthesis problem is formulated as a convex optimization subject to LMIs that require no approximation or discretization. A scalable algorithm is presented to synthesize the estimator. The algorithm is illustrated by suitable examples. 
    more » « less