A dynamic global-coefficient mixed subgrid-scale model for large-eddy simulation of turbulent flows

A dynamic global-coefficient mixed subgrid-scale eddy-viscosity model for large-eddy simulation of turbulent flows in complex geometries is developed. In the present model, the subgrid-scale stress is decomposed into the modified Leonard stress, cross stress, and subgrid-scale Reynolds stress. The modified Leonard stress is explicitly computed assuming a scale similarity, while the cross stress and the subgrid-scale Reynolds stress are modeled using the global-coefficient eddy-viscosity model. The model coefficient is determined by a dynamic procedure based on the global-equilibrium between the subgrid-scale dissipation and the viscous dissipation. The new model relieves some of the difficulties associated with an eddy-viscosity closure, such as the nonalignment of the principal axes of the subgrid-scale stress tensor and the strain rate tensor and the anisotropy of turbulent flow fields, while, like other dynamic global-coefficient models, it does not require averaging or clipping of the model coefficient for numerical stabilization. The combination of the global-coefficient eddy-viscosity model and a scale-similarity model is demonstrated to produce improved predictions in a number of turbulent flow simulations.     

Structural subgrid-scale modeling for large-eddy simulation: A review

Accurately modeling nonlinear interactions in turbulence is one of the key challenges for large-eddy simu-lation (LES) of turbulence. In this article, we review recent studies on structural subgrid scale modeling, focusing on evaluating how well these models predict the effects of small scales. The article discusses a priori and a posteriori test results. Other nonlinear models are briefly discussed, and future prospects are noted.     

Direct numerical simulation of a turbulent axisymmetric jet with buoyancy induced acceleration

Direct numerical simulations of an axisymmetric jet with off-source volumetric heat addition are presented in this paper. The system solved here involves a three-way coupling between velocity, concentration and temperature. The computations are performed on a spherical coordinate system, and application of a traction free boundary condition at the lateral edges allows physical entrainment into the computational domain. The Reynolds and Richardson numbers based on local scales employed in the simulations are 1000 and 12 respectively. A strong effect of heat addition on the jet is apparent. Heating causes acceleration of the jet, and an increased dilution due to an increase in entrainment. Further, the streamwise velocity profile is distorted, and the cross-stream velocity is inward for all radial locations for the heated jet. Interestingly, the maximum temperature is realized off-axis and a short distance upstream of the exit of the heat injection zone (HIZ). The temperature width is intermediate between the scalar and velocity widths in the HIZ. Normalized rms of the concentration and temperature increases in the HIZ, whereas that of streamwise, cross-stream and tangential velocities increases rapidly after decreasing. Both mass flux and entrainment are larger for the heated jet as compared to their unheated counterparts. The buoyancy flux increases monotonically in the HIZ, and subsequently remains constant.     

Evaluation of subgrid-scale models in large-eddy simulation of flow past a two-dimensional block

Large-eddy simulations of flow past a two-dimensional (2D) block were performed to evaluate four subgrid-scale (SGS) models: (i) the traditional Smagorinsky model, (ii) the Lagrangian dynamic model, (iii) the Lagrangian scale-dependent dynamic model, and (iv) the modulated gradient model. An immersed boundary method was employed to simulate the 2D block boundaries on a uniform Cartesian grid. The sensitivity of the simulation results to grid refinement was investigated by using four different grid resolutions. The velocity streamlines and the vertical profiles of the mean velocities and variances were compared with experimental results. The modulated gradient model shows the best overall agreement with the experimental results among the four SGS models. In particular, the flow recirculation, the reattachment position and the vertical profiles are accurately reproduced with a relative coarse grid resolution of (N_x × N_y × N_z=) 160 × 40 × 160 (n_x × n_z = 13 × 16 covering the block). Besides the modulated gradient model, the Lagrangian scale-dependent dynamic model is also able to give reasonable prediction of the flow statistics with some discrepancies compared with the experimental results. Relatively poor performance by the Lagrangian dynamic model and the Smagorinsky model is observed, with simulated recirculating patterns that differ from the measured ones. Analysis of the turbulence kinetic energy (TKE) budget in this flow shows evidence of a strong production of TKE in the shear layer that forms as the flow is deflected around the block.     

Transient behavior of an axisymmetric turbulent jet

The prediction of the response of unsteady flows submitted to external excitation is a real challenge for the optimization of industrial processes. As the jet flow is a very basic turbulent flow related to mixing and entrainment phenomena via turbulent structure dynamic, we investigate the transient behavior of an axisymmetric jet submitted to a large and sudden decrease of its ejection velocity. The non stationary flow evolution is studied experimentally. Measurements along the jet axis based on pure ensemble averaging show clearly the convective motion of the perturbation and the adaptation of the local interaction to the local jet time scale. A transverse investigation in the non stationary region show that the mean flow and its turbulence is deeply affected during the local velocity decrease.     

Numerical simulation of axisymmetric turbulent jets

