Horizontal turbulent diffusion in a convective mixed layer

2014 ◽  
Vol 758 ◽  
pp. 553-564 ◽  
Author(s):  
Junshi Ito ◽  
Hiroshi Niino ◽  
Mikio Nakanishi

AbstractA large eddy simulation (LES) is used to estimate a reliable horizontal turbulent diffusion coefficient, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}K_{{h}}$, in a convective mixed layer (CML). The introduction of a passive scalar field with a fixed horizontal gradient at a given time enables $K_{{h}}$ estimation as a function of height, based on the simulated turbulent horizontal scalar flux. Here $K_{{h}}$ is found to be of the order of $100\ {\mathrm{m}}^2\ {\mathrm{s}}^{-1}$ for a typical terrestrial atmospheric CML. It is shown to scale by the product of the CML convective velocity, $w_{*}$, and its depth, $h$. Here $K_{{h}}$ is characterized by a vertical profile in the CML: it is large near both the bottom and top of the CML, where horizontal flows associated with convection are large. The equation pertaining to the temporal rate of change of a horizontal scalar flux suggests that $K_{{h}}$ is determined by a balance between production and pressure correlation at a fully developed stage. Pressure correlation near the bottom of the CML is localized in convergence zones near the boundaries of convective cells and becomes large within an eddy turnover time, $h/w_{*}$, after the introduction of the passive scalar field.

1991 ◽  
Vol 231 ◽  
pp. 665-688 ◽  
Author(s):  
James B. Young ◽  
Thomas J. Hanratty

An extension of an axial viewing optical technique (first used by Lee, Adrian & Hanratty) is described which allows the determination of the turbulence characteristics of solid particles being transported by water in a pipe. Measurements are presented of the mean radial velocity, the mean rate of change radial velocity, the mean-square of the radial and circumferential fluctuations, the Eulerian turbulent diffusion coefficient, and the Lagrangian turbulent diffusion coefficient. A particular focus is to explore the influence of slip velocity for particles which have small time constants. It is found that with increasing slip velocity the magnitude of the turbulent velocity fluctuations remains unchanged but that the turbulent diffusivity decreases. The measurements of the average rate of change of particle velocity are consistent with the notion that particles move from regions of high fluid turbulence to regions of low fluid turbulence. Measurements of the root-mean-square of the fluctuations of the rate of change of particle velocity allow an estimation of the average magnitude of the particle slip in a highly turbulent flow, which needs to be known to analyse the motion of particles not experiencing a Stokes drag.


2014 ◽  
Vol 742 ◽  
pp. 701-719
Author(s):  
Daan D. J. A. van Sommeren ◽  
C. P. Caulfield ◽  
Andrew W. Woods

AbstractWe perform experiments to study the mixing of passive scalar by a buoyancy-induced turbulent flow in a long narrow vertical tank. The turbulent flow is associated with the downward mixing of a small flux of dense aqueous saline solution into a relatively large upward flux of fresh water. In steady state, the mixing region is of finite extent, and the intensity of the buoyancy-driven mixing is described by a spatially varying turbulent diffusion coefficient $\kappa _v(z)$ which decreases linearly with distance $z$ from the top of the tank. We release a pulse of passive scalar into either the fresh water at the base of the tank, or the saline solution at the top of the tank, and we measure the subsequent mixing of the passive scalar by the flow using image analysis. In both cases, the mixing of the passive scalar (the dye) is well-described by an advection–diffusion equation, using the same turbulent diffusion coefficient $\kappa _v(z)$ associated with the buoyancy-driven mixing of the dynamic scalar. Using this advection–diffusion equation with spatially varying turbulent diffusion coefficient $\kappa _v(z)$, we calculate the residence time distribution (RTD) of a unit mass of passive scalar released as a pulse at the bottom of the tank. The variance in this RTD is equivalent to that produced by a uniform eddy diffusion coefficient with value $\kappa _e= 0.88 \langle \kappa _v \rangle $, where $\langle \kappa _v \rangle $ is the vertically averaged eddy diffusivity. The structure of the RTD is also qualitatively different from that produced by a flow with uniform eddy diffusion coefficient. The RTD using $\kappa _v$ has a larger peak value and smaller values at early times, associated with the reduced diffusivity at the bottom of the tank, and manifested mathematically by a skewness $\gamma _1\approx 1.60$ and an excess kurtosis $\gamma _2\approx 4.19 $ compared to the skewness and excess kurtosis of $\gamma _1\approx 1.46$, $\gamma _2 \approx 3.50$ of the RTD produced by a constant eddy diffusion coefficient with the same variance.


1994 ◽  
Vol 29 (2-3) ◽  
pp. 221-232
Author(s):  
M.J. McCormick

Abstract Four one-dimensional models which have been used to characterize surface mixed layer (ML) processes and the thermal structure are described. Although most any model can be calibrated to mimic surface water temperatures, it does not imply that the corresponding mixing processes are well described. Eddy diffusion or "K" models can exhibit this problem. If a ML model is to be useful for water quality applications, then it must be able to resolve storm events and, therefore, be able to simulate the ML depth, h, and its time rate of change, dh/dt. A general water quality model is derived from mass conservation principles to demonstrate how ML models can be used in a physically meaningful way to address water quality issues.


