BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Charles Meneveau

Charles Meneveau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynnlin.CDcond-mat.softphysics.comp-ph

Catalog footprint

What is connected

26works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Resolved and subgrid-scale crossing trajectory effects in Eulerian large eddy simulations of polydisperse droplet transport

We study the dispersion characteristics of slightly buoyant droplets in a turbulent jet using large eddy simulations (LES). The droplet number density fields are represented using an Eulerian approach with the dispersed phase modelled using the Fast-Eulerian method (Ferry & Balachandar 2001) that includes the droplet rise velocity. Radial concentration profiles and turbulent concentration fluxes for droplets of different sizes are analyzed to quantify the "trajectory crossing effect", when relative motions between particles and turbulent eddies tend to reduce turbulent diffusion. For finer LES grid resolutions, the model captures the differential, size based dispersion characteristics of the droplets with the transverse dispersion of the larger droplet sizes suppressed, since trajectory crossing effects are explicitly resolved in LES. We examine a similarity solution model for the size dependent radial concentration profiles based on a modified Schmidt number derived from the theory of turbulent diffusion of particles in the atmosphere proposed by Csanady 1963. The results are validated with the high resolution LES data and show good agreement. Then, the size dependent Schmidt number model is reformulated as a model for unresolved subgrid-scale trajectory crossing effects and used to calculate the subgrid concentration flux in a coarse LES of a turbulent jet, with slightly buoyant droplets injected at the centerline in the self-similar region of the jet. The results are compared to a simulation with higher grid resolution and a coarse simulation with a constant Schmidt number SGS model. We find that the subgrid model enhances the prediction accuracy of the concentration profiles and turbulent concentration flux for the coarse LES.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Lagrangian relaxation towards equilibrium wall model for large eddy simulation

A large eddy simulation wall model is developed based on a formal interpretation of quasi-equilibrium that governs the momentum balance integrated in the wall-normal direction. The model substitutes the law-of-the-wall velocity profile for a smooth surface into the wall-normal integrated momentum balance, leading to a Lagrangian relaxation towards equilibrium (LaRTE) transport equation for the friction velocity vector $\mathbf{u}_τ(x,z,t)$. This partial differential equation includes a relaxation timescale governing the rate at which the wall stress can respond to imposed fluctuations due to the inertia of the fluid layer from the wall to the wall-model height. A-priori tests based on channel flow direct numerical simulation (DNS) data show that the identified relaxation timescale ensures self-consistency with assumed quasi-equilibrium conditions. The new approach enables us to formally distinguish quasi-equilibrium from additional, non-equilibrium contributions to the wall stress. A particular model for non-equilibrium contributions is derived, motivated by laminar Stokes layer dynamics in the viscous sublayer when applying fast varying pressure gradients. The new wall modeling approach is first tested in standard equilibrium channel flow in order to document various properties of the approach. The model is then applied in LES of channel flow with a suddenly applied spanwise pressure gradient (SSPG). The resulting mean wall stress evolution is compared to DNS with good agreement. At the onset of the SSPG, the laminar Stokes layer develops rapidly while the LaRTE portion of the stress has a delayed response due to its inherent relaxation dynamics. Results also highlight open challenges such as modeling the response of near-wall turbulence occurring above the viscous sublayer and at timescales faster than quasi-equilibrium conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Data compression for turbulence databases using spatio-temporal sub-sampling and local re-simulation

Motivated by specific data and accuracy requirements for building numerical databases of turbulent flows, data compression using spatio-temporal sub-sampling and local re-simulation is proposed. Numerical re-simulation experiments for decaying isotropic turbulence based on sub-sampled data are undertaken. The results and error analyses are used to establish parameter choices for sufficiently accurate sub-sampling and sub-domain re-simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Filtered actuator disks: Theory and application to wind turbine models in large eddy simulation

