Fibonacci wavelets operational matrix approach for solving chemistry problems

The key idea and contribution of this study are to present the innovative functional matrix approach for solving the two chemical, mathematical problems such as the absorption of CO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C{O}_{2}$$\end{document} into phenyl glycidyl ether (PGE) and the chemical kinetic problem. These chemical problems are characterized as a system of nonlinear ODEs with initial and boundary conditions. Numerical outcomes are obtained to establish the simplicity and efficacy of the developed scheme. Graphs and tables show how consistently and effectively the developed strategy works. Results obtained show that the newly adopted technique is more precise and effective than other methods such as Adam–Bashforth–Moultan method (ABM), Runge–Kutta method (RK4), Adomian Decomposition method (ADM), and residual method (RM). Mathematical software called Mathematica 11.3 has been used to perform all calculations. Theorems explain the convergence of this approach.


Introduction
In real life, many variables and parameters are connected under certain circumstances.When we mathematically represent the relationship between these variables and parameters, we usually arrive at a mathematical model of the problem: a differential equation, an integral equation, a system of differential or integral equations, etc. Mathematical modelling is the best way to formulate real-world problems.It is familiar that numerous mathematical characterizations of enormous growth in physical and chemical sciences are represented by SODEs.In the field of chemistry, the above-considered two models that are nonlinear SODEs can describe different kinds of Dirichlet and Neumann-type boundary conditions.Carbon dioxide is vital in plant photosynthesis and the industrial of carbonated soft drinks, removing caffeine from coffee, used in fire extinguishers, and energizing pneumatic models in robotics [1,  2].CO 2 is a beneficial gas that forms one carbon atom and two oxygen atoms.In recent times the chemical fascination of CO 2 is a significant research focus since the hazard postured by global warming and the conversion of CO 2 into valuable substances is a tremendously elegant result.
Consider the general form of the SODEs, which represents the CO 2 absorbed into the PGE problem [3]: Vol:.(1234567890)
The mathematical model of the chemical kinetic problem is one of the well-described nonlinear SODEs that Robertson introduced in 1966.There are three spaces in the model of chemical procedure, which are indicated by P, Q, and R. We can define the three reactions as: The concentration of P, Q, and R can be denoted by u 1 , u 2 and u 3 respectively.It is worth considering that these are the combinations of three concentrations in one.Let b 1 denote the reaction rate of Eq. (2), and this means that the rate at which u 2 rises and at which u 1 reduces, due of the above reaction, it will be equivalent to b 1 u 1 .R acts as a catalyst in the production of P from Q and b 2 represents reaction rate in the reaction Eq. ( 3), meaning that the decrease of u 3 and the increase of u 1 this reaction will be equivalent to b 2 u 2 u 3 .At last, the rate at which this reaction will be equal to b 3 u 2 3 .Due to the production of R from Q will have a constant rate equivalent to b 3 .
Consider the general form of the nonlinear SODEs, which represents the chemical kinetic problem: [4,5] where ( ′ ) represents the derivative of the dependent variable u i concerning the independent variable x and corre- sponding initial conditions; u i (x) x=0 = i , where i 's are the specified constant vector and u i (x) is the solution vector and i = 1, 2, 3, … q .Different published works were carried out on the above models, such as the Chemistry problem by Abbasbandy and Shirzadi [6], the chemistry problem by Matinfar et al. [7], An analytical approximation to the solution of the chemical kinetic problem by Hossein A (1) [4], chemistry problems by Jawary and Raham [8], Kumar et al. [9] implemented a comparative study for fractional chemical kinetics and carbon dioxide CO 2 absorbed into phenyl glycidyl ether problems, CO 2 absorbed into the PGE problem by Robertson [10], Ganji et al. [5] proposed He's Homotopy Perturbation method for the chemistry kinetic problem, Jawary et al. implemented an effective iteration method and VIM (Variational Iteration Method) for the chemical kinetic problem, and the absorption of carbon dioxide into phenyl glycidyl ether problem respectively [11].The ADM method was used for the simple steadystate condensations of CO 2 and PGE, and it is presented in [3], Duan et al. [12] successfully attempted the simple steady-state condensations of CO 2 and PGE by the ADM Method, Kaya [13] attempted the Adomian decomposition method to the chemical kinetic problem, Khader [14] applied the Laplace-Pade approximation technique to improve the performance of the Chemical Kinetic Problem (CKP) model, MSH Chowdary implemented a novel iterative method for the CK System [15], Boundary domain integral method and HAM for nonlinear BVP's were proposed by Jawary et al. [16], Singha and Wazwaz [17] proposed an optimal Homotopy analysis method for steady-state concentrations of absorption of carbon dioxide into phenyl glycidyl ether problem.The wavelet theory is a relatively renewed and emerging theory in mathematical research.Wavelet analysis is a great tool that significantly impacts study and engineering.The primary criteria that draw researchers to the wavelet are orthogonality and effective localization.With the help of wavelets that are mathematical operations, data may be divided into several frequency components, each of which can then be analyzed at a different resolution.It is possible to employ wavelets as a mathematical tool to extract information from various data formats.Wavelet acquired the attention of many mathematics researchers for solving various types of highly nonlinear typical problems.Different wavelet collocation schemes applied for some of the common mathematical problems are Hermite, Legendre, and Laguerre wavelets [18][19][20][21][22][23][24][25][26], Bernoulli wavelets [28], Fractional delay differential equations [27], Fractional Riccati differential equation [29].Few wavelet collocation methods generally used for solving various engineering, physical and biological models include Legendre wavelets [30], Haar wavelets [31][32][33], Chebyshev wavelets [34], Bernoulli wavelets [27,[43][44][45], and Ultraspherical wavelets [35,36].
Fibonacci wavelets produced by Fibonacci polynomials are a new accumulation to the field of wavelet families.It has an added advantage in contrast with the other wavelet methods as follows [39]: • The number of terms of the Fibonacci polynomials P m (x) is less than the number of the terms of the Legendre polynomials L m (x).It helps to reduce CPU time.• Error components in the OMI representing Fibonacci polynomials are less than that of Legendre polynomials.• Fibonacci polynomials P m (x) have less coefficient of indi- vidual terms than in Legendre polynomials L m (x).Com- putational errors can be reduced using this property.• Using the Mathematica command Fibonacci [m, x], the coefficients of the Fibonacci polynomials can be easily obtained in computer programs.
Due to its superior properties and advantages over other wavelets, it grabs the attention of many researchers.As a consequence, researchers start using this package to solve mathematical problems such as the Nonlinear Hunter-Saxton Equation [48], fractional order optimal control problems [37], time-varying delay problems [36], nonlinear Stratonovich Volterra integral equations [38], time-fractional telegraph equations [39], time-fractional bioheat transfer model [40], nonlinear Volterra integral equations [41], Dual-phase lag heat transfer model [42].Also, some articles referred to the improvisation of concept [49][50][51].
This paper is organized as follows.Section 2, named "Preliminaries," defines Fibonacci wavelets.OMI of Fibonacci wavelets carried out in Sect.3. The method of solution and application of the proposed scheme is explained in Sects.4 and 5, respectively.Finally, Sect.6 gives the conclusion of the article.

