Source author record

Robert L. Pego

Robert L. Pego 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

18works
12topics
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

18 published item(s)

preprint2022arXiv

Global Dynamics and Photon Loss in the Kompaneets Equation

The Kompaneets equation governs dynamics of the photon energy spectrum in certain high temperature (or low density) plasmas. We prove several results concerning the long-time convergence of solutions to Bose--Einstein equilibria and the failure of photon conservation. In particular, we show the total photon number can decrease with time via an outflux of photons at the zero-energy boundary. The ensuing accumulation of photons at zero energy is analogous to Bose--Einstein condensation. We provide two conditions that guarantee that photon loss occurs, and show that once loss is initiated then it persists forever. We prove that as $t\to \infty$, solutions necessarily converge to equilibrium and we characterize the limit in terms of the total photon loss. Additionally, we provide a few results concerning the behavior of the solution near the zero-energy boundary, an Oleinik inequality, a comparison principle, and show that the solution operator is a contraction in $L^1$. None of these results impose a boundary condition at the zero-energy boundary.

preprint2019arXiv

A Markov jump process modelling animal group size statistics

We translate a coagulation-framentation model, describing the dynamics of animal group size distributions, into a model for the population distribution and associate the \blue{nonlinear} evolution equation with a Markov jump process of a type introduced in classic work of H.~McKean. In particular this formalizes a model suggested by H.-S. Niwa [J.~Theo.~Biol.~224 (2003)] with simple coagulation and fragmentation rates. Based on the jump process, we develop a numerical scheme that allows us to approximate the equilibrium for the Niwa model, validated by comparison to analytical results by Degond et al. [J.~Nonlinear Sci.~27 (2017)], and study the population and size distributions for more complicated rates. Furthermore, the simulations are used to describe statistical properties of the underlying jump process. We additionally discuss the relation of the jump process to models expressed in stochastic differential equations and demonstrate that such a connection is justified in the case of nearest-neighbour interactions, as opposed to global interactions as in the Niwa model.

preprint2019arXiv

Temporal oscillations in Becker-Doering equations with atomization

We prove that time-periodic solutions arise via Hopf bifurcation in a finite closed system of coagulation-fragmentation equations. The system we treat is a variant of the Becker-Doering equations, in which clusters grow or shrink by addition or deletion of monomers. To this is added a linear atomization reaction for clusters of maximum size. The structure of the system is motivated by models of gas evolution oscillators in physical chemistry, which exhibit temporal oscillations under certain input/output conditions.

preprint2019arXiv

Well-posedness and derivative blow-up for a dispersionless regularized shallow water system

We study local-time well-posedness and breakdown for solutions of regularized Saint-Venant equations (regularized classical shallow water equations) recently introduced by Clamond and Dutykh. The system is linearly non-dispersive, and smooth solutions conserve an $H^1$-equivalent energy. No shock discontinuities can occur, but the system is known to admit weakly singular shock-profile solutions that dissipate energy. We identify a class of small-energy smooth solutions that develop singularities in the first derivatives in finite time.

preprint2016arXiv

Bose-Einstein condensation in a hyperbolic model for the Kompaneets equation

In low-density or high-temperature plasmas, Compton scattering is the dominant process responsible for energy transport. Kompaneets in 1957 derived a non-linear degenerate parabolic equation for the photon energy distribution. In this paper we consider a simplified model obtained by neglecting diffusion of the photon number density in a particular way. We obtain a non-linear hyperbolic PDE with a position-dependent flux, which permits a one-parameter family of stationary entropy solutions to exist. We completely describe the long-time dynamics of each non-zero solution, showing that it approaches some non-zero stationary solution. While the total number of photons is formally conserved, if initially large enough it necessarily decreases after finite time through an out-flux of photons with zero energy. This corresponds to formation of a Bose-Einstein condensate, whose mass we show can only increase with time.

preprint2015arXiv

Coagulation-fragmentation model for animal group-size statistics

We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.

preprint2015arXiv

Limit theorems for Smoluchowski dynamics associated with critical continuous-state branching processes

