Source author record

Jacob Bedrossian

Jacob Bedrossian appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Kinetic shock profiles for the Landau equation

The physical quantities in a gas should vary continuously across a shock. However, the physics inherent in the compressible Euler equations is insufficient to describe the width or structure of the shock. We demonstrate the existence of weak shock profiles to the kinetic Landau equation, that is, traveling wave solutions with Maxwellian asymptotic states whose hydrodynamic quantities satisfy the Rankine-Hugoniot conditions. These solutions serve to capture the structure of weak shocks at the kinetic level. Previous works considered only the Boltzmann equation with hard sphere and angular cut-off potentials.

preprint2022arXiv

A regularity method for lower bounds on the Lyapunov exponent for stochastic differential equations

We put forward a new method for obtaining quantitative lower bounds on the top Lyapunov exponent of stochastic differential equations (SDEs). Our method combines (i) an (apparently new) identity connecting the top Lyapunov exponent to a Fisher information-like functional of the stationary density of the Markov process tracking tangent directions with (ii) a novel, quantitative version of Hörmander's hypoelliptic regularity theory in an $L^1$ framework which estimates this (degenerate) Fisher information from below by a $W^{s,1}_{\mathrm{loc}}$ Sobolev norm. This method is applicable to a wide range of systems beyond the reach of currently existing mathematically rigorous methods. As an initial application, we prove the positivity of the top Lyapunov exponent for a class of weakly-dissipative, weakly forced SDE; in this paper we prove that this class includes the Lorenz 96 model in any dimension, provided the additive stochastic driving is applied to any consecutive pair of modes.

preprint2022arXiv

Chaos in stochastic 2d Galerkin-Navier-Stokes

We prove that all Galerkin truncations of the 2d stochastic Navier-Stokes equations in vorticity form on any rectangular torus subjected to hypoelliptic, additive stochastic forcing are chaotic at sufficiently small viscosity, provided the frequency truncation satisfies $N\geq 392$. By "chaotic" we mean having a strictly positive Lyapunov exponent, i.e. almost-sure asymptotic exponential growth of the derivative with respect to generic initial conditions. A sufficient condition for such results was derived in previous joint work with Alex Blumenthal which reduces the question to the non-degeneracy of a matrix Lie algebra implying Hörmander's condition for the Markov process lifted to the sphere bundle (projective hypoellipticity). The purpose of this work is to reformulate this condition to be more amenable for Galerkin truncations of PDEs and then to verify this condition using a) a reduction to genericity properties of a diagonal sub-algebra inspired by the root space decomposition of semi-simple Lie algebras and b) computational algebraic geometry executed by Maple in exact rational arithmetic. Note that even though we use a computer assisted proof, the result is valid for all aspect ratios and all sufficiently high dimensional truncations; in fact, certain steps simplify in the formal infinite dimensional limit.

preprint2022arXiv

Stationary measures for stochastic differential equations with degenerate damping

A variety of physical phenomena involve the nonlinear transfer of energy from weakly damped modes subjected to external forcing to other modes which are more heavily damped. In this work we explore this in (finite-dimensional) stochastic differential equations in $\mathbb R^n$ with a quadratic, conservative nonlinearity $B(x,x)$ and a linear damping term $-Ax$ which is degenerate in the sense that $\mathrm{ker} A \neq \emptyset$. We investigate sufficient conditions to deduce the existence of a stationary measure for the associated Markov semigroups. Existence of such measures is straightforward if $A$ is full rank, but otherwise, energy could potentially accumulate in $\mathrm{ker} A$ and lead to almost-surely unbounded trajectories, making the existence of stationary measures impossible. We give a relatively simple and general sufficient condition based on time-averaged coercivity estimates along trajectories in neighborhoods of $\mathrm{ker} A$ and many examples where such estimates can be made.

preprint2022arXiv

Vortex filament solutions of the Navier-Stokes equations

