Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Journal of Non-Equilibrium Thermodynamics

Founded by Keller, Jürgen U.

Editor-in-Chief: Hoffmann, Karl Heinz

Managing Editor: Prehl, Janett / Schwalbe, Karsten

Ed. by Michaelides, Efstathios E. / Rubi, J. Miguel

4 Issues per year

IMPACT FACTOR 2017: 1.633
5-year IMPACT FACTOR: 1.642

CiteScore 2017: 1.70

SCImago Journal Rank (SJR) 2017: 0.591
Source Normalized Impact per Paper (SNIP) 2017: 1.160

See all formats and pricing
More options …
Volume 43, Issue 1


Fractional Effects on the Light Scattering Properties of a Simple Binary Mixture

Rosalío F. Rodríguez
  • Corresponding author
  • Instituto de Física, Universidad Nacional Autónoma de México, Apdo. Postal 20-364, 01000 CDMX, México
  • Fellow of SNI, México. Also at FENOMEC, UNAM, México
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Jorge Fujioka
  • Instituto de Física, Universidad Nacional Autónoma de México, Apdo. Postal 20-364, 01000 CDMX, México
  • Fellow of SNI, México. Also at FENOMEC, UNAM, México
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Elizabeth Salinas-Rodríguez
  • Departamento I. P. H., Universidad Autónoma Metropolitana, Iztapalapa, Apdo. Postal 55-534, 09340 CDMX, México; Fellow of SNI, México
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2017-12-14 | DOI: https://doi.org/10.1515/jnet-2017-0036


A fractional generalized hydrodynamic (GH) model of the concentration fluctuations correlation function is analyzed. Our analysis is based on a previously proposed irreversible thermodynamics non-fractional model that describes the first GH deviations in a finite frequency and wavelength regime where the local equilibrium assumption is not valid. Our basic purpose is to generalize this theory to investigate the time fractional effects on the intermediate scattering function (ISF) of a model of a binary mixture. We discuss two different forms of introducing the fractional derivatives into the hydrodynamic time evolution equation of the ISF and examine their consistency in the non-fractional limit. We calculate analytically the fractional intermediate scattering function (FISF) and compare it with the non-fractional one. We find that although the FISFs are in general less than ISF, the FISF may be a significant fraction of the ISF (∼36%) and might be measurable by light or neutron scattering techniques. We show that fractional time derivatives provide a consistent description of this correlation function in the GH regime and reduce to its well-known behavior in the Navier–Stokes limit (NS) where local equilibrium is restored. We also suggest an experimental verification of our results for near equilibrium states. Finally, we summarize the main results of our work and make some further physical remarks.

Keywords: fractional derivatives; fluctuations; generalized hydrodynamics; irreversible thermodynamics

1 Introduction

A successful account of many irreversible processes has been provided over the years by linear (classical) irreversible thermodynamics (LIT), a theory built upon the local equilibrium hypothesis and restricted to the linear regime [1, 2, 3, 4, 5, 6]. However, to describe more general processes where the local equilibrium assumption if no longer valid, a generalization of LIT is called for. This situation arises in the generalized hydrodynamic (GH) regime of intermediate wave vectors k and frequencies ω, where local equilibrium breaks down [7]. These states may arise for systems like fractals or self-organizing structures where the microscopic conditions for the molecular chaos assumption may not be met and, therefore, the local equilibrium hypothesis is no longer valid [8, 9].

Among the several phenomenological approaches used for this purpose is extended irreversible thermodynamics (EIT). The central idea in EIT is to enlarge the space of state variables by assuming that the entropy depends not only on the conserved hydrodynamic densities but also on the dissipative fluxes or non-conserved variables as well. In this process the assumption of local equilibrium breaks down, but when the subspace spanned by the non-conserved variables reduces to that of the conserved ones, the local equilibrium assumption is restored and EIT reduces to the well-known LIT [10, 11, 12, 13]. In this way EIT provides for relaxation type equations which play the role of constitutive equations for the fast variables of the system in terms of the slow ones, and associates to each of them a specific time scale. In this way a set of characteristic relaxation times is established for nonequilibrium states beyond local equilibrium.

Distinctive features of the nature of the spectrum of characteristic relaxation times of nonequilibrium states may be such that a well-defined separation between the time scales associated with the macroscopic transport processes and the microscopic ones giving rise to them, no longer exist. Under these conditions the fluctuations of the state variables may exhibit strange dynamics where there is no preferred time scale and the central limit theorem (CLT) may not be applicable [14, 15, 16, 17]. Then the usual smoothing process of statistical mechanics is no longer feasible to describe the dynamics of fluctuations and the usual formalism based on ordinary calculus would no longer be adequate to describe the transport processes. One way to consider these effects of non-localizability in time (memory) and non-localizability in spatial coordinates (spatial correlations) is to develop new methodology based on fractional calculus. This leads to a slow relaxation of the correlation functions associated with many transport properties, such as the concentration fluctuations correlation function producing the scattered spectrum of light from complex fluids [18, 19]. In this sense, the violation of the principle of local equilibrium makes the traditional methods in transport theory to become unsuitable. Examples of macroscopic systems which may exhibit these features are the large variety of relaxation processes that occur in viscoelastic fluids, glassy materials, synthetic polymers or biopolymers, all of which have in common that their relaxation functions are non-exponential, due to the large number of highly coupled elementary units responsible for the relaxation. It is this requirement of high cooperation between these elements that calls for a fractional description [20, 21].

