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.     

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.     

Control parameter estimation in a semi-linear parabolic inverse problem using a high accurate method

Parabolic inverse problems have an important role in many branches of science and technology. The aim of this research work is to solve these classes of equations using a high order compact finite difference scheme. We consider the following inverse problem for finding u(x, t) and p(t) governed by u_t = u_xx + p(t)u + φ(x, t) with an over specified condition inside the domain. Spatial derivatives are approximated using central difference scheme. The time advancement of the simulation is performed using a “third order compact Runge-Kutta method”. The convergence orders for the approximation of both u and p are of o(k³ + h²) which improves the results obtained in the literature. An exact test case is used to evaluate the validity of our numerical analysis. We found that the accuracy of the results is better than that of previous works in the literature.     

An inverse problem for computing a leading coefficient in the Sturm-Liouville operator by using the boundary data

We consider an inverse problem for identifying a leading coefficient α(x) in −(α(x)y′(x))′ + q(x)y(x) = H(x), which is known as an inverse coefficient problem for the Sturm-Liouville operator. We transform y(x) to u(x, t) =  (1 + t)y(x) and derive a parabolic type PDE in a fictitious time domain of t. Then we develop a Lie-group adaptive method (LGAM) to find the coefficient function α(x). When α(x) is a continuous function of x, we can identify it very well, by giving boundary data of y, y′ and α. The efficiency of LGAM is confirmed by comparing the numerical results with exact solutions. Although the data used in the identification are limited, we can provide a rather accurate solution of α(x).     

Classification of type I and type II behaviors for a supercritical nonlinear heat equation

We study blow-up of radially symmetric solutions of the nonlinear heat equation ut=Δu+|u|^p−1u either on RN or on a finite ball under the Dirichlet boundary conditions. We assume and that the initial data is bounded, possibly sign-changing. Our first goal is to establish various characterizations of type I and type II blow-ups. Among many other things we show that the following conditions are equivalent: (a) the blow-up is of type II; (b) the rescaled solution w(y,s) converges to either φ^∗(y) or −φ^∗(y) as s→∞, where φ^∗ denotes the singular stationary solution; (c) u(x,T)/φ^∗(x) tends to ±1 as x→0, where T is the blow-up time.Our second goal is to study continuation beyond blow-up. Among other things we show that if a blow-up is of type I and incomplete, then its limit L¹ continuation becomes smooth immediately after blow-up, and that type I blow-up implies “type I regularization,” that is, (t−T)^1/(p−1)‖u(⋅,t)L^∞_‖ is bounded as t↘T. We also give various criteria for complete and incomplete blow-ups.     

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 variational approach to dislocation problems for periodic Schrödinger operators

As a simple model for lattice defects like grain boundaries in solid state physics we consider potentials which are obtained from a periodic potential V=V(x,y) on R² with period lattice Z² by setting Wt(x,y)=V(x+t,y) for x<0 and Wt(x,y)=V(x,y) for x?0, for t∈[0,1]. For Lipschitz-continuous V it is shown that the Schrödinger operators Ht=−Δ+Wt have spectrum (surface states) in the spectral gaps of H₀, for suitable t∈(0,1). We also discuss the density of these surface states as compared to the density of the bulk. Our approach is variational and it is first applied to the well-known dislocation problem (Korotyaev (2000, 2005) [15] and [16]) on the real line. We then proceed to the dislocation problem for an infinite strip and for the plane. In Appendix A, we discuss regularity properties of the eigenvalue branches in the one-dimensional dislocation problem for suitable classes of potentials.     

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.     

Stability interval for explicit difference schemes for multi-dimensional second-order hyperbolic equations with significant first-order space derivative terms

In this piece of work, we introduce a new idea and obtain stability interval for explicit difference schemes of O(k²+h²) for one, two and three space dimensional second-order hyperbolic equations u_tt=a(x,t)u_xx+α(x,t)u_x-2η²(x,t)u,u_tt=a(x,y,t)u_xx+b(x,y,t)u_yy+α(x,y,t)u_x+β(x,y,t)u_y-2η²(x,y,t)u, and u_tt=a(x,y,z,t)u_xx+b(x,y,z,t)u_yy+c(x,y,z,t)u_zz+α(x,y,z,t)u_x+β(x,y,z,t)u_y+γ(x,y,z,t)u_z-2η²(x,y,z,t)u,0<x,y,z<1,t>0 subject to appropriate initial and Dirichlet boundary conditions, where h>0 and k>0 are grid sizes in space and time coordinates, respectively. A new idea is also introduced to obtain explicit difference schemes of O(k²) in order to obtain numerical solution of u at first time step in a different manner.     

Automatic Controls for Nonlinear Dynamical Systems with Lipschitzian Trajectories

We would like to investigate on the solution to the automatic control problem given by the differential equation y′(t) = f(t, y(t), w(t)) for a given initial function x in the initial domain D(x, ω, Y) for almost all t in the interval I, with controls given by w(t) = g(t, y(t), T(y)(t)), where T is a nonanticipating and Lipschitzian operator. The result will be generalized for a dynamical system y′(t) = f(t, y(t), T(y), u(t)).     

Eigenvalue problems for second-order nonlinear dynamic equations on time scales