The flow in axisymmetric turbulent jets is numerically simulated with the use of a semi-empirical second-order turbulence model including differential transport equations for the normal Reynolds stresses. Calculated results are demonstrated to agree with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 55–60, September–October, 2008.     

Nonisothermal axisymmetric turbulent jet curved by gravity forces

Characteristic of a free turbulent jet which has a material density differing from the density of the surrounding medium is a variation of the momentum along the stream, caused by the action of the Archimedes lift force. The present study makes an attempt to find the relations which define the variation of the axial velocities and the differential temperatures for the most general case of nonvertical discharge of a free nonisothermal axisymmetric jet of air, with commensurate values of the inertial and gravitational forces. In solving this problem we shall make use of the same jet scheme used in [1] with the corresponding basic assumptions.     

The entrainment behavior of a turbulent axisymmetric jet in a viscous host fluid

Although turbulent jets have been studied extensively, one configuration that has not received much attention is the viscosity-stratified jet, wherein a turbulent jet of lower viscosity issues into a density-matched host liquid of higher viscosity. We present experimental data for scalar dispersion and two-dimensional velocity measurements in the axial plane of a turbulent axisymmetric jet with a Reynolds number (Re) of 2,000 issuing into a viscous host liquid at viscosity ratios (m) ranging from 1 to 55. The presence of a strong viscosity discontinuity across the jet edge results in a significant decrease in the scalar spread rate. We attribute this to the rapid reduction in turbulence intensity and the suppression of large engulfing eddies at the jet edge. The velocity profile, on the other hand, indicates that the velocity width and mass flux reduce with increasing m up to about 20, but then increase for higher values of m. This non-monotonic variation is explained by the growing influence of viscous stress for m>20. The scalar spread rate, the velocity spread rate, the centerline velocity decay rate, and the jet mass flux are all minimized for m20 for Re=2,000.

---

The development of an axisymmetric curved turbulent wall jet

An experimental study has been carried out of the low speed Coanda wall jet with both streamwise and axisymmetric curvature. A single component laser Doppler technique was used, and by taking several orientations at a given point, values of the three mean velocities and five of the six Reynolds stresses were obtained. The lateral divergence and convex streamwise curvature both enhanced the turbulence in the outer part of the jet compared with a plane two-dimensional wall jet. The inner layer exhibited a large separation of the positions of maximum velocity and zero shear stress. It was found that the streamwise mean velocity profile became established very rapidly downstream of the slot exit. The profile appeared fairly similar at later downstream positions, but the mean radial velocity and turbulence parameters showed the expected nonself preservation of the flow. Removal of the streamwise curvature resulted in a general return of the jet conditions toward those expected of a plane wall jet. The range and accuracy of the data may be used for developing turbulence models and computational techniques for this type of flow.     

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.     

Constrained large-eddy simulation of turbulent flow over inhomogeneous rough surfaces

In this work we extend the method of the constrained large-eddy simulation(CLES) to simulate the turbulent flow over inhomogeneous rough walls. In the original concept of CLES, the subgrid-scale(SGS) stress is constrained so that the mean part and the fluctuation part of the SGS stress can be modelled separately to improve the accuracy of the simulation result. Here in the simulation of the rough-wall flows, we propose to interpret the extra stress terms in the CLES formulation as the roughness-induced stress so that the roughness inhomogeneity can be incorporated by modifying the formulation of the constrained SGS stress. This is examined with the simulations of the channel flow with the spanwise alternating high/low roughness strips. Then the CLES method is employed to investigate the temporal response of the turbulence to the change of the wall condition from rough to smooth. We demonstrate that the temporal development of the internal boundary layer is just similar to that in a spatial rough-tosmooth transition process, and the spanwise roughness inhomogeneity has little impact on the transition process.     

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

We analyze the error of large-eddy simulation(LES) in wall pressure fluctuation of a turbulent channel flow. To separate different sources of the error, we conduct both direct numerical simulations(DNS)and LES, and apply an explicit filter on DNS data to obtain filtered DNS(FDNS) data. The error of LES is consequently decomposed into two parts: The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation. The second part is the difference between LES and FDNS induced by the error of LES in velocity field. By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation, it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.     

Temporal large-eddy simulation of unstratified and stably stratified turbulent channel flows

Recently, Pruett et al. [Pruett, C.D., Gatski, T.B., Grosch, C.E., Thacker, W.D., 2003. The temporally filtered Navier–Stokes equations: properties of the residual stress. Phys. Fluids 15, 2127–2140] proposed an approach to large-eddy simulation (LES) based on time-domain filtering; their approach was termed temporal large-eddy simulation or TLES. In a continuation of their work, Pruett and collaborators tested their methodology by successfully performing TLES of unstratified turbulent channel flow up to Reynolds number of 590 (based on channel half-height and friction velocity) [Pruett, C.D., Thomas, B.C., Grosch, C.E., Gatski, T.B., 2006. A temporal approximate deconvolution model for LES. Phys. Fluids 18, 028104, 4p]. Here, we carefully analyze the TLES methodology in order to understand the role of its key components and in the process compare TLES to more traditional approaches of spatial LES. Furthermore, we extend the methodology to stably stratified turbulent channel flow.     

