scholarly journals Stratification Effects in the Turbulent Boundary Layer beneath a Melting Ice Shelf: Insights from Resolved Large-Eddy Simulations

2019 ◽  
Vol 49 (7) ◽  
pp. 1905-1925 ◽  
Author(s):  
Catherine A. Vreugdenhil ◽  
John R. Taylor
Keyword(s):  
Boundary Layer ◽  
Large Scale ◽  
Rectangular Domain ◽  
Large Eddy Simulations ◽  
Intermittent Flow ◽  
Transfer Coefficients ◽  
Ice Shelf ◽  
Obukhov Length ◽  
Mixed Regime ◽  
Large Eddy

AbstractOcean turbulence contributes to the basal melting and dissolution of ice shelves by transporting heat and salt toward the ice. The meltwater causes a stable salinity stratification to form beneath the ice that suppresses turbulence. Here we use large-eddy simulations motivated by the ice shelf–ocean boundary layer (ISOBL) to examine the inherently linked processes of turbulence and stratification, and their influence on the melt rate. Our rectangular domain is bounded from above by the ice base where a dynamic melt condition is imposed. By varying the speed of the flow and the ambient temperature, we identify a fully turbulent, well-mixed regime and an intermittently turbulent, strongly stratified regime. The transition between regimes can be characterized by comparing the Obukhov length, which provides a measure of the distance away from the ice base where stratification begins to dominate the flow, to the viscous length scale of the interfacial sublayer. Upper limits on simulated turbulent transfer coefficients are used to predict the transition from fully to intermittently turbulent flow. The predicted melt rate is sensitive to the choice of the heat and salt transfer coefficients and the drag coefficient. For example, when coefficients characteristic of fully developed turbulence are applied to intermittent flow, the parameterized three-equation model overestimates the basal melt rate by almost a factor of 10. These insights may help to guide when existing parameterizations of ice melt are appropriate for use in regional or large-scale ocean models, and may also have implications for other ice–ocean interactions such as fast ice or drifting ice.

scholarly journals Ice-shelf ocean boundary layer dynamics from large-eddy simulations

2021 ◽  
Author(s):  
Carolyn Branecky Begeman ◽  
Xylar Asay-Davis ◽  
Luke Van Roekel
Keyword(s):  
Boundary Layer ◽  
Ocean Circulation ◽  
Large Eddy Simulations ◽  
Global Ocean ◽  
Small Scale ◽  
Parameter Study ◽  
Ice Shelf ◽  
Ice Shelves ◽  
Large Eddy ◽  
Ocean Boundary

Abstract. Small scale, turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small scale processes, which regulate heat transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

scholarly journals An intercomparison of Large-Eddy Simulations of the Martian daytime convective boundary layer

2016 ◽  
Author(s):  
Tanguy Bertrand ◽  
Aymeric Spiga ◽  
Scot Rafkin ◽  
Arnaud Colaitis ◽  
François Forget ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Radiative Forcing ◽  
Large Scale ◽  
Numerical Models ◽  
Horizontal Resolution ◽  
Large Eddy Simulations ◽  
Horizontal Wind ◽  
Convective Boundary ◽  
Airborne Dust ◽  
Large Eddy

Abstract. Large-Eddy Simulations (LES) for Mars resolve the Planetary Boundary Layer (PBL) turbulent dynamics by using a very fine horizontal resolution of a few tens of meters. LES modeling is becoming a more and more useful tool to prepare the robotic exploration of Mars by providing means to evaluate the intensity of convective plumes and vortices, horizontal wind gustiness, and turbulent fluctuations of temperature in the Martian PBL. In such context, and given the relative paucity of turbulence-related measurements on Mars, an intercomparison of LES models is a fruitful way to evaluate the models' predictions and to indicate possible areas of improvement. Thus, to prepare the landing of the ExoMars Schiaparelli lander (also named ExoMars Demonstrator Module, EDM), scheduled for October 2016, the results of the Laboratoire de Météorologie Dynamique (LMD) and South-West Research Institute (SwRI) LES models have been compared. The objective of this study is to determine the range of uncertainties, and dispersions, of the two numerical models' predictions, for the critical phase of the spacecraft's descent in the Martian daytime turbulent PBL. First, a strategy is defined to ensure similar radiative forcing in both the LMD and SwRI models. Then, LES are performed over a flat terrain with and without large-scale ambient horizontal wind. The LMD and SwRI Martian LES models predict similar temporal evolution of the PBL and organization in the horizontal and vertical wind fields. However, the convective motions in the daytime PBL are more vigorous by a factor 1.5–2 in SwRI results than in LMD results, independently of the presence or not of ambient horizontal wind. This discrepancy is further investigated through sensitivity studies to surface conditions, ambient wind, and airborne dust loading.

