**Open Journal of Fluid Dynamics**

Vol.04 No.02(2014), Article ID:46731,43 pages

10.4236/ojfd.2014.42015

Multimoment Hydrodynamics in Problem on Flow around a Sphere: Entropy Interpretation of the Appearance and Development of Instability

Igor V. Lebed

Zhukovsky Central Institute of Aerohydrodynamics, Moscow, Russia

Email: lebed-i@newmail.ru

Copyright © 2014 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 31 March 2014; revised 30 April 2014; accepted 8 May 2014

ABSTRACT

Multimoment hydrodynamics equations are applied to investigate the phenomena of appearance and development of instability in problem on a flow around a solid sphere at rest. The simplest solution to the multimoment hydrodynamics equations coincides with the Stokes solution to the classic hydrodynamics equations in the limit of small Reynolds number values,. Solution to the multimoment hydrodynamics equations reproduces recirculating zone in the wake behind the sphere having the form of an axisymmetric toroidal vortex ring. The solution remains stable while the entropy production in the system exceeds the entropy outflow through the surface confining the system. The passage of the first critical value is accompanied by the solution stability loss. The solution, when loses its stability, reproduces periodic pulsations of the periphery of the recirculating zone in the wake behind the sphere. The and solutions to the multimoment hydrodynamics equations interpret a vortex shedding. After the second critical value is reached, the solution at the periphery of the recirculating zone and in the far wake is replaced by the solution. In accordance with the solution, the periphery of the recirculating zone periodically detached from the core and moves downstream in the form of a vortex ring. After the attainment of the third critical value, the solution at the periphery of the recirculating zone and in the far wake is replaced by the solution. In accordance with the solution, vortex rings penetrate into each other and form the continuous vortex sheet in the wake behind the sphere. The replacement of one unstable flow regime by another unstable regime is governed the tendency of the system to discover the fastest path to depart from the state of statistical equilibrium. Having lost the stability, the system does not reach a new stable position. Such a scenario differs from the ideas of classic hydrodynamics, which interprets the development of instability in terms of bifurcations from one stable state to another stable state. Solutions to the multimoment hydrodynamics equations indicate the direction of instability development, which qualitatively reproduces the experimental data in a wide range of Re values. The problems encountered by classic hydrodynamics when interpreting the observed instability development process are solved on the way toward an increase in the number of principle hydrodynamic values.

**Keywords:**

Multimoment Hydrodynamics, Instability, Entropy, Vortex Shedding

1. Introduction

Detailed evaluating the results of the direct numerical integration of the Navier-Stokes equations against experiment in problem on flow past a hard sphere at rest is carried out in [1] . Experiment records three stable medium states. Each of these three stable flows begins to develop in its own direction, qualitatively different from other flows when it loses stability. Calculation satisfactory reproduces all three stable medium states observed experimentally. However, the calculation is incapable of producing anything that corresponds to seven unstable regimes observed along the three directions of instability development. The analysis of numerous divergences between the results of numerical integration of the Navier-Stokes equations and the experiment [1] -[3] led to the following conclusion. Solutions to the classic hydrodynamics equations successfully reach the border of the instability field represented by the dashed slanting line in Figure 1 from [3] . As Reynolds number grows, these solutions move along the border of the field. However, the solutions to the classic hydrodynamics equations are unable to cross this border and to pass into the instability field.

The Navier-Stokes equations themselves are called as the most probable reason for discrepancies between calculation and experiment [1] - [4] . It may be likely that the Navier-Stokes equations become inapplicable to unstable phenomena. The responsibility for the failure of the classic hydrodynamics was laid on the approximation made in deriving the Boltzmann equation, namely, the hypothesis of molecular chaos “Stosszahlansatz”. The Boltzmann hypothesis closes the kinetic equation, allowing classic hydrodynamics to be constructed for only three lower principle hydrodynamic values. It turns out that the neglect of higher principle hydrodynamic values does not introduce visually observable changes into the picture of stable flows. This error, however, grows very rapidly after the loss of stability.

In present paper the multimoment hydrodynamics equations [5] are applied to solve a problem on flow past a hard sphere at rest at a wide range of values. The direction of evolution of the ground axisymmetric flow [1] after losing its stability is studied. In Section 2, the problem is formulated. In Section 3, the problem is solved for the Stokes flow around a sphere, i.e., at. The Section 4 is devoted to construction of the solution to the multimoment hydrodynamics equations suitable of interpreting the stationary axisymmetric flow at a moderately high values. The Section 5 examines the behavior of the solution that loses its stability after the passage of the first critical Reynolds number value. In Section 6, the characteristic features of the appearance of instability are interpreted in terms of entropy. The principle of retention and loss of the open system stability is formulated. In Section 7, the characteristic features of the development of instability are interpreted in terms of entropy. The evolution criterion is formulated. The Section 8 is devoted to finding the solutions to the multimoment hydrodynamics equations capable of reproducing a vortex shedding. The Section 9 provides an algorithm to select one unstable solution of many found unstable solutions. The selected solution indicates the direction of system evolution. In Section 10, the results are compared with the experimental data.

2. Problem Statement

Consider a space filled with thermodynamically equilibrium gas. Suppose that a solid sphere of radius moves in gas at constant velocity along the axis of the immobile Cartesian frame of reference. Let us now pass from to the Cartesian frame of reference with the axes parallel to those of and the origin made coincident with the center of the moving sphere. In the frame of reference, the sphere is at rest, the inflowing-gas velocity at an infinite distance from the sphere, , is aligned with the positive direction of the axis, and the flow problem is stationary.

The pair distribution function corresponding to the equilibrium state is

(2.1)

where and are the local density of the number of particles and the temperature of unperturbed gas, , , is the mass of a gas particle, and and are the velocity of the center of mass and relative velocity of the pair of particles. In, the pair distribution function can be written as

(2.2)

Let, and be the Cartesian coordinates of point in the space;, and are its spherical coor- dinates, and, Figure 1. Bring the origin of the Cartesian frame of

Figure 1. The frame with its origin at the center of the sphere. The frame with its origin at some point x in the gas.

reference into coincidence with point so that the axis is directed along vector. Mark the projections of vector onto the axes by, and, its projections onto the axes of, by, and, and its spherical coordinates, by , and, Figure 1.

The basic property of pair distribution functions (17) from [3] is that is conserved with time along the trajectory of the center of mass of a pair of particles. In the stationary case, the function remains unchanged along a straight line parallel to vector

(2.3)

The, , and functions are the first integrals of Equation (2.3). These functions were termed in [6] trajectory invariants. In this study, consideration is restricted to gas flows around a sphere, which are invariant under rotation through an arbitrary angle about the axis. Let us compose combinations of trajectory invariants, , and, invariant with respect to this rotation

(2.4)

By virtue of the flow geometry, there are two independent regions of integration with respect to at each point in the gas. The first region, denote it by, is a spherical cone, Figure 2. This region incorporates the trajectories of centers of mass of pairs of particles that originate and terminate at the surface of the sphere. The integration limits for this group of trajectories are. The second region, , embraces all possible trajectories of centers of mass of particle pairs converging to point. The integration limits for the second group of trajectories are unbounded.

Figure 2. Regions of integration with respect to G.

According to the general approach to solving the multimoment hydrodynamic equations in terms of pair functions, outlined in [6] , solution is sought for in the form of an infinite series of the products of trajectory invariants (2.4)

(2.5)

The spatial dependence of (2.5) is controlled by trajectory invariants (2.4), and coefficients are independent of. In keeping with definitions given in [5] , the expressions for the principal hydrodynamic values―the moments of (2.2)―can be written as

,

,

(2.6)

,

, ,

,

,

,

Here, is the local density of the number of particles; is the hydrodynamic velocity;, , , and are the pressure, temperature, stress tensor, and heat flux corresponding to motion of the centers of mass of pairs of particles;, , and correspond to relative motion of particles within pairs;.

Nonprincipal hydrodynamic values and are given by Equation (21) from [3] with the Navier- Stokes accuracy. In terms of spherical coordinates, these values for cylindrically symmetric flows considered in this study (i.e., for flows invariant with respect to rotation through an arbitrary angle about the axis) assume the form

(2.7)

Here, is the coefficient of the dynamic viscosity. From the definitions of gas pressure, temperature, viscous stress tensor, and heat flux [5] it follows that they are linear combinations of moments (2.6) and (2.7)

(2.8)

According to [5] , the overall equations of conservation of the number of particles, momentum, and energy assume the form

(2.9)

(2.10)

(2.11)

The basic property of pair distribution functions (17) from [3] subdivides Equations (2.10) and (2.11). As recast in terms of the spherical coordinates, the and components assume the following form for the cylindrically symmetric flow in question

(2.12)

(2.13)

The boundary conditions must allow for the fact that the contribution of the function given by Equation (2.5) to the hydrodynamic values vanishes at an infinite distance from sphere

(2.14)

At the surface of the sphere, i.e., at, one has to impose the no-leak and no-slip conditions

(2.15)

Heat flux traversing each element of the sphere surface must be counterbalanced by the heat flux within the sphere and the thermal radiation

(2.16)

Here, is the coefficient of the thermal conductivity of the sphere, and is the coefficient of the sphere emissivity.

3. Flow around a Sphere at Small Re

Let us formulate dimensionless parameters, where denotes the Mach number, and the Reynolds number,. Multimoment hydrodynamics Equations (2.9)-(2.13) will be brought to a dimensionless form in a conventional manner [7] . It turns out that the dimensionless equations contain and. Hydrodynamic values can be represented as parametric series

(3.1)

Here, sequentially assumes values;;, , and;

and, and sequentially assumes values;;;. Symbols and appear

over dimensionless quantities. Multimoment hydrodynamics Equations (2.12), (2.13) are written with the Navier-Stokes accuracy, by which token expansion (3.1) contains no terms with.

Consider the flow regime at small Reynolds number. At the hydrodynamic stage of description characterized by small values of the Knudsen number,. Terms of the order of dominate expansion (3.1) of the second-order moments and in the flow regime in question. In expansions (3.1) of hydrodynamic values, , and, it will suffice to retain only, and, respectively. Expansion (3.1) of the particle density flux is reasonably be re-

stricted to, which provides for the desired order of the terms in Equation (2.7) for

. As noted in [6] , the Navier-Stokes approximation is not accurate enough to calculate temperature components and. Thus, at, the spatial dependence of temperature is neglected.

Let us truncate expansion (2.5) to the three lowest-order terms

(3.2)

Each of three terms of series (3.2) can contribute to each principal hydrodynamic value (, , , , , , and). Considering that the principal hydrodynamic values are linearly independent, coefficients, of expansion (3.2) can be written as linear combinations of six arbitrarily chosen components

(3.3)

To calculate the hydrodynamic values (2.6), the first component of series (3.2) has to be integrated with respect to with the appropriate weight function of velocities and within region, and the second and third components, within region. Why is approximated by Equation (3.2), and its components are integrated in such a way is explained below. The integration yields

(3.4)

Here,

(3.5)

The order of coefficients and appearing in Equations (3.4, 3.5) depends on the only density flux component, , retained in expansion (3.1) at. The order of coefficients, , and de-

pends on components and of expansion (3.1). Using one of six linear combinations of

coefficients, we eliminate the contribution to the particle density, proportional to, which is missed in expansion (3.1). Using two linear combinations of coefficients, we eliminate the contribution

to and that is proportional to and is missed in expansion (3.1). Two other linear combina-

tions of enable us to eliminate the contribution to and that is proportional to. The reason why and contain no terms proportional to will be elucidated in Section 4. The last, sixth linear combination of yields in Equation (3.5). Three linear combinations of enable us