For these reasons, it appears worthwhile to investigate these issues on a simple, but not unrealistic model, which lends itself to explicit analytic mathematical manipulation. In previous work we have developed such a model using EIT [12]. In [12] we discuss how the behavior of the concentration correlation function is modified for nonequilibrium states in the GH regime, where additional time scales influence the dynamics of state variables fluctuations. In the present work, we show that the inclusion of fractional time derivatives in the transport equations provides a consistent description to include these scales and examine their effect on time correlation functions. For this purpose, we propose two different novel methods to introduce fractional time derivatives into the macroscopic equations for the fluctuations of a simple model for an isothermal-isobaric, heat insulating, inviscid and chemically inert binary mixture. Although at first sight this model system might seem to be somewhat restrictive, mixtures of this sort arise in real permeation systems where, for instance, an inert binary mixture of He-N₂ is produced in different regions of a tubular separator involving porous membranes used for gas–gas separation processes [22, 23] and may comply with the requirements of the present model. This type of mixtures also appears in other systems of engineering and physico-chemical interest like the steady-state mass transfer of non-reactive species [24].

More specifically, we calculate analytically the fractional intermediate scattering function FISF in the GH regime and show how when the fast variables relax and local equilibrium is restored, this correlation function recovers its NS behavior. We find that although the FISFs are in general less than the ISF, the FISF may be a significant fraction (∼36 %) of ISF and might be measurable by light or neutron scattering techniques. This model calculations show that fractional time derivatives provide a consistent description of this correlation function in the GH regime and reduces to its well-known behavior in the NS limit. Although we are not aware of the existence of light or neutron scattering experimental results for these systems, we suggest light scattering experiments which could shed some light on the confrontation of our theoretical model predictions and experimental measurements. This even partial verification could hopefully serve to test the nature of these processes in a fractional and non-classical hydrodynamic regime [7].

2 Model

Since the formal theory behind our hydrodynamic model has been developed in our previous work [12] using EIT, in this section we shall only briefly review the underlying theoretical ideas which lead to the final dynamic equations that we shall use later. The state variables for the binary mixture were chosen as one conserved variable, the relative concentration of one of the components c(r,τ), and two non-conserved variables, the local mass flux J(r,τ) and a second-order tensor (r,τ), which accounts for the spatial inhomogeneities in the space of state variables. The mass conservation equation reads dcdτ=J(r,τ)(1)

where d/ denotes the hydrodynamic derivative; the equation of motion is given by ρdudτ=p(2)

where ρ is the local mass density, p is the pressure and u is the macroscopic flow velocity. If as explained in [12], we apply the usual procedure of EIT, restricting the calculations to lowest order in the non-conserved variables, we end up with the following coupled set of relaxation equations for the non-conserved variables (J,) (see eqs. (2.8)–(2.9) in [12])


where μ10, l1 and l2 are phenomenological coefficients. The relaxation times τJ and τ of J and are defined explicitly in [12] in terms of phenomenological coefficients (see eqs. (2.8)–(2.11) in [12]) which can only be determined either from experiment or from a microscopic theory. On intuitive grounds one would be tempted to expect that will vary on a shorter time scale than that of J, so we restrict our description to a time scale in which has already relaxed and τ is negligible so that eq. (3) becomes (see eq. (2.14) in [12]) τJdJdτ=J+Dcl1l22J.(5)

Since in most experimental situations involving mass diffusion the pressure is kept constant, D denotes the isothermal-isobaric diffusion coefficient D=(Tμ10)1(μ/c)T,p, where μ is the chemical potential. Note that although this relaxing equation (generalized constitutive equation) is valid in the GH limit where local equilibrium is lost, in the NS regime with τJ=l1l2=0, it reduces to Fick’s law restoring the validity of local equilibrium.

Equations (1)–(5) constitute our basic closed set of hydrodynamic equations and to derive transport properties from them, they should be linearized in the small deviations with respect to an initial reference state identified with the subscript 0, δc(r,t)=c(r,t)c0(r), δJ(r,t)=J(r,t)J0(r). Since these deviations are random quantities, they will manifest themselves as the hydrodynamic fluctuations. The relevant correlation function to calculate the light scattering properties of the mixture is the equilibrium correlation function of concentration fluctuations [25]


where the angular brackets denote an average over an equilibrium ensemble. The intensity distribution of the (isotropic) component of the Rayleigh–Brillouin scattering spectrum of a molecular liquid, i. e., the intermediate scattering function I(k,τ) (ISF) or van Hove’s self-correlation function [7], is proportional to the spatial Fourier transform of δc(r,τ) within the scattering volume I(k,τ)=δc(k,τ)δc(k,0)〉;,(7)

where the asterisk (*) indicates complex conjugate. Substitution of eq. (1) into eq. (5) yields the time evolution equation for I(k,τ)


