The impact of memory effect on space fractional strong quantum couplers with tunable decay behavior and its numerical simulation

The nontrivial behavior of wave packets in the space fractional coupled nonlinear Schrödinger equation has received considerable theoretical attention. The difficulty comes from the fact that the Riesz fractional derivative is inherently a prehistorical operator. In contrast, nonlinear Schrödinger equation with both time and space nonlocal operators, which is the cornerstone in the modeling of a new type of fractional quantum couplers, is still in high demand of attention. This paper is devoted to numerically study the propagation of solitons through a new type of quantum couplers which can be called time-space fractional quantum couplers. The numerical methodology is based on the finite-difference/Galerkin Legendre spectral method with an easy to implement numerical algorithm. The time-fractional derivative is considered to describe the decay behavior and the nonlocal memory of the model. We conduct numerical simulations to observe the performance of the tunable decay and the sharpness behavior of the time-space fractional strongly coupled nonlinear Schrödinger model as well as the performance of the numerical algorithm. Numerical simulations show that the time and space fractional-order operators control the decay behavior or the memory and the sharpness of the interface and undergo a seamless transition of the fractional-order parameters

Effect of spatial configuration of an extended nonlinear Kierstead-Slobodkin reaction-transport model with adaptive numerical scheme

In this paper, we consider the numerical simulations of an extended nonlinear form of Kierstead-Slobodkin reaction-transport system in one and two dimensions. We employ the popular fourth-order exponential time differencing Runge-Kutta (ETDRK4) schemes proposed by Cox and Matthew (J Comput Phys 176:430-455, 2002), that was modified by Kassam and Trefethen (SIAM J Sci Comput 26:1214-1233, 2005), for the time integration of spatially discretized partial differential equations. We demonstrate the supremacy of ETDRK4 over the existing exponential time differencing integrators that are of standard approaches and provide timings and error comparison. Numerical results obtained in this paper have granted further insight to the question "What is the minimal size of the spatial domain so that the population persists?" posed by Kierstead and Slobodkin (J Mar Res 12:141-147, 1953 ), with a conclusive remark that the popula- tion size increases with the size of the domain. In attempt to examine the biological wave phenomena of the solutions, we present the numerical results in both one- and two-dimensional space, which have interesting ecological implications. Initial data and parameter values were chosen to mimic some existing patternsScopus 201

Exponential Integrator Methods for Nonlinear Fractional Reaction-diffusion Models

Nonlocality and spatial heterogeneity of many practical systems have made fractional differential equations very useful tools in Science and Engineering. However, solving these type of models is computationally demanding. In this work, we propose an exponential integrator method for nonlinear fractional reaction-diffusion equations. This scheme is based on using a real distinct poles discretization for the underlying matrix exponentials. Due to these real distinct poles, the algorithm could be easily implemented in parallel to take advantage of multiple processors for increased computational efficiency. The method is established to be second-order convergent; and proven to be robust for problems involving non-smooth/mismatched initial and boundary conditions and steep solution gradients. We examine the stability of the scheme through its amplification factor and plot the boundaries of the stability regions comparative to other second-order FETD schemes. This numerical scheme combined with fractional centered differencing is used for simulating many important nonlinear fractional models in applications. We demonstrate the superiority of our method over competing second order FETD schemes, BDF2 scheme, and IMEX schemes. Our experiments show that the proposed scheme is computationally more efficient (in terms of cpu time). Furthermore, we investigate the trade-off between using fractional centered differencing and matrix transfer technique in discretization of Riesz fractional derivatives. The generalized Mittag-Leffler function and its inverse is very useful in solving fractional differential equations and structural derivatives, respectively. However, their computational complexities have made them difficult to deal with numerically. We propose a real distinct pole rational approximation of the generalized Mittag-Leffler function. Under some mild conditions, this approximation is proven and empirically shown to be L-Acceptable. Due to the complete monotonicity property of the Mittag-Leffler function, we derive a rational approximation for the inverse generalized Mittag-Leffler function. These approximations are especially useful in developing efficient and accurate numerical schemes for partial differential equations of fractional order. Several applications are presented such as complementary error function, solution of fractional differential equations, and the ultraslow diffusion model using the structural derivative. Furthermore, we present a preliminary result of the application of the M-L RDP approximation to develop a generalized exponetial integrator scheme for time-fractional nonlinear reaction-diffusion equation

The impact of memory effect on space fractional strong quantum couplers with tunable decay behavior and its numerical simulation

The nontrivial behavior of wave packets in the space fractional coupled nonlinear Schrödinger equation has received considerable theoretical attention. The difficulty comes from the fact that the Riesz fractional derivative is inherently a prehistorical operator. In contrast, nonlinear Schrödinger equation with both time and space nonlocal operators, which is the cornerstone in the modeling of a new type of fractional quantum couplers, is still in high demand of attention. This paper is devoted to numerically study the propagation of solitons through a new type of quantum couplers which can be called time-space fractional quantum couplers. The numerical methodology is based on the finite-difference/Galerkin Legendre spectral method with an easy to implement numerical algorithm. The time-fractional derivative is considered to describe the decay behavior and the nonlocal memory of the model. We conduct numerical simulations to observe the performance of the tunable decay and the sharpness behavior of the time-space fractional strongly coupled nonlinear Schrödinger model as well as the performance of the numerical algorithm. Numerical simulations show that the time and space fractional-order operators control the decay behavior or the memory and the sharpness of the interface and undergo a seamless transition of the fractional-order parameters. © 2021, The Author(s).This study was supported financially by RFBR Grant (19-01-00019), the National Research Centre of Egypt (NRC) and Ghent university

