On the role of quadratic oscillations in nonlinear Schrödinger equations

We consider a nonlinear semi-classical Schrödinger equation for which it is known that quadratic oscillations lead to focusing at one point, described by a nonlinear scattering operator. If the initial data is an energy bounded sequence, we prove that the nonlinear term has an effect at leading order only if the initial data have quadratic oscillations; the proof relies on a linearizability condition (which can be expressed in terms of Wigner measures). When the initial data is a sum of such quadratic oscillations, we prove that the associate solution is the superposition of the nonlinear evolution of each of them, up to a small remainder term. In an appendix, we transpose those results to the case of the nonlinear Schrödinger equation with harmonic potential.     

An instability theorem for a certain sixth order nonlinear delay differential equation

The aim of this paper is to establish an instability theorem for a certain sixth order nonlinear delay differential equation. The proof of the theorem is based on the use of Lyapunov–Krasovskii functional approach. By this work, we improve an instability result obtained in the literature for a certain sixth order nonlinear differential equation without delay to the instability of the zero solution of a certain sixth order nonlinear delay differential equation.     

A General Approach to Obtain Series Solutions of Nonlinear Differential Equations

Based on homotopy, which is a basic concept in topology, a general analytic method (namely the homotopy analysis method) is proposed to obtain series solutions of nonlinear differential equations. Different from perturbation techniques, this approach is independent of small/large physical parameters. Besides, different from all previous analytic methods, it provides us with a simple way to adjust and control the convergence of solution series. Especially, it provides us with great freedom to replace a nonlinear differential equation of order n into an infinite number of linear differential equations of order k , where the order k is even unnecessary to be equal to the order n . In this paper, a nonlinear oscillation problem is used as example to describe the basic ideas of the homotopy analysis method. We illustrate that the second-order nonlinear oscillation equation can be replaced by an infinite number of (2κ)th-order linear differential equations, where κ≥ 1 can be any a positive integer. Then, the homotopy analysis method is further applied to solve a high-dimensional nonlinear differential equation with strong nonlinearity, i.e., the Gelfand equation. We illustrate that the second-order two or three-dimensional nonlinear Gelfand equation can be replaced by an infinite number of the fourth or sixth-order linear differential equations, respectively. In this way, it might be greatly simplified to solve some nonlinear problems, as illustrated in this paper. All of our series solutions agree well with numerical results. This paper illustrates that we might have much larger freedom and flexibility to solve nonlinear problems than we thought traditionally. It may keep us an open mind when solving nonlinear problems, and might bring forward some new and interesting mathematical problems to study.     

Approximation order from stability for nonlinear subdivision schemes

This paper proves approximation order properties of various nonlinear subdivision schemes. Building on some recent results on the stability of nonlinear multiscale transformations, we are able to give very short and concise proofs. In particular we point out an interesting connection between stability properties and approximation order for nonlinear subdivision schemes.     

On surjective second order non-linear Markov operators and associated nonlinear integral equations

It was known that orthogonality preserving property and surjectivity of nonlinear Markov operators, acting on finite dimensional simpleces, are equivalent. It turns out that these notions are no longer equivalent when such kind of operators are considered over on infinite dimensional spaces. In the present paper, we find necessary and sufficient condition to be equivalent of these notions, for the second order nonlinear Markov operators. To do this, we fully describe all surjective second order nonlinear Markov operators acting on infinite dimensional simplex. As an application of this result, we provided some sufficient conditions for the existence of positive solutions of nonlinear integral equations whose domain are not compact.     

Synchronization of integer order and fractional order Chua’s systems using robust observer

In this paper we present a fractional order Chua’s circuit that behaves chaotically based on the use of a fractional order low pass filter. Next, an integer order robust observer will be designed to synchronize the fractional order Chua’s circuit as well as integer order Chua’s circuit with unknown nonlinearity. This method consists in designing a Luenberger like observer appended with an estimator of the unknown nonlinear function. The estimator assumes that the nonlinear function is slowly varying and that the observer converges quickly and uses the backward difference formula to approximate the state derivative. The efficiency of the proposed method is confirmed using numerical simulations.     