Calculation of an axisymmetric turbulent wall jet over a surface of convex curvature

One of the many applications of curved wall jets of engineering importance is the Coanda Flare, which is used for burning waste gases in the petroleum industry and which gave rise to this work. The gas jet flows over an axisymmetric tulip shaped body, entraining ambient air and so promoting clean combustion. The object of this work was to calculate the development of the jet with the extra rates of strain imposed by both longitudinal curvature and divergence. A differential ‘partially-parabolic’ technique was used with uncoupling of the streamwise and cross-stream momentum equations, leading to an efficient computer program. The extra rates of strain were modelled by corrections to a mixing length model with the two effects being assumed to be additive. The calculation method was compared with seven test ccases of experimental data. The first five were from published literature, and included the plane wall jet and axisymmetric free jet, and the separate effects of longitudinal curvature and divergence. The lost two cases were measurements of the wall jet flow over a model Coanda Flare. The calculation method gave generally good results for the main features of the flow such as growth rate and velocity decay. Details of the flow were not so well predicted, particularly the turbulent shear stress, as a result of the relatively simple turbulence model employed. The calculation method should provide a useful engineering tool, but some profitable developments could be made, particularly in the area of turbulence modelling.     

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45° expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k–? model were solved in a non-orthogonal curvillinear co-ordinate system. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite differencing methods—hybrid, quadratic upwind and skew upwind—were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.     

Implicit large-eddy simulation applied to turbulent channel flow with periodic constrictions

The subgrid-scale (SGS) model in a large-eddy simulation (LES) operates on a range of scales which is marginally resolved by discretization schemes. Accordingly, the discretization scheme and the subgrid-scale model are linked. One can exploit this link by developing discretization methods from subgrid-scale models, or the converse. Approaches where SGS models and numerical discretizations are fully merged are called implicit LES (ILES). Recently, we have proposed a systematic framework for the design, analysis, and optimization of nonlinear discretization schemes for implicit LES. In this framework parameters inherent to the discretization scheme are determined in such a way that the numerical truncation error acts as a physically motivated SGS model. The resulting so-called adaptive local deconvolution method (ALDM) for implicit LES allows for reliable predictions of isotropic forced and decaying turbulence and of unbounded transitional flows for a wide range of Reynolds numbers. In the present paper, ALDM is evaluated for the separated flow through a channel with streamwise-periodic constrictions at two Reynolds numbers Re = 2,808 and Re = 10,595. We demonstrate that, although model parameters of ALDM have been determined for isotropic turbulence at infinite Reynolds number, it successfully predicts mean flow and turbulence statistics in the considered physically complex, anisotropic, and inhomogeneous flow regime. It is shown that the implicit model performs at least as well as an established explicit model.      

Numerical simulation of turbulent impinging jet on a rotating disk

The calculations of quasi‐three‐dimensional momentum equations were carried out to study the influence of wall rotation on the characteristics of an impinging jet. The pressure coefficient, the mean velocity distributions and the components of Reynolds stress are calculated. The flow is assumed to be steady, incompressible and turbulent. The finite volume scheme is used to solve the continuity equation, momentum equations and k–ε model equations. The flow characteristics were studied by varying rotation speed ω for 0?ω?167.6 rad/s, the distance from nozzle to disk (H/d) was (3, 5, 8 and 10) and the Reynolds number Re base on V_J and d was 1.45 × 10⁴. The results showed that, the radial velocity and turbulence intensity increase by increasing the rotation speed and decrease in the impingement zone as nozzle to disk spacing increases. When the centrifugal force increases, the radial normal stresses and shear stresses increase. The location of maximum radial velocity decreases as the local velocity ratio (α) increases. The pressure coefficient depends on the centrifugal force and it decreases as the distance from nozzle to plate increases. In impingement zone and radial wall jet, the spread of flow increases as the angular velocity decreases The numerical results give good agreement with the experiment data of Minagawa and Obi (Int. J. of Heat and Fluid Flow 2004; 25 :759–766). Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation and visualization of fiber suspension in a turbulent round jet

The turbulent properties of the fiber suspension in a turbulent round jet are numerically simulated and visualized, and some of the results are compared with the experimental data. The effects of the Reynolds number, fiber volume fraction, and aspect ratio are analyzed. The results show that the fiber injection in the flow has a delay effect on the streamwise velocity decay along the jet axis, and such an effect becomes more obvious with the increases in the fiber volume fraction and aspect ratio and the decrease in the Reynolds number. The flow with fibers shows an increase in the streamwise velocity along the radial direction, and the increase magnitude is directly proportional to the fiber volume fraction and aspect ratio and inversely proportional to the Reynolds number. The presence of fibers makes the turbulent kinetic energy and Reynolds stress increase, and the extent increases with the fiber volume fraction, Reynolds number, and fiber aspect ratio.