Universal self-similarity of porous media equation with absorption: the critical exponent case

In this paper we study the large time behavior of non-negative solutions to the Cauchy problem of u_t=Δu^m−u^q in R^N×(0,∞), where m>1 and q=q_c≡m+2/N is a critical exponent. For non-negative initial value u(x,0)=u₀(x)∈L¹(R^N), we show that the solution converges, if u₀(x)(1+|x|)^k is bounded for some k>N, to a unique fundamental solution of u_t=Δu^m, independent of the initial value, with additional logarithmic anomalous decay exponent in time as t→∞.     

Properties of non-simultaneous blow-up in heat equations coupled via different localized sources

This paper deals with u_t = Δu + u^m(x, t)e^pv(0,t), v_t = Δv + u^q(0, t)e^nv(x,t), subject to homogeneous Dirichlet boundary conditions. The complete classification on non-simultaneous and simultaneous blow-up is obtained by four sufficient and necessary conditions. It is interesting that, in some exponent region, large initial data u₀(v₀) leads to the blow-up of u(v), and in some betweenness, simultaneous blow-up occurs. For all of the nonnegative exponents, we find that u(v) blows up only at a single point if m > 1(n > 0), while u(v) blows up everywhere for 0 ? m ? 1 (n = 0). Moreover, blow-up rates are considered for both non-simultaneous and simultaneous blow-up solutions.     

On critical Fujita exponents for heat equations with nonlinear flux conditions on the boundary

We consider nonnegative solutions of initial-boundary value problems for parabolic equationsu _t=u_xx, u_t=(u^m)_xxand (m>1) forx>0,t>0 with nonlinear boundary conditions−u _x=u^p,−(u ^m)_x=u^pand forx=0,t>0, wherep>0. The initial function is assumed to be bounded, smooth and to have, in the latter two cases, compact support. We prove that for each problem there exist positive critical valuesp ₀,p_c(withp ₀<p_c)such that forp∃(0,p ₀],all solutions are global while forp∃(p₀,p_c] any solutionu≢0 blows up in a finite time and forp>p _csmall data solutions exist globally in time while large data solutions are nonglobal. We havep _c=2,p _c=m+1 andp _c=2m for each problem, whilep ₀=1,p ₀=1/2(m+1) andp ₀=2m/(m+1) respectively. This work was done during visits of the first author to Iowa State University and the Institute for Mathematics and its Applications at the University of Minnesota. The second author was supported in part by NSF Grant DMS-9102210.     

Blow-up problem for semilinear heat equation with absorption and a nonlocal boundary condition

In this paper we consider a semilinear parabolic equation u_t=Δu−c(x,t)u^p for (x,t)∈Ω×(0,∞) with nonlinear and nonlocal boundary condition u∣_{∂Ω×(0,∞)}=∫_Ωk(x,y,t)u^ldy and nonnegative initial data where p>0 and l>0. We prove some global existence results. Criteria on this problem which determine whether the solutions blow up in finite time for large or for all nontrivial initial data are also given.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

Blow-up results for a class of first-order nonlinear evolution inequalities

We investigate the non-existence of solutions to a class of evolution inequalities; in this case, as it happens in a relatively small number of blow-up studies, nonlinearities depend also on time-variable t and spatial derivatives of the unknown. The present results, which in great part do not require any assumption on the regularity of data, are completely new and shown with various applications. Some of these results referring to the problem u_t=Δu+a(x)|u|^p+λf(x) in R^N, t>0 include the non-existence results of positive global solutions obtained by Fujita and others when a≡1 and f≡0, Bandle-Levine and Levine-Meier when a≡|x|^m and f≡0, Pinsky when either f≡0 or f?0 and λ>0, Zhang and Bandle-Levine-Zhang when a≡1 and λ=1.     

Locally Lipschitz solutions of a single conservation law in several space variables

We consider the Cauchy problem for a single conservation law in several space variables. Letting u(x, t) denote the solution with initial data u₀, we state necessary and sufficient conditions on u₀ so that u(x, t) is locally Lipschitz continuous in the half space {t > 0}. These conditions allow for the preservation of smoothness of u₀ as well as for the smooth resolution of discontinuities in u₀. One consequence of our result is that u(x, t) cannot be locally Lipschitz unless u₀ has locally bounded variation. Another is that solutions which are bounded and locally Lipschitz continuous in {t > 0} automatically have boundary values u₀ at t = 0 in the sense that u(·, t) → u₀ in L_loc¹. Finally, we give an elementary proof that locally Lipschitz solutions satisfy Kruzkov's uniqueness condition.     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

Problem with periodicity conditions for a degenerating equation of mixed type

