High‐order detached eddy simulation,zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth‐order discretization scheme is used and is validated using the test case of a two‐dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS–NLES approach does not suffer from this trait. For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton–Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS–NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. Copyright © 2013 John Wiley & Sons, Ltd.     

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.     

An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview

The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

Explicit algebraic stress modelling of homogeneous and inhomogeneous flows

Capability of the explicit algebraic stress models to predict homogeneous and inhomogeneous shear flows are examined. The importance of the explicit solution of the production to dissipation ratio is first highlighted by examining the algebraic stress models performance at purely irrotational strain conditions. Turbulent recirculating flows within sudden expanding pipes are further simulated with explicit algebraic stress model and anisotropic eddy viscosity model. Both models predict better stress–strain interactions, showing reasonable shear layer developments. The anisotropic stress field are also accurately predicted by the models, though the anisotropic eddy viscosity model of Craft et al. returns marginally better results. Copyright © 2005 John Wiley & Sons, Ltd.     

Modifications for an explicit algebraic stress model

An extension of the explicit algebraic stress model, developed by Gatski and Speziale [Gatski TB, Speziale CG. On the explicit algebraic stress models for complex turbulent flows. Journal of Fluid Mechanics 1993; 254: 59–78] is proposed. The extension implicates some essential characteristics of second‐order closure models. The strain‐dependent coefficients are modified, resulting in an alleviation of the numerical instabilities involved in the model. A new near‐wall damping function f_μ in the eddy viscosity relation is introduced. To enhance dissipation in near‐wall regions, the model constant C_ϵ1 is modified and an extra positive source term is included in the dissipation equation. In addition, a realizable time scale is incorporated to remove the wall singularity. Computed results show that the modified Gatski–Speziale (MGS) model predictions are in better agreement with the direct numerical simulation (DNS) and experimental data than those of the original Gatski–Speziale (OGS) model. Copyright © 2001 John Wiley & Sons, Ltd.     

On the Realizability of Nonlinear Stress–Strain Relationships for Reynolds Stress Closures

A class of recently developed explicit algebraic stress models based on tensorially quadratic stress--strain relations [7] is subjected to a systematical realizability analysis. It is found that these models, which are of particular interest for their rigorous derivation from linear second-moment closure models, tend to produce inappropriate unrealizable results like negative turbulence energy components, even in simple shear flows. The cause of the defect is identified in conjunction with a set of realizability-furnishing constraints on the model coefficients. With the exception of the silent normal stress component in accelerated flow, the nature and rationale of the explicit algebraic stress model suggested by Gatski and Speziale [7] can be extended to maintain the realizability principle. Results obtained from the corresponding quasi-realizable quadratic eddy-viscosity model are reported in comparison with other nonlinear modelling practices.     

Numerical simulation of thin airfoil stall by using a modified DES approach

The detached‐Eddy simulation (DES) method was applied to calculate pre‐ and post‐stall aerodynamic characteristics of airfoil stall. A discrepancy between numerical and experimental data was observed near the stall regime for the airfoil NACA64_A‐ 006 which is a thin airfoil stall type. The reason of this discrepancy and one possible way for improvement of the numerical model are discussed here. It is shown that the use of the Baldwin–Lomax model in the RANS region improves the DES results in this case. If the relevant factors (grid density, time step, turbulence model, etc.) are appropriately taken into account, the DES approach could reliably predict stall aerodynamical characteristics. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparative studies of numerical methods for evaluating aerodynamic characteristics of two-dimensional airfoil at low Reynolds numbers

This study investigates the predictability of the aerodynamic performance of some numerical methods at low Reynolds numbers and their dependency on the geometric shape of airfoil. We conducted three-dimensional large-eddy simulations (3-D LES), two-dimensional laminar simulations (2-D Lam), and Reynolds-averaged Navier–Stokes simulations with Baldwin–Lomax (2-D RANS(BL)) and Spalart–Allmaras (2-D RANS(SA)) turbulence models. Although there is little discrepancy between the 3-D LES, 2-D Lam, and 2-D RANS(SA) results in terms of the lift and drag characteristics, significant differences are observed in the predictability of the separation and reattachment points. The predicted lift, separation, and reattachment points of the 2-D Lam are qualitatively similar to those of the 3-D LES, except for high angles of attack at which a massive separation occurs. The 2-D RANS(SA) shows good predictability of the lift and separation points, but it does not estimate reattachment points accurately. The 2-D RANS(BL) fails to predict the precise separation points, which results in a poor lift predictability. These characteristics appear regardless of the airfoil geometry shapes. The results suggest that a 2-D Lam without any turbulence models can be used to estimate qualitative airfoil aerodynamic characteristics at the low Reynolds numbers.     

Detached eddy simulation of flow around two wall-mounted cubes in tandem

