Generalized continued fractions and ergodic theory

NASA Astrophysics Data System (ADS)

Pustyl'nikov, L. D.

2003-02-01

In this paper a new theory of generalized continued fractions is constructed and applied to numbers, multidimensional vectors belonging to a real space, and infinite-dimensional vectors with integral coordinates. The theory is based on a concept generalizing the procedure for constructing the classical continued fractions and substantially using ergodic theory. One of the versions of the theory is related to differential equations. In the finite-dimensional case the constructions thus introduced are used to solve problems posed by Weyl in analysis and number theory concerning estimates of trigonometric sums and of the remainder in the distribution law for the fractional parts of the values of a polynomial, and also the problem of characterizing algebraic and transcendental numbers with the use of generalized continued fractions. Infinite-dimensional generalized continued fractions are applied to estimate sums of Legendre symbols and to obtain new results in the classical problem of the distribution of quadratic residues and non-residues modulo a prime. In the course of constructing these continued fractions, an investigation is carried out of the ergodic properties of a class of infinite-dimensional dynamical systems which are also of independent interest.

A structure-preserving method for a class of nonlinear dissipative wave equations with Riesz space-fractional derivatives

NASA Astrophysics Data System (ADS)

Macías-Díaz, J. E.

2017-12-01

In this manuscript, we consider an initial-boundary-value problem governed by a (1 + 1)-dimensional hyperbolic partial differential equation with constant damping that generalizes many nonlinear wave equations from mathematical physics. The model considers the presence of a spatial Laplacian of fractional order which is defined in terms of Riesz fractional derivatives, as well as the inclusion of a generic continuously differentiable potential. It is known that the undamped regime has an associated positive energy functional, and we show here that it is preserved throughout time under suitable boundary conditions. To approximate the solutions of this model, we propose a finite-difference discretization based on fractional centered differences. Some discrete quantities are proposed in this work to estimate the energy functional, and we show that the numerical method is capable of conserving the discrete energy under the same boundary conditions for which the continuous model is conservative. Moreover, we establish suitable computational constraints under which the discrete energy of the system is positive. The method is consistent of second order, and is both stable and convergent. The numerical simulations shown here illustrate the most important features of our numerical methodology.

Continued-fraction representation of the Kraus map for non-Markovian reservoir damping

NASA Astrophysics Data System (ADS)

van Wonderen, A. J.; Suttorp, L. G.

2018-04-01

Quantum dissipation is studied for a discrete system that linearly interacts with a reservoir of harmonic oscillators at thermal equilibrium. Initial correlations between system and reservoir are assumed to be absent. The dissipative dynamics as determined by the unitary evolution of system and reservoir is described by a Kraus map consisting of an infinite number of matrices. For all Laplace-transformed Kraus matrices exact solutions are constructed in terms of continued fractions that depend on the pair correlation functions of the reservoir. By performing factorizations in the Kraus map a perturbation theory is set up that conserves in arbitrary perturbative order both positivity and probability of the density matrix. The latter is determined by an integral equation for a bitemporal matrix and a finite hierarchy for Kraus matrices. In the lowest perturbative order this hierarchy reduces to one equation for one Kraus matrix. Its solution is given by a continued fraction of a much simpler structure as compared to the non-perturbative case. In the lowest perturbative order our non-Markovian evolution equations are applied to the damped Jaynes–Cummings model. From the solution for the atomic density matrix it is found that the atom may remain in the state of maximum entropy for a significant time span that depends on the initial energy of the radiation field.

General PFG signal attenuation expressions for anisotropic anomalous diffusion by modified-Bloch equations

NASA Astrophysics Data System (ADS)

Lin, Guoxing

2018-05-01

Anomalous diffusion exists widely in polymer and biological systems. Pulsed-field gradient (PFG) anomalous diffusion is complicated, especially in the anisotropic case where limited research has been reported. A general PFG signal attenuation expression, including the finite gradient pulse (FGPW) effect for free general anisotropic fractional diffusion { 0 < α , β ≤ 2 } based on the fractional derivative, has not been obtained, where α and β are time and space derivative orders. It is essential to derive a general PFG signal attenuation expression including the FGPW effect for PFG anisotropic anomalous diffusion research. In this paper, two recently developed modified-Bloch equations, the fractal differential modified-Bloch equation and the fractional integral modified-Bloch equation, were extended to obtain general PFG signal attenuation expressions for anisotropic anomalous diffusion. Various cases of PFG anisotropic anomalous diffusion were investigated, including coupled and uncoupled anisotropic anomalous diffusion. The continuous-time random walk (CTRW) simulation was also carried out to support the theoretical results. The theory and the CTRW simulation agree with each other. The obtained signal attenuation expressions and the three-dimensional fractional modified-Bloch equations are important for analyzing PFG anisotropic anomalous diffusion in NMR and MRI.

Linking Models: Reasoning from Patterns to Tables and Equations

ERIC Educational Resources Information Center

Switzer, J. Matt

2013-01-01

