In this paper, we consider the existence and stability of singular patterns in a fractional Ginzburg–Landau equation with a mean field. We prove the existence of three types of singular steady-state patterns (double fronts, single spikes, and double spikes) by solving their respective consistency conditions. In the case of single spikes, we prove the stability of single small spike solution for sufficiently large spatial period by studying an explicit non-local eigenvalue problem which is equivalent to the original eigenvalue problem. For the other solutions, we prove the instability by using the variational characterisation of eigenvalues. Finally, we present the results of some numerical computations of spike solutions based on the finite difference methods of Crank–Nicolson and Adams–Bashforth.

Suspended fibres significantly alter fluid rheology, as exhibited in for example solutions of DNA, RNA and synthetic biological nanofibres. It is of interest to determine how this altered rheology affects flow stability. Motivated by the fact thermal gradients may occur in biomolecular analytic devices, and recent stability results, we examine the problem of Rayleigh–Bénard convection of the transversely isotropic fluid of Ericksen. A transversely isotropic fluid treats these suspensions as a continuum with an evolving preferred direction, through a modified stress tensor incorporating four viscosity-like parameters. We consider the linear stability of a stationary, passive, transversely isotropic fluid contained between two parallel boundaries, with the lower boundary at a higher temperature than the upper. To determine the marginal stability curves the Chebyshev collocation method is applied, and we consider a range of initially uniform preferred directions, from horizontal to vertical, and three orders of magnitude in the viscosity-like anisotropic parameters. Determining the critical wave and Rayleigh numbers, we find that transversely isotropic effects delay the onset of instability; this effect is felt most strongly through the incorporation of the anisotropic shear viscosity, although the anisotropic extensional viscosity also contributes. Our analysis confirms the importance of anisotropic rheology in the setting of convection.

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I–ΔtL which maps the identity (no preconditioner) for Δt≪1 to Laplacian preconditioning for Δt≫1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for , away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

In this paper we study Bénard-Marangoni convection in confined containers where a thin fluid layer is heated from below. We consider containers with circular, square and hexagonal cross-sections. For Marangoni numbers close to the critical Marangoni number, the flow patterns are dominated by the appearance of the well-known hexagonal convection cells. The main purpose of this computational study is to explore the possible patterns the system may end up in for a given set of parameters. In a series of numerical experiments, the coupled fluid-thermal system is started with a zero initial condition for the velocity and a random initial condition for the temperature. For a given set of parameters we demonstrate that the system can end up in more than one state. For example, the final state of the system may be dominated by a steady convection pattern with a fixed number of cells, however, the same system may occasionally end up in a steady pattern involving a slightly different number of cells, or it may end up in a state where most of the cells are stationary, while one or more cells end up in an oscillatory state. For larger aspect ratio containers, we are also able to reproduce dislocations in the convection pattern, which have also been observed experimentally. It has been conjectured that such imperfections (e.g., a localized star-like pattern) are due to small irregularities in the experimental setup (e.g., the geometry of the container). However, we show, through controlled numerical experiments, that such phenomena may appear under otherwise ideal conditions. By repeating the numerical experiments for the same non-dimensional numbers, using a different random initial condition for the temperature in each case, we are able to get an indication of how rare such events are. Next, we study the effect of symmetrizing the initial conditions. Finally, we study the effect of selected geometry deformations on the resulting convection patterns.

A mathematical model for the steady, mixed convection heat and mass transfer along a semi-infinite vertical plate embedded in a micropolar fluid in the presence of Soret and Dufour effects is presented. The non-linear governing equations and their associated boundary conditions are initially cast into dimensionless forms using local similarity transformations. The resulting system of equations is then solved numerically using the Keller-box method. The numerical results are compared and found to be in good agreement with previously published results as special cases of the present investigation. The non-dimensional velocity, microrotation, temperature and concentration profiles are displayed graphically for different values of coupling number, Soret and Dufour numbers. In addition, the skin-friction coefficient, the Nusselt number and Sherwood number are shown in a tabular form.