We investigate the performance of unsteady Reynolds-averaged Navier–Stokes (URANS) computation and various versions of detached eddy simulation (DES) in resolving coherent structures in turbulent flow around two cubes mounted in tandem on a flat plate at Reynolds number (Re) of 22,000 and for a thin incoming boundary layer. Calculations are carried out using four different coherent structure resolving turbulence models: (1) URANS with the Spalart–Allmaras model; (2) the standard DES [Spalart, P.R., Jou, W.H., Strelets, M., Allmaras, S.R., 1997. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In: Liu, C., Liu, Z., (Eds.), Advances in DNS/LES. Greyden Press, Columbus, OH]; (3) the Delayed DES (DDES); and (4) the DES with a low-Re modification (DES-LR) [Spalart, P., Deck, S., Shur, M., Squires, K., Strelets, M., Travin, A., 2006. A new version of detached eddy simulation, resistant to ambiguous grid densities. Theor. Comput. Fluid Dyn. 20 (3), 181–195]. The grid sensitivity of the computed solutions is examined by carrying out simulations on two successively refined grids. The computed results for all cases are compared with the experimental measurements of Martinuzzi and Havel [Martinuzzi, R., Havel, B., 2000. Turbulent flow around two interfering surface-mounted cubic obstacles in tandem arrangement. ASME J. Fluids Eng. 122, 24–31] for two different cube spacings. All turbulence models reproduce essentially identical separation of the approach thin boundary layer and yield an unsteady horseshoe vortex system consisting of multiple vortices in the leading edge region of the upstream cube. Significant discrepancies between the URANS and all DES solutions are observed, however, in other regions of interest such as the shear layers emanating from the cubes, the inter-cube gap and the downstream wake. Regardless of the grid refinement, URANS fails to capture key features of the mean flow, including the second horseshoe vortex in the upstream junction and recirculating flow on the top surface of the downstream cube for the large cube spacing, and underestimates significantly turbulence statistics in most regions of the flow for both cases. On the coarse mesh, all three DES approaches appear to yield very similar results and fail to reproduce the second horseshoe vortex. The standard DES and DDES solutions obtained on the fine meshes are essentially identical and both suffer from premature switching to unresolved DNS, due to the mis-interpretation of grid refinement as wall proximity, which leads to spurious vortices in the inter-cube region. Numerical solutions show that the low-Re modification (DES-LR) is critical prerequisite in DES on the ambiguously fine – not fine enough for full LES – mesh to prevent excessive nonlinear drop of the subgrid eddy viscosity in low cell-Re regions like in the inter-obstacle gap. Mean flow quantities and turbulence statistics obtained with DES-LR on the fine mesh are in good overall agreement with the measurements in most regions of interest for both cases.     

Accounting for Buoyancy Effects in the Explicit Algebraic Stress Model: Homogeneous Turbulent Shear Flows

Most explicit algebraic stress models are formulated for turbulent shear flows without accounting for external body force effects, such as the buoyant force. These models yield fairly good predictions of the turbulence field generated by mean shear. As for thermal turbulence generated by the buoyant force, the models fail to give satisfactory results. The reason is that the models do not explicitly account for buoyancy effects, which interact with the mean shear to enhance or suppress turbulent mixing. Since applicable, coupled differential equations have been developed describing these thermal turbulent fields, it is possible to develop corresponding explicit algebraic stress models using tensor representation theory. While the procedure to be followed has been employed previously, unique challenges arise in extending the procedure for developing the algebraic representations to turbulent buoyant flows. In this paper the development of an explicit algebraic stress model (EASM) is confined to the homogeneous buoyant shear flow case to illustrate the methodology needed to develop the proper polynomial representations. The derivation is based on the implicit formulation of the Reynolds stress anisotropy at buoyant equilibrium. A five-term representation is found to be necessary to account properly for the effect of the thermal field. Thus derived, external buoyancy effects are represented in the scalar coefficients of the basis tensors, and structural buoyancy effects are accounted for in additional terms in the stress anisotropy tensor. These terms will not vanish even in the absence of mean shear. The performance of the new EASM, together with a two-equation (2-Eq) model, the non-buoyant EASM of Gatski and Speziale (1993) and a full second-order model, is assessed against direct numerical simulations of homogeneous, buoyant shear flows at two different Richardson numbers representing weak and strong buoyancy effects. The calculations show that this five-term representation yields better results than the 2-Eq model and the EASM of Gatski and Speziale where buoyancy effects are not explicitly accounted for. Received 5 March 2001 and accepted 15 January 2002     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Unsteady RANS and detached-eddy simulations of flow around a circular cylinder in ground effect

Unsteady Reynolds-averaged Navier–Stokes (URANS) simulations and detached-eddy simulations (DES) were performed of flow around a circular cylinder placed near and parallel to a moving ground, on which substantially no boundary layer developed to interfere with the cylinder. The results were compared with experiments previously reported by the authors to examine how accurately the URANS and DES can predict the cessation of von Kármán-type vortex shedding and the attendant critical drag reduction of the cylinder in ground effect. The DES, which were performed in a three-dimensional domain with spanwise periodicity imposed, correctly captured the cessation of the vortex shedding, whereas both two- and three-dimensional URANS also predicted it but at a much smaller gap-to-diameter ratio compared with the experiments. The wake structures of the cylinder predicted by the DES were in good agreement with the experiments in both large- and small-gap regimes, and also in the intermediate-gap regime, where the DES captured the intermittence of the vortex shedding in the near-wake region. Based on the results obtained, further discussions are also given to the reason why the von Kármán-type vortices in the URANS solutions incorrectly ‘survived’ until the cylinder came much closer to the ground.     