We investigate the well-posedness and asymptotic self-similarity of solutions to a generalized Smoluchowski coagulation equation recently introduced by Bertoin and Le Gall in the context of continuous-state branching theory. In particular, this equation governs the evolution of the Lévy measure of a critical continuous-state branching process which becomes extinct (i.e., is absorbed at zero) almost surely. We show that a nondegenerate scaling limit of the Lévy measure (and the process) exists if and only if the branching mechanism is regularly varying at 0. When the branching mechanism is regularly varying, we characterize nondegenerate scaling limits of arbitrary finite-measure solutions in terms of generalized Mittag-Leffler series.

preprint2015arXiv

Polynomial decay to equilibrium for the Becker-Döring equations

This paper studies rates of decay to equilibrium for the Becker-Döring equations with subcritical initial data. In particular, polynomial rates of decay are established when initial perturbations of equilibrium have polynomial moments. This is proved by using new dissipation estimates in polynomially weighted $\ell^1$ spaces, operator decomposition techniques from kinetic theory, and interpolation estimates from the study of travelling waves.

preprint2013arXiv

Global existence for two extended Navier-Stokes systems

We prove global existence of weak solutions to two systems of equations which extend the dynamics of the Navier-Stokes equations for incompressible viscous flow with no-slip boundary condition. The systems of equations we consider arise as formal limits of time discrete pressure-Poisson schemes introduced by Johnston and Liu (J. Comp. Phys. 199 (2004) 221-259) and by Shirokoff and Rosales (J Comp. Phys 230 (2011) 8619-8646) when the initial data does not satisfy the required compatibility condition. Unlike the results of Iyer et al (J. Math. Phys. 53 (2012) 115605), our approach proves existence of weak solutions in domains with less than C^1 regularity. Our approach also addresses uniqueness in 2D and higher regularity.

preprint2012arXiv

Asymptotic stability of solitary waves in the Benney-Luke model of water waves

We study asymptotic stability of solitary wave solutions in the one-dimensional Benney-Luke equation, a formally valid approximation for describing two-way water wave propagation. For this equation, as for the full water wave problem, the classic variational method for proving orbital stability of solitary waves fails dramatically due to the fact that the second variation of the energy-momentum functional is infinitely indefinite. We establish nonlinear stability in energy norm under the spectral stability hypothesis that the linearization admits no non-zero eigenvalues of non-negative real part. We also verify this hypothesis for waves of small energy.

preprint2012arXiv

Coercivity and stability results for an extended Navier-Stokes system

In this article we study a system of equations that is known to {\em extend} Navier-Stokes dynamics in a well-posed manner to velocity fields that are not necessarily divergence-free. Our aim is to contribute to an understanding of the role of divergence and pressure in developing energy estimates capable of controlling the nonlinear terms. We address questions of global existence and stability in bounded domains with no-slip boundary conditions. Even in two space dimensions, global existence is open in general, and remains so, primarily due to the lack of a self-contained $L^2$ energy estimate. However, through use of new $H^1$ coercivity estimates for the linear equations, we establish a number of global existence and stability results, including results for small divergence and a time-discrete scheme. We also prove global existence in 2D for any initial data, provided sufficient divergence damping is included.

preprint2011arXiv

Spectral stability of vortices in two-dimensional Bose-Einstein condensates via the Evans function and Krein signature

We investigate spectral stability of vortex solutions of the Gross-Pitaevskii equation, a mean-field approximation for Bose-Einstein condensates (BEC) in an effectively two-dimensional axisymmetric harmonic trap. We study eigenvalues of the linearization both rigorously and through computation of the Evans function, a sensitive and robust technique whose use we justify mathematically. The absence of unstable eigenvalues is justified a posteriori through use of the Krein signature of purely imaginary eigenvalues, which also can be used to significantly reduce computational effort. In particular, we prove general basic continuation results on Krein signature for finite systems of eigenvalues in infinite-dimensional problems.

preprint2010arXiv

Asymptotic linear stability of solitary water waves

