Shifted Jacobi spectral collocation method with convergence analysis for solving integro-differential equations and system of integro-differential equations

This article addresses the solution of multi-dimensional integro-differential equations (IDEs) by means of the spectral collocation method and taking the advantage of the properties of shifted Jacobi polynomials. The applicability and accuracy of the present technique have been examined by the given numerical examples in this paper. By means of these numerical examples, we ensure that the present technique is simple and very accurate. Furthermore, an error analysis is performed to verify the correctness and feasibility of the proposed method when solving IDE

Post-Processing Techniques and Wavelet Applications for Hammerstein Integral Equations

This dissertation is focused on the varieties of numerical solutions of nonlinear Hammerstein integral equations. In the first part of this dissertation, several acceleration techniques for post-processed solutions of the Hammerstein equation are discussed. The post-processing techniques are implemented based on interpolation and extrapolation. In this connection, we generalize the results in [29] and [28] to nonlinear integral equations of the Hammerstein type. Post-processed collocation solutions are shown to exhibit better accuracy. Moreover, an extrapolation technique for the Galerkin solution of Hammerstein equation is also obtained. This result appears new even in the setting of the linear Fredholm equation. In the second half of this dissertation, the wavelet-collocation technique of solving nonlinear Hammerstein integral equation is discussed. The main objective is to establish a fast wavelet-collocation method for Hammerstein equation by using a \u27linearization\u27 technique. The sparsity in the Jacobian matrix takes place in the fast wavelet-collocation method for Hammerstein equation with smooth as well as weakly singular kernels. A fast algorithm is based upon the block truncation strategy which was recently proposed in [10]. A multilevel augmentation method for the linearized Hammerstein equation is subsequently proposed which further accelerates the solution process while maintaining the order of convergence. Numerical examples are given throughout this dissertation

Numerical Approximate Methods for Solving Linear and Nonlinear Integral Equations

Integral equation has been one of the essential tools for various area of applied mathematics. In this work, we employed different numerical methods for solving both linear and nonlinear Fredholm integral equations. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study on numerical methods for solving integral equations. Integral equations can be viewed as equations which are results of transformation of points in a given vector spaces of integrable functions by the use of certain specific integral operators to points in the same space. If, in particular, one is concerned with function spaces spanned by polynomials for which the kernel of the corresponding transforming integral operator is separable being comprised of polynomial functions only, then several approximate methods of solution of integral equations can be developed. This work, specially, deals with the development of different wavelet methods for solving integral and intgro-differential equations. Wavelets theory is a relatively new and emerging area in mathematical research. It has been applied in a wide range of engineering disciplines; particularly, wavelets are very successfully used in signal analysis for waveform representations and segmentations, time frequency analysis, and fast algorithms for easy implementation. Wavelets permit the accurate representation of a variety of functions and operators. Moreover, wavelets establish a connection with fast numerical algorithms. Wavelets can be separated into two distinct types, orthogonal and semi-orthogonal. The preliminary concept of integral equations and wavelets are first presented in Chapter 1. Classification of integral equations, construction of wavelets and multi-resolution analysis (MRA) have been briefly discussed and provided in this chapter. In Chapter 2, different wavelet methods are constructed and function approximation by these methods with convergence analysis have been presented. In Chapter 3, linear semi-orthogonal compactly supported B-spline wavelets together with their dual wavelets have been applied to approximate the solutions of Fredholm integral equations (both linear and nonlinear) of the second kind and their systems. Properties of these wavelets are first presented; these properties are then utilized to reduce the computation of integral equations to some algebraic equations. Convergence analysis of B-spline method has been discussed in this chapter. Again, in Chapter 4, system of nonlinear Fredholm integral equations have been solved by using hybrid Legendre Block-Pulse functions and xiii Bernstein collocation method. In Chapter 5, two practical problems arising from chemical phenomenon, have been modeled as Fredholm- Hammerstein integral equations and solved numerically by different numerical techniques. First, COSMO-RS model has been solved by Bernstein collocation method, Haar wavelet method and Sinc collocation method. Second, Hammerstein integral equation arising from chemical reactor theory has been solved by B-spline wavelet method. Comparison of results have been demonstrated through illustrative examples. In Chapter 6, Legendre wavelet method and Bernoulli wavelet method have been developed to solve system of integro-differential equations. Legendre wavelets along with their operational matrices are developed to approximate the solutions of system of nonlinear Volterra integro-differential equations. Also, nonlinear Volterra weakly singular integro-differential equations system has been solved by Bernoulli wavelet method. The properties of these wavelets are used to reduce the system of integral equations to a system of algebraic equations which can be solved numerically by Newton's method. Rigorous convergence analysis has been done for these wavelet methods. Illustrative examples have been included to demonstrate the validity and applicability of the proposed techniques. In Chapter 7, we have solved the second order Lane-Emden type singular differential equation. First, the second order differential equation is transformed into integro-differential equation and then solved by Legendre multi-wavelet method and Chebyshev wavelet method. Convergence of these wavelet methods have been discussed in this chapter. In Chapter 8, we have developed a efficient collocation technique called Legendre spectral collocation method to solve the Fredholm integro-differential-difference equations with variable coefficients and system of two nonlinear integro-differential equations which arise in biological model. The proposed method is based on the Gauss-Legendre points with the basis functions of Lagrange polynomials. The present method reduces this model to a system of nonlinear algebraic equations and again this algebraic system has been solved numerically by Newton's method. The study of fuzzy integral equations and fuzzy differential equations is an emerging area of research for many authors. In Chapter 9, we have proposed some numerical techniques for solving fuzzy integral equations and fuzzy integro-differential equations. Fundamentals of fuzzy calculus have been discussed in this chapter. Nonlinear fuzzy Hammerstein integral equation has been solved by Bernstein polynomials and Legendre wavelets, and then compared with homotopy analysis method. We have solved nonlinear fuzzy Hammerstein Volterra integral equations with constant delay by Bernoulli wavelet method and then compared with B-spline wavelet method. Finally, fuzzy integro-differential equation has been solved by Legendre wavelet method and compared with homotopy analysis method. In fuzzy case, we have applied two-dimensional numerical methods which are discussed in chapter 2. Convergence analysis and error estimate have been also provided for Bernoulli wavelet method. xiv The study of fractional calculus, fractional differential equations and fractional integral equations has a great importance in the field of science and engineering. Most of the physical phenomenon can be best modeled by using fractional calculus. Applications of fractional differential equations and fractional integral equations create a wide area of research for many researchers. This motivates to work on fractional integral equations, which results in the form of Chapter 10. First, the preliminary definitions and theorems of fractional calculus have been presented in this chapter. The nonlinear fractional mixed Volterra-Fredholm integro-differential equations along with mixed boundary conditions have been solved by Legendre wavelet method. A numerical scheme has been developed by using Petrov-Galerkin method where the trial and test functions are Legendre wavelets basis functions. Also, this method has been applied to solve fractional Volterra integro-differential equations. Uniqueness and existence of the problem have been discussed and the error estimate of the proposed method has been presented in this work. Sinc Galerkin method is developed to approximate the solution of fractional Volterra-Fredholm integro-differential equations with weakly singular kernels. The proposed method is based on the Sinc function approximation. Uniqueness and existence of the problem have been discussed and the error analysis of the proposed method have been presented in this chapte

