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.

The aim of this paper is to describe a technique based on matched asymptotic expansions that allows us to derive the variation of the stress intensity factors in a homogeneous isotropic elastic medium under plane strain deformation. The case of antiplane shearing is also considered.

Using the two-scale convergence approach, we derive equations which govern transversal time-harmonic waves through a layered medium taking the form of a poroelastic composite saturated with a viscous fluid. To improve convergence, we construct a corrector. We study how wave speed and attenuation time depend on porosity and frequency. We prove that the Darcy permeability and the acoustic permeability in the Biot equations do not coincide.

We investigate self-similar solutions of the thin film equation in the case of zero contact angle boundary conditions on a finite domain. We prove existence and uniqueness of such a solution and determine the asymptotic behaviour as the exponent in the equation approaches the critical value at which zero contact angle boundary conditions become untenable. Numerical and power-series solutions are also presented.

For the Cauchy problem for the nonlinear infiltration equation $$\left\{\begin{array}{@{}l@{\qquad}l} u_{t}=\frac{1}{m}(u^{m})_{xx},&x\in{\mathbb{R}}, t>0,m\geq{}1,\\[3pt] u|_{t=0}=u_{0}(x),&x\in{\mathbb{R}}, \end{array} \right.$$ we use its linear solution $u(x,t,1)$ to approach the nonlinear solution $u(x,t,m)$, and obtain the explicit estimate: $$\int_{0}^{T}\int_{\mathbb{R}}|u(x,t,m)-u(x,t,1)|^{2}\,dx\,dt{} \leq{}(C^{\ast}(m-1))^{2},$$ where $C^{\ast}=O(T^{\gamma})$ and $\gamma=\frac{1+m-\alpha}{2(1+m)}$ for any $0<\alpha<1$.

Two-dimensional flow is considered in a fluid bead located in the gap between a pair of contra-rotating cylinders and bounded by two curved menisci. The stability of such bead flows with two inlet films, and hence no contact line, are analysed as the roll speed ratio S is increased. One of the inlet films can be regarded as an ‘input flux’ whilst the other is a ‘returning film’ whose thickness is specified as a fraction ζ of the outlet film on that roll. The flow is modelled via lubrication theory and for Ca [Lt ] 1, where Ca represents the capillary number, boundary conditions are formally developed that account for S ≠ 1 and the non-constant gap. It is shown that there is a qualitative difference in the results between the single and double inlet film models unless small correction terms to the pressure drops at the interfaces are taken into account. Futhermore, it is shown that the inclusion of these small terms produces an O(1) effect on the prediction of the critical value of S at which bead break occurs. When the limits of the returning film fraction are examined it is found that as ζ → 0 results are in good agreement with those for the single inlet film. Further it is shown for a fixed input flux that as ζ → 1 a transition from bead break to upstream flooding of the nip can occur and multiple two-dimensionally stable solutions exist. For a varying input flux and fixed and ‘sufficiently large’ values of ζ there is a critical input flux &λmacr;(ζ) such that as S is increased from zero:

(i) bead break occurs for λ < &λmacr;;

(ii) upstream flooding occurs for λ > &λmacr;;

(iii) when λ = &λmacr; the flow becomes neutrally stable at a specific value of S beyond which there exist two steady solutions (two-dimensionally stable) leading to bead break and upstream flooding, respectively.

A two-dimensional computational study is performed on the texturing of fibre-filled nematic liquid crystals using the Landau-de Gennes model describing the spatio-temporal evolution of the second moment of the orientation distribution function or quadrupolar tensor order parameter. The investigation is performed on a consistent computational domain comprising a square array of four circular fibres embedded within a unit square containing a uniaxial low molar mass calamitic liquid crystal. Interest is focused on the role of temperature, boundary conditions and their effect on the nucleation and evolution of defect structures. Thermal effects are characterised below and above the temperature at which the nematic state is stable. Simulations in the stable nematic state serves as a scenario for investigating the effect of imposing different external boundary conditions, namely periodic and Dirichlet; the former describes a square lattice array of fibres embedded in a nematic liquid crystal, and the latterdescribes a four-fibre arrangement in an aligned nematic material. In each case, the influence of temperature is characterised, with defect structures forming and either remaining or splitting into lower strength defects. For fibre lattices, splitting transitions of defects at the centre of the domain occur at a critical temperature, but for the four-fibre arrangement, defect transitions occur continuously over a temperature range. The discontinuous defect splitting transition in fibre arrays occurs at lower temperatures than the continuous defect transformation in the four-fibre arrangement. At sufficiently low temperatures, the four-fibre arrangement and the fibre lattice give the same texture consisting of two disclination lines close to each fibre. The evolution of the texture with respect to temperature can be characterised as a change from single-fibre mode at low temperature to a collective mode with a centre-located heterogeneity at higher temperature. At higher temperatures, in the stable isotropic state, it is shown that surface-induced ordering arising from the fibre/liquid crystal interaction propagates into the bulk forming thin disclination lattices around the four-fibre configuration.

Steady-state diffusion in long axisymmetric structures is considered. The goal is to assess one-dimensional approximations by comparing them with axisymmetric eigenfunction expansions. Two problems are considered in detail: a finite tube with one end that is partly absorbing and partly reflecting; and two finite coaxial tubes with different cross-sectional radii joined together abruptly. Both problems may be modelled using effective boundary conditions, containing a parameter known as the trapping rate. We show that trapping rates depend on the lengths of the finite tubes (and that they decay slowly as these lengths increase) and we show how trapping rates are related to blockage coefficients, which are well known in the context of potential flow along tubes of infinite length.

Consider a deterministic dynamical system in a domain containing a stable equilibrium, e.g., a particle in a potential well. The particle, independent of initial conditions, eventually reaches the bottom of the well. If however, the particle is subjected to white noise, due, e.g., to collisions with a population of smaller, lighter particles comprising the medium through which the particle travels, a dramatic difference in the behaviour of the Brownian particle occurs. The particle will exit the well. The natural questions then are how long will it take for it to exit and from where on the boundary of the domain of attraction of the deterministic equilibrium (the rim of the well) will it exit. We compute the mean first passage time to the boundary and the mean probabilities of the exit positions. When the noise is small each quantity satisfies a singularly perturbed deterministic boundary value problem. We treat the problem by the method of matched asymptotic expansions (MAE) and generalizations thereof. MAE has been used successfully to solve problems in many applications. However, there exist problems for which MAE does not suffice. Among these are problems exhibiting boundary layer resonance, i.e., the problem of ‘spurious solutions’, which led some to conclude that this was ‘the failure of MAE’. We present a physical argument and four mathematical arguments to modify or augment MAE to make it successful. Finally, we discuss applications of the theory.

The connection between Residual Neural Networks (ResNets) and continuous-time control systems (known as NeurODEs) has led to a mathematical analysis of neural networks, which has provided interesting results of both theoretical and practical significance. However, by construction, NeurODEs have been limited to describing constant-width layers, making them unsuitable for modelling deep learning architectures with layers of variable width. In this paper, we propose a continuous-time Autoencoder, which we call AutoencODE, based on a modification of the controlled field that drives the dynamics. This adaptation enables the extension of the mean-field control framework originally devised for conventional NeurODEs. In this setting, we tackle the case of low Tikhonov regularisation, resulting in potentially non-convex cost landscapes. While the global results obtained for high Tikhonov regularisation may not hold globally, we show that many of them can be recovered in regions where the loss function is locally convex. Inspired by our theoretical findings, we develop a training method tailored to this specific type of Autoencoders with residual connections, and we validate our approach through numerical experiments conducted on various examples.

An analysis is made of the damping of sound and structural vibrations by vorticity production in the apertures of a bias flow, perforated elastic plate. Unsteady motion causes vorticity to be generated at the aperture edges; the vorticity and its energy are swept away by the bias flow and result in a net loss of acoustic and vibrational energy. In this paper we investigate the interaction of an arbitrary fluid-structure disturbance with a small circular aperture in the presence of a high Reynolds number, low Mach number bias flow. By considering the limit in which the aperture is small compared to the length scale of the impinging disturbance, it is shown that the effect of the interaction can be represented by a concentrated source in the plate bending wave equation consisting of a delta function and two of its axisymmetric derivatives. A generalized bending wave equation is then formulated for a plate perforated with an homogeneous distribution of small, bias flow circular apertures. This equation is used to predict the attenuation of sound and resonant bending waves by vorticity production. Acoustic damping is found to be significant provided the fluid loading is sufficiently small for the plate to be regarded as rigid (e.g. for an aluminium plate in air when the frequency is not too small). On the other hand, a bending wave is effectively damped only when the fluid loading is large enough for the wave to produce a substantial pressure drop across the plate; when this occurs the predicted attenuations are comparable with those usually achieved by the application of elastomeric damping materials. Numerical predictions are presented for steel and aluminium plates in air and water.

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.

A spectral method is developed to study the steady and unsteady flow of fluid into a line sink from a horizontally confined aquifer, and the results are compared to solutions obtained implementing the finite element package COMSOLTM. The aquifer or drain is considered to be confined below so that the solutions are fundamentally unsteady. Comparison is made between the two methods in determining the drawdown of the surface.

The present work is part of a series of studies conducted by the authors on analytical models of avascular tumour growth that exhibit both geometrical anisotropy and physical inhomogeneity. In particular, we consider a tumour structure formed in distinct ellipsoidal regions occupied by cell populations at a certain stage of their biological cycle. The cancer cells receive nutrient by diffusion from an inhomogeneous supply and they are subject to also an inhomogeneous pressure field imposed by the tumour microenvironment. It is proved that the lack of symmetry is strongly connected to a special condition that should hold between the data imposed by the tumour’s surrounding medium, in order for the ellipsoidal growth to be realizable, a feature already present in other non-symmetrical yet more degenerate models. The nutrient and the inhibitor concentration, as well as the pressure field, are provided in analytical fashion via closed-form series solutions in terms of ellipsoidal eigenfunctions, while their behaviour is demonstrated by indicative plots. The evolution equation of all the tumour’s ellipsoidal interfaces is postulated in ellipsoidal terms and a numerical implementation is provided in view of its solution. From the mathematical point of view, the ellipsoidal system is the most general coordinate system that the Laplace operator, which dominates the mathematical models of avascular growth, enjoys spectral decomposition. Therefore, we consider the ellipsoidal model presented in this work, as the most general analytic model describing the avascular growth in inhomogeneous environment. Additionally, due to the intrinsic degrees of freedom inherited to the model by the ellipsoidal geometry, the ellipsoidal model presented can be adapted to a very populous class of avascular tumours, varying in figure and in orientation.

This paper addresses a free boundary problem for a steady, uniform patch of vorticity surrounding a single flat plate of zero thickness and finite length. Exact solutions to this problem have previously been found in terms of conformal maps represented by Cauchy-type integrals. Here, however, it is demonstrated how, by considering an associated circular-arc polygon and using ideas from automorphic function theory, these maps can be expressed in a simple non-integral form.

It is proved that c-cyclical monotonicity is a sufficient condition for optimality in the multi-marginal optimal transport problem with Coulomb repulsive cost. The notion of c-splitting set and some mild regularity property are the tools. The result may be extended to Coulomb like costs.

We study a problem related to the spin-coating process in which a fluid coats a rotating surface. Our interest lies in the contact-line region for which we propose a simplified travelling wave approximation. We construct solutions to this problem by a shooting method that matches solution branches in the contact-line region and in the interior of the droplet. Furthermore, we prove uniqueness and qualitative properties of the solution connected to the fourth-order nature of the equation, such as a global maximum in the film height close to the contact line, elevated from the average height of the film.

Finite element methods developed for unfitted meshes have been widely applied to various interface problems. However, many of them resort to non-conforming spaces for approximation, which is a critical obstacle for the extension to $\textbf{H}(\text{curl})$ equations. This essential issue stems from the underlying Sobolev space $\textbf{H}^s(\text{curl};\,\Omega)$, and even the widely used penalty methodology may not yield the optimal convergence rate. One promising approach to circumvent this issue is to use a conforming test function space, which motivates us to develop a Petrov–Galerkin immersed finite element (PG-IFE) method for $\textbf{H}(\text{curl})$-elliptic interface problems. We establish the Nédélec-type IFE spaces and develop some important properties including their edge degrees of freedom, an exact sequence relating to the $H^1$ IFE space and optimal approximation capabilities. We analyse the inf-sup condition under certain assumptions and show the optimal convergence rate, which is also validated by numerical experiments.

We study plane shear waves in a continuum model for a saturated soil, assuming a hypoplastic flow rule. Following an earlier paper on velocity-controlled waves [1], we now specify a time- periodic stress disturbance on the boundary of a semi-infinite medium. As before, the stress saturates, tending at large times to a position-independent value, away from the boundary. By contrast, the dependence of this asymptotic state on the boundary data is strikingly different: in the previous case, this dependence was wildly discontinuous, while here it is very simple and smooth.

Short-time existence and uniqueness results in Sobolev spaces are proved for Hele-Shaw flow with kinetic undercooling and for Stokes flow without surface tension. In both cases, the flow is driven by arbitrarily distributed sources and sinks in the interior of the liquid domain. The proofs are based on a general approach consisting of the reformulation of the problem as a Cauchy problem for a nonlinear, nonlocal evolution equation on the unit sphere, quasilinearization by equivariance, investigation of the linearization, and Galerkin approximations. In the situation discussed here, the linearized evolution operator is a first-order differential operator, and thus the evolution equation is of hyperbolic type. Finally, a brief survey of the properties of the evolution equations that arise from Hele-Shaw flow and Stokes flow with and without regularization is given.