 Award ID(s):
 1935453
 Publication Date:
 NSFPAR ID:
 10113773
 Journal Name:
 Proceedings of the IEEE Conference on Decision & Control
 ISSN:
 25762370
 Sponsoring Org:
 National Science Foundation
More Like this

We present a framework for stability analysis of systems of coupled linear PartialDifferential 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 wellposed 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.

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 semiseparable class. We show that by redefining the state, the PDEs can be represented using operators that embed the boundary conditions and inputoutput 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.

Abstract This paper examines a class of involutionconstrained PDEs where some part of the PDE system evolves a vector field whose curl remains zero or grows in proportion to specified source terms. Such PDEs are referred to as curlfree or curlpreserving, respectively. They arise very frequently in equations for hyperelasticity and compressible multiphase flow, in certain formulations of general relativity and in the numerical solution of Schrödinger’s equation. Experience has shown that if nothing special is done to account for the curlpreserving vector field, it can blow up in a finite amount of simulation time. In this paper, we catalogue a class of DGlike schemes for such PDEs. To retain the globally curlfree or curlpreserving constraints, the components of the vector field, as well as their higher moments, must be collocated at the edges of the mesh. They are updated using potentials collocated at the vertices of the mesh. The resulting schemes: (i) do not blow up even after very long integration times, (ii) do not need any special cleaning treatment, (iii) can operate with large explicit timesteps, (iv) do not require the solution of an elliptic system and (v) can be extended to higher orders using DGlike methods. Themore »

Abstract Several important PDE systems, like magnetohydrodynamics and computational electrodynamics, are known to support involutions where the divergence of a vector field evolves in divergencefree or divergence constraintpreserving fashion. Recently, new classes of PDE systems have emerged for hyperelasticity, compressible multiphase flows, socalled firstorder reductions of the Einstein field equations, or a novel firstorder hyperbolic reformulation of Schrödinger’s equation, to name a few, where the involution in the PDE supports curlfree or curl constraintpreserving evolution of a vector field. We study the problem of curl constraintpreserving reconstruction as it pertains to the design of mimetic finite volume (FV) WENOlike schemes for PDEs that support a curlpreserving involution. (Some insights into discontinuous Galerkin (DG) schemes are also drawn, though that is not the prime focus of this paper.) This is done for two and threedimensional structured mesh problems where we deliver closed form expressions for the reconstruction. The importance of multidimensional Riemann solvers in facilitating the design of such schemes is also documented. In two dimensions, a von Neumann analysis of structurepreserving WENOlike schemes that mimetically satisfy the curl constraints, is also presented. It shows the tremendous value of higher order WENOlike schemes in minimizing dissipation and dispersion for this classmore »

Embedding properties of network realizations of dissipative reduced order models Jörn Zimmerling, Mikhail Zaslavsky,Rob Remis, Shasri Moskow, Alexander Mamonov, Murthy Guddati, Vladimir Druskin, and Liliana Borcea Mathematical Sciences Department, Worcester Polytechnic Institute https://www.wpi.edu/people/vdruskin Abstract Realizations of reduced order models of passive SISO or MIMO LTI problems can be transformed to tridiagonal and blocktridiagonal forms, respectively, via dierent modications of the Lanczos algorithm. Generally, such realizations can be interpreted as ladder resistorcapacitorinductor (RCL) networks. They gave rise to network syntheses in the rst half of the 20th century that was at the base of modern electronics design and consecutively to MOR that tremendously impacted many areas of engineering (electrical, mechanical, aerospace, etc.) by enabling ecient compression of the underlining dynamical systems. In his seminal 1950s works Krein realized that in addition to their compressing properties, network realizations can be used to embed the data back into the state space of the underlying continuum problems. In more recent works of the authors Krein's ideas gave rise to socalled nitedierence Gaussian quadrature rules (FDGQR), allowing to approximately map the ROM statespace representation to its full order continuum counterpart on a judicially chosen grid. Thus, the state variables can be accessed directly from themore »