A note on the behavior of blow-up solutions for one-phase Stefan problems

ＡＮＯＴＥＯＮＴＨＥＢＥＨＡＶＩＯＲＯＦＢＬＯＷ┐ＵＰＳＯＬＵＴＩＯＮＳＦＯＲＯＮＥ┐ＰＨＡＳＥＳＴＥＦＡＮＰＲＯＢＬＥＭＳＺＨＵＮＩＮＧＡｂｓｔｒａｃｔ．Ｉｎｔｈｉｓｐａｐｅｒ,ｔｈｅｆｏｌｏｗｉｎｇｏｎｅ－ｐｈａｓｅＳｔｅｆａｎｐｒｏｂｌｅｍｉｓ...     

Reaction diffusion equations with non-linear boundary conditions,blowup and steady states

In this paper we study the following problem: u_t−Δu=−f(u) in Ω×(0, T)≡Q_T, ∂u ∂n=g(u) on ∂Ω×(0, T)≡S_T, u(x, 0)=u₀(x) in Ω , where Ω⊂ℝ^N is a smooth bounded domain, f and g are smooth functions which are positive when the argument is positive, and u₀(x)>0 satisfies some smooth and compatibility conditions to guarantee the classical solution u(x, t) exists. We first obtain some existence and non-existence results for the corresponding elliptic problems. Then, we establish certain conditions for a finite time blow-up and global boundedness of the solutions of the time-dependent problem. Further, we analyse systems with same kind of boundary conditions and find some blow-up results. In the last section, we study the corresponding elliptic problems in one-dimensional domain. Our main method is the comparison principle and the construction of special forms of upper–lower solutions using related equations. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

The Cauchy problem for a linear parabolic partial differential equation

Numerical approximation of the solution of the Cauchy problem for the linear parabolic partial differential equation is considered. The problem: (p(x)u_x)_x ? q(x)u = p(x)u_t, 0 < x < 1,0 < t? T; $\text{u(0, t) = ?}{}_1\text{(t), 0 < t ? T}$; $\text{u(1,t) = ?}{}_2\text{(t), 0 < t ? T}$; p(0) u_x(0, t) = g(t), 0 < t₀ ? t ? T, is ill-posed in the sense of Hadamard. Complex variable and Dirichlet series techniques are used to establish Hölder continuous dependence of the solution upon the data under the additional assumption of a known uniform bound for ¦ u(x, t)¦ when 0 ? x ? 1 and 0 ? t ? T. Numerical results are obtained for the problem where the data ?₁, ?₂ and g are known only approximately.     

A Threshold Result for a Non-local Parabolic Equation

In this paper, we consider the Cauchy problem: (ECP) u_t−Δu+p(x)u=u(x,t)∫u²(y,t)/∣x−y∣dy; x∈ℝ³, t>0, u(x, 0)=u₀(x)⩾0 x∈ℝ³, (0.2) The stationary problem for (ECP) is the famous Choquard–Pekar problem, and it has a unique positive solution ū(x) as long as p(x) is radial, continuous in ℝ³, p(x)⩾ā>0, and lim_{∣x∣→∞}p(x)=p¯>0. In this paper, we prove that if the initial data 0⩽u₀(x)⩽(≢)ū(x), then the corresponding solution u(x, t) exists globally and it tends to the zero steady-state solution as t→∞, if u₀(x)⩾(≢)ū(x), then the solution u(x,t) blows up in finite time. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Removable singularities and non-uniqueness of solutions of a singular diffusion equation

 We prove that the solution u of the equation u _t =Δlog u, u>0, in (Ω\{x ₀})×(0,T), Ω⊂ℝ², has removable singularities at {x ₀}×(0,T) if and only if for any 0<α<1, 0<a<b<T, there exist constants ρ₀, C ₁, C ₂>0, such that C ₁ |x−x ₀|^α≤u(x,t)≤C ₂|x−x ₀|^−α holds for all 0<|x−x ₀|≤ρ₀ and a≤t≤b. As a consequence we obtain a sufficient condition for removable singularities at {∞}×(0,T) for solutions of the above equation in ℝ²×(0,T) and we prove the existence of infinitely many finite mass solutions for the equation in ℝ²×(0,T) when 0≤u ₀∉L ¹ (ℝ²) is radially symmetric and u ₀L _loc ¹(ℝ²). Received: 16 December 2001 / Revised version: 20 May 2002 / Published online: 10 February 2003 Mathematics Subject Classification (1991): 35B40, 35B25, 35K55, 35K65     