to eliminate the contribution to being proportional to and the contribution to and being proportional to, which are missed in Equation (3.1). The first term of series (3.2) con-

tributes only to the first and the third moments of the function. Two linear combinations of offer means of eliminating the contribution to and being proportional to.

Using boundary conditions (2.15) we obtain,. The components of tensor (2.7) can be determined from Equation (3.4). Substituting and from (3.4) into any of the equations of momentum conservation (2.12), we obtain. When analyzing the multimoment hydrody-

namics equations in the order of [6] , it was inferred that, , hence, ,. Substituting the foregoing expressions for the coefficients, , , , into Equation (3.4) and using definitions (2.8), we obtain

(3.6)

(3.7)

Distributions (3.6), (3.7) follow from the classic hydrodynamics equations [7] for the Stokes flow, which are identical to the multimoment hydrodynamics equations in the limit [6] .

The hydrodynamic values (3.6), (3.7) assume the form of a linear combination of the products of by trigonometric functions. To calculate (3.6), (3.7), the first component of Equation (3.2) was integrated with respect to within the region, and the second and third components, within the region. If, on the contrary, the first component is integrated within, and the others, within, Equation (3.4) contain the

products of irrational and trigonometric functions. However, with such products in Equation

(3.4), it is impossible to satisfy boundary conditions (2.15) and equations of momentum conservation (2.12) simultaneously. When the higher-order terms are being retained in expansion (2.5), the functions will appear in Equation (3.4) which are nonlinear in cosθ and sinθ and the Equations (2.15) and (2.12) cannot, as previously, be satisfied simultaneously.

4. Stationary Flow around a Sphere at Moderately High Re

Consider a flow at a moderately high and. In this case, will be retained together with

two other components of expansion (3.1) of the flow density: and, In ad-

dition to already obtained components of expansions (3.1) of the zero- and second-order moments of the pair distribution function, , , and with и, will be taken into account. Expansions (3.1) of the energy fluxes and will be limited to and with и

Expressions for hydrodynamic values (2.6) will be constructed by the same principle as in Section 3―as linear combinations of the products of and trigonometric functions. This principle governs both the structure of the retained terms of expansion (2.5) and the region of integration with respect to for each term of expansion (2.5). In addition to three terms given by Equation (3.2), the following components of expansion (2.5) will be included

(4.1)

(4.2)

Generally, each component of expansions (4.1, 4.2) contributes to each principle hydrodynamic value, irrespective of the order of terms retained in expansion (3.1). Considering that the principle hydrodynamic values are linearly independent, each coefficient, should be represented as series (3.3) of six terms for any

considered order of expansion (3.1). Function given by (4.1) is used to calculate the zero- and the second-order moments. The contribution of is eliminated from the first- and third-order moments

through the use of Equation (3.3), as in Section 3. Function (4.2) is used exclusively to calculate the first- and the third-order moments. The terms of Equation (4.1) proportional to, , and the terms of Equation (4.2) proportional to, are integrated with respect to within region. The terms of Equation (4.1) proportional to and the terms of Equation (4.2) proportional to are integrated with respect to within region.

Distributions of the principle hydrodynamics values (2.6) are given in [8] by Equations (A.1)-(A.7). The expressions (A.1)-(A.7) are written in terms of coefficients. Coefficients are expressed in terms of which appear in Equations (4.1), (4.2). Each coefficient is expanded into series (3.3) of six terms for each considered order of expansion (3.1). In what follows, the relationships between and required to derive explicit expressions for and, will be omitted. Each coefficient, contains dimensionless multiplier.

In order to calculate, one needs, as in Section 3, to satisfy boundary conditions (2.14)-(2.16) and equations of conservation (2.12), (2.13). Imposing boundary conditions (2.15) upon the particle density flux (A.5) from [8] , we express 16 coefficients, in terms of 34 remained coefficients [8] . As revealed in [6] , the -order approximation to the total energy-flux vector is

(4.3)

Here is the coefficient of the thermal conductivity. Substituting the expressions for hydrodynamic values (A.1)-(A.4) from [8] into Equation (4.3), it is possible to express 8 coefficients, , in terms of 26 remained coefficients. Substitute Equations (A.2), (A.5), and (A.6) from [8] into equation of mo-

mentum conservation (2.12) and equate the coefficients of the products for either of the

orders, and. Then 6 coefficients, , can be expressed in terms of 20 remained coefficients, , [8] . Eventually, the expressions for the principle hydrodynamic values (2.6) can be written as

(4.4)

(4.5)

(4.6)

(4.7)

(4.8)

(4.9)

(4.10)

In (4.4)-(4.10), for, for.

The terms proportional to and, were dropped from Equations (4.4)-(4.7) for the zero-and the second-order moments, respectively, because for Equation (2.13) to be satisfied in the order of, all the coefficients of these terms must be zero. For Equation (2.13) to be good in the order of

, one has to put expression (4.9) for From the conditions for no heat flux at the surface of the sphere (at) and at infinite distance from the sphere the expression follows (4.10) for Let us next substitute Equations (4.4)-(4.10), and (3.6), (3.7) into Equation (2.13) and equate the coefficients of products of the order of. It turns out that all but the

and products satisfy the law of energy conservation. These two products satisfy Equation (2.13) at

(4.11)

As previously, Equation (2.13) must be satisfied for each product of the order of individually. Eventually we obtained a nonlinear set of eighteen algebraic equations (A.8-A.10) from [8] . In [8] , of eighteen equations we retained sixteen ones.

The total energy flux equals the heat flux at the surface of the sphere,. Using the solutions of

the internal boundary problem for the Laplace equation, the temperature distribution inside the sphere can be expanded in terms of Legendre polynomials [9] . The coefficients of this expansion can be calculated by matching the temperature distribution (4.4)-(4.7) at the surface of the sphere with the distribution derived from the solution of the internal boundary problem. The resulting balance of heat fluxes at the surface of the sphere (2.16) in order of becomes

(4.12)

Here,.

Sixteen equations (П.8, П.9, П.13) from [8] supplemented with three Equation (4.12), and Equation (4.11) form closed nonlinear set of twenty algebraic equations for twenty coefficients,. Four coefficients, namely, , , and, appear in the distribution of particle density (4.4). Seven coefficients, , , , , , and, are responsible for the distributions of pressure and stress (4.6), (4.7). Seven coefficients, , , , , , , and, appear in the distributions of energy flux (4.9), (4.10), two coefficients, and , govern the distribution of particle-density flux (4.8).

When constructing the distributions of hydrodynamic values in Section 3, the series (2.5) was truncated to the terms that contribute linearly in and to these distributions. Going beyond the limits of the case, we retained seven terms in Equation (4.1) and eleven terms in Equation (4.2). The expansion (4.1) is used to calculate the zero- and the second-order moments and makes the contributions to Equations (4.4)-(4.7) that are proportional to 1, , and. The expansion (4.2) is used to calculate the first- and the third- order moments and gives rise to the -components of Equations (4.8)-(4.10) proportional to and. With higher-order terms of expansions (4.1), (4.2) taken into account, the contributions to Equations (4.4)-(4.10) are proportional to more high powers of. We were compelled to make one exception associated with retention of the trajectory invariants proportional to, , and in Equation (4.2). Owing to these invariants, the distribution (4.8) involves the component proportional to, lacking in the expansions of the other hydrodynamic values (2.6). The terms of expansion (4.2) proportional to, , and enabled us to allow for the terms of the order of, , in the distribution of particle-density flux (4.8). The terms of this order dominate at.

Numerical integration of the nonlinear set of 20 algebraic equations was carried out at, and. Calculations have revealed a great many roots. However, at, the set has only one invariably stable root, , displayed in Figure 4 in [8] . According to the, , solution, an axisymmetric recirculating zone is formed in the wake behind the sphere at Re~20. This recirculating zone has the shape of an axisymmetric toroidal ring. It expands as grows but its shape remains unchanged (Figure 3(а)). At, the system becomes unstable. At, the Barnett corrections [7] are known to be commensurate with the Navier-Stokes terms. For this reason the calculations [8] were interrupted at.

5. First Unstable Flow Regime

As the flow around sphere becomes unstable, the problem becomes nonstationary. As in the stationary case, the pair distribution function is represented as

(5.1)

Here, is given by Equation (2.2), and is given by Equation (2.5) with time-de- pendent coefficients. In going to a stationary flow, pair distribution functions lose their main property (17) from [3] . Retaining, as previously, 21 trajectory invariants in expansion (2.5), we arrive at distributions (A.1)-(A.7) from [8] of hydrodynamic values with time-dependent coefficients,. Coefficients are derived from the condition that hydrodynamic values satisfy boundary conditions (2.14)-(2.16)

(a) (b)

Figure 3. (а) Schematic representation of recirculating zone in near wake behind the sphere; (b) Schematic representation of vortex ring in far wake.

and nonstationary Equations (2.9)-(2.11). Executing the sequence of transformations, which in Section 4 led us to the closed algebraic set of twenty equations, we obtain a closed nonstationary set of twenty equations denoted hereafter by. Sixteen algebraic Equations (A.8), (A.9), and (A.13) from [8] supplemented with three algebraic Equation (4.12), remain unchanged when passing from the stationary problem to the nonstationary one. However, Equation (4.11) changes to

(5.2)

where.

Expressions (4.4)-(4.10) with time-dependent coefficients, , define distributions of principle hydrodynamic values. The set is applied to study the time history of coefficients

. The stationary solution is described in Section 4. Functions give deviations from the stationary solution.

As revealed by numerical integration of the set, solution, , remains asymptotically stable up to certain critical value,. This means that the small deviations from the

stationary solution, , are damped at. The passage of is accompanied by the

, , solution stability loss. Starting with, the, , small axisymmetric deviations begin to increase exponentially at. The solid curve in Figure 4 is the time history of at. The, , deviations grow up to time. At, the solution is cut off. Why the solution to the set terminates at is best explained in terms of trajectories on the phase plane, (Figure 6 in [8] ). It turned out that at, the set becomes unsuitable for modeling evolution of the physical system.

Time history of the coefficient at creates time dependence at the distribution of particle- density flux (4.8). The distribution (4.8) corresponds to observed evolution of the periphery of the recirculating zone in the wake behind a sphere [10] [11] . In accordance with experiment, starting with, and up to, the periphery of the recirculating zone in the wake behind a sphere moves translationally, receding from the sphere. Starting with the periphery of the recirculating zone moves back towards the sphere [10] [11] .

Figure 4. Time evolution of and, Re = 180. Solid curve corresponds to -ensemble, dashed curve corresponds to -en- semble. = 2.333 for -ensemble, = 1.333 for -ensemble.

Suppose that at time, velocities of all gas particles flowing around the sphere reversed their direction, and the boundary conditions were also reversed. Then the periphery of recirculating zone in the wake behind the sphere starts moving back towards the sphere at. As proven in [6] , macroscopic motion of the gas with reversed velocities of its constituent particles is governed by the reverse equations for pair functions. So, observed reverse motion of the periphery of recirculating zone should be described by means of the reverse equations.

However, upon reversing velocities of all gas particles and the boundary conditions, certain details of macroscopic motion of the gas disagree with experiment [10] [11] . For example, upon reversal of gas-particle velocities, the gas involved in vortex motion starts circulating in the opposite direction. However, all gas particles cannot be set in reverse motion at. Nevertheless macroscopic reversal of all gas-particle velocities and boundary conditions is by no means the only way of initiating reverse motion in the wake as a whole. At, reverse motion in the physical system starts spontaneously [10] [11] , i.e., without any interference from outside. Insofar as the reverse equations apply to a vortex propagating backward upon reversal of all gas-particle velocities, it is reasonably assumed that the same equations apply to a vortex set in reverse motion in some other way.

