Search

This manuscript deals with the three-dimensional version of a flux-limited Keller-Segel system coupled to the incompressible Stokes equations through transport and buoyancy. The main goal consists in verifying that within a certain parameter regime, known as being optimal therefor in some fluid-free simplification, a feature of blow-up prevention by suitably strong flux limitation persists also in the framework of the considered full chemotaxis-fluid system. To achieve this, as a secondary objective of possibly independent interest the manuscript separately establishes some conditional bounds for fluid fields and taxis gradients in a fairly general setting. The application thereof to the specific problem under consideration thereafter facilitates the derivation of a result on global existence of bounded classical solutions for widely arbitrary initial data actually, indeed within the entire and essentially optimal parameter range, through an argument which appears to be signficantly condensed when compared to reasonings pursued in previous related works.

We study the analyticity of the semigroups generated by a class of degenerate second order differential operators in the space $C(S_d)$, where $S_d$ is the canonical simplex of $\R^d$. The semigroups arise from the theory of Fleming--Viot processes in population genetics., Comment: 32 pages

We study the large-time behaviour of the solutions of the evolution equation involving nonlinear diffusion and gradient absorption, $$ \partial_t u - \Delta_p u + |\nabla u|^q=0 . $$ We consider the problem posed for $x\in \real^N$ and t>0 with nonnegative and compactly supported initial data. We take the exponent p>2 which corresponds to slow p-Laplacian diffusion. The main feature of the paper is that the exponent q takes the critical value q=p-1 which leads to interesting asymptotics. This is due to the fact that in this case both the Hamilton-Jacobi term $ |\nabla u|^q$ and the diffusive term $\Delta_p u$ have a similar size for large times. The study performed in this paper shows that a delicate asymptotic equilibrium happens, so that the large-time behaviour of the solutions is described by a rescaled version of a suitable self-similar solution of the Hamilton-Jacobi equation $|\nabla W|^{p-1}=W$, with logarithmic time corrections. The asymptotic rescaled profile is a kind of sandpile with a cusp on top, and it is independent of the space dimension., Comment: 26 pages, 2 figures

In this paper we consider a two-phase flow problem in porous media and study its singular limit as the viscosity of the air tends to zero; more precisely, we prove the convergence of subsequences to solutions of a generalized Richards model.

The Perona-Malik equation is a celebrated example of forward-backward parabolic equation. The forward behavior takes place in the so-called subcritical region, in which the gradient of the solution is smaller than a fixed threshold. In this paper we show that this subcritical region evolves in a different way in the following three cases: dimension one, radial solutions in dimension greater than one, general solutions in dimension greater than one. In the first case subcritical regions increase, but there is no estimate on the expansion rate. In the second case they expand with a positive rate and always spread over the whole domain after a finite time, depending only on the (outer) radius of the domain. As a by-product, we obtain a non-existence result for global-in-time classical radial solutions with large enough gradient. In the third case we show an example where subcritical regions do not expand. Our proofs exploit comparison principles for suitable degenerate and non-smooth free boundary problems., Comment: 23 pages, 2 figures

We establish the longtime existence and convergence results of the mean curvature flow of entire Lagrangian graphs in Pseudo-Euclidean space which is related to Logarithmic gradient flow., Comment: 13pages

By making use of the approximation method, we obtain the existence and regularity of the viscosity solutions for the generalized mean curvature flow. The asymptotic behavior of the flow is also considered. In particular, the Dirichlet problem of the degenerate elliptic equation $$ -|\nabla v|(\mathrm{div}(\frac{\nabla v}{|\nabla v|})+\nu)=0 $$ is solvable in viscosity sense., Comment: 13pages

It is shown that the curvature function satisfies a nonlinear evolution equation under the general curve shortening flow and a detailed asymptotic behavior of the closed curves is presented when they contract to a point in finite time., Comment: 23pages