Exact simulation of final, minimal and maximal values of Brownian motion and jump-diffusions with applications to option pricing

We introduce a method for generating (W_x,T^(m,s),m_x,T^(m,s),M_x,T^(m,s))(W_{x,T}^{(\mu,\sigma)},m_{x,T}^{(\mu,\sigma)},M_{x,T}^{(\mu,\sigma)}) , where W_x,T^(m,s)W_{x,T}^{(\mu,\sigma)} denotes the final value of a Brownian motion starting in x with drift μ and volatility σ at some final time T, m_x,T^(m,s) = inf_{0 ￡ t ￡ T}W_x,t^(m,s)m_{x,T}^{(\mu,\sigma)} = {\rm inf}_{0\leq t \leq T}W_{x,t}^{(\mu,\sigma)} and M_x,T^(m,s) = sup_{0 ￡ t ￡ T} W_x,t^(m,s)M_{x,T}^{(\mu,\sigma)} = {\rm sup}_{0\leq t \leq T} W_{x,t}^{(\mu,\sigma)} . By using the trivariate distribution of (W_x,T^(m,s),m_x,T^(m,s),M_x,T^(m,s))(W_{x,T}^{(\mu,\sigma)},m_{x,T}^{(\mu,\sigma)},M_{x,T}^{(\mu,\sigma)}) , we obtain a fast method which is unaffected by the well-known random walk approximation errors. The method is extended to jump-diffusion models. As sample applications we include Monte Carlo pricing methods for European double barrier knock-out calls with continuous reset conditions under both models. The proposed methods feature simple importance sampling techniques for variance reduction.     

The zero dispersion limit of the korteweg-de vries equation for initial potentials with non-trivial reflection coefficient

The inverse scattering method is used to determine the distribution limit as ? → 0 of the solution u(x, t, ?) of the initial value problem. U_t ? 6uu_x + ?²u_xxx = 0, u(x, 0) = v(x), where v(x) is a positive bump which decays sufficiently fast as x x→±α. The case v(x) ? 0 has been solved by Peter D. Lax and C. David Levermore [8], [9], [10]. The computation of the distribution limit of u(x, t, ?) as ? → 0 is reduced to a quadratic maximization problem, which is then solved.     

Inverse Problems for Time-Dependent Singular Heat Conductivities: Multi-Dimensional Case

We consider an inverse boundary value problem for the heat equation ? _t u = div (γ? _x u) in (0, T) × Ω, u = f on (0, T) × ?Ω, u| _t=0 = u ₀, in a bounded domain Ω ? ? ⁿ , n ≥ 2, where the heat conductivity γ(t, x) is piecewise constant and the surface of discontinuity depends on time: γ(t, x) = k ² (x ∈ D(t)), γ(t, x) = 1 (x ∈ Ω?D(t)). Fix a direction e* ∈ 𝕊 ^n?1 arbitrarily. Assuming that ?D(t) is strictly convex for 0 ≤ t ≤ T, we show that k and sup {e*·x; x ∈ D(t)} (0 ≤ t ≤ T), in particular D(t) itself, are determined from the Dirichlet-to-Neumann map : f → ?_ν u(t, x)|_{(0, T)×?Ω}. The knowledge of the initial data u ₀ is not used in the proof. If we know min_0≤t≤T (sup _x∈D(t) x·e*), we have the same conclusion from the local Dirichlet-to-Neumann map. Numerical examples of stationary and moving circles inside the unit disk are shown. The results have applications to nondestructive testing. Consider a physical body consisting of homogeneous material with constant heat conductivity except for a moving inclusion with different conductivity. Then the location and shape of the inclusion can be monitored from temperature and heat flux measurements performed at the boundary of the body. Such a situation appears for example in blast furnaces used in ironmaking.     

Global stabilization of nonlinear systems by quadratic Lyapunov functions