Both the direct and the reverse multimoment hydrodynamics equations are specified by Equations (54)-(56) in [5] . However, expressions for reverse nonprinciple hydrodynamic values and are specified by Equation (2.7), in which each term of the right-hand side has opposite sign [6] .

Executing the sequence of transformations, which led us to the closed nonstationary set, we obtain a closed nonstationary set of twenty equations for coefficients, , denoted hereafter by. Sixteen algebraic Equations (A.8), (A.9), and (A.13) from [8] supplemented with three algebraic equations (4.12) remain unchanged in going from the direct to the reverse nonstationary problem. However, Equation (5.2) changes to

(5.3)

It turned out that the, , solution to the set existed in the neighborhood of the cut-off point. According to [12] [13] , the law of large numbers is violated near singular points (bifurcations, regions of the coexistence of several stable solutions, etc.), and large spontaneous fluctuations may appear in the system. In conformity with [12] [13] , a large spontaneous fluctuation causes the transfer of the system from the cut-off point to the, , solution at time. Starting with, axi-

symmetric deviations, , begin to decay exponentially,.

At time, the solution reaches the neighborhood of the, , stationary solution. At time further evolution of the, , solution finishes. The, , solution exists at, the, , solution exists at,. At any time instant:

(5.4)

A large spontaneous fluctuation causes the transfer of the system from the, , solution

to the, , solution at time. The, , solution is unstable. It follows that

small axisymmetric deviations, , begin to grow starting with time. This pro-

cess is repeated periodically. So, we obtained the solution, , , referred hereinaf-

ter as.

Numerical integration of the, set draws the time history of coefficients,. The

time history of and shown in Figure 4 by a solid line corresponds to initial deviation

. Deviations are dependent quantities and can be calculated from the set, in which equation (5.2) is substituted by equation. Calculation gives,. The equations of the, set were integrated with the accuracy of the order of. Accuracy is a measure of fluctuations of hydrodynamic values at thermodynamic equilibrium. The time history of the coefficients and from Figure 4 corresponds to axisymmetric pulsations of the periphery of recirculating zone in the wake behind a sphere. The half-period of pulsations was estimated at approximately one minute at and.

The dashed curve in Figure 4 is the, dependence calculated at

with accuracy. Thus, the, nonlinear set “draws apart” initially close trajectories. This sensitivity to initial conditions was called the Loretz “butterfly effect” [14] .

The multimoment hydrodynamics equations [5] , as well as the classic hydrodynamic equations, govern space and time evolution of the whole ensemble of systems (Gibbs ensemble) rather than of some individual system. All the microscopic parameters of each individual system are compatible with the initial macroscopic parameters which are present not in the form of particular values but in the form of intervals allowing for their possible fluctuations [15] . Thus, each statistical coefficient, , , is a linear combination of a great many dynamic coefficients, denote their number by K which can be infinitely large

(5.5)

Dynamic coefficients are calculated within the classic mechanics. Fluctuation at any

time is defined as a difference between the dynamic and statistical coefficients,

, or which is the same, as a difference between two arbitrary dynamic coefficients,.

Consider the range of the system parameters within which the equations for statistical coeffi-

cients have stable solution,. Within this range, the overwhelming majority of dynamic coefficients passing in the immediate vicinity of at will remain in the vicinity of indefinitely long. Statistical solution, , describes most of the dynamic trajectories, , with the accuracy.

The situation in the unstable range is radically different. Consider two ensembles of systems. The first -ensemble at incorporates systems with coefficients within

small vicinity near -solution, here,. The se- cond -ensemble at incorporates systems with coefficients within small vicinity near -solution, here,. Let, at time. Time history of the coefficients and is

drawn in Figure 4 (solid and dashed curves). Initially close trajectories diverge. It immediately follows that each dynamic -coefficient, , strictly speaking, behaves in its own way. There is no unique, , which describes any dynamic coefficient from the set, , for the whole ensemble with the accuracy. The Gibbs ensemble disintegrates.

Disintegration of the Gibbs ensemble in the unstable region suggests the follows. The multimoment hydrodynamics Equations (2.9)-(2.11) governing the Gibbs ensemble as a whole are invalid in the region where solutions to these equations become unstable. The and dissipation moments (2.7) are also true for the ensem-

ble as a whole. Deviations, , are nothing but fluctuations of stationary coefficients,

.

Fluctuation of some coefficient, namely, , is preset at, for example,

. Fluctuations, , of other coefficients are calculated from the set, in which Equation (5.2) is substituted by equation. Time evolution of the so for-

mulated fluctuation is approximated by solution of the, set. Fluctuation of any hydrodynamic value at any instant and any point is calculated from Equations (4.4)-(4.10). The fluctuations are interrelated in time and space. In [8] , these interrelated fluctuations have been termed regular. However, there is a factor always present in real physical systems. It is the spontaneous chaotic fluctuation [16] . Spontaneous fluctuations are random independent events. Thus, fluctuation of the statistical coefficient is

(5.6)

After the attainment of the first critical value, the solution to the multimoment hydrodynamics equations loses its stability. The conservation laws (2.9)-(2.11) governing the Gibbs ensemble as a whole become invalid. Regular fluctuations alone cannot provide for fulfillment of Equations (2.9)-(2.11). Strictly speaking, to solve the unstable problem accurately, one needs to switch from the statistical to the dynamic level of description and apply the equations of classical mechanics modeling the dynamics of each individual gas particle. However, numerical integration of the classical mechanics equations for a tremendous number of particles (which is sometimes infinite) is an extremely arduous problem. This line of attack seems ill advised.

In [17] , when modeling an individual system, each hydrodynamic value in the equations of conservation was supplemented with its spontaneous fluctuation component. As was done in [17] , let us attract spontaneous fluctuations. With spontaneous fluctuations taken into account, equations of conservation (2.9)-(2.11) are satisfied in the case of both direct and reverse multimoment hydrodynamics equations. To satisfy the Equations (2.9)-(2.11), the contribution of the time derivative of regular fluctuations must be counterbalanced by the con-

tribution of time and space derivatives of spontaneous fluctuations. Namely

(5.7)

The functions incorporate space derivatives of spontaneous fluctuations. In the case of reverse equations of conservation, besides of Equation (5.7), the change in third equation of set (4.12) (namely, the replacement of by) must be counterbalanced by the contribution of spontaneous fluctuations.

Let us now transform Equation (5.7) to the dimensional form and assess the order of spontaneous fluctuations varying on time and space scales and. Possible large-scale spontaneous fluctuations with and are of the order of large-scale regular fluctuations. Realistic small scale spontaneous fluctuations with and contribute very little, if at all, to the distributions of hydrodynamic values. However, their time and space gradients are of the basic order of magnitude.

6. Interpretation of System Stability Loss in Terms of Pair Entropy

Let us mentally circumscribe a sphere of radius around the sphere of radius. Let us term the gas confined between the surfaces of the coaxial spheres the physical system or simply the system. Let us now turn the radius of the larger sphere to infinity. In pursuance of the condition, the number of particles clinging to the surface of the sphere of radius at any instant is negligibly small compared to the total number of particles in the system. The condition makes it possible to form quasi-isolated system. The processes occurring in a quasi-isolated system are studied without regard for interaction of this system with the ambient medium [16] . According to [16] , let us neglect fluctuations of hydrodynamic values characterizing the system as a whole, which arise due to permeability of the external sphere of radius. In pursuance of the, these fluctuations cannot alter the general physical pattern of the processes.

When deriving equations of entropy conservation (А.6), we reasoned from the concept of a Gibbs ensemble of systems. When modeling an individual system, each hydrodynamic value in the equations of conservation should be supplemented with its fluctuation component [17] . Let us correct Equation (А.9) for the

fluctuation of the pair entropy,. Here, superscript (0) corres-

ponds to the 0-solution to the multimoment hydrodynamics equations. Integrating the resulting equations with respect to over the volume of the system yields

(6.1)

In accordance with (6.1), evolution of the pair entropy is defined by two factors, by the

entropy production in the system and the entropy outflow through the surface confining the system.

Balance Equation (6.1) is valid for and individually. The function characterizes the ensemble of systems as a whole. The second term of the left-hand side of Equation (6.1) is integrated over surface confining volume (over the surface of the spheres of radii and). Note that the physical meaning of the Boltzmann entropy was established without resorting to the concept of a Gibbs ensemble [18] . Inequalities (А.3) were also derived without invoking this concept. Thus, inequalities (А.3) are also valid for, , and, and entropies, , and also convey the meaning of volume occupied by the system in the Г-space.

In terms of Equation (5.6), the fluctuation can be written as

(6.2)

Here, superscripts “R” and “S” mark the contribution of regular and spontaneous fluctuations to the entropy.

Explicit analytical expressions for space distributions of principle hydrodynamic values (4.4)-(4.10) contain dimensionless parameters (and). The -dependence is contained within the, , coefficients. In terms of Equations (4.4-4.10), the components of Equation (6.1) can be written in the form of an infinite series in

(6.3)

Here, , ,. Substituting distributions of hydrodynamic values (3.6), (3.7), (4.4)-(4.10) into (А.7), (А.10) and integrating the resulting expressions with respect to over the volume yield

(6.4)

(6.5)

(6.6)

(6.7)

(6.8)

Here, , , ,. Coefficients

, are given in Section 4. The terms of distributions (4.4)-(4.10) proportional to give rise to terms proportional to and in Equation (6.6) upon integration with respect to over volume. That is why, in deriving (6.6), the integration limit for indefinitely increasing terms is changed by putting.

The function (6.5) incorporates the total number of particles in the system and the total thermal energy. Owing to aforesaid condition and the boundary conditions (2.16) adapted to an individual system, fluctuations of these quantities (and) can be excluded from consideration. Thus,. So, pair entropy accurate to the order of does not change at time. Hence, the system stability is independent of terms of this order. The study of Equation (6.1) undertaken in [19] has revealed that entropy

production in the system (6.7) compensates the entropy outflow through the surface confining the system, then, because the term is

omitted for the reason to be explained later in this section.

In problem with time independent boundary conditions, entropy balance Equation (6.1) accurate to the order of assume the form

(6.9)

As is noted in Section 5, spontaneous fluctuations in large systems tend to be as small as possible. Further, the nearest vicinities of the point at which the solution breaks will be omitted from consideration. In this case, in accordance with the law of large numbers [12] [13] , large spontaneous fluctuations can be disregarded. The contempt for large spontaneous fluctuations makes it possible to omit the, , terms allowing for small spontaneous fluctuations. However, varying with space and time on the scales and, small spontaneous fluctuations keep equations of conservation (6.9) valid. Thus, space and time gradients of small spontaneous fluctuations must be retained on the left hand side of Equation (6.9).

Derivative is calculated from Equation (6.6). One of ways of allowing for the contribution of spontaneous fluctuations is to introduce a random source [17] . However, numerical integration of the, set is as much in error as the random source method. Thus, the contribution of spontaneous fluctuations is taken into account involuntary, owing to the error in numerical calculation.

The study of the spatial variation of the pair entropy [19] has revealed that the

term models entropy removal from the system exclusively through the surface of the

solid sphere of radius. By virtue of the condition, the term models entropy removal by spontaneous fluctuations also through the surface of the sphere. The entropy balance in the solid sphere is maintained by thermal radiation.

