A hybrid collocation method for a nonlinear Volterra integral equation with weakly singular kernel

This work is concerned with the numerical solution of a nonlinear weakly singular Volterra integral equation. Owing to the singular behavior of the solution near the origin, the global convergence order of product integration and collocation methods is not optimal. In order to recover the optimal orders a hybrid collocation method is used which combines a non-polynomial approximation on the first subinterval followed by piecewise polynomial collocation on a graded mesh. Some numerical examples are presented which illustrate the theoretical results and the performance of the method. A comparison is made with the standard graded collocation method.     

On the stability of the one-step exact collocation methods for the numerical solution of the second kind Volterra integral equation

The purpose of this paper is to analyze the stability properties of one-step collocation methods for the second kind Volterra integral equation through application to the basic test and the convolution test equation.Stability regions are determined when the collocation parameters are symmetric and when they are zeros of ultraspherical polynomials.     

Biorthogonal systems for solving Volterra integral equation systems of the second kind

With the aid of biorthogonal systems in adequate Banach spaces, the problem of approximating the solution of a system of nonlinear Volterra integral equations of the second kind is turned into a numerical method that allows it to be solved numerically.     

The Galerkin scheme for Lavrentiev’sm-times iterated method to solve linear accretive Volterra integral equations of the first kind

In this paper, for the numerical solution of linear accretive Volterra integral equations of the first kind in Hilbert spaces we consider the Galerkin scheme for Lavrentiev’sm-times iterated method, i.e., for each parameter choice for Lavrentiev’sm-times iterated method the arisingm stabilized equations are discretized by the Galerkin scheme. An associated discrepancy principle as parameter choice strategy for this finite-dimensional version of Lavrentiev’sm-times iterated method is proposed, and corresponding convergence results are provided.     

Legendre spectral collocation method for volterra-hammerstein integral equation of the second kind

This paper is concerned with obtaining the approximate solution for VolterraHammerstein integral equation with a regular kernel. We choose the Gauss points associated with the Legendre weight function ω(x) = 1 as the collocation points. The Legendre collocation discretization is proposed for Volterra-Hammerstein integral equation. We provide an error analysis which justifies that the errors of approximate solution decay exponentially in L~2 norm and L~∞ norm. We give two numerical examples in order to illustrate the validity of the proposed Legendre spectral collocation method.     

Dependence on the parameters of the set of trajectories of the control system described by a nonlinear Volterra integral equation

In this paper the control system with limited control resources is studied, where the behavior of the system is described by a nonlinear Volterra integral equation. The admissible control functions are chosen from the closed ball centered at the origin with radius µ in L _p (p > 1). It is proved that the set of trajectories generated by all admissible control functions is Lipschitz continuous with respect to µ for each fixed p, and is continuous with respect to p for each fixed µ. An upper estimate for the diameter of the set of trajectories is given.     

Local controllability to stationary trajectories of a Burgers equation with nonlocal viscosity

This article studies the local controllability to trajectories of a Burgers equation with nonlocal viscosity. By linearization we are led to an equation with a non local term whose controllability properties are analyzed by using Fourier decomposition and biorthogonal techniques. Once the existence of controls is proved and the dependence of their norms with respect to the time is established for the linearized model, a fixed point method allows us to deduce the result for the nonlinear initial problem.     

On the nonnegative self-adjoint solutions of the operator Riccati equation for infinite dimensional systems

The nonnegative self-adjoint solutions of the operator Riccati equation (ORE) are studied for stabilizable semigroup Hilbert state space systems with bounded sensing and control. Basic properties of the maximal solution of the ORE are investigated: stability of the corresponding closed loop system, structure of the kernel, Hilbert-Schmidt property. Similar properties are obtained for the nonnegative self-adjoint solutions of the ORE. The analysis leads to a complete classification of all nonnegative self-adjoint solutions, which is based on a bijection between these solutions and finite dimensional semigroup invariant subspaces contained in the antistable unobservable subspace.     

Approximate the Fokker-Planck equation by a class of nonlocal dispersal problems

This paper is concerned with an inhomogeneous nonlocal dispersal equation. We study the limit of the re-scaled problem of this nonlocal operator and prove that the solutions of the re-scaled equation converge to a solution of the Fokker-Planck equation uniformly. We then analyze the nonlocal dispersal equation of an inhomogeneous diffusion kernel and find that the heterogeneity in the classical diffusion term coincides with the inhomogeneous kernel when the scaling parameter goes to zero.     

