A stable linear time-invariant classical digital control system with several widely different small coefficients multiplying the lowest functions is considered. It is formulated as a multi-parameter singularly perturbed system. Perturbation methods are developed for both initial and boundary value problems based on asymptotic expansions of the perturbation parameters. The approximate solution consists of an outer solution and a number of boundary layer correction solutions equal to the number of initial conditions lost in the process of degeneration. An example is provided for illustration.

A theory is provided for the natural seiching frequencies and radiative decay rates for a shallow-water basin whose connection to open water is restricted by a submerged wall or reef. The transition from an essentially-open basin to a closed basin, as the aperture reduces to zero, is discussed using a matching procedure. Graphs of frequencies and damping factors as functions of aperture size are obtained for idealized two-dimensional shelf configurations, involving a constant-depth shallow basin connected to constant-depth, but not necessarily shallow, open water.

Scattering of membrane coupled gravity waves in deep water by partial vertical barriers is investigated by the recently developed expansion formulae for wave structure interaction problems. The horizontal thin membrane is considered to be under uniform tension and is covering the free surface. The analysis is based on the linearized theory of water waves, and by combining the kinematic and dynamic conditions at the membrane covered surface, one may derive a not so well-posed mixed boundary value problem for Laplace’s equation with third-order boundary condition. The flexible membrane is attached by a spring to the surface piercing barrier, giving suitable edge conditions for the unique solution. The boundary value problem has been converted into dual integral equations with kernels composed of trigonometric functions, which are then solved analytically. The important physical quantities such as reflection and transmission coefficients for both cases of submerged and surface piercing barriers are obtained analytically in terms of modified Bessel functions. It is found that complete reflection or transmission is possible at certain resonant frequencies for the incident membrane coupled waves. Numerical results are plotted and discussed for different values of the nondimensional membrane tension parameter.

We study the stability of inviscid, incompressible swirling flows of variable density with respect to azimuthal, normal mode disturbances. We prove that the wave velocity of neutral modes is bounded. A further refinement of Fung’s semi-elliptical instability region is given. This new instability region depends not only on the minimum Richardson number, and the lower and upper bounds for the angular velocity like Fung’s semi-ellipse, but also on the azimuthal wave number and the radii of the inner and outer cylinders. An estimation for the growth rate of unstable disturbances is obtained and it is compared to some of the recent asymptotic results.

We interpret a boundary-value problem arising in a cell growth model as a singular Sturm–Liouville problem that involves a functional differential equation of the pantograph type. We show that the probability density function of the cell growth model corresponds to the first eigenvalue and that there is a family of rapidly decaying eigenfunctions.

Bounds are presented for the modulus of the complex growth rate p of an arbitrary oscillatory perturbation, neutral or unstable, in some double-diffusive problems of relevance in oceanography, astrophysics and non-Newtonian fluid mechanics.

Gamma-Weibull variates with five parameters are defined by multiplication of gamma and Weibull densities and renormalising. Sums of independent such variates are distributed as combinations of products of gammas and confluent hypergeometric functions and are explicitly determined. Sums of independent non-identical Weibulls arise as a special case. These variates can be used to model moderately extreme scenarios between gamma and Weibull that occur in many natural applications. All results are exact.