scholarly journals Large-Eddy Simulations of Stratified Atmospheric Boundary Layers: Comparison of Different Subgrid Models

2020 ◽  
Author(s):  
Srinidhi N. Gadde ◽  
Anja Stieren ◽  
Richard J. A. M. Stevens
Keyword(s):  
Boundary Layer ◽  
Boundary Layers ◽  
Large Scale ◽  
Wind Farm ◽  
Research Area ◽  
Large Eddy Simulations ◽  
Turbulence Energy ◽  
Atmospheric Boundary Layers ◽  
Smagorinsky Model ◽  
Large Eddy

Abstract The development and assessment of subgrid-scale (SGS) models for large-eddy simulations of the atmospheric boundary layer is an active research area. In this study, we compare the performance of the classical Smagorinsky model, the Lagrangian-averaged scale-dependent (LASD) model, and the anisotropic minimum dissipation (AMD) model. The LASD model has been widely used in the literature for 15 years, while the AMD model was recently developed. Both the AMD and the LASD models allow three-dimensional variation of SGS coefficients and are therefore suitable to model heterogeneous flows over complex terrain or around a wind farm. We perform a one-to-one comparison of these SGS models for neutral, stable, and unstable atmospheric boundary layers. We find that the LASD and the AMD models capture the logarithmic velocity profile and the turbulence energy spectra better than the Smagorinsky model. In stable and unstable boundary-layer simulations, the AMD and LASD model results agree equally well with results from a high-resolution reference simulation. The performance analysis of the models reveals that the computational overhead of the AMD model and the LASD model compared to the Smagorinsky model is approximately 10% and 30% respectively. The LASD model has a higher computational and memory overhead because of the global filtering operations and Lagrangian tracking procedure, which can result in bottlenecks when the model is used in extensive simulations. These bottlenecks are absent in the AMD model, which makes it an attractive SGS model for large-scale simulations of turbulent boundary layers.

scholarly journals Comments on “Langmuir Turbulence and Surface Heating in the Ocean Surface Boundary Layer”

2018 ◽  
Vol 48 (2) ◽  
pp. 455-458 ◽  
Author(s):  
Yign Noh ◽  
Yeonju Choi
Keyword(s):  
Boundary Layer ◽  
Ocean Surface ◽  
Large Eddy Simulations ◽  
Surface Heating ◽  
Length Scale ◽  
Surface Boundary ◽  
Langmuir Turbulence ◽  
Surface Boundary Layer ◽  
Obukhov Length ◽  
Large Eddy

AbstractUsing large-eddy simulations (LES) it is shown that the depth of a diurnal thermocline h should be scaled by the Zilitinkevich scale LZ, not by the Monin–Obukhov length scale LMO, contrary to the proposition by Pearson et al. Their argument to explain the slower increase of h than LMO using the effect of the preexisting thermocline is also invalid.

Toward a new ice-shelf melt rate parameterization with large-eddy simulations

2020 ◽  
Author(s):  
Carolyn Branecky Begeman ◽  
Xylar Asay-Davis ◽  
Luke Van Roekel
Keyword(s):  
Large Scale ◽  
Ocean Model ◽  
Large Eddy Simulations ◽  
Stratified Turbulence ◽  
Ice Shelf ◽  
Ice Shelves ◽  
Ocean Coupling ◽  
Large Eddy ◽  
Dynamic Melting ◽  
Ocean Boundary

