The hybrid Euler–Lagrange procedure using an extension of Moffatt's method

2010 ◽  
Vol 661 ◽  
pp. 45-72 ◽  
Author(s):  
A. M. SOWARD ◽  
P. H. ROBERTS
Keyword(s):  
Fluid Dynamic ◽  
Fluid Motion ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Constant Density ◽  
Mean Flow ◽  
Tensor Calculus ◽  
Inviscid Incompressible Fluid ◽  
Large Eddy ◽  
The Mean

The hybrid Euler–Lagrange (HEL) description of fluid mechanics, pioneered largely by Andrews & McIntyre (J. Fluid Mech., vol. 89, 1978, pp. 609–646), has had to face the fact, in common with all Lagrangian descriptions of fluid motion, that the variables used do not describe conditions at the coordinate x, upon which they depend, but conditions elsewhere at some displaced position xL(x, t) = x + ξ(x, t), generally dependent on time t. To address this issue, we employ ‘Lie dragging’ techniques of general tensor calculus to extend a method introduced by Moffatt (J. Fluid Mech., vol. 166, 1986, pp. 359–378) in the fluid dynamic context, whereby the point x is dragged to xL(x, t) by a ‘fictitious steady flow’ η*(x, t) in a unit of ‘fictitious time’. Whereas ξ(x, t) is a Lagrangian concept intimately linked to the location xL(x, t), the ‘dragging velocity’ η*(x, t) has an essentially Eulerian character, because it describes the fictitious velocity at x itself. For the case of constant-density fluids, we show, using solenoidal η*(x, t) instead of solenoidal ξ(x, t), how the HEL theory can be cast into Eulerian form. A useful aspect of this Eulerian development is that the mean flow itself remains solenoidal, a feature that traditional HEL theories lack. Our method realizes the objective sought by Holm (Physica D, vol. 170, 2002, pp. 253–286) in his derivation of the Navier–Stokes–α equation, which is the basis of one of the methods currently employed to represent the sub-grid scales in large-eddy simulations. His derivation, based on expansion to second order in ξ, contained an error which, when corrected, implied a violation of Kelvin's theorem on the constancy of circulation in inviscid incompressible fluid. We show that this is rectified when the expansion is in η* rather than ξ, Kelvin's theorem then being satisfied to all orders for which the expansion converges. We discuss the implications of our approach using η* for the Navier–Stokes–α theory.

scholarly journals Three-Dimensional Turbulent Vortex Shedding From a Surface-Mounted Square Cylinder: Predictions With Large-Eddy Simulations and URANS

2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi
Keyword(s):  
Experimental Data ◽  
Turbulence Model ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
Large Eddy

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.

Large eddy simulations of variable property turbulent flow in horizontal and vertical channels with buoyancy and heat transfer effects

2007 ◽  
Vol 221 (4) ◽  
pp. 429-441 ◽  
Author(s):  
J S Lee ◽  
R H Pletcher
Keyword(s):  
Heat Transfer ◽  
Large Scale ◽  
Vertical Channel ◽  
Large Eddy Simulations ◽  
Mean Flow ◽  
Mean Velocity ◽  
Velocity Fluctuations ◽  
Combined Flow ◽  
Large Eddy ◽  
The Mean

Turbulent combined flow of forced and natural convection was investigated using large eddy simulations for horizontal and vertical channels with significant heat transfer. The walls were maintained at constant temperatures, one heated and the other cooled, at temperature ratios of 1.01, 1.99, and 3.00, respectively. Results showed that with increasing the Grashof number, large-scale turbulent motions emerged near the wall, resulting in significant changes in turbulent intensities for the horizontal channel flow case. Aiding and opposing flows, however, for the vertical channel, emerge near the heated and cooled walls, respectively, while the pressure gradient drives the mean flow upwards. Buoyancy effects on the mean velocity, temperature, and turbulent intensities were observed near the walls. In the aiding flow, the turbulent intensities and heat transfer were suppressed and the flow was relaminarized at large values of the Grash of number. In the opposing flow, however, turbulence was enhanced with increasing velocity fluctuations.

scholarly journals Large-Eddy Simulations of a Drizzling, Stratocumulus-Topped Marine Boundary Layer

