On Nonexistence of type II blowup for a supercritical nonlinear heat equation

In this paper we study blowup of radially symmetric solutions of the nonlinear heat equation u_t = Δu + |u|^p?1u either on ?^N or on a finite ball under the Dirichlet boundary conditions. We assume that the exponent p is supercritical in the Sobolev sense, that is, We prove that if p_s < p < p^*, then blowup is always of type I, where p^* is a certain (explicitly given) positive number. More precisely, the rate of blowup in the L^∞ norm is always the same as that for the corresponding ODE dv/dt = |v|^p?1v. Because it is known that “type II” blowup (or, equivalently, “fast blowup”) can occur if p > p^*, the above range of exponent p is optimal. We will also derive various fundamental estimates for blowup that hold for any p > p_s and regardless of type of blowup. Among other things we classify local profiles of type I and type II blowups in the rescaled coordinates. We then establish useful estimates for the so‐called incomplete blowup, which reveal that incomplete blowup solutions belong to nice function spaces even after the blowup time. © 2004 Wiley Periodicals, Inc.     

Continuation of blowup solutions of nonlinear heat equations in several space dimensions

The possible continuation of solutions of the nonlinear heat equation in R^N × R₊ u_t = Δu^m + u^p with m > 0, p > 1, after the blowup time is studied and the different continuation modes are discussed in terms of the exponents m and p. Thus, for m + p ≤ 2 we find a phenomenon of nontrivial continuation where the region {x : u(x, t) = ∞} is bounded and propagates with finite speed. This we call incomplete blowup. For N ≥ 3 and p > m(N + 2)/(N − 2) we find solutions that blow up at finite t = T and then become bounded again for t > T. Otherwise, we find that blowup is complete for a wide class of initial data. In the analysis of the behavior for large p, a list of critical exponents appears whose role is described. We also discuss a number of related problems and equations. We apply the same technique of analysis to the problem of continuation after the onset of extinction, for example, for the equation u_t = Δu^m − u^p, m > 0. We find that no continuation exists if p + m ≤ 0 (complete extinction), and there exists a nontrivial continuation if p + m > 0 (incomplete extinction). © 1997 John Wiley & Sons, Inc.     

Stability and localization of unbounded solutions of a nonlinear heat equation in a plane

The article investigates unbounded solutions of the equation u _t = div (u ^σgrad u) + u ^β in a plane. We numerically analyze the stability of two-dimensional self-similar solutions (structures) that increase with blowup. We confirm structural stability of the simple structure with a single maximum and metastability of complex structures. We prove structural stability of the radially symmetrical structure with a zero region at the center and investigate its attraction region. We study the effect of various perturbations of the initial function on the evolution of self-similar solutions. We further investigate how arbitrary compact-support initial distributions attain the self-similar mode, including distributions whose support is different from a disk. We show that the self-similar mode described by a simple radially symmetrical structure is achieved only in the central region, while the entire localization region does not have enough time to transform into a disk during blowup. We show for the first time that simple structures may merge into a complex structure, which evolves for a long time according to self-similar law.     

A note on gradient blowup rate of the inhomogeneous hamilton-jacobi equations

The gradient blowup of the equation u_t = Δu + a(x)|∇u|^p + h(x), where p > 2, is studied. It is shown that the gradient blowup rate will never match that of the self-similar variables. The exact blowup rate for radial solutions is established under the assumptions on the initial data so that the solution is monotonically increasing in time.     

Blowup rate for heat equation in Lipschitz domains with nonlinear heat source terms on the boundary

This paper establishes the blowup rate estimate near the blowup time for the heat equation u_t=Δu with the nonlinear boundary condition u_n=u^p on the boundary ∂Ω×[0,T) for a bounded Lipschitz domain and for the nonlinear boundary condition u_n=exp(u) on the boundary for a two-dimensional disc with nonradially symmetric data.     

Solving the equation −uxx − ϵuyy = f(x,y, u) by an O(h4) finite difference method

The semi‐linear equation −u_xx − ϵu_yy = f(x, y, u) with Dirichlet boundary conditions is solved by an O(h⁴) finite difference method, which has local truncation error O(h²) at the mesh points neighboring the boundary and O(h⁴) at most interior mesh points. It is proved that the finite difference method is O(h⁴) uniformly convergent as h → 0. The method is considered in the form of a system of algebraic equations with a nine diagonal sparse matrix. The system of algebraic equations is solved by an implicit iterative method combined with Gauss elimination. A Mathematica module is designed for the purpose of testing and using the method. To illustrate the method, the equation of twisting a springy rod is solved. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 395–407, 2000     

Thermal avalanche for blowup solutions of semilinear heat equations

