A Recent paper by Ma et al.,claims to estimate the state of charge of Lithium-ion batteries with a fractionalorder impedance model including a Warburg and a constant phase element(CPE)with a maximum error of 0.5%[1].T...A Recent paper by Ma et al.,claims to estimate the state of charge of Lithium-ion batteries with a fractionalorder impedance model including a Warburg and a constant phase element(CPE)with a maximum error of 0.5%[1].The proposed equivalent circuit model from[1]is reproduced in Fig.1.展开更多
Recently, a conformable fractional derivative has been proposed to calculate the derivative of non-integer order of time functions. It has been shown that this new fractional derivative definition obeys many advantage...Recently, a conformable fractional derivative has been proposed to calculate the derivative of non-integer order of time functions. It has been shown that this new fractional derivative definition obeys many advantages over the preceding definitions. For mathematical models in applied sciences and to preserve the dimensionality of the physical quantities, an auxiliary parameter (~r) which has the dimension of seconds should be implemented in the fractional derivative definition. We obtain analytic solutions for the resulting conformable fractional differential equations describing the vertical velocity and the height of the falling body. It is shown that the dimensions of velocity and height are always correct without any restrictions on the auxiliary parameter cr which contradicts with the results in the literature when applying the Caputo definition to the same problem. This may open the door for many future works either to describe the role of such an auxiliary parameter or to derive a more suitable definition for the fractional derivative.展开更多
Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensiona...Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.展开更多
We propose a theoretical study investigating the spread of the novel coronavirus(COVID-19)reported inWuhan City of China in 2019.We develop a mathematical model based on the novel corona virus’s characteristics and t...We propose a theoretical study investigating the spread of the novel coronavirus(COVID-19)reported inWuhan City of China in 2019.We develop a mathematical model based on the novel corona virus’s characteristics and then use fractional calculus to fractionalize it.Various fractional order epidemicmodels have been formulated and analyzed using a number of iterative and numerical approacheswhile the complications arise due to singular kernel.We use the well-known Caputo-Fabrizio operator for the purposes of fictionalization because this operator is based on the non-singular kernel.Moreover,to analyze the existence and uniqueness,we will use the well-known fixed point theory.We also prove that the considered model has positive and bounded solutions.We also draw some numerical simulations to verify the theoretical work via graphical representations.We believe that the proposed epidemic model will be helpful for health officials to take some positive steps to control contagious diseases.展开更多
The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced.The flow near a wall suddenly set in mo- tion is studied for a non-Newtonian viscoelastic fluid with the frac...The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced.The flow near a wall suddenly set in mo- tion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model.Exact solutions of velocity and stress are obtained by using the discrete in- verse Laplace transform of the sequential fractional derivatives.It is found that the effect of the fractional orders in the constitutive relationship on the flow field is signif- icant.The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate,for large times the viscoelastic effects become weak.展开更多
The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases ...The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.展开更多
State of charge(SOC)estimation for lithium ion batteries plays a critical role in battery management systems for electric vehicles.Battery fractional order models(FOMs)which come from frequency-domain modelling have p...State of charge(SOC)estimation for lithium ion batteries plays a critical role in battery management systems for electric vehicles.Battery fractional order models(FOMs)which come from frequency-domain modelling have provided a distinct insight into SOC estimation.In this article,we compare five state-of-the-art FOMs in terms of SOC estimation.To this end,firstly,characterisation tests on lithium ion batteries are conducted,and the experimental results are used to identify FOM parameters.Parameter identification results show that increasing the complexity of FOMs cannot always improve accuracy.The model R(RQ)W shows superior identification accuracy than the other four FOMs.Secondly,the SOC estimation based on a fractional order unscented Kalman filter is conducted to compare model accuracy and computational burden under different profiles,memory lengths,ambient temperatures,cells and voltage/current drifts.The evaluation results reveal that the SOC estimation accuracy does not necessarily positively correlate to the complexity of FOMs.Although more complex models can have better robustness against temperature variation,R(RQ),the simplest FOM,can overall provide satisfactory accuracy.Validation results on different cells demonstrate the generalisation ability of FOMs,and R(RQ)outperforms other models.Moreover,R(RQ)shows better robustness against truncation error and can maintain high accuracy even under the occurrence of current or voltage sensor drift.展开更多
To fully display the modeling mechanism of the novelfractional order grey model (FGM (q,1)), this paper decomposesthe data matrix of the model into the mean generation matrix, theaccumulative generation matrix and...To fully display the modeling mechanism of the novelfractional order grey model (FGM (q,1)), this paper decomposesthe data matrix of the model into the mean generation matrix, theaccumulative generation matrix and the raw data matrix, whichare consistent with the fractional order accumulative grey model(FAGM (1,1)). Following this, this paper decomposes the accumulativedata difference matrix into the accumulative generationmatrix, the q-order reductive accumulative matrix and the rawdata matrix, and then combines the least square method, findingthat the differential order affects the model parameters only byaffecting the formation of differential sequences. This paper thensummarizes matrix decomposition of some special sequences,such as the sequence generated by the strengthening and weakeningoperators, the jumping sequence, and the non-equidistancesequence. Finally, this paper expresses the influences of the rawdata transformation, the accumulation sequence transformation,and the differential matrix transformation on the model parametersas matrices, and takes the non-equidistance sequence as an exampleto show the modeling mechanism.展开更多
In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-d...In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.展开更多
Because charge carriers of many organic semiconductors(OSCs)exhibit fractional drift diffusion(Fr-DD)transport properties,the need to develop a Fr-DD model solver becomes more apparent.However,the current research on ...Because charge carriers of many organic semiconductors(OSCs)exhibit fractional drift diffusion(Fr-DD)transport properties,the need to develop a Fr-DD model solver becomes more apparent.However,the current research on solving the governing equations of the Fr-DD model is practically nonexistent.In this paper,an iterative solver with high precision is developed to solve both the transient and steady-state Fr-DD model for organic semiconductor devices.The Fr-DD model is composed of two fractionalorder carriers(i.e.,electrons and holes)continuity equations coupled with Poisson’s equation.By treating the current density as constants within each pair of consecutive grid nodes,a linear Caputo’s fractional-order ordinary differential equation(FrODE)can be produced,and its analytic solution gives an approximation to the carrier concentration.The convergence of the solver is guaranteed by implementing a successive over-relaxation(SOR)mechanism on each loop of Gummel’s iteration.Based on our derivations,it can be shown that the Scharfetter–Gummel discretization method is essentially a special case of our scheme.In addition,the consistency and convergence of the two core algorithms are proved,with three numerical examples designed to demonstrate the accuracy and computational performance of this solver.Finally,we validate the Fr-DD model for a steady-state organic field effect transistor(OFET)by fitting the simulated transconductance and output curves to the experimental data.展开更多
The COVID-19 pandemic is a curse and a threat to global health, development, the economy, and peaceful society because of its massive transmission and high rates of mutation. More than 220 countries have been affected...The COVID-19 pandemic is a curse and a threat to global health, development, the economy, and peaceful society because of its massive transmission and high rates of mutation. More than 220 countries have been affected by COVID-19. The world is now facing a drastic situation because of this ongoing virus. Bangladesh is also dealing with this issue, and due to its dense population, it is particularly vulnerable to the spread of COVID-19. Recently, many non-linear systems have been proposed to solve the SIR (Susceptible, Infected, and Recovered) model for predicting Coronavirus cases. In this paper, we have discussed the fractional order SIR epidemic model of a non-fatal disease in a population of a constant size. Using the Laplace Adomian Decomposition method, we get an approximate solution to the model. To predict the dynamic transmission of COVID-19 in Bangladesh, we provide a numerical argument based on real data. We also conducted a comparative analysis among susceptible, infected, and recovered people. Furthermore, the most sensitive parameters for the Basic Reproduction Number (<em>R</em><sub>0</sub>) are graphically presented, and the impact of the compartments on the transmission dynamics of the COVID-19 pandemic is thoroughly investigated.展开更多
In this paper,the fractional natural decomposition method(FNDM)is employed to find the solution for the Kundu-Eckhaus equation and coupled fractional differential equations describing the massive Thirring model.Themas...In this paper,the fractional natural decomposition method(FNDM)is employed to find the solution for the Kundu-Eckhaus equation and coupled fractional differential equations describing the massive Thirring model.Themassive Thirring model consists of a system of two nonlinear complex differential equations,and it plays a dynamic role in quantum field theory.The fractional derivative is considered in the Caputo sense,and the projected algorithm is a graceful mixture of Adomian decomposition scheme with natural transform technique.In order to illustrate and validate the efficiency of the future technique,we analyzed projected phenomena in terms of fractional order.Moreover,the behaviour of the obtained solution has been captured for diverse fractional order.The obtained results elucidate that the projected technique is easy to implement and very effective to analyze the behaviour of complex nonlinear differential equations of fractional order arising in the connected areas of science and engineering.展开更多
We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gor...We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gordon equation expansion(ESh GEE) schemes are utilized. The solutions obtained include dark, bright, dark-bright, periodic and other kinds of solitons. These analytical wave solutions are gained and verified with the use of Mathematica software. These solutions do not exist in literature. Some of the solutions are demonstrated by 2D, 3D and contour graphs. This model is mostly used in circuit theory, transmission of nerve impulses, and population genetics. Finally, both the schemes are more applicable, reliable and significant to deal with the fractional nonlinear partial differential equations.展开更多
In this paper, we provide a new approach to solve approximately a system of fractional differential equations (FDEs). We extend this approach for approximately solving a fractional-order differential equation model of...In this paper, we provide a new approach to solve approximately a system of fractional differential equations (FDEs). We extend this approach for approximately solving a fractional-order differential equation model of HIV infection of CD4<sup>+</sup>T cells with therapy effect. The fractional derivative in our approach is in the sense of Riemann-Liouville. To solve the problem, we reduce the system of FDE to a discrete optimization problem. By obtaining the optimal solutions of new problem by minimization the total errors, we obtain the approximate solution of the original problem. The numerical solutions obtained from the proposed approach indicate that our approximation is easy to implement and accurate when it is applied to a systems of FDEs.展开更多
The work is devoted to the fractional characterization of time-dependent coupled convection-diffusion systems arising in magnetohydrodynamics(MHD)flows.The time derivative is expressed by means of Caputo’s fractional...The work is devoted to the fractional characterization of time-dependent coupled convection-diffusion systems arising in magnetohydrodynamics(MHD)flows.The time derivative is expressed by means of Caputo’s fractional derivative concept,while the model is solved via the full-spectral method(FSM)and the semi-spectral scheme(SSS).The FSM is based on the operational matrices of derivatives constructed by using higher-order orthogonal polynomials and collocation techniques.The SSS is developed by discretizing the time variable,and the space domain is collocated by using equal points.A detailed comparative analysis is made through graphs for various parameters and tables with existing literature.The contour graphs are made to show the behaviors of the velocity and magnetic fields.The proposed methods are reasonably efficient in examining the behavior of convection-diffusion equations arising in MHD flows,and the concept may be extended for variable order models arising in MHD flows.展开更多
This paper presents an original theoretical framework to model steel material properties in continuous casting line process. Specific properties arising from non-Newtonian dynamics are herein used to indicate the natu...This paper presents an original theoretical framework to model steel material properties in continuous casting line process. Specific properties arising from non-Newtonian dynamics are herein used to indicate the natural convergence of distributed parameter systems to fractional order transfer function models. Data driven identification from a real continuous casting line is used to identify model of the electromagnetic actuator device to control flow velocity of liquid steel. To ensure product specifications, a fractional order control is designed and validated on the system. A projection of the closed loop performance onto the quality assessment at end production line is also given in this paper.展开更多
The motion of pore water directly influences mechanical properties of soils, which are variable during creep. Accurate description of the evolution of mechanical properties of soils can help to reveal the internal beh...The motion of pore water directly influences mechanical properties of soils, which are variable during creep. Accurate description of the evolution of mechanical properties of soils can help to reveal the internal behavior of pore water. Based on the idea of using the fractional order to reflect mechanical properties of soils, a fractional creep model is proposed by introducing a variable-order fractional operator, and realized on a series of creep responses in soft soils. A comparative analysis illustrates that the evolution of mechanical properties, shown through the simulated results, exactly corresponds to the motion of pore water and the solid skeleton. This demonstrates that the proposed variable-order fractional model can be employed to characterize the evolution of mechanical properties of and the pore water motion in soft soils during creep. It is observed that the fractional order from the proposed model is related to the dissipation rate of pore water pressure.展开更多
In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a h...In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.展开更多
Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operation...Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.展开更多
Understanding the variations in microscopic pore-fracture structures(MPFS) during coal creep under pore pressure and stress coupling is crucial for coal mining and effective gas treatment. In this manuscript, a triaxi...Understanding the variations in microscopic pore-fracture structures(MPFS) during coal creep under pore pressure and stress coupling is crucial for coal mining and effective gas treatment. In this manuscript, a triaxial creep test on deep coal at various pore pressures using a test system that combines in-situ mechanical loading with real-time nuclear magnetic resonance(NMR) detection was conducted.Full-scale quantitative characterization, online real-time detection, and visualization of MPFS during coal creep influenced by pore pressure and stress coupling were performed using NMR and NMR imaging(NMRI) techniques. The results revealed that seepage pores and microfractures(SPM) undergo the most significant changes during coal creep, with creep failure gradually expanding from dense primary pore fractures. Pore pressure presence promotes MPFS development primarily by inhibiting SPM compression and encouraging adsorption pores(AP) to evolve into SPM. Coal enters the accelerated creep stage earlier at lower stress levels, resulting in more pronounced creep deformation. The connection between the micro and macro values was established, demonstrating that increased porosity at different pore pressures leads to a negative exponential decay of the viscosity coefficient. The Newton dashpot in the ideal viscoplastic body and the Burgers model was improved using NMR experimental results, and a creep model that considers pore pressure and stress coupling using variable-order fractional operators was developed. The model’s reasonableness was confirmed using creep experimental data. The damagestate adjustment factors ω and β were identified through a parameter sensitivity analysis to characterize the effect of pore pressure and stress coupling on the creep damage characteristics(size and degree of difficulty) of coal.展开更多
文摘A Recent paper by Ma et al.,claims to estimate the state of charge of Lithium-ion batteries with a fractionalorder impedance model including a Warburg and a constant phase element(CPE)with a maximum error of 0.5%[1].The proposed equivalent circuit model from[1]is reproduced in Fig.1.
文摘Recently, a conformable fractional derivative has been proposed to calculate the derivative of non-integer order of time functions. It has been shown that this new fractional derivative definition obeys many advantages over the preceding definitions. For mathematical models in applied sciences and to preserve the dimensionality of the physical quantities, an auxiliary parameter (~r) which has the dimension of seconds should be implemented in the fractional derivative definition. We obtain analytic solutions for the resulting conformable fractional differential equations describing the vertical velocity and the height of the falling body. It is shown that the dimensions of velocity and height are always correct without any restrictions on the auxiliary parameter cr which contradicts with the results in the literature when applying the Caputo definition to the same problem. This may open the door for many future works either to describe the role of such an auxiliary parameter or to derive a more suitable definition for the fractional derivative.
基金the financial support from the National Natural Science Foundation of China(No.52109119)the Guangxi Natural Science Foundation(No.2021GXNSFBA075030)+2 种基金the Guangxi Science and Technology Project(No.Guike AD20325002)the Chinese Postdoctoral Science Fund Project(No.2022 M723408)the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin(China Institute of Water Resources and Hydropower Research)(No.IWHR-SKL-202202).
文摘Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.
基金supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project No. (PNURSP2022R14),Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia。
文摘We propose a theoretical study investigating the spread of the novel coronavirus(COVID-19)reported inWuhan City of China in 2019.We develop a mathematical model based on the novel corona virus’s characteristics and then use fractional calculus to fractionalize it.Various fractional order epidemicmodels have been formulated and analyzed using a number of iterative and numerical approacheswhile the complications arise due to singular kernel.We use the well-known Caputo-Fabrizio operator for the purposes of fictionalization because this operator is based on the non-singular kernel.Moreover,to analyze the existence and uniqueness,we will use the well-known fixed point theory.We also prove that the considered model has positive and bounded solutions.We also draw some numerical simulations to verify the theoretical work via graphical representations.We believe that the proposed epidemic model will be helpful for health officials to take some positive steps to control contagious diseases.
基金The project supported by the National Natural Science Foundation of China (10002003)Foundation for University Key Teacher by the Ministry of EducationResearch Fund for the Doctoral Program of Higher Education
文摘The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced.The flow near a wall suddenly set in mo- tion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model.Exact solutions of velocity and stress are obtained by using the discrete in- verse Laplace transform of the sequential fractional derivatives.It is found that the effect of the fractional orders in the constitutive relationship on the flow field is signif- icant.The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate,for large times the viscoelastic effects become weak.
基金The project supported by the National Natural Science Foundation of China (10272067, 10426024)the Doctoral Program Foundation of the Education Ministry of China (20030422046)the Natural Science Foundation of Shandong University at Weihai.
文摘The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.
基金Beijing Municipal Natural Science Foundation of China(Grant No.3182035)National Natural Science Foundation of China(Grant No.51877009).
文摘State of charge(SOC)estimation for lithium ion batteries plays a critical role in battery management systems for electric vehicles.Battery fractional order models(FOMs)which come from frequency-domain modelling have provided a distinct insight into SOC estimation.In this article,we compare five state-of-the-art FOMs in terms of SOC estimation.To this end,firstly,characterisation tests on lithium ion batteries are conducted,and the experimental results are used to identify FOM parameters.Parameter identification results show that increasing the complexity of FOMs cannot always improve accuracy.The model R(RQ)W shows superior identification accuracy than the other four FOMs.Secondly,the SOC estimation based on a fractional order unscented Kalman filter is conducted to compare model accuracy and computational burden under different profiles,memory lengths,ambient temperatures,cells and voltage/current drifts.The evaluation results reveal that the SOC estimation accuracy does not necessarily positively correlate to the complexity of FOMs.Although more complex models can have better robustness against temperature variation,R(RQ),the simplest FOM,can overall provide satisfactory accuracy.Validation results on different cells demonstrate the generalisation ability of FOMs,and R(RQ)outperforms other models.Moreover,R(RQ)shows better robustness against truncation error and can maintain high accuracy even under the occurrence of current or voltage sensor drift.
基金supported by the National Natural Science Foundation of China(5147915151279149+2 种基金71540027)the China Postdoctoral Science Foundation Special Foundation Project(2013T607552012M521487)
文摘To fully display the modeling mechanism of the novelfractional order grey model (FGM (q,1)), this paper decomposesthe data matrix of the model into the mean generation matrix, theaccumulative generation matrix and the raw data matrix, whichare consistent with the fractional order accumulative grey model(FAGM (1,1)). Following this, this paper decomposes the accumulativedata difference matrix into the accumulative generationmatrix, the q-order reductive accumulative matrix and the rawdata matrix, and then combines the least square method, findingthat the differential order affects the model parameters only byaffecting the formation of differential sequences. This paper thensummarizes matrix decomposition of some special sequences,such as the sequence generated by the strengthening and weakeningoperators, the jumping sequence, and the non-equidistancesequence. Finally, this paper expresses the influences of the rawdata transformation, the accumulation sequence transformation,and the differential matrix transformation on the model parametersas matrices, and takes the non-equidistance sequence as an exampleto show the modeling mechanism.
基金the National Natural Science Fund(11661058,11761053)Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2017MS0107)+1 种基金Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07)National Undergraduate Innovative Training Project of Inner Mongolia University(201710126026).
文摘In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.
基金This work was supported in part by the National Science Foundation through Grant CNS-1726865by the USDA under Grant 2019-67021-28990.
文摘Because charge carriers of many organic semiconductors(OSCs)exhibit fractional drift diffusion(Fr-DD)transport properties,the need to develop a Fr-DD model solver becomes more apparent.However,the current research on solving the governing equations of the Fr-DD model is practically nonexistent.In this paper,an iterative solver with high precision is developed to solve both the transient and steady-state Fr-DD model for organic semiconductor devices.The Fr-DD model is composed of two fractionalorder carriers(i.e.,electrons and holes)continuity equations coupled with Poisson’s equation.By treating the current density as constants within each pair of consecutive grid nodes,a linear Caputo’s fractional-order ordinary differential equation(FrODE)can be produced,and its analytic solution gives an approximation to the carrier concentration.The convergence of the solver is guaranteed by implementing a successive over-relaxation(SOR)mechanism on each loop of Gummel’s iteration.Based on our derivations,it can be shown that the Scharfetter–Gummel discretization method is essentially a special case of our scheme.In addition,the consistency and convergence of the two core algorithms are proved,with three numerical examples designed to demonstrate the accuracy and computational performance of this solver.Finally,we validate the Fr-DD model for a steady-state organic field effect transistor(OFET)by fitting the simulated transconductance and output curves to the experimental data.
文摘The COVID-19 pandemic is a curse and a threat to global health, development, the economy, and peaceful society because of its massive transmission and high rates of mutation. More than 220 countries have been affected by COVID-19. The world is now facing a drastic situation because of this ongoing virus. Bangladesh is also dealing with this issue, and due to its dense population, it is particularly vulnerable to the spread of COVID-19. Recently, many non-linear systems have been proposed to solve the SIR (Susceptible, Infected, and Recovered) model for predicting Coronavirus cases. In this paper, we have discussed the fractional order SIR epidemic model of a non-fatal disease in a population of a constant size. Using the Laplace Adomian Decomposition method, we get an approximate solution to the model. To predict the dynamic transmission of COVID-19 in Bangladesh, we provide a numerical argument based on real data. We also conducted a comparative analysis among susceptible, infected, and recovered people. Furthermore, the most sensitive parameters for the Basic Reproduction Number (<em>R</em><sub>0</sub>) are graphically presented, and the impact of the compartments on the transmission dynamics of the COVID-19 pandemic is thoroughly investigated.
文摘In this paper,the fractional natural decomposition method(FNDM)is employed to find the solution for the Kundu-Eckhaus equation and coupled fractional differential equations describing the massive Thirring model.Themassive Thirring model consists of a system of two nonlinear complex differential equations,and it plays a dynamic role in quantum field theory.The fractional derivative is considered in the Caputo sense,and the projected algorithm is a graceful mixture of Adomian decomposition scheme with natural transform technique.In order to illustrate and validate the efficiency of the future technique,we analyzed projected phenomena in terms of fractional order.Moreover,the behaviour of the obtained solution has been captured for diverse fractional order.The obtained results elucidate that the projected technique is easy to implement and very effective to analyze the behaviour of complex nonlinear differential equations of fractional order arising in the connected areas of science and engineering.
文摘We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gordon equation expansion(ESh GEE) schemes are utilized. The solutions obtained include dark, bright, dark-bright, periodic and other kinds of solitons. These analytical wave solutions are gained and verified with the use of Mathematica software. These solutions do not exist in literature. Some of the solutions are demonstrated by 2D, 3D and contour graphs. This model is mostly used in circuit theory, transmission of nerve impulses, and population genetics. Finally, both the schemes are more applicable, reliable and significant to deal with the fractional nonlinear partial differential equations.
文摘In this paper, we provide a new approach to solve approximately a system of fractional differential equations (FDEs). We extend this approach for approximately solving a fractional-order differential equation model of HIV infection of CD4<sup>+</sup>T cells with therapy effect. The fractional derivative in our approach is in the sense of Riemann-Liouville. To solve the problem, we reduce the system of FDE to a discrete optimization problem. By obtaining the optimal solutions of new problem by minimization the total errors, we obtain the approximate solution of the original problem. The numerical solutions obtained from the proposed approach indicate that our approximation is easy to implement and accurate when it is applied to a systems of FDEs.
基金Project supported by the National Natural Science Foundation of China(Nos.12250410244,11872151)the Jiangsu Province Education Development Special Project-2022 for Double First-ClassSchool Talent Start-up Fund of China(No.2022r109)the Longshan Scholar Program of Jiangsu Province of China。
文摘The work is devoted to the fractional characterization of time-dependent coupled convection-diffusion systems arising in magnetohydrodynamics(MHD)flows.The time derivative is expressed by means of Caputo’s fractional derivative concept,while the model is solved via the full-spectral method(FSM)and the semi-spectral scheme(SSS).The FSM is based on the operational matrices of derivatives constructed by using higher-order orthogonal polynomials and collocation techniques.The SSS is developed by discretizing the time variable,and the space domain is collocated by using equal points.A detailed comparative analysis is made through graphs for various parameters and tables with existing literature.The contour graphs are made to show the behaviors of the velocity and magnetic fields.The proposed methods are reasonably efficient in examining the behavior of convection-diffusion equations arising in MHD flows,and the concept may be extended for variable order models arising in MHD flows.
基金supported by Research Foundation Flanders(FWO)(1S04719N,12X6819N)partially supported by a grant of the Ministry of Research+2 种基金Innovation and DigitizationCNCS-UEFISCDIproject number PN-Ⅲ-P1-1.1-PD-2021-0204,within PNCDIⅢ。
文摘This paper presents an original theoretical framework to model steel material properties in continuous casting line process. Specific properties arising from non-Newtonian dynamics are herein used to indicate the natural convergence of distributed parameter systems to fractional order transfer function models. Data driven identification from a real continuous casting line is used to identify model of the electromagnetic actuator device to control flow velocity of liquid steel. To ensure product specifications, a fractional order control is designed and validated on the system. A projection of the closed loop performance onto the quality assessment at end production line is also given in this paper.
基金supported by the Natural Science Foundation of Jiangsu Province of China(Grant No.BK2012810)the Fundamental Research Funds for the Central Universities(Grant No.2009B15114)
文摘The motion of pore water directly influences mechanical properties of soils, which are variable during creep. Accurate description of the evolution of mechanical properties of soils can help to reveal the internal behavior of pore water. Based on the idea of using the fractional order to reflect mechanical properties of soils, a fractional creep model is proposed by introducing a variable-order fractional operator, and realized on a series of creep responses in soft soils. A comparative analysis illustrates that the evolution of mechanical properties, shown through the simulated results, exactly corresponds to the motion of pore water and the solid skeleton. This demonstrates that the proposed variable-order fractional model can be employed to characterize the evolution of mechanical properties of and the pore water motion in soft soils during creep. It is observed that the fractional order from the proposed model is related to the dissipation rate of pore water pressure.
基金supported by the Project of the National Key R&D Program(Grant No.2021YFA1000202)National Natural Science Foundation of China(Grant Nos.12120101001,12001326 and 12171283)+2 种基金Natural Science Foundation of Shandong Province(Grant Nos.ZR2021ZD03,ZR2020QA032 and ZR2019ZD42)China Postdoctoral Science Foundation(Grant Nos.BX20190191 and 2020M672038)the Startup Fund from Shandong University(Grant No.11140082063130)。
文摘In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.
基金funded by the joint fund of the National Key Research and Development Program of China(No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金Open Foundation of National Energy shale gas R&D(experiment) center(2022-KFKT-12)the 111 Project(B17009)。
文摘Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.
基金supported by the National Natural Science Foundation of China(Nos.52121003,51827901 and 52204110)China Postdoctoral Science Foundation(No.2022M722346)+1 种基金the 111 Project(No.B14006)the Yueqi Outstanding Scholar Program of CUMTB(No.2017A03).
文摘Understanding the variations in microscopic pore-fracture structures(MPFS) during coal creep under pore pressure and stress coupling is crucial for coal mining and effective gas treatment. In this manuscript, a triaxial creep test on deep coal at various pore pressures using a test system that combines in-situ mechanical loading with real-time nuclear magnetic resonance(NMR) detection was conducted.Full-scale quantitative characterization, online real-time detection, and visualization of MPFS during coal creep influenced by pore pressure and stress coupling were performed using NMR and NMR imaging(NMRI) techniques. The results revealed that seepage pores and microfractures(SPM) undergo the most significant changes during coal creep, with creep failure gradually expanding from dense primary pore fractures. Pore pressure presence promotes MPFS development primarily by inhibiting SPM compression and encouraging adsorption pores(AP) to evolve into SPM. Coal enters the accelerated creep stage earlier at lower stress levels, resulting in more pronounced creep deformation. The connection between the micro and macro values was established, demonstrating that increased porosity at different pore pressures leads to a negative exponential decay of the viscosity coefficient. The Newton dashpot in the ideal viscoplastic body and the Burgers model was improved using NMR experimental results, and a creep model that considers pore pressure and stress coupling using variable-order fractional operators was developed. The model’s reasonableness was confirmed using creep experimental data. The damagestate adjustment factors ω and β were identified through a parameter sensitivity analysis to characterize the effect of pore pressure and stress coupling on the creep damage characteristics(size and degree of difficulty) of coal.