We consider an extension of the classical Fisher–Kolmogorov equation, called the “Fisher–Stefan” model, which is a moving boundary problem on $0<x<L(t)$. A key property of the Fisher–Stefan model is the “spreading–vanishing dichotomy”, where solutions with $L(t)>L_{\text{c}}$ will eventually spread as $t\rightarrow \infty$, whereas solutions where $L(t)\ngtr L_{\text{c}}$ will vanish as $t\rightarrow \infty$. In one dimension it is well known that the critical length is $L_{\text{c}}=\unicode[STIX]{x1D70B}/2$. In this work, we re-formulate the Fisher–Stefan model in higher dimensions and calculate $L_{\text{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\text{c}}$ depends upon the dimension of the problem, and numerical solutions of the governing partial differential equation are consistent with our calculations.

Rumours have become part of our daily lives, and their spread has a negative impact on a variety of human affairs. Therefore, how to control the spread of rumours is an important topic. In this paper, we extend the classic Maki–Thompson model from a deterministic framework to a stochastic framework with a forgetting mechanism, because real-world person-to-person communications are inevitably affected by random factors. By constructing suitable stochastic Lyapunov functions, we show that the asymptotic behaviour of the stochastic rumour model is governed by the basic reproductive number. If this number is less than one, then the solution of the stochastic rumour model oscillates around the rumour-free equilibrium under extra mild conditions, indicating the extinction of the rumour with a probability of one. Otherwise, the solution always fluctuates around the endemic equilibrium under certain parametric restrictions, implying that the rumour will continually persist. In addition, we discuss a possible intervention strategy that stops the spread of rumours by strengthening the intensity of white noise, which is very different from the deterministic rumour model without white noise. Also, numerical simulations are conducted to support our analytical results.

A numerical method for calculation of the generalized Chakalov-Popoviciu quadrature formulae of Radau and Lobatto type, using the results given for the generalized Chakalov-Popoviciu quadrature formula, is given. Numerical results are included. As an application we discuss the problem of approximating a function f on the finite interval I = [a, b] by a spline function of degree m and variable defects dv, with n (variable) knots, matching as many of the initial moments of f as possible. An analytic formula for the coefficients in the generalized Chakalov-Popoviciu quadrature formula is given.

There is much interest within the mathematical biology and statistical physics community in converting stochastic agent-based models for random walkers into a partial differential equation description for the average agent density. Here a collection of noninteracting biased random walkers on a one-dimensional lattice is considered. The usual master equation approach requires that two continuum limits, involving three parameters, namely step length, time step and the random walk bias, approach zero in a specific way. We are interested in the case where the two limits are not consistent. New results are obtained using a Fokker–Planck equation and the results are highly dependent on the simulation update schemes. The theoretical results are confirmed with examples. These findings provide insight into the importance of updating schemes to an accurate macroscopic description of stochastic local movement rules in agent-based models when the lattice spacing represents a physical object such as cell diameter.

Fluid motion established by an oscillatory pressure gradient superimposed on a mean, in a tube of slowly varying section, is studied when the temperature of the tube wall varies with axial distance. Particular attention is focussed on the mean flow and steady streaming components of the oscillatory flow of higher approximation. For the velocity components, the axial component takes the pride of place, since this component is responsible for convection of nutrients to various parts of the body of a mammal in systematic circulation. A salient point in the paper concerns consequences of free convection currents at a constriction (stenosis).

The square-root fixed-interval discrete-time smoother has been used extensively in discrete recursive estimation since it was first developed by Rauch, Tung and Streibel [10]. Various people, for example Bierman [2], [3], have recognized the inherent instability in employing this kind of smoother in its original form; they have investigated implementing the recursion more stably. Bierman's paper [3] is one such contribution. In this paper we plan to present a more comprehensive development of Bierman's approach, and to show that this algorithm can be implemented more stably as a square-root smoother. Throughout this paper the fixed-interval discrete-time smoother will be referred to as the RTS smoother. Numerical results are given for the usual form of the RTS smoother, Bierman's algorithm and our square-root formulation of his algorithm. These confirm that the square-root formulation is more desirable than Bierman's algorithm, which performs better than the usual implementation of the RTS smoother.

This paper investigates the Cauchy problem for two classes of parabolic systems with localised sources. We first give the blowup criterion, and then deal with the possibilities of simultaneous blowup or non-simultaneous blowup under some suitable assumptions. Moreover, when simultaneous blowup occurs, we also establish precise blowup rate estimates. Finally, using similar ideas and methods, we shall consider several nonlocal problems with homogeneous Neumann boundary conditions.

In this paper we present a computational method for solving a class of time-lag optimal control problems with restricted phase coordinates.

A new formula for Adomian polynomials is introduced and applied to obtain truncated series solutions for fractional initial value problems with nondifferentiable functions. These kinds of equations contain a fractional single term which is examined using Jumarie fractional derivatives and fractional Taylor series for nondifferentiable functions. The property of nonlocality of these equations is examined, and the existence and uniqueness of solutions are discussed. Convergence and error analysis for the Adomian series solution are also studied. Numerical examples show the accuracy and efficiency of this formula for solving initial value problems for high-order fractional differential equations.

A theoretical investigation of the unsteady flow of a Newtonian fluid through a channel is presented using an alternative boundary condition to the standard no-slip condition, namely the Navier boundary condition, independently proposed over a hundred years ago by both Navier and Maxwell. This boundary condition contains an extra parameter called the slip length, and the most general case of a constant but different slip length on each channel wall is studied. An analytical solution for the velocity distribution through the channel is obtained via a Fourier series, and is used as a benchmark for numerical simulations performed utilizing a finite element analysis modified with a penalty method to implement the slip boundary condition. Comparison between the analytical and numerical solution shows excellent agreement for all combinations of slip lengths considered.

In this paper, we develop a discretisation algorithm with an adaptive scheme for solving a class of combined semi-infinite and semi-definite programming problems. We show that any sequence of points generated by the algorithm contains a convergent subsequence; and furthermore, each accumulation point is a local optimal solution of the combined semi-infinite and semi-definite programming problem. To illustrate the effectiveness of the algorithm, two specific classes of problems are solved. They are relaxations of quadratically constrained semi-infinite quadratic programming problems and semi-infinite eigenvalue problems.

A three-dimensional barotropic and baroclinic model is developed to simulate currents and temperature changes induced by tropical cyclones traversing the continental shelf and slope region of the Australian North West Shelf. The model is based on a layered, explicit, finite difference formulation using the nonlinear primitive equations with an embedded entrainment scheme; a mixed-surface-layer interface is defined, which is allowed to shift from one interface to another, depending on the strength of a storm. The model has been tested by simulating the currents and temperature changes induced by tropical cyclones Orson and Ian. The modelled currents and temperatures agreed well with the available measured records except near the seabed. It has been found that the pre-storm currents have very little influence on the peak of the storm-induced currents and the currents in the wake of a tropical cyclone. The model contained no coefficients which must be calibrated for a particular application and clearly illustrated the importance of the baroclinic effects on the storm-induced response over the North West Shelf of Australia.

Necessary conditions for optimal controls are given for non-linear systems with time delayed effects in both control and state variables.