No Arabic abstract
Motivated by engineering applications of subsea installation by deepwater construction vessels in oil drilling, and of aid delivery by unmanned aerial vehicles in disaster relief, we develop output-feedback boundary control of heterodirectional coupled hyperbolic PDEs sandwiched between two ODEs, where the measurement is the output state of one ODE and suffers a time delay. After rewriting the time-delay dynamics as a transport PDE of which the left boundary connects with the sandwiched system, a state observer is built to estimate the states of the overall system of ODE-heterodirectional coupled hyperbolic PDEs-ODE-transport PDE using the right boundary state of the last transport PDE. An observer-based output-feedback controller acting at the first ODE is designed to stabilize the overall system using backstepping transformations and frequency-domain designs. The exponential stability results of the closed-loop system, boundedness and exponential convergence of the control input are proved. The obtained theoretical result is applied to control of a deepwater oil drilling construction vessel as a simulation case, where the simulation results show the proposed control design reduces cable oscillations and places the oil drilling equipment to be installed in the target area on the sea floor. Performance deterioration under extreme and unmodeled disturbances is also illustrated.
We analyze a class of cell-bulk coupled PDE-ODE models, motivated by quorum and diffusion sensing phenomena in microbial systems, that characterize communication between localized spatially segregated dynamically active signaling compartments that have a permeable boundary. Each cell secretes a signaling chemical into the bulk region at a constant rate and receives a feedback of the bulk chemical from the entire collection of cells. This global feedback, which activates signaling pathways within the cells, modifies the intracellular dynamics according to the external environment. The cell secretion and global feedback are regulated by permeability parameters across the cell membrane. For arbitrary reaction-kinetics within each cell, the method of matched asymptotic expansions is used in the limit of small cell radius to construct steady-state solutions of the PDE-ODE model, and to derive a globally coupled nonlinear matrix eigenvalue problem (GCEP) that characterizes the linear stability properties of the steady-states. In the limit of large bulk diffusivity an asymptotic analysis of the PDE-ODE model leads to a limiting ODE system for the spatial average of the concentration in the bulk region that is coupled to the intracellular dynamics within the cells. Results from the linear stability theory and ODE dynamics are illustrated for Selkov reaction-kinetics, where the kinetic parameters are chosen so that each cell is quiescent when uncoupled from the bulk medium. For various specific spatial configurations of cells, the linear stability theory is used to construct phase diagrams in parameter space characterizing where a switch-like emergence of intracellular oscillations can occur through a Hopf bifurcation.
When mathematical/computational problems reach infinity, extending analysis and/or numerical computation beyond it becomes a notorious challenge. We suggest that, upon suitable singular transformations (that can in principle be computationally detected on the fly) it becomes possible to go beyond infinity to the other side, with the solution becoming again well behaved and the computations continuing normally. In our lumped, Ordinary Differential Equation (ODE) examples this infinity crossing can happen instantaneously; at the spatially distributed, Partial Differential Equation (PDE) level the crossing of infinity may even persist for finite time, necessitating the introduction of conceptual (and computational) buffer zones in which an appropriate singular transformation is continuously (locally) detected and performed. These observations (and associated tools) could set the stage for a systematic approach to bypassing infinity (and thus going beyond it) in a broader range of evolution equations; they also hold the promise of meaningfully and seamlessly performing the relevant computations. Along the path of our analysis, we present a regularization process via complexification and explore its impact on the dynamics; we also discuss a set of compactification transformations and their intuitive implications.
We present a new paradigm for Neural ODE algorithms, called ODEtoODE, where time-dependent parameters of the main flow evolve according to a matrix flow on the orthogonal group O(d). This nested system of two flows, where the parameter-flow is constrained to lie on the compact manifold, provides stability and effectiveness of training and provably solves the gradient vanishing-explosion problem which is intrinsically related to training deep neural network architectures such as Neural ODEs. Consequently, it leads to better downstream models, as we show on the example of training reinforcement learning policies with evolution strategies, and in the supervised learning setting, by comparing with previous SOTA baselines. We provide strong convergence results for our proposed mechanism that are independent of the depth of the network, supporting our empirical studies. Our results show an intriguing connection between the theory of deep neural networks and the field of matrix flows on compact manifolds.
We present designs for exponential stabilization of an ODE-heat PDE-ODE coupled system where the control actuation only acts in one ODE. The combination of PDE backstepping and ODE backstepping is employed in a state-feedback control law and in an observer that estimates PDE and two ODE states only using one PDE boundary measurement. Based on the state-feedback control law and the observer, the output-feedback control law is then proposed. The exponential stability of the closed-loop system and the boundedness and exponential convergence of the control law are proved via Lyapunov analysis. Finally, numerical simulations validate the effectiveness of this method for the `sandwiched system.
Quadratization problem is, given a system of ODEs with polynomial right-hand side, transform the system to a system with quadratic right-hand side by introducing new variables. Such transformations have been used, for example, as a preprocessing step by model order reduction methods and for transforming chemical reaction networks. We present an algorithm that, given a system of polynomial ODEs, finds a transformation into a quadratic ODE system by introducing new variables which are monomials in the original variables. The algorithm is guaranteed to produce an optimal transformation of this form (that is, the number of new variables is as small as possible), and it is the first algorithm with such a guarantee we are aware of. Its performance compares favorably with the existing software, and it is capable to tackle problems that were out of reach before.