Source author record

Lenya Ryzhik

Lenya Ryzhik 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

22works
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

22 published item(s)

preprint2023arXiv

Pushed, pulled and pushmi-pullyu fronts of the Burgers-FKPP equation

We consider the long time behavior of the solutions to the Burgers-FKPP equation with advection of a strength $β\in\mathbb{R}$. This equation exhibits a transition from pulled to pushed front behavior at $β_c=2$. We prove convergence of the solutions to a traveling wave in a reference frame centered at a position $m_β(t)$ and study the asymptotics of the front location $m_β(t)$. When $β< 2$, it has the same form as for the standard Fisher-KPP equation established by Bramson \cite{Bramson1,Bramson2}: $m_β(t) = 2t - (3/2)\log(t) + x_\infty + o(1)$ as $t\to+\infty$. This form is typical of pulled fronts. When $β> 2$, the front is located at the position $m_β(t)=c_*(β)t+x_\infty+o(1)$ with $c_*(β)=β/2+2/β$, which is the typical form of pushed fronts. However, at the critical value $β_c = 2$, the expansion changes to $m_β(t) = 2t - (1/2)\log(t) + x_\infty + o(1)$, reflecting the "pushmi-pullyu" nature of the front. The arguments for $β<2$ rely on a new weighted Hopf-Cole transform that allows to control the advection term, when combined with additional steepness comparison arguments. The case $β>2$ relies on standard pushed front techniques. The proof in the case $β=β_c$ is much more intricate and involves arguments not usually encountered in the study of the Bramson correction. It relies on a somewhat hidden viscous conservation law structure of the Burgers-FKPP equation at $β_c=2$ and utilizes a dissipation inequality, which comes from a relative entropy type computation, together with a weighted Nash inequality involving dynamically changing weights.

preprint2022arXiv

Voting models and semilinear parabolic equations

We present probabilistic interpretations of solutions to semi-linear parabolic equations with polynomial nonlinearities in terms of the voting models on the genealogical trees of branching Brownian motion (BBM). These extend the connection between the Fisher-KPP equation and BBM discovered by McKean in~\cite{McK}. In particular, we present ``random outcome'' and ``random threshold'' voting models that yield any polynomial nonlinearity $f$ satisfying $f(0)=f(1)=0$ and a ``recursive up the tree'' model that allows to go beyond this restriction on $f$. We compute a few examples of particular interest; for example, we obtain a curious interpretation of the heat equation in terms of a nontrivial voting model and a ``group-based'' voting rule that leads to a probabilistic view of the pushmi-pullyu transition for a class of nonlinearities introduced by Ebert and van Saarloos.

preprint2020arXiv

Chemotaxis and Reactions in Biology

Chemotaxis plays a crucial role in a variety of processes in biology and ecology. Quite often it acts to improve efficiency of biological reactions. One example is the immune system signalling, where infected tissues release chemokines attracting monocytes to fight invading bacteria. Another example is reproduction, where eggs release pheromones that attract sperm. A macro scale example is flower scent appealing to pollinators. In this paper we consider a system of PDE designed to model such processes. Our interest is to quantify the effect of chemotaxis on reaction rates compared to pure reaction-diffusion. We limit consideration to surface chemotaxis, which is well motivated from the point of view of many applications. Our results provide the first insight into situations where chemotaxis can be crucial for reaction success, and where its effect is likely to be limited. The proofs are based on new analytical tools; a significant part of the paper is dedicated to building up the linear machinery that can be useful in more general settings. In particular we establish precise estimates on the rates of convergence to ground state for a class of Fokker-Planck operators with potentials that grow at a logarithmic rate at infinity. These estimates are made possible by a new sharp weak weighted Poincaré inequality improving in particular a result of Bobkov and Ledoux.

preprint2019arXiv

Current Distribution on Capacitive Electrode-Electrolyte Interfaces

The distribution of electric current on an electrode surface in electrolyte varies with time due to charge accumulation at a capacitive interface, as well as due to electrode kinetics and concentration polarization in the medium. Initially, the potential at the electrode-electrolyte interface is uniform, resulting in a non-uniform current distribution due to the uneven ohmic drop of the potential in the medium. Over time, however, the non-uniform current density causes spatially varying rate of the charge accumulation at the interface, breaking down its equipotentiality. We developed an analytical model to describe such transition at a capacitive interface when the current is below the mass-transfer limitation, and demonstrated that the steady distribution of the current is achieved when the current density is proportional to the capacitance per unit area, which leads to linear voltage ramp at the electrode. More specific results regarding the dynamics of this transition are provided for a disk electrode, along with an experimental validation of the theoretical result. These findings are important for many electrochemical applications, and in particular, for proper design of the electro-neural interfaces.