where A1(k)(1+l1l2k2)/τJ and A2(k)Dk2/τJ are two free parameters which are combinations of phenomenological coefficients τJ, l1, l2. For given I(k,0) with vanishing time derivative I(k,0)=0 its solution reads I(k,τ)I(k,0)=exp(12A1τ)(cosRτ+A12RsinRτ)(9)

with R(k)=[A2(k)A12(k)/4]1/2>0, which imposes the condition A2(k)>A12(k)/4. It must be emphasized once again that the inclusion of spatial inhomogeneities through the non-conserved variable directly took us out of the NS regime and drove us into the domain of GH. However, the local equilibrium regime is easily regained by taking the limit τ0, l1l20, in eq. (9) which yields the classical exponential decay of I(k,τ)


Furthermore, I(k,τ)must satisfy certain restrictions, the so-called sum rules [7], which yield a relation between the static structure factor S(k) and A2(k), namely (see eq. (3.11) in [12]) S(k)=kBTA2(k)k2=kBTDτJ.(11)

Since S(k) may be determined from a microscopic theory in terms of the measurable quantities T and D at least in principle, this equation allows us to determine τJ. This would leave us with one free parameter A1(k), or equivalently l1l2. Usually, phenomenological approaches like EIT introduce a wealth of phenomenological coefficients in the equations of motion or the transport equations, which can only be determined either experimentally or from microscopic theories [10]. However, in the present model neither the values of A1 or A2 are known for the binary mixture, but since the condition A2(k)>A12(k)/4 must be fulfilled, it is convenient to reduce the number of phenomenological coefficients in eq. (8) by rescaling the equation through the following change of variables IaJ, τbt. These yields ab2A2Jtt+aA1bA2Jt+aJ=0(12)

where JtJ(k,t)/t, Jtt2J(k,t)/t2. If we now choose a and b in such a way that the coefficients of both derivatives equal unity, then a=A2/A12,b=A11 and eq. (12) is rewritten in the rescaled simpler form Jtt+Jt+aJ=0(13)

with only one phenomenological coefficient a. For given J(k,0), Jt(k,0)=0 its solution becomes Jk,tJk,0=et/2cosMt+12MsinMt,(14)

where M=a1/41/2. By taking the Laplace transform of Jk,t with respect to time and setting s=iω gives the structure factor


where we have defined Sk=Jk,0 and used the abbreviations


3 Time fractional derivatives

Physical models where first-order conventional time derivatives have been formally replaced by fractional differential operators have been proposed for a variety of physical systems [17, 26]. It should be remarked, however, that there are several different, but not necessarily equivalent, ways of defining the time fractional derivative (TFD) of a function f(t). Among these, the definitions of Riemann–Liouville [27], Grünwald–Letnikov [28], Caputo [29, 30], Hadamard [31], Erdelyi [32] or Ortigueira [33] can be found in the literature; however, despite their differences, all the TFDs share the common feature of being non-local operators. In other words, the TFD of f(t) at a point to not only depends on the behavior of f(t) in the neighborhood of t0, but also on the values of f(t) over finite intervals of the form t1tt0 or t0tt2. Therefore, each of the above possible definitions of a FTD (except for Ortigueira’s centered derivative) provides an expression for the left-hand sided derivative, and a different expression for the right-handed derivative.

The selection of the appropriate FTD to be used in our analysis is a delicate matter which should be performed by considering criteria consistent with the physical situation at hand; we have used the following criteria: First, the FTDs should be left-handed because they can only be calculated in terms of their past values. Secondly, the solution of a fractional differential equation requires to know appropriate initial conditions to obtain a particular solution of the fractional version of eq. (13). We choose the left-handed Caputo time fractional derivative (LHCD) because then this equation can be solved by Laplace transformation, which incorporates the natural initial conditions, i. e., the integer order derivatives at t=0. Since these are the initial conditions that can be specified or controlled in the relaxation processes described in the present model, it is consistent to use a LHCD in eq. (13). In contrast, the Laplace transform of the Riemann–Liouville derivative demands to know fractional derivatives at t=0, which are unknown for the usual type of relaxation processes toward equilibrium in dissipative systems.

For an arbitrary function f(k,t) the LHCD is defined by [28, 29], 0CDtαfk,t=1Γnα0tf(n)k,τtτα+1ndτ(17)

Here f(n)(k,t)dnf/dtn is the conventional n order time-derivative, being n the smallest integer greater than α, i. e., n-1 ≤ α<n, and where Γ(nα) is the Euler gamma function. An intuitive picture of the action of 0CDtαf(k,t) implies that the LHCD considers the past behavior of f(n) by adding to its initial values, the successive weighted increments over time. These increments per unit time are represented by f(n)(k,t) in eq. (17), while the weights are indicated by the factors 1/(tτ)α+1n, whose value decreases with increasing time separation from time t.

To introduce TFD in eq. (8) we consider two different approaches in the two following subsections. The first one (Case A) is motivated by the description of the dispersion of ultra-short pulses in fiber optics and consists in replacing the first and second derivatives which appear in eq. (13) by a fractional derivative of order 1<α<2 [34]. In our second approach (Case B) instead of the rather abrupt replacement of the first and second integer derivatives by a unique fractional derivative, we shall gradually approximate both integer derivatives toward an intermediate value. That is, we shall replace the first (integer) order derivative by a fractional derivative of order 1+ε and the second one by a fractional derivative of order 2ε. Let us now examine in detail these two generalizations of eq. (8).

