In this paper we revisit the problem of finding an orthogonal similarity transformation that puts an $n\times n$ matrix $A$ in a block upper-triangular form that reveals its Jordan structure at a particular eigenvalue $\lambda_0$. The obtained form in fact reveals the dimensions of the null spaces of $(A-\lambda_0 I)^i$ at that eigenvalue via the sizes of the leading diagonal blocks, and from this the Jordan structure at $\lambda_0$ is then easily recovered. The method starts from a Hessenberg form that already reveals several properties of the Jordan structure of $A$.

The largest part of the Earth's atmosphere ozone is located in the stratosphere, forming the so-called ozone layer. This layer played a key role in the development of life on Earth and still protects the planet from the most Dangerous ultraviolet radiation. After the discovery of the high ozone depletion potential of some anthropogenic origin substances (e.g. chlorofluorocarbons), some limitations in the production of the major ozone-depleting substances (ODS) have been applied with the Montreal Protocol in 1987.

In this paper we consider drug binding in the arterialwall following delivery by a drug-eluting stent. Whilst it is now generally accepted that a non-linear saturable reversible binding model is required to properly describe the binding process, the precise form of the binding model varies between authors. Our particular interest in this manuscript is in assessing to what extent modelling specific and non-specific binding in the arterial wall as separate phases is important.

We present a method for applying a class of velocity-dependent forces within a multicomponent lattice Boltzmann equation simulation that is designed to recover continuum regime incompressible hydrodynamics. This method is applied to the problem, in two dimensions, of constraining to uniformity the tangential velocity of a vesicle membrane implemented within a recent multicomponent lattice Boltzmann simulation method, which avoids the use of Lagrangian boundary tracers.

Thanks to innovative sample-preparation and sequencing technologies, gene expression in individual cells can now be measured for thousands of cells in a single experiment. Since its introduction, single-cell RNA sequencing (scRNA-seq) approaches have revolutionized the genomics field as they created unprecedented opportunities for resolving cell heterogeneity by exploring gene expression profiles at a single-cell resolution. However, the rapidly evolving field of scRNA-seq invoked the emergence of various analytics approaches aimed to maximize the full potential of this novel strategy.

We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is well known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux.

We propose a numerical method for approximating integro-differential equations arising in age-of-infection epidemic models. The method is based on a non-standard finite differences approximation of the integral term appearing in the equation. The study of convergence properties and the analysis of the qualitative behavior of the numerical solution show that it preserves all the basic properties of the continuous model with no restrictive conditions on the step-length h of integration and that it recovers the continuous dynamic as h tends to zero.

A (2+2)-dimensional kinetic equation, directly inspired by the run-and-tumble modeling of chemotaxis dynamics is studied so as to derive a both ''2D well-balanced'' and ''asymptotic-preserving'' numerical approximation. To this end, exact stationary regimes are expressed by means of Laplace transforms of Fourier-Bessel solutions of associated elliptic equations. This yields a scattering S-matrix which permits to formulate a timemarching scheme in the form of a convex combination in kinetic scaling.