We consider radially symmetric solutions of the degenerate Keller–Segel system

\begin{align*}\begin{cases}\partial_t u=\nabla\cdot (u^{m-1}\nabla u - u\nabla v),\\0=\Delta v -\mu +u,\quad\mu =\frac{1}{|\Omega|}\int_\Omega u,\end{cases}\end{align*}

in balls $\Omega\subset\mathbb R^n$, $n\ge 1$, where m > 1 is arbitrary. Our main result states that the initial evolution of the positivity set of u is essentially determined by the shape of the (nonnegative, radially symmetric, Hölder continuous) initial data u0 near the boundary of its support $\overline{B_{r_1}(0)}\subsetneq\Omega$: It shrinks for sufficiently flat and expands for sufficiently steep u0. More precisely, there exists an explicit constant $A_{\mathrm{crit}} \in (0, \infty)$ (depending only on $m, n, R, r_1$ and $\int_\Omega u_0$) such that if $u_0(x)\le A(r_1-|x|)^\frac{1}{m-1}$ for all $|x|\in(r_0, r_1)$ and some $r_0\in(0,r_1)$ and $A \lt A_{\mathrm{crit}}$ then there are T > 0 and ζ > 0 such that $\sup\{\, |x| \mid x \in \operatorname{supp} u(\cdot, t)\,\}\le r_1 -\zeta t$ for all $t\in(0, T)$, while if $u_0(x)\ge A(r_1-|x|)^\frac{1}{m-1}$ for all $|x|\in(r_0, r_1)$ and some $r_0 \in (0, r_1)$ and $A \gt A_{\mathrm{crit}}$ then we can find T > 0 and ζ > 0 such that $\sup\{\, |x| \mid x \in \operatorname{supp} u(\cdot, t)\,\}\ge r_1 +\zeta t$ for all $t\in(0, T)$.

In this paper, we consider the following PDE-ODE system modelling cancer invasion with slow diffusion and ECM remodelling,

\[ \begin{cases} u_t=\Delta u^m-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), \\ v_t=\Delta v+u-v, \\ \omega_t={-}v\omega+\eta \omega(1-u-\omega). \end{cases} \]

$\eta =0$

$\eta >0$

$\eta =0$

$\eta >0$

Existence of specific eternal solutions in exponential self-similar form to the following quasilinear diffusion equation with strong absorption

\[ \partial_t u=\Delta u^m-|x|^{\sigma}u^q, \]

$(t,\,x)\in (0,\,\infty )\times \mathbb {R}^N$

$m>1$

$q\in (0,\,1)$

$\sigma =\sigma _c:=2(1-q)/ (m-1)$

\[ u(t,x)={\rm e}^{-\alpha t}f(|x|\,{\rm e}^{\beta t}), \quad \alpha=\frac{2}{m-1}\beta, \]

$\beta ^*\in (0,\,\infty )$

$(u_A)_{A>0}$

$(f_A)_{A>0}$

$f_A(0)=A$

$f_A'(0)=0$

$\sigma \in (0,\,\sigma _c)$

In this paper, we mainly investigate the well-posedness of the four-order degenerate differential equation ($P_4$): $(Mu)''''(t) + \alpha (Lu)'''(t) + (Lu)''(t)$ $=\beta Au(t) + \gamma Bu'(t) + Gu'_t + Fu_t + f(t),\,( t\in [0,\,2\pi ])$ in periodic Lebesgue–Bochner spaces $L^p(\mathbb {T}; X)$ and periodic Besov spaces $B_{p,q}^s\;(\mathbb {T}; X)$, where $A$, $B$, $L$ and $M$ are closed linear operators on a Banach space $X$ such that $D(A)\cap D(B)\subset D(M)\cap D(L)$ and $\alpha,\,\beta,\,\gamma \in \mathbb {C}$, $G$ and $F$ are bounded linear operators from $L^p([-2\pi,\,0];X)$ (respectively $B_{p,q}^s([-2\pi,\,0];X)$) into $X$, $u_t(\cdot ) = u(t+\cdot )$ and $u'_t(\cdot ) = u'(t+\cdot )$ are defined on $[-2\pi,\,0]$ for $t\in [0,\, 2\pi ]$. We completely characterize the well-posedness of ($P_4$) in the above two function spaces by using known operator-valued Fourier multiplier theorems.

