The nonlinear and weakly dispersive Serre equations contain higher-order dispersive terms. These include mixed spatial and temporal derivative flux terms which are difficult to handle numerically. These terms can be replaced by an alternative combination of equivalent temporal and spatial terms, so that the Serre equations can be written in conservation law form. The water depth and new conserved quantities are evolved using a second-order finite-volume scheme. The remaining primitive variable, the depth-averaged horizontal velocity, is obtained by solving a second-order elliptic equation using simple finite differences. Using an analytical solution and simulating the dam-break problem, the proposed scheme is shown to be accurate, simple to implement and stable for a range of problems, including flows with steep gradients. It is only slightly more computationally expensive than solving the shallow water wave equations.

Large-scale low-pressure systems in the atmosphere are occasionally observed to possess Kelvin–Helmholtz fingers, and an example is shown in this paper. However, these structures are hundreds of kilometres long, so that they are necessarily affected strongly by nonlinearity. They are evidently unstable and are observed to dissipate after a few days.

A model for this phenomenon is presented here, based on the usual $f$-plane equations of meteorology, assuming an atmosphere governed by the ideal gas law. Large-amplitude perturbations are accounted for, by retaining the equations in their nonlinear forms, and these are then solved numerically using a spectral method. Finger formation is modelled as an initial perturbation to the $n$th Fourier mode, and the numerical results show that the fingers grow in time, developing structures that depend on the particular mode. Results are presented and discussed, and are also compared with the predictions of the ${\it\beta}$-plane theory, in which changes of the Coriolis acceleration with latitude are included. An idealized vortex in the northern hemisphere is considered, but the results are at least in qualitative agreement with an observation of such systems in the southern hemisphere.

The steady, axisymmetric flow induced by a point sink (or source) submerged in an inviscid fluid of infinite depth is computed and the resulting deformation of the free surface is obtained. The effect of surface tension on the free surface is determined and is the new component of this work. The maximum Froude numbers at which steady solutions exist are computed. It is found that the determining factor in reaching the critical flow changes as more surface tension is included. If there is zero or a very small amount of surface tension, the limiting factor appears to be the formation of small wavelets on the free surface; but, as the surface tension increases, this is replaced by a tendency for the lowest point on the free surface to descend sharply as the Froude number is increased.

We investigate the dynamics of a susceptible infected recovered (SIR) epidemic model on small networks with different topologies, as a stepping stone to determining how the structure of a contact network impacts the transmission of infection through a population. For an SIR model on a network of $N$ nodes, there are $3^{N}$ configurations that the network can be in. To simplify the analysis, we group the states together based on the number of nodes in each infection state and the symmetries of the network. We derive analytical expressions for the final epidemic size of an SIR model on small networks composed of three or four nodes with different topological structures. Differential equations which describe the transition of the network between states are also derived and solved numerically to confirm our analysis. A stochastic SIR model is numerically simulated on each of the small networks with the same initial conditions and infection parameters to confirm our results independently. We show that the structure of the network, degree of the initial infectious node, number of initial infectious nodes and the transmission rate all significantly impact the final epidemic size of an SIR model on small networks.

We consider a coupled, logistic predator–prey system with delay. Mainly, by choosing the delay time ${\it\tau}$ as a bifurcation parameter, we show that Hopf bifurcation can occur as the delay time ${\it\tau}$ passes some critical values. Based on the normal-form theory and the centre manifold theorem, we also derive formulae to obtain the direction, stability and the period of the bifurcating periodic solution at critical values of ${\it\tau}$. Finally, numerical simulations are investigated to support our theoretical results.

We consider an approximation scheme using Haar wavelets for solving optimal path planning problems. The problem is first expressed as an optimal control problem. A computational method based on Haar wavelets in the time domain is then proposed for solving the obtained optimal control problem. A Haar wavelets integral operational matrix and a direct collocation method are used to find an approximate optimal trajectory of the original problem. Numerical results are also presented for several examples to demonstrate the applicability and efficiency of the proposed method.

We present a high-order upwind finite volume element method to solve optimal control problems governed by first-order hyperbolic equations. The method is efficient and easy for implementation. Both the semi-discrete error estimates and the fully discrete error estimates are derived. Optimal order error estimates in the sense of $L^{2}$-norm are obtained. Numerical examples are provided to confirm the effectiveness of the method and the theoretical results.