Approximate solution of the system of nonlinear integral equation by Newton-Kantorovich method

The Newton-Kantorovich method is developed for solving the system of nonlinear integral equations. The existence and uniqueness of the solution are proved, and the rate of convergence of the approximate solution is established. Finally, numerical examples are provided to show the validity and the efficiency of the method presented.     

Numerical exploitation of symmetry in integral equations

Linear integral operators describing physical problems on symmetric domains often are equivariant, which means that they commute with certain symmetries, i.e., with a group of orthogonal transformations leaving the domain invariant. Under suitable discretizations the resulting system matrices are also equivariant. A method for exploiting this equivariance in the numerical solution of linear equations and eigenvalue problems via symmetry reduction is described. A very significant reduction in the computational expense in both the assembling of the system matrix and in solving linear systems can be obtained in this way. This reduction is particularly important because the system matrices are typically full. The basic ideas underlying our method and its analysis involve group representation theory. We concentrate here on the use of symmetry adapted bases and their automated generation. In this paper symmetry reduction is studied in connection with quadrature formulae and the Nyström method. A software package has been posted on the Internet.     

Second kind integral equation formulation for the mode calculation of optical waveguides

We present a second kind integral equation (SKIE) formulation for calculating the electromagnetic modes of optical waveguides, where the unknowns are only on material interfaces. The resulting numerical algorithm can handle optical waveguides with a large number of inclusions of arbitrary irregular cross section. It is capable of finding the bound, leaky, and complex modes for optical fibers and waveguides including photonic crystal fibers (PCF), dielectric fibers and waveguides. Most importantly, the formulation is well conditioned even in the case of nonsmooth geometries. Our method is highly accurate and thus can be used to calculate the propagation loss of the electromagnetic modes accurately, which provides the photonics industry a reliable tool for the design of more compact and efficient photonic devices. We illustrate and validate the performance of our method through extensive numerical studies and by comparison with semi-analytical results and previously published results.     

Solvability of an implicit fractional integral equation via a measure of noncompactness argument

In this paper, we study the existence of solutions to an implicit functional equation involving a fractional integral with respect to a certain function, which generalizes the Riemann-Liouville fractional integral and the Hadamard fractional integral. We establish an existence result to such kind of equations using a generalized version of Darbo's theorem associated to a certain measure of noncompactness. Some examples are presented.     

On the solution of a mixed nonlinear integral equation

In this paper, we consider a mixed nonlinear integral equation of the second kind in position and time. The existence of a unique solution of this equation is discussed and proved. A numerical method is used to obtain a system of Harmmerstein integral equations of the second kind in position. Then the modified Toeplitz matrix method, as a numerical method, is used to obtain a nonlinear algebraic system. Many important theorems related to the existence and uniqueness solution to the produced nonlinear algebraic system are derived. The rate of convergence of the total error is discussed. Finally, numerical examples when the kernel of position takes a logarithmic and Carleman forms, are presented and the error estimate, in each case, is calculated.     

An efficient parallel processing optimal control scheme for a class of nonlinear composite systems

This article presents an efficient parallel processing approach for solving the optimal control problem of nonlinear composite systems. In this approach, the original high-order coupled nonlinear two-point boundary value problem (TPBVP) derived from the Pontryagin's maximum principle is first transformed into a sequence of lower-order decoupled linear time-invariant TPBVPs. Then, an optimal control law which consists of both feedback and forward terms is achieved by using the modal series method for the derived sequence. The feedback term specified by local states of each subsystem is determined by solving a matrix Riccati differential equation. The forward term for each subsystem derived from its local information is an infinite sum of adjoint vectors. The convergence analysis and parallel processing capability of the proposed approach are also provided. To achieve an accurate feedforward-feedback suboptimal control, we apply a fast iterative algorithm with low computational effort. Finally, some comparative results are included to illustrate the effectiveness of the proposed approach.     

A nonstandard multigrid method with flexible multiple semicoarsening for the numerical solution of the pressure equation in a Navier-Stokes solver

