In ice‐covered lakes, penetrative radiation warms fluid beneath a diffusive boundary layer, thereby increasing its density and providing energy for convection in a diurnally active, deepening mixed layer. Shallow regions are differentially heated to warmer temperatures, driving turbulent gravity currents that transport warm water downslope and into the basin interior. We examine the energetics of these processes, focusing on the rate at which penetrative radiation supplies energy that is available to drive fluid motion. Using numerical simulations that resolve convective plumes, gravity currents, and the secondary instabilities leading to entrainment, we show that advective fluxes due to differential heating contribute to the evolution of the mixed layer in waterbodies with significant shallow areas. A heat balance is used to assess the relative importance of differential heating to the one‐dimensional effects of radiative heating and diffusive cooling at the ice‐water interface in lakes of varying morphologies.
- Award ID(s):
- 1723400
- NSF-PAR ID:
- 10109398
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 873
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 942 to 976
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.more » « less
-
Abstract The presence of supercooled water in polar regions causes anchor ice to grow on submerged objects, generating costly problems for engineered materials and life-endangering risks for benthic communities. The factors driving underwater ice accretion are poorly understood, and passive prevention mechanisms remain unknown. Here we report that the Antarctic scallop Adamussium colbecki appears to remain ice-free in shallow Antarctic marine environments where underwater ice growth is prevalent. In contrast, scallops colonized by bush sponges in the same microhabitat grow ice and are removed from the population. Characterization of the Antarctic scallop shells revealed a hierarchical micro-ridge structure with sub-micron nano-ridges which promotes directed icing. This concentrates the formation of ice on the growth rings while leaving the regions in between free of ice, and appears to reduce ice-to-shell adhesion when compared to temperate species that do not possess highly ordered surface structures. The ability to control the formation of ice may enable passive underwater anti-icing protection, with the removal of ice possibly facilitated by ocean currents or scallop movements. We term this behavior cryofouling avoidance. We posit that the evolution of natural anti-icing structures is a key trait for the survival of Antarctic scallops in anchor ice zones.more » « less
-
Abstract Evaporation from undulating soil surfaces is rarely studied due to limited modeling theory and inadequate experimental data linking dynamic soil and atmospheric interactions. The goal of this paper is to provide exploratory insights into evaporation behavior from undulating soil surfaces under turbulent conditions through numerical and experimental approaches. A previously developed and verified coupled free flow and porous media flow model was extended by incorporating turbulent airflow through Reynolds‐averaged Navier–Stokes equations. The model explicitly describes the relevant physical processes and the key properties in the free flow, porous media, and at the interface, allowing for the analysis of coupled exchange fluxes. An experiment aiming to simultaneously collect the data of the boundary layer and soil evaporation in the area around the soil–atmosphere interface was conducted using a wind tunnel integrated with a soil tank. The turbulent boundary layer above the undulating soil surface was captured using high‐resolution hot‐wire anemometry, confirming the presence of recirculation zones in the undulating valleys and locally low evaporative flux. Experimental data were used to validate the extended model, and modeling results demonstrate that turbulent airflow enhances evaporation and shortens the duration of Stage I. The surface geometry significantly affects the local evaporative flux by influencing the vapor distribution, concentration gradient, and water availability at the soil surface, especially when recirculation zones form in the valleys. As a joint result of turbulence and surface undulations, the influence of wind speed on both the local and system‐level evaporation rate is restricted.
-
An incoming canonical spatially developing turbulent boundary layer (SDTBL) over a 2-D curved hill is numerically investigated via the Reynolds-averaged Navier–Stokes (RANS) equations plus two eddy-viscosity models: the K−ω SST (henceforth SST) and the Spalart–Allmaras (henceforth SA) turbulence models. A spatially evolving thermal boundary layer has also been included, assuming temperature as a passive scalar (Pr = 0.71) and a turbulent Prandtl number, Prt, of 0.90 for wall-normal turbulent heat flux modeling. The complex flow with a combined strong adverse/favorable streamline curvature-driven pressure gradient caused by concave/convex surface curvatures has been replicated from wind-tunnel experiments from the literature, and the measured velocity and pressure fields have been used for validation purposes (the thermal field was not experimentally measured). Furthermore, direct numerical simulation (DNS) databases from the literature were also employed for the incoming turbulent flow assessment. Concerning first-order statistics, the SA model demonstrated a better agreement with experiments where the turbulent boundary layer remained attached, for instance, in Cp, Cf, and Us predictions. Conversely, the SST model has shown a slightly better match with experiments over the flow separation zone (in terms of Cp and Cf) and in Us profiles just upstream of the bubble. The Reynolds analogy, based on the St/(Cf/2) ratio, holds in zero-pressure gradient (ZPG) zones; however, it is significantly deteriorated by the presence of streamline curvature-driven pressure gradient, particularly due to concave wall curvature or adverse-pressure gradient (APG). In terms of second-order statistics, the SST model has better captured the positively correlated characteristics of u′ and v′ or positive Reynolds shear stresses ( > 0) inside the recirculating zone. Very strong APG induced outer secondary peaks in and turbulence production as well as an evident negative slope on the constant shear layer.more » « less