Laboratory experiments were performed on a geometrically scaled vertical-axis wind turbine model over an unprecedented range of Reynolds numbers, including and exceeding those of the full-scale turbine. The study was performed in the high-pressure environment of the Princeton High Reynolds number Test Facility (HRTF). Utilizing highly compressed air as the working fluid enabled extremely high Reynolds numbers while still maintaining dynamic similarity by matching the tip speed ratio (defined as the ratio of tip velocity to free stream, $\unicode[STIX]{x1D706}=\unicode[STIX]{x1D714}R/U$ ) and Mach number (defined at the turbine tip, $Ma=\unicode[STIX]{x1D714}R/a$ ). Preliminary comparisons are made with measurements from the full-scale field turbine. Peak power for both the field data and experiments resides around $\unicode[STIX]{x1D706}=1$ . In addition, a systematic investigation of trends with Reynolds number was performed in the laboratory, which revealed details about the asymptotic behaviour. It was shown that the parameter that characterizes invariance in the power coefficient was the Reynolds number based on blade chord conditions ( $Re_{c}$ ). The power coefficient reaches its asymptotic value when $Re_{c}>1.5\times 10^{6}$ , which is higher than what the field turbine experiences. The asymptotic power curve is found, which is invariant to further increases in Reynolds number.
more »
« less
Stochastic theory and direct numerical simulations of the relative motion of high-inertia particle pairs in isotropic turbulence
The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$ is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.
more »
« less
- Award ID(s):
- 1435953
- NSF-PAR ID:
- 10025836
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 813
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 205 to 249
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ( $0.1\leqslant \unicode[STIX]{x1D6F7}\leqslant 0.4$ ) and particle Reynolds numbers ( $0.1\leqslant Re_{p}\leqslant 10$ ). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal $\unicode[STIX]{x1D6F7}$ . This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. ( J. Rheol. , vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the $Re_{p}$ . The wide spectrum of the local shear rate and its dependency on $\unicode[STIX]{x1D6F7}$ and $Re_{p}$ point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, $O(Re_{p})\sim 10$ .more » « less
-
The dynamics of air bubbles in turbulent Rayleigh–Bénard (RB) convection is described for the first time using laboratory experiments and complementary numerical simulations. We performed experiments at $Ra=5.5\times 10^{9}$ and $1.1\times 10^{10}$ , where streams of 1 mm bubbles were released at various locations from the bottom of the tank along the path of the roll structure. Using three-dimensional particle tracking velocimetry, we simultaneously tracked a large number of bubbles to inspect the pair dispersion, $R^{2}(t)$ , for a range of initial separations, $r$ , spanning one order of magnitude, namely $25\unicode[STIX]{x1D702}\leqslant r\leqslant 225\unicode[STIX]{x1D702}$ ; here $\unicode[STIX]{x1D702}$ is the local Kolmogorov length scale. Pair dispersion, $R^{2}(t)$ , of the bubbles within a quiescent medium was also determined to assess the effect of inhomogeneity and anisotropy induced by the RB convection. Results show that $R^{2}(t)$ underwent a transition phase similar to the ballistic-to-diffusive ( $t^{2}$ -to- $t^{1}$ ) regime in the vicinity of the cell centre; it approached a bulk behavior $t^{3/2}$ in the diffusive regime as the distance away from the cell centre increased. At small $r$ , $R^{2}(t)\propto t^{1}$ is shown in the diffusive regime with a lower magnitude compared to the quiescent case, indicating that the convective turbulence reduced the amplitude of the bubble’s fluctuations. This phenomenon associated to the bubble path instability was further explored by the autocorrelation of the bubble’s horizontal velocity. At large initial separations, $R^{2}(t)\propto t^{2}$ was observed, showing the effect of the roll structure.more » « less
-
more » « less
Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.
-
We study the variation of $\unicode[STIX]{x1D707}$ -invariants in Hida families with residually reducible Galois representations. We prove a lower bound for these invariants which is often expressible in terms of the $p$ -adic zeta function. This lower bound forces these $\unicode[STIX]{x1D707}$ -invariants to be unbounded along the family, and we conjecture that this lower bound is an equality. When $U_{p}-1$ generates the cuspidal Eisenstein ideal, we establish this conjecture and further prove that the $p$ -adic $L$ -function is simply a power of $p$ up to a unit (i.e. $\unicode[STIX]{x1D706}=0$ ). On the algebraic side, we prove analogous statements for the associated Selmer groups which, in particular, establishes the main conjecture for such forms.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2016.859