Numerical methods for the incompressible Reynolds-averaged Navier-Stokes equations discretized by finite difference techniques on collocated cell-centered structured grids are considered in this paper. A widespread solution method to solve the pressure-velocity coupling problem is to use a segregated approach, in which the computational work is deeply controlled by the solution of the pressure problem. This pressure equation is an elliptic partial differential equation with possibly discontinuous or anisotropic coeffficients. The resulting singular linear system needs efficient solution strategies especially for 3-dimensional applications. A robust method (close to MG-S [22,34]) combining multiple cell-centered semicoarsening strategies, matrix-independent transfer operators, Galerkin coarse grid approximation is therefore designed. This strategy is both evaluated as a solver or as a preconditioner for Krylov subspace methods on various 2- or 3-dimensional fluid flow problems. The robustness of this method is shown. This revised version was published online in June 2006 with corrections to the Cover Date.     

Regularity results and optimal control problems for the perturbation of boussinesq equations of the ocean

We study in this article a method which computes the variability of current, density and pressure in an oceanic domain. The equations are of Navier-Stokes type for the velocity and pressure, of transport-diffusion type for the density. They are linearized around a given mean circulation and modified by the Boussinesq approximation: density variations are neglected except in the terms of gravity acceleration. The existence and uniqueness of a solution are proved for two sets of equations: first the three-dimensional problem and then the two-dimensional cyclic problem derived by assuming a sinusoidal x-dependence for the perturbation of mean flow. The latter corresponds to a modelization of tropical instability waves which are illustrated by the El Nino phenomenon.

The value of the pressure p on the surface of ocean is of great interest for physical interpretation. To define that quantity, it is necessary to have the regularity p ? H ¹. We have proved that the perturbation (u,ρ,p) of mean circulation is such that: u ? L ²(0T,H ²), ρ ? L ²(0,T H ²) and p ? L ² L ²(0,T H ¹), provided the perturbation of the windstress is sufficiently regular and satisfies compatibility relations. It is proved by means of an extension method, with even-odd reflection. We then develop a problem of control. The observation is the Variability of pressure on the surface of ocean. The control is the variability of windstress f, which acts as to forcing of the perturbation. We prove the existence and uniqueness of an optimal control, which is characterized by a set of equations including the direct problem and the adjoint problem. These results are valid for the three-dimensional problem and the two-dimensional cyclic problem.     

The projected newton method for solving inverse eigenvalue problems - the case of multiple eigenvaues

This paper deals with the optimal solution of ill-posed linear problems, i.e..linear problems for which the solution operator is unbounded. We consider worst-case ar,and averagecase settings. Our main result is that algorithms having finite error (for a given setting) exist if and only if the solution operator is bounded (in that setting). In the worst-case setting, this means that there is no algorithm for solving ill-posed problems having finite error. In the average-case setting, this means that algorithms having finite error exist if and only lf the solution operator is bounded on the average. If the solution operator is bounded on the average, we find average-case optimal information of cardinality n and optimal algorithms using this information, and show that the average error of these algorithms tends to zero as n→∞. These results are then used to determine the [euro]-complexity, i.e., the minimal costof finding an [euro]-accurate approximation. In the worst-case setting, the [euro]comp1exity of an illposed problem is infinite for all [euro]>0; that is, we cannot find an approximation having finite error and finite cost. In the average-case setting, the [euro]-complexity of an ill-posed problem is infinite for all [euro]>0 iff the solution operator is not bounded on the average, moreover, if the the solutionoperator is bounded on the average, then the [euro]-complexity is finite for all [euro]>0.     

Direct pseudo‐spectral method for optimal control of obstacle problem – an optimal control problem governed by elliptic variational inequality

In this paper, a computational technique based on the pseudo‐spectral method is presented for the solution of the optimal control problem constrained with elliptic variational inequality. In fact, our aim in this paper is to present a direct approach for this class of optimal control problems. By using the pseudo‐spectral method, the infinite dimensional mathematical programming with equilibrium constraint, which can be an equivalent form of the considered problem, is converted to a finite dimensional mathematical programming with complementarity constraint. Then, the finite dimensional problem can be solved by the well‐developed methods. Finally, numerical examples are presented to show the validity and efficiency of the technique. Copyright © 2017 John Wiley & Sons, Ltd.