Next Article in Journal
Hyperspectral Imaging System with Rotation Platform for Investigation of Jujube Skin Defects
Previous Article in Journal
Automated Data Acquisition in Construction with Remote Sensing Technologies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation for Fractional-Order Bloch Equation Arising in Nuclear Magnetic Resonance by Using the Jacobi Polynomials

by
Harendra Singh
1 and
H. M. Srivastava
2,3,4,*
1
Department of Mathematics, Post-Graduate College, Ghazipur 233001, Uttar Pradesh, India
2
Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8W 3R4, Canada
3
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
4
Department of Mathematics and Informatics, Azerbaijan University, 71 Jeyhun Hajibeyli Street, Baku AZ1007, Azerbaijan
*
Author to whom correspondence should be addressed.
Submission received: 17 March 2020 / Revised: 10 April 2020 / Accepted: 13 April 2020 / Published: 20 April 2020

Abstract

:
In the present paper, we numerically simulate fractional-order model of the Bloch equation by using the Jacobi polynomials. It arises in chemistry, physics and nuclear magnetic resonance (NMR). It also arises in magnetic resonance imaging (MRI) and electron spin resonance (ESR). It is used for purity determination, provided that the molecular weight and structure of the compound is known. It can also be used for structural determination. By the study of NMR, chemists can determine the structure of many compounds. The obtained numerical results are compared and simulated with the known solutions. Accuracy of the proposed method is shown by providing tables for absolute errors and root mean square errors. Different orders of the time-fractional derivatives results are illustrated by using figures.

1. Introduction

The Bloch equation is a system of differential equations. It is mainly valuable for studying expensive biological samples like RNA, DNA, proteins and nucleic acids. It has many real-life applications like process control, liquid media, petrochemical plants and process optimization in oil refineries. Surface magnetic resonance is based on the principle of NMR, and the measurements can be used to indirectly estimate the water content of the saturated and unsaturated zones. The standard system of Bloch equations is given as follows:
d N x ( t ) d t = w 0 N y ( t ) N x ( t ) T 2 d N y ( t ) d t = w 0 N x ( t ) N y ( t ) T 2
d N z ( t ) d t = N 0 N z ( t ) T 1 ,
with the initial conditions N x ( 0 ) = a 1 ,   N y ( 0 ) = a 2 and N z ( 0 ) = a 3 .
Here N x ( t ) ,   N y ( t ) and N z ( t ) denote the system magnetisation in x ,   y   and   z components, respectively; w 0 is the resonant frequency given by the Larmor relationship w 0 = γ M 0 , where M 0 is the static magnetic field in z component, N 0 is the equilibrium magnetisation, T 1 and T 2 are the spin-lattice and spin-spin relaxation time, respectively, and a 1 ,   a 2 and a 3 are real constants. The analytical solution of Equation (1) is given by
N x ( t ) = e t T 2 ( N x ( 0 ) cos w 0 t + N y ( 0 ) sin w 0 t ) N y ( t ) = e t T 2 ( N y ( 0 ) cos w 0 t N x ( 0 ) sin w 0 t )
N z ( t ) = N z ( 0 ) e t T 1 + N 0 ( 1 e t T 1 )
Real-life applications of fractional calculus are in subjects like biology [1], viscoelasticity [2,3,4], signal processing [5], control theory [6], fluid dynamics [7]. For more details, the reader should refer to [8]. There are many magnetic resonance systems which are modelled by the fractional-order Bloch equation, and it is well known that the fractional-order derivatives are non-local in nature. Therefore, we will replace the integer-order Bloch equation by the fractional-order Bloch equation with a view to further understand the resulting magnetic resonance systems. Therefore, we replace the integer-order time-derivative by the non-integer- order time-derivative:
d α N x ( t ) d t α = w 0 N y ( t ) N x ( t ) T 2 d β N y ( t ) d t β = w 0 N x ( t ) N y ( t ) T 2
d γ N z ( t ) d t γ = N 0 N z ( t ) T 1 ,
where 0 < α , β , γ 1 .
The non-integer-order derivative is in the Liouville–Caputo (LC) sense. The LC non-integer-order derivative of order β is defined as follows [8]:
D β f ( x ) = I l β D l f ( x ) = 1 Γ ( l β ) 0 x ( x t ) l β 1 d l d t l f ( t ) d t ,   l 1 < β < l ,   x > 0 .
In this paper, we are considering that β ( 0 ,   1 ) ; therefore, we will take l = 1 . The time-fractional derivatives play a key role upsetting the spin dynamics defined by the Bloch equations in Equation (3) (see [9,10]). The magnetic resonance components of the magnetisation are identified in the initial state of the system, and hence, these should be visibly predictable. The physical meaning of the non-integer order Bloch equations can be understood in the basic preparation of the non-integer-order Schrödinger equation.
Bloch equations in NMR can be simulated numerically and analytically (see, for details, [11,12,13,14,15,16]). The time-fractional order Bloch equation having fractional derivative in Caputo sense is solved in [17]. Recently, Kumar et al. [18] solved fractional-order Bloch equation by using homotopy perturbation method (HPM). Use has been made of operational matrix method with Legendre polynomials in [19] and with the Laguerre polynomials in [20] for the solution of this equation. In [21], this equation was solved numerically by using the iterative method. Furthermore, in [22], by using numerical approximation, a special class of this equation, namely the fuzzy time-fractional Bloch equation, was solved. In this paper we propose to solve the fractional-order Bloch equation by using the Jacobi polynomials. Some developments on orthogonal approximations can be found in [23,24,25,26,27,28,29,30]. Some introductory overview and recent development of fractional calculus can be seen in [31]. In this method, we get unknown coefficients for the approximated parameter in the model and, by the use of these coefficients, we obtain approximate solutions of the fractional-order Bloch equation in NMR.