Patterns are commonly used in middle years mathematics classrooms to teach students about functions and modelling with tables, graphs, and equations. Grade 6 students are expected to, "continue and create sequences involving whole numbers, fractions and decimals," and "describe the rule used to create the sequence." (Australian…

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

NASA Technical Reports Server (NTRS)

Rosenfeld, Moshe; Kwak, Dochan; Vinokur, Marcel

1992-01-01

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases.

Theory of a ring laser. [electromagnetic field and wave equations

NASA Technical Reports Server (NTRS)

Menegozzi, L. N.; Lamb, W. E., Jr.

1973-01-01

Development of a systematic formulation of the theory of a ring laser which is based on first principles and uses a well-known model for laser operation. A simple physical derivation of the electromagnetic field equations for a noninertial reference frame in uniform rotation is presented, and an attempt is made to clarify the nature of the Fox-Li modes for an open polygonal resonator. The polarization of the active medium is obtained by using a Fourier-series method which permits the formulation of a strong-signal theory, and solutions are given in terms of continued fractions. It is shown that when such a continued fraction is expanded to third order in the fields, the familiar small-signal ring-laser theory is obtained.

Analytical pricing of geometric Asian power options on an underlying driven by a mixed fractional Brownian motion

NASA Astrophysics Data System (ADS)

Zhang, Wei-Guo; Li, Zhe; Liu, Yong-Jun

2018-01-01

In this paper, we study the pricing problem of the continuously monitored fixed and floating strike geometric Asian power options in a mixed fractional Brownian motion environment. First, we derive both closed-form solutions and mixed fractional partial differential equations for fixed and floating strike geometric Asian power options based on delta-hedging strategy and partial differential equation method. Second, we present the lower and upper bounds of the prices of fixed and floating strike geometric Asian power options under the assumption that both risk-free interest rate and volatility are interval numbers. Finally, numerical studies are performed to illustrate the performance of our proposed pricing model.

Martian tidal pressure and wind fields obtained from the Mariner 9 infrared spectroscopy experiment

NASA Technical Reports Server (NTRS)

Pirraglia, J. A.; Conrath, B. J.

1973-01-01

Using temperature fields derived from the Mariner 9 infrared spectroscopy experiment, the Martian atmospheric tidal pressure and wind fields are calculated. Temperature as a function of local time, latitude, and atmospheric pressure level is obtained by secular and longitudinal averaging of the data. The resulting temperature field is approximated by a spherical harmonic expansion, retaining one symmetric and one asymmetric term for wavenumber zero and wavenumber one. Vertical averaging of the linearized momentum and continuity equations results in an inhomogeneous tidal equation for surface pressure fluctuations with the driving function related to the temperature field through the geopotential function and the hydrostatic equation. Solutions of the tidal equation show a diurnal fractional pressure amplitude approximately equal to one half of the vertically averaged diurnal fractional temperature amplitude.

Martian tidal pressure and wind fields obtained from the Mariner 9 infrared spectroscopy experiment

NASA Technical Reports Server (NTRS)

Pirraglia, J. A.; Conrath, B. J.

1974-01-01

Using temperature fields derived from the Mariner 9 infrared spectroscopy experiment, the Martian atmospheric tidal pressure and wind fields are calculated. Temperature as a function of local time, latitude, and atmospheric pressure level is obtained by secular and longitudinal averaging of the data. The resulting temperature field is approximated by a spherical harmonic expansion, retaining one symmetric and one asymmetric term each for wavenumber zero and wavenumber one. Vertical averaging of the linearized momentum and continuity equations results in an inhomogeneous tidal equation for surface pressure fluctuations with the driving function related to the temperature field through the geopotential function and the hydrostatic equation. Solutions of the tidal equation show a diurnal fractional pressure amplitude approximately equal to one-half the vertically averaged diurnal fractional temperature amplitude.

Optimal Control of Stochastic Systems Driven by Fractional Brownian Motions

DTIC Science & Technology

2014-10-09

problems for stochastic partial differential equations driven by fractional Brownian motions are explicitly solved. For the control of a continuous time...linear systems with Brownian motion or a discrete time linear system with a white Gaussian noise and costs 1. REPORT DATE (DD-MM-YYYY) 4. TITLE AND...Army Research Office P.O. Box 12211 Research Triangle Park, NC 27709-2211 stochastic optimal control, fractional Brownian motion , stochastic

Initial value problem of space dynamics in universal Stumpff anomaly

NASA Astrophysics Data System (ADS)

Sharaf, M. A.; Dwidar, H. R.

2018-05-01

In this paper, the initial value problem of space dynamics in universal Stumpff anomaly ψ is set up and developed in analytical and computational approach. For the analytical expansions, the linear independence of the functions U_{j} (ψ;σ); {j=0,1,2,3} are proved. The differential and recurrence equations satisfied by them and their relations with the elementary functions are given. The universal Kepler equation and its validations for different conic orbits are established together with the Lagrangian coefficients. Efficient representations of these functions are developed in terms of the continued fractions. For the computational developments we consider the following items: 1. Top-down algorithm for continued fraction evaluation. 2. One-point iteration formulae. 3. Determination of the coefficients of Kepler's equation. 4. Derivatives of Kepler's equation of any integer order. 5. Determination of the initial guess for the solution of the universal Kepler equation. Finally we give summary on the computational design for the initial value problem of space dynamics in universal Stumpff anomaly. This design based on the solution of the universal Kepler's equation by an iterative schemes of quadratic up to any desired order ℓ.

Nonlinear Conservation Laws and Finite Volume Methods

NASA Astrophysics Data System (ADS)

Leveque, Randall J.

Introduction Software Notation Classification of Differential Equations Derivation of Conservation Laws The Euler Equations of Gas Dynamics Dissipative Fluxes Source Terms Radiative Transfer and Isothermal Equations Multi-dimensional Conservation Laws The Shock Tube Problem Mathematical Theory of Hyperbolic Systems Scalar Equations Linear Hyperbolic Systems Nonlinear Systems The Riemann Problem for the Euler Equations Numerical Methods in One Dimension Finite Difference Theory Finite Volume Methods Importance of Conservation Form - Incorrect Shock Speeds Numerical Flux Functions Godunov's Method Approximate Riemann Solvers High-Resolution Methods Other Approaches Boundary Conditions Source Terms and Fractional Steps Unsplit Methods Fractional Step Methods General Formulation of Fractional Step Methods Stiff Source Terms Quasi-stationary Flow and Gravity Multi-dimensional Problems Dimensional Splitting Multi-dimensional Finite Volume Methods Grids and Adaptive Refinement Computational Difficulties Low-Density Flows Discrete Shocks and Viscous Profiles Start-Up Errors Wall Heating Slow-Moving Shocks Grid Orientation Effects Grid-Aligned Shocks Magnetohydrodynamics The MHD Equations One-Dimensional MHD Solving the Riemann Problem Nonstrict Hyperbolicity Stiffness The Divergence of B Riemann Problems in Multi-dimensional MHD Staggered Grids The 8-Wave Riemann Solver Relativistic Hydrodynamics Conservation Laws in Spacetime The Continuity Equation The 4-Momentum of a Particle The Stress-Energy Tensor Finite Volume Methods Multi-dimensional Relativistic Flow Gravitation and General Relativity References

Unsteady free convection flow of viscous fluids with analytical results by employing time-fractional Caputo-Fabrizio derivative (without singular kernel)

NASA Astrophysics Data System (ADS)

Ali Shah, Nehad; Mahsud, Yasir; Ali Zafar, Azhar

2017-10-01

This article introduces a theoretical study for unsteady free convection flow of an incompressible viscous fluid. The fluid flows near an isothermal vertical plate. The plate has a translational motion with time-dependent velocity. The equations governing the fluid flow are expressed in fractional differential equations by using a newly defined time-fractional Caputo-Fabrizio derivative without singular kernel. Explicit solutions for velocity, temperature and solute concentration are obtained by applying the Laplace transform technique. As the fractional parameter approaches to one, solutions for the ordinary fluid model are extracted from the general solutions of the fractional model. The results showed that, for the fractional model, the obtained solutions for velocity, temperature and concentration exhibit stationary jumps discontinuity across the plane at t=0 , while the solutions are continuous functions in the case of the ordinary model. Finally, numerical results for flow features at small-time are illustrated through graphs for various pertinent parameters.

A fast semi-discrete Kansa method to solve the two-dimensional spatiotemporal fractional diffusion equation

NASA Astrophysics Data System (ADS)

Sun, HongGuang; Liu, Xiaoting; Zhang, Yong; Pang, Guofei; Garrard, Rhiannon

2017-09-01

Fractional-order diffusion equations (FDEs) extend classical diffusion equations by quantifying anomalous diffusion frequently observed in heterogeneous media. Real-world diffusion can be multi-dimensional, requiring efficient numerical solvers that can handle long-term memory embedded in mass transport. To address this challenge, a semi-discrete Kansa method is developed to approximate the two-dimensional spatiotemporal FDE, where the Kansa approach first discretizes the FDE, then the Gauss-Jacobi quadrature rule solves the corresponding matrix, and finally the Mittag-Leffler function provides an analytical solution for the resultant time-fractional ordinary differential equation. Numerical experiments are then conducted to check how the accuracy and convergence rate of the numerical solution are affected by the distribution mode and number of spatial discretization nodes. Applications further show that the numerical method can efficiently solve two-dimensional spatiotemporal FDE models with either a continuous or discrete mixing measure. Hence this study provides an efficient and fast computational method for modeling super-diffusive, sub-diffusive, and mixed diffusive processes in large, two-dimensional domains with irregular shapes.

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Diethelm, Kai; Gray, Hugh R. (Technical Monitor)

2002-01-01

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

A procedure to construct exact solutions of nonlinear fractional differential equations.

PubMed

Güner, Özkan; Cevikel, Adem C

2014-01-01

We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions.

Non-linear continuous time random walk models★

NASA Astrophysics Data System (ADS)

Stage, Helena; Fedotov, Sergei

2017-11-01

A standard assumption of continuous time random walk (CTRW) processes is that there are no interactions between the random walkers, such that we obtain the celebrated linear fractional equation either for the probability density function of the walker at a certain position and time, or the mean number of walkers. The question arises how one can extend this equation to the non-linear case, where the random walkers interact. The aim of this work is to take into account this interaction under a mean-field approximation where the statistical properties of the random walker depend on the mean number of walkers. The implementation of these non-linear effects within the CTRW integral equations or fractional equations poses difficulties, leading to the alternative methodology we present in this work. We are concerned with non-linear effects which may either inhibit anomalous effects or induce them where they otherwise would not arise. Inhibition of these effects corresponds to a decrease in the waiting times of the random walkers, be this due to overcrowding, competition between walkers or an inherent carrying capacity of the system. Conversely, induced anomalous effects present longer waiting times and are consistent with symbiotic, collaborative or social walkers, or indirect pinpointing of favourable regions by their attractiveness. Contribution to the Topical Issue "Continuous Time Random Walk Still Trendy: Fifty-year History, Current State and Outlook", edited by Ryszard Kutner and Jaume Masoliver.

A Procedure to Construct Exact Solutions of Nonlinear Fractional Differential Equations

PubMed Central

Güner, Özkan; Cevikel, Adem C.

2014-01-01

We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions. PMID:24737972

Discrete adjoint of fractional step Navier-Stokes solver in generalized coordinates

NASA Astrophysics Data System (ADS)

Wang, Mengze; Mons, Vincent; Zaki, Tamer

2017-11-01

Optimization and control in transitional and turbulent flows require evaluation of gradients of the flow state with respect to the problem parameters. Using adjoint approaches, these high-dimensional gradients can be evaluated with a similar computational cost as the forward Navier-Stokes simulations. The adjoint algorithm can be obtained by discretizing the continuous adjoint Navier-Stokes equations or by deriving the adjoint to the discretized Navier-Stokes equations directly. The latter algorithm is necessary when the forward-adjoint relations must be satisfied to machine precision. In this work, our forward model is the fractional step solution to the Navier-Stokes equations in generalized coordinates, proposed by Rosenfeld, Kwak & Vinokur. We derive the corresponding discrete adjoint equations. We also demonstrate the accuracy of the combined forward-adjoint model, and its application to unsteady wall-bounded flows. This work has been partially funded by the Office of Naval Research (Grant N00014-16-1-2542).

August median streamflow on ungaged streams in Eastern Coastal Maine

USGS Publications Warehouse

Lombard, Pamela J.

2004-01-01

Methods for estimating August median streamflow were developed for ungaged, unregulated streams in eastern coastal Maine. The methods apply to streams with drainage areas ranging in size from 0.04 to 73.2 square miles and fraction of basin underlain by a sand and gravel aquifer ranging from 0 to 71 percent. The equations were developed with data from three long-term (greater than or equal to 10 years of record) continuous-record streamflow-gaging stations, 23 partial-record streamflow- gaging stations, and 5 short-term (less than 10 years of record) continuous-record streamflow-gaging stations. A mathematical technique for estimating a standard low-flow statistic, August median streamflow, at partial-record streamflow-gaging stations and short-term continuous-record streamflow-gaging stations was applied by relating base-flow measurements at these stations to concurrent daily streamflows at nearby long-term continuous-record streamflow-gaging stations (index stations). Generalized least-squares regression analysis (GLS) was used to relate estimates of August median streamflow at streamflow-gaging stations to basin characteristics at these same stations to develop equations that can be applied to estimate August median streamflow on ungaged streams. GLS accounts for different periods of record at the gaging stations and the cross correlation of concurrent streamflows among gaging stations. Thirty-one stations were used for the final regression equations. Two basin characteristics?drainage area and fraction of basin underlain by a sand and gravel aquifer?are used in the calculated regression equation to estimate August median streamflow for ungaged streams. The equation has an average standard error of prediction from -27 to 38 percent. A one-variable equation uses only drainage area to estimate August median streamflow when less accuracy is acceptable. This equation has an average standard error of prediction from -30 to 43 percent. Model error is larger than sampling error for both equations, indicating that additional or improved estimates of basin characteristics could be important to improved estimates of low-flow statistics. Weighted estimates of August median streamflow at partial- record or continuous-record gaging stations range from 0.003 to 31.0 cubic feet per second or from 0.1 to 0.6 cubic feet per second per square mile. Estimates of August median streamflow on ungaged streams in eastern coastal Maine, within the range of acceptable explanatory variables, range from 0.003 to 45 cubic feet per second or 0.1 to 0.6 cubic feet per second per square mile. Estimates of August median streamflow per square mile of drainage area generally increase as drainage area and fraction of basin underlain by a sand and gravel aquifer increase.

Comment on "Fractional quantum mechanics" and "Fractional Schrödinger equation"

NASA Astrophysics Data System (ADS)

Wei, Yuchuan

2016-06-01

In this Comment we point out some shortcomings in two papers [N. Laskin, Phys. Rev. E 62, 3135 (2000), 10.1103/PhysRevE.62.3135; N. Laskin, Phys. Rev. E 66, 056108 (2002), 10.1103/PhysRevE.66.056108]. We prove that the fractional uncertainty relation does not hold generally. The probability continuity equation in fractional quantum mechanics has a missing source term, which leads to particle teleportation, i.e., a particle can teleport from a place to another. Since the relativistic kinetic energy can be viewed as an approximate realization of the fractional kinetic energy, the particle teleportation should be an observable relativistic effect in quantum mechanics. With the help of this concept, superconductivity could be viewed as the teleportation of electrons from one side of a superconductor to another and superfluidity could be viewed as the teleportation of helium atoms from one end of a capillary tube to the other. We also point out how to teleport a particle to an arbitrary destination.

Comment on "Fractional quantum mechanics" and "Fractional Schrödinger equation".

PubMed

Wei, Yuchuan

2016-06-01

In this Comment we point out some shortcomings in two papers [N. Laskin, Phys. Rev. E 62, 3135 (2000)10.1103/PhysRevE.62.3135; N. Laskin, Phys. Rev. E 66, 056108 (2002)10.1103/PhysRevE.66.056108]. We prove that the fractional uncertainty relation does not hold generally. The probability continuity equation in fractional quantum mechanics has a missing source term, which leads to particle teleportation, i.e., a particle can teleport from a place to another. Since the relativistic kinetic energy can be viewed as an approximate realization of the fractional kinetic energy, the particle teleportation should be an observable relativistic effect in quantum mechanics. With the help of this concept, superconductivity could be viewed as the teleportation of electrons from one side of a superconductor to another and superfluidity could be viewed as the teleportation of helium atoms from one end of a capillary tube to the other. We also point out how to teleport a particle to an arbitrary destination.

Nonlinear vibration of viscoelastic beams described using fractional order derivatives

NASA Astrophysics Data System (ADS)

Lewandowski, Roman; Wielentejczyk, Przemysław

2017-07-01

The problem of non-linear, steady state vibration of beams, harmonically excited by harmonic forces is investigated in the paper. The viscoelastic material of the beams is described using the Zener rheological model with fractional derivatives. The constitutive equation, which contains derivatives of both stress and strain, significantly complicates the solution to the problem. The von Karman theory is applied to take into account geometric nonlinearities. Amplitude equations are obtained using the finite element method together with the harmonic balance method, and solved using the continuation method. The tangent matrix of the amplitude equations is determined in an explicit form. The stability of the steady-state solution is also examined. A parametric study is carried out to determine the influence of viscoelastic properties of the material on the beam's responses.

Fractional Diffusion Equations and Anomalous Diffusion

NASA Astrophysics Data System (ADS)

Evangelista, Luiz Roberto; Kaminski Lenzi, Ervin

2018-01-01

Preface; 1. Mathematical preliminaries; 2. A survey of the fractional calculus; 3. From normal to anomalous diffusion; 4. Fractional diffusion equations: elementary applications; 5. Fractional diffusion equations: surface effects; 6. Fractional nonlinear diffusion equation; 7. Anomalous diffusion: anisotropic case; 8. Fractional Schrödinger equations; 9. Anomalous diffusion and impedance spectroscopy; 10. The Poisson–Nernst–Planck anomalous (PNPA) models; References; Index.

Existence of solution for a general fractional advection-dispersion equation

NASA Astrophysics Data System (ADS)

Torres Ledesma, César E.

2018-05-01

In this work, we consider the existence of solution to the following fractional advection-dispersion equation -d/dt ( p {_{-∞}}It^{β }(u'(t)) + q {t}I_{∞}^{β }(u'(t))) + b(t)u = f(t, u(t)),t\\in R where β \\in (0,1) , _{-∞}It^{β } and tI_{∞}^{β } denote left and right Liouville-Weyl fractional integrals of order β respectively, 0

Fractional calculus and morphogen gradient formation

NASA Astrophysics Data System (ADS)

Yuste, Santos Bravo; Abad, Enrique; Lindenberg, Katja

2012-12-01

Some microscopic models for reactive systems where the reaction kinetics is limited by subdiffusion are described by means of reaction-subdiffusion equations where fractional derivatives play a key role. In particular, we consider subdiffusive particles described by means of a Continuous Time Random Walk (CTRW) model subject to a linear (first-order) death process. The resulting fractional equation is employed to study the developmental biology key problem of morphogen gradient formation for the case in which the morphogens are subdiffusive. If the morphogen degradation rate (reactivity) is constant, we find exponentially decreasing stationary concentration profiles, which are similar to the profiles found when the morphogens diffuse normally. However, for the case in which the degradation rate decays exponentially with the distance to the morphogen source, we find that the morphogen profiles are qualitatively different from the profiles obtained when the morphogens diffuse normally.

General pulsed-field gradient signal attenuation expression based on a fractional integral modified-Bloch equation

NASA Astrophysics Data System (ADS)

Lin, Guoxing

2018-10-01

Anomalous diffusion has been investigated in many polymer and biological systems. The analysis of PFG anomalous diffusion relies on the ability to obtain the signal attenuation expression. However, the general analytical PFG signal attenuation expression based on the fractional derivative has not been previously reported. Additionally, the reported modified-Bloch equations for PFG anomalous diffusion in the literature yielded different results due to their different forms. Here, a new integral type modified-Bloch equation based on the fractional derivative for PFG anomalous diffusion is proposed, which is significantly different from the conventional differential type modified-Bloch equation. The merit of the integral type modified-Bloch equation is that the original properties of the contributions from linear or nonlinear processes remain unchanged at the instant of the combination. From the modified-Bloch equation, the general solutions are derived, which includes the finite gradient pulse width (FGPW) effect. The numerical evaluation of these PFG signal attenuation expressions can be obtained either by the Adomian decomposition, or a direct integration method that is fast and practicable. The theoretical results agree with the continuous-time random walk (CTRW) simulations performed in this paper. Additionally, the relaxation effect in PFG anomalous diffusion is found to be different from that in PFG normal diffusion. The new modified-Bloch equations and their solutions provide a fundamental tool to analyze PFG anomalous diffusion in nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI).

Exact solutions for STO and (3+1)-dimensional KdV-ZK equations using (G‧/G2) -expansion method

NASA Astrophysics Data System (ADS)

Bibi, Sadaf; Mohyud-Din, Syed Tauseef; Ullah, Rahmat; Ahmed, Naveed; Khan, Umar

This article deals with finding some exact solutions of nonlinear fractional differential equations (NLFDEs) by applying a relatively new method known as (G‧/G2) -expansion method. Solutions of space-time fractional Sharma-Tasso-Olever (STO) equation of fractional order and (3+1)-dimensional KdV-Zakharov Kuznetsov (KdV-ZK) equation of fractional order are reckoned to demonstrate the validity of this method. The fractional derivative version of modified Riemann-Liouville, linked with Fractional complex transform is employed to transform fractional differential equations into the corresponding ordinary differential equations.

Lie symmetries and conservation laws for the time fractional Derrida-Lebowitz-Speer-Spohn equation

NASA Astrophysics Data System (ADS)

Rui, Wenjuan; Zhang, Xiangzhi

2016-05-01

This paper investigates the invariance properties of the time fractional Derrida-Lebowitz-Speer-Spohn (FDLSS) equation with Riemann-Liouville derivative. By using the Lie group analysis method of fractional differential equations, we derive Lie symmetries for the FDLSS equation. In a particular case of scaling transformations, we transform the FDLSS equation into a nonlinear ordinary fractional differential equation. Conservation laws for this equation are obtained with the aid of the new conservation theorem and the fractional generalization of the Noether operators.

Series expansion solutions for the multi-term time and space fractional partial differential equations in two- and three-dimensions

NASA Astrophysics Data System (ADS)

Ye, H.; Liu, F.; Turner, I.; Anh, V.; Burrage, K.

2013-09-01

Fractional partial differential equations with more than one fractional derivative in time describe some important physical phenomena, such as the telegraph equation, the power law wave equation, or the Szabo wave equation. In this paper, we consider two- and three-dimensional multi-term time and space fractional partial differential equations. The multi-term time-fractional derivative is defined in the Caputo sense, whose order belongs to the interval (1,2],(2,3],(3,4] or (0, m], and the space-fractional derivative is referred to as the fractional Laplacian form. We derive series expansion solutions based on a spectral representation of the Laplacian operator on a bounded region. Some applications are given for the two- and three-dimensional telegraph equation, power law wave equation and Szabo wave equation.

Periodicity computation of generalized mathematical biology problems involving delay differential equations.

PubMed

Jasim Mohammed, M; Ibrahim, Rabha W; Ahmad, M Z

2017-03-01

In this paper, we consider a low initial population model. Our aim is to study the periodicity computation of this model by using neutral differential equations, which are recognized in various studies including biology. We generalize the neutral Rayleigh equation for the third-order by exploiting the model of fractional calculus, in particular the Riemann-Liouville differential operator. We establish the existence and uniqueness of a periodic computational outcome. The technique depends on the continuation theorem of the coincidence degree theory. Besides, an example is presented to demonstrate the finding.

Exact solutions to the time-fractional differential equations via local fractional derivatives

NASA Astrophysics Data System (ADS)

Guner, Ozkan; Bekir, Ahmet

2018-01-01

This article utilizes the local fractional derivative and the exp-function method to construct the exact solutions of nonlinear time-fractional differential equations (FDEs). For illustrating the validity of the method, it is applied to the time-fractional Camassa-Holm equation and the time-fractional-generalized fifth-order KdV equation. Moreover, the exact solutions are obtained for the equations which are formed by different parameter values related to the time-fractional-generalized fifth-order KdV equation. This method is an reliable and efficient mathematical tool for solving FDEs and it can be applied to other non-linear FDEs.

Analytical solutions of the space-time fractional Telegraph and advection-diffusion equations

NASA Astrophysics Data System (ADS)

Tawfik, Ashraf M.; Fichtner, Horst; Schlickeiser, Reinhard; Elhanbaly, A.

2018-02-01

The aim of this paper is to develop a fractional derivative model of energetic particle transport for both uniform and non-uniform large-scale magnetic field by studying the fractional Telegraph equation and the fractional advection-diffusion equation. Analytical solutions of the space-time fractional Telegraph equation and space-time fractional advection-diffusion equation are obtained by use of the Caputo fractional derivative and the Laplace-Fourier technique. The solutions are given in terms of Fox's H function. As an illustration they are applied to the case of solar energetic particles.

Symmetry classification of time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Naeem, I.; Khan, M. D.

2017-01-01

In this article, a new approach is proposed to construct the symmetry groups for a class of fractional differential equations which are expressed in the modified Riemann-Liouville fractional derivative. We perform a complete group classification of a nonlinear fractional diffusion equation which arises in fractals, acoustics, control theory, signal processing and many other applications. Introducing the suitable transformations, the fractional derivatives are converted to integer order derivatives and in consequence the nonlinear fractional diffusion equation transforms to a partial differential equation (PDE). Then the Lie symmetries are computed for resulting PDE and using inverse transformations, we derive the symmetries for fractional diffusion equation. All cases are discussed in detail and results for symmetry properties are compared for different values of α. This study provides a new way of computing symmetries for a class of fractional differential equations.

Analytic solution for the space-time fractional Klein-Gordon and coupled conformable Boussinesq equations

NASA Astrophysics Data System (ADS)

Shallal, Muhannad A.; Jabbar, Hawraz N.; Ali, Khalid K.

2018-03-01

In this paper, we constructed a travelling wave solution for space-time fractional nonlinear partial differential equations by using the modified extended Tanh method with Riccati equation. The method is used to obtain analytic solutions for the space-time fractional Klein-Gordon and coupled conformable space-time fractional Boussinesq equations. The fractional complex transforms and the properties of modified Riemann-Liouville derivative have been used to convert these equations into nonlinear ordinary differential equations.

An approach for modeling thermal destruction of hazardous wastes in circulating fluidized bed incinerator.

PubMed

Patil, M P; Sonolikar, R L

2008-10-01

This paper presents a detailed computational fluid dynamics (CFD) based approach for modeling thermal destruction of hazardous wastes in a circulating fluidized bed (CFB) incinerator. The model is based on Eular - Lagrangian approach in which gas phase (continuous phase) is treated in a Eularian reference frame, whereas the waste particulate (dispersed phase) is treated in a Lagrangian reference frame. The reaction chemistry hasbeen modeled through a mixture fraction/ PDF approach. The conservation equations for mass, momentum, energy, mixture fraction and other closure equations have been solved using a general purpose CFD code FLUENT4.5. Afinite volume method on a structured grid has been used for solution of governing equations. The model provides detailed information on the hydrodynamics (gas velocity, particulate trajectories), gas composition (CO, CO2, O2) and temperature inside the riser. The model also allows different operating scenarios to be examined in an efficient manner.

Lie symmetry analysis, exact solutions and conservation laws for the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation

NASA Astrophysics Data System (ADS)

Baleanu, Dumitru; Inc, Mustafa; Yusuf, Abdullahi; Aliyu, Aliyu Isa

2018-06-01

In this work, we investigate the Lie symmetry analysis, exact solutions and conservation laws (Cls) to the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera (CDGDK) equation with Riemann-Liouville (RL) derivative. The time fractional CDGDK is reduced to nonlinear ordinary differential equation (ODE) of fractional order. New exact traveling wave solutions for the time fractional CDGDK are obtained by fractional sub-equation method. In the reduced equation, the derivative is in Erdelyi-Kober (EK) sense. Ibragimov's nonlocal conservation method is applied to construct Cls for time fractional CDGDK.

High-accuracy power series solutions with arbitrarily large radius of convergence for the fractional nonlinear Schrödinger-type equations

NASA Astrophysics Data System (ADS)

Khawaja, U. Al; Al-Refai, M.; Shchedrin, Gavriil; Carr, Lincoln D.

2018-06-01

Fractional nonlinear differential equations present an interplay between two common and important effective descriptions used to simplify high dimensional or more complicated theories: nonlinearity and fractional derivatives. These effective descriptions thus appear commonly in physical and mathematical modeling. We present a new series method providing systematic controlled accuracy for solutions of fractional nonlinear differential equations, including the fractional nonlinear Schrödinger equation and the fractional nonlinear diffusion equation. The method relies on spatially iterative use of power series expansions. Our approach permits an arbitrarily large radius of convergence and thus solves the typical divergence problem endemic to power series approaches. In the specific case of the fractional nonlinear Schrödinger equation we find fractional generalizations of cnoidal waves of Jacobi elliptic functions as well as a fractional bright soliton. For the fractional nonlinear diffusion equation we find the combination of fractional and nonlinear effects results in a more strongly localized solution which nevertheless still exhibits power law tails, albeit at a much lower density.

Multiscale functions, scale dynamics, and applications to partial differential equations

NASA Astrophysics Data System (ADS)

Cresson, Jacky; Pierret, Frédéric

2016-05-01

Modeling phenomena from experimental data always begins with a choice of hypothesis on the observed dynamics such as determinism, randomness, and differentiability. Depending on these choices, different behaviors can be observed. The natural question associated to the modeling problem is the following: "With a finite set of data concerning a phenomenon, can we recover its underlying nature? From this problem, we introduce in this paper the definition of multi-scale functions, scale calculus, and scale dynamics based on the time scale calculus [see Bohner, M. and Peterson, A., Dynamic Equations on Time Scales: An Introduction with Applications (Springer Science & Business Media, 2001)] which is used to introduce the notion of scale equations. These definitions will be illustrated on the multi-scale Okamoto's functions. Scale equations are analysed using scale regimes and the notion of asymptotic model for a scale equation under a particular scale regime. The introduced formalism explains why a single scale equation can produce distinct continuous models even if the equation is scale invariant. Typical examples of such equations are given by the scale Euler-Lagrange equation. We illustrate our results using the scale Newton's equation which gives rise to a non-linear diffusion equation or a non-linear Schrödinger equation as asymptotic continuous models depending on the particular fractional scale regime which is considered.

Analyzing Lie symmetry and constructing conservation laws for time-fractional Benny-Lin equation

NASA Astrophysics Data System (ADS)

Rashidi, Saeede; Hejazi, S. Reza

This paper investigates the invariance properties of the time fractional Benny-Lin equation with Riemann-Liouville and Caputo derivatives. This equation can be reduced to the Kawahara equation, fifth-order Kdv equation, the Kuramoto-Sivashinsky equation and Navier-Stokes equation. By using the Lie group analysis method of fractional differential equations (FDEs), we derive Lie symmetries for the Benny-Lin equation. Conservation laws for this equation are obtained with the aid of the concept of nonlinear self-adjointness and the fractional generalization of the Noether’s operators. Furthermore, by means of the invariant subspace method, exact solutions of the equation are also constructed.

Long-Term Dynamics of Autonomous Fractional Differential Equations

NASA Astrophysics Data System (ADS)

Liu, Tao; Xu, Wei; Xu, Yong; Han, Qun

This paper aims to investigate long-term dynamic behaviors of autonomous fractional differential equations with effective numerical method. The long-term dynamic behaviors predict where systems are heading after long-term evolution. We make some modification and transplant cell mapping methods to autonomous fractional differential equations. The mapping time duration of cell mapping is enlarged to deal with the long memory effect. Three illustrative examples, i.e. fractional Lotka-Volterra equation, fractional van der Pol oscillator and fractional Duffing equation, are studied with our revised generalized cell mapping method. We obtain long-term dynamics, such as attractors, basins of attraction, and saddles. Compared with some existing stability and numerical results, the validity of our method is verified. Furthermore, we find that the fractional order has its effect on the long-term dynamics of autonomous fractional differential equations.

Exp-function method for solving fractional partial differential equations.

PubMed

Zheng, Bin

2013-01-01

We extend the Exp-function method to fractional partial differential equations in the sense of modified Riemann-Liouville derivative based on nonlinear fractional complex transformation. For illustrating the validity of this method, we apply it to the space-time fractional Fokas equation and the nonlinear fractional Sharma-Tasso-Olver (STO) equation. As a result, some new exact solutions for them are successfully established.

Fractional-order difference equations for physical lattices and some applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tarasov, Vasily E., E-mail: tarasov@theory.sinp.msu.ru

2015-10-15

Fractional-order operators for physical lattice models based on the Grünwald-Letnikov fractional differences are suggested. We use an approach based on the models of lattices with long-range particle interactions. The fractional-order operators of differentiation and integration on physical lattices are represented by kernels of lattice long-range interactions. In continuum limit, these discrete operators of non-integer orders give the fractional-order derivatives and integrals with respect to coordinates of the Grünwald-Letnikov types. As examples of the fractional-order difference equations for physical lattices, we give difference analogs of the fractional nonlocal Navier-Stokes equations and the fractional nonlocal Maxwell equations for lattices with long-range interactions.more » Continuum limits of these fractional-order difference equations are also suggested.« less

Lie Symmetry Analysis and Explicit Solutions of the Time Fractional Fifth-Order KdV Equation

PubMed Central

Wang, Gang wei; Xu, Tian zhou; Feng, Tao

2014-01-01

In this paper, using the Lie group analysis method, we study the invariance properties of the time fractional fifth-order KdV equation. A systematic research to derive Lie point symmetries to time fractional fifth-order KdV equation is performed. In the sense of point symmetry, all of the vector fields and the symmetry reductions of the fractional fifth-order KdV equation are obtained. At last, by virtue of the sub-equation method, some exact solutions to the fractional fifth-order KdV equation are provided. PMID:24523885

THE FUNDAMENTAL SOLUTIONS FOR MULTI-TERM MODIFIED POWER LAW WAVE EQUATIONS IN A FINITE DOMAIN.

PubMed

Jiang, H; Liu, F; Meerschaert, M M; McGough, R J

2013-01-01

Fractional partial differential equations with more than one fractional derivative term in time, such as the Szabo wave equation, or the power law wave equation, describe important physical phenomena. However, studies of these multi-term time-space or time fractional wave equations are still under development. In this paper, multi-term modified power law wave equations in a finite domain are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals (1, 2], [2, 3), [2, 4) or (0, n ) ( n > 2), respectively. Analytical solutions of the multi-term modified power law wave equations are derived. These new techniques are based on Luchko's Theorem, a spectral representation of the Laplacian operator, a method of separating variables and fractional derivative techniques. Then these general methods are applied to the special cases of the Szabo wave equation and the power law wave equation. These methods and techniques can also be extended to other kinds of the multi-term time-space fractional models including fractional Laplacian.

Oscillation of a class of fractional differential equations with damping term.

PubMed

Qin, Huizeng; Zheng, Bin

2013-01-01

We investigate the oscillation of a class of fractional differential equations with damping term. Based on a certain variable transformation, the fractional differential equations are converted into another differential equations of integer order with respect to the new variable. Then, using Riccati transformation, inequality, and integration average technique, some new oscillatory criteria for the equations are established. As for applications, oscillation for two certain fractional differential equations with damping term is investigated by the use of the presented results.

Toward lattice fractional vector calculus

NASA Astrophysics Data System (ADS)

Tarasov, Vasily E.

2014-09-01

An analog of fractional vector calculus for physical lattice models is suggested. We use an approach based on the models of three-dimensional lattices with long-range inter-particle interactions. The lattice analogs of fractional partial derivatives are represented by kernels of lattice long-range interactions, where the Fourier series transformations of these kernels have a power-law form with respect to wave vector components. In the continuum limit, these lattice partial derivatives give derivatives of non-integer order with respect to coordinates. In the three-dimensional description of the non-local continuum, the fractional differential operators have the form of fractional partial derivatives of the Riesz type. As examples of the applications of the suggested lattice fractional vector calculus, we give lattice models with long-range interactions for the fractional Maxwell equations of non-local continuous media and for the fractional generalization of the Mindlin and Aifantis continuum models of gradient elasticity.

Dynamical property analysis of fractionally damped van der pol oscillator and its application

NASA Astrophysics Data System (ADS)

Zhong, Qiuhui; Zhang, Chunrui

2012-01-01

In this paper, the fractionally damped van der pol equation was studied. Firstly, the fractionally damped van der pol equation was transformed into a set of integer order equations. Then the Lyapunov exponents diagram was given. Secondly, it was transformed into a set of fractional integral equations and solved by a predictor-corrector method. The time domain diagrams and phase trajectory were used to describe the dynamic behavior. Finally, the fractionally damped van der pol equation was used to detect a weak signal.

The fractional dynamics of quantum systems

NASA Astrophysics Data System (ADS)

Lu, Longzhao; Yu, Xiangyang

2018-05-01

The fractional dynamic process of a quantum system is a novel and complicated problem. The establishment of a fractional dynamic model is a significant attempt that is expected to reveal the mechanism of fractional quantum system. In this paper, a generalized time fractional Schrödinger equation is proposed. To study the fractional dynamics of quantum systems, we take the two-level system as an example and derive the time fractional equations of motion. The basic properties of the system are investigated by solving this set of equations in the absence of light field analytically. Then, when the system is subject to the light field, the equations are solved numerically. It shows that the two-level system described by the time fractional Schrödinger equation we proposed is a confirmable system.

Quasinormal modes of Reissner-Nordstrom black holes

NASA Technical Reports Server (NTRS)

Leaver, Edward W.

1990-01-01

A matrix-eigenvalue algorithm is presented for accurately computing the quasi-normal frequencies and modes of charged static blackholes. The method is then refined through the introduction of a continued-fraction step. The approach should generalize to a variety of nonseparable wave equations, including the Kerr-Newman case of charged rotating blackholes.

Analytical study of fractional equations describing anomalous diffusion of energetic particles

NASA Astrophysics Data System (ADS)

Tawfik, A. M.; Fichtner, H.; Schlickeiser, R.; Elhanbaly, A.

2017-06-01

To present the main influence of anomalous diffusion on the energetic particle propagation, the fractional derivative model of transport is developed by deriving the fractional modified Telegraph and Rayleigh equations. Analytical solutions of the fractional modified Telegraph and the fractional Rayleigh equations, which are defined in terms of Caputo fractional derivatives, are obtained by using the Laplace transform and the Mittag-Leffler function method. The solutions of these fractional equations are given in terms of special functions like Fox’s H, Mittag-Leffler, Hermite and Hyper-geometric functions. The predicted travelling pulse solutions are discussed in each case for different values of fractional order.

Nonlinear acoustic wave equations with fractional loss operators.

PubMed

Prieur, Fabrice; Holm, Sverre

2011-09-01

Fractional derivatives are well suited to describe wave propagation in complex media. When introduced in classical wave equations, they allow a modeling of attenuation and dispersion that better describes sound propagation in biological tissues. Traditional constitutive equations from solid mechanics and heat conduction are modified using fractional derivatives. They are used to derive a nonlinear wave equation which describes attenuation and dispersion laws that match observations. This wave equation is a generalization of the Westervelt equation, and also leads to a fractional version of the Khokhlov-Zabolotskaya-Kuznetsov and Burgers' equations. © 2011 Acoustical Society of America

A more fundamental approach to the derivation of nonlinear acoustic wave equations with fractional loss operators (L).

PubMed

Prieur, Fabrice; Vilenskiy, Gregory; Holm, Sverre

2012-10-01

A corrected derivation of nonlinear wave propagation equations with fractional loss operators is presented. The fundamental approach is based on fractional formulations of the stress-strain and heat flux definitions but uses the energy equation and thermodynamic identities to link density and pressure instead of an erroneous fractional form of the entropy equation as done in Prieur and Holm ["Nonlinear acoustic wave equations with fractional loss operators," J. Acoust. Soc. Am. 130(3), 1125-1132 (2011)]. The loss operator of the obtained nonlinear wave equations differs from the previous derivations as well as the dispersion equation, but when approximating for low frequencies the expressions for the frequency dependent attenuation and velocity dispersion remain unchanged.

Formulation, Implementation and Validation of a Two-Fluid model in a Fuel Cell CFD Code

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jain, Kunal; Cole, J. Vernon; Kumar, Sanjiv

2008-12-01

Water management is one of the main challenges in PEM Fuel Cells. While water is essential for membrane electrical conductivity, excess liquid water leads to flooding of catalyst layers. Despite the fact that accurate prediction of two-phase transport is key for optimal water management, understanding of the two-phase transport in fuel cells is relatively poor. Wang et. al. have studied the two-phase transport in the channel and diffusion layer separately using a multiphase mixture model. The model fails to accurately predict saturation values for high humidity inlet streams. Nguyen et. al. developed a two-dimensional, two-phase, isothermal, isobaric, steady state modelmore » of the catalyst and gas diffusion layers. The model neglects any liquid in the channel. Djilali et. al. developed a three-dimensional two-phase multicomponent model. The model is an improvement over previous models, but neglects drag between the liquid and the gas phases in the channel. In this work, we present a comprehensive two-fluid model relevant to fuel cells. Models for two-phase transport through Channel, Gas Diffusion Layer (GDL) and Channel-GDL interface, are discussed. In the channel, the gas and liquid pressures are assumed to be same. The surface tension effects in the channel are incorporated using the continuum surface force (CSF) model. The force at the surface is expressed as a volumetric body force and added as a source to the momentum equation. In the GDL, the gas and liquid are assumed to be at different pressures. The difference in the pressures (capillary pressure) is calculated using an empirical correlations. At the Channel-GDL interface, the wall adhesion affects need to be taken into account. SIMPLE-type methods recast the continuity equation into a pressure-correction equation, the solution of which then provides corrections for velocities and pressures. However, in the two-fluid model, the presence of two phasic continuity equations gives more freedom and more complications. A general approach would be to form a mixture continuity equation by linearly combining the phasic continuity equations using appropriate weighting factors. Analogous to mixture equation for pressure correction, a difference equation is used for the volume/phase fraction by taking the difference between the phasic continuity equations. The relative advantages of the above mentioned algorithmic variants for computing pressure correction and volume fractions are discussed and quantitatively assessed. Preliminary model validation is done for each component of the fuel cell. The two-phase transport in the channel is validated using empirical correlations. Transport in the GDL is validated against results obtained from LBM and VOF simulation techniques. The Channel-GDL interface transport will be validated against experiment and empirical correlation of droplet detachment at the interface.« less

Analytical approach for the fractional differential equations by using the extended tanh method

NASA Astrophysics Data System (ADS)

Pandir, Yusuf; Yildirim, Ayse

2018-07-01

In this study, we consider analytical solutions of space-time fractional derivative foam drainage equation, the nonlinear Korteweg-de Vries equation with time and space-fractional derivatives and time-fractional reaction-diffusion equation by using the extended tanh method. The fractional derivatives are defined in the modified Riemann-Liouville context. As a result, various exact analytical solutions consisting of trigonometric function solutions, kink-shaped soliton solutions and new exact solitary wave solutions are obtained.

Fractional Klein-Gordon equation composed of Jumarie fractional derivative and its interpretation by a smoothness parameter

NASA Astrophysics Data System (ADS)

Ghosh, Uttam; Banerjee, Joydip; Sarkar, Susmita; Das, Shantanu

2018-06-01

Klein-Gordon equation is one of the basic steps towards relativistic quantum mechanics. In this paper, we have formulated fractional Klein-Gordon equation via Jumarie fractional derivative and found two types of solutions. Zero-mass solution satisfies photon criteria and non-zero mass satisfies general theory of relativity. Further, we have developed rest mass condition which leads us to the concept of hidden wave. Classical Klein-Gordon equation fails to explain a chargeless system as well as a single-particle system. Using the fractional Klein-Gordon equation, we can overcome the problem. The fractional Klein-Gordon equation also leads to the smoothness parameter which is the measurement of the bumpiness of space. Here, by using this smoothness parameter, we have defined and interpreted the various cases.

Equation-of-State Scaling Factors

DOE Office of Scientific and Technical Information (OSTI.GOV)

Scannapieco, Anthony J.

2016-06-28

Equation-of-State scaling factors are needed when using a tabular EOS in which the user de ned material isotopic fractions di er from the actual isotopic fractions used by the table. Additionally, if a material is dynamically changing its isotopic structure, then an EOS scaling will again be needed, and will vary in time and location. The procedure that allows use of a table to obtain information about a similar material with average atomic mass Ms and average atomic number Zs is described below. The procedure is exact for a fully ionized ideal gas. However, if the atomic number is replacemore » by the e ective ionization state the procedure can be applied to partially ionized material as well, which extends the applicability of the scaling approximation continuously from low to high temperatures.« less

The convergence of the order sequence and the solution function sequence on fractional partial differential equation

NASA Astrophysics Data System (ADS)

Rusyaman, E.; Parmikanti, K.; Chaerani, D.; Asefan; Irianingsih, I.

2018-03-01

One of the application of fractional ordinary differential equation is related to the viscoelasticity, i.e., a correlation between the viscosity of fluids and the elasticity of solids. If the solution function develops into function with two or more variables, then its differential equation must be changed into fractional partial differential equation. As the preliminary study for two variables viscoelasticity problem, this paper discusses about convergence analysis of function sequence which is the solution of the homogenous fractional partial differential equation. The method used to solve the problem is Homotopy Analysis Method. The results show that if given two real number sequences (αn) and (βn) which converge to α and β respectively, then the solution function sequences of fractional partial differential equation with order (αn, βn) will also converge to the solution function of fractional partial differential equation with order (α, β).

Correcting the initialization of models with fractional derivatives via history-dependent conditions

NASA Astrophysics Data System (ADS)

Du, Maolin; Wang, Zaihua

2016-04-01

Fractional differential equations are more and more used in modeling memory (history-dependent, non-local, or hereditary) phenomena. Conventional initial values of fractional differential equations are defined at a point, while recent works define initial conditions over histories. We prove that the conventional initialization of fractional differential equations with a Riemann-Liouville derivative is wrong with a simple counter-example. The initial values were assumed to be arbitrarily given for a typical fractional differential equation, but we find one of these values can only be zero. We show that fractional differential equations are of infinite dimensions, and the initial conditions, initial histories, are defined as functions over intervals. We obtain the equivalent integral equation for Caputo case. With a simple fractional model of materials, we illustrate that the recovery behavior is correct with the initial creep history, but is wrong with initial values at the starting point of the recovery. We demonstrate the application of initial history by solving a forced fractional Lorenz system numerically.

Fractional Number Operator and Associated Fractional Diffusion Equations

NASA Astrophysics Data System (ADS)

Rguigui, Hafedh

2018-03-01

In this paper, we study the fractional number operator as an analog of the finite-dimensional fractional Laplacian. An important relation with the Ornstein-Uhlenbeck process is given. Using a semigroup approach, the solution of the Cauchy problem associated to the fractional number operator is presented. By means of the Mittag-Leffler function and the Laplace transform, we give the solution of the Caputo time fractional diffusion equation and Riemann-Liouville time fractional diffusion equation in infinite dimensions associated to the fractional number operator.

Continuous time random walk with local particle-particle interaction

NASA Astrophysics Data System (ADS)

Xu, Jianping; Jiang, Guancheng

2018-05-01

The continuous time random walk (CTRW) is often applied to the study of particle motion in disordered media. Yet most such applications do not allow for particle-particle (walker-walker) interaction. In this paper, we consider a CTRW with particle-particle interaction; however, for simplicity, we restrain the interaction to be local. The generalized Chapman-Kolmogorov equation is modified by introducing a perturbation function that fluctuates around 1, which models the effect of interaction. Subsequently, a time-fractional nonlinear advection-diffusion equation is derived from this walking system. Under the initial condition of condensed particles at the origin and the free-boundary condition, we numerically solve this equation with both attractive and repulsive particle-particle interactions. Moreover, a Monte Carlo simulation is devised to verify the results of the above numerical work. The equation and the simulation unanimously predict that this walking system converges to the conventional one in the long-time limit. However, for systems where the free-boundary condition and long-time limit are not simultaneously satisfied, this convergence does not hold.

Periodicity and positivity of a class of fractional differential equations.

PubMed

Ibrahim, Rabha W; Ahmad, M Z; Mohammed, M Jasim

2016-01-01

Fractional differential equations have been discussed in this study. We utilize the Riemann-Liouville fractional calculus to implement it within the generalization of the well known class of differential equations. The Rayleigh differential equation has been generalized of fractional second order. The existence of periodic and positive outcome is established in a new method. The solution is described in a fractional periodic Sobolev space. Positivity of outcomes is considered under certain requirements. We develop and extend some recent works. An example is constructed.

An enriched finite element method to fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam

2017-08-01

In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.

A new approach to the Schrödinger equation with rational potentials

NASA Astrophysics Data System (ADS)

Dong, Ming-de; Chu, Jue-Hui

1984-04-01

A new analytic theory is established for the Schrödinger equation with a rational potential, including a complete classification of the regular eigenfunctions into three different types, an exact method of obtaining wavefunctions, an explicit formulation of the spectral equation (3 x 3 determinant) etc. All representations are exhibited in a unifying way via function-theoretic methods and therefore given in explicit form, in contrast to the prevailing discussion appealing to perturbation or variation methods or continued-fraction techniques. The irregular eigenfunctions at infinity can be obtained analogously and will be discussed separately as another solvable case for singular potentials.

Systematic derivation of reaction-diffusion equations with distributed delays and relations to fractional reaction-diffusion equations and hyperbolic transport equations: application to the theory of Neolithic transition.

PubMed

Vlad, Marcel Ovidiu; Ross, John

2002-12-01

We introduce a general method for the systematic derivation of nonlinear reaction-diffusion equations with distributed delays. We study the interactions among different types of moving individuals (atoms, molecules, quasiparticles, biological organisms, etc). The motion of each species is described by the continuous time random walk theory, analyzed in the literature for transport problems, whereas the interactions among the species are described by a set of transformation rates, which are nonlinear functions of the local concentrations of the different types of individuals. We use the time interval between two jumps (the transition time) as an additional state variable and obtain a set of evolution equations, which are local in time. In order to make a connection with the transport models used in the literature, we make transformations which eliminate the transition time and derive a set of nonlocal equations which are nonlinear generalizations of the so-called generalized master equations. The method leads under different specified conditions to various types of nonlocal transport equations including a nonlinear generalization of fractional diffusion equations, hyperbolic reaction-diffusion equations, and delay-differential reaction-diffusion equations. Thus in the analysis of a given problem we can fit to the data the type of reaction-diffusion equation and the corresponding physical and kinetic parameters. The method is illustrated, as a test case, by the study of the neolithic transition. We introduce a set of assumptions which makes it possible to describe the transition from hunting and gathering to agriculture economics by a differential delay reaction-diffusion equation for the population density. We derive a delay evolution equation for the rate of advance of agriculture, which illustrates an application of our analysis.

Solution of a modified fractional diffusion equation

NASA Astrophysics Data System (ADS)

Langlands, T. A. M.

2006-07-01

Recently, a modified fractional diffusion equation has been proposed [I. Sokolov, J. Klafter, From diffusion to anomalous diffusion: a century after Einstein's brownian motion, Chaos 15 (2005) 026103; A.V. Chechkin, R. Gorenflo, I.M. Sokolov, V.Yu. Gonchar, Distributed order time fractional diffusion equation, Frac. Calc. Appl. Anal. 6 (3) (2003) 259279; I.M. Sokolov, A.V. Checkin, J. Klafter, Distributed-order fractional kinetics, Acta. Phys. Pol. B 35 (2004) 1323.] for describing processes that become less anomalous as time progresses by the inclusion of a second fractional time derivative acting on the diffusion term. In this letter we give the solution of the modified equation on an infinite domain. In contrast to the solution of the traditional fractional diffusion equation, the solution of the modified equation requires an infinite series of Fox functions instead of a single Fox function.

Exact solutions of fractional mBBM equation and coupled system of fractional Boussinesq-Burgers

NASA Astrophysics Data System (ADS)

Javeed, Shumaila; Saif, Summaya; Waheed, Asif; Baleanu, Dumitru

2018-06-01

The new exact solutions of nonlinear fractional partial differential equations (FPDEs) are established by adopting first integral method (FIM). The Riemann-Liouville (R-L) derivative and the local conformable derivative definitions are used to deal with the fractional order derivatives. The proposed method is applied to get exact solutions for space-time fractional modified Benjamin-Bona-Mahony (mBBM) equation and coupled time-fractional Boussinesq-Burgers equation. The suggested technique is easily applicable and effectual which can be implemented successfully to obtain the solutions for different types of nonlinear FPDEs.

The impact of fraction magnitude knowledge on algebra performance and learning.

PubMed

Booth, Julie L; Newton, Kristie J; Twiss-Garrity, Laura K

2014-02-01

Knowledge of fractions is thought to be crucial for success with algebra, but empirical evidence supporting this conjecture is just beginning to emerge. In the current study, Algebra 1 students completed magnitude estimation tasks on three scales (0-1 [fractions], 0-1,000,000, and 0-62,571) just before beginning their unit on equation solving. Results indicated that fraction magnitude knowledge, and not whole number knowledge, was especially related to students' pretest knowledge of equation solving and encoding of equation features. Pretest fraction knowledge was also predictive of students' improvement in equation solving and equation encoding skills. Students' placement of unit fractions (e.g., those with a numerator of 1) was not especially useful for predicting algebra performance and learning in this population. Placement of non-unit fractions was more predictive, suggesting that proportional reasoning skills might be an important link between fraction knowledge and learning algebra. Copyright © 2013 Elsevier Inc. All rights reserved.

Analytical solutions for coupling fractional partial differential equations with Dirichlet boundary conditions

NASA Astrophysics Data System (ADS)

Ding, Xiao-Li; Nieto, Juan J.

2017-11-01

In this paper, we consider the analytical solutions of coupling fractional partial differential equations (FPDEs) with Dirichlet boundary conditions on a finite domain. Firstly, the method of successive approximations is used to obtain the analytical solutions of coupling multi-term time fractional ordinary differential equations. Then, the technique of spectral representation of the fractional Laplacian operator is used to convert the coupling FPDEs to the coupling multi-term time fractional ordinary differential equations. By applying the obtained analytical solutions to the resulting multi-term time fractional ordinary differential equations, the desired analytical solutions of the coupling FPDEs are given. Our results are applied to derive the analytical solutions of some special cases to demonstrate their applicability.

Number Theoretic Background

NASA Astrophysics Data System (ADS)

Rudnick, Z.

Contents: 1. Introduction 2. Divisibility 2.1. Basics on Divisibility 2.2. The Greatest Common Divisor 2.3. The Euclidean Algorithm 2.4. The Diophantine Equation ax+by=c 3. Prime Numbers 3.1. The Fundamental Theorem of Arithmetic 3.2. There Are Infinitely Many Primes 3.3. The Density of Primes 3.4. Primes in Arithmetic Progressions 4. Continued Fractions 5. Modular Arithmetic 5.1. Congruences 5.2. Modular Inverses 5.3. The Chinese Remainder Theorem 5.4. The Structure of the Multiplicative Group (Z/NZ)^* 5.5. Primitive Roots 6. Quadratic Congruences 6.1. Euler's Criterion 6.2. The Legendre Symbol and Quadratic Reciprocity 7. Pell's Equation 7.1. The Group Law 7.2. Integer Solutions 7.3. Finding the Fundamental Solution 8. The Riemann Zeta Function 8.1 Analytic Continuation and Functinal Equation of ζ(s) 8.2 Connecting the Primes and the Zeros of ζ(s) 8.3 The Riemann Hypothesis References

Fractional calculus in hydrologic modeling: A numerical perspective

PubMed Central

Benson, David A.; Meerschaert, Mark M.; Revielle, Jordan

2013-01-01

Fractional derivatives can be viewed either as handy extensions of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Lévy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Lévy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus. PMID:23524449

Time-fractional Cahn-Allen and time-fractional Klein-Gordon equations: Lie symmetry analysis, explicit solutions and convergence analysis

NASA Astrophysics Data System (ADS)

Inc, Mustafa; Yusuf, Abdullahi; Isa Aliyu, Aliyu; Baleanu, Dumitru

2018-03-01

This research analyzes the symmetry analysis, explicit solutions and convergence analysis to the time fractional Cahn-Allen (CA) and time-fractional Klein-Gordon (KG) equations with Riemann-Liouville (RL) derivative. The time fractional CA and time fractional KG are reduced to respective nonlinear ordinary differential equation of fractional order. We solve the reduced fractional ODEs using an explicit power series method. The convergence analysis for the obtained explicit solutions are investigated. Some figures for the obtained explicit solutions are also presented.

New analytical exact solutions of time fractional KdV-KZK equation by Kudryashov methods

NASA Astrophysics Data System (ADS)

S Saha, Ray

2016-04-01

In this paper, new exact solutions of the time fractional KdV-Khokhlov-Zabolotskaya-Kuznetsov (KdV-KZK) equation are obtained by the classical Kudryashov method and modified Kudryashov method respectively. For this purpose, the modified Riemann-Liouville derivative is used to convert the nonlinear time fractional KdV-KZK equation into the nonlinear ordinary differential equation. In the present analysis, the classical Kudryashov method and modified Kudryashov method are both used successively to compute the analytical solutions of the time fractional KdV-KZK equation. As a result, new exact solutions involving the symmetrical Fibonacci function, hyperbolic function and exponential function are obtained for the first time. The methods under consideration are reliable and efficient, and can be used as an alternative to establish new exact solutions of different types of fractional differential equations arising from mathematical physics. The obtained results are exhibited graphically in order to demonstrate the efficiencies and applicabilities of these proposed methods of solving the nonlinear time fractional KdV-KZK equation.

Mean first passage times of Brownian rotators from differential recurrence relations

NASA Astrophysics Data System (ADS)

Coffey, W. T.

1999-11-01

An exact method of calculation of mean first passage times (analogous to that previously used [W. T. Coffey, Yu. P. Kalmykov, E. S. Massawe, and J. T. Waldron, J. Chem. Phys. 99, 4011 (1993)] for the correlation time) is developed in terms of continued fractions from the zero frequency limit of the Laplace transform of the set of differential recurrence relations generated by the Fokker-Planck or Langevin equations. The method because it is based on a Floquet representation avoids the use of quadratures and so may be easily generalized to multidegree of freedom systems by the use of matrix continued fractions. The procedure is illustrated by considering the mean first passage time of a fixed axis rotator with two equivalent sites.

Numerical approximations for fractional diffusion equations via a Chebyshev spectral-tau method

NASA Astrophysics Data System (ADS)

Doha, Eid H.; Bhrawy, Ali H.; Ezz-Eldien, Samer S.

2013-10-01

In this paper, a class of fractional diffusion equations with variable coefficients is considered. An accurate and efficient spectral tau technique for solving the fractional diffusion equations numerically is proposed. This method is based upon Chebyshev tau approximation together with Chebyshev operational matrix of Caputo fractional differentiation. Such approach has the advantage of reducing the problem to the solution of a system of algebraic equations, which may then be solved by any standard numerical technique. We apply this general method to solve four specific examples. In each of the examples considered, the numerical results show that the proposed method is of high accuracy and is efficient for solving the time-dependent fractional diffusion equations.

Power-law spatial dispersion from fractional Liouville equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tarasov, Vasily E.

2013-10-15

A microscopic model in the framework of fractional kinetics to describe spatial dispersion of power-law type is suggested. The Liouville equation with the Caputo fractional derivatives is used to obtain the power-law dependence of the absolute permittivity on the wave vector. The fractional differential equations for electrostatic potential in the media with power-law spatial dispersion are derived. The particular solutions of these equations for the electric potential of point charge in this media are considered.

Design of distributed PID-type dynamic matrix controller for fractional-order systems

NASA Astrophysics Data System (ADS)

Wang, Dawei; Zhang, Ridong

2018-01-01

With the continuous requirements for product quality and safety operation in industrial production, it is difficult to describe the complex large-scale processes with integer-order differential equations. However, the fractional differential equations may precisely represent the intrinsic characteristics of such systems. In this paper, a distributed PID-type dynamic matrix control method based on fractional-order systems is proposed. First, the high-order approximate model of integer order is obtained by utilising the Oustaloup method. Then, the step response model vectors of the plant is obtained on the basis of the high-order model, and the online optimisation for multivariable processes is transformed into the optimisation of each small-scale subsystem that is regarded as a sub-plant controlled in the distributed framework. Furthermore, the PID operator is introduced into the performance index of each subsystem and the fractional-order PID-type dynamic matrix controller is designed based on Nash optimisation strategy. The information exchange among the subsystems is realised through the distributed control structure so as to complete the optimisation task of the whole large-scale system. Finally, the control performance of the designed controller in this paper is verified by an example.

On the singular perturbations for fractional differential equation.

PubMed

Atangana, Abdon

2014-01-01

The goal of this paper is to examine the possible extension of the singular perturbation differential equation to the concept of fractional order derivative. To achieve this, we presented a review of the concept of fractional calculus. We make use of the Laplace transform operator to derive exact solution of singular perturbation fractional linear differential equations. We make use of the methodology of three analytical methods to present exact and approximate solution of the singular perturbation fractional, nonlinear, nonhomogeneous differential equation. These methods are including the regular perturbation method, the new development of the variational iteration method, and the homotopy decomposition method.

A comparison of artificial compressibility and fractional step methods for incompressible flow computations

NASA Technical Reports Server (NTRS)

Chan, Daniel C.; Darian, Armen; Sindir, Munir

1992-01-01

We have applied and compared the efficiency and accuracy of two commonly used numerical methods for the solution of Navier-Stokes equations. The artificial compressibility method augments the continuity equation with a transient pressure term and allows one to solve the modified equations as a coupled system. Due to its implicit nature, one can have the luxury of taking a large temporal integration step at the expense of higher memory requirement and larger operation counts per step. Meanwhile, the fractional step method splits the Navier-Stokes equations into a sequence of differential operators and integrates them in multiple steps. The memory requirement and operation count per time step are low, however, the restriction on the size of time marching step is more severe. To explore the strengths and weaknesses of these two methods, we used them for the computation of a two-dimensional driven cavity flow with Reynolds number of 100 and 1000, respectively. Three grid sizes, 41 x 41, 81 x 81, and 161 x 161 were used. The computations were considered after the L2-norm of the change of the dependent variables in two consecutive time steps has fallen below 10(exp -5).

Stochastic oscillations in models of epidemics on a network of cities

NASA Astrophysics Data System (ADS)

Rozhnova, G.; Nunes, A.; McKane, A. J.

2011-11-01

We carry out an analytic investigation of stochastic oscillations in a susceptible-infected-recovered model of disease spread on a network of n cities. In the model a fraction fjk of individuals from city k commute to city j, where they may infect, or be infected by, others. Starting from a continuous-time Markov description of the model the deterministic equations, which are valid in the limit when the population of each city is infinite, are recovered. The stochastic fluctuations about the fixed point of these equations are derived by use of the van Kampen system-size expansion. The fixed point structure of the deterministic equations is remarkably simple: A unique nontrivial fixed point always exists and has the feature that the fraction of susceptible, infected, and recovered individuals is the same for each city irrespective of its size. We find that the stochastic fluctuations have an analogously simple dynamics: All oscillations have a single frequency, equal to that found in the one-city case. We interpret this phenomenon in terms of the properties of the spectrum of the matrix of the linear approximation of the deterministic equations at the fixed point.

Bright and singular soliton solutions of the conformable time-fractional Klein-Gordon equations with different nonlinearities

NASA Astrophysics Data System (ADS)

Hosseini, Kamyar; Mayeli, Peyman; Ansari, Reza

2018-07-01

Finding the exact solutions of nonlinear fractional differential equations has gained considerable attention, during the past two decades. In this paper, the conformable time-fractional Klein-Gordon equations with quadratic and cubic nonlinearities are studied. Several exact soliton solutions, including the bright (non-topological) and singular soliton solutions are formally extracted by making use of the ansatz method. Results demonstrate that the method can efficiently handle the time-fractional Klein-Gordon equations with different nonlinearities.

A novel approach for solitary wave solutions of the generalized fractional Zakharov-Kuznetsov equation

NASA Astrophysics Data System (ADS)

Batool, Fiza; Akram, Ghazala

2018-01-01

In this article the solitary wave solutions of generalized fractional Zakharov-Kuznetsov (GZK) equation which appear in the electrical transmission line model are investigated. The (G'/G)-expansion method is used to obtain the solitary solutions of fractional GZK equation via local fractional derivative. Three classes of solutions, hyperbolic, trigonometric and rational wave solutions of the associated equation are characterized with some free parameters. The obtained solutions reveal that the proposed technique is effective and powerful.

Time-fractional characterization of brine reaction and precipitation in porous media

NASA Astrophysics Data System (ADS)

Xu, Jianping; Jiang, Guancheng

2018-04-01

Brine reaction and precipitation in porous media sometimes occur in the presence of a strong fluid flowing field, which induces the mobilization of the precipitated salts and distorts their spatial distribution. It is interesting to investigate how the distribution responds to such mobilization. We view these precipitates as random walkers in the complex inner space of the porous media, where they make stochastic jumps among locations and possibly wait between successive transitions. In consideration of related experimental results, the waiting time of the precipitates at a particular position is allowed to range widely from short sojourn to permanent residence. Through the model of a continuous-time random walk, a class of time-fractional equations for the precipitate's concentration profile is derived, including that in the Riemann-Liouville formalism and the Prabhakar formalism. The solutions to these equations show the general pattern of the precipitate's spatiotemporal evolution: a coupling of mass accumulation and mass transport. And the degree to which the mass is mobilized turns out to be monotonically correlated to the fractional exponent α . Moreover, to keep the completeness of the model, we further discuss how the interaction among the precipitates influences the precipitation process. In doing so, a time-fractional non-linear Fokker-Planck equation with source term is introduced and solved. It is shown that the interaction among the precipitates slightly perturbs their spatial distribution. This distribution is largely dominated by the brine reaction itself and the interaction between the precipitates and the porous media.

Application of the principal fractional meta-trigonometric functions for the solution of linear commensurate-order time-invariant fractional differential equations.

PubMed

Lorenzo, C F; Hartley, T T; Malti, R

2013-05-13

A new and simplified method for the solution of linear constant coefficient fractional differential equations of any commensurate order is presented. The solutions are based on the R-function and on specialized Laplace transform pairs derived from the principal fractional meta-trigonometric functions. The new method simplifies the solution of such fractional differential equations and presents the solutions in the form of real functions as opposed to fractional complex exponential functions, and thus is directly applicable to real-world physics.

On the Singular Perturbations for Fractional Differential Equation

PubMed Central

Atangana, Abdon

2014-01-01

The goal of this paper is to examine the possible extension of the singular perturbation differential equation to the concept of fractional order derivative. To achieve this, we presented a review of the concept of fractional calculus. We make use of the Laplace transform operator to derive exact solution of singular perturbation fractional linear differential equations. We make use of the methodology of three analytical methods to present exact and approximate solution of the singular perturbation fractional, nonlinear, nonhomogeneous differential equation. These methods are including the regular perturbation method, the new development of the variational iteration method, and the homotopy decomposition method. PMID:24683357

Discrete fractional solutions of a Legendre equation

NASA Astrophysics Data System (ADS)

Yılmazer, Resat

2018-01-01

One of the most popular research interests of science and engineering is the fractional calculus theory in recent times. Discrete fractional calculus has also an important position in fractional calculus. In this work, we acquire new discrete fractional solutions of the homogeneous and non homogeneous Legendre differential equation by using discrete fractional nabla operator.

Generalized Lie symmetry approach for fractional order systems of differential equations. III

NASA Astrophysics Data System (ADS)

Singla, Komal; Gupta, R. K.

2017-06-01

The generalized Lie symmetry technique is proposed for the derivation of point symmetries for systems of fractional differential equations with an arbitrary number of independent as well as dependent variables. The efficiency of the method is illustrated by its application to three higher dimensional nonlinear systems of fractional order partial differential equations consisting of the (2 + 1)-dimensional asymmetric Nizhnik-Novikov-Veselov system, (3 + 1)-dimensional Burgers system, and (3 + 1)-dimensional Navier-Stokes equations. With the help of derived Lie point symmetries, the corresponding invariant solutions transform each of the considered systems into a system of lower-dimensional fractional partial differential equations.

Dark Soliton Solutions of Space-Time Fractional Sharma-Tasso-Olver and Potential Kadomtsev-Petviashvili Equations

NASA Astrophysics Data System (ADS)

Guner, Ozkan; Korkmaz, Alper; Bekir, Ahmet

2017-02-01

Dark soliton solutions for space-time fractional Sharma-Tasso-Olver and space-time fractional potential Kadomtsev-Petviashvili equations are determined by using the properties of modified Riemann-Liouville derivative and fractional complex transform. After reducing both equations to nonlinear ODEs with constant coefficients, the \\tanh ansatz is substituted into the resultant nonlinear ODEs. The coefficients of the solutions in the ansatz are calculated by algebraic computer computations. Two different solutions are obtained for the Sharma-Tasso-Olver equation as only one solution for the potential Kadomtsev-Petviashvili equation. The solution profiles are demonstrated in 3D plots in finite domains of time and space.

Lie symmetry analysis, explicit solutions and conservation laws for the space-time fractional nonlinear evolution equations

NASA Astrophysics Data System (ADS)

Inc, Mustafa; Yusuf, Abdullahi; Aliyu, Aliyu Isa; Baleanu, Dumitru

2018-04-01

This paper studies the symmetry analysis, explicit solutions, convergence analysis, and conservation laws (Cls) for two different space-time fractional nonlinear evolution equations with Riemann-Liouville (RL) derivative. The governing equations are reduced to nonlinear ordinary differential equation (ODE) of fractional order using their Lie point symmetries. In the reduced equations, the derivative is in Erdelyi-Kober (EK) sense, power series technique is applied to derive an explicit solutions for the reduced fractional ODEs. The convergence of the obtained power series solutions is also presented. Moreover, the new conservation theorem and the generalization of the Noether operators are developed to construct the nonlocal Cls for the equations . Some interesting figures for the obtained explicit solutions are presented.

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Angstmann, C.N.; Donnelly, I.C.; Henry, B.I., E-mail: B.Henry@unsw.edu.au

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also showmore » that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.« less

General solution of a fractional Parker diffusion-convection equation describing the superdiffusive transport of energetic particles

NASA Astrophysics Data System (ADS)

Tawfik, Ashraf M.; Fichtner, Horst; Elhanbaly, A.; Schlickeiser, Reinhard

2018-06-01

Anomalous diffusion models of energetic particles in space plasmas are developed by introducing the fractional Parker diffusion-convection equation. Analytical solution of the space-time fractional equation is obtained by use of the Caputo and Riesz-Feller fractional derivatives with the Laplace-Fourier transforms. The solution is given in terms of the Fox H-function. Profiles of particle densities are illustrated for different values of the space fractional order and the so-called skewness parameter.

Characterization of non-diffusive transport in plasma turbulence by means of flux-gradient integro-differential kernels

NASA Astrophysics Data System (ADS)

Alcuson, J. A.; Reynolds-Barredo, J. M.; Mier, J. A.; Sanchez, Raul; Del-Castillo-Negrete, Diego; Newman, David E.; Tribaldos, V.

2015-11-01

A method to determine fractional transport exponents in systems dominated by fluid or plasma turbulence is proposed. The method is based on the estimation of the integro-differential kernel that relates values of the fluxes and gradients of the transported field, and its comparison with the family of analytical kernels of the linear fractional transport equation. Although use of this type of kernels has been explored before in this context, the methodology proposed here is rather unique since the connection with specific fractional equations is exploited from the start. The procedure has been designed to be particularly well-suited for application in experimental setups, taking advantage of the fact that kernel determination only requires temporal data of the transported field measured on an Eulerian grid. The simplicity and robustness of the method is tested first by using fabricated data from continuous-time random walk models built with prescribed transport characteristics. Its strengths are then illustrated on numerical Eulerian data gathered from simulations of a magnetically confined turbulent plasma in a near-critical regime, that is known to exhibit superdiffusive radial transport

NUMERICAL METHODS FOR SOLVING THE MULTI-TERM TIME-FRACTIONAL WAVE-DIFFUSION EQUATION.

PubMed

Liu, F; Meerschaert, M M; McGough, R J; Zhuang, P; Liu, Q

2013-03-01

In this paper, the multi-term time-fractional wave-diffusion equations are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], [1,2), [0,2), [0,3), [2,3) and [2,4), respectively. Some computationally effective numerical methods are proposed for simulating the multi-term time-fractional wave-diffusion equations. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and techniques can also be extended to other kinds of the multi-term fractional time-space models with fractional Laplacian.

NUMERICAL METHODS FOR SOLVING THE MULTI-TERM TIME-FRACTIONAL WAVE-DIFFUSION EQUATION

PubMed Central

Liu, F.; Meerschaert, M.M.; McGough, R.J.; Zhuang, P.; Liu, Q.

2013-01-01

In this paper, the multi-term time-fractional wave-diffusion equations are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], [1,2), [0,2), [0,3), [2,3) and [2,4), respectively. Some computationally effective numerical methods are proposed for simulating the multi-term time-fractional wave-diffusion equations. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and techniques can also be extended to other kinds of the multi-term fractional time-space models with fractional Laplacian. PMID:23772179

Heavy-tailed fractional Pearson diffusions.

PubMed

Leonenko, N N; Papić, I; Sikorskii, A; Šuvak, N

2017-11-01

We define heavy-tailed fractional reciprocal gamma and Fisher-Snedecor diffusions by a non-Markovian time change in the corresponding Pearson diffusions. Pearson diffusions are governed by the backward Kolmogorov equations with space-varying polynomial coefficients and are widely used in applications. The corresponding fractional reciprocal gamma and Fisher-Snedecor diffusions are governed by the fractional backward Kolmogorov equations and have heavy-tailed marginal distributions in the steady state. We derive the explicit expressions for the transition densities of the fractional reciprocal gamma and Fisher-Snedecor diffusions and strong solutions of the associated Cauchy problems for the fractional backward Kolmogorov equation.

A class of traveling wave solutions for space-time fractional biological population model in mathematical physics

NASA Astrophysics Data System (ADS)

Akram, Ghazala; Batool, Fiza

2017-10-01

The (G'/G)-expansion method is utilized for a reliable treatment of space-time fractional biological population model. The method has been applied in the sense of the Jumarie's modified Riemann-Liouville derivative. Three classes of exact traveling wave solutions, hyperbolic, trigonometric and rational solutions of the associated equation are characterized with some free parameters. A generalized fractional complex transform is applied to convert the fractional equations to ordinary differential equations which subsequently resulted in number of exact solutions. It should be mentioned that the (G'/G)-expansion method is very effective and convenient for solving nonlinear partial differential equations of fractional order whose balancing number is a negative integer.

Density-Dependent Conformable Space-time Fractional Diffusion-Reaction Equation and Its Exact Solutions

NASA Astrophysics Data System (ADS)

Hosseini, Kamyar; Mayeli, Peyman; Bekir, Ahmet; Guner, Ozkan

2018-01-01

In this article, a special type of fractional differential equations (FDEs) named the density-dependent conformable fractional diffusion-reaction (DDCFDR) equation is studied. Aforementioned equation has a significant role in the modelling of some phenomena arising in the applied science. The well-organized methods, including the \\exp (-φ (\\varepsilon )) -expansion and modified Kudryashov methods are exerted to generate the exact solutions of this equation such that some of the solutions are new and have been reported for the first time. Results illustrate that both methods have a great performance in handling the DDCFDR equation.

The numerical solution of linear multi-term fractional differential equations: systems of equations

NASA Astrophysics Data System (ADS)

Edwards, John T.; Ford, Neville J.; Simpson, A. Charles

2002-11-01

In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.

Numerical solution of distributed order fractional differential equations

NASA Astrophysics Data System (ADS)

Katsikadelis, John T.

2014-02-01

In this paper a method for the numerical solution of distributed order FDEs (fractional differential equations) of a general form is presented. The method applies to both linear and nonlinear equations. The Caputo type fractional derivative is employed. The distributed order FDE is approximated with a multi-term FDE, which is then solved by adjusting appropriately the numerical method developed for multi-term FDEs by Katsikadelis. Several example equations are solved and the response of mechanical systems described by such equations is studied. The convergence and the accuracy of the method for linear and nonlinear equations are demonstrated through well corroborated numerical results.

A multiphase approach to model ultrafiltration of deformable colloids

NASA Astrophysics Data System (ADS)

Haribabu, Malavika; Dunstan, Dave; Davidson, Malcolm; Harvie, Dalton

2017-11-01

Ultrafiltration (UF) is widely used in the dairy industry to fractionate and concentrate proteins, during the manufacture of milk protein concentrate and cheese. The protein build-up, comprising casein micelles (CM) and whey proteins, at the membrane surface during UF increases the resistance of the membrane system, thereby decreasing the performance of the process unit. CM have a complex structure that hydrodynamically behaves as a hard-sphere when dilute, but deforms beyond the random packing limit, forming a shear-thinning gel. This study employs a mixture model, based on the mixture phase continuity, Navier-Stokes equations, and solids continuity equation, to predict the solid concentration and velocity distribution during UF of CM. Micelle deformation is modelled as a function of volume fraction and dependent on its elastic modulus and particle size. The effect of deformation on gel permeability is implemented via Happel's permeability for hard spheres. Under crossflow conditions, the gel thickness is observed to increase along the membrane length, followed by a decrease towards the end of the membrane, resulting in an increase in flux at the latter section of the membrane. This study demonstrates that the membrane end-effects are important in determining UF performance.

On the solutions of fractional order of evolution equations

NASA Astrophysics Data System (ADS)

Morales-Delgado, V. F.; Taneco-Hernández, M. A.; Gómez-Aguilar, J. F.

2017-01-01

In this paper we present a discussion of generalized Cauchy problems in a diffusion wave process, we consider bi-fractional-order evolution equations in the Riemann-Liouville, Liouville-Caputo, and Caputo-Fabrizio sense. Through Fourier transforms and Laplace transform we derive closed-form solutions to the Cauchy problems mentioned above. Similarly, we establish fundamental solutions. Finally, we give an application of the above results to the determination of decompositions of Dirac type for bi-fractional-order equations and write a formula for the moments for the fractional vibration of a beam equation. This type of decomposition allows us to speak of internal degrees of freedom in the vibration of a beam equation.

Nonlocal electrical diffusion equation

NASA Astrophysics Data System (ADS)

Gómez-Aguilar, J. F.; Escobar-Jiménez, R. F.; Olivares-Peregrino, V. H.; Benavides-Cruz, M.; Calderón-Ramón, C.

2016-07-01

In this paper, we present an analysis and modeling of the electrical diffusion equation using the fractional calculus approach. This alternative representation for the current density is expressed in terms of the Caputo derivatives, the order for the space domain is 0<β≤1 and for the time domain is 0<γ≤2. We present solutions for the full fractional equation involving space and time fractional derivatives using numerical methods based on Fourier variable separation. The case with spatial fractional derivatives leads to Levy flight type phenomena, while the time fractional equation is related to sub- or super diffusion. We show that the mathematical concept of fractional derivatives can be useful to understand the behavior of semiconductors, the design of solar panels, electrochemical phenomena and the description of anomalous complex processes.

Stable multi-domain spectral penalty methods for fractional partial differential equations

NASA Astrophysics Data System (ADS)

Xu, Qinwu; Hesthaven, Jan S.

2014-01-01

We propose stable multi-domain spectral penalty methods suitable for solving fractional partial differential equations with fractional derivatives of any order. First, a high order discretization is proposed to approximate fractional derivatives of any order on any given grids based on orthogonal polynomials. The approximation order is analyzed and verified through numerical examples. Based on the discrete fractional derivative, we introduce stable multi-domain spectral penalty methods for solving fractional advection and diffusion equations. The equations are discretized in each sub-domain separately and the global schemes are obtained by weakly imposed boundary and interface conditions through a penalty term. Stability of the schemes are analyzed and numerical examples based on both uniform and nonuniform grids are considered to highlight the flexibility and high accuracy of the proposed schemes.

The effect of magnetohydrodynamic nano fluid flow through porous cylinder

NASA Astrophysics Data System (ADS)

Widodo, Basuki; Arif, Didik Khusnul; Aryany, Deviana; Asiyah, Nur; Widjajati, Farida Agustini; Kamiran

2017-08-01

This paper concerns about the analysis of the effect of magnetohydrodynamic nano fluid flow through horizontal porous cylinder on steady and incompressible condition. Fluid flow is assumed opposite gravity and induced by magnet field. Porous cylinder is assumed had the same depth of porous and was not absorptive. The First thing to do in this research is to build the model of fluid flow to obtain dimentional governing equations. The dimentional governing equations are consist of continuity equation, momentum equation, and energy equation. Furthermore, the dimensional governing equations are converted to non-dimensional governing equation by using non-dimensional parameters and variables. Then, the non-dimensional governing equations are transformed into similarity equations using stream function and solved using Keller-Box method. The result of numerical solution further is obtained by taking variation of magnetic parameter, Prandtl number, porosity parameter, and volume fraction. The numerical results show that velocity profiles increase and temperature profiles decrease when both of the magnetic and the porosity parameter increase. However, the velocity profiles decrease and the temperature profiles increase when both of the magnetic and the porosity parameter increase.

Two-time scale subordination in physical processes with long-term memory

NASA Astrophysics Data System (ADS)

Stanislavsky, Aleksander; Weron, Karina

2008-03-01

We describe dynamical processes in continuous media with a long-term memory. Our consideration is based on a stochastic subordination idea and concerns two physical examples in detail. First we study a temporal evolution of the species concentration in a trapping reaction in which a diffusing reactant is surrounded by a sea of randomly moving traps. The analysis uses the random-variable formalism of anomalous diffusive processes. We find that the empirical trapping-reaction law, according to which the reactant concentration decreases in time as a product of an exponential and a stretched exponential function, can be explained by a two-time scale subordination of random processes. Another example is connected with a state equation for continuous media with memory. If the pressure and the density of a medium are subordinated in two different random processes, then the ordinary state equation becomes fractional with two-time scales. This allows one to arrive at the Bagley-Torvik type of state equation.

Systems-based decomposition schemes for the approximate solution of multi-term fractional differential equations

NASA Astrophysics Data System (ADS)

Ford, Neville J.; Connolly, Joseph A.

2009-07-01

We give a comparison of the efficiency of three alternative decomposition schemes for the approximate solution of multi-term fractional differential equations using the Caputo form of the fractional derivative. The schemes we compare are based on conversion of the original problem into a system of equations. We review alternative approaches and consider how the most appropriate numerical scheme may be chosen to solve a particular equation.

Fractional diffusion on bounded domains

DOE PAGES

Defterli, Ozlem; D'Elia, Marta; Du, Qiang; ...

2015-03-13

We found that the mathematically correct specification of a fractional differential equation on a bounded domain requires specification of appropriate boundary conditions, or their fractional analogue. In this paper we discuss the application of nonlocal diffusion theory to specify well-posed fractional diffusion equations on bounded domains.

Boundary value problems for multi-term fractional differential equations

NASA Astrophysics Data System (ADS)

Daftardar-Gejji, Varsha; Bhalekar, Sachin

2008-09-01

Multi-term fractional diffusion-wave equation along with the homogeneous/non-homogeneous boundary conditions has been solved using the method of separation of variables. It is observed that, unlike in the one term case, solution of multi-term fractional diffusion-wave equation is not necessarily non-negative, and hence does not represent anomalous diffusion of any kind.

A new Jacobi spectral collocation method for solving 1+1 fractional Schrödinger equations and fractional coupled Schrödinger systems

NASA Astrophysics Data System (ADS)

Bhrawy, A. H.; Doha, E. H.; Ezz-Eldien, S. S.; Van Gorder, Robert A.

2014-12-01

The Jacobi spectral collocation method (JSCM) is constructed and used in combination with the operational matrix of fractional derivatives (described in the Caputo sense) for the numerical solution of the time-fractional Schrödinger equation (T-FSE) and the space-fractional Schrödinger equation (S-FSE). The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations, which greatly simplifies the solution process. In addition, the presented approach is also applied to solve the time-fractional coupled Schrödinger system (T-FCSS). In order to demonstrate the validity and accuracy of the numerical scheme proposed, several numerical examples with their approximate solutions are presented with comparisons between our numerical results and those obtained by other methods.

Investigation of the Dirac Equation by Using the Conformable Fractional Derivative

NASA Astrophysics Data System (ADS)

Mozaffari, F. S.; Hassanabadi, H.; Sobhani, H.; Chung, W. S.

2018-05-01

In this paper,the Dirac equation is constructed using the conformable fractional derivative so that in its limit for the fractional parameter, the normal version is recovered. Then, the Cornell potential is considered as the interaction of the system. In this case, the wave function and the energy eigenvalue equation are derived with the aim of the bi-confluent Heun functions. use of the conformable fractional derivative is proven to lead to a branching treatment for the energy of the system. Such a treatment is obvious for small values of the fractional parameter, and a united value as the fractional parameter approaches unity.

Fractional Stochastic Field Theory

NASA Astrophysics Data System (ADS)

Honkonen, Juha

2018-02-01

Models describing evolution of physical, chemical, biological, social and financial processes are often formulated as differential equations with the understanding that they are large-scale equations for averages of quantities describing intrinsically random processes. Explicit account of randomness may lead to significant changes in the asymptotic behaviour (anomalous scaling) in such models especially in low spatial dimensions, which in many cases may be captured with the use of the renormalization group. Anomalous scaling and memory effects may also be introduced with the use of fractional derivatives and fractional noise. Construction of renormalized stochastic field theory with fractional derivatives and fractional noise in the underlying stochastic differential equations and master equations and the interplay between fluctuation-induced and built-in anomalous scaling behaviour is reviewed and discussed.

Parallel Cartesian grid refinement for 3D complex flow simulations

NASA Astrophysics Data System (ADS)

Angelidis, Dionysios; Sotiropoulos, Fotis

2013-11-01

A second order accurate method for discretizing the Navier-Stokes equations on 3D unstructured Cartesian grids is presented. Although the grid generator is based on the oct-tree hierarchical method, fully unstructured data-structure is adopted enabling robust calculations for incompressible flows, avoiding both the need of synchronization of the solution between different levels of refinement and usage of prolongation/restriction operators. The current solver implements a hybrid staggered/non-staggered grid layout, employing the implicit fractional step method to satisfy the continuity equation. The pressure-Poisson equation is discretized by using a novel second order fully implicit scheme for unstructured Cartesian grids and solved using an efficient Krylov subspace solver. The momentum equation is also discretized with second order accuracy and the high performance Newton-Krylov method is used for integrating them in time. Neumann and Dirichlet conditions are used to validate the Poisson solver against analytical functions and grid refinement results to a significant reduction of the solution error. The effectiveness of the fractional step method results in the stability of the overall algorithm and enables the performance of accurate multi-resolution real life simulations. This material is based upon work supported by the Department of Energy under Award Number DE-EE0005482.

Fractional Poisson Fields and Martingales

NASA Astrophysics Data System (ADS)

Aletti, Giacomo; Leonenko, Nikolai; Merzbach, Ely

2018-02-01

We present new properties for the Fractional Poisson process (FPP) and the Fractional Poisson field on the plane. A martingale characterization for FPPs is given. We extend this result to Fractional Poisson fields, obtaining some other characterizations. The fractional differential equations are studied. We consider a more general Mixed-Fractional Poisson process and show that this process is the stochastic solution of a system of fractional differential-difference equations. Finally, we give some simulations of the Fractional Poisson field on the plane.

Fractional Diffusion Analysis of the Electromagnetic Field In Fractured Media Part II: 2.5-D Approach

NASA Astrophysics Data System (ADS)

Ge, J.; Everett, M. E.; Weiss, C. J.

2012-12-01

A 2.5D finite difference (FD) frequency-domain modeling algorithm based on the theory of fractional diffusion of electromagnetic (EM) fields generated by a loop source lying above a fractured geological medium is addressed in this paper. The presence of fractures in the subsurface, usually containing highly conductive pore fluids, gives rise to spatially hierarchical flow paths of induced EM eddy currents. The diffusion of EM eddy currents in such formations is anomalous, generalizing the classical Gaussian process described by the conventional Maxwell equations. Based on the continuous time random walk (CTRW) theory, the diffusion of EM eddy currents in a rough medium is governed by the fractional Maxwell equations. Here, we model the EM response of a 2D subsurface containing fractured zones, with a 3D loop source, which results the so-called 2.5D model geometry. The governing equation in the frequency domain is converted using Fourier transform into k domain along the strike direction (along which the model conductivity doesn't vary). The resulting equation system is solved by the multifrontal massively parallel solver (MUMPS). The data obtained is then converted back to spatial domain and the time domain. We find excellent agreement between the FD and analytic solutions for a rough halfspace model. Then FD solutions are calculated for a 2D fault zone model with variable conductivity and roughness. We compare the results with responses from several classical models and explore the relationship between the roughness and the spatial density of the fracture distribution.

Optimal Variational Asymptotic Method for Nonlinear Fractional Partial Differential Equations.

PubMed

Baranwal, Vipul K; Pandey, Ram K; Singh, Om P

2014-01-01

We propose optimal variational asymptotic method to solve time fractional nonlinear partial differential equations. In the proposed method, an arbitrary number of auxiliary parameters γ 0, γ 1, γ 2,… and auxiliary functions H 0(x), H 1(x), H 2(x),… are introduced in the correction functional of the standard variational iteration method. The optimal values of these parameters are obtained by minimizing the square residual error. To test the method, we apply it to solve two important classes of nonlinear partial differential equations: (1) the fractional advection-diffusion equation with nonlinear source term and (2) the fractional Swift-Hohenberg equation. Only few iterations are required to achieve fairly accurate solutions of both the first and second problems.

Fractional corresponding operator in quantum mechanics and applications: A uniform fractional Schrödinger equation in form and fractional quantization methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Xiao; Science and Technology on Electronic Information Control Laboratory, 610036, Chengdu, Sichuan; Wei, Chaozhen

2014-11-15

In this paper we use Dirac function to construct a fractional operator called fractional corresponding operator, which is the general form of momentum corresponding operator. Then we give a judging theorem for this operator and with this judging theorem we prove that R–L, G–L, Caputo, Riesz fractional derivative operator and fractional derivative operator based on generalized functions, which are the most popular ones, coincide with the fractional corresponding operator. As a typical application, we use the fractional corresponding operator to construct a new fractional quantization scheme and then derive a uniform fractional Schrödinger equation in form. Additionally, we find thatmore » the five forms of fractional Schrödinger equation belong to the particular cases. As another main result of this paper, we use fractional corresponding operator to generalize fractional quantization scheme by using Lévy path integral and use it to derive the corresponding general form of fractional Schrödinger equation, which consequently proves that these two quantization schemes are equivalent. Meanwhile, relations between the theory in fractional quantum mechanics and that in classic quantum mechanics are also discussed. As a physical example, we consider a particle in an infinite potential well. We give its wave functions and energy spectrums in two ways and find that both results are the same.« less

Two-dimensional simulation of GaAsSb/GaAs quantum dot solar cells

NASA Astrophysics Data System (ADS)

Kunrugsa, Maetee

2018-06-01

Two-dimensional (2D) simulation of GaAsSb/GaAs quantum dot (QD) solar cells is presented. The effects of As mole fraction in GaAsSb QDs on the performance of the solar cell are investigated. The solar cell is designed as a p-i-n GaAs structure where a single layer of GaAsSb QDs is introduced into the intrinsic region. The current density–voltage characteristics of QD solar cells are derived from Poisson’s equation, continuity equations, and the drift-diffusion transport equations, which are numerically solved by a finite element method. Furthermore, the transition energy of a single GaAsSb QD and its corresponding wavelength for each As mole fraction are calculated by a six-band k · p model to validate the position of the absorption edge in the external quantum efficiency curve. A GaAsSb/GaAs QD solar cell with an As mole fraction of 0.4 provides the best power conversion efficiency. The overlap between electron and hole wave functions becomes larger as the As mole fraction increases, leading to a higher optical absorption probability which is confirmed by the enhanced photogeneration rates within and around the QDs. However, further increasing the As mole fraction results in a reduction in the efficiency because the absorption edge moves towards shorter wavelengths, lowering the short-circuit current density. The influences of the QD size and density on the efficiency are also examined. For the GaAsSb/GaAs QD solar cell with an As mole fraction of 0.4, the efficiency can be improved to 26.2% by utilizing the optimum QD size and density. A decrease in the efficiency is observed at high QD densities, which is attributed to the increased carrier recombination and strain-modified band structures affecting the absorption edges.

Continuous time anomalous diffusion in a composite medium.

PubMed

Stickler, B A; Schachinger, E

2011-08-01

The one-dimensional continuous time anomalous diffusion in composite media consisting of a finite number of layers in immediate contact is investigated. The diffusion process itself is described with the help of two probability density functions (PDFs), one of which is an arbitrary jump-length PDF, and the other is a long-tailed waiting-time PDF characterized by the waiting-time index β∈(0,1). The former is assumed to be a function of the space coordinate x and the time coordinate t while the latter is a function of x and the time interval. For such an environment a very general form of the diffusion equation is derived which describes the continuous time anomalous diffusion in a composite medium. This result is then specialized to two particular forms of the jump-length PDF, namely the continuous time Lévy flight PDF and the continuous time truncated Lévy flight PDF. In both cases the PDFs are characterized by the Lévy index α∈(0,2) which is regarded to be a function of x and t. It is possible to demonstrate that for particular choices of the indices α and β other equations for anomalous diffusion, well known from the literature, follow immediately. This demonstrates the very general applicability of the derivation and of the resulting fractional differential equation discussed here.

Approximate series solution of multi-dimensional, time fractional-order (heat-like) diffusion equations using FRDTM.

PubMed

Singh, Brajesh K; Srivastava, Vineet K

2015-04-01

The main goal of this paper is to present a new approximate series solution of the multi-dimensional (heat-like) diffusion equation with time-fractional derivative in Caputo form using a semi-analytical approach: fractional-order reduced differential transform method (FRDTM). The efficiency of FRDTM is confirmed by considering four test problems of the multi-dimensional time fractional-order diffusion equation. FRDTM is a very efficient, effective and powerful mathematical tool which provides exact or very close approximate solutions for a wide range of real-world problems arising in engineering and natural sciences, modelled in terms of differential equations.

Approximate series solution of multi-dimensional, time fractional-order (heat-like) diffusion equations using FRDTM

PubMed Central

Singh, Brajesh K.; Srivastava, Vineet K.

2015-01-01

The main goal of this paper is to present a new approximate series solution of the multi-dimensional (heat-like) diffusion equation with time-fractional derivative in Caputo form using a semi-analytical approach: fractional-order reduced differential transform method (FRDTM). The efficiency of FRDTM is confirmed by considering four test problems of the multi-dimensional time fractional-order diffusion equation. FRDTM is a very efficient, effective and powerful mathematical tool which provides exact or very close approximate solutions for a wide range of real-world problems arising in engineering and natural sciences, modelled in terms of differential equations. PMID:26064639

Group iterative methods for the solution of two-dimensional time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Balasim, Alla Tareq; Ali, Norhashidah Hj. Mohd.

2016-06-01

Variety of problems in science and engineering may be described by fractional partial differential equations (FPDE) in relation to space and/or time fractional derivatives. The difference between time fractional diffusion equations and standard diffusion equations lies primarily in the time derivative. Over the last few years, iterative schemes derived from the rotated finite difference approximation have been proven to work well in solving standard diffusion equations. However, its application on time fractional diffusion counterpart is still yet to be investigated. In this paper, we will present a preliminary study on the formulation and analysis of new explicit group iterative methods in solving a two-dimensional time fractional diffusion equation. These methods were derived from the standard and rotated Crank-Nicolson difference approximation formula. Several numerical experiments were conducted to show the efficiency of the developed schemes in terms of CPU time and iteration number. At the request of all authors of the paper an updated version of this article was published on 7 July 2016. The original version supplied to AIP Publishing contained an error in Table 1 and References 15 and 16 were incomplete. These errors have been corrected in the updated and republished article.

Exact soliton of (2 + 1)-dimensional fractional Schrödinger equation

NASA Astrophysics Data System (ADS)

Rizvi, S. T. R.; Ali, K.; Bashir, S.; Younis, M.; Ashraf, R.; Ahmad, M. O.

2017-07-01

The nonlinear fractional Schrödinger equation is the basic equation of fractional quantum mechanics introduced by Nick Laskin in 2002. We apply three tools to solve this mathematical-physical model. First, we find the solitary wave solutions including the trigonometric traveling wave solutions, bell and kink shape solitons using the F-expansion and Improve F-expansion method. We also obtain the soliton solution, singular soliton solutions, rational function solution and elliptic integral function solutions, with the help of the extended trial equation method.

Calculation of spontaneous emission from a V-type three-level atom in photonic crystals using fractional calculus

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Chih-Hsien; Hsieh, Wen-Feng; Institute of Electro-Optical Science and Engineering, National Cheng Kung University, 1 Dahsueh Rd., Tainan 701, Taiwan

2011-07-15

Fractional time derivative, an abstract mathematical operator of fractional calculus, is used to describe the real optical system of a V-type three-level atom embedded in a photonic crystal. A fractional kinetic equation governing the dynamics of the spontaneous emission from this optical system is obtained as a fractional Langevin equation. Solving this fractional kinetic equation by fractional calculus leads to the analytical solutions expressed in terms of fractional exponential functions. The accuracy of the obtained solutions is verified through reducing the system into the special cases whose results are consistent with the experimental observation. With accurate physical results and avoidingmore » the complex integration for solving this optical system, we propose fractional calculus with fractional time derivative as a better mathematical method to study spontaneous emission dynamics from the optical system with non-Markovian dynamics.« less

Couple of the Variational Iteration Method and Fractional-Order Legendre Functions Method for Fractional Differential Equations

PubMed Central

Song, Junqiang; Leng, Hongze; Lu, Fengshun

2014-01-01

We present a new numerical method to get the approximate solutions of fractional differential equations. A new operational matrix of integration for fractional-order Legendre functions (FLFs) is first derived. Then a modified variational iteration formula which can avoid “noise terms” is constructed. Finally a numerical method based on variational iteration method (VIM) and FLFs is developed for fractional differential equations (FDEs). Block-pulse functions (BPFs) are used to calculate the FLFs coefficient matrices of the nonlinear terms. Five examples are discussed to demonstrate the validity and applicability of the technique. PMID:24511303

Analytical solutions for the motion of a charged particle in electric and magnetic fields via non-singular fractional derivatives

NASA Astrophysics Data System (ADS)

Morales-Delgado, V. F.; Gómez-Aguilar, J. F.; Taneco-Hernandez, M. A.

2017-12-01

In this work we propose fractional differential equations for the motion of a charged particle in electric, magnetic and electromagnetic fields. Exact solutions are obtained for the fractional differential equations by employing the Laplace transform method. The temporal fractional differential equations are considered in the Caputo-Fabrizio-Caputo and Atangana-Baleanu-Caputo sense. Application examples consider constant, ramp and harmonic fields. In addition, we present numerical results for different values of the fractional order. In all cases, when α = 1, we recover the standard electrodynamics.

Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus

NASA Astrophysics Data System (ADS)

Chen, Shanzhen; Jiang, Xiaoyun

2012-08-01

In this paper, analytical solutions to time-fractional partial differential equations in a multi-layer annulus are presented. The final solutions are obtained in terms of Mittag-Leffler function by using the finite integral transform technique and Laplace transform technique. In addition, the classical diffusion equation (α=1), the Helmholtz equation (α→0) and the wave equation (α=2) are discussed as special cases. Finally, an illustrative example problem for the three-layer semi-circular annular region is solved and numerical results are presented graphically for various kind of order of fractional derivative.

Boundary particle method for Laplace transformed time fractional diffusion equations

NASA Astrophysics Data System (ADS)

Fu, Zhuo-Jia; Chen, Wen; Yang, Hai-Tian

2013-02-01

This paper develops a novel boundary meshless approach, Laplace transformed boundary particle method (LTBPM), for numerical modeling of time fractional diffusion equations. It implements Laplace transform technique to obtain the corresponding time-independent inhomogeneous equation in Laplace space and then employs a truly boundary-only meshless boundary particle method (BPM) to solve this Laplace-transformed problem. Unlike the other boundary discretization methods, the BPM does not require any inner nodes, since the recursive composite multiple reciprocity technique (RC-MRM) is used to convert the inhomogeneous problem into the higher-order homogeneous problem. Finally, the Stehfest numerical inverse Laplace transform (NILT) is implemented to retrieve the numerical solutions of time fractional diffusion equations from the corresponding BPM solutions. In comparison with finite difference discretization, the LTBPM introduces Laplace transform and Stehfest NILT algorithm to deal with time fractional derivative term, which evades costly convolution integral calculation in time fractional derivation approximation and avoids the effect of time step on numerical accuracy and stability. Consequently, it can effectively simulate long time-history fractional diffusion systems. Error analysis and numerical experiments demonstrate that the present LTBPM is highly accurate and computationally efficient for 2D and 3D time fractional diffusion equations.

A remark on fractional differential equation involving I-function

NASA Astrophysics Data System (ADS)

Mishra, Jyoti

2018-02-01

The present paper deals with the solution of the fractional differential equation using the Laplace transform operator and its corresponding properties in the fractional calculus; we derive an exact solution of a complex fractional differential equation involving a special function known as I-function. The analysis of the some fractional integral with two parameters is presented using the suggested Theorem 1. In addition, some very useful corollaries are established and their proofs presented in detail. Some obtained exact solutions are depicted to see the effect of each fractional order. Owing to the wider applicability of the I-function, we can conclude that, the obtained results in our work generalize numerous well-known results obtained by specializing the parameters.

A finite element formulation preserving symmetric and banded diffusion stiffness matrix characteristics for fractional differential equations

NASA Astrophysics Data System (ADS)

Lin, Zeng; Wang, Dongdong

2017-10-01

Due to the nonlocal property of the fractional derivative, the finite element analysis of fractional diffusion equation often leads to a dense and non-symmetric stiffness matrix, in contrast to the conventional finite element formulation with a particularly desirable symmetric and banded stiffness matrix structure for the typical diffusion equation. This work first proposes a finite element formulation that preserves the symmetry and banded stiffness matrix characteristics for the fractional diffusion equation. The key point of the proposed formulation is the symmetric weak form construction through introducing a fractional weight function. It turns out that the stiffness part of the present formulation is identical to its counterpart of the finite element method for the conventional diffusion equation and thus the stiffness matrix formulation becomes trivial. Meanwhile, the fractional derivative effect in the discrete formulation is completely transferred to the force vector, which is obviously much easier and efficient to compute than the dense fractional derivative stiffness matrix. Subsequently, it is further shown that for the general fractional advection-diffusion-reaction equation, the symmetric and banded structure can also be maintained for the diffusion stiffness matrix, although the total stiffness matrix is not symmetric in this case. More importantly, it is demonstrated that under certain conditions this symmetric diffusion stiffness matrix formulation is capable of producing very favorable numerical solutions in comparison with the conventional non-symmetric diffusion stiffness matrix finite element formulation. The effectiveness of the proposed methodology is illustrated through a series of numerical examples.

Variational iteration method — a promising technique for constructing equivalent integral equations of fractional order

NASA Astrophysics Data System (ADS)

Wang, Yi-Hong; Wu, Guo-Cheng; Baleanu, Dumitru

2013-10-01

The variational iteration method is newly used to construct various integral equations of fractional order. Some iterative schemes are proposed which fully use the method and the predictor-corrector approach. The fractional Bagley-Torvik equation is then illustrated as an example of multi-order and the results show the efficiency of the variational iteration method's new role.

A 2D multi-term time and space fractional Bloch-Torrey model based on bilinear rectangular finite elements

NASA Astrophysics Data System (ADS)

Qin, Shanlin; Liu, Fawang; Turner, Ian W.

2018-03-01

The consideration of diffusion processes in magnetic resonance imaging (MRI) signal attenuation is classically described by the Bloch-Torrey equation. However, many recent works highlight the distinct deviation in MRI signal decay due to anomalous diffusion, which motivates the fractional order generalization of the Bloch-Torrey equation. In this work, we study the two-dimensional multi-term time and space fractional diffusion equation generalized from the time and space fractional Bloch-Torrey equation. By using the Galerkin finite element method with a structured mesh consisting of rectangular elements to discretize in space and the L1 approximation of the Caputo fractional derivative in time, a fully discrete numerical scheme is derived. A rigorous analysis of stability and error estimation is provided. Numerical experiments in the square and L-shaped domains are performed to give an insight into the efficiency and reliability of our method. Then the scheme is applied to solve the multi-term time and space fractional Bloch-Torrey equation, which shows that the extra time derivative terms impact the relaxation process.

Material point method modeling in oil and gas reservoirs

DOEpatents

Vanderheyden, William Brian; Zhang, Duan

2016-06-28

A computer system and method of simulating the behavior of an oil and gas reservoir including changes in the margins of frangible solids. A system of equations including state equations such as momentum, and conservation laws such as mass conservation and volume fraction continuity, are defined and discretized for at least two phases in a modeled volume, one of which corresponds to frangible material. A material point model technique for numerically solving the system of discretized equations, to derive fluid flow at each of a plurality of mesh nodes in the modeled volume, and the velocity of at each of a plurality of particles representing the frangible material in the modeled volume. A time-splitting technique improves the computational efficiency of the simulation while maintaining accuracy on the deformation scale. The method can be applied to derive accurate upscaled model equations for larger volume scale simulations.

Properties of wavelet discretization of Black-Scholes equation

NASA Astrophysics Data System (ADS)

Finěk, Václav

2017-07-01

Using wavelet methods, the continuous problem is transformed into a well-conditioned discrete problem. And once a non-symmetric problem is given, squaring yields a symmetric positive definite formulation. However squaring usually makes the condition number of discrete problems substantially worse. This note is concerned with a wavelet based numerical solution of the Black-Scholes equation for pricing European options. We show here that in wavelet coordinates a symmetric part of the discretized equation dominates over an unsymmetric part in the standard economic environment with low interest rates. It provides some justification for using a fractional step method with implicit treatment of the symmetric part of the weak form of the Black-Scholes operator and with explicit treatment of its unsymmetric part. Then a well-conditioned discrete problem is obtained.

Fractional quantum mechanics on networks: Long-range dynamics and quantum transport

NASA Astrophysics Data System (ADS)

Riascos, A. P.; Mateos, José L.

2015-11-01

In this paper we study the quantum transport on networks with a temporal evolution governed by the fractional Schrödinger equation. We generalize the dynamics based on continuous-time quantum walks, with transitions to nearest neighbors on the network, to the fractional case that allows long-range displacements. By using the fractional Laplacian matrix of a network, we establish a formalism that combines a long-range dynamics with the quantum superposition of states; this general approach applies to any type of connected undirected networks, including regular, random, and complex networks, and can be implemented from the spectral properties of the Laplacian matrix. We study the fractional dynamics and its capacity to explore the network by means of the transition probability, the average probability of return, and global quantities that characterize the efficiency of this quantum process. As a particular case, we explore analytically these quantities for circulant networks such as rings, interacting cycles, and complete graphs.

Fractional Poisson-Nernst-Planck Model for Ion Channels I: Basic Formulations and Algorithms.

PubMed

Chen, Duan

2017-11-01

In this work, we propose a fractional Poisson-Nernst-Planck model to describe ion permeation in gated ion channels. Due to the intrinsic conformational changes, crowdedness in narrow channel pores, binding and trapping introduced by functioning units of channel proteins, ionic transport in the channel exhibits a power-law-like anomalous diffusion dynamics. We start from continuous-time random walk model for a single ion and use a long-tailed density distribution function for the particle jump waiting time, to derive the fractional Fokker-Planck equation. Then, it is generalized to the macroscopic fractional Poisson-Nernst-Planck model for ionic concentrations. Necessary computational algorithms are designed to implement numerical simulations for the proposed model, and the dynamics of gating current is investigated. Numerical simulations show that the fractional PNP model provides a more qualitatively reasonable match to the profile of gating currents from experimental observations. Meanwhile, the proposed model motivates new challenges in terms of mathematical modeling and computations.

Fractional quantum mechanics on networks: Long-range dynamics and quantum transport.

PubMed

Riascos, A P; Mateos, José L

2015-11-01

In this paper we study the quantum transport on networks with a temporal evolution governed by the fractional Schrödinger equation. We generalize the dynamics based on continuous-time quantum walks, with transitions to nearest neighbors on the network, to the fractional case that allows long-range displacements. By using the fractional Laplacian matrix of a network, we establish a formalism that combines a long-range dynamics with the quantum superposition of states; this general approach applies to any type of connected undirected networks, including regular, random, and complex networks, and can be implemented from the spectral properties of the Laplacian matrix. We study the fractional dynamics and its capacity to explore the network by means of the transition probability, the average probability of return, and global quantities that characterize the efficiency of this quantum process. As a particular case, we explore analytically these quantities for circulant networks such as rings, interacting cycles, and complete graphs.

Numerical analysis for the fractional diffusion and fractional Buckmaster equation by the two-step Laplace Adam-Bashforth method

NASA Astrophysics Data System (ADS)

Jain, Sonal

2018-01-01

In this paper, we aim to use the alternative numerical scheme given by Gnitchogna and Atangana for solving partial differential equations with integer and non-integer differential operators. We applied this method to fractional diffusion model and fractional Buckmaster models with non-local fading memory. The method yields a powerful numerical algorithm for fractional order derivative to implement. Also we present in detail the stability analysis of the numerical method for solving the diffusion equation. This proof shows that this method is very stable and also converges very quickly to exact solution and finally some numerical simulation is presented.

Non-standard finite difference and Chebyshev collocation methods for solving fractional diffusion equation

NASA Astrophysics Data System (ADS)

Agarwal, P.; El-Sayed, A. A.

2018-06-01

In this paper, a new numerical technique for solving the fractional order diffusion equation is introduced. This technique basically depends on the Non-Standard finite difference method (NSFD) and Chebyshev collocation method, where the fractional derivatives are described in terms of the Caputo sense. The Chebyshev collocation method with the (NSFD) method is used to convert the problem into a system of algebraic equations. These equations solved numerically using Newton's iteration method. The applicability, reliability, and efficiency of the presented technique are demonstrated through some given numerical examples.

Imaging of the internal structure of comet 67P/Churyumov-Gerasimenko from radiotomography CONSERT Data (Rosetta Mission) through a full 3D regularized inversion of the Helmholtz equations on functional spaces

NASA Astrophysics Data System (ADS)

Barriot, Jean-Pierre; Serafini, Jonathan; Sichoix, Lydie; Benna, Mehdi; Kofman, Wlodek; Herique, Alain

We investigate the inverse problem of imaging the internal structure of comet 67P/ Churyumov-Gerasimenko from radiotomography CONSERT data by using a coupled regularized inversion of the Helmholtz equations. A first set of Helmholtz equations, written w.r.t a basis of 3D Hankel functions describes the wave propagation outside the comet at large distances, a second set of Helmholtz equations, written w.r.t. a basis of 3D Zernike functions describes the wave propagation throughout the comet with avariable permittivity. Both sets are connected by continuity equations over a sphere that surrounds the comet. This approach, derived from GPS water vapor tomography of the atmosphere,will permit a full 3D inversion of the internal structure of the comet, contrary to traditional approaches that use a discretization of space at a fraction of the radiowave wavelength.

FRACTIONAL PEARSON DIFFUSIONS.

PubMed

Leonenko, Nikolai N; Meerschaert, Mark M; Sikorskii, Alla

2013-07-15

Pearson diffusions are governed by diffusion equations with polynomial coefficients. Fractional Pearson diffusions are governed by the corresponding time-fractional diffusion equation. They are useful for modeling sub-diffusive phenomena, caused by particle sticking and trapping. This paper provides explicit strong solutions for fractional Pearson diffusions, using spectral methods. It also presents stochastic solutions, using a non-Markovian inverse stable time change.

New solitary wave solutions of the time-fractional Cahn-Allen equation via the improved (G'/G)-expansion method

NASA Astrophysics Data System (ADS)

Batool, Fiza; Akram, Ghazala

2018-05-01

An improved (G'/G)-expansion method is proposed for extracting more general solitary wave solutions of the nonlinear fractional Cahn-Allen equation. The temporal fractional derivative is taken in the sense of Jumarie's fractional derivative. The results of this article are generalized and extended version of previously reported solutions.

Analytical Solutions, Moments, and Their Asymptotic Behaviors for the Time-Space Fractional Cable Equation

NASA Astrophysics Data System (ADS)

Li, Can; Deng, Wei-Hua

2014-07-01

Following the fractional cable equation established in the letter [B.I. Henry, T.A.M. Langlands, and S.L. Wearne, Phys. Rev. Lett. 100 (2008) 128103], we present the time-space fractional cable equation which describes the anomalous transport of electrodiffusion in nerve cells. The derivation is based on the generalized fractional Ohm's law; and the temporal memory effects and spatial-nonlocality are involved in the time-space fractional model. With the help of integral transform method we derive the analytical solutions expressed by the Green's function; the corresponding fractional moments are calculated; and their asymptotic behaviors are discussed. In addition, the explicit solutions of the considered model with two different external current injections are also presented.

Plasmon dispersion in a multilayer solid torus in terms of three-term vector recurrence relations and matrix continued fractions

DOE PAGES

Garapati, K. V.; Bagherian, M.; Passian, A.; ...

2018-01-03

Toroidal confinement, which has played a crucial role in magnetized plasmas and Tokamak physics, is emerging as an effective means to obtain useful electronic and optical response in solids. In particular, excitation of surface plasmons in metal nanorings by photons or electrons finds important applications due to the engendered field distribution and electromagnetic energy confinement. However, in contrast to the case of a plasma, often the solid nanorings are multilayered and/or embedded in a medium. The non-simply connected geometry of the torus results in surface modes that are not linearly independent. A three-term difference equation was recently shown to arisemore » when seeking the nonretarded plasmon dispersion relations for a stratified solid torus (Garapati et al 2017 Phys. Rev. B 95 165422). The reported generalized plasmon dispersion relations are here investigated in terms of the involved matrix continued fractions and their convergence properties including the determinant forms of the dispersion relations obtained for computing the plasmon eigenmodes. We also present the intricacies of the derivation and properties of the Green's function employed to solve the three term amplitude equation that determines the response of the toroidal structure to arbitrary external excitations.« less

Plasmon dispersion in a multilayer solid torus in terms of three-term vector recurrence relations and matrix continued fractions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Garapati, K. V.; Bagherian, M.; Passian, A.

Toroidal confinement, which has played a crucial role in magnetized plasmas and Tokamak physics, is emerging as an effective means to obtain useful electronic and optical response in solids. In particular, excitation of surface plasmons in metal nanorings by photons or electrons finds important applications due to the engendered field distribution and electromagnetic energy confinement. However, in contrast to the case of a plasma, often the solid nanorings are multilayered and/or embedded in a medium. The non-simply connected geometry of the torus results in surface modes that are not linearly independent. A three-term difference equation was recently shown to arisemore » when seeking the nonretarded plasmon dispersion relations for a stratified solid torus (Garapati et al 2017 Phys. Rev. B 95 165422). The reported generalized plasmon dispersion relations are here investigated in terms of the involved matrix continued fractions and their convergence properties including the determinant forms of the dispersion relations obtained for computing the plasmon eigenmodes. We also present the intricacies of the derivation and properties of the Green's function employed to solve the three term amplitude equation that determines the response of the toroidal structure to arbitrary external excitations.« less

Gaussian solitary waves and compactons in Fermi–Pasta–Ulam lattices with Hertzian potentials

PubMed Central

James, Guillaume; Pelinovsky, Dmitry

2014-01-01

We consider a class of fully nonlinear Fermi–Pasta–Ulam (FPU) lattices, consisting of a chain of particles coupled by fractional power nonlinearities of order α>1. This class of systems incorporates a classical Hertzian model describing acoustic wave propagation in chains of touching beads in the absence of precompression. We analyse the propagation of localized waves when α is close to unity. Solutions varying slowly in space and time are searched with an appropriate scaling, and two asymptotic models of the chain of particles are derived consistently. The first one is a logarithmic Korteweg–de Vries (KdV) equation and possesses linearly orbitally stable Gaussian solitary wave solutions. The second model consists of a generalized KdV equation with Hölder-continuous fractional power nonlinearity and admits compacton solutions, i.e. solitary waves with compact support. When , we numerically establish the asymptotically Gaussian shape of exact FPU solitary waves with near-sonic speed and analytically check the pointwise convergence of compactons towards the limiting Gaussian profile. PMID:24808748

Relational Priming Based on a Multiplicative Schema for Whole Numbers and Fractions.

PubMed

DeWolf, Melissa; Son, Ji Y; Bassok, Miriam; Holyoak, Keith J

2017-11-01

Why might it be (at least sometimes) beneficial for adults to process fractions componentially? Recent research has shown that college-educated adults can capitalize on the bipartite structure of the fraction notation, performing more successfully with fractions than with decimals in relational tasks, notably analogical reasoning. This study examined patterns of relational priming for problems with fractions in a task that required arithmetic computations. College students were asked to judge whether or not multiplication equations involving fractions were correct. Some equations served as structurally inverse primes for the equation that immediately followed it (e.g., 4 × 3/4 = 3 followed by 3 × 8/6 = 4). Students with relatively high math ability showed relational priming (speeded solution times to the second of two successive relationally related fraction equations) both with and without high perceptual similarity (Experiment 2). Students with relatively low math ability also showed priming, but only when the structurally inverse equation pairs were supported by high perceptual similarity between numbers (e.g., 4 × 3/4 = 3 followed by 3 × 4/3 = 4). Several additional experiments established boundary conditions on relational priming with fractions. These findings are interpreted in terms of componential processing of fractions in a relational multiplication context that takes advantage of their inherent connections to a multiplicative schema for whole numbers. Copyright © 2017 Cognitive Science Society, Inc.

Accessible solitons of fractional dimension

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhong, Wei-Ping, E-mail: zhongwp6@126.com; Texas A&M University at Qatar, P.O. Box 23874, Doha; Belić, Milivoj

We demonstrate that accessible solitons described by an extended Schrödinger equation with the Laplacian of fractional dimension can exist in strongly nonlocal nonlinear media. The soliton solutions of the model are constructed by two special functions, the associated Legendre polynomials and the Laguerre polynomials in the fraction-dimensional space. Our results show that these fractional accessible solitons form a soliton family which includes crescent solitons, and asymmetric single-layer and multi-layer necklace solitons. -- Highlights: •Analytic solutions of a fractional Schrödinger equation are obtained. •The solutions are produced by means of self-similar method applied to the fractional Schrödinger equation with parabolic potential.more » •The fractional accessible solitons form crescent, asymmetric single-layer and multilayer necklace profiles. •The model applies to the propagation of optical pulses in strongly nonlocal nonlinear media.« less

Fluorescence correlation spectroscopy: the case of subdiffusion.

PubMed

Lubelski, Ariel; Klafter, Joseph

2009-03-18

The theory of fluorescence correlation spectroscopy is revisited here for the case of subdiffusing molecules. Subdiffusion is assumed to stem from a continuous-time random walk process with a fat-tailed distribution of waiting times and can therefore be formulated in terms of a fractional diffusion equation (FDE). The FDE plays the central role in developing the fluorescence correlation spectroscopy expressions, analogous to the role played by the simple diffusion equation for regular systems. Due to the nonstationary nature of the continuous-time random walk/FDE, some interesting properties emerge that are amenable to experimental verification and may help in discriminating among subdiffusion mechanisms. In particular, the current approach predicts 1), a strong dependence of correlation functions on the initial time (aging); 2), sensitivity of correlation functions to the averaging procedure, ensemble versus time averaging (ergodicity breaking); and 3), that the basic mean-squared displacement observable depends on how the mean is taken.

Transient aging in fractional Brownian and Langevin-equation motion.

PubMed

Kursawe, Jochen; Schulz, Johannes; Metzler, Ralf

2013-12-01

Stochastic processes driven by stationary fractional Gaussian noise, that is, fractional Brownian motion and fractional Langevin-equation motion, are usually considered to be ergodic in the sense that, after an algebraic relaxation, time and ensemble averages of physical observables coincide. Recently it was demonstrated that fractional Brownian motion and fractional Langevin-equation motion under external confinement are transiently nonergodic-time and ensemble averages behave differently-from the moment when the particle starts to sense the confinement. Here we show that these processes also exhibit transient aging, that is, physical observables such as the time-averaged mean-squared displacement depend on the time lag between the initiation of the system at time t=0 and the start of the measurement at the aging time t(a). In particular, it turns out that for fractional Langevin-equation motion the aging dependence on t(a) is different between the cases of free and confined motion. We obtain explicit analytical expressions for the aged moments of the particle position as well as the time-averaged mean-squared displacement and present a numerical analysis of this transient aging phenomenon.

Fractional analysis for nonlinear electrical transmission line and nonlinear Schroedinger equations with incomplete sub-equation

NASA Astrophysics Data System (ADS)

Fendzi-Donfack, Emmanuel; Nguenang, Jean Pierre; Nana, Laurent

2018-02-01

We use the fractional complex transform with the modified Riemann-Liouville derivative operator to establish the exact and generalized solutions of two fractional partial differential equations. We determine the solutions of fractional nonlinear electrical transmission lines (NETL) and the perturbed nonlinear Schroedinger (NLS) equation with the Kerr law nonlinearity term. The solutions are obtained for the parameters in the range (0<α≤1) of the derivative operator and we found the traditional solutions for the limiting case of α =1. We show that according to the modified Riemann-Liouville derivative, the solutions found can describe physical systems with memory effect, transient effects in electrical systems and nonlinear transmission lines, and other systems such as optical fiber.

Dynamic stiffness of chemically and physically ageing rubber vibration isolators in the audible frequency range. Part 1: constitutive equations

NASA Astrophysics Data System (ADS)

Kari, Leif

2017-09-01

The constitutive equations of chemically and physically ageing rubber in the audible frequency range are modelled as a function of ageing temperature, ageing time, actual temperature, time and frequency. The constitutive equations are derived by assuming nearly incompressible material with elastic spherical response and viscoelastic deviatoric response, using Mittag-Leffler relaxation function of fractional derivative type, the main advantage being the minimum material parameters needed to successfully fit experimental data over a broad frequency range. The material is furthermore assumed essentially entropic and thermo-mechanically simple while using a modified William-Landel-Ferry shift function to take into account temperature dependence and physical ageing, with fractional free volume evolution modelled by a nonlinear, fractional differential equation with relaxation time identical to that of the stress response and related to the fractional free volume by Doolittle equation. Physical ageing is a reversible ageing process, including trapping and freeing of polymer chain ends, polymer chain reorganizations and free volume changes. In contrast, chemical ageing is an irreversible process, mainly attributed to oxygen reaction with polymer network either damaging the network by scission or reformation of new polymer links. The chemical ageing is modelled by inner variables that are determined by inner fractional evolution equations. Finally, the model parameters are fitted to measurements results of natural rubber over a broad audible frequency range, and various parameter studies are performed including comparison with results obtained by ordinary, non-fractional ageing evolution differential equations.

Khater method for nonlinear Sharma Tasso-Olever (STO) equation of fractional order

NASA Astrophysics Data System (ADS)

Bibi, Sadaf; Mohyud-Din, Syed Tauseef; Khan, Umar; Ahmed, Naveed

In this work, we have implemented a direct method, known as Khater method to establish exact solutions of nonlinear partial differential equations of fractional order. Number of solutions provided by this method is greater than other traditional methods. Exact solutions of nonlinear fractional order Sharma Tasso-Olever (STO) equation are expressed in terms of kink, travelling wave, periodic and solitary wave solutions. Modified Riemann-Liouville derivative and Fractional complex transform have been used for compatibility with fractional order sense. Solutions have been graphically simulated for understanding the physical aspects and importance of the method. A comparative discussion between our established results and the results obtained by existing ones is also presented. Our results clearly reveal that the proposed method is an effective, powerful and straightforward technique to work out new solutions of various types of differential equations of non-integer order in the fields of applied sciences and engineering.

Stochastic Evolution Equations Driven by Fractional Noises

DTIC Science & Technology

2016-11-28

rate of convergence to zero or the error and the limit in distribution of the error fluctuations. We have studied time discrete numerical schemes...error fluctuations. We have studied time discrete numerical schemes based on Taylor expansions for rough differential equations and for stochastic...variations of the time discrete Taylor schemes for rough differential equations and for stochastic differential equations driven by fractional Brownian

Solving Nonlinear Fractional Differential Equation by Generalized Mittag-Leffler Function Method

NASA Astrophysics Data System (ADS)

Arafa, A. A. M.; Rida, S. Z.; Mohammadein, A. A.; Ali, H. M.

2013-06-01

In this paper, we use Mittag—Leffler function method for solving some nonlinear fractional differential equations. A new solution is constructed in power series. The fractional derivatives are described by Caputo's sense. To illustrate the reliability of the method, some examples are provided.

Size-dependent geometrically nonlinear free vibration analysis of fractional viscoelastic nanobeams based on the nonlocal elasticity theory

NASA Astrophysics Data System (ADS)

Ansari, R.; Faraji Oskouie, M.; Gholami, R.

2016-01-01

In recent decades, mathematical modeling and engineering applications of fractional-order calculus have been extensively utilized to provide efficient simulation tools in the field of solid mechanics. In this paper, a nonlinear fractional nonlocal Euler-Bernoulli beam model is established using the concept of fractional derivative and nonlocal elasticity theory to investigate the size-dependent geometrically nonlinear free vibration of fractional viscoelastic nanobeams. The non-classical fractional integro-differential Euler-Bernoulli beam model contains the nonlocal parameter, viscoelasticity coefficient and order of the fractional derivative to interpret the size effect, viscoelastic material and fractional behavior in the nanoscale fractional viscoelastic structures, respectively. In the solution procedure, the Galerkin method is employed to reduce the fractional integro-partial differential governing equation to a fractional ordinary differential equation in the time domain. Afterwards, the predictor-corrector method is used to solve the nonlinear fractional time-dependent equation. Finally, the influences of nonlocal parameter, order of fractional derivative and viscoelasticity coefficient on the nonlinear time response of fractional viscoelastic nanobeams are discussed in detail. Moreover, comparisons are made between the time responses of linear and nonlinear models.

Analytical solution of the nonlinear diffusion equation

NASA Astrophysics Data System (ADS)

Shanker Dubey, Ravi; Goswami, Pranay

2018-05-01

In the present paper, we derive the solution of the nonlinear fractional partial differential equations using an efficient approach based on the q -homotopy analysis transform method ( q -HATM). The fractional diffusion equations derivatives are considered in Caputo sense. The derived results are graphically demonstrated as well.

Modeling transport across the running-sandpile cellular automaton by means of fractional transport equations

NASA Astrophysics Data System (ADS)

Sánchez, R.; Newman, D. E.; Mier, J. A.

2018-05-01

Fractional transport equations are used to build an effective model for transport across the running sandpile cellular automaton [Hwa et al., Phys. Rev. A 45, 7002 (1992), 10.1103/PhysRevA.45.7002]. It is shown that both temporal and spatial fractional derivatives must be considered to properly reproduce the sandpile transport features, which are governed by self-organized criticality, at least over sufficiently long or large scales. In contrast to previous applications of fractional transport equations to other systems, the specifics of sand motion require in this case that the spatial fractional derivatives used for the running sandpile must be of the completely asymmetrical Riesz-Feller type. Appropriate values for the fractional exponents that define these derivatives in the case of the running sandpile are obtained numerically.

Exact traveling wave solutions of fractional order Boussinesq-like equations by applying Exp-function method

NASA Astrophysics Data System (ADS)

Rahmatullah; Ellahi, Rahmat; Mohyud-Din, Syed Tauseef; Khan, Umar

2018-03-01

We have computed new exact traveling wave solutions, including complex solutions of fractional order Boussinesq-Like equations, occurring in physical sciences and engineering, by applying Exp-function method. The method is blended with fractional complex transformation and modified Riemann-Liouville fractional order operator. Our obtained solutions are verified by substituting back into their corresponding equations. To the best of our knowledge, no other technique has been reported to cope with the said fractional order nonlinear problems combined with variety of exact solutions. Graphically, fractional order solution curves are shown to be strongly related to each other and most importantly, tend to fixate on their integer order solution curve. Our solutions comprise high frequencies and very small amplitude of the wave responses.

Fractional dynamics pharmacokinetics–pharmacodynamic models

PubMed Central

2010-01-01

While an increasing number of fractional order integrals and differential equations applications have been reported in the physics, signal processing, engineering and bioengineering literatures, little attention has been paid to this class of models in the pharmacokinetics–pharmacodynamic (PKPD) literature. One of the reasons is computational: while the analytical solution of fractional differential equations is available in special cases, it this turns out that even the simplest PKPD models that can be constructed using fractional calculus do not allow an analytical solution. In this paper, we first introduce new families of PKPD models incorporating fractional order integrals and differential equations, and, second, exemplify and investigate their qualitative behavior. The families represent extensions of frequently used PK link and PD direct and indirect action models, using the tools of fractional calculus. In addition the PD models can be a function of a variable, the active drug, which can smoothly transition from concentration to exposure, to hyper-exposure, according to a fractional integral transformation. To investigate the behavior of the models we propose, we implement numerical algorithms for fractional integration and for the numerical solution of a system of fractional differential equations. For simplicity, in our investigation we concentrate on the pharmacodynamic side of the models, assuming standard (integer order) pharmacokinetics. PMID:20455076

A study of fractional Schrödinger equation composed of Jumarie fractional derivative

NASA Astrophysics Data System (ADS)

Banerjee, Joydip; Ghosh, Uttam; Sarkar, Susmita; Das, Shantanu

2017-04-01

In this paper we have derived the fractional-order Schrödinger equation composed of Jumarie fractional derivative. The solution of this fractional-order Schrödinger equation is obtained in terms of Mittag-Leffler function with complex arguments, and fractional trigonometric functions. A few important properties of the fractional Schrödinger equation are then described for the case of particles in one-dimensional infinite potential well. One of the motivations for using fractional calculus in physical systems is that the space and time variables, which we often deal with, exhibit coarse-grained phenomena. This means infinitesimal quantities cannot be arbitrarily taken to zero - rather they are non-zero with a minimum spread. This type of non-zero spread arises in the microscopic to mesoscopic levels of system dynamics, which means that, if we denote x as the point in space and t as the point in time, then limit of the differentials d x (and d t) cannot be taken as zero. To take the concept of coarse graining into account, use the infinitesimal quantities as (Δ x) α (and (Δ t) α ) with 0 < α < 1; called as `fractional differentials'. For arbitrarily small Δ x and Δ t (tending towards zero), these `fractional' differentials are greater than Δ x (and Δ t), i.e. (Δ x) α > Δ x and (Δ t) α > Δ t. This way of defining the fractional differentials helps us to use fractional derivatives in the study of dynamic systems.

On exact traveling-wave solutions for local fractional Korteweg-de Vries equation.

PubMed

Yang, Xiao-Jun; Tenreiro Machado, J A; Baleanu, Dumitru; Cattani, Carlo

2016-08-01

This paper investigates the Korteweg-de Vries equation within the scope of the local fractional derivative formulation. The exact traveling wave solutions of non-differentiable type with the generalized functions defined on Cantor sets are analyzed. The results for the non-differentiable solutions when fractal dimension is 1 are also discussed. It is shown that the exact solutions for the local fractional Korteweg-de Vries equation characterize the fractal wave on shallow water surfaces.

A Predictor-Corrector Approach for the Numerical Solution of Fractional Differential Equations

NASA Technical Reports Server (NTRS)

Diethelm, Kai; Ford, Neville J.; Freed, Alan D.; Gray, Hugh R. (Technical Monitor)

2002-01-01

We discuss an Adams-type predictor-corrector method for the numerical solution of fractional differential equations. The method may be used both for linear and for nonlinear problems, and it may be extended to multi-term equations (involving more than one differential operator) too.

On the comparison of perturbation-iteration algorithm and residual power series method to solve fractional Zakharov-Kuznetsov equation

NASA Astrophysics Data System (ADS)

Şenol, Mehmet; Alquran, Marwan; Kasmaei, Hamed Daei

2018-06-01

In this paper, we present analytic-approximate solution of time-fractional Zakharov-Kuznetsov equation. This model demonstrates the behavior of weakly nonlinear ion acoustic waves in a plasma bearing cold ions and hot isothermal electrons in the presence of a uniform magnetic field. Basic definitions of fractional derivatives are described in the Caputo sense. Perturbation-iteration algorithm (PIA) and residual power series method (RPSM) are applied to solve this equation with success. The convergence analysis is also presented for both methods. Numerical results are given and then they are compared with the exact solutions. Comparison of the results reveal that both methods are competitive, powerful, reliable, simple to use and ready to apply to wide range of fractional partial differential equations.

A spectral tau algorithm based on Jacobi operational matrix for numerical solution of time fractional diffusion-wave equations

NASA Astrophysics Data System (ADS)

Bhrawy, A. H.; Doha, E. H.; Baleanu, D.; Ezz-Eldien, S. S.

2015-07-01

In this paper, an efficient and accurate spectral numerical method is presented for solving second-, fourth-order fractional diffusion-wave equations and fractional wave equations with damping. The proposed method is based on Jacobi tau spectral procedure together with the Jacobi operational matrix for fractional integrals, described in the Riemann-Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations in the unknown expansion coefficients of the sought-for spectral approximations. The validity and effectiveness of the method are demonstrated by solving five numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions more easier.

Solution of the Bagley Torvik equation by fractional DTM

NASA Astrophysics Data System (ADS)

Arora, Geeta; Pratiksha

2017-07-01

In this paper, fractional differential transform method(DTM) is implemented on the Bagley Torvik equation. This equation models the viscoelastic behavior of geological strata, metals, glasses etc. It explains the motion of a rigid plate immersed in a Newtonian fluid. DTM is a simple, reliable and efficient method that gives a series solution. Caputo fractional derivative is considered throughout this work. Two examples are given to demonstrate the validity and applicability of the method and comparison is made with the existing results.

Fractional vector calculus and fluid mechanics

NASA Astrophysics Data System (ADS)

Lazopoulos, Konstantinos A.; Lazopoulos, Anastasios K.

2017-04-01

Basic fluid mechanics equations are studied and revised under the prism of fractional continuum mechanics (FCM), a very promising research field that satisfies both experimental and theoretical demands. The geometry of the fractional differential has been clarified corrected and the geometry of the fractional tangent spaces of a manifold has been studied in Lazopoulos and Lazopoulos (Lazopoulos KA, Lazopoulos AK. Progr. Fract. Differ. Appl. 2016, 2, 85-104), providing the bases of the missing fractional differential geometry. Therefore, a lot can be contributed to fractional hydrodynamics: the basic fractional fluid equations (Navier Stokes, Euler and Bernoulli) are derived and fractional Darcy's flow in porous media is studied.

Blood flow problem in the presence of magnetic particles through a circular cylinder using Caputo-Fabrizio fractional derivative

NASA Astrophysics Data System (ADS)

Uddin, Salah; Mohamad, Mahathir; Khalid, Kamil; Abdulhammed, Mohammed; Saifullah Rusiman, Mohd; Che – Him, Norziha; Roslan, Rozaini

2018-04-01

In this paper, the flow of blood mixed with magnetic particles subjected to uniform transverse magnetic field and pressure gradient in an axisymmetric circular cylinder is studied by using a new trend of fractional derivative without singular kernel. The governing equations are fractional partial differential equations derived based on the Caputo-Fabrizio time-fractional derivatives NFDt. The current result agrees considerably well with that of the previous Caputo fractional derivatives UFDt.

Efficient numerical simulation of non-integer-order space-fractional reaction-diffusion equation via the Riemann-Liouville operator

NASA Astrophysics Data System (ADS)

Owolabi, Kolade M.

2018-03-01

In this work, we are concerned with the solution of non-integer space-fractional reaction-diffusion equations with the Riemann-Liouville space-fractional derivative in high dimensions. We approximate the Riemann-Liouville derivative with the Fourier transform method and advance the resulting system in time with any time-stepping solver. In the numerical experiments, we expect the travelling wave to arise from the given initial condition on the computational domain (-∞, ∞), which we terminate in the numerical experiments with a large but truncated value of L. It is necessary to choose L large enough to allow the waves to have enough space to distribute. Experimental results in high dimensions on the space-fractional reaction-diffusion models with applications to biological models (Fisher and Allen-Cahn equations) are considered. Simulation results reveal that fractional reaction-diffusion equations can give rise to a range of physical phenomena when compared to non-integer-order cases. As a result, most meaningful and practical situations are found to be modelled with the concept of fractional calculus.

Fractional Bateman—Feshbach Tikochinsky Oscillator

NASA Astrophysics Data System (ADS)

Dumitru, Baleanu; Jihad, H. Asad; Ivo, Petras

2014-02-01

In the last few years the numerical methods for solving the fractional differential equations started to be applied intensively to real world phenomena. Having these things in mind in this manuscript we focus on the fractional Lagrangian and Hamiltonian of the complex Bateman—Feshbach Tikochinsky oscillator. The numerical analysis of the corresponding fractional Euler-Lagrange equations is given within the Grünwald—Letnikov approach, which is power series expansion of the generating function.

Solution of time fractional Black-Scholes European option pricing equation arising in financial market

NASA Astrophysics Data System (ADS)

Ravi Kanth, A. S. V.; Aruna, K.

2016-12-01

In this paper, we present fractional differential transform method (FDTM) and modified fractional differential transform method (MFDTM) for the solution of time fractional Black-Scholes European option pricing equation. The method finds the solution without any discretization, transformation, or restrictive assumptions with the use of appropriate initial or boundary conditions. The efficiency and exactitude of the proposed methods are tested by means of three examples.

Maximum Principle for General Controlled Systems Driven by Fractional Brownian Motions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Han Yuecai; Hu Yaozhong; Song Jian, E-mail: jsong2@math.rutgers.edu

2013-04-15

We obtain a maximum principle for stochastic control problem of general controlled stochastic differential systems driven by fractional Brownian motions (of Hurst parameter H>1/2). This maximum principle specifies a system of equations that the optimal control must satisfy (necessary condition for the optimal control). This system of equations consists of a backward stochastic differential equation driven by both fractional Brownian motions and the corresponding underlying standard Brownian motions. In addition to this backward equation, the maximum principle also involves the Malliavin derivatives. Our approach is to use conditioning and Malliavin calculus. To arrive at our maximum principle we need tomore » develop some new results of stochastic analysis of the controlled systems driven by fractional Brownian motions via fractional calculus. Our approach of conditioning and Malliavin calculus is also applied to classical system driven by standard Brownian motions while the controller has only partial information. As a straightforward consequence, the classical maximum principle is also deduced in this more natural and simpler way.« less

A new analysis of the Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler-type kernel

NASA Astrophysics Data System (ADS)

Kumar, Devendra; Singh, Jagdev; Baleanu, Dumitru

2018-02-01

The mathematical model of breaking of non-linear dispersive water waves with memory effect is very important in mathematical physics. In the present article, we examine a novel fractional extension of the non-linear Fornberg-Whitham equation occurring in wave breaking. We consider the most recent theory of differentiation involving the non-singular kernel based on the extended Mittag-Leffler-type function to modify the Fornberg-Whitham equation. We examine the existence of the solution of the non-linear Fornberg-Whitham equation of fractional order. Further, we show the uniqueness of the solution. We obtain the numerical solution of the new arbitrary order model of the non-linear Fornberg-Whitham equation with the aid of the Laplace decomposition technique. The numerical outcomes are displayed in the form of graphs and tables. The results indicate that the Laplace decomposition algorithm is a very user-friendly and reliable scheme for handling such type of non-linear problems of fractional order.

A set for relational reasoning: Facilitation of algebraic modeling by a fraction task.

PubMed

DeWolf, Melissa; Bassok, Miriam; Holyoak, Keith J

2016-12-01

Recent work has identified correlations between early mastery of fractions and later math achievement, especially in algebra. However, causal connections between aspects of reasoning with fractions and improved algebra performance have yet to be established. The current study investigated whether relational reasoning with fractions facilitates subsequent algebraic reasoning using both pre-algebra students and adult college students. Participants were first given either a relational reasoning fractions task or a fraction algebra procedures control task. Then, all participants solved word problems and constructed algebraic equations in either multiplication or division format. The word problems and the equation construction tasks involved simple multiplicative comparison statements such as "There are 4 times as many students as teachers in a classroom." Performance on the algebraic equation construction task was enhanced for participants who had previously completed the relational fractions task compared with those who completed the fraction algebra procedures task. This finding suggests that relational reasoning with fractions can establish a relational set that promotes students' tendency to model relations using algebraic expressions. Copyright © 2016 Elsevier Inc. All rights reserved.

Anomalous dielectric relaxation with linear reaction dynamics in space-dependent force fields.

PubMed

Hong, Tao; Tang, Zhengming; Zhu, Huacheng

2016-12-28

The anomalous dielectric relaxation of disordered reaction with linear reaction dynamics is studied via the continuous time random walk model in the presence of space-dependent electric field. Two kinds of modified reaction-subdiffusion equations are derived for different linear reaction processes by the master equation, including the instantaneous annihilation reaction and the noninstantaneous annihilation reaction. If a constant proportion of walkers is added or removed instantaneously at the end of each step, there will be a modified reaction-subdiffusion equation with a fractional order temporal derivative operating on both the standard diffusion term and a linear reaction kinetics term. If the walkers are added or removed at a constant per capita rate during the waiting time between steps, there will be a standard linear reaction kinetics term but a fractional order temporal derivative operating on an anomalous diffusion term. The dielectric polarization is analyzed based on the Legendre polynomials and the dielectric properties of both reactions can be expressed by the effective rotational diffusion function and component concentration function, which is similar to the standard reaction-diffusion process. The results show that the effective permittivity can be used to describe the dielectric properties in these reactions if the chemical reaction time is much longer than the relaxation time.

Closed form solutions of two time fractional nonlinear wave equations

NASA Astrophysics Data System (ADS)

Akbar, M. Ali; Ali, Norhashidah Hj. Mohd.; Roy, Ripan

2018-06-01

In this article, we investigate the exact traveling wave solutions of two nonlinear time fractional wave equations. The fractional derivatives are described in the sense of conformable fractional derivatives. In addition, the traveling wave solutions are accomplished in the form of hyperbolic, trigonometric, and rational functions involving free parameters. To investigate such types of solutions, we implement the new generalized (G‧ / G) -expansion method. The extracted solutions are reliable, useful and suitable to comprehend the optimal control problems, chaotic vibrations, global and local bifurcations and resonances, furthermore, fission and fusion phenomena occur in solitons, the relativistic energy-momentum relation, scalar electrodynamics, quantum relativistic one-particle theory, electromagnetic interactions etc. The results reveal that the method is very fruitful and convenient for exploring nonlinear differential equations of fractional order treated in theoretical physics.

Parareal algorithms with local time-integrators for time fractional differential equations

NASA Astrophysics Data System (ADS)

Wu, Shu-Lin; Zhou, Tao

2018-04-01

It is challenge work to design parareal algorithms for time-fractional differential equations due to the historical effect of the fractional operator. A direct extension of the classical parareal method to such equations will lead to unbalance computational time in each process. In this work, we present an efficient parareal iteration scheme to overcome this issue, by adopting two recently developed local time-integrators for time fractional operators. In both approaches, one introduces auxiliary variables to localized the fractional operator. To this end, we propose a new strategy to perform the coarse grid correction so that the auxiliary variables and the solution variable are corrected separately in a mixed pattern. It is shown that the proposed parareal algorithm admits robust rate of convergence. Numerical examples are presented to support our conclusions.

Existence of a coupled system of fractional differential equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ibrahim, Rabha W.; Siri, Zailan

2015-10-22

We manage the existence and uniqueness of a fractional coupled system containing Schrödinger equations. Such a system appears in quantum mechanics. We confirm that the fractional system under consideration admits a global solution in appropriate functional spaces. The solution is shown to be unique. The method is based on analytic technique of the fixed point theory. The fractional differential operator is considered from the virtue of the Riemann-Liouville differential operator.

Approximate Solution of Time-Fractional Advection-Dispersion Equation via Fractional Variational Iteration Method

PubMed Central

İbiş, Birol

2014-01-01

This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662

Kappa Distribution in a Homogeneous Medium: Adiabatic Limit of a Super-diffusive Process?

NASA Astrophysics Data System (ADS)

Roth, I.

2015-12-01

The classical statistical theory predicts that an ergodic, weakly interacting system like charged particles in the presence of electromagnetic fields, performing Brownian motions (characterized by small range deviations in phase space and short-term microscopic memory), converges into the Gibbs-Boltzmann statistics. Observation of distributions with a kappa-power-law tails in homogeneous systems contradicts this prediction and necessitates a renewed analysis of the basic axioms of the diffusion process: characteristics of the transition probability density function (pdf) for a single interaction, with a possibility of non-Markovian process and non-local interaction. The non-local, Levy walk deviation is related to the non-extensive statistical framework. Particles bouncing along (solar) magnetic field with evolving pitch angles, phases and velocities, as they interact resonantly with waves, undergo energy changes at undetermined time intervals, satisfying these postulates. The dynamic evolution of a general continuous time random walk is determined by pdf of jumps and waiting times resulting in a fractional Fokker-Planck equation with non-integer derivatives whose solution is given by a Fox H-function. The resulting procedure involves the known, although not frequently used in physics fractional calculus, while the local, Markovian process recasts the evolution into the standard Fokker-Planck equation. Solution of the fractional Fokker-Planck equation with the help of Mellin transform and evaluation of its residues at the poles of its Gamma functions results in a slowly converging sum with power laws. It is suggested that these tails form the Kappa function. Gradual vs impulsive solar electron distributions serve as prototypes of this description.

On the joint bimodality of temperature and moisture near stratocumulus cloud tops

NASA Technical Reports Server (NTRS)

Randall, D. A.

1983-01-01

The observed distributions of the thermodynamic variables near stratocumulus top are highly bimodal. Two simple models of sub-grid fractional cloudiness motivated by this observed bimodality are examined. In both models, certain low order moments of two independent, moist-conservative thermodynamic variables are assumed to be known. The first model is based on the assumption of two discrete populations of parcels: a warm-day population and a cool-moist population. If only the first and second moments are assumed to be known, the number of unknowns exceeds the number of independent equations. If the third moments are assumed to be known as well, the number of independent equations exceeds the number of unknowns. The second model is based on the assumption of a continuous joint bimodal distribution of parcels, obtained as the weighted sum of two binormal distributions. For this model, the third moments are used to obtain 9 independent nonlinear algebraic equations in 11 unknowns. Two additional equations are needed to determine the covariance within the two subpopulations. In case these two internal covariance vanish, the system of equations can be solved analytically.

Lax Integrability and the Peakon Problem for the Modified Camassa-Holm Equation

NASA Astrophysics Data System (ADS)

Chang, Xiangke; Szmigielski, Jacek

2018-02-01

Peakons are special weak solutions of a class of nonlinear partial differential equations modelling non-linear phenomena such as the breakdown of regularity and the onset of shocks. We show that the natural concept of weak solutions in the case of the modified Camassa-Holm equation studied in this paper is dictated by the distributional compatibility of its Lax pair and, as a result, it differs from the one proposed and used in the literature based on the concept of weak solutions used for equations of the Burgers type. Subsequently, we give a complete construction of peakon solutions satisfying the modified Camassa-Holm equation in the sense of distributions; our approach is based on solving certain inverse boundary value problem, the solution of which hinges on a combination of classical techniques of analysis involving Stieltjes' continued fractions and multi-point Padé approximations. We propose sufficient conditions needed to ensure the global existence of peakon solutions and analyze the large time asymptotic behaviour whose special features include a formation of pairs of peakons that share asymptotic speeds, as well as Toda-like sorting property.

Multiple positive solutions to a coupled systems of nonlinear fractional differential equations.

PubMed

Shah, Kamal; Khan, Rahmat Ali

2016-01-01

In this article, we study existence, uniqueness and nonexistence of positive solution to a highly nonlinear coupled system of fractional order differential equations. Necessary and sufficient conditions for the existence and uniqueness of positive solution are developed by using Perov's fixed point theorem for the considered problem. Further, we also established sufficient conditions for existence of multiplicity results for positive solutions. Also, we developed some conditions under which the considered coupled system of fractional order differential equations has no positive solution. Appropriate examples are also provided which demonstrate our results.

Maximal regularity in lp spaces for discrete time fractional shifted equations

NASA Astrophysics Data System (ADS)

Lizama, Carlos; Murillo-Arcila, Marina

2017-09-01

In this paper, we are presenting a new method based on operator-valued Fourier multipliers to characterize the existence and uniqueness of ℓp-solutions for discrete time fractional models in the form where A is a closed linear operator defined on a Banach space X and Δα denotes the Grünwald-Letnikov fractional derivative of order α > 0. If X is a UMD space, we provide this characterization only in terms of the R-boundedness of the operator-valued symbol associated to the abstract model. To illustrate our results, we derive new qualitative properties of nonlinear difference equations with shiftings, including fractional versions of the logistic and Nagumo equations.

Reply to "Comment on 'Fractional quantum mechanics' and 'Fractional Schrödinger equation' ".

PubMed

Laskin, Nick

2016-06-01

The fractional uncertainty relation is a mathematical formulation of Heisenberg's uncertainty principle in the framework of fractional quantum mechanics. Two mistaken statements presented in the Comment have been revealed. The origin of each mistaken statement has been clarified and corrected statements have been made. A map between standard quantum mechanics and fractional quantum mechanics has been presented to emphasize the features of fractional quantum mechanics and to avoid misinterpretations of the fractional uncertainty relation. It has been shown that the fractional probability current equation is correct in the area of its applicability. Further studies have to be done to find meaningful quantum physics problems with involvement of the fractional probability current density vector and the extra term emerging in the framework of fractional quantum mechanics.

A pseudo energy-invariant method for relativistic wave equations with Riesz space-fractional derivatives

NASA Astrophysics Data System (ADS)

Macías-Díaz, J. E.; Hendy, A. S.; De Staelen, R. H.

2018-03-01

In this work, we investigate a general nonlinear wave equation with Riesz space-fractional derivatives that generalizes various classical hyperbolic models, including the sine-Gordon and the Klein-Gordon equations from relativistic quantum mechanics. A finite-difference discretization of the model is provided using fractional centered differences. The method is a technique that is capable of preserving an energy-like quantity at each iteration. Some computational comparisons against solutions available in the literature are performed in order to assess the capability of the method to preserve the invariant. Our experiments confirm that the technique yields good approximations to the solutions considered. As an application of our scheme, we provide simulations that confirm, for the first time in the literature, the presence of the phenomenon of nonlinear supratransmission in Riesz space-fractional Klein-Gordon equations driven by a harmonic perturbation at the boundary.

A multi-domain spectral method for time-fractional differential equations

NASA Astrophysics Data System (ADS)

Chen, Feng; Xu, Qinwu; Hesthaven, Jan S.

2015-07-01

This paper proposes an approach for high-order time integration within a multi-domain setting for time-fractional differential equations. Since the kernel is singular or nearly singular, two main difficulties arise after the domain decomposition: how to properly account for the history/memory part and how to perform the integration accurately. To address these issues, we propose a novel hybrid approach for the numerical integration based on the combination of three-term-recurrence relations of Jacobi polynomials and high-order Gauss quadrature. The different approximations used in the hybrid approach are justified theoretically and through numerical examples. Based on this, we propose a new multi-domain spectral method for high-order accurate time integrations and study its stability properties by identifying the method as a generalized linear method. Numerical experiments confirm hp-convergence for both time-fractional differential equations and time-fractional partial differential equations.

On absolutely continuous weakly mixing cocycles over irrational rotations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rozhdestvenskii, A V

2003-06-30

A weakly mixing cocycle over a rotation {alpha} is a measurable function {phi}:S{sup 1}{yields}S{sup 1}, where S{sup 1}={l_brace}z element of C:|z|=1{r_brace}, such that the equation {phi}{sup n}(z)=c (h(exp(2{pi}i{alpha})z))/(h(z))for almost all z; (*) has no measurable solutions h( {center_dot} ):S{sup 1}{yields}S{sup 1} for any n element of Z{l_brace}0{r_brace} and c element of C, |c|=1. If the irrational number {alpha} has bounded convergents in its continued fraction expansion and a function M(y) increases more slowly than y ln{sup 1/2}y, then it is proved that there exists a weakly mixing cocycle of the form {phi}(exp(2{pi}ix))=exp(2{pi}i{phi}-tilde(x)), where {phi}-tilde:T{yields}R belongs to the class W{sup 1}(M(L)(T)).more » In addition, it is shown that equation (*) (and also the corresponding additive cohomological equation) is soluble for {phi}-tilde element of W{sup 1}(L log{sub +}{sup 1/2}L(T))« less

Time-domain comparisons of power law attenuation in causal and noncausal time-fractional wave equations

PubMed Central

Zhao, Xiaofeng; McGough, Robert J.

2016-01-01

The attenuation of ultrasound propagating in human tissue follows a power law with respect to frequency that is modeled by several different causal and noncausal fractional partial differential equations. To demonstrate some of the similarities and differences that are observed in three related time-fractional partial differential equations, time-domain Green's functions are calculated numerically for the power law wave equation, the Szabo wave equation, and for the Caputo wave equation. These Green's functions are evaluated for water with a power law exponent of y = 2, breast with a power law exponent of y = 1.5, and liver with a power law exponent of y = 1.139. Simulation results show that the noncausal features of the numerically calculated time-domain response are only evident very close to the source and that these causal and noncausal time-domain Green's functions converge to the same result away from the source. When noncausal time-domain Green's functions are convolved with a short pulse, no evidence of noncausal behavior remains in the time-domain, which suggests that these causal and noncausal time-fractional models are equally effective for these numerical calculations. PMID:27250193

Comment on “Maxwell's equations and electromagnetic Lagrangian density in fractional form” [J. Math. Phys. 53, 033505 (2012)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rabei, Eqab M.; Al-Jamel, A.; Widyan, H.

In a recent paper, Jaradat et al. [J. Math. Phys. 53, 033505 (2012)] have presented the fractional form of the electromagnetic Lagrangian density within the Riemann-Liouville fractional derivative. They claimed that the Agrawal procedure [O. P. Agrawal, J. Math. Anal. Appl. 272, 368 (2002)] is used to obtain Maxwell's equations in the fractional form, and the Hamilton's equations of motion together with the conserved quantities obtained from fractional Noether's theorem are reported. In this comment, we draw the attention that there are some serious steps of the procedure used in their work are not applicable even though their final resultsmore » are correct. Their work should have been done based on a formulation as reported by Baleanu and Muslih [Phys. Scr. 72, 119 (2005)].« less

General solution of the Bagley-Torvik equation with fractional-order derivative

NASA Astrophysics Data System (ADS)

Wang, Z. H.; Wang, X.

2010-05-01

This paper investigates the general solution of the Bagley-Torvik equation with 1/2-order derivative or 3/2-order derivative. This fractional-order differential equation is changed into a sequential fractional-order differential equation (SFDE) with constant coefficients. Then the general solution of the SFDE is expressed as the linear combination of fundamental solutions that are in terms of α-exponential functions, a kind of functions that play the same role of the classical exponential function. Because the number of fundamental solutions of the SFDE is greater than 2, the general solution of the SFDE depends on more than two free (independent) constants. This paper shows that the general solution of the Bagley-Torvik equation involves actually two free constants only, and it can be determined fully by the initial displacement and initial velocity.

A higher order numerical method for time fractional partial differential equations with nonsmooth data

NASA Astrophysics Data System (ADS)

Xing, Yanyuan; Yan, Yubin

2018-03-01

Gao et al. [11] (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate O (k 3 - α), 0 < α < 1 by directly approximating the integer-order derivative with some finite difference quotients in the definition of the Caputo fractional derivative, see also Lv and Xu [20] (2016), where k is the time step size. Under the assumption that the solution of the time fractional partial differential equation is sufficiently smooth, Lv and Xu [20] (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate O (k 3 - α), 0 < α < 1 uniformly with respect to the time variable t. However, in general the solution of the time fractional partial differential equation has low regularity and in this case the numerical method fails to have the convergence rate O (k 3 - α), 0 < α < 1 uniformly with respect to the time variable t. In this paper, we first obtain a similar approximation scheme to the Riemann-Liouville fractional derivative with the convergence rate O (k 3 - α), 0 < α < 1 as in Gao et al. [11] (2014) by approximating the Hadamard finite-part integral with the piecewise quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate O (k 3 - α), 0 < α < 1 for any fixed tn > 0 for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

Numerical simulation for solution of space-time fractional telegraphs equations with local fractional derivatives via HAFSTM

NASA Astrophysics Data System (ADS)

Pandey, Rishi Kumar; Mishra, Hradyesh Kumar

2017-11-01

In this paper, the semi-analytic numerical technique for the solution of time-space fractional telegraph equation is applied. This numerical technique is based on coupling of the homotopy analysis method and sumudu transform. It shows the clear advantage with mess methods like finite difference method and also with polynomial methods similar to perturbation and Adomian decomposition methods. It is easily transform the complex fractional order derivatives in simple time domain and interpret the results in same meaning.

Homotopy decomposition method for solving one-dimensional time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Abuasad, Salah; Hashim, Ishak

2018-04-01

In this paper, we present the homotopy decomposition method with a modified definition of beta fractional derivative for the first time to find exact solution of one-dimensional time-fractional diffusion equation. In this method, the solution takes the form of a convergent series with easily computable terms. The exact solution obtained by the proposed method is compared with the exact solution obtained by using fractional variational homotopy perturbation iteration method via a modified Riemann-Liouville derivative.

The meshless local Petrov-Galerkin method based on moving Kriging interpolation for solving the time fractional Navier-Stokes equations.

PubMed

Thamareerat, N; Luadsong, A; Aschariyaphotha, N

2016-01-01

In this paper, we present a numerical scheme used to solve the nonlinear time fractional Navier-Stokes equations in two dimensions. We first employ the meshless local Petrov-Galerkin (MLPG) method based on a local weak formulation to form the system of discretized equations and then we will approximate the time fractional derivative interpreted in the sense of Caputo by a simple quadrature formula. The moving Kriging interpolation which possesses the Kronecker delta property is applied to construct shape functions. This research aims to extend and develop further the applicability of the truly MLPG method to the generalized incompressible Navier-Stokes equations. Two numerical examples are provided to illustrate the accuracy and efficiency of the proposed algorithm. Very good agreement between the numerically and analytically computed solutions can be observed in the verification. The present MLPG method has proved its efficiency and reliability for solving the two-dimensional time fractional Navier-Stokes equations arising in fluid dynamics as well as several other problems in science and engineering.

Hipergeometric solutions to some nonhomogeneous equations of fractional order

NASA Astrophysics Data System (ADS)

Olivares, Jorge; Martin, Pablo; Maass, Fernando

2017-12-01

In this paper a study is performed to the solution of the linear non homogeneous fractional order alpha differential equation equal to I 0(x), where I 0(x) is the modified Bessel function of order zero, the initial condition is f(0)=0 and 0 < alpha < 1. Caputo definition for the fractional derivatives is considered. Fractional derivatives have become important in physical and chemical phenomena as visco-elasticity and visco-plasticity, anomalous diffusion and electric circuits. In particular in this work the values of alpha=1/2, 1/4 and 3/4. are explicitly considered . In these cases Laplace transform is applied, and later the inverse Laplace transform leads to the solutions of the differential equation, which become hypergeometric functions.

A note on the evolution equations from the area fraction and the thickness of a floating ice cover

NASA Astrophysics Data System (ADS)

Schulkes, R. M. S. M.

1995-03-01

In this paper, two sets of evolution equations for the area fraction and the ice thickness are investigated. First of all, a simplified alternative derivation of the evolution equations as presented by Gray and Morland (1994) is given. In addition, it is shown that with proper identification of ridging functions, there is a close connection between the derived equations and the thickness distribution model introduced by Thorndike et al. (1975).

An efficient technique for higher order fractional differential equation.

PubMed

Ali, Ayyaz; Iqbal, Muhammad Asad; Ul-Hassan, Qazi Mahmood; Ahmad, Jamshad; Mohyud-Din, Syed Tauseef

2016-01-01

In this study, we establish exact solutions of fractional Kawahara equation by using the idea of [Formula: see text]-expansion method. The results of different studies show that the method is very effective and can be used as an alternative for finding exact solutions of nonlinear evolution equations (NLEEs) in mathematical physics. The solitary wave solutions are expressed by the hyperbolic, trigonometric, exponential and rational functions. Graphical representations along with the numerical data reinforce the efficacy of the used procedure. The specified idea is very effective, expedient for fractional PDEs, and could be extended to other physical problems.

The Use of Generalized Laguerre Polynomials in Spectral Methods for Solving Fractional Delay Differential Equations.

PubMed

Khader, M M

2013-10-01

In this paper, an efficient numerical method for solving the fractional delay differential equations (FDDEs) is considered. The fractional derivative is described in the Caputo sense. The proposed method is based on the derived approximate formula of the Laguerre polynomials. The properties of Laguerre polynomials are utilized to reduce FDDEs to a linear or nonlinear system of algebraic equations. Special attention is given to study the error and the convergence analysis of the proposed method. Several numerical examples are provided to confirm that the proposed method is in excellent agreement with the exact solution.

Solving Fuzzy Fractional Differential Equations Using Zadeh's Extension Principle

PubMed Central

Ahmad, M. Z.; Hasan, M. K.; Abbasbandy, S.

2013-01-01

We study a fuzzy fractional differential equation (FFDE) and present its solution using Zadeh's extension principle. The proposed study extends the case of fuzzy differential equations of integer order. We also propose a numerical method to approximate the solution of FFDEs. To solve nonlinear problems, the proposed numerical method is then incorporated into an unconstrained optimisation technique. Several numerical examples are provided. PMID:24082853

Tempered fractional calculus

NASA Astrophysics Data System (ADS)

Sabzikar, Farzad; Meerschaert, Mark M.; Chen, Jinghua

2015-07-01

Fractional derivatives and integrals are convolutions with a power law. Multiplying by an exponential factor leads to tempered fractional derivatives and integrals. Tempered fractional diffusion equations, where the usual second derivative in space is replaced by a tempered fractional derivative, govern the limits of random walk models with an exponentially tempered power law jump distribution. The limiting tempered stable probability densities exhibit semi-heavy tails, which are commonly observed in finance. Tempered power law waiting times lead to tempered fractional time derivatives, which have proven useful in geophysics. The tempered fractional derivative or integral of a Brownian motion, called a tempered fractional Brownian motion, can exhibit semi-long range dependence. The increments of this process, called tempered fractional Gaussian noise, provide a useful new stochastic model for wind speed data. A tempered fractional difference forms the basis for numerical methods to solve tempered fractional diffusion equations, and it also provides a useful new correlation model in time series.

Tempered fractional calculus

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sabzikar, Farzad, E-mail: sabzika2@stt.msu.edu; Meerschaert, Mark M., E-mail: mcubed@stt.msu.edu; Chen, Jinghua, E-mail: cjhdzdz@163.com

2015-07-15

Fractional derivatives and integrals are convolutions with a power law. Multiplying by an exponential factor leads to tempered fractional derivatives and integrals. Tempered fractional diffusion equations, where the usual second derivative in space is replaced by a tempered fractional derivative, govern the limits of random walk models with an exponentially tempered power law jump distribution. The limiting tempered stable probability densities exhibit semi-heavy tails, which are commonly observed in finance. Tempered power law waiting times lead to tempered fractional time derivatives, which have proven useful in geophysics. The tempered fractional derivative or integral of a Brownian motion, called a temperedmore » fractional Brownian motion, can exhibit semi-long range dependence. The increments of this process, called tempered fractional Gaussian noise, provide a useful new stochastic model for wind speed data. A tempered fractional difference forms the basis for numerical methods to solve tempered fractional diffusion equations, and it also provides a useful new correlation model in time series.« less

Constructing and predicting solitary pattern solutions for nonlinear time-fractional dispersive partial differential equations

NASA Astrophysics Data System (ADS)

Arqub, Omar Abu; El-Ajou, Ahmad; Momani, Shaher

2015-07-01

Building fractional mathematical models for specific phenomena and developing numerical or analytical solutions for these fractional mathematical models are crucial issues in mathematics, physics, and engineering. In this work, a new analytical technique for constructing and predicting solitary pattern solutions of time-fractional dispersive partial differential equations is proposed based on the generalized Taylor series formula and residual error function. The new approach provides solutions in the form of a rapidly convergent series with easily computable components using symbolic computation software. For method evaluation and validation, the proposed technique was applied to three different models and compared with some of the well-known methods. The resultant simulations clearly demonstrate the superiority and potentiality of the proposed technique in terms of the quality performance and accuracy of substructure preservation in the construct, as well as the prediction of solitary pattern solutions for time-fractional dispersive partial differential equations.

Fractional-order in a macroeconomic dynamic model

NASA Astrophysics Data System (ADS)

David, S. A.; Quintino, D. D.; Soliani, J.

2013-10-01

In this paper, we applied the Riemann-Liouville approach in order to realize the numerical simulations to a set of equations that represent a fractional-order macroeconomic dynamic model. It is a generalization of a dynamic model recently reported in the literature. The aforementioned equations have been simulated for several cases involving integer and non-integer order analysis, with some different values to fractional order. The time histories and the phase diagrams have been plotted to visualize the effect of fractional order approach. The new contribution of this work arises from the fact that the macroeconomic dynamic model proposed here involves the public sector deficit equation, which renders the model more realistic and complete when compared with the ones encountered in the literature. The results reveal that the fractional-order macroeconomic model can exhibit a real reasonable behavior to macroeconomics systems and might offer greater insights towards the understanding of these complex dynamic systems.

Competition between plant functional types in the Canadian Terrestrial Ecosystem Model (CTEM) v. 2.0

NASA Astrophysics Data System (ADS)

Melton, J. R.; Arora, V. K.

2015-06-01

The Canadian Terrestrial Ecosystem Model (CTEM) is the interactive vegetation component in the Earth system model of the Canadian Centre for Climate Modelling and Analysis. CTEM models land-atmosphere exchange of CO2 through the response of carbon in living vegetation, and dead litter and soil pools, to changes in weather and climate at timescales of days to centuries. Version 1.0 of CTEM uses prescribed fractional coverage of plant functional types (PFTs) although, in reality, vegetation cover continually adapts to changes in climate, atmospheric composition, and anthropogenic forcing. Changes in the spatial distribution of vegetation occur on timescales of years to centuries as vegetation distributions inherently have inertia. Here, we present version 2.0 of CTEM which includes a representation of competition between PFTs based on a modified version of the Lotka-Volterra (L-V) predator-prey equations. Our approach is used to dynamically simulate the fractional coverage of CTEM's seven natural, non-crop PFTs which are then compared with available observation-based estimates. Results from CTEM v. 2.0 show the model is able to represent the broad spatial distributions of its seven PFTs at the global scale. However, differences remain between modelled and observation-based fractional coverages of PFTs since representing the multitude of plant species globally, with just seven non-crop PFTs, only captures the large scale climatic controls on PFT distributions. As expected, PFTs that exist in climate niches are difficult to represent either due to the coarse spatial resolution of the model, and the corresponding driving climate, or the limited number of PFTs used. We also simulate the fractional coverages of PFTs using unmodified L-V equations to illustrate its limitations. The geographic and zonal distributions of primary terrestrial carbon pools and fluxes from the versions of CTEM that use prescribed and dynamically simulated fractional coverage of PFTs compare reasonably well with each other and observation-based estimates. The parametrization of competition between PFTs in CTEM v. 2.0 based on the modified L-V equations behaves in a reasonably realistic manner and yields a tool with which to investigate the changes in spatial distribution of vegetation in response to future changes in climate.

Competition between plant functional types in the Canadian Terrestrial Ecosystem Model (CTEM) v. 2.0

NASA Astrophysics Data System (ADS)

Melton, J. R.; Arora, V. K.

2016-01-01

The Canadian Terrestrial Ecosystem Model (CTEM) is the interactive vegetation component in the Earth system model of the Canadian Centre for Climate Modelling and Analysis. CTEM models land-atmosphere exchange of CO2 through the response of carbon in living vegetation, and dead litter and soil pools, to changes in weather and climate at timescales of days to centuries. Version 1.0 of CTEM uses prescribed fractional coverage of plant functional types (PFTs) although, in reality, vegetation cover continually adapts to changes in climate, atmospheric composition and anthropogenic forcing. Changes in the spatial distribution of vegetation occur on timescales of years to centuries as vegetation distributions inherently have inertia. Here, we present version 2.0 of CTEM, which includes a representation of competition between PFTs based on a modified version of the Lotka-Volterra (L-V) predator-prey equations. Our approach is used to dynamically simulate the fractional coverage of CTEM's seven natural, non-crop PFTs, which are then compared with available observation-based estimates. Results from CTEM v. 2.0 show the model is able to represent the broad spatial distributions of its seven PFTs at the global scale. However, differences remain between modelled and observation-based fractional coverage of PFTs since representing the multitude of plant species globally, with just seven non-crop PFTs, only captures the large-scale climatic controls on PFT distributions. As expected, PFTs that exist in climate niches are difficult to represent either due to the coarse spatial resolution of the model, and the corresponding driving climate, or the limited number of PFTs used. We also simulate the fractional coverage of PFTs using unmodified L-V equations to illustrate its limitations. The geographic and zonal distributions of primary terrestrial carbon pools and fluxes from the versions of CTEM that use prescribed and dynamically simulated fractional coverage of PFTs compare reasonably well with each other and observation-based estimates. The parametrization of competition between PFTs in CTEM v. 2.0 based on the modified L-V equations behaves in a reasonably realistic manner and yields a tool with which to investigate the changes in spatial distribution of vegetation in response to future changes in climate.

Kinetic theory of oxygen isotopic exchange between minerals and water

USGS Publications Warehouse

Criss, R.E.; Gregory, R.T.; Taylor, H.P.

1987-01-01

Kinetic and mass conservation equations are used to describe oxygen isotopic exchange between minerals and water in "closed" and open hydrothermal systems. In cases where n coexisting mineral phases having different reaction rates are present, the exchange process is described by a system of n + 1 simultaneous differential equations consisting of n pseudo first-order rate equations and a conservation of mass equation. The simultaneous solutions to these equations generate curved exchange trajectories on ??-?? plots. Families of such trajectories generated under conditions allowing for different fluid mole fractions, different fluid isotopic compositions, or different fluid flow rates are connected by positive-sloped isochronous lines. These isochrons reproduce the effects observed in hydrothermally exchanged mineral pairs including 1) steep positive slopes, 2) common reversals in the measured fractionation factors (??), and 3) measured fractionations that are highly variable over short distances where no thermal gradient can be geologically demonstrated. ?? 1987.

Black Hole and Galaxy Coevolution from Continuity Equation and Abundance Matching

NASA Astrophysics Data System (ADS)

Aversa, R.; Lapi, A.; de Zotti, G.; Shankar, F.; Danese, L.

2015-09-01

We investigate the coevolution of galaxies and hosted supermassive black holes (BHs) throughout the history of the universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source terms to reconstruct the supermassive BH mass function from the active galactic nucleus (AGN) luminosity functions. Such an approach includes physically motivated AGN light curves tested on independent data sets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disk conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the ultraviolet and far-IR luminosity functions associated with star formation in galaxies. These results imply that the build-up of stars and BHs in galaxies occurs via in situ processes, with dry mergers playing a marginal role at least for stellar masses ≲ 3× {10}11 {M}⊙ and BH masses ≲ {10}9 {M}⊙ , where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter (DM) component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. In addition, they may be operationally implemented in numerical simulations to populate DM halos or to gauge subgrid physics. Moreover, they may be exploited to investigate the galaxy/AGN clustering as a function of redshift, mass, and/or luminosity. In fact, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, thus further validating our results from the continuity equation. Finally, our analysis highlights that (i) the fraction of AGNs observed in the slim-disk regime, where most of the BH mass is accreted, increases with redshift; and (ii) already at z≳ 6 a substantial amount of dust must have formed over timescales ≲ {10}8 yr in strongly star-forming galaxies, making these sources well within the reach of ALMA surveys in (sub)millimeter bands.

TEMPERED FRACTIONAL CALCULUS.

PubMed

Meerschaert, Mark M; Sabzikar, Farzad; Chen, Jinghua

2015-07-15

Fractional derivatives and integrals are convolutions with a power law. Multiplying by an exponential factor leads to tempered fractional derivatives and integrals. Tempered fractional diffusion equations, where the usual second derivative in space is replaced by a tempered fractional derivative, govern the limits of random walk models with an exponentially tempered power law jump distribution. The limiting tempered stable probability densities exhibit semi-heavy tails, which are commonly observed in finance. Tempered power law waiting times lead to tempered fractional time derivatives, which have proven useful in geophysics. The tempered fractional derivative or integral of a Brownian motion, called a tempered fractional Brownian motion, can exhibit semi-long range dependence. The increments of this process, called tempered fractional Gaussian noise, provide a useful new stochastic model for wind speed data. A tempered difference forms the basis for numerical methods to solve tempered fractional diffusion equations, and it also provides a useful new correlation model in time series.

TEMPERED FRACTIONAL CALCULUS

PubMed Central

MEERSCHAERT, MARK M.; SABZIKAR, FARZAD; CHEN, JINGHUA

2014-01-01

Fractional derivatives and integrals are convolutions with a power law. Multiplying by an exponential factor leads to tempered fractional derivatives and integrals. Tempered fractional diffusion equations, where the usual second derivative in space is replaced by a tempered fractional derivative, govern the limits of random walk models with an exponentially tempered power law jump distribution. The limiting tempered stable probability densities exhibit semi-heavy tails, which are commonly observed in finance. Tempered power law waiting times lead to tempered fractional time derivatives, which have proven useful in geophysics. The tempered fractional derivative or integral of a Brownian motion, called a tempered fractional Brownian motion, can exhibit semi-long range dependence. The increments of this process, called tempered fractional Gaussian noise, provide a useful new stochastic model for wind speed data. A tempered difference forms the basis for numerical methods to solve tempered fractional diffusion equations, and it also provides a useful new correlation model in time series. PMID:26085690

Generalized fractional diffusion equations for accelerating subdiffusion and truncated Lévy flights

NASA Astrophysics Data System (ADS)

Chechkin, A. V.; Gonchar, V. Yu.; Gorenflo, R.; Korabel, N.; Sokolov, I. M.

2008-08-01

Fractional diffusion equations are widely used to describe anomalous diffusion processes where the characteristic displacement scales as a power of time. For processes lacking such scaling the corresponding description may be given by diffusion equations with fractional derivatives of distributed order. Such equations were introduced in A. V. Chechkin, R. Gorenflo, and I. Sokolov [Phys. Rev. E 66, 046129 (2002)] for the description of the processes getting more anomalous in the course of time (decelerating subdiffusion and accelerating superdiffusion). Here we discuss the properties of diffusion equations with fractional derivatives of the distributed order for the description of anomalous relaxation and diffusion phenomena getting less anomalous in the course of time, which we call, respectively, accelerating subdiffusion and decelerating superdiffusion. For the former process, by taking a relatively simple particular example with two fixed anomalous diffusion exponents we show that the proposed equation effectively describes the subdiffusion phenomenon with diffusion exponent varying in time. For the latter process we demonstrate by a particular example how the power-law truncated Lévy stable distribution evolves in time to the distribution with power-law asymptotics and Gaussian shape in the central part. The special case of two different orders is characteristic for the general situation in which the extreme orders dominate the asymptotics.

Direct computation of thermodynamic properties of chemically reacting air with consideration to CFD

NASA Astrophysics Data System (ADS)

Iannelli, Joe

2003-10-01

This paper details a two-equation procedure to calculate exactly mass and mole fractions, pressure, temperature, specific heats, speed of sound and the thermodynamic and jacobian partial derivatives of pressure and temperature for a five-species chemically reacting equilibrium air. The procedure generates these thermodynamic properties using as independent variables either pressure and temperature or density and internal energy, for CFD applications. An original element in this procedure consists in the exact physically meaningful solution of the mass-fraction and mass-action equations. Air-equivalent molecular masses for oxygen and nitrogen are then developed to account, within a mixture of only oxygen and nitrogen, for the presence of carbon dioxide, argon and the other noble gases within atmospheric air. The mathematical formulation also introduces a versatile system non-dimensionalization that makes the procedure uniformly applicable to flows ranging from shock-tube flows with zero initial velocity to aerothermodynamic flows with supersonic/hypersonic free-stream Mach numbers. Over a temperature range of more than 10000 K and pressure and density ranges corresponding to an increase in altitude in standard atmosphere of 30000 m above sea level, the predicted distributions of mole fractions, constant-volume specific heat, and speed of sound for the model five species agree with independently published results, and all the calculated thermodynamic properties, including their partial derivatives, remain continuous, smooth, and physically meaningful.

A new modification in the exponential rational function method for nonlinear fractional differential equations

NASA Astrophysics Data System (ADS)

Ahmed, Naveed; Bibi, Sadaf; Khan, Umar; Mohyud-Din, Syed Tauseef

2018-02-01

We have modified the traditional exponential rational function method (ERFM) and have used it to find the exact solutions of two different fractional partial differential equations, one is the time fractional Boussinesq equation and the other is the (2+1)-dimensional time fractional Zoomeron equation. In both the cases it is observed that the modified scheme provides more types of solutions than the traditional one. Moreover, a comparison of the recent solutions is made with some already existing solutions. We can confidently conclude that the modified scheme works better and provides more types of solutions with almost similar computational cost. Our generalized solutions include periodic, soliton-like, singular soliton and kink solutions. A graphical simulation of all types of solutions is provided and the correctness of the solution is verified by direct substitution. The extended version of the solutions is expected to provide more flexibility to scientists working in the relevant field to test their simulation data.

Well-posedness of the Einstein-Euler system in asymptotically flat spacetimes: The constraint equations

NASA Astrophysics Data System (ADS)

Brauer, Uwe; Karp, Lavi

This paper deals with the construction of initial data for the coupled Einstein-Euler system. We consider the condition where the energy density might vanish or tend to zero at infinity, and where the pressure is a fractional power of the energy density. In order to achieve our goals we use a type of weighted Sobolev space of fractional order. The common Lichnerowicz-York scaling method (Choquet-Bruhat and York, 1980 [9]; Cantor, 1979 [7]) for solving the constraint equations cannot be applied here directly. The basic problem is that the matter sources are scaled conformally and the fluid variables have to be recovered from the conformally transformed matter sources. This problem has been addressed, although in a different context, by Dain and Nagy (2002) [11]. We show that if the matter variables are restricted to a certain region, then the Einstein constraint equations have a unique solution in the weighted Sobolev spaces of fractional order. The regularity depends upon the fractional power of the equation of state.

Approximate solution of space and time fractional higher order phase field equation

NASA Astrophysics Data System (ADS)

Shamseldeen, S.

2018-03-01

This paper is concerned with a class of space and time fractional partial differential equation (STFDE) with Riesz derivative in space and Caputo in time. The proposed STFDE is considered as a generalization of a sixth-order partial phase field equation. We describe the application of the optimal homotopy analysis method (OHAM) to obtain an approximate solution for the suggested fractional initial value problem. An averaged-squared residual error function is defined and used to determine the optimal convergence control parameter. Two numerical examples are studied, considering periodic and non-periodic initial conditions, to justify the efficiency and the accuracy of the adopted iterative approach. The dependence of the solution on the order of the fractional derivative in space and time and model parameters is investigated.

Application of fractional derivative with exponential law to bi-fractional-order wave equation with frictional memory kernel

NASA Astrophysics Data System (ADS)

Cuahutenango-Barro, B.; Taneco-Hernández, M. A.; Gómez-Aguilar, J. F.

2017-12-01

Analytical solutions of the wave equation with bi-fractional-order and frictional memory kernel of Mittag-Leffler type are obtained via Caputo-Fabrizio fractional derivative in the Liouville-Caputo sense. Through the method of separation of variables and Laplace transform method we derive closed-form solutions and establish fundamental solutions. Special cases with homogeneous Dirichlet boundary conditions and nonhomogeneous initial conditions, as well as for the external force are considered. Numerical simulations of the special solutions were done and novel behaviors are obtained.

A modification of \\mathsf {WKB} method for fractional differential operators of Schrödinger's type

NASA Astrophysics Data System (ADS)

Sayevand, K.; Pichaghchi, K.

2017-09-01

In this paper, we were concerned with the description of the singularly perturbed differential equations within the scope of fractional calculus. However, we shall note that one of the main methods used to solve these problems is the so-called WKB method. We should mention that this was not achievable via the existing fractional derivative definitions, because they do not obey the chain rule. In order to accommodate the WKB to the scope of fractional derivative, we proposed a relatively new derivative called the local fractional derivative. By use of properties of local fractional derivative, we extend the WKB method in the scope of the fractional differential equation. By means of this extension, the WKB analysis based on the Borel resummation, for fractional differential operators of WKB type are investigated. The convergence and the Mittag-Leffler stability of the proposed approach is proven. The obtained results are in excellent agreement with the existing ones in open literature and it is shown that the present approach is very effective and accurate. Furthermore, we are mainly interested to construct the solution of fractional Schrödinger equation in the Mittag-Leffler form and how it leads naturally to this semi-classical approximation namely modified WKB.

A Simple Mathematical Model Inspired by the Purkinje Cells: From Delayed Travelling Waves to Fractional Diffusion.

PubMed

Dipierro, Serena; Valdinoci, Enrico

2018-07-01

Recently, several experiments have demonstrated the existence of fractional diffusion in the neuronal transmission occurring in the Purkinje cells, whose malfunctioning is known to be related to the lack of voluntary coordination and the appearance of tremors. Also, a classical mathematical feature is that (fractional) parabolic equations possess smoothing effects, in contrast with the case of hyperbolic equations, which typically exhibit shocks and discontinuities. In this paper, we show how a simple toy-model of a highly ramified structure, somehow inspired by that of the Purkinje cells, may produce a fractional diffusion via the superposition of travelling waves that solve a hyperbolic equation. This could suggest that the high ramification of the Purkinje cells might have provided an evolutionary advantage of "smoothing" the transmission of signals and avoiding shock propagations (at the price of slowing a bit such transmission). Although an experimental confirmation of the possibility of such evolutionary advantage goes well beyond the goals of this paper, we think that it is intriguing, as a mathematical counterpart, to consider the time fractional diffusion as arising from the superposition of delayed travelling waves in highly ramified transmission media. The case of a travelling concave parabola with sufficiently small curvature is explicitly computed. The new link that we propose between time fractional diffusion and hyperbolic equation also provides a novelty with respect to the usual paradigm relating time fractional diffusion with parabolic equations in the limit. This paper is written in such a way as to be of interest to both biologists and mathematician alike. In order to accomplish this aim, both complete explanations of the objects considered and detailed lists of references are provided.

Numerical study of fractional nonlinear Schrödinger equations.

PubMed

Klein, Christian; Sparber, Christof; Markowich, Peter

2014-12-08

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation.

Numerical study of fractional nonlinear Schrödinger equations

PubMed Central

Klein, Christian; Sparber, Christof; Markowich, Peter

2014-01-01

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation. PMID:25484604

Predicting fractional bed load transport rates: Application of the Wilcock‐Crowe equations to a regulated gravel bed river

USGS Publications Warehouse

Gaeuman, David; Andrews, E.D.; Krause, Andreas; Smith, Wes

2009-01-01

Bed load samples from four locations in the Trinity River of northern California are analyzed to evaluate the performance of the Wilcock‐Crowe bed load transport equations for predicting fractional bed load transport rates. Bed surface particles become smaller and the fraction of sand on the bed increases with distance downstream from Lewiston Dam. The dimensionless reference shear stress for the mean bed particle size (τ*rm) is largest near the dam, but varies relatively little between the more downstream locations. The relation between τ*rm and the reference shear stresses for other size fractions is constant across all locations. Total bed load transport rates predicted with the Wilcock‐Crowe equations are within a factor of 2 of sampled transport rates for 68% of all samples. The Wilcock‐Crowe equations nonetheless consistently under‐predict the transport of particles larger than 128 mm, frequently by more than an order of magnitude. Accurate prediction of the transport rates of the largest particles is important for models in which the evolution of the surface grain size distribution determines subsequent bed load transport rates. Values of τ*rm estimated from bed load samples are up to 50% larger than those predicted with the Wilcock‐Crowe equations, and sampled bed load transport approximates equal mobility across a wider range of grain sizes than is implied by the equations. Modifications to the Wilcock‐Crowe equation for determining τ*rm and the hiding function used to scale τ*rm to other grain size fractions are proposed to achieve the best fit to observed bed load transport in the Trinity River.

Fractional Solitons in Excitonic Josephson Junctions.

PubMed

Hsu, Ya-Fen; Su, Jung-Jung

2015-10-29

The Josephson effect is especially appealing to physicists because it reveals macroscopically the quantum order and phase. In excitonic bilayers the effect is even subtler due to the counterflow of supercurrent as well as the tunneling between layers (interlayer tunneling). Here we study, in a quantum Hall bilayer, the excitonic Josephson junction: a conjunct of two exciton condensates with a relative phase ϕ0 applied. The system is mapped into a pseudospin ferromagnet then described numerically by the Landau-Lifshitz-Gilbert equation. In the presence of interlayer tunneling, we identify a family of fractional sine-Gordon solitons which resemble the static fractional Josephson vortices in the extended superconducting Josephson junctions. Each fractional soliton carries a topological charge Q that is not necessarily a half/full integer but can vary continuously. The calculated current-phase relation (CPR) shows that solitons with Q = ϕ0/2π is the lowest energy state starting from zero ϕ0 - until ϕ0 > π - then the alternative group of solitons with Q = ϕ0/2π - 1 takes place and switches the polarity of CPR.

Fractional Solitons in Excitonic Josephson Junctions

NASA Astrophysics Data System (ADS)

Su, Jung-Jung; Hsu, Ya-Fen

The Josephson effect is especially appealing because it reveals macroscopically the quantum order and phase. Here we study this effect in an excitonic Josephson junction: a conjunct of two exciton condensates with a relative phase ϕ0 applied. Such a junction is proposed to take place in the quantum Hall bilayer (QHB) that makes it subtler than in superconductor because of the counterflow of excitonic supercurrent and the interlayer tunneling in QHB. We treat the system theoretically by first mapping it into a pseudospin ferromagnet then describing it by the Landau-Lifshitz-Gilbert equation. In the presence of interlayer tunneling, the excitonic Josephson junction can possess a family of fractional sine-Gordon solitons that resemble the static fractional Josephson vortices in the extended superconducting Josephson junctions. Interestingly, each fractional soliton carries a topological charge Q which is not necessarily a half/full integer but can vary continuously. The resultant current-phase relation (CPR) shows that solitons with Q =ϕ0 / 2 π are the lowest energy states for small ϕ0. When ϕ0 > π , solitons with Q =ϕ0 / 2 π - 1 take place - the polarity of CPR is then switched.

Fractional Solitons in Excitonic Josephson Junctions

NASA Astrophysics Data System (ADS)

Hsu, Ya-Fen; Su, Jung-Jung

2015-10-01

The Josephson effect is especially appealing to physicists because it reveals macroscopically the quantum order and phase. In excitonic bilayers the effect is even subtler due to the counterflow of supercurrent as well as the tunneling between layers (interlayer tunneling). Here we study, in a quantum Hall bilayer, the excitonic Josephson junction: a conjunct of two exciton condensates with a relative phase ϕ0 applied. The system is mapped into a pseudospin ferromagnet then described numerically by the Landau-Lifshitz-Gilbert equation. In the presence of interlayer tunneling, we identify a family of fractional sine-Gordon solitons which resemble the static fractional Josephson vortices in the extended superconducting Josephson junctions. Each fractional soliton carries a topological charge Q that is not necessarily a half/full integer but can vary continuously. The calculated current-phase relation (CPR) shows that solitons with Q = ϕ0/2π is the lowest energy state starting from zero ϕ0 - until ϕ0 > π - then the alternative group of solitons with Q = ϕ0/2π - 1 takes place and switches the polarity of CPR.

Efficient Application of Continuous Fractional Component Monte Carlo in the Reaction Ensemble

PubMed Central

2017-01-01

A new formulation of the Reaction Ensemble Monte Carlo technique (RxMC) combined with the Continuous Fractional Component Monte Carlo method is presented. This method is denoted by serial Rx/CFC. The key ingredient is that fractional molecules of either reactants or reaction products are present and that chemical reactions always involve fractional molecules. Serial Rx/CFC has the following advantages compared to other approaches: (1) One directly obtains chemical potentials of all reactants and reaction products. Obtained chemical potentials can be used directly as an independent check to ensure that chemical equilibrium is achieved. (2) Independent biasing is applied to the fractional molecules of reactants and reaction products. Therefore, the efficiency of the algorithm is significantly increased, compared to the other approaches. (3) Changes in the maximum scaling parameter of intermolecular interactions can be chosen differently for reactants and reaction products. (4) The number of fractional molecules is reduced. As a proof of principle, our method is tested for Lennard-Jones systems at various pressures and for various chemical reactions. Excellent agreement was found both for average densities and equilibrium mixture compositions computed using serial Rx/CFC, RxMC/CFCMC previously introduced by Rosch and Maginn (Journal of Chemical Theory and Computation, 2011, 7, 269–279), and the conventional RxMC approach. The serial Rx/CFC approach is also tested for the reaction of ammonia synthesis at various temperatures and pressures. Excellent agreement was found between results obtained from serial Rx/CFC, experimental results from literature, and thermodynamic modeling using the Peng–Robinson equation of state. The efficiency of reaction trial moves is improved by a factor of 2 to 3 (depending on the system) compared to the RxMC/CFCMC formulation by Rosch and Maginn. PMID:28737933

Partial Fractions via Calculus

ERIC Educational Resources Information Center

Bauldry, William C.

2018-01-01

The standard technique taught in calculus courses for partial fraction expansions uses undetermined coefficients to generate a system of linear equations; we present a derivative-based technique that calculus and differential equations instructors can use to reinforce connections to calculus. Simple algebra shows that we can use the derivative to…

New optical solitons of space-time conformable fractional perturbed Gerdjikov-Ivanov equation by sine-Gordon equation method

NASA Astrophysics Data System (ADS)

Yaşar, Elif; Yıldırım, Yakup; Yaşar, Emrullah

2018-06-01

This paper devotes to conformable fractional space-time perturbed Gerdjikov-Ivanov (GI) equation which appears in nonlinear fiber optics and photonic crystal fibers (PCF). We consider the model with full nonlinearity in order to give a generalized flavor. The sine-Gordon equation approach is carried out to model equation for retrieving the dark, bright, dark-bright, singular and combined singular optical solitons. The constraint conditions are also reported for guaranteeing the existence of these solitons. We also present some graphical simulations of the solutions for better understanding the physical phenomena of the behind the considered model.

Fractional diffusion: recovering the distributed fractional derivative from overposed data

NASA Astrophysics Data System (ADS)

Rundell, W.; Zhang, Z.

2017-03-01

There has been considerable recent study in ‘subdiffusion’ models that replace the standard parabolic equation model by a one with a fractional derivative in the time variable. There are many ways to look at this newer approach and one such is to realize that the order of the fractional derivative is related to the time scales of the underlying diffusion process. This raises the question of what order α of derivative should be taken and if a single value actually suffices. This has led to models that combine a finite number of these derivatives each with a different fractional exponent {αk} and different weighting value c k to better model a greater possible range of time scales. Ultimately, one wants to look at a situation that combines derivatives in a continuous way—the so-called distributional model with parameter μ ≤ft(α \\right) . However all of this begs the question of how one determines this ‘order’ of differentiation. Recovering a single fractional value has been an active part of the process from the beginning of fractional diffusion modeling and if this is the only unknown then the markers left by the fractional order derivative are relatively straightforward to determine. In the case of a finite combination of derivatives this becomes much more complex due to the more limited analytic tools available for such equations, but recent progress in this direction has been made, (Li et al 2015 Appl. Math. Comput. 257 381-97, Li and Yamamoto 2015 Appl. Anal. 94 570-9). This paper considers the full distributional model where the order is viewed as a function μ ≤ft(α \\right) on the interval (0, 1]. We show existence, uniqueness and regularity for an initial-boundary value problem including an important representation theorem in the case of a single spatial variable. This is then used in the inverse problem of recovering the distributional coefficient μ ≤ft(α \\right) from a time trace of the solution and a uniqueness result is proven.

New Operational Matrices for Solving Fractional Differential Equations on the Half-Line

PubMed Central

2015-01-01

In this paper, the fractional-order generalized Laguerre operational matrices (FGLOM) of fractional derivatives and fractional integration are derived. These operational matrices are used together with spectral tau method for solving linear fractional differential equations (FDEs) of order ν (0 < ν < 1) on the half line. An upper bound of the absolute errors is obtained for the approximate and exact solutions. Fractional-order generalized Laguerre pseudo-spectral approximation is investigated for solving nonlinear initial value problem of fractional order ν. The extension of the fractional-order generalized Laguerre pseudo-spectral method is given to solve systems of FDEs. We present the advantages of using the spectral schemes based on fractional-order generalized Laguerre functions and compare them with other methods. Several numerical examples are implemented for FDEs and systems of FDEs including linear and nonlinear terms. We demonstrate the high accuracy and the efficiency of the proposed techniques. PMID:25996369

New operational matrices for solving fractional differential equations on the half-line.

PubMed

Bhrawy, Ali H; Taha, Taha M; Alzahrani, Ebraheem O; Alzahrani, Ebrahim O; Baleanu, Dumitru; Alzahrani, Abdulrahim A

2015-01-01

In this paper, the fractional-order generalized Laguerre operational matrices (FGLOM) of fractional derivatives and fractional integration are derived. These operational matrices are used together with spectral tau method for solving linear fractional differential equations (FDEs) of order ν (0 < ν < 1) on the half line. An upper bound of the absolute errors is obtained for the approximate and exact solutions. Fractional-order generalized Laguerre pseudo-spectral approximation is investigated for solving nonlinear initial value problem of fractional order ν. The extension of the fractional-order generalized Laguerre pseudo-spectral method is given to solve systems of FDEs. We present the advantages of using the spectral schemes based on fractional-order generalized Laguerre functions and compare them with other methods. Several numerical examples are implemented for FDEs and systems of FDEs including linear and nonlinear terms. We demonstrate the high accuracy and the efficiency of the proposed techniques.

Generalized hydrodynamic correlations and fractional memory functions

NASA Astrophysics Data System (ADS)

Rodríguez, Rosalio F.; Fujioka, Jorge

2015-12-01

A fractional generalized hydrodynamic (GH) model of the longitudinal velocity fluctuations correlation, and its associated memory function, for a complex fluid is analyzed. The adiabatic elimination of fast variables introduces memory effects in the transport equations, and the dynamic of the fluctuations is described by a generalized Langevin equation with long-range noise correlations. These features motivate the introduction of Caputo time fractional derivatives and allows us to calculate analytic expressions for the fractional longitudinal velocity correlation function and its associated memory function. Our analysis eliminates a spurious constant term in the non-fractional memory function found in the non-fractional description. It also produces a significantly slower power-law decay of the memory function in the GH regime that reduces to the well-known exponential decay in the non-fractional Navier-Stokes limit.

Evaporative fractionation of marine water isotopes in the Arctic Ocean help understand a changing Arctic water cycle

NASA Astrophysics Data System (ADS)

Klein, E. S.; Welker, J. M.

2017-12-01

Most of the global hydrologic cycle occurs in oceanic waters. This oceanic derived moisture is critical to the precipitation and evapotranspiration regimes that influence terrestrial Earth systems. Thus understanding oceanic water processes has important global implications for our knowledge of modern and past hydrologic cycles. As they are influenced by environmental variables such as sea surface temperature and atmospheric humidity, water isotope ratios (e.g., δ18O, δ2H) can help understand the patterns driving the water cycle. However, our knowledge of marine isotopes is relatively limited. In particular, the fractionation of water isotopes during evaporation of oceanic water, essentially the start of the hydrologic cycle, is largely based on theoretical relationships derived from spatially and temporally limited data sets. This constrained understanding of oceanic evaporation fractionation patterns is especially pronounced in the rapidly changing Arctic Ocean. These changes are associated with reduced sea ice coverage, which is increasing the amount of local Artic Ocean sourced moisture in atmospheric and terrestrial systems and amplifying the Arctic hydrologic cycle. Here we present new data revealing the nuances of evaporative fractionation of Arctic Ocean water isotopes with the first collection of continuous, contemporaneous sea water and vapor isotopes. These data, collected in situ aboard the icebreaker Healy, show that the difference between actual ocean vapor isotope values and vapor values estimated by the closure equation increases progressively with latitude (especially beyond 70°) and varies between δ18O and δ2H. These differences are likely due to more isotopic mixing in the troposphere and/or closure equation assumptions inapplicable to Arctic regions. Moreover, we find: 1) a positive relationship between fractionation magnitude and latitude; and 2) the influence of evaporative fractionation from environmental variables such as wind and relative humidity reverses with the presence of sea ice. These new data increase our understanding of the patterns and processes governing past, present, and future changes to the Arctic hydrologic cycle.

