The article “Mathematical Foundations of Physical Systems Modeling Languages”, co-authored by Albert Benveniste, Benoît Caillaud and Mathias Malandain, has been accepted by the IFAC journal Annual Reviews in Control. It can be read and downloaded until January 20th, 2021 for free by following this link:
The author version was made available in August 2020.
Modern modeling languages for general physical systems, such as Modelica, Amesim, or Simscape, rely on Differential Algebraic Equations (DAEs), i.e., constraints of the form ? (?′ , ?, ?) = 0. This drastically facilitates modeling from first principles of the physics, as well as the reuse of models. In this paper, we develop the mathematical theory needed to establish the development of compilers and tools for DAE-based physical modeling languages on solid mathematical bases.
Unlike Ordinary Differential Equations (ODEs, of the form ?′ = ?(?, ?)), DAEs exhibit subtle issues because of the notion of differentiation index and related latent equations—ODEs are DAEs of index zero, for which no latent equation needs to be considered. Prior to generating execution code and calling solvers, the compilation of such languages requires a nontrivial structural analysis step that reduces the differentiation index to a level acceptable by DAE solvers.
The models supported by tools of the Modelica class involve multiple modes, with mode-dependent DAE- based dynamics and state-dependent mode switching. However, multimode DAEs are much more difficult to handle than DAEs, especially because of the events of mode change. Unfortunately, the large literature devoted to the mathematical analysis of DAEs does not cover the multimode case, typically saying nothing about mode changes. This lack of foundations causes numerous difficulties to the existing modeling tools. Some models are well handled, others are not, with no clear boundary between the two classes.
In this paper, we develop a comprehensive mathematical approach supporting compilation and code generation for this class of languages. Its core is the structural analysis of multimode DAE systems. As a byproduct of this structural analysis, we propose sound criteria for accepting or rejecting multimode models. Our mathematical development relies on nonstandard analysis, which allows us to cast hybrid system dynamics to discrete-time dynamics with infinitesimal step size, thus providing a uniform framework for handling both continuous dynamics and mode change events.