<p>Predictions of ice shelf melting depend on dynamical insights into ocean boundary layers below ice shelves. Fundamental questions regarding the nature of stratified turbulence below the sloped and ablating ice shelf base remain. Laboratory experiments, direct numerical simulations, and observations have yielded important insights, but have yet to produce a robust relationship between ice shelf melt rates and shear- and buoyancy-driven mixing. This relationship is the target of our Large-Eddy Simulations (LES) of the ice-shelf ocean boundary layer. Several new developments were applied to the LES code PALM to produce dynamic melting as well as tides. In this presentation, we demonstrate these new model capabilities. We contrast profiles of vertical turbulent fluxes of heat, salt and momentum across different simulated ice shelf settings: cold, shear-dominated settings vs. warm, buoyancy-dominated settings. We also discuss our recent work toward a new ice-shelf melt parameterization for use in large-scale ocean models on the basis of these simulations. A new melt parameterization is a critical component of ongoing ice-ocean coupling efforts, both to place melt rate predictions on a more physical footing and to achieve convergence with vertical ocean model resolution, on which current parameterizations fail.</p>

Turbulent dynamics and ice-shelf basal melt rates from large-eddy simulations

2021 ◽  
Author(s):  
Carolyn Branecky Begeman ◽  
Xylar Asay-Davis ◽  
Luke Van Roekel
Keyword(s):  
Boundary Layer ◽  
Ocean Circulation ◽  
Large Eddy Simulations ◽  
Far Field ◽  
Ice Shelf ◽  
Thermal Exchange ◽  
Boundary Layer Turbulence ◽  
Shelf Slope ◽  
Linear Relationships ◽  
Large Eddy

<p>Large-eddy simulations are used to investigate boundary layer turbulence and its control on ice-shelf basal melt rates in Antarctic settings. We present simulations at relatively low thermal driving and low ice-shelf basal slopes, resulting in simulated melt rates from 10s cm/yr to several m/yr. Our results are broadly consistent with the linear relationships between far-field thermal driving and melt rate and between ice-shelf slope and melt rate reported by previous studies. The simulated thermal exchange coefficient is lower than recommended values; thermal exchange becomes less efficient as stratification increases.  In our simulations, shear production of turbulent kinetic energy outweighs buoyant production, as found below Larsen C Ice Shelf through recent microstructure measurements. We also find that turbulent intensity and melt rate vary significantly with the orientation between the ice-shelf slope and the far-field flow, even at low ice-shelf slopes. Our results suggest that numerical ocean models employing the standard ice-shelf melt parameterization will underestimate slope effects on ice-shelf melt rates even if they capture the mean buoyancy effects on boundary layer flow. The proposed slope effects would modify feedbacks between ocean circulation and ice-shelf geometry and tidal variability in ice-shelf melt rates.</p>

scholarly journals Using the Sensitivity of Large-Eddy Simulations to Evaluate Atmospheric Boundary Layer Models

2012 ◽  
Vol 69 (5) ◽  
pp. 1582-1601 ◽  
Author(s):  
Gilles Bellon ◽  
Bjorn Stevens
Keyword(s):  
Boundary Layer ◽  
Atmospheric Boundary Layer ◽  
Stationary State ◽  
Large Scale ◽  
Large Eddy Simulations ◽  
Tropospheric Temperature ◽  
Unstable Regime ◽  
Shallow Cumulus ◽  
Large Eddy ◽  
Bulk Model

