We consider the minimizers of the Gibbs free energy which couples a complex Ginzburg–Landau order parameter with a magnetic potential. It is established that the set on which the complex order parameter equals zero consists only of isolated points. Some estimates concerning the set on which the absolute value of the order parameter is small are also given. Numerical simulations are presented for the problem without a magnetic potential.

The Stefan problem with surface tension in the three-dimensional case with spherical symmetry is considered. We first establish the existence and uniqueness of the classical solution with surface tension and kinetic undercooling effects for all time, and then pass to the limit as the kinetic undercooling tends to zero. The limiting solution is the global-in-time weak solution and the local-in-time classical solution for the Stefan problem with surface tension. This solution cannot be the global-in-time classical solution. If S(t) is the radius of a solid ball in a supercooled liquid, then (1) there exists at least one point t* of discontinuity of the function S(t):

or (2) the continuous function S(t) cannot be absolutely continuous, and it maps some zero-measure set of (t*, T*) onto some set of Ω with a strictly positive measure.

We consider travelling wave solutions of a reaction–diffusion system arising in a model for infiltration with changing porosity due to reaction. We show that the travelling wave solution exists, and is unique modulo translations. A small parameter ε appears in this problem. The formal limit as ε → 0 is a free boundary problem. We show that the solution for ε > 0 tends, in a suitable norm, to the solution of the formal limit.

The bifurcation from a normally conducting to a superconducting state as an external magnetic field is lowered is examined using the Ginzburg-Landau theory. The results for three specific examples are reviewed, extended and unified in the framework of a systematic perturbation theory introduced in [1].

The stability theory for rolls in stress-free convection at finite Prandtl number is affected by coupling with low wavenumber two-dimensional mean-flow modes. In this work, a set of modified Ginzburg–Landau equations describing the onset of convection is derived which accounts for these additional modes. These equations can be used to extend the modulation equations of Zippelius & Siggia describing the breakup of rolls, bringing their stability theory into agreement with the results of Busse & Bolton.

We consider the classical Hele-Shaw situation with two parallel planes separated by a narrow gap. A blob of Newtonian fluid is sandwiched between the planes, and we suppose its plan-view to occupy a bounded, multiply-connected domain; physically, we have a viscous fluid with the holes giving rise to the multiple connectivity occupied by relatively inviscid air. The relevant free boundary condition is taken to be one of constant pressure, but we allow different pressures to act within the different holes, and at the outer boundary. The motion is driven either by injection of further fluid into the blob at certain points, or by injection of air into the holes to change their area, or by a combination of these; suction, instead of injection, is also contemplated. A general mathematical theory of the above class of problems is developed, and applied to the particular situation that arises when fluid is injected into an initially empty gap bounded by two straight, semi-infinite barriers meeting at right-angles: injection into a quarterplane. For a range of positions of the injection point, air is trapped in the corner and, invoking images, the problem is equivalent to one involving a doubly-connected blob. When there is an air vent in the corner, so that the pressure is the same on the two free boundaries in these circumstances, the air hole rapidly disappears, as might be expected. If, however, there is no air vent and we suppose the air to be incompressible, so that the area of the region occupied by the air in the plan-view remains constant, we find there to be no solution within the framework of our model. Other scenarios within this same geometry, involving both suction and injection of fluid at the injection point, and air at the corner, are also examined.

This paper considers the stability of melt-solid interfaces to eigenfunction perturbations for a system of equations which describe the melting and freezing of helium. The analysis is carried out in both planar and spherical geometries. The principal results are that when the melt is freezing, under certain far-field conditions, the interface is stable in the sense of Mullins and Sekerka. On the other hand, when the solid is melting (at least when the melting is sufficiently fast), the interface is unstable. In some circumstances these instabilities are oscillatory, with amplitude and growth rate increasing with surface tension and frequency. The last section considers the original problem of Mullins and Sekerka in the present notation.

