This paper proves the energy equality for distributional solutions to fractional Navier-Stokes equations, which gives a new proof and covers the classical result of Galdi [Proc. Amer. Math. Soc. 147 (2019), 785–792].

Evaporation within porous media is both a multiscale and interface-driven process, since the phase change at the evaporating interfaces within the pores generates a vapour flow and depends on the transport of vapour through the porous medium. While homogenised models of flow and chemical transport in porous media allow multiscale processes to be modelled efficiently, it is not clear how the multiscale effects impact the interface conditions required for these homogenised models. In this paper, we derive a homogenised model, including effective interface conditions, for the motion of an evaporation front through a porous medium, using a combined homogenisation and boundary layer analysis. This analysis extends previous work for a purely diffusive problem to include both gas flow and the advective–diffusive transport of material. We investigate the effect that different microscale models describing the chemistry of the evaporation have on the homogenised interface conditions. In particular, we identify a new effective parameter, $\mathcal{L}$, the average microscale interface length, which modifies the effective evaporation rate in the homogenised model. Like the effective diffusivity and permeability of a porous medium, $\mathcal{L}$ may be found by solving a periodic cell problem on the microscale. We also show that the different microscale models of the interface chemistry result in fundamentally different fine-scale behaviour at, and near, the interface.

The invariant Galton–Watson (IGW) tree measures are a one-parameter family of critical Galton–Watson measures invariant with respect to a large class of tree reduction operations. Such operations include the generalized dynamical pruning (also known as hereditary reduction in a real tree setting) that eliminates descendant subtrees according to the value of an arbitrary subtree function that is monotone nondecreasing with respect to an isometry-induced partial tree order. We show that, under a mild regularity condition, the IGW measures are attractors of arbitrary critical Galton–Watson measures with respect to the generalized dynamical pruning. We also derive the distributions of height, length, and size of the IGW trees.

We consider the global existence for the following fully parabolic chemotaxis system with two populations

