Blow-up and pattern formation in hyperbolic models for chemotaxis in 1-D

In this paper we study finite time blow-up of solutions of a hyperbolic model for chemotaxis. Using appropriate scaling this hyperbolic model leads to a parabolic model as studied by Othmer and Stevens (1997) and Levine and Sleeman (1997). In the latter paper, explicit solutions which blow-up in finite time were constructed. Here, we adapt their method to construct a corresponding blow-up solution of the hyperbolic model. This construction enables us to compare the blow-up times of the corresponding models. We find that the hyperbolic blow-up is always later than the parabolic blow-up. Moreover, we show that solutions of the hyperbolic problem become negative near blow-up. We calculate the zero-turning-rate time explicitly and we show that this time can be either larger or smaller than the parabolic blow-up time. The blow-up models as discussed here and elsewhere are limiting cases of more realistic models for chemotaxis. At the end of the paper we discuss the relevance to biology and exhibit numerical solutions of more realistic models.     

Global Solutions of Some Chemotaxis and Angiogenesis Systems in High Space Dimensions

We consider two simple conservative systems of parabolic-elliptic and parabolic-degenerate type arising in modeling chemotaxis and angiogenesis. Both systems share the same property that when the norm of initial data is small enough, where d 2 is the space dimension, then there is a global (in time) weak solution that stays in all the L^p spaces with max This result is already known for the parabolic-elliptic system of chemotaxis, moreover blow-up can occur in finite time for large initial data and Dirac concentrations can occur. For the parabolic-degenerate system of angiogenesis in two dimensions, we also prove that weak solutions (which are equi-integrable in L¹) exist even for large initial data. But break-down of regularity or propagation of smoothness is an open problem.Lecture by B. Perthame held at the Presentation of MJM, Milano, October 18, 2002Received: March, 2003     

Numerical approximations and regularizations of Riccati equations arising in hyperbolic dynamics with unbounded control operators

This paper provides an approximation theory for numerical computations of the solutions to algebraic Riccati equations arising in hyperbolic, boundary control problems. One of the difficulties in the approximation theory for Riccati equations is that many attractive numerical methods (such as standard finite elements) do not satisfy a uniform stabilizability condition, which is necessary for the stability of the approximate Riccati solutions. To deal with these problems, a regularizationapproximation technique, based on the introduction of special artificial terms to the dynamics of the original model, is proposed. The need for this regularization appears to be a distinct feature of hyperbolic (hyperbolic-like) equations, rather than parabolic (parabolic-like) problems where the smoothing effect of the dynamics is beneficial for the convergence and stability properties of approximate solutions to the associated Riccati equations (see [14]). The ultimate result demonstrates that the regularized, finite-dimensional feedback control yields near optimal performance and that the corresponding Riccati solution satisfies all the desired convergence properties. The general theory is illustrated by an example of a boundary control problem associated with the Kirchoff plate model. Some numerical results are provided for the given example.     

On the computation of modules of long quadrilaterals

We study quadrilateralsQ which are given by two intervals on {:Im = 0} and {:Im = 1}, and two Jordan arcs ₁, ₂, in {:0 Im 1} connecting these two intervals. Many practical problems require the determination of the modulem(Q) ofQ, but ifQ is long, i.e., if

Compact disjointness preserving operators

Scientific, Technical, and Designing Union Lensistemotekhnika. Translated from Funktsional'nyi Analiz i Ego Prilozheniya, Vol. 26, No. 2, pp. 55–57, April–June, 1992.     

A note on treating a second order cone program as a special case of a semidefinite program

It is well known that a vector is in a second order cone if and only if its arrow matrix is positive semidefinite. But much less well-known is about the relation between a second order cone program (SOCP) and its corresponding semidefinite program (SDP). The correspondence between the dual problem of SOCP and SDP is quite direct and the correspondence between the primal problems is much more complicated. Given a SDP primal optimal solution which is not necessarily arrow-shaped, we can construct a SOCP primal optimal solution. The mapping from the primal optimal solution of SDP to the primal optimal solution of SOCP can be shown to be unique. Conversely, given a SOCP primal optimal solution, we can construct a SDP primal optimal solution which is not an arrow matrix. Indeed, in general no primal optimal solutions of the SOCP-related SDP can be an arrow matrix.Mathematics Subject Classification (2000): 20E28, 20G40, 20C20     

