Molecular Origins of Turbulence

: The twin problems of closure and dissipation have been barriers to the analytical solution of the Navier–Stokes equation for fluid flow by top-down methods for two centuries. Here, the statistical multifractal analysis of airborne observations is used to argue that bottom-up approaches based on the dynamic behaviour of the basic constituent particles are necessary. Contrasts among differing systems will yield scale invariant turbulence, but not with universal analytical solutions to the Navier– Stokes equation. The small number of publications regarding a molecular origin for turbulence are briefly considered. Research approaches using suitable observations are recommended.


Introduction
Turbulence has been used to describe motion in many systems where a fluid flow is evident.Water and air were historically the most apparent, but it has been applied to systems as large as giant astronomical molecular clouds, to galaxies, to stars, to magnetohydrodynamic systems such as solar winds, to the atmosphere, to the ocean, to rivers and to smaller systems such as those that exist in industrial facilities and in laboratories, in superconducting fluids and in biological cells.
The basic argument used here is that rather than taking a top-down approach, a bottom-up approach should be adopted.That implies knowing what the basic constituent particles are and using their dynamics to transition upscale to a flow.
Although Leonardo da Vinci employed the word 'turbulenza' to describe the flow of water in a stream, the initial successful attempt at the mathematical formulation of fluid flow was made by Navier [1].Stokes later amended the formulation to produce the well-known Navier-Stokes equation [2].At the time, the existence and behaviour of atoms and molecules was not known, so a top-down approach was unavoidable.Since we are concerned with meteorology, we concentrate on air, a gas composed of molecules, noting that the hydrogen-bonded structures in liquid water make the ocean different than the atmosphere, although both are turbulent.
The format of this article will be a text with embedded references, but introducing a small number of key equations, equivalences and figures located in the Discussion section below.This perspective can be short because few authors have considered a bottom-up approach to turbulence, particularly in the atmosphere.

Brief Historical Progression
The experiments of Reynolds [3] showed a transition from a laminar to a turbulent flow in viscous liquids forced down a pipe by a pressure gradient.That transition was the subject of much research, including the thesis work of Heisenberg [4].Heisenberg's theory supervisor, Sommerfeld, gave him the highest grade, while his experimental supervisor Wien gave him the lowest.A compromise gave Heisenberg a mediocre pass, whereupon he left to work for Born in Göttingen, remarking that turbulence was too difficult, a fortunate circumstance for the development of quantum theory.Chapman and Enskog approached the transport properties of gases from an atomic and molecular standpoint [5,6].Their work did not extend to fluid flow.The fourth edition of Lamb's book [7] introduced the term 'vorticity' to specify the spin of a local fluid element about its axis, a concept originated by Helmholtz [8].Research starting with the Navier-Stokes equation continued for half a century, but even with numerical simulation, twin difficulties arose, from dissipation and from the closure problem, and the fact that the expression for the n th moment of a variable involves its (n + 1) th power.Chapman and Cowling [9] utilized the solution of the Maxwell-Boltzmann equations to derive transport functions for non-uniform gases.Alder and Wainwright [10] initiated the study of fluids via the computer integration of the motion of hard, elastic spheres, molecular dynamics.In the J. G. Kirkwood memorial issue of The Journal of Chemical Physics, Grad [11] considered the need for a molecular approach to gas dynamics.Alder and Wainwright continued their work [12,13] and showed, in a crucial development, that hydrodynamic behaviour emerged, as 'ring currents', in a thermalized population subject to a symmetry-breaking flux of energetic 'molecules' represented as hard elastic spheres.Vorticity was generated nonlinearly at the smallest possible scale and could propagate upscale in a self-sustaining manner.
The development of statistical multifractality for application to the atmosphere by Schertzer and Lovejoy was a significant event [14][15][16].Their formulation of generalized scale invariance demonstrated that the formal equivalence between equilibrium statistical thermodynamic variables and their scaling equivalents in the Hamiltonian flux dynamics by a Legendre transform was more than a mathematical coincidence, it was a mapping [16].When applied to the largely 'horizontal' airborne data collected in the lower stratosphere between 1987 and 2004, and the 'vertical' data collected via a dropsonde through the winter troposphere over the northeastern Pacific Ocean 2004Ocean , 2005Ocean and 2006 [17] [17], it vindicated the theory and yielded new results about atmospheric dimensionality and scaling [18,19].
Muriel [20,21] argued for a molecular origin of turbulence, even claiming an exact solution for the Navier-Stokes equation [22], a claim that did not find favour with the Clay Institute.The final sentence of the abstract reads "No turbulence is obtained from the solution".Abramov [23] has argued that the rotational motion of molecules can cause turbulence in gases via the effect on intermolecular potentials.There have been calls to measure the rotational state of O 2 and N 2 in the atmosphere [17][18][19].

