Large-eddy simulation (LES) is used to model turbulent winds in a nominally neutral atmospheric boundary layer at varying mesh resolutions. The boundary layer is driven by wind shear with zero surface heat flux and is capped by a stable inversion. Because of entrainment the boundary layer is in a weakly stably stratified regime. The simulations use meshes varying from 1282× 64 to 10242× 512 grid points in a fixed computational domain of size (2560, 2560, 896) m. The subgrid-scale (SGS) parameterizations used in the LES vary with the mesh spacing. Low-order statistics, spectra, and structure functions are compared on the different meshes and are used to assess grid convergence in the simulations. As expected, grid convergence is primarily achieved in the middle of the boundary layer where there is scale separation between the energy-containing and dissipative eddies. Near the surface second-order statistics do not converge on the meshes studied. The analysis also highlights differences between one-dimensional and two-dimensional velocity spectra; differences are attributed to sampling errors associated with aligning the horizontal coordinates with the vertically veering mean wind direction. Higher-order structure functions reveal non-Gaussian statistics on all scales, but are highly dependent on the mesh resolution. A generalized logarithmic law and a k−1spectral scaling regime are identified with mesh-dependent parameters in agreement with previously published results.
- Award ID(s):
- 2009418
- NSF-PAR ID:
- 10411233
- Date Published:
- Journal Name:
- Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Volume:
- 380
- Issue:
- 2219
- ISSN:
- 1364-503X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
null (Ed.)We study bubble motion in a vertical capillary tube under an external flow. Bretherton ( J. Fluid Mech. , vol. 10, issue 2, 1961, pp. 166–188) has shown that, without external flow, a bubble can spontaneously rise when the Bond number ( ${Bo} \equiv \rho g R^2 / \gamma$ ) is above the critical value ${Bo}_{cr}=0.842$ , where $\rho$ is the liquid density, $g$ the gravitational acceleration, $R$ the tube radius and $\gamma$ the surface tension. It was then shown by Magnini et al. ( Phys. Rev. Fluids , vol. 4, issue 2, 2019, 023601) that the presence of an imposed liquid flow, in the same (upward) direction as buoyancy, accelerates the bubble and thickens the liquid film around it. In this work we carry out a systematic study of the bubble motion under a wide range of upward and downward external flows, focusing on the inertialess regime with Bond numbers above the critical value. We show that a rich variety of bubble dynamics occurs when an external downward flow is applied, opposing the buoyancy-driven rise of the bubble. We reveal the existence of a critical capillary number of the external downward flow ( ${Ca}_l \equiv \mu U_l/\gamma$ , where $\mu$ is the fluid viscosity and $U_l$ is the mean liquid speed) at which the bubble arrests and changes its translational direction. Depending on the relative direction of gravity and the external flow, the thickness of the film separating the bubble surface and the tube inner wall follows two distinct solution branches. The results from theory, experiments and numerical simulations confirm the existence of the two solution branches and reveal that the two branches overlap over a finite range of ${Ca}_l$ , thus suggesting non-unique, history-dependent solutions for the steady-state film thickness under the same external flow conditions. Furthermore, inertialess symmetry-breaking shape profiles at steady state are found as the bubble transits near the tipping points of the solution branches, which are shown in both experiments and three-dimensional numerical simulations.more » « less
-
more » « less
Nearly, all dense suspensions undergo dramatic and abrupt thickening transitions in their flow behavior when sheared at high stresses. Such transitions occur when the dominant interactions between the suspended particles shift from hydrodynamic to frictional. Here, we interpret abrupt shear thickening as a precursor to a rigidity transition and give a complete theory of the viscosity in terms of a universal crossover scaling function from the frictionless jamming point to a rigidity transition associated with friction, anisotropy, and shear. Strikingly, we find experimentally that for two different systems—cornstarch in glycerol and silica spheres in glycerol—the viscosity can be collapsed onto a single universal curve over a wide range of stresses and volume fractions. The collapse reveals two separate scaling regimes due to a crossover between frictionless isotropic jamming and frictional shear jamming, with different critical exponents. The material-specific behavior due to the microscale particle interactions is incorporated into a scaling variable governing the proximity to shear jamming, that depends on both stress and volume fraction. This reformulation opens the door to importing the vast theoretical machinery developed to understand equilibrium critical phenomena to elucidate fundamental physical aspects of the shear thickening transition.
-
more » « less
Many of the cell membrane's vital functions are achieved by the self-organization of the proteins and biopolymers embedded in it. The protein dynamics is in part determined by its drag. A large number of these proteins can polymerize to form filaments.
In vitro studies of protein–membrane interactions often involve using rigid beads coated with lipid bilayers, as a model for the cell membrane. Motivated by this, we use slender-body theory to compute the translational and rotational resistance of a single filamentous protein embedded in the outer layer of a supported bilayer membrane and surrounded on the exterior by a Newtonian fluid. We first consider the regime where the two layers are strongly coupled through their inter-leaflet friction. We find that the drag along the parallel direction grows linearly with the filament's length and quadratically with the length for the perpendicular and rotational drag coefficients. These findings are explained using scaling arguments and by analysing the velocity fields around the moving filament. We then present and discuss the qualitative differences between the drag of a filament moving in a freely suspended bilayer and a supported membrane as a function of the membrane's inter-leaflet friction. Finally, we briefly discuss how these findings can be used in experiments to determine membrane rheology. In summary, we present a formulation that allows computation of the effects of membrane properties (its curvature, viscosity and inter-leaflet friction), and the exterior and interior three-dimensional fluids’ depth and viscosity on the drag of a rod-like/filamentous protein, all in a unified theoretical framework. -
more » « less
Abstract We provide a first-principles analysis of the energy fluxes in the oceanic internal wave field. The resulting formula is remarkably similar to the renowned phenomenological formula for the turbulent dissipation rate in the ocean, which is known as the finescale parameterization. The prediction is based on the wave turbulence theory of internal gravity waves and on a new methodology devised for the computation of the associated energy fluxes. In the standard spectral representation of the wave energy density, in the two-dimensional vertical wavenumber–frequency (
m –ω ) domain, the energy fluxes associated with the steady state are found to be directed downscale in both coordinates, closely matching the finescale parameterization formula in functional form and in magnitude. These energy transfers are composed of a “local” and a “scale-separated” contributions; while the former is quantified numerically, the latter is dominated by the induced diffusion process and is amenable to analytical treatment. Contrary to previous results indicating an inverse energy cascade from high frequency to low, at odds with observations, our analysis of all nonzero coefficients of the diffusion tensor predicts a direct energy cascade. Moreover, by the same analysis fundamental spectra that had been deemed “no-flux” solutions are reinstated to the status of “constant-downscale-flux” solutions. This is consequential for an understanding of energy fluxes, sources, and sinks that fits in the observational paradigm of the finescale parameterization, solving at once two long-standing paradoxes that had earned the name of “oceanic ultraviolet catastrophe.”Significance Statement The global circulation models cannot resolve the scales of the oceanic internal waves. The finescale parameterization of turbulent dissipation, a formula grounded in observations, is the standard tool by which the energy transfers due to internal waves are incorporated in the global models. Here, we provide an interpretation of this parameterization formula building on the first-principles statistical theory describing energy transfers between waves at different scales. Our result is in agreement with the finescale parameterization and points out a large contribution to the energy fluxes due to a type of wave interactions (local) usually disregarded. Moreover, the theory on which the traditional understanding of the parameterization is mainly built, a “diffusion approximation,” is known to be partly in contradiction with observations. We put forward a solution to this problem, visualized by means of “streamlines” that improve the intuition of the direction of the energy cascade.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1098/rsta.2021.0094
