First integrals and exact solutions of the SIRI and tuberculosis models

The first integrals and exact solutions of mathematical models of epidemiology: a susceptible‐infected‐recovered‐infected (SIRI) model and a tuberculosis model with demographic growth are analyzed. These models are represented by systems of first‐order nonlinear ordinary differential equations, and this system is replaced by one which contains a second‐order ordinary differential equation. The partial Lagrangian approach is then utilized to derive the first integrals of these models. Several cases arise. Then, we utilize the derived first integrals to construct exact solutions for the models under investigation and determine new solutions. The dynamic properties of these models are studied too. Copyright © 2016 John Wiley & Sons, Ltd.     

Nonautonomous first integrals and general solutions of the KdV‐Burgers and mKdV‐Burgers equations with the source

The method for constructing first integrals and general solutions of nonlinear ordinary differential equations is presented. The method is based on index accounting of the Fuchs indices, which appeared during the Painlevé test of a nonlinear differential equation. The Fuchs indices indicate us the leading members of the first integrals for the origin differential equation. Taking into account the values of the Fuchs indices, we can construct the auxiliary equation, which allows to look for the first integrals of nonlinear differential equations. The method is used to obtain the first integrals and general solutions of the KdV‐Burgers and the mKdV‐Burgers equations with a source. The nonautonomous first integrals in the polynomials form are found. The general solutions of these nonlinear differential equations under at some additional conditions on the parameters of differential equations are also obtained. Illustrations of some solutions of the KdV‐Burgers and the mKdV‐Burgers are given.     

On the Existence and Stability of Spatially Structured Solutions of the Reaction-Diffusion Equations

The reaction-diffusion equations for the well-known ‘Brusselator’chemical kinetic model are investigated when the model is madeconsistent with the principle of detailed balance. In contrastto the original model, the corrected one does not show solutionswith ‘spatial structure’ i.e. solutions with multipleinternal maxima and multiple internal global minima in bothdependent variables. Sufficient conditions for stability ofthe solutions are given in terms of a Rayleigh quotient forgeneral boundary conditions for nonlinear reaction-diffusionequations in general. For the particular case of the ‘Brusselator’it is shown that conditions for a change of stability are muchmore unlikely to be attained as a result of the restrictionsof detailed balancing. The detailed sufficiency condition forthe stability of any steady-state solution and for no branchingfrom the ‘equilibrium’ branch solution depends onwhether the solution has global extrema inside the region, onits boundary, or both     

Spatial discretization of partial differential equations with integrals

We consider the problem of constructing spatial finite-differenceapproximations on an arbitrary fixed grid which preserve anynumber of integrals of the partial differential equation andpreserve some of its symmetries. A basis for the space of suchfinite-difference operators is constructed; most cases of interestinvolve a single such basis element. (The ‘Arakawa’Jacobian is such an element, as are discretizations satisfying‘summation by parts’ identities.) We show how thegrid, its symmetries, and the differential operator interactto affect the complexity of the finite difference.     

Construction of Euclidian Monopoles

This paper describes a procedure for the construction of monopoleson three-dimensional Euclidean space, starting from their rationalmaps. A companion paper, ‘Euclidean monopoles and rationalmaps’, to appear in the same journal, describes the assignmentto a monopole of a rational map, from CP¹ to a suitable flagmanifold. In describing the reverse direction, this paper completesthe proof of the main theorem therein. A construction of monopoles from solutions to Nahm's equations(a system of ordinary differential equations) has been well-knownfor certain gauge groups for some time. These solutions arehard to construct however, and the equations themselves becomeincreasingly unwieldy when the gauge group is not SU(2). Here, in contrast, a rational map is the only initial data.But whereas one can be reasonably explicit in moving from Nahmdata to a monopole, here the monopole is only obtained fromthe rational map after solving a partial differential equation. A non-linear flow equation, essentially just the path of steepestdescent down the Yang-Mills-Higgs functional, is set up. Itis shown that, starting from an ‘approximate monopole’- constructed explicitly from the rational map - a solutionto the flow must exist, and converge to an exact monopole havingthe desired rational map. 1991 Mathematics Subject Classification:53C07, 53C80, 58D27, 58E15, 58G11.     