2. Preliminaries

The Jacobi polynomial of degree i on [0, 1] is given by [28]
σ i ( t ) = k = 0 i ( 1 ) i k Γ ( i + b + 1 ) Γ ( i + k + a + b + 1 ) Γ ( k + b + 1 ) Γ ( i + a + b + 1 ) ( i k ) ! k ! t k
The orthogonal property of the Jacobi polynomials with respect to the weight function w ( a , b ) ( t ) = ( 1 t ) a t b is given by
0 1 σ n ( t ) σ m ( t ) w ( a , b ) ( t ) d t = v n a , b δ m n
where δ m n   is Kronecker delta function and
v n a , b = Γ ( n + a + 1 ) Γ ( n + b + 1 ) ( 2 n + a + b + 1 ) n ! Γ ( n + a + b + 1 )
A function f L 2 [ 0 ,   1 ] , with | f ( t ) | Q , can be expanded as follows:
f ( t ) = lim n i = 0 n c i σ i ( t ) ,
where c i = 1 v i a ,   b 0 1 σ i ( t ) f ( t ) w ( a , b ) ( t ) d t .
Equation (7), for finite dimensional approximation, is written in the following form:
f i = 0 m c i σ i ( t ) = C T q m ( t ) ,
where C and q m ( t ) are ( m + 1 ) × 1 matrices given by C = [ c 0 , c 1 , . , c m ] T and q m ( t ) = [ σ 0 , σ 1 , . , σ m ] T .
Theorem 1. 
If q n ( t ) = [ σ 0 , σ 1 , . , σ n ] T denotes the shifted Jacobi vector and if v > 0 , then I v σ i ( t ) = I ( v ) q n ( t ) , where I ( v ) = ( u ( i , j ) ) , is the ( n + 1 ) × ( n + 1 ) operational matrix of fractional integral of order v , and its ( i , j ) th entry is given by
u ( i , j ) = k = 0 i l = 0 j ( 1 ) i + j k l Γ ( a + 1 ) Γ ( i + b + 1 ) Γ ( i + k + a + b + 1 ) Γ ( j + l + a + b + 1 ) Γ ( v + k + l + a + b + 1 ) ( 2 j + a + b + 1 ) j ! ( i k ) ! ( j l ) ! ( l ) !   Γ ( k + b + 1 ) Γ ( i + a + b + 1 ) Γ ( v + k + 1 ) Γ ( j + a + 1 ) Γ ( l + b + 1 ) Γ ( k + l + v + a + b + 1 )
Proof. 
Please see [28]. □

