Upwind difference approximations for degenerate parabolic convection-diffusion equations with a discontinuous coefficient

We analyse approximate solutions generated by an upwind differencescheme (of Engquist–Osher type) for nonlinear degenerateparabolic convection–diffusion equations where the nonlinearconvective flux function has a discontinuous coefficient (x)and the diffusion function A(u) is allowed to be strongly degenerate(the pure hyperbolic case is included in our setup). The mainproblem is obtaining a uniform bound on the total variationof the difference approximation u, which is a manifestationof resonance. To circumvent this analytical problem, we constructa singular mapping (, ·) such that the total variationof the transformed variable z = (, u) can be bounded uniformlyin . This establishes strong L¹ compactness of z and, since(, ·) is invertible, also u. Our singular mapping isnovel in that it incorporates a contribution from the diffusionfunction A(u). We then show that the limit of a converging sequenceof difference approximations is a weak solution as well as satisfyinga Krukov-type entropy inequality. We prove that the diffusionfunction A(u) is Hölder continuous, implying that the constructedweak solution u is continuous in those regions where the diffusionis nondegenerate. Finally, some numerical experiments are presentedand discussed.     

Long-Time Behavior of Solutions of the Fast Diffusion Equations with Critical Absorption Terms

This paper is devoted to the long-time behavior of solutionsto the Cauchy problem of the porous medium equation u_t = (u^m)– u^p in Rⁿ x (0,) with (1 – 2/n)₊ < m < 1and the critical exponent p = m + 2/n. For the strictly positiveinitial data u(x,0) = O(1 + |x|)^–k with n + mn(2 –n + nm)/(2[2 – m + mn(1 – m)]) k < 2/(1 –m), we prove that the solution of the above Cauchy problem convergesto a fundamental solution of u_t = (u^m) with an additional logarithmicanomalous decay exponent in time as t .     

Cohen-Macaulay Properties of Thom-Boardman Strata I: Morin's Ideal

Thom–Boardman strata ^I are fundamental tools in studyingsingularities of maps. The Zariski closures of the strata ^Iare components of the set of zeros of the ideals ^I defined by B. Morin using iterated jacobian extensions in his paper‘Calcul jacobien’ (Ann. Sci. École Norm.Sup.} 8 (1975) 1–98). In this paper, we consider the questionof when the Morin ideals ^I define Cohen–Macaulay spaces.We determine all I=(i₁...,i_k) such that ^I defines a Cohen–Macaulayspace alongthe stratum. 1991 Mathematics Subject Classification: 13D25, 14B05, 14M12, 58C25.     

On Lp boundedness of wave operators for 4-dimensional Schrodinger operators with threshold singularities

Let H=–+V(x) be a Schrödinger operator on L²(R⁴),H₀=–. Assume that |V(x)|+| V(x)|C x ^– for some>8. Let be the wave operators. It is known that W_± extend to bounded operators in L^p(R⁴)for all 1p, if 0 is neither an eigenvalue nor a resonance ofH. We show that if 0 is an eigenvalue, but not a resonance ofH, then the W_± are still bounded in L^p(R⁴) for all psuch that 4/3<p<4.     

On the existence and profile of nodal solutions for a two-dimensional elliptic problem with large exponent in nonlinearity

We study the existence of nodal solutions to the boundary valueproblem – u = |u|^{p – 1}u in a bounded, smooth domain in ², with homogeneous Dirichlet boundary condition, when pis a large exponent. We prove that, for p large enough, thereexist at least two pairs of solutions which change sign exactlyonce and whose nodal lines intersect the boundary of .     

Nodal Solutions of a p-Laplacian Equation

We prove that the p-Laplacian problem –_p u = f(x, u),with u on a bounded domain R^N, with p > 1 arbitrary, has a nodal solution providedthat f : x R R is subcritical, and f(x, t) / |t|^p –² is superlinear. Infinitely many nodal solutions are obtainedif, in addition, f(x, –t) = –f(x, t). 2000 MathematicsSubject Classification 35J20, 35J65, 58E05.     

On Hardy, BMO and Lipschitz spaces of invariantly harmonic functions in the unit ball