3.1 Case A: One centered fractional derivative

According to the first approach, we replace the two integer derivatives in eq. (13) by a single LHCD of order ε, 1<ε<2. Since in our numerical tests we shall use ε=1.5, eq. (13) becomes D1.5JFA+aJFA=0,(18)

whose solution is


where E1.5 is the one-parameter Mittag–Leffler function. To simplify the notation here and in what follows we shall not write explicitly the vector character of k.

3.2 Case B: Two fractional derivatives

In our second approach, we replace the first and second derivatives by fractional derivatives of order 1+ε and 2ε, respectively, thus obtaining an equation of the form D2εJF+D1+εJF+aJF=0.(20)

It should be noted that ε could be positive or negative. However, we will only consider the case ε>0 so that eq. (20) approaches eq. (18) when ε0.5. In this equation the parameter ε only affects the first two terms of the equation, but has no effect whatsoever on the third term. This feature is unsatisfactory, but it can be improved by first introducing the fractional derivatives into eq. (8) and afterwards performing the scaling. Following this procedure, we obtain a fractional equation of the form c2ε,kD2εJFB+c1ε,kD1+εJFB+A2kJFB=0(21)

where c1ε,k=A1k(1ε)+ε,(22) c2ε,k=1ε+εA1k.(23)

With this choice of c1 and c2, when ε=1/2 the two fractional derivatives in eq. (21) coalesce into a single derivative of order 1.5, as in eq. (18). If now we perform the change of variable JFB=aˆJFB in eq. (21) and if we also make the change of variables T=bˆt in eq. (13), we arrive at 0CDtαJFBk,T=1Γnα0tbˆnα1tYnα11bˆndnJFBdYnbˆdY,(24)

where n1α<n. If we now set y=bˆdY, eq. (24) takes de form 0CDtαJFBk,T=bˆαΓnα0ttYnα1dnJFBdYndY=1bˆα0CDtαJFBk,T.(25)

On the other hand, by setting JF=aˆJFB in eq. (21) and then by introducing the change T=bˆt in the resulting equation, we obtain


To determine the values of aˆ and bˆ we require that the coefficients of the two derivatives in this equation be equal to unity with the result




Finally, substituting eqs. (22) and (23) in these equations yields 0Dt2εJFBt+0Dt1+εJFBt+aˆJFBt=0,(29)

where now aˆ depends on A1, A2 and ε,


instead of a=A2/A12 as used before regarding eq. (12).

The solution of eq. (29) for 0<ε<0.5 is obtained as follows. We first take the Laplace transform of eq. (29) with respect to t, recalling that the Laplace transform of the Caputo’s derivative is given by L0Dtαf=1sδLf(n),(31)

where n1α<n and δ=nα. In this way, we arrive at the Laplace transform of the solution of eq. (29), JˆFBk,sJ0FB=Uˆs+Vˆs,(32)





By inserting these expressions into eq. (32) and by taking the inverse Laplace transform of the resulting equation we obtain JˆFBk,tJ0FB=t12εm=0t12εmm!E2ε,2mm+2(m)aˆt2ε+m=0t12εmm!E2ε,1m1+ε(m)aˆ1εt2ε,(35)

Here Eα,βmz denote the derivatives of the generalized Mittag–Leffler function (also known as Erdely’s function) which for negative and large values of z have the following asymptotic behavior [28] Eα,βmzm=01m+1m!zm+1Γαβ.(36)

For future use we use this expression to obtain the following asymptotic (long-time) behavior of JˆFBk,t/J0FB [19], JˆFBk,tJ0FB=m=01tm2ε2α1Γ1t2ε+1Γ2,(37)

where we have defined Γ1=Γ2εm+ε1/2,(38) Γ2=Γ2εm+1.(39)

Taking the leading term of the sum in eq. (37) yields the power law behavior


4 Results

For the purpose of having a reference case to compare with, let us first examine the behavior of the non-fractional normalized correlation function Mk,t;A2Jk,t/Jk,0 defined by eq. (14) as a function of the dimensionless time t, for different values of A2(k)Dk2/τJ and setting A1=1. Figure 1 shows its damped oscillating behavior. For all the considered values A2=5, 10, 20, the correlation Mt;A2 shows the dominant oscillatory damped exponential decay leading to the exponential decay behavior of the NS regime given by eq. (10). Small values of A2 correspond to large relaxation time τJ, indicating that the fast variable J has not yet relaxed and prevents the validity of local equilibrium. Thus, A2 controls the degree of departure of the state from the NS regime toward the GH regime. However, as A2 decreases the amplitudes of the GH correlations also decrease, but their range (wavelength) increases.

Plot of the non-fractional correlation Mt;A2=Jt;A2/J0;A2$\mathcal {M}\left( {t; {A_2}} \right) = \mathcal {J}\left( {t; {A_2}} \right)/\mathcal {J}\left( {0;{A_2}} \right)$, as given by eq. (17), for A1=1 and A2=5 (_ _), 10( ….), 20 (____).
Figure 1:

