Numerical Solutions of Fractional Variable Order Differential Equations via Using Shifted Legendre Polynomials

In this manuscript,an algorithm for the computation of numerical solutions to some variable order fractional differential equations(FDEs)subject to the boundary and initial conditions is developed.We use shifted Legendre polynomials for the required numerical algorithm to develop some operational matrices.Further,operational matrices are constructed using variable order differentiation and integration.We are finding the operationalmatrices of variable order differentiation and integration by omitting the discretization of data.With the help of aforesaid matrices,considered FDEs are converted to algebraic equations of Sylvester type.Finally,the algebraic equations we get are solved with the help of mathematical software like Matlab or Mathematica to compute numerical solutions.Some examples are given to check the proposed method’s accuracy and graphical representations.Exact and numerical solutions are also compared in the paper for some examples.The efficiency of the method can be enhanced further by increasing the scale level.

APPROXIMATE SOLUTION OF INTEGRAL EQUATIONS AND CONVOLUTION INTEGRALS USING LEGENDRE POLYNOMIALS

It has been argued that Chebyshev polynomials are ideal to use as approximating functions to obtain solutions of integral equations and convolution integrals on account of their fast convergence. Using the standard deviation as a measure of the accuracy of the approximation and the CPU time as a measure of the speed, we find that for reasonable accuracy Legendre polynomials are more efficient. ’

On THE NUMERICAL INVERSION OF THE LAPLACE TRANSFORM BY THE USE OF AN OPTIMIZED LEGENDRE POLYNOMIALS

A method for inverting the Laplace transform is presented, using a finite series of the classical Legendre polynomials. The method recovers a real valued function f(t) in a finite interval of the positive real axis when f(t) belongs to a certain class βand requires the knowledge of its Laplace transform F(s) only at a finite number of discrete points on the real axis s>0. The choice of these points will be carefully considered so as to improve the approximation error as well as to minimize the number of steps needed in the evaluations. The method is tested on few examples, with particular emphasis on the estimation of the error bounds involved.

Characterization of Optical Aberrations Induced by Thermal Gradients and Vibrations via Zernike and Legendre Polynomials

For every astronomical instrument, the operating conditions are undoubtedly different from those defined in a setup experiment. Besides environmental conditions, the drives, the electronic cabinets containing heaters and fans introduce disturbances that must be taken into account already in the preliminary design phase. Such disturbances can be identified as being mostly of two types: heat sources/sinks or cooling systems responsible for heat transfer via conduction, radiation, free and forced convection on one side and random and periodic vibrations on the other. For this reason, a key role already from the very beginning of the design process is played by integrated model merging the outcomes based on a Finite Element Model from thermo-structural and modal analysis into the optical model to estimate the aberrations. The current paper presents the status of such model, capable of analyzing the deformed surfaces deriving from both thermo-structural and vibrational analyses and measuring their effect in terms of optical aberrations by fitting them by Zernike and Legendre polynomial fitting respectively for circular and rectangular apertures. The independent contribution of each aberration is satisfied by the orthogonality of the polynomials and mesh uniformity.

Legendre Polynomial Kernel: Application in SVM

In machines learning problems, Support Vector Machine is a method of classification. For non-linearly separable data, kernel functions are a basic ingredient in the SVM technic. In this paper, we briefly recall some useful results on decomposition of RKHS. Based on orthogonal polynomial theory and Mercer theorem, we construct the high power Legendre polynomial kernel on the cube [-1,1]d. Following presentation of the theoretical background of SVM, we evaluate the performance of this kernel on some illustrative examples in comparison with Rbf, linear and polynomial kernels.

Galerkin approximation with Legendre polynomials for a continuous-time nonlinear optimal control problem

We investigate the use of an approximation method for obtaining near-optimal solutions to a kind of nonlinear continuous-time(CT) system. The approach derived from the Galerkin approximation is used to solve the generalized Hamilton-Jacobi-Bellman(GHJB) equations. The Galerkin approximation with Legendre polynomials(GALP) for GHJB equations has not been applied to nonlinear CT systems. The proposed GALP method solves the GHJB equations in CT systems on some well-defined region of attraction. The integrals that need to be computed are much fewer due to the orthogonal properties of Legendre polynomials, which is a significant advantage of this approach. The stabilization and convergence properties with regard to the iterative variable have been proved.Numerical examples show that the update control laws converge to the optimal control for nonlinear CT systems.