We consider the semilinear heat equation u_t = Δu + u^p both in ?^N and in a bounded domain with homogeneous Dirichlet boundary conditions, with 1 < p < p_s where p_s is the Sobolev exponent. This problem has solutions with finite‐time blowup; that is, for large enough initial data there exists T < ∞ such that u is a classical solution for 0 < t < T, while it becomes unbounded as t ↗ T. In order to understand the situation for t > T, we consider a natural approximation by reaction problems with global solution and pass to the limit. As is well‐known, the limit solution undergoes complete blowup: after it blows up at t = T, the continuation is identically infinite for all t > T. We perform here a deeper analysis of how complete blowup occurs. Actually, the singularity set of a solution that blows up as t ↗ T propagates instantaneously at time t = T to cover the whole space, producing a catastrophic discontinuity between t = T^? and t = T⁺. This is called the “avalanche.” We describe its formation as a layer appearing in the limit of the natural approximate problems. After a suitable scaling, the initial structure of the layer is given by the solution of a limit problem, described by a simple ordinary differential equation. As t proceeds past T, the solutions of the approximate problems have a traveling wave behavior with a speed that we compute. The situation in the inner region depends on the type of approximation: a conical pattern is formed with time when we approximate the power u^p by a flat truncation at level n, while for tangent truncations we get an exponential increase in time and a diffusive spatial pattern. © 2003 Wiley Periodicals, Inc.     

Stability of Cahn‐Hilliard fronts

We prove stability of the kink solution of the Cahn‐Hilliard equation ∂_tu = ∂( ∂u − u/2 + u³/2), x ∈ ℝ. The proof is based on an inductive renormalization group method, and we obtain detailed asymptotics of the solution as t → ∞. We prove stability of the kink solution of the Cahn‐Hilliard equation ∂_tu = ∂( ∂u − u/2 + u³/2), x ∈ ℝ. The proof is based on an inductive renormalization group method, and we obtain detailed asymptotics of the solution as t → ∞. © 1999 John Wiley & Sons, Inc.     

Blowup for nonlinear wave equations describing boson stars

We consider the nonlinear wave equation modeling the dynamics of (pseudorelativistic) boson stars. For spherically symmetric initial data, u₀(x) ∈ C (?³), with negative energy, we prove blowup of u(t, x) in the H^1/2‐norm within a finite time. Physically this phenomenon describes the onset of “gravitational collapse” of a boson star. We also study blowup in external, spherically symmetric potentials, and we consider more general Hartree‐type nonlinearities. As an application, we exhibit instability of ground state solitary waves at rest if m = 0. © 2007 Wiley Periodicals, Inc.     

Blowup of solutions for the “bad” Boussinesq-type equation

The paper studies the blowup of solutions to the initial boundary value problem for the “bad” Boussinesq-type equation u_tt−u_xx−bu_xxxx=σ(u)_xx, where b>0 is a real number and σ(s) is a given nonlinear function. By virtue of the energy method and the Fourier transform method, respectively, it proves that under certain assumptions on σ(s) and initial data, the generalized solutions of the above-mentioned problem blow up in finite time. And a few examples are shown, especially for the “bad” Boussinesq equation, two examples of blowup of solutions are obtained numerically.     

Stable Blowup Dynamics for the 1‐Corotational Energy Critical Harmonic Heat Flow

We exhibit a stable finite time blowup regime for the 1‐corotational energy critical harmonic heat flow from ?² into a smooth compact revolution surface of ?³ that reduces to the semilinear parabolic problem for a suitable class of functions f. The corresponding initial data can be chosen smooth, well localized, and arbitrarily close to the ground state harmonic map in the energy‐critical topology. We give sharp asymptotics on the corresponding singularity formation that occurs through the concentration of a universal bubble of energy at the speed predicted by van den Berg, Hulshof, and King. Our approach lies in the continuation of the study of the 1‐equivariant energy critical wave map and Schrödinger map with ??² target by Merle, Raphaël, and Rodnianski. © 2012 Wiley Periodicals, Inc.     

The focusing problem for the Eikonal equation

We study the focusing problem for the eikonal equation¶¶ ?_tu=| ?u| ², \partial _{t}u=\left| \nabla u\right| ^{2}, ¶¶i.e., the initial value problem in which the support of the initial datum is outside some compact set in R^d \mathbf{R}^{d} . The hole in the support will be filled in finite time and we are interested in the asymptotics of the hole as it closes. We show that in the radially symmetric case there are self-similar asymptotics, while in the absence of radial symmetry essentially any convex final shape is possible. However in R² \mathbf{R}^2 , for generic initial data the asymptotic shape will be either a vanishing triangle or the region between two parabolas moving in opposite directions (a closing eye). We compare these results with the known results for the porous medium pressure equation which approaches the eikonal equation in the limit as m? 1 m\rightarrow 1 .     

