Canards from Chua's circuit

The aim of this work is to extend Beno\^it's theorem for the generic existence of "canards" solutions in singularly perturbed dynamical systems of dimension three with one fast variable to those of dimension four. Then, it is established that this result can be found according to the Flow Curvature Method. Applications to Chua's cubic model of dimension three and four enable to state the existence of "canards" solutions in such systems.Comment: arXiv admin note: text overlap with arXiv:1408.489

Landscapes of Non-gradient Dynamics Without Detailed Balance: Stable Limit Cycles and Multiple Attractors

Landscape is one of the key notions in literature on biological processes and physics of complex systems with both deterministic and stochastic dynamics. The large deviation theory (LDT) provides a possible mathematical basis for the scientists' intuition. In terms of Freidlin-Wentzell's LDT, we discuss explicitly two issues in singularly perturbed stationary diffusion processes arisen from nonlinear differential equations: (1) For a process whose corresponding ordinary differential equation has a stable limit cycle, the stationary solution exhibits a clear separation of time scales: an exponential terms and an algebraic prefactor. The large deviation rate function attains its minimum zero on the entire stable limit cycle, while the leading term of the prefactor is inversely proportional to the velocity of the non-uniform periodic oscillation on the cycle. (2) For dynamics with multiple stable fixed points and saddles, there is in general a breakdown of detailed balance among the corresponding attractors. Two landscapes, a local and a global, arise in LDT, and a Markov jumping process with cycle flux emerges in the low-noise limit. A local landscape is pertinent to the transition rates between neighboring stable fixed points; and the global landscape defines a nonequilibrium steady state. There would be nondifferentiable points in the latter for a stationary dynamics with cycle flux. LDT serving as the mathematical foundation for emergent landscapes deserves further investigations.Comment: 4 figur

Nonlinear zero-sum differential game analysis by singular perturbation methods

A class of nonlinear, zero-sum differential games, exhibiting time-scale separation properties, can be analyzed by singular-perturbation techniques. The merits of such an analysis, leading to an approximate game solution, as well as the 'well-posedness' of the formulation, are discussed. This approach is shown to be attractive for investigating pursuit-evasion problems; the original multidimensional differential game is decomposed to a 'simple pursuit' (free-stream) game and two independent (boundary-layer) optimal-control problems. Using multiple time-scale boundary-layer models results in a pair of uniformly valid zero-order composite feedback strategies. The dependence of suboptimal strategies on relative geometry and own-state measurements is demonstrated by a three dimensional, constant-speed example. For game analysis with realistic vehicle dynamics, the technique of forced singular perturbations and a variable modeling approach is proposed. Accuracy of the analysis is evaluated by comparison with the numerical solution of a time-optimal, variable-speed 'game of two cars' in the horizontal plane

Singularly Perturbed Monotone Systems and an Application to Double Phosphorylation Cycles

The theory of monotone dynamical systems has been found very useful in the modeling of some gene, protein, and signaling networks. In monotone systems, every net feedback loop is positive. On the other hand, negative feedback loops are important features of many systems, since they are required for adaptation and precision. This paper shows that, provided that these negative loops act at a comparatively fast time scale, the main dynamical property of (strongly) monotone systems, convergence to steady states, is still valid. An application is worked out to a double-phosphorylation ``futile cycle'' motif which plays a central role in eukaryotic cell signaling.Comment: 21 pages, 3 figures, corrected typos, references remove

Rapid near-optimal aerospace plane trajectory generation and guidance

Effort was directed toward the problems of the real time trajectory optimization and guidance law development for the National Aerospace Plane (NASP) applications. In particular, singular perturbation methods were used to develop guidance algorithms suitable for onboard, real time implementation. The progress made in this research effort is reported

Delayed loss of stability in singularly perturbed finite-dimensional gradient flows

In this paper we study the singular vanishing-viscosity limit of a gradient flow in a finite dimensional Hilbert space, focusing on the so-called delayed loss of stability of stationary solutions. We find a class of time-dependent energy functionals and initial conditions for which we can explicitly calculate the first discontinuity time $t^*$ of the limit. For our class of functionals, $t^*$ coincides with the blow-up time of the solutions of the linearized system around the equilibrium, and is in particular strictly greater than the time $t_c$ where strict local minimality with respect to the driving energy gets lost. Moreover, we show that, in a right neighborhood of $t^*$, rescaled solutions of the singularly perturbed problem converge to heteroclinic solutions of the gradient flow. Our results complement the previous ones by Zanini, where the situation we consider was excluded by assuming the so-called transversality conditions, and the limit evolution consisted of strict local minimizers of the energy up to a negligible set of times

Rapid near-optimal trajectory generation and guidance law development for single-stage-to-orbit airbreathing vehicles

General problems associated with on-board trajectory optimization, propulsion system cycle selection, and with the synthesis of guidance laws were addressed for an ascent to low-earth-orbit of an air-breathing single-stage-to-orbit vehicle. The NASA Generic Hypersonic Aerodynamic Model Example and the Langley Accelerator aerodynamic sets were acquired and implemented. Work related to the development of purely analytic aerodynamic models was also performed at a low level. A generic model of a multi-mode propulsion system was developed that includes turbojet, ramjet, scramjet, and rocket engine cycles. Provisions were made in the dynamic model for a component of thrust normal to the flight path. Computational results, which characterize the nonlinear sensitivity of scramjet performance to changes in vehicle angle of attack, were obtained and incorporated into the engine model. Additional trajectory constraints were introduced: maximum dynamic pressure; maximum aerodynamic heating rate per unit area; angle of attack and lift limits; and limits on acceleration both along and normal to the flight path. The remainder of the effort focused on required modifications to a previously derived algorithm when the model complexity cited above was added. In particular, analytic switching conditions were derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another for two cases: the case in which engine cycle operations can overlap, and the case in which engine cycle operations are mutually exclusive. The resulting guidance algorithm was implemented in software and exercised extensively. It was found that the approximations associated with the assumed time scale separation employed in this work are reasonable except over the Mach range from roughly 5 to 8. This phenomenon is due to the very large thrust capability of scramjets in this Mach regime when sized to meet the requirement for ascent to orbit. By accounting for flight path angle and flight path angle rate in construction of the flight path over this Mach range, the resulting algorithm provides the means for rapid near-optimal trajectory generation and propulsion cycle selection over the entire Mach range from take-off to orbit

Adaptive Pseudo-Transient-Continuation-Galerkin Methods for Semilinear Elliptic Partial Differential Equations

In this paper we investigate the application of pseudo-transient-continuation (PTC) schemes for the numerical solution of semilinear elliptic partial differential equations, with possible singular perturbations. We will outline a residual reduction analysis within the framework of general Hilbert spaces, and, subsequently, employ the PTC-methodology in the context of finite element discretizations of semilinear boundary value problems. Our approach combines both a prediction-type PTC-method (for infinite dimensional problems) and an adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully adaptive PTC-Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for different examples.Comment: arXiv admin note: text overlap with arXiv:1408.522