The actuator disk model (ADM) continues to be a popular wind turbine representation in large eddy simulations (LES) of large wind farms. Computational restrictions typically limit the number of grid points across the rotor of each actuator disk and require spatial filtering to smoothly distribute the applied force distribution on discrete grid points. At typical grid resolutions, simulations cannot capture all of the vorticity shed behind the disk and subsequently over-predict power by upwards of 10%. To correct these modeling errors, we propose a vortex cylinder model to quantify the shed vorticity when a filtered force distribution is applied at the actuator disk. This model is then used to derive a correction factor for numerical simulations that collapses the power curve for simulations at various filter widths and grid resolutions onto the curve obtained using axial momentum theory. The proposed correction therefore facilitates accurate power measurements in LES without resorting to highly refined numerical grids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A closure for Lagrangian velocity gradient evolution in turbulence using recent deformation mapping of initially Gaussian fields

The statistics of the velocity gradient tensor in turbulent flows are of both theoretical and practical importance. The Lagrangian view provides a privileged perspective for studying the dynamics of turbulence in general, and of the velocity gradient tensor in particular. Stochastic models for the Lagrangian evolution of velocity gradients in isotropic turbulence, with closure models for the pressure Hesssian and viscous Laplacian, have been shown to reproduce important features such as non-Gaussian probability distributions, skewness and vorticity strain-rate alignments. The Recent Fluid Deformation (RFD) closure introduced the idea of mapping an isotropic Lagrangian pressure Hessian as upstream initial condition using the fluid deformation tensor. Recent work on a Gaussian fields closure, however, has shown that even Gaussian isotropic velocity fields contain significant anisotropy for the conditional pressure Hessian tensor due to the inherent velocity-pressure couplings, and that assuming an isotropic pressure Hessian as upstream condition may not be realistic. In this paper, Gaussian isotropic field statistics are used to generate more physical upstream conditions for the recent fluid deformation mapping. In this new framework, known isotropy relations can be satisfied {\it a priori} and no DNS-tuned coefficients are necessary. A detailed comparison of results from the new model, referred to as the recent deformation of Gaussian fields (RDGF) closure, with existing models and DNS shows the improvements gained, especially in various single-time statistics of the velocity gradient tensor at moderate Reynolds numbers. Application to arbitrarily high Reynolds numbers remains an open challenge for this type of model, however.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extended self-similarity in moment-generating-functions in wall-bounded turbulence at high Reynolds number

In wall-bounded turbulence, the moment generating functions (MGFs) of the streamwise velocity fluctuations $\left<\exp(qu_z^+)\right>$ develop power-law scaling as a function of the wall normal distance $z/δ$. Here $u$ is the streamwise velocity fluctuation, $+$ indicates normalization in wall units (averaged friction velocity), $z$ is the distance from the wall, $q$ is an independent variable and $δ$ is the boundary layer thickness. Previous work has shown that this power-law scaling exists in the log-region {\small $3Re_τ^{0.5}\lesssim z^+$, $z\lesssim 0.15δ$}, where $Re_τ$ is the friction velocity-based Reynolds numbers. Here we present empirical evidence that this self-similar scaling can be extended, including bulk and viscosity-affected regions $30<z^+$, $z<δ$, provided the data are interpreted with the Extended-Self-Similarity (ESS), i.e. self-scaling of the MGFs as a function of one reference value, $q_o$. ESS also improves the scaling properties, leading to more precise measurements of the scaling exponents. The analysis is based on hot-wire measurements from boundary layers at $Re_τ$ ranging from $2700$ to $13000$ from the Melbourne High-Reynolds-Number-Turbulent-Boundary-Layer-Wind-Tunnel. Furthermore, we investigate the scalings of the filtered, large-scale velocity fluctuations $u^L_z$ and of the remaining small-scale component, $u^S_z=u_z-u^L_z$. The scaling of $u^L_z$ falls within the conventionally defined log region and depends on a scale that is proportional to {\small $l^+\sim Re_τ^{1/2}$}; the scaling of $u^{S}_z$ extends over a much wider range from $z^+\approx 30$ to $z\approx 0.5δ$. Last, we present a theoretical construction of two multiplicative processes for $u^L_z$ and $u^S_z$ that reproduce the empirical findings concerning the scalings properties as functions of $z^+$ and in the ESS sense.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimal smoothing length scale for actuator line models of wind turbine blades