A Mechanism
We pursue the view that hydrodynamics has to be formulated from the bottom scales up, which means that the basic constituent particles, which are the photons, atoms, molecules, electrons, ions, or nucleons, are the starting point.In turn, that may imply there is no unique formulation of turbulence.One size does not fit all, the basic particles and their interaction with boundary conditions offer differing possibilities for different systems.Scale invariance is a feature common to many differing turbulent systems.
The application of Lovejoy and Schertzer's statistical multifractal approach [16] to airborne observations [24,25] led to the conclusion that the correlation of the intermittency of temperature with the ozone photodissociation rate implied that local thermodynamic equilibrium was not attained in the lower stratosphere.An explanation was sought in Alder and Wainwright's molecular dynamics calculations that hydrodynamic behaviours ('ring currents') emerged on very short timescales in a thermalized population of Maxwellian hard spheres subjected to a symmetry-breaking flux of such high-speed particles [13].The statistical multifractal approach also offers scaling versions of entropy and Gibbs free energy [25], and when combined with the Langevin equation, points to a mechanism for turbulence in the air [18,19].The persistence of molecular velocity after collision prevents instantaneous thermalization and breaks the continuous translational symmetry of a Maxwell-Boltzmann equilibrated gas.Rotational symmetries in real molecules [23,26], N 2 and O 2 , will also break the continuous translational symmetry of any thermalized air.This view of the generation of turbulence implies effects on radiative transfer via spectral line shapes and, hence, dissipation via radiation to space.It implies that turbulence is associated with the operation of the Second Law of Thermodynamics.The Gibbs free energy enables the work that allows the emergent turbulence to propagate nonlinearly upscale.It will possess the directional variation that is observed in the turbulent structure and that is expected from the nature of molecular rotation and collisions.Laminar flow is unlikely to emerge from such conditions.The real atmosphere has more energy than numerical models which assume local thermodynamic equilibrium, as illustrated in Figure 3 of reference [19].
An unexpected result arising from the application of statistical multifractality to airborne data was the correlation between the intermittency of temperature with the ozone photodissociation rate and with temperature itself [24,25].The excited photofragments were not thermalized immediately, but instead, the persistence of their velocity after collision [9] induced vorticity via the Alder-Wainwright mechanism [13], the nonlinearity of which enabled its propagation upscale in the form of turbulence.
The importance of the boundary conditions is evidenced by comparing the atmosphere's turbulence with that in giant molecular clouds, which are vastly larger, much colder, much more dilute and in which gravity plays a different role [27,28].Nevertheless, fractality is evident and characterizes the observed turbulence.Its intermittency, the concentration of energy in local structures, is evident both astronomically and in Earth's atmosphere [29,30].The scaling exponents in the two milieux show that the contrasting physical processes produce different manifestations of turbulence.It is unrealistic to expect the integration of the Navier-Stokes equation for fluid motion to produce a universal solution for turbulence, given the differing mechanisms of dissipation.
The Langevin equation has been used to interpret the molecular dynamical approach to atmospheric turbulence [18].The emergent hydrodynamic behaviour exhibited by Alder and Wainwright [13] is viewed as organization sustained by the most energetic molecules, while the maintenance of an operational temperature by the majority of the near-average molecules is dissipation.That is the reverse of the meteorological convention that means are organized while eddies are dissipative.The most energetic molecules in the probability distribution carry the Gibbs free energy that provides the work that sustains the flow, while the more numerous, average molecules define an operational temperature, which does not, however, equal the theoretical temperature of the equilibrium of the Maxwell-Boltzmann distribution [18,19].
Symmetry breaking is an important factor in the mechanism of turbulence and its generation [18,19,31].Symmetry is also involved in the effect of boundary conditions.When studying laminar flow in a cylindrical pipe, Couette flow in a channel, the flow over an aerofoil, or over the planetary surface the solid boundaries impose a symmetry on the fluid that it does not naturally have.The transition from a laminar to a turbulent flow is a manifestation of that, and of the lack of sensitivity and resolution in the examining instruments.All angles occur in the scattering from molecular collisions, but are rarely considered in experimental or theoretical studies of turbulence.The numerical modelling and diagnostic approaches that employ such techniques as Fourier analysis and large eddy simulation seek to impose symmetries on the air that it does not have, respectively, sine waves and cubes.

