Large eddy simulation of a stratified boundary layer under an oscillatory current

2009 ◽  
Vol 643 ◽  
pp. 233-266 ◽  
Author(s):  
BISHAKHDATTA GAYEN ◽  
SUTANU SARKAR ◽  
JOHN R. TAYLOR

A numerical study based on large eddy simulation is performed to investigate a bottom boundary layer under an oscillating tidal current. The focus is on the boundary layer response to an external stratification. The thermal field shows a mixed layer that is separated from the external stratified fluid by a thermocline. The mixed layer grows slowly in time with an oscillatory modulation by the tidal flow. Stratification strongly affects the mean velocity profiles, boundary layer thickness and turbulence levels in the outer region although the effect on the near-bottom unstratified fluid is relatively mild. The turbulence is asymmetric between the accelerating and decelerating stages. The asymmetry is more pronounced with increasing stratification. There is an overshoot of the mean velocity in the outer layer; this jet is linked to the phase asymmetry of the Reynolds shear stress gradient by using the simulation data to examine the mean momentum equation. Depending on the height above the bottom, there is a lag of the maximum turbulent kinetic energy, dissipation and production with respect to the peak external velocity and the value of the lag is found to be influenced by the stratification. Flow instabilities and turbulence in the bottom boundary layer excite internal gravity waves that propagate away into the ambient. Unlike the steady case, the phase lines of the internal waves change direction during the tidal cycle and also from near to far field. The frequency spectrum of the propagating wave field is analysed and found to span a narrow band of frequencies clustered around 45°.

2013 ◽  
Vol 13 (1) ◽  
pp. 1489-1526 ◽  
Author(s):  
D. Jarecka ◽  
H. Pawlowska ◽  
W. W. Grabowski ◽  
A. A. Wyszogrodzki

Abstract. This paper discusses aircraft observations and large-eddy simulation (LES) of the 15 May 2008, North Sea boundary-layer clouds from the EUCAARI-IMPACT field campaign. These clouds were advected from the north-east by the prevailing lower-tropspheric winds, and featured stratocumulus-over-cumulus cloud formations. Almost-solid stratocumulus deck in the upper part of the relatively deep weakly decoupled marine boundary layer overlaid a field of small cumuli with a cloud fraction of ~10%. The two cloud formations featured distinct microphysical characteristics that were in general agreement with numerous past observations of strongly-diluted shallow cumuli on the one hand and solid marine boundary-layer stratocumulus on the other. Macrophysical and microphysical cloud properties were reproduced well by the double-moment warm-rain microphysics large-eddy simulation. A novel feature of the model is its capability to locally predict homogeneity of the subgrid-scale mixing between the cloud and its cloud-free environment. In the double-moment warm-rain microphysics scheme, the homogeneity is controlled by a single parameter α, that ranges from 0 to 1 and limiting values representing the homogeneous and the extremely inhomogeneous mixing scenarios, respectively. Parameter α depends on the characteristic time scales of the droplet evaporation and of the turbulent homogenization. In the model, these scales are derived locally based on the subgrid-scale turbulent kinetic energy, spatial scale of cloudy filaments, the mean cloud droplet radius, and the humidity of the cloud-free air entrained into the cloud. Simulated mixing is on average quite inhomogeneous, with the mean parameter α around 0.7 across the entire depth of the cloud field, but with local variations across almost the entire range, especially near the base and the top of the cloud field.


1996 ◽  
Vol 118 (2) ◽  
pp. 248-254 ◽  
Author(s):  
B. J. Boersma ◽  
F. T. M. Nieuwstadt

In this paper, we use Large-Eddy Simulation (LES) to compute a fully-developed turbulent flow in a curved pipe. The results allow us to study how the curvature influences the mean velocity profile and also various turbulent statistics. We find reasonable agreement with the few experiments that are available. Our simulation also allows a detailed study of secondary motion in the cross section of the pipe which are caused by the centrifugal acceleration due to the pipe curvature. It is known that this secondary motion may consist of one, two, or three circulation cells. In our simulation results we find one circulation cell.


2009 ◽  
Vol 66 (6) ◽  
pp. 1465-1487 ◽  
Author(s):  
Roel A. J. Neggers ◽  
Martin Köhler ◽  
Anton C. M. Beljaars

Abstract This study considers the question of what is the least complex bulk mass flux framework that can still conceptually reproduce the smoothly varying coupling between the shallow convective cloud layer and the subcloud mixed layer. To this end, the model complexity of the classic single bulk mass flux scheme is enhanced. Inspired by recent large-eddy simulation results, the authors argue that two relatively minor but key conceptual modifications are already sufficient to achieve this goal: (i) retaining a dry transporting updraft in the moist limit and (ii) applying continuous updraft area partitioning to this dual mass flux (DualM) framework. The dry updraft represents all internal mixed layer updrafts that terminate near the mixed layer top, whereas the moist updraft represents all updrafts that condense and rise out of the mixed layer as buoyant cumulus clouds. The continuous area partitioning between the dry and moist updraft is a function of moist convective inhibition above the mixed layer top. Updraft initialization is a function of the updraft area fraction and is therefore consistent with the updraft definition. It is argued that the model complexity thus enhanced is sufficient to allow reproduction of various phenomena involved in the cloud–subcloud coupling, namely (i) dry countergradient transport within the mixed layer that is independent of the moist updraft, (ii) soft triggering of moist convective flux throughout the boundary layer, and (iii) a smooth response to smoothly varying forcings, including the reproduction of gradual transitions to and from shallow cumulus convection. The DualM framework is evaluated by implementing in the Eddy Diffusivity Mass Flux (EDMF) boundary layer scheme of the ECMWF’s Integrated Forecasting System. Single column model experiments are evaluated against large-eddy simulation results for a range of different cases that span a broad parameter space of cloud–subcloud coupling intensities. The results illustrate that also in numerical practice the DualM framework can reproduce gradual transitions to and from shallow cumulus convection. Model behavior is further explored through experiments in which model complexity is purposely reduced, thus mimicking a single bulk updraft setup. This gives more insight into the new model-internal interactions and explains the obtained case results.


2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Bobin Wang ◽  
Guixiang Cui

Thermal instability induced by solar radiation is the most common condition of urban atmosphere in daytime. Compared to researches under neutral conditions, only a few numerical works studied the unstable urban boundary layer and the effect of buoyancy force is unclear. In this paper, unstably stratified turbulent boundary layer flow over three-dimensional urban-like building arrays with ground heating is simulated. Large eddy simulation is applied to capture main turbulence structures and the effect of buoyancy force on turbulence can be investigated. Lagrangian dynamic subgrid scale model is used for complex flow together with a wall function, taking into account the large pressure gradient near buildings. The numerical model and method are verified with the results measured in wind tunnel experiment. The simulated results satisfy well with the experiment in mean velocity and temperature, as well as turbulent intensities. Mean flow structure inside canopy layer varies with thermal instability, while no large secondary vortex is observed. Turbulent intensities are enhanced, as buoyancy force contributes to the production of turbulent kinetic energy.


Sign in / Sign up

Export Citation Format

Share Document