Singularity analysis of a variant of the Painlevé–Ince equation

We examine by singularity analysis an equation derived by reduction using Lie point symmetries from the Euler–Bernoulli Beam equation which is the Painlevé–Ince Equation with additional terms. The equation possesses the same leading-order behaviour and resonances as the Painlevé–Ince Equation and has a Right Painlevé Series. However, it has no Left Painlevé Series. A conjecture for the existence of Left Painlevé Series for ordinary differential equations is given.     

The Interaction of Shocks with Dispersive Waves II. Incompressible-Integrable Limit

This is the second in a two-part series of articles in which we analyze a system similar in structure to the well-known Zakharov equations from weak plasma turbulence theory, but with a nonlinear conservation equation allowing finite time shock formation. In this article we analyze the incompressible limit in which the shock speed is large compared to the underlying group velocity of the dispersive wave (a situation typically encountered in applications). After presenting some exact solutions of the full system, a multiscale perturbation method is used to resolve several basic wave interactions. The analysis breaks down into two categories: the nonlinear limit and the linear limit, corresponding to the form of the equations when the group velocity to shock speed ratio, denoted by ε, is zero. The former case is an integrable limit in which the model reduces to the cubic nonlinear Schrödinger equation governing the dispersive wave envelope. We focus on the interaction of a “fast” shock wave and a single hump soliton. In the latter case, the ε=0 problem reduces to the linear Schrödinger equation, and the focus is on a fast shock interacting with a dispersive wave whose amplitude is cusped and exponentially decaying. To motivate the time scales and structure of the shock-dispersive wave interactions at lowest orders, we first analyze a simpler system of ordinary differential equations structurally similar to the original system. Then we return to the fully coupled partial differential equations and develop a multiscale asymptotic method to derive the effective leading-order shock equations and the leading-order modulation equations governing the phase and amplitude of the dispersive wave envelope. The leading-order interaction equations admit a fairly complete analysis based on characteristic methods. Conditions are derived in which: (a) the shock passes through the soliton, (b) the shock is completely blocked by the soliton, or (c) the shock reverses direction. In the linear limit, a phenomenon is described in which the dispersive wave induces the formation of a second, transient shock front in the rapidly moving hyperbolic wave. In all cases, we can characterize the long-time dynamics of the shock. The influence of the shock on the dispersive wave is manifested, to leading order, in the generalized frequency of the dispersive wave: the fast-time part of the frequency is the shock wave itself. Hence, the frequency undergoes a sudden jump across the shock layer.In the last section, a sequence of numerical experiments depicting some of the interesting interactions predicted by the analysis is performed on the leading-order shock equations.     

Asymptotic analysis of compound liquid jets at low Reynolds numbers

Asymptotic methods based on the slenderness ratio are used to obtain the leading-order equations which govern the fluid dynamics of both hollow and solid, compound jets such as those employed in the manufacture of textile fibers, composite fibers and optical fibers. These fibers consist of an inner material which may be a round jet or an annular one and which, in turn, is surrounded by an annular jet in contact with ambient air. It is shown that the leading-order equations are one-dimensional and that it is possible to obtain analytical solutions for several flow regimes for steady, compound jets. These analytical solutions are presented here. The one-dimensional leading-order equations for the fluid dynamics of annular liquid jets at low Reynolds numbers are also derived here.     

Invariant solutions as internal singularities of nonlinear differential equations and their use for qualitative analysis of implicit and numerical solutions

Lie group analysis of nonlinear differential equations reveals existence of singularities provided by invariant solutions and invisible from the form of the equation in question. We call them internal singularities in contrast with external singularities manifested by the form of the equation. It is illustrated by way of examples that internal singularities are useful for analyzing a behaviour of solutions of nonlinear differential equations near external singularities.     

Asymptotic formulae of Liouville–Green type for higher even-order equations

Asymptotic formulae of Liouville–Green type for general linear ordinary differential equations of an arbitrary even-order 2m are investigated. A theorem on asymptotic behaviour at the infinity of 2m linearly independent solutions is proved. It is shown that numerous results known in the literature are contained in this theorem as particular cases.     

Modeling and simulation of nitrogen regulation in Corynebacterium glutamicum

