Comparison between the sea-breeze circulation day and normal day Reynolds stress anisotropy in the lower atmospheric region

Author(s):  
Sayahnya Roy ◽  
Alexei Sentchev ◽  
François G. Schmitt ◽  
Patrick Augustin ◽  
Marc Fourmentin

<p>This study shows the comparison between the sea-breeze circulation (SBC) day and normal day turbulent characteristics and the Reynolds stress anisotropy in the lower atmospheric region. The Reynolds stress tensor is responsible for the dissipation and transport of mean kinetic energy. The variability of the turbulent kinetic energy due to the Reynolds stress anisotropy modulates the air quality. A 20 Hz Ultrasonic anemometer was deployed in the coastal area of northern France to measure the temporal wind variability for the duration of one year five months.</p><p>The SBC was detected by a change in wind direction from the West to the East during the day time. We found that the axial component of the turbulent kinetic energy is higher than the other two through an axisymmetric expansion, and energy ellipsoid has a cigar shape due to SBC. During this time the dominance of small scale zonal turbulent motions was observed. Also, the probability of a higher degree of wind anisotropy due to SBC generates large mean kinetic energy within the lower troposphere. Moreover, the production of larger negative turbulent kinetic energy due to SBC was evident.</p>

2012 ◽  
Vol 7 (1) ◽  
pp. 53-69
Author(s):  
Vladimir Dulin ◽  
Yuriy Kozorezov ◽  
Dmitriy Markovich

The present paper reports PIV (Particle Image Velocimetry) measurements of turbulent velocity fluctuations statistics in development region of an axisymmetric free jet (Re = 28 000). To minimize measurement uncertainty, adaptive calibration, image processing and data post-processing algorithms were utilized. On the basis of theoretical analysis and direct measurements, the paper discusses effect of PIV spatial resolution on measured statistical characteristics of turbulent fluctuations. Underestimation of the second-order moments of velocity derivatives and of the turbulent kinetic energy dissipation rate due to a finite size of PIV interrogation area and finite thickness of laser sheet was analyzed from model spectra of turbulent velocity fluctuations. The results are in a good agreement with the measured experimental data. The paper also describes performance of possible ways to account for unresolved small-scale velocity fluctuations in PIV measurements of the dissipation rate. In particular, a turbulent viscosity model can be efficiently used to account for the unresolved pulsations in a free turbulent flow


2009 ◽  
Vol 9 (7) ◽  
pp. 2335-2353 ◽  
Author(s):  
W. W. Grabowski ◽  
L.-P. Wang

Abstract. A large set of rising adiabatic parcel simulations is executed to investigate the combined diffusional and accretional growth of cloud droplets in maritime and continental conditions, and to assess the impact of enhanced droplet collisions due to small-scale cloud turbulence. The microphysical model applies the droplet number density function to represent spectral evolution of cloud and rain/drizzle drops, and various numbers of bins in the numerical implementation, ranging from 40 to 320. Simulations are performed applying two traditional gravitational collection kernels and two kernels representing collisions of cloud droplets in the turbulent environment, with turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3. The overall result is that the rain initiation time significantly depends on the number of bins used, with earlier initiation of rain when the number of bins is low. This is explained as a combination of the increase of the width of activated droplet spectrum and enhanced numerical spreading of the spectrum during diffusional and collisional growth when the number of model bins is low. Simulations applying around 300 bins seem to produce rain at times which no longer depend on the number of bins, but the activation spectra are unrealistically narrow. These results call for an improved representation of droplet activation in numerical models of the type used in this study. Despite the numerical effects that impact the rain initiation time in different simulations, the turbulent speedup factor, the ratio of the rain initiation time for the turbulent collection kernel and the corresponding time for the gravitational kernel, is approximately independent of aerosol characteristics, parcel vertical velocity, and the number of bins used in the numerical model. The turbulent speedup factor is in the range 0.75–0.85 and 0.60–0.75 for the turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3, respectively.


2020 ◽  
Author(s):  
João Bettencourt ◽  
Carlos Guedes Soares