In a ball $\Omega \subset \mathbb {R}^{n}$ with $n\ge 2$, the chemotaxis system

\[ \left\{ \begin{array}{@{}l} u_t = \nabla \cdot \big( D(u)\nabla u\big) + \nabla\cdot \big(\dfrac{u}{v} \nabla v\big), \\ 0=\Delta v - uv \end{array} \right. \]

$u$

$v$

$D\in C^{3}([0,\infty ))$

$0< D(\xi ) \le {K_D} (\xi +1)^{-\alpha }$

$\xi >0$

${K_D}>0$

$\alpha >0$

$\Omega$

Structural changes of the pore space and clogging phenomena are inherent to many porous media applications. However, related analytical investigations remain challenging due to potentially vanishing coefficients in the respective systems of partial differential equations. In this research, we apply an appropriate scaling of the unknowns and work with porosity-weighted function spaces. This enables us to prove existence, uniqueness and non-negativity of weak solutions to a combined flow and transport problem with vanishing, but prescribed porosity field, permeability and diffusion.

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.

In this article, we consider diffusive transport of a reactive substance in a saturated porous medium including variable porosity. Thereby, the evolution of the microstructure is caused by precipitation of the transported substance. We are particularly interested in analysing the model when the equations degenerate due to clogging. Introducing an appropriate weighted function space, we are able to handle the degeneracy and obtain analytical results for the transport equation. Also the decay behaviour of this solution with respect to the porosity is investigated. There a restriction on the decay order is assumed, that is, besides low initial concentration also dense precipitation leads to possible high decay. We obtain nonnegativity and boundedness for the weak solution to the transport equation. Moreover, we study an ordinary differential equation (ODE) describing the change of porosity. Thereby, the control of an appropriate weighted norm of the gradient of the porosity is crucial for the analysis of the transport equation. In order to obtain global in time solutions to the overall coupled system, we apply a fixed point argument. The problem is solved for substantially degenerating hydrodynamic parameters.

We prove the existence of global minimisers for a class of attractive–repulsive interaction potentials that are in general not radially symmetric. The global minimisers have compact support. For potentials including degenerate power-law diffusion, the interaction potential can be unbounded from below. Further, a formal calculation indicates that for non-symmetric potentials global minimisers may neither be radial symmetric nor unique.

Let n ⩾ 3 and 0 < m < (n − 2)/n. We extend the results of Vazquez and Winkler (2011, J. Evol. Equ. 11, no. 3, 725–742) and prove the uniqueness of finite points blow-up solutions of the fast diffusion equation ut = Δum in both bounded domains and ℝn × (0, ∞). We also construct initial data such that the corresponding solution of the fast diffusion equation in bounded domain oscillates between infinity and some positive constant as t → ∞.

We present a full classification of the short-time behaviour of the interfaces and local solutions to the nonlinear parabolic p-Laplacian-type reaction-diffusion equation of non-Newtonian elastic filtration

$$u_t-\Big(|u_x|^{p-2}u_x\Big)_x+bu^{\beta}=0, \ 1 \lt p \lt 2, \beta \gt 0.$$

This paper deals with nonnegative solutions of the one-dimensional degenerate parabolic equations with zero homogeneous Dirichlet boundary condition. To obtain an existence result, we prove a sharp estimate for |ux|. Besides, we investigate the qualitative behaviours of nonnegative solutions such as the quenching phenomenon, and the finite speed of propagation. Our results of the Dirichlet problem are also extended to the associated Cauchy problem on the whole domain ℝ. In addition, we also consider the instantaneous shrinking of compact support of nonnegative solutions.

The Wasserstein gradient flow structure of the partial differential equation system governing multiphase flows in porous media was recently highlighted in Cancès et al. [Anal. PDE10(8), 1845–1876]. The model can thus be approximated by means of the minimising movement (or JKO after Jordan, Kinderlehrer and Otto [SIAM J. Math. Anal.29(1), 1–17]) scheme that we solve thanks to the ALG2-JKO scheme proposed in Benamou et al. [ESAIM Proc. Surv.57, 1–17]. The numerical results are compared to a classical upstream mobility finite volume scheme, for which strong stability properties can be established.

We propose a numerical method for the simulation of a quasi-linear parabolic biofilm model that exhibits three non-linear diffusion effects: (i) a power law degeneracy, (ii) a super diffusion singularity and (iii) non-linear cross-diffusion. The method is based on a spatial Finite Volume discretisation in which cross-diffusion terms are formally treated as convection terms. Time-integration of the resulting semi-discretised system is carried out using an error-controlled, time-adaptive, embedded Rosenbrock–Wanner method. We compare several variants of the method and two variants of the model to investigate how details such as the choice cross-diffusion coefficients, and specific variants of the time integrator affect simulation time.

We consider a scalar conservation law with zero-flux boundary conditions imposed on the boundary of a rectangular multidimensional domain. We study monotone schemes applied to this problem. For the Godunov version of the scheme, we simply set the boundary flux equal to zero. For other monotone schemes, we additionally apply a simple modification to the numerical flux. We show that the approximate solutions produced by these schemes converge to the unique entropy solution, in the sense of [7], of the conservation law. Our convergence result relies on a BV bound on the approximate numerical solution. In addition, we show that a certain functional that is closely related to the total variation is nonincreasing from one time level to the next. We extend our scheme to handle degenerate convection-diffusion equations and for the one-dimensional case we prove convergence to the unique entropy solution.

We study a system of two coupled parabolic equations that models the spreading of a drop of an insoluble surfactant on a thin liquid film. Unlike the previously known results, the surface diffusion coefficient is not assumed constant and depends on the surfactant concentration. We obtain sufficient conditions for finite speed of support propagation and for waiting-time phenomenon by application of an extension of Stampacchia's lemma for a system of functional equations.

We study the asymptotic behaviour as t → T– , near a finite blow-up time T > 0, of decreasing-in-x solutions to the following semilinear heat equation with a non-local term:

with Neumann boundary conditions and strictly decreasing initial function u0(x) with zero mass. We prove sharp estimates for u(x, t) as t → T– , revealing a non-uniform global blow-up:

uniformly on any compact set [δ, 1], δ ∈ (0, 1).

Multiclass Lighthill-Whitham-Richards traffic models [Benzoni-Gavage and Colombo, Euro. J. Appl. Math., 14 (2003), pp. 587–612; Wong and Wong, Transp. Res. A, 36 (2002), pp. 827–841] give rise to first-order systems of conservation laws that are hyperbolic under usual conditions, so that their associated Cauchy problems are well-posed. Anticipation lengths and reaction times can be incorporated into these models by adding certain conservative second-order terms to these first-order conservation laws. These terms can be diffusive under certain circumstances, thus, in principle, ensuring the stability of the solutions. The purpose of this paper is to analyze the stability of these diffusively corrected models under varying reaction times and anticipation lengths. It is demonstrated that instabilities may develop for high reaction times and short anticipation lengths, and that these instabilities may have controlled frequencies and amplitudes due to their nonlinear nature.

Finite difference computations that involve spatial adaptation commonly employ an equidistribution principle. In these cases, a new mesh is constructed such that a given monitor function is equidistributed in some sense. Typical choices of the monitor function involve the solution or one of its many derivatives. This straightforward concept has proven to be extremely effective and practical. However, selections of core monitoring functions are often challenging and crucial to the computational success. This paper concerns six different designs of the monitoring function that targets a highly nonlinear partial differential equation that exhibits both quenching-type and degeneracy singularities. While the first four monitoring strategies are within the so-called primitive regime, the rest belong to a later category of the modified type, which requires the priori knowledge of certain important quenching solution characteristics. Simulated examples are given to illustrate our study and conclusions.

Considered here is the pullback attractor of the process associated with the first initial boundary value problem for the non-autonomous semilinear degenerate parabolic equation

\begin{linenomath}u_t-\text{div}(\sigma(x)\nabla u)+f(u)=g(x,t)\end{linenomath}

$D_0^2(\Omega,\sigma)$