3. Construction of Algorithm

In this section, we construct an algorithm to get the approximate solution of the Bloch equation. Using this algorithm, we can then obtain magnetisation in each direction.
Let us take the following approximations:
d α N x ( t ) d t α = C 1 T q ( t ) ,   d β N y ( t ) d t β = C 2 T q ( t ) ,   d γ N z ( t ) d t γ = C 3 T q ( t ) ,
and
N x ( 0 ) = L T q ( t ) ,   N y ( 0 ) = M T q ( t ) ,   N z ( 0 ) = N T q ( t ) ,   N 0 T 1 = O T q ( t ) .
From Equations (10) and (11), we can write
N x ( t ) = C 1 T I ( α ) q ( t ) + L T q ( t ) ,
N y ( t ) = C 2 T I ( β ) q ( t ) + M T q ( t ) ,
N z ( t ) = C 3 T I ( γ ) q ( t ) + N T q ( t )
Using Equations (10), (12), (13) and (14) in Equation (3), we get
C 1 T ( I + 1 T 2 I ( α ) ) w 0 C 2 T I ( β ) = w 0 M T 1 T 2 L T
w 0 C 1 T I ( α ) + C 2 T ( I + 1 T 2 I ( β ) ) = w 0 L T 1 T 2 M T
C 3 T ( I + 1 T 1 I ( γ ) ) = O T 1 T 1 N T
where I ( α ) ,   I ( β )   and   I ( γ ) are the operational matrices of non-integer-order integration of order α , β   and   γ , respectively. Here I is an identity matrix.
The simpler form for Equations (15)–(17) is given as follows:
C 1 T W 1 C 2 T W 5 = G 1 ,
C 1 T W 4 + C 2 T W 2 = G 2 ,
C 3 T W 3 = G 3 ,
where
W 1 = I + 1 T 2 I ( α ) ,
W 2 = I + 1 T 2 I ( β ) ,
W 3 = I + 1 T 1 I ( γ ) ,
W 4 = w 0 I ( α ) ,
W 5 = w 0 I ( β ) ,
G 1 = w 0 M T 1 T 2 L T ,
G 2 = w 0 L T 1 T 2 M T ,
G 3 = O T 1 T 1 N T
The matrices W 1 ,   W 2 ,   W 3 , W 4 ,   W 5 ,   G 1 ,   G 2   and G 3 are given in terms of known values, so these matrices are known matrices.
On solving Equations (18)–(20), we get
C 1 T = ( G 1 W 5 1 + G 2 W 2 1 ) ( W 1 W 5 1 + W 4 W 2 1 ) 1 ,
C 2 T = { ( G 1 W 5 1 + G 2 W 2 1 ) ( W 1 W 5 1 + W 4 W 2 1 ) 1 W 1 G 1 } W 5 1 ,
C 3 T = G 3 W 3 1
Using Equations (29)–(31) in Equations (12)–(14), respectively, we get the system magnetisation N x ( t ) ,   N y ( t ) and N z ( t ) for Bloch equations in NMR.

4. Convergence Analysis