preprint2019arXiv

Global well-posedness for the Euler alignment system with mildly singular interactions

We consider the Euler alignment system with mildly singular interaction kernels. When the local repulsion term is of the fractional type, global in time existence of smooth solutions was proved in\cite{do2018global,shvydkoy2017eulerian1,shvydkoy2017eulerian2,shvydkoy2017eulerian3}. Here, we consider a class of less singular interaction kernels and establish the global regularity of solutions as long as the interaction kernels are not integrable. The proof relies on modulus of continuity estimates for a class of parabolic integro-differential equations with a drift and mildly singular kernels.

preprint2016arXiv

Convergence to a single wave in the Fisher-KPP equation

We study the large time asymptotics of a solution of the Fisher-KPP reaction-diffusion equation, with an initial condition that is a compact perturbation of a step function. A well-known result of Bramson states that, in the reference frame moving as $2t - ({3}/{2}) \log t +x_\infty$, the solution of the equation converges as $t\to+\infty$ to a translate of the traveling wave corresponding to the minimal speed~$c_*=2$. The constant $x_\infty$ depends on the initial condition $u(0,x)$. The proof is elaborate, and based on probabilistic arguments. The purpose of this paper is to provide a simple proof based on PDE arguments.

preprint2016arXiv

The Bramson logarithmic delay in the cane toads equations

We study a nonlocal reaction-diffusion-mutation equation modeling the spreading of a cane toads population structured by a phenotypical trait responsible for the spatial diffusion rate. When the trait space is bounded, the cane toads equation admits traveling wave solutions [7]. Here, we prove a Bramson type spreading result: the lag between the position of solutions with localized initial data and that of the traveling waves grows as (3/(2$λ$ *)) log t. This result relies on a present-time Harnack inequality which allows to compare solutions of the cane toads equation to those of a Fisher-KPP type equation that is local in the trait variable.

preprint2015arXiv

Finite time singularity formation for the modified SQG patch equation

It is well known that the incompressible Euler equations in two dimensions have globally regular solutions. The inviscid surface quasi-geostrophic (SQG) equation has a Biot-Savart law which is one derivative less regular than in the Euler case, and the question of global regularity for its solutions is still open. We study here the patch dynamics in the half-plane for a family of active scalars which interpolates between these two equations, via a parameter $α\in[0,\frac 12]$ appearing in the kernels of their Biot-Savart laws. The values $α=0$ and $α=\frac 12$ correspond to the 2D Euler and SQG cases, respectively. We prove global in time regularity for the 2D Euler patch model, even if the patches initially touch the boundary of the half-plane. On the other hand, for any sufficiently small $α>0$, we exhibit initial data which lead to a singularity in finite time. Thus, these results show a phase transition in the behavior of solutions to these equations, and provide a rigorous foundation for classifying the 2D Euler equations as critical.

preprint2015arXiv

Radiative transfer with long-range interactions: regularity and asymptotics

This work is devoted to radiative transfer equations with long-range interactions. Such equations arise in the modeling of high frequency wave propagation in random media with long-range dependence. In the regime we consider, the singular collision operator modeling the interaction between the wave and the medium is conservative, and as a consequence wavenumbers take values on the unit sphere. Our goals are to investigate the regularizing effects of grazing collisions, the diffusion limit, and the peaked forward limit. As in the case where wavenumbers take values in R^{d+1}, we show that the transport operator is hypoelliptic, so that the solutions are infinitely differerentiable in all variables. Using probabilistic techniques, we show as well that the diffusion limit can be carried on as in the case of a regular collision operator, and as a consequence that the diffusion coefficient is non-zero and finite. We finally consider the regime where grazing collisions are dominant.

preprint2015arXiv

The random Schrödinger equation: homogenization in time-dependent potentials

We analyze the solutions of the Schrödinger equation with the low frequency initial data and a time-dependent weakly random potential. We prove a homogenization result for the low frequency component of the wave field. We also show that the dynamics generates a non-trivial energy in the high frequencies, which do not homogenize -- the high frequency component of the wave field remains random and the evolution of its energy is described by a kinetic equation. The transition from the homogenization of the low frequencies to the random limit of the high frequencies is illustrated by understanding the size of the small random fluctuations of the low frequency component.

