Wall-resolved large-eddy simulation of turbulent channel flows with rough walls

Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.

LARGE EDDY SIMULATION OF FULLY DEVELOPED TURBULENT FLOW IN A CURVED CHANNEL

In this paper large eddy simulation of the fully developed turbulent flow in a curved channel is carried out.The computational results are presented and compared with the experimental results of Eskinazi and Yeh.It is shown that the numerical results of the present LES are reliable and the influence of the curvature on the turbulence feature is correctly revealed.

Large-eddy simulation of circular cylinder flow at subcritical Reynolds number:Turbulent wake and sound radiation

The flows past a number 3900 are simulated circular cylinder at Reynolds using large-eddy simulation （LES） and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier- Stokes equations. The dynamic global coefficient version of the Vreman＇s subgrid scale （SGS） model is used to com- pute the sub-grid stresses. Curie＇s integral of Lighthill＇s acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctua- tions obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the -5/3 power law. It is shown that Vreman＇s SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.

Large Eddy Simulation of Stratified Turbulent Channel Flow with a Dynamic Subgrid Scale Model

In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.

Large Eddy Simulation of Particle Wake Effect and RANS Simulation of Turbulence Modulation in Gas-Particle Flows

The turbulence enhancement by particle wake effect is studied by large eddy simulation （LES） of turbulent gas flows passing a single particle. The predicted time-averaged and root-mean-square fluctuation velocities behind the particle are in agreement with the Reynolds-averaged Navier-Stokes modeling results and experimental results. A semi-empirical turbulence enhancement model is proposed by the present-authors based on the LES resuits. This model is incorporated into the second-order moment two-phase turbulence model for simulating vertical gas-particle pipe flows and horizontal gas-particle channel flows. The simulation results show that compared with the model not accounting for the particle wake effect, the present model gives simulation results for the gas turbulence modulation in much better agreement with the experimental results.

Large eddy simulation of vertical turbulent jets under JONSWAP waves

The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.

Investigation of Wall Pressure Fluctuations in a Turbulent Boundary Layer by Large Eddy Simulation

Large eddy simulation (LES) was used to investigate the space-time field of the low Mach number, fully developed turbulent boundary layer on a smooth, rigid flat plate. The wall-pressure field simulated by LES was analyzed to obtain the pressure statistics, including the wall-pressure root-mean square, skewness and flatness factors, which show the wall pressure distribution was not Gaussian. The profile of the auto-power spectral density and the contour of the streamwise wavenumber-frequency spectral density of wall-pressure were plotted. The "convection ridge" can be observed clearly and the convection velocity can be calculated from the location of the convection peak.

Large Eddy Simulations of Rayleigh-Bénard Convection in a Cubic Cavity

This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.

Large eddy simulations of zinc ions transfer to turbulent flows from hyporheic zone

Metal contaminants from surface water pollution events often enter hyporheic zones,under certain conditions,they may be released back into streams,causing secondary pollution to the water quality.The present study investigated the effects of adsorption,permeability,and anisotropy of sediment beds on the release of zinc ions(Zn2+)from the hyporheic zone into overlying turbulent flows using large-eddy simulations(LES).The volume-averaged Navier-Stokes equations and advection-diffusion equation with adsorption term were used to describe the sediment in-flow,adsorption,and convective diffusion of Zn2+within the sediment layer.The effects of sediment permeability on the Zn2+concentration distribution and mass transfer processes were investigated by time-averaged statistics of flow and concentration fields.The results show that adsorption becomes stronger as the pH value increases,leading to a slow increase in Zn2+concentration in the overlying water layer and reaching a lower steady-state concentration.Higher overall permeability of the sediment layer can enhance mass and momentum exchange near the sediment-water interface(SWI),and intensify the release of Zn2+from the sediment layer into the overlying water.As the wall-normal permeability of the sediment layer increases,the normal turbulent intensity strengthens,momentum transport enhances,the wall-normal Zn2+concentration flux increases,the effective diffusion coefficient increases,and the concentration in the overlying water increases.

Investigation on the width-to-depth ratio effect on turbulent flows in a sharp meandering channel with periodic boundaries using large eddy simulations

As one of the most common river patterns in nature,meandering river has very complex flow structures in its curved channel bends,including secondary flow structure and primary flow velocity redistributions.To date,most of the studies have been carried out on the flow structures in channel bends with unavoidable influences from inlet and outlet boundaries,while a streamwise periodic boundary can overcome this shortcoming elegantly.In this paper,large eddy simulations(LES)are employed to investigate the complex flow structures in periodically continuous sharp sine-generated bends.The influence of width-to-depth ratios and dimensionless curvature radiuses are studied.The results highlight two additional vortex structures beyond the commonly known secondary currents:The recirculation zone(RZ)and the inner bank cell(IBC).The width-to-depth ratio shows the determining effect on the recirculation zone.The size of recirculation zone is usually bigger in sine-generated-curve(SGC)channel with large width-to-depth ratios.The biggest recirculation zones appear between the zero-curvature section and the apex section.The inner bank cell only forms in SGC channels with small width-to-depth ratios and low curvature.For SGC channel with large width-to-depth ratios,only one circulation cell is observed near the inner bank.The spatial variations of turbulent features are also revealed by statistical analysis based on the LES sampling data.Results highlight remarkable effect of width-to-depth ratio and dimensionless curvature radius on the turbulent kinetic energy(TKE)and bed shear stress in SGC channels.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