Exponential integrators: tensor structured problems and applications

The solution of stiff systems of Ordinary Differential Equations (ODEs), that typically arise after spatial discretization of many important evolutionary Partial Differential Equations (PDEs), constitutes a topic of wide interest in numerical analysis. A prominent way to numerically integrate such systems involves using exponential integrators. In general, these kinds of schemes do not require the solution of (non)linear systems but rather the action of the matrix exponential and of some specific exponential-like functions (known in the literature as phi-functions). In this PhD thesis we aim at presenting efficient tensor-based tools to approximate such actions, both from a theoretical and from a practical point of view, when the problem has an underlying Kronecker sum structure. Moreover, we investigate the application of exponential integrators to compute numerical solutions of important equations in various fields, such as plasma physics, mean-field optimal control and computational chemistry. In any case, we provide several numerical examples and we perform extensive simulations, eventually exploiting modern hardware architectures such as multi-core Central Processing Units (CPUs) and Graphic Processing Units (GPUs). The results globally show the effectiveness and the superiority of the different approaches proposed

Trends, directions for further research, and some open problems of fractional calculus

The area of fractional calculus (FC) has been fast developing and is presently being applied in all scientific fields. Therefore, it is of key relevance to assess the present state of development and to foresee, if possible, the future evolution, or, at least, the challenges identified in the scope of advanced research works. This paper gives a vision about the directions for further research as well as some open problems of FC. A number of topics in mathematics, numerical algorithms and physics are analyzed, giving a systematic perspective for future research.info:eu-repo/semantics/publishedVersio

Error analysis of splitting methods for wave type equations

In dieser Doktorarbeit analysieren wir Splittingverfahren für zwei wellenartige Gleichungen. Wir untersuchen das Lie- und das Strangsplitting für die kubisch nichtlineare Schrödingergleichung auf dem Ganzraum und dem Torus in bis zu drei Raumdimensionen. Wir beweisen, dass das Strangsplitting in $L^2$ mit Ordnung $1+\theta$ für Anfangsfunktionen in $H^{2+2\theta}$ mit $\theta\in (0,1)$ konvergiert und dass beide Splittingverfahren mit Ordnung eins für Anfangsfunktionen in $H^2$ konvergieren. Wir bestätigen die theoretischen Konvergenzordnungen durch numerische Experimente. Außerdem analysieren wir ein "Alternating direction implicit"-Zeitsplittingverfahren für die Maxwellgleichungen mit Quellen, Strömen und Leitfähigkeit. Wir zeigen, dass es effizient ist, dass es mit Ordnung zwei in $L^2$ und in einem schwachen Sinne konvergiert, und dass es die Divergenzbedingungen bis auf Ordnung eins in $L^2$ und in einem schwachen Sinne erhält. Wir bestätigen die $L^2$-Resultate durch numerische Experimente

Special Functions: Fractional Calculus and the Pathway for Entropy

Historically, the notion of entropy emerged in conceptually very distinct contexts. This book deals with the connection between entropy, probability, and fractional dynamics as they appeared, for example, in solar neutrino astrophysics since the 1970's (Mathai and Rathie 1975, Mathai and Pederzoli 1977, Mathai and Saxena 1978, Mathai, Saxena, and Haubold 2010). The original solar neutrino problem, experimentally and theoretically, was resolved through the discovery of neutrino oscillations and was recently enriched by neutrino entanglement entropy. To reconsider possible new physics of solar neutrinos, diffusion entropy analysis, utilizing Boltzmann entropy, and standard deviation analysis was undertaken with Super-Kamiokande solar neutrino data. This analysis revealed a non-Gaussian signal with harmonic content. The Hurst exponent is different from the scaling exponent of the probability density function and both Hurst exponent and scaling exponent of the Super-Kamiokande data deviate considerably from the value of ½, which indicates that the statistics of the underlying phenomenon is anomalous. Here experiment may provide guidance about the generalization of theory of Boltzmann statistical mechanics. Arguments in the so-called Boltzmann-Planck-Einstein discussion related to Planck's discovery of the black-body radiation law are recapitulated mathematically and statistically and emphasize from this discussion is pursued that a meaningful implementation of the complex ‘entropy-probability-dynamics’ may offer two ways for explaining the results of diffusion entropy analysis and standard deviation analysis. One way is to consider an anomalous diffusion process that needs to use the fractional space-time diffusion equation (Gorenflo and Mainardi) and the other way is to consider a generalized Boltzmann entropy by assuming a power law probability density function. Here new mathematical framework, invented by sheer thought, may provide guidance for the generalization of Boltzmann statistical mechanics. In this book Boltzmann entropy, generalized by Tsallis and Mathai, is considered. The second one contains a varying parameter that is used to construct an entropic pathway covering generalized type-1 beta, type-2 beta, and gamma families of densities. Similarly, pathways for respective distributions and differential equations can be developed. Mathai's entropy is optimized under various conditions reproducing the well-known Boltzmann distribution, Raleigh distribution, and other distributions used in physics. Properties of the entropy measure for the generalized entropy are examined. In this process the role of special functions of mathematical physics, particularly the H-function, is highlighted