The actuator line model (ALM) is a commonly used method to represent lifting surfaces such as wind turbine blades within Large-Eddy Simulations (LES). In the ALM the lift and drag forces are replaced by an imposed body force which is typically smoothed over several grid points using a Gaussian kernel with some prescribed smoothing width $ε$. To date, the choice of $ε$ has most often been based on numerical considerations related to the grid spacing used in LES. However, especially for finely resolved LES with grid spacings on the order of or smaller than the chord-length of the blade, the best choice of $ε$ is not known. In this work, a theoretical approach is followed to determine the most suitable value of $ε$. Firstly, we develop an analytical solution to the linearized flow response to a Gaussian lift and drag force and use the results to establish a relationship between the local and far-field velocity required to specify lift and drag forces. Then, focusing first on the lift force, we find $ε$ and the force center location that minimize the square difference between the velocity fields induced by the Gaussian force and 2D potential flow over Joukowski airfoils. We find that the optimal smoothing width $ε^{\rm opt}$ is on the order of 14-25\% of the chord length of the blade, and the center of force is located at about 13-26\% downstream of the leading edge of the blade, for the cases considered. These optimal values do not depend on angle of attack and depend only weakly on the type of lifting surface. To represent the drag force, the optimal width of the circular Gaussian drag force field is shown to be equal to the momentum thickness of the wake.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wake Structure of Wind Turbines in Yaw under Uniform Inflow Conditions

Reducing wake losses in wind farms by deflecting the wakes through turbine yawing has been shown to be a feasible wind farm controls approach. Nonetheless, the effectiveness of yawing depends not only on the degree of wake deflection but also on the resulting shape of the wake. In this work, the deflection and morphology of wakes behind a wind turbine operating in yawed conditions are studied using wind tunnel experiments of a wind turbine modeled as a porous disk in a uniform inflow. First, by measuring velocity distributions at various downstream positions and comparing with prior studies, we confirm that the non-rotating wind turbine model in yaw generates realistic wake deflections. Second, we characterize the wake shape and make first observations of what is termed a curled wake, displaying significant spanwise asymmetry. The wake curling observed in the experiments is also reproduced qualitatively in large eddy simulations using both actuator disk and actuator line models. When a wind turbine is yawed for the benefit of downstream turbines, the asymmetric shape of the wake must be taken into account since it affects how much of it intersects the downstream turbines.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wind Turbine Large-Eddy Simulations on Very Coarse Grid Resolutions using an Actuator Line Model

In this work the accuracy of the Actuator Line Model (ALM) in Large Eddy Simulations of wind turbine flow is studied under the specific conditions of very coarse spatial resolutions. For finely-resolved conditions, it is known that ALM provides better accuracy compared to the standard Actuator Disk Model (ADM) without rotation. However, we show here that on very coarse resolutions, flow induction occurring at rotor scales can affect the predicted inflow angle and can adversely affect the ALM predictions. We first provide an illustration of coarse LES to reproduce wind tunnel measurements. The resulting flow predictions are good, but the challenges in predicting power outputs from the detailed ALM motivate more detailed analysis on a case with uniform inflow. We present a theoretical framework to compare the filtered quantities that enter the Large-Eddy Simulation equations as body forces with a scaling relation between the filtered and unfiltered quantities. The study aims to apply the theoretical derivation to the simulation framework and improve the current results for an ALM, especially in the near wake where the largest differences are observed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupled wake boundary layer model of wind-farms

