Recently Biscari, DiCarlo and Turzi [Anisotropic wave propagation in nematic liquid crystals, Soft Matter 10, 8296–8307.] proposed a theory for nematoacustics which models nematic liquid crystals as nematic elastomers with molecular relaxation. Here, we extend the analysis of this theory to account for the director motion possibly induced by the propagation of a sound wave. We find that the director vibration is related to the - usually small - anisotropy of the molecular distribution, thus providing a justification to the relative high ultrasonic intensities required to observe non-negligible acousto-optic responses.

In this paper we study some aspects of the integrability problem for polynomial vector fields $\dot{x}=P(x,y)$, $\skew1\dot{y}=Q(x,y)$. We analyze the possible existence of first integrals of the form $I(x,y)=e^{ h_1(x) \prod_{k=1}^r (y-a_k(x))/ \prod_{j=1}^s(y-f_j(x))} h_2(x)$$\prod_{i=1}^{\ell} (y-g_i(x))^{\alpha_i}$, where $g_i(x)$ and $f_j(x)$ are unknown particular solutions of $dy/dx=Q(x,y)/P(x,y)$, $\alpha_i \in \mathbb{C}$ are unknown constants, and $a_k(x)$, $h_1(x)$ and $h_2(x)$ are unknown functions. We give an algorithmic method to determine if the polynomial vector field has a first integral of the form above described. In the case when some of the particular solutions remain arbitrary and the other ones are explicitly determined or are functionally related to the arbitrary particular solutions, we will obtain a generalized nonlinear superposition principle, see [6]. In the case when all the particular solutions $g_i(x)$ and $f_j(x)$ are determined, they are algebraic functions and our algorithm gives an alternative method for determining such type of solutions.

We consider a non-stationary incompressible non-Newtonian Stokes system in a porous medium with characteristic size of the pores ϵ and containing a thin fissure of width ηϵ. The viscosity is supposed to obey the power law with flow index $\frac{5}{3}\leq q\leq 2$. The limit when size of the pores tends to zero gives the homogenized behaviour of the flow. We obtain three different models depending on the magnitude ηϵ with respect to ϵ: if ηϵ ≪ $\varepsilon^{q\over 2q-1}$ the homogenized fluid flow is governed by a time-dependent non-linear Darcy law, while if ηϵ ≫ $\varepsilon^{q\over 2q-1}$ is governed by a time-dependent non-linear Reynolds problem. In the critical case, ηϵ ≈ $\varepsilon^{q\over 2q-1}$, the flow is described by a time-dependent non-linear Darcy law coupled with a time-dependent non-linear Reynolds problem.

Most of the results available on the inverse problem of determining loads acting on elastic beams or plates under transverse vibration refer to single beam or single plate. In this paper, we consider the determination of sources in multi-span systems obtained by connecting either two Euler–Bernoulli elastic beams or two rectangular Kirchhoff–Love elastic plates. The material of the structure is assumed to be homogeneous and isotropic. The transverse load is of the form g(t)f(x), where g(t) is a known function of time and f(x) is the unknown term depending on the position variable x. Under slight a priori assumptions, we prove a uniqueness result for f(x) in terms of observations of the dynamic response taken at interior points of the structure in an arbitrary small interval of time. A numerical implementation of the method is included to show the possible application of the results in the practical identification of the source term.

In both historical and modern conflicts, space plays a critical role in how interactions occur over time. Despite its importance, the spatial distribution of adversaries has often been neglected in mathematical models of conflict. In this paper, we propose an entropy-maximising spatial interaction method for disaggregating the impact of space, employing a general notion of ‘threat’ between two adversaries. This approach addresses a number of limitations that are associated with partial differential equation approaches to spatial disaggregation. We use this method to spatially disaggregate the Richardson model of conflict escalation, and then explore the resulting model with both analytical and numerical treatments. A bifurcation is identified that dramatically influences the resulting spatial distribution of conflict and is shown to persist under a range of model specifications. Implications of this finding for real-world conflicts are discussed.