<p>The Azores Current-Front system coincides with the northern limit of the subtropical gyre in  the Eastern North Atlantic. The mean zonal jet is positioned south of the Azores archipelago  and extends from west of the mid-atlantic ridge to the Gulf of Cadiz, where it partially  turns south. North of the main jet, a sub-surface counter-current is found, flowing westwards. The associated thermal front separates the warm subtropical waters from the colder subpolar waters. The instantaneous flow in the Azores Current/Front system is characterized by the presence of meandering currents with length scales of 200 km that regularly shed anticyclonic warm water and cyclonic cold water eddies to the north and south of the mean jet axis, respectively, due to vortex stretching and the planetary beta effect. The time scale of eddy shedding is 100-200 days. On the meandering arms of the current, downwelling <br>and upwelling cells are found and sharp thermal gradients are formed and a residual poleward heat transport is observed. The instability cycle that originates the mesoscale meanders and the eddies is well-known from quasi-geostrophic and primitive equation models initialized from a basic baroclinic state: a first phase of baroclinic instability feeds on available potential energy to raise eddy kinetic energy levels, that, in a second phase feed the mean kinetic energy by Reynolds stress convergence. The cycle repeats itself as long as the APE reservoir is filled at the end of each cycle.</p><p>However, seasonal variability of the zonal jet dynamics has not been addressed before and it can provide valuable insights in to the variations of the Eastern North Atlantic between the subtropical and subpolar gyres. We use a primitive equation regional ocean model of the Eastern Central North Atlantic with realistic climatological wind and thermal forcing to study the yearly cycle of meandering, eddy shedding and restoration of the mean jet in the Azores/Current system. We observe an semi-annual cycle in the jet's kinetic energy with maxima in Summer/Winter and minima in early Spring/Autumn. Potential energy conversion by baroclinic instability occurs throughout the year but is predominant in the first half of the year. The mean kinetic energy draws from the turbulent kinetic energy through Reynolds stress convergence in periods of 50 - 100 days, that are followed by short barotropic instability periods. During Winter, Reynolds stress convergence, and thus mean jet reinforcement from the mesoscale eddy field, occurs along the jet meridional extent, in the top 500 m of the water column, but from Spring to Autumn it is observed only in the southern flank of the mean jet axis.</p>


2014 ◽  
Vol 760 ◽  
pp. 304-312 ◽  
Author(s):  
Farid Karimpour ◽  
Subhas K. Venayagamoorthy

AbstractIn this study, we revisit the consequence of assuming equilibrium between the rates of production ($P$) and dissipation $({\it\epsilon})$ of the turbulent kinetic energy $(k)$ in the highly anisotropic and inhomogeneous near-wall region. Analytical and dimensional arguments are made to determine the relevant scales inherent in the turbulent viscosity (${\it\nu}_{t}$) formulation of the standard $k{-}{\it\epsilon}$ model, which is one of the most widely used turbulence closure schemes. This turbulent viscosity formulation is developed by assuming equilibrium and use of the turbulent kinetic energy $(k)$ to infer the relevant velocity scale. We show that such turbulent viscosity formulations are not suitable for modelling near-wall turbulence. Furthermore, we use the turbulent viscosity $({\it\nu}_{t})$ formulation suggested by Durbin (Theor. Comput. Fluid Dyn., vol. 3, 1991, pp. 1–13) to highlight the appropriate scales that correctly capture the characteristic scales and behaviour of $P/{\it\epsilon}$ in the near-wall region. We also show that the anisotropic Reynolds stress ($\overline{u^{\prime }v^{\prime }}$) is correlated with the wall-normal, isotropic Reynolds stress ($\overline{v^{\prime 2}}$) as $-\overline{u^{\prime }v^{\prime }}=c_{{\it\mu}}^{\prime }(ST_{L})(\overline{v^{\prime 2}})$, where $S$ is the mean shear rate, $T_{L}=k/{\it\epsilon}$ is the turbulence (decay) time scale and $c_{{\it\mu}}^{\prime }$ is a universal constant. ‘A priori’ tests are performed to assess the validity of the propositions using the direct numerical simulation (DNS) data of unstratified channel flow of Hoyas & Jiménez (Phys. Fluids, vol. 18, 2006, 011702). The comparisons with the data are excellent and confirm our findings.