\[\left\{ \begin{array}{@{}ll} \partial_tu_i=\kappa_i\Delta u_i-\chi_i\nabla\cdot(u_i\nabla v),\quad i\in\{1,2\}, & x\in\Omega,\ t>0, \\ v_t=\Delta v-v+u_1+u_2, & x\in\Omega,\ t>0,\\ u_i(x,t=0)=u_{i0}(x),\quad v(x,t=0)=v_0(x), & x\in\Omega, \end{array} \right. \]

$\Omega =\mathbb {R}^2$

$\Omega =B_R(0)\subset \mathbb {R}^2$

$\kappa _i,\chi _i>0,$

$i=1,2$

$0=\Delta v-v+u_1+u_2$

$0=\Delta v+u_1+u_2$

The frameworks of Baikov–Gazizov–Ibragimov (BGI) and Fushchich–Shtelen (FS) approximate symmetries are used to study symmetry properties of partial differential equations with a small parameter. In general, it is shown that unlike the case of ordinary differential equations (ODEs), unstable BGI point symmetries of unperturbed partial differential equations (PDEs) do not necessarily yield local approximate symmetries for the perturbed model. While some relations between the BGI and FS approaches can be established, the two methods yield different approximate symmetry classifications. Detailed classifications are presented for two nonlinear PDE families. The second family includes a one-dimensional wave equation describing the wave motion in a hyperelastic material with a single family of fibers. For this model, approximate symmetries can be used to compute approximate closed-form solutions. Wave breaking times are found numerically and using the approximate solutions, which yield comparable results.

Let $\sigma \in (0,\,2)$, $\chi ^{(\sigma )}(y):={\mathbf 1}_{\sigma \in (1,2)}+{\mathbf 1}_{\sigma =1} {\mathbf 1}_{y\in B(\mathbf {0},\,1)}$, where $\mathbf {0}$ denotes the origin of $\mathbb {R}^n$, and $a$ be a non-negative and bounded measurable function on $\mathbb {R}^n$. In this paper, we obtain the boundedness of the non-local elliptic operator

\[ Lu(x):=\int_{\mathbb{R}^n}\left[u(x+y)-u(x)-\chi^{(\sigma)}(y)y\cdot\nabla u(x)\right]a(y)\,\frac{{\rm d}y}{|y|^{n+\sigma}} \]

$\mathrm {BMO}(\mathbb {R}^n)\cap (\bigcup _{p\in (1,\infty )}L^p(\mathbb {R}^n))$

$\mathrm {BMO}(\mathbb {R}^n)$

$H^1(\mathbb {R}^n)$

$H^1(\mathbb {R}^n)$

$\lambda \in (0,\,\infty )$

$Lu-\lambda u=f$

$\mathbb {R}^n$

$f\in \mathrm {BMO}(\mathbb {R}^n)\cap (\bigcup _{p\in (1,\infty )}L^p(\mathbb {R}^n))$

$H^1(\mathbb {R}^n)$

$\mathrm {BMO}(\mathbb {R}^n)$

$H^1(\mathbb {R}^n)$

$L^p(\mathbb {R}^n)$

$p\in (1,\,\infty )$

$p=1$

$p=\infty$

We consider positive solutions to a class of quasilinear elliptic problems involving the Hardy potential under zero Dirichlet boundary condition. Via moving plane method, proving a weak comparison principle, we prove symmetry and monotonicity properties for the solutions defined on strictly convex symmetric domains.

Narrow escape and narrow capture problems which describe the average times required to stop the motion of a randomly travelling particle within a domain have applications in various areas of science. While for general domains, it is known how the escape time decreases with the increase of the trap sizes, for some specific 2D and 3D domains, higher-order asymptotic formulas have been established, providing the dependence of the escape time on the sizes and locations of the traps. Such results allow the use of global optimisation to seek trap arrangements that minimise average escape times. In a recent paper (Iyaniwura (2021) SIAM Rev. 63(3), 525–555), an explicit size- and trap location-dependent expansion of the average mean first passage time (MFPT) in a 2D elliptic domain was derived. The goal of this work is to systematically seek global minima of MFPT for $1\leq N\leq 50$ traps in elliptic domains using global optimisation techniques and compare the corresponding putative optimal trap arrangements for different values of the domain eccentricity. Further, an asymptotic formula for the average MFPT in elliptic domains with N circular traps of arbitrary sizes is derived, and sample optimal configurations involving non-equal traps are computed.

New classes of conditionally integrable systems of nonlinear reaction–diffusion equations are introduced. They are obtained by extending a well-known nonclassical symmetry of a scalar partial differential equation to a vector equation. New exact solutions of nonlinear predator–prey systems with cross-diffusion are constructed. Infinite dimensional classes of exact solutions are made available for such nonlinear systems. Some of these solutions decay towards extinction and some oscillate or spiral around an interior fixed point. It is shown that the conditionally integrable systems are closely related to the standard diffusive Lotka–Volterra system, but they have additional features.

We study the asymptotic behaviour of the periodically mixed Zaremba problem. We cover the part of the boundary by a chess board with a small period (square size) $\varepsilon$ and impose the Dirichlet condition on black and the Neumann condition on white squares. As $\varepsilon \to 0$, we get the effective boundary condition which is always of the Dirichlet type. The Dirichlet data on the boundary, however, depend on the ratio between the magnitudes of the two boundary values.

Symmetries and adjoint-symmetries are two fundamental (coordinate-free) structures of PDE systems. Recent work has developed several new algebraic aspects of adjoint-symmetries: three fundamental actions of symmetries on adjoint-symmetries; a Lie bracket on the set of adjoint-symmetries given by the range of a symmetry action; a generalised Noether (pre-symplectic) operator constructed from any non-variational adjoint-symmetry. These results are illustrated here by considering five examples of physically interesting nonlinear PDE systems – nonlinear reaction-diffusion equations, Navier-Stokes equations for compressible viscous fluid flow, surface-gravity water wave equations, coupled solitary wave equations and a nonlinear acoustic equation.

We devote this paper to study semi-stable nonconstant radial solutions of $S_k(D^2u)=w(\left \vert x \right \vert )g(u)$ on the Euclidean space $\mathbb {R}^n$. We establish pointwise estimates and necessary conditions for the existence of such solutions (not necessarily bounded) for this equation. For bounded solutions we estimate their asymptotic behaviour at infinity. All the estimates are given in terms of the spatial dimension $n$, the values of $k$ and the behaviour at infinity of the growth rate function of $w$.

In this paper, we consider the existence and stability of singular patterns in a fractional Ginzburg–Landau equation with a mean field. We prove the existence of three types of singular steady-state patterns (double fronts, single spikes, and double spikes) by solving their respective consistency conditions. In the case of single spikes, we prove the stability of single small spike solution for sufficiently large spatial period by studying an explicit non-local eigenvalue problem which is equivalent to the original eigenvalue problem. For the other solutions, we prove the instability by using the variational characterisation of eigenvalues. Finally, we present the results of some numerical computations of spike solutions based on the finite difference methods of Crank–Nicolson and Adams–Bashforth.

This note studies local integral gradient bounds for distributional solutions of a large class of partial differential inequalities with diffusion in divergence form and power-like first-order terms. The applications of these estimates are two-fold. First, we show the (sharp) global Hölder regularity of distributional semi-solutions to this class of diffusive PDEs with first-order terms having supernatural growth and right-hand side in a suitable Morrey class posed on a bounded and regular open set $\Omega$. Second, we provide a new proof of entire Liouville properties for inequalities with superlinear first-order terms without assuming any one-side bound on the solution for the corresponding homogeneous partial differential inequalities. We also discuss some extensions of the previous properties to problems arising in sub-Riemannian geometry and also to partial differential inequalities posed on noncompact complete Riemannian manifolds under appropriate area-growth conditions of the geodesic spheres, providing new results in both these directions. The methods rely on integral arguments and do not exploit maximum and comparison principles.

In this paper, we consider the following Hardy–Littlewood–Sobolev inequality with extended kernel(0.1)

\begin{equation} \int_{\mathbb{R}_+^{n}}\int_{\partial\mathbb{R}^{n}_+} \frac{x_n^{\beta}}{|x-y|^{n-\alpha}}f(y)g(x) {\rm d}y{\rm d}x\leq C_{n,\alpha,\beta,p}\|f\|_{L^{p}(\partial\mathbb{R}_+^{n})} \|g\|_{L^{q'}(\mathbb{R}_+^{n})}, \end{equation}

$f\in L^{p}(\partial \mathbb {R}_+^{n})$

$g\in L^{q'}(\mathbb {R}_+^{n})$

$p,\,\ q'\in (1,\,\infty )$

$\beta \geq 0$

$\alpha +\beta >1$

$\frac {n-1}{n}\frac {1}{p}+\frac {1}{q'}-\frac {\alpha +\beta -1}{n}=1$

We prove the existence of all extremal functions for (0.1). We show that if $f$ and $g$ are extremal functions for (0.1) then both of $f$ and $g$ are radially decreasing. Moreover, we apply the regularity lifting method to obtain the smoothness of extremal functions. Finally, we derive the sufficient and necessary condition of the existence of any nonnegative nontrivial solutions for the Euler–Lagrange equations by using Pohozaev identity.

We establish multiplicity results for the following class of quasilinear problemsP

\begin{equation*} \left\{ \begin{array}{@{}l} -\Delta_{\Phi}u=f(x,u) \quad \mbox{in} \quad \Omega, \\ u=0 \quad \mbox{on} \quad \partial \Omega, \end{array} \right. \end{equation*}

$\Delta _{\Phi }u=\text {div}(\varphi (x,|\nabla u|)\nabla u)$

$\Phi (x,t)=\int _{0}^{|t|}\varphi (x,s)s\,ds$

$\Omega \subset \mathbb {R}^{N}$

$\Omega _N$

$\Omega _p$

$\overline {\Omega _N}\cap \overline {\Omega _p}=\emptyset$

$(P)$

$-\Delta _{\Phi }$

$-\Delta _N$

$\Omega _N$

$-\Delta _p$

$\Omega _p$

$f:\Omega \times \mathbb {R}\to \mathbb {R}$

$e^{\alpha |t|^{\frac {N}{N-1}}}$

$\Omega _N$

$|t|^{p^{*}-2}t$

$\Omega _p$

$|t|$

$\alpha >0$

$p^{*}=\frac {Np}{N-p}$

$1< p< N$

$f$

$f$

$\Omega$

$(P)$

Infinitesimal symmetries of a partial differential equation (PDE) can be defined algebraically as the solutions of the linearisation (Frechet derivative) equation holding on the space of solutions to the PDE, and they are well-known to comprise a linear space having the structure of a Lie algebra. Solutions of the adjoint linearisation equation holding on the space of solutions to the PDE are called adjoint-symmetries. Their algebraic structure for general PDE systems is studied herein. This is motivated by the correspondence between variational symmetries and conservation laws arising from Noether’s theorem, which has a modern generalisation to non-variational PDEs, where infinitesimal symmetries are replaced by adjoint-symmetries, and variational symmetries are replaced by multipliers (adjoint-symmetries satisfying a certain Euler-Lagrange condition). Several main results are obtained. Symmetries are shown to have three different linear actions on the linear space of adjoint-symmetries. These linear actions are used to construct bilinear adjoint-symmetry brackets, one of which is a pull-back of the symmetry commutator bracket and has the properties of a Lie bracket. The brackets do not use or require the existence of any local variational structure (Hamiltonian or Lagrangian) and thus apply to general PDE systems. One of the symmetry actions is shown to encode a pre-symplectic (Noether) operator, which leads to the construction of symplectic 2-form and Poisson bracket for evolution systems. The generalised KdV equation in potential form is used to illustrate all of the results.

In this paper, we establish gradient continuity for solutions to

\[ (\partial_t - \operatorname{div}(A(x) \nabla ))^{s} u =f,\quad s \in (1/2, 1), \]

$f$

$L\left (\frac {n+2}{2s-1}, 1\right )$

$f$

In this work we study the Cauchy problem in Gevrey spaces for a generalized class of equations that contains the case $b=0$ of the $b$-equation. For the generalized equation, we prove that it is locally well-posed for initial data in Gevrey spaces. Moreover, as we move to global well-posedness, we show that for a particular choice of the parameter in the equation the local solution is global analytic in both time and spatial variables.

For $N\geq 2$, a bounded smooth domain $\Omega$ in $\mathbb {R}^{N}$, and $g_0,\, V_0 \in L^{1}_{loc}(\Omega )$, we study the optimization of the first eigenvalue for the following weighted eigenvalue problem:

\[ -\Delta_p \phi + V |\phi|^{p-2}\phi = \lambda g |\phi|^{p-2}\phi \text{ in } \Omega, \quad \phi=0 \text{ on } \partial \Omega, \]

$g$

$V$

$g_0$

$V_0$

$(\underline {g},\,\underline {V})$

$(\overline {g},\,\overline {V})$

$g_0$

$V_0$

$p=2$