We consider solutions of the Navier-Stokes equations in $3d$ with vortex filament initial data of arbitrary circulation, that is, initial vorticity given by a divergence-free vector-valued measure of arbitrary mass supported on a smooth curve. First, we prove global well-posedness for perturbations of the Oseen vortex column in scaling-critical spaces. Second, we prove local well-posedness (in a sense to be made precise) when the filament is a smooth, closed, non-self-intersecting curve. Besides their physical interest, these results are the first to give well-posedness in a neighborhood of large self-similar solutions of $3d$ Navier-Stokes, as well as solutions which are locally approximately self-similar.

preprint2020arXiv

Linearized wave-damping structure of Vlasov-Poisson in $\mathbb R^3$

In this paper we study the linearized Vlasov-Poisson equation for localized disturbances of an infinite, homogeneous Maxwellian background distribution in $\mathbb R^3_x \times \mathbb R^3_v$. In contrast with the confined case $\mathbb T^d _x \times \mathbb R_v ^d$, or the unconfined case $\mathbb R^d_x \times \mathbb R^d_v$ with screening, the dynamics of the disturbance are not scattering towards free transport as $t \to \pm \infty$: we show that the electric field decomposes into a very weakly-damped Klein-Gordon-type evolution for long waves and a Landau-damped evolution. The Klein-Gordon-type waves solve, to leading order, the compressible Euler-Poisson equations linearized about a constant density state, despite the fact that our model is collisionless, i.e. there is no trend to local or global thermalization of the distribution function in strong topologies. We prove dispersive estimates on the Klein-Gordon part of the dynamics. The Landau damping part of the electric field decays faster than free transport at low frequencies and damps as in the confined case at high frequencies; in fact, it decays at the same rate as in the screened case. As such, neither contribution to the electric field behaves as in the vacuum case.

preprint2020arXiv

Suppression of plasma echoes and Landau damping in Sobolev spaces by weak collisions in a Vlasov-Fokker-Planck equation

In this paper, we study Landau damping in the weakly collisional limit of a Vlasov-Fokker-Planck equation with nonlinear collisions in the phase-space $(x,v) \in \mathbb T_x^n \times \mathbb R^n_v$. The goal is four-fold: (A) to understand how collisions suppress plasma echoes and enable Landau damping in agreement with linearized theory in Sobolev spaces, (B) to understand how phase mixing accelerates collisional relaxation, (C) to understand better how the plasma returns to global equilibrium during Landau damping, and (D) to rule out that collision-driven nonlinear instabilities dominate. We give an estimate for the scaling law between Knudsen number and the maximal size of the perturbation necessary for linear theory to be accurate in Sobolev regularity. We conjecture this scaling to be sharp (up to logarithmic corrections) due to potential nonlinear echoes in the collisionless model.

preprint2019arXiv

Sufficient conditions for dual cascade flux laws in the stochastic 2d Navier-Stokes equations

We provide sufficient conditions for mathematically rigorous proofs of the third order universal laws capturing the energy flux to large scales and enstrophy flux to small scales for statistically stationary, forced-dissipated 2d Navier-Stokes equations in the large-box limit. These laws should be regarded as 2d turbulence analogues of the $4/5$ law in 3d turbulence, predicting a constant flux of energy and enstrophy (respectively) through the two inertial ranges in the dual cascade of 2d turbulence. Conditions implying only one of the two cascades are also obtained, as well as compactness criteria which show that the provided sufficient conditions are not far from being necessary. The specific goal of the work is to provide the weakest characterizations of the '0-th laws' of 2d turbulence in order to make mathematically rigorous predictions consistent with experimental evidence.

preprint2018arXiv

The linearized Vlasov and Vlasov-Fokker-Planck equations in a uniform magnetic field

We study the linearized Vlasov equations and the linearized Vlasov-Fokker-Planck equations in the weakly collisional limit in a uniform magnetic field. In both cases, we consider periodic confinement and Maxwellian (or close to Maxwellian) backgrounds. In the collisionless case, for modes transverse to the magnetic field, we provide a precise decomposition into a countably infinite family of standing waves for each spatial mode. These are known as Bernstein modes in the physics literature, though the decomposition is not an obvious consequence of any existing arguments that we are aware of. We show that other modes undergo Landau damping. In the presence of collisions with collision frequency $ν\ll 1$, we show that these modes undergo uniform-in-$ν$ Landau damping and enhanced collisional relaxation at the time-scale $O(ν^{-1/3})$. The modes transverse to the field are uniformly stable and exponentially thermalize on the time-scale $O(ν^{-1})$. Most of the results are proved using Laplace transform analysis of the associated Volterra equations, whereas a simple case of Yan Guo's energy method for hypocoercivity of collision operators is applied for stability in the collisional case.

