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     

The Generalized Dirichlet to Neumann Map for the KdV Equation on the Half-Line

For the two versions of the KdV equation on the positive half-line an initial-boundary value problem is well posed if one prescribes an initial condition plus either one boundary condition if q _t and q _xxx have the same sign (KdVI) or two boundary conditions if q _t and q _xxx have opposite sign (KdVII). Constructing the generalized Dirichlet to Neumann map for the above problems means characterizing the unknown boundary values in terms of the given initial and boundary conditions. For example, if {q(x,0),q(0,t)} and {q(x,0),q(0,t),q _x (0,t)} are given for the KdVI and KdVII equations, respectively, then one must construct the unknown boundary values {q _x (0,t),q _xx (0,t)} and {q _xx (0,t)}, respectively. We show that this can be achieved without solving for q(x,t) by analysing a certain “global relation” which couples the given initial and boundary conditions with the unknown boundary values, as well as with the function Φ ^(t)(t,k), where Φ ^(t) satisfies the t-part of the associated Lax pair evaluated at x=0. The analysis of the global relation requires the construction of the so-called Gelfand–Levitan–Marchenko triangular representation for Φ ^(t). In spite of the efforts of several investigators, this problem has remained open. In this paper, we construct the representation for Φ ^(t) for the first time and then, by employing this representation, we solve explicitly the global relation for the unknown boundary values in terms of the given initial and boundary conditions and the function Φ ^(t). This yields the unknown boundary values in terms of a nonlinear Volterra integral equation. We also discuss the implications of this result for the analysis of the long t-asymptotics, as well as for the numerical integration of the KdV equation.     

Boundary Control of the Wave Process with Elastic Clamping

The article investigates the boundary control problem for the wave process described by the equation u _tt(x,t) − u _xx(x,t) = 0 in the time interval 0 < t ≤ T with elastic clamping at the point x = l. For 0 < T ≤ 2l necessary and sufficient conditions are obtained ensuring the existence of a unique boundary control and its analytical form is determined. For 2l < T ≤ 3l we derive an explicit analytical expression for this boundary control that contains two arbitrary functions (of class ) defined on a segment of length T − 2l.     

A transform method for linear evolution PDEs on a finite interval

We study initial boundary value problems for linear scalar evolutionpartial differential equations, with spatial derivatives ofarbitrary order, posed on the domain {t > 0, 0 < x <L}. We show that the solution can be expressed as an integralin the complex k-plane. This integral is defined in terms ofan x-transform of the initial condition and a t-transform ofthe boundary conditions. The derivation of this integral representationrelies on the analysis of the global relation, which is an algebraicrelation defined in the complex k-plane coupling all boundaryvalues of the solution. For particular cases, such as the case of periodic boundaryconditions, or the case of boundary value problems for even-orderPDEs, it is possible to obtain directly from the global relationan alternative representation for the solution, in the formof an infinite series. We stress, however, that there existinitial boundary value problems for which the only representationis an integral which cannot be written as an infinite series.An example of such a problem is provided by the linearized versionof the KdV equation. Similarly, in general the solution of odd-orderlinear initial boundary value problems on a finite intervalcannot be expressed in terms of an infinite series.     

S‐asymptotically ω‐periodic solutions for semilinear Volterra equations

We study S‐asymptotically ω‐periodic mild solutions of the semilinear Volterra equation u′(t)=(a* Au)(t)+f(t, u(t)), considered in a Banach space X, where A is the generator of an (exponentially) stable resolvent family. In particular, we extend the recent results for semilinear fractional integro‐differential equations considered in (Appl. Math. Lett. 2009; 22:865–870) and for semilinear Cauchy problems of first order given in (J. Math. Anal. Appl. 2008; 343(2): 1119–1130). Applications to integral equations arising in viscoelasticity theory are shown. Copyright © 2010 John Wiley & Sons, Ltd.     

Current work and open problems in the numerical analysis of Volterra functional equations with vanishing delays

The aims of this paper are (i) to present a survey of recent advances in the analysis of superconvergence of collocation solutions for linear Volterra-type functional integral and integro-differential equations with delay functions θ(t) vanishing at the initial point of the interval of integration (with ia(t) = qt (0 < q < 1, t ⩾ 0) being an important special case), and (ii) to point, by means of a list of open problems, to areas in the numerical analysis of such Volterra functional equations where more research needs to be carried out.      

Asymptotic analysis of a volterra integral equation

A uniformly valid aymptotic solution is obtained for a class of perturbed Volterra integral equations, in which a naive expansion breaks down as t → ∞. The procedure used is an adaption of the formal methodology presented in [1] for the construction of a uniform asymptotic solution to Volterra equations which possess a boundary layer near t = 0.     

Operator time‐splitting techniques combined with quintic B‐spline collocation method for the generalized Rosenau–KdV equation

In this article, the generalized Rosenau–KdV equation is split into two subequations such that one is linear and the other is nonlinear. The resulting subequations with the prescribed initial and boundary conditions are numerically solved by the first order Lie–Trotter and the second‐order Strang time‐splitting techniques combined with the quintic B‐spline collocation by the help of the fourth order Runge–Kutta (RK‐4) method. To show the accuracy and reliability of the proposed techniques, two test problems having exact solutions are considered. The computed error norms L₂ and L_∞ with the conservative properties of the discrete mass Q(t) and energy E(t) are compared with those available in the literature. The convergence orders of both techniques have also been calculated. Moreover, the stability analyses of the numerical schemes are investigated.     

Solution of nonlinear weakly singular Volterra integral equations using the fractional‐order Legendre functions and pseudospectral method

