Discrete and continuous coupled nonlinear integrable systems via the dressing method

A discrete analog of the dressing method is presented and used to derive integrable nonlinear evolution equations, including two infinite families of novel continuous and discrete coupled integrable systems of equations of nonlinear Schrödinger type. First, a demonstration is given of how discrete nonlinear integrable equations can be derived starting from their linear counterparts. Then, starting from two uncoupled, discrete one‐directional linear wave equations, an appropriate matrix Riemann‐Hilbert problem is constructed, and a discrete matrix nonlinear Schrödinger system of equations is derived, together with its Lax pair. The corresponding compatible vector reductions admitted by these systems are also discussed, as well as their continuum limits. Finally, by increasing the size of the problem, three‐component discrete and continuous integrable discrete systems are derived, as well as their generalizations to systems with an arbitrary number of components.     

Exact solutions of some nonlinear systems of partial differential equations by using the first integral method

In recent years, many approaches have been utilized for finding the exact solutions of nonlinear systems of partial differential equations. In this paper, the first integral method introduced by Feng is adopted for solving some important nonlinear systems of partial differential equations, including, KdV, Kaup–Boussinesq and Wu–Zhang systems, analytically. By means of this method, some exact solutions for these systems of equations are formally obtained. The results obtained confirm that the proposed method is an efficient technique for analytic treatment of a wide variety of nonlinear systems of partial differential equations.     

Iterative Solution of SPARK Methods Applied to DAEs

In this article a broad class of systems of implicit differential–algebraic equations (DAEs) is considered, including the equations of mechanical systems with holonomic and nonholonomic constraints. Solutions to these DAEs can be approximated numerically by applying a class of super partitioned additive Runge–Kutta (SPARK) methods. Several properties of the SPARK coefficients, satisfied by the family of Lobatto IIIA-B-C-C^*-D coefficients, are crucial to deal properly with the presence of constraints and algebraic variables. A main difficulty for an efficient implementation of these methods lies in the numerical solution of the resulting systems of nonlinear equations. Inexact modified Newton iterations can be used to solve these systems. Linear systems of the modified Newton method can be solved approximately with a preconditioned linear iterative method. Preconditioners can be obtained after certain transformations to the systems of nonlinear and linear equations. These transformations rely heavily on specific properties of the SPARK coefficients. A new truly parallelizable preconditioner is presented.     

Universal theory of dynamical chaos in nonlinear dissipative systems of differential equations

A new universal theory of dynamical chaos in nonlinear dissipative systems of differential equations including ordinary and partial, autonomous and non-autonomous differential equations and differential equations with delay arguments is presented in this paper. Four corner-stones lie in the foundation of this theory: the Feigenbaum’s theory of period doubling bifurcations in one-dimensional mappings, the Sharkovskii’s theory of bifurcations of cycles of an arbitrary period up to the cycle of period three in one-dimensional mappings, the Magnitskii’s theory of rotor type singular points of two-dimensional non-autonomous systems of differential equations as a bridge between one-dimensional mappings and differential equations and the theory of homoclinic cascade of bifurcations of stable cycles in nonlinear differential equations. All propositions of the theory are strictly proved and illustrated by numerous analytical and computing examples.     

Multidimensional nonlinear wave equations with multivalued solutions

We present the theory of breaking waves in nonlinear systems whose dynamics and spatial structure are described by multidimensional nonlinear hyperbolic wave equations. We obtain a general relation between systems of first-order quasilinear equations and nonlinear hyperbolic equations of higher orders, which, in particular, describe electromagnetic waves in a medium with nonlinear polarization of an arbitrary form. We use this approach to construct exact multivalued solutions of such equations and to study their spatial structure and dynamics. The results are generalized to a wide class of multidimensional equations such as d’Alembert equations, nonlinear Klein-Gordon equations, and nonlinear telegraph equations.     

MIP-based instantaneous control of mixed-integer PDE-constrained gas transport problems

We study the transient optimization of gas transport networks including both discrete controls due to switching of controllable elements and nonlinear fluid dynamics described by the system of isothermal Euler equations, which are partial differential equations in time and 1-dimensional space. This combination leads to mixed-integer optimization problems subject to nonlinear hyperbolic partial differential equations on a graph. We propose an instantaneous control approach in which suitable Euler discretizations yield systems of ordinary differential equations on a graph. This networked system of ordinary differential equations is shown to be well-posed and affine-linear solutions of these systems are derived analytically. As a consequence, finite-dimensional mixed-integer linear optimization problems are obtained for every time step that can be solved to global optimality using general-purpose solvers. We illustrate our approach in practice by presenting numerical results on a realistic gas transport network.     