We analyze the diffusion processes associated to equations of Wright-Fisher type in one spatial dimension. These are defined by a degenerate second order operator on the interval [0, 1], where the coefficient of the second order term vanishes simply at the endpoints, and the first order term is an inward-pointing vector field. We consider various aspects of this problem, motivated by applications in population genetics, including a sharp regularity theory for the zero flux boundary conditions, as well as a derivation of the precise asymptotics for solutions of this equation, both as t goes to 0 and infinity, and as x goes to 0, 1., Comment: 55 pages

In this paper we analyze a mathematical model focusing on key events of the cells invasion process. Global well-possedness and asymptotic behaviour of nonnegative solutions to the corresponding coupled system of three nonlinear partial differential equations are studied., Comment: 29 pages

We consider the Cauchy problem for the Perona-Malik equation in an open subset of R^{n}, with Neumann boundary conditions. It is well known that in the one-dimensional case this problem does not admit any global C^{1} solution if the initial condition is transcritical, namely when the norm of the gradient of the initial condition is smaller than 1 in some region, and larger than 1 in some other region In this paper we show that this result cannot be extended to higher dimension. We show indeed that for n >= 2 the problem admits radial solutions of class C^{2,1} with a transcritical initial condition., Comment: 38 pages, 3 figures

We give sufficient conditions under which the convergence of finite difference approximations in the space variable of possibly degenerate second order parabolic and elliptic equations can be accelerated to any given order of convergence by Richardson's method., Comment: 30 pages

In this review paper we present a stable Lagrangian numerical method for computing plane curves evolution driven by the fourth order geometric equation. The numerical scheme and computational examples are presented., Comment: submitted to: Proceedings of Equadiff 2007 Conference

We study qualitative and quantitative properties of local weak solutions of the fast $p$-Laplacian equation, $\partial_t u=\Delta_{p}u$, with $1

This paper deals with homogenization of second order divergence form parabolic operators with locally stationary coefficients. Roughly speaking, locally stationary coefficients have two evolution scales: both an almost constant microscopic one and a smoothly varying macroscopic one. The homogenization procedure aims to give a macroscopic approximation that takes into account the microscopic heterogeneities. This paper follows "Diffusion in a locally stationary random environment" (published in Probability Theory and Related Fields) and improves this latter work by considering possibly degenerate diffusion matrices. The geometry of the homogenized equation shows that the particle is trapped in subspace of R^d., Comment: To appear in Annales de l'institut henri poincare (link of the journal: http://www.imstat.org/aihp/)

The parabolic N-membranes problem for the p-Laplacian and the complete order constraint on the components of the solution is studied in what concerns the approximation, the regularity and the stability of the variational solutions. We extend to the evolutionary case the characterization of the Lagrange multipliers associated with the ordering constraint in terms of the characteristic functions of the coincidence sets. We give continuous dependence results, and study the asymptotic behavior as $t \to \infty$ of the solution and the coincidence sets, showing that they converge to their stationary counterparts., Comment: 16 pages

A new simple Lagrangian method with favorable stability and efficiency properties for computing general plane curve evolutions is presented. The method is based on the flowing finite volume discretization of the intrinsic partial differential equation for updating the position vector of evolving family of plane curves. A curve can be evolved in the normal direction by a combination of fourth order terms related to the intrinsic Laplacian of the curvature, second order terms related to the curvature, first order terms related to anisotropy and by a given external velocity field. The evolution is numerically stabilized by an asymptotically uniform tangential redistribution of grid points yielding the first order intrinsic advective terms in the governing system of equations. By using a semi-implicit in time discretization it can be numerically approximated by a solution to linear penta-diagonal systems of equations (in presence of the fourth order terms) or tri-diagonal systems (in the case of the second order terms). Various numerical experiments of plane curve evolutions, including, in particular, nonlinear, anisotropic and regularized backward curvature flows, surface diffusion and Willmore flows, are presented and discussed., Comment: 15 pages