This paper is concerned with the existence and nonexistence of positive solutions of the second-order nonlinear dynamic equation uΔΔ(t)+λa(t)f(u(σ(t)))=0, t∈[0,1], satisfying either the conjugate boundary conditions u(0)=u(σ(1))=0 or the right focal boundary conditions u(0)=uΔ(σ(1))=0, where a and f are positive. We show that there exists a λ^∗>0 such that the above boundary value problem has at least two, one and no positive solutions for 0<λ<λ^∗, λ=λ^∗ and λ>λ^∗, respectively. Furthermore, by using the semiorder method on cones of the Banach space, we establish an existence and uniqueness criterion for positive solution of the problem. In particular, such a positive solution uλ(t) of the problem depends continuously on the parameter λ, i.e., uλ(t) is nondecreasing in λ, limλ→⁰⁺‖uλ‖=0 and limλ→+∞‖uλ‖=+∞.     

On solvability of two-dimensional Lp-Minkowski problem

Given p≠0 and a positive continuous function g, with g(x+T)=g(x), for some 0<T<1 and all real x, it is shown that for suitable choice of a constant C>0 the functional has a minimizer in the class of positive functions u∈C¹(R) for which u(x+T)=u(x) for all x∈R. This minimizer is used to prove the existence of a positive periodic solution y∈C²(R) of two-dimensional L_p-Minkowski problem y^1−p(x)(y″(x)+y(x))=g(x), where p∉{0,2}.     

Study of weak solutions for parabolic equations with nonstandard growth conditions

The authors of this paper study the existence and uniqueness of weak solutions of the initial and boundary value problem for ut=div((uσ+d₀)|∇u|^p(x,t)−2∇u)+f(x,t). Localization property of weak solutions is also discussed.     

Duality for a Class of Multiobjective Control Problems

In this paper, a class of multiobjective control problems is considered, where the objective and constraint functions involved are f(t, x(t), ?(t), y(t), z(t)) with x(t) ∈ Rⁿ, y(t) ∈ Rⁿ, and z(t) ∈ R^m, where x(t) and z(t) are the control variables and y(t) is the state variable. Under the assumption of invexity and its generalization, duality theorems are proved through a parametric approach to related properly efficient solutions of the primal and dual problems.     

Convergence to strong nonlinear diffusion waves for solutions to p-system with damping on quadrant

In this paper, we consider the so-called p-system with linear damping on quadrant. We show that for a certain class of given large initial data (v₀(x),u₀(x)), the corresponding initial-boundary value problem admits a unique global smooth solution (v(x,t),u(x,t)) and such a solution tends time-asymptotically, at the Lp (2?p?∞) optimal decay rates, to the corresponding nonlinear diffusion wave which satisfies (1.9) provided the corresponding prescribed initial error function (V₀(x),U₀(x)) lies in (H³(R⁺)∩L¹(R⁺))×(H²(R⁺)∩L¹(R⁺)).     

Global solutions of the fast diffusion equation with gradient absorption terms

We investigate the existence of nonnegative weak solutions to the problem ut=Δ(um)−p|∇u| in Rn×(0,∞) with ⁺(1−2/n)<m<1. It will be proved that: (i) When 1<p<2, if the initial datum u₀∈D(Rn) then there exists a solution; (ii) When 1<p<(2+mn)/(n+1), if the initial datum u₀(x) is a bounded and nonnegative measure then the solution exists; (iii) When (2+mn)/(n+1)?p<2, if the initial datum is a Dirac mass then the solution does not exist. We also study the large time behavior of the L¹-norm of solutions for 1<p?(2+mn)/(n+1), and the large time behavior of t1/β‖u(⋅,t)−Ec(⋅,t)L^∞_‖ for (2+mn)/(n+1)<p<2.     

Critical Exponents of Quasilinear Parabolic Equations

In this paper we study the critical exponents of the Cauchy problem in Rⁿ of the quasilinear singular parabolic equations: u_t = div(|∇u|^m − 1∇u) + t^s|x|^σu^p, with non-negative initial data. Here s ≥ 0, (n − 1)/(n + 1) < m < 1, p > 1 and σ > n(1 − m) − (1 + m + 2s). We prove that p_c ≡ m + (1 + m + 2s + σ)/n > 1 is the critical exponent. That is, if 1 < p ≤ p_c then every non-trivial solution blows up in finite time, but for p > p_c, a small positive global solution exists.     

Fractal oscillations of self-adjoint and damped linear differential equations of second-order

For a prescribed real number s ∈ [1, 2), we give some sufficient conditions on the coefficients p(x) and q(x) such that every solution y = y(x), y ∈ C²((0, T]) of the linear differential equation (p(x)y′)′ + q(x)y = 0 on (0, T], is bounded and fractal oscillatory near x = 0 with the fractal dimension equal to s. This means that y oscillates near x = 0 and the fractal (box-counting) dimension of the graph Γ(y) of y is equal to s as well as the s dimensional upper Minkowski content (generalized length) of Γ(y) is finite and strictly positive. It verifies that y admits similar kind of the fractal geometric asymptotic behaviour near x = 0 like the chirp function y_ch(x) = a(x)S(φ(x)), which often occurs in the time-frequency analysis and its various applications. Furthermore, this kind of oscillations is established for the Bessel, chirp and other types of damped linear differential equations given in the form y″ + (μ/x)y′ + g(x)y = 0, x ∈ (0, T]. In order to prove the main results, we state a new criterion for fractal oscillations near x = 0 of real continuous functions which essentially improves related one presented in [1].