Modeling electro-magneto-hydrodynamic thermo-fluidic transport of biofluids with new trend of fractional derivative without singular kernel

NASA Astrophysics Data System (ADS)

Abdulhameed, M.; Vieru, D.; Roslan, R.

2017-10-01

This paper investigates the electro-magneto-hydrodynamic flow of the non-Newtonian behavior of biofluids, with heat transfer, through a cylindrical microchannel. The fluid is acted by an arbitrary time-dependent pressure gradient, an external electric field and an external magnetic field. The governing equations are considered as fractional partial differential equations based on the Caputo-Fabrizio time-fractional derivatives without singular kernel. The usefulness of fractional calculus to study fluid flows or heat and mass transfer phenomena was proven. Several experimental measurements led to conclusion that, in such problems, the models described by fractional differential equations are more suitable. The most common time-fractional derivative used in Continuum Mechanics is Caputo derivative. However, two disadvantages appear when this derivative is used. First, the definition kernel is a singular function and, secondly, the analytical expressions of the problem solutions are expressed by generalized functions (Mittag-Leffler, Lorenzo-Hartley, Robotnov, etc.) which, generally, are not adequate to numerical calculations. The new time-fractional derivative Caputo-Fabrizio, without singular kernel, is more suitable to solve various theoretical and practical problems which involve fractional differential equations. Using the Caputo-Fabrizio derivative, calculations are simpler and, the obtained solutions are expressed by elementary functions. Analytical solutions of the biofluid velocity and thermal transport are obtained by means of the Laplace and finite Hankel transforms. The influence of the fractional parameter, Eckert number and Joule heating parameter on the biofluid velocity and thermal transport are numerically analyzed and graphic presented. This fact can be an important in Biochip technology, thus making it possible to use this analysis technique extremely effective to control bioliquid samples of nanovolumes in microfluidic devices used for biological analysis and medical diagnosis.

