A one-dimensional model of clarifier-thickener units in engineering applications can be expressed as a conservation law with a flux that is discontinuous with respect to the spatial variable. This model also includes a singular feed source. In this paper, the clarifier-thickener model studied in a previous paper (Numer. Math.97 (2004) 25–65) is extended by a singular sink through which material is extracted from the unit. A difficulty is that in contrast to the singular source, the sink term cannot be incorporated into the flux function; rather, the sink is represented by a new non-conservative transport term. To focus on the new analytical difficulties arising due to this non-conservative term, a reduced problem is formulated, which contains the new sink term of the extended clarifier-thickener model, but not the source term and flux discontinuities. The paper is concerned with numerical methods for both models (extended and reduced) and with the well-posedness analysis for the reduced problem. For the reduced problem, a definition of entropy solutions, based on Kružkov-type entropy functions and fluxes, is provided. Jump conditions are derived and uniqueness of the entropy solution is shown. Existence of an entropy solution is shown by proving convergence of a monotone difference scheme. Two variants of the numerical scheme are introduced. Numerical examples illustrate that all three variants converge to the entropy solution, but introduce different amounts of numerical diffusion.

Rational solutions and rational-oscillatory solutions of the defocusing nonlinear Schrödinger equation are expressed in terms of special polynomials associated with rational solutions of the fourth Painlevé equation. The roots of these special polynomials have a regular, symmetric structure in the complex plane. The rational solutions verify results of Nakamura and Hirota [J. Phys. Soc. Japan, 54 (1985) 491–499] whilst the rational-oscillatory solutions appear to be new solutions of the defocusing nonlinear Schrödinger equation.

We discuss two different approaches for the analysis of the Poisson and of the non-homogeneous biharmonic equations in two dimensions. The first approach yields the solution as an integral in the complex $z$-plane (the physical plane), involving explicitly the given boundary conditions. The second approach yields an integral in the complex $k$-plane (the Fourier plane), involving the Fourier transforms of the given boundary conditions. For simple boundary value problems, such as certain problems formulated in the half complex plane, the first approach is easier. However, for more complicated problems, such as those formulated in the interior of an equilateral triangle, it appears that only the second approach can be used. Furthermore, the second approach also seems more efficient for numerical computations.

In this paper we investigate the model for the motion of a contact line over a smooth solid surface developed by Shikhmurzaev, [24]. We show that the formulation is incomplete as it stands, since the mathematical structure of the model indicates that an additional condition is required at the contact line. Recent work by Bedeaux, [4], provides this missing condition, and we examine the consequences of this for the relationship between the contact angle and contact line speed for Stokes flow, using asymptotic methods to investigate the case of small capillary number, and a boundary integral method to find the solution for general capillary number, which allows us to include the effect of viscous bending. We compare the theory with experimental data from a plunging tape experiment with water/glycerol mixtures of varying viscosities [11]. We find that we are able to obtain a reasonable fit using Shikhmurzaev's model, but that it remains unclear whether the linearized surface thermodynamics that underlies the theory provide an adequate description for the motion of a contact line.