Suppose that at time, the system produces an entropy fluctuation. The

values at close to are plotted at Figure 5. As revealed by calculations, at and at. Thus, at, any fluctuation with generated by the system tends to fade out, whereas at, any fluctuation with tends to grow. Time history of is illustrated in Figure 6. From Figure 6 follows that at and. Hence, entropy fluctuations generated by the system die

down monotonically, and the entropy of the system tends to its stationary value. At and,. Thus, entropy fluctuations produced by the system build up

Figure 5. Pair entropy derivative as a function of Re,. The curves correspond to at time. Pair entropy deviation at is taken to be independent of Re:.

Figure 6. Time history of the pair entropy derivative. Curve 1 corresponds to at a subcritical Re value of Re = 129. Curve 2 corresponds to at a supercritical Re value of Re = 130.

monotonically up to point where the solution breaks, and the entropy of the system deviates increasingly farther from its stationary value with increasing. As revealed by calculations, fluctuations always fade out at, and build up at. The function is given by Equation (6.6).

From balance Equation (6.1) follows that at, the pair entropy formed in the system per unit time due to binary collisions,

exceeds the pair entropy removed through the surface confining the system per unit time,

with the result that fluctuations produced by the system die out. At, conversely, the pair entropy outflow exceeds its production, and fluctuations build up. Thus, the case of instability onset in a flow around a sphere is a prevalence of pair entropy outflow over its production at the instant the Reynolds number reaches its critical value.

As revealed by calculations [19] , the Boltzmann entropy behaves exactly as the pair en-

tropy does, when the system passes through. Hence, the cause of instability onset in the system, regardless of whether it is formulated in terms of pair entropy or in terms of Boltzmann entropy, remains the same.

Earlier in the analysis we restricted consideration exclusively to entropy fluctuations with. The reason for this “asymmetry” is as follows. Under the second principle of thermodynamics, the probability of

fluctuations in an isolated system is so much greater than that of

fluctuations that the latter rarely, if ever, occur in nature [16] . As revealed by study [19] , entro-

py fluctuations with direct the system along extremely unlikely, i.e., impractic-

able path. Nevertheless, the probability of fluctuations in an open system is not to be ruled out [19] .

Evolution of fluctuations generated by system depends on two factors, on entropy production and removal through the surface confining the system. The fact that these factors were analyzed without resorting to any kinds of approximations encourages us to believe in the universal nature of the established cause of instability onset. Therefore, the principle according to which an open system retains (or loses) its stability can be formulated as follows.

An open system with time-independent boundary conditions has a stable stationary a-state with entropy while entropy production in it exceeds entropy outflow through the surface confining the system for and does not exceed entropy outflow for

(6.10a)

(6.10b)

As soon as the parameters characterizing the system reach the values, at which at least one of inequalities (6.10а) and (6.10b) fails, the stationary a-state of the open system becomes unstable.

The principle originally formulated for open system with time-independent boundary conditions can be expected to the case of open systems with time-dependent boundary conditions. Generally entropy corresponding to an ensemble of systems may not be reckoned as stationary value. If so, Equation (6.1) no longer implies that the entropy production equals its outflow. That is why, generally, the stability principle is formulated in terms of excess of the entropy production and excess of the entropy outflow.

The a-state with entropy of an open system remains stable while the excess of entropy production generated in the system exceeds its excess of outflow through the surface confining the system for and does not exceed the excess of outflow for

(6.11a)

(6.11b)

As soon as the parameters characterizing the system reach the values, at which at least one of inequalities (6.11а) and (6.11b) fails, the a-state of the open system becomes unstable.

Inequalities (6.10) for systems with time-independent boundary conditions are reduced to inequalities (6.11). However, the stability principle for stationary states (6.10) seems to be more “transparent”. The above formulated stability principle remains invariant in going from pair to Boltzmann entropy [19] .

7. Interpretation of System Evolution in Terms of Pair Entropy

In accordance with the principle of retention and loss of stability (6.11), in an open unstable system, any entropy fluctuation begins to grow. In particular, for a system with time-independent boundary conditions

(7.1)

Based on the expression (7.1), we can formulate the criterion of evolution of an open system with lost stability.

An open unstable system with time-independent boundary conditions, takes a direction of evolution that provides the most rapid decrease in entropy. Namely, of the two directions of development of the instability, having the same values of the entropy and entropy derivative at the time, fluctuations find such a direction that is characterized by lower value of the second derivative of entropy:

(7.2)

wherein

That is, at the time, the system takes the α-direction, for which the second derivative of the entropy with respect to time has a lower value compared to the respectively derivative for the β-direction.

In constructing the approximate solution to the equation for the pair function, only a limited number of terms is retain in the series of products of trajectory invariants (2.5). Different approximate solutions compatible with the boundary conditions of the problem differ in the number of terms retained in expression (2.5). To select the optimal approximate solutions, it is necessary to introduce an additional criterion. The logic of selecting one of the set of approximate solutions can be seen in the formulation of the criterion of evolution (7.2).

Let be the set of parameters characterizing the jth stable stationary solution, , , to the multimoment hydrodynamic equations, and let an increase in be accompanied by the de-

parture of the system state from the state of statistical equilibrium. In the simplest case, the entropy is

calculated in the entire space from one jth solution. Suppose that, at some value of the parameter, the

pair entropy, calculated from the solution, , , occurs in a small vicinity

of the pair entropy. The entropy is calculated based on the solution to the clas-

sic hydrodynamics equations, which are valid in the limit. Based on expression (7.2), let us formulate a criterion for selecting the approximate stable solution for open system.

In interpreting the behavior of open system with time-independent boundary conditions, a solution that provides the fastest drop in entropy should be chosen from the set of stable approximate solutions to the multimoment hydrodynamics equations. Namely, at a certain value of the parameter, an approximate solution with the lowest value of the entropy derivative should be chosen among a few approximate solutions with the same entropy values in a small vicinity of the pair entropy:

(7.3)

This means that, at the given, the priority lies with the kth approximate stable solution to the multimoment hydrodynamics equations for which the derivative of the pair entropy with respect to has the lowest value among the respective derivatives provided by other similar solutions. In more complex cases, it may be necessary to compare the second derivatives of the entropy, as is done in (7.2).

Criteria (7.2), (7.3) are suitable for interpreting open systems with time-dependent boundary conditions.

8. Vortex Shedding Regimes

According to the solution, the recirculating zone is formed in the wake behind a sphere. After the attainment of, the periphery of the recirculating zone begins to pulsate periodically. Pulsating periphery demonstrates the absence of slightest indications of detachment from the core of the recirculating zone. Consequently, there is no vortex street in the far wake behind the sphere. Therefore, the solution does not describe the vortex shedding.

Let the distribution of the particle-density flux be:

(8.1)

To obtain distribution (8.1), the second term on the right side of Equation (4.8) is written in -variables. Distribution (8.1) defines a vortex ring at distance from the sphere center, Figure 3(b). The vortex ring center is located at the axis, which forms the -angle with the -axis, Figure 7. The inflowing gas velocity is aligned with the positive direction of the axis of reference frame XYZ. The XYZ frame has its origin at the center of the sphere.

Let us retain all the terms proportional to in the distributions of principal hydrodynamic values (4.4-4.10). These distributions have the form of linear combinations of the products of the inverse power functions of and trigonometric functions of the polar angle. To bring distribution (8.1) to this form, we change in Equation (8.1) from the variables to the variables, expand these expressions in a Taylor series in powers of, and retain in this series only the terms independent of and linear in:

(8.2)

Then, we change from the spherical coordinates to the spherical coordinates. The variables and are, respectively, the polar and azimuthal angles of the vector in the Cartesian reference frame with the axis defining the direction of the flow incoming on the sphere and in the Cartesian reference frame with the axis forming an angle with the axis, Figure 7. The transformations performed give rise to azimuthal angle dependent terms in the distributions of the particle-density flux:

(8.3)

Here,

Along with the coefficient, the transformed expression (8.3) for the distribution of the particle-density

flux contains two additional coefficients: and. The coefficient characteriz-

es the distance from the vortex ring to the sphere surface, whereas the coefficient characterizes the deviation of the vortex ring center from the axis. Trajectory invariants for expansion (4.2), which lead to the distribution (8.3), are not written in the present Section.

Figure 7. The Z axis of the frame with its origin at the center of the sphere; and are the spherical coordinates of some point x.

The distributions of the principal hydrodynamic values (4.4)-(4.10) associated with the coefficients should be supplemented with distribution (8.3). The presence of -dependent terms in Equation (8.3) for the distribution of the particle-density flux leads to the appearance of these terms in the hydrodynamics equations. Taking into account the dependence on is beyond the accuracy limits of the approximation under consideration.

Expand expression (8.3) for the distribution of the particle-density flux in a Fourier series in the azimuthal angle and substitute the zeroth-order Fourier expansion term in the multimoment hydrodynamics Equations (2.9)-(2.11). Eighteen nonlinear algebraic equations are presented by relationships (A8)-(A10) in [8] . Right- hand sides of Equations (A8)-(A10) from [8] contain twenty coefficients,. In accordance with the algorithm presented in [8] , we supplement these equations by terms containing the coefficients and. Let us supplement eighteen algebraic Equations (A8-A10) from [8] with three algebraic Equation (4.12) and differential Equation (4.2). As a result, we obtain a closed set of nonlinear equations of twenty-second order S22 for the coefficients,. It turned out that, in the investigated range of, of great many solutions to the system S22, only two solutions correspond to such an entropy value that allows these solutions to compete with the solution. We denote these solutions as and.

In the rearrangement of the distribution (8.1), the terms nonlinear in were omitted. Calculations showed that the omitted terms have no significant influence on the solutions and. Substitute the full Fourier expansion of expression (8.3) for the distribution of the particle-density flux into the multimoment hydrodynamics Equations (2.9)-(2.11) and integrate the resulting nonlinear system of differential equations over. In this case, we arrive at system S22 supplemented by several terms. Calculations have shown that these supplements produce no significant influence on the solutions and.

9. Selecting the Direction of Instability Development

Figure 8 shows the time dependence of the dimensionless pair entropy (6.3) calculated from the so-

Figure 8. Time behavior of the pair entropy, calculated within the layer at, Re = 400,.

lution at Re = 400. The pair entropy is calculated in the dominant order,. In

calculating the pair entropy, the spatial integration (A.4) is performed over a semispherical concentric layer, Figure 9. The subscript in brackets corresponds to the limits of spatial integration. The solution describes periodic pulsations of the recirculating zone in the wake behind the sphere. The movement of the representative point over the curve (Figure 8) from to corresponds to the expansion or, in other words, the excitation of the recirculating zone,. The pair entropy will be used as a measure of excitation degree of recirculating zone. Beginning from the time up to the time, the entropy decreases permanently. This behavior of the entropy corresponds to the movement away of the system state that lost its stability from the state of statistical equilibrium. By the time, the degree of excitation of

the recirculating zone reaches a maximum, which corresponds to the minimum value of the function.

At the lettime, the solution to the multimoment hydrodynamics equations breaks. The movement of the representative point over the curve (Figure 8) from to corresponds to the return of the most excited recirculating zone to its original position, i.e., the position corresponding to the time. Since the time, the movement of the representative point is described by the reverse equations of multimoment hydrodynamics [6] . The reverse equations are solved with regressive timing along the time axis, from the past to the future. This timing order is represented on the axis beneath the abscissa in Figure 3 (within

). The pair entropy exists in the range, whereas the pair entropy