Preliminaries of the Fibonacci wavelets
On the interval [0, 1] , Fibonacci wavelets are defined as [36,  39], with where P m (x) is the Fibonacci polynomial of degree m = 0, 1, 2 … , M − 1 , t r a n s l a t i o n p a r a m e t e r n = 1, 2, … , 2 k−1 and k represents the level of resolution k = 1, 2, … , respectively.The quantity 1 √ S m is a normalization factor.The Fibonacci polynomials are defined as follows in the form of the recurrence relation for every x ∈ R + : with initial conditions P 0 (x) = 1, P 1 (x) = x.Fibonacci wavelets are compactly supported wavelets formed by Fibonacci polynomials over the interval [0,1].
Basic Definitions [46,47]: Convergence: A real sequence {a n } converges to a limit a if for given  > 0 , there exists Uniform Convergence: Let E be a nonempty subset of ℝ .A sequence of functions f n ∶ E → ℝ is said to converge uniformly on E to a function f if and only if for every  > 0 there is an N ∈ ℕ such that Absolute Continuity: A real-valued function f defined on [a, b] is said to be absolutely continuous on [a, b] if, given  > 0, there is a  > 0 such that Outer Measure: Let A be any set of real numbers, then the outer measure of A written as m * A is defined by , where infimum is taken over the countable covering of A by open intervals.
Theorem 1 [27] Let L 2 [0, 1] be the Hilbert space generated by the Fibonacci wavelet basis.Let (x) be the continuous bounded function in L 2 [0, 1] .Then the Fibonacci wavelet expansion of (x) converges with it.
Proof Let ∶ [0, 1] → R be a continuous function and | (x)| ≤ , where be any real number.Then Fibonacci wavelet dilation of y(x) can be expressed as, a n,m − = ⟨ (x), n,m (x)⟩ denotes inner product: Since n,m are the orthogonal basis: Then substitute 2 k−1 x − n + 1 = y then we get, a n,m is absolutely convergent.Hence the Fibonacci wavelet series expansion (x) converges uniformly to it.
Theorem 2 [27] Let I ⊂ R be a finite interval with length m(I) .Furthermore, f (x) is an integrable function defined on I and j=1 a i,j i,j (x) be a good Fibonacci wavelet approximation of f on I with for some  > 0 ,