Theorem 2. 
If d α N x d t α   ϵ C ( m + 1 ) [ 0 , 1 ] and R m ( d α N x d t α ) are the m t h approximations of d α N x d t α by using P m ( t ) = s p a n { σ 0 ( t ) , σ 1 ( t ) , . , σ m ( t ) } , then lim m d α N x d t α R m ( d α N x d t α ) w ( a , b ) 0 .
Proof. 
Since d α N x d t α   ϵ C ( m + 1 ) [ 0 , 1 ] , so the Taylor polynomial of d α N x d t α at t = 0 is given as follows:
M 1 ( t ) = ( d α N x d t α ) t = 0 + ( d α N x d t α ) t = 0 t + + ( d α N x d t α ) t = 0 m t m m ! .
The upper bound of the error of the Taylor polynomial is given by
| d α N x d t α M 1 ( t ) | K t m + 1 ( m + 1 ) ! ,
where
K = max t [ 0 , 1 ] | ( d α N x d t α ) m + 1 ( t ) |
Since R m ( d α N x d t α ) and M 1 P m , we have
d α N x d t α R m ( d α N x d t α ) w ( a , b ) 2 d α N x d t α M 1 w ( a , b ) 2 ,   ( K ( m + 1 ) ! ) 2 0 1 t 2 m + 2 + b ( 1 t ) a d t , = ( K ( m + 1 ) ! ) 2 ( 1 + a ) ( 3 + 2 m + b ) ( 4 + 2 m + a + b ) , d α N x d t α R m ( d α N x d t α ) w ( a , b ) K ( m + 1 ) ! ( 1 + a ) ( 3 + 2 m + b ) ( 4 + 2 m + a + b )
Taking m in Equation (35), we get
lim m d α N x d t α R m ( d α N x d t α ) w ( a , b ) 0 .

5. Numerical Results and Discussion

In this section, we will numerically simulate our results with known results. For each numerical simulation, we will consider i. c. N x ( 0 ) = 0 ,   N y ( 0 ) = 100 and N z ( 0 ) = 0 . In Figure 1 and Figure 2, we have presented 3D and 2D plots of the numerical solutions of the Bloch equation for integer order, respectively. These figures show the dynamics of N x ,   N y   and   N z for integer-order relaxation. In Figure 1, the entire trajectory of magnetisation is shown in 3D for integer order starting at i. c. ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 2, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 3 and Figure 4, we have presented 3D and 2D plots of the numerical solutions of the Bloch equation for fractional order ( α = β = γ = 0.9 ) , respectively. These figures show the dynamics of N x ,   N y   and   N z for fractional-order relaxation. In Figure 3, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = β = γ = 0.9 ) starting at i. c. ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 4, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 5 and Figure 6, we have presented the 3D and 2D plots of the numerical solutions of the Bloch equation for fractional order ( α = β = γ = 0.8 ) , respectively. These figures show the dynamics of N x ,     N y   and   N z for the fractional-order relaxation. In Figure 5, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = β = γ = 0.8 ) starting at i. c. ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 6, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 7 and Figure 8, we have presented the 3D and 2D plots of the numerical solutions of the Bloch equation for fractional order ( α = β = γ = 0.7 ) , respectively. These figures show the dynamics of N x ,   N y   and   N z for the fractional-order relaxation. In Figure 7, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = β = γ = 0.7 ) with the starting initially ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 8, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 9 and Figure 10, we have presented the 3D and 2D plots of the numerical solutions of the Bloch equation for fractional order ( α = 1.0 ,     β = 0.9 ,     γ = 0.8 ) , respectively. These figures show the dynamics of N x ,   N y   and   N z for the fractional-order relaxation. In Figure 9, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = 1.0 ,     β = 0.9 ,     γ = 0.8 ) starting at i. c. ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 10, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 11 and Figure 12 we have presented 3D and 2D plots of numerical solutions of Bloch equation for fractional order ( α = 0.9 ,     β = 0.9 ,     γ = 1.0 ) , respectively. These figures show the dynamic of N x ,     N y   and   N z for fractional order relaxation. In Figure 11, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = 0.9 ,     β = 0.9 ,     γ = 1.0 ) starting at i.c. ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 12, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 13 and Figure 14, we have presented the 3D and 2D plots of the numerical solutions of the Bloch equation for fractional order ( α = 1.0 ,     β = 1.0 ,     γ = 0.9 ) , respectively. These figures show the dynamics of N x ,     N y   and   N z for the fractional-order relaxation. In Figure 13, the entire trajectory of magnetisation is shown in 3D for fractional order ( α = 1.0 ,     β = 1.0 ,     γ = 0.9 ) starting at the initial level ( N x ( 0 ) , N y ( 0 ) , N z ( 0 ) ) and returning to N 0 . From Figure 14, it is clear that the magnetisation N x in x direction increases with time, and the magnetisation N y in y direction decreases with time.
In Figure 15 and Figure 16, we have shown the numerical simulation of analytical and numerical solutions of the Bloch equation for N x ( t ) and N y ( t ) for integer order, respectively. From Figure 16, it is clear that the solution has periodic behaviour at low frequency. This solution varies periodically for N x ( t ) and N y ( t ) at low frequency.
In Figure 17, Figure 18 and Figure 19, we have shown the absolute errors for N x ( t ) ,   N y ( t )   and   N z ( t ) , respectively, at different values of m = 3 ,   6   and   9 . In Figure 17, Figure 18 and Figure 19, the absolute errors are denoted by E 1 , E 2 and E 3 for m = 3 ,   6 and 9, respectively. In all these figures, E 2 and E 3 are multiplied by 10 4 and 10 5 , respectively.
From these figures, we can see that the errors decrease with the increase of m . In Figure 20, Figure 21 and Figure 22, we have shown the behaviour of the solutions of N x ( t ) ,     N y ( t )   and   N z ( t ) at different values of α ,     β   and   γ ,   respectively. In Figure 20, Figure 21 and Figure 22, exact solution means the analytical solution for integer order ( α = β = γ = 1 ) Bloch equation as given by Equation (2).
From these figures it is clear that the solution varies consistently from non-integer order to integer order.
In Table 1, we have listed the maximum absolute errors ( l ) and the root-mean-square errors ( l 2 ) for two different values of m = 4 and 8 at α = β = γ = 1 ,   w 0 = 1 ,   T 1 = 1 ,   T 2 = 20 ,   a = b = 0.8 . We have calculated these errors for integer order by taking the exact solution as given by Equation (2).
From Table 1, it is detected that the errors decrease with the increase of m .

