BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
W. D. McComb contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
From a critical review of Onsager's pioneering paper which was published in 1949 we conclude that his analysis was at odds with other workers in the field because instead of taking the infinite Reynolds number limit he simply set the viscosity equal to zero which is not the same thing.
A new form of the well known Lin equation for the energy balance in isotropic turbulence is given which takes account of the filtered-partitioned transfer spectra which were used to resolve the scale paradox. It is argued that this new formulation provides a new starting point for both theoretical and numerical methods which may lead to a better understanding of the underlying energy processes in turbulence.
A model for the Reynolds number dependence of the dimensionless dissipation rate $C_{\varepsilon}$ was derived from the dimensionless Kármán-Howarth equation, resulting in $C_{\varepsilon}=C_{\varepsilon, \infty} + C/R_L + O(1/R_L^2)$, where $R_L$ is the integral scale Reynolds number. The coefficients $C$ and $C_{\varepsilon,\infty}$ arise from asymptotic expansions of the dimensionless second- and third-order structure functions. This theoretical work was supplemented by direct numerical simulations (DNSs) of forced isotropic turbulence for integral scale Reynolds numbers up to $R_L=5875$ ($R_λ=435$), which were used to establish that the decay of dimensionless dissipation with increasing Reynolds number took the form of a power law $R_L^n$ with exponent value $n = -1.000\pm 0.009$, and that this decay of $C_{\varepsilon}$ was actually due to the increase in the Taylor surrogate $U^3/L$. The model equation was fitted to data from the DNS which resulted in the value $C=18.9\pm 1.3$ and in an asymptotic value for $C_\varepsilon$ in the infinite Reynolds number limit of $C_{\varepsilon,\infty} = 0.468 \pm 0.006$.
We observe a symmetry-breaking transition from a turbulent to a self-organized state in direct numerical simulation of the Navier-Stokes equation at very low Reynolds number. In this self-organized state the kinetic energy is contained only in modes at the lowest resolved wavenumber, the skewness vanishes, and visualization of the flows shows a lack of small-scale structure, with the vorticity and velocity vectors becoming aligned (a Beltrami flow).
The pseudospectral method, in conjunction with a new technique for obtaining scaling exponents $ζ_n$ from the structure functions $S_n(r)$, is presented as an alternative to the extended self-similarity (ESS) method and the use of generalized structure functions. We propose plotting the ratio $|S_n(r)/S_3(r)|$ against the separation $r$ in accordance with a standard technique for analysing experimental data. This method differs from the ESS technique, which plots $S_n(r)$ against $S_3(r)$, with the assumption $S_3(r) \sim r$. Using our method for the particular case of $S_2(r)$ we obtain the new result that the exponent $ζ_2$ decreases as the Taylor-Reynolds number increases, with $ζ_2 \to 0.679 \pm 0.013$ as $R_λ \to \infty$. This supports the idea of finite-viscosity corrections to the K41 prediction for $S_2$, and is the opposite of the result obtained by ESS. The pseudospectral method also permits the forcing to be taken into account exactly through the calculation of the energy input in real space from the work spectrum of the stirring forces.
The low wavenumber expansion of the energy spectrum takes the well known form: $ E(k,t) = E_2(t) k^2 + E_4(t) k^4 + ... $, where the coefficients are weighted integrals against the correlation function $C(r,t)$. We show that expressing $E(k,t)$ in terms of the longitudinal correlation function $f(r,t)$ immediately yields $E_2(t)=0$ by cancellation. We verify that the same result is obtained using the correlation function $C(r,t)$, provided only that $f(r,t)$ falls off faster than $r^{-3}$ at large values of $r$. As power-law forms are widely studied for the purpose of establishing bounds, we consider the family of model correlations $f(r,t)=α_n(t)r^{-n}$, for positive integer $n$, at large values of the separation $r$. We find that for the special case $n=3$, the relationship connecting $f(r,t)$ and $C(r,t)$ becomes indeterminate, and (exceptionally) $E_2 \neq 0$, but that this solution is unphysical in that the viscous term in the Kármán-Howarth equation vanishes. Lastly, we show that $E_4(t)$ is independent of time, without needing to assume the exponential decrease of correlation functions at large distances.
Procedures for time-ordering the covariance function, as given in a previous paper (K. Kiyani and W.D. McComb Phys. Rev. E 70, 066303 (2004)), are extended and used to show that the response function associated at second order with the Kraichnan-Wyld perturbation series can be determined by a local (in wavenumber) energy balance. These time-ordering procedures also allow the two-time formulation to be reduced to time-independent form by means of exponential approximations and it is verified that the response equation does not have an infra-red divergence at infinite Reynolds number. Lastly, single-time Markovianised closure equations (stated in the previous paper above) are derived and shown to be compatible with the Kolmogorov distribution without the need to introduce an ad hoc constant.