We study the self-similar solutions of any sign of the equation u_{t}-div(|∇u|^{p-2}∇u)=|u|^{q-1}u, in R^{N}, where p,q>1. We extend the results of Haraux-Weissler obtained for p=2 to the case q>p-1>0. In particular we study the existence of slow or fast decaying solutions. For given t>0, the fast solutions u(t,.) have a compact support in R^{N} when p>2, and |x|^{p/(2-p)}u(t,x) is bounded at infinity when p<2. We describe the behaviour for large |x| of all the solutions. According to the position of q with respect to the first critical exponent p-1+p/N and the critical Sobolev exponent q^{∗}, we study the existence of positive solutions, or the number of the zeros of u(t,.). We prove that any solution u(t,.) is oscillatory when p<2 and q is closed to 1.

We study the Cauchy problem for the parabolic infinity Laplace equation. We prove a new comparison principle and obtain uniqueness of viscosity solutions in the class of functions with a polinomial growth at infinity, improving previous results obtained assuming a linear growth., Comment: This paper has been withdrawn by the author due to a gap in the proof of Lemma 2.3. In fact, the condition $\Phi (x,y,0,0) \leq 0$ does not directly follow from (6) and (7)

We prove a global existence result for a model describing the swarming phenomenon of the bacterium Proteus mirabilis. The model consists of an ordinary differential equation coupled with an age-structured equation involving nonlinear degenerate diffusion and an additional drift term.

This paper addresses the existence and regularity of weak solutions for a fully parabolic model of chemotaxis, with prevention of overcrowding, that degenerates in a two-sided fashion, including an extra nonlinearity represented by a $p$-Laplacian diffusion term. To prove the existence of weak solutions, a Schauder fixed-point argument is applied to a regularized problem and the compactness method is used to pass to the limit. The local H\"older regularity of weak solutions is established using the method of intrinsic scaling. The results are a contribution to showing, qualitatively, to what extent the properties of the classical Keller-Segel chemotaxis models are preserved in a more general setting. Some numerical examples illustrate the model.

The convergence to non-diffusive self-similar solutions is investigated for non-negative solutions to the Cauchy problem $\partial_t u = \Delta_p u + |\nabla u|^q$ when the initial data converge to zero at infinity. Sufficient conditions on the exponents $p>2$ and $q>1$ are given that guarantee that the diffusion becomes negligible for large times and the $L^\infty$-norm of $u(t)$ converges to a positive value as $t\to\infty$.

We present a fully adaptive multiresolution scheme for spatially two-dimensional, possibly degenerate reaction-diffusion systems, focusing on combustion models and models of pattern formation and chemotaxis in mathematical biology. Solutions of these equations in these applications exhibit steep gradients, and in the degenerate case, sharp fronts and discontinuities. The multiresolution scheme is based on finite volume discretizations with explicit time stepping. The multiresolution representation of the solution is stored in a graded tree. By a thresholding procedure, namely the elimination of leaves that are smaller than a threshold value, substantial data compression and CPU time reduction is attained. The threshold value is chosen optimally, in the sense that the total error of the adaptive scheme is of the same slope as that of the reference finite volume scheme. Since chemical reactions involve a large range of temporal scales, but are spatially well localized (especially in the combustion model), a locally varying adaptive time stepping strategy is applied. It turns out that local time stepping accelerates the adaptive multiresolution method by a factor of two, while the error remains controlled., Comment: 27 pages, 14 figures

This work deals with the numerical solution of the monodomain and bidomain models of electrical activity of myocardial tissue. The bidomain model is a system consisting of a possibly degenerate parabolic PDE coupled with an elliptic PDE for the transmembrane and extracellular potentials, respectively. This system of two scalar PDEs is supplemented by a time-dependent ODE modeling the evolution of the so-called gating variable. In the simpler sub-case of the monodomain model, the elliptic PDE reduces to an algebraic equation. Two simple models for the membrane and ionic currents are considered, the Mitchell-Schaeffer model and the simpler FitzHugh-Nagumo model. Since typical solutions of the bidomain and monodomain models exhibit wavefronts with steep gradients, we propose a finite volume scheme enriched by a fully adaptive multiresolution method, whose basic purpose is to concentrate computational effort on zones of strong variation of the solution. Time adaptivity is achieved by two alternative devices, namely locally varying time stepping and a Runge-Kutta-Fehlberg-type adaptive time integration. A series of numerical examples demonstrates thatthese methods are efficient and sufficiently accurate to simulate the electrical activity in myocardial tissue with affordable effort. In addition, an optimalthreshold for discarding non-significant information in the multiresolution representation of the solution is derived, and the numerical efficiency and accuracy of the method is measured in terms of CPU time speed-up, memory compression, and errors in different norms., Comment: 25 pages, 41 figures

