Two-parameter Sturm-Liouville problems

This paper deals with the computation of the eigenvalues of two-parameter Sturm- Liouville (SL) problems using the Regularized Sampling Method, a method which has been effective in computing the eigenvalues of broad classes of SL problems (Singular, Non-Self-Adjoint, Non-Local, Impulsive,...). We have shown, in this work that it can tackle two-parameter SL problems with equal ease. An example was provided to illustrate the effectiveness of the method.Comment: 9 page

Differential-Algebraic Equations and Beyond: From Smooth to Nonsmooth Constrained Dynamical Systems

The present article presents a summarizing view at differential-algebraic equations (DAEs) and analyzes how new application fields and corresponding mathematical models lead to innovations both in theory and in numerical analysis for this problem class. Recent numerical methods for nonsmooth dynamical systems subject to unilateral contact and friction illustrate the topicality of this development.Comment: Preprint of Book Chapte

Startup process in the Richtmyer-Meshkov instability

An analytical model for the initial growth period of the planar Richtmyer–Meshkov instability is presented for the case of a reflected shock, which corresponds in general to light-to-heavy interactions. The model captures the main features of the interfacial perturbation growth before the regime with linear growth in time is attained. The analysis provides a characteristic time scale τ for the startup phase of the instability, expressed explicitly as a function of the perturbation wavenumber k, the algebraic transmitted and reflected shock speeds U_(S1) 0 (defined in the frame of the accelerated interface), and the postshock Atwood number A^+: τ=[(1-A^+)/U_(S2)+(1+A^+)/(-U_(s1))]/(2k). Results are compared with computations obtained from two-dimensional highly resolved numerical simulations over a wide range of incident shock strengths S and preshock Atwood ratios A. An interesting observation shows that, within this model, the amplitude of small perturbations across a light-to-heavy interface evolves quadratically in time (and not linearly) in the limit A→1^−

A Numerical Slow Manifold Approach to Model Reduction for Optimal Control of Multiple Time Scale ODE

Time scale separation is a natural property of many control systems that can be ex- ploited, theoretically and numerically. We present a numerical scheme to solve optimal control problems with considerable time scale separation that is based on a model reduction approach that does not need the system to be explicitly stated in singularly perturbed form. We present examples that highlight the advantages and disadvantages of the method

Mixed finite difference method for singularly perturbed differential difference equations with mixed shifts via domain decomposition

AbstractIn this paper, a mixed finite difference method is proposed to solve singularly perturbed differential difference equations with mixed shifts, solutions of which exhibit boundary layer behaviour at the left end of the interval using domain decomposition. A terminal boundary point is introduced into the domain, to decompose it into inner and outer regions. The original problem is reduced to an asymptotically equivalent singular perturbation problem and with the terminal point the singular perturbation problem is treated as inner region and outer region problems separately. The outer region and the modified inner region problems are solved by mixed finite difference method. The method is repeated for various choices of the terminal point. To validate the computational efficiency of the method model examples have been solved for different values of perturbation, delay and advanced parameters. Convergence of the proposed scheme has also been investigated

Fifty Years of Stiffness

The notion of stiffness, which originated in several applications of a different nature, has dominated the activities related to the numerical treatment of differential problems for the last fifty years. Contrary to what usually happens in Mathematics, its definition has been, for a long time, not formally precise (actually, there are too many of them). Again, the needs of applications, especially those arising in the construction of robust and general purpose codes, require nowadays a formally precise definition. In this paper, we review the evolution of such a notion and we also provide a precise definition which encompasses all the previous ones.Comment: 24 pages, 11 figure

A neighboring extremal solution for an optimal switched impulsive control problem

This paper presents a neighboring extremal solution for a class of optimal switched impulsive control problems with perturbations in the initial state, terminal condition and system's parameters. The sequence of mode's switching is pre-specified, and the decision variables, i.e. the switching times and parameters of the system involved, have inequality constraints. It is assumed that the active status of these constraints is unchanged with the perturbations. We derive this solution by expanding the necessary conditions for optimality to first-order and then solving the resulting multiple-point boundary-value problem by the backward sweep technique. Numerical simulations are presented to illustrate this solution method