Existence, Uniqueness, and Construction of the Solution of a System of Ordinary Functional Differential Equations, with Application to the Design of Perfectly Focusing Symmetric Lenses

In the design of perfectly focusing symmetric lenses, one isled, in a natural way, to a set offunctional differential equations;that is, differential equations involving composites of unknownfunctions, with initial conditions prescribed on the lens axis.This paper concentrateson those features of the equations whichmake them uniquely solvable. They are: (i) a contractivenessproperty of the equations near the axis; (ii) a uniform retardationin the arguments of thecomposite functions away from the axis.The second and third sections of this paper generalize and formalizethese properties and provide proofs of existence, uniqueness,and continuous dependence on the data for solutions of suchgeneralized systems of functional differential equations. Becauseof the lens context which motivates our study, the problem inwhich the contractiveness property (i) above holds is calledthe ‘local’ problem, and the problem in which thearguments of composite functions are uniformly retarded is calledthe ‘global’ problem. In the final section of thepaper we apply the general results of the preceding sectionsto prove existence and uniqueness of perfectly focusing symmetriclenses up to distances from the lens axis at which various typesof breakdown, discussed in the text, may occur.     

Exact master equations describing reduced dynamics of the Wigner function

Master equations of different types describe the evolution (reduced dynamics) of a subsystem of a larger system generated by the dynamic of the latter system. Since, in some cases, the (exact) master equations are relatively complicated, there exist numerous approximations for such equations, which are also called master equations. In the paper, we develop an exact master equation describing the reduced dynamics of the Wigner function for quantum systems obtained by a quantization of a Hamiltonian system with a quadratic Hamilton function. First, we consider an exact master equation for first integrals of ordinary differential equations in infinite-dimensional locally convex spaces. After this, we apply the results obtained to develop an exact master equation corresponding to a Liouville-type equation (which is the equation for first integrals of the (system of) Hamilton equation(s)); the latter master equation is called the master Liouville equation; it is a linear first-order differential equation with respect to a function of real variables taking values in a space of functions on the phase space. If the Hamilton equation generating the Liouville equation is linear, then the vector fields that define the first-order linear differential operators in the master Liouville equations are also linear, which in turn implies that for a Gaussian reference state the Fourier transform of a solution of the master Liouville equation also satisfies a linear differential equation. __________ Translated from Fundamentalnaya i Prikladnaya Matematika, Vol. 12, No. 5, pp. 203–219, 2005.     

Tollmien-Schlichting/vortex interactions in compressible boundary-layer flows

A weakly nonlinear interaction of oblique Tollmien-Schlichtingwaves and longitudinal vortices in compressible, high Reynoldsnumber, boundary-layer flow over a flat plate is consideredfor all ranges of the Mach number. The interaction equationsconsist of equations for the vortex which is indirectly forcedby the waves via a boundary condition, whereas a vortex termappears in the amplitude equation for the wave pressure. Thedownstream solution properties of interaction equations arefound to depend on the sign of an interaction coefficient. Thisparticular type of weakly nonlinear interaction was first proposedby Hall & Smith (1989), who considered incompressible flows;however, there are some errors in their formulation. Correctedresults for the incompressible regime are presented for comparisonwith those calculated for compressible flows. Compressibilityis found to have a significant effect on the interaction properties,principally through its impact on the waves and their governingmechanism, the ‘triple-deck’ structure. It is foundthat, in general, the flow quantities will grow slowly withincreasing downstream coordinate. However, for flows with Machnumber values below 2, there exists a small band of wave obliquenessangles for which the solutions terminate in abrupt, finite-distance‘break-ups’.     

Reduced differential transform method for nonlinear integral member of Kadomtsev–Petviashvili hierarchy differential equations