2016 ◽  
Vol 799 ◽  
pp. 159-199 ◽  
Author(s):  
A. Briard ◽  
T. Gomez ◽  
C. Cambon

The present work aims at developing a spectral model for a passive scalar field and its associated scalar flux in homogeneous anisotropic turbulence. This is achieved using the paradigm of eddy-damped quasi-normal Markovian (EDQNM) closure extended to anisotropic flows. In order to assess the validity of this approach, the model is compared to several detailed direct numerical simulations (DNS) and experiments of shear-driven flows and isotropic turbulence with a mean scalar gradient at moderate Reynolds numbers. This anisotropic modelling is then used to investigate the passive scalar dynamics at very high Reynolds numbers. In the framework of homogeneous isotropic turbulence submitted to a mean scalar gradient, decay and growth exponents for the cospectrum and scalar energies are obtained analytically and assessed numerically thanks to EDQNM closure. With the additional presence of a mean shear, the scaling of the scalar flux and passive scalar spectra in the inertial range are investigated and confirm recent theoretical predictions. Finally, it is found that, in shear-driven flows, the small scales of the scalar second-order moments progressively return to isotropy when the Reynolds number increases.


2021 ◽  
Author(s):  
Shengyu Shi ◽  
Jiale Chen ◽  
Clarisse Bourdelle ◽  
Xiang Jian ◽  
Tomas Odstrcil ◽  
...  

Abstract The behavior of heavy/high-Z impurity tungsten (W) in an improved high-performance fully non-inductive discharge on EAST with ITER-like divertor (ILD) is analyzed. It is found that W could be well controlled. The causes of no W accumulation are clarified by analyzing the background plasma parameters and modeling the W transport. It turns out that the electron temperature (T_e) and its gradient are usually high while the toroidal rotation and density peaking of the bulk plasma are small. In this condition, the modeled W turbulent diffusion coefficient is big enough to offset the total turbulent and neoclassical pinch, so that W density profile for zero particle flux will not be very peaked. Combining NEO and TGLF for the W transport coefficient and the impurity transport code STRAHL, not only the core W density profile is predicted but also the radiated information mainly produced by W in the experiment can be closely reconstructed. At last, the physics of controlling W accumulation by electron cyclotron resonance heating (ECRH) is illustrated considering the effects of changed T_e by ECRH on ionization balance and transport of W. It shows that the change of ionization and recombination balance by changed T_e is not enough to explain the experimental observation of W behavior, which should be attributed to the changed W transport. By comparing the W transport coefficients in two kinds of plasmas with different T_e profiles, it is shown that high T_e and its gradient play a key role to generate large turbulent diffusion through increasing the growth rate of linear instability so that W accumulation is prevented.


2013 ◽  
Vol 721 ◽  
pp. 541-577 ◽  
Author(s):  
Amin Rasam ◽  
Geert Brethouwer ◽  
Arne V. Johansson

AbstractIn Marstorpet al. (J. Fluid Mech., vol. 639, 2009, pp. 403–432), an explicit algebraic subgrid stress model (EASSM) for large-eddy simulation (LES) was proposed, which was shown to considerably improve LES predictions of rotating and non-rotating turbulent channel flow. In this paper, we extend that work and present a new explicit algebraic subgrid scalar flux model (EASSFM) for LES, based on the modelled transport equation of the subgrid-scale (SGS) scalar flux. The new model is derived using the same kind of methodology that leads to the explicit algebraic scalar flux model of Wikströmet al. (Phys. Fluids, vol. 12, 2000, pp. 688–702). The algebraic form is based on a weak equilibrium assumption and leads to a model that depends on the resolved strain-rate and rotation-rate tensors, the resolved scalar-gradient vector and, importantly, the SGS stress tensor. An accurate prediction of the SGS scalar flux is consequently strongly dependent on an accurate description of the SGS stresses. The new EASSFM is therefore primarily used in connection with the EASSM, since this model can accurately predict SGS stresses. The resulting SGS scalar flux is not necessarily aligned with the resolved scalar gradient, and the inherent dependence on the resolved rotation-rate tensor makes the model suitable for LES of rotating flow applications. The new EASSFM (together with the EASSM) is validated for the case of passive scalar transport in a fully developed turbulent channel flow with and without system rotation. LES results with the new model show good agreement with direct numerical simulation data for both cases. The new model predictions are also compared to those of the dynamic eddy diffusivity model (DEDM) and improvements are observed in the prediction of the resolved and SGS scalar quantities. In the non-rotating case, the model performance is studied at all relevant resolutions, showing that its predictions of the Nusselt number are much less dependent on the grid resolution and are more accurate. In channel flow with wall-normal rotation, where all the SGS stresses and fluxes are non-zero, the new model shows significant improvements over the DEDM predictions of the resolved and SGS quantities.


1999 ◽  
Vol 90 (1) ◽  
pp. 47-82 ◽  
Author(s):  
Margaret a. Lemone ◽  
Mingyu Zhou ◽  
Chin-Hoh Moeng ◽  
Donald H. Lenschow ◽  
L. Jay Miller ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document