For the equation K(t)u _xx + u _tt − b ² K(t)u = 0 in the rectangular domain D = “(x, t)‖ 0 < x < 1, −α < t < β”, where K(t) = (sgnt)|t| ^m , m > 0, and b > 0, α > 0, and β > 0 are given real numbers, we use the spectral method to obtain necessary and sufficient conditions for the unique solvability of the boundary value problem u(0, t) = u(1, t), u _x (0, t) = u _x (1, t), −α ≤ t ≤ β, u(x, β) = φ(x), u(x,−α) = ψ(x), 0 ≤ x ≤ 1.     

Existence and decay rates of solutions to the generalized Burgers equation

In this paper we study the generalized Burgers equation u_t+(u²/2)_x=f(t)u_xx, where f(t)>0 for t>0. We show the existence and uniqueness of classical solutions to the initial value problem of the generalized Burgers equation with rough initial data belonging to , as well it is obtained the decay rates of u in L^p norm are algebra order for p∈[1,∞[.     

On the existence of shock fronts for a Burgers equation with a singular source term

In this paper we consider the Cauchy problem for the equation ∂u/∂t + u ∂u/∂x + u/x = 0 for x > 0, t ⩾ 0, with u(x, 0) = u⁰₋(x) for x < x₀, u(x, 0) = u⁰₊(x) for x > x₀, u⁰₋(x₀) > u⁰₊(x₀). Following the ideas of Majda, 1984 and Lax, 1973, we construct, for smooth u⁰₋ and u⁰₊, a global shock front weak solution u(x, t) = u₋(x, t) for x < ϕ(t), u(x, t) = u₊(x, t) for x > ϕ(t), where u₋ and u₊ are the strong solutions corresponding (respectively) to u⁰₋ and u⁰₊ and the curve t → ϕ(t) is defined by dϕ/dt (t) = 1/2[u₋(ϕ(t), t) + u₊(ϕ(t), t)], t ⩾ 0 and ϕ(0) = x₀. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

Complete quenching for singular parabolic problems

We prove finite time extinction of the solution of the equation ut−Δu+χ{u>0}(u−β−λf(u))=0 in Ω×(0,∞) with boundary data u(x,t)=0 on ∂Ω×(0,∞) and initial condition u(x,0)=u₀(x) in Ω, where Ω⊂RN is a bounded smooth domain, 0<β<1 and λ>0 is a parameter. For every small enough λ>0 there exists a time t₀>0 such that the solution is identically equal to zero.     

Existence and asymptotic behaviour of solutions of the very fast diffusion equation

Let n ≥ 3, 0 < m ≤ (n ? 2)/n, p > max(1, (1 ? m)n/2), and ${0 \le u_0 \in L_{loc}^p(\mathbb{R}^n)}$ satisfy ${{\rm lim \, inf}_{R\to\infty}R^{-n+\frac{2}{1-m}} \int_{|x|\le R}u_0\,dx = \infty}$ . We prove the existence of unique global classical solution of u _t = Δu ^m , u > 0, in ${\mathbb{R}^n \times (0, \infty), u(x, 0) = u_0(x)}$ in ${\mathbb{R}^n}$ . If in addition 0 < m < (n ? 2)/n and u ₀(x) ≈ A|x|^?q as |x| → ∞ for some constants A > 0, q < n/p, we prove that there exist constants α, β, such that the function v(x, t) = t ^α u(t ^β x, t) converges uniformly on every compact subset of ${\mathbb{R}^n}$ to the self-similar solution ψ(x, 1) of the equation with ψ(x, 0) = A|x|^?q as t → ∞. Note that when m = (n ? 2)/(n + 2), n ≥ 3, if ${g_{ij} = u^{\frac{4}{n+2}}\delta_{ij}}$ is a metric on ${\mathbb{R}^n}$ that evolves by the Yamabe flow ?g _ij /?t = ?Rg _ij with u(x, 0) = u ₀(x) in ${\mathbb{R}^n}$ where R is the scalar curvature, then u(x, t) is a global solution of the above fast diffusion equation.     

Convergence of the Dirichlet solutions of the very fast diffusion equation

For any −1<m<0, positive functions f, g and u₀≥0, we prove that under some mild conditions on f, g and u₀ as R→∞ the solution u^R of the Dirichlet problem u_t=(u^m/m)_xx in (−R,R)×(0,∞), u(R,t)=(f(t)|m|R)^1/m, u(−R,t)=(g(t)|m|R)^1/m for all t>0, u(x,0)=u₀(x) in (−R,R), converges uniformly on every compact subset of R×(0,T) to the solution of the equation u_t=(u^m/m)_xx in R×(0,T), u(x,0)=u₀(x) in R, which satisfies some mass loss formula on (0,T) where T is the maximal time such that the solution u is positive. We also prove that the solution constructed is equal to the solution constructed in Hui (2007) [15] using approximation by solutions of the corresponding Neumann problem in bounded cylindrical domains.     

Lower dimensional invariant tori with prescribed frequency for nonlinear wave equation

In this paper, one-dimensional (1D) nonlinear wave equation utt−uxx+mu+u³=0, subject to Dirichlet boundary conditions is considered. We show that for each given m>0, and each prescribed integer b>1, the above equation admits a Whitney smooth family of small-amplitude quasi-periodic solutions with b-dimensional Diophantine frequencies, which correspond to b-dimensional invariant tori of an associated infinite-dimensional dynamical system. In particular, these Diophantine frequencies are the small dilation of a prescribed Diophantine vector. The proof is based on a partial Birkhoff normal form reduction and an improved KAM method.     

Plane structures in thermal runaway

We consider the problem

1. u _t=u _xx+e ^u whenx ∈ ?,t > 0,
2. u(x, 0) =u ₀(x) whenx ∈ ?,

whereu ₀(x) is continuous, nonnegative and bounded. Equation (1) appears as a limit case in the analysis of combustion of a one-dimensional solid fuel. It is known that solutions of (1), (2) blow-up in a finite timeT, a phenomenon often referred to as thermal runaway. In this paper we prove the existence of blow-up profiles which are flatter than those previously observed. We also derive the asymptotic profile ofu(x, T) near its blow-up points, which are shown to be isolated.     

Invariant tori of nonlinear Schrödinger equation

In this paper, we consider one-dimensional nonlinear Schrödinger equation iut−uxx+V(x)u+f(²|u|)u=0 on [0,π]×R under the boundary conditions a₁u(t,0)−b₁ux(t,0)=0, a₂u(t,π)+b₂ux(t,π)=0, , for i=1,2. It is proved that for a prescribed and analytic positive potential V(x), the above equation admits small-amplitude quasi-periodic solutions corresponding to d-dimensional invariant tori of the associated infinite-dimensional dynamical system.     

Analysis of a semigroup approach in the inverse problem of identifying an unknown parameters

This article presents a semigroup approach to the mathematical analysis of the inverse parameter problems of identifying the unknown parameters p(t) and q in the linear parabolic equation u_t(x, t)  = u_xx + qu_x(x, t) + p(t)u(x, t), with Dirichlet boundary conditions u(0, t) = ψ₀, u(1, t) = ψ₁. The main purpose of this paper is to investigate the distinguishability of the input-output mapping Φ[·]:P→H^1,2[0,T], via semigroup theory. In this paper, it is shown that if the nullspace of the semigroup T(t) consists of only zero function, then the input-output mapping Φ[·] has the distinguishability property. It is also shown that the types of the boundary conditions and the region on which the problem is defined play an important role in the distinguishability property of the mapping. Moreover, under the light of the measured output data u_x(0, t) = f(t) the unknown parameter p(t) at (x, t) = (0, 0) and the unknown coefficient q are determined via the input data. Furthermore, it is shown that measured output data f(t) can be determined analytically by an integral representation. Hence the input-output mapping Φ[·]:P→H^1,2[0,T] is given explicitly interms of the semigroup.     

Numerical solutions with a priori error bounds for coupled self-adjoint time dependent partial differential systems

This paper is concerned with the construction of accurate continuous numerical solutions for partial self-adjoint differential systems of the type (P(t) u_t)_t = Q(t)u_xx, u(0, t) = u(d, t) = 0, u(x, 0) = f(x), u_t(x, 0) = g(x), 0 ≤ x ≤ d, t >- 0, where P(t), Q(t) are positive definite oR^r×r-valued functions such that P′(t) and Q′(t) are simultaneously semidefinite (positive or negative) for all t ≥ 0. First, an exact theoretical series solution of the problem is obtained using a separation of variables technique. After appropriate truncation strategy and the numerical solution of certain matrix differential initial value problems the following question is addressed. Given T > 0 and an admissible error ϵ > 0 how to construct a continuous numerical solution whose error with respect to the exact series solution is smaller than ϵ, uniformly in D(T) = {(x, t); 0 ≤ x ≤ d, 0 ≤ t ≤ T}. Uniqueness of solutions is also studied.     

On singular diffusion equations with applications to self-organized criticality

We consider solutions of the singular diffusion equation _t, = (u^m?1 u_x)_x, m ≦ 0, associated with the flux boundary condition lim_x→?∞ (u^m?1u_x)_x = λ > 0. The evolutions defined by this problem depend on both m and λ. We prove existence and stability of traveling wave solutions, parameterized by λ. Each traveling wave is stable in its appropriate evolution. These traveling waves are in L¹ for ?1 < m ≦ 0, but have non-integrable tails for m ≦ ?1. We also show that these traveling waves are the same as those for the scalar conservation law u_t = ?[f(u)]_x + u_xx for a particular nonlinear convection term f(u) = f(u;m, λ). © 1993 John Wiley & Sons, Inc.