Semilinear wave equation with time dependent potential

We consider the following semilinear wave equation: (1) for (t,x) ∈ ?_t × ?. We prove that if the potential V(t,x) is a measurable function that satisfies the following decay assumption: ∣V(t,x)∣?C(1+t)(1+∣x∣) for a.e. (t,x) ∈ ?_t × ? where C, σ₀>0 are real constants, then for any real number λ that satisfies there exists a real number ρ(f,g,λ)>0 such that the equation has a global solution provided that 0<ρ?ρ(f,g,λ). Copyright © 2004 John Wiley & Sons, Ltd.     

The Cauchy problem for quasilinear SG‐hyperbolic systems

We study the Cauchy problem for a class of quasilinear hyperbolic systems with coefficients depending on (t, x) ∈ [0, T ] × ?ⁿ and presenting a linear growth for |x | → ∞. We prove well‐posedness in the Schwartz space ?? (?ⁿ ). The result is obtained by deriving an energy estimate for the solution of the linearized problem in some weighted Sobolev spaces and applying a fixed point argument. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Gerdjikov‐Ivanov–type derivative nonlinear Schrödinger equation: Long‐time dynamics of nonzero boundary conditions

We consider the Gerdjikov‐Ivanov–type derivative nonlinear Schrödinger equation on the line. The initial value q(x,0) is given and satisfies the symmetric, nonzero boundary conditions at infinity, that is, q(x,0)→q_± as x→±∞, and |q_±|=q₀>0. The goal of this paper is to study the asymptotic behavior of the solution of this initial value problem as t→∞. The main tool is the asymptotic analysis of an associated matrix Riemann‐Hilbert problem by using the steepest descent method and the so‐called g‐function mechanism. We show that the solution q(x,t) of this initial value problem has a different asymptotic behavior in different regions of the xt‐plane. In the regions and , the solution takes the form of a plane wave. In the region , the solution takes the form of a modulated elliptic wave.     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

A best approximation for the solution of one‐dimensional variable‐coefficient Burgers' equation

In this article, an iterative method for the approximate solution to one‐dimensional variable‐coefficient Burgers' equation is proposed in the reproducing kernel space W_(3,2). It is proved that the approximation w_n(x,t) converges to the exact solution u(x,t) for any initial function w₀(x,t) ε W_(3,2), and the approximate solution is the best approximation under a complete normal orthogonal system . Moreover the derivatives of w_n(x,t) are also uniformly convergent to the derivatives of u(x,t).© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Determination of a control function in three‐dimensional parabolic equations by Legendre pseudospectral method

A Legendre pseudospectral method is proposed for solving approximately an inverse problem of determining an unknown control parameter p(t) which is the coefficient of the solution u(x, y, z, t) in a diffusion equation in a three‐dimensional region. The diffusion equation is to be solved subject to suitably prescribed initial‐boundary conditions. The presence of the unknown coefficient p(t) requires an extra condition. This extra condition considered as the integral overspecification over the spacial domain. For discretizing the problem, after homogenization of the boundary conditions, we apply the Legendre pseudospectral method in a matrix based manner. As a results a system of nonlinear differential algebraic equations is generated. Then by using suitable transformation, the problem will be converted to a homogeneous time varying system of linear ordinary differential equations. Also a pseudospectral method for efficient solving of the resulted system of ordinary differential equations is proposed. The solution of this system gives the approximation to values of u and p. The matrix based structure of the present method makes it easy to implement. Numerical experiments are presented to demonstrate the accuracy and the efficiency of the proposed computational procedure. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 74‐93, 2012     

On computing minimal realizable spectral radii of non‐negative matrices

For decades considerable efforts have been exerted to resolve the inverse eigenvalue problem for non‐negative matrices. Yet fundamental issues such as the theory of existence and the practice of computation remain open. Recently, it has been proved that, given an arbitrary (n–1)‐tuple ?? = (λ₂,…,λ_n) ∈ ?^n–1 whose components are closed under complex conjugation, there exists a unique positive real number ?(??), called the minimal realizable spectral radius of ??, such that the set {λ₁,…,λ_n} is precisely the spectrum of a certain n × n non‐negative matrix with λ₁ as its spectral radius if and only if λ₁ ? ?(??). Employing any existing necessary conditions as a mode of checking criteria, this paper proposes a simple bisection procedure to approximate the location of ?(??). As an immediate application, it offers a quick numerical way to check whether a given n‐tuple could be the spectrum of a certain non‐negative matrix. Copyright © 2004 John Wiley & Sons, Ltd.     