2009 ◽  
Vol 137 (3) ◽  
pp. 1083-1110 ◽  
Author(s):  
Andrew S. Ackerman ◽  
Margreet C. vanZanten ◽  
Bjorn Stevens ◽  
Verica Savic-Jovcic ◽  
Christopher S. Bretherton ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Marine Boundary Layer ◽  
Cloud Water ◽  
Large Eddy Simulations ◽  
Cloud Base ◽  
Entrainment Rate ◽  
Liquid Water Path ◽  
Average Maximum ◽  
Large Eddy ◽  
The Mean

Abstract Cloud water sedimentation and drizzle in a stratocumulus-topped boundary layer are the focus of an intercomparison of large-eddy simulations. The context is an idealized case study of nocturnal stratocumulus under a dry inversion, with embedded pockets of heavily drizzling open cellular convection. Results from 11 groups are used. Two models resolve the size distributions of cloud particles, and the others parameterize cloud water sedimentation and drizzle. For the ensemble of simulations with drizzle and cloud water sedimentation, the mean liquid water path (LWP) is remarkably steady and consistent with the measurements, the mean entrainment rate is at the low end of the measured range, and the ensemble-average maximum vertical wind variance is roughly half that measured. On average, precipitation at the surface and at cloud base is smaller, and the rate of precipitation evaporation greater, than measured. Including drizzle in the simulations reduces convective intensity, increases boundary layer stratification, and decreases LWP for nearly all models. Including cloud water sedimentation substantially decreases entrainment, decreases convective intensity, and increases LWP for most models. In nearly all cases, LWP responds more strongly to cloud water sedimentation than to drizzle. The omission of cloud water sedimentation in simulations is strongly discouraged, regardless of whether or not precipitation is present below cloud base.

Investigation on hydrodynamic forces leading to coarse particle entrainment using Large-Eddy Simulations

2021 ◽  
Author(s):  
Gaston Latessa ◽  
Angela Busse ◽  
Manousos Valyrakis
Keyword(s):  
Experimental Data ◽  
Flow Field ◽  
Large Scale ◽  
Large Eddy Simulations ◽  
Flow Conditions ◽  
Mean Flow ◽  
Protection Measures ◽  
Practical Applications ◽  
Particle Entrainment ◽  
Large Eddy

<p>The prediction of particle motion in a fluid flow environment presents several challenges from the quantification of the forces exerted by the fluid onto the solids -normally with fluctuating behaviour due to turbulence- and the definition of the potential particle entrainment from these actions. An accurate description of these phenomena has many practical applications in local scour definition and to the design of protection measures.</p><p>In the present work, the actions of different flow conditions on sediment particles is investigated with the aim to translate these effects into particle entrainment identification through analytical solid dynamic equations.</p><p>Large Eddy Simulations (LES) are an increasingly practical tool that provide an accurate representation of both the mean flow field and the large-scale turbulent fluctuations. For the present case, the forces exerted by the flow are integrated over the surface of a stationary particle in the streamwise (drag) and vertical (lift) directions, together with the torques around the particle’s centre of mass. These forces are validated against experimental data under the same bed and flow conditions.</p><p>The forces are then compared against threshold values, obtained through theoretical equations of simple motions such as rolling without sliding. Thus, the frequency of entrainment is related to the different flow conditions in good agreement with results from experimental sediment entrainment research.</p><p>A thorough monitoring of the velocity flow field on several locations is carried out to determine the relationships between velocity time series at several locations around the particle and the forces acting on its surface. These results a relevant to determine ideal locations for flow investigation both in numerical and physical experiments.</p><p>Through numerical experiments, a large number of flow conditions were simulated obtaining a full set of actions over a fixed particle sitting on a smooth bed. These actions were translated into potential particle entrainment events and validated against experimental data. Future work will present the coupling of these LES models with Discrete Element Method (DEM) models to verify the entrainment phenomena entirely from a numerical perspective.</p>

scholarly journals Stability and Sensitivity Analysis of Hydrodynamic Instabilities in Industrial Swirled Injection Systems

2017 ◽  
Author(s):  
Thomas L. Kaiser ◽  
Thierry Poinsot ◽  
Kilian Oberleithner
Keyword(s):  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Vortex Core ◽  
Large Eddy Simulations ◽  
Mode Shapes ◽  
Mean Flow ◽  
Precessing Vortex Core ◽  
Large Eddy ◽  
Good Agreement