(А.15), within. The distribution of the hydrodynamics values involved in the expressions for

and are conjugated with the coefficients and,. Relationships (5.4) between the coefficients and yield for any time within,

. Figure 8 shows that most of the pulsation half-period, approximately, the entropy of the system decreases rather weakly. A sharp decrease in the entropy occurs within a short time of.

The function is a portion of the pair entropy, which specifies the behavior of the en-

tropy in the neighborhood, , of the point of solution break,

(9.1)

Represent the function as a sum of functions:

(9.2)

The entropy is calculated within a semispherical concentric layer (Figure 9). The

Figure 9. Semispherical concentric layer:; semispherical concentric layer:; semis- pherical concentric layer:.

entropy is calculated within a semispherical concentric layer. In calculating

and, the parameters and were set constant.

Let the representative point move upwards along the curve in Figure 8. By the time, the recirculating

zone reaches a degree of excitation characterized by the function. The time corresponds to

the point of intersection of the curves 1 and 2 in Figure 10. At the time, the relations

(9.3a)

(9.3b)

are valid.

The left-hand side of Equation (9.3a) defines the function calculated within the layer.

The function on the right-hand side of (9.3a) is calculated within the layer. In accordance

with definition

(9.4)

The function on the right-hand side of Equation (9.3a) was calculated within the layer

from the solution at Re = 400. The solution describes the motion of a single vortex structure in the

wake behind the sphere. The coefficient in the solution highly accurately defines the position of the vortex structure at the Z axis, since (Figure 7). In accordance with the solution,

the vortex structure monotonically recedes from the surface of the sphere, being accompanied by dissipation of

the vortex structure. This behavior of the vortex structure causes the increase in the pair entropy as grows. In calculating, the initial time for the movement of a single vortex structure, , is

chosen so that, at the time, the vortex structure moving in space would reach the position

where. That is, at the time, the position of a single vortex structure at the

Z axis is determined by, which defines the inner boundary of the layer.

The left-hand side of Equation (9.3) corresponds to the state of the system at which the layer is completely occupied by the recirculating zone, given by the solution. By the time, the recirculating

zone reaches the degree of excitation that is characterized by the function. According to Equa-

tion (9.2), in the state, the core of the recirculating zone has a degree of excitation characterized by the

function, whereas the periphery recirculating zone has a degree of excitation characterized by the function. The right-hand side of Equation (9.3) corresponds to a qualitatively different state

of the system. Namely, in the state, the periphery of the recirculating zone with the degree of excita-

tion, corresponding to the state of the system, is replaced by a single vortex structure in the layer, described by the solution. This vortex structure is characterized by the entropy. In the state, the core of recirculating zone with the degree of excitation corresponding to the state of the system, is replaced in the layer by a recirculating zone core with the different degree of

Figure 10. Time dependence of the entropy derivative, ,. The function calculated within the layer at, Re = 400, is repre- sented by curve 1. Curve 2 is at Re = 400,. The function is calculated in the layer at; the function is calculated in the layer at,. The time of restructuring is. The time count origin for the movement of a single vortex structure is for. The function calculated within the layer at, Re = 400, is represented by curve 3. Curve 4 is at Re = 400,. The function is calculated in the layer at, whereas the function in the layer at,. The restructuring time is. The initial point of time count for the movement of a single vortex structure is at. The function calculated within the layer at , Re = 260, is represented by curve 5. Curve 6 is at Re = 260,. The function calculated within the layer at; the function is calculated in the layer at,. The time of restructuring is. The initial point of time count for the movement of a single vortex structure is at. For curves 1, 2, 3, 4, the abscissa is given at the bottom, and the ordinate, on the left. For curves 5, 6, the abscissa is given at the top, and the ordinate, on the right.

excitation. The degree of excitation of the recirculating zone core in the state is higher as compared to the degree of excitation of the recirculating zone in the state

(9.5)

Thus, expressions (9.3) describe the instantaneous rearrangement of the flow at the time.

The curve 1 in Figure 10 specifies the time evolution of the entropy derivative which is

calculated within the layer. The curve 2 in Figure 10 specifies the time evolution of the derivative

. The function is calculated in the layer. The function is calculated in the layer. In calculating, the variable initial time for the movement of a single vortex structure is chosen at each moment of time so that, by the time, the vortex structure moving in space would reach the position,

where specifies the inner boundary of the layer. For each moment of time, not too close to the time of recovery, there is a time interval, , which ensures the validity of equality (9.3a). However, as follows from Figure 10, only at the time of restructuring, the validity of Equation (9.3a) entails the validity of Equation (9.3b). Then, in accordance with the criterion (7.2), the system takes the direction of evolution set by the flow restructured at the time. The state of the restructured flow is characterized by the terms on the right-hand side of Equation (9.3). Information required to build the curves in Figure 10 and Figure 11 is given in Figures 4- 7 in [20] .

The evolution of the reconstructed flow proceeds as follows. A single vortex structure located at the time in the layer, begins to move downstream in strict accordance with the solution. The inner

boundary of the layer, which is simultaneously the outer boundary of the layer, is bound to the position of a single vortex structure:; i.e., this boundary moves together with the vortex structure at a fixed corresponding to. The time behavior of the pair entropy is represented by the curve 3 in Figure 11. The function is calculated from the solution within the layer with the moving outer boundary. At the restructuring time, a single vortex structure is positioned at. By the time, , the vortex structure reaches the position, in accordance with the solution. Defined by the solution, the outer boundary of the layer at the time sets the value of:. Since the time of recovery, the movement of the system is described by the reverse multimoment hydrodynamics equations [6] , which enable to calculate the time dependence of the reverse pair entropy. The reverse equations are solved with regressive timing direction

along the time axis, from the past to the future. This timing order, , is shown on the axis beneath the abscissa in Figure 11 in the range. The time for progressive timing corresponds to the time for the regressive timing. The curves 3 and 2 in Figure 11 are intersected at a time of. Then,

(9.6)

In Figure 11, the functions and are represented by curves 1 and 2.

Figure 11. Time behavior of the pair entropy, Re = 400. The restructuring time is; the time of detachment is, ,. The pair entropy, , exist in the range. The pair entropy, , exist in the range,. Curve 1 represents the functions and calculated within the layer at. Curve 2 represents the functions and calculated within the layer at. Curve 3 represents the functions and calculated within the layer. The movement of the outer boundary of the layer is given by,.

In accordance with Equations (9.4b) and (5.4), and the definition of the pair entropy, at the time (equivalent to), we have

(9.7а)

(9.7b)

Thus, by the time, the system reaches conditions identical to those the system had at the time; i.e., the single vortex structure leaves the layer. The layer is entirely occupied by the recirculating zone described by the solution. By the time, the single vortex structure crosses the external border of the periphery of the recirculating zone,. We assume that the time is the moment of separation of the single vortex structure from the recirculating zone.

Further, as previously, the periphery of the recirculating zone located in the layer is replaced by a single vortex structure given by the solution. At the same time, in the layer, the recirculating zone core with the degree of excitation corresponding to the state is replaced by the core with the degree of excitation corresponding to the state. At the time, if the equality (9.3a) is held, so is equality (9.3b). Then, in accordance with the criterion (7.2), the direction taken after the replacement is more preferable for the system. As at the time, at, a single vortex structure, in accordance with the solution, starts to move from the outer boundary of the recirculating zone core, , to the outer boundary of the recirculating zone periphery,. The periodic shedding of a single vortex structure occurs with a period of . The positions of shedding vortex structures are schematically indicated by vertical lines in the lower branch in Figure 12.

It turns out that, for equality (9.3b) to be held, the degree of excitation of the recirculating zone core

should increase with decreasing a parameter that defines the position of the vortex structure at the periphery of the recirculating zone at the time of restructuring. At, the degree of excitation of the recirculating zone core reaches its maximum,

since at, the time coincides with the time of recovery,. A further decrease in the parameter excludes the possibility of flow restructuring, i.e., the possibility of emergence of a vortex structure at the periphery of the recirculating zone.

The time behavior of the pair entropy derivatives and

is shown in Figure 10. Curve 3 corresponds to the

state of the system before restructuring, , curve 4 corresponds to the state of the restructured system,. The curve 3 in Figure 10 is positioned slightly below curve 1. However, in the scale of Figure 10, curve 1 and 3 are indistinguishable. The time of restructuring corresponds to the point of intersection of curves 3 and 4 in Figure 10. At the time of restructuring, equalities (9.3a) and (9.3b) are held simultaneously. In accordance with the criterion (7.2), the system takes the direction of evolution that is predetermined by the flow restructured at the time.

Figure 12. Positions z of periodically shedding vortex structures at the Z axis. The leftmost value is the position of the vortex in the recirculating zone at the restructuring time. The lower branch:, shedding period is, Re = 400; the middle branch:, , Re = 400; the upper branch:, , Re = 260, the next vortex structure at the upper branch occupies the position z = 10.3.

The use of the extremely low value gives extremely low values of the entropy derivatives

and at the time of restructuring. The values of these derivatives at correspond to the intersection of curves 3 and 4 in Figure 10. The values of and at the point of intersection of curves 3 and 4 are lower than the respective values at the point of intersection of curves 1 and 2, for which. As a result, the derivative at and

has an extremely low value. Therefore, in accordance with the criterion (7.2), the extremely low value of which defines the position of a single vortex structure in the recirculating zone at the time of restructuring is preferable.

The use of extremely low value of leads to a situation where the representative point takes a minimum time to move from position of break to the position of separation. The choice of an extremely low value of leads to an extremely low value of the period of separation of a vortex structure from the recirculating zone,. In Figure 12, the vertical lines on the middle branch indicate the spatial position of periodically shedding vortex structures for an extremely low value of. Shedding vortex structures are located extremely close to each other (middle branch).

The time behavior of the pair entropy derivatives is shown in Figure 10. Curve 5 corresponds to the function, which is calculated from the solution at Re = 260. The curve 6 corresponds to the function, in which the de-

rivative is calculated from the solution at Re = 260. The values of the parameters used in the calculations, , are specified in the Figure captions.

The curve 5 corresponds to the state of the system before restructuring, , while the curve 6 corresponds to the state of the restructured system,. The time of restructuring corresponds to the point of intersection of curves 5 and 6 in Figure 10. At the time of restructuring, the equalities (9.3a) and (9.3b) are fulfilled simultaneously. In accordance with the criterion (7.2), the system takes the direction of evolution that is set by the flow restructured at the time.

It is found that, with increasing parameter, defining the position of the vortex structure in the recirculating zone at the time of restructuring, the entropy of the recirculating zone periphery grows much faster than the entropy of single vortex structure. At an extremely high value of the parameter, the difference is not so great as to make the equali-

ties (9.3a) and (9.3b) simultaneously valid. The extremely high value of the parameter at Re = 260 is given in the Figure captions. A further increase in the parameter excludes the possibility of flow restructuring, i.e., the possibility of emergence of a vortex structure in the recirculating zone. The extremely high value of leads to the extremely low value of the period of shedding of single vortex structures from the recirculating zone core. The extremely low value of provides extremely low values of the entropy derivatives at the time of restructuring, which are represented by the point of intersection of curves 5 and 6 in Figure 10. Therefore, in accordance with criterion (7.2), the extremely high value of the parameter is selected from the spectrum of its values. The upper branch in Figure 12 shows the spatial positions of periodically shedding vortex structures for the extremely high values of, marked by vertical lines.

The calculations have revealed the dependence of the domain of existence of the solution on Re. At Re = 260, the solution exists at, here the variable highly accurately defines the posi-