Operational matrix of integration (OMI)
At k = 1 and M = 10, the Fibonacci wavelet basis is examined as follows: where, Integrating the first ten bases described above for the range of x limits from 0 to x in the form of a linear combina- tion of Fibonacci wavelet basis.We obtain, x Again, integrating the above ten basis, we obtain, Hence, where In the same way, we can generate matrices of different sizes for our handiness.

Fibonacci wavelet method
Consider the two non-linear reactions in normalized form and boundary conditions will become: where 1 , 2 , 1 , 2 are normalized arguments, z(x) denotes the condensation of PGE and y(x) denotes the condensa- tion of CO 2 .Here, we would like to bring the solution of system (8) that is extracted from the CO 2 absorbed into the PGE problem (1) in the Fibonacci wavelet space.
Next, we will consider the system of ODEs, which represents the Chemical kinetic problem as follows: We want to bring the solution of the SODEs in the Fibonacci wavelet space that is extracted from a chemical kinetic problem (5).
Assume that, where, and the initial conditions are We obtain the following set of equations on integrating the Eqs.( 18)-( 20) with respect ' x ' from ' 0 ' to ' x' .Using Eq. ( 6) and physical constraints in Eq. ( 17) formulated in terms of (x) .We obtain, where F, D, and E are the known vectors.On Substituting Eqs. ( 18), ( 19), (20), and ( 21) in ( 16), we obtain,   On collocating all the equations of ( 22) by using the grid points x i = 2i−1 2 k M , i = 1, 2 … M , we arrive the system with 2 k M algebraic equations as below: .   Utilizing the Newton-Raphson/Secant method, solve the aforementioned equations to determine the values of the unknown Fibonacci wavelet coefficients.The numerical solution of the system (16) using the Fibonacci wavelet is obtained by substituting these coefficient values into (21).( 23) .