preprint2016arXiv

Invariant measures for passive scalars in the small noise inviscid limit

We consider a class of invariant measures for a passive scalar $f$ driven by an incompressible velocity field $\boldsymbol{u}$, on a $d$-dimensional periodic domain, satisfying $$ \partial_t f + \boldsymbol{u} \cdot \nabla f = 0, \qquad f(0)=f_0. $$ The measures are obtained as limits of stochastic viscous perturbations. We prove that the span of the $H^1$ eigenfunctions of the operator $\boldsymbol{u}\cdot\nabla$ contains the support of these measures. We also analyze several explicit examples: when $\boldsymbol{u}$ is a shear flow or a relaxation enhancing flow (a generalization of weakly mixing), we can characterize the limiting measure uniquely and compute its covariance structure. We also consider the case of two-dimensional cellular flows, for which further regularity properties of the functions in the support of the measure can be deduced. The main results are proved with the use of spectral theory results, in particular the RAGE theorem, which are used to characterize large classes of orbits of the inviscid problem that are growing in $H^1$.

preprint2016arXiv

Landau damping in finite regularity for unconfined systems with screened interactions

We prove Landau damping for the collisionless Vlasov equation with a class of $L^1$ interaction potentials (including the physical case of screened Coulomb interactions) on $\mathbb R^3_x \times \mathbb R^3_v$ for localized disturbances of an infinite, homogeneous background. Unlike the confined case $\mathbb T^3_x \times \mathbb R_v^3$, results are obtained for initial data in Sobolev spaces (as well as Gevrey and analytic classes). For spatial frequencies bounded away from zero, the Landau damping of the density is similar to the confined case. The finite regularity is possible due to an additional dispersive mechanism available on $\mathbb R_x^3$ which reduces the strength of the plasma echo resonance.

preprint2016arXiv

Suppression of blow-up in Patlak-Keller-Segel via shear flows

In this paper we consider the parabolic-elliptic Patlak-Keller-Segel models in $\mathbb T^d$ with $d=2,3$ with the additional effect of advection by a large shear flow. Without the shear flow, the model is $L^1$ critical in two dimensions with critical mass $8π$; solutions with mass less than $8π$ are global and solutions with mass larger than $8 π$ with finite second moment, all blow up in finite time. In three dimensions, the model is $L^{3/2}$ critical and $L^1$ supercritical; there exists solutions with arbitrarily small mass which blow up in finite time arbitrarily fast. We show that the additional shear flow, if it is chosen sufficiently large, suppresses one dimension of the dynamics and hence can suppress blow-up. In two dimensions, the problem becomes effectively $L^1$ subcritical and so all solutions are global in time (if the shear flow is chosen large). In three dimensions, the problem is effectively $L^1$ critical, and solutions with mass less than $8π$ are global in time (and for all mass larger than $8π$, there exists solutions which blow up in finite time).

preprint2016arXiv

The Sobolev stability threshold for 2D shear flows near Couette

We consider the 2D Navier-Stokes equation on $\mathbb T \times \mathbb R$, with initial datum that is $\varepsilon$-close in $H^N$ to a shear flow $(U(y),0)$, where $\| U(y) - y\|_{H^{N+4}} \ll 1$ and $N>1$. We prove that if $\varepsilon \ll ν^{1/2}$, where $ν$ denotes the inverse Reynolds number, then the solution of the Navier-Stokes equation remains $\varepsilon$-close in $H^1$ to $(e^{t ν\partial_{yy}}U(y),0)$ for all $t>0$. Moreover, the solution converges to a decaying shear flow for times $t \gg ν^{-1/3}$ by a mixing-enhanced dissipation effect, and experiences a transient growth of gradients. In particular, this shows that the stability threshold in finite regularity scales no worse than $ν^{1/2}$ for 2D shear flows close to the Couette flow.

preprint2015arXiv

Dynamics near the subcritical transition of the 3D Couette flow I: Below threshold case