tion of the vortex structure at the Z axis. In this range of values, the pair entropy is not

small enough to ensure that the relations (2.3) will be held. So, at Re = 260, the solution is the only solution allowing the vortex shedding. At Re = 400, each of the solutions and allows the vortex shedding. Each of the solutions and is characterized by its own time of flow restructuring,. Calculated from the solution, the period of shedding of single vortex structures from the recirculating zone core is significantly smaller than the corresponding period calculated by the solution. The entropy

derivatives, and

at the point of restructuring for the solution are consi-

derably higher than the respective values represented by curves 3 and 4 in Figure 10 at their point of intersection. Thus, in accordance with the criterion (7.2), the solution is more preferable than the solution at Re = 400.

10. Comparison with Experiment

All the experiments with flows past a sphere available at present are divided into two groups in the review [1] . The first group includes experiments in which a body sinks freely in a medium under the action of gravitation. Besides, the first group includes the experiments in which a hard sphere is fixed rigidly in the air or water tunnel. The second group experiments were made by towing a body through a bath filled with a liquid. The reason for such a division is as follows. The velocity of a free-stream flow can be considered as strictly uniform in none of the first group experiments. A strictly uniform profile of the flow running against a body can only be achieved in experiments of the second group [10] [11] by using a rigid support system.

According to the measurement data reported in [10] [11] , the stationary recirculating zone formed in the near wake behind a sphere at Re~10 - 20 increases as grows. After the attainment of the first critical Reynolds number, the periphery of the recirculating zone in the near wake behind a sphere begins to pulsate periodically. The recirculating zone remains toroidal during pulsations, and its forefront is firmly fixed on the sphere. The pulsating flow remains axisymmetric.

The pulsations become increasingly well defined as Re increases, and their amplitude grows. After the passage of, the periphery of the recirculating zone begins to be periodically shed from its core and moves downstream. The shed vortex structure has the shape of ring. Separate vortex rings depart from a sphere downstream and move along the spiral path [11] . The attainment of is accompanied by the change in the regime of vortex shedding from sphere. The frequency, which characterizes the vortex shedding, monotonically increases as Re grows starting with. The increase in the frequency causes the complete disappearance of intervals between vortex rings. Vortex rings penetrate into each other and form the continuous spiral sheet in the wake behind a sphere [11] . According to the measurement data reported in [10] ,. In [11] , the, the and values were obtained. The continuous spiral vortex sheet was observed in [11] over the whole range of Reynolds numbers studied, up to Re = 30,000.

The simplest solution to the multimoment hydrodynamic equations gives the distribution (3.6-3.7). This distribution coincides with the Stokes solution to the classic hydrodynamics equations in the limit. The stationary axisymmetric solution to the S20 set of the multimoment hydrodynamic equations satisfactorily reproduces the stable flow around a sphere. The attainment of the first critical Reynolds number is accompanied by the solution stability loss. The cause for stability loss is as follows. After the passage of, the entropy outflow through the surface confining the system begins to exceed the entropy production in the system. Such interpretation corresponds to the principle of retention and loss of the open system stability (6.11). The unstable axisymmetric solution reproduces satisfactorily the first unstable flow regime.

After the attainment of the second critical Reynolds number value, the solution for the periphery of the recirculating zone and in the far wake is replaced by the solution to the S22 set, which describes a vortex ring moving downstream. The cause for the replacement is that the combination of the solutions and becomes more preferable in comparison with the solution. The combination of the solutions and provides a sharper drop in the entropy in the course of evolution than the solution does (Figure 10). The criterion (7.2) dictates the choice of the direction of development that is given by combination of the and solutions. The positions of detached vortices at the Z axis at given by the solution, are shown in the upper branch in Figure 12. According to the solution, there are significant intervals between neighboring vortex rings. The combination of the solutions and describes the vortex street.

After the attainment of the third critical Reynolds number value, the solution at the periphery of the recirculating zone and in the far wake behind the sphere is replaced by the solution to the S22 set, which also describes a vortex ring moving downstream. The cause for the replacement is that the combination of the solutions and provides a sharper decrease in the entropy in the course of evolution than the combination of the solutions and does (Figure 10). Criterion (7.2) dictates the choice of the direction of development that is given by combination of the and solutions. The positions of detached vortices at the Z axis at, given by the solution, are shown in the middle branch of Figure 12. The combination of the solutions and describes the vortex sheet.

The observed process of formation and separation of the vortex structure from the recirculating zone in the wake behind the sphere is shown in photographs 40 and 41 in [21] . A vortex structure arising near the surface of the sphere moves inside the core zone of the recirculating zone to its periphery, increasing in size and acquiring a definite shape. At the periphery of the recirculating zone, the vortex structure is separated from this zone and moves downstream.

In the calculated flow pattern, a formed vortex structure appears in the recirculating zone instantly at the time. This simplification of the observed complex process of formation and separation of the vortex structure is motivated by a lack of computational resources for simulating the process of vortex shedding.

Eighteen trajectory invariants are used to generate the functions (5.1), (5.2). These invariants give the distributions of the principle hydrodynamic values (4.4)-(4.10) coupled with coefficients,. Only two components coupled with the coefficients and form the distribution of the particle-density flux (4.8). Moreover, only one of these components coupled with the coefficient dominates the particle density flux (4.8), reproducing recirculating zone in the near wake behind the sphere. The drag coefficient of the sphere also depends only on one component coupled with:. Moreover, the same coefficient controls the pressure drop at a surface of the sphere at points and:

.

Further, the dominant component in the distribution (4.8) that is coupled with the coefficient is used to reproduce a single downstream-moving vortex structure unrelated to the sphere, expression (8.1). Thus, a detailed description of the observed pattern of nucleation and separation of the vortex structure can be achieved by increasing the number of trajectory invariants involved in calculating the particle-density flux distribution. This increase is necessary both in the formulation of the set S20 describing the pulsations of the recirculating zone in the near wake and in the formulation of the set S22 describing the motion of single vortex structure outside the recirculating zone. Moreover, constructing the set of equations suitable for describing the far wake only (system S22) is barely an approximation of the calculation procedure. A further improvement necessitates combining sets S20 and S22. The solutions to the combined set must describe the movement of the single vortex structure not only in the far wake, but also in the recirculating zone. It is these solutions to the combined set of equations that must compete with the solutions to the set S20.

Forming the set of equations S20 and S22 was conducted so that these sets would be able to reproduce the time behavior of the recirculating zone and the evolution of the vortex structure separated from it in the wake of the sphere. The terms of expansion (4.2) proportional to, , and enabled us to describe the wake behind a sphere. By virtue of the proportional to, , and terms, however, yet another physically improbable vortex configurations arise in front of the sphere at its surface at a sufficiently high Re [8] . The appearance of physically improbable vortex configurations in the distribution of particle density flux (4.8) generally does not allow of using this distribution outside the wake. To obtain the distributions of the hydrodynamic values outside the wake, it is necessary to formulate a different set of the multimoment hydrodynamics equations. The boundary separating solutions to different sets should be sought based on the principle (7.3). In other words, the boundary separating the solutions should be positioned in space so as to ensure the most rapid decrease of the entropy.

In the Section 6, in calculating the pair entropy, the solution to the set S20 was used both in the wake behind the sphere and outside the wake. At close to, physically improbable vortex configurations in front of the sphere are very weak. That is why, these configurations distort the calculation very little at. Another result of going with the solution trace beyond the wake produced divergent terms in the course of spatial integration over the volume V. The mergence of divergent terms forced us to restrict the integration limits in the Section 6 to a spherical concentric layer directly adjacent to the surface of the sphere. Namely, in deriving Equation (6.6), the integration limit for indefinitely increasing terms is changed by putting. As revealed by calculations, the appropriate numerical coefficients in Equation (6.6) change under changes in the integration limit, however, characteristic features of variation of

and remains unchanged.

The approximate calculation of the pair entropy in Section 9 was performed within hemispherical concentric layers, , and (Figure 9). This approximation was used because the dominant contribution to the deviation of the entropy from its equilibrium value comes from processes in the near wake behind the sphere. These are the processes of excitation of the recirculating zone and of separation of the vortex structure that appeared in this zone. The outer boundary of the layer is specified by the position of the separating vortex structure at the time of its detachment, , from the recirculating zone. Estimates show that the approximation of the boundaries of the near wake behind the sphere by the boundaries of the hemispherical layer does not lead to a qualitative distortion of the results.

11. Discussion

Let us substitute the inequality (6.11a) into Equation (6.1) that is written in terms. We obtain that the - state of the open system remains stable if

(11.1a)

or, in terms of Boltzmann entropy

(11.1b)

The negative definite value serves as Lyapunov function, and then, the first of inequalities (11.1b) is a sufficient condition for the system stability in re to fluctuations. According to the Lyapunov theorem, the condition (11.1b) is compatible with the fundamental Lyapunov definition of stability. Inequality (11.1b) makes up a formal mathematical stability criterion. A formal mathematical criterion, however, gives no way of revealing the physical cause of stability or instability of an open system.

Behind the modern concepts of open system stability is the Glansdorff-Prigogine theory [12] . This theory, in turn, is formulated in terms of nonequilibrium thermodynamics and, therefore, allows for the possibility of constructing the universal function with non-positive fluctuation at any fluctuations of hydrodynamic values

The fluctuation of was used as Lyapunov function. The set of stability conditions derived in [12] in terms of this fluctuation involves inequality (1.11b) and

(11.2)

Based on inequalities (11.1b) and (11.2), P. Glansdorff and I. Prigogine formulated the condition for stability of open system state in terms of entropy production and entropy outflow [12] . They argue that the main quality controlling the stability and evolution of the system is excess of entropy production. The Glansdorff-Prigogine stability condition states that for the system state to be stable, the entropy excess must be generated in the system rather than being absorbed

(11.3)

In accordance with Glansdorff-Prigogine theory [12] , for the -state to be stable, the effect of excess of entropy outflow on the system must be equivalent to the effect of excess of entropy production. In other words, for the state to be stable, entropy excess must flow into the system through the confining surface, thus augmenting the entropy excess generated in the system

(11.4)

Inequalities (11.3) and (11.4) supplemented with similar inequalities for the fluctuation of the function (11.2) incorporating the terms of the order responsible for system stability make up the Glansdorff- Prigogine stability criterion [12] . Therefore, in [12] , the additional constraints (11.2) and (11.4) are imposed on stability condition (11.3).

The methods of kinetic theory and hydrodynamics amplified in this study offer a much more effective means of constructing Lyapunov function than nonequilibrium thermodynamics. In particular, there was no need in imposing constraint (11.2) on fluctuation of the function. Having made sure that a positive definite fluctuation set the system on a very unlike unrealistic path, we bounded the domain of definition of by condition. In the bounded domain of definition, , as revealed by calculation, is negative definite and, therefore, can be used as a Lyapunov function. In the Glansdorff -Prigogine theory, excess of entropy outflow through the surface confining the system is treated as a quantity of secondary importance. The Glansdorff-Prigogine stability criterion was formulated, assuming that excess of entropy outflow exerted a strictly limited effect on the system, identical to the effect of excess of entropy production. The central idea of novel concept of stability is that the excess of entropy production and excess of entropy outflow are equally important. According to principle proposed in this study, the factor controlling the stability of the system is a balance of excess of entropy production and excess of entropy outflow rather than one of these quantities alone.