Discussion
The discussion will centre around two equations, a table and three key figures which have been generated in analysing atmospheric observations [9,13,18,24,25].
Chapman and Cowling [5] formulated an equation describing the persistence of molecular velocity after collision: The persistence ratio ϖ 12 is the ratio of the mean velocity after collision to what it was prior to the collision between molecules of masses m 1 and m 2 .If m 1 = m 2 , then ϖ 12 = 0.406, but in the general case for m 1 ̸ = m 2 , the more massive molecule will be decelerated less than the lighter one.
This process, when applied to the Alder and Wainwright mechanism [13], will break the continuous translational symmetry of a sample of an equilibrated Maxwell-Boltzmann gas, especially that of an inhomogeneous gas like air, which is never isotropic or at equilibrium.Alder and Wainwright discovered, via a molecular dynamics calculation, that hydrodynamic behaviour was induced in an equilibrated population of 'billiards' on very short time and space scales when subjected to a symmetry-breaking energetic flux of such elastic spheres.See Figure 1.These 'ring currents' are what a meteorologist would call vortices; the generation of vorticity thus occurs on the smallest scales and exhibits scaling behaviour [32]-it is emergent turbulence.By virtue of the nonlinearity of the mechanism, it can propagate upscale.The closure problem is automatically avoided.Dissipation is automatically included via the action of the energetically near-average molecules in defining an operational temperature, which, in the atmosphere, provide dissipation via radiation to space.It is why large-scale global models can successfully integrate the Navier-Stokes equation for weather forecasting-provided they are continually fed with observations.In free-running mode, they will have difficulties [19].The absence of local thermodynamic equilibrium also affects, in principle, the line shapes of the infrared active molecules-carbon dioxide, water vapour and its dimer, nitrous oxide, methane, ozone, CFCs, HCFCs-via collisional effects.
Meteorology 2024, 3, FOR PEER REVIEW 4 The persistence ratio 12 is the ratio of the mean velocity after collision to what it was prior to the collision between molecules of masses m1 and m2.If m1 = m2, then 12 = 0.406, but in the general case for m1  m2, the more massive molecule will be decelerated less than the lighter one.
This process, when applied to the Alder and Wainwright mechanism [13], will break the continuous translational symmetry of a sample of an equilibrated Maxwell-Boltzmann gas, especially that of an inhomogeneous gas like air, which is never isotropic or at equilibrium.Alder and Wainwright discovered, via a molecular dynamics calculation, that hydrodynamic behaviour was induced in an equilibrated population of 'billiards' on very short time and space scales when subjected to a symmetry-breaking energetic flux of such elastic spheres.See Figure 1.These 'ring currents' are what a meteorologist would call vortices; the generation of vorticity thus occurs on the smallest scales and exhibits scaling behaviour [32]-it is emergent turbulence.By virtue of the nonlinearity of the mechanism, it can propagate upscale.The closure problem is automatically avoided.Dissipation is automatically included via the action of the energetically near-average molecules in defining an operational temperature, which, in the atmosphere, provide dissipation via radiation to space.It is why large-scale global models can successfully integrate the Navier-Stokes equation for weather forecasting-provided they are continually fed with observations.In free-running mode, they will have difficulties [19].The absence of local thermodynamic equilibrium also affects, in principle, the line shapes of the infrared active molecules-carbon dioxide, water vapour and its dimer, nitrous oxide, methane, ozone, CFCs, HCFCs-via collisional effects.The fluctuating variability in atmospheric observations [17][18][19] led to the application of the generalized scale invariance pioneered by Schertzer and Lovejoy [14,15] and the statistical multifractality of Lovejoy and Schertzer's work [16].The results are summarized in Table 1.The equivalences between statistical thermodynamical and scale invariant quantities are not merely formal but are mappings via a Legendre transform between the energy in the former and its Hamiltonian flux dynamics in the latter.
Table 1.The equivalence between statistical thermodynamic and scaling variables.