Plot of the non-fractional correlation Mt;A2=Jt;A2/J0;A2, as given by eq. (17), for A1=1 and A2=5 (_ _), 10( ….), 20 (____).

Fractional effects are considered in Figure 2 which shows the time behavior of the fractional correlation Nt;A2JFAt;A2/JFA0;A2 for A2=1, 5, 10, 20, as calculated from eq. (19) for Case A. The curves do not oscillate in the same time interval considered in Figure 1, and their ranges and amplitudes in the GH regime (A2=1) increase with respect to the non-fractional correlation in Figure 1. The amplitudes become smaller as the NS limit is approached (A2=10, 20).

Plot of the fractional correlation Nt;A2≡JFAt;A2/JFA0;A2$\mathcal {N}\left( {t; {A_2}} \right) \equiv {\mathcal {J}^{FA}}\left( {t; {A_2}} \right)/{\mathcal {J}^{FA}}\left( {0;{A_2}} \right)$ as defined by eq. (18) for A1=1 and A2=1 (-.-.-), 5 (_ _), 10 (…), 20 (____). In this case the curves do not oscillate and they vanish as t grows.
Figure 2:

Plot of the fractional correlation Nt;A2JFAt;A2/JFA0;A2 as defined by eq. (18) for A1=1 and A2=1 (-.-.-), 5 (_ _), 10 (…), 20 (____). In this case the curves do not oscillate and they vanish as t grows.

A more quantitative comparison between both behaviors is more clearly seen in Figure 3 where both, the non-fractional Mt;A2 and the fractional Nt;A2 correlation functions, are plotted vs. t for A2=1, 5, 10, 50. In the NS regime (A2=50) the fractional effects are almost negligible, however, as A2 decreases (τJ increases) and moves toward the GH limit, the amplitude of Nt;A2 increases. These variations can be quantified; for instance, if A2=1, in the interval 0.2<t<1.78 the amplitude of Nt;A2 varies between 0.69 and 0.25, which amounts to a maximum per cent difference (ΔC)MN=[(M(t;A2)-N(t;A2))/M(t;A2)]×100=29.6 %.

The fractional correlation function Nt;A2≡JFAt;A2/JFA0;A2$\mathcal {N}\left( {t; {A_2}} \right) \equiv {\mathcal {J}^{FA}}\left( {t; {A_2}} \right)/{\mathcal {J}^{FA}}\left( {0;{A_2}} \right)$ plotted vs. t for A2=1, 5, 10, 50. The solid black line is the non-fractional correlation Mt;A2$\mathcal {M}\left( {t; {A_2}} \right)$ and it is plotted as a reference.
Figure 3:

The fractional correlation function Nt;A2JFAt;A2/JFA0;A2 plotted vs. t for A2=1, 5, 10, 50. The solid black line is the non-fractional correlation Mt;A2 and it is plotted as a reference.

An important feature of these behaviors should be pointed out. Notice that since EIT gradually departs from the NS limit, it is to be expected that the behavior of Nt;A2 should reproduce some of the oscillations observed in the non-fractional solution Mt;A2. However, this is not observed in Figure 3. This suggests that the fractional description of Case A might not be a sufficiently good fractional extension of the standard non-fractional model. This was the motivation to introduce the second fractional model (Case B) described in subsection 3.2. Indeed, if the method of Case B is used, Figure 4 shows the time behavior of the fractional correlation function Pt;A2JFBt;A2/JFB0;A2, as calculated from eq. (35). Notice that the most striking feature is that now it reproduces part of the oscillatory behavior of the non-fractional case for A2=10, thus correcting the inability of the first fractional model (Case A). Moreover, for A2=20 both, the fractional and the non-fractional correlations are very similar, in agreement with the analytic result given by eq. (35), since it can be proven that this equation reduces to the non-fractional solution given by eq. (14) when ε=0. Although for A2=10 this similarity is less satisfactory, still exist.

The fractional intermediate scattering function Pt;A2≡JFBt;A2/JFB0;A2$\mathcal{P}\left( {t; {A_2}} \right) \equiv {J^{FB}}\left( {t; {A_2}} \right)/{J^{FB}}\left( {0;{A_2}} \right)$ for case B is plotted as a function of t and for A2=5, 10 and for ε=0.1. For reference the non-fractional (black solid) curve Mt;A2$\mathcal {M}\left( {t; {A_2}} \right)$ for A2=20 is also shown.
Figure 4:

The fractional intermediate scattering function Pt;A2JFBt;A2/JFB0;A2 for case B is plotted as a function of t and for A2=5, 10 and for ε=0.1. For reference the non-fractional (black solid) curve Mt;A2 for A2=20 is also shown.

The corresponding per cent differences (ΔC)MP=[(M(t;A2)-P(t;A2))/M(t;A2)]×100 can be quantified from these curves. For instance, for A2=10, t=1.8 and ε=0.4 we find that (ΔC)MP=−21.43 %, which means that the magnitude of the fractional correlation P(1.8;10) is greater than the non-fractional M (1.8;10) by a factor of ∼22 %. Thus, the fractional effect is not small. Likewise, for A2=5, t=2.8 and ε=0.4, this difference is even larger and amounts to ∼36 %.