Fully explicit and self-consistent algebraic Reynolds stress model

A fully explicit, self-consistent algebraic expression (for Reynolds stress) which is the exact solution to the Reynolds stress transport equation in the weak-equilibrium limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the weak-equilibrium assumption treat the production-to-dissipation (P/) ratio as a constant, resulting in an effective viscosity that can be singular away from the equilibrium limit. In this paper the set of simultaneous algebraic Reynolds stress equations in the weak-equilibrium limit are solved in the full nonlinear form and the eddy viscosity is found to be nonsingular. Preliminary tests indicate that the model performs adequately, even for three-dimensional mean-flow cases. Due to the explicit and nonsingular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-19480.     

An efficient reduced-order modelling approach based on fluid eigenmodes and boundary element method

This paper presents an efficient reduced-order modelling approach based on the boundary element method. In this approach, the eigenvalue problem of the unsteady flows is defined based on the unknown wake singularities. By constructing this reduced-order model, the body quasi-static eigenmodes are removed from the eigensystem and it is possible to obtain satisfactory results without using the static correction technique when enough eigenmodes are used. In addition to the conventional method, eigenanalysis and reduced-order modelling of unsteady flows over a NACA 0012 airfoil, a wing with NACA 0012 section and a wing–body combination are performed using the proposed reduced order modelling (ROM) method. Numerical examples are presented that demonstrate the accuracy and computational efficiency of the present method.     

Implicit large eddy simulation using the high‐order correction procedure via reconstruction scheme

We investigate implicit large eddy simulation of the Taylor–Green vortex, Comte‐Bellot–Corrsin experiment, turbulent channel flow and transitional and turbulent flow over an SD7003 airfoil using the high‐order unstructured correction procedure via reconstruction (CPR) scheme, also known as the flux reconstruction scheme. We employ P1 (second‐order) to P5 (sixth‐order) spatial discretizations. Results show that the CPR scheme can accurately predict turbulent flows without the addition of a sub‐grid scale model. Numerical dissipation, concentrated at the smallest resolved scales, is found to filter high‐frequency content from the solution. In addition, the high‐order schemes are found to be more accurate than the low‐order schemes on a per degree of freedom basis for the canonical test cases we consider. These results motivate the further investigation and use of the CPR scheme for simulating turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.     

A novel method for automated grid generation of ice shapes for local‐flow analysis

Modelling a complex geometry, such as ice roughness, plays a key role for the computational flow analysis over rough surfaces. This paper presents two enhancement ideas in modelling roughness geometry for local flow analysis over an aerodynamic surface. The first enhancement is use of the leading‐edge region of an airfoil as a perturbation to the parabola surface. The reasons for using a parabola as the base geometry are: it resembles the airfoil leading edge in the vicinity of its apex and it allows the use of a lower apparent Reynolds number. The second enhancement makes use of the Fourier analysis for modelling complex ice roughness on the leading edge of airfoils. This method of modelling provides an analytical expression, which describes the roughness geometry and the corresponding derivatives. The factors affecting the performance of the Fourier analysis were also investigated. It was shown that the number of sine–cosine terms and the number of control points are of importance. Finally, these enhancements are incorporated into an automated grid generation method over the airfoil ice accretion surface. The validations for both enhancements demonstrate that they can improve the current capability of grid generation and computational flow field analysis around airfoils with ice roughness. Copyright © 2004 John Wiley & Sons, Ltd.     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part II Buoyant flow

This paper presents a derivation of an explicit algebraic model for two-dimensional (2-D) buoyant flows. It is an extension of the work reported in Part I (So et al. [27]). The tensor representation method of Jongen and Gatski [14] is extended to derive an explicit algebraic Reynolds stress model (EASM) for 2-D buoyant flow invoking the Boussinesq approximation. The projection methodology is further extended to treat the heat flux transport equation in the derivation of an explicit algebraic heat flux model (EAHFM) for buoyant flow. Again, the weak equilibrium assumption is invoked for the scaled Reynolds stress and scaled heat flux equation. An explicit algebraic model for buoyant flows is then formed with the EASM and EAHFM. From the derived EAHFM, an expression for the thermal diffusivity tensor in buoyant shear flows is deduced. Furthermore, a turbulent Prandtl number (Pr_T) for each of the three heat flux directions is determined. These directional Pr_T are found to be a function of the gradient Richardson number. Alternatively, a scalar Pr_T can be derived; its value is compared with the directional Pr_T. The EASM and EAHFM are used to calculate 2-D homogeneous buoyant shear flows and the results are compared with direct numerical simulation data and other model predictions, where good agreement is obtained. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University