We prove an asymptotic stability result for the water wave equations linearized around small solitary waves. The equations we consider govern irrotational flow of a fluid with constant density bounded below by a rigid horizontal bottom and above by a free surface under the influence of gravity neglecting surface tension. For sufficiently small amplitude waves, with waveform well-approximated by the well-known sech-squared shape of the KdV soliton, solutions of the linearized equations decay at an exponential rate in an energy norm with exponential weight translated with the wave profile. This holds for all solutions with no component in (i.e., symplectically orthogonal to) the two-dimensional neutral-mode space arising from infinitesimal translational and wave-speed variation of solitary waves. We also obtain spectral stability in an unweighted energy norm.

preprint2010arXiv

On Optimal Estimates for the Laplace-Leray Commutator in Planar Domains with Corners

For smooth domains, Liu et al. (Comm. Pure Appl. Math. 60: 1443-1487, 2007) used optimal estimates for the commutator of the Laplacian and the Leray projection operator to establish well-posedness of an extended Navier-Stokes dynamics. In their work, the pressure is not determined by incompressibility, but rather by a certain formula involving the Laplace-Leray commutator. A key estimate of Liu et al. controls the commutator strictly by the Laplacian in energy norm at leading order. In this paper we show that this strict control fails in a large family of bounded planar domains with corners. However, when the domain is an infinite cone, we find that strict control may be recovered in certain power-law weighted norms.

preprint2008arXiv

Dynamics and self-similarity in min-driven clustering

We study a mean-field model for a clustering process that may be described informally as follows. At each step a random integer $k$ is chosen with probability $p_k$, and the smallest cluster merges with $k$ randomly chosen clusters. We prove that the model determines a continuous dynamical system on the space of probability measures supported in $(0,\infty)$, and we establish necessary and sufficient conditions for approach to self-similar form. We also characterize eternal solutions for this model via a Levy-Khintchine formula. The analysis is based on an explicit solution formula discovered by Gallay and Mielke, extended using a careful choice of time scale.

preprint2006arXiv

The scaling attractor and ultimate dynamics for Smoluchowski's coagulation equations

We describe a basic framework for studying dynamic scaling that has roots in dynamical systems and probability theory. Within this framework, we study Smoluchowski's coagulation equation for the three simplest rate kernels $K(x,y)=2$, $x+y$ and $xy$. In another work, we classified all self-similar solutions and all universality classes (domains of attraction) for scaling limits under weak convergence (Comm. Pure Appl. Math 57 (2004)1197-1232). Here we add to this a complete description of the set of all limit points of solutions modulo scaling (the scaling attractor) and the dynamics on this limit set (the ultimate dynamics). The main tool is Bertoin's Lévy-Khintchine representation formula for eternal solutions of Smoluchowski's equation (Adv. Appl. Prob. 12 (2002) 547--64). This representation linearizes the dynamics on the scaling attractor, revealing these dynamics to be conjugate to a continuous dilation, and chaotic in a classical sense. Furthermore, our study of scaling limits explains how Smoluchowski dynamics ``compactifies'' in a natural way that accounts for clusters of zero and infinite size (dust and gel).

preprint2006arXiv

Universality classes in Burgers turbulence

We establish necessary and sufficient conditions for the shock statistics to approach self-similar form in Burgers turbulence with Lévy process initial data. The proof relies upon an elegant closure theorem of Bertoin and Carraro and Duchon that reduces the study of shock statistics to Smoluchowski's coagulation equation with additive kernel, and upon our previous characterization of the domains of attraction of self-similar solutions for this equation.

preprint1999arXiv

Non-self-similar behavior in the LSW theory of Ostwald ripening

The classical Lifshitz-Slyozov-Wagner theory of domain coarsening predicts asymptotically self-similar behavior for the size distribution of a dilute system of particles that evolve by diffusional mass transfer with a common mean field. Here we consider the long-time behavior of measure-valued solutions for systems in which particle size is uniformly bounded, i.e., for initial measures of compact support. We prove that the long-time behavior of the size distribution depends sensitively on the initial distribution of the largest particles in the system. Convergence to the classically predicted smooth similarity solution is impossible if the initial distribution function is comparable to any finite power of distance to the end of the support. We give a necessary criterion for convergence to other self-similar solutions, and conditional stability theorems for some such solutions. For a dense set of initial data, convergence to any self-similar solution is impossible.