We present and test the coupled wake boundary layer (CWBL) model that describes the distribution of the power output in a wind-farm. The model couples the traditional, industry-standard wake model approach with a "top-down" model for the overall wind-farm boundary layer structure. This wake model captures the effect of turbine positioning, while the "top-down" portion of the model adds the interactions between the wind-turbine wakes and the atmospheric boundary layer. Each portion of the model requires specification of a parameter that is not known a-priori. For the wake model, the wake expansion coefficient is required, while the "top-down" model requires an effective spanwise turbine spacing within which the model's momentum balance is relevant. The wake expansion coefficient is obtained by matching the predicted mean velocity at the turbine from both approaches, while the effective spanwise turbine spacing depends on turbine positioning and thus can be determined from the wake model. Coupling of the constitutive components of the CWBL model is achieved by iterating these parameters until convergence is reached. We illustrate the performance of the model by applying it to both developing wind-farms including entrance effects and to fully developed (deep-array) conditions. Comparisons of the CWBL model predictions with results from a suite of large eddy simulations (LES) shows that the model closely represents the results obtained in these high-fidelity numerical simulations. A comparison with measured power degradation at the Horns Rev and Nysted wind-farms shows that the model can also be successfully applied to real wind-farms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Deformation statistics of sub-Kolmogorov-scale ellipsoidal neutrally buoyant drops in isotropic turbulence

Small droplets in turbulent flows can undergo highly variable deformations and orientational dynamics. For neutrally buoyant droplets smaller than the Kolmogorov scale, the dominant effects from the surrounding turbulent flow arise through Lagrangian time histories of the velocity gradient tensor. Here we study the evolution of representative droplets using a model that includes rotation and stretching effects from the surrounding fluid, and restoration effects from surface tension including a constant droplet volume constraint, while assuming that the droplets maintain an ellipsoidal shape. The model is combined with Lagrangian time histories of the velocity gradient tensor extracted from DNS of turbulence to obtain simulated droplet evolutions. These are used to characterize the size, shape and orientation statistics of small droplets in turbulence. A critical capillary number, $Ca_c$ is identified associated with unbounded growth of one or two of the droplet's semi-axes. Exploiting analogies with dynamics of polymers in turbulence, the $Ca_c$ number can be predicted based on the large deviation theory for the largest Finite Time Lyapunov exponent. Also, for sub-critical $Ca$ the theory enables predictions of the slope of the power-law tails of droplet size distributions in turbulence. For cases when the viscosities of droplet and outer fluid differ in a way that enables vorticity to decorrelate the shape from the straining directions, the large deviation formalism based on the stretching properties of the velocity gradient tensor loses validity and its predictions fail. Even considering the limitations of the assumed ellipsoidal droplet shape, the results highlight the complex coupling between droplet deformation, orientation and the local fluid velocity gradient tensor to be expected when small viscous drops interact with turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spatio-temporal spectra in the logarithmic layer of wall turbulence: large-eddy simulations and simple models

Motivated by the need to characterize the spatio-temporal structure of turbulence in wall-bounded flows, we study wavenumber-frequency spectra of the streamwise velocity component based on large-eddy simulation (LES) data. The LES data are used to measure spectra as a function of the two wall-parallel wavenumbers and the frequency in the equilibrium (logarithmic) layer. We then reformulate one of the simplest models that is able to reproduce the observations: the random sweeping model with a Gaussian large-scale fluctuating velocity and with additional mean flow. Comparison with LES data shows that the model captures the observed temporal decorrelation, which is related to the Doppler broadening of frequencies. We furthermore introduce a parameterization for the entire wavenumber-frequency spectrum $E_{11}(k_1,k_2,ω;z)$, where $k_1$, $k_2$ are the streamwise and spanwise wavenumbers, $ω$ is the frequency and $z$ is the distance to the wall. The results are found to be in good agreement with LES data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A concurrent precursor inflow method for Large Eddy Simulations and applications to finite length wind farms

In order to enable simulations of developing wind turbine array boundary layers with highly realistic inflow conditions a concurrent precursor method for Large Eddy Simulations is proposed. In this method we consider two domains simultaneously, i.e. in one domain a turbulent Atmospheric Boundary Layer (ABL) without wind turbines is simulated in order to generate the turbulent inflow conditions for a second domain in which the wind turbines are placed. The benefit of this approach is that a) it avoids the need for large databases in which the turbulent inflow conditions are stored and the correspondingly slow I/O operations and b) we are sure that the simulations are not negatively affected by statically swept fixed inflow fields or synthetic fields lacking the proper ABL coherent structures. Sample applications are presented, in which, in agreement with field data a strong decrease of the power output of downstream wind-turbines with respect to the first row of wind-turbines is observed for perfectly aligned inflow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Large Eddy Simulation studies of the effects of alignment and wind farm length

