An official website of the United States government
Abstract We study higher uniformity properties of the Möbius function$$\mu $$, the von Mangoldt function$$\Lambda $$, and the divisor functions$$d_k$$on short intervals$$(X,X+H]$$with$$X^{\theta +\varepsilon } \leq H \leq X^{1-\varepsilon }$$for a fixed constant$$0 \leq \theta < 1$$and any$$\varepsilon>0$$. More precisely, letting$$\Lambda ^\sharp $$and$$d_k^\sharp $$be suitable approximants of$$\Lambda $$and$$d_k$$and$$\mu ^\sharp = 0$$, we show for instance that, for any nilsequence$$F(g(n)\Gamma )$$, we have$$\begin{align*}\sum_{X < n \leq X+H} (f(n)-f^\sharp(n)) F(g(n) \Gamma) \ll H \log^{-A} X \end{align*}$$ when$$\theta = 5/8$$and$$f \in \{\Lambda , \mu , d_k\}$$or$$\theta = 1/3$$and$$f = d_2$$. As a consequence, we show that the short interval Gowers norms$$\|f-f^\sharp \|_{U^s(X,X+H]}$$are also asymptotically small for any fixedsfor these choices of$$f,\theta $$. As applications, we prove an asymptotic formula for the number of solutions to linear equations in primes in short intervals and show that multiple ergodic averages along primes in short intervals converge in$$L^2$$. Our innovations include the use of multiparameter nilsequence equidistribution theorems to control type$$II$$sums and an elementary decomposition of the neighborhood of a hyperbola into arithmetic progressions to control type$$I_2$$sums.
more »
« less
Lagrangian model for passive scalar gradients in turbulence
The equation for the fluid velocity gradient along a Lagrangian trajectory immediately follows from the Navier–Stokes equation. However, such an equation involves two terms that cannot be determined from the velocity gradient along the chosen Lagrangian path: the pressure Hessian and the viscous Laplacian. A recent model handles these unclosed terms using a multi-level version of the recent deformation of Gaussian fields (RDGF) closure (Johnson & Meneveau,Phys. Rev. Fluids, vol. 2 (7), 2017, 072601). This model is in remarkable agreement with direct numerical simulations (DNS) data and works for arbitrary Taylor Reynolds numbers$$\textit {Re}_\lambda$$. Inspired by this, we develop a Lagrangian model for passive scalar gradients in isotropic turbulence. The equation for passive scalar gradients also involves an unclosed term in the Lagrangian frame, namely the scalar gradient diffusion term, which we model using the RDGF approach. However, comparisons of the statistics obtained from this model with DNS data reveal substantial errors due to erroneously large fluctuations generated by the model. We address this defect by incorporating into the closure approximation information regarding the scalar gradient production along the local trajectory history of the particle. This modified model makes predictions for the scalar gradients, their production rates, and alignments with the strain-rate eigenvectors that are in very good agreement with DNS data. However, while the model yields valid predictions up to$$\textit {Re}_\lambda \approx 500$$, beyond this, the model breaks down.
more »
« less
- Award ID(s):
- 2042346
- PAR ID:
- 10512557
- Publisher / Repository:
- Journal of Fluid Mechanics
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 964
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Except in the trivial case of spatially uniform flow, the advection–diffusion operator of a passive scalar tracer is linear and non-self-adjoint. In this study, we exploit the linearity of the governing equation and present an analytical eigenfunction approach for computing solutions to the advection–diffusion equation in two dimensions given arbitrary initial conditions, and when the advecting flow field at any given time is a plane parallel shear flow. Our analysis illuminates the specific role that the non-self-adjointness of the linear operator plays in the solution behaviour, and highlights the multiscale nature of the scalar mixing problem given the explicit dependence of the eigenvalue–eigenfunction pairs on a multiscale parameter$$q=2{\rm i}k\,{\textit {Pe}}$$, where$$k$$is the non-dimensional wavenumber of the tracer in the streamwise direction, and$${\textit {Pe}}$$is the Péclet number. We complement our theoretical discussion on the spectra of the operator by computing solutions and analysing the effect of shear flow width on the scale-dependent scalar decay of tracer variance, and characterize the distinct self-similar dispersive processes that arise from the shear flow dispersion of an arbitrarily compact tracer concentration. Finally, we discuss limitations of the present approach and future directions.more » « less
-
Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ($$Pr$$) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behaviour and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing$$Pr$$. Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as$$Pr$$is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when$$Pr$$deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.more » « less
-
Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius$$R$$and axial length$$10{\rm \pi} R$$at friction Reynolds numbers up to$$Re_\tau =5200$$using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor$$\lambda$$overshoots by$$2\,\%$$and undershoots by$$0.6\,\%$$the Prandtl friction law at low and high$$Re$$ranges, respectively. In addition,$$\lambda$$in our results is slightly higher than in Pirozzoliet al.(J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichiet al.(Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to$$y^+=250$$, has not yet developed a plateau farther away from the wall in the pipes even for the$$Re_\tau =5200$$cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with$$Re_\tau$$– are lower in the pipe than in the channel, but the difference decreases with increasing$$Re_\tau$$. While the wall value is slightly lower in the channel than in the pipe at the same$$Re_\tau$$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a$$k^{-1}$$dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.more » « less
-
The motion of a disk in a Langmuir film bounding a liquid substrate is a classical hydrodynamic problem, dating back to Saffman (J. Fluid Mech., vol. 73, 1976, p. 593) who focused upon the singular problem of translation at large Boussinesq number,$${\textit {Bq}}\gg 1$$. A semianalytic solution of the dual integral equations governing the flow at arbitrary$${\textit {Bq}}$$was devised by Hugheset al.(J. Fluid Mech., vol. 110, 1981, p. 349). When degenerated to the inviscid-film limit$${\textit {Bq}}\to 0$$, it produces the value$$8$$for the dimensionless translational drag, which is$$50\,\%$$larger than the classical$$16/3$$-value corresponding to a free surface. While that enhancement has been attributed to surface incompressibility, the mathematical reasoning underlying the anomaly has never been fully elucidated. Here we address the inviscid limit$${\textit {Bq}}\to 0$$from the outset, revealing a singular mechanism where half of the drag is contributed by the surface pressure. We proceed beyond that limit, considering a nearly inviscid film. A naïve attempt to calculate the drag correction using the reciprocal theorem fails due to an edge singularity of the leading-order flow. We identify the formation of a boundary layer about the edge of the disk, where the flow is primarily in the azimuthal direction with surface and substrate stresses being asymptotically comparable. Utilising the reciprocal theorem in a fluid domain tailored to the asymptotic topology of the problem produces the drag correction$$(8\,{\textit {Bq}}/{\rm \pi} ) [ \ln (2/{\textit {Bq}}) + \gamma _E+1]$$,$$\gamma _E$$being the Euler–Mascheroni constant.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.375