Variable
Statistical Thermodynamics Scaling Equivalent The variables are obtained as follows.q defines the qth order structure function of the observed quantity.The scaling exponent K(q) is derived from the slope of a log-log plot [25].Equation ( 2) expresses the relation of the Hurst exponent H to the Gibbs free energy equivalent.
The examination of energy E in terms of a scale ratio produces an expression for the fractal co-dimension c(γ).C 1 is the co-dimension of the mean, characterizing the intensity of the intermittency.The Lévy exponent α characterizes the generator of the intermittency, which is the logarithm of the turbulent flux; for a real system, its value may not be confined to the theoretical range.Means converge but the variance does not, a result expressed as 2.0 > α >1.5.
The calculation of the intermittency of temperature C 1 (T) is exemplified in Figure 2. The implication is that local thermodynamic equilibrium is not obtained in the lower stratosphere, possibly extending to the entire atmosphere above and below [19].The correlation between the intermittency of temperature C 1 (T) with both ozone photodissociation rate J[O 3 ] and with the temperature itself is shown in Figure 3.The scale invariance associated with turbulence is revealed and enables that condition [14][15][16][17][18][19].Note that any change in dynamical phase, say, between the lower stratosphere and the troposphere below or the mesosphere above, would show up as a change in the slope of the log-log plot.
The way this mechanism is envisaged in the stratosphere is that the translationally hot molecules created by ozone photodissociation pile up a higher-number density ahead of themselves, leaving lower-velocity molecules at lower number density behind.The average molecules move to eradicate the gradient, creating ring currents as they do so.The highernumber density ahead of the fast molecules will experience higher ozone photodissociation rates, so enabling the acceleration and propagation of the emergent hydrodynamic flow.This accounts for the correlation seen in Figure 3.At the same time, the action of the more numerous, average molecules ensures thermodynamic behaviour by establishing an operational temperature.Those average molecules will emit infrared radiation appropriate to their state and, in doing so, ensure that the Second Law of Thermodynamics is obeyed by the atmosphere.
The statistical multifractal analysis also applies to atmospheric aerosols as well as to air molecules [33].Such particles span the mesoscale gap between microscopic and macroscopic scales [34].Turbulent flow is probably manifested both two-dimensionally in the surface film and three-dimensionally in the interior.The graph from which C1(T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).The graph from which C 1 (T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).The graph from which C1(T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).Future research should address the translational and rotational state of the N 2 and O 2 molecules, from the surface to the mesosphere.In situ observations over the same range would assess whether statistical multifractal analysis would provide scaling exponents capable of addressing the cold bias evident in many free-running models of the atmosphere.If N 2 and O 2 are not instantly equilibrated, their rotation will ensure that the directionality of the emergent hydrodynamic flow will vary through the whole 360 o , as frequently observed in turbulence.The real atmosphere will have more energy than models that assume local thermodynamic equilibrium.
The thermodynamic form of multifractality is discussed in [16], Chapter 5.1, and the scaling from atomic scale-up is examined in [32].A recent perspective can be found in [35], and for hydromagnetic systems, in [36].An examination of Stokes's Law at different molecular scales has appeared very recently [37].