1997 ◽  
Vol 334 ◽  
pp. 61-86 ◽  
Author(s):  
PAUL PICCIRILLO ◽  
CHARLES W. VAN ATTA

Experiments were carried out in a new type of stratified flow facility to study the evolution of turbulence in a mean flow possessing both uniform stable stratification and uniform mean shear.The new facility is a thermally stratified wind tunnel consisting of ten independent supply layers, each with its own blower and heaters, and is capable of producing arbitrary temperature and velocity profiles in the test section. In the experiments, four different sized turbulence-generating grids were used to study the effect of different initial conditions. All three components of the velocity were measured, along with the temperature. Root-mean-square quantities and correlations were measured, along with their corresponding power and cross-spectra.As the gradient Richardson number Ri = N2/(dU/dz)2 was increased, the downstream spatial evolution of the turbulent kinetic energy changed from increasing, to stationary, to decreasing. The stationary value of the Richardson number, Ricr, was found to be an increasing function of the dimensionless shear parameter Sq2/∈ (where S = dU/dz is the mean velocity shear, q2 is the turbulent kinetic energy, and ∈ is the viscous dissipation).The turbulence was found to be highly anisotropic, both at the small scales and at the large scales, and anisotropy was found to increase with increasing Ri. The evolution of the velocity power spectra for Ri [les ] Ricr, in which the energy of the large scales increases while the energy in the small scales decreases, suggests that the small-scale anisotropy is caused, or at least amplified, by buoyancy forces which reduce the amount of spectral energy transfer from large to small scales. For the largest values of Ri, countergradient buoyancy flux occurred for the small scales of the turbulence, an effect noted earlier in the numerical results of Holt et al. (1992), Gerz et al. (1989), and Gerz & Schumann (1991).


2016 ◽  
Vol 33 (4) ◽  
pp. 817-837 ◽  
Author(s):  
Justine M. McMillan ◽  
Alex E. Hay ◽  
Rolf G. Lueck ◽  
Fabian Wolk

AbstractThe ability to estimate the rate of dissipation (ε) of turbulent kinetic energy at middepth in a high-speed tidal channel using broadband acoustic Doppler current profilers (ADCPs) is assessed by making comparisons to direct measurements of ε obtained using shear probes mounted on a streamlined underwater buoy. The investigation was carried out in Grand Passage, Nova Scotia, Canada, where the depth-averaged flow speed reached 2 m s−1 and the Reynolds number was 8 × 107. The speed bin–averaged dissipation rates estimated from the ADCP data agree with the shear probe data to within a factor of 2. Both the ADCP and the shear probe measurements indicate a linear dependence of ε on the cube of the flow speed during flood and much lower dissipation rates during ebb. The ebb–flood asymmetry and the small-scale intermittency in ε are also apparent in the lognormal distributions of the shear probe data. Possible sources of bias and error in the ε estimates are investigated, and the most likely causes of the discrepancy between the ADCP and shear probe estimates are the cross-channel separation of the instruments and the high degree of spatial variability that exists in the channel.


2010 ◽  
Vol 37 (4) ◽  
pp. 648-656 ◽  
Author(s):  
Ahmad Sana ◽  
Hitoshi Tanaka

A total of seven versions of two-equation turbulence models (four versions of low Reynolds number k–ε model, one k–ω model and two versions of k–ε / k–ω blended models) are tested against the direct numerical simulation (DNS) data of a one-dimensional oscillatory boundary layer with flat crested free-stream velocity that results from a steep pressure gradient. A detailed comparison has been made for cross-stream velocity, turbulent kinetic energy (TKE), Reynolds stress, and ratio of Reynolds stress and turbulent kinetic energy. It is observed that the newer versions of k–ε model perform very well in predicting the velocity, turbulent kinetic energy, and Reynolds stress. The k–ω model and blended models underestimate the peak value of turbulent kinetic energy that may be explained by the Reynolds stress to TKE ratio in the logarithmic zone. The maximum bottom shear stress is well predicted by the k–ε model proposed by Sana et al. and the original k–ω model.