Analysis of a non-linear difference scheme in reaction-diffusion

Summary A nonlinear partial difference equation resulting from discretising in space and time the parabolic reaction diffusion equation, which models the spruce budworm problem, is analysed and accuracy estimates obtained for solutions over afinite time range and ast. Although the analysis is restricted to the logistic model in one space dimension, the techniques and comparison principles developed in the paper should prove useful in assessing the merits of numerical solutions of other nonlinear parabolic difference equations.During the period of this research Professor Guo Ben Yu was supported by a Science and Engineering Research Council visiting fellowship     

The classification of compact punctally cohesive Desarguesian projective Klingenberg planes

It is shown that for a compact Desarguesian projective Klingenberg plane P with incidence structure P=#x2119;, , I and neighbour relation , where two distinct points always lie on some line, exactly one of the following holds: P is a non-discrete connected or totally disconnected ordinary projective plane with =id, P is a finite projective plane with =id, P is a finite projective Hjelmslev plane with id, or P is a non-discrete totally disconnected ordinary projective plane with id.Dedicated to H. Salzmann on his 60th birthdayThe author wishes to thank the Natural Sciences and Engineering Research Council of Canada (NSERC) for their financial assistance in the writing of this paper.     

Singular Rank One Perturbations of Self-Adjoint Operators and Krein Theory of Self-Adjoint Extensions

Gesztesy and Simon recently have proven the existence of the strong resolvent limit A_, for A_, = A + (·), where A is a self-adjoint positive operator, being the A-scale). In the present note it is remarked that the operator A_, also appears directly as the Friedrichs extension of the symmetric operator :=A \{f (A)| f,=0\}. It is also shown that Krein's resolvents formula: (A_{_b,}-z)^-1 =(A-z)^-1+ (·, ) _z, with b=b-(1+z) (_z,_-1),_z= (A-z)^-1 defines a self-adjoint operator A_b, for each and b R¹. Moreover it is proven that for any sequence _n which goes to in there exists a sequence _n0 such that A_b, in the strong resolvent sense.     

A degenerate and strongly coupled quasilinear parabolic system not in divergence form

This paper deals with positive solutions of degenerate and strongly coupled quasi-linear parabolic system not in divergence form: u_t=v^p(u+au), v_t=u^q (v+bv) with null Dirichlet boundary condition and positive initial condition, where p, q, a and b are all positive constants, and p, q 1. The local existence of positive classical solution is proved. Moreover, it will be proved that: (i) When min {a, b} ₁ then there exists global positive classical solution, and all positive classical solutions can not blow up in finite time in the meaning of maximum norm (we can not prove the uniqueness result in general); (ii) When min {a, b} > ₁, there is no global positive classical solution (we can not still prove the uniqueness result), if in addition the initial datum (u₀v₀) satisfies u₀ + au₀ 0, v₀+bv₀ 0 in , then the positive classical solution is unique and blows up in finite time, where ₁ is the first eigenvalue of – in with homogeneous Dirichlet boundary condition.     

Sur les immeubles hyperboliques (On hyperbolic buildings)

The aim of this paper is to give a geometric approach to Tits' amalgam method to construct buildings and to initiate a study of hyperbolic buildings, i.e. whose types are reflexion systems of the real hyperbolic space. We construct lots of examples and study their cohomology at infinity. We construct CAT(–1) polyhedral complexes having big discrete parabolic groups of isometries.     

Two-Dimensional Discrete Hardy Inequalities

A sufficient condition on nonnegative double-sequences

On the relation between dynamic relaxation and semi-iterative matrix methods

Summary The well known relation between some matrix iterative methods for the solution of elliptic partial differential equations and time dependent parabolic equations is extended to establish a correspondence between the Chebyshev semi-iterative method and a time dependent hyperbolic equation. The method of dynamic relaxation based upon finding the transient solution to the corresponding dynamical system of equations is shown to be an exact analogue of the Chebyshev method provided the appropriate choice of inertia and damping coefficients is made. The possibility of developing further techniques based upon different physical models and leading to improved convergence is noted.     