Spontaneous potential well-logging is an important technique in petroleum exploitation. To make the corresponding log interpretation chart, it is supposed that the geometrical structure of the formation, the resistivity in each subdomain, and the spontaneous potential difference on each interface are all known; then in the direct problem, the spontaneous potential u = u(r, z) satisfies an elliptic boundary value problem with jump conditions on interfaces. At the joint points A and B of the interfaces (figure 2), the jumps of the spontaneous potential do not, in general, satisfy the compatibility condition. It turns out that it is impossible to find a piecewise H1 solution to the problem, and the standard finite element method cannot be applied to get an approximate solution. In this paper, by means of a method of removing the singularities at A and B, it is proved that the problem admits a unique weak solution that is piecewise W1,p for any fixed p with 1 ≤ p <2. Moreover, based on this method a numerical scheme is suggested, and some numerical examples and some conclusions of practical interest are given. The techniques used in this paper will find a wider applicability in other problems.

Subsonic potential flow around the leading edge of a thin two-dimensional general aerofoil with a parabolic nose at small angles of attack is analysed. Asymptotic expansions of the velocity potential function are constructed at a fixed Mach number in terms of the thickness ratio of the aerofoil in an outer region around the aerofoil and in an inner region near the nose. These expansions are matched asymptotically. The outer expansion consists of linearized aerofoil theory and its second-order problem, where the leading-edge singularity appears. The inner expansion accounts for the flow around the nose, where a stagnation point exists. A boundary-value problem is formulated in the inner region for the solution of a compressible uniform flow with the same Mach number, past a semi-infinite two-dimensional parabola. The analysis shows that the drag of any general aerofoil with a round nose is zero for every subsonic Mach number below the critical Mach number, where the drag component that is calculated by the integration of the outer pressure distribution over the aerofoil is cancelled exactly by the drag of the parabolic nose. The numerical solution of the flow field in the inner region results in the pressure distribution on the parabolic nose. The stagnation point tends to shift toward the leading edge due to the compressibility effects as the Mach number is increased at a fixed angle of attack and camber of the aerofoil. Good agreement is found in the leading edge region between the present asymptotic solution and numerical solutions of the Euler equations. Using the outer linearized solution and the nose solution a uniformly valid pressure distribution on the entire aerofoil surface in a subsonic flow is derived.

The bifurcation from a normally conducting to a superconducting state as an external magnetic field is lowered is examined using the Ginzburg–Landau theory. Linear and weakly nonlinear stability analyses are performed near the bifurcation point, and the implications of the results for each of three examples is considered.

Nonlinear evolution equations including hysteresis functionals are studied. It is the purpose of this paper to investigate the asymptotic stability of the solution as time t → + ∞. In the case when the forcing term of the equation tends to a time-independent function as t → + ∞, we shall show that the solution stabilizes at + ∞ and the system is asymptotically stable (equilibrium stability). In the case when the forcing term is periodic in time, we shall show that there is at least one periodic solution and that in some restricted cases the periodic solution is unique.