We study small disturbances to the periodic, plane Couette flow in the 3D incompressible Navier-Stokes equations at high Reynolds number $\textbf{Re}$. We prove that for sufficiently regular initial data of size $ε\leq c_0\textbf{Re}^{-1}$ for some universal $c_0 > 0$, the solution is global, remains within $O(c_0)$ of the Couette flow in $L^2$, and returns to the Couette flow as $t \rightarrow \infty$. For times $t \gtrsim \textbf{Re}^{1/3}$, the streamwise dependence is damped by a mixing-enhanced dissipation effect and the solution is rapidly attracted to the class of "2.5 dimensional" streamwise-independent solutions referred to as streaks. Our analysis contains perturbations that experience a transient growth of kinetic energy from $O(\textbf{Re}^{-1})$ to $O(c_0)$ due to the algebraic linear instability known as the lift-up effect. Furthermore, solutions can exhibit a direct cascade of energy to small scales. The behavior is very different from the 2D Couette flow, in which stability is independent of $\textbf{Re}$, enstrophy experiences a direct cascade, and inviscid damping is dominant (resulting in a kind of inverse energy cascade). In 3D, inviscid damping will play a role on one component of the velocity, but the primary stability mechanism is the mixing-enhanced dissipation. Central to the proof is a detailed analysis of the interplay between the stabilizing effects of the mixing and enhanced dissipation and the destabilizing effects of the lift-up effect, vortex stretching, and weakly nonlinear instabilities connected to the non-normal nature of the linearization.

preprint2015arXiv

Dynamics near the subcritical transition of the 3D Couette flow II: Above threshold case

This is the second in a pair of works which study small disturbances to the plane, periodic 3D Couette flow in the incompressible Navier-Stokes equations at high Reynolds number $\textbf{Re}$. In this work, we show that there is constant $0 < c_0 \ll 1$, independent of $\textbf{Re}$, such that sufficiently regular disturbances of size $ε\lesssim \textbf{Re}^{-2/3-δ}$ for any $δ> 0$ exist at least until $t = c_0ε^{-1}$ and in general evolve to be $O(c_0)$ due to the lift-up effect. Further, after times $t \gtrsim \textbf{Re}^{1/3}$, the streamwise dependence of the solution is rapidly diminished by a mixing-enhanced dissipation effect and the solution is attracted back to the class of "2.5 dimensional" streamwise-independent solutions (sometimes referred to as "streaks"). The largest of these streaks are expected to eventually undergo a secondary instability at $t \approx ε^{-1}$. Hence, our work strongly suggests, for all (sufficiently regular) initial data, the genericity of the "lift-up effect $\Rightarrow$ streak growth $\Rightarrow$ streak breakdown" scenario for turbulent transition of the 3D Couette flow near the threshold of stability forwarded in the applied mathematics and physics literature.

preprint2015arXiv

On the stability threshold for the 3D Couette flow in Sobolev regularity

We study Sobolev regularity disturbances to the periodic, plane Couette flow in the 3D incompressible Navier-Stokes equations at high Reynolds number $\textbf{Re}$. Our goal is to estimate how the stability threshold scales in $\textbf{Re}$: the largest the initial perturbation can be while still resulting in a solution that does not transition away from Couette flow. In this work we prove that initial data which satisfies $\| u_{in} \|_{H^σ} \leq δ\textbf{Re}^{-3/2}$ for any $σ> 9/2$ and some $δ= δ(σ) > 0$ depending only on $σ$, is global in time, remains within $O(\textbf{Re}^{-1/2})$ of the Couette flow in $L^2$ for all time, and converges to the class of "2.5 dimensional" streamwise-independent solutions referred to as streaks for times $t \gtrsim \textbf{Re}^{1/3}$. Numerical experiments performed by Reddy et. al. with "rough" initial data estimated a threshold of $\sim \textbf{Re}^{-31/20}$, which shows very close agreement with our estimate.

preprint2014arXiv

Enhanced dissipation and inviscid damping in the inviscid limit of the Navier-Stokes equations near the 2D Couette flow

