Constraint reasoning for differential models

The basic motivation of this work was the integration of biophysical models within the interval constraints framework for decision support. Comparing the major features of biophysical models with the expressive power of the existing interval constraints framework, it was clear that the most important inadequacy was related with the representation of differential equations. System dynamics is often modelled through differential equations but there was no way of expressing a differential equation as a constraint and integrate it within the constraints framework. Consequently, the goal of this work is focussed on the integration of ordinary differential equations within the interval constraints framework, which for this purpose is extended with the new formalism of Constraint Satisfaction Differential Problems. Such framework allows the specification of ordinary differential equations, together with related information, by means of constraints, and provides efficient propagation techniques for pruning the domains of their variables. This enabled the integration of all such information in a single constraint whose variables may subsequently be used in other constraints of the model. The specific method used for pruning its variable domains can then be combined with the pruning methods associated with the other constraints in an overall propagation algorithm for reducing the bounds of all model variables. The application of the constraint propagation algorithm for pruning the variable domains, that is, the enforcement of local-consistency, turned out to be insufficient to support decision in practical problems that include differential equations. The domain pruning achieved is not, in general, sufficient to allow safe decisions and the main reason derives from the non-linearity of the differential equations. Consequently, a complementary goal of this work proposes a new strong consistency criterion, Global Hull-consistency, particularly suited to decision support with differential models, by presenting an adequate trade-of between domain pruning and computational effort. Several alternative algorithms are proposed for enforcing Global Hull-consistency and, due to their complexity, an effort was made to provide implementations able to supply any-time pruning results. Since the consistency criterion is dependent on the existence of canonical solutions, it is proposed a local search approach that can be integrated with constraint propagation in continuous domains and, in particular, with the enforcing algorithms for anticipating the finding of canonical solutions. The last goal of this work is the validation of the approach as an important contribution for the integration of biophysical models within decision support. Consequently, a prototype application that integrated all the proposed extensions to the interval constraints framework is developed and used for solving problems in different biophysical domains

An Interval Constraint Programming Approach for Quasi Capture Tube Validation

Proving that the state of a controlled nonlinear system always stays inside a time moving bubble (or capture tube) amounts to proving the inconsistency of a set of nonlinear inequalities in the time-state space. In practice however, even with a good intuition, it is difficult for a human to find such a capture tube except for simple examples. In 2014, Jaulin et al. established properties that support a new interval approach for validating a quasi capture tube, i.e. a candidate tube (with a simple form) from which the mobile system can escape, but into which it enters again before a given time. A quasi capture tube is easy to find in practice for a controlled system. Merging the trajectories originated from the candidate tube yields the smallest capture tube enclosing it. This paper proposes an interval constraint programming solver dedicated to the quasi capture tube validation. The problem is viewed as a differential CSP where the functional variables correspond to the state variables of the system and the constraints define system trajectories that escape from the candidate tube "for ever". The solver performs a branch and contract procedure for computing the trajectories that escape from the candidate tube. If no solution is found, the quasi capture tube is validated and, as a side effect, a corrected smallest capture tube enclosing the quasi one is computed. The approach is experimentally validated on several examples having 2 to 5 degrees of freedom

Advances in Reachability Analysis for Nonlinear Dynamic Systems

Systems of nonlinear ordinary differential equations (ODEs) are used to model an incredible variety of dynamic phenomena in chemical, oil and gas, and pharmaceutical industries. In reality, such models are nearly always subject to significant uncertainties in their initial conditions, parameters, and inputs. This dissertation provides new theoretical and numerical techniques for rigorously enclosing the set of solutions reachable by a given systems of nonlinear ODEs subject to uncertain initial conditions, parameters, and time-varying inputs. Such sets are often referred to as reachable sets, and methods for enclosing them are critical for designing systems that are passively robust to uncertainty, as well as for optimal real-time decision-making. Such enclosure methods are used extensively for uncertainty propagation, robust control, system verification, and optimization of dynamic systems arising in a wide variety of applications. Unfortunately, existing methods for computing such enclosures often provide an unworkable compromise between cost and accuracy. For example, interval methods based on differential inequalities (DI) can produce bounds very efficiently but are often too conservative to be of any practical use. In contrast, methods based on more complex sets can achieve sharp bounds, but are far too expensive for real-time decision-making and scale poorly with problem size. Recently, it has been shown that bounds computed via differential inequalities can often be made much less conservative while maintaining high efficiency by exploiting redundant model equations that are known to hold for all trajectories of interest (e.g., linear relationships among chemical species in a reaction network that hold due to the conservation of mass or elements). These linear relationships are implied by the governing ODEs, and can thus be considered redundant. However, these advances are only applicable to a limited class of system in which pre-existing linear redundant model equations are available. Moreover, the theoretical results underlying these algorithms do not apply to redundant equations that depend on time-varying inputs and rely on assumptions that prove to be very restrictive for nonlinear redundant equations, etc. This dissertation continues a line of research that has recently achieved very promising bounding results using methods based on differential inequalities. In brief, the major contributions can be divided into three categories: (1) In regard to algorithms, this dissertation significantly improves existing algorithms that exploit linear redundant model equations to achieve more accurate and efficient enclosures. It also develops new fast and accurate bounding algorithms that can exploit nonlinear redundant model equations. (2) Considering theoretical contributions, it develops a novel theoretical framework for the introduction of redundant model equations into arbitrary dynamic models to effectively reduce conservatism. The newly developed theories have more generality in terms of application. For example, complex nonlinear constraints that involve states, time derivatives of the system states, and time- varying inputs are allowed to be exploited. (3) A new differential inequalities method called Mean Value Differential Inequalities (MVDI) is developed that can automatically introduce redundant model equations for arbitrary dynamic systems and has a second-order convergence rate reported the first time among DI-based methods

