A fractional nonlinear system for release oscillation of silver ions from hollow fibers

Introduction

A nonlinear oscillator can be generally expressed as ¹

y ″ + y + N ( y , y ′ , y ″ ) = 0

(1)

N ( y , x , x ′ ) y > 0

(2)

The nonlinear oscillator given in equation (1) is equivalent to the following nonlinear system

{ y ′ = x x ′ = − y − N ( u , x , x ′ )

(3)

Analytical methods for nonlinear oscillators include the homotopy perturbation method,^1–9 the variational iteration method,^1,10,11 frequency formulation^12–16 and others.

In this paper we will consider the following generalized nonlinear system

{ x ′ = a x − b x y y ′ = c x y − d y

(4)

Equation (4) is the well-known Lotka–Volterra predator–prey model, and it can also describe a damped oscillation of a bridge or release of silver ions from hollow fibers.^17,18 As an example, the well-known van der Pol oscillator can be expressed as

{ x ′ = μ ( x − 1 3 x 3 − y ) y ′ = 1 μ x

(5)

{ x ′ = y y ′ = μ ( 1 − x 2 ) y − x

(6)

Fractional modification of the model

In this paper we study the fractional modification of equation (4)

{ d α x d t α = a x − b x y d α x d t α = c x y − d y

(7)

d α x d t α = 1 Γ ( n − α ) d n d t n ∫ t 0 t ( s − x ) n − α − 1 [ x 0 ( s ) − x ( s ) ] d s

(8)

The physical meaning of the fractional derivative was discussed in literature, ^{3 , 6 , 19} and it can be widely applied to discontinuous problems.^20–23

All predator–prey models assume that time is continuous; however, the number of either prey or predator changes not smoothly with time. The discontinuous growth rate is of their intrinsic property, and a fractional model can best model its property of periodic property.

The release of ions from hollow fibers can be also described by the fractional calculus, where the water can be considered as a porous medium.

By the fractional transform ⁶

s = t a Γ ( 1 + a )

(9)

Equation (6) can be converted into the following one

{ d x d s = f ( x , y ) = a x − b x y d x d s = g ( x , y ) = c x y − d y

(10)

The equilibrium point locates at

{ x 0 = d c y 0 = a b

(11)

The equilibrium point for van der Pol oscillator is (0,0). In this paper we search for the periodicity near the equilibrium point, so the system, equation (10), can be linearized as

{ f ( x , y ) = f ( x 0 , y 0 ) + ( x − x 0 ) ∂ f ∂ x ( x 0 , y 0 ) + ( y − x 0 ) ∂ f ∂ y ( x 0 , y 0 ) g ( x , y ) = g ( x 0 , y 0 ) + ( x − x 0 ) ∂ g ∂ x ( x 0 , y 0 ) + ( y − x 0 ) ∂ g ∂ y ( x 0 , y 0 )

(12)

Equation (10) becomes

{ x ′ = a x − b x y ≈ a x 0 − b x 0 y 0 + ( a − b y 0 ) ( x − x 0 ) − b x 0 ( y − y 0 ) y ′ = c x y − d y ≈ c x 0 y 0 − d y 0 + c y 0 ( x − x 0 ) − ( c x 0 − d ) ( y − y 0 )

(13)

We assume that the system has the periodic solution

{ x = A + B sin ⁡ ( ω s + α ) y = m + n sin ⁡ ( ω s + β )

(14)

In view of equation (14), we write equation (13) in the form

{ B ω cos ⁡ ( ω s + α ) = x 0 − b x 0 y 0 + ( a − b y 0 ) ( A + B sin ⁡ ( ω s + α ) − x 0 ) − b x 0 ( m + n sin ⁡ ( ω s + β ) − y 0 ) n ω cos ⁡ ( ω s + β ) = c x 0 y 0 − d y 0 + c y 0 ( A + B sin ⁡ ( ω s + α ) − x 0 ) − ( c x 0 − d ) ( m + n sin ⁡ ( ω s + β ) − y 0 )

(15)

In order to guarantee a periodic solution, the following criterion has to be satisfied

{ cos ⁡ ( ω s + α ) = sin ⁡ ( ω s + β ) cos ⁡ ( ω s + β ) = sin ⁡ ( ω s + α ) a − b y 0 = 0 c x 0 − d = 0 1 − b m = 0 B ω = − b x 0 n − d + c A = 0 n ω = c y 0 B

(16)

Equation (16) implies that

{ β = π 2 − α d c = a b m = 1 b A = d c a B 2 + b n 2 = 0

(17)

The frequency is obtained as

ω = − b d n c B

(18)

As an example, considering the case a = d = −1, b = c = 1, α = 0 , β = π / 2 , m = 1, A = −1, B = 1, n = 1, and ω = 1 , we obtain the solution

{ x = − 1 + sin ⁡ ( s ) y = 1 + sin ⁡ ( s + π 2 )

(19)

or

{ x = − 1 + sin ⁡ ( t a Γ ( 1 + a ) ) y = 1 + sin ⁡ ( t a Γ ( 1 + a ) + π 2 )

(20)

Figure 1 shows that the periodicity criterion is satisfied, the frequency is not changed for α  = 1, and the frequency increases with time for α  < 1.

OPEN IN VIEWER

Figure 1.

Frequency change with time.

Conclusions

This short paper elucidates the periodic property of a fractional nonlinear system; the predator–prey model is used to illustrate the solution process. As the predator–prey model can describe the release of silver ions from hollow fibers or nanofiber membrane,^17,18 Figure 1 can be explained as the release oscillation, which is greatly affected by the fractional order. Generally the frequency changes with time, and when time tends to infinity, the frequency becomes infinite large.