Large eddy simulations of wind farms are performed to study the effects of wind turbine row alignment with respect to the incoming flow direction. Various wind farms with fixed stream-wise spacing (7.85 rotor diameters) and varying lateral displacements and span-wise turbine spacings are considered, for a fixed inflow direction. Simulations show that, contrary to common belief, a perfectly staggered (checker-board) configuration does not necessarily give the highest average power output. Instead, the highest mean wind farm power output is found to depend on several factors, the most important one being the alignment that leads to minimization of wake effects from turbines in several upstream rows. This alignment typically occurs at significantly smaller angles than those corresponding to perfect staggering. The observed trends have implications for wind farm designs, especially in sites with a well-defined prevailing wind direction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Large-eddy simulation study of the logarithmic law for second and higher-order moments in turbulent wall-bounded flow

The logarithmic law for the mean velocity in turbulent boundary layers has long provided a valuable and robust reference for comparison with theories, models, and large-eddy simulations (LES) of wall-bounded turbulence. More recently, analysis of high-Reynolds number experimental boundary layer data has shown that also the variance and higher-order moments of the streamwise velocity fluctuations $u'^{+}$ display logarithmic laws. Such experimental observations motivate the question whether LES can accurately reproduce the variance and the higher-order moments, in particular their logarithmic dependency on distance to the wall. In this study we perform LES of very high Reynolds number wall-modeled channel flow and focus on profiles of variance and higher-order moments of the streamwise velocity fluctuations. In agreement with the experimental data, we observe an approximately logarithmic law for the variance in the LES, with a `Townsend-Perry' constant of $A_1\approx 1.25$. The LES also yields approximate logarithmic laws for the higher-order moments of the streamwise velocity. Good agreement is found between $A_p$, the generalized `Townsend-Perry' constants for moments of order $2p$, from experiments and simulations. Both are indicative of sub-Gaussian behavior of the streamwise velocity fluctuations. The near-wall behavior of the variance, the ranges of validity of the logarithmic law and in particular possible dependencies on characteristic length scales such as the roughness scale $z_0$, the LES grid scale $Δ$, and sub-grid scale (SGS) mixing length $C_sΔ$ are examined. We also present LES results on moments of spanwise and wall-normal fluctuations of velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pressure Hessian and viscous contributions to velocity gradient statistics based on Gaussian random fields

Understanding the non-local pressure contributions and viscous effects on the small-scale statistics remains one of the central challenges in the study of homogeneous isotropic turbulence. Here we address this issue by studying the impact of the pressure Hessian as well as viscous diffusion on the statistics of the velocity gradient tensor in the framework of an exact statistical evolution equation. This evolution equation shares similarities with earlier phenomenological models for the Lagrangian velocity gradient tensor evolution, yet constitutes the starting point for a systematic study of the unclosed pressure Hessian and viscous diffusion terms. Based on the assumption of incompressible Gaussian velocity fields, closed expressions are obtained as the results of an evaluation of the characteristic functionals. The benefits and shortcomings of this Gaussian closure are discussed, and a generalization is proposed based on results from direct numerical simulations. This enhanced Gaussian closure yields, for example, insights on how the pressure Hessian prevents the finite-time singularity induced by the local self-amplification and how its interaction with viscous effects leads to the characteristic strain skewness phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Standard logarithmic mean velocity distribution in a band-limited restricted nonlinear model of turbulent flow in a half-channel

Numerical simulations of wall-turbulence using the restricted nonlinear (RNL) model generate realistic mean velocity profiles in plane Couette and channel flow at low Reynolds numbers. The results are less accurate at higher Re, and while a logarithmic region is observed, its von-Kármán constant is not consistent with the standard logarithmic law. In half-channel flow we show that limiting the streamwise-varying wavenumber support of RNL turbulence to one or few empirically determined modes improves its predictions considerably. In particular, the mean velocity profiles obtained with the band-limited RNL model follow standard logarithmic behavior for the higher Reynolds numbers in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temporal structure of aggregate power fluctuations in large-eddy simulations of extended wind-farms

Fluctuations represent a major challenge for the incorporation of electric power from large wind-farms into power grids. Wind farm power output fluctuates strongly in time, over various time scales. Understanding these fluctuations, especially their spatio-temporal characteristics, is particularly important for the design of backup power systems that must be readily available in conjunction with wind-farms. In this work we analyze the power fluctuations associated with the wind-input variability at scales between minutes to several hours, using large eddy simulations (LES) of extended wind-parks, interacting with the atmospheric boundary layer. LES studies enable careful control of parameters and availability of wind-velocities simultaneously across the entire wind-farm. The present study focuses on neutral atmospheric conditions and flat terrain, using actuator-disk representations of the individual wind-turbines. We consider power from various aggregates of wind-turbines such as the total average power signal, or signals from sub-averages within the wind-farm. Non-trivial correlations are observed due to the complex interactions between turbines placed downstream of each other, and they lead to noticeable spectral peaks at frequencies associated with the inter-turbine spacings when the wind-direction is completely fixed. In that case we observe that the frequency spectra of the total wind-farm output show a decay that follows approximately a $-5/3$ power-law scaling regime, qualitatively consistent with some observations made in field-scale operational wind-parks (Apt, 2007). We find that these features are still observed when the wind-speed varies in magnitude. However, significant changes in the wind-direction over time tend to smooth out the observed spectral peak and reduce the extent of the observed $-5/3$ power-law.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orientation dynamics of small, triaxial-ellipsoidal particles in isotropic turbulence

The orientation dynamics of small anisotropic tracer particles in turbulent flows is studied using direct numerical simulation (DNS) and results are compared with Lagrangian stochastic models. Generalizing earlier analysis for axisymmetric ellipsoidal particles (Parsa et al. 2012), we measure the orientation statistics and rotation rates of general, triaxial ellipsoidal tracer particles using Lagrangian tracking in DNS of isotropic turbulence. Triaxial ellipsoids that are very long in one direction, very thin in another, and of intermediate size in the third direction exhibit reduced rotation rates that are similar to those of rods in the ellipsoid's longest direction, while exhibiting increased rotation rates that are similar to those of axisymmetric discs in the thinnest direction. DNS results differ significantly from the case when the particle orientations are assumed to be statistically independent from the velocity gradient tensor. They are also different from predictions of a Gaussian process for the velocity gradient tensor, which does not provide realistic preferred vorticity-strain-rate tensor alignments. DNS results are also compared with a stochastic model for the velocity gradient tensor based on the recent fluid deformation approximation (RFDA). Unlike the Gaussian model, the stochastic model accurately predicts the reduction in rotation rate in the longest direction of triaxial ellipsoids since this direction aligns with the flow's vorticity, with its rotation perpendicular to the vorticity being reduced. For disc-like particles, or in directions perpendicular to the longest direction in triaxial particles, the model predicts {noticeably} smaller rotation rates than those observed in DNS, a behavior that can be understood based on the probability of vorticity orientation with the most contracting strain-rate eigen-direction in the model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flow visualization using momentum and energy transport tubes and applications to turbulent flow in wind farms

As a generalization of the mass-flux based classical stream-tube, the concept of momentum and energy transport tubes is discussed as a flow visualization tool. These transport tubes have the property, respectively, that no fluxes of momentum or energy exist over their respective tube mantles. As an example application using data from large-eddy simulation, such tubes are visualized for the mean-flow structure of turbulent flow in large wind farms, in fully developed wind-turbine-array boundary layers. The three-dimensional organization of energy transport tubes changes considerably when turbine spacings are varied, enabling the visualization of the path taken by the kinetic energy flux that is ultimately available at any given turbine within the array.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lagrangian time correlations of vorticity alignments in isotropic turbulence: observations and model predictions

Motivated by results from recent particle tracking experiments in turbulence (Xu et al., Nat. Phys. 7, 709 (2011)), we study the Lagrangian time correlations of vorticity alignments with the three eigenvectors of the deformation-rate tensor. We use data from direct numerical simulations (DNS), and explore the predictions of a Lagrangian model for the velocity gradient tensor. We find that the initial increase of correlation of vorticity direction with the most extensive eigen-direction observed by Xu et al. is reproduced accurately using the Lagrangian model, as well as the evolution of correlation with the other two eigendirections. Conversely, time correlations of vorticity direction with the eigen-frame of the pressure Hessian tensor show differences with the model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Local and nonlocal pressure Hessian effects in real and synthetic fluid turbulence

The Lagrangian dynamics of the velocity gradient tensor A in isotropic and homogeneous turbulence depend on the joint action of the self-streching term and the pressure Hessian. Existing closures for pressure effects in terms of A are unable to reproduce one important statistical role played by the anisotropic part of the pressure Hessian, namely the redistribution of the probabilities towards enstrophy production dominated regions. As a step towards elucidating the required properties of closures, we study several synthetic velocity fields and how well they reproduce anisotropic pressure effects. It is found that synthetic (i) Gaussian, (ii) Multifractal and (iii) Minimal Turnover Lagrangian Map (MTLM) incompressible velocity fields reproduce many features of real pressure fields that are obtained from numerical simulations of the Navier Stokes equations, including the redistribution towards enstrophy-production regions. The synthetic fields include both spatially local, and nonlocal, anisotropic pressure effects. However, we show that the local effects appear to be the most important ones: by assuming that the pressure Hessian is local in space, an expression in terms of the Hessian of the second invariant Q of the velocity gradient tensor can be obtained. This term is found to be well correlated with the true pressure Hessian both in terms of eigenvalue magnitudes and eigenvector alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lagrangian Refined Kolmogorov Similarity Hypothesis for Gradient Time-evolution in Turbulent Flows

We study the time evolution of velocity and pressure gradients in isotropic turbulence, by quantifying their decorrelation time scales as one follows fluid particles in the flow. The Lagrangian analysis uses data in a public database generated using direct numerical simulation of the Naiver-Stokes equations, at a Reynolds number 430. It is confirmed that when averaging over the entire domain, correlation functions decay on timescales on the order of the mean Kolmogorov turnover time scale, computed from the globally averaged rate of dissipation and viscosity. However, when performing the analysis in different subregions of the flow, turbulence intermittency leads to large spatial variability in the decay time scales. Remarkably, excellent collapse of the auto-correlation functions is recovered when using the `local Kolmogorov time-scale' defined using the locally averaged, rather than the global, dissipation-rate. This provides new evidence for the validity of Kolmogorov's Refined Similarity Hypothesis, but from a Lagrangian viewpoint that provides a natural frame to describe the dynamical time evolution of turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Recent Fluid Deformation closure for velocity gradient tensor dynamics in turbulence: time-scale effects and expansions