Let us verify that stability condition (6.11) is more general than the Glansdorff-Prigogine criterion. Indeed, if inequalities (11.3) and (11.4) are true, condition (11.1b) is surely met. If the reverse is the case (inequalities (11.3) and (11.4) fail), the system is not necessarily unstable. The excesses of entropy production and calculated with due regard for the terms of the order of responsible for the system stability are plotted against Re in Figure 13. The function is given by Equa-

tion (6.8), the function is defined in [19] . From Figure 13 follows that,

at. Thus, although the Glansdorff-Prigogine criterion breaks, the system is stable by virtue of relationship (6.11a) between entropy production and entropy outflow. Therefore, the Glansdorff-Prigogine stability criterion generally is incapable of interpreting hydrodynamic instability. More detailed comparison of the criterion (6.11) with the Glansdorff-Prigogine criterion was carried out in [19] .

To derive the evolution criterion, the authors of [12] set up a balance equation for certain part of the total entropy of the system and differentiated it with respect to time. The equation takes into account variation of the entropy production solely under variation of generalized forces with time. The formulation of the evolution criterion was made possible by the neglect of the time differential of the excess of entropy outflow through the confining surface. There is generally no reason to disregard the excess of entropy outflow, and therefore, the Glansdorff-Prigogine evolution criterion is not universally true. The criterion states that evolution of an open system with time independent boundary conditions follows the path along which the entropy production in the system has no tendency for an increase with time

(11.5)

The time history of of the dominating order of is shown in Figure 14. The excess of pair entropy is given by Equation (6.7). As revealed by calculation, the system

develops in accord with inequality (11.5) at Re = 129 and. At Re = 130 and (the covered Re values did not exceed 400) the Glansdorff-Prigogine evolution criterion is invalid. So, the property (11.5) is valid when system relaxes to the stable state. However, system evolution after stability loss is not accompanied by fulfillment of the property (11.5). The above formulated conclusion remains invariant in going from terms of pair entropy to terms of Boltzmann entropy.

The need to use the Glansdorff-Prigogine criterion of evolution on a search of solutions to the classic hydrodynamics equations is missing. The criterion demonstrates the fulfillment (or absence of the fulfillment) of the property (11.5) for the found solution to the hydrodynamics equations. The function of the criterion (7.2) is very significant. Criterion (7.2) chooses a direction of evolution for the system. i.e., after stability loss, the calculation of hydrodynamic values is impossible without the attraction of the criterion (7.2).

Figure 13. The excess entropy production component of the order of responsible for instability as a function of Re. Curve 1 cor- responds to the excess of the pair entropy production, and curve 2, the excess of Boltzmann entropy production . Entropy fluctuation at time is assumed to be independent of Re:; .

Figure 14. Time history of the excess entropy production component of the dominating order of: curve 1 corresponds to a subcritical Re value of Re = 129, and curve 2, to supercritical Re value of Re = 130.

Since the time of Boltzmann, the responsibility for directing the evolution of the system rests with the initial conditions realized in the system, namely the set of initial values of the coordinates and velocities of all particles [22] [23] . For a given mutual arrangement of the particles, the system evolves in the direction that we see everywhere and every second. However, there are such arrangements of particles that direct the system in an extremely unlikely, rarely realized direction. Even weak disturbance of the configuration of particles can change the direction of evolution (is characteristic size of particles) [24] .

The local pair entropy corresponding to the direct equations for pair distribution functions (1.12) from [6] and the multimoment hydrodynamics equations they yield can only be produced in the system due to binary collisions at any space point and at any instant,. Thus, at any instant binary collisions merely rise the pair entropy of the system,. Such behavior of the entropy is in full accordance with the second law of thermodynamics [25] . The solutions to the direct multimoment hydrodynamics equations describe the direction of evolution of the system that is everywhere and every second is found in nature.

The local pair entropy corresponding to the reverse equations for the pair distribution functions (1.15) from [6] and the reverse multimoment hydrodynamics equations they yield can only be absorbed in the system due to binary collisions at any space point and at any instant,. Thus, at any instant binary collisions absorb the pair entropy of the system,. The solutions to the reverse multimoment hydrodynamics equations describe the evolution of the system in the opposite direction, which, as is commonly believed, is extremely rare in nature.

The direct multimoment hydrodynamics equations are valid for the progressive direction of timing on the time axis pointing from the past to the future. The reverse multimoment hydrodynamics equations are valid for multimoment regressive direction of timing on the same time axis [6] . Processes occurring in nature are objective events, while the choice of the direction of timing on the time axis is subjective process. Time is counted by an observer, while processes occurring in nature absolutely insensitive to the direction in which the observer counts the time. Let two observers, agreeing upon the origin, began to observe some phenomenon. Let the first observer counts time in the regressive direction at the time axis, directed from the past to the future and let the first observer managed to describe the phenomenon he observed using the reverse equations. Let now the second observer counts time in the progressive direction at the same time axis. Then, the second observer can claim that the reverse equations described phenomena for progressive direction of timing.

This means the following. On finding the solution to the inverse equations, the first observer obtained the distribution of all hydrodynamic variables and their spatial and temporal derivatives. The first observer found agreement between the calculated and measured values. The second observer used the solution to the reverse equations to compare with his observations and, on establishing the appropriate relationship between the two time scales, found that the reverse distributions agree with the observed values in their time scale. Changing the sign of the time derivatives, the second observer also obtained an agreement with his observations. However, the second observer failed to derive direct equations that would satisfy altered distribution.

This example returns us to the Boltzmann time, when there were heated debates about the correctness of the Boltzmann equation and the H-theorem. Boltzmann opponents, E. Zermelo and J. Loschmidt, gave examples of processes that are not described by the Boltzmann equation [23] [26] . Basically, L. Boltzmann acknowledged the objections of opponents, agreeing with the existence of processes that can not be governed by his theory. However, L. Boltzmann argued that such processes was extremely unlikely, i.e., impracticable [22] . Probably, at a weak deviation of the system state from the state of thermodynamic equilibrium, as predicted L. Boltzmann, conditions guided a system in the unlikely direction are extremely rare. However, as the degree of nonequilibrium increases the probability of occurrence of such conditions may grow. The penetration into the instability field confirmed this assumption. It turned out that the motion of unstable system along the unlikely direction becomes an event that repeats periodically.

12. Conclusions

The experiment records two stable stationary medium states represented by the and velocity distributions, and a stable state of the central type with the velocity distribution. Each of these three stable flows begins to develop in its own direction qualitatively different from other flows when it loses stability. The development occurs through a sequence of regular nonstationary periodic states schematically shown in Figure 1 from [3] . Each of these three experimentally observed directions inevitably reaches the periodic vortex shedding mode. Vortex shedding along each of three directions is characterized by its own characteristic features intrinsic in it. However, irrespective of the direction selected experimentally, periodic vortex shedding is obligatory, well defined, and fairly prolonged along Reynolds numbers mode of the development of a turbulent process. Experiment records six vortex shedding modes, and one pulsation mode along the three directions of instability development. The development of the ground flow occurs through the pulsating regime, and further, through the vortex shedding regimes and.

The direct numerical integration of Navier-Stokes equations in the problem on a flow around a solid sphere at rest was performed by various numerical methods. Nevertheless, the results of all these numerical experiments were absolutely identical (see review [1] ). Calculations find two stable stationary solutions, and and a nonstationary stable solution,. Apart from the, , and solutions, the Navier-Stokes equations only have a multiperiodic, that is, essentially chaotic, solution. The, , and solutions satisfactorily reproduce, , and stable flows.

According to calculations, the development of instability occurs in strict correspondence to the classic Landau-Hopf scenario [14] . After first critical Reynolds number value is reached, the solution loses stability and bifurcates to the solution (regular bifurcation). The solution, after it loses stability, bifurcates to the limiting cycle (the Hopf bifurcation). After next critical Reynolds number value is reached, the limiting cycle loses stability and is substituted by the chaotic solution. It follows that, according to the Landau-Hopf scenario, the system, after it loses stability, inevitably reaches a new stable position and exercises either periodic or chaotic motion about it [27] . So, calculations cannot put anything in correspondence to seven of ten experimentally observed modes schematically shown in Figure 1 from [3] . In particular, to periodic pulsation of the recirculating zone not accompanied by vortex shedding, , and two periodic vortex shedding modes, and.

The ideas of experiment on the nature of the phenomenon of vortex shedding fundamentally disagree with the concepts of classic hydrodynamics [14] . In accordance with experiment, a system, when loses stability after the attainment of a certain critical Reynolds number value, remains further in unstable state. One unstable regime is replaced by another unstable regime as Re grows. Experiment interprets the observed phenomenon of vortex shedding as the von Karman instability, or as the Kelvin-Helmholtz instability. By definition [28] , the instability occurs only during the flow transformation, i.e., the time intervals when one solution loses its stability and another stable solution is saturated. The process of saturation of new stable state is extremely limited in time and non-periodic. i.e., the bulk of time the system is at a stable state. Time-limited non-periodic process of saturation can in no way correlate with the strictly periodic vortex shedding phenomenon, which, generally, is not limited in time. Thus, the solutions to the classic hydrodynamics equations are incapable of corresponding to ideas of experiment on unstable nature of the vortex shedding phenomenon. Classic hydrodynamics can put only stable solutions in correspondence with observed vortex shedding modes. In particular, in problem on flow around a sphere the limiting cycle is the only case where it is possible to correlate the observed vortex shedding from a sphere to calculation result. However, numerous attempts to find a vortex street in the far wake behind a body failed [1] [2] [21] .

The multimoment hydrodynamics confirms the ideas of experiment on unstable nature of the phenomenon of vortex shedding. The crossing of the first critical Reynolds number value is accompanied by the stability loss. The system loses its stability when entropy produced in the system can not compensate entropy outflow through the surface confining the system. Such interpretation follows directly from the principle of retention and loss of the open system stability formulated in Section 6. In accordance with solutions to the multimoment hydrodynamics equations, the system, when loses its stability, remains further unstable. One unstable flow is replaced by another unstable flow as Re grows. The replacement of one unstable flow regime by another unstable regime is governed the tendency of the system to discover the fastest path to depart from the state of statistical equilibrium. This striving follows directly from the evolution criterion formulated in Section 7.

Figure 8 demonstrates the behavior of pair entropy of the system losing its stability. According to the solution, pair entropy describes periodic pulsations of the recirculating zone in the wake behind a sphere. Starting with, pair entropy begins to decay monotonously. The monotonous decay is replaced by surge. However, the sharp decrease of entropy is not unlimited. At, the decrease of entropy finishes. Starting with, pair entropy begins to increase up to time. At, the system returns to its original position. Pair entropy begins to decay again starting with. The process is repeated with the period .

A similar behavior demonstrates pair entropy, which is associated with periodic vortex shedding, Figure 9. The sharp decrease of entropy is completed at. Starting with, pair entropy begins to increase up to time. At, the system reaches the position that it occupied at. Pair entropy begins to decay again starting with. The process is repeated with the period.

The cut-off of the solution is the reason for the completion of the sharp decrease of pair entropy at time. It turned out that at the multimoment hydrodynamics equations [5] become unsuitable for modeling evolution of the system. However, in the neighborhood of the cut-off point there is a solution to the reverse multimoment hydrodynamics equations [6] . The solution to the reverse equations changes the direction of evolution, striving the system to the state of statistical equilibrium. The pair entropy corresponding to the direct multimoment hydrodynamics equations can only be produced in the system due to binary collisions at each time point. Such behavior of the entropy is in full accordance with the second law of thermodynamics. The pair entropy corresponding to the reverse multimoment hydrodynamics equations can only be absorbed in the system due to binary collisions at each time point. The solutions to the reverse multimoment hydrodynamics equations describe the evolution of the system in the opposite direction, which, as is commonly believed, is extremely rare in nature.