New directions in Markov processes and research on master equations are showcased by example. The utility of Magnus expansions for handling time-varying rates is demonstrated. The useful notion in applied mathematics often turns out to be the pseudospectra and not simply the eigenvalues. We highlight that general principle with our own examples of Markov processes where exact eigenvalues are found and contrasted with the large errors produced by standard numerical methods. As a motivating application, isomerisation provides a running example and an illustration of our approaches to chemical kinetics. We also present a brief example of a totally asymmetric exclusion process.

We consider two-dimensional bubbles in a corner flow in a Hele–Shaw cell of a viscous incompressible fluid that occupies the complement to a bubble. We discuss the governing equations, some basic properties of the free interface of the bubbles, their geometry, and construct explicit solutions that present asymmetric long bubbles analogous to the famous Saffman–Taylor fingers in a wedge of arbitrary angle $\alpha\in (0,2\pi)$.

New modulation equations for hexagonal patterns in reaction–diffusion systems are derived for parameter régimes corresponding to the onset of patterns. These systems include additional nonlinearities which are not present in Rayleigh–Bénard convection or Swift–Hohenberg type models. The dynamics of hexagonal and roll patterns are studied using a combination of analytical and computational approaches which exploit the hexagonal structure of the modulation equations. The investigation demonstrates instabilities and new phenomena not found in other systems, and is applied to patterns of flame fronts in a certain model of burner stabilized flames.

A mathematical model of the development of an oil and gas field is presented. The field development process is treated as sequential in nature. Completion of a well and its production are considered to be random processes. The model uses results from renewal theory where the completion of a well and failure to produce economical amount of oil or gas are analogous to the failure of a component. In principle, the theory described can give the complete probability distribution associated with a field development. Explicit expressions are given for the expected value and variance of the number of completed wells.

We consider an autonomous Hamiltonian system $\ddot {q}+V_q(q)=0$ in ${\bf R}^2$, where the potential $V$ has a global maximum at the origin and singularities at some points $\xi_1$, $\xi_2 \in {\bf R}^2 \setminus \{0\}$. Under some compactness conditions on $V$ at infinity and assuming a strong force type condition at the singularities, we study, using variational arguments, the existence of various types of heteroclinic and homoclinic solutions of the system.

Online parameter identification in time-dependent differential equations from time course observations related to the physical state can be understood as a non-linear inverse and ill-posed problem and appears in a variety of applications in science and engineering. The feature as well as the challenge of online identification is that sensor data have to be continuously processed during the operation of the real dynamic process in order to support simulation-based decision making. In this paper we present an online parameter identification method that is based on a non-linear parameter-to-output operator and, as opposed to methods available so far, works both for finite- and infinite-dimensional dynamical systems, e.g., both for ordinary differential equations and time-dependent partial differential equations. A further advantage of the method suggested is that it renders typical restrictive assumptions such as full state observability, linear parametrisation of the underlying model and data differentiation or filtering unnecessary. Assuming existence of a solution for exact data, a convergence analysis based on Lyapunov theory is presented. Numerical illustrations given are by means of online identification both of aerodynamic coefficients in a 3DoF-longitudinal aircraft model and of a (distributed) conductivity coefficient in a heat equation.

We consider the equation of turbulent burst propagation

with m > 1 and k ≥ 0, or, in an equivalent formulation

We consider initial data u(x, 0) which are radially symmetric and have compact support; x = l(t) = r(T) describes the interface of the solution.

This paper investigates global dynamics of an infection age-space structured cholera model. The model describes the vibrio cholerae transmission in human population, where infection-age structure of vibrio cholerae and infectious individuals are incorporated to measure the infectivity during the different stage of disease transmission. The model is described by reaction–diffusion models involving the spatial dispersal of vibrios and the mobility of human populations in the same domain Ω ⊂ ℝn. We first give the well-posedness of the model by converting the model to a reaction–diffusion model and two Volterra integral equations and obtain two constant equilibria. Our result suggest that the basic reproduction number determines the dichotomy of disease persistence and extinction, which is achieved by studying the local stability of equilibria, disease persistence and global attractivity of equilibria.