In order to model pressure and viscous terms in the equation for the Lagrangian dynamics of the velocity gradient tensor in turbulent flows, Chevillard & Meneveau (Phys. Rev. Lett. 97, 174501, 2006) introduced the Recent Fluid Deformation closure. Using matrix exponentials, the closure allows to overcome the unphysical finite-time blow-up of the well-known Restricted Euler model. However, it also requires the specification of a decorrelation time scale of the velocity gradient along the Lagrangian evolution, and when the latter is chosen too short (or, equivalently, the Reynolds number is too high), the model leads to unphysical statistics. In the present paper, we explore the limitations of this closure by means of numerical experiments and analytical considerations. We also study the possible effects of using time-correlated stochastic forcing instead of the previously employed white-noise forcing. Numerical experiments show that reducing the correlation time scale specified in the closure and in the forcing does not lead to a commensurate reduction of the autocorrelation time scale of the predicted evolution of the velocity gradient tensor. This observed inconsistency could explain the unrealistic predictions at increasing Reynolds numbers.We perform a series expansion of the matrix exponentials in powers of the decorrelation time scale, and we compare the full original model with a linearized version. The latter is not able to extend the limits of applicability of the former but allows the model to be cast in terms of a damping term whose sign gives additional information about the stability of the model as function of the second invariant of the velocity gradient tensor.

0 citations0 reviewsNo code linkedOpen access