Entropy Generation through Non-Equilibrium Ordered ... - MDPI

entropy Article

Entropy Generation through Non-Equilibrium Ordered Structures in Corner Flows with Sidewall Mass Injection LaVar King Isaacson Mechanical Engineering, University of Utah, 2067 Browning Avenue, Salt Lake City, UT 84108, USA; [email protected]; Tel.: +1-801-583-1756 Academic Editor: Brian Agnew Received: 22 June 2016; Accepted: 25 July 2016; Published: 28 July 2016

Abstract: Additional entropy generation rates through non-equilibrium ordered structures are predicted for corner flows with sidewall mass injection. Well-defined non-equilibrium ordered structures are predicted at a normalized vertical station of approximately eighteen percent of the boundary-layer thickness. These structures are in addition to the ordered structures previously reported at approximately thirty-eight percent of the boundary layer thickness. The computational procedure is used to determine the entropy generation rate for each spectral velocity component at each of several stream wise stations and for each of several injection velocity values. Application of the procedure to possible thermal system processes is discussed. These results indicate that cooling sidewall mass injection into a horizontal laminar boundary layer may actually increase the heat transfer to the horizontal surface. Keywords: corner flow mass injection; non-equilibrium ordered structures; empirical entropies; entropy generation rates; laminar-turbulent transition

1. Introduction Entropy generation rates for corner flows with sidewall mass injection were reported by the author in a previous article [1]. These computations indicated the generation of non-equilibrium ordered structures at a normalized vertical station of approximately thirty-eight percent of the normalized laminar boundary layer thickness. The present article presents additional computations for entropy generation rates predicted at a normalized vertical station of approximately eighteen percent of the normalized laminar boundary layer thickness for corner flows with a variety of injection mass velocities. The flow configuration considered consists of a corner flow environment, with sidewall mass injection, with the mass injection uniform in the steam wise direction, along the vertical sidewall at the intersection of the sidewall with the horizontal surface. A schematic diagram for this configuration is shown in Figure 1. The first computational component is the calculation of the steady flow laminar velocity profiles along the horizontal surface in the stream wise-vertical (x-y) plane and the steady flow orthogonal boundary layer profile in the span wise-vertical (z-y) plane. These profiles are computed using the program listings provided by Cebeci and Bradshaw [2] and Cebeci and Cousteix [3]. The similarity characteristics for these orthogonal profiles proved by Hansen [4] justify the simultaneous use of these profiles in our three-dimensional flow configuration. The next computational component casts the time-dependent spectral equations of Townsend [5] and Hellberg et al. [6] into a Lorenz format (Sparrow [7]) to yield the nonlinear time series solutions for the fluctuating components of the spectral velocity field.

Entropy 2016, 18, 279; doi:10.3390/e18080279

www.mdpi.com/journal/entropy

Entropy 2016, 18, 279

2 of 23

Entropy 2016, 18, 279

2 of 2

y-axis Sidewall surface

Horizontal surface x-axis

Horizontal surface boundary layer we

ue

Main stream flow velocity, ue z-axis Crossflow velocity, we

Figure Figure 1. 1. Schematic Schematic of of the the corner corner flow flow boundary boundary layer layer configuration configuration is is shown shown [1]. [1].

The nonlinear time series solutions of the modified Lorenz equations yield both smooth The nonlinear time series solutions of the modified Lorenz equations yield both smooth trajectories trajectories and trajectories with fluctuating spectral velocity components within the series. Statistical and trajectories with fluctuating spectral velocity components within the series. Statistical analysis analysis of the fluctuating spectral time-series solutions yields the entropy generation rates for the of the fluctuating spectral time-series solutions yields the entropy generation rates for the given given flow configurations. The correspondence of the peaks of the spectral power density analysis flow configurations. The correspondence of the peaks of the spectral power density analysis and the empirical modes of the singular value decomposition analysis is brought about by the and the empirical modes of the singular value decomposition analysis is brought about by the Weiner–Khintchine theorem relating the power density spectrum and the autocorrelation function Weiner–Khintchine theorem relating the power density spectrum and the autocorrelation function for for the nonlinear time series data [8] (pp. 354–355). the nonlinear time series data [8] (pp. 354–355). The power spectral densities within the nonlinear time series solutions for the fluctuating The power spectral densities within the nonlinear time series solutions for the fluctuating spectral velocity components are computed using Burg’s method [9], providing the spectral peaks of spectral velocity components are computed using Burg’s method [9], providing the spectral peaks the ordered structures within the time series solutions [10] (pp. 572–574). The singular value of the ordered structures within the time series solutions [10] (pp. 572–574). The singular value decomposition method [10] (pp. 59–70) is also used to compute the empirical modes within the same decomposition method [10] (pp. 59–70) is also used to compute the empirical modes within the time series, corresponding to the power spectral density peaks found from Burg’s method. These same time series, corresponding to the power spectral density peaks found from Burg’s method. results yield the values for the empirical entropies for the ordered structures as defined by These results yield the values for the empirical entropies for the ordered structures as defined by Rissanen [11]. The empirical entropy for each of the modes from the singular value decomposition Rissanen [11]. The empirical entropy for each of the modes from the singular value decomposition process yields a corresponding entropic index of the format postulated by Tsallis [12]. Then, from process yields a corresponding entropic index of the format postulated by Tsallis [12]. Then, from each of these entropic indices, the method of Arimitsu and Arimitsu [13] yields a corresponding value each of these entropic indices, the method of Arimitsu and Arimitsu [13] yields a corresponding for the intermittency exponent. Combining these computational components into an overall value for the intermittency exponent. Combining these computational components into an overall computational procedure, we then compute the entropy generation rate for each of the ordered computational procedure, we then compute the entropy generation rate for each of the ordered structures identified in the nonlinear time series solutions of the non-equilibrium spectral equations. structures identified in the nonlinear time series solutions of the non-equilibrium spectral equations. The specific computational components are more fully described in Isaacson [1,14]. This article The specific computational components are more fully described in Isaacson [1,14]. This article includes includes the following sections: the following sections: In Section 2, the thermodynamic and transport processes of the working substance required for In Section 2, the thermodynamic and transport processes of the working substance required for the computational procedures are discussed. In this study, the working substance is the flow of air at the computational procedures are discussed. In this study, the working substance is the flow of air at a a specified temperature and pressure. In Section 3, the mathematical and computational bases for the specified temperature and pressure. In Section 3, the mathematical and computational bases for the computation of the steady three-dimensional boundary layer environment are reviewed. In computation of the steady three-dimensional boundary layer environment are reviewed. In Section 4, Section 4, the fluctuation equations of Townsend [5] are transformed into the spectral plane and the fluctuation equations of Townsend [5] are transformed into the spectral plane and written in the written in the Lorenz format. Section 5 presents computational results for the time-dependent Lorenz format. Section 5 presents computational results for the time-dependent spectral velocity spectral velocity components for various sidewall mass injection velocities at various stream wise components for various sidewall mass injection velocities at various stream wise locations. Section 6 locations. Section 6 discusses the extraction of empirical entropies, empirical entropic indices, and discusses the extraction of empirical entropies, empirical entropic indices, and intermittency exponents intermittency exponents from the nonlinear time series solutions of the modified Lorenz equations. from the nonlinear time series solutions of the modified Lorenz equations. Section 7 presents the Section 7 presents the results for the computation of the entropy generation rates through ordered structures in the boundary layer environment for various sidewall injection velocities and various stream wise locations. Section 8 includes a comparison of the entropy generation rates for two vertical

Entropy 2016, 18, 279

3 of 23