Comparison of stochastic lung deposition fractions with experimental data.

PubMed

Majid, Hussain; Hofmann, Werner; Winkler-Heil, Renate

2012-04-01

Deposition fractions of inhaled particles predicted by different computational models vary with respect to physical and biological factors and mathematical modeling techniques. These models must be validated by comparison with available experimental data. Experimental data supplied by different deposition studies with surrogate airway models or lung casts were used in this study to evaluate the stochastic deposition model Inhalation, Deposition and Exhalation of Aerosols in the Lung at the airway generation level. Furthermore, different analytical equations derived for the three major deposition mechanisms, diffusion, impaction, and sedimentation, were applied to different cast or airway models to quantify their effect on calculated particle deposition fractions. The experimental results for ultrafine particles (0.00175 and 0.01) were found to be in close agreement with the stochastic model predictions; however, for coarse particles (3 and 8 μm), experimental deposition fractions became higher with increasing flow rate. An overall fair agreement among the calculated deposition fractions for the different cast geometries was found. However, alternative deposition equations resulted in up to 300% variation in predicted deposition fractions, although all equations predicted the same trends as functions of particle diameter and breathing conditions. From this comparative study, it can be concluded that structural differences in lung morphologies among different individuals are responsible for the apparent variability in particle deposition in each generation. The use of different deposition equations yields varying deposition results caused primarily by (i) different lung morphometries employed in their derivation and the choice of the central bifurcation zone geometry, (ii) the assumption of specific flow profiles, and (iii) different methods used in the derivation of these equations.