Generalized extrapolation methods for solving nonlinear Fredholm integral equations

In this paper we develop a class of generalized extrapolation methods for numerical solution of nonlinear Fredholm integral equations of the second kind. The direct representation of this class allows us to simply discretize the nonlinear integral equations with smooth kernels. This approach enjoys several outstanding features of numerical methods such as: economized computational cost, the high order accuracy, direct implementation, discrretization on arbitrary nodes and applying the methods with positive weights. The comparison results demonstrate the superior results of the new class of methods versus the classical and recent approaches

Applied Mathematics and Fractional Calculus

In the last three decades, fractional calculus has broken into the field of mathematical analysis, both at the theoretical level and at the level of its applications. In essence, the fractional calculus theory is a mathematical analysis tool applied to the study of integrals and derivatives of arbitrary order, which unifies and generalizes the classical notions of differentiation and integration. These fractional and derivative integrals, which until not many years ago had been used in purely mathematical contexts, have been revealed as instruments with great potential to model problems in various scientific fields, such as: fluid mechanics, viscoelasticity, physics, biology, chemistry, dynamical systems, signal processing or entropy theory. Since the differential and integral operators of fractional order are nonlinear operators, fractional calculus theory provides a tool for modeling physical processes, which in many cases is more useful than classical formulations. This is why the application of fractional calculus theory has become a focus of international academic research. This Special Issue "Applied Mathematics and Fractional Calculus" has published excellent research studies in the field of applied mathematics and fractional calculus, authored by many well-known mathematicians and scientists from diverse countries worldwide such as China, USA, Canada, Germany, Mexico, Spain, Poland, Portugal, Iran, Tunisia, South Africa, Albania, Thailand, Iraq, Egypt, Italy, India, Russia, Pakistan, Taiwan, Korea, Turkey, and Saudi Arabia

Identification of low order models for large scale processes