At his time, Boltzmann, defending his point of view in disputes with opponents, suggested that the exclusive conditions that guide the system in a highly unlikely direction arise very rarely. Apparently, Boltzmann’s assumption is correct for a weak deviation of the system from the state of statistical equilibrium. However, after crossing the border of the instability field, exclusive conditions arise with periodic regularity. This regularity manifests itself in each of three unstable modes that represent the flows, , and.

The tendency of an unstable physical system to find the fastest path to recede from the state of statistical equilibrium does not lead the system to infinity. At time, the entropy stops decreasing. The multimoment hydrodynamics equations are unable to provide solutions that would continue to divert the system from the state of statistical equilibrium. Time intervals during which the system moves away from the state of statistical equilibrium are periodically followed by intervals within which the system tends to equilibrium. It is this periodicity that permits to interpret vortex shedding, a graphic example of periodic unstable phenomena.

References

- Lebed, I.V. and Umanskii, S.Y. (2007) The Appearance and Development of Turbulence in a Flow Past a Sphere: Problems and the Existing Approaches to Their Solution. Russian Journal of Physical Chemistry B, 1, 52-73.
- Lebed, I.V. and Umanskii, S.Y. (2012) On the Possibility of Improving Classic Hydrodynamics Equations by an In- crease in the Number of Hydrodynamic Values. Russian Journal of Physical Chemistry B, 6, 149-162.
- Lebed, I.V. (2013) About the Prospects for Passage to Instability. Open Journal of Fluid Dynamics, 3, 214-229. http://dx.doi.org/10.4236/ojfd.2013.33027
- Lebed, I.V. (2002) On the Inapplicability of Navier-Stokes Equations to Interpreting the Turbulence. Physica A, 315, 228-235. http://dx.doi.org/10.1016/S0378-4371(02)01254-2
- Lebed, I.V. (1996) Method of Two-Particle Distribution Functions. Hydrodynamic Equations. Chemical Physics Re- ports, 15, 861-883.
- Lebed, I.V. (1996) Hydrodynamic Equations Stemming from Two Particle Distributions in the Limit of Weak Non- Equilibrium. Analysis of Invertibility of Equations. Chemical Physics Reports, 15, 1725-1750.
- Loitsyanskii, L.G. (1966) Mechanics of Liquids and Gases. Pergamon, Oxford.
- Lebed, I.V. (1997) The Method of Pair Functions as Applied to the Problem of a Flow around a Quiescent Solid Sphere. Chemical Physics Reports, 16, 1263-1301.
- Tikhonov, A.N. and Samarskii, A.A. (1953) Equations of Mathematical Physics. Gostechizdat, Moscow.
- Taneda, S. (1956) Experimental Investigation of the Wake behind a Sphere at Low Reynolds Numbers. Journal of the Physical Society of Japan, 11, 1104-1108. http://dx.doi.org/10.1143/JPSJ.11.1104
- Chomaz, J.M., Bonneton, P. and Hopfinger, E.J. (1993) The Structure of the Near Wake of a Sphere Moving Horizon- tally in a Stratified Fluid. Journal of Fluid Mechanics, 254, 1-21. http://dx.doi.org/10.1017/S0022112093002009
- Glansdorff, P. and Prigogine, I. (1971) Thermodynamic Theory of Structure, Stability and Fluctuations. Wiley-Inter- science, A Division of John Wiley & Sons, Ltd., London.
- Prigogine, I. (1980) From Being to Becoming: Time and Complexity in the Physical Sciences. W.H. Freeman and Company, San Francisco.
- Schuster, H.G. (1984) Deterministic Chaos. Physik_Verlag, Weinheim.
- Zubarev, D.N. (1971) Nonequilibrium Statistical Thremodynamics. Nauka, Moscow.
- Landau, L.D. and Lifshitz, E.M. (1976) Course of Theoretical Physics. Vol. 5, Statistical Physics, Part 1, Nauka, Mos- cow.
- Lifshitz, E.M. and Pitaevskii, L.P. (1980) Course of Theoretical Physics. Vol. 9, Statistical Physics, Part 2, Pergamon, New York.
- Ferziger, J.H. and Kaper, H.G. (1972) Mathematical Theory of Transport Processes in Gases. North-Holland Publish- ing Company, Amsterdam.
- Lebed, I.V. (1998) About the Behavior of the Entropy of a Gas Flow Losing Its Stability. Chemical Physics Reports, 17, 411-439.
- Lebed, I.V. (2014) Development of Instability in the Problem of Flow around a Sphere. Russian Journal of Physical Chemistry B, 8, 240-253.
- Johnson, T.A. and Patel, V.S. (1999) Flow Past a Sphere Up to a Reynolds Number of 300. Journal of Fluid Mechan- ics, 378, 19-70. http://dx.doi.org/10.1017/S0022112098003206
- Boltzmann, L. (1896) Entgegnung auf die Wärmetheoretischen Betrachtungen des Hrn. Zermelo. Ann. Phys. Chem., 57, 773-784.
- Zermelo, E. (1896) Uber einen Satz der Dynamik and die Mechanische Wärmetheorie. Ann. Phys. Chem., 57, 485-494.
- Orban, J. and Bellemance, A. (1967) Velocity-Inversion and Irreversibility in a Dilute Gas of Hard Discs. Physics Let- ters A, 24, 620-621. http://dx.doi.org/10.1016/0375-9601(67)90651-2
- Liboff, R.L. (1969) Introduction to the Theory of Kinetic Equations. Willey, New York.
- Loschmidt, J. (1876) Uber den Zustand des Wärmegleichgewichtes eines Systems von Körpern mit Rücksicht auf die Schwerkraft. Wien. Ber., 73, 128-142.
- Tomboulides, A.G. and Orszag, S.A. (2000) Numerical Investigation of Transitional and Weak Turbulent Flow Past a Sphere. Journal of Fluid Mechanics, 416, 45-73. http://dx.doi.org/10.1017/S0022112000008880
- Hannemann, K. and Oertel Jr., H. (1989) Numerical Simulation of the Absolutely and Convectively Unstable Wake. Journal of Fluid Mechanics, 199, 55-88. http://dx.doi.org/10.1017/S0022112089000297
- Isihara, A. (1971) Statistical Physics. Academic Press, New York.
- Grad, H. (1949) About Kinetic Theory of Rarefied Gases. Communications on Pure and Applied Mathematics, 2, 331- 407. http://dx.doi.org/10.1002/cpa.3160020403

Appendix

Suppose that and are the pair and paired-particle distribution functions corresponding to approaching particles [5] . The correlation between and a one-particle distribution function was established in [5] using

(А.1)

Then from the Gibbs inequality [29] we obtain

(А.2)

Let us set up the and functions as follows

(A.3)

Integrals with respect to velocity in expressions (А.1, А.2) are taken over a three-dimensional volume.

(A.4)

Spatial integration (А.4) is carried out over the physical volume of the system.

(А.5)

Here, is the Boltzmann entropy. By virtue of inequalities (А.2) we obtain

(A.6)

The term in inequality (А.6) brings and to the same dimensionality, and is a total number of particles in the system. Thus, is the lower bound of both and. It is well known that Boltzmann entropy conveys the meaning of a volume occupied by the system in the Г-space

[18] . The possibility of comparing, , and using relationships (A.6) implies that

and have the same physical meaning as the Boltzmann entropy. The physical meaning of and can also be elucidated in a conventional way [18] . Functions, , and approximate the total entropy of the system to different accuracies. If is termed the pair entropy of the system, is its local pair entropy.

When calculating, is expanded in terms of a double set of Hermitian polynomials (see Equation (29) in [5] ). Furthermore, let be equal to approximations (31, 34) with coefficients (57) from [5] , composed on seven principle hydrodynamic values (, , , , , , and). The so composed function is then substituted into Equation (А.3), and the resulting expression for is truncated to the first few non-vanishing terms

(А.7)

Following the Grad’s procedure for calculating [30] , let us expend into a series in terms of Hermitian polynomials. The hydrodynamic values appearing as coefficients in this expansion are reasonably presented in the form of linear combinations of G- and v-components, as was done in [5] . Substituting the so formulated function into Equation (А.5) yields the expression for, in which only the first non-vanishing terms are retained

(А.8)

In analyzing the evolution of the pair entropy with space and time we shall proceed from Equation

(28b) given in [5] . Multiplying this equation by and integrating the result with respect to velocities yield

(А.9)

The expressions for the components of the left-hand side of Equation (A.9) are written in [19] . The expression for the pair entropy production in binary collisions assumes the form

(А.10а)

(А.10b)

Inequality (А.10b) implies that is non-negative at any space point and at any instant. Thus, binary collisions merely raise the pair entropy of the system. Since the Boltzmann’s time, this result constitutes the essence of the Н-theorem.

The space and time evolution of the Boltzmann entropy will be inferred from Equation (A.2b) derived in [5] . Integrating Equation (A.2b) from [5] with respect to and taking advantage of Equations (18b) and (19.b) from [5] , we obtain the equation for with the collision integral unrelated to the Boltz- mann hypothesis of molecular chaos. Let us multiply the resulting equation by and integrate it with respect to. We obtain

(А.11)

The expressions for the components of the left-hand side of Equation (A.11) are written in [19] . The expression for the Boltzmann entropy production in binary collisions assumes the form

(А.12)

The study undertaken in [6] has revealed that in some special cases it is possible to establish linear relations between and and between and. In particular, in the limit [6] . In the limit, the second and the third terms of the right-hand side of Equation (А.12) are negligibly small compared to the first term. In this case

(А.13)

It does not generally follow from Equation (А.12) that is non-negative value. The property

(A.10) of pair entropy is the main reason why the pair entropy should be preferred to Boltzmann en-

tropy when interpreting non-equilibrium phenomena. The tendency of the Boltzmann entropy of a gas system to build up in binary collisions was inferred from the Boltzmann equation based on the molecular chaos hypothesis [18] . If the molecular chaos hypothesis is abandoned, this tendency generally disappears. From Equation (А.13) follows that binary particle collisions produce Boltzmann entropy at any space point and at any instant, only if gas is close to statistical equilibrium.

Integrating Equation (А.9) with respect to over the volume of the system yields

(А.14)

The second term of the left-hand side of Equation (А.14) is integrated over surface confining volume.

In virtue of inequality (А.10b), the right-hand side of Equation (А.14) is non-negative. Pair entropy is defined in terms of pair distribution function labeled by app superscript. In the derivation of Equation (А.9), the multimoment hydrodynamics Equations (54)-(56) from [5] supplemented with direct expressions (53) from [5] , are used. Equation (А.14) is valid only in the direction of increasing time. The positive direction of the time axis runs from the past to the future [6] .

Pair entropy is defined in terms of pair distribution function labeled by div superscript. In analyzing the evolution of the pair entropy with space and time we shall proceed from Equation (28a) given in [5] . We obtain

(А.15)

In the derivation of Equation (А.15), the multimoment hydrodynamics Equations (54)-(56) from [5] supplemented with reverse expressions (1.17) from [6] , are used. The expression for is given by Equation (А.7). The expression for differs from its counterpart (the expression (А.10a) for) by the sign. That is why, at any space point and at any instant the function is non-posi- tive. i.e., binary particle collisions absorb the pair entropy at any space point and at any instant, and right-hand side of Equation (А.15) is non-positive. Equation (А.14) is valid only in the direction of decreasing time. The positive direction of the time axis runs from the past to the future [6] .