Stability and convergence of the spectral Galerkin method for the Cahn‐Hilliard equation

A spectral Galerkin method in the spatial discretization is analyzed to solve the Cahn‐Hilliard equation. Existence, uniqueness, and stabilities for both the exact solution and the approximate solution are given. Using the theory and technique of a priori estimate for the partial differential equation, we obtained the convergence of the spectral Galerkin method and the error estimate between the approximate solution u_N(t) and the exact solution u(t). © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

A Riemann-Hilbert Approach to the Chen-Lee-Liu Equation on the Half Line

In this paper, the Fokas unified method is used to analyze the initial-boundary value for the Chen- Lee-Liu equation

$i{\partial _t}u + {\partial_{xx}u - i |u{|^2}{\partial _x}u = 0}$

on the half line (?∞, 0] with decaying initial value. Assuming that the solution u(x, t) exists, we show that it can be represented in terms of the solution of a matrix Riemann-Hilbert problem formulated in the plane of the complex spectral parameter λ. The jump matrix has explicit (x, t) dependence and is given in terms of the spectral functions {a(λ), b(λ)} and {A(λ), B(λ)}, which are obtained from the initial data u₀(x) = u(x, 0) and the boundary data g₀(t) = u(0, t), g₁(t) = u_x(0, t), respectively. The spectral functions are not independent, but satisfy a so-called global relation.

     

Asymptotic analysis and estimates of blow‐up time for the radial symmetric semilinear heat equation in the open‐spectrum case

We estimate the blow‐up time for the reaction diffusion equation u_t=Δu+ λf(u), for the radial symmetric case, where f is a positive, increasing and convex function growing fast enough at infinity. Here λ>λ^*, where λ^* is the ‘extremal’ (critical) value for λ, such that there exists an ‘extremal’ weak but not a classical steady‐state solution at λ=λ^* with ∥w(?, λ)∥_∞→∞ as 0<λ→λ^*?. Estimates of the blow‐up time are obtained by using comparison methods. Also an asymptotic analysis is applied when f(s)=e^s, for λ?λ^*?1, regarding the form of the solution during blow‐up and an asymptotic estimate of blow‐up time is obtained. Finally, some numerical results are also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Computable analysis of a boundary‐value problem for the generalized KdV–Burgers equation

In this paper, we investigate the computability of the solution operator of the generalized KdV‐Burgers equation with initial‐boundary value problem. Here, the solution operator is a nonlinear map H^{3m ? 1}(R⁺) × H^m(0,T)→C([0,T];H^{3m ? 1}(R⁺)) from the initial‐boundary value data to the solution of the equation. By a technique that is widely used for the study of nonlinear dispersive equation, and using the type 2 theory of effectivity as computable model, we prove that the solution map is Turing computable, for any integer m ≥ 2, and computable real number T > 0. Copyright © 2014 John Wiley & Sons, Ltd.     

On the Complexity of Bifurcation in a Family of Equations Involving Lipschitz Operators

We investigate branches of eigenfunctions starting from a bifurcation point of Equation (1.2), x — F(x, λ) = 0, in a Banach space X for real λ. Elements of the total derivative of F(·, λ) in zero may create different directions of bifurcating branches depending in a complicated way on neighbouring λ. In special situatons we are able to calculate curves of eigenfunctions leading to segments consisting entirely of bifurcation points λ_t. In this investigation the finite dimensional Ljapunov‐Schmidt branching equations are used, simplified by trunkation.     

Continuity of the scattering function and Levinson‐type formula of Klein–Gordon equation

We considered the inverse problem of scattering theory for a boundary value problem on the half line generated by Klein–Gordon differential equation with a nonlinear spectral parameter‐dependent boundary condition. We defined the scattering data, and we proved the continuity of the scattering function S(λ); in a special case, the relation for the difference of the logarithm of the scattering function, which is called the Levinson‐type formula, was obtained. Copyright © 2012 John Wiley & Sons, Ltd.     

A rank‐exploiting infinite Arnoldi algorithm for nonlinear eigenvalue problems