Duality between Bessel Functions and Chebyshev Polynomials in Expansions of Functions

In expansions of arbitrary functions in Bessel functions or Spherical Bessel functions, a dual partner set of polynomials play a role. For the Bessel functions, these are the Chebyshev polynomials of first kind and for the Spherical Bessel functions the Legendre polynomials. These two sets of functions appear in many formulas of the expansion and in the completeness and (bi)-orthogonality relations. The analogy to expansions of functions in Taylor series and in moment series and to expansions in Hermite functions is elaborated. Besides other special expansion, we find the expansion of Bessel functions in Spherical Bessel functions and their inversion and of Chebyshev polynomials of first kind in Legendre polynomials and their inversion. For the operators which generate the Spherical Bessel functions from a basic Spherical Bessel function, the normally ordered (or disentangled) form is found.

A Chebyshev/Legendre polynomial interpolation approach for fingerprint orientation estimation smoothing and prediction

We introduce a novel coarse ridge orientation smoothing algorithm based on orthogonal polynomials, which can be used to estimate the orientation field (OF) for fingerprint areas of no ridge information. This method does not need any base information of singular points (SPs). The algorithm uses a consecutive application of filtering-and model-based orientation smoothing methods. A Gaussian filter has been employed for the former. The latter conditionally employs one of the orthogonal polynomials such as Legendre and Chebyshev type I or II, based on the results obtained at the filtering-based stage. To evaluate our proposed method, a variety of exclusive fingerprint classification and minutiae-based matching experiments have been conducted on the fingerprint images of FVC2000 DB2, FVC2004 DB3 and DB4 databases. Results showed that our proposed method has achieved higher SP detection, classification, and verification performance as compared to competing methods.

Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

Geometric Covariance Modeling for Surface Variation of Compliant Parts Based on Hybrid Polynomial Approximation and Spectrum Analysis

Part variation characterization is essential to analyze the variation propagation in flexible assemblies. Aiming at two governing types of surface variation,warping and waviness,a comprehensive approach of geometric covariance modeling based on hybrid polynomial approximation and spectrum analysis is proposed,which can formulate the level and the correlation of surface variations accurately. Firstly,the form error data of compliant part is acquired by CMM. Thereafter,a Fourier-Legendre polynomial decomposition is conducted and the error data are approximated by a Legendre polynomial series. The weighting coefficient of each component is decided by least square method for extracting the warping from the surface variation. Consequently,a geometrical covariance expression for warping deformation is established. Secondly,a Fourier-sinusoidal decomposition is utilized to approximate the waviness from the residual error data. The spectrum is analyzed is to identify the frequency and the amplitude of error data. Thus,a geometrical covariance expression for the waviness is deduced. Thirdly,a comprehensive geometric covariance model for surface variation is developed by the combination the Legendre polynomials with the sinusoidal polynomials. Finally,a group of L-shape sheet metals is measured along a specific contour,and the covariance of the profile errors is modeled by the proposed method. Thereafter,the result is compared with the covariance from two other methods and the real data. The result shows that the proposed covariance model can match the real surface error effectively and represents a tighter approximation error compared with the referred methods.

Operator Methods and SU(1,1) Symmetry in the Theory of Jacobi and of Ultraspherical Polynomials

Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting properties. It leads to an alternative definition of the Ultraspherical polynomials by a fixed integral operator in application to powers of the variable u in an analogous way as it is possible for Hermite polynomials. From this follows a generating function which is apparently known only for the Legendre and Chebyshev polynomials as their special case. Furthermore, we show that the Ultraspherical polynomials form a realization of the SU(1,1) Lie algebra with lowering and raising operators which we explicitly determine. By reordering of multiplication and differentiation operators we derive new operator identities for the whole set of Jacobi polynomials which may be applied to arbitrary functions and provide then function identities. In this way we derive a new “convolution identity” for Jacobi polynomials and compare it with a known convolution identity of different structure for Gegenbauer polynomials. In short form we establish the connection of Jacobi polynomials and their related orthonormalized functions to the eigensolution of the Schrödinger equation to Pöschl-Teller potentials.