Conclusions
Turbulence is a state arising from the dynamics of the basic constituent particles of a fluid medium.In the atmosphere, that is nitrogen and oxygen molecules subject to an influx of solar photons.In other systems composed differently and where gravity and temperature interact differently, scale invariant turbulence is still manifested, but not with the same scaling exponents.The turbulence characterized in this way is argued to be consistent with the statistical multifractal formulation of energy, entropy and Gibbs free energy and, hence, is involved in the action of the Second Law of Thermodynamics.In principle, the Navier-Stokes equation can be solved from the bottom up, but not from the top down.The problems of closure and dissipation are automatically avoided.

Figure 1 .
Figure 1.The emergence of ring currents, vortices, in a thermalized population of 220 Maxwellian elastic atomic spheres subjected to a symmetry-breaking flux of energetic atoms.The molecular dynamics simulation was originally performed by Alder and Wainwright [13].The blue arrows are

Figure 1 .
Figure 1.The emergence of ring currents, vortices, in a thermalized population of 220 Maxwellian elastic atomic spheres subjected to a symmetry-breaking flux of energetic atoms.The molecular dynamics simulation was originally performed by Alder and Wainwright[13].The blue arrows are averages of the atom vectors after 9.9 collisions, the black arrows are from a Navier-Stokes equation integration.Later simulations showed disagreements between the two approaches.

Figure 2 .
Figure 2. (Upper): The temperature record observed from the ER-2 during POLARIS at approximately 55 mbar, 19,970,506 (yyyymmdd format) in low wind speed conditions.(Lower left): The log-log graph from which H(T) was calculated, corresponding with the 5/9 theoretical value of statistical multifractality.(Lower right):The graph from which C1(T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).

Figure 3 .
Figure 3.These results are from Arctic flights [19].J[O3] is averaged along the flight segment concerned and is shown on the ordinate of the left plot, with T averaged in the same way and displayed on the right plot.The abscissa on both shows the intermittency of T, calculated from the relevant lower right graph in Figure 2. The vertical and horizontal bars are standard deviations.The nonconvergence of variance and the changes of air mass within the time periods account for the scatter; nevertheless, the effect is clear and makes physical sense [17-19,24].

Figure 2 .
Figure 2. (Upper): The temperature record observed from the ER-2 during POLARIS at approximately 55 mbar, 19970506 (yyyymmdd format) in low wind speed conditions.(Lower left): The log-log graph from which H(T) was calculated, corresponding with the 5/9 theoretical value of statistical multifractality.(Lower right):The graph from which C 1 (T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).

Meteorology 2024, 3 , 6 Figure 2 .
Figure 2. (Upper): The temperature record observed from the ER-2 during POLARIS at approximately 55 mbar, 19,970,506 (yyyymmdd format) in low wind speed conditions.(Lower left): The log-log graph from which H(T) was calculated, corresponding with the 5/9 theoretical value of statistical multifractality.(Lower right):The graph from which C1(T), the intermittency of temperature, was calculated using scaling variables K(q) and q, see Equation (2).

Figure 3 .
Figure 3.These results are from Arctic flights [19].J[O3] is averaged along the flight segment concerned and is shown on the ordinate of the left plot, with T averaged in the same way and displayed on the right plot.The abscissa on both shows the intermittency of T, calculated from the relevant lower right graph in Figure 2. The vertical and horizontal bars are standard deviations.The nonconvergence of variance and the changes of air mass within the time periods account for the scatter; nevertheless, the effect is clear and makes physical sense [17-19,24].

Figure 3 .
Figure 3.These results are from Arctic flights [19].J[O 3 ] is averaged along the flight segment concerned and is shown on the ordinate of the left plot, with T averaged in the same way and displayed on the right plot.The abscissa on both shows the intermittency of T, calculated from the relevant lower right graph in Figure 2. The vertical and horizontal bars are standard deviations.The nonconvergence of variance and the changes of air mass within the time periods account for the scatter; nevertheless, the effect is clear and makes physical sense [17-19,24].