preprint2015arXiv

The random Schrödinger equation: slowly decorrelating time-dependent potentials

We analyze the weak-coupling limit of the random Schrödinger equation with low frequency initial data and a slowly decorrelating random potential. For the probing signal with a sufficiently long wavelength, we prove a homogenization result, that is, the properly compensated wave field admits a deterministic limit in the "very low" frequency regime. The limit is "anomalous" in the sense that the solution behaves as $\exp(-Dt^{s})$ with $s>1$ rather than the "usual"~$\exp(-Dt)$ homogenized behavior when the random potential is rapidly decorrelating. Unlike in rapidly decorrelating potentials, as we decrease the wavelength of the probing signal, stochasticity appears in the asymptotic limit -- there exists a critical scale depending on the random potential which separates the deterministic and stochastic regimes.

preprint2014arXiv

Anisotropic interactions in a first-order aggregation model: a proof of concept

We extend a well-studied ODE model for collective behaviour by considering anisotropic interactions among individuals. Anisotropy is modelled by limited sensorial perception of individuals, that depends on their current direction of motion. Consequently, the first-order model becomes implicit, and new analytical issues, such as non-uniqueness and jump discontinuities in velocities, are being raised. We study the well-posedness of the anisotropic model and discuss its modes of breakdown. To extend solutions beyond breakdown we propose a relaxation system containing a small parameter $\varepsilon$, which can be interpreted as a small amount of inertia or response time. We show that the limit $\varepsilon \to 0$ can be used as a jump criterion to select the physically correct velocities. In smooth regimes, the convergence of the relaxation system as $\varepsilon \to 0$ is guaranteed by a theorem due to Tikhonov. We illustrate the results with numerical simulations in two dimensions.

preprint2013arXiv

On the nonlocal Fisher-KPP equation: steady states, spreading speed and global bounds

We consider the Fisher-KPP equation with a non-local interaction term. We establish a condition on the interaction that allows for existence of non-constant periodic solutions, and prove uniform upper bounds for the solutions of the Cauchy problem, as well as upper and lower bounds on the spreading rate of the solutions with compactly supported initial data.

preprint2013arXiv

Traveling Wave Solutions in a Reaction-Diffusion Model for Criminal Activity

We study a reaction-diffusion system of partial differential equations, which can be taken to be a basic model for criminal activity. We show that the assumption of a populations natural tendency towards crime significantly changes the long-time behavior of criminal activity patterns. Under the right assumptions on these natural tendencies we first show that there exists traveling wave solutions connecting zones with no criminal activity and zones with high criminal activity, known as hotspots. This corresponds to an invasion of criminal activity onto all space. Second, we study the problem of preventing such invasions by employing a finite number of resources that reduce the payoff committing a crime in a finite region. We make the concept of wave propagation mathematically rigorous in this situation by proving the existence of entire solutions that approach traveling waves as time approaches negative infinity. Furthermore, we characterize the minimum amount of resources necessary to prevent the invasion in the case when prevention is possible. Finally, we apply our theory to what is commonly known as the gap problem in the excitable media literature, proving existing conjectures in the literature.

preprint2012arXiv

A differential equations approach to $l_1$-minimization with applications to array imaging

We present an ordinary differential equations approach to the analysis of algorithms for constructing $l_1$ minimizing solutions to underdetermined linear systems of full rank. It involves a relaxed minimization problem whose minimum is independent of the relaxation parameter. An advantage of using the ordinary differential equations is that energy methods can be used to prove convergence. The connection to the discrete algorithms is provided by the Crandall-Liggett theory of monotone nonlinear semigroups. We illustrate the effectiveness of the discrete optimization algorithm in some sparse array imaging problems.

preprint2012arXiv

Asymptotics of the solutions of the stochastic lattice wave equation

We consider the long time limit theorems for the solutions of a discrete wave equation with a weak stochastic forcing. The multiplicative noise conserves the energy and the momentum. We obtain a time-inhomogeneous Ornstein-Uhlenbeck equation for the limit wave function that holds both for square integrable and statistically homogeneous initial data. The limit is understood in the point-wise sense in the former case, and in the weak sense in the latter. On the other hand, the weak limit for square integrable initial data is deterministic.

preprint2012arXiv

Biomixing by chemotaxis and efficiency of biological reactions: the critical reaction case