results for the computation of the entropy generation rates through ordered structures in the boundary layer environment for various sidewall injection velocities and various stream wise locations. Section 8 includes a comparison of the entropy generation rates for two vertical locations within the laminar boundary layer and compares these values with the entropy generated in a corresponding turbulent boundary layer. The article closes with a discussion of the results and final conclusions. 2. Selection of Heated Air as the Working Substance The flow configuration we are simulating in this exploratory study is shown in Figure 1. The figure is intended to represent the corner flow in the entrance to a square channel in, for example, an industrial heating unit. In this section, we describe the characteristics of the working gas we use in our computational procedure. The working gas is assumed to be heated air with the following thermodynamic and transport properties: Air temperature, T : 1068.0 K (1) Static pressure, p : 0.912 ˆ 105 N{m2

(2)

Kinematic viscosity, ν : 1.51634 ˆ 10-4 m2 {s

(3)

It should be noted that the value for the kinematic viscosity used in this study is very similar to the values used in previous studies. The explanation for this is that a series of coupled, nonlinear time-dependent differential equations are included in the computational procedure. The solutions for such equations are very sensitive to the initial conditions imposed on the equations and are also sensitive to the control parameters applied in the solution. The solution of the steady-state boundary-layer equations that provides these control parameters is dependent on the particular value of the kinematic viscosity applied in the calculations. We have found that only a narrow range of kinematic viscosities yields the prediction of ordered structures from the computational procedure [1]. The choice of air at the given temperature and pressure yields an appropriate value for the kinematic viscosity. The range of appropriate values of kinematic viscosities has not been delineated and remains as work to be done. A fundamental feature of the computational procedure discussed in this article is the embedding of the time-dependent Lorenz-format equations in the reservoir of steady-state boundary layer velocity gradients obtained from the time independent boundary layer equations. To clarify the essential nature of this feature, the next section presents a summary of the mathematical and computational methods used for the determination of the reservoir steady-state laminar boundary-layer velocity gradients. These boundary layer mean velocity gradients are time independent but vary with the stream wise distance x. The mean velocity gradients serve as control parameters for the time-dependent initiation of instabilities within the boundary layer for each stream wise station, as summarized in Section 4 of the article. 3. Steady-Flow Laminar Boundary-Layer Environment The boundary-layer configuration considered in this article consists of a laminar boundary layer in the x-y plane produced by the stream wise velocity along the horizontal surface and a laminar boundary layer in the z-y plane produced by mass injection from the sidewall in the z-direction. The momentum equation for the thin-shear boundary layer approximation may be written ([2], pp. 39–45) as: „  Bu Bu 1 dp 1 B Bu 1 1 u `v “´ ` µ ´ ρu v . (4) Bx By ρ dx ρ By By The boundary conditions for Equation (4) are: y “ 0, u “ v “ 0,

(5)

Entropy 2016, 18, 279

4 of 23

y “ δpxq, u “ ue pxq.

(6)

The Reynolds shear stress for the computation of turbulent boundary layers is modeled with the “eddy viscosity”, εm , having the dimensions of (viscosity)/(density), by: ´ ρu1 v1 “ ρεm

Bu . By

(7)

The computer program we have chosen to implement for the solution of the boundary layer equation (Equation (4)) is based on the Keller–Cebeci box method presented by Cebeci and Bradshaw [2] and Cebeci and Cousteix [3]. One of the basic aspects of this method is to transform Equation (4) into a system of first-order ordinary differential equations. The Falkner–Skan transformation, in the form: ´ u ¯1{2 e η“ y, (8) νx is introduced into the transformation process. The dimensionless stream function, f px, ηq, is defined by: Ψ px, yq “ pue νxq1{2 f px, ηq . (9) These definitions yield the results for the mean boundary layer velocities u and v as: u “ ue f 1 , v “ ´

B η ´ ue ν ¯1{2 1 rpue ν xq1{2 f s ` f . Bx 2 x

(10)

Differentiation with respect to η is indicated by the prime in these expressions. dp e From Bernoulli’s equation, the pressure gradient term is given by dx “ ´ρue du dx . To simplify the resulting equations, the parameter m is defined as: m“

x due . ue dx

(11)

Through the application of the transformation process, (we refer the reader to Cebeci and Bradshaw [2] for the details) the basic partial differential equation for the boundary layer (Equation (4)) is replaced with three first-order partial nonlinear differential equations in the following form: f 1 “ u,

(12)

u1 “ v, ˆ

˙

(13) ˆ

˙

m`1 Bu Bf f v ` mp1 ´ u2 q “ x u ´v . (14) 2 Bx Bx ˘ εm ` In Equation (14), b “ p1 ` ε` m and εm “ ν . The corresponding boundary conditions for Equations (12)–(14) are: fpx, 0q “ 0, upx, 0q “ 0, upx, η8 q “ 1 (15) pbvq1 `

Note that in Equation (13), v is not the y-component velocity. Cebeci and Bradshaw [2] present computer program listings for the numerical solutions for both laminar and turbulent boundary layers over flat plate surfaces. The program listings used in the study reported here are those presented in [2]. Hansen [4] has indicated that orthogonal laminar boundary layer profiles in a three-dimensional coordinate system possess the characteristic of similarity. We therefore use the boundary layer computations for both the profiles in the x-y plane and in the z-y plane. The steady state boundary layer equations and the corresponding velocity gradients serve as the thermodynamic steady state reservoir that provides the control parameters for the time-dependent development of the spectral fluctuations within the boundary layer environment [8] (pp. 26–31).

Entropy 2016, 18, 279

5 of 23

4. Modified Time-Dependent Lorenz Equations in the Spectral Plane 4.1. Transformation of the Townsend Equations to the Modified Lorenz Format The crosswind velocity of the mass injected in the z-direction from the sidewall produces a laminar boundary layer along the horizontal surface in the z-y plane. This boundary layer is orthogonal to the laminar boundary layer in the x-y plane. The nonlinear interactions of this z-y boundary layer within the laminar boundary layer in the x-y plane generate instabilities within the three-dimensional laminar boundary layer flow. The methods of Townsend [5] are used to obtain a set of non-equilibrium, time-dependent equations in the spectral domain for these non-equilibrium instabilities, separate from the equations for the steady boundary layer flow. Separating the equations of motion into steady and unsteady equations, the equations for the velocity fluctuations may then be written as [5,6]: Bui Bu BU Bu 1 Bp B 2 ui ` Uj i ` u j i ` u j i “ ´ `ν . Bt Bx j Bx j Bx j ρ Bxi Bu j Bu j

(16)

In these equations, ρ is the density and ν is the kinematic viscosity, Ui represent the mean boundary layer velocity components with i = 1, 2, 3 indicating the x, y, and z components, and xj , with j = 1, 2, 3, designate the x, y and z directions. The pressure term is eliminated by taking the divergence of Equation (16) and invoking incompressibility, yielding: ´

BU Bum Bul Bum 1 B2 p “2 l ` . ρ Bxl2 Bxm Bxl Bxm Bxl

(17)

We wish to transform these equations into the spectral plane for the computational solutions for the time-dependent fluctuating spectral components. The solutions of the spectral equations yield the spectral wave vector components and the spectral velocity wave components. Through Parseval’s theorem, the products of the spectral velocity components represent the products of the fluctuating velocity components in the physical plane. The statistical analysis of the spectral velocity component time series solutions then yields the entropy generation rates within the boundary layer non-equilibrium ordered structures. The fluctuating velocity and pressure fields of Equations (16) and (17) may be expanded in terms of the Fourier components (Mathieu and Scott [15]) as: ui px, tq “

ÿ

ai pk, tqexppik ‚ xq

(18)

ppx, tq ÿ “ bpk, tqexppik ‚ xq. ρ

(19)

k

and:

k

The pressure component in Equation (16) is transformed into a function of spectral velocity components and boundary layer velocity gradients through Equations (17) and (19). Substituting the resulting equations and Equation (18) into Equation (16) yields an expression for the fluctuations of the spectral components with time. The equations for the time-dependent three spectral velocity wave components, ai (k), are then given as: dai pkq dt

“

´νk2 ai pkq ´ BBUx i ai pkq ` 2 kki k2 l BBxUml ai pkq i

`i

ř k 1 `k 2 “k

pk l kikk2m ´ δim k l qal pk1 qam pk2 q

.

(20)

Entropy 2016, 18, 279

6 of 23

The general equations for the balance of transferable properties give the equations for the spectral wave numbers, ki as: BU dk i “ ´ l kl . (21) dt Bxi The set of equations for the time-dependent wave number components, including the gradients of the mean velocities in the x-y and z-y boundary layers, may be written: dk x BU BVx “´ kx ´ ky , dt Bx Bx

(22)

dk y BU BVx BW “´ kx ´ ky ´ kz , dt By By By

(23)

dk z BVz BW “´ ky ´ kz . dt Bz Bz

(24)

The nonlinear products of the spectral velocity components in Equations (20) are retained in our series of equations by characterizing the coefficients: k l pδim ´

ki km q k2

(25)

as a projection matrix [15]. This coefficient represents the projection of a given velocity wave vector component, ai , normal to the direction of the corresponding wave number component, ki . A model equation for this expression in the form: 1 ´ K¨ cospkptqq

(26)

is introduced to retain the effect of the projection matrix on the nonlinear interactive terms in our equations. K is an empirical weighting amplitude factor [14] and k(t) is given by: kptq “

b pk2x q.

(27)

With F = Kcos(k(t)), the equations for the spectral velocity components, Equations (20), are written in Lorenz format as [14]: da x “ σy a y ´ σ x a x , (28) dt day “ ´p1 ´ Fqa x az ` r1 a x ´ s1 ay , dt

(29)

daz “ p1 ´ Fqa x ay ´ b1 az . dt

(30)

From Equations (20), the coefficients of the velocity wave component terms have the following forms [14]: „ˆ ˙  2k x k y BU 2k x k y BVx 2k x k z BW σy “ ´ 1 ` ` , (31) By k2 k2 By k2 By " „ˆ ˙ * 2k x k x BU 2k x k y BVx 2k x k z BW σx “ υk2 ´ ´ 1 ` ` , (32) Bx k2 k2 Bx k2 Bx ˆ  „ ˙ 2k y k x BU 2k y k y BVx 2k y k z BW r1 “ ` ´ 1 ` , (33) Bx k2 Bx k2 k2 Bx " „ ˆ ˙ * 2k y k x BU 2k y k y BVz 2k y k z BW s1 “ υk2 ´ ` ´ 1 ` , (34) By k2 By k2 k2 By " „ ˆ ˙ * 2k z k z BW 2k z k x BU 2k z k y BVz ` ` ´ 1 b1 “ υk2 ´ . (35) Bz k2 Bz k2 Bz k2

Entropy 2016, 18, 279

7 of 23

The nonlinear time series solutions for the spectral velocity wave number components in the spectral domain for each of the seven downstream stations make up the non-equilibrium thermodynamic system. Solutions of the overall set of equations are obtained at each of the seven stations along the x-axis. The application of these equations at the initial station requires that additional assumptions be applied to the modified Lorenz equations. Isaacson [14], Mathieu and Scott [15], and Manneville [16] have discussed both the format and justification for the particular choice of these assumptions for the initial station. These solutions contain the spectral velocity component fluctuations that we wish to study. We apply our computational procedure to the six stations following the initial station at x = 0.08 along the stream wise direction. However, as Attard [8] (pp. 42–45) points out, we must take into account that the nonlinear time series solutions obtained for the second and subsequent stations will be influenced by the fluctuations produced in the first and following stations. To accomplish this, we use the synchronization properties of the modified Lorenz set of equations describing the nonlinear solutions for the spectral components [1,14]. 4.2. Synchronization Properties of the Modified Lorenz Equations Pecora and Carroll [17], Pérez and Cerdeiral [18], and Cuomo and Oppenheim [19] have used the synchronization properties of systems of Lorenz-type equations to extract messages masked by chaotic signals. The synchronization properties of Lorenz-type equations are adapted here to extract ordered signals from the nonlinear time series generated for each of the spectral components in the solutions of the modified Lorenz equations. As indicated previously, the solutions are obtained for a sequence of stations along the x-axis. We apply the transformation of the projection matrix (Equation (26)) to the initial station at x = 0.080. We then apply the synchronization properties to each of the following downstream stations. The various boundary layer coefficients at each of these stations are computed in the same manner as in the initial station. Following the results in [14], the time-dependent output for the x-direction spectral velocity component from the initial station is used as input to the nonlinear coupled terms in the modified Lorenz equations at the next station, which we denote as the second station in the x-direction. Then, the input to the nonlinear-coupled terms at the next downstream station is made up of the sum of the stream wise velocity wave component output from the initial station plus the x-direction spectral velocity wave component output from the next downstream station. This process is repeated for each of the six downstream stations. With this method, the memory of the initial velocity fluctuations and the influence of subsequent fluctuations are retained in the overall computational procedure. 4.3. Deterministic Results for the Modified Lorenz Equations The stream wise velocity has been chosen as ue = 1.00 for all of the cases considered in this study. We have chosen to present the results obtained at the station of x = 0.140 since the results for this station indicate the clearest power spectral density peaks in the overall set of statistical results. However, it turns out that the results for this station also indicate the lowest level of entropy generation rates at the normalized boundary layer location of η = 1.40 for the crosswind velocity of we = 0.0800. This choice for the presentation of the results should therefore be kept in mind in the comparison of entropy generation rates for the overall set of crosswind velocities. Figure 2 presents the deterministic three-dimensional spectral velocity wave trajectory for the downstream station of x = 0.140, for a sidewall mass injection velocity of we = 0.0800. Figure 3 shows the corresponding phase diagram for az4 –ay4 , where ay4 is the normal spectral velocity wave component and az4 is the span wise spectral velocity wave component, again at the station x = 0.140.

However, turns out normalized that the results for this station also the the lowest level of entropy generation itrates at the boundary layer location of indicate η = 1.40 for crosswind velocity of generation the normalized boundary of layer of η = 1.40 for thebecrosswind velocity of we = 0.0800.rates This at choice for the presentation the location results should therefore kept in mind in the w e = 0.0800. This choicegeneration for the presentation ofoverall the results therefore be kept in mind in the comparison of entropy rates for the set ofshould crosswind velocities. comparison entropy the generation rates for the overall set ofspectral crosswind velocities. Figure 2ofpresents deterministic three-dimensional velocity wave trajectory for the Figure 2 station presents deterministic three-dimensional spectral velocity trajectory the downstream of the x = 0.140, for a sidewall mass injection velocity of we =wave 0.0800. Figure 3 for shows Entropy 2016, 18, 279 station phase downstream of x = 0.140, for afor sidewall injection of we =spectral 0.0800. Figure 3 shows the corresponding diagram az4–ay4mass , where ay4 isvelocity the normal velocity wave 8 of 23 the corresponding diagram for az4–a y4, where is the normal velocity wave component and az4 isphase the span wise spectral velocity waveay4component, againspectral at the station x = 0.140. component and az4 is the span wise spectral velocity wave component, again at the station x = 0.140.

Figure 2. Shown is a three-dimensional deterministic trajectory of the spectral velocity components,

Figure 2. Shown is a three-dimensional deterministic trajectory of the spectral velocity components, az4, at is x =a 0.140 for η = 1.40 and we = 0.0800. trajectory of the spectral velocity components, ax4, ay4, and Figure 2. Shown three-dimensional deterministic ax4 , ay4 , and az4 , at x = 0.140 for η = 1.40 and we = 0.0800. ax4, ay4, and az4, at x = 0.140 for η = 1.40 and we = 0.0800.

y4

Normal spectral velocity component, a a Normal spectral velocity component,

y4

15.0 15.0 10.0 10.0 5.0 5.0 0.0 0.0 -5.0 -5.0 -10.0 -10.0 -15.0 -15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 -15.0 -15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 Span wise spectral velocity component, a z4

Entropy 2016, 18, 279

9 of 9

Span wise spectral velocity component, a z4

Figure 3. The phase diagram of the span wise and normal spectral velocity components, az4–ay4, is Figure 5 shows corresponding phase diagram fornormal az4–ay4, where avelocity y4 is thecomponents, normal spectral shown x the = 0.140, ηdiagram = 1.40ofand e = span 0.0800. Figure 3.phase The phase of w the wiseand and az4–ay4velocity is –a , is Figure 3. Thefor diagram the span wise normalspectral spectral velocity components, a, z4 y4 waveshown component and a z4 is the span wise spectral velocity wave component, again at the station x = 0.140 for x = 0.140, η = 1.40 and we = 0.0800.

shown for x = 0.140, η = 1.40 and we = 0.0800. and crosswind velocity the of wdeterministic e = 0.0925. Figure 4 presents three-dimensional trajectory of the spectral velocity wave Figure 4for presents the deterministic three-dimensional trajectory of injection the spectral wave components the downstream station of x = 0.140, for a sidewall mass rate velocity of we = 0.0925. components for the downstream station of x = 0.140, for a sidewall mass injection rate of we = 0.0925.

Figure 4. A three-dimensional representation of the deterministic trajectories of the spectral velocity

Figure 4. A three-dimensional representation of the deterministic trajectories of the spectral velocity components, ax4, ay4, and az4, is shown for x = 0.140, η = 1.40 and we = 0.0925. components, ax4 , ay4 , and az4 , is shown for x = 0.140, η = 1.40 and we = 0.0925. y4

velocity compoment, a

80.0 60.0 40.0 20.0 0.0

Entropy 2016, 18, 279

9 of 23

Figure 4 presents the deterministic three-dimensional trajectory of the spectral velocity wave components for the downstream station of x = 0.140, for a sidewall mass injection rate of we = 0.0925. Figure 5 shows the corresponding phase diagram for az4 –ay4 , where ay4 is the normal spectral velocity waveFigure component and az4 is the representation span wise spectral velocity wave component, the station 4. A three-dimensional of the deterministic trajectories of the again spectralatvelocity x = 0.140 and crosswind velocity of w = 0.0925. e for x = 0.140, η = 1.40 and we = 0.0925. components, ax4, ay4, and az4, is shown

y4

Normal spectral velocity compoment, a

80.0 60.0 40.0 20.0 0.0 -20.0 -40.0 -60.0 -80.0 -80.0

-60.0

-40.0

-20.0

0.0

20.0

40.0

60.0

80.0

Span wise spectral velocity component, a z4

Figure 5. The The phase phase diagram diagram of of the the span span wise wise and andnormal normalspectral spectralvelocity velocitycomponents, components,aaz4z4–a –ay4 y4,, is is Figure 5. shown for x = 0.140, η = 1.40 and w e = 0.0925. shown for x = 0.140, η = 1.40 and we = 0.0925.

5. 5. Power PowerSpectral SpectralDensity Densitywithin withinthe theNon-Equilibrium Non-EquilibriumSpectral SpectralVelocity Velocity Fluctuations Fluctuations 5.1. Structures 5.1. Power Power Spectral Spectral Densities Densities Indicating Indicating Ordered Ordered Structures The Lorenz synchronization synchronization procedure the The Lorenz procedure yields yields the the nonlinear nonlinear time time series series solution solution for for each each of of the six stations in the x-direction, following the initial station. These solutions indicate the generation of six stations in the x-direction, following the initial station. These solutions indicate the generation deterministic spiral structures, with significant non-equilibrium velocity fluctuations embedded of deterministic spiral structures, with significant non-equilibrium velocity fluctuations embedded within some some of of the the solutions. solutions. Entropy these within Entropy generation generation comes comes about about through through the the dissipation dissipation of of these fluctuations into equilibrium thermodynamic states. We have found that Burg’s method [9] for the fluctuations into equilibrium thermodynamic states. We have found that Burg’s method [9] for the power spectral spectral density the time-series time-series data is an an effective effective method method for for extracting extracting the the underlying underlying power density of of the data is structural characteristics of of the the velocity velocity fluctuations fluctuations within within the the nonlinear nonlinear time time series series solutions. solutions. structural characteristics The resulting power spectral density results for the normal spectral velocity wave component, ay4 at the fourth station at x = 0.140, are presented in Figure 6. For each of the power spectral density results, we have assigned empirical mode numbers to these peaks, starting with mode j = 1 representing the highest peak in the distribution, continuing to mode j = 16, representing the corresponding lowest peak among the sixteen peaks. The peaks of the spectral power density analysis and the empirical modes of the singular value decomposition analysis are computed from the same set of time series data. Therefore, we can relate each power spectral density peak with a corresponding empirical mode through the Weiner´Khintchine theorem relating the power density spectrum and the autocorrelation function for the nonlinear time series data [8] (pp. 354–355). The resulting power spectral density results for the span wise spectral velocity wave component, ay4 at the fourth station at x = 0.140, are presented in Figure 7.

Entropy 2016, 18, 279

4.50 10

10 of 10

0

y4

0

4.50 10 0 3.00 10

y4

Power spectral density for a spectral density for a Power

The resulting power spectral density results for the normal spectral velocity wave component, ay4 at the fourth station at 4.00 x = 10 0.140, are presented in Figure 6. For each of the power spectral density 0 results, we have assigned empirical mode numbers to these peaks, starting with mode j = 1 representing Entropy 2016, 18, 279 the highest peak in the distribution, continuing to mode j = 16, representing the 10 of 23 0 3.50 10 corresponding lowest peak among the sixteen peaks.

0

4.00 10 0 2.50 10

0

3.50 10 0 2.00 10

3.00 10

0

20

40

0

60

80

100

120

140

Frequency step number

Figure 6. The power spectral density for the normal spectral velocity component is shown for 0 e = 0.0800. x = 0.140, η = 1.40 and w2.50 10

The peaks of the spectral power density analysis and the empirical modes of the singular value 0 2.00computed 10 decomposition analysis are from the same set of time series data. Therefore, we can relate 0 20 40 60 80 100 120 140 each power spectral density peak with a corresponding empirical mode through the Frequency step number and the autocorrelation function Weiner−Khintchine theorem relating the power density spectrum for the nonlinear time series data [8] (pp. 354–355). 6. The power spectral density for the normal spectral velocity component is shown for Figure Figure 6. The powerpower spectral density for theresults normal component is wave showncomponent, for x = 0.140, The resulting spectral density forspectral the spanvelocity wise spectral velocity x = 0.140, η = 1.40 and we = 0.0800. η =ay41.40 andfourth we = 0.0800. at the station at x = 0.140, are presented in Figure 7.

0

9.00 101 1.20 10

0

8.00 101 1.10 10

z4

Powerspectral density for a Powerspectral density for a

z4

The peaks of the spectral power density analysis and the empirical modes of the singular value 1 10 decomposition analysis are 1.20 computed from the same set of time series data. Therefore, we can relate each power spectral density peak with a corresponding empirical mode through the 1 1.10 10 Weiner−Khintchine theorem relating the power density spectrum and the autocorrelation function for the nonlinear time series data [8] (pp. 354–355). The resulting power spectral 1.00 101 density results for the span wise spectral velocity wave component, ay4 at the fourth station at x = 0.140, are presented in Figure 7.

Entropy 2016, 18, 279

0

7.00 10 1.00 101

0

6.00 100 9.00 10 0

50

100

150

200

250

11 of 11

300

Frequency numbervelocity wave component, ay4 at the The power spectral density0 results for the normalstep spectral 8.00 10 fourth station at x = 0.140 for a crosswind velocity of we = 0.0925, are presented in Figure 8. The Figure 7. The power spectral density for the span wise spectral velocity component is shown for Figure 7. The power spectral density forfor the component is shown resulting power spectral density results thespan spanwise wisespectral spectral velocity velocity wave component, ay4 atfor x = 0.140, η = 1.40 and we = 0.0800. 0 10 x =the 0.140, η =at1.40 wefor = 0.0800. station x = and 0.140, a7.00 crosswind velocity of we = 0.0925, are presented in Figure 9. 8.00 10

3

7.00 10

3

6.00 10

3

4.00 10

3

3.00 10

3

2.00 10

3

1.00 10

3

y4

Power spectral density for a

6.00 10

0

0

50

100

150

200

250

300

Frequency step number

Figure 7. The power spectral density for the span wise spectral velocity component is shown for x = 0.140, η = 1.40 and w5.00 e = 0.0800. 3 10

0

20

40

60

80

100

120

140

Frequency step number

Figure 8. The power spectral density for the normal spectral velocity component is shown for

2.00 10

5

1.50 10

5

z4

sity for a

Figure 8. The power spectral density for the normal spectral velocity component is shown for x = 0.140, x = 0.140, η = 1.40 and we = 0.0925. η = 1.40 and we = 0.0925.

Power spectra

Entropy 2016, 18, 279

4.00 10

3

3.00 10

3

2.00 10

3

1.00 10

3

11 of 23 0

20

40

60

80

100

120

140

2.00 10

5

1.50 10

5

1.00 10

5

5.00 10

4

z4

Power spectral density for a

The power spectral density results for the normal spectral velocity wave component, ay4 at the Frequency step number fourth station at x = 0.140 for a crosswind velocity of we = 0.0925, are presented in Figure 8. The resulting power spectral wisespectral spectralvelocity velocitycomponent wave component, Figure 8. The power density spectralresults densityfor forthe thespan normal is shownay4 forat the station x = η0.140, aw crosswind x = at 0.140, = 1.40for and e = 0.0925. velocity of we = 0.0925, are presented in Figure 9.

0.00

0

50

100

150

200

250

300

Frequency step number

Figure forfor thethe span wise spectral velocity component is shown for Figure 9. 9. The Thepower powerspectral spectraldensity density span wise spectral velocity component is shown 0.0925. xfor = 0.140, η = 1.40 and and we = w x = 0.140, η = 1.40 e = 0.0925.

A significant advantage of Burg’s method (the maximum entropy method) is the enhancement A significant advantage of Burg’s method (the maximum entropy method) is the enhancement of of the spectral peaks in the power spectral density distribution. Press et al. [10] (pp. 572–575) present the spectral peaks in the power spectral density distribution. Press et al. [10] (pp. 572–575) present computer program listings for the prediction of the power spectral density using the maximum computer program listings for the prediction of the power spectral density using the maximum entropy method. The results indicate that the power spectral density is distributed over a selection of entropy method. The results indicate that the power spectral density is distributed over a selection of well-defined spectral peaks. These spectral peaks represent a collection of ordered structures well-defined spectral peaks. These spectral peaks represent a collection of ordered structures embedded embedded within the fluctuating time series solutions of the modified Lorenz equations, each within the fluctuating time series solutions of the modified Lorenz equations, each identified with a identified with a corresponding mode number. corresponding mode number. These spectral peaks have been organized in descending order into sixteen empirical modes. These spectral peaks have been organized in descending order into sixteen empirical modes. Simpson’s integration rule is used to obtain the kinetic energy available for dissipation for each of the Simpson’s integration rule is used to obtain the kinetic energy available for dissipation for each of the power spectral density modes. The total spectral energy available for dissipation is then obtained as power spectral density modes. The total spectral energy available for dissipation is then obtained as the sum of the individual contributions across the modes. This value is then used to get the fraction the sum of the individual contributions across the modes. This value is then used to get the fraction of dissipation energy in each mode that is eventually available for dissipation into the internal energy of the final equilibrium thermodynamic state. 5.2. Empirical Entropies from Singular Value Decomposition The singular value decomposition procedure [20] can also serve as a means of identifying fundamental characteristics of the nonlinear time series solutions of the coupled nonlinear spectral equations. We have incorporated into our numerical procedure the singular value decomposition computer program listings presented by Press et al. [10] (pp. 59–65). The computational procedure is made up of two parts, the computation of the autocorrelation matrix and the singular value decomposition of that matrix [10]. The overall computational procedure yields the empirical eigenvalues for each of the empirical eigenfunctions for the given nonlinear time series data segment. The application of the singular value decomposition procedure to a specified segment of the nonlinear time-series solution for each of the spectral velocity wave components yields the distribution of the component eigenvalues λ j across the empirical modes, j, for the flow conditions

Entropy 2016, 18, 279

12 of 23

listed in Equations (1)–(3). The empirical entropy, Sempj , is defined from these eigenvalues by the expression [11]: Semp j “ ´lnpλ j q. (36) Here, λj is the empirical eigenvalue computed from the singular value decomposition procedure applied to the nonlinear time-series solution. The distribution of the empirical entropy across the decomposition empirical modes has been shown in [14]. We have applied the singular value decomposition procedure to the nonlinear time-series solutions of the modified Lorenz non-equilibrium equations. Again, invoking Parseval’s formula [14], the empirical eigenvalues produced by the singular value decomposition procedure represent twice the kinetic energy associated with the fluctuating velocity fields [20]. 5.3. Empirical Entropic Indices for the Ordered Structures The empirical entropies for the fluctuating spectral velocity wave component time series indicate different characteristics for the various non-equilibrium collections embedded within the time series. These results indicate that the majority of the kinetic energy is contained within the first six empirical modes of the singular value decompositions, with relatively low empirical entropy. The following three empirical modes indicate a transition into non-equilibrium ordered structures. These structures have been classified as coherent [20] with well-defined structural boundaries. We wish to find a way to follow these structures through the process of transition from ordered structures into equilibrium thermodynamic states. For this, we turn to the concept of the Tsallis entropic format (Tsallis [12]). The Tsallis entropic format is applicable to a variety of microscopic and macroscopic sub-systems. We are working with a collection of statistical sub-systems spread over a limited number of empirical modes, j. Note that this is not a sequence over time, but is a sequence in phase space [8] (pp. 409–412). The empirical entropy, Sempj is introduced to describe the entropy of an ordered structure described by the empirical eigenvalue, λ j , for the singular value decomposition empirical mode, j. Hence, we simply adopt, in an ad hoc fashion, an expression from which we may extract an empirical entropic index, qj , from the empirical entropy, Sempj . This expression is written as [14]: Semp j “ ´lnpλ j q “

1 ´ pλ j qq j . pq j ´ 1q

(37)

This expression includes the effects of the nonlinear, non-equilibrium nature of the ordered structures we are following. It has the format of an entropic index; hence, we simply call qj the empirical entropic index or simply entropic index. The empirical entropic indices for the normal spectral velocity wave component, ay4 at the fourth station at x = 0.140, with a crosswind velocity of we = 0.08000 are presented as a function of the empirical mode j in Figure 10. The empirical entropic indices for the normal spectral velocity wave component, ay4 at the fourth stream wise station at x = 0.140, and a crosswind velocity of 0.0925 are presented as a function of the empirical mode, j in Figure 11. 5.4. Empirical Intermittency Exponents for the Ordered Structures In this section, we introduce a heuristic method to connect the non-equilibrium results for the entropic indices with the final phase of the dissipation of fluctuating kinetic energy into thermodynamic internal energy. We explore this computational connection through the concept of intermittency exponents and a relaxation process into the final thermodynamic entropy state. The non-equilibrium ordered structures discussed in previous sections are of a macroscopic nature embedded within the nonlinear time series solutions of the nonlinear equations for the fluctuating spectral velocity wave field. Singular value decomposition of the time series solutions provides empirical entropies for these non-equilibrium structures. The empirical entropic indices of the Tsallis form extracted from the empirical entropies using Equation (37) have been used to

Entropy 2016, 18, 279

13 of 23

Entropy 2016, 18, 279 Entropy 2016, 18, 279

13 of 13 13 of 13

obtain the intermittency exponents for the ordered structures. We heuristically apply a relationship, indices for the velocity wave aay4 at the fourth stream station at found andspectral Arimitsu [13], connecting the entropic index of Tsallis to thewise intermittency indicesby forArimitsu the normal normal spectral velocity wave component, component, y4 at the fourth stream wise station at xexponent, aa crosswind velocity of are as of empirical mode, jj in ζj . This intermittency exponent describes the fraction of fluctuating energy within x == 0.140, 0.140, and and crosswind velocity of 0.0925 0.0925 are presented presented as aa function function of the thekinetic empirical mode, in Figure 11. the non-equilibrium ordered structure that is dissipated into thermodynamic internal energy [13]. Figure 11. 0.00 0.00 -0.50 -0.50

j

Entropyindex, index,qq Entropy j

-1.00 -1.00 -1.50 -1.50 -2.00 -2.00 -2.50 -2.50 -3.00 -3.00 -3.50 -3.50 -4.00 -4.00 0 0

2 2

4 4

6 6

8 8

10 10

12 12

14 14

16 16

Empirical mode, j Empirical mode, j

Figure Figure 10. 10. The normal spectral spectral velocity velocity component as aa The empirical empirical entropic entropic index index for for the the normal component is is shown shown as function function of of the the empirical empirical mode, mode, j,j, for for xx == 0.140, 0.140, ηη η == 1.40 1.40 and and w weee===0.0800. 0.0800. 0.0800. 0.00 0.00

j

Entropicindex, index,qq Entropic j

-1.00 -1.00

-2.00 -2.00

-3.00 -3.00

-4.00 -4.00

-5.00 -5.00 0 0

2 2

4 4

6 6

8 8

10 10

12 12

14 14

16 16

Empirical mode, j Empirical mode, j

Figure empirical entropic index for the normal is shown as Figure 11. 11. The The normal spectral spectral velocity velocity component component as aaa Figure 11. The empirical empirical entropic entropic index index for for the the normal spectral velocity component is is shown shown as function η == 1.40 and w e = 0.0925. function of of the the empirical empirical mode, mode, j,j,j, for for xxx == = 0.140, 0.140, function of the empirical mode, for 0.140, η η = 1.40 1.40 and and w wee ==0.0925. 0.0925.

5.4. 5.4. Empirical Empirical Intermittency Intermittency Exponents Exponents for for the the Ordered Ordered Structures Structures We will substitute the absolute value of the empirical entropic index discussed in the previous In this we introduce aa heuristic method to the non-equilibrium results for In into this section, section, we equation introducederived heuristic method and to connect connect the[13]. non-equilibrium results for the the section the original by Arimitsu Arimitsu This expression is written as: entropic indices with the final phase of the dissipation of fluctuating kinetic energy into entropic indices with the final phase of the dissipation of fluctuating kinetic energy into a a ζj thermodynamic this the ˇ ˇ energy. 1 `We ζ j ´explore log2 p1 ` ´ 2´ζj q ¨ log2 p1connection ´ 1 ´ 2´through q thermodynamic internal internal energy. We explore this 1computational computational connection through the concept concept of of ˇ ˇ a a q “ 1 ´ (38) intermittency exponents the final thermodynamic entropy state. j and a relaxation process into ´ζ j the final thermodynamic intermittency exponents and a relaxation log2 p1 `process 1 ´ 2into q ´ log2 p1 ´ 1 ´ 2´ζj q entropy state. The The non-equilibrium non-equilibrium ordered ordered structures structures discussed discussed in in previous previous sections sections are are of of aa macroscopic macroscopic nature embedded within the nonlinear time series solutions of the nonlinear equations nature embedded within the nonlinear time series solutions of the nonlinear equations for for the the fluctuating spectral velocity wave field. Singular value decomposition of the time series solutions fluctuating spectral velocity wave field. Singular value decomposition of the time series solutions provides provides empirical empirical entropies entropies for for these these non-equilibrium non-equilibrium structures. structures. The The empirical empirical entropic entropic indices indices of of the Tsallis form extracted from the empirical entropies using Equation (37) have been used to obtain the Tsallis form extracted from the empirical entropies using Equation (37) have been used to obtain

equilibrium ordered structure that is dissipated into thermodynamic internal energy [13]. We will substitute the absolute value of the empirical entropic index discussed in the previous We will substitute the absolute value of the empirical entropic index discussed in the previous section into the original equation derived by Arimitsu and Arimitsu [13]. This expression is section into the original equation derived by Arimitsu and Arimitsu [13]. This expression is written as: written as: −ζ

−ζ

1 + ζ − log (1 + 1 − 2 −ζ j ) ⋅ log (1 − 1 − 2 −ζ j ) 14 of 23 q j = 1 − 1 + ζ jj − log 22 (1 + 1 − 2 j ) ⋅ log 22 (1 − 1 − 2 j ) (38) −ζ j −ζ j qj = 1− (38) log 2 (1 + 1 − 2−ζ j ) − log 2 (1 − 1 − 2 −ζ j ) log 2 (1 + 1 − 2 ) − log 2 (1 − 1 − 2 ) Figure 12 shows the results for the intermittency exponents obtained from Equation (38) for the Figure 12 shows the results for the intermittency exponents obtained from Equation (38) for the stream wise wise location location of of xx == 0.140 0.140 with with aa crosswind crosswindvelocity velocityof ofwwe e==0.0800. 0.0800.These These results indicate stream results indicate a stream wise location of x = 0.140 with a crosswind velocity of we = 0.0800. These results indicate a a dominant intermittency exponentatatempirical empiricalmode modej j==33and and intermittency intermittency exponents exponents for for the dominant intermittency exponent dominant intermittency exponent at empirical mode j = 3 and intermittency exponents for the remaining empirical empiricalmodes modesclose closetotothe theobserved observedvalue value 0.24 [13]. intermittency exponents remaining ofof 0.24 [13]. TheThe intermittency exponents for remaining empirical modes close to the observed value of 0.24 [13]. The intermittency exponents for for the stream wise location of x = 0.140 with a crosswind velocity of w = 0.0925 are shown in the stream wise location of x = 0.140 with a crosswind velocity of we = 0.0925e are shown in Figure 13. the stream wise location of x = 0.140 with a crosswind velocity of we = 0.0925 are shown in Figure 13. Figure results 13. These results again indicate a dominant intermittency exponent at empirical mode and These again indicate a dominant intermittency exponent at empirical mode j =j =33 and These results again indicate a dominant intermittency exponent at empirical mode j = 3 and intermittency exponents for the remaining empirical modes close to the observed value of 0.24 [13]. intermittency exponents for the remaining empirical modes close to the observed value of 0.24 [13]. Entropy 2016, 18, 279

0.90 0.90

Intermittency exponent Intermittency exponent

0.80 0.80 0.70 0.70 0.60 0.60 0.50 0.50 0.40 0.40 0.30 0.30 0.20 0.20 0.10 0.10 0 0

2 2

4 4

6 6

8 8

10 10

12 12

14 14

16 16

Empirical mode, j Empirical mode, j

Figure 12. The empirical intermittency exponent for the normal spectral velocity component is shown Figure 12. 12. The The empirical empirical intermittency intermittency exponent exponent for for the the normal normal spectral spectral velocity velocity component component is is shown shown Figure e = 0.0800. as a function of the empirical mode, j, for x = 0.140, η = 1.40 and w as aa function function of of the the empirical empirical mode, mode, j, j, for for xx == 0.140, 0.140, ηη == 1.40 1.40 and and w wee==0.0800. 0.0800. as 0.90 0.90

Intermittency exponent Intermittency exponent

0.80 0.80 0.70 0.70 0.60 0.60 0.50 0.50 0.40 0.40 0.30 0.30 0.20 0.20 0.10 0.10 0 0

2 2

4 4

6 6

8 8

10 10

12 12

14 14

16 16

Empirical mode, j Empirical mode, j

Figure 13. The empirical intermittency exponent for the normal spectral velocity component is shown Figure 13. empirical intermittency exponent for normal spectral velocity Figure 13. The The intermittency exponent forηthe the normal velocity component component is is shown shown as a function ofempirical the empirical mode, j, for x = 0.140, = 1.40 and spectral we = 0.0925. e = 0.0925. as a function of the empirical mode, j, for x = 0.140, η = 1.40 and w as a function of the empirical mode, j, for x = 0.140, η = 1.40 and w = 0.0925. e

6. Entropy Generation Rates through the Ordered Structures 6.1. Kinetic Energy Available for Dissipation The local mean flow kinetic energy, u2 /2, at the normalized vertical distance, η = 1.40 in the x-y plane boundary layer, is considered as the source of kinetic energy to be dissipated through the spiral structures. This available kinetic energy is distributed over the stream wise component, the normal component and the span wise component. The fraction of kinetic energy in the x-direction velocity

Entropy 2016, 18, 279

15 of 23

component is denoted as κx , the fraction of kinetic energy in the y-direction velocity component is denoted as κy in the z-direction velocity component denoted as κz . The fraction of dissipation kinetic energy within each empirical mode of the power spectral energy distribution is denoted as ξj . Then the total rate of dissipation of the available fluctuating kinetic energy for the stream wise, normal and span wise velocity components is the summation, over the empirical modes, j, of the product of the kinetic energy fraction of each mode times the intermittency exponent for that mode, ζj [1]. The empirical intermittency exponent for each of the empirical modes within the ordered structures has been extracted from the empirical expression (Equation (38)) given by Arimitsu and Arimitsu [13]. At this point in the computational procedure, values are available for the input energy source for the non-equilibrium ordered structures, the fraction of the fluctuation kinetic energy available in each of the empirical modes within the non-equilibrium ordered structures, and the fraction of the energy in each of the empirical modes that dissipates into background thermal energy, thus increasing the thermodynamic entropy. We consider the dissipation process for the ordered structures as a general relaxation process and use concepts from non-equilibrium thermodynamics to describe this dissipation process. 6.2. Entropy Generation Rates for Relaxation Processes From the concepts of non-equilibrium thermodynamics, de Groot and Mazur [21] (pp. 221–230) write the equation for the entropy generation rate in an internal relaxation process as: Bs Jpxq Bµpxq “´ . Bt ρT Bx

(39)

Here, s is the entropy per unit mass, ρ is the density, T is the temperature, µ is the mechanical potential for the transport of the ordered structures in an external context and J(x) is the flux of dissipation energy through the ordered structures available for dissipation into thermal internal energy. The dissipation of the ordered structures into background thermal energy may be considered as a two-stage process from the transition of the ordered structures into equilibrium thermodynamic states and a relaxation process of the downstream velocity in the initial state to the final equilibrium state of the velocity over the internal distance x. At the final equilibrium state, the dissipated ordered structure vanishes into thermal equilibrium with the reservoir. The local boundary layer steady state velocity is written as u = ue f’, where f” is the derivative of the Falkner–Skan stream function f with respect to the normalized distance η. The expression for the entropy generation rate (in J/(m3 ¨ K¨ s)) through the non-equilibrium ordered structures is then written as [1]: » ˆ

u2e

‚ — 1 S gen “ ρ – 2 T

˙

$ ¨ fi ˛ , ¨ ˛ ¨ ˛ / ’ 16 16 16 . & ´ ¯ ÿ ÿ ÿ ` 1 ˘2 ue ffi f κ x ˝ ξ j ζ j ‚ ` κy ˝ ξ j ζ j ‚ ` κz ˝ ξ j ζ j ‚ fl . / ’ - x % j “1 j “1 j“1 x

y

(40)

z

In this expression, ρ is the density of the working substance, in this case air at the given pressure and temperature. The first two terms within the brackets in Equation (40) represent the source of total kinetic energy available from the local steady boundary layer velocity in the stream wise direction. The first factor in each of the three terms within the braces represents the fraction within each of the velocity components of the total kinetic energy from the boundary layer. The first factor in each of the summation terms represents the fraction of kinetic energy within the given spectral mode from the power spectral density distribution. The second factor in each summation is the corresponding intermittency exponent obtained from the empirical entropic index for that mode. The final factor within the brackets is the overall mechanical turnover rate that transports the ordered structures in the given flow environment. The kinetic energy in each spectral mode available for final dissipation into equilibrium internal energy is computed for each of the spectral peaks. The empirical entropy for each of the structures

3

Entropy generation rate, joules / (m K s)

the power spectral density distribution. The second factor in each summation is the corresponding intermittency exponent obtained from the empirical entropic index for that mode. The final factor within the brackets is the overall mechanical turnover rate that transports the ordered structures in the given flow environment. The kinetic Entropy 2016, 18, 279 energy in each spectral mode available for final dissipation into equilibrium internal 16 of 23 energy is computed for each of the spectral peaks. The empirical entropy for each of the structures indicated by the spectral peaks is found from the singular value decomposition process applied to indicated the series spectral is found from the singular value process the givenbytime datapeaks segment. The connecting parameter, the decomposition empirical entropic index,applied is then to theextracted given time data segment. parameter, the empirical entropic index, is then fromseries the resulting value ofThe the connecting empirical entropy. Glandsdorff and Prigogine findempirical that for the general evolution criterion for non-equilibrium extracted from the resulting value[22] of the entropy. processes, dSempand j/dt Prigogine < 0. When[22] thefind Tsallis index isevolution negative, criterion Mariz [23] that the Glandsdorff thatentropic for the general forfound non-equilibrium empiricaldSemp entropy change is also negative, dSemp j/dt < 0. The results presented in Figures 10 and 11 processes, /dt < 0. When the Tsallis entropic index is negative, Mariz [23] found that the empirical j indicate that significant non-equilibrium structures exist within the specified time frame of the entropy change is also negative, dSempj /dt < 0. The results presented in Figures 10 and 11 indicate that particularnon-equilibrium nonlinear time structures series solution. These the regions may time therefore as ordered, significant exist within specified framebeof classified the particular nonlinear non-equilibrium structures. Therefore, the significant negative nature for the extracted empirical time series solution. These regions may therefore be classified as ordered, non-equilibrium structures. entropic indices at the fourth stationnature at x = 0.140 is in agreement with both the Prigogine criterion Therefore, the significant negative for the extracted empirical entropic indices at the and fourth the Mariz results for the Tsallis entropic index. The ad hoc introduction of an empirical entropy station at x = 0.140 is in agreement with both the Prigogine criterion and the Mariz results for the Tsallis index thus provides a representation of the nonlinear, non-equilibrium ordered structures in a entropic index. The ad hoc introduction of an empirical entropy index thus provides a representation significant way. of the nonlinear, non-equilibrium ordered structures in a significant way. Given the absolute value of the empirical entropic index, qj, the intermittency exponent, ζj for Given the absolute value of the empirical entropic index, qj , the intermittency exponent, ζj for the the mode, j, is extracted from Equation (38) [13] by the use of Brent’s method [10] (pp. 397–405). mode, We j, is extracted from Equation (38) [13] by the use of Brent’s method [10] (pp. 397–405). consider the entropy generation rate to be a phase transition of the dissipation of the We consider ordered the entropy generation rate to be a phase transition of the dissipation non-equilibrium structure into background thermal energy and a relaxation process of of thethe non-equilibrium ordered structure into background thermal energy and a relaxation of the global mechanical transport process of the stream wise velocity, ue over the distance x. process At the final global mechanical transport process of the stream wise velocity, u over the distance x. At the final e equilibrium state, the stream wise velocity of the dissipated structure vanishes. The entropy equilibrium state, the stream wise velocity of the dissipated structure vanishes. The entropy generation generation rates through the non-equilibrium ordered structures are then determined from rates through Equation (40)the [1].non-equilibrium ordered structures are then determined from Equation (40) [1]. The entropy of the thedownstream downstreamstations stations are shown Figure The entropygeneration generationrates rates as as a function function of are shown in in Figure 14, 14, with fourvalues valuesofofsidewall sidewallinjection injection velocities velocities as Note that with thethe four as the theimposed imposedcrosswind crosswindparameters. parameters. Note that each crosswindvelocity velocityisisconstant constantin in the the x-direction. x-direction. each crosswind 1.60 10

-4

1.40 10

-4

1.20 10

-4

1.00 10

-4

8.00 10

-5

6.00 10

-5

4.00 10

-5

2.00 10

-5

we = 0.0800 we = 0.0830 we = 0.0921 we = 0.0925

0.00 0.080

0.120

0.160

0.200

0.240

0.280

0.320

Stream wise distance, x

Figure Theentropy entropygeneration generation rates rates for for the the non-equilibrium at at a normalized Figure 14.14.The non-equilibriumordered orderedstructures structures a normalized η = 1.40 are shown as a function of the x-direction stations for several applied vertical station vertical station η = 1.40 are shown as a function of the x-direction stations for several applied crosswind velocities. crosswind velocities.

Figure 14 shows only the entropy generation rates at various stream wise stations, excluding the initial station at x = 0.080. It is interesting to note that the rate of entropy generation for the stream wise station at x = 0.140 is considerably lower than the rates for the other stations. Also, as indicated in Figure 14, at x = 0.010, entropy generation rates were obtained for only the two lowest crosswind velocity values. These results indicate that the deterministic spiral structures for the larger crosswind velocity values are primarily smooth functions, with much lower levels of fluctuation energy than for the following stream wise stations.

2.00 10

-4

1.60 10

-4

1.20 10

-4

8.00 10

-5

3

Entropy generation rate, joules / (m K s)

initial station at x = 0.080. It is interesting to note that the rate of entropy generation for the stream wise station at x = 0.140 is considerably lower than the rates for the other stations. Also, as indicated in Figure 14, at x = 0.010, entropy generation rates were obtained for only the two lowest crosswind velocity values. These results indicate that the deterministic spiral structures for the larger crosswind Entropy 2016, 18, 279 of 23 velocity values are primarily smooth functions, with much lower levels of fluctuation energy17than for the following stream wise stations. Figure Figure15 15 shows shows the the entropy entropy generation generationrates ratesat at various various stream stream wise wise stations, stations, for for aa normalized normalized boundary layer distance of η = 3.00 [1]. Also, as indicated in Figure 15, at x = 0.080 boundary layer distance of η = 3.00 [1]. Also, as indicated in Figure 15, at x = 0.080and andxx ==0.100, 0.100, entropy generation rates were obtained for only two of the crosswind velocity values. Again, these entropy generation rates were obtained for only two of the crosswind velocity values. Again, these results resultsindicate indicatethat thatthe the deterministic deterministicspiral spiralstructures structuresfor forthe theother othertwo two crosswind crosswindvelocities velocitiesat at this this station are primarily smooth functions, with much lower levels of fluctuation energy than for the station are primarily smooth functions, with much lower levels of fluctuation energy than for the other other stream wise stations. stream wise stations.

we = 0.0860 we = 0.086375 we = 0.08675

4.00 10

we = 0.087125

-5

0.00 0.040

0.080

0.120

0.160

0.200

0.240

0.280

Stream wise distancw, x

Figure Figure 15. 15. The The entropy entropy generation generationrates ratesfor for the the non-equilibrium non-equilibriumordered orderedstructures structuresat at aa normalized normalized vertical verticalstation stationηη=3.00 =3.00are areshown shownas asaafunction functionofofthe thex-direction x-directionstations stationsfor forseveral severalapplied appliedcrosswind crosswind velocities [1]. velocities [1].

The abrupt increases in the entropy generation rates at the third station at x = 0.120 for the The abrupt increases in the entropy generation rates at the third station at x = 0.120 for the normalized distances of η = 1.40 and η = 3.00 indicate a burst of fluctuation energy in the nonlinear normalized distances of η = 1.40 and η = 3.00 indicate a burst of fluctuation energy in the nonlinear time series solutions for the spectral velocity wave components at this station. time series solutions for the spectral velocity wave components at this station. 6.3. 6.3.Entropy EntropyGeneration GenerationRates Ratesfor foraa Turbulent Turbulent Boundary Boundary Layer Layer For Foraa comparison comparisonof ofthese thesevalues valuesfor for the the entropy entropy generation generationrates, rates, the the entropy entropy generation generationrates rates within a turbulent boundary layer are computed for each given stream wise location. Moore within a turbulent boundary layer are computed for each given stream wise location. Moore and and Moore Moore[24] [24] give give the the entropy entropy generation generation rate rate near near the the wall wall in in aa turbulent turbulent boundary boundary layer layer as: as: • τ ˆ ∂u˙ ‚ S tbl =τw w Bu . Stbl “ ρTe  ∂y  w . ρTe By w

(41) (41)

Introducing the skin friction coefficient as: Introducing the skin friction coefficient as: Cf “

τw ρu2e 2

(42)

and applying the Falkner–Skan transformation (Equation (8)) to the velocity gradient, we may write: ‚

Stbl “ ρ

ˆ

1 u2e 2 T

˙´

¯ u 1 e C f p q pRex q 2 p f 2 q. x

(43)

2

and applying the Falkner–Skan transformation (Equation (8)) to the velocity gradient, we may write: 1 •  1 u2  u S tbl = ρ  e  ( C f ) ( e ) ( Re x ) 2 ( f '') . x 2 T 

Entropy 2016, 18, 279

(43)

18 of 23

The often-quoted expression for the skin friction coefficient for a turbulent boundary layer on a The may often-quoted expression flat plate be written as [25]: for the skin friction coefficient for a turbulent boundary layer on a flat plate may be written as [25]: 0.0592( Reqx´)1−{1/55 (44) C f C“f = 0.0592pRe (44) x Substituting into Equation Equation (43) (43) yields: yields: Substituting this this expression expression for for the the skin skin friction friction coefficient coefficient into 2 ˆ  1 2u˙ • ‚ .30 ¯ 1 u e  uuee ´ .0592( Re )00.30 ( fp ''). x q f 2 q. Stbl S“tbl ρ= ρ  2 eT p( x q) 00.0592pRe x 2 Te e  x

(

)

(45) (45)

2.00 10

-4

S_dot_tbl

3

Entropy generation rate, joules / (m K s)

This This expression expression is is used used to to compute compute the the entropy entropy generation generation rates rates across across aa hypothetical hypothetical turbulent turbulent boundary layer as a function of the normalized distance η along the horizontal boundary layer as a function of the normalized distance η along the horizontal surface surface of of the the corner corner flow the turbulent boundary layer begins at the initial station x= flow configuration. configuration.The Thecomputation computationofof the turbulent boundary layer begins at the initial station 0.02 with transition enforced at this location. Hence, the turbulent boundary layer for our calculations x = 0.02 with transition enforced at this location. Hence, the turbulent boundary layer for our at the stream at wise = 0.140 is much for smaller a naturally occurring transition further calculations thelocation stream xwise location x =smaller 0.140 isthan much than for a naturally occurring along the x-direction. The distribution of the entropy generation rates across a flat plate turbulent transition further along the x-direction. The distribution of the entropy generation rates across a flat boundary layer boundary is shown in Figure 16. in Figure 16. plate turbulent layer is shown

vw=0.0800

1.50 10

-4

vw=0.0830 vw=0.0921 vw=0.0925 vw=0.0860

1.00 10

-4

vw=0.086375 vw=0.086750 vw=0.087125

5.00 10

-5

0.00 0.00

1.00

2.00

3.00

4.00

5.00

6.00

7.00

8.00

Normalized distance from the surface, η

Figure Figure 16. 16. The The entropy entropy generation generation rates rates for for various various crosswind crosswind velocities velocities are are shown shown at at normalized normalized distances of η = 1.40 and η = 3.00 at a distance of x = 0.140. Also shown is the distribution of the distances of η = 1.40 and η = 3.00 at a distance of x = 0.140. Also shown is the distribution of theentropy entropy generation generation rate rate at at aa distance distance of of xx == 0.140 0.140 for for aa turbulent turbulent boundary boundary layer layer initiated initiated at at xx == 0.02 0.02 on on the the horizontal surface. horizontal surface.

The contrast in the fundamental meaning of entropy generation rates between the two sets of The contrast in the fundamental meaning of entropy generation rates between the two sets of data shown in Figure 16 should be clarified. The entropy generation rates computed for the turbulent data shown in Figure 16 should be clarified. The entropy generation rates computed for the turbulent boundary layer use a well-established empirical relation for the turbulent skin friction coefficient boundary layer use a well-established empirical relation for the turbulent skin friction coefficient (Equation (44)) and an empirical eddy viscosity expression for the computation of the mean turbulent (Equation (44)) and an empirical eddy viscosity expression for the computation of the mean turbulent velocity profile across the boundary layer. These types of results are useful to the designer of thermal power equipment for the estimation of the loss in stagnation pressure due to turbulent irreversible processes within the system. On the other hand, the computation of the entropy generation rates for the ordered structures proceeds from the identification of time-dependent instabilities within the laminar boundary layer due to nonlinear interactions with an imposed crosswind velocity, through to the evaluation of the dissipation of energy within non-equilibrium ordered structures. With the availability of

Entropy 2016, 18, 279

19 of 23

an appropriate turbulent boundary layer velocity profile, the enhanced skin friction coefficient produced by these ordered structures may then be obtained from the resulting entropy generation rate values. This approach thus introduces basic physical processes into the evaluation of time-dependent non-equilibrium entropy generation processes occurring in this three-dimensional laminar flow environment. 7. Discussion It has previously been found that a crosswind velocity in a corner boundary layer flow could trigger instabilities in the three-dimensional laminar boundary layer along the horizontal surface [1]. This flow into the boundary layer provides the development of an additional boundary layer profile in the z–y plane along the horizontal surface, thus providing a three-dimensional orthogonal boundary layer environment. Application of coupled, nonlinear modified Lorenz equations in the spectral plane of the flow environment indicates the generation of non-equilibrium spiral structures in the three-dimensional environment. Fluctuating spectral velocity components are found within the three spectral velocity component time-series solutions for these equations. Statistical processing of the solutions indicates the presence of ordered structures embedded within the fluctuating time-series environment. The transition of these ordered structures into equilibrium thermodynamic states yields the entropy generation rates for the overall process. For the results reported here, the effects of four sidewall injection velocity values are evaluated, we = 0.0800, 0.0830, 0.0921, 0.925, at seven stream wise stations along the x-axis. Each value of crosswind velocity is applied uniformly through the sidewall along the edge of the stream wise boundary layer. The synchronization properties of the modified Lorenz equations are applied from the second to the seventh x-direction stations. Computational results for the initial x-direction station and for the subsequent synchronized stations indicate significant non-equilibrium ordered structure generation at each station at the normalized vertical distance of η = 1.40 in the laminar boundary layer. The intensity of the power spectral densities within these non-equilibrium structures varies from station to station, indicating the triggering of significantly enhanced structures for several stations. It should also be noted that not all of the spectral velocity component time series solutions exhibit fluctuating content. Only the contributions from fluctuating time series solutions from the synchronized stream wise stations are included in the computation of the respective rates of entropy generation. For example, only the fluctuating spectral velocity components for the two lowest injection velocities contribute to the entropy generation at the second downstream station. At the second stream wise station, x = 0.100, only fluctuations from the lowest injection velocity and the highest injection velocity contribute to the entropy generation, yielding a significantly lower rate of entropy generation. At the third stream wise station, x = 0.120, the four crosswind velocities indicate high rates of entropy production. The following station, x = 0.140, indicates a significant decrease in entropy generation rates. However, at this station, the power spectral density distribution indicates sharper spectral peaks. We have presented in this article the sequence of the statistical analysis yielding the entropy generation rates for the ordered structures predicted at the stream wise station x = 0.140 at the normalized boundary layer distance of η = 1.40. To gain a perspective on the magnitude of the predicted rates of entropy generation through the transition of ordered structures, a comparison is made with the rates of entropy generation in a flat plate turbulent boundary layer for the same given flow conditions. The distribution of the entropy generation rates across a turbulent boundary layer at the distance of x = 0.140 is computed for a turbulent boundary layer initiated at x = 0.02 from the leading edge of the horizontal surface. The entropy generation rates through the ordered structures at normalized vertical distances of η = 1.40 and η = 3.0 are then compared with the turbulent boundary layer distribution. At the normalized vertical distance of η = 1.40, the entropy generation rates through the ordered structures are of the same magnitude as the values generated in the turbulent boundary layer. However, at the vertical

Entropy 2016, 18, 279

20 of 23

distance of η = 3.00, the entropy generation rates for the ordered structures considerably exceed those predicted for the turbulent boundary layer. Thus, the influence of the crosswind velocity on the horizontal laminar boundary layer is to initiate instabilities within the laminar boundary layer and to considerably enhance the resulting skin friction coefficient and resulting heat transfer to the horizontal surface [26]. Although the results presented in this article are obtained for the flow configuration of corner flows with sidewall mass injection, the computational procedure is quite generic in its application to general dynamic flow environments. For example, consider the first of three components in the formation of turbulent spots in the transition of laminar to turbulent flows [27]. Initially, a stream wise counter-rotating vortex is formed from a point-wise initiating obstacle. A stream wise laminar boundary layer is initiated on the right side of this vortex, along the floor of the cavity. The counter-rotating vortex produces a crosswind surface boundary layer orthogonal to the stream wise boundary layer [28]. The computational procedure used in the present study appears to be directly applicable to this dynamic flow configuration. The remaining components in the dynamic processes of the transition of laminar to turbulent flows should also be amenable to analysis with the present computational procedure. However, there are several significant issues concerning the computational procedure discussed in this article. First, the range of kinematic viscosity values for which instabilities have been observed in three-dimensional boundary layer configurations is very narrow [1]. The range of kinematic viscosities that may yield instabilities should be explored to delineate the applicability of the procedure. Second, there has been no validation of the computational procedure by independent computational results. Inclusion of the time-dependent modified Lorenz equations in the computational procedure should yield additional insight into the process of the transition of laminar to turbulent flow. An independent computational validation would be welcome. A third issue is that experimental verification of the predicted results is sparse. Results presented in Boiko et al. [29] (pp.120–121) indicate that small amplitude vortices of approximately ten percent of the free-stream velocity in a flat plate boundary layer have been observed. In our calculations, the amplitude of each of the injection crosswind velocities is less than ten percent of the free-stream velocity. In an important experimental measurement, Boiko et al. [29] (pp.115–118) reported the observation of turbulent spots at two vertical locations within a laminar boundary layer, similar to our prediction of the generation of two regions of non-equilibrium ordered structures within the three-dimensional laminar boundary layer in a corner flow with sidewall mass injection. The vertical locations within the boundary layer of the turbulent spots reported in [29] correspond to the vertical locations of the ordered structures reported in this study. This is an important experimental consistency for the results obtained through our computational procedure. 8. Conclusions An innovative computational procedure for the three-dimensional laminar boundary layer in a corner flow with sidewall mass injection indicates the generation of two regions of non-equilibrium ordered structures within the horizontal laminar boundary layer configuration. The entropy generation rates through the dissipation of these ordered structures are compared. The vertical locations of these ordered structures in the laminar boundary layer are compared with reported experimental observations. Detailed results for the entropy generation rates generated through the dissipation of the ordered structures near the wall are reported. Comparison of the entropy generation rates through the ordered structures with a simulated turbulent boundary layer indicates that the ordered structures near the wall are of similar magnitude, but the dissipation rates for the structures further out are quite strong. These results indicate that transpiration cooling from sidewall injection into a laminar boundary layer may actually increase the heat transfer to the horizontal surface through the generation of these ordered structures. The inclusion of time-dependent non-equilibrium equations

Entropy 2016, 18, 279

21 of 23

into a reservoir of steady laminar boundary layer profiles may provide additional insight into the transition of laminar to turbulent flow. Acknowledgments: The author would like to thank the editor and referees for their timely and helpful suggestions. Incorporation of these suggestions has significantly strengthened the presentation. Conflicts of Interest: The author declares no conflict of interest.

Nomenclature ai b b b1 Cf F j J k k ki K m p qj r1 s s1 Sempj ‚ S gen ‚ Stbl t u u’ ue ui Ui v v’ we x xi xj y z

Fluctuating i-th component of velocity wave vector Coefficient in Equation (14) Fluctuating Fourier component of the static pressure Coefficient in modified Lorenz equations defined by Equation (35) Skin-friction coefficient Time-dependent feedback factor Mode number empirical eigenvalue Net source of kinetic energy dissipation rate, Equation (39) Time-dependent wave number magnitude Dimensional constant, Equation (27) Fluctuating i-th wave number of Fourier expansion Adjustable weighting factor Pressure gradient parameter, Equation (11) Local static pressure Empirical entropic index for the empirical entropy of mode, j Coefficient in modified Lorenz equations defined by Equation (33) Entropy per unit mass Coefficient in modified Lorenz equations defined by Equation (34) Empirical entropy for empirical mode, j Entropy generation rate through kinetic energy dissipation Entropy generation rate in a turbulent boundary layer Time Mean stream wise velocity in the x-direction in Equation (4) Fluctuating stream wise velocity in Equation (4) Stream wise velocity at the outer edge of the x-y plane boundary layer The i-th component of the fluctuating velocity Mean velocity in the i-th direction in the modified Lorenz equations Mean normal velocity in Equation (4) Fluctuating normal velocity in Equation (4) Span wise velocity at the outer edge of the z-y plane boundary layer Stream wise distance i-th direction j-th direction Normal distance Span wise distance

Greek Letters δ δlm εm

Boundary layer thickness Kronecker delta Eddy viscosity

Entropy 2016, 18, 279

ε` m ζj η λj µ ν ξj ρ σy σx τw

22 of 23

Normalized eddy viscosity Intermittency exponent for the j-th mode in Equation (38) Transformed normal distance parameter Eigenvalue for the empirical mode, j Mechanical potential in Equation (39) Kinematic viscosity of the gas mixture Kinetic energy dissipation rate in the j-th empirical mode Density Coefficient in modified Lorenz equations defined by Equation (31) Coefficient in modified Lorenz equations defined by Equation (32) Wall shear stress

Subscripts E i, j, l, m x y z

Outer edge of the laminar boundary layer Tensor indices Component in the x-direction Component in the y-direction Component in the z-direction

References 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17.

Isaacson, L.K. Entropy Generation through Non-equilibrium Spiral Structures in Corner Flows with Sidewall Surface Mass Injection. Entropy 2016, 18, 47. [CrossRef] Cebeci, T.; Bradshaw, P. Momentum Transfer in Boundary Layers; Hemisphere: Washington, DC, USA, 1977; pp. 319–321. Cebeci, T.; Cousteix, J. Modeling and Computation of Boundary-Layer Flows; Horizons: Long Beach, CA, USA, 2005. Hansen, A.G. Similarity Analyses of Boundary Value Problems in Engineering; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964; pp. 86–92. Townsend, A.A. The Structure of Turbulent Shear Flow, 2nd ed.; Cambridge University Press: Cambridge, UK, 1976. Hellberg, C.S.; Orszag, S.A. Chaotic behavior of interacting elliptical instability modes. Phys. Fluids 1988, 31, 6–8. [CrossRef] Sparrow, C. The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors; Springer: New York, NY, USA, 1982. Attard, P. Non-Equilibrium Thermodynamics and Statistical Mechanics: Foundation and Applications; Oxford University Press: Oxford, UK, 2012. Chen, C.H. Digital Waveform Processing and Recognition; CRC Press: Boca Raton, FL, USA, 1982; pp. 131–158. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes in C: The Art of Scientific Computing, 2nd ed.; Cambridge University Press: Cambridge, UK, 1992. Rissanen, J. Information and Complexity in Statistical Modeling; Springer: New York, NY, USA, 2007. Tsallis, C. Introduction to Nonextensive Statistical Mechanics; Springer: New York, NY, USA, 2009; pp. 37–43. Arimitsu, T.; Arimitsu, N. Analysis of fully developed turbulence in terms of Tsallis statistics. Phys. Rev. E 2000, 61, 3237–3240. [CrossRef] Isaacson, L.K. Transitional Intermittency Exponents through Non-equilibrium Boundary-Layer Structures and Empirical Entropic Indices. Entropy 2014, 16, 2729–2755. [CrossRef] Mathieu, J.; Scott, J. An Introduction to Turbulent Flow; Cambridge University Press: Cambridge, UK, 2000; pp. 251–261. Manneville, P. Non-equilibrium Structures and Weak Turbulence; Academic Press: San Diego, CA, USA, 1990. Pecora, L.M.; Carroll, T.L. Synchronization in chaotic systems. In Controlling Chaos: Theoretical and Practical Methods in Nonlinear Dynamics; Kapitaniak, T., Ed.; Academic Press: San Diego, CA, USA, 1996; pp. 142–145.

Entropy 2016, 18, 279

18. 19.

20. 21. 22. 23. 24.

25. 26. 27. 28. 29.

23 of 23

Pérez, G.; Cerdeiral, H.A. Extracting messages masked by chaos. In Controlling Chaos: Theoretical and Practical Methods in Nonlinear Dynamics; Kapitaniak, T., Ed.; Academic Press: San Diego, CA, USA, 1996; pp. 157–160. Cuomo, K.M.; Oppenheim, A.V. Circuit implementation of synchronized chaos with applications to communications. In Controlling Chaos: Theoretical and Practical Methods in Nonlinear Dynamics; Kapitaniak, T., Ed.; Academic Press: San Diego, CA, USA, 1996; pp. 153–156. Holmes, P.; Lumley, J.L.; Berkooz, G.; Rowley, C.W. Turbulence, Coherent Structures, Dynamical Systems and Symmetry, 2nd ed.; Cambridge University Press: Cambridge, UK, 2012. De Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics; North-Holland: Amsterdam, Holland, 1962. Glansdorff, P.; Prigogine, I. Thermodynamic Theory of Structure, Stability and Fluctuations; Wiley: London, UK, 1971. Mariz, A.M. On the irreversible nature of the Tsallis and Renyi entropies. Phys. Lett. A 1992, 165, 409–411. [CrossRef] Moore, J.; Moore, J.G. Entropy Production Rates from Viscous Flow Calculations: Part I—A Turbulent Boundary Layer Flow. In Proceedings of the 28th ASME International Gas Turbine Conference and Exhibit, Phoenix, AZ, USA, 27–31 March 1983. White, F.M. Viscous Fluid Flow; McGraw-Hill: New York, NY, USA, 1974. Bayley, F.J.; Turner, A.B. The Transpiration-Cooled Gas Turbine. J. Eng. Power 1970, 92, 351–358. [CrossRef] Cantwell, B.J. Organized Motion in Turbulent Flow. Ann. Rev. Fluid Mech. 1981, 13, 457–515. [CrossRef] Ersoy, S.; Walker, J.D.A. Viscous flow induced by counter-rotating vortices. Phys. Fluids 1985, 28, 2687–2698. [CrossRef] Boiko, A.V.; Grek, G.R.; Dovgal, A.V.; Kozlov, V.V. The Origin of Turbulence in Near-Wall Flows; Springer: Berlin/Heidelberg, Germany, 2010. © 2016 by the author; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).