Abstract A simple framework to study the sensitivity of atmospheric boundary layer (ABL) models to the large-scale conditions and forcings is introduced. This framework minimizes the number of parameters necessary to describe the large-scale conditions, subsidence, and radiation. Using this framework, the sensitivity of the stationary ABL to the large-scale boundary conditions [underlying sea surface temperature (SST) and overlying humidity and temperature in the free troposphere] is investigated in large-eddy simulations (LESs). For increasing SST or decreasing free-tropospheric temperature, the LES exhibits a transition from a cloud-free, well-mixed ABL stationary state, through a cloudy, well-mixed stationary state and a stable shallow cumulus stationary state, to an unstable regime with a deepening shallow cumulus layer. For a warm SST, when increasing free-tropospheric humidity, the LES exhibits a transition from a stable shallow cumulus stationary state, through a stable cumulus-under-stratus stationary state, to an unstable regime with a deepening, cumulus-under-stratus layer. For a cool SST, when increasing the free-tropospheric humidity, the LES stationary state exhibits a transition from a cloud-free, well-mixed ABL regime, through a well-mixed cumulus-capped regime, to a stratus-capped regime with a decoupling between the subcloud and cloud layers. This dataset can be used to evaluate other ABL models. As an example, the sensitivity of a bulk model based on the mixing-line model is presented. This bulk model reproduces the LES sensitivity to SST and free-tropospheric temperature for the stable and unstable shallow cumulus regimes, but it is less successful at reproducing the LES sensitivity to free-tropospheric humidity for both shallow cumulus and well-mixed regimes.

scholarly journals Coherent structure of the convective boundary layer derived from large-eddy simulations

1989 ◽  
Vol 200 ◽  
pp. 511-562 ◽  
Author(s):  
Helmut Schmidt ◽  
Ulrich Schumann
Keyword(s):  
Boundary Layer ◽  
Boundary Condition ◽  
Large Scale ◽  
Three Dimensional ◽  
Large Eddy Simulations ◽  
Stable Layer ◽  
Convective Boundary ◽  
Velocity Fluctuations ◽  
Large Eddy ◽  
The Mean

Turbulence in the convective boundary layer (CBL) uniformly heated from below and topped by a layer of uniformly stratified fluid is investigated for zero mean horizontal flow using large-eddy simulations (LES). The Rayleigh number is effectively infinite, the Froude number of the stable layer is 0.09 and the surface roughness height relative to the height of the convective layer is varied between 10−6 and 10−2. The LES uses a finite-difference method to integrate the three-dimensional grid-volume-averaged Navier–Stokes equations for a Boussinesq fluid. Subgrid-scale (SGS) fluxes are determined from algebraically approximated second-order closure (SOC) transport equations for which all essential coefficients are determined from the inertial-range theory. The surface boundary condition uses the Monin–Obukhov relationships. A radiation boundary condition at the top of the computational domain prevents spurious reflections of gravity waves. The simulation uses 160 × 160 × 48 grid cells. In the asymptotic state, the results in terms of vertical mean profiles of turbulence statistics generally agree very well with results available from laboratory and atmospheric field experiments. We found less agreement with respect to horizontal velocity fluctuations, pressure fluctuations and dissipation rates, which previous investigations tend to overestimate. Horizontal spectra exhibit an inertial subrange. The entrainment heat flux at the top of the CBL is carried by cold updraughts and warm downdraughts in the form of wisps at scales comparable with the height of the boundary layer. Plots of instantaneous flow fields show a spoke pattern in the lower quarter of the CBL which feeds large-scale updraughts penetrating into the stable layer aloft. The spoke pattern has also been found in a few previous investigations. Small-scale plumes near the surface and remote from strong updraughts do not merge together but decay while rising through large-scale downdraughts. The structure of updraughts and downdraughts is identified by three-dimensional correlation functions and conditionally averaged fields. The mean circulation extends vertically over the whole boundary layer. We find that updraughts are composed of quasi-steady large-scale plumes together with transient rising thermals which grow in size by lateral entrainment. The skewness of the vertical velocity fluctuations is generally positive but becomes negative in the lowest mesh cells when the dissipation rate exceeds the production rate due to buoyancy near the surface, as is the case for very rough surfaces. The LES results are used to determine the root-mean-square value of the surface friction velocity and the mean temperature difference between the surface and the mixed layer as a function of the roughness height. The results corroborate a simple model of the heat transfer in the surface layer.

Sign in / Sign up

Export Citation Format

Share Document