A non-linear reaction–diffusion system with cross-diffusion describing the COVID-19 outbreak is studied using the Lie symmetry method. A complete Lie symmetry classification is derived and it is shown that the system with correctly specified parameters admits highly non-trivial Lie symmetry operators, which do not occur for all known reaction–diffusion systems. The symmetries obtained are also applied for finding exact solutions of the system in the most interesting case from applicability point of view. It is shown that the exact solutions derived possess typical properties for describing the pandemic spread under 1D approximation in space and lead to the distributions, which qualitatively correspond to the measured data of the COVID-19 spread in Ukraine.

In this paper, we study a mathematical model for an infectious disease caused by a virus such as Cholera without lifetime immunity. Due to the different mobility for susceptible, infected human and recovered human hosts, the diffusion coefficients are assumed to be different. The resulting system is governed by a strongly coupled reaction–diffusion system with different diffusion coefficients. Global existence and uniqueness are established under certain assumptions on known data. Moreover, global asymptotic behaviour of the solution is obtained when some parameters satisfy certain conditions. These results extend the existing results in the literature. The main tool used in this paper comes from the delicate theory of elliptic and parabolic equations. Moreover, the energy method and Sobolev embedding are used in deriving a priori estimates. The analysis developed in this paper can be employed to study other epidemic models in biological, ecological and health sciences.

Silicon is produced in submerged arc furnaces which are heated by electric currents passing through the furnace. It is important to understand the distribution of heating within the furnace in order to accurately model the silicon production process, yet many existing studies neglect aspects of this current flow. In the present paper, we formulate a model that couples the electrical current to thermal, material flow and chemical processes in the furnace. We then exploit disparate timescales to homogenise the model over the timescale of the alternating current, deriving averaged equations for the slow evolution of the system. Our numerical simulations predict a minimum applied current that is required in order to obtain steady-state solutions of the homogenised model and show that for high enough applied currents, two spatially heterogeneous steady-state solutions exist, with distinct crater sizes. We show that the system evolves to the steady state with a larger crater radius and explain this behaviour in terms of the overall power balance typically found within a furnace. We find that the industrial practice of stoking furnaces increases the overall rate of material consumption in the furnace, thereby improving the efficiency of silicon production.

Existence of non-negative weak solutions is shown for a full curvature thin-film model of a liquid thin film flowing down a vertical fibre. The proof is based on the application of a priori estimates derived for energy-entropy functionals. Long-time behaviour of these weak solutions is analysed and, under some additional constraints for the model parameters and initial values, convergence towards a travelling wave solution is obtained. Numerical studies of energy minimisers and travelling waves are presented to illustrate analytical results.

We consider patterns formed in a two-dimensional thin film on a planar substrate with a Derjaguin disjoining pressure and periodic wettability stripes. We rigorously clarify some of the results obtained numerically by Honisch et al. [Langmuir 31: 10618–10631, 2015] and embed them in the general theory of thin-film equations. For the case of constant wettability, we elucidate the change in the global structure of branches of steady-state solutions as the average film thickness and the surface tension are varied. Specifically we find, by using methods of local bifurcation theory and the continuation software package AUTO, both nucleation and metastable regimes. We discuss admissible forms of spatially non-homogeneous disjoining pressure, arguing for a form that differs from the one used by Honisch et al., and study the dependence of the steady-state solutions on the wettability contrast in that case.

We consider the no-flux initial-boundary value problem for the cross-diffusive evolution system:

\begin{eqnarray*} \left\{ \begin{array}{ll} u_t = u_{xx} - \chi \big(\frac{u}{v} \partial_x v \big)_x - uv +B_1(x,t), \qquad & x\in \Omega, \ t>0, \\[1mm] v_t = v_{xx} +uv - v + B_2(x,t), \qquad & x\in \Omega, \ t>0, \end{array} \right. \end{eqnarray*}

$\chi=2$

In bounded intervals $\Omega\subset\mathbb{R}$ and with prescribed suitably regular non-negative functions $B_1$ and $B_2$ , we first prove the existence of global classical solutions for any choice of $\chi>0$ and all reasonably regular non-negative initial data.

We next address the issue of determining the qualitative behaviour of solutions under appropriate assumptions on the asymptotic properties of $B_1$ and $B_2$ . Indeed, for arbitrary $\chi>0$ , we obtain boundedness of the solutions given strict positivity of the average of $B_2$ over the domain; moreover, it is seen that imposing a mild decay assumption on $B_1$ implies that u must decay to zero in the long-term limit. Our final result, valid for all $\chi\in\left(0,\frac{\sqrt{6\sqrt{3}+9}}{2}\right),$ which contains the relevant value $\chi=2$ , states that under the above decay assumption on $B_1$ , if furthermore $B_2$ appropriately stabilises to a non-trivial function $B_{2,\infty}$ , then (u,v) approaches the limit $(0,v_\infty)$ , where $v_\infty$ denotes the solution of

\begin{eqnarray*} \left\{ \begin{array}{l} -\partial_{xx}v_\infty + v_\infty = B_{2,\infty}, \qquad x\in \Omega, \\[1mm] \partial_x v_{\infty}=0, \qquad x\in\partial\Omega. \end{array} \right. \end{eqnarray*}

$\chi$

$\chi$

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.