Chaos control in passive walking dynamics of a compass-gait model

The compass-gait walker is a two-degree-of-freedom biped that can walk passively and steadily down an incline without any actuation. The mathematical model of the walking dynamics is represented by an impulsive hybrid nonlinear model. It is capable of displaying cyclic motions and chaos. In this paper, we propose a new approach to controlling chaos cropped up from the passive dynamic walking of the compass-gait model. The proposed technique is to linearize the nonlinear model around a desired passive hybrid limit cycle. Then, we show that the nonlinear model is transformed to an impulsive hybrid linear model with a controlled jump. Basing on the linearized model, we derive an analytical expression of a constrained controlled Poincaré map. We present a method for the numerical simulation of this constrained map where bifurcation diagrams are plotted. Relying on these diagrams, we show that the linear model is fairly close to the nonlinear one. Using the linearized controlled Poincaré map, we design a state feedback controller in order to stabilize the fixed point of the Poincaré map. We show that this controller is very efficient for the control of chaos for the original nonlinear model.     

Coherent States for Systems of L 2-Supercritical Nonlinear Schrödinger Equations

We consider the propagation of wave packets for a nonlinear Schrödinger equation, with a matrix-valued potential, in the semi-classical limit. For a matrix-valued potential, Strichartz estimates are available under long range assumptions. Under these assumptions, for an initial coherent state polarized along an eigenvector, we prove that the wave function remains in the same eigenspace, at leading order, in a scaling such that nonlinear effects cannot be neglected. We also prove a nonlinear superposition principle for these nonlinear wave packets.     

A penalty function method for solving inverse optimal value problem

In order to consider the inverse optimal value problem under more general conditions, we transform the inverse optimal value problem into a corresponding nonlinear bilevel programming problem equivalently. Using the Kuhn–Tucker optimality condition of the lower level problem, we transform the nonlinear bilevel programming into a normal nonlinear programming. The complementary and slackness condition of the lower level problem is appended to the upper level objective with a penalty. Then we give via an exact penalty method an existence theorem of solutions and propose an algorithm for the inverse optimal value problem, also analysis the convergence of the proposed algorithm. The numerical result shows that the algorithm can solve a wider class of inverse optimal value problem.     

High-order splitting schemes for semilinear evolution equations

We first derive necessary and sufficient stiff order conditions, up to order four, for exponential splitting schemes applied to semilinear evolution equations. The main idea is to identify the local splitting error as a sum of quadrature errors. The order conditions of the quadrature rules then yield the stiff order conditions in an explicit fashion, similarly to that of Runge–Kutta schemes. Furthermore, the derived stiff conditions coincide with the classical non-stiff conditions. Secondly, we propose an abstract convergence analysis, where the linear part of the vector field is assumed to generate a group or a semigroup and the nonlinear part is assumed to be smooth and to satisfy a set of compatibility requirements. Concrete applications include nonlinear wave equations and diffusion-reaction processes. The convergence analysis also extends to the case where the nonlinear flows in the exponential splitting scheme are approximated by a sufficiently accurate one-step method.     

A second‐order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection

In this article, we study a new second‐order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for higher‐order‐in‐time temporal discretizations is how to ensure an unconditional energy stability without compromising numerical efficiency or accuracy. We propose a framework for designing a second‐order numerical scheme with unconditional energy stability using the BDF method with constant coefficient stabilizing terms. Based on the unconditional energy stability property that we establish, we derive an stability for the numerical solution and provide an optimal convergence analysis. To deal with the highly nonlinear four‐Laplacian term at each time step, we apply efficient preconditioned steepest descent and preconditioned nonlinear conjugate gradient algorithms to solve the corresponding nonlinear system. Various numerical simulations are presented to demonstrate the stability and efficiency of the proposed schemes and solvers. Comparisons with other second‐order schemes are presented.     

