Abstract
Studying and analyzing the random motion of a particle immersed in a liquid represented in the Langevin fractional model by Caputo’s independent derivative is one of the aims of applied physics. In this article, we will attend to a new, accurate, and comprehensive numerical solution to the aforementioned model using the reproducing kernel Hilbert approach. Basically, numerical and exact solutions of the fractional Langevin model are represented using an infinite/finite sum, simultaneously, in the space. The proof has been sketched for many mathematical theorems such as independence, convergence, error behavior, and completeness of the solution. A sufficient set of tabular results and two-dimensional graphs are shown, and absolute/relative error graphs that express the dynamic behavior of the fractional parameters are utilized as well. From an analytical and practical point of view, we noticed that the simulation process and the iterative approach are appropriate, easy, and highly efficient tools for solving the studied model. In conclusion, what we have carried out is presented with a set of recommendations and an outlook on the most important literature used.
1 Introduction
The Langevin equation, in its fractional issue, is a mathematical dynamic model fundamental in Brownian motion applications to characterize the emergence of physical episodes in oscillating mediums. It is a popularization of the conventional model that utilizes a fractional Gaussian procedure formalized by two vertices, which is much more adaptable to the parameterization of the fractal evolution processes [1–4]. The applications of the FLM can be seen in the stock market, motor control system modeling, photoelectron counting, fluid suspensions, deuteron-cluster dynamics, protein dynamics, evacuation process modeling, financial markets, single-file diffusion, and anomalous transport [5–8]. Commonly, FDMs have wide applications in the formation of many engineering and fluid physics phenomena. In the second half of the last century, the search for more serious numerical algorithms began to control many of the nonlinear problems that appeared with the emergence of fractional derivatives and their emerging applications. However, many studies came into existence, as in the following literature [9–13]. To date, effective analytical and numerical schemes have been developed and successfully applied to deal with different classes of FDMs [14–20].
Although finding accurate solutions to FDMs of different orders in Brownian particle motion placements is an important problem for understanding the dynamic attitudes of oscillatory environments in trifling media, in this article, we contemplate creating an accurate numerical solution to FLMs utilizing the CFD with appropriate boundary data using a new renewal in the RKHA. Here, we will generate effective and straightforward numerical solutions without imposing any restrictions on the nature of the proposed FLM and obtain sufficient convergence while reducing the computation time by exposure to the following model [1–8]:
equipped with the posterior boundary condition:
Typically, the FDM (Eqs
1,
2) consists of the posterior parameters, functions, and variable effects as attached:
1) stands for the time-coordinate independent domain.
2) stands for the rank of fractional derivatives applied.
3) stands for the particle position.
4) is the CFD of the rank of and is given as follows:
6) is a bounded variational map that symbolizes the imposed acting on Brownian particles.
7) is a nontrivial parameter that represents the damping or viscosity term.
8) are nontrivial parameters that represent the initial and terminal positions of the particle simultaneously.
Ordinarily, no conventional schemes produce an accurate prototype solution for FLMs. Thus, there is a great need for the RKHA, which satisfies the purpose, as usual, and has reached the desired and satisfactory numerical results with fully stochastic properties. Here, the relevant theories and facts have been confirmed by numeric emulations, graphical representation, and scale tables for three types of FLMs. However, the tenor of the paper is arranged as follows: Section 1: Presentation: FDMs and CFDs. Section 2: RKHA: Preliminaries and definitions. Section 3: RKHA: Construction and properties. Section 4: RKHA: Solutions and convergence. Section 5: Error: Estimation and bound. Section 6: Justifications: Algorithms, applications, and analyses. Section 7: Outline: Conclusion and outlook.
2 Reproducing kernel: Preliminaries and definitions
The approach of the reproducing kernel is a novel solver built to find solutions to FDMs emerging in physics, waves, statistics, and engineering [21–23]. This technique is based on the Gram–Schmidt process and Fourier expansion approach for an arbitrary order and is used for optimizing an orthogonal basis to detect unknown compounds. The RKHA has many motivational aspects and a great ability to handle complicated problems without imposing any restrictions on the style of the models. Therefore, it has been gaining a lot of solicitude and examination lately [24–32].
First, a reproducing kernel induced from a given Hilbert space is called the RKHA. Here, is a set of maps that are continuous absolutely on . At the outset, some requirements that are necessary to go further in our RKHA scheme will be sought.
Remark 1[24] The frameworks of are as follows:
Remark 2[24] The frameworks of are as follows:with .
Definition 1 The framework of is as follows:Here, one can find that is a subset in and is closed. Next, to generate the kernel function of , we put .
If , then the framework function of is as follows:
Proof: Because , . However, for each in , one can obtainor is the reproducing kernel of .
If , then , , and .
Proof: Since , applying successive integration from to for , , and , one can obtainTaking values and using the fact , we can obtainTo complete, one can get , , , and .
3 Reproducing kernel: Construction and properties
Herein, the boundaries in Eq. 2 will first be homogenized to zero to obtain easy-to-access modeling in the proposed space. The form of the operator that formulates the required solution will also be determined in addition to some of what is needed for the scheme.
To achieve this, one must first carry out the following:
equipped with the posterior boundary condition:For the briefing, the normalizing modified version in Eqs 13, 14 was obtained from the posterior underlying conversion, taking into account that all the extra terms transformed into as
The conversions in Eqs
15,
16are needful to insert the equipped boundaries in Eq.
2inside
. Indeed, we will denote
to new and old solutions always.
1) Define the map such that
2) Build the operator as
3) Reframe the FLM problem to solve such that
is a bounded linear operator.
Proof: From Remark 1, one can obtainUsing the reproducing property of , one can obtainWith the use of the Schwarz inequality, one can obtainSo or . By picking out a countable dense subset in , defining , and setting , one can fit the orthogonal function system of . Furthermore, by the Gram–Schmidt process, one can fit the orthonormal function systems on as
is the complete function system of with
Proof: First, indicates that applies to a function of . Certainly,So can be written as in . To illustrate more effective properties in , the set is linearly independent as agrees well with , and for , taking in as , thenSo is linearly independent for each .
4 Reproducing kernel: Solutions and convergence
This section aims to construct exact and RKHA numeric solutions of Eq. 19, together with some convergence theories, to ensure this analysis is more efficient. Here, and indicate a set of continuous maps on the group inside parentheses.
If and are orthonormal, then , are the Fourier maps of for and its Fourier expansion , wherein is dense throughout .
Suppose that
are orthogonalization coefficients of
and a unique solution of Eq.
19exists, then the posteriors are accomplished:
Proof: For the first side: Through Theorem 4, is an orthonormal basis in , which is complete. Using as the Fourier expansion concerning , one can obtain in . SoFor the second side: Because and are the Hilbert space, one can truncate Eq. 27 using the -idiom RKHA numeric solution of to generate Eq. 28.
Remark 3Assume that and are dense throughout . Then, according to Eqs 27–28, the portrayal effect can be determined as follows:
Assume that , as , , and . Then, as .
Proof: It all started where as . Evidently, one can obtainwith in between and . Applying Theorem 2, one can obtain or equivalently as and . Using and , one can obtain equivalently as . By , it gives a glimpse that as .
One gains as .
Proof: Because and the orthogonality of , one can obtainSo and such that or . To getit is adequate to have for thatwith . However, on the other line, as . Using completeness, such that as .
5 Error: Estimation and bound
As a matter of fact, the exact solution of FDMs (Eqs 1, 2) restricted with the CFD depends on and . In most cases, this real problem is complicated to be solved traditionally because it lacks an exact solution. This induces us to develop numerical schemes like RKHA that generate approximated realizations of the exact solution depending on the values of and .
Here, we will fix with , , , and . Furthermore, we will fix to the residual truncated error at .
Lemma 1One gains ,.Proof: Set . SoIn other words, .
Lemma 2Assume that , for a fixed , and at with . So and a parameter exist withProof: Clearly, and , , one can obtainOn , applying Roll’s theorem on yields with , . So for fixed , such that . Similarly, one can writeSo . Sequentially, a parameter exists with and . By compiling the previous results, one can obtainwith .
A parameter exists with
Proof: Utilizing Lemma 3, one can obtainHowever, since , then and a parameter exists with
Applying Theorem 2, one can obtainIn another mode, one can obtain , , where .
Ultimately, one can see that is decreasing for a large as
However, as , one can obtain as .
6 Justifications: Algorithms, applications, and analyses
To clarify the portability and effectiveness of the presented numeric approach, we need some of the steps, the first of which is to provide sufficient algorithms to demonstrate the mechanism of the solution, the second of which is to present several tangible applications, and then finally to provide several tables, figures, and numeric explanations of the solution procedures. However, all of this is the content of the following sections.
6.1 Algorithms
Next, three used algorithms in our RKHA implementation are given. These algorithms are problem initialization, the Gram–Schmidt process, and RKHA solution steps, simultaneously. However, an expert in the Mathematica platform can interpret these steps in the form of programs.
Stride 1: Set the assumptions.
Output: Homogenous FLM.
Stride 2: Define a suitable operator.
Output: Homogenous FLM in the functional form.
Algorithm 1
Stride 1: At and , evaluate
Output: parameters.
Stride 2: At , evaluate
Output: system.
Algorithm 2
Stride I: At fix on evaluating and in .
Output: system.
Stride II: Wherein and , evaluate Gram–Schmidt.
Output: parameters.
Stride III: Set in .
Output: system.
Stride IV: Set in and evaluate
Output: -term numeric approximation of .
Algorithm 3
6.2 Applications
Next, three test applications that coincide with the FLM platform are utilized on the basis of the CFD: the first is , the second is , and the third is . Actually, these applications are solved and analyzed using the presented RKHA as utilized in Algorithm 1, Algorithm 2, Algorithm 3.
Application 1: Theorize the posterior and evaluate :
equipped with the posterior boundary condition:Here, is the exact smooth solution of Eqs 58, 59 on .
Application 2: Theorize the posterior and evaluate :
equipped with the posterior boundary condition:Here, is the exact smooth solution of Eqs 60, 61 on .
Application 3: Theorize the posterior and evaluate :
equipped with the posterior boundary condition:
Here, is the exact smooth solution of (Eqs 62, 63) on .
6.3 Analyses
Eventually, to clarify the portability and effectiveness of the presented RKHA, we present and provide several RKHA solution tables, RKHA numeric solution figures, RKHA absolute error figures, and RKHA relative error figures. However, in this section, we used with in on and then executed Algorithm 3 throughout and its related steps.
The tabulated data presenting
concerning the features of the memory inherited, are utilized in detail by the RKHA performance of the applications addressed previously, as shown in the included tables (Tables 1–3).
TABLE 1
Associated RKHA scores for application 1 with , when .
TABLE 2
Associated RKHA scores for application 2 with , when .
TABLE 3
Associated RKHA scores for application 3 with , when .
The - Cartesian plots presenting , concerning the features of the memory inherited, are utilized in detail by the RKHA performance of the applications addressed previously, as shown in the included graphs (Figures 1A–C).
FIGURE 1
The - Cartesian plots presenting , concerning the features of the memory inherited, are utilized in detail by the RKHA performance of the applications addressed previously, as shown in the included graphs (Figures 2A–C).
FIGURE 2
Ultimately, the - Cartesian plots presenting , concerning the features of the memory inherited, are utilized in detail by the RKHA performance of the applications addressed previously, as shown in the included graphs (Figures 3A–C).
FIGURE 3
7 Outline: Conclusion and outlook
This research embraced the RKHA to handle a type of well-known FDM called FLM on the basis of the CFD by including three test applications. A detailed presentation of the theories related to the establishment of the solution and the formulation of the approximate was interspersed with the construction of the necessary spaces and the associations’ form of Green’s functions used, with many new results that centered on convergence, error, and independence. Based on the RKHA, , , , and have been sketched in - and tabulated for various value parameters of and . Conclusively, obtaining analytical solutions for many types of FDMs is not a simple procedure, which motivates us to conduct more studies and scientific research to obtain innovative approximations of FDMs subject to influence the CFD. The RKHA has diverse feasible and favorable benefits. First, the RKHA is appropriate and delicate since the approximation is very closest to the required solution. Second, by utilizing small terms, we can obtain high accuracy. Third, it is an easy, simple, and soft method to be applied since it does not require sophisticated mathematical tools or an adept professional programmer. Fourth, it is global since it may be utilized to handle different types of fractional complex models. Fifth, the main characteristic of the RKHA is that it may be used with other orthogonal basis sequences. Our outlook study will focus on solving the FLM concerning fuzzy boundaries.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
MA: data curation, investigation, software, methodology, validation, roles/writing—original draft, and writing—review and editing. OA: funding acquisition, investigation, resources, supervision, visualization, and roles/writing—original draft. BM: conceptualization, formal analysis, investigation, project administration, software, and writing—review and editing.
Acknowledgments
The authors are grateful to the Middle East University, Amman, Jordan for the financial support granted to cover the publication fee of this research article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
FLM, fractional Langevin model; RKHA, reproducing kernel Hilbert approach; FDM, fractional differential model; CFD, Caputo’s fractional derivative.
References
1.
BaghaniO. On fractional Langevin equation involving two fractional orders. Commun Nonlinear Sci Numer Simulation (2017) 42:675–81. 10.1016/j.cnsns.2016.05.023
2.
SalemAMsharyN. On the existence and uniqueness of solution to fractional-order Langevin equation. Adv Math Phys (2020) 2020:–11. 10.1155/2020/8890575
3.
SalemAAlneggaM. Fractional Langevin equations with multi-point and non-local integral boundary conditions. Cogent Math Stat (2020) 7:1758361. 10.1080/25742558.2020.1758361
4.
SudsutadWNtouyasSKTariboonJ. Systems of fractional Langevin equations of Riemann-Liouville and Hadamard types. Adv Differ Equ (2015) 2015:235. 10.1186/s13662-015-0566-8
5.
YadavSKumarDNisarKS. A reliable numerical method for solving fractional reaction-diffusion equations. J King Saud Univ – Sci (2021) 33:101320. 10.1016/j.jksus.2020.101320
6.
FangDLiL. Numerical approximation and fast evaluation of the overdamped generalized Langevin equation with fractional noise. Math Model Numer Anal (2020) 54:431–63. 10.1051/m2an/2019067
7.
FazliHSunHGNietoJJ. Fractional Langevin equation involving two fractional orders: Existence and uniqueness revisited. Mathematics (2020) 8:743. 10.3390/math8050743
8.
GuoPZengCLiCChenYQ. Numerics for the fractional Langevin equation driven by the fractional Brownian motion. Fract Calc Appl Anal (2013) 16:123–41. 10.2478/s13540-013-0009-8
9.
MainardiF. Fractional calculus and waves in linear viscoelasticity. UK: Imperial College Press (2010).
10.
ZaslavskyGM. Hamiltonian chaos and fractional dynamics. UK: Oxford University Press (2005).
11.
PodlubnyI. Fractional differential equations. USA: Academic Press (1999).
12.
SamkoSGKilbasAAMarichevOI. Fractional integrals and derivatives theory and applications. USA: Gordon & Breach (1993).
13.
KilbasASrivastavaHTrujilloJ. Theory and applications of fractional differential equations. Amsterdam, Netherlands: Elsevier (2006).
14.
AtanganaABaleanuD. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm Sci (2016) 20:763–9. 10.2298/tsci160111018a
15.
PindzaKMOE. Dynamics of fractional chaotic systems with Chebyshev spectral approximation method. Int J Appl Comput Math (2022) 8:140. 10.1007/s40819-022-01340-2
16.
Avalos-RuizLFGomez-AguilarJFAtanganaAOwolabiKM. On the dynamics of fractional map with power-law, exponential decay and Mittag-Leffler memory. Chaos, Solitons and Fractals (2019) 127:364–88.
17.
OwolabiKMAtanganaA. Computational study of multi-species fractional reaction-diffusion system with ABC operator. Chaos Solitons Fractals (2019) 128:280–9. 10.1016/j.chaos.2019.07.050
18.
AlzabutJBolatYAbdeljawadT. Almost periodic dynamics of a discrete Nicholson’s blowflies model involving a linear harvesting term. Adv Differ Equ (2012) 2012:158. 10.1186/1687-1847-2012-158
19.
AlzabutJTyagiSAbbasS. Discrete fractional-order BAM neural networks with leakage delay: Existence and stability results. Asian J Control (2020) 22:143–55. 10.1002/asjc.1918
20.
ArfanMLashinMMASunthrayuthPShahKUllahAIskakovaKet alOn nonlinear dynamics of COVID-19 disease model corresponding to nonsingular fractional order derivative. Med Biol Eng Comput (2022) 60:3169–85. 10.1007/s11517-022-02661-6
21.
CuiMLinY. Nonlinear numerical analysis in the reproducing kernel space. USA: Nova Science (2009).
22.
BerlinetAAgnanCT. Reproducing kernel Hilbert space in probability and statistics. USA: Kluwer Academic Publishers (2004).
23.
DanielA. Reproducing kernel spaces and applications. Switzerland: Springer (2003).
24.
Abu ArqubO. The reproducing kernel algorithm for handling differential algebraic systems of ordinary differential equations. Math Methods Appl Sci (2016) 39:4549–62. 10.1002/mma.3884
25.
KaraagacBOwolabiKMNisarKS. Analysis and dynamics of illicit drug use described by fractional derivative with Mittag-Leffler kernel. Comput Mater Contin (2020) 65:1905–24. 10.32604/cmc.2020.011623
26.
JiangWChenZ. A collocation method based on reproducing kernel for a modified anomalous subdiffusion equation. Numer Methods Partial Differ Equ (2014) 30:289–300. 10.1002/num.21809
27.
GengFZQianSPLiS. A numerical method for singularly perturbed turning point problems with an interior layer. J Comput Appl Math (2014) 255:97–105. 10.1016/j.cam.2013.04.040
28.
LinYCuiMYangL. Representation of the exact solution for a kind of nonlinear partial differential equation. Appl Math Lett (2006) 19:808–13. 10.1016/j.aml.2005.10.010
29.
AkgülA. A novel method for a fractional derivative with non-local and non-singular kernel. Chaos Solitons Fractals (2018) 114:478–82. 10.1016/j.chaos.2018.07.032
30.
AhmadSUllahAAkgülABaleanuD. Analysis of the fractional tumour-immune-vitamins model with Mittag–Leffler kernel. Results Phys (2020) 19:103559. 10.1016/j.rinp.2020.103559
31.
SiddiqueIAkgülA. Analysis of MHD generalized first problem of Stokes’ in view of local and non-local fractal fractional differential operators. Chaos Solitons Fractals (2020) 140:110161. 10.1016/j.chaos.2020.110161
32.
ZhouaYCuiMLinY. Numerical algorithm for parabolic problems with non-classical conditions. J Comput Appl Math (2009) 230:770–80. 10.1016/j.cam.2009.01.012
Summary
Keywords
fractional Langevin model, reproducing kernel Hilbert approach, fractional differential model, Caputo fractional derivative MSC2020, 65L10, fluid dynamics
Citation
Aal MA, Arqub OA and Maayah B (2022) Hilbert solution, iterative algorithms, convergence theoretical results, and error bound for the fractional Langevin model arising in fluids with Caputo’s independent derivative. Front. Phys. 10:1072746. doi: 10.3389/fphy.2022.1072746
Received
17 October 2022
Accepted
07 November 2022
Published
14 December 2022
Volume
10 - 2022
Edited by
Jordan Yankov Hristov, University of Chemical Technology and Metallurgy, Bulgaria
Reviewed by
Kamal Shah, University of Malakand, Pakistan
Ndolane Sene, Cheikh Anta Diop University, Senegal
Kolade Matthew Owolabi, Federal University of Technology, Nigeria
Updates
Copyright
© 2022 Aal, Arqub and Maayah.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Omar Abu Arqub, o.abuarqub@bau.edu.jo
This article was submitted to Statistical and Computational Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.