The system x = A (t, x)x + B(t, x)u, where A(t, x) and B(t, x) are, respectively, n × n and n × m (m<n) continuous matrices whose elements are uniformly bounded for t ≽ t ₀ and x ∈ ℝ ⁿ , is considered. It is assumed that the system has relative degree q = n - m + 1, and the determinant of the matrix composed of the last m rows of the matrix B(t, x) is bounded away from zero for t ≽ t ₀ and x ∈ ℝ ⁿ . A special quadratic Lyapunov function with constant positive definite coefficient matrix H depending only on the range of variation of the coefficients in the matrices A(t, x) and B(t, x) is constructed and applied to obtain a control u(t, x) =7n ~B⋆ (t, x)H depending on a scalar parameter 7n under which the system is globally asymptotically stable provided that it is closed. Here, ~B (t, x) is the scalar matrix obtained from the matrix B(t, x) by setting the first n - m rows to zero.     

On a smooth solution of a nonlinear periodic boundary-value problem

We establish conditions for the existence of a smooth solution of a quasilinear hyperbolic equationu _tt - u_xx = ƒ(x, t, u, u, u _x),u (0,t) = u (π,t) = 0,u (x, t+ T) = u (x, t), (x, t) ∈ [0, π] ×R, and prove a theorem on the existence and uniqueness of a solution. Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 51, No. 11, pp. 1574–1576, November, 1999.     

An asymptotic of the negative discrete spectrum of the Schrödinger operator

The Schrödinger operator Hu = -Δu + V(x)u, where V(x) → 0 as ¦x¦ → ∞, is considered in L₂(R^m) for m?3. The asymptotic formula $$N(\lambda ,V) \sim \Upsilon _m \int {(\lambda - V(x))_ + ^{{m \mathord{\left/ {\vphantom {m {2_{dx} }}} \right. \kern-\nulldelimiterspace} {2_{dx} }}} ,} \lambda \to ---0,$$ is established for the number N(λ, V) of the characteristic values of the operator H which are less than λ. It is assumed about the potential V that V = V_o + V₁; V_o < 0, ¦V_o =o (¦V_o¦^3/2) as ¦x¦ → ∞; σ (t/2, V_o) ?cσ (t. V_o) and V₁∈L_m/2,loc, σ(t, V₁) =o (σ (t, V_o)), where σ (t,f)= mes {x:¦f (x) ¦ > t).     

Quenching phenomena for a non-local diffusion equation with a singular absorption

In this paper we study the quenching problem for the non-local diffusion equation

The radial basis functions method for identifying an unknown parameter in a parabolic equation with overspecified data

Parabolic partial differential equations with overspecified data play a crucial role in applied mathematics and engineering, as they appear in various engineering models. In this work, the radial basis functions method is used for finding an unknown parameter p(t) in the inverse linear parabolic partial differential equation u_t = u_xx + p(t)u + φ, in [0,1] × (0,T], where u is unknown while the initial condition and boundary conditions are given. Also an additional condition ∫₀¹k(x)u(x,t)dx = E(t), 0 ≤ t ≤ T, for known functions E(t), k(x), is given as the integral overspecification over the spatial domain. The main approach is using the radial basis functions method. In this technique the exact solution is found without any mesh generation on the domain of the problem. We also discuss on the case that the overspecified condition is in the form ∫₀^s(t) u(x,t)dx = E(t), 0 < t ≤ T, 0 < s(t) < 1, where s and E are known functions. Some illustrative examples are presented to show efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On Osgood theorem in Banach spaces

Let X be a real Banach space, ω : [0, +∞) → ? be an increasing continuous function such that ω(0) = 0 and ω(t + s) ≤ ω(t) + ω(s) for all t, s ∈ [0, +∞). According to the infinite dimensional analog of the Osgood theorem if ∫¹₀ (ω(t))^?1 dt = ∞, then for any (t₀, x₀) ∈ ?×X and any continuous map f : ?×X → X such that ∥f(t, x) – f(t, y)∥ ≤ ω(∥x – y∥) for all t ∈ ?, x, y ∈ X, the Cauchy problem (t) = f(t, x(t)), x(t₀) = x₀ has a unique solution in a neighborhood of t₀. We prove that if X has a complemented subspace with an unconditional Schauder basis and ∫¹₀ (ω(t))^?1 dt < ∞ then there exists a continuous map f : ? × X → X such that ∥f(t, x) – f(t, y)∥ ≤ ω(∥x – y∥) for all (t, x, y) ∈ ? × X × X and the Cauchy problem (t) = f(t, x(t)), x(t₀) = x₀ has no solutions in any interval of the real line.     