Modeling the dynamic behaviour of biochemical systems at a molecular level aims at understanding and predicting the interactions of macromolecules inside the cell. Models of small subsystems based on differential equations not only prepare the way for the long-term goal of understanding a whole cell, but are inherently valuable due to their ability to predict the behaviour of the subsystem for varying external conditions or parameters. Nitrogen supply is essential for prokaryotes, thus the nitrogen uptake is an interesting target for model building. The goal is to provide new information about the interactions of the relevant proteins by performing various simulations.A model based on piecewise linear differential equations is formulated for the nitrogen uptake in Corynebacterium glutamicum. We theoretically derive a model for biochemical networks and introduce a general method for the parameter estimation which is also applicable in the case of very short time series. This approach is applied to a special system concerning the nitrogen uptake using Western blot experiments. The equations are set up for the main components of this system, the optimization problem for parameter estimation is formulated and solved, and simulations for the evaluation of the model as well as for predictions are carried out.We show that model building based on differential equations can also, when only a few measurements are performed, lead to a satisfactory model which provides valuable insights into the way it’s network components function. For example, we are able to make predictions about the maximal value of the time course as well as the steady-state level of the signal transduction protein GlnK in case of restricted activity of the proteases when considering the transition of nitrogen starvation to nitrogen excess or vice versa.     

Asymptotic stability of impulsive stochastic partial differential equations with infinite delays

In this paper, we study the existence and asymptotic stability in pth moment of mild solutions to nonlinear impulsive stochastic partial differential equations with infinite delay. By employing a fixed point approach, sufficient conditions are derived for achieving the required result. These conditions do not require the monotone decreasing behaviour of the delays.     

Stability and error analysis of one-leg methods for nonlinear delay differential equations

This paper is concerned with the numerical solution of delay differential equations (DDEs). We focus on the stability behaviour and error analysis of one-leg methods with respect to nonlinear DDEs. The new concepts of GR-stability, GAR-stability and weak GAR-stability are introduced. It is proved that a strongly A-stable one-leg method with linear interpolation is GAR-stable, and that an A-stable one-leg method with linear interpolation is GR-stable, weakly GAR-stable and D-convergent of order s, if it is consistent of order s in the classical sense.     

Verticum-type systems applied to ecological monitoring

In the paper ecological interaction chains of the type resource - producer - primary user - secondary consumer are considered. The dynamic behaviour of these four-level chains is modelled by a system of differential equations, the linearization of which is a verticum-type system introduced for the study of industrial verticums. Applying the technique of such systems, for the monitoring of the considered ecological system, an observer system is constructed, which makes it possible to recover the whole state process from the partial observation of the ecological interaction chain.     

Ostrowski’s fourth-order iterative method speedily solves cubic equations of state

Pressure-volume-temperature (P-V-T) data are required in simulating chemical plants because the latter usually involve production, separation, transportation, and storage of fluids. In the absence of actual experimental data, the pertinent mathematical model must rely on phase behaviour prediction by the so-called equations of state (EOS). When the plant model is a combination of differential and algebraic equations, simulation generally relies on numerical integration which proceeds in a piecewise fashion unless an approximate solution is needed at a single point. Needless to say, the constituent algebraic equations must be efficiently re-solved before each update of derivatives. Now, Ostrowski’s fourth-order iterative technique is a partial substitution variant of Newton’s popular second-order method. Although simple and powerful, this two-point variant has been utilised very little since its publication over forty years ago. After a brief introduction to cubic equations of state and their solution, this paper solves five of them. The results clearly demonstrate the superiority of Ostrowski’s method over Newton’s, Halley’s, and Chebyshev’s solvers.     

Invariant tensors and partial differential equations

We consider tensors with coefficients in a commutative differential algebra A. Using the Lie derivative, we introduce the notion of a tensor invariant under a derivation on an ideal of A. Each system of partial differential equations generates an ideal in some differential algebra. This makes it possible to study invariant tensors on such an ideal. As examples we consider the equations of gas dynamics and magnetohydrodynamics.     

A change of variables in the asymptotic theory of differential equations with unbounded delays

In this paper, the author investigates the asymptotic properties of solutions of the nonhomogeneous linear differential equation$\text{x}\text{̇}\text{(t)=ax(τ(t))+bx(t)+f(t)}$with nonzero real scalars a,b and the unbounded lag. Using the change of the independent and dependent variable he relates the asymptotic behaviour of solutions of this equation to the asymptotic behaviour of solutions of auxiliary functional (nondifferential) equations.     