We consider a family of multi-phase Stefan problems for a certain one-dimensional model of cell-to-cell adhesion and diffusion, which takes the form of a non-linear forward–backward parabolic equation. In each material phase the cell density stays either high or low, and phases are connected by jumps across an ‘unstable’ interval. We develop an existence theory for such problems which allows for the annihilation of phases and the subsequent continuation of solutions. Stability results for the long-time behaviour of solutions are also obtained, and, where necessary, the analysis is complemented by numerical simulations.

We consider the dynamic process of an elastic body in unilateral frictional contact with a rigid foundation. Friction is modelled with the Coulomb law with a coefficient that depends on the slip velocity. To allow for velocity discontinuities we use the elastodynamic (hyperbolic) framework. Nevertheless, this does not lead to a well-posed problem. To remedy this, we perturb the solution of the elastodynamic problem in a thin layer next to the contact boundary. This is a generalisation of an approach previously studied in a one-dimensional case. We establish existence and uniqueness results for the perturbed and regularised problem and provide an interpretation of this perturbation.

Following the derivation of amplitude equations through a new two-time-scale method [O'Malley, R. E., Jr. & Kirkinis, E (2010) A combined renormalization group-multiple scale method for singularly perturbed problems. Stud. Appl. Math. 124, 383–410], we show that a multi-scale method may often be preferable for solving singularly perturbed problems than the method of matched asymptotic expansions. We illustrate this approach with 10 singularly perturbed ordinary and partial differential equations.

We present a method of solution of a class of fracture problems in the theory of elasticity. The method can be applied to any problem reducible to Poisson's equation, e.g. heat conduction and mass diffusion in solids, theory of consolidation and the like. The novelty of the paper is that we address regions of layered composites with notches, or, in a particular case, with a crack. Within the framework of classical analysis, we apply Fourier and Mellin transforms, 'fit' them together, and reduce the problem to solving a singular integral equation with fixed singularities on a semi-axis. We show the existence and uniqueness of solutions of the equations under consideration, and justify the asymptotics necessary for applications. We show the practical usefulness of the method on the examples of an antiplane problem of fracture mechanics. From our solution, we are able to find the stress intensity factor in the case when a crack tip penetrates a layered composite consisting of 60 layers, and show the limits of applicability of the anisotropic model of such composites.

This paper investigates the stability of a fully parabolic–parabolic-fluid (PP-fluid) system of the Keller–Segel–Navier–Stokes type in a bounded planar domain under the natural volume-filling hypothesis. In the limit of fast signal diffusion, we first show that the global classical solutions of the PP-fluid system will converge to the solution of the corresponding parabolic–elliptic-fluid (PE-fluid) system. As a by-product, we obtain the global well-posedness of the PE-fluid system for general large initial data. We also establish some new exponential time decay estimates for suitable small initial cell mass, which in particular ensure an improvement of convergence rate on time. To further explore the stability property, we carry out three numerical examples of different types: the nontrivial and trivial equilibriums, and the rotating aggregation. The simulation results illustrate the possibility to achieve the optimal convergence and show the vanishment of the deviation between the PP-fluid system and PE-fluid system for the equilibriums and the drastic fluctuation of error for the rotating solution.

We consider a non-local version of the Cahn–Hilliard equation characterized by the presence of a fractional diffusion operator, and which is subject to fractional dynamic boundary conditions. Our system generalizes the classical system in which the dynamic boundary condition was used to describe any relaxation dynamics of the order-parameter at the walls. The proposed fractional dynamic boundary condition appears to be more general in the sense that it incorporates non-local effects which were completely ignored in the classical approach. We aim to deduce well-posedness and regularity results as well as to establish the existence of finite-dimensional attractors for this system.

An analysis is given of brine transport through a porous medium, which incorporates the effect of volume changes due to variations in the salt concentration. Two specific situations are investigated which lead to self-similarity. We develop the existence and uniqueness theory for the corresponding ordinary differential equations, and give a number of qualitative properties of the solutions. In particular, we present an asymptotic expression for the solution in terms of the relative density difference (ρs−ρf)/ρf. Finally, we show some numerical results. It is found that the volume changes have a noticeable effect on the mass transport only when salt concentrations are large.