亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Wave solutions and numerical validation for the coupled reaction-advection-diffusion dynamical model in a porous medium

        2022-12-12 00:08:44AliMubarakiHwajoonKimNuruddeenUroojAkramandYasirAkbar
        Communications in Theoretical Physics 2022年12期

        Ali M Mubaraki ,Hwajoon Kim ,R I Nuruddeen ,Urooj Akram and Yasir Akbar

        1 Department of Mathematics and Statistics,Collage of Science,Taif University,P O Box 11099,Taif 21944,Saudi Arabia

        2 Kyungdong University,Yangju,Republic of Korea

        3 Department of Mathematics,Faculty of Science,Federal University Dutse,P O Box 7156 Dutse,Jigawa State,Nigeria

        4 Department of Mathematics,COMSATS University Islamabad,Lahore,Pakistan

        5 Department of Mathematics,COMSATS University Islamabad,44000 Islamabad,Pakistan

        Abstract The current study examines the special class of a generalized reaction-advection-diffusion dynamical model that is called the system of coupled Burger’s equations.This system plays a vital role in the essential areas of physics,including fluid dynamics and acoustics.Moreover,two promising analytical integration schemes are employed for the study;in addition to the deployment of an efficient variant of the eminent Adomian decomposition method.Three sets of analytical wave solutions are revealed,including exponential,periodic,and dark-singular wave solutions;while an amazed rapidly convergent approximate solution is acquired on the other hand.At the end,certain graphical illustrations and tables are provided to support the reported analytical and numerical results.No doubt,the present study is set to bridge the existing gap between the analytical and numerical approaches with regard to the solution validity of various models of mathematical physics.

        Keywords: reaction-advection-diffusion model,Burger’s equations,METEM,KM,NLDM

        1.Introduction

        Physical systems involving more than one process are greatly admired in the essential areas of science and engineering.One of these models that describes physical phenomena in a variety of processes is the generalized coupled system of reactionadvection-diffusion dynamical equations in porous medium that follows from the modification of [1] as (integer-order)

        wherew1=w1(x,t),andw2=w2(x,t),are the respective concentrations,α1,α2are the respective diffusion coefficients,γ1,γ2are the respective advection coefficients,θ1,θ2are the respective reaction coefficients,η1,λ1are real constants,m1,m2are non-negative integers,β1,β2,η2,λ2are coupling parameters;whilef1(x,t) andf2(x,t) are the respective source functions,otherwise known as inhomogeneous terms Additionally,whilem1andm2are non-negative integers.This model takes into consideration three important phenomena involving the motion of particles in a given medium,read [2] and the references therein.In particular,the first process,which is a reaction,describes possible processes like the adsorption,decay and reaction of the substances with other components.More so,the advection process refers to the transport of a substance or heat by the flow of a liquid;while diffusion is the movement of a substance from a region of high concentration to a region of low concentration throughout the physical domain of the problem’[3].Moreover,to highlight some universal or rather significant features of the governing model given in the above equation,several important coupled dynamical models may be extracted from this very model including,for instance,the reaction-advection model,reaction-diffusion model,advection-diffusion model,the system of Burger’s equation,and many other linear and nonlinear evolution equations to mention a few.

        Furthermore,having stated some universal nature of the governing model,and on the other hand,its application in a variety of areas of physics like fluid dynamics,environmental engineering,groundwater,and biological model;the present paper aims further to extract one of the celebrated models from the governing model that is called the system of coupled Burger’s equations [4].This model plays a vital role in modelling fluid flow problems,among others.It is pertinent to recall that the classical Burger’s equation plays an eminent part in modeling various phenomena in fluid dynamics,traffic flow,gas dynamics and nonlinear acoustics,to state a few.Additionally,different mathematical methods consisting of analytical and computations have in the past literature been employed to examine the coupled system of Burger’s equation.One may,however,read the works of[5,6]and the references therein with regard to various methods.Equally,one may find different extensions and modifications of this very important extract,including,for example,the coupled inhomogeneous singular Burger’s equations [7] among others.Besides,this model is an important coupled nonlinear evolution equation [8–11] that admits various forms of solutions.More so,it is important to recall that various integration schemes have been devised and successfully utilized on nonlinear evolution equations,see [12–17] and the references therewith.

        However,the current study examines the system of coupled Burger’s equations;being one of the astonishing members of the generalized reaction-advection-diffusion dynamical model with vast applications.In light of this,the study would employ two integration schemes,together with a modification of the famous Adomian decomposition method for the numerical study.More so,the modified extended tanh expansion method (METEM) [18–20],and the Kudryashov method(KM)[21–23]would be considered as the analytical methods;while a version of the Adomian’s decomposition method [24] that we call in this study the numerical Laplace decomposition method(NLDM) [25,26]will be adopted as the numerical approach.In addition,a number of plots and tables will be provided in the end,to portray the results to be acquired.Also,as for the limitations,some interesting cases the inhomogeneous variant of the coupled system of Burger’s equation would be analyzed in the present study (see section 3);besides,all the computational simulations and graphical illustrations would be implemented and produced on theWolfram Mathematica 9,sequentially.Lastly,the paper is organized as follows:section 2 gives the methods of solution.Section 3 presents the application of the methods given in section 2;while section 4 gives the application of the NLDM approach with examples.Section 5 discusses the obtained results;while section 6 gives some concluding comments.

        2.Methodology

        This section under consideration gives the outline of the adopted methodology of the present investigation.Basically,two modes of techniques will be employed,consisting of analytical and numerical approaches.For the analytical method,the modified extended tanh expansion method(METEM) [18–20] and the Kudryashov method (KM)[21–23] will be utilized;while the numerical Laplace decomposition method(NLDM)[25,26]will be adopted for the numerical purpose.The METEM and KM are elegant analytical methods that construct exact solutions to diverse evolution and Schrodinger equations,with great relevance in nonlinear sciences like optics and plasma.On the other hand,KM mainly gives exponential solution structures.For the applications of these solutions,kindly read section 5.Moreover,the two analytical methods give a series solution that depends onM(?N)to be computed based on the homogeneous balancing way—this way depends on the orders of the highest linear operator and the nonlinear term present in the governing equation.

        2.1.Analytical methods

        For the analytical methods,let us make consideration to the generalized nonlinear partial differential equation of the following pattern

        More so,upon making use of the following wave transformation

        in equation (2),one gets an ordinary differential equation of the form

        where in equation (3),kandcare real constants (nonzero);whileDin equation (4) is a differential operator.

        Therefore,we now make consideration to two of the aiming analytical methods as follows

        (i)METEM

        In accordance with METEM,the following finite series is offered [18–20]

        for equation (4),whereM(?N)is a natural number to be obtained by making a balance between the orders of the highest linear operator,and that of the nonlinear term present;whilea0,aj,bj,forj=1,2,…,Mare unknown constants that would be determined little later,ensuring thataM≠0,andbM≠0 for allM.

        Additionally,the application of METEM requires the function Ψ(ξ) in equation (5) to satisfy a Riccati differential equation of the form

        that admits the following solution sets

        Remarkably,METEM reveals three different types of solutions,including hyperbolic,periodic and irrational solutions.For the hyperbolic and periodic solutions,equation (7) demonstrated the modality of how these structures are constructed for any given evolution and Schrodinger equation;while we intentionally ignore the irrational solution possibilities atz=0.Complete sets of METEM solutions and applications can be found in[18–20] and the references therewith.

        (ii)KM

        According to KM,equation (4) is presumed to have a finite series solution as follows [21–23]

        whereM(?)Nis equally a natural number to be obtained by making a balance between the orders of the highest linear operator,and that of the nonlinear term present;whilea0,aj,forj=1,2,…,Mare unknown constants that would be determined,such thataM≠0,for allM.

        What’s more,KM requires the function Ψ(ξ) in equation(5)to satisfy the following differential equation

        such that the following exponential solution is satisfied by the latter equation is given by the function

        wheredis an arbitrary constant.This form of solution expressed in the above equation is basically the only form of solution posed by KM;in comparison with the METEM which gives three types of solutions.

        Hence,with the application of the METEM and KM as rightly presented above,the assumed solution in the respective cases given in equations (5) and (8) is thus substituted into equation (4) to yield a polynomial equation in Ψ(ξ).More so,upon collecting the coefficients of the resulting polynomial,and subsequently setting each coefficient to zero,one gets an algebraic system of equations.Of course,this tedious computation is done via the application of mathematical software.Therefore,we,state here that theWolfram Mathematica 9software will be utilized in the present study for all computational as well as graphical purposes.

        2.2.NLDM

        The current section uses the known Laplace integral transform in conjunction with the famous Adomian decomposition method [24] to present the NLDM procedure.In short,this section derives an iterative closed-form solution to a generalized nonlinear equation in mathematical physics.Thus,to present this approach,we consider the following one-dimensional nonlinear inhomogeneous partial differential equation

        with the following prescribed initial data

        Also,from equations (11)–(12),Lis a linear differential operator,Nis a nonlinear operator,h(x,t)is an inhomogeneous term;whileg(x) is the prescribed nice initial profile.

        Thus,as the celebrated Laplace integral transform will be utilized,let us define this integral transform and its corresponding inverse transform here as follows [27]

        Therefore,taking the Laplace transform of equation (11) intalongside using equation (12),the equation transforms to the following

        Next,we take the inverse Laplace transformL-1of the above equation ins,which gives

        More so,via the Adomian’s approach,we decompose the unknown solutionw(x,t) using the following infinite sum[24–26]

        and the nonlinear termed operatorN(w(x,t))using the followinginfniite su m of Adomian’s polynomials as follows

        whereRn?s are the polynomials devised by Adomian [24–26],and to be computed via the following recurrent formula

        So,by rewriting equation(15)in terms of the summations given in equations (16) and (17),one gets

        Lastly,identifying the terms arising from the prescribed initial data and the inhomogeneous function with the first componentw0(x,t),and the rest of the terms recurrently follow as suggested by the approach,we thus get the following recurrent scheme

        such that upon taking the net sum of them-term approximations yields

        where

        3.Application

        This section gives the application of the analytical and numerical approaches outlined in the above section.Additionally,the section makes consideration to yet another important nonlinear model that emanates from the generalized coupled reaction-advection-diffusion dynamical model earlier given in equation(1).Precisely,we make consideration of the coupled inhomogeneous Burger’s equation upon setting the respective diffusion coefficients,advection coefficients,and the reaction coefficients to zero,that is,α1=0=α2,γ1=0=γ2,and θ1=0=θ2.More so,the coupled inhomogeneous Burger’s equations of interest reads

        wherew1=w1(x,t),andw2=w2(x,t)are the respective wave profiles,η1,λ1are real constants,η2,λ2are coupling parameters;whilef1(x,t) andf2(x,t) are the respective source functions.Additionally,Burger’s equation has been used to describe a variety of phenomena,including a mathematical model of turbulence and an approximate theory of shock wave flow in viscous fluid.It is a simple model of sedimentation or evolution of scaled volume concentrations of two types of particles in fluid suspensions or colloids under gravity’s influence.This equation is interesting from a numerical viewpoint because analytical solutions are not generally available.Furthermore,we demonstrate the applicability of the presented analytical and numerical approaches on the governing coupled inhomogeneous Burger’s equation given above in what follows.

        3.1.Analytical methods

        To start off,let us consider the coupled homogeneous Burger’s equation from the above model,that isf1(x,t)=0=f(x,t).Therefore,we make use of the wave transformation for the coupled model as follows

        Therefore,the governing model in equation (23) becomes

        such that upon making a balance between the orders of the highest linear operator,and that of the nonlinear term present in the respective equation gives that

        Thus,we proceed as follows:

        (i)METEM

        WithM1=1,andM2=1,equation (25) via the application of the METEM admits the following solution

        wherea0,a1,b1,c0,c1andd1are unknown constants to be determined.Therefore,substituting the above solution form into equation (25) reveals the following algebraic system of equations,after setting each coefficient of Ψ(ξ) to zero as follows

        Now,upon solving the above algebraic system of equations,the unknown constantsa0,a1,b1,c0,c1andd1are determined as follows

        Figure 1.3D plots for the solution of coupled Burger’s equation given in equation (33) via METEM.

        Figure 2.3D plots for the solution of coupled Burger’s equation given in equation (35) via METEM.

        Therefore,fork=0.61,c=0.25,d1=1,andd1=-1,we give the three-dimensional (3D) plots for the solution given in equations (33) and (35),respectively,of the coupled homogeneous Burger’s equation via the application of the METEM in figures 1 and 2,sequentially.Specifically,the solutions given in equations(33)and (32) are dark-singular solutions;while those reported in equations(34)and(35)are singular periodic solutions.

        (ii)KM

        WithM1=1,andM2=1,equation (25) via the application of the KM admits the following solution

        wherea0,a1,b0andb1are unknown constants to be determined.Therefore,substituting the above solution form into equation (25) reveals the following algebraic system of equations,after setting each coefficient of Ψ(ξ) to zero as follows

        Then,on solving the above algebraic system of equations,the unknown constantsa0,a1,b0,andb1are determined as follows

        which yields only one wave solution set which is as follows

        Figure 3.3D plots for the solution of coupled Burger’s equation given in equation (39) via KM.

        wheredis an arbitrary constant.Thus,ford=0.1,k=1.5,andc=0.25,we give the 3D plots for the above solution of the coupled Burger’s equation via the application of the KM in figure 3.Again,the solution given in equation(39)through the KM is an exponential wave solution.As can be observed in figure 3,the plot is a kink-shaped profile.

        3.2.NLDM

        To successfully apply the NLDM,let us prescribe the following generalized initial conditions

        Therefore,without further delay,the coupled inhomogeneous Burger’s equation given in equation(23)admits the following recursive solution

        withAk?s,Bk?s andCk?s,respectively denote the polynomials by Adomian in favour of the following nonlinear terms

        where we express a few terms from these nonlinear terms using the application of equation (18) as follows

        and

        Lastly,upon taking the net sum ofm-term components,one obtains

        where

        4.Numerical validation

        In relation to the derived numerical scheme that is associated with the NLDM,this section further makes consideration to certain numerical examples that will give more light to the scheme.More so,we try to present some comparisons of results,where necessary.

        Example 4.1.Considering the coupled homogeneous Burger’s equation given in equation (23) with [1]

        together with the following prescribed initial data

        This problem further satisfies the following exact solution

        What’s more,the obtained recursive solution in equations(41)becomes

        where the Adomian polynomialsAk,BkandCkare determined from equation (42),accordingly.More so,we evaluate a few solution components from the above scheme as follows

        Therefore,taking the net sum of the components,one gets the following expressions

        that obviously converge to the following exact analytical solution

        Figure 4.3D plot for the approximate solution of coupled Burger’s equation given in Example 4.1.

        Figure 5.2D comparison between the exact and approximate solution of coupled Burger’s equation given in Example 4.1.

        Hence,we give the 3D and 2D plots for the obtained approximate solution of coupled homogeneous Burger’s equation in figures 4 and 5,respectively.Additionally,the absolute error difference between the exact solution(w1(x,t),w2(x,t))=(e-tsin (x),e-tsin(x))and that of the approximate solution via the sum of the first 15 componentsare reported in table 1.What’s more,a very good agreement between the exact solution and the reported numerical solution has been graphically shown in figure 5,just for using the sum of the first 15 components.In fact,this further proves the efficacy of the used method.

        Example 4.2.Considering the coupled homogeneous Burger’s equation given in equation (23) with [1]

        Table 1.Error difference between the exact solution of (w1(x,t),w2(x,t)) and the obtained numerical solutionat t=1.

        Table 1.Error difference between the exact solution of (w1(x,t),w2(x,t)) and the obtained numerical solutionat t=1.

        with the following initial data

        The problem further satisfies the following exact solution

        Accordingly,the recursive scheme in this particular example takes the following form

        where the Adomian polynomialsAk,BkandCkare equally obtained from equation (42).Further,we simulate this example and the absolute error difference between the exact solution (w1(x,t),w2(x,t)) and that of the approximate solution via the sum of the first 20 componentsare reported in table 2.Furthermore,we give the 3D and 2D plots for the obtained approximate solution in figures 6 and 7,respectively.What’s more,a very good agreement between the exact solution and the reported numerical solution has been graphically shown in figure 7,just after using the sum of the first 20 components.In fact,this further proves the efficacy of the used method.

        Example 4.3.Considering the coupled singular inhomogeneous Burger’s equation with [4,7]

        together with the following prescribed initial data

        Table 2.Error difference between the exact solution of (w1(x,t),w2(x,t)) and the obtained numerical solutionat t=1..

        Table 2.Error difference between the exact solution of (w1(x,t),w2(x,t)) and the obtained numerical solutionat t=1..

        This problem further satisfies the following exact solution

        Accordingly,the following recursive solution is thus obtained

        where the Adomian polynomialsAk,BkandCkare equally obtained from equation (42).More so,we evaluate a few solution components from the above scheme as follows

        Therefore,taking the net sum of the components,one gets the following expressions

        which is,in fact,the obvious available exact analytical solution.

        5.Discussion of results

        Figure 6.3D plots for the approximate solution of coupled Burger’s equation given in Example 4.2.

        Figure 7.2D comparisons between the exact and approximate solution of coupled Burger’s equation given in Example 4.2.

        The current study examines the system of coupled Burger’s equations[4–7];being one of the astonishing members of the generalized reaction-advection-diffusion dynamical model with vast applications.In doing so,two promising analytical integrations are employed;in addition to the deployment of an efficient variant of the eminent Adomian decomposition method.More specifically,the analytical methods involved are the METEM [18–20],and the KM [21–23];while the numerical method is a version of the Adomian’s method called NLDM [25,26].Three sets of analytical wave solutions are revealed including exponential,periodic,and darksingular wave solutions;while an amazed rapidly convergent approximate solution is acquired on the other hand via NLDM.Equally,self-explanatory graphical illustrations and tables are provided to support the reported analytical and numerical results.The reported exponential wave solution graphically turned out to be a kink-shaped profile.One would recall the application of kink solitons in polarization switches among dissimilar optical logic domains;while,a dark soliton solution is a localized wave that gives rise to a transitory decrease in wave amplitude (think also of the bright soliton solutions which are known to amplify the ocean waves).Also,periodic wave solutions are used in self-oscillatory and excitable systems,respectively,among others.More so,the numerical test examples considered reveal amazing closed form solutions,which get hold of their respective exact solutions when the net sums of the solution components are taken.In addition,the present study affirms the findings of the available literature.More comparatively,the present study generalizes many studies,including the integer-order derivative versions of the models considered in[4,6,7].Moreover,the employed numerical scheme which doubly serves as a semi-analytical method competes with so many methods in the literature,including [28–37];while some known robust integration schemes that happen to contend with the two analytical methods of interest in tackling different class of evolution equations include [38–43].

        6.Conclusion

        In conclusion,the current study examined the class of a reaction-advection-diffusion dynamical model that plays a vital role in essential areas of physics like fluid dynamics and acoustics.More specifically,we comprehensively studied the coupled system of Burger’s equations;being a special case of the above-mentioned governing model after setting the diffusion,advection,and reaction coefficients,respectively,to zero.Two known analytical integration schemes by the name METEM and KM were employed to scrutinize the model;in addition to the deployment of an efficient variant of the decomposition method,NLDM.Three sets of exact analytical wave solutions were revealed by the analytical methods,including exponential,periodic,and dark-singular wave solutions.Additionally,the derived numerical schemes revealed amazing approximate solutions that rapidly converged to the available exact solutions.Furthermore,certain graphical illustrations and tables were provided to support the presented analytical and numerical results;a perfect agreement between the exact solutions and the acquired numerical results,on the other hand,has been achieved with low number of iterations-see figures 5 and 7,for instance.We also state here that theWolfram Mathematica 9software was used for the simulation of numerical results,and their graphical depictions.No doubt,the present study is set to bridge the existing gap between the analytical and numerical approaches with regard to the solution validity of various models of mathematical physics;also,the employed methods can equally be used to tackle various higher-order dynamical systems of evolution equations.Moreover,in the future,we intended to make use of other reliable numerical approaches like the finite difference and finite element methods to further validate the exactitude of the reported exact analytical solutions in this study.

        ORCID iDs

        最新福利姬在线视频国产观看| 日韩精品久久久久久免费| 国产精品_国产精品_k频道w| 在线综合网| 国产精品福利视频一区| 1717国产精品久久| 亚洲产国偷v产偷v自拍色戒| 999国产精品视频| 国产性一交一乱一伦一色一情| 日韩精品有码在线视频| 亚洲码无人客一区二区三区| 男女深夜视频网站入口| 免费视频无打码一区二区三区| 特级精品毛片免费观看| 精品丝袜人妻久久久久久| 少妇的诱惑免费在线观看| 免费在线观看亚洲视频| 亚洲视频免费在线观看| 欧美精品videosex极品| 亚洲精品92内射| 国产精品女同一区二区| 亚洲精品白浆高清久久| 国产一区二区美女主播| 亚洲天堂亚洲天堂亚洲色图| 欧美午夜理伦三级在线观看| 日韩精品专区av无码| 一边吃奶一边摸做爽视频| 中文字幕av日韩精品一区二区| AV成人午夜无码一区二区| 亚洲啪啪AⅤ一区二区三区| 国产精品美女一区二区av| 成人自拍小视频在线看| 亚洲一区二区三区四区五区黄| 国产无套护士在线观看| 少妇高潮惨叫久久久久电影| 偷拍美女上厕所一区二区三区| 欧美大屁股xxxx高潮喷水| 国产伦久视频免费观看视频| 国产综合精品久久亚洲| 加勒比东京热久久综合| 亚洲中文字幕乱码在线观看|