The hydrodynamic instability in an industrial, two-staged, counter-rotative, swirled injector of highly complex geometry is under investigation. Large eddy simulations show that the complicated and strongly nonparallel flow field in the injector is superimposed by a strong precessing vortex core. Mean flow fields of large eddy simulations, validated by experimental particle image velocimetry measurements are used as input for both local and global linear stability analysis. It is shown that the origin of the instability is located at the exit plane of the primary injector. Mode shapes of both global and local linear stability analysis are compared to a dynamic mode decomposition based on large eddy simulation snapshots, showing good agreement. The estimated frequencies for the instability are in good agreement with both the experiment and the simulation. Furthermore, the adjoint mode shapes retrieved by the global approach are used to find the best location for periodic forcing in order to control the precessing vortex core.

Disruption of the bottom log layer in large-eddy simulations of full-depth Langmuir circulation

2012 ◽  
Vol 699 ◽  
pp. 79-93 ◽  
Author(s):  
A. E. Tejada-Martínez ◽  
C. E. Grosch ◽  
N. Sinha ◽  
C. Akan ◽  
G. Martinat
Keyword(s):  
Gravity Waves ◽  
Large Eddy Simulations ◽  
Stokes Drift ◽  
Navier Stokes ◽  
Wake Region ◽  
Langmuir Circulation ◽  
Mean Velocity ◽  
Wave Forcing ◽  
Large Eddy ◽  
Shear Current

AbstractWe report on disruption of the log layer in the resolved bottom boundary layer in large-eddy simulations (LES) of full-depth Langmuir circulation (LC) in a wind-driven shear current in neutrally-stratified shallow water. LC consists of parallel counter-rotating vortices that are aligned roughly in the direction of the wind and are generated by the interaction of the wind-driven shear with the Stokes drift velocity induced by surface gravity waves. The disruption is analysed in terms of mean velocity, budgets of turbulent kinetic energy (TKE) and budgets of TKE components. For example, in terms of mean velocity, the mixing due to LC induces a large wake region eroding the classical log-law profile within the range $90\lt { x}_{3}^{+ } \lt 200$. The dependence of this disruption on wind and wave forcing conditions is investigated. Results indicate that the amount of disruption is primarily determined by the wavelength of the surface waves generating LC. These results have important implications for turbulence parameterizations for Reynolds-averaged Navier–Stokes simulations of the coastal ocean.

scholarly journals Flat vs. curved rigid-lid LES computations of an open-channel confluence

2019 ◽  
Vol 21 (2) ◽  
pp. 318-334 ◽  
Author(s):  
Pedro Xavier Ramos ◽  
Laurent Schindfessel ◽  
João Pedro Pêgo ◽  
Tom De Mulder
Keyword(s):  
Free Surface ◽  
Secondary Flow ◽  
Open Channel ◽  
Pressure Field ◽  
Large Eddy Simulations ◽  
Numerical Treatment ◽  
Mean Flow ◽  
Large Eddy ◽  
Flow Case ◽  
Channel Confluence

Abstract This paper describes the application of four Large Eddy Simulations (LES) to an open-channel confluence flow, making use of a frictionless rigid-lid to treat the free-surface. Three simulations are conducted with a flat rigid-lid, at different elevations. A fourth simulation is carried out with a curved rigid-lid which is a closer approximation to the real free-surface of the flow. The curved rigid-lid is obtained from the time-averaged pressure field on the flat rigid-lid from one of the initial three simulations. The aim is to investigate the limitations of the free-surface treatment by means of a rigid-lid in the simulation of an asymmetric confluence, showing the differences that both approaches produce in terms of mean flow, secondary flow and turbulence. After validation with experimental data, the predictions are used to understand the differences between adopting a flat and a curved rigid-lid onto the confluence hydrodynamics. For the present flow case, although it was characterized by a moderately low downstream Froude number (Fr ≈ 0.37), it was found that an oversimplification of the numerical treatment of the free-surface leads to a decreased accuracy of the predictions of the secondary flow and turbulent kinetic energy.

Sign in / Sign up

Export Citation Format

Share Document