The main objective of this paper is to use the reduced differential to transform method (RDTM) for finding the analytical approximate solutions of two integral members of nonlinear Kadomtsev–Petviashvili (KP) hierarchy equations. Comparing the approximate solutions which obtained by RDTM with the exact solutions to show that the RDTM is quite accurate, reliable and can be applied for many other nonlinear partial differential equations. The RDTM produces a solution with few and easy computation. This method is a simple and efficient method for solving the nonlinear partial differential equations. The analysis shows that our analytical approximate solutions converge very rapidly to the exact solutions.     

Modeling weakly nonlinear acoustic wave propagation

Three weakly nonlinear models of lossless, compressible fluidflow—a straightforward weakly nonlinear equation (WNE),the inviscid Kuznetsov equation (IKE) and the Lighthill–Westerveltequation (LWE)—are derived from first principles and theirrelationship to each other is established. Through a numericalstudy of the blow-up of acceleration waves, the weakly nonlinearequations are compared to the ‘exact’ Euler equations,and the ranges of applicability of the approximate models areassessed. By reformulating these equations as hyperbolic systemsof conservation laws, we are able to employ a Godunov-type finite-differencescheme to obtain numerical solutions of the approximate modelsfor times beyond the instant of blow-up (that is, shock formation),for both density and velocity boundary conditions. Our studyreveals that the straightforward WNE gives the best results,followed by the IKE, with the LWE's performance being the poorestoverall.     

The Exact Traveling Wave Solutions to Two Integrable KdV6 Equations

The exact explicit traveling solutions to the two completely integrable sixthorder nonlinear equations KdV6 are given by using the method of dynamical systems and Cosgrove's work.It is proved that thes...     

Bounds and Self-Consistent Estimates for the Overall Properties of Nonlinear Composites

Variational ‘self-consistent’ estimates for nonlinearproblems are formulated, building on a variational formulationpreviously developed by the authors. The formulation employsa linear ‘comparison medium’ for whose propertiessome ‘self-consistent’ choice is made. In contrastto linear problems, three possible self-consistent choices presentthemselves. The results that they give are analysed for twoparticular systems (a nonlinear dielectric and a nonlinear lossycomposite) for which bounds are already available. Estimatesbased on self-consistent embedding of a single inclusion ina homogeneous matrix composed of ‘comparison material’are also developed.     

Two-dimensional nonlinear advection-diffusion in a model of surfactant spreading on a thin liquid film

The spreading of a localized monolayer of dilute, insoluble surfactant, discharged from a point source that moves at constant speed over a thin liquid film coating a planar substrate, is described according to lubrication theory by a pair of coupled nonlinear evolution equations for the monolayer concentration and the film depth h. Numerical and asymptotic techniquesare here used to show that the extent and structure of sucha spreading asymmetric monolayer can be well approximated bya single nonlinear advection–diffusion equation involving alone. At large times the solution is composed of three, spatiallydistinct, asymptotic regions: (i) a quasi-steady ‘nose’region (containing the source), in which there is a dominantbalance between two-dimensional nonlinear diffusion and advection;(ii) an ‘advective’ region, in which longitudinaladvection balances transverse diffusion; and (iii) a ‘tail’region, in which unsteady diffusion is dominant. In each region,local similarity solutions are obtained either exactly (inthe advective region) or approximately (elsewhere) by rescalingnumerical solutions of the initial-value problem. If the sourceconcentration decreases with time, it is demonstrated that the monolayer’s width is greatest in the tail region, whereasfor a source of increasing concentration the monolayer is widestin the advective region. For the simpler one-dimensional problemof a monolayer spreading from a line source, the same balanceshold but with transverse diffusion eliminated; here self-similarsolutions are found in all three regions that agree closelywith numerical solutions of the initial-value problem. Received 7 October, 1998. Revised 11 April, 2000. ⁺ antoine@mip.ups-tlse.fr Present address: Division of Theoretical Mechanics, Schoolof Mathematical Sciences, University of Nottingham , UniversityPark, Nottingham NG7 2RD, UK. Oliver.Jensen@nottingham.ac.uk.     