Many phenomena in biology involve both reactions and chemotaxis. These processes can clearly influence each other, and chemotaxis can play an important role in sustaining and speeding up the reaction. In continuation of our earlier work, we consider a model with a single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The model is motivated, in particular, by the studies of coral broadcast spawning, where experimental observations of the efficiency of fertilization rates significantly exceed the data obtained from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical secreted by egg gametes) into account. We consider the case of the weakly coupled quadratic reaction term, which is the most natural from the biological point of view and was left open. The result is that similarly to higher power coupling, the chemotaxis plays a crucial role in ensuring efficiency of reaction. However, mathematically, the picture is quite different in the quadratic reaction case and is more subtle. The reaction is now complete even in the absence of chemotaxis, but the timescales are very different. Without chemotaxis, the reaction is very slow, especially for the weak reaction coupling coefficient. With chemotaxis, the timescale and efficiency of reaction are independent of the coupling parameter.

preprint2011arXiv

Biomixing by chemotaxis and enhancement of biological reactions

Many processes in biology involve both reactions and chemotaxis. However, to the best of our knowledge, the question of interaction between chemotaxis and reactions has not yet been addressed either analytically or numerically. We consider a model with a single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The model is motivated, in particular, by studies of coral broadcast spawning, where experimental observations of the efficiency of fertilization rates significantly exceed the data obtained from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical secreted by egg gametes) into account. We prove that in the framework of our model, chemotaxis plays a crucial role. There is a rigid limit to how much the fertilization efficiency can be enhanced if there is no chemotaxis but only advection and diffusion. On the other hand, when chemotaxis is present, the fertilization rate can be arbitrarily close to being complete provided that the chemotactic attraction is sufficiently strong. Moreover, an interesting feature of the estimates on fertilization rate and timescales in the chemotactic case is that they do not depend on the amplitude of the reaction term.

preprint2011arXiv

From homogenization to averaging in cellular flows

We consider an elliptic eigenvalue problem with a fast cellular flow of amplitude $A$, in a two-dimensional domain with $L^2$ cells. For fixed $A$, and $L \to \infty$, the problem homogenizes, and has been well studied. Also well studied is the limit when $L$ is fixed, and $A \to \infty$. In this case the solution equilibrates along stream lines. In this paper, we show that if \textit{both} $A \to \infty$ and $L \to \infty$, then a transition between the homogenization and averaging regimes occurs at $A \approx L^4$. When $A\gg L^4$, the principal Dirichlet eigenvalue is approximately constant. On the other hand, when $A\ll L^4$, the principal eigenvalue behaves like ${\bar σ(A)}/L^2$, where $\bar σ(A) \approx \sqrt{A} I$ is the effective diffusion matrix. A similar transition is observed for the solution of the exit time problem. The proof in the homogenization regime involves bounds on the second correctors. Miraculously, if the slow profile is quadratic, these estimates can be obtained using drift independent $L^p \to L^\infty$ estimates for elliptic equations with an incompressible drift. This provides effective sub and super-solutions for our problem.

preprint2010arXiv

A simple model for asset price bubble formation and collapse

We consider a simple stochastic differential equation for modeling bubbles in social context. A prime example is bubbles in asset pricing, but similar mechanisms may control a range of social phenomena driven by psychological factors (for example, popularity of rock groups, or a number of students pursuing a given major). Our goal is to study the simplest possible model in which every term has a clear meaning and which demonstrates several key behaviors. The main factors that enter are tendency of mean reversion to a stable value, speculative social response triggered by trend following and random fluctuations. The interplay of these three forces may lead to bubble formation and collapse. Numerical simulations show that the equation has distinct regimes depending on the values of the parameters. We perform rigorous analysis of the weakly random regime, and study the role of change in fundamentals in igniting the bubble.

preprint2010arXiv

Existence and Non-existence of Fisher-KPP Transition Fronts

We consider Fisher-KPP-type reaction-diffusion equations with spatially inhomogeneous reaction rates. We show that a sufficiently strong localized inhomogeneity may prevent existence of transition-front-type global in time solutions while creating a global in time bump-like solution. This is the first example of a medium in which no reaction-diffusion transition front exists. A weaker localized inhomogeneity leads to existence of transition fronts but only in a finite range of speeds. These results are in contrast with both Fisher-KPP reactions in homogeneous media as well as ignition-type reactions in inhomogeneous media.