Unveiling the relationships among the viscosity equations of glass liquids and colloidal suspensions for obtaining universal equations with the generic free volume concept.

PubMed

Hao, Tian

2015-09-14

The underlying relationships among viscosity equations of glass liquids and colloidal suspensions are explored with the aid of free volume concept. Viscosity equations of glass liquids available in literature are focused and found to have a same physical basis but different mathematical expressions for the free volume. The glass transitions induced by temperatures in glass liquids and the percolation transition induced by particle volume fractions in colloidal suspensions essentially are a second order phase transition: both those two transitions could induce the free volume changes, which in turn determines how the viscosities are going to change with temperatures and/or particle volume fractions. Unified correlations of the free volume to both temperatures and particle volume fractions are thus proposed. The resulted viscosity equations are reducible to many popular viscosity equations currently widely used in literature; those equations should be able to cover many different types of materials over a wide temperature range. For demonstration purpose, one of the simplified versions of those newly developed equations is compared with popular viscosity equations and the experimental data: it can well fit the experimental data over a wide temperature range. The current work reveals common physical grounds among various viscosity equations, deepening our understanding on viscosity and unifying the free volume theory across many different systems.

Boundedness of the solutions for certain classes of fractional differential equations with application to adaptive systems.

PubMed

Aguila-Camacho, Norelys; Duarte-Mermoud, Manuel A