The invariantly harmonic functions in the unit ball Bⁿ in Cⁿare those annihilated by the Bergman Laplacian . The Poisson-Szegökernel P(z,) solves the Dirichlet problem for : if f C(Sⁿ),the Poisson-Szegö transform of f, where d is the normalized Lebesgue measure on Sⁿ,is the unique invariantly harmonic function u in Bⁿ, continuousup to the boundary, such that u=f on Sⁿ. The Poisson-Szegötransform establishes, loosely speaking, a one-to-one correspondencebetween function theory in Sⁿ and invariantly harmonic functiontheory in Bⁿ. When n 2, it is natural to consider on Sⁿ functionspaces related to its natural non-isotropic metric, for theseare the spaces arising from complex analysis. In the paper,different characterizations of such spaces of smooth functionsare given in terms of their invariantly harmonic extensions,using maximal functions and area integrals, as in the correspondingEuclidean theory. Particular attention is given to characterizationin terms of purely radial or purely tangential derivatives.The smoothness is measured in two different scales: that ofSobolev spaces and that of Lipschitz spaces, including BMO andBesov spaces. 1991 Mathematics Subject Classification: 32A35,32A37, 32M15, 42B25.     

Error Analysis of the Enthalpy Method for the Stefan Problem

In this paper an error bound is derived for a practical piecewiselinear finite-element approximation of an enthalpy formulationof the multidimensional Stefan problem with an implicit timediscretization. It is shown that if the time step t is O(h),then the error in the temperature measured in the L² norm isO(h).     

On Extensions of Myers' Theorem