Many industrial chemical processes are complex, multi-phase and large scale in nature. These processes are characterized by various nonlinear physiochemical effects and fluid flows. Such processes often show coexistence of fast and slow dynamics during their time evolutions. The increasing demand for a flexible operation of a complex process, a pressing need to improve the product quality, an increasing energy cost and tightening environmental regulations make it rewarding to automate a large scale manufacturing process. Mathematical tools used for process modeling, simulation and control are useful to meet these challenges. Towards this purpose, development of process models, either from the first principles (conservation laws) i.e. the rigorous models or the input-output data based models constitute an important step. Both types of models have their own advantages and pitfalls. Rigorous process models can approximate the process behavior reasonably well. The ability to extrapolate the rigorous process models and the physical interpretation of their states make them more attractive for the automation purpose over the input-output data based identified models. Therefore, the use of rigorous process models and rigorous model based predictive control (R-MPC) for the purpose of online control and optimization of a process is very promising. However, due to several limitations e.g. slow computation speed and the high modeling efforts, it becomes difficult to employ the rigorous models in practise. This thesis work aims to develop a methodology which will result in smaller, less complex and computationally efficient process models from the rigorous process models which can be used in real time for online control and dynamic optimization of the industrial processes. Such methodology is commonly referred to as a methodology of Model (order) Reduction. Model order reduction aims at removing the model redundancy from the rigorous process models. The model order reduction methods that are investigated in this thesis, are applied to two benchmark examples, an industrial glass manufacturing process and a tubular reactor. The complex, nonlinear, multi-phase fluid flow that is observed in a glass manufacturing process offers multiple challenges to any model reduction technique. Often, the rigorous first principle models of these benchmark examples are implemented in a discretized form of partial differential equations and their solutions are computed using the Computational Fluid Dynamics (CFD) numerical tools. Although these models are reliable representations of the underlying process, computation of their dynamic solutions require a significant computation efforts in the form of CPU power and simulation time. The glass manufacturing process involves a large furnace whose walls wear out due to the high process temperature and aggressive nature of the molten glass. It is shown here that the wearing of a glass furnace walls result in change of flow patterns of the molten glass inside the furnace. Therefore it is also desired from the reduced order model to approximate the process behavior under the influence of changes in the process parameters. In this thesis the problem of change in flow patterns as result of changes in the geometric parameter is treated as a bifurcation phenomenon. Such bifurcations exhibited by the full order model are detected using a novel framework of reduced order models and hybrid detection mechanisms. The reduced order models are obtained using the methods explained in the subsequent paragraphs. The model reduction techniques investigated in this thesis are based on the concept of Proper Orthogonal Decompositions (POD) of the process measurements or the simulation data. The POD method of model reduction involves spectral decomposition of system solutions and results into arranging the spatio-temporal data in an order of increasing importance. The spectral decomposition results into spatial and temporal patterns. Spatial patterns are often known as POD basis while the temporal patterns are known as the POD modal coefficients. Dominant spatio-temporal patterns are then chosen to construct the most relevant lower dimensional subspace. The subsequent step involves a Galerkin projection of the governing equations of a full order first principle model on the resulting lower dimensional subspace. This thesis can be viewed as a contribution towards developing the databased nonlinear model reduction technique for large scale processes. The major contribution of this thesis is presented in the form of two novel identification based approaches to model order reduction. The methods proposed here are based on the state information of a full order model and result into linear and nonlinear reduced order models. Similar to the POD method explained in the previous paragraph, the first step of the proposed identification based methods involve spectral decomposition. The second step is different and does not involve the Galerkin projection of the equation residuals. Instead, the second step involves identification of reduced order models to approximate the evolution of POD modal coefficients. Towards this purpose, two different methods are presented. The first method involves identification of locally valid linear models to represent the dynamic behavior of the modal coefficients. Global behavior is then represented by ‘blending’ the local models. The second method involves direct identification of the nonlinear models to represent dynamic evolution of the model coefficients. In the first proposed model reduction method, the POD modal coefficients, are treated as outputs of an unknown reduced order model that is to be identified. Using the tools from the field of system identification, a blackbox reduced order model is then identified as a linear map between the plant inputs and the modal coefficients. Using this method, multiple local reduced LTI models corresponding to various working points of the process are identified. The working points cover the nonlinear operation range of the process which describes the global process behavior. These reduced LTI models are then blended into a single Reduced Order-Linear Parameter Varying (ROLPV) model. The weighted blending is based on nonlinear splines whose coefficients are estimated using the state information of the full order model. Along with the process nonlinearity, the nonlinearity arising due to the wear of the furnace wall is also approximated using the RO-LPV modeling framework. The second model reduction method that is proposed in this thesis allows approximation of a full order nonlinear model by various (linear or nonlinear) model structures. It is observed in this thesis, that, for certain class of full order models, the POD modal coefficients can be viewed as the states of the reduced order model. This knowledge is further used to approximate the dynamic behavior of the POD modal coefficients. In particular, reduced order nonlinear models in the form of tensorial (multi-variable polynomial) systems are identified. In the view of these nonlinear tensorial models, the stability and dissipativity of these models is investigated. During the identification of the reduced order models, the physical interpretation of the states of the full order rigorous model is preserved. Due to the smaller dimension and the reduced complexity, the reduced order models are computationally very efficient. The smaller computation time allows them to be used for online control and optimization of the process plant. The possibility of inferring reduced order models from the state information of a full order model alone i.e. the possibility to infer the reduced order models in the absence of access to the governing equations of a full order model (as observed for many commercial software packages) make the methods presented here attractive. The resulting reduced order models need further system theoretic analysis in order to estimate the model quality with respect to their usage in an online controller setting