The basic constitutive relations for elastoplasticity and elastoviscoplasticity are shown to form a typical boundary layer-type stiff system of ordinary differential equations. Three numerical algorithms are discussed: (i) The singular perturbation method (O'Malley, 1971a, b; Hoppensteadt, 1971; Miranker, 1981; Smith, 1985), which yields accurate results for both the rate-independent and rate-dependent cases, where in the former case, the algorithm is explicit, whereas in the latter case, it is implicit and requires the solution of a nonlinear equation; therefore it is impractical as a constitutive algorithm for large-scale finite-element applications, where the constitutive algorithm is used a great number of times at each finite-element node. (ii) The new constitutive algorithm (Nemat-Nasser, 1991; Nemat-Nasser & Chung, 1989, 1992) which is explicit and accurate for both the rate-independent and rate-dependent cases; the underlying mathematical feature of this new method is investigated, and it is shown that it can be classified as a simplified perturbation method; computable error bounds for this algorithm are obtained, and when the flow rule is given by the commonly used power law, it is shown that the errors are very small, (iii) A modified outer-solution method, which combines the above two techniques, and is simple, explicit, and accurate.

The internal layer behaviour, in one spatial dimension, associated with two classes of Ginzbug–Landau equation with double-well nonlinearities and small diffusivities is investigated. The problems that are examined are the Ginzburg–Landau equation with and without a constant mass constraint. For the constrained problem, steady-state internal layer solutions are constructed using a formal projection method. This method is also used to derive a differential-algebraic system describing the slow dynamics of the constrained internal layer motion. The dynamics of a two-layer evolution is studied in detail. For the unconstrained problem, a nonlinear WKB-type transformation is introduced that magnifies exponentially weak layer interactions and leads to well-conditioned steady problems. A conventional singular perturbation method, without the need for exponential asymptotics, is used on the resulting transformed problem as an alternative method to construct equilibrium solutions and metastable patterns. Exponentially sensitive steady-state internal layer solutions as well as a one-layer evolution are computed accurately using the transformed problem.

This paper deals with Maxwell's equations in a quasi-stationary electromagnetic field subject to the effects of temperature. This model is encountered in the penetration of a magnetic field in substances where the electrical conductivity depends on the temperature. Similar phenomena also occur in some industrial problems such as the thermistor. Taking the effect of Joule heating into the consideration, we obtain a strongly coupled nonlinear system. Global solvability is established for this system.

A new integral representation for the Barnes double gamma function is derived. This is canonical in the sense that solutions to a class of functional difference equations of first order with trigonometrical coefficients can be expressed in terms of the Barnes function. The integral representation given here makes these solutions very simple to compute. Several well-known difference equations are solved by this method, and a treatment of the linear theory for moving contact line flow in an inviscid fluid wedge is given.

We consider the behaviour of a premixed flame anchored on a flat burner. For Lewis numbers L < L* < 1, stationary spatially periodic solutions corresponding to stationary cellular flames bifurcate from the basic solution which corresponds to a steady planar flame. We study the existence and stability of two-dimensional patterns which correspond to certain imposed symmetries by considering the evolution of N pairs of wave vectors, each of which is separated from the next by angle π/N. In the neighbourhood of the critical Lewis number L*, we derive evolution equations for the amplitudes corresponding to N = 2, which corresponds to square patterns, and N = 3, which corresponds to triangular or hexagonal patterns. We determine existence and stability results in terms of m∈(0, 1), the flow rate of the fuel, and K > 2/e, the scaled heat loss to the burner. Square patterns exist for L < L* and are stable for values of m and K above a stability boundary in the m−K plane, which has a maximum at K = K* ∼ 4.77, so that for K > K* square patterns are stable for all m. The stability of the square patterns does not vary with L. Hexagonal patterns exist for L < LH, where 1 > LH > L*. The size of the stability region increases with decreasing L < L*. For a range of values of L there is bistability, that is, for given parameter values rolls and hexagons are simultaneously stable, each with its own domain of attraction.

In this note we develop and analyse a system of equations describing a molten linear polymer being extruded in a capillary rheometer which is operating under the controlled condition that the mean velocity at the capillary inlet is maintained at a constant value. Slipping of the melt at the pipe wall is permitted and an evolution equation for the boundary slip parameter is postulated. The combined system governing the mean cross-sectional velocity and slip parameter is shown to exhibit relaxation oscillations similar to those observed in actual rheometers.

A simple method of reducing a parabolic partial differential equation to canonical form if it has only one term involving second derivatives is the following: find the general solution of the first-order equation obtained by ignoring that term and then seek a solution of the original equation which is a function of one more independent variable. Special cases of the method have been given before, but are not well known. Applications occur in fluid mechanics and the theory of finance, where the Black-Scholes equation yields to the method, and where the variable corresponding to time appears to run backwards, but there is an information-theoretic reason why it should.

A mathematical model is presented for the washing of a band, or tow, of fibres which is pulled past a submerged obstacle. A lubrication flow between the obstacle and the tow is coupled to a flow through the tow, which is modelled as a porous medium. Coupled first-order ordinary differential equations for the fluid pressure and volume flux are solved numerically and the results are analysed. The washing of the tow induced by the reversal of the viscous boundary layers lying above and below the tow is discussed briefly. In Part 2, the model is extending to allow for fibre bending. The numerical solution of the resulting differential equation shows that at a critical value of the tow porosity, the tow touches the obstacle. An existence proof for a family of touching solutions is obtained. Using the existence of these solutions, a method for constructing a solution with a rubbing region, where the gap width is zero, is outlined. The numerical solutions for these solutions are also given.