Ulam-type stability for differential equations driven by measures

Based on a notion of Stieltjes derivative of a function with respect to another function, we provide Ulam–Hyers type stability results for nonlinear differential equations driven by measures on compact or on unbounded intervals, in the lack of Lipschitz continuity assumptions. In particular, one can deduce stability results for generalized differential equations, dynamic equations on time scales or impulsive differential equations (including the case of an infinite number of impulses that accumulate in the considered interval, thus allowing the study of Zeno hybrid systems).     

Decomposition of a hierarchy of nonlinear evolution equations

The generalized Hamiltonian structures for a hierarchy of nonlinear evolution equations are established with the aid of the trace identity. Using the nonlinearization approach, the hierarchy of nonlinear evolution equations is decomposed into a class of new finite-dimensional Hamiltonian systems. The generating function of integrals and their generator are presented, based on which the finite-dimensional Hamiltonian systems are proved to be completely integrable in the Liouville sense. As an application, solutions for the hierarchy of nonlinear evolution equations are reduced to solving the compatible Hamiltonian systems of ordinary differential equations.     

Universal theory of dynamical chaos in nonlinear dissipative systems of differential equations

The article presents a new universal theory of dynamical chaos in nonlinear dissipative systems of differential equations, including autonomous and nonautonomous ordinary differential equations (ODE), partial differential equations, and delay differential equations. The theory relies on four remarkable results: Feigenbaum’s period doubling theory for cycles of one-dimensional unimodal maps, Sharkovskii’s theory of birth of cycles of arbitrary period up to cycle of period three in one-dimensional unimodal maps, Magnitskii’s theory of rotor singular point in two-dimensional nonautonomous ODE systems, acting as a bridge between one-dimensional maps and differential equations, and Magnitskii’s theory of homoclinic bifurcation cascade that follows the Sharkovskii cascade. All the theoretical propositions are rigorously proved and illustrated with numerous analytical examples and numerical computations, which are presented for all classical chaotic nonlinear dissipative systems of differential equations.     

成组Broyden修正矩阵的紧凑形式与成组记忆修正算法

1 引言 成组型线性方程组 其中,p是适中的数值,由于其有相当的实际应用背景,人们一直在研究有效的数值方法,特别是近年来,实际问题中归结出来的成组型方程组,其规模越来越大,又具有稀疏结构,因而使用迭代法是一种有效的途径,目前使用比较多的是Krylov子空间方法中的Lanczos方法,CG方法,GMRES方法等等。这种成组型算法的建立,其基本出发点是使算法具有较少的计算量和存储量,具体体现在: 1)成组型算法在应用于问题(1.1)的求解时,也具有有限终止性性质,而其终止步数一般要比单个型算法的步数减少了户倍,由于成组型算法每迭代一步的计算量基本上等同于单个型算法使用户次的计算量,如此,算法的计算量会有明显的改善。 2)当A存储在二级(secondary)内存时,在迭代计算时需要不断地进行存取交换,由于成组型算法的迭代步数减少了户倍,如此,用在这种交换的时间也要减少户倍,相当有效。 3)由于在成组型算法中,出现的多是AX的形式,其中,故成组型算法便于计算并行化。 4)即使用于求解单个方程组,当A的少数几个极端特征值分离甚远时,这种成组型算法也有可能改善其收敛速度,如成组型的CG方法。 目前,这种成组型算法已体现出很大的实用计算价值,然而其进一步的理论分析还有待深入研究。     

Finite difference reaction diffusion equations with nonlinear boundary conditions

This article is concerned with numerical solutions of finite difference systems of reaction diffusion equations with nonlinear internal and boundary reaction functions. The nonlinear reaction functions are of general form and the finite difference systems are for both time-dependent and steady-state problems. For each problem a unified system of nonlinear equations is treated by the method of upper and lower solutions and its associated monotone iterations. This method leads to a monotone iterative scheme for the computation of numerical solutions as well as an existence-comparison theorem for the corresponding finite difference system. Special attention is given to the dynamical property of the time-dependent solution in relation to the steady-state solutions. Application is given to a heat-conduction problem where a nonlinear radiation boundary condition obeying the Boltzmann law of cooling is considered. This application demonstrates a bifurcation property of two steady-state solutions, and determines the dynamic behavior of the time-dependent solution. Numerical results for the heat-conduction problem, including a test problem with known analytical solution, are presented to illustrate the various theoretical conclusions. © 1995 John Wiley & Sons, Inc.     