At this point it is convenient to explain how the infinite series appearing in eq. (35) were approximated to calculate the fractional solutions shown in Figure 4. When we take eight terms in each of the infinite series in eq. (35), we obtain the dashed curve shown in Figure 5 for A2=20 and ε=0.01.

Plot of Mt;A2=Jt;A2/J0;A2$\mathcal {M}\left( {t; {A_2}} \right) = \mathcal {J}\left( {t; {A_2}} \right)/\mathcal {J}\left( {0;{A_2}} \right)$ and Pt;A2≡JFBt;A2/JFB0;A2$\mathcal {P}\left( {t; {A_2}} \right) \equiv {J^{FB}}\left( {t; {A_2}} \right)/{J^{FB}}\left( {0;{A_2}} \right)$ for A2=20 and ε=0.01. P$\mathcal {P}$8 results from the first eight terms of the series in eq. (33), whereas P$\mathcal {P}$10 arises from the first ten terms.
Figure 5:

Plot of Mt;A2=Jt;A2/J0;A2 and Pt;A2JFBt;A2/JFB0;A2 for A2=20 and ε=0.01. P8 results from the first eight terms of the series in eq. (33), whereas P10 arises from the first ten terms.

In this figure, we can observe a spurious growth of the fractional solution beyond t∼8 and beyond t∼10. It can be observed that as we increase the number of terms in the sums of eq. (31), the decaying oscillations of the fractional solution do not change appreciably and show a similar structure. The observed difference is that the time where the spurious growth starts is shifted toward larger times. Therefore, only eight terms are sufficient if we are only interested in calculating the solution in the interval 0<t<8. It is worth mentioning that if the values of A2 were different, more terms in the sums in eq. (31) might be needed. Thus, the approximation of considering only a few terms in the infinite series seems to be consistent.

5 Discussion

In this work, we reported model calculations to analyze the effect of including fractional time derivatives in the hydrodynamic equation for the intermediate scattering function of an isothermal-isobaric inert, inviscid and incompressible binary mixture in the GH regime of intermediate wave vectors k and frequencies ω. We summarize the results obtained and add some further comments to put them into a more proper context. We also suggest alternative ways to confront our theoretical results with experiment.

The expression for the non-fractional dynamics of the ISF was calculated analytically from an extended irreversible thermodynamic model. The effect produced by fractional derivatives on the dynamics of the ISF was also calculated analytically by considering two different ways of introducing the fractional derivatives. In the first approach, we replaced the two integer-order derivatives which appear in the standard model described by eq. (9), by a single fractional derivative of order 1<ε<2. In the second approach, the two derivatives appearing in eq. (9) were replaced by fractional derivatives of orders 2-ε and 1-ε, respectively. One outcome is that the second approach describes better the propagative character of the FISF because it reproduces the oscillatory behavior of ISF near the NS limit. We find that although the FISFs are in general less than the ISF, FISF may be a significant fraction (∼36 %) of ISF and might be measurable by light or neutron scattering techniques. As shown in Figure 5, this approach also shows quantitatively that the inclusion of fractional time scales into the dynamics of fluctuations may produce a FISF greater than the ISF showing that the FISF may be a measurable effect. These qualitative and quantitative behaviors for the model considered, are consistent with the ones previously reported for the behavior of other fractional correlation functions in complex fluids [19].

We have interpreted the fractional fluctuations as a superposition of conventional thermal fluctuations with successive weighted increments over time, but a deeper interpretation of the nature of fractional fluctuation is still being formulated [35].

Our results and predictions have been obtained for a realistic, although in several aspects a rather idealized model of a binary mixture in a nonequilibrium state. For this reason, it is difficult to know real experimental values for the phenomenological coefficients A1 and A2. We avoided this necessity by restricting our analysis to the time scale of the slowest non-conserved variable (A2) and by reformulating the transport equation only in terms of aˆ (eq. (28)) which depends on the ratio of both coefficients. In this way by setting A1=1 and by changing the value of A2 we could vary the nonequilibrium states in the GH regime, where the local equilibrium assumption is not strictly valid, to near-equilibrium states where the system is not that far from the NS limit.

There are or course many open questions and further possibilities regarding the predicted behavior of the calculated FISF, but its eventual validity can only come from experiments. Yet, to our knowledge there are no available experimental measurements of the ISF, the FISF or the (non-fractional) dynamic structure factor Sk,ω defined by eq. (15). Usually ISF or Sk,ω are measured by neutron scattering [7], but other experimental techniques such as forced Rayleigh scattering has been also used to measure the ISF in simple fluids [36, 37, 38]. If such experiments were available, they would provide alternative ways of comparing our theoretical predictions with experiments. From our theoretical results in this work we now suggest how a quantitative and analytic estimation of the FISF can be made when the initial sate of the system starts in the GH regime and approaches and gets very close to the NS limit. In this regime the system will eventually be accessible to measurements of ISF and FISF by using conventional or forced Rayleigh scattering, but still containing information of the GH states.