Solution of Initial Value Problems with Monogenic Initial Functions in Banach Spaces with L p -Norm

This paper deals with the initial value problem of the type

New Logarithmic Sobolev Inequalities and an ɛ‐Regularity Theorem for the Ricci Flow

In this note, we prove an ?‐regularity theorem for the Ricci flow. Let (Mⁿ,g(t)) with t ? [?T,0] be a Ricci flow, and let H_x0(y,s) be the conjugate heat kernel centered at some point (x₀,0) in the final time slice. By substituting H_x0(?,s) into Perelman's W‐functional, we obtain a monotone quantity W_x0(s) that we refer to as the pointed entropy. This satisfies W_x0(s) ≤ 0, and W_x0(s) = 0 if and only if (Mⁿ,g(t)) is isometric to the trivial flow on Rⁿ. Then our main theorem asserts the following: There exists ? > 0, depending only on T and on lower scalar curvature and μ‐entropy bounds for the initial slice (Mⁿ,g(?T)) such that W_x0(s) ≥ ?? implies |Rm| ≤ r^?2 on P_{? r}(x₀,0), where r² ≡ |s| and P_ρ(x,t) ≡ B_ρ(x,t) × (t?ρ²,t] is our notation for parabolic balls. The main technical challenge of the theorem is to prove an effective Lipschitz bound in x for the s‐average of W_x(s). To accomplish this, we require a new log‐Sobolev inequality. Perelman's work implies that the metric measure spaces (Mⁿ,g(t),dvol_g(t)) satisfy a log‐Sobolev; we show that this is also true for the heat kernel weighted spaces (Mⁿ,g(t),H_x0(?,t)dvol_g(t)). Our log‐Sobolev constants for these weighted spaces are in fact universal and sharp. The weighted log‐Sobolev has other consequences as well, including certain average Gaussian upper bounds on the conjugate heat kernel. © 2014 Wiley Periodicals, Inc.     

Identification of the unknown diffusion coefficient in a quasi‐linear parabolic equation by semigroup approach with mixed boundary conditions

In this article, a semigroup approach is presented for the mathematical analysis of the inverse coefficient problems of identifying the unknown diffusion coefficient k(u(x, t)) in the quasi‐linear parabolic equation u_t(x, t)=(k(u(x, t))u_x(x, t))_x, with Dirichlet boundary conditions u_x(0, t)=ψ₀, u(1, t)=ψ₁. The main purpose of this work is to analyze the distinguishability of the input–output mappings Φ[·] : ??→C¹[0, T], Ψ[·] : ??→C¹[0, T] using semigroup theory. In this article, it is shown that if the null space of semigroups T(t) and S(t) consists of only a zero function, then the input–output mappings Φ[·] and Ψ[·] have the distinguishability property. Copyright © 2008 John Wiley & Sons, Ltd.     

Compactness of support of solutions to nonlinear second-order elliptic and parabolic equations in a half-cylinder

We study equations of the form $$\begin{gathered} u_{tt} + Lu + b(x,t)u_t = a(x,t)\left| u \right|^{\sigma - 1} u, \hfill \\ - u_t + Lu = a(x,t)\left| u \right|^{\sigma - 1} u \hfill \\ \end{gathered}$$ , whereL is a uniformly elliptic operator and 0<σ<1. In the half-cylinder II_0,∞={(x, t):x= (x ₁,...,x _n )∈ ω,t>0}, where ? ? ? _n is a bounded domain, we consider solutions satisfying the homogeneous Neumann condition forx∈?ω andt>0. We find conditions under which these solutions have compact support and prove statements of the following type: ifu(x, t)=o(t ^γ) ast→∞, then there exists aT such thatu(x, t)≡0 fort>T. In this case γ depends on the coefficients of the equation and on the exponent σ.