On universality of blow-up profile for <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> critical nonlinear Schrödinger equation

We consider finite time blow-up solutions to the critical nonlinear Schrödinger equation iu_t=-u-|u|^4/Nu with initial condition u₀H¹. Existence of such solutions is known, but the complete blow-up dynamic is not understood so far. For a specific set of initial data, finite time blow-up with a universal sharp upper bound on the blow-up rate has been proved in [22], [23].We establish in this paper the existence of a universal blow-up profile which attracts blow-up solutions in the vicinity of blow-up time. Such a property relies on classification results of a new type for solutions to critical NLS. In particular, a new characterization of soliton solutions is given, and a refined study of dispersive effects of (NLS) in L² will remove the possibility of self similar blow-up in energy space H¹.     

Relativistic Euler equations for isentropic fluids: Stability of Riemann solutions with large oscillation

We analyze global entropy solutions of the 2 × 2 relativistic Euler equations for isentropic fluids in special relativity and establish the uniqueness of Riemann solutions in the class of entropy solutions in L     BV_loc with arbitrarily large oscillation. The uniqueness result does not require specific reference to any particular method for constructing the entropy solutions. Then the uniqueness of Riemann solutions implies their inviscid time-asymptotic stability under arbitrarily large L¹     L     BV_loc perturbation of the Riemann initial data, provided that the corresponding solutions are in L and have local bounded total variation that allows the linear growth in time. This approach is also extended to deal with the stability of Riemann solutions containing vacuum in the class of entropy solutions in L with arbitrarily large oscillation.     

Strong coupling limit for a certain class of polaron models

Summary It will be shown, that a Polaron-like asymptotic is given for a class of functions which are called of decreasing interaction and nicely approximable with respect to the Wiener process. Using Large Deviation techniques one can see, that the asymptotics for thestrong coupling limit for those Polaron models is given by a meanfield model and can be described by a variational problem.     

On Ovals of Riemann Surfaces of Even Genera

Let X be a Riemann surface of genus g2. A symmetry of of X is an antiholomorphic involution acting of X. A classical theorem of Harnack states that the set Fix () of fixed points of is either emplty or it consists of g+1 disjoint simple closed curves called, following Hilberts terminology, the ovals of . A Riemann surface admitting a symmetry corresponds to a real algebraic curve and nonconjugate symmetries correspond to different real models of the curve. The number of ovals of the symmetry equals the number of connected components of the corresponding real model. It is well known that two symmetries of a Riemann surface of genus g have at most 2g+2 ovals, and the bound is attained for every genus and just for commuting symmetries. Natanzon showed that three and four nonconjugate symmetries of a Riemann surface of genus g have at most 2g+4 and 2g+8 ovals respectively, and these bounds are attained for every odd genus and for commuting symmetries. Natanzon found that a Riemann surface of genus g has at most 2( +1) nonconjugate symmetries and, again, this bound is attained for infinitely many of g. Recently we have showed that a Riemann surface of even genus g admits at most four symmetries. Our aim here is to show, using NEC groups and combinatorial methods, that three nonconjugate symmetries of a surface of even genus g has at most 2g+3 ovals and, surprisingly, if such a surface admits four nonconjugate symmetries then its total number of ovals does not exceed 2g+2. Furthermore, we show that this last bound is sharp for every even genus g and for surfaces with automorphism group D _n × Z₂, for each n dividing 2g.     

Uniqueness theorem for measures in C(K) and its application in the theory of random processes

Let (, , P) be a probability space, let K be a separable Hausdorff topological space, and let :×K be a random process with continuous realizations on K. By the deviation of the random process from the function aC(K) we mean the random variable . It is proved in the paper that the random process is uniquely determined by the p-th moments of its deviations from continuous functions for a fixed p, p0, 2, 4, 6, .... The basis of the proof is the uniqueness theorem for measures in C(K), generalizing the uniqueness theorem for measures of order p on the line and refining the Hoffman-Jørgensen result on measures that coincide on the balls of the space C(K).Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 119, pp. 144–153, 1982.