Evaluation of the effectiveness of the interval computation method to simulate the dynamic behavior of subdefinite system: application on an active suspension system

International audienceA new design approach based on methods by intervals adapted to the integration of the simulation step at the earliest stage of preliminary design for dynamic systems is proposed in this study. The main idea consists on using the interval computation method to make a simulation by intervals in order to minimize the number of simulations which allow obtaining a set of solutions instead of a single one. These intervals represent the domains of possible values for the design parameters of the subdefinite system. So the parameterized model of the system is solved by interval. This avoids launching n simulations with n values for each design parameter. The proposed method is evaluated by several tests on a scalable numerical example. It has been applied to solve parameterized differential equations of a Macpher-son suspension system and to study its dynamic behavior in its passive and active form. The dynamic model of the active suspension is nonlinear but linearisable. It is transformed into a parameterized state equation by intervals. The solution to this state equation is given in the form of a matrix exponential. Three digital implementations of exponential have been tested to obtain convergent results. Simulations results are presented and discussed

Global optimisation for dynamic systems using novel overestimation reduction techniques

The optimisation of dynamic systems is of high relevance in chemical engineering as many practical systems can be described by ordinary differential equations (ODEs) or differential algebraic equations (DAEs). The current techniques for solving these problems rigorously to global optimality rely mainly on sequential approaches in which a branch and bound framework is used for solving the global optimisation part of the problem and a verified simulator (in which rounding errors are accounted for in the computations) is used for solving the dynamic constraints. The verified simulation part is the main bottleneck since tight bounds are difficult to obtain for high dimensional dynamic systems. Additionally, uncertainty in the form of, for example, intervals is introduced in the parameters of the dynamic constraints which are also the decision variables of the optimisation problem. Nevertheless, in the verified simulation the accumulation of trajectories that do not belong to the exact solution (overestimation) makes the state bounds overconservative and in the worst case they blow up and tend towards ±∞. In this thesis, methods for verified simulation in global optimisation for dynamic systems were investigated. A novel algorithm that uses an interval Taylor series (ITS) method with enhanced overestimation reduction capabilities was developed. These enhancements for the reduction of the overestimation rely on interval contractors (Krawczyk, Newton, ForwardBackward) and model reformulation based on pattern substitution and input scaling. The method with interval contractors was also extended to Taylor Models (TM) for comparison purposes. The two algorithms were tested on several case studies to demonstrate the effectiveness of the methods. The case studies have a different number of state variables and system parameters and they use uncertain amounts in some of the system parameters and initial conditions. Both of the methods were also used in a sequential approach to address the global optimisation for dynamic systems problem subject to uncertainty. The simulation results demonstrated that the ITS method with overestimation reduction techniques provided tighter state bounds with less computational expense than the traditional method. In the case of the forward-backward contractor additional constraints can be introduced that can potentially contribute significantly to the reduction of the overestimation. Similarly, the novel TM method with enhanced overestimation reduction capabilities provided tighter bounds than the TM method alone. On the other hand, the optimisation results showed that the global optimisation algorithm with the novel ITS method with overestimation reduction techniques converged faster to a rigorous solution due to the improved state bounds

Chebyshev model arithmetic for factorable functions

This article presents an arithmetic for the computation of Chebyshev models for factorable functions and an analysis of their convergence properties. Similar to Taylor models, Chebyshev models consist of a pair of a multivariate polynomial approximating the factorable function and an interval remainder term bounding the actual gap with this polynomial approximant. Propagation rules and local convergence bounds are established for the addition, multiplication and composition operations with Chebyshev models. The global convergence of this arithmetic as the polynomial expansion order increases is also discussed. A generic implementation of Chebyshev model arithmetic is available in the library MC++. It is shown through several numerical case studies that Chebyshev models provide tighter bounds than their Taylor model counterparts, but this comes at the price of extra computational burden

Model Checking Delay Differential Equations Against Metric Interval Temporal Logic

Delay differential equations (DDEs) play an important role in the modeling of dynamic processes. Delays arise in contemporary control schemes like networked distributed control and can cause deterioration of control performance, invalidating both stability and safety properties. This induces an interest in DDE especially in the area of modeling and verification of embedded control. In this article, we present an approach aiming at automatic safety verification of a simple class of DDEs against requirements expressed in a linear-time temporal logic. As requirements specification language, we exploit metric interval temporal logic (MITL) with a continuous-time semantics evaluating signals over metric spaces. We employ an over-approximation method based on interval Taylor series to enclose the solution of the DDE and thereby reduce the continuous-time verification problem for MITL formulae to a discrete-time problem over sequences of Taylor coefficients. We encode sufficient conditions for satisfaction as SMT formulae over polynomial arithmetic and use the iSAT3 SMT solver in its bounded model-checking mode for discharging the resulting proof obligations, thus proving satisfaction of time-bounded MITL specifications by the trajectories induced by a DDE. In contrast to our preliminary work in [44], we can verify arbitrary time-bounded MITL formulae, including nesting of modalities, rather than just invariance properties