A fully adaptive finite volume multiresolution scheme for one-dimensional strongly degenerate parabolic equations with discontinuous flux is presented. The numerical scheme is based on a finite volume discretization using the Engquist--Osher approximation for the flux and explicit time--stepping. An adaptivemultiresolution scheme with cell averages is then used to speed up CPU time and meet memory requirements. A particular feature of our scheme is the storage of the multiresolution representation of the solution in a dynamic graded tree, for the sake of data compression and to facilitate navigation. Applications to traffic flow with driver reaction and a clarifier--thickener model illustrate the efficiency of this method.

This text is a survey of recent results obtained by the author and collaborators on different problems for non-self-adjoint operators. The topics are: Kramers-Fokker-Planck type operators, spectral asymptotics in two dimensions and Weyl asymptotics for the eigenvalues of non-self-adjoint operators with small random perturbations. In the introduction we also review the notion of pseudo-spectrum and its relation to non-self-adjoint spectral problems.

We proved the existence and uniqueness of a traveling wave solution to the thin film equation with a Navier slip condition at the liquid-solid interface. We obtain explicit lower and upper bounds for the solution and an absolute error estimate of approximation of a solution to the thin films equation by the traveling-wave solution., Comment: 14 pages, 2 figure

For a class of degenerate diffusion processes of rank 2, i.e. when only Poisson brackets of order one are needed to span the whole space, we obtain a parametrix representation of the density from which we derive some explicit Gaussian controls that characterize the additional singularity induced by the degeneracy. We then give a local limit theorem with the usual convergence rate for an associated Markov chain approximation. The key point is that the "weak" degeneracy allows to exploit the techniques first introduced by Konakov and Molchanov and then developed by Konakov and Mammen that rely on Gaussian approximations., Comment: 33 pages