A two-dimensional lithospheric magnetic anomaly field model of Egypt using the measurements from Swarm satellites

We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct,for the first time,a two-dimensional(2 D)lithospheric magnetic anomaly model of Egypt and its surrounding area.Nighttime data during quiet geomagnetic conditions has been expanded in terms of the Legendre polynomial in harmonic terms N=6-50.The damped least square method has been used to estimate the model coefficients based on the lithospheric magnetic data.Modeled data at two different altitudes(438-448 km and 503-511 km)were compared with the CHAOS model.Results exhibit that the 2 D model is superior to the CHAOS model in the capability of extracting more information about small-scale crustal anomaly field.At low altitudes(438-448 km),the strength of the anomaly field increases,but the noise of the external fields has greatly reduced at high altitudes(503-511 km).Besides,the magnetic anomaly field at low altitudes has illuminated short-scale anomalies that didn’t appear at high altitudes.Both the total and vertical magnetic anomaly vectors showed their ability to reveal tectonic structures compared with Moho depth map and the geological maps.

A Unified Petrov-Galerkin Spectral Method and Fast Solver for Distributed-Order Partial Differential Equations

Fractional calculus and fractional-order modeling provide effective tools for modeling and simulation of anomalous diffusion with power-law scalings.In complex multi-fractal anomalous transport phenomena,distributed-order partial differential equations appear as tractable mathematical models,where the underlying derivative orders are distributed over a range of values,hence taking into account a wide range of multi-physics from ultraslow-to-standard-to-superdiffusion/wave dynamics.We develop a unified,fast,and stable Petrov–Galerkin spectral method for such models by employing Jacobi poly-fractonomials and Legendre polynomials as temporal and spatial basis/test functions,respectively.By defining the proper underlying distributed Sobolev spaces and their equivalent norms,we rigorously prove the well-posedness of the weak formulation,and thereby,we carry out the corresponding stability and error analysis.We finally provide several numerical simulations to study the performance and convergence of proposed scheme.

Model of the Nerve Impulse with Account of Mechanosensory Processes: Stationary Solutions.

Mechanotransduction refers to a physiological process by which mechanical forces, such as pressures exerted by ionized fluids on cell membranes and tissues, can trigger excitations of electrical natures that play important role in the control of various sensory (i.e. stimuli-responsive) organs and homeostasis of living organisms. In this work, the influence of mechanotransduction processes on the generic mechanism of the action potential is investigated analytically, by considering a mathematical model that consists of two coupled nonlinear partial differential equations. One of these two equations is the Korteweg-de Vries equation governing the spatio-temporal evolution of the density difference between intracellular and extracellular fluids across the nerve membrane, and the other is Hodgkin-Huxley cable equation for the transmembrane voltage with a self-regulatory (i.e. diode-type) membrane capacitance. The self-regulatory feature here refers to the assumption that membrane capacitance varies with the difference in density of ion-carrying intracellular and extracellular fluids, thus ensuring an electromechanical feedback mechanism and consequently an effective coupling of the two nonlinear equations. The exact one-soliton solution to the density-difference equation is obtained in terms of a pulse excitation. With the help of this exact pulse solution the Hodgkin-Huxley cable equation is shown to transform, in steady state, to a linear eigenvalue problem some bound states of which can be obtained exactly. Few of such bound-state solutions are found analytically.

On Polynomial Maximum EntropyMethod for ClassicalMoment Problem

The maximum entropy method for the Hausdorff moment problem suffers from ill conditioning as it uses monomial basis{1,x,x2,···,xn}.Themaximum entropy method for the Chebyshev moment probelm was studied to overcome this drawback in[4].In this paper we review and modify the maximum entropy method for the Hausdorff and Chebyshev moment problems studied in[4]and present the maximum entropy method for the Legendre moment problem.We also give the algorithms of converting the Hausdorff moments into the Chebyshev and Lengendre moments,respectively,and utilizing the corresponding maximum entropy method.