2015 ◽  
Vol 770 ◽  
pp. 210-246 ◽  
Author(s):  
M. Mehrabadi ◽  
S. Tenneti ◽  
R. Garg ◽  
S. Subramaniam

Gas-phase velocity fluctuations due to mean slip velocity between the gas and solid phases are quantified using particle-resolved direct numerical simulation. These fluctuations are termed pseudo-turbulent because they arise from the interaction of particles with the mean slip even in ‘laminar’ gas–solid flows. The contribution of turbulent and pseudo-turbulent fluctuations to the level of gas-phase velocity fluctuations is quantified in initially ‘laminar’ and turbulent flow past fixed random particle assemblies of monodisperse spheres. The pseudo-turbulent kinetic energy $k^{(f)}$ in steady flow is then characterized as a function of solid volume fraction ${\it\phi}$ and the Reynolds number based on the mean slip velocity $\mathit{Re}_{m}$. Anisotropy in the Reynolds stress is quantified by decomposing it into isotropic and deviatoric parts, and its dependence on ${\it\phi}$ and $Re_{m}$ is explained. An algebraic stress model is proposed that captures the dependence of the Reynolds stress on ${\it\phi}$ and $Re_{m}$. Gas-phase velocity fluctuations in freely evolving suspensions undergoing elastic and inelastic particle collisions are also quantified. The flow corresponds to homogeneous gas–solid systems, with high solid-to-gas density ratio and particle diameter greater than dissipative length scales. It is found that for the parameter values considered here, the level of pseudo-turbulence differs by only 15 % from the values for equivalent fixed beds. The principle of conservation of interphase turbulent kinetic energy transfer is validated by quantifying the interphase transfer terms in the evolution equations of kinetic energy for the gas-phase and solid-phase fluctuating velocity. It is found that the collisional dissipation is negligible compared with the viscous dissipation for the cases considered in this study where the freely evolving suspensions attain a steady state starting from an initial condition where the particles are at rest.


2008 ◽  
Vol 8 (4) ◽  
pp. 14717-14763 ◽  
Author(s):  
W. W. Grabowski ◽  
L.-P. Wang

Abstract. A large set of rising adiabatic parcel simulations is executed to investigate the combined diffusional and accretional growth of cloud droplets in maritime and continental conditions, and to assess the impact of enhanced droplet collisions due to small-scale cloud turbulence. The microphysical model applies the droplet number density function to represent spectral evolution of cloud and rain/drizzle drops, and various numbers of bins in the numerical implementation, ranging from 40 to 320. Simulations are performed applying two traditional gravitational collection kernels and two kernels representing collisions of cloud droplets in the turbulent environment, with turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3. The overall result is that the rain initiation time significantly depends on the number of bins used, with earlier initiation of rain when the number of bins is low. This is explained as a combination of the increase of the width of activated droplet spectrum and enhanced numerical spreading of the spectrum during diffusional and collisional growth when the number of model bins is low. Simulations applying around 300 bins seem to produce rain at times which no longer depend on the number of bins, but the activation spectra are unrealistically narrow. These results call for an improved representation of droplet activation in numerical models of the type used in this study. Despite the numerical effects that impact the rain initiation time in different simulations, the turbulent speedup factor, the ratio of the rain initiation time for the turbulent collection kernel and the corresponding time for the gravitational kernel, is approximately independent of aerosol characteristics, parcel vertical velocity, and the number of bins used in the numerical model. The turbulent speedup factor is in the range 0.75–0.85 and 0.60–0.75 for the turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3, respectively.


Sign in / Sign up

Export Citation Format

Share Document