Convexity and all‐time C∞‐regularity of the interface in flame propagation

We consider the following one‐phase free boundary problem: Find (u, Ω) such that Ω = {u > 0} and with Q_T = ?ⁿ × (0, T). Under the condition that Ω_o is convex and log u_o is concave, we show that the convexity of Ω(t) and the concavity of log u(·, t) are preserved under the flow for 0 ≤ t ≤ T as long as ?Ω(t) and u on Ω(t) are smooth. As a consequence, we show the existence of a smooth‐up‐to‐the‐interface solution, on 0 < t < T_c, with T_c denoting the extinction time of Ω(t). We also provide a new proof of a short‐time existence with C^2,α initial data on the general domain. © 2002 John Wiley & Sons, Inc.     

Solutions of the Differential Equation u‴+λ2zu+(α−1)λ2u=0

This paper deals with the solutions of the differential equation u?+λ²zu+(α?1)λ²u=0, in which λ is a complex parameter of large absolute value and α is an arbitrary constant, real or complex. After a discussion of the structure of the solutions of the differential equation, an integral representation of the solution is given, from which the series solutions and their asymptotic representations are derived. A third independent solution is needed for the special case when α?1 is a positive integer, and two derivations for this are given. Finally, a comparison is made with the results obtained by R. E. Langer.     

On –δu = K(x)u5 in ℝ3

In this paper, we study the equation –Δu = K(x)u⁵ in ?³ and provide a large class of positive functions K(x) for which we obtain infinitely many positive solutions which decay at infinity at the rate of |x|^?1. © 1993 John Wiley & Sons, Inc.     

Finite‐difference approximation for the u(k)‐derivative with O(hM−k+1) accuracy: An analytical expression

An approximation of function u(x) as a Taylor series expansion about a point x₀ at M points x_i, ～ i = 1,2,…,M is used where x_i are arbitrary‐spaced. This approximation is a linear system for the derivatives u^(k) with an arbitrary accuracy. An analytical expression for the inverse matrix A ^?1 where A = [A_ik] = (x_i ? x₀)^k is found. A finite‐difference approximation of derivatives u^(k) of a given function u(x) at point x₀ is derived in terms of the values u(x_i). © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

The finite-time blowup of the solution of an initial boundary-value problem for the nonlinear equation of ion sound waves

We obtain blowup conditions for the solutions of initial boundary-value problems for the nonlinear equation of ion sound waves in a hydrogen plasma in the approximation of “hot” electrons and “heavy” ions. A specific characteristic of this nonlinear equation is the noncoercive nonlinearity of the form ?_t|?u|², which complicates its study by any energy method. We solve this problem by the Mitidieri–Pohozaev method of nonlinear capacity.     

Certain classes of potentials for p ‐Laplacian to be non‐degenerate

Given a positive integer n and an exponent 1 ≤ α ≤ ∞. We will find explicitly the optimal bound r_n such that if the L^α norm of a potential q (t ) satisfies ‖q ‖equation/tex2gif-inf-2.gif < r_n then the n ^th Dirichlet eigenvalue of the onedimensional p ‐Laplacian with the potential q (t ): (|u ′|^{p –2} u ′)′ + (λ + q (t )) |u |^{p –2}u = 0 (1 < p < ∞) will be positive. Using these bounds, we will construct, for the Dirichlet, the Neumann, the periodic or the antiperiodic boundary conditions, certain classes of potentials q (t ) so that the p ‐Laplacian with the potential q (t ) is non‐degenerate, which means that the above equation with λ = 0 has only the trivial solution verifying the corresponding boundary condition. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic behavior at infinity for solutions of Emden-Fowler type equations

The semilinear equation Δu = |u|^σ?1 u is considered in the exterior of a ball in ? ⁿ , _n ≥ 3. It is shown that if the exponent σ is greater than a “critical” value (= n/n?2), then for x → ∞ the leading term of the asymptotics of any solution is a linear combination of derivatives of the fundamental solution. It is shown that there exist solutions with the indicated leading term of an asymptotics of such a type.     

Principal Eigenvalues and Anti – Maximum Principle for Some Quasilinear Elliptic Equations on ℝN

We improve some previous existence and nonexistence results for positive principal eigenvalues of the problem —Δ_pu = λg(x)ψ^p(u), x ∈ ℝ^N, lim_{‖x‖⇒+∞}u(x) = 0. Also we discuss existence, nonexistence and antimaximum principle questions concerning the perturbed problem —Δ_pu = λg(x)ψ^p(u) + f(x), x∈ ℝ^N.