Discrete unified gas kinetic scheme for multiscale anisotropic radiative heat transfer

In this work,a discrete unified gas kinetic scheme(DUGKS)is developed for radiative transfer in anisotropic scattering media.The method is an extension of a previous one for isotropic radiation problems[1].The present scheme is a finite-volume discretization of the anisotropic gray radiation equation,where the anisotropic scattering phase function is approximated by the Legendre polynomial expansion.With the coupling of free transport and scattering processes in the reconstruction of the flux at cell interfaces,the present DUGKS has the nice unified preserving properties such that the cell size is not limited by the photon mean free path even in the optical thick regime.Several one-and two-dimensional numerical tests are conducted to validate the performance of the present DUGKS,and the numerical results demonstrate that the scheme is a reliable method for anisotropic radiative heat transfer problems.

A Robust Hybrid Spectral Method for Nonlocal Problems with Weakly Singular Kernels

In this paper,we propose a hybrid spectral method for a type of nonlocal problems,nonlinear Volterra integral equations(VIEs)of the second kind.The main idea is to use the shifted generalized Log orthogonal functions(GLOFs)as the basis for the first interval and employ the classical shifted Legendre polynomials for other subintervals.This method is robust for VIEs with weakly singular kernel due to the GLOFs can efficiently approximate one-point singular functions as well as smooth functions.The well-posedness and the related error estimates will be provided.Abundant numerical experiments will verify the theoretical results and show the high-efficiency of the new hybrid spectral method.

On the solution of the Human Immunodeficiency Virus(HIV)infection model using spectral collocation method

In this research work,we study the Human Immunodeficiency Virus(HIV)infection on helper T cells governed by a mathematical model consisting of a system of three first-order nonlinear differential equations.The objective of the analysis is to present an approximate mathematical solution to the model that gives the count of the numbers of uninfected and infected helper T cells and the number of free virus particles present at a given instant of time.The system of nonlinear ODEs is converted into a system of nonlinear algebraic equations using spectral collocation method with three different basis functions such as Chebyshev,Legendre and Jacobi polynomials.Some factors such as the production of helper T cells and infection of these cells play a vital role in infected and uninfected cell counts.Detailed error analysis is done to compare our results with the existing methods.It is shown that the spectral collocation method is a very reliable,efficient and robust method of solution compared to many other solution procedures available in the literature.All these results are presented in the form of tables and figures.

Navier-Stokes Spectral Solver in a Finite Cylinder

A primitive variable spectral method for simulating incompressible viscous flows inside a finite cylinder is presented.One element of originality of the proposed method is that the radial discretization of the Fourier coefficients depends on the Fourier mode,its dimension decreasing with the increase of the azimuthal modal number.This principle was introduced independently by Matsushima and Marcus and by Verkley in polar coordinates and is adopted here for the first time to formulate a 3D cylindrical Galerkin projection method.A second element of originality is the use of a special basis of Jacobi polynomials introduced recently for the radial dependence in the solution of Dirichlet problems.In this basis the radial operators are represented by matrices of minimal sparsity-diagonal stiffness and tridiagonal mass-provided here in closed form for the first time,and lead to a Helmholtz operator characterized by a favorable condition number.Finally,a new method is presented for eliminating the singular behaviour of the solution originated by the rotation of the lid with respect to the cylindrical wall.Thanks to these elements,the resulting Navier-Stokes spectral solver guarantees the differentiability to any order of the solution in the entire computational domain and does not suffer from the time-step stability restriction occurring in spectral methods with a point clustering close to the axis.Several test examples are offered that demonstrate the spectral accuracy of the solution method under different representative conditions.

THE NUMERICAL SOLUTION FOR A PARTIAL INTEGRO-DIFFERENTIAL EQUATION WITH A WEAKLY SINGULAR KERNEL

In this paper, a first order semi-discrete method of a partial integro-differential equation with a weakly singular kernel is considered. We apply Galerkin spectral method in one direction, and the inversion technique for the Laplace transform in another direction, the result of the numerical experiment proves the accuracy of this method.