We consider the nonlinear eigenvalue problem M(λ)x = 0, where M(λ) is a large parameter‐dependent matrix. In several applications, M(λ) has a structure where the higher‐order terms of its Taylor expansion have a particular low‐rank structure. We propose a new Arnoldi‐based algorithm that can exploit this structure. More precisely, the proposed algorithm is equivalent to Arnoldi's method applied to an operator whose reciprocal eigenvalues are solutions to the nonlinear eigenvalue problem. The iterates in the algorithm are functions represented in a particular structured vector‐valued polynomial basis similar to the construction in the infinite Arnoldi method [Jarlebring, Michiels, and Meerbergen, Numer. Math., 122 (2012), pp. 169–195]. In this paper, the low‐rank structure is exploited by applying an additional operator and by using a more compact representation of the functions. This reduces the computational cost associated with orthogonalization, as well as the required memory resources. The structure exploitation also provides a natural way in carrying out implicit restarting and locking without the need to impose structure in every restart. The efficiency and properties of the algorithm are illustrated with two large‐scale problems. Copyright © 2016 John Wiley & Sons, Ltd.     

Recovering a time‐dependent coefficient in a parabolic equation from overspecified boundary data using the pseudospectral Legendre method

The aim of this article is to discuss the problem of finding the unknown function u(x,t) and the time‐dependent coefficient a(t) in a parabolic partial differential equation. The pseudospectral Legendre method is employed to solve this problem. The results of numerical experiments are given. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

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.     

On pointwise convergence of the family of Urysohn‐type integral operators

In this paper, we consider a family of nonlinear integral operators of Urysohn‐type and study the pointwise convergence of the family at characteristic points of L₁?function. The kernel K_λ(x,t,u(t)) depends on the positive parameter λ changing on the set of numbers with the accumulation point at infinity and K_λ(x,t,u(t)) is an entire analytic function of variable u, which is a bounded function belonging to L₁(R).     

Fundamental solution for linear two-point boundary value problem

The quadratic functional minimization with differential restrictions represented by the command linear systems is considered. The optimal solution determination implies the solving of a linear problem with two points boundary values. The proposed method implies the construction of a fundamental solution S(t)—a n×n matrix—and of a vector h(t) defining an adjoint variable λ(t) depending of the state variable x(t). From the extremum necessary conditions it is obtained the Riccati matrix differential equation having the S(t) as unknown fundamental solution is obtained. The paper analyzes the existence of the Riccati equation solution S(t) and establishes as the optimal solution of the proposed optimum problem. Also a superior limit of the minimum for the considered quadratic functionals class are evaluated.     

An Application of a Mountain Pass Theorem

We are concerned with the following Dirichlet problem: −Δu(x) = f(x, u), x∈Ω, u∈H ¹ ₀(Ω), (P) where f(x, t) ∈C (×ℝ), f(x, t)/t is nondecreasing in t∈ℝ and tends to an L ^∞-function q(x) uniformly in x∈Ω as t→ + ∞ (i.e., f(x, t) is asymptotically linear in t at infinity). In this case, an Ambrosetti-Rabinowitz-type condition, that is, for some θ > 2, M > 0, 0 > θF(x, s) ≤f(x, s)s, for all |s|≥M and x∈Ω, (AR) is no longer true, where F(x, s) = ∫ ^s ₀ f(x, t)dt. As is well known, (AR) is an important technical condition in applying Mountain Pass Theorem. In this paper, without assuming (AR) we prove, by using a variant version of Mountain Pass Theorem, that problem (P) has a positive solution under suitable conditions on f(x, t) and q(x). Our methods also work for the case where f(x, t) is superlinear in t at infinity, i.e., q(x) ≡ +∞. Received June 24, 1998, Accepted January 14, 2000.     

Identification of a time‐dependent coefficient in a partial differential equation subject to an extra measurement

The problem of recovering a time‐dependent coefficient in a parabolic partial differential equation has attracted considerable recent attention. Several finite difference schemes are presented for identifying the function u(x, t) and the unknown coefficient a(t) in a one‐dimensional partial differential equation. These schemes are developed to determine the unknown properties in a region by measuring only data on the boundary. Our goal has been focused on coefficients that presents physical quantities, for example, the conductivity of a medium. For the convenience of discussion, we will present the results of numerical experiment on several test problems. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005