Exponential stability of impulsive stochastic functional differential equations

In this paper, we investigate the pth moment and almost sure exponential stability of impulsive stochastic functional differential equations with finite delay by using Lyapunov method. Several stability theorems of impulsive stochastic functional differential equations with finite delay are derived. These new results are employed to impulsive stochastic equations with bounded time-varying delays and stochastically perturbed equations. Meanwhile, an example and simulations are given to show that impulses play an important role in pth moment and almost sure exponential stability of stochastic functional differential equations with finite delay.     

The generalized Bochner condition about classical orthogonal polynomials revisited

We bring a new proof for showing that an orthogonal polynomial sequence is classical if and only if any of its polynomial fulfils a certain differential equation of order 2k, for some k?1. So, we build those differential equations explicitly. If k=1, we get the Bochner's characterization of classical polynomials. With help of the formal computations made in Mathematica, we explicitly give those differential equations for k=1,2 and 3 for each family of the classical polynomials. Higher order differential equations can be obtained similarly.     

Nonlinear Transition Layers—The Second Painleve Transcendent

We investigate a large class of weakly nonlinear second-order ordinary differential equations with slowly varying coefficients. We show that the standard two-timing perturbation solution is not valid during the transition from oscillatory to exponentially decaying behavior. In all cases this difficulty is remedied by a nonlinear transition layer, whose leading-order character is described by one special nonlinear differential equation known as the second Painlevé transcendent (in essence a nonlinear Airy equation). The method of matched asymptotic expansions yields the desired connection formula. The second Painlevé transcendent also provides two other types of transitions: (1) between weakly nonlinear solutions (either oscillatory or exponentially decaying) and special fully nonlinear solutions, and (2) between two of these special nonlinear solutions. These special solutions are of three: different kinds: (a) slowly varying stable equilibrium solutions, (b) “exploding” solutions, and (c) solutions depending on both the fast and slow scales (which emerge from the unstable zero equilibrium solution).     

A Finite-Dimensional Dynamical Model for Gelation in Coagulation Processes

Summary. We study a finite-dimensional system of ordinary differential equations derived from Smoluchowski's coagulation equations and whose solutions mimic the behaviour of the nondensity-conserving (geling) solutions in those equations. The analytic and numerical studies of the finite-dimensional system reveals an interesting dynamic behaviour in several respects: Firstly, it suggests that some special geling solutions to Smoluchowski's equations discovered by Leyvraz can have an important dynamic role in gelation studies, and, secondly, the dynamics is interesting in its own right with an attractor possessing an unexpected structure of equilibria and connecting orbits. Received April 26, 1997; revised October 3, 1997; accepted October 23, 1997     

Small solutions of nonlinear differential equations near branching points

We construct parametric families of small branching solutions to nonlinear differential equations of the nth order near branching points. We use methods of the analytical theory of branching solutions of nonlinear equations and the theory of differential equations with a regular singular point. We illustrate the general existence theorems with an example of a nonlinear differential equation in a certain magnetic insulation problem.     

Tropical differential equations

Tropical differential equations are introduced and an algorithm is designed which tests solvability of a system of tropical linear differential equations within the complexity polynomial in the size of the system and in the absolute values of its coefficients. Moreover, we show that there exists a minimal solution, and the algorithm constructs it (in case of solvability). This extends a similar complexity bound established for tropical linear systems. In case of tropical linear differential systems in one variable a polynomial complexity algorithm for testing its solvability is designed.We prove also that the problem of solvability of a system of tropical non-linear differential equations in one variable is NP-hard, and this problem for arbitrary number of variables belongs to NP. Similar to tropical algebraic equations, a tropical differential equation expresses the (necessary) condition on the dominant term in the issue of solvability of a differential equation in power series.     

Mixed monotone iterative technique for a class of semilinear impulsive evolution equations in Banach spaces

This paper deals with the existence of mild L-quasi-solutions to the initial value problem for a class of semilinear impulsive evolution equations in an ordered Banach space E. Under a new concept of upper and lower solutions, a new monotone iterative technique on the initial value problem of impulsive evolution equations has been established. The results improve and extend some relevant results in ordinary differential equations and partial differential equations. An example is also given.