We study the asymptotic behavior of nonnegative solutions of the semilinear parabolic problem {u_t=\Delta u + u^{p}, x\in\mathbb{R}^{N}, t>0 u(0)=u_{0}, x\in\mathbb{R}^{N}, t=0. It is known that the nonnegative solution $u(t)$ of this problem blows up in finite time for $1 1+ 2/N$ and the norm of $u_{0}$ is small enough, the problem admits global solution. In this work, we use the entropy method to obtain the decay rate of the global solution $u(t)$., Comment: 15

This paper is devoted to the analysis of non-negative solutions for a generalisation of the classical parabolic-elliptic Patlak-Keller-Segel system with $d\ge3$ and porous medium-like non-linear diffusion. Here, the non-linear diffusion is chosen in such a way that its scaling and the one of the Poisson term coincide. We exhibit that the qualitative behaviour of solutions is decided by the initial mass of the system. Actually, there is a sharp critical mass $M_c$ such that if $M \in (0,M_c]$ solutions exist globally in time, whereas there are blowing-up solutions otherwise. We also show the existence of self-similar solutions for $M \in (0,M_c)$. While characterising the eventual infinite time blowing-up profile for $M=M_c$, we observe that the long time asymptotics are much more complicated than in the classical Patlak-Keller-Segel system in dimension two.

Proteus mirabilis are bacteria that make strikingly regular spatial-temporal patterns on agar surfaces. In this paper we investigate a mathematical model that has been shown to display these structures when solved numerically. The model consists of an ordinary differential equation coupled with a partial differential equation involving a first-order hyperbolic aging term together with nonlinear degenerate diffusion. The system is shown to admit global weak solutions.

Motivated by emerging applications from imaging processing, the heat flow of a generalized $p$-harmonic map into spheres is studied for the whole spectrum, $1\leq p<\infty$, in a unified framework. The existence of global weak solutions is established for the flow using the energy method together with a regularization and a penalization technique. In particular, a $BV$-solution concept is introduced and the existence of such a solution is proved for the 1-harmonic map heat flow. The main idea used to develop such a theory is to exploit the properties of measures of the forms $\cA\cdot\nab\bv$ and $\cA\wedge\nab\bv$; which pair a divergence-$L^1$, or a divergence-measure, tensor field $\cA$, and a $BV$-vector field $\bv$. Based on these analytical results, a practical fully discrete finite element method is then proposed for approximating weak solutions of the $p$-harmonic map heat flow, and the convergence of the proposed numerical method is also established., Comment: 27 pages

The goal of this paper is to give a non-local sufficient condition for generalized Poincar\'e inequalities, which extends the well-known Bakry-Emery condition. Such generalized Poincar\'e inequalities have been introduced by W. Beckner in the gaussian case and provide, along the Ornstein-Uhlenbeck flow, the exponential decay of some generalized entropies which interpolate between the $L^2$ norm and the usual entropy. Our criterion improves on results which, for instance, can be deduced from the Bakry-Emery criterion and Holley-Stroock type perturbation results. In a second step, we apply the same strategy to non-linear equations of porous media type. This provides new interpolation inequalities and decay estimates for the solutions of the evolution problem. The criterion is again a non-local condition based on the positivity of the lowest eigenvalue of a Schr\"odinger operator. In both cases, we relate the Fisher information with its time derivative. Since the resulting criterion is non-local, it is better adapted to potentials with, for instance, a non-quadratic growth at infinity, or to unbounded perturbations of the potential.

In [5] the authors obtained Mean-Field backward stochastic differential equations (BSDE) associated with a Mean-field stochastic differential equation (SDE) in a natural way as limit of some highly dimensional system of forward and backward SDEs, corresponding to a large number of ``particles'' (or ``agents''). The objective of the present paper is to deepen the investigation of such Mean-Field BSDEs by studying them in a more general framework, with general driver, and to discuss comparison results for them. In a second step we are interested in partial differential equations (PDE) whose solutions can be stochastically interpreted in terms of Mean-Field BSDEs. For this we study a Mean-Field BSDE in a Markovian framework, associated with a Mean-Field forward equation. By combining classical BSDE methods, in particular that of ``backward semigroups" introduced by Peng [14], with specific arguments for Mean-Field BSDEs we prove that this Mean-Field BSDE describes the viscosity solution of a nonlocal PDE. The uniqueness of this viscosity solution is obtained for the space of continuous functions with polynomial growth. With the help of an example it is shown that for the nonlocal PDEs associated to Mean-Field BSDEs one cannot expect to have uniqueness in a larger space of continuous functions., Comment: The results were presented by Rainer Buckdahn at the Mittag-Leffler Institute (Stockholm, Sweden) in November 8th, 2007. The paper was submitted

We study the large-time behaviour of the solutions $u$ of the evolution equation involving nonlinear diffusion and gradient absorption $\partial_t u - \Delta_p u + |\nabla u|^q=0$. We consider the problem posed for $x\in {\mathbb R}^N $ and $t>0$ with non-negative and compactly supported initial data. We take the exponent $p>2$ which corresponds to slow $p$-Laplacian diffusion, and the exponent $q$ in the superlinear range $1

We demonstrate the equivalence between the two popular models of nonlinear diffusion, the porous medium equation and the p-Laplacian equation. The equivalence is shown at the level of selfsimilar solutions., Comment: 6 pages

This memoir attempts at a systematic study of convergence to stationary state for certain classes of degenerate diffusive equations, by means of well-chosen Lyapunov functionals. Typical examples are the kinetic Fokker--Planck and Boltzmann equation. Many open problems and possible directions for future research are discussed.

It is shown that an elliptic scattering operator $A$ on a compact manifold with boundary with coefficients in the bounded operators of a bundle of Banach spaces of class (HT) and Pisier's property $(\alpha)$ has maximal regularity (up to a spectral shift), provided that the spectrum of the principal symbol of $A$ on the scattering cotangent bundle of the manifold avoids the right half-plane. This is deduced directly from a Seeley theorem, i.e. the resolvent is represented in terms of pseudodifferential operators with R-bounded symbols, thus showing by an iteration argument the R-boundedness of $\lambda(A-\lambda)^{-1}$ for $\Re(\lambda) \geq 0$. To this end, elements of a symbolic and operator calculus of pseudodifferential operators with R-bounded symbols are introduced. The significance of this method for proving maximal regularity results for partial differential operators is underscored by considering also a more elementary situation of anisotropic elliptic operators on $R^d$ with operator valued coefficients., Comment: 21 pages

We prove a sharp version of the Hopf boundary point lemma for Black-Scholes type equations. We also investigate the existence and the regularity of the spatial derivative of the solutions at the spatial boundary., Comment: 10 pages

We consider local solutions of the two-phase Stefan problem with a "mushy" region. We show that if a (distributional) solution u is locally square integrable then the temperature is continuous., Comment: 11 pages

We study a Hele-Shaw problem with a mushy region obtained as a Mesa type limit of one phase Stefan problems in exterior domains. We study the convergence, determine some of the qualitative properties and regularity of the unique limiting solution, and prove regularity of the free boundary of this limit under very general conditions on the initial data. Indeed, our results handle changes in topology and multiple injection slots., Comment: 28 pages, 4 figures, latex file - Corrections to the sixth section of the last version

We introduce a variation of the classical Ricci flow equation that modifies the unit volume constraint of that equation to a scalar curvature constraint. The resulting equations are named the Conformal Ricci Flow Equations because of the role that conformal geometry plays in constraining the scalar curvature. These equations are analogous to the incompressible Navier-Stokes equations of fluid mechanics inasmuch as a conformal pressure arises as a Lagrange multiplier to conformally deform the metric flow so as to maintain the scalar curvature constraint. The equilibrium points are Einstein metrics with a negative Einstein constant and the conformal pressue is shown to be zero at an equilibrium point and strictly positive otherwise. The geometry of the conformal Ricci flow is discussed as well as the remarkable analytic fact that the constraint force does not lose derivatives and thus analytically the conformal Ricci equation is a bounded perturbation of the classical unnormalized Ricci equation. That the constraint force does not lose derivatives is exactly analogous to the fact that the real physical pressure force that occurs in the Navier-Stokes equations is a bounded function of the velocity. Using a nonlinear Trotter product formula, existence and uniqueness of solutions to the conformal Ricci flow equations is proven. Lastly, we discuss potential applications to Perelman's proposed implementation of Hamilton's program to prove Thurston's 3-manifold geometrization conjectures., Comment: 52 pages, 1 figure

The theory of degenerate parabolic equations of the forms \[ u_t=(\Phi(u_x))_{x} \quad {\rm and} \quad v_{t}=(\Phi(v))_{xx} \] is used to analyze the process of contour enhancement in image processing, based on the evolution model of Sethian and Malladi. The problem is studied in the framework of nonlinear diffusion equations. It turns out that the standard initial-value problem solved in this theory does not fit the present application since it it does not produce image concentration. Due to the degenerate character of the diffusivity at high gradient values, a new free boundary problem with singular boundary data can be introduced, and it can be solved by means of a non-trivial problem transformation. The asymptotic convergence to a sharp contour is established and rates calculated., Comment: 29 pages, includes 6 figures

We prove estimates for the stationary state $n$-point functions at zero molecular diffusivity in the Kraichnan model. This is done by proving upper bounds for the heat kernels and Green's functions of the degenerate elliptic operators $M_n$ that occur in the Hopf equations for the $n$-point functions., Comment: 52 pages