6. Conclusions and Future Scope

In this paper, we have presented the numerical solution and the simulation for fractional-order and integer-order Bloch equations. Mathematical model for NMR allows us to explore and define magnetisation for spin dynamics at resonance frequency in a static magnetic field. Implementation of our proposed technique is easy in comparison to the existing methods because the operational matrices are easy to construct. The numerical section shows how the solution given by the used technique varies consistently at different values of non-integer-order time-derivatives. Moreover, for integer order, the solution by the used technique is identical to the exact solution for the Bloch equation. The error table shows the accuracy of the proposed method. For future work, we can construct operational matrices for different polynomials in order to attain better exactness.

Author Contributions

Both authors have equal contribution. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bagley, R.L.; Torvik, P.J. A theoretical basis for the application of fractional calculus to viscoelasticity. J. Rheol. 1983, 27, 201–210. [Google Scholar] [CrossRef]
  2. Srivastava, H.M.; Shah, F.A.; Abass, R. An application of the Gegenbauer wavelet method for the numerical solution of the fractional Bagley-Torvik equation. Russ. J. Math. Phys. 2019, 26, 77–93. [Google Scholar] [CrossRef]
  3. Bagley, R.L.; Torvik, P.J. Fractional calculus in the transient analysis of viscoelasticity damped structures. AIAA J. 1985, 23, 918–925. [Google Scholar] [CrossRef]
  4. Robinson, A.D. The use of control systems analysis in neurophysiology of eye movements. Annu. Rev. Neurosci. 1981, 4, 462–503. [Google Scholar] [CrossRef] [PubMed]
  5. Singh, H. A new stable algorithm for fractional Navier-Stokes equation in polar coordinate. Int. J. Appl. Comput. Math. 2017, 3, 3705–3722. [Google Scholar] [CrossRef]
  6. Bohannan, G.W. Analog fractional order controller in temperature and motor control applications. J. Vib. Control 2008, 14, 1487–1498. [Google Scholar] [CrossRef]
  7. Panda, R.; Dash, M. Fractional generalized splines and signal processing. Signal Process 2006, 86, 2340–2350. [Google Scholar] [CrossRef]
  8. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, North-Holland Mathematical Studies; Elsevier (North-Holland) Science Publishers: Amsterdam, The Netherlands, 2006. [Google Scholar]
  9. Li, X. Numerical solution of fractional partial differential equations using cubic B-spline wavelet collocation method. Adv. Comput. Math. Appl. 2012, 1, 159–164. [Google Scholar]
  10. Daftardar-Gejji, V.; Bhalekar, S. Solving fractional diffusion-wave equations using a new iterative method. Fract. Calc. Appl. Anal. 2008, 11, 193–202. [Google Scholar]
  11. Awojoyogber, O.B. Analytical solution of the time dependent Bloch NMR, flow equations: A translational mechanical analysis. Phys. A Stat. Mech. Appl. 2004, 339, 437–460. [Google Scholar] [CrossRef]
  12. Murase, K.; Tanki, N. Numerial solution to the time dependent Bloch equations revisited. Magn. Reson. Imaging 2011, 29, 126–131. [Google Scholar] [CrossRef] [PubMed]
  13. Leyte, J.C. Some solutions of the Bloch equations. Chem. Phys. Lett. 1990, 165, 231–240. [Google Scholar] [CrossRef]
  14. Yan, H.; Chen, B.; Gore, J.C. Approximate solutions of the Bloch equations for selective excitation. J. Magn. Reson. 1987, 75, 83–95. [Google Scholar] [CrossRef]
  15. Hoult, D.I. The solution of the Bloch equation in presence of varying B 1 field—An approach to selective pulse analysis. J. Magn. Reson. 1979, 35, 69–86. [Google Scholar] [CrossRef]
  16. Xu, Z.; Chan, A.K. A Near-Resonance solution to the Bloch equations and its application to RF pulse design. J. Magn. Reson. 1999, 138, 225–231. [Google Scholar] [CrossRef]
  17. Magin, R.; Feng, X.; Baleanu, D. Solving the fractional order Bloch Equation. Concepts Magn. Reson. Part A 2009, 34, 16–23. [Google Scholar] [CrossRef]
  18. Kumar, S.; Faraz, N.; Sayevand, K. A fractional model of Bloch equation in NMR and its analytic approximate solution. Walailak J. Sci. Technol. 2014, 11, 273–285. [Google Scholar]
  19. Singh, H. A New Numerical Algorithm for Fractional Model of Bloch equation in nuclear magnetic resonance. Alex. Eng. J. 2016, 55, 2863–2869. [Google Scholar] [CrossRef] [Green Version]
  20. Singh, H. Operational matrix approach for approximate solution of fractional model of Bloch equation. J. King Saud Univ. Sci. 2017, 29, 235–240. [Google Scholar] [CrossRef]
  21. Petráš, L. Modeling and numerical analysis of fractional-order Bloch equations. Comput. Math. Appl. 2011, 61, 341–356. [Google Scholar] [CrossRef]
  22. Ahmadian, A.; Chan, C.S.; Salahshour, S.; Vaitheeswaran, V. FTFBE: A numerical approximation for fuzzy time-fractional Bloch equation. In Proceedings of the International conference on fuzzy systems (FUZZ-IEEE), World Congress on Computational Intelligence, Beijing, China, 6–11 July 2014; pp. 418–423. [Google Scholar]
  23. Wu, J.L. A wavelet operational method for solving fractional partial differential equations numerically. Appl. Math. Comput. 2009, 214, 31–40. [Google Scholar] [CrossRef]
  24. Singh, H.; Srivastava, H.M.; Kumar, D. A reliable numerical algorithm for the fractional vibration equation. Chaos Solitons Fractals 2017, 103, 131–138. [Google Scholar] [CrossRef]
  25. Tohidi, E.; Bhrawy, A.H.; Erfani, K. A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation. Appl. Math. Model. 2013, 37, 4283–4294. [Google Scholar] [CrossRef]
  26. Kazem, S.; Abbasbandy, S.; Kumar, S. Fractional order Legendre functions for solving fractional-order differential equations. Appl. Math. Model. 2013, 37, 5498–5510. [Google Scholar] [CrossRef]
  27. Singh, C.S.; Singh, H.; Singh, V.K.; Singh, O.P. Fractional order operational matrix methods for fractional singular integro-differential equation. Appl. Math. Model. 2016, 40, 10705–10718. [Google Scholar] [CrossRef]
  28. Singh, H.; Srivastava, H.M. Jacobi collocation method for the approximate solution of some fractional-order Riccati differential equations with variable coefficients. Phys. A Stat. Mech. Appl. 2019, 523, 1130–1149. [Google Scholar] [CrossRef]
  29. Singh, C.S.; Singh, H.; Singh, S.; Kumar, D. An efficient computational method for solving system of nonlinear generalized Abel integral equations arising in astrophysics. Phys. A Stat. Mech. Appl. 2019, 525, 1440–1448. [Google Scholar] [CrossRef]
  30. Singh, H.; Pandey, R.K.; Baleanu, D. Stable numerical approach for fractional delay differential equations. Few-Body Syst. 2017, 58, 156. [Google Scholar] [CrossRef]
  31. Srivastava, H.M. Fractional-order derivatives and integrals: Introductory overview and recent developments. Kyungpook Math. J. 2020, 60, 73–116. [Google Scholar]