In this article, our main goal is to render an idea to convert a nonlinear weakly singular Volterra integral equation to a non‐singular one by new fractional‐order Legendre functions. The fractional‐order Legendre functions are generated by change of variable on well‐known shifted Legendre polynomials. We consider a general form of singular Volterra integral equation of the second kind. Then the fractional Legendre–Gauss–Lobatto quadratures formula eliminates the singularity of the kernel of the integral equation. Finally, the Legendre pseudospectral method reduces the solution of this problem to the solution of a system of algebraic equations. This method also can be utilized on fractional differential equations as well. The comparison of results of the presented method and other numerical solutions shows the efficiency and accuracy of this method. Also, the obtained maximum error between the results and exact solutions shows that using the present method leads to accurate results and fast convergence for solving nonlinear weakly singular Volterra integral equations. Copyright © 2015 John Wiley & Sons, Ltd.     

A numerical study of parabolic problems with nonlinear boundary conditions

In this paper, we consider an initial‐boundary value problem for a parabolic equation with nonlinear boundary conditions. The solution to the problem can be expressed as a convolution integral of a Green's function and two unknown functions. We change the problem to a system of two nonlinear Volterra integral equations of convolution type. By using an explicit procedure on the basis of Sinc‐function properties, the resulting integral equations are replaced by a system of nonlinear algebraic equations, whose solution yields an accurate approximate solution to the parabolic problem. Some examples are considered to illustrate the ability of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Damping of Periodic Waves in Physically Significant Wave Systems

Damping of periodic waves in the classically important nonlinear wave systems—nonlinear Schrödinger, Korteweg–deVries (KdV), and modified KdV—is considered here. For small damping, asymptotic analysis is used to find an explicit equation that governs the temporal evolution of the solution. These results are then confirmed by direct numerical simulations. The undamped periodic solutions are given in terms of Jacobi elliptic functions. The damping structure is found as a function of the elliptic function modulus, m=m(t) . The damping rate of the maximum amplitude is ascertained and is found to vary smoothly from the linear solution when m= 0 to soliton waves when m= 1 .     

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.     

On the solutions of the linear integral equations of Volterra type

Some boundaries about the solution of the linear Volterra integral equations of the form f(t)=1?K*f were obtained as |f(t)|?1, |f(t)|?2 and |f(t)|?4 in (J. Math. Anal. Appl. 1978; 64 :381–397; Int. J. Math. Math. Sci. 1982; 5 (1):123–131). The boundary of the solution function of an equation in this type was found as |f(t)|?2ⁿ in (Integr. Equ. Oper. Theory 2002; 43 :466–479), where t∈[0, ∞) and n is a natural number such that n?2. In (Math. Comp. 2006; 75 :1175–1199), it is shown that the boundary of the solution function of an equation in the same form can also be derived as that of (Integr. Equ. Oper. Theory 2002; 43 :466–479) under different conditions than those of (Integr. Equ. Oper. Theory 2002; 43 :466–479). In the present paper, the sufficient conditions for the boundedness of functions f, f′, f′′, …, f⁽ⁿ⁺³⁾, (n∈?) defined on the infinite interval [0, ∞) are given by our method, where f is the solution of the equation f(t)=1?K*f. Copyright © 2007 John Wiley & Sons, Ltd.     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

Optimal boundary control of forced vibrations by the displacement at one end of the string with the other end fixed

For a string vibration process described by an inhomogeneous wave equation, we consider the problem of boundary control at one end of the string with the other end being fixed. For any time interval T > 2l, where l is the string length, we find a function u(0, t) = µ(t) bringing the vibration system from a given initial state into a given terminal state and minimizing the boundary energy integral.     

Resolvents and solutions of weakly singular linear Volterra integral equations

The structure of the resolvent R(t,s) for a weakly singular matrix function B(t,s) is determined, where B(t,s) is the kernel of the linear Volterra vector integral equation

(E )     

An Initial-Boundary Value Problem for the Korteweg–de Vries Equation with Dominant Surface Tension

We consider the initial-boundary value problem (IBVP) for the Korteweg–de Vries equation with zero boundary conditions at x=0 and arbitrary smooth decreasing initial data. We prove that the solution of this IBVP can be found by solving two linear inverse scattering problems (SPs) on two different spectral planes. The first SP is associated with the KdV equation. The second SP is self-conjugate and its scattering function is found in terms of entries of the scattering matrix s(k) for the first SP. Knowing the scattering function, we solve the second inverse SP for finding the potential self-conjugate matrix. Consequently, the unknown object entering coefficients in the system of evolution equations for s(k,t) is found. Then, the time-dependent scattering matrix s(k,t) is expressed in terms of s(k)=s(k,0) and of solutions of the self-conjugate SP. Knowing s(k,t), we find the solution of the IBVP in terms of the solution of the Gelfand–Levitan–Marchenko equation in the first inverse SP.     

Analysis and numerical solution of a Riemann-Liouville fractional derivative two-point boundary value problem

A two-point boundary value problem is considered on the interval [0, 1], where the leading term in the differential operator is a Riemann-Liouville fractional derivative of order 2 ? δ with 0 < δ < 1. It is shown that any solution of such a problem can be expressed in terms of solutions to two associated weakly singular Volterra integral equations of the second kind. As a consequence, existence and uniqueness of a solution to the boundary value problem are proved, the structure of this solution is elucidated, and sharp bounds on its derivatives (in terms of the parameter δ) are derived. These results show that in general the first-order derivative of the solution will blow up at x = 0, so accurate numerical solution of this class of problems is not straightforward. The reformulation of the boundary problem in terms of Volterra integral equations enables the construction of two distinct numerical methods for its solution, both based on piecewise-polynomial collocation. Convergence rates for these methods are proved and numerical results are presented to demonstrate their performance.     

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