Large Time Behavior via the Method of ?-Trajectories

The method of ?-trajectories is presented in a general setting as an alternative approach to the study of the large-time behavior of nonlinear evolutionary systems. It can be successfully applied to the problems where solutions suffer from lack of regularity or when the leading elliptic operator is nonlinear. Here we concentrate on systems of a parabolic type and apply the method to an abstract nonlinear dissipative equation of the first order and to a class of equations pertinent to nonlinear fluid mechanics. In both cases we prove the existence of a finite-dimensional global attractor and the existence of an exponential attractor.     

A tau method for nonlinear dynamical systems

In this work we present a new Tau method for the solution of nonlinear systems of differential equations which are linear in the derivative of highest order and polynomial in the remaining. We avoid the linearization of the problem by associating to it a nonlinear algebraic system and combine a forward substitution with the Tau method. We develop an adaptive step by step version of this alternative nonlinear tau method and we apply it to several nonlinear dynamical systems.     

Homoclinic solutions for nonlinear general fourth‐order differential equations

This work provides sufficient conditions for the existence of homoclinic solutions of fourth‐order nonlinear ordinary differential equations. Using Green's functions, we formulate a new modified integral equation that is equivalent to the original nonlinear equation. In an adequate function space, the corresponding nonlinear integral operator is compact, and it is proved an existence result by Schauder's fixed point theorem. Copyright © 2017 John Wiley & Sons, Ltd.     

A second‐order finite difference approximation for a mathematical model of erythropoiesis

We present a second‐order finite difference scheme for approximating solutions of a mathematical model of erythropoiesis, which consists of two nonlinear partial differential equations and one nonlinear ordinary differential equation. We show that the scheme achieves second‐order accuracy for smooth solutions. We compare this scheme to a previously developed first‐order method and show that the first order method requires significantly more computational time to provide solutions with similar accuracy. We also compare this numerical scheme with other well‐known second‐order methods and show that it has better capability in approximating discontinuous solutions. Finally, we present an application to recovery after blood loss. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Blow-up of solutions of wave equation with a nonlinear boundary condition and interior focusing source of variable order of growth

In this paper, being investigated an initial-boundary value problem for a one-dimensional wave equation with a nonlinear source of variable order and nonlinear dissipation at the boundary. The existence of a local solution of the problem under consideration is proved. Then the question of the absence of global solutions is investigated. Depending on the relationship between the order of growth of the nonlinear source and the nonlinear boundary dissipation, different results are obtained on the blow-up of weak solutions in a finite time interval.     

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation (PME)is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13–32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on $f$log $f$ as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete $W$^1,∞norm and a refined estimate is applied to derive the optimal error order.     

The convergence ball and error analysis of the two-step Secant method

Under the assumption that the nonlinear operator has Lipschitz continuous divided differences for the first order, we obtain an estimate of the radius of the convergence ball for the two-step secant method. Moreover, we also provide an error estimate that matches the convergence order of the two-step secant method. At last, we give an application of the proposed theorem.     

On the control of a flock by a leader

The collective and individual behaviors of the elements of a set of identical units are investigated from the point of view of their response to a coordination signal broadcast by a leader. In terms of the model used here, we show that for such systems nonlinearity plays a critical role. A new method for establishing controllability of nonlinear, replicated systems is given. Similarly, stabilization also depends on high order nonlinear effects. For the purpose of stabilization, a distinction between odd order nonlinear terms and even order nonlinear terms plays a role.     

A Robin–Dirichlet problem for a nonlinear wave equation with the source term containing a nonlinear integral

In this paper, we consider the Robin–Dirichlet problem for a nonlinear wave equation with the source term containing a nonlinear integral. Using the Faedo–Galerkin method and the linearization method for nonlinear terms, we prove the existence and uniqueness of a weak solution. We also discuss an asymptotic expansion of high order in a small parameter of a weak solution.