Figure 1. Numerical solutions of the Bloch equation with parameters: α = β = γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,     b = 1 .
Figure 1. Numerical solutions of the Bloch equation with parameters: α = β = γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,     b = 1 .
Applsci 10 02850 g001
Figure 2. Numerical solutions of the Bloch equation in plane ( N x   vs .   N y ) with parameters: α = β = γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 2. Numerical solutions of the Bloch equation in plane ( N x   vs .   N y ) with parameters: α = β = γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g002
Figure 3. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 3. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g003
Figure 4. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = β = γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 4. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = β = γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g004
Figure 5. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 5. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g005
Figure 6. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = β = γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 6. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = β = γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g006
Figure 7. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.7 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 7. Numerical solutions of the Bloch equation with parameters: α = β = γ = 0.7 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g007
Figure 8. Numerical solutions of the Bloch equation in the plane ( N x   vs .     N y ) with parameters: α = β = γ = 0.7 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 8. Numerical solutions of the Bloch equation in the plane ( N x   vs .     N y ) with parameters: α = β = γ = 0.7 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g008
Figure 9. Numerical solutions of the Bloch equation with parameters: α = 1.0 ,   β = 0.9 ,   γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Figure 9. Numerical solutions of the Bloch equation with parameters: α = 1.0 ,   β = 0.9 ,   γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Applsci 10 02850 g009
Figure 10. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 1.0 ,   β = 0.9 ,   γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Figure 10. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 1.0 ,   β = 0.9 ,   γ = 0.8 , w 0 = 12 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Applsci 10 02850 g010
Figure 11. Numerical solutions of the Bloch equation with parameters: α = 0.9 ,     β = 0.9 ,     γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 11. Numerical solutions of the Bloch equation with parameters: α = 0.9 ,     β = 0.9 ,     γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g011
Figure 12. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 0.9 ,   β = 0.9 ,   γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 12. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 0.9 ,   β = 0.9 ,   γ = 1 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g012
Figure 13. Numerical solutions of the Bloch equation with parameters: α = 1 ,     β = 1 ,   γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,     b = 1 .
Figure 13. Numerical solutions of the Bloch equation with parameters: α = 1 ,     β = 1 ,   γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,     b = 1 .
Applsci 10 02850 g013
Figure 14. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 1 ,   β = 1 ,   γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Figure 14. Numerical solutions of the Bloch equation in the plane ( N x   vs .   N y ) with parameters: α = 1 ,   β = 1 ,   γ = 0.9 , w 0 = 12 , T 1 = 1 , T 2 = 20 , a = 1 ,   b = 1 .
Applsci 10 02850 g014
Figure 15. Numerical simulation of solutions of the Bloch equation for N x ( t ) with parameters: α = β = γ = 1 , w 0 = 14 , T 1 = 1 , T 2 = 20 ,   a = 0.9 ,   b = 0.9 .
Figure 15. Numerical simulation of solutions of the Bloch equation for N x ( t ) with parameters: α = β = γ = 1 , w 0 = 14 , T 1 = 1 , T 2 = 20 ,   a = 0.9 ,   b = 0.9 .
Applsci 10 02850 g015
Figure 16. Numerical simulation of solutions of the Bloch equation for N y ( t ) with parameters: α = β = γ = 1 , w 0 = 14 , T 1 = 1 , T 2 = 20 ,   a = 0.9 ,   b = 0.9 .
Figure 16. Numerical simulation of solutions of the Bloch equation for N y ( t ) with parameters: α = β = γ = 1 , w 0 = 14 , T 1 = 1 , T 2 = 20 ,   a = 0.9 ,   b = 0.9 .
Applsci 10 02850 g016
Figure 17. Errors for N x ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Figure 17. Errors for N x ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Applsci 10 02850 g017
Figure 18. Errors for N y ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Figure 18. Errors for N y ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Applsci 10 02850 g018
Figure 19. Errors for N z ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Figure 19. Errors for N z ( t ) at m = 3 ,   6   and   9 , with parameters: α = β = γ = 1 , w 0 = 1 , T 1 = 1 , T 2 = 20 , a = 1 , b = 1 .
Applsci 10 02850 g019
Figure 20. Behaviour of the approximate solution of N x ( t ) at α = 0.6 ,   0.7 ,   0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Figure 20. Behaviour of the approximate solution of N x ( t ) at α = 0.6 ,   0.7 ,   0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Applsci 10 02850 g020
Figure 21. Behaviour of the approximate solution of N y ( t ) at β = 0.6 ,   0.7   , 0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Figure 21. Behaviour of the approximate solution of N y ( t ) at β = 0.6 ,   0.7   , 0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Applsci 10 02850 g021
Figure 22. Behaviour of the approximate solution of N z ( t ) at γ = 0.6 ,   0.7 ,   0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Figure 22. Behaviour of the approximate solution of N z ( t ) at γ = 0.6 ,   0.7 ,   0.8 ,   0.9   and   1 , with parameters: w 0 = 1 , T 1 = 1 , T 2 = 20 ,   a = 1 ,   b = 1 .
Applsci 10 02850 g022
Table 1. Comparison of ( l ) and ( l 2 ) errors at a = b = 0.8 , m = 5 ,   8 for integer order solution.
Table 1. Comparison of ( l ) and ( l 2 ) errors at a = b = 0.8 , m = 5 ,   8 for integer order solution.
N i ( t ) m = 5
l - Error
l 2 - Error
m = 8
l - Error
l 2 - Error
N x ( t ) 2.0349 × 10−44.2101 × 10−8
5.7209 × 10−69.2050 × 10−10
N y ( t ) 1.9648 × 10−46.2953 × 10−10
5.5559 × 10−66.2953 × 10−10
N z ( t ) 1.7733 × 10−65.7426 × 10−10
5.0010 × 10−86.8075 × 10−12

Share and Cite

MDPI and ACS Style

Singh, H.; Srivastava, H.M. Numerical Simulation for Fractional-Order Bloch Equation Arising in Nuclear Magnetic Resonance by Using the Jacobi Polynomials. Appl. Sci. 2020, 10, 2850. https://0-doi-org.brum.beds.ac.uk/10.3390/app10082850

AMA Style

Singh H, Srivastava HM. Numerical Simulation for Fractional-Order Bloch Equation Arising in Nuclear Magnetic Resonance by Using the Jacobi Polynomials. Applied Sciences. 2020; 10(8):2850. https://0-doi-org.brum.beds.ac.uk/10.3390/app10082850

Chicago/Turabian Style

Singh, Harendra, and H. M. Srivastava. 2020. "Numerical Simulation for Fractional-Order Bloch Equation Arising in Nuclear Magnetic Resonance by Using the Jacobi Polynomials" Applied Sciences 10, no. 8: 2850. https://0-doi-org.brum.beds.ac.uk/10.3390/app10082850

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop