The pth ($p\geq 1$) moment exponential stability, almost surely exponential stability and stability in distribution for stochastic McKean–Vlasov equation are derived based on some distribution-dependent Lyapunov function techniques.

We derive and analyse well-posed boundary conditions for the linear shallow water wave equation. The analysis is based on the energy method and it identifies the number, location and form of the boundary conditions so that the initial boundary value problem is well-posed. A finite-volume method is developed based on the summation-by-parts framework with the boundary conditions implemented weakly using penalties. Stability is proven by deriving a discrete energy estimate analogous to the continuous estimate. The continuous and discrete analysis covers all flow regimes. Numerical experiments are presented verifying the analysis.

Assuming linear theory, the phenomenon of scattering of waves by a circular arc shaped barrier with nonuniform porosity is studied. The water region is considered to be of infinite or finite depth. Based on a judicious application of Green’s integral theorem, the corresponding boundary value problem is reduced to a hypersingular integral equation of second kind. The boundary element method and the collocation method are adopted to solve the hypersingular integral equation, and we ensure a good matching of the solutions obtained by the two methods. The reflection coefficient and energy dissipation are evaluated by using the solution of the integral equation which is then studied graphically. Different choices of distributions of pores on the barrier are considered, and we observe that the nonuniform porosity of the barrier has significant effect on the reflected wave and the energy dissipation.

The behaviour of an axisymmetric bubble in a pure liquid forced by an acoustic pressure field is analysed. The bubble is assumed to have a sharp deformable interface, which is subject both to surface tension and to Rayleigh viscosity damping. Two modelling regimes are considered. The first is a linearized solution, based on the assumption of small axisymmetric deformations to an otherwise spherical bubble. The second involves a semi-numerical solution of the fully nonlinear problem, using a novel spectral method of high accuracy. For large-amplitude nonspherical bubble oscillations, the fully nonlinear solutions show that a complicated resonance structure is possible and that curvature singularities may occur at the interface, even in the presence of surface tension. Rayleigh viscosity at the interface prevents singularity formation, but eventually causes the bubble to become purely spherical unless shape-mode resonances occur. An extended analysis is also presented for purely spherical bubbles, which allows for a more detailed study of the effects of resonance and the Rayleigh viscosity at the bubble surface.

We study the influence of a low-frequency harmonic vibration on the formation of the two-dimensional rolling solitary waves in vertically co-flowing two-layer liquid films. The system consists of two adjacent layers of immiscible fluids with the first layer being sandwiched between a vertical solid plate and the second fluid layer. The solid plate oscillates harmonically in the horizontal direction inducing Faraday waves at the liquid–liquid and liquid–air interfaces. We use a reduced hydrodynamic model derived from the Navier–Stokes equations in the long-wave approximation. Linear stability of the base flow in a flat two-layer film is determined semi-analytically using Floquet theory. We consider sub-millimetre-thick films and focus on the competition between the long-wavelength gravity-driven and finite wavelength Faraday instabilities. In the linear regime, the range of unstable wave vectors associated with the gravity-driven instability broadens at low and shrinks at high vibration frequencies. In nonlinear regimes, we find multiple metastable states characterized by solitary-like travelling waves and short pulsating waves. In particular, we find the range of the vibration parameters at which the system is multistable. In this regime, depending on the initial conditions, the long-time dynamics is dominated either by the fully developed solitary-like waves or by the shorter pulsating Faraday waves.

We propose a novel time-asymptotically stable, implicit–explicit, adaptive, time integration method (denoted by the $\theta $-method) for the solution of the fractional advection–diffusion-reaction (FADR) equations. The spectral analysis of the method (involving the group velocity and the phase speed) indicates a region of favourable dispersion for a limited range of Péclet number. The numerical inversion of the coefficient matrix is avoided by exploiting the sparse structure of the matrix in the iterative solver for the Poisson equation. The accuracy and the efficacy of the method is benchmarked using (a) the two-dimensional fractional diffusion equation, originally proposed by researchers earlier, and (b) the incompressible, subdiffusive dynamics of a planar viscoelastic channel flow of the Rouse chain melts (FADR equation with fractional time-derivative of order ) and the Zimm chain solution (). Numerical simulations of the viscoelastic channel flow effectively capture the nonhomogeneous regions of high viscosity at low fluid inertia (or the so-called “spatiotemporal macrostructures”), experimentally observed in the flow-instability transition of subdiffusive flows.