2016-01-01

This paper presents the analysis of three classes of fractional differential equations appearing in the field of fractional adaptive systems, for the case when the fractional order is in the interval α ∈(0,1] and the Caputo definition for fractional derivatives is used. The boundedness of the solutions is proved for all three cases, and the convergence to zero of the mean value of one of the variables is also proved. Applications of the obtained results to fractional adaptive schemes in the context of identification and control problems are presented at the end of the paper, including numerical simulations which support the analytical results. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Convective flows of generalized time-nonlocal nanofluids through a vertical rectangular channel

NASA Astrophysics Data System (ADS)

Ahmed, Najma; Vieru, Dumitru; Fetecau, Constantin; Shah, Nehad Ali

2018-05-01

Time-nonlocal generalized model of the natural convection heat transfer and nanofluid flows through a rectangular vertical channel with wall conditions of the Robin type are studied. The generalized mathematical model with time-nonlocality is developed by considering the fractional constitutive equations for the shear stress and thermal flux defined with the time-fractional Caputo derivative. The Caputo power-law non-local kernel provides the damping to the velocity and temperature gradient; therefore, transport processes are influenced by the histories at all past and present times. Analytical solutions for dimensionless velocity and temperature fields are obtained by using the Laplace transform coupled with the finite sine-cosine Fourier transform which is suitable to problems with boundary conditions of the Robin type. Particularizing the fractional thermal and velocity parameters, solutions for three simplified models are obtained (classical linear momentum equation with damped thermal flux; fractional shear stress constitutive equation with classical Fourier's law for thermal flux; classical shear stress and thermal flux constitutive equations). It is found that the thermal histories strongly influence the thermal transport for small values of time t. Also, the thermal transport can be enhanced if the thermal fractional parameter decreases or by increasing the nanoparticles' volume fraction. The velocity field is influenced on the one hand by the temperature of the fluid and on the other by the damping of the velocity gradient introduced by the fractional derivative. Also, the transport motions of the channel walls influence the motion of the fluid layers located near them.

Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions.

PubMed

Langlands, T A M; Henry, B I; Wearne, S L

2009-12-01

We introduce fractional Nernst-Planck equations and derive fractional cable equations as macroscopic models for electrodiffusion of ions in nerve cells when molecular diffusion is anomalous subdiffusion due to binding, crowding or trapping. The anomalous subdiffusion is modelled by replacing diffusion constants with time dependent operators parameterized by fractional order exponents. Solutions are obtained as functions of the scaling parameters for infinite cables and semi-infinite cables with instantaneous current injections. Voltage attenuation along dendrites in response to alpha function synaptic inputs is computed. Action potential firing rates are also derived based on simple integrate and fire versions of the models. Our results show that electrotonic properties and firing rates of nerve cells are altered by anomalous subdiffusion in these models. We have suggested electrophysiological experiments to calibrate and validate the models.

Multidimensional fractional Schrödinger equation

NASA Astrophysics Data System (ADS)

Rodrigues, M. M.; Vieira, N.

2012-11-01

This work is intended to investigate the multi-dimensional space-time fractional Schrödinger equation of the form (CDt0+αu)(t,x) = iħ/2m(C∇βu)(t,x), with ħ the Planck's constant divided by 2π, m is the mass and u(t,x) is a wave function of the particle. Here (CDt0+α,C∇β are operators of the Caputo fractional derivatives, where α ∈]0,1] and β ∈]1,2]. The wave function is obtained using Laplace and Fourier transforms methods and a symbolic operational form of solutions in terms of the Mittag-Leffler functions is exhibited. It is presented an expression for the wave function and for the quantum mechanical probability density. Using Banach fixed point theorem, the existence and uniqueness of solutions is studied for this kind of fractional differential equations.

Global uniqueness in an inverse problem for time fractional diffusion equations

NASA Astrophysics Data System (ADS)

Kian, Y.; Oksanen, L.; Soccorsi, E.; Yamamoto, M.

2018-01-01

Given (M , g), a compact connected Riemannian manifold of dimension d ⩾ 2, with boundary ∂M, we consider an initial boundary value problem for a fractional diffusion equation on (0 , T) × M, T > 0, with time-fractional Caputo derivative of order α ∈ (0 , 1) ∪ (1 , 2). We prove uniqueness in the inverse problem of determining the smooth manifold (M , g) (up to an isometry), and various time-independent smooth coefficients appearing in this equation, from measurements of the solutions on a subset of ∂M at fixed time. In the "flat" case where M is a compact subset of Rd, two out the three coefficients ρ (density), a (conductivity) and q (potential) appearing in the equation ρ ∂tα u - div (a∇u) + qu = 0 on (0 , T) × M are recovered simultaneously.

A unifying fractional wave equation for compressional and shear waves.

PubMed

Holm, Sverre; Sinkus, Ralph

2010-01-01

This study has been motivated by the observed difference in the range of the power-law attenuation exponent for compressional and shear waves. Usually compressional attenuation increases with frequency to a power between 1 and 2, while shear wave attenuation often is described with powers less than 1. Another motivation is the apparent lack of partial differential equations with desirable properties such as causality that describe such wave propagation. Starting with a constitutive equation which is a generalized Hooke's law with a loss term containing a fractional derivative, one can derive a causal fractional wave equation previously given by Caputo [Geophys J. R. Astron. Soc. 13, 529-539 (1967)] and Wismer [J. Acoust. Soc. Am. 120, 3493-3502 (2006)]. In the low omegatau (low-frequency) case, this equation has an attenuation with a power-law in the range from 1 to 2. This is consistent with, e.g., attenuation in tissue. In the often neglected high omegatau (high-frequency) case, it describes attenuation with a power-law between 0 and 1, consistent with what is observed in, e.g., dynamic elastography. Thus a unifying wave equation derived properly from constitutive equations can describe both cases.

Fractional vector calculus for fractional advection dispersion

NASA Astrophysics Data System (ADS)

Meerschaert, Mark M.; Mortensen, Jeff; Wheatcraft, Stephen W.

2006-07-01

We develop the basic tools of fractional vector calculus including a fractional derivative version of the gradient, divergence, and curl, and a fractional divergence theorem and Stokes theorem. These basic tools are then applied to provide a physical explanation for the fractional advection-dispersion equation for flow in heterogeneous porous media.

Fractional cable equation for general geometry: A model of axons with swellings and anomalous diffusion.

PubMed

López-Sánchez, Erick J; Romero, Juan M; Yépez-Martínez, Huitzilin

2017-09-01

Different experimental studies have reported anomalous diffusion in brain tissues and notably this anomalous diffusion is expressed through fractional derivatives. Axons are important to understand neurodegenerative diseases such as multiple sclerosis, Alzheimer's disease, and Parkinson's disease. Indeed, abnormal accumulation of proteins and organelles in axons is a hallmark of these diseases. The diffusion in the axons can become anomalous as a result of this abnormality. In this case the voltage propagation in axons is affected. Another hallmark of different neurodegenerative diseases is given by discrete swellings along the axon. In order to model the voltage propagation in axons with anomalous diffusion and swellings, in this paper we propose a fractional cable equation for a general geometry. This generalized equation depends on fractional parameters and geometric quantities such as the curvature and torsion of the cable. For a cable with a constant radius we show that the voltage decreases when the fractional effect increases. In cables with swellings we find that when the fractional effect or the swelling radius increases, the voltage decreases. Similar behavior is obtained when the number of swellings and the fractional effect increase. Moreover, we find that when the radius swelling (or the number of swellings) and the fractional effect increase at the same time, the voltage dramatically decreases.

Fractional cable equation for general geometry: A model of axons with swellings and anomalous diffusion

NASA Astrophysics Data System (ADS)

López-Sánchez, Erick J.; Romero, Juan M.; Yépez-Martínez, Huitzilin

2017-09-01

Different experimental studies have reported anomalous diffusion in brain tissues and notably this anomalous diffusion is expressed through fractional derivatives. Axons are important to understand neurodegenerative diseases such as multiple sclerosis, Alzheimer's disease, and Parkinson's disease. Indeed, abnormal accumulation of proteins and organelles in axons is a hallmark of these diseases. The diffusion in the axons can become anomalous as a result of this abnormality. In this case the voltage propagation in axons is affected. Another hallmark of different neurodegenerative diseases is given by discrete swellings along the axon. In order to model the voltage propagation in axons with anomalous diffusion and swellings, in this paper we propose a fractional cable equation for a general geometry. This generalized equation depends on fractional parameters and geometric quantities such as the curvature and torsion of the cable. For a cable with a constant radius we show that the voltage decreases when the fractional effect increases. In cables with swellings we find that when the fractional effect or the swelling radius increases, the voltage decreases. Similar behavior is obtained when the number of swellings and the fractional effect increase. Moreover, we find that when the radius swelling (or the number of swellings) and the fractional effect increase at the same time, the voltage dramatically decreases.

Novel Numerical Methods for Optimal Control Problems Involving Fractional-Order Differential Equations

DTIC Science & Technology

2018-03-14

pricing, Appl. Math . Comp. Vol.305, 174-187 (2017) 5. W. Li, S. Wang, Pricing European options with proportional transaction costs and stochastic...for fractional differential equation. Numer. Math . Theor. Methods Appl. 5, 229–241, 2012. [23] Kilbas A.A. and Marzan, S.A., Cauchy problem for...numerical technique for solving fractional optimal control problems, Comput. Math . Appl., 62, Issue 3, 1055–1067, 2011. [26] Lotfi A., Yousefi SA., Dehghan M

The Vertical Linear Fractional Initialization Problem

NASA Technical Reports Server (NTRS)

Lorenzo, Carl F.; Hartley, Tom T.

1999-01-01

This paper presents a solution to the initialization problem for a system of linear fractional-order differential equations. The scalar problem is considered first, and solutions are obtained both generally and for a specific initialization. Next the vector fractional order differential equation is considered. In this case, the solution is obtained in the form of matrix F-functions. Some control implications of the vector case are discussed. The suggested method of problem solution is shown via an example.

Mathematical treatment of isotopologue and isotopomer speciation and fractionation in biochemical kinetics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Maggi, F.M.; Riley, W.J.

2009-11-01

We present a mathematical treatment of the kinetic equations that describe isotopologue and isotopomer speciation and fractionation during enzyme-catalyzed biochemical reactions. These equations, presented here with the name GEBIK (general equations for biochemical isotope kinetics) and GEBIF (general equations for biochemical isotope fractionation), take into account microbial biomass and enzyme dynamics, reaction stoichiometry, isotope substitution number, and isotope location within each isotopologue and isotopomer. In addition to solving the complete GEBIK and GEBIF, we also present and discuss two approximations to the full solutions under the assumption of biomass-free and enzyme steady-state, and under the quasi-steady-state assumption as applied tomore » the complexation rate. The complete and approximate approaches are applied to observations of biological denitrification in soils. Our analysis highlights that the full GEBIK and GEBIF provide a more accurate description of concentrations and isotopic compositions of substrates and products throughout the reaction than do the approximate forms. We demonstrate that the isotopic effects of a biochemical reaction depend, in the most general case, on substrate and complex concentrations and, therefore, the fractionation factor is a function of time. We also demonstrate that inverse isotopic effects can occur for values of the fractionation factor smaller than 1, and that reactions that do not discriminate isotopes do not necessarily imply a fractionation factor equal to 1.« less

Fractional Stochastic Differential Equations Satisfying Fluctuation-Dissipation Theorem

NASA Astrophysics Data System (ADS)

Li, Lei; Liu, Jian-Guo; Lu, Jianfeng

2017-10-01

We propose in this work a fractional stochastic differential equation (FSDE) model consistent with the over-damped limit of the generalized Langevin equation model. As a result of the `fluctuation-dissipation theorem', the differential equations driven by fractional Brownian noise to model memory effects should be paired with Caputo derivatives, and this FSDE model should be understood in an integral form. We establish the existence of strong solutions for such equations and discuss the ergodicity and convergence to Gibbs measure. In the linear forcing regime, we show rigorously the algebraic convergence to Gibbs measure when the `fluctuation-dissipation theorem' is satisfied, and this verifies that satisfying `fluctuation-dissipation theorem' indeed leads to the correct physical behavior. We further discuss possible approaches to analyze the ergodicity and convergence to Gibbs measure in the nonlinear forcing regime, while leave the rigorous analysis for future works. The FSDE model proposed is suitable for systems in contact with heat bath with power-law kernel and subdiffusion behaviors.

Ground state sign-changing solutions for fractional Kirchhoff equations in bounded domains

NASA Astrophysics Data System (ADS)

Luo, Huxiao; Tang, Xianhua; Gao, Zu

2018-03-01

We study the existence of ground state sign-changing solutions for the fractional Kirchhoff problem. Under mild assumptions on the nonlinearity, by using some new analytical skills and the non-Nehari manifold method, we prove that the fractional Kirchhoff problem possesses a ground state sign-changing solution ub. Moreover, we show that the energy of ub is strictly larger than twice that of the ground state solutions of Nehari-type. Finally, we establish the convergence property of ub as the parameter b ↘ 0. Our results generalize some results obtained by Shuai [J. Differ. Equations 259, 1256 (2015)] and Tang and Cheng [J. Differ. Equations 261, 2384 (2016)].

Multiplicity and asymptotic behavior of solutions to a class of Kirchhoff-type equations involving the fractional p-Laplacian.

PubMed

Shen, Liejun

2018-01-01

The present study is concerned with the following fractional p -Laplacian equation involving a critical Sobolev exponent of Kirchhoff type: [Formula: see text] where [Formula: see text], [Formula: see text] and [Formula: see text] are constants, and [Formula: see text] is the fractional p -Laplacian operator with [Formula: see text] and [Formula: see text]. For suitable [Formula: see text], the above equation possesses at least two nontrivial solutions by variational method for any [Formula: see text]. Moreover, we regard [Formula: see text] and [Formula: see text] as parameters to obtain convergent properties of solutions for the given problem as [Formula: see text] and [Formula: see text], respectively.

(N+1)-dimensional fractional reduced differential transform method for fractional order partial differential equations

NASA Astrophysics Data System (ADS)

Arshad, Muhammad; Lu, Dianchen; Wang, Jun

2017-07-01

In this paper, we pursue the general form of the fractional reduced differential transform method (DTM) to (N+1)-dimensional case, so that fractional order partial differential equations (PDEs) can be resolved effectively. The most distinct aspect of this method is that no prescribed assumptions are required, and the huge computational exertion is reduced and round-off errors are also evaded. We utilize the proposed scheme on some initial value problems and approximate numerical solutions of linear and nonlinear time fractional PDEs are obtained, which shows that the method is highly accurate and simple to apply. The proposed technique is thus an influential technique for solving the fractional PDEs and fractional order problems occurring in the field of engineering, physics etc. Numerical results are obtained for verification and demonstration purpose by using Mathematica software.

Societal Statistics by virtue of the Statistical Drake Equation

NASA Astrophysics Data System (ADS)

Maccone, Claudio

2012-09-01

The Drake equation, first proposed by Frank D. Drake in 1961, is the foundational equation of SETI. It yields an estimate of the number N of extraterrestrial communicating civilizations in the Galaxy given by the product N=Ns×fp×ne×fl×fi×fc×fL, where: Ns is the number of stars in the Milky Way Galaxy; fp is the fraction of stars that have planetary systems; ne is the number of planets in a given system that are ecologically suitable for life; fl is the fraction of otherwise suitable planets on which life actually arises; fi is the fraction of inhabited planets on which an intelligent form of life evolves; fc is the fraction of planets inhabited by intelligent beings on which a communicative technical civilization develops; and fL is the fraction of planetary lifetime graced by a technical civilization. The first three terms may be called "the astrophysical terms" in the Drake equation since their numerical value is provided by astrophysical considerations. The fourth term, fl, may be called "the origin-of-life term" and entails biology. The last three terms may be called "the societal terms" inasmuch as their respective numerical values are provided by anthropology, telecommunication science and "futuristic science", respectively. In this paper, we seek to provide a statistical estimate of the three societal terms in the Drake equation basing our calculations on the Statistical Drake Equation first proposed by this author at the 2008 IAC. In that paper the author extended the simple 7-factor product so as to embody Statistics. He proved that, no matter which probability distribution may be assigned to each factor, if the number of factors tends to infinity, then the random variable N follows the lognormal distribution (central limit theorem of Statistics). This author also proved at the 2009 IAC that the Dole (1964) [7] equation, yielding the number of Habitable Planets for Man in the Galaxy, has the same mathematical structure as the Drake equation. So the number of Habitable Planets follows the lognormal distribution as well. But the Dole equation is described by the first FOUR factors of the Drake equation. Thus, we may "divide" the 7-factor Drake equation by the 4-factor Dole equation getting the probability distribution of the last-3-factor Drake equation, i.e. the probability distribution of the SOCIETAL TERMS ONLY. These we study in detail in this paper, achieving new statistical results about the SOCIETAL ASPECTS OF SETI.

Fractional bosonic strings

NASA Astrophysics Data System (ADS)

Diaz, Victor Alfonzo; Giusti, Andrea

2018-03-01

The aim of this paper is to present a simple generalization of bosonic string theory in the framework of the theory of fractional variational problems. Specifically, we present a fractional extension of the Polyakov action, for which we compute the general form of the equations of motion and discuss the connection between the new fractional action and a generalization the Nambu-Goto action. Consequently, we analyze the symmetries of the modified Polyakov action and try to fix the gauge, following the classical procedures. Then we solve the equations of motion in a simplified setting. Finally, we present a Hamiltonian description of the classical fractional bosonic string and introduce the fractional light-cone gauge. It is important to remark that, throughout the whole paper, we thoroughly discuss how to recover the known results as an "integer" limit of the presented model.

R-Function Relationships for Application in the Fractional Calculus

NASA Technical Reports Server (NTRS)

Lorenzo, Carl F.; Hartley, Tom T.

2000-01-01

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, e(t), and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, el, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

Spectral iterative method and convergence analysis for solving nonlinear fractional differential equation

NASA Astrophysics Data System (ADS)

Yarmohammadi, M.; Javadi, S.; Babolian, E.

2018-04-01

In this study a new spectral iterative method (SIM) based on fractional interpolation is presented for solving nonlinear fractional differential equations (FDEs) involving Caputo derivative. This method is equipped with a pre-algorithm to find the singularity index of solution of the problem. This pre-algorithm gives us a real parameter as the index of the fractional interpolation basis, for which the SIM achieves the highest order of convergence. In comparison with some recent results about the error estimates for fractional approximations, a more accurate convergence rate has been attained. We have also proposed the order of convergence for fractional interpolation error under the L2-norm. Finally, general error analysis of SIM has been considered. The numerical results clearly demonstrate the capability of the proposed method.

R-function relationships for application in the fractional calculus.

PubMed

Lorenzo, Carl F; Hartley, Tom T

2008-01-01

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, et, and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, et, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

Extending Maxwell's equations for dielectric materials using analytical principles from viscoelasticity based on the fractional calculus

NASA Astrophysics Data System (ADS)

Wharmby, Andrew William

Existing fractional calculus models having a non-empirical basis used to describe constitutive relationships between stress and strain in viscoelastic materials are modified to employ all orders of fractional derivatives between zero and one. Parallels between viscoelastic and dielectric theory are drawn so that these modified fractional calculus based models for viscoelastic materials may be used to describe relationships between electric flux density and electric field intensity in dielectric materials. The resulting fractional calculus based dielectric relaxation model is tested using existing complex permittivity data in the radio-frequency bandwidth of a wide variety of homogeneous materials. The consequences that the application of this newly developed fractional calculus based dielectric relaxation model has on Maxwell's equations are also examined through the effects of dielectric dissipation and dispersion.

A Tutorial Review on Fractal Spacetime and Fractional Calculus

NASA Astrophysics Data System (ADS)

He, Ji-Huan

2014-11-01

This tutorial review of fractal-Cantorian spacetime and fractional calculus begins with Leibniz's notation for derivative without limits which can be generalized to discontinuous media like fractal derivative and q-derivative of quantum calculus. Fractal spacetime is used to elucidate some basic properties of fractal which is the foundation of fractional calculus, and El Naschie's mass-energy equation for the dark energy. The variational iteration method is used to introduce the definition of fractional derivatives. Fractal derivative is explained geometrically and q-derivative is motivated by quantum mechanics. Some effective analytical approaches to fractional differential equations, e.g., the variational iteration method, the homotopy perturbation method, the exp-function method, the fractional complex transform, and Yang-Laplace transform, are outlined and the main solution processes are given.

Localization and Ballistic Diffusion for the Tempered Fractional Brownian-Langevin Motion

NASA Astrophysics Data System (ADS)

Chen, Yao; Wang, Xudong; Deng, Weihua

2017-10-01

This paper discusses the tempered fractional Brownian motion (tfBm), its ergodicity, and the derivation of the corresponding Fokker-Planck equation. Then we introduce the generalized Langevin equation with the tempered fractional Gaussian noise for a free particle, called tempered fractional Langevin equation (tfLe). While the tfBm displays localization diffusion for the long time limit and for the short time its mean squared displacement (MSD) has the asymptotic form t^{2H}, we show that the asymptotic form of the MSD of the tfLe transits from t^2 (ballistic diffusion for short time) to t^{2-2H}, and then to t^2 (again ballistic diffusion for long time). On the other hand, the overdamped tfLe has the transition of the diffusion type from t^{2-2H} to t^2 (ballistic diffusion). The tfLe with harmonic potential is also considered.

A Solution to the Fundamental Linear Fractional Order Differential Equation

NASA Technical Reports Server (NTRS)

Hartley, Tom T.; Lorenzo, Carl F.

1998-01-01

This paper provides a solution to the fundamental linear fractional order differential equation, namely, (sub c)d(sup q, sub t) + ax(t) = bu(t). The impulse response solution is shown to be a series, named the F-function, which generalizes the normal exponential function. The F-function provides the basis for a qth order "fractional pole". Complex plane behavior is elucidated and a simple example, the inductor terminated semi- infinite lossy line, is used to demonstrate the theory.

On Chaotic Behavior of Temperature Distribution in a Heat Exchanger

NASA Astrophysics Data System (ADS)

Bagyalakshmi, Morachan; Gangadharan, Saisundarakrishnan; Ganesh, Madhu

The objective of this paper is to introduce the notion of fractional derivatives in the energy equations and to study the chaotic nature of the temperature distribution in a heat exchanger with variation of temperature dependent transport properties. The governing fractional partial differential equations are transformed to a set of recurrence relations using fractional differential transform method and solved using inverse transform. The approximate analytical solution obtained by the proposed method has good agreement with the existing results.

Fractional calculus phenomenology in two-dimensional plasma models

NASA Astrophysics Data System (ADS)

Gustafson, Kyle; Del Castillo Negrete, Diego; Dorland, Bill

2006-10-01

Transport processes in confined plasmas for fusion experiments, such as ITER, are not well-understood at the basic level of fully nonlinear, three-dimensional kinetic physics. Turbulent transport is invoked to describe the observed levels in tokamaks, which are orders of magnitude greater than the theoretical predictions. Recent results show the ability of a non-diffusive transport model to describe numerical observations of turbulent transport. For example, resistive MHD modeling of tracer particle transport in pressure-gradient driven turbulence for a three-dimensional plasma reveals that the superdiffusive (2̂˜t^α where α> 1) radial transport in this system is described quantitatively by a fractional diffusion equation Fractional calculus is a generalization involving integro-differential operators, which naturally describe non-local behaviors. Our previous work showed the quantitative agreement of special fractional diffusion equation solutions with numerical tracer particle flows in time-dependent linearized dynamics of the Hasegawa-Mima equation (for poloidal transport in a two-dimensional cold-ion plasma). In pursuit of a fractional diffusion model for transport in a gyrokinetic plasma, we now present numerical results from tracer particle transport in the nonlinear Hasegawa-Mima equation and a planar gyrokinetic model. Finite Larmor radius effects will be discussed. D. del Castillo Negrete, et al, Phys. Rev. Lett. 94, 065003 (2005).

Applications and Implications of Fractional Dynamics for Dielectric Relaxation

NASA Astrophysics Data System (ADS)

Hilfer, R.

This article summarizes briefly the presentation given by the author at the NATO Advanced Research Workshop on "Broadband Dielectric Spectroscopy and its Advanced Technological Applications", held in Perpignan, France, in September 2011. The purpose of the invited presentation at the workshop was to review and summarize the basic theory of fractional dynamics (Hilfer, Phys Rev E 48:2466, 1993; Hilfer and Anton, Phys Rev E Rapid Commun 51:R848, 1995; Hilfer, Fractals 3(1):211, 1995; Hilfer, Chaos Solitons Fractals 5:1475, 1995; Hilfer, Fractals 3:549, 1995; Hilfer, Physica A 221:89, 1995; Hilfer, On fractional diffusion and its relation with continuous time random walks. In: Pekalski et al. (eds) Anomalous diffusion: from basis to applications. Springer, Berlin, p 77, 1999; Hilfer, Fractional evolution equations and irreversibility. In: Helbing et al. (eds) Traffic and granular flow'99. Springer, Berlin, p 215, 2000; Hilfer, Fractional time evolution. In: Hilfer (ed) Applications of fractional calculus in physics. World Scientific, Singapore, p 87, 2000; Hilfer, Remarks on fractional time. In: Castell and Ischebeck (eds) Time, quantum and information. Springer, Berlin, p 235, 2003; Hilfer, Physica A 329:35, 2003; Hilfer, Threefold introduction to fractional derivatives. In: Klages et al. (eds) Anomalous transport: foundations and applications. Wiley-VCH, Weinheim, pp 17-74, 2008; Hilfer, Foundations of fractional dynamics: a short account. In: Klafter et al. (eds) Fractional dynamics: recent advances. World Scientific, Singapore, p 207, 2011) and demonstrate its relevance and application to broadband dielectric spectroscopy (Hilfer, J Phys Condens Matter 14:2297, 2002; Hilfer, Chem Phys 284:399, 2002; Hilfer, Fractals 11:251, 2003; Hilfer et al., Fractional Calc Appl Anal 12:299, 2009). It was argued, that broadband dielectric spectroscopy might be useful to test effective field theories based on fractional dynamics.

Optimal temperature control of tissue embedded with gold nanoparticles for enhanced thermal therapy based on two-energy equation model.

PubMed

Wang, Shen-Ling; Qi, Hong; Ren, Ya-Tao; Chen, Qin; Ruan, Li-Ming

2018-05-01

Thermal therapy is a very promising method for cancer treatment, which can be combined with chemotherapy, radiotherapy and other programs for enhanced cancer treatment. In order to get a better effect of thermal therapy in clinical applications, optimal internal temperature distribution of the tissue embedded with gold nanoparticles (GNPs) for enhanced thermal therapy was investigated in present research. The Monte Carlo method was applied to calculate the heat generation of the tissue embedded with GNPs irradiated by continuous laser. To have a better insight into the physical problem of heat transfer in tissues, the two-energy equation was employed to calculate the temperature distribution of the tissue in the process of GNPs enhanced therapy. The Arrhenius equation was applied to evaluate the degree of permanent thermal damage. A parametric study was performed to investigate the influence factors on the tissue internal temperature distribution, such as incident light intensity, the GNPs volume fraction, the periodic heating and cooling time, and the incident light position. It was found that period heating and cooling strategy can effectively avoid overheating of skin surface and heat damage of healthy tissue. Lower GNPs volume fraction will be better for the heat source distribution. Furthermore, the ring heating strategy is superior to the central heating strategy in the treatment effect. All the analysis provides theoretical guidance for optimal temperature control of tissue embedded with GNP for enhanced thermal therapy. Copyright © 2018 Elsevier Ltd. All rights reserved.

Transverse mode coupling instability threshold with space charge and different wakefields

DOE Office of Scientific and Technical Information (OSTI.GOV)

Balbekov, V.

Transverse mode coupling instability of a bunch with space charge and wake field is considered in frameworks of the boxcar model. Eigenfunctions of the bunch without wake are used as the basis for solution of the equations with the wake field included. Dispersion equation for the bunch eigentunes is obtained in the form of an infinite continued fraction. It is shown that influence of space charge on the instability essentially depends on the wake sign. In particular, threshold of the negative wake increases in absolute value until the space charge tune shift is rather small, and goes to zero atmore » higher space charge. The explanation of this behavior is developed by analysis of the bunch spectrum. As a result, a comparison of the results with published articles is represented.« less

Transverse mode coupling instability threshold with space charge and different wakefields

DOE PAGES

Balbekov, V.

2017-03-10

Transverse mode coupling instability of a bunch with space charge and wake field is considered in frameworks of the boxcar model. Eigenfunctions of the bunch without wake are used as the basis for solution of the equations with the wake field included. Dispersion equation for the bunch eigentunes is obtained in the form of an infinite continued fraction. It is shown that influence of space charge on the instability essentially depends on the wake sign. In particular, threshold of the negative wake increases in absolute value until the space charge tune shift is rather small, and goes to zero atmore » higher space charge. The explanation of this behavior is developed by analysis of the bunch spectrum. As a result, a comparison of the results with published articles is represented.« less

Numerical analysis of fractional MHD Maxwell fluid with the effects of convection heat transfer condition and viscous dissipation

NASA Astrophysics Data System (ADS)

Bai, Yu; Jiang, Yuehua; Liu, Fawang; Zhang, Yan

2017-12-01

This paper investigates the incompressible fractional MHD Maxwell fluid due to a power function accelerating plate with the first order slip, and the numerical analysis on the flow and heat transfer of fractional Maxwell fluid has been done. Moreover the deformation motion of fluid micelle is simply analyzed. Nonlinear velocity equation are formulated with multi-term time fractional derivatives in the boundary layer governing equations, and convective heat transfer boundary condition and viscous dissipation are both taken into consideration. A newly finite difference scheme with L1-algorithm of governing equations are constructed, whose convergence is confirmed by the comparison with analytical solution. Numerical solutions for velocity and temperature show the effects of pertinent parameters on flow and heat transfer of fractional Maxwell fluid. It reveals that the fractional derivative weakens the effects of motion and heat conduction. The larger the Nusselt number is, the greater the heat transfer capacity of fluid becomes, and the temperature gradient at the wall becomes more significantly. The lower Reynolds number enhances the viscosity of the fluid because it is the ratio of the viscous force and the inertia force, which resists the flow and heat transfer.

Derivation of kinetic equations from non-Wiener stochastic differential equations

NASA Astrophysics Data System (ADS)

Basharov, A. M.

2013-12-01

Kinetic differential-difference equations containing terms with fractional derivatives and describing α -stable Levy processes with 0 < α < 1 have been derived in a unified manner in terms of one-dimensional stochastic differential equations controlled merely by the Poisson processes.

Representing Numbers as Continued Fractions and an N-Spire. tns Document to Do Some Basic Continued Fraction Arithmetic

ERIC Educational Resources Information Center

Leinbach, L. Carl

2015-01-01

This paper illustrates a TI N-Spire .tns file created by the author for generating continued fraction representations of real numbers and doing arithmetic with them. The continued fraction representation provides an alternative to the decimal representation. The .tns file can be used as tool for studying continued fractions and their properties as…

Numerical Solution of Time-Dependent Problems with a Fractional-Power Elliptic Operator

NASA Astrophysics Data System (ADS)

Vabishchevich, P. N.

2018-03-01

A time-dependent problem in a bounded domain for a fractional diffusion equation is considered. The first-order evolution equation involves a fractional-power second-order elliptic operator with Robin boundary conditions. A finite-element spatial approximation with an additive approximation of the operator of the problem is used. The time approximation is based on a vector scheme. The transition to a new time level is ensured by solving a sequence of standard elliptic boundary value problems. Numerical results obtained for a two-dimensional model problem are presented.

Evaluation of the pentose phosphate pathway from 14CO2 data. Fallibility of a classic equation when applied to non-homogeneous tissues.

PubMed

Larrabee, M G

1990-11-15

A classic equation that has frequently been used to estimate the fraction of glucose metabolized by the pentose phosphate pathway, using 14CO2 data, is more simply re-derived with careful consideration of the assumptions involved and the conditions under which it is applicable. The equation is shown to be unreliable for non-homogeneous tissues, depending on the fraction of triose phosphate converted to CO2. The formula in question is as follows: ([1]CO2/G-[6]CO2/G)/(1-[6]CO2/G) = 3Fmet./(1 + 2Fmet.) where [1]CO2 and [6]CO2 are output rates of carbons 1 and 6 of glucose respectively to CO2, G is the rate of glucose uptake and Fmet. is the fraction of the glucose that is metabolized to CO2 and triose phosphate by the pentose phosphate pathway, allowing for recycling of an appropriate fraction of the fructose-6-phosphate produced by the pathway. This analysis illustrates the importance of suitably testing any equation that assumes homogeneity before application to non-homogeneous tissues.

A semi-analytical method for the computation of the Lyapunov exponents of fractional-order systems

NASA Astrophysics Data System (ADS)

Caponetto, Riccardo; Fazzino, Stefano

2013-01-01

Fractional-order differential equations are interesting for their applications in the construction of mathematical models in finance, materials science or diffusion. In this paper, an application of a well known transformation technique, Differential Transform Method (DTM), to the area of fractional differential equation is employed for calculating Lyapunov exponents of fractional order systems. It is known that the Lyapunov exponents, first introduced by Oseledec, play a crucial role in characterizing the behaviour of dynamical systems. They can be used to analyze the sensitive dependence on initial conditions and the presence of chaotic attractors. The results reveal that the proposed method is very effective and simple and leads to accurate, approximately convergent solutions.

Abundant closed form solutions of the conformable time fractional Sawada-Kotera-Ito equation using (G‧ / G) -expansion method

NASA Astrophysics Data System (ADS)

Al-Shawba, Altaf Abdulkarem; Gepreel, K. A.; Abdullah, F. A.; Azmi, A.

2018-06-01

In current study, we use the (G‧ / G) -expansion method to construct the closed form solutions of the seventh order time fractional Sawada-Kotera-Ito (TFSKI) equation based on conformable fractional derivative. As a result, trigonometric, hyperbolic and rational functions solutions with arbitrary constants are obtained. When the arbitrary constants are taken some special values, the periodic and soliton solutions are obtained from the travelling wave solutions. The obtained solutions are new and not found elsewhere. The effect of the fractional order on some of these solutions are represented graphically to illustrate the behavior of the exact solutions when the parameter take some special choose.

Analytical solution to the fractional polytropic gas spheres

NASA Astrophysics Data System (ADS)

Nouh, Mohamed I.; Abdel-Salam, Emad A.-B.

2018-04-01

The Lane-Emden equation can be used to model stellar interiors, star clusters and many configurations in astrophysics. Unfortunately, there is an exact solution only for the polytropic indices n = 0, 1 and 5. In the present paper, a series solution for the fractional Lane-Emden equation is presented. The solution is performed in the frame of modified Rienmann Liouville derivatives. The obtained results recover the well-known series solutions when α =1. The fractional model of n = 3 is calculated and the mass-radius relation, density ratio, pressure ratio and temperature ratio are investigated. The fractional star appears much different than the integer star, as it is denser, more stressed and hotter than the integer star.

Dynamics of the magnetization of single domain particles having triaxial anisotropy subjected to a uniform dc magnetic field

NASA Astrophysics Data System (ADS)

Ouari, Bachir; Kalmykov, Yury P.

2006-12-01

Thermally induced relaxation of the magnetization of single domain ferromagnetic particles with triaxial (orthorhombic) anisotropy in the presence of a uniform external magnetic field H0 is considered in the context of Brown's continuous diffusion model. Simple analytic equations, which allow one to describe qualitatively the field effects in the relaxation behavior of the system for wide ranges of the field strength and damping parameters are derived. It is shown that these formulas are in complete agreement with the exact matrix continued fraction solution of the infinite hierarchy of linear differential-recurrence equations for the statistical moments, which governs the magnetization dynamics of an individual particle (this hierarchy is derived by averaging the underlying stochastic Landau-Lifshitz-Gilbert equation over its realizations). It is also demonstrated that in strong fields the longitudinal relaxation of the magnetization is essentially modified by the contribution of the high-frequency "intrawell" modes to the relaxation process. This effect discovered for uniaxial particles by Coffey et al. [Phys. Rev. B 51, 15947 (1995)] is the natural consequence of the depletion of population of the shallow potential well. However, in contrast to uniaxial anisotropy, for orthorhombic crystals there is an inherent geometric dependence of the complex magnetic susceptibility and the relaxation time on the damping parameter α arising from the coupling of longitudinal and transverse relaxation modes.

Discovering variable fractional orders of advection-dispersion equations from field data using multi-fidelity Bayesian optimization

NASA Astrophysics Data System (ADS)

Pang, Guofei; Perdikaris, Paris; Cai, Wei; Karniadakis, George Em

2017-11-01

The fractional advection-dispersion equation (FADE) can describe accurately the solute transport in groundwater but its fractional order has to be determined a priori. Here, we employ multi-fidelity Bayesian optimization to obtain the fractional order under various conditions, and we obtain more accurate results compared to previously published data. Moreover, the present method is very efficient as we use different levels of resolution to construct a stochastic surrogate model and quantify its uncertainty. We consider two different problem set ups. In the first set up, we obtain variable fractional orders of one-dimensional FADE, considering both synthetic and field data. In the second set up, we identify constant fractional orders of two-dimensional FADE using synthetic data. We employ multi-resolution simulations using two-level and three-level Gaussian process regression models to construct the surrogates.

A parallel algorithm for the two-dimensional time fractional diffusion equation with implicit difference method.

PubMed

Gong, Chunye; Bao, Weimin; Tang, Guojian; Jiang, Yuewen; Liu, Jie

2014-01-01

It is very time consuming to solve fractional differential equations. The computational complexity of two-dimensional fractional differential equation (2D-TFDE) with iterative implicit finite difference method is O(M(x)M(y)N(2)). In this paper, we present a parallel algorithm for 2D-TFDE and give an in-depth discussion about this algorithm. A task distribution model and data layout with virtual boundary are designed for this parallel algorithm. The experimental results show that the parallel algorithm compares well with the exact solution. The parallel algorithm on single Intel Xeon X5540 CPU runs 3.16-4.17 times faster than the serial algorithm on single CPU core. The parallel efficiency of 81 processes is up to 88.24% compared with 9 processes on a distributed memory cluster system. We do think that the parallel computing technology will become a very basic method for the computational intensive fractional applications in the near future.

Nonclassical point of view of the Brownian motion generation via fractional deterministic model

NASA Astrophysics Data System (ADS)

Gilardi-Velázquez, H. E.; Campos-Cantón, E.

In this paper, we present a dynamical system based on the Langevin equation without stochastic term and using fractional derivatives that exhibit properties of Brownian motion, i.e. a deterministic model to generate Brownian motion is proposed. The stochastic process is replaced by considering an additional degree of freedom in the second-order Langevin equation. Thus, it is transformed into a system of three first-order linear differential equations, additionally α-fractional derivative are considered which allow us to obtain better statistical properties. Switching surfaces are established as a part of fluctuating acceleration. The final system of three α-order linear differential equations does not contain a stochastic term, so the system generates motion in a deterministic way. Nevertheless, from the time series analysis, we found that the behavior of the system exhibits statistics properties of Brownian motion, such as, a linear growth in time of mean square displacement, a Gaussian distribution. Furthermore, we use the detrended fluctuation analysis to prove the Brownian character of this motion.

Maximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation

NASA Astrophysics Data System (ADS)

Yun, Dongfang; Protas, Bartosz

2018-02-01

This investigation is a part of a research program aiming to characterize the extreme behavior possible in hydrodynamic models by analyzing the maximum growth of certain fundamental quantities. We consider here the rate of growth of the classical and fractional enstrophy in the fractional Burgers equation in the subcritical and supercritical regimes. Since solutions to this equation exhibit, respectively, globally well-posed behavior and finite-time blowup in these two regimes, this makes it a useful model to study the maximum instantaneous growth of enstrophy possible in these two distinct situations. First, we obtain estimates on the rates of growth and then show that these estimates are sharp up to numerical prefactors. This is done by numerically solving suitably defined constrained maximization problems and then demonstrating that for different values of the fractional dissipation exponent the obtained maximizers saturate the upper bounds in the estimates as the enstrophy increases. We conclude that the power-law dependence of the enstrophy rate of growth on the fractional dissipation exponent has the same global form in the subcritical, critical and parts of the supercritical regime. This indicates that the maximum enstrophy rate of growth changes smoothly as global well-posedness is lost when the fractional dissipation exponent attains supercritical values. In addition, nontrivial behavior is revealed for the maximum rate of growth of the fractional enstrophy obtained for small values of the fractional dissipation exponents. We also characterize the structure of the maximizers in different cases.

Fractal ladder models and power law wave equations

PubMed Central

Kelly, James F.; McGough, Robert J.

2009-01-01

The ultrasonic attenuation coefficient in mammalian tissue is approximated by a frequency-dependent power law for frequencies less than 100 MHz. To describe this power law behavior in soft tissue, a hierarchical fractal network model is proposed. The viscoelastic and self-similar properties of tissue are captured by a constitutive equation based on a lumped parameter infinite-ladder topology involving alternating springs and dashpots. In the low-frequency limit, this ladder network yields a stress-strain constitutive equation with a time-fractional derivative. By combining this constitutive equation with linearized conservation principles and an adiabatic equation of state, a fractional partial differential equation that describes power law attenuation is derived. The resulting attenuation coefficient is a power law with exponent ranging between 1 and 2, while the phase velocity is in agreement with the Kramers–Kronig relations. The fractal ladder model is compared to published attenuation coefficient data, thus providing equivalent lumped parameters. PMID:19813816

BRYNTRN: A baryon transport model

NASA Technical Reports Server (NTRS)

Wilson, John W.; Townsend, Lawrence W.; Nealy, John E.; Chun, Sang Y.; Hong, B. S.; Buck, Warren W.; Lamkin, S. L.; Ganapol, Barry D.; Khan, Ferdous; Cucinotta, Francis A.

1989-01-01

The development of an interaction data base and a numerical solution to the transport of baryons through an arbitrary shield material based on a straight ahead approximation of the Boltzmann equation are described. The code is most accurate for continuous energy boundary values, but gives reasonable results for discrete spectra at the boundary using even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O). The resulting computer code is self-contained, efficient and ready to use. The code requires only a very small fraction of the computer resources required for Monte Carlo codes.

A new approach to exact optical soliton solutions for the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Morales-Delgado, V. F.; Gómez-Aguilar, J. F.; Baleanu, Dumitru

2018-05-01

By using the modified homotopy analysis transform method, we construct the analytical solutions of the space-time generalized nonlinear Schrödinger equation involving a new fractional conformable derivative in the Liouville-Caputo sense and the fractional-order derivative with the Mittag-Leffler law. Employing theoretical parameters, we present some numerical simulations and compare the solutions obtained.

Lie Symmetry Analysis, Conservation Laws and Exact Power Series Solutions for Time-Fractional Fordy-Gibbons Equation

NASA Astrophysics Data System (ADS)

Feng, Lian-Li; Tian, Shou-Fu; Wang, Xiu-Bin; Zhang, Tian-Tian

2016-09-01

In this paper, the time fractional Fordy-Gibbons equation is investigated with Riemann-Liouville derivative. The equation can be reduced to the Caudrey-Dodd-Gibbon equation, Savada-Kotera equation and the Kaup-Kupershmidt equation, etc. By means of the Lie group analysis method, the invariance properties and symmetry reductions of the equation are derived. Furthermore, by means of the power series theory, its exact power series solutions of the equation are also constructed. Finally, two kinds of conservation laws of the equation are well obtained with aid of the self-adjoint method. Supported by the Fundamental Research Funds for Key Discipline Construction under Grant No. XZD201602, the Fundamental Research Funds for the Central Universities under Grant Nos. 2015QNA53 and 2015XKQY14, the Fundamental Research Funds for Postdoctoral at the Key Laboratory of Gas and Fire Control for Coal Mines, the General Financial Grant from the China Postdoctoral Science Foundation under Grant No. 2015M570498, and Natural Sciences Foundation of China under Grant No. 11301527

Computational Algorithms or Identification of Distributed Parameter Systems

DTIC Science & Technology

1993-04-24

delay-differential equations, Volterra integral equations, and partial differential equations with memory terms . In particular we investigated a...tested for estimating parameters in a Volterra integral equation arising from a viscoelastic model of a flexible structure with Boltzmann damping. In...particular, one of the parameters identified was the order of the derivative in Volterra integro-differential equations containing fractional

Pharmacokinetics and Pharmacodynamics of Piperacillin-Tazobactam in 42 Patients Treated with Concomitant CRRT

PubMed Central

Bauer, Seth R.; Salem, Charbel; Connor, Michael J.; Groszek, Joseph; Taylor, Maria E.; Wei, Peilin; Tolwani, Ashita J.

2012-01-01

Summary Background and objectives Current recommendations for piperacillin-tazobactam dosing in patients receiving continuous renal replacement therapy originate from studies with relatively few patients and lower continuous renal replacement therapy doses than commonly used today. This study measured the pharmacokinetic and pharmacodynamic characteristics of piperacillin-tazobactam in patients treated with continuous renal replacement therapy using contemporary equipment and prescriptions. Design, setting, participants, & measurements A multicenter prospective observational study in the intensive care units of two academic medical centers was performed, enrolling patients with AKI or ESRD receiving piperacillin-tazobactam while being treated with continuous renal replacement therapy. Pregnant women, children, and patients with end stage liver disease were excluded from enrollment. Plasma and continuous renal replacement therapy effluent samples were analyzed for piperacillin and tazobactam levels using HPLC. Pharmacokinetic and pharmacodynamic parameters were calculated using standard equations. Multivariate analyses were used to examine the association of patient and continuous renal replacement therapy characteristics with piperacillin pharmacokinetic parameters. Results Forty-two of fifty-five subjects enrolled had complete sampling. Volume of distribution (median=0.38 L/kg, intraquartile range=0.20 L/kg) and elimination rate constants (median=0.104 h−1, intraquartile range=0.052 h−1) were highly variable, and clinical parameters could explain only a small fraction of the large variability in pharmacokinetic parameters. Probability of target attainment for piperacillin was 83% for total drug but only 77% when the unbound fraction was considered. Conclusions There is significant patient to patient variability in pharmacokinetic/pharmacodynamic parameters in patients receiving continuous renal replacement therapy. Many patients did not achieve pharmacodynamic targets, suggesting that therapeutic drug monitoring might optimize therapy. PMID:22282479

Elements of an improved model of debris‐flow motion

USGS Publications Warehouse

Iverson, Richard M.

2009-01-01

A new depth‐averaged model of debris‐flow motion describes simultaneous evolution of flow velocity and depth, solid and fluid volume fractions, and pore‐fluid pressure. Non‐hydrostatic pore‐fluid pressure is produced by dilatancy, a state‐dependent property that links the depth‐averaged shear rate and volumetric strain rate of the granular phase. Pore‐pressure changes caused by shearing allow the model to exhibit rate‐dependent flow resistance, despite the fact that the basal shear traction involves only rate‐independent Coulomb friction. An analytical solution of simplified model equations shows that the onset of downslope motion can be accelerated or retarded by pore‐pressure change, contingent on whether dilatancy is positive or negative. A different analytical solution shows that such effects will likely be muted if downslope motion continues long enough, because dilatancy then evolves toward zero, and volume fractions and pore pressure concurrently evolve toward steady states.

Elements of an improved model of debris-flow motion

USGS Publications Warehouse

Iverson, R.M.

2009-01-01

A new depth-averaged model of debris-flow motion describes simultaneous evolution of flow velocity and depth, solid and fluid volume fractions, and pore-fluid pressure. Non-hydrostatic pore-fluid pressure is produced by dilatancy, a state-dependent property that links the depth-averaged shear rate and volumetric strain rate of the granular phase. Pore-pressure changes caused by shearing allow the model to exhibit rate-dependent flow resistance, despite the fact that the basal shear traction involves only rate-independent Coulomb friction. An analytical solution of simplified model equations shows that the onset of downslope motion can be accelerated or retarded by pore-pressure change, contingent on whether dilatancy is positive or negative. A different analytical solution shows that such effects will likely be muted if downslope motion continues long enough, because dilatancy then evolves toward zero, and volume fractions and pore pressure concurrently evolve toward steady states. ?? 2009 American Institute of Physics.

A comparative mathematical analysis of RL and RC electrical circuits via Atangana-Baleanu and Caputo-Fabrizio fractional derivatives

NASA Astrophysics Data System (ADS)

Abro, Kashif Ali; Memon, Anwar Ahmed; Uqaili, Muhammad Aslam

2018-03-01

This research article is analyzed for the comparative study of RL and RC electrical circuits by employing newly presented Atangana-Baleanu and Caputo-Fabrizio fractional derivatives. The governing ordinary differential equations of RL and RC electrical circuits have been fractionalized in terms of fractional operators in the range of 0 ≤ ξ ≤ 1 and 0 ≤ η ≤ 1. The analytic solutions of fractional differential equations for RL and RC electrical circuits have been solved by using the Laplace transform with its inversions. General solutions have been investigated for periodic and exponential sources by implementing the Atangana-Baleanu and Caputo-Fabrizio fractional operators separately. The investigated solutions have been expressed in terms of simple elementary functions with convolution product. On the basis of newly fractional derivatives with and without singular kernel, the voltage and current have interesting behavior with several similarities and differences for the periodic and exponential sources.

New Insights into the Fractional Order Diffusion Equation Using Entropy and Kurtosis.

PubMed

Ingo, Carson; Magin, Richard L; Parrish, Todd B

2014-11-01

Fractional order derivative operators offer a concise description to model multi-scale, heterogeneous and non-local systems. Specifically, in magnetic resonance imaging, there has been recent work to apply fractional order derivatives to model the non-Gaussian diffusion signal, which is ubiquitous in the movement of water protons within biological tissue. To provide a new perspective for establishing the utility of fractional order models, we apply entropy for the case of anomalous diffusion governed by a fractional order diffusion equation generalized in space and in time. This fractional order representation, in the form of the Mittag-Leffler function, gives an entropy minimum for the integer case of Gaussian diffusion and greater values of spectral entropy for non-integer values of the space and time derivatives. Furthermore, we consider kurtosis, defined as the normalized fourth moment, as another probabilistic description of the fractional time derivative. Finally, we demonstrate the implementation of anomalous diffusion, entropy and kurtosis measurements in diffusion weighted magnetic resonance imaging in the brain of a chronic ischemic stroke patient.

Partitioning of photosynthate within and distal to the growth zone of tall fescue leaf blades

DOE Office of Scientific and Technical Information (OSTI.GOV)

Allard, G.; Nelson, C.J.

The growth zone of developing tall fescue leaf blades, consisting of zones of cell division, cell elongation and cell maturation, are a strong sink for photosynthate. Distribution of {sup 14}C along the growth zone and partitioning between water-soluble carbohydrates (WSC) and WSC-free dry matter (SDM) were observed for up to 64 h after labelling the youngest fully expanded leaf or the exposed tip of the elongating leaf with {sup 14}CO{sub 2}. Deposition rates of {sup 14}C were estimated using the continuity equation. The cell elongation zone, especially the WSC fraction, was the strongest sink for {sup 14}C. In the proximalmore » end of the maturation zone, partitioning of {sup 14}C shifted from being allocated mainly to the WSC fraction after 2 h to an equal distribution between WSC and SDM at the distal end of the zone after 64 h. A significant proportion of {sup 14}C in the SDM fraction of the maturation zone could be attributed to redistribution from WSC.« less

The random continued fraction transformation

NASA Astrophysics Data System (ADS)

Kalle, Charlene; Kempton, Tom; Verbitskiy, Evgeny

2017-03-01

We introduce a random dynamical system related to continued fraction expansions. It uses random combinations of the Gauss map and the Rényi (or backwards) continued fraction map. We explore the continued fraction expansions that this system produces, as well as the dynamical properties of the system.

The second hyperpolarizability of systems described by the space-fractional Schrödinger equation

NASA Astrophysics Data System (ADS)

Dawson, Nathan J.; Nottage, Onassis; Kounta, Moussa

2018-01-01

The static second hyperpolarizability is derived from the space-fractional Schrödinger equation in the particle-centric view. The Thomas-Reiche-Kuhn sum rule matrix elements and the three-level ansatz determines the maximum second hyperpolarizability for a space-fractional quantum system. The total oscillator strength is shown to decrease as the space-fractional parameter α decreases, which reduces the optical response of a quantum system in the presence of an external field. This damped response is caused by the wavefunction dependent position and momentum commutation relation. Although the maximum response is damped, we show that the one-dimensional quantum harmonic oscillator is no longer a linear system for α ≠ 1, where the second hyperpolarizability becomes negative before ultimately damping to zero at the lower fractional limit of α → 1 / 2.

Analysis of the cable equation with non-local and non-singular kernel fractional derivative

NASA Astrophysics Data System (ADS)

Karaagac, Berat

2018-02-01

Recently a new concept of differentiation was introduced in the literature where the kernel was converted from non-local singular to non-local and non-singular. One of the great advantages of this new kernel is its ability to portray fading memory and also well defined memory of the system under investigation. In this paper the cable equation which is used to develop mathematical models of signal decay in submarine or underwater telegraphic cables will be analysed using the Atangana-Baleanu fractional derivative due to the ability of the new fractional derivative to describe non-local fading memory. The existence and uniqueness of the more generalized model is presented in detail via the fixed point theorem. A new numerical scheme is used to solve the new equation. In addition, stability, convergence and numerical simulations are presented.

On mixed derivatives type high dimensional multi-term fractional partial differential equations approximate solutions

NASA Astrophysics Data System (ADS)

Talib, Imran; Belgacem, Fethi Bin Muhammad; Asif, Naseer Ahmad; Khalil, Hammad

2017-01-01

In this research article, we derive and analyze an efficient spectral method based on the operational matrices of three dimensional orthogonal Jacobi polynomials to solve numerically the mixed partial derivatives type multi-terms high dimensions generalized class of fractional order partial differential equations. We transform the considered fractional order problem to an easily solvable algebraic equations with the aid of the operational matrices. Being easily solvable, the associated algebraic system leads to finding the solution of the problem. Some test problems are considered to confirm the accuracy and validity of the proposed numerical method. The convergence of the method is ensured by comparing our Matlab software simulations based obtained results with the exact solutions in the literature, yielding negligible errors. Moreover, comparative results discussed in the literature are extended and improved in this study.

Optimal control problem for linear fractional-order systems, described by equations with Hadamard-type derivative

NASA Astrophysics Data System (ADS)

Postnov, Sergey

2017-11-01

Two kinds of optimal control problem are investigated for linear time-invariant fractional-order systems with lumped parameters which dynamics described by equations with Hadamard-type derivative: the problem of control with minimal norm and the problem of control with minimal time at given restriction on control norm. The problem setting with nonlocal initial conditions studied. Admissible controls allowed to be the p-integrable functions (p > 1) at half-interval. The optimal control problem studied by moment method. The correctness and solvability conditions for the corresponding moment problem are derived. For several special cases the optimal control problems stated are solved analytically. Some analogies pointed for results obtained with the results which are known for integer-order systems and fractional-order systems describing by equations with Caputo- and Riemann-Liouville-type derivatives.

The first boundary-value problem for a fractional diffusion-wave equation in a non-cylindrical domain

NASA Astrophysics Data System (ADS)

Pskhu, A. V.

2017-12-01

We solve the first boundary-value problem in a non-cylindrical domain for a diffusion-wave equation with the Dzhrbashyan- Nersesyan operator of fractional differentiation with respect to the time variable. We prove an existence and uniqueness theorem for this problem, and construct a representation of the solution. We show that a sufficient condition for unique solubility is the condition of Hölder smoothness for the lateral boundary of the domain. The corresponding results for equations with Riemann- Liouville and Caputo derivatives are particular cases of results obtained here.

Sediment composition for the assessment of water erosion and nonpoint source pollution in natural and fire-affected landscapes.

PubMed

Carkovic, Athena B; Pastén, Pablo A; Bonilla, Carlos A

2015-04-15

Water erosion is a leading cause of soil degradation and a major nonpoint source pollution problem. Many efforts have been undertaken to estimate the amount and size distribution of the sediment leaving the field. Multi-size class water erosion models subdivide eroded soil into different sizes and estimate the aggregate's composition based on empirical equations derived from agricultural soils. The objective of this study was to evaluate these equations on soil samples collected from natural landscapes (uncultivated) and fire-affected soils. Chemical, physical, and soil fractions and aggregate composition analyses were performed on samples collected in the Chilean Patagonia and later compared with the equations' estimates. The results showed that the empirical equations were not suitable for predicting the sediment fractions. Fine particles, including primary clay, primary silt, and small aggregates (<53 μm) were over-estimated, and large aggregates (>53 μm) and primary sand were under-estimated. The uncultivated and fire-affected soils showed a reduced fraction of fine particles in the sediment, as clay and silt were mostly in the form of large aggregates. Thus, a new set of equations was developed for these soils, where small aggregates were defined as particles with sizes between 53 μm and 250 μm and large aggregates as particles>250 μm. With r(2) values between 0.47 and 0.98, the new equations provided better estimates for primary sand and large aggregates. The aggregate's composition was also well predicted, especially the silt and clay fractions in the large aggregates from uncultivated soils (r(2)=0.63 and 0.83, respectively) and the fractions of silt in the small aggregates (r(2)=0.84) and clay in the large aggregates (r(2)=0.78) from fire-affected soils. Overall, these new equations proved to be better predictors for the sediment and aggregate's composition in uncultivated and fire-affected soils, and they reduce the error when estimating soil loss in natural landscapes. Copyright © 2015 Elsevier B.V. All rights reserved.

Analytical and numerical treatment of the heat conduction equation obtained via time-fractional distributed-order heat conduction law

NASA Astrophysics Data System (ADS)

Želi, Velibor; Zorica, Dušan

2018-02-01

Generalization of the heat conduction equation is obtained by considering the system of equations consisting of the energy balance equation and fractional-order constitutive heat conduction law, assumed in the form of the distributed-order Cattaneo type. The Cauchy problem for system of energy balance equation and constitutive heat conduction law is treated analytically through Fourier and Laplace integral transform methods, as well as numerically by the method of finite differences through Adams-Bashforth and Grünwald-Letnikov schemes for approximation derivatives in temporal domain and leap frog scheme for spatial derivatives. Numerical examples, showing time evolution of temperature and heat flux spatial profiles, demonstrate applicability and good agreement of both methods in cases of multi-term and power-type distributed-order heat conduction laws.

Fractional populations in sex-linked inheritance

NASA Astrophysics Data System (ADS)

Pyo Lee, Seung; Chung, Myung-Hoon; Koo Kim, Chul; Nahm, Kyun

2001-03-01

We study the fractional populations in chromosome inherited diseases. The governing equations for the fractional populations are found and solved in the presence of mutation and selection. The physical fixed points obtained are used to discuss the cases of color blindness and hemophilia.

Smoothed particle hydrodynamics method for evaporating multiphase flows.

PubMed

Yang, Xiufeng; Kong, Song-Charng

2017-09-01

The smoothed particle hydrodynamics (SPH) method has been increasingly used for simulating fluid flows; however, its ability to simulate evaporating flow requires significant improvements. This paper proposes an SPH method for evaporating multiphase flows. The present SPH method can simulate the heat and mass transfers across the liquid-gas interfaces. The conservation equations of mass, momentum, and energy were reformulated based on SPH, then were used to govern the fluid flow and heat transfer in both the liquid and gas phases. The continuity equation of the vapor species was employed to simulate the vapor mass fraction in the gas phase. The vapor mass fraction at the interface was predicted by the Clausius-Clapeyron correlation. An evaporation rate was derived to predict the mass transfer from the liquid phase to the gas phase at the interface. Because of the mass transfer across the liquid-gas interface, the mass of an SPH particle was allowed to change. Alternative particle splitting and merging techniques were developed to avoid large mass difference between SPH particles of the same phase. The proposed method was tested by simulating three problems, including the Stefan problem, evaporation of a static drop, and evaporation of a drop impacting a hot surface. For the Stefan problem, the SPH results of the evaporation rate at the interface agreed well with the analytical solution. For drop evaporation, the SPH result was compared with the result predicted by a level-set method from the literature. In the case of drop impact on a hot surface, the evolution of the shape of the drop, temperature, and vapor mass fraction were predicted.

40 CFR 98.213 - Calculating GHG emissions.

Code of Federal Regulations, 2010 CFR

2010-07-01

... calcination fractions with Equation U-1 of this section. ER30OC09.077 Where: ECO2 = Annual CO2 mass emissions... ton carbonate consumed. Fi = Fraction calcination achieved for each particular carbonate type i (decimal fraction). As an alternative to measuring the calcination fraction, a value of 1.0 can be used. n...

A numerical framework for bubble transport in a subcooled fluid flow

NASA Astrophysics Data System (ADS)

Jareteg, Klas; Sasic, Srdjan; Vinai, Paolo; Demazière, Christophe

2017-09-01

In this paper we present a framework for the simulation of dispersed bubbly two-phase flows, with the specific aim of describing vapor-liquid systems with condensation. We formulate and implement a framework that consists of a population balance equation (PBE) for the bubble size distribution and an Eulerian-Eulerian two-fluid solver. The PBE is discretized using the Direct Quadrature Method of Moments (DQMOM) in which we include the condensation of the bubbles as an internal phase space convection. We investigate the robustness of the DQMOM formulation and the numerical issues arising from the rapid shrinkage of the vapor bubbles. In contrast to a PBE method based on the multiple-size-group (MUSIG) method, the DQMOM formulation allows us to compute a distribution with dynamic bubble sizes. Such a property is advantageous to capture the wide range of bubble sizes associated with the condensation process. Furthermore, we compare the computational performance of the DQMOM-based framework with the MUSIG method. The results demonstrate that DQMOM is able to retrieve the bubble size distribution with a good numerical precision in only a small fraction of the computational time required by MUSIG. For the two-fluid solver, we examine the implementation of the mass, momentum and enthalpy conservation equations in relation to the coupling to the PBE. In particular, we propose a formulation of the pressure and liquid continuity equations, that was shown to correctly preserve mass when computing the vapor fraction with DQMOM. In addition, the conservation of enthalpy was also proven. Therefore a consistent overall framework that couples the PBE and two-fluid solvers is achieved.

From Fractals to Fractional Vector Calculus: Measurement in the Correct Metric

NASA Astrophysics Data System (ADS)

Wheatcraft, S. W.; Meerschaert, M. M.; Mortensen, J.

2005-12-01

Traditional (stationary) stochastic theories have been fairly successful in reproducing transport behavior at relatively homogeneous field sites such as the Borden and Cape Code sites. However, the highly heterogeneous MADE site has produced tracer data that can not be adequately explained with traditional stochastic theories. In recent years, considerable attention has been focused on developing more sophisticated theories that can predict or reproduce the behavior of complex sites such as the MADE site. People began to realize that the model for geologic complexity may in many cases be very different than the model required for stochastic theory. Fractal approaches were useful in conceptualizing scale-invariant heterogeneity by demonstrating that scale dependant transport was just an artifact of our measurement system. Fractal media have dimensions larger than the dimension that measurement is taking place in, thus assuring the scale-dependence of parameters such as dispersivity. What was needed was a rigorous way to develop a theory that was consistent with the fractal dimension of the heterogeneity. The fractional advection-dispersion equation (FADE) was developed with this idea in mind. The second derivative in the dispersion term of the advection-dispersion equation is replaced with a fractional derivative. The order of differentiation, α, is fractional. Values of α in the range: 1 < α < 2 produce super-Fickian dispersion; in essence, the dispersion scaling is controlled by the value of α. When α = 2, the traditional advection-dispersion equation is recovered. The 1-D version of the FADE has been used successfully to back-predict tracer test behavior at several heterogeneous field sites, including the MADE site. It has been hypothesized that the order of differentiation in the FADE is equivalent to (or at least related to) the fractal dimension of the particle tracks (or geologic heterogeneity). With this way of thinking, one can think of the FADE as a governing equation written for the correct dimension, thus eliminating scale-dependent behavior. Before a generalized multi-dimensional form of the FADE can be developed, it has been necessary to develop a generalized fractional vector calculus. The authors have recently developed generalized canonical fractional forms of the gradient, divergence and curl. This fractional vector calculus will be useful in developing fractional forms of many governing equations in physics.

Variable order fractional Fokker-Planck equations derived from Continuous Time Random Walks

NASA Astrophysics Data System (ADS)

Straka, Peter

2018-08-01

Continuous Time Random Walk models (CTRW) of anomalous diffusion are studied, where the anomalous exponent β(x) ∈(0 , 1) varies in space. This type of situation occurs e.g. in biophysics, where the density of the intracellular matrix varies throughout a cell. Scaling limits of CTRWs are known to have probability distributions which solve fractional Fokker-Planck type equations (FFPE). This correspondence between stochastic processes and FFPE solutions has many useful extensions e.g. to nonlinear particle interactions and reactions, but has not yet been sufficiently developed for FFPEs of the "variable order" type with non-constant β(x) . In this article, variable order FFPEs (VOFFPE) are derived from scaling limits of CTRWs. The key mathematical tool is the 1-1 correspondence of a CTRW scaling limit to a bivariate Langevin process, which tracks the cumulative sum of jumps in one component and the cumulative sum of waiting times in the other. The spatially varying anomalous exponent is modelled by spatially varying β(x) -stable Lévy noise in the waiting time component. The VOFFPE displays a spatially heterogeneous temporal scaling behaviour, with generalized diffusivity and drift coefficients whose units are length2/timeβ(x) resp. length/timeβ(x). A global change of the time scale results in a spatially varying change in diffusivity and drift. A consequence of the mathematical derivation of a VOFFPE from CTRW limits in this article is that a solution of a VOFFPE can be approximated via Monte Carlo simulations. Based on such simulations, we are able to confirm that the VOFFPE is consistent under a change of the global time scale.

Troposphere mapping functions for GPS and very long baseline interferometry from European Centre for Medium-Range Weather Forecasts operational analysis data

NASA Astrophysics Data System (ADS)

Boehm, Johannes; Werl, Birgit; Schuh, Harald

2006-02-01

In the analyses of geodetic very long baseline interferometry (VLBI) and GPS data the analytic form used for mapping of the atmosphere delay from zenith to the line of site is most often a three-parameter continued fraction in 1/sin(elevation). Using the 40 years reanalysis (ERA-40) data of the European Centre for Medium-Range Weather Forecasts for the year 2001, the b and c coefficients of the continued fraction form for the hydrostatic mapping functions have been redetermined. Unlike previous mapping functions based on data from numerical weather models (isobaric mapping functions (Niell, 2000) and Vienna mapping functions (VMF) (Boehm and Schuh, 2004)), the new c coefficients are dependent on the day of the year, and unlike the Niell mapping functions (Niell, 1996) they are no longer symmetric with respect to the equator (apart from the opposite phase for the two hemispheres). Compared to VMF, this causes an effect on the VLBI or GPS station heights that is constant and as large as 2 mm at the equator and that varies seasonally between 4 mm and 0 mm at the poles. The updated VMF, based on these new coefficients and called VMF1 hereinafter, yields slightly better baseline length repeatabilities for VLBI data. The hydrostatic and wet mapping functions are applied in various combinations with different kinds of a priori zenith delays in the analyses of all VLBI International VLBI Service for Geodesy and Astrometry (IVS)-R1 and IVS-R4 24-hour sessions of 2002 and 2003; the investigations concentrate on baseline length repeatabilities, as well as on absolute changes of station heights.

Plasma protein osmotic pressure equations and nomogram for sheep.

PubMed

Yamada, S; Grady, M K; Licko, V; Staub, N C

1991-08-01

The equations developed by Landis and Pappenheimer (Handbook of Physiology. Circulation, 1963, p. 961-1034) for calculating the protein osmotic pressure of human plasma proteins have been frequently used for other animal species without regard to the fractional albumin concentration or correction for protein-protein interaction. Using an electronic osmometer, we remeasured the protein osmotic pressure of purified sheep albumin and sheep plasma partially depleted of albumin. We measured protein osmotic pressures of serial dilutions over the concentration range 0-180 g/l for albumin and 0-100 g/l for the albumin-depleted proteins at room temperature (26 degrees C). Using a nonlinear least squares parameter-fitting computer program, we obtained the equation of best fit for purified albumin, and then we used that equation together with the measured albumin fraction to obtain the best-fit equation for the nonalbumin proteins. The equation for albumin is IIcmH2O,39 degrees C = 0.382C + 0.0028C2 + 0.000013C3, where C is albumin concentration in g/l. The equation for the nonalbumin fraction is IIcmH2O,39 degrees C = 0.119C + 0.0016C2. Up to 200- and 100-g/l protein concentration, respectively, these equations give the least standard error of the estimate for each of the virial coefficients. The computed number-average molecular weight for the nonalbumin proteins is 222,000. Using the new equations, we constructed a nomogram, based on the one of Nitta and co-workers (Tohoku J. Exp. Med. 135: 43-49, 1981). We tested the nomogram using 144 random samples of sheep plasma and lymph from 31 sheep. We obtained a correlation coefficient of 0.99 between the measured and nomogram estimates of protein osmotic pressure.

Multiplicative noise removal through fractional order tv-based model and fast numerical schemes for its approximation

NASA Astrophysics Data System (ADS)

Ullah, Asmat; Chen, Wen; Khan, Mushtaq Ahmad

2017-07-01

This paper introduces a fractional order total variation (FOTV) based model with three different weights in the fractional order derivative definition for multiplicative noise removal purpose. The fractional-order Euler Lagrange equation which is a highly non-linear partial differential equation (PDE) is obtained by the minimization of the energy functional for image restoration. Two numerical schemes namely an iterative scheme based on the dual theory and majorization- minimization algorithm (MMA) are used. To improve the restoration results, we opt for an adaptive parameter selection procedure for the proposed model by applying the trial and error method. We report numerical simulations which show the validity and state of the art performance of the fractional-order model in visual improvement as well as an increase in the peak signal to noise ratio comparing to corresponding methods. Numerical experiments also demonstrate that MMAbased methodology is slightly better than that of an iterative scheme.

A Kronecker product splitting preconditioner for two-dimensional space-fractional diffusion equations

NASA Astrophysics Data System (ADS)

Chen, Hao; Lv, Wen; Zhang, Tongtong

2018-05-01

We study preconditioned iterative methods for the linear system arising in the numerical discretization of a two-dimensional space-fractional diffusion equation. Our approach is based on a formulation of the discrete problem that is shown to be the sum of two Kronecker products. By making use of an alternating Kronecker product splitting iteration technique we establish a class of fixed-point iteration methods. Theoretical analysis shows that the new method converges to the unique solution of the linear system. Moreover, the optimal choice of the involved iteration parameters and the corresponding asymptotic convergence rate are computed exactly when the eigenvalues of the system matrix are all real. The basic iteration is accelerated by a Krylov subspace method like GMRES. The corresponding preconditioner is in a form of a Kronecker product structure and requires at each iteration the solution of a set of discrete one-dimensional fractional diffusion equations. We use structure preserving approximations to the discrete one-dimensional fractional diffusion operators in the action of the preconditioning matrix. Numerical examples are presented to illustrate the effectiveness of this approach.

Representation of solution for fully nonlocal diffusion equations with deviation time variable

NASA Astrophysics Data System (ADS)

Drin, I. I.; Drin, S. S.; Drin, Ya. M.

2018-01-01

We prove the solvability of the Cauchy problem for a nonlocal heat equation which is of fractional order both in space and time. The representation formula for classical solutions for time- and space- fractional partial differential operator Dat + a2 (-Δ) γ/2 (0 <= α <= 1, γ ɛ (0, 2]) and deviation time variable is given in terms of the Fox H-function, using the step by step method.

Existence and uniqueness theorems for impulsive fractional differential equations with the two-point and integral boundary conditions.

PubMed

Mardanov, M J; Mahmudov, N I; Sharifov, Y A

2014-01-01

We study a boundary value problem for the system of nonlinear impulsive fractional differential equations of order α (0 < α ≤ 1) involving the two-point and integral boundary conditions. Some new results on existence and uniqueness of a solution are established by using fixed point theorems. Some illustrative examples are also presented. We extend previous results even in the integer case α = 1.

Transport equations for subdiffusion with nonlinear particle interaction.

PubMed

Straka, P; Fedotov, S

2015-02-07

We show how the nonlinear interaction effects 'volume filling' and 'adhesion' can be incorporated into the fractional subdiffusive transport of cells and individual organisms. To this end, we use microscopic random walk models with anomalous trapping and systematically derive generic non-Markovian and nonlinear governing equations for the mean concentrations of the subdiffusive cells or organisms. We uncover an interesting interaction between the nonlinearities and the non-Markovian nature of the transport. In the subdiffusive case, this interaction manifests itself in a nontrivial combination of nonlinear terms with fractional derivatives. In the long time limit, however, these equations simplify to a form without fractional operators. This provides an easy method for the study of aggregation phenomena. In particular, this enables us to show that volume filling can prevent "anomalous aggregation," which occurs in subdiffusive systems with a spatially varying anomalous exponent. Copyright © 2014 Elsevier Ltd. All rights reserved.

Modeling of heat conduction via fractional derivatives

NASA Astrophysics Data System (ADS)

Fabrizio, Mauro; Giorgi, Claudio; Morro, Angelo

2017-09-01

The modeling of heat conduction is considered by letting the time derivative, in the Cattaneo-Maxwell equation, be replaced by a derivative of fractional order. The purpose of this new approach is to overcome some drawbacks of the Cattaneo-Maxwell equation, for instance possible fluctuations which violate the non-negativity of the absolute temperature. Consistency with thermodynamics is shown to hold for a suitable free energy potential, that is in fact a functional of the summed history of the heat flux, subject to a suitable restriction on the set of admissible histories. Compatibility with wave propagation at a finite speed is investigated in connection with temperature-rate waves. It follows that though, as expected, this is the case for the Cattaneo-Maxwell equation, the model involving the fractional derivative does not allow the propagation at a finite speed. Nevertheless, this new model provides a good description of wave-like profiles in thermal propagation phenomena, whereas Fourier's law does not.

A Legendre tau-spectral method for solving time-fractional heat equation with nonlocal conditions.

PubMed

Bhrawy, A H; Alghamdi, M A

2014-01-01

We develop the tau-spectral method to solve the time-fractional heat equation (T-FHE) with nonlocal condition. In order to achieve highly accurate solution of this problem, the operational matrix of fractional integration (described in the Riemann-Liouville sense) for shifted Legendre polynomials is investigated in conjunction with tau-spectral scheme and the Legendre operational polynomials are used as the base function. The main advantage in using the presented scheme is that it converts the T-FHE with nonlocal condition to a system of algebraic equations that simplifies the problem. For demonstrating the validity and applicability of the developed spectral scheme, two numerical examples are presented. The logarithmic graphs of the maximum absolute errors is presented to achieve the exponential convergence of the proposed method. Comparing between our spectral method and other methods ensures that our method is more accurate than those solved similar problem.

A Legendre tau-Spectral Method for Solving Time-Fractional Heat Equation with Nonlocal Conditions

PubMed Central

Bhrawy, A. H.; Alghamdi, M. A.

2014-01-01

We develop the tau-spectral method to solve the time-fractional heat equation (T-FHE) with nonlocal condition. In order to achieve highly accurate solution of this problem, the operational matrix of fractional integration (described in the Riemann-Liouville sense) for shifted Legendre polynomials is investigated in conjunction with tau-spectral scheme and the Legendre operational polynomials are used as the base function. The main advantage in using the presented scheme is that it converts the T-FHE with nonlocal condition to a system of algebraic equations that simplifies the problem. For demonstrating the validity and applicability of the developed spectral scheme, two numerical examples are presented. The logarithmic graphs of the maximum absolute errors is presented to achieve the exponential convergence of the proposed method. Comparing between our spectral method and other methods ensures that our method is more accurate than those solved similar problem. PMID:25057507

Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

NASA Astrophysics Data System (ADS)

Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em

2017-12-01

Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

Numerical Analysis of an H 1-Galerkin Mixed Finite Element Method for Time Fractional Telegraph Equation

PubMed Central

Wang, Jinfeng; Zhao, Meng; Zhang, Min; Liu, Yang; Li, Hong

2014-01-01

We discuss and analyze an H 1-Galerkin mixed finite element (H 1-GMFE) method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H 1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H 1-GMFE method. Based on the discussion on the theoretical error analysis in L 2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H 1-norm. Moreover, we derive and analyze the stability of H 1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure. PMID:25184148

A meshless method for solving two-dimensional variable-order time fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Tayebi, A.; Shekari, Y.; Heydari, M. H.

2017-07-01

Several physical phenomena such as transformation of pollutants, energy, particles and many others can be described by the well-known convection-diffusion equation which is a combination of the diffusion and advection equations. In this paper, this equation is generalized with the concept of variable-order fractional derivatives. The generalized equation is called variable-order time fractional advection-diffusion equation (V-OTFA-DE). An accurate and robust meshless method based on the moving least squares (MLS) approximation and the finite difference scheme is proposed for its numerical solution on two-dimensional (2-D) arbitrary domains. In the time domain, the finite difference technique with a θ-weighted scheme and in the space domain, the MLS approximation are employed to obtain appropriate semi-discrete solutions. Since the newly developed method is a meshless approach, it does not require any background mesh structure to obtain semi-discrete solutions of the problem under consideration, and the numerical solutions are constructed entirely based on a set of scattered nodes. The proposed method is validated in solving three different examples including two benchmark problems and an applied problem of pollutant distribution in the atmosphere. In all such cases, the obtained results show that the proposed method is very accurate and robust. Moreover, a remarkable property so-called positive scheme for the proposed method is observed in solving concentration transport phenomena.

Lord Brouncker's Forgotten Sequence of Continued Fractions for Pi

ERIC Educational Resources Information Center

Osler, Thomas J.

2010-01-01

"Lord Brouncker's continued fraction for pi" is a well-known result. In this article, we show that Brouncker found not only this one continued fraction, but an entire infinite sequence of related continued fractions for pi. These were recorded in the "Arithmetica Infinitorum" by John Wallis, but appear to have been ignored and forgotten by modern…

Green function of the double-fractional Fokker-Planck equation: path integral and stochastic differential equations.

PubMed

Kleinert, H; Zatloukal, V

2013-11-01

The statistics of rare events, the so-called black-swan events, is governed by non-Gaussian distributions with heavy power-like tails. We calculate the Green functions of the associated Fokker-Planck equations and solve the related stochastic differential equations. We also discuss the subject in the framework of path integration.

Enriched reproducing kernel particle method for fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Ying, Yuping; Lian, Yanping; Tang, Shaoqiang; Liu, Wing Kam

2018-06-01

The reproducing kernel particle method (RKPM) has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advection-diffusion equation (FADE), which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM, leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.

Fractional power-law spatial dispersion in electrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tarasov, Vasily E., E-mail: tarasov@theory.sinp.msu.ru; Departamento de Análisis Matemático, Universidad de La Laguna, 38271 La Laguna, Tenerife; Trujillo, Juan J., E-mail: jtrujill@ullmat.es

2013-07-15

Electric fields in non-local media with power-law spatial dispersion are discussed. Equations involving a fractional Laplacian in the Riesz form that describe the electric fields in such non-local media are studied. The generalizations of Coulomb’s law and Debye’s screening for power-law non-local media are characterized. We consider simple models with anomalous behavior of plasma-like media with power-law spatial dispersions. The suggested fractional differential models for these plasma-like media are discussed to describe non-local properties of power-law type. -- Highlights: •Plasma-like non-local media with power-law spatial dispersion. •Fractional differential equations for electric fields in the media. •The generalizations of Coulomb’s lawmore » and Debye’s screening for the media.« less

Tissue responses to fractional transient heating with sinusoidal heat flux condition on skin surface.

PubMed

Ezzat, Magdy A; El-Bary, Alaa A; Al-Sowayan, Noorah S

2016-10-01

A fractional model of Bioheat equation for describing quantitatively the thermal responses of skin tissue under sinusoidal heat flux conditions on skin surface is given. Laplace transform technique is used to obtain the solution in a closed form. The resulting formulation is applied to one-dimensional application to investigate the temperature distribution in skin with instantaneous surface heating for different cases. According to the numerical results and its graphs, conclusion about the fractional bioheat transfer equation has been constructed. Sensitivity analysis is performed to explore the thermal effects of various control parameters on tissue temperature. The comparisons are made with the results obtained in the case of the absence of time-fractional order. © 2016 Japanese Society of Animal Science. © 2016 Japanese Society of Animal Science.

The effect of multi-directional nanocomposite materials on the vibrational response of thick shell panels with finite length and rested on two-parameter elastic foundations

NASA Astrophysics Data System (ADS)

Tahouneh, Vahid; Naei, Mohammad Hasan

2016-03-01

The main purpose of this paper is to investigate the effect of bidirectional continuously graded nanocomposite materials on free vibration of thick shell panels rested on elastic foundations. The elastic foundation is considered as a Pasternak model after adding a shear layer to the Winkler model. The panels reinforced by randomly oriented straight single-walled carbon nanotubes are considered. The volume fractions of SWCNTs are assumed to be graded not only in the radial direction, but also in axial direction of the curved panel. This study presents a 2-D six-parameter power-law distribution for CNTs volume fraction of 2-D continuously graded nanocomposite that gives designers a powerful tool for flexible designing of structures under multi-functional requirements. The benefit of using generalized power-law distribution is to illustrate and present useful results arising from symmetric, asymmetric and classic profiles. The material properties are determined in terms of local volume fractions and material properties by Mori-Tanaka scheme. The 2-D differential quadrature method as an efficient numerical tool is used to discretize governing equations and to implement boundary conditions. The fast rate of convergence of the method is shown and results are compared against existing results in literature. Some new results for natural frequencies of the shell are prepared, which include the effects of elastic coefficients of foundation, boundary conditions, material and geometrical parameters. The interesting results indicate that a graded nanocomposite volume fraction in two directions has a higher capability to reduce the natural frequency than conventional 1-D functionally graded nanocomposite materials.

Zubarev's Nonequilibrium Statistical Operator Method in the Generalized Statistics of Multiparticle Systems

NASA Astrophysics Data System (ADS)

Glushak, P. A.; Markiv, B. B.; Tokarchuk, M. V.

2018-01-01

We present a generalization of Zubarev's nonequilibrium statistical operator method based on the principle of maximum Renyi entropy. In the framework of this approach, we obtain transport equations for the basic set of parameters of the reduced description of nonequilibrium processes in a classical system of interacting particles using Liouville equations with fractional derivatives. For a classical systems of particles in a medium with a fractal structure, we obtain a non-Markovian diffusion equation with fractional spatial derivatives. For a concrete model of the frequency dependence of a memory function, we obtain generalized Kettano-type diffusion equation with the spatial and temporal fractality taken into account. We present a generalization of nonequilibrium thermofield dynamics in Zubarev's nonequilibrium statistical operator method in the framework of Renyi statistics.

Noninvasive prediction of left ventricular end-diastolic pressure in patients with coronary artery disease and preserved ejection fraction.

PubMed

Abd-El-Aziz, Tarek A

2012-01-01

The aim of this study was to compare 3 different available methods for estimating left ventricular end-diastolic pressure (LVEDP) noninvasively in patients with coronary artery disease and preserved left ventricular ejection fraction (EF). We used 3 equations for noninvasive estimation of LVEDP: The equation of Mulvagh et al., LVEDP(1) = 46 - 0.22 (IVRT) - 0.10 (AFF) - 0.03 (DT) - (2 ÷ E/A) + 0.05 MAR; the equation of Stork et al., LVEDP(2) = 1.06 + 15.15 × Ai/Ei; and the equation of Abd-El-Aziz, LVEDP(3) = [0.54 (MABP) × (1 - EF)] - 2.23. ( A, A-wave velocity; AFF, atrial filling fraction; Ai, time velocity integral of A wave; DT, deceleration time; E, E-wave velocity; Ei, time velocity integral of E wave; IVRT, isovolumic relaxation time; MABP, mean arterial blood pressure; MAR, time from termination of mitral flow to the electrocardiographic R wave; Ti, time velocity integral of total wave.) LVEDP measured by catheterization was correlated with LVEDP(1) (r = 0.52, P < 0.001), LVEDP(2) (r = 0.31, P < 0.05), and LVEDP(3) (r = 0.81, P < 0.001). The equation described by Abd-El-Aziz, LVEDP = [0.54 MABP × (1 - EF)] - 2.23, appears to be the most accurate, reliable, and easily applied method for estimating LVEDP noninvasively in patients with preserved left ventricular ejection fraction and an LVEDP < 20 mm Hg. Copyright © 2012 Canadian Cardiovascular Society. Published by Elsevier Inc. All rights reserved.

Rogers-Schur-Ramanujan Type Identities for the M(p,p') Minimal Models of Conformal Field Theory

NASA Astrophysics Data System (ADS)

Berkovich, Alexander; McCoy, Barry M.; Schilling, Anne

We present and prove Rogers-Schur-Ramanujan (Bose/Fermi) type identities for the Virasoro characters of the minimal model M(p,p'). The proof uses the continued fraction decomposition of p'/p introduced by Takahashi and Suzuki for the study of the Bethe's Ansatz equations of the XXZ model and gives a general method to construct polynomial generalizations of the fermionic form of the characters which satisfy the same recursion relations as the bosonic polynomials of Forrester and Baxter. We use this method to get fermionic representations of the characters for many classes of r and s.

Existence and discrete approximation for optimization problems governed by fractional differential equations

NASA Astrophysics Data System (ADS)

Bai, Yunru; Baleanu, Dumitru; Wu, Guo-Cheng

2018-06-01

We investigate a class of generalized differential optimization problems driven by the Caputo derivative. Existence of weak Carathe ´odory solution is proved by using Weierstrass existence theorem, fixed point theorem and Filippov implicit function lemma etc. Then a numerical approximation algorithm is introduced, and a convergence theorem is established. Finally, a nonlinear programming problem constrained by the fractional differential equation is illustrated and the results verify the validity of the algorithm.

Fractions: A New Slant on Slope

ERIC Educational Resources Information Center

Cheng, Ivan

2010-01-01

An understanding of fractions is foundational to proportional reasoning and algebra. For example, a direct-variation relationship may include a constant of proportionality that is rational, and linear equations may have solutions that are fractions. Beyond that, a grasp of the concept of proportionality and of rational numbers is necessary for…

On Partial Fraction Decompositions by Repeated Polynomial Divisions

ERIC Educational Resources Information Center

Man, Yiu-Kwong

2017-01-01

We present a method for finding partial fraction decompositions of rational functions with linear or quadratic factors in the denominators by means of repeated polynomial divisions. This method does not involve differentiation or solving linear equations for obtaining the unknown partial fraction coefficients, which is very suitable for either…

How to Obtain the Covariant Form of Maxwell's Equations from the Continuity Equation

ERIC Educational Resources Information Center

Heras, Jose A.

2009-01-01

The covariant Maxwell equations are derived from the continuity equation for the electric charge. This result provides an axiomatic approach to Maxwell's equations in which charge conservation is emphasized as the fundamental axiom underlying these equations.

Numerical Solutions of the Nonlinear Fractional-Order Brusselator System by Bernstein Polynomials

PubMed Central

Khan, Rahmat Ali; Tajadodi, Haleh; Johnston, Sarah Jane

2014-01-01

In this paper we propose the Bernstein polynomials to achieve the numerical solutions of nonlinear fractional-order chaotic system known by fractional-order Brusselator system. We use operational matrices of fractional integration and multiplication of Bernstein polynomials, which turns the nonlinear fractional-order Brusselator system to a system of algebraic equations. Two illustrative examples are given in order to demonstrate the accuracy and simplicity of the proposed techniques. PMID:25485293

High-order fractional partial differential equation transform for molecular surface construction.

PubMed

Hu, Langhua; Chen, Duan; Wei, Guo-Wei

2013-01-01

Fractional derivative or fractional calculus plays a significant role in theoretical modeling of scientific and engineering problems. However, only relatively low order fractional derivatives are used at present. In general, it is not obvious what role a high fractional derivative can play and how to make use of arbitrarily high-order fractional derivatives. This work introduces arbitrarily high-order fractional partial differential equations (PDEs) to describe fractional hyperdiffusions. The fractional PDEs are constructed via fractional variational principle. A fast fractional Fourier transform (FFFT) is proposed to numerically integrate the high-order fractional PDEs so as to avoid stringent stability constraints in solving high-order evolution PDEs. The proposed high-order fractional PDEs are applied to the surface generation of proteins. We first validate the proposed method with a variety of test examples in two and three-dimensional settings. The impact of high-order fractional derivatives to surface analysis is examined. We also construct fractional PDE transform based on arbitrarily high-order fractional PDEs. We demonstrate that the use of arbitrarily high-order derivatives gives rise to time-frequency localization, the control of the spectral distribution, and the regulation of the spatial resolution in the fractional PDE transform. Consequently, the fractional PDE transform enables the mode decomposition of images, signals, and surfaces. The effect of the propagation time on the quality of resulting molecular surfaces is also studied. Computational efficiency of the present surface generation method is compared with the MSMS approach in Cartesian representation. We further validate the present method by examining some benchmark indicators of macromolecular surfaces, i.e., surface area, surface enclosed volume, surface electrostatic potential and solvation free energy. Extensive numerical experiments and comparison with an established surface model indicate that the proposed high-order fractional PDEs are robust, stable and efficient for biomolecular surface generation.

Numerical simulation of hydrogen-air reacting flows in rectangular channels with catalytic surface reactions

NASA Astrophysics Data System (ADS)

Amano, Ryoichi S.; Abou-Ellail, Mohsen M.; Elhaw, Samer; Saeed Ibrahim, Mohamed

2013-09-01

In this work a prediction was numerically modeled for a catalytically stabilized thermal combustion of a lean homogeneous mixture of air and hydrogen. The mixture flows in a narrow rectangular channel lined with a thin coating of platinum catalyst. The solution using an in-house code is based on the steady state partial differential continuity, momentum and energy conservation equations for the mixture and species involved in the reactions. A marching technique is used along the streamwise direction to solve the 2-D plane-symmetric laminar flow of the gas. Two chemical kinetic reaction mechanisms were included; one for the gas phase reactions consisting of 17 elementary reactions; of which 7 are forward-backward reactions while the other mechanism is for the surface reactions—which are the prime mover of the combustion under a lean mixture condition—consisting of 16 elementary reactions. The results were compared with a former congruent experimental work where temperature was measured using thermocouples, while using PLIF laser for measuring water and hydrogen mole fractions. The comparison showed good agreement. More results for the velocities, mole fractions of other species were carried out across the transverse and along the streamwise directions providing a complete picture of overall mechanism—gas and surface—and on the production, consumptions and travel of the different species. The variations of the average OH mole fraction with the streamwise direction showed a sudden increase in the region where the ignition occurred. Also the rate of reactions of the entire surface species were calculated along the streamwise direction and a surface water production flux equation was derived by calculating the law of mass action's constants from the concentrations of hydrogen, oxygen and the rate of formation of water near the surface.

Extraction of Peace River bitumen using supercritical ethane

NASA Astrophysics Data System (ADS)

Rose, Jeffrey Lawrence

2000-10-01

As the supply of conventional crude oil continues to decline, petroleum companies are looking for alternative hydrocarbon sources. The vast reserves of heavy oil and bitumen located in northern Alberta are among the alternatives. The challenge facing engineers is to develop a process for recovering this oil which is economic, efficient and environmentally acceptable. Supercritical fluid extraction is one method being investigated which could potentially meet all of these criteria. In this study, Peace River bitumen was extracted using supercritical ethane. The bitumen was mixed with sand and packed into a semi-batch extractor. Ethane contacted the bitumen/sand mixture and the fraction of the bitumen soluble in the ethane was removed and subsequently collected in a two phase separator. The flow of ethane was such that the experiments were governed by equilibrium and not mass transfer. Experimental temperatures and pressures were varied in order to observe the effect of these parameters on the mass and composition of the extracted material. The extraction yields increased as the temperature decreased and pressure increased. Samples were collected at various time intervals to measure changes in the properties of the extracted bitumen over the duration of the process. As the extraction proceeded, the samples became heavier and more viscous. The bitumen feed was characterised and the experimental data was then modelled using the Peng-Robinson equation of state. The characterisation process involved the distillation of the bitumen into five fractions. The distillation curve and density of each fraction was measured and this data was used in conjunction with correlations to determine the critical properties of the bitumen. Interaction parameters in the equation of state were then optimised until the predicted composition of extracted bitumen matched the experimental results.

Suspension concentration distribution in turbulent flows: An analytical study using fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Kundu, Snehasis

2018-09-01

In this study vertical distribution of sediment particles in steady uniform turbulent open channel flow over erodible bed is investigated using fractional advection-diffusion equation (fADE). Unlike previous investigations on fADE to investigate the suspension distribution, in this study the modified Atangana-Baleanu-Caputo fractional derivative with a non-singular and non-local kernel is employed. The proposed fADE is solved and an analytical model for finding vertical suspension distribution is obtained. The model is validated against experimental as well as field measurements of Missouri River, Mississippi River and Rio Grande conveyance channel and is compared with the Rouse equation and other fractional model found in literature. A quantitative error analysis shows that the proposed model is able to predict the vertical distribution of particles more appropriately than previous models. The validation results shows that the fractional model can be equally applied to all size of particles with an appropriate choice of the order of the fractional derivative α. It is also found that besides particle diameter, parameter α depends on the mass density of particle and shear velocity of the flow. To predict this parameter, a multivariate regression is carried out and a relation is proposed for easy application of the model. From the results for sand and plastic particles, it is found that the parameter α is more sensitive to mass density than the particle diameter. The rationality of the dependence of α on particle and flow characteristics has been justified physically.

Analysis of the acceleration region in a circulating fluidized bed riser operating above fast fluidization velocities

DOE Office of Scientific and Technical Information (OSTI.GOV)

Monazam, E.R.; Shadle, L.J.

2008-11-05

In commercial circulating fluidized bed (CFB) processes the acceleration zone greatly contributes to solids mixing, gas and solids dispersion, and particle residence times. A new analysis was developed to describe the relative gas-solids concentration in the acceleration region of a transport system with air as the fluidizing agent for Geldart-type B particles. A theoretical expression was derived from a drag relationship and momentum and continuity equations to describe the evolution of the gas-solids profile along the axial direction. The acceleration zone was characterized using nondimensional analysis of the continuum equations (balances of masses and momenta) that described multiphase flows. Inmore » addition to acceleration length, the boundary condition for the solids fraction at the bottom of the riser and the fully developed regions were measured using an industrial scale CFB of 0.3 m diameter and 15 m tall. The operating factors affecting the flow development in the acceleration region were determined for three materials of various sizes and densities in core annular and dilute regimes of the riser. Performance data were taken from statistically designed experiments over a wide range of Fr (0.5-39), Re (8-600), Ar (29-3600), load ratio (0.2-28), riser to particle diameter ratio (375-5000), and gas to solids density ratio (138-1381). In this one-dimensional system of equations, velocities and solid fractions were assumed to be constant over any cross section. The model and engineering correlations were compared with literature expressions to assess their validity and range of applicability. These expressions can be used as tools for simulation and design of a CFB riser and can also be easily coupled to a kinetics model for process simulation.« less

A mass-energy preserving Galerkin FEM for the coupled nonlinear fractional Schrödinger equations

NASA Astrophysics Data System (ADS)

Zhang, Guoyu; Huang, Chengming; Li, Meng

2018-04-01

We consider the numerical simulation of the coupled nonlinear space fractional Schrödinger equations. Based on the Galerkin finite element method in space and the Crank-Nicolson (CN) difference method in time, a fully discrete scheme is constructed. Firstly, we focus on a rigorous analysis of conservation laws for the discrete system. The definitions of discrete mass and energy here correspond with the original ones in physics. Then, we prove that the fully discrete system is uniquely solvable. Moreover, we consider the unconditionally convergent properties (that is to say, we complete the error estimates without any mesh ratio restriction). We derive L2-norm error estimates for the nonlinear equations and L^{∞}-norm error estimates for the linear equations. Finally, some numerical experiments are included showing results in agreement with the theoretical predictions.

A Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger-Boussinesq Equations Involving Riesz Fractional Derivative

NASA Astrophysics Data System (ADS)

Saha Ray, S.

2017-09-01

In the present paper the Riesz fractional coupled Schrödinger-Boussinesq (S-B) equations have been solved by the time-splitting Fourier spectral (TSFS) method. This proposed technique is utilized for discretizing the Schrödinger like equation and further, a pseudospectral discretization has been employed for the Boussinesq-like equation. Apart from that an implicit finite difference approach has also been proposed to compare the results with the solutions obtained from the time-splitting technique. Furthermore, the time-splitting method is proved to be unconditionally stable. The error norms along with the graphical solutions have also been presented here. Supported by NBHM, Mumbai, under Department of Atomic Energy, Government of India vide Grant No. 2/48(7)/2015/NBHM (R.P.)/R&D II/11403

Apparatus and method for determining solids circulation rate

DOEpatents

Ludlow, J Christopher [Morgantown, WV; Spenik, James L [Morgantown, WV

2012-02-14

The invention relates to a method of determining bed velocity and solids circulation rate in a standpipe experiencing a moving packed bed flow, such as the in the standpipe section of a circulating bed fluidized reactor The method utilizes in-situ measurement of differential pressure over known axial lengths of the standpipe in conjunction with in-situ gas velocity measurement for a novel application of Ergun equations allowing determination of standpipe void fraction and moving packed bed velocity. The method takes advantage of the moving packed bed property of constant void fraction in order to integrate measured parameters into simultaneous solution of Ergun-based equations and conservation of mass equations across multiple sections of the standpipe.

Combinatorial interpretation of Haldane-Wu fractional exclusion statistics.

PubMed

Aringazin, A K; Mazhitov, M I

2002-08-01

Assuming that the maximal allowed number of identical particles in a state is an integer parameter, q, we derive the statistical weight and analyze the associated equation that defines the statistical distribution. The derived distribution covers Fermi-Dirac and Bose-Einstein ones in the particular cases q=1 and q--> infinity (n(i)/q-->1), respectively. We show that the derived statistical weight provides a natural combinatorial interpretation of Haldane-Wu fractional exclusion statistics, and present exact solutions of the distribution equation.

Application of the enhanced homotopy perturbation method to solve the fractional-order Bagley-Torvik differential equation

NASA Astrophysics Data System (ADS)

Zolfaghari, M.; Ghaderi, R.; Sheikhol Eslami, A.; Ranjbar, A.; Hosseinnia, S. H.; Momani, S.; Sadati, J.

2009-10-01

The enhanced homotopy perturbation method (EHPM) is applied for finding improved approximate solutions of the well-known Bagley-Torvik equation for three different cases. The main characteristic of the EHPM is using a stabilized linear part, which guarantees the stability and convergence of the overall solution. The results are finally compared with the Adams-Bashforth-Moulton numerical method, the Adomian decomposition method (ADM) and the fractional differential transform method (FDTM) to verify the performance of the EHPM.

Adomian decomposition

NASA Astrophysics Data System (ADS)

Daftardar-Gejji, Varsha; Jafari, Hossein

2005-01-01

Adomian decomposition method has been employed to obtain solutions of a system of fractional differential equations. Convergence of the method has been discussed with some illustrative examples. In particular, for the initial value problem: where A=[aij] is a real square matrix, the solution turns out to be , where E([alpha]1,...,[alpha]n),1 denotes multivariate Mittag-Leffler function defined for matrix arguments and Ai is the matrix having ith row as [ai1...ain], and all other entries are zero. Fractional oscillation and Bagley-Torvik equations are solved as illustrative examples.

Stability of standing wave for the fractional nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Peng, Congming; Shi, Qihong

2018-01-01

In this paper, we study the stability and instability of standing waves for the fractional nonlinear Schrödinger equation i∂tu = (-Δ)su - |u|2σu, where (t ,x ) ∈R × RN, 1/2