Numerical diagnostics of solution blowup in differential equations

New simple and robust methods have been proposed for detecting poles, logarithmic poles, and mixed-type singularities in systems of ordinary differential equations. The methods produce characteristics of these singularities with a posteriori asymptotically precise error estimates. This approach is applicable to an arbitrary parametrization of integral curves, including the arc length parametrization, which is optimal for stiff and ill-conditioned problems. The method can be used to detect solution blowup for a broad class of important nonlinear partial differential equations, since they can be reduced to huge-order systems of ordinary differential equations by applying the method of lines. The method is superior in robustness and simplicity to previously known methods.     

Stochastic multiobjective problems with complementarity constraints and applications in healthcare management

We consider a class of stochastic multiobjective problems with complementarity constraints (SMOPCCs) in this paper. We derive the first-order optimality conditions including the Clarke/Mordukhovich/strong-type stationarity in the Pareto sense for the SMOPCC. Since these first-order optimality systems involve some unknown index sets, we reformulate them as nonlinear equations with simple constraints. Then, we introduce an asymptotic method to solve these constrained equations. Furthermore, we apply this methodology results to a patient allocation problem in healthcare management.     

Numerical detection and study of singularities in solutions of differential equations

New simple and robust methods are proposed for detecting singularities, such as poles, logarithmic poles, and mixed singularities, in systems of ordinary differential equations. The methods produce characteristics of these singularities with an a posteriori asymptotically precise error estimate. They are applicable in the case of an arbitrary parametrization of integral curves, including one in terms of the arc length, which is optimal for stiff and ill-conditioned problems. Following this approach, blowup solutions can be detected for a broad class of important nonlinear partial differential equations, since they are reducible by the method of lines to systems of ordinary differential equations of huge orders. The simplicity and reliability of the approach are superior to those of previously known methods.     

Numerical methods for simulation of large systems

The numerical solution of large initial value problems, including those that are derived as approximations to systems of partial differential equations, may encounter difficulties using conventional numerical methods because of stiffness (large range of eigenvalues of the associated linear system). In a nonlinear system, the eigenvalues may change greatly during the solution and a system that is initially well behaved may become stiff, yielding increased computer cost or inaccuracies. This paper contains a discussion of various definitions of stiffness, and several methods for overcoming it, including a new method for identifying and partitioning a two-time-scale system into fast and slow sub-systems. Also included are some experiences using the DARE continuous system simulation language for systems as large as 200 coupled nonlinear ordinary differential equations.     

Blow-up of the solution of a nonlinear system of equations with positive energy

We consider the Dirichlet problem for a nonlinear system of equations, continuing our study of nonlinear hyperbolic equations and systems of equations with an arbitrarily large positive energy. We use a modified Levine method to prove the blow-up.     

On a Class of Systems of Linear and Nonlinear Fredholm Integral Equations of the Third Kind with Multipoint Singularities

Based on a new approach, we show that finding solutions for a class of systems of linear (respectively, nonlinear) Fredholm integral equations of the third kind with multipoint singularities is equivalent to finding solutions of systems of linear (respectively, nonlinear) Fredholm integral equations of the second kind with additional conditions. We study the existence, nonexistence, uniqueness, and nonuniqueness of solutions for this class of systems of Fredholm integral equations of the third kind with multipoint singularities.     

AN ENCLOSURE METHOD FOR FREE BOUNDARY PROBLEMS BASED ON A LINEAR COMPLEMENTARITY PROBLEM WITH INTERVAL DATA*

This paper presents the linear complementarity problem with interval data and emphasizes its application in including the solution of an ordinary free boundary problem. This study is based on Taylor's formula with remainder and the reformulation of linear complementarity problems as nonsmooth nonlinear systems of equations.

     

Bifurcations and exact traveling wave solutions of the equivalent complex short-pulse equations

In this paper, we study the traveling wave solutions for a complex short-pulse equation of both focusing and defocusing types, which governs the propagation of ultrashort pulses in nonlinear optical fibers. It can be viewed as an analog of the nonlinear Schrodinger (NLS) equation in the ultrashort-pulse regime. The corresponding traveling wave systems of the equivalent complex short-pulse equations are two singular planar dynamical systems with four singular straight lines. By using the method of dynamical systems, bifurcation diagrams and explicit exact parametric representations of the solutions are given, including solitary wave solution, periodic wave solution, peakon solution, periodic peakon solution and compacton solution under different parameter conditions.