In this work we study the long time, inviscid limit of the 2D Navier-Stokes equations near the periodic Couette flow, and in particular, we confirm at the nonlinear level the qualitative behavior predicted by Kelvin's 1887 linear analysis. At high Reynolds number Re, we prove that the solution behaves qualitatively like 2D Euler for times t \lesssim Re^(1/3), and in particular exhibits inviscid damping (e.g. the vorticity weakly approaches a shear flow). For times t \gtrsim Re^(1/3), which is sooner than the natural dissipative time scale O(Re), the viscosity becomes dominant and the streamwise dependence of the vorticity is rapidly eliminated by an enhanced dissipation effect. Afterward, the remaining shear flow decays on very long time scales t \gtrsim Re back to the Couette flow. When properly defined, the dissipative length-scale in this setting is L_D \sim Re^(-1/3), larger than the scale L_D \sim Re^(-1/2) predicted in classical Batchelor-Kraichnan 2D turbulence theory. The class of initial data we study is the sum of a sufficiently smooth function and a small (with respect to Re^(-1)) $L^2$ function.

preprint2014arXiv

Inviscid damping and the asymptotic stability of planar shear flows in the 2D Euler equations

We prove asymptotic stability of shear flows close to the planar Couette flow in the 2D inviscid Euler equations on $\Torus \times \Real$. That is, given an initial perturbation of the Couette flow small in a suitable regularity class, specifically Gevrey space of class smaller than 2, the velocity converges strongly in L^2 to a shear flow which is also close to the Couette flow. The vorticity is asymptotically driven to small scales by a linear evolution and weakly converges as $t \rightarrow \pm\infty$. The strong convergence of the velocity field is sometimes referred to as inviscid damping, due to the relationship with Landau damping in the Vlasov equations. This convergence was formally derived at the linear level by Kelvin in 1887 and it occurs at an algebraic rate first computed by Orr in 1907; our work appears to be the first rigorous confirmation of this behavior on the nonlinear level.

preprint2014arXiv

Large mass global solutions for a class of L1-critical nonlocal aggregation equations and parabolic-elliptic Patlak-Keller-Segel models