Let M be a compact Riemannian manifold, and let h be a smoothfunction on M. Let p^h(x) = inf||–₁(Ric_x(,)–2Hess(h_x(,)).Here Ric_x denotes the Ricci curvature at x and Hess(h) is theHessian of h. Then M has finite fundamental group if ^h–p^h<0. Here ^h =:+2L_h is the Bismut-Witten Laplacian. This leadsto a quick proof of recent results on extension of Myers' theoremto manifolds with mostly positive curvature. There is also asimilar result for noncompact manifolds.     

The P-version Penalty Finite Element Method

A p-version penalty finite element method is used to solve themodel problem –u=f in , u=g on . Error estimates are derivedin H¹-norm. The p-version penalty method with extrapolationyields an approximate solution which converges at the optimalrate. Numerical results show the effectiveness of the p-versionpenalty method with extrapolation.     

Asymptotic behavior of positive solutions of some quasilinear elliptic problems

We discuss the asymptotic behavior of positive solutions ofthe quasilinear elliptic problem –_pu = a u^p–1–b(x)u^q, u| = 0, as q p – 1 + 0 and as q , via a scale argument.Here _p is the p-Laplacian with 1 < p and q > p –1. If p = 2, such problems arise in population dynamics. Ourmain results generalize the results for p = 2, but some technicaldifficulties arising from the nonlinear degenerate operator–_p are successfully overcome. As a by-product, we cansolve a free boundary problem for a nonlinear p-Laplacian equation.     

The Model Plasma Problem: Existence of a Solution Branch

We are interested in the model plasma problem –u = u⁺in ,u = –d on , _au⁺ dx=j where is a bounded domain in with boundary ; here, j isa given positive number, the function u and the positive number are the unknowns of the problem, and d is a real parameter.Using a variant of the implicit function theorem, we can provethe existence of a global solution branch parametrized by d.The method has the advantage that it can be used for analysingthe approximation of the above problem by a finite-element method.     

Kato Class Potentials for Higher Order Elliptic Operators

Our goal in this paper is to determine conditions on a potentialV which ensure that an operator such as H:=(–)^m+V (1) acting on L²(R^N) defines a semigroup in L^p(R^N) for various valuesof p including p=1. The operator is defined as a quadratic formsum. That is, we put for (all integrals are on R^N and are with respect to Lebesgue measure), and note thatthe closure of the form is non-negative and has domain equalto the Sobolev space W^m,2. We then assume that the potentialhas quadratic form bound less than 1 with respect to Q₀, anddefine This form is closed and is associated with a semibounded self-adjointoperator H in L² (see [17, p. 348; 5, Theorem 4.23]). One canthen ask whether the semigroup e^–Ht defined on L² fort0 is extendable to a strongly continuous one-parameter semigroupon L^p for other values of p, and if so whether one can describethe domain and spectrum of its generator.     

Logistic Type Equations on RN by a Squeezing Method Involving Boundary Blow-Up Solutions

We study, on the entire space R^N(N 1), the diffusive logisticequation u_t–u=u–u^p, u0 (1.1) and its generalizations. Here p > 1 is a constant. Problem(1.1) plays an important role in understanding various populationmodels and some other problems in applied mathematics. When = 1 and p = 2, it is also known as the Fisher equation andKPP equation, due to the pioneering works of Fisher [8] andKolmogoroff, Petrovsky and Piscounoff [18].     

Using random sets as oracles

Let R be a notion of algorithmic randomness for individual subsetsof . A set B is a base for R randomness if there is a Z _T Bsuch that Z is R random relative to B. We show that the basesfor 1-randomness are exactly the K-trivial sets, and discussseveral consequences of this result. On the other hand, thebases for computable randomness include every ₂⁰ set that isnot diagonally noncomputable, but no set of PA-degree. As aconsequence, an n-c.e. set is a base for computable randomnessif and only if it is Turing incomplete.     

Trace Inequalities of Sobolev Type in the Upper Triangle Case

We give a new non-capacitary characterization of positive Borelmeasures µ on Rⁿ such that the potential space I*L^p isimbedded in L^q(dµ) for $1qp+, that is, the trace inequality holds, for Riesz potentials I = (- )². A weak-type trace inequality is also characterized. A non-isotropic version on the unit sphere Sⁿ is studied,as well as the holomorphic case for Hardy–Sobolev spaces in the ball. 1991 MathematicsSubject Classification: primary 31C15, 42B20; secondary 32A35.     

Convergent and spurious solutions of nonlinear elliptic equations

In this paper we investigate finite element approximations ofnonlinear elliptic equations in three dimensions. By applyingand extending the results of Lopez-Marcos and Sanz-Serna, weprove that the finite element approximation on a mesh of sizeh, has a solution U_k which converges to an exact solution ofthe differential equation as h0. This solution is unique withina suitably defined stability ball Bh. For the particular nonlinearequation u + (u + u^p) we show that the size of B_h depends uponh only if p > 5 when it tends to zero as h 0. In this casewe prove the existence of spurious solutions V_h of the Galerkinapproximation which become unbounded in the maximum norm ash0. The stability ball B_h then acts to separate the convergentand the spurious solutions. We present the results of some numericalexperiments to substantiate our claims.     

Geometry of Critical Loci

Let :(Z,z)(U,0) be the germ of a finite (that is, proper with finite fibres)complex analytic morphism from a complex analytic normal surfaceonto an open neighbourhood U of the origin 0 in the complexplane C². Let u and v be coordinates of C² defined on U. Weshall call the triple (, u, v) the initial data. Let stand for the discriminant locus of the germ , that is,the image by of the critical locus of . Let ()_A be the branches of the discriminant locus at O whichare not the coordinate axes. For each A, we define a rational number d by where I(–, –) denotes the intersection number at0 of complex analytic curves in C². The set of rational numbersd, for A, is a finite subset D of the set of rational numbersQ. We shall call D the set of discriminantal ratios of the initialdata (, u, v). The interesting situation is when one of thetwo coordinates (u, v) is tangent to some branch of , otherwiseD = {1}. The definition of D depends not only on the choiceof the two coordinates, but also on their ordering. In this paper we prove that the set D is a topological invariantof the initial data (, u, v) (in a sense that we shall definebelow) and we give several ways to compute it. These resultsare first steps in the understanding of the geometry of thediscriminant locus. We shall also see the relation with thegeometry of the critical locus.     

On Positive Multipeak Solutions of a Nonlinear Elliptic Problem

In this paper we continue our investigation in [5, 7, 8] onmultipeak solutions to the problem –²u+u=Q(x)|u|^q–2u, xR^N, uH¹(R^N) (1.1) where = ^N_i=1²/x²_i is the Laplace operator in R^N, 2 < q < for N = 1, 2, 2 < q < 2N/(N–2) for N3, and Q(x)is a bounded positive continuous function on R^N satisfying thefollowing conditions. (Q₁) Q has a strict local minimum at some point x₀R^N, that is,for some > 0 Q(x)>Q(x₀) for all 0 < |x–x₀| < . (Q₂) There are constants C, > 0 such that |Q(x)–Q(y)|C|x–y| for all |x–x₀| , |y–y₀| . Our aim here is to show that corresponding to each strict localminimum point x₀ of Q(x) in R^N, and for each positive integerk, (1.1) has a positive solution with k-peaks concentratingnear x₀, provided is sufficiently small, that is, a solutionwith k-maximum points converging to x₀, while vanishing as 0 everywhere else in R^N.     

Error estimates for the approximate solutions of a nonlinear hyperbolic equation given by finite volume schemes

This paper is devoted to the study of an error estimate of thefinite volume, approximation to the solution u L(R_N x R) ofthe equation u_t + div(Vf(u)) = 0, where v is a vector functiondepending on time and space. A 'h' error estimate for an initialvalue in BV(R^N) is shown for a large variety of finite volumemonotonous flux schemes, with an explicit or implicit time discretization.For this purpose, the error estimate is given for the generalsetting of approximate entropy solutions, where the error isexpressed in terms of measures in R^N and R^N x R. The study ofthe implicit schemes involves the study of the existence anduniqueness of the approximate solution. The cases where an 'h'error estimate can be achieved are also discussed.