The general expression for FISF is given by eq. (35), and in the asymptotic limit t it reduces to eq. (40). As mentioned in the previous section, A2 controls the degree of departure of the state from the NS regime. Since A2(k)Dk2/τJ, for given D and k, small values of A2 correspond to a large relaxation time τJ, indicating that the fast variable J has not yet relaxed and prevent the validity of local equilibrium. In contrast, large values of A2 correspond to near local equilibrium states. Therefore, for given A2, k may be large or small depending on the value of D, which is unknown. Of course, since k is related to the scattering angle θ in a light scattering experiment, k=2sinθ/2, a choice of A2 may correspond to very small θ. In this case the above-mentioned forced Rayleigh scattering technique is suitable. In Figure 6 we illustrate this approximated but quantitative estimation of FISF. Note that all curves decay to the exponential amplitude Ik,τ as given by eq. (10). Since we have approximated eq. (37) by its leading term, the oscillations of JˆFBk,t are no longer perceptible. To observe them would require a more complicated approximation to eq. (37) involving more terms than the leading one. It should be pointed out that the decay rate of the different curves depends on the value of ε, which controls the order of fractional derivatives, and in this sense, reflects the extent of the memory or influence of other time scales.

Plot of the FISF JˆFBk,t→∞${\hat J^{FB}}\left( {k, t \to \infty } \right)$ vs t as given by eq. (37) for A2=20 and ε=0.01, 0.1, 0.5.
Figure 6:

Plot of the FISF JˆFBk,t vs t as given by eq. (37) for A2=20 and ε=0.01, 0.1, 0.5.