Large Eddy Simulation of Liquid Flow in a Stirred Tank with Improved Inner-Outer Iterative Algorithm

In this study, the large eddy simulation technique was applied on the flow in a baffled stirred tank driven by a Rushton turbine at Re=29000. The interaction between the rotating impeller and the static baffles was accounted for by means of the improved inner-outer iterative algorithm. The sub-grid scale model was a conventional Smagorinsky model. The numerical solution of the governing equations was conducted in a cylindrical staggered grid. The momentum and the continuity equations were discretized using the finite difference method, with a third-order QUICK scheme used for convective terms. The phase-resolved predictions were compared with the experimental data of Wu and Patterson and good agreement was observed for both the mean and the turbulence quantities. They were much better than the Reynolds-averaged Navier-Stokes model including the Reynolds Stress Model for simulating the turbulence. The study also suggests the feasibility of LES in combination with the improved inner-outer iterative algorithm for the simulation of turbulent flow in stirred tanks.

Numerical Study on Flow-induced Noise for a Steam Stop-valve Using Large Eddy Simulation

The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.

LARGE-EDDY SIMULATION OF STRATIFIED CHANNEL FLOW

An investigation of large-eddy simulation(LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smagorinsky eddy-viscosity model and Yoshizawa eddy-viscosity model were used to describe the unresolved subgrid scale (SGS) fluctuations respectively After some numerical testing, the latter was further simplified so that if can be used in the dynamic model closure. A LES code was developed for parallel computations by using the parallel technique, and aas run on the Dawn-1000 parallel computer. To demonstrate the viability and accuracy of the code, our results are compared with and found in good agreement with available LES results.

Large Eddy Simulations of Flow Instabilities in a Stirred Tank Gen-erated by a Rushton Turbine

The aim of this work is to investigate the flow instabilities in a baffled, stirred tank generated by a single Rushton turbine by means of large eddy simulation （LES）. The sliding mesh method was used for the coupling between the rotating and the stationary frame of references. The calculations were carried out on the ＂Shengcao-21C＂ supercomputer using a computational fluid dynamics （CFD） code CFX5. The flow fields predicted by the LES simulation and the simulation using standard κ-ε model were compared to the results from particle image velocimetry （PIV） measurements. It is shown that the CFD simulations using the LES approach and the standard κ-ε model agree well with the PIV measurements. Fluctuations of the radial and axial velocity are predicted at different frequencies by the LES simulation. Velocity fluctuations of high frequencies are seen in the impeller region, while low frequencies velocity fluctuations are observed in the bulk flow. A low frequency velocity fluctuation with a nondimensional frequency of 0.027Hz is predicted by the LES simulation, which agrees with experimental investigations in the literature. Flow circulation patterns predicted by the LES simulation are asymmetric, stochastic and complex, spanning a large portion of the tanks and varying with time, while circulation patterns calculated by the simulation using the standard κ-ε model are symmetric. The results of the present work give better understanding to the flow instabilities in the mechanically agitated tank. However, further analysis of the LES calculated velocity series by means of fast Fourier transform （FFT） and/or spectra analysis are recommended in future work in order to gain more knowledge of the complicated flow phenomena.

Large eddy simulation of turbulence in ocean surface boundary layer at Zhangzi Island offshore

This study uses a large eddy simulation （LES） model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island （39°N, 122°E） during July 2009 are used to drive the LES model. The LES results capture a clear diurnal cycle in the oceanic turbulence boundary layer. The process of the heat penetration and heat distribution characteristics are analyzed through the heat flux results from the LES and their differences between two diurnal cycles are discussed as well. Energy balance and other dynamics are investigated which show that the tide-induced shear production is the main source of the turbulence energy that balanced dissipation. Momentum flux near the surface shows better agreement with atmospheric data computed by the eddy correlation method than those computed by bulk formula.

Fourier neural operator approach to large eddy simulation of three-dimensional turbulence

Fourier neural operator(FNO)model is developed for large eddy simulation(LES)of three-dimensional(3D)turbulence.Velocity fields of isotropic turbulence generated by direct numerical simulation(DNS)are used for training the FNO model to predict the filtered velocity field at a given time.The input of the FNO model is the filtered velocity fields at the previous several time-nodes with large time lag.In the a posteriori study of LES,the FNO model performs better than the dynamic Smagorinsky model(DSM)and the dynamic mixed model(DMM)in the prediction of the velocity spectrum,probability density functions(PDFs)of vorticity and velocity increments,and the instantaneous flow structures.Moreover,the proposed model can significantly reduce the computational cost,and can be well generalized to LES of turbulence at higher Taylor-Reynolds numbers.