Timeout control scheme for overloaded M/E2/1 queue

A timeout scheme is considered for controlling an infinite ‘firstcome, first served’ overloaded single-server queue. Inthe overload situation, a customer-rejection mechanism is usedfor timing out ‘older’ customers in the queue, i.e.excluding those who have waited longer than a certain time.Applying ‘level-crossing analysis’ to an M/E₂/1queue, exact analytic expressions of performance such as thedensity and distribution functions of waiting time of the customerswho get served, the mean delay of customers, successful throughput,and ‘goodput’ are determined for this queue.     

New stress and velocity fields for highly frictional granular materials

The idealized theory for the quasi-static flow of granular materialswhich satisfy the Coulomb–Mohr hypothesis is considered.This theory arises in the limit as the angle of internal frictionapproaches /2, and accordingly these materials may be referredto as being ‘highly frictional’. In this limit,the stress field for both two-dimensional and axially symmetricflows may be formulated in terms of a single nonlinear second-orderpartial differential equation for the stress angle. To obtainan accompanying velocity field, a flow rule must be employed.Assuming the non-dilatant double-shearing flow rule, a furtherpartial differential equation may be derived in each case, thistime for the streamfunction. Using Lie symmetry methods, a completeset of group-invariant solutions is derived for both systems,and through this process new exact solutions are constructed.Only a limited number of exact solutions for gravity-drivengranular flows are known, so these results are potentially importantin many practical applications. The problem of mass flow througha two-dimensional wedge hopper is examined as an illustration.     

Exact Solutions of a Mathematical Model for Fluid Transport in Peritoneal Dialysis

A mathematical model for fluid transport in peritoneal dialysis is constructed. The model is based on a nonlinear system of two-dimensional partial differential equations with corresponding boundary and initial conditions. Using the classical Lie scheme, we establish that the base system of partial differential equations (under some restrictions on coefficients) is invariant under the infinite-dimensional Lie algebra, which enables us to construct families of exact solutions. Moreover, exact solutions with a more general structure are found using another (non-Lie) technique. Finally, it is shown that some of the solutions obtained describe the hydrostatic pressure and the glucose concentration in peritoneal dialysis. __________ Published in Ukrains'kyi Matematychnyi Zhurnal, Vol. 57, No. 8, pp. 1112–1119, August, 2005.     

The first integral method for constructing exact and explicit solutions to nonlinear evolution equations

Problems that are modeled by nonlinear evolution equations occur in many areas of applied sciences. In the present study, we deal with the negative order KdV equation and the generalized Zakharov system and derive some further results using the so‐called first integral method. By means of the established first integrals, some exact traveling wave solutions are obtained in a concise manner. Copyright © 2012 John Wiley & Sons, Ltd.     

Analytic solutions of the (2 + 1)-dimensional nonlinear evolution equations using the sine-cosine method

In this paper, we establish exact solutions for (2 + 1)-dimensional nonlinear evolution equations. The sine-cosine method is used to construct exact periodic and soliton solutions of (2 + 1)-dimensional nonlinear evolution equations. Many new families of exact traveling wave solutions of the (2 + 1)-dimensional Boussinesq, breaking soliton and BKP equations are successfully obtained. These solutions may be important of significance for the explanation of some practical physical problems. It is shown that the sine-cosine method provides a powerful mathematical tool for solving a great many nonlinear partial differential equations in mathematical physics.     

Exact travelling wave solutions of the discrete sine–Gordon equation obtained via the exp-function method

In this paper, we generalize the exp-function method, which was used to find new exact travelling wave solutions of nonlinear partial differential equations (NPDEs) or coupled nonlinear partial differential equations, to nonlinear differential–difference equations (NDDEs). As an illustration, two series of exact travelling wave solutions of the discrete sine–Gordon equation are obtained by means of the exp-function method. As some special examples, these new exact travelling wave solutions can degenerate into the kink-type solitary wave solutions reported in the open literature.