Abstract
The Lie symmetry method is applied, and exact homotopic solutions of a non-linear double-diffusion problem are obtained. Additionally, we derived Lie point symmetries and corresponding transformations for equations representing heat and mass transfer in a thin liquid film over an unsteady stretching surface, using MAPLE. We used these symmetries to construct new (Lie) similarity transformations that are different from those that are commonly used for flow and mass transfer problems. These new (Lie) similarity transformations map the partial differential equations of a mathematical model under consideration to ordinary differential equations along with boundary conditions. Lie similarity transformations are shown to lead to new solutions for the considered flow problem. These solutions are obtained using the homotopy analysis method to analytically solve the ordinary differential equations that resulted from the reduction of considered flow equations through Lie similarity transformations. With the aid of these solutions, effects of various parameters on the flow and heat transfer are discussed and presented graphically in this study.
1 Introduction
Fluid flow and heat transfer phenomena have a wide range of applications in engineering. By varying these transporters and enforcing various physical conditions, it is possible to produce a variety of industrial products at their best. As a result, it has drawn a significant amount of attention during the past several decades. The Navier–Stokes equations are used to quantitatively represent these heat and flow exchanges, with the appropriate circumstances. These are extremely non-linear partial differential equations (PDEs) of order two or higher. Such non-linearities lessen the likelihood of obtaining precise results. As a result, flow studies are related to approximation techniques and analytical solution schemes, and heat transfer techniques are frequently used.
The Runge–Kutta and shot method are combined for the derivation of the former type of solutions, whereas homotopy analysis and perturbation techniques are frequently used for the latter.
These problem-solving methods are not directly related to the PDEs that describe the flow problems. The system of ordinary differential equations (ODEs) relating to these flow issues is, nevertheless, solved using these methods. The similarity transformation is the technology that makes this kind of reduction possible. The dependent and/or independent variables of flow equations are reduced using these adjustments.
First, the fact that there are more established and diverse solution methods for ODEs than PDEs accounts for this reduction. Second, running ODEs through mathematical symbolic and numeric software requires less time and equipment compared to other approaches. Following the reduction of flow equations to ODEs via similarity transformations, one finds several applications of such solution algorithms in the literature.
With this procedure, the flow and heat transfers have been studied under different sets of conditions, for example, in a liquid film on an unsteady stretching surface [1, 2], under the effects of variable fluid properties and thermo capillarity [3], with Soret and Dufour effects on a viscoelastic fluid in three dimensions [4], in a rotating channel three-dimensional squeezing flow [5], in a three-dimensional flow of a nanofluid over a non-linearly stretching sheet [6], and for an Oldroyd-B nanofluid thin film over an unsteady stretching sheet [7]. Likewise, magnetohydrodynamic (MHD) flow and heat transfer have been studied for the following: thermosolutal Marangoni convection with heat generation [8], viscoelastic fluid flow over a vertical stretching sheet under the effects of Soret and Dufour [9], Jeffrey fluid over a stretching sheet considering the chemical reaction and thermal radiation [10], three-dimensional flow of an Oldroyd-B nanofluid on a radiative surface [11], thermally radiative flow in three dimensions of a Jeffrey nanofluid under internal heat generation [12], a shrinking sheet with thermal slip [13], a vertical stretching sheet under the effects of heat sink or source [14], mixed convection on the inclined stretching plate in the Darcy porous medium with a Soret effect considering variable surface conditions [15], and mixed convective flow of a Maxwell nanofluid past a porous vertical stretching sheet with a chemical reaction [16].
There are countless studies through an area of research known as the Lie symmetry method, which helps to accurately derive the analytical or approximate solutions for flow and heat transfer equations. For instance, Lie group theory has been employed to study the flow and heat transfer in a non-Newtonian fluid over a stretching surface with thermal radiation [17], MHD boundary layer flow over a stretching sheet with viscous dissipation and uniform heat source/sink [18], MHD mixed flow of unsteady convection on a vertical porous plate with radiation [19], MHD double-diffusion convection of a Casson nanofluid on a vertical stretching/shrinking surface under the effects of thermal radiation and chemical reaction [20], heat flux effect on MHD second slip flow past a stretching sheet along with heat generation [21], MHD Casson fluid flow near a stagnation point on a linearly stretching sheet taking variable viscosity and thermal conductivity into account [22], thermophysical properties of a magnetized Williamson fluid subject to porous/non-porous surfaces [23], two-parameter Lie scaling approach on an unsteady MHD Casson fluid over a porous rigid plate with a stagnation point flow [24], double-diffusive MHD tangent hyperbolic fluid flow on a stretching sheet [25], MHD thermally slip Carreau fluid subject to multiple flow regimes [26], and for a liquid film on an unsteady stretching sheet using Lie point symmetries [27].
The governing equations in the aforementioned flow models are non-linear. Therefore, numerous approaches are adopted to deal with the non-linearity of the governing equations. The Lie symmetry method is one of those that provide a systematic procedure to construct similarity transformation that is a pivotal component of solution schemes employed on fluid flows mentioned previously. Non-linear phenomena impose constraints on the studies conducted to analyze physical models appearing in numerous applications due to the availability of few techniques that are employed to deal with it. As far as the Lie approach is concerned, one may linearize the governing equations (28)–(31). There are many non-Lie procedures that are also available in the literature, for example, effective treatments of the non-linearity of differential equations have been reported in [32–34].
A Lie point symmetry transformation can be associated with a differential or an algebraic equation if it leaves it form invariant. It implies that a heat equation remains a heat equation after mapping it under its Lie point transformation. Every Lie point transformation possesses a Lie symmetry generator. For basic theory and the algebraic computations of the Lie symmetry generators and transformation, readers are referred to [35, 36]. MAPLE contains all these procedures to build symmetry transformations in the “PDEtools” package, which, on applying “Infinitesimals” on differential equations, reveals their symmetries. MAPLE is used to find out symmetry generators and corresponding transformations for flow problems that are being taken into consideration in this study.
We deduce Lie point symmetries for the momentum, energy, and concentration equations representing the flow problem under consideration. There exist nine Lie symmetries, and by using them, Lie similarity transformations are obtained. However, we employ only those symmetries which leave the associated boundary conditions in a particular form. Based on these constraints, we consider three linear combinations (that are also Lie point symmetries) of the derived Lie symmetries. In one of these, we combine two symmetries, while the remaining two consist of three symmetries. These three combinations provided a different type of similarity transformation which transformed flow equations into three different types of ODE systems. Arbitrary constants are used in the linear combinations of the Lie point symmetries, and these constants also appear in the resulting system of ODEs due to their presence in the Lie similarity transformations we construct. We use them to control the convergence of solutions of the flow model we are considering.
The outline of the paper is as follows. The second section is about flow equations and their Lie symmetries. The subsequent section is on similarity transformations and mapping of flow PDEs to ODEs. In the fourth section, analytical solutions are constructed and presented with graphs and tables. The last section is the conclusion.
2 Flow equations
The flow of heat and mass in a thin liquid film has been studied [37] on an unsteady stretching surface with thermosolutal capillarity and variable magnetic field. Here, we are considering the flow model without the magnetic field and thermosolutal capillarity. The governing equations for the flow of heat and mass transfer in a thin liquid film over an unsteady surface are given by the following system of PDEs:subject to boundary conditions as follows:
The Lie point symmetries of the flow mathematical model (Eq. 1) are derived by using the MAPLE “PDEtools” package and the built-in command “Infinitesimals.”
However, for a detailed algebraic procedure to obtain symmetries of system (Eq. 1), the reader is referred to [27]. The Lie symmetry transformations corresponding to symmetry generators (Eq. 3) leave equations of system (Eq. 1) form invariant. These Lie transformations are given in Table 1. Furthermore, all the associated conditions (Eq. 2) should also remain invariant. For this purpose, we employ eachwhere denotes the extension of the symmetry generator; here, we require the first extension of , for m , and denotes the conditions (Eq. 2) for , e.g., , and vice versa.
Table 1
| Generator | Transformation |
|---|---|
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , | |
| , , , , , , |
Lie symmetry generators and transformations.
The symmetry generator from (3) is applied to each of these conditions through the invariance criterion.
3 Lie similarity transformations of flow equations
We construct the Lie similarity transformations corresponding to a few linear combinations for the derived Lie point symmetries . These combinations are based on the unknown functions they determine for , , , and . In this work, only those cases are of interest in which all these functions remain dependent on their arguments. Hence, we consider the combination of Lie symmetries in Case-I, in Case-II, and in Case-III, where , , and are any non-zero real numbers. All other symmetries from the list (Eq. 3) are not suitable in any form to construct the similarity transformations due to stretching sheet velocity and temperature obtained for these symmetries and their combinations. Hence, we consider only those linear combinations that are mentioned previously. These three linear combinations of symmetries leave both and in the stretching sheet velocity and temperature ; i.e., we want to keep them as functions of time and space variable . Moreover, is also left as a function of .
In the study conducted earlier on this type of fluid and heat transports [38], both the said quantities are set to be dependent on both and .
3.1 Case-I: Similarity transformations for
These symmetry generators provided the similarity transformationswhich map the system of PDEs (Eq. 1) into the following system of ODEs:where is the new independent variable. The associated boundary conditions are
3.2 Case-II: Similarity transformations for
In this case, the following similarity transformations are obtained:
These similarity transformations map the system of PDEs (Eq. 1) into the following system of ODEs:and the associated boundary conditions are given as follows:
3.3 Case-III: Similarity transformations for
Here, we obtain the following similarity transformations:
These similarity transformations map the system of PDEs (Eq. 1) into the following system of ODEs:
The associated boundary conditions map to
4 Analytic solution by the homotopy analysis method
In this section, the velocity and temperature profiles are constructed with the aid of the analytical solution of order ten derived through the HAM. It has been observed that the first equation in all three cases that are under consideration here is the same. First, we draw -curves that are presented graphically for in Figure 1. The reason to consider this range is the dimensionless film thickness which remains negative or zero for . Hence, all the velocity, temperature, and concentration profiles are presented here for . The dimensionless film thickness increases with an increase in , under the conditions provided by Lie similarity conditions. This situation changes and opposite trends have been found in [39] using Lie similarity transformations with an introduction of a magnetic term. Figure 2 shows the velocity profiles for the same range of an unsteadiness parameter, which shows an increase in the velocity with this parameter. The temperature and concentration profiles are expected to be different in all three cases as, apparently, the second and third equations in the systems of ODEs (Eq. 6), (Eq. 9), and (Eq. 12) are different. Hence, they are written separately in the following cases to present the trends that are followed by these quantities under the influence of , , and . Moreover, the constants , , and that are used in forming the linear combinations of the Lie symmetry generators (Eq. 3) also affect the temperature and concentration profiles. These are all present in the second and third equations of the systems in Case (3.1)–(3.3).
FIGURE 1
FIGURE 2
4.1 Velocity and concentration profiles for Case-I
For system (Eq. 6), we draw the -curves in Figures 3A–C, for and for three different values of . The -curves show a decline for From these curves, we select a value for to construct the temperature profiles in Figures 4A–C, which also exhibit a decreasing trend with a decrease in the values of Likewise, we draw curves in Figures 5A–C for and for multiple values of . These figures show a decrease in the curves with a decrease in and an increase in . The concentration profiles behave in a similar manner as curves. Here, we present these profiles for and a variation in and The temperature and concentration profiles follow the same trends as system (Eq. 6) equations for both are the same; however, we are presenting them here separately. In both the mentioned set of figures, we considered different values of the unsteadiness parameter . It can be observed from these figures that the unsteadiness parameter and concentration are inversely proportional, i.e., .
FIGURE 3
FIGURE 4
FIGURE 5
4.2 Velocity and concentration profiles for Case-II
System (Eq. 9) involves three arbitrary constants , and , which appear here due to the linear combination of Lie point symmetries we used to construct the corresponding Lie similarity transformation. We draw common curves for and as -curves for this system in Figures 6A–C. These curves are drawn for and a variation in the unsteadiness parameter and , and , and a range of and . These curves and corresponding set of graphs for temperature and concentration show an increase when the unsteadiness parameter decreases from to Similar is the case when goes from to , as shown in Figures 7A–C and Figures 8A–8C.
FIGURE 6
FIGURE 7
FIGURE 8
4.3 Velocity and concentration profiles for Case-III
System (Eq. 12) involves three arbitrary constants , and that are also part of the associated Lie similarity transformation. Figures 9A–C show the -curves for both and These curves are constructed with the same values of as in the previous case and for a different value of the unsteadiness parameter . When the unsteadiness parameter decreases from to , the - and -curves are decreasing. Similar behavior is shown by temperature and concentration profiles in Figures 10A–C; that is, for , the temperature and concentration are increasing. However, for , a decrease in the temperature and concentration is evident from these figures.
FIGURE 9
FIGURE 10
5 Conclusion
Lie point symmetries for heat and mass transfer in a thin liquid film on an unsteady stretching sheet are derived. These symmetries are used to construct Lie similarity transformations which map the PDEs representing the heat and flow model to ODE systems. We showed that there exist three different types of such reductions of the considered flow equations. In the Lie similarity transformation derivation, linear combinations of Lie symmetry generators are utilized. These linear combinations are derived with the help of arbitrary constants, which gives rise to multiple solutions of the flow and heat equations. We use the HAM to analytically solve the obtained non-linear ODEs with a -order of approximation. Velocity, temperature, and concentration profiles are drawn with the aid of these -order HAM solutions. These profiles are presented graphically with variations in the unsteadiness parameter , Prandtl number , Schmidt number , and the arbitrary constants used in the linear combinations of the Lie point symmetries.
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
RK designed and analyzed the results; ST prepared the figures and was involved in discussion; SA was responsible for software and coding; IK analyzed and wrote the manuscript and discussed the results; and SE was responsible for software and coding, supervised the research, and acquired funding.
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.
References
1.
AnderssonHIAarsethJBDandapatBS. Heat transfer in a liquid film on an unsteady stretching surface. Int J Heat Mass Transfer (2000) 43(1):69–74. 10.1016/s0017-9310(99)00123-4
2.
WangC. Analytic solutions for a liquid film on an unsteady stretching surface. Heat Mass Transfer (2006) 42(8):759–66. 10.1007/s00231-005-0027-0
3.
DandapatBSantraBVajraveluK. The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet. Int J Heat Mass Transfer (2007) 50(5-6):991–6. 10.1016/j.ijheatmasstransfer.2006.08.007
4.
HayatTSafdarAAwaisMMesloubS. Soret and Dufour effects for three-dimensional flow in a viscoelastic fluid over a stretching surface. Int J Heat Mass Transfer (2012) 55(7-8):2129–36. 10.1016/j.ijheatmasstransfer.2011.12.016
5.
MunawarSMehmoodAAliA. Three-dimensional squeezing flow in a rotating channel of lower stretching porous wall. Comput Math Appl (2012) 64(6):1575–86. 10.1016/j.camwa.2012.01.003
6.
KhanJAMustafaMHayatTAlsaediA. Three-dimensional flow of nanofluid over a non-linearly stretching sheet: An application to solar energy. Int J Heat Mass Transfer (2015) 86:158–64. 10.1016/j.ijheatmasstransfer.2015.02.078
7.
ZhangYZhangMBaiY. Flow and heat transfer of an Oldroyd-B nanofluid thin film over an unsteady stretching sheet. J Mol Liquids (2016) 220:665–70. 10.1016/j.molliq.2016.04.108
8.
ZhangYZhengL. Analysis of MHD thermosolutal Marangoni convection with the heat generation and a first-order chemical reaction. Chem Eng Sci (2012) 69(1):449–55. 10.1016/j.ces.2011.10.069
9.
RashidiMMAliMRostamiBRostamiPXieGN. Heat and mass transfer for MHD viscoelastic fluid flow over a vertical stretching sheet with considering Soret and Dufour effects. Math Probl Eng (2015) 2015:1–12. 10.1155/2015/861065
10.
NarayanaPSBabuDH. Numerical study of MHD heat and mass transfer of a Jeffrey fluid over a stretching sheet with chemical reaction and thermal radiation. J Taiwan Inst Chem Eng (2016) 59:18–25. 10.1016/j.jtice.2015.07.014
11.
ShehzadSAbdullahZAbbasiFHayatTAlsaediA. Magnetic field effect in three-dimensional flow of an Oldroyd-B nanofluid over a radiative surface. J Magnetism Magn Mater (2016) 399:97–108. 10.1016/j.jmmm.2015.09.001
12.
ShehzadSAbdullahZAlsaediAAbbasiFHayatT. Thermally radiative three-dimensional flow of Jeffrey nanofluid with internal heat generation and magnetic field. J Magnetism Magn Mater (2016) 397:108–14. 10.1016/j.jmmm.2015.07.057
13.
AhmadSYousafMKhanAZamanG. Magnetohydrodynamic fluid flow and heat transfer over a shrinking sheet under the influence of thermal slip. Heliyon (2018) 4(10):e00828. 10.1016/j.heliyon.2018.e00828
14.
AlarifiIMAbokhalilAOsmanMLundLAyedMBelmabroukHet alMHD flow and heat transfer over vertical stretching sheet with heat sink or source effect. Symmetry (2019) 11(3):297. 10.3390/sym11030297
15.
MandalBBhattacharyyaKBanerjeeAKumar VermaAKumar GautamA. MHD mixed convection on an inclined stretching plate in Darcy porous medium with Soret effect and variable surface conditions. Nonlinear Eng (2020) 9(1):457–69. 10.1515/nleng-2020-0029
16.
DessieHFisshaD. MHD mixed convective flow of Maxwell nanofluid past a porous vertical stretching sheet in presence of chemical reaction. Appl Appl Math Int J (2020) 15(1):31.
17.
TufailMN. Group theoretical analysis of non-Newtonian fluid flow, heat and mass transfer over a stretching surface in the presence of thermal radiation. Appl Fluid Mech (2016) 9:1515. 10.18869/acadpub.jafm.68.228.24069
18.
MetriPGGuarigliaESilvestrovS. Lie group analysis for MHD boundary layer flow and heat transfer over stretching sheet in presence of viscous dissipation and uniform heat source/sink. AIP Conf Proc (2017) 1798:020096. 10.1063/1.4972688
19.
RajputGRPatilVJadhavB. MHD mixed flow of unsteady convection with radiation over a vertical porous plate: Lie group symmetry analysis. Appl Comput Mech (2017) 11(2):1–12. 10.24132/acm.2017.380
20.
PalDRoyN. Lie group transformation on MHD double-diffusion convection of a Casson nanofluid over a vertical stretching/shrinking surface with thermal radiation and chemical reaction. Int J Appl Comput Math (2018) 4(1):13–23. 10.1007/s40819-017-0449-7
21.
MegahedAMAbdel-RahmanRG. Lie group analysis of heat flux effect on MHD second slip flow for a slightly rarefied gas past a stretching sheet with heat generation. Tech Sci (2019) 22(1):45–59. 10.31648/ts.4347
22.
Nazim TufailMZaibF. Symmetry analysis of MHD Casson fluid flow for heat and mass transfer near a stagnation point over a linearly stretching sheet with variable viscosity and thermal conductivity. Heat Transfer (2021) 50(6):5418–38. 10.1002/htj.22131
23.
RehmanKUShatanawiWAbodayehK. Thermophysical aspects of magnetized Williamson fluid flow subject to both porous and non-porous surfaces: A lie symmetry analysis. Case Stud Therm Eng (2021) 28:101688. 10.1016/j.csite.2021.101688
24.
SaleemMTufailMNChaudhryQA. Unsteady MHD Casson fluid flow with heat transfer passed over a porous rigid plate with stagnation point flow: Two-parameter Lie scaling approach. Pramana (2021) 95(1):28–9. 10.1007/s12043-020-02054-0
25.
ZebSKhanSUllahZYousafMKhanIAlshammariNet alLie group analysis of double diffusive MHD tangent hyperbolic fluid flow over a stretching sheet. Math Probl Eng (2022) 2022:1–14. 10.1155/2022/9919073
26.
RehmanKUShatanawiWAbodayehK. A group theoretic analysis on heat transfer in MHD thermally slip Carreau fluid subject to multiple flow regimes (MFRs). Case Stud Therm Eng (2022) 30:101787. 10.1016/j.csite.2022.101787
27.
SafdarMIjaz KhanMTajSMalikMShiQH. Construction of similarity transformations and analytic solutions for a liquid film on an unsteady stretching sheet using lie point symmetries. Chaos, Solitons and Fractals (2021) 150:111115. 10.1016/j.chaos.2021.111115
28.
SohCWMahomedF. Linearization criteria for a system of second-order ordinary differential equations. J Int J Non-Linear Mech (2001) 36(4):671–7. 10.1016/s0020-7462(00)00032-9
29.
MaharajALeachP. The method of reduction of order and linearization of the two‐dimensional Ermakov system. J Math Methods Appl Sci (2007) 30(16):2125–45. 10.1002/mma.919
30.
AyubMKhanMMahomedF. Algebraic linearization criteria for systems of ordinary differential equations. J Nonlinear Dyn (2012) 67(3):2053–62. 10.1007/s11071-011-0128-x
31.
DuttHMSafdarMQadirA. Linearization criteria for two-dimensional systems of third-order ordinary differential equations by complex approach. J Arabian J Math (2019) 8(3):163–70. 10.1007/s40065-019-0238-8
32.
SongL-MYangZJLiXLZhangSM. Coherent superposition propagation of Laguerre–Gaussian and Hermite–Gaussian solitons. J Appl Math Lett (2020) 102:106114. 10.1016/j.aml.2019.106114
33.
ShenSYangZLiXZhangS. Periodic propagation of complex-valued hyperbolic-cosine-Gaussian solitons and breathers with complicated light field structure in strongly nonlocal nonlinear media. J Commun Nonlinear Sci Numer Simulation (2021) 103:106005. 10.1016/j.cnsns.2021.106005
34.
ShenSYangZJPangZGGeYR. The complex-valued astigmatic cosine-Gaussian soliton solution of the nonlocal nonlinear Schrödinger equation and its transmission characteristics. J Appl Math Lett (2022) 125:107755. 10.1016/j.aml.2021.107755
35.
StephaniH. Differential equations: Their solution using symmetries. Cambridge, United Kingdom: Cambridge University Press (1989).
36.
IbragimovNH. CRC handbook of Lie group analysis of differential equations. Florida, United States: CRC Press (1995).
37.
ZhangYZhangMQiS. Heat and mass transfer in a thin liquid film over an unsteady stretching surface in the presence of thermosolutal capillarity and variable magnetic field. Math Probl Eng (2016) 2016:8521580. 10.1155/2016/8521580
38.
WuYZhangXZhangX. Simplified analysis of heat and mass transfer model in droplet evaporation process. Appl Therm Eng (2016) 99:938–43. 10.1016/j.applthermaleng.2016.01.020
39.
SafdarMKhanMIKhanRATajSAbbasFElattarSet alAnalytic solutions for the MHD flow and heat transfer in a thin liquid film over an unsteady stretching surface with Lie symmetry and homotopy analysis method. Waves in Random Complex Media (2022) 2022:1–19. 10.1080/17455030.2022.2073402
Summary
Keywords
Lie similarity transformations, homotopy analysis method, symmetry, exact solutions, thin-film flow
Citation
Khan RA, Taj S, Ahmed S, Khan I and Eldin SM (2023) Lie symmetry and exact homotopic solutions of a non-linear double-diffusion problem. Front. Phys. 11:1150176. doi: 10.3389/fphy.2023.1150176
Received
23 January 2023
Accepted
08 February 2023
Published
28 February 2023
Volume
11 - 2023
Edited by
Samir A. El-Tantawy, Port Said University, Egypt
Reviewed by
Nadeem Sheikh, City University of Science and Information Technology, Pakistan
Muhammad Khan, University of Technology Malaysia, Malaysia
Updates
Copyright
© 2023 Khan, Taj, Ahmed, Khan and Eldin.
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: Ilyas Khan, i.said@mu.edu.sa
This article was submitted to Mathematical 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.