We consider a class of $L^1$ critical nonlocal aggregation equations with linear or nonlinear porous media-type diffusion which are characterized by a long-range interaction potential that decays faster than the Newtonian potential at infinity. The fast decay breaks the $L^1$ scaling symmetry and we prove that `sufficiently spread out' initial data, regardless of the mass, result in global spreading solutions. This is in contrast to the classical parabolic-elliptic PKS for which essentially all solutions with more than critical mass are known to blow up in finite time. In all cases, the long-time asymptotics are given by the self-similar solution to the linear heat equation or by the Barenblatt solutions of the porous media equation. The results with linear diffusion are proved using properties of the Fokker-Planck semi-group whereas the results with nonlinear diffusion are proved using a more interesting bootstrap argument coupling the entropy-entropy dissipation methods of the porous media equation together with higher $L^p$ estimates similar to those used in small-data and local theory for PKS-type equations.

preprint2013arXiv

Asymptotic stability for the Couette flow in the 2D Euler equations

In this expository note we discuss our recent work [arXiv:1306.5028] on the nonlinear asymptotic stability of shear flows in the 2D Euler equations of ideal, incompressible flow. In that work it is proved that perturbations to the Couette flow which are small in a suitable regularity class converge strongly in $L^2$ to a shear flow which is close to the Couette flow. Enstrophy is mixed to small scales by an almost linear evolution and is generally lost in the weak limit as t -> +/- infinity. In this note we discuss the most important physical and mathematical aspects of the result and the key ideas of the proof.

preprint2013arXiv

Blister patterns and energy minimization in compressed thin films on compliant substrates

This paper is motivated by the complex blister patterns sometimes seen in thin elastic films on thick, compliant substrates. These patterns are often induced by an elastic misfit which compresses the film. Blistering permits the film to expand locally, reducing the elastic energy of the system. It is natural to ask: what is the minimum elastic energy achievable by blistering on a fixed area fraction of the substrate? This is a variational problem involving both the {\it elastic deformation} of the film and substrate and the {\it geometry} of the blistered region. It involves three small parameters: the {\it nondimensionalized thickness} of the film, the {\it compliance ratio} of the film/substrate pair and the {\it mismatch strain}. In formulating the problem, we use a small-slope (Föppl-von Kármán) approximation for the elastic energy of the film, and a local approximation for the elastic energy of the substrate. For a 1D version of the problem, we obtain "matching" upper and lower bounds on the minimum energy, in the sense that both bounds have the same scaling behavior with respect to the small parameters. For a 2D version of the problem, our results are less complete. Our upper and lower bounds only "match" in their scaling with respect to the nondimensionalized thickness, not in the dependence on the compliance ratio and the mismatch strain. The upper bound considers a 2D lattice of blisters, and uses ideas from the literature on the folding or "crumpling" of a confined elastic sheet. Our main 2D result is that in a certain parameter regime, the elastic energy of this lattice is significantly lower than that of a few large blisters.

preprint2013arXiv

Existence, Uniqueness and Lipschitz Dependence for Patlak-Keller-Segel and Navier-Stokes in R2 with Measure-valued Initial Data

We establish a new local well-posedness result in the space of finite Borel measures for mild solutions of the parabolic-elliptic Patlak-Keller-Segel (PKS) model of chemotactic aggregation in two dimensions. Our result only requires that the initial measure satisfy the necessary assumption $\max_{x \in \Real^2}μ(\set{x}) < 8π$. This work improves the small-data results of Biler and the existence results of Senba and Suzuki. Our work is based on that of Gallagher and Gallay, who prove the uniqueness and log-Lipschitz continuity of the solution map for the 2D Navier-Stokes equations (NSE) with measure-valued initial vorticity. We refine their techniques and present an alternative version of their proof which yields existence, uniqueness and \emph{Lipschitz} continuity of the solution maps of both PKS and NSE. Many steps are more difficult for PKS than for NSE, particularly on the level of the linear estimates related to the self-similar spreading solutions.

preprint2013arXiv

Landau damping: paraproducts and Gevrey regularity

We give a new, simpler, proof of nonlinear Landau damping on T^d in Gevrey-1/s regularity (s > 1/3) which matches the regularity requirement predicted by the formal analysis of Mouhot and Villani in the original proof of Landau damping [Acta Mathematica 2011]. Our proof combines in a novel way ideas from the original proof of Landau damping and the proof of inviscid damping in 2D Euler [arXiv:1306.5028]. As in the work on 2D Euler, we use paraproduct decompositions and controlled regularity loss to replace the Newton iteration scheme employed in the original proof. We perform time-response estimates adapted from the original proof to control the plasma echoes and couple them to energy estimates on the distribution function in the style of the work on 2D Euler.

preprint2012arXiv

Bounded Mean Oscillation and the Uniqueness of Active Scalar Equations

We consider a number of uniqueness questions for several wide classes of active scalar equations, unifying and generalizing the techniques of several authors. As special cases of our results, we provide a significantly simplified proof to the known uniqueness result for the 2D Euler equations in $L^1 \cap BMO$ and provide a mild improvement to the recent results of Rusin for the 2D inviscid surface quasi-geostrophic (SQG) equations, which are now to our knowledge, the best results known for this model. We also obtain what are (to our knowledge) the strongest known uniqueness results for the Patlak-Keller-Segel models. We obtain these results via technical refinements of energy methods which are well-known in the $L^2$ setting but are less well-known in the $\dot{H}^{-1}$ setting. The $\dot{H}^{-1}$ method can be considered a generalization of Yudovich's classical method and is naturally applied to equations such as the Patlak-Keller-Segel models with nonlinear diffusion, and other variants. Important points of our analysis are an $L^p$-$BMO$ interpolation lemma and a Sobolev embedding lemma which shows that velocity fields $v$ with $\grad v \in BMO$ are locally log-Lipschitz; the latter is known in harmonic analysis but does not seem to have been connected to this setting.

preprint2012arXiv

Global Existence and Finite Time Blow-Up for Critical Patlak-Keller-Segel Models with Inhomogeneous Diffusion

The L^1-critical parabolic-elliptic Patlak-Keller-Segel system is a classical model of chemotactic aggregation in micro-organisms well-known to have critical mass phenomena. In this paper we study this critical mass phenomenon in the context of Patlak-Keller-Segel models with spatially varying diffusivity and decay rate of the chemo-attractant. The primary tool for the proof of global existence below the critical mass is the use of pseudo-differential operators to precisely evaluate the leading order quadratic portion of the potential energy (interaction energy). Under the assumption of radial symmetry, blow-up is proved above critical mass using a maximum-principle type argument based on comparing the mass distribution of solutions to a barrier consisting of the unique stationary solutions of the scale-invariant PKS. Although effective where standard Virial methods do not apply, this method seems to be dependent on the assumption of radial symmetry. For technical reasons we work in dimensions three and higher where L^1-critical variants of the PKS have porous media-type nonlinear diffusion on the organism density, resulting in finite speed of propagation and simplified functional inequalities.

preprint2011arXiv

Inhomogeneous Patlak-Keller-Segel models and Aggregation Equations with Nonlinear Diffusion in $\Real^d$

Aggregation equations and Patlak-Keller-Segel (PKS) models for chemotaxis with nonlinear diffusion are popular models for nonlocal aggregation phenomenon and are a source of a number of interesting mathematical problems in nonlinear PDE. The purpose of this work is twofold. First, we continue our previous work, which focused on nonlocal aggregation, modeled with a convolution. The goal was to unify the local and global theory of these convolution-type models, including the identification of a sharp critical mass; however, some cases involving unbounded domains were left open. In particular, the biologically relevant case $\Real^2$ was not treated. In this paper, we present an alternative proof of local existence, which now applies to $\Real^d$ for all $d \geq 2$ and give global results that were left open. The proof departs from previous work in that it uses a more direct and intuitive regularization that constructs approximate solutions on $\Real^d$ instead of on sequences of bounded domains. Second, this work develops the local, subcritical, and small data critical theory for a variety of Patlak-Keller-Segel models with spatially varying diffusion and decay rate of the chemo-attractant.

preprint2011arXiv

Intermediate asymptotics for critical and supercritical aggregation equations and Patlak-Keller-Segel models

We examine the long-term asymptotic behavior of dissipating solutions to aggregation equations and Patlak-Keller-Segel models with degenerate power-law and linear diffusion. The purpose of this work is to identify when solutions decay to the self-similar spreading solutions of the homogeneous diffusion equations. Combined with strong decay estimates, entropy-entropy dissipation methods provide a natural solution to this question and make it possible to derive quantitative convergence rates in $L^1$. The estimated rate depends only on the nonlinearity of the diffusion and the strength of the interaction kernel at long range.

preprint2011arXiv

Local and Global Well-Posedness for Aggregation Equations and Patlak-Keller-Segel Models with Degenerate Diffusion

Recently, there has been a wide interest in the study of aggregation equations and Patlak-Keller-Segel (PKS) models for chemotaxis with degenerate diffusion. The focus of this paper is the unification and generalization of the well-posedness theory of these models. We prove local well-posedness on bounded domains for dimensions $d\geq 2$ and in all of space for $d\geq 3$, the uniqueness being a result previously not known for PKS with degenerate diffusion. We generalize the notion of criticality for PKS and show that subcritical problems are globally well-posed. For a fairly general class of problems, we prove the existence of a critical mass which sharply divides the possibility of finite time blow up and global existence. Moreover, we compute the critical mass for fully general problems and show that solutions with smaller mass exists globally. For a class of supercritical problems we prove finite time blow up is possible for initial data of arbitrary mass.

preprint2010arXiv

Global Minimizers for Free Energies of Subcritical Aggregation Equations with Degenerate Diffusion

We prove the existence of non-trivial global minimizers of a class of free energies related to aggregation equations with degenerate diffusion on $\Real^d$. Such equations arise in mathematical biology as models for organism group dynamics which account for competition between the tendency to aggregate into groups and nonlinear diffusion to avoid over-crowding. The existence of non-zero optimal free energy stationary solutions representing coherent groups in $\Real^d$ is therefore of interest. The primary contribution is the investigation of a notion of criticality associated with the global minimizer problem. The notion arises from the scaling of diffusive and aggregative forces as mass spreads and is shown to dictate the existence, and sometimes non-existence, of global minimizers.