A novel difference schemes for analyzing the fractional Navier- Stokes equations

DOI: 0.55/auom-207-005 An. Şt. Univ. Ovidius Constanţa Vol. 25(),207, 95 206 A novel difference schemes for analyzing the fractional Navier- Stokes equations Khosro Sayevand, Dumitru Baleanu, Fatemeh Sahsavand Abstract In this report, a novel difference scheme is used to analyzing the Navier - Stokes problems of fractional order. Existence and uniqueness of the suggested approach with a Lipschitz condition and Picard theorem are proved. Furthermore, we find a discrete analogue of the derivative and then stability and convergence of our strategy in multi dimensional domain are proved. Introduction The theory of fractional calculus deals with derivatives and integrals of any arbitrary orders. Fractional calculus have gained importance, mainly due to their demonstrated applications in many area of physics, economics and engineering. The fractional calculus has been occurring in many physical problems such as damping law, perfusion processes and also motion of a large thin plate in a Newtonian fluid. For more details on the scientific applications of fractional calculus, see [, 2]. Difference schemes [3 5] (grid methods) is explained at numerical approximation of various problems in real world application. Under such an approach the solution of differential equations amounts for solving the systems of algebraic equations. Difference scheme is based on composition of discrete (difference) approximations to equations of mathematics and verifying a priori Key Words: Fractional calculus, Difference scheme, Navier - Stokes equations, Riemann Liouville fractional derivative 200 Mathematics Subject Classification: Primary 34A08; Secondary 49S05. Received: 2.0.205 Accepted:.03.206 95 Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 96 quality characteristics of these approximations, mainly stability, convergence, error of approximation and accuracy of the difference schemes obtained. We were interested in approximate substitutions of difference operators for differential ones. However, many problems of physics [6, 7] involve not only differential equations, but also the supplementary conditions such as boundary and initial which guide a proper choice of a unique solution from the collection of possible solutions. This work deals with introducing a novel difference schemes for Navier- Stokes fractional differential equations in multi dimensional domains [8, 9]. A fractional Navier-Stokes equations in cylindrical coordinates is in the following form: D α 0tu(r, t) = p + v(d 2 r + r D ru), 0 < α, () u is the velocity, r is the distance along an x-axis, p is the pressure, v is the kinematics viscosity. Here the fractional derivative (D) is described in the Riemann Liouville form. As we know the fractional Navier - Stokes equations are nonlinear. Therefore there is no known methodology to solve these equations and there are very few cases in which the exact solution of the time fractional Navier - Stokes equations can be achieved. 2 Preliminaries and basic concepts 2. Cauchy problem for fractional ordinary differential equation Consider the following Cauchy problem in the form of fractional order α 0tu = D α 0tu α 0tu = f(u, t), u(0) = u 0, (2) u(0) Γ( α)t α = t u ()d Γ( α) 0 (t ) α, (3) is a regularized derivative of order α, 0 < α <. Furthermore, D α 0tu is the Riemann Liouville fractional derivative of order α, 0 < α <. It is to be noted that (Da α +u)(t) = ( d ) t n u d, n < α < n, (4) Γ(n α) dt 0 (t ) α n+ and (D α a +()β ) = Γ(β) (β α) tβ α, (5) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 97 and D α α +(t) = Γ() Γ( α)t α = u(0), if β =. (6) Γ( α)tα In (2), the function f(u, t) is defined on a rectangle D = {0 t T, u u 0 U} and satisfies the Lipschitz condition. Also, on 0 t T, we introduce the grid ω = {t j = j, j = 0,..., j 0 }. Denote a grid function by y j = y(t j ). In calculus, in the study of differential equations, the Picard- Lindelof theorem, Picard s existence theorem or Cauchy-Lipschitz theorem is an fundemental theorem on existence and uniqueness of solutions with given initial and boundary conditions. 2.2 Picard-Lindelof theorem Theorem. Consider the following initial value problem y (t) = f(t, y(t)), y(t 0 ) = y 0, t [t 0 ɛ, t 0 + ɛ]. (7) If f be a uniformly Lipschitz continuous in y and also ontinuous in t. Then, for some value ɛ > 0, there exists a unique solution y(t) to the initial value problem on the interval [t 0 ɛ, t 0 + ɛ]. Theorem 2. Let U R n be an open and connected set. Also, assume (0, b) R + and define D = (0, b) U. Furthermore, supposed that f(x, y) be a real valued function on D. Now consider the following equation: under the conditions,. f is continuous in D, ( D is closure of D), 2. f satisfies a Lipschitz condition in D, y α = f(x, y), 0 < α, (8) then (x 0, y 0 ) D, a positive number β can be found such that the closed interval I = [x 0 β, x 0 + β] is contained in (0, b) and there exists a unique continuous function y : I U such that y α = f(x, y), x I, y(x 0 ) = y 0. (9) Proposition. Theorem 2, implies that the Cauchy problem (2), has a unique solution. 3 Discretization process of differential equations of fractional order Hereunder, we find a discrete methodology of the fractional derivative of order α. It is assumed that the solutions have the required smoothness. Thus, Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 98 we have = α 0tu t=tj = tj Γ( α) 0 ((j s + )) α ((j s)) α u(η)dη (t j η) α u t,s + O() 2, (0) ξ is an intermediate point between η and t s and ü = 2 u u, u = 2 t2 t and u t,s = u(t s+) u(t s ) and big O notation describes the limiting behavior of a function when the argument tends towards a particular value or infinity. If t j = t 0 + j, t s = (s ), t s = s, and u i = u(t i ), then by using Taylor series expansion we will have u(t s ) = u(t s ) + u(t 2 s )(t s t 2 s ) +..., 2 u(t s+ ) = u(t s ) + u(t 2 s+ )(t s t 2 s ) +..., () 2 u(t s+ ) u(t s ) = u(t s+ ) + (t s t 2 s ) + O( 2 ). 2 We say that the expression s= u(t s ) = u s+ u s = u(t s+) u(t s ) + O(), (2) 2 approximates the first derivative u = du. Now, to find an upper bound on dt the second sum in expression (0) one will set M Γ( α) s ts Γ( α) Assume that ü M and Consequently, we have s= t s (t j η) α dη α 0t j u = ts 2 ) ü(ξ)(η t s t s (t j η) α dη M Γ( α) (t α j s+ t α j s s ) + O(). (3) D α 0t j u = α 0t j u + O(). (4) (t α j s+ t α s= j s )u t,s. (5) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 99 which is a discrete methodology of the fractional derivative. Hence (t α j s+ t α s= y(0) = y 0 = u 0, y t,s = (y s+ y s ) Now, for the error z = y u, we have (t α j s+ t α s= j s )y t,s = f(y j, t j ), (6) = y(t s+) y(t s ). (7) j s )z t,s = f u z j + Ψ j, (8) z(0) = z 0 = 0, Ψ j = O(). (9) It is to be noted that f u denotes the derivative at ξ and u j ξ u j + z j. Let y = z + u, consequently Let A = (t α j s+ t α j s In the other words ((t α j s+ t α s= s= ), hence A(z t,s ) + j s )(z t,s + u t,s)) = f(z j + u j, t j ), (20) Au t,s = f(z j + u j, t j ). (2) s= A(z t,s )+f(z j +u j, t j ) s= Consequently and As we know Au t,s = f(z j +u j, t j ). (22) s= D α 0t j u = α 0t j u + O(), (23) Az t = f(z j + u j, t j ) f(u j, t j ) + O(). (24) s= ξ (u j, u j + z j ) s.t f(u j + z j, t j ) f(u j, t j ) = f u (ξ, t j)z j = f u z j. (25) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 200 Consequently s= Now, consider the following equation Az t,s = f u z j + O(). (26) (t α j s+ t α s=0 j s )z s+ z s = f (ξ, t)z j + Ψ j, (27) z t,s = z s+ z s. (28) As a direct consequence of (27) we obtain [ ] z j+ = ( (t α 2 t α )) α + α f (ξ, t) z j (29) Assume that + α [ (t α j+ t α j z 0 + ( t α j+ + 2t α j t α j )z +... ] +( t α 3 2t α 2 t α z j ) + α Ψ j. x, Let us first show the validity of the inequality or Define the function 2x α (x + ) α + (x ) α 2α. (30) A j = t α j+ 2t α j + t α j < 0, j, (3) 2j α > (j + ) α + (j ) α, j. (32) f(t) = 2t α (t + ) α, t. (33) + (t ) α Since we have f() = 2 α, lim f(t) =, f (t) < 0, (34) t < f(t) 2 α, x, t [, ), t α j+ + 2t α j t α j > 0, t. (35) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 20 Assume that f(ξ, t) < a when a < 0. Consequently (29) implies the estimate: or z j+ (2 α ) α a max 0 s j z s (36) +(2 α ) max 0 s j z s + α Ψ j, z j+ α a max 0 s j z s + α Ψ j. (37) Whence, (37) implies that z j+ j =0 α max Ψ 0 s j j. (38) In the other words, expression (38) implies that scheme (6), (7) converges at a rate of O( α ). 4 Navier - Stokes flow with measures as initial vorticity As we know, () can be written in the following form [0] D α 0tu(x, t) = q(x, t)u(x, t) k(x, t)u(x, t) + f(x, t), (39) u(0, t) = u(l, t) = 0, D0t α u t=0 = u 0 (x), k(x, t) c > 0, q(x, t) 0. Now, in the cylinder Q T = G [0 < t T ], consider the following problem α 0tu = Lu + f(x, t), (x, t) Q T, (40) u(x, t) Γ = µ(x, t), t 0, u(x, 0) = u 0 (x), x Ḡ, Lu = L k u, L k u = 2 u x 2, k =, 2, p. k k= Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 202 Furthermore, 0tu α = t u (x,) Γ( α) 0 (t ) d, 0 < α <, is the Riemann Liouville α fractional derivative of order α. Introduce a grid w h that is uniform in each direction x k : w h = {x i = (i h,, i p h p ) G, i k = 0,,, N k, h k = l k }, (4) N k γ h is boundary nodes. w = {t j = j, j =, 2,, j 0 }, Problem (40) is approximated by the difference scheme: Hence one will set: Λy = α 0ty = Λy (σ) + ϕ, x w h, t w, (42) y γh = µ(x, t), x γ h, t w, (43) y(x, 0) = u 0 (x), x w h, u xx = Λ. (44) (Λ k )y, Λ k = y xk x k, k =, 2,, p. (45) k= Now, we get a one-parameter family of difference schemes y xx = y xx,i = (y i+ 2y i y i+ )/h 2, (46) y (σ) = σŷ + ( σ)y, 0 σ. (47) Sometimes, scheme (47) will be treated as a scheme with weights. We denote by y j i the value at the node (x i, t j ) of the grid function y given on w j. Now, consider the following equation ŷ = y j+, y = y j ϕ = f(x, t ), t = t j+ 2, (48) α 0t j y = (t α j s+ t α s=0 j s )ys t. (49) we define sum y = ν + ν and then estimate the solution of difference problem (42), (43), ν is the solution to problem (42),(43) with ϕ = 0 and ν is the solution to the same problem with µ = 0 and u 0 (x) = 0. First, we rewire ν to the canonical form (see [3, 4]) ( α + σ h 2 k= k ) ν j+ i k = σ h 2 k= k ( ν j+ i k + + ν j+ i k )+ (50) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 203 Since ( 2 2 α α k= ( t α + and B(P, Q) > 0 if α < 2( σ) h 2 k j+ t α j ) ν j i k + k= ν 0 i k + t α 2( σ) h 2 ( ν j i k + + ν j i k ) k j+ + 2t α j t α 3 + 2t α 2 t α A(p) = ν j i k ). α + 2σ h 2 k 2 2 α ( 2( σ) k= k= t α j ν i k + + > 0, (5) ), 0 < σ <, (52) h 2 k t α j+ + 2t α j t α j > 0, j, (53) by using maximum principle we have (see [3]) ν 0 C µ C, (54) and µ = { µ, x γ h, t w, u 0 (x), x w h, t = 0, (55) µ C = max x γ h, t w h µ(x, t). (56) To estimate ν rewriting problem for ν ( α + Φ(P (j+) ) = k= ) 2σ h 2 k ν j+ i k = ( 2 2 α α ( t α + k= j+ t α j σ h 2 k= k ) 2( σ) h 2 k ν 0 i k + t α ( ν j+ i k + + ν j+ i k ) + Φ(P (j+) ), (57) ν ik + k= j+ + 2t α j σ h 2 ( ν ik ++ ν ik ) (58) k t α j ν i k + Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 204 Consequently Φ(P (j+) ) = Q S j + t α 3 + 2t α 2 t α ν j i k )+ϕ j i k. B(P, Q) ν(q) + F (P (j+) ), P (j+) = P (x, t j+ ). (59) Thus, canonical form (57) can be rewritten to yield the equation A(P (j+) ) ν(p (j+) ) = B(P, Q) ν(q) + φ(p (j+) ), (60) Q S j+ S j+ is the set of nodes Q(ξ, t j+) S (P (x, t j+ )) and S j of nodes Q(ξ, t j ),, Q(ξ, t ) S (P (x, t j+ ). Introduce the notation Note that for all Q S j D (P (j+) ) Consequently D (P (j+) ) = A(P (j+) ) + Q S j+ is the set B(P (j+), Q). (6) D (P (j+) ) = α, A(P (j+)) > 0, (62) B(P (j+), Q) = and Q S j+ if Q S j α < σ h 2 k= k 2 2 α ( 2Γ( α)( σ) k= > 0, 0 < σ <, (63) ), (64) h 2 k B(P (j+), Q) = α ν j+ Ch y 0 Ch + Combining (57) and (66) gives j =0 Q S B(P (j+), Q) =. (65) α y j+ Ch ν 0 Ch + max 0<k j+ µ(t k) + max 0 s j ϕs Ch. (66) j =0 α max 0 s j ϕs Ch. (67) Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 205 The error z = y u satisfies the estimate z j+ Ch j =0 α max 0 s j ψs Ch. (68) Since ψ = O( h 2 + ɛ), h 2 = h 2 + h 2 p, it follows from (68) that For α, (69) implies the well-known result z j+ Ch = O( α + h2 ). (69) α z j+ Ch = O( h 2 + ). (70) References [] D. Baleanu, K. Diethelm, E. Scalas, J. J. Trujillo, Fractional calculus models and numerical methods, World Scientific, Singapore, 202. [2] K. Sayevand, K. Pichaghchi, Successive approximation: A survey on stable manifold of fractional differential systems, Fract. Calculus Appl. Anal., 8 (3), 2554-2556 (205). [3] A. A. Samarskii, A. V. Gulin, Stability of difference schemes, Nauka, Moscow, 973. [4] K. Shkhanukov, On the convergence of difference schemes for differential equations with a fractional derivative, Dokl. Akad. Nauk, 348, 2554-2556 (996). [5] C. Tadjeran, M. Meerschaert, H.P. Scheffler, A second-order accurate numerical approximation for the fractional diffusion equation, J. Comput. Phys., 23 (), 205-23 (2006). [6] M. Marin, An evolutionary equation in thermoelasticity of dipolar bodies J. Math.l Phys., A.I.P., 40(3), 39-399 (999). [7] M. Marin, R. Agarwal, S. Mahmoud, Nonsimple material problems addressed by the Lagrange s identity, Bound. Value Prob., 203: doi: 0.86/687-2770-203-3 Download Date 3/6/8 5:03 PM

A NOVEL DIFFERENCE SCHEME... 206 [8] S. Momani, Z. Odibat, Analytical solution of a time-fractional Navier - Stokes equation by Adomian decomposition method, Appl. Math. Comput., 77, 488-494 (2006). [9] M. El-Shahed, A. Salem, On the generalized Navier - Stokes equations, Appl. Math. Comput., 56, 287-293 (2004). [0] H. Jafari, H. Tajadodi, D. Baleanu, A modified variational iteration method for solving fractional Riccati differential equation by Adomian polynomials, Frac. Cal. Appli. Anal., 6 (), 09-22 (203). Khosro SAYEVAND, Department of Mathematics, Malayer University Malayer, Iran Email: ksayehvand@malayeru.ac.ir Dumitru BALEANU, Çankaya University, Faculty of Art and Sciences, Department of Mathematics, Ankara, Turkey and Institute of Space Sciences, MG-23, R 76900, Magurele-Bucharest, Romania Email: dumitru@cankaya.edu.tr Fatemeh SAHSAVAND, Department of Mathematics, Malayer University Malayer, Iran Email: f.sahsavand@gmail.com Download Date 3/6/8 5:03 PM