Numerical results and discussion
Example 5.1 The non-linear SODEs of the CO 2 absorbed into phenyl glycidyl ether problem is of the form (24)     and the boundary conditions will be ∶ y(0) = , y(1) = , z � (0) = , z(1) = numerical simulations and graphical results of this particular case.We have compared the proposed method with the optimal homotopy analysis method (OHAM) [17], the Adomian decomposition method (ADM) [17], and the ND Solve solution due to the non-availability of the exact solution.It can be observed that the FWCM solution is very much close to the NDSolve solution.The CPU time taken is 1.017 s for the two tables of values together.
Case 2: 3 and 4 and Figs. 5, 6, 7 and 8 depict this case's numerical simulations and graphical results.We have compared the proposed method (FWCM) with the Adomian Decomposition method (ADM), Residual method (RM), and ND solve solution due to the non-availability of the exact solution.It can be observed that the FWCM solution is very close to the NDSolve solution.The CPU time taken is 1.827 s for the two tables of values together.A n d t h e p h y s i c a l c o n s t r a i n t s a r e : u 1 (0) = 1, u 2 (0) = 0, u 3 (0) = 0.
The Exact solutions of the above system are: w 1 (x) = x 3 − x, w 2 (x) = x 2 − x .We compared the solutions obtained from the FWCM and ND Solve methods with the Exact solutions in Tables 11 and 12. Also, we brought the Absolute Errors of FWCM and NDSolve method with the exact solutions in the same tables.

Conclusion
In this article, we have implemented a new collocation approach called the Fibonacci wavelet collocation method (FWCM).The functional matrix of integration is obtained by using the Fibonacci wavelets.We solved two nonlinear SODEs representing the two familiar chemistry models using this matrix.The accuracy and proficiency of this novel technique were demonstrated by comparing the obtained solutions with other standard methods such as RK4 and ND Solve.Tables and Graphs reveal that the FWCM converges rapidly compared to other existing techniques (ADM, RM, HPM, and VIM) in the literature.Also, numerical illustrations support the claim that only a few Fibonacci wavelets are sufficient to obtain satisfactory results.Hence the proposed method is a very attractive and efficient numerical technique to solve different chemical, biological, and physical problems.From Tables 8, 9, and 10 and Figs.15,16, and 17, we can see that the values of k and M incresases the accuracy of the solution.The future recommendation of the proposed scheme is that it is suitable for solutions with sharp edge/jump discontinuities.By slightly modifying the method, the Fibonacci wavelet method can solve the higher-order system of ordinary differential equations and extend to PDEs and other mathematical models with different physical conditions.It is used to obtain the solution of the differential equation in the universal domain by taking the suitable transformation (Figs. 18,19,20,21).

Fig. 1 Fig. 2 AFig. 3
Fig. 1 The plot of FWCM solution at k = 1 for y(x) with different methods of solution

Fig. 5 Fig. 6
Fig.5 The plot of FWCM solution at k = 1 for y(x) with different methods of solution

Fig. 7 Fig. 8
Fig. 7 Absolute error (AE) comparison of y(x) at k = 1 and M = 6 with different methods in the literature

Fig. 9 Fig. 10
Fig.9 The plot of the FWCM solution at k = 1 for u 1 (x) with different methods of solution

Example 5 . 2
Next, Consider the system of ODEs of Chemical Kinetic Problem as follows:

Fig. 11 Fig. 12 Fig. 13 Fig. 14
Fig.11The plot of the FWCM solution at k = 1 for u 3 (x) with different methods of solution observed that the FWCM solution is very much close to the ND Solve solution.It is clear from the Tables and Graphs that the approximate solution obtained from the FWCM is close to the ND Solve solution, which shows the efficiency of the proposed method.Also, we got a better solution by increasing the value of M and k, which can be seen in Tables 8, 9 and 10.Comparison of AE with distinct values

Fig. 15 Example 5 . 3
Fig. 15 AE comparison for u 1 (x) at various values of M and k Fig. 16 AE comparison for u 2 (x) at various values of M and k

Table 1
FWCM solution compared with the various methods for y(x)

Table 2
FWCM solution compared with the various methods for z(x)

Table 3
FWCM solution compared with the various methods for y(x)

Table 6
FWCM solution compared with the various approaches for u 2(x)

Table 7
FWCM solution compared with the various approaches for u 3 (x)