If such experiments would exist for our system, the comparison of our theoretical curves against experimental measurements would provide an additional confrontation between theory and experiment. However, we are not aware of any experiments of this sort and the experimental verification of our model predictions remains to be assessed.


  • [1]

    L. Onsager, Reciprocal relations in irreversible processes. I., Phys. Rev. 37 (1931), 405, 2265. Google Scholar

  • [2]

    S. R. De Groot and P. Mazur, Nonequilibrium Thermodynamic, North Holland, Amsterdam, 1962. Google Scholar

  • [3]

    L. D. Landau and E. M. Lifshitz, Fluid Mechanics, 2nd ed., Pergamon, Oxford, 1987. Google Scholar

  • [4]

    H. J. M. Hanley editor, Transport Phenomena in Fluids, M. Dekker, New York, 1971. Google Scholar

  • [5]

    E. L. Cussler, Multicomponent Diffusion, Elsevier, New York, 1976. Google Scholar

  • [6]

    R. D. Mountain and J. M. Deutch, Light scattering from binary solutions, J. Chem. Phys. 50 (1969), 1103–1108. CrossrefGoogle Scholar

  • [7]

    J. P. Boon and S. Yip, Molecular Hydrodynamics, McGraw-Hill, New York, 1980. Google Scholar

  • [8]

    R. P. Meilanov and R. A. Magomedov, Thermodynamics in fractional calculus, J. Eng. Phys. Thermophysics 87 (2014), 1521–1531. CrossrefGoogle Scholar

  • [9]

    S. L. Sobolev, Locally nonequilibrium transfer processes, Usp. Fiz. Nauk. 167 (1997), 1095–1106. Google Scholar

  • [10]

    D. Jou, C. Pérez-García, L. S. García-Colín, M. López De Haro and R. F. Rodríguez, Generalized hydrodynamics and extended irreversible thermodynamics, Phys. Rev. A 31 (1985), 2502–2508. CrossrefGoogle Scholar

  • [11]

    R. F. Rodríguez, L. S. García-Colín, M. López De Haro, D. Jou and C. Pérez-García, The underlying thermodynamic aspects of generalized hydrodynamics, Phys. Lett. A 107 (1985), 17–20. CrossrefGoogle Scholar

  • [12]

    R. F. Rodríguez, L. S. García-Colín and M. López De Haro, Extended thermodynamic description of diffusion in an inert binary mixture, J. Chem. Phys. 83 (1985), 4099–4102. CrossrefGoogle Scholar

  • [13]

    R. F. Rodríguez, M. López De Haro and L. S. García-Colín, Mutual diffusion in a binary mixture, Lect. Notes Phys. 253 (1986), 343–348. CrossrefGoogle Scholar

  • [14]

    B. J. West and S. Picozzi, Fractional Langevin model of memory in financial time series, Phys. Rev. E 65 (2002), Article 037106. CrossrefGoogle Scholar

  • [15]

    P. Pramukkul, A. Svenkeson, P. Grigolini, M. Bologna and B. J. West, Complexity and the fractional calculus, Adv. Math. Phys. 2013 (2013), Article 498789. Google Scholar

  • [16]

    B. J. West, Physiology, Promiscuity and Prophecy at the Millennium: A Tale of Tails, Studies of Nonlinear Phenomena in the Life Sciences, Vol. 7, World Scientific, Singapore, 1999. Google Scholar

  • [17]

    R. Metzler and J. Klafter, The random walk’s guide to anomalous diffusion: A fractional dynamics approach, Phys. Rep. 339 (2000), 1–77. CrossrefGoogle Scholar

  • [18]

    R. F. Rodríguez, J. Fujioka and E. Salinas-Rodríguez, Fractional fluctuation effects on the light scattered by a viscoelastic suspension, Phys. Rev. E 88 (2013), Article 022154. CrossrefWeb of ScienceGoogle Scholar

  • [19]

    R. F. Rodríguez, J. Fujioka and E. Salinas-Rodríguez, Fractional correlation functions in simple viscoelastic liquids, Physica A 427 (2015), 326–340. CrossrefWeb of ScienceGoogle Scholar

  • [20]

    G. Stolovitzky and K. R. Sreenivasan, Kolmogorov’s refined similarity hypotheses for turbulence and general stochastic processes, Rev. Mod. Phys. 66 (1994), 229–240. CrossrefGoogle Scholar

  • [21]

    R. Gorenflo, F. Mainardi and A. Vivoli, Continuous-time random walk and parametric subordination in fractional diffusion, Chaos Solitons & Fractals 34 (2007), 87–103. CrossrefWeb of ScienceGoogle Scholar

  • [22]

    M. G. Hernández, E. Salinas-Rodríguez, S. A. Gómez, J. A. E. Roa-Neri, S. Alfaro and F. J. Valdés-Parada, Helium permeation through a silicalite-1 tubular membrane, Heat Mass Transfer 51 (2015), 847–857. Web of ScienceCrossrefGoogle Scholar

  • [23]

    F. J. Valdés-Parada, J. A. Ochoa-Tapia, E. Salinas-Rodríguez, S. Gómez-Torres and M. G. Hernández, Upscaled model for dispersive mass transfer in a tubular porous membrane separator, Rev. Mex. Ingen. Quím. 13 (2014), 237–257. Google Scholar

  • [24]

    J. M. Van De Graaf, F. Kapteijn and J. A. Moulijn, Permeation of weakly adsorbing components through a silicalite-1 membrane, Chem. Eng. Sci. 54 (1999), 1081–1092. CrossrefGoogle Scholar

  • [25]

    B. J. Berne and R. Pecora, Dynamic Light Scattering, J. Wiley, New York, 1976. Google Scholar

  • [26]

    B. J. West, Colloquium: Fractional calculus view of complexity: A tutorial, Rev. Mod. Phys. 86 (2014), 1169–1189. CrossrefWeb of ScienceGoogle Scholar

  • [27]

    L. Podlubny, Fractional Differential Equations, Academic Press, New York, 1999. Google Scholar

  • [28]

    A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, 2006. Google Scholar

  • [29]

    M. Caputo and C. Cametti, Diffusion with memory in two cases of biological interest, J. Theoret. Biol. 254 (2008), 697–703. CrossrefGoogle Scholar

  • [30]

    S. G. Samko, A. A. Kilbas and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Beach, Amsterdam, 1993. Google Scholar

  • [31]

    K. S. Miller and B. Ross, An Introduction to Fractional Calculus and Fractional Differential Equations, Wiley, New York, 1993. Google Scholar

  • [32]

    K. B. Oldham and J. Spanier, The Fractional Calculus, Vol. 11, Academic Press, New York, 1974. Google Scholar

  • [33]

    M. D. Ortigueira, J. A. T. Machado and J. S. Da Costa, Which differ integration? IEE proceedings-vision, Image Signal Process 152 (2005), 846–850. Google Scholar

  • [34]

    J. Fujioka, A. Espinosa and R. F. Rodríguez, Fractional optical solitons, Phys. Lett. A 374 (2010), 1126–1134. CrossrefGoogle Scholar

  • [35]

    P. Pramukkul, A. Svenkeson, P. Grigolini, M. Bologna and B. West, Complexity and the fractional calculus, Adv. Math. Phys. 2013 (2013), Article 498789. Google Scholar

  • [36]

    J. P. Boon, C. Allain and P. Lallemand, Propagating thermal modes in a fluid under thermal constraint, Phys. Rev. Lett. 43 (1979), 199–203. CrossrefGoogle Scholar

  • [37]

    P. N. Segrè, R. W. Gammon, J. V. Sengers and B. M. Law, Rayleigh scattering in a liquid far from thermal equilibrium, Phys. Rev. A 45 (1992), 714–724. CrossrefGoogle Scholar

  • [38]

    B. M. Law, P. N. Segrè, R. W. Gammon and J. V. Sengers, Light-scattering measurements of entropy and viscous fluctuations in a liquid far from thermal equilibrium, Phys. Rev. A 41 (1990), 816–824. CrossrefGoogle Scholar

About the article

Received: 2017-07-11

Accepted: 2017-10-27

Revised: 2017-10-20

Published Online: 2017-12-14

Published in Print: 2018-01-26

Citation Information: Journal of Non-Equilibrium Thermodynamics, Volume 43, Issue 1, Pages 43–55, ISSN (Online) 1437-4358, ISSN (Print) 0340-0204, DOI: https://doi.org/10.1515/jnet-2017-0036.

Export Citation

© 2018 Walter de Gruyter GmbH, Berlin/Boston.Get Permission

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Rosalío Rodríguez, Elizabeth Salinas-Rodríguez, and Jorge Fujioka
Entropy, 2018, Volume 20, Number 1, Page 28

Comments (0)

Please log in or register to comment.
Log in