Nonlinear Versus Linear Stress-Strain Relations in Engine Turbulence Modeling Under Swirl and Squish Flow Conditions

Author(s):  
Mirko Baratta ◽  
Andrea E. Catania ◽  
Stefano d’Ambrosio

A general form of the stress-strain constitutive relation was introduced for the application of two nonlinear k-ε turbulence models, namely, the algebraic Reynolds stress model of Gatski and Speziale (1993, “On Explicit Algebraic Stress Models for Complex Turbulent Flows,” J. Fluid Mech., 254, pp. 59–78) and the cubic model of Lien et al. (1996, “Low Reynolds Number Eddy-Viscosity Modeling Based on Non-Linear Stress-Strain/Vorticity Relations,” Proceedings of Third Symposium on Engineering Turbulence Modeling and Measurements, Crete, Greece), to the numerical analysis of flow fields in a test engine with flat-piston and bowl-in-piston arrangements, under swirling and no-swirling flow motored conditions. The model capabilities in capturing turbulent flow features were compared to those of the upgraded linear RNG k-ε model, which was previously indicated as a good compromise between accuracy and computational cost. Evaluations were made on the basis of the predicted flow evolution throughout the whole engine cycle, as well as of the comparison between the numerical and experimental results. Furthermore, the effect of the stress-strain relationship on the predicted averaged turbulence quantities and anisotropy-invariant values were examined, in addition to the sensitivity of the nonlinear models to the mesh quality. Finally, prospects concerning possible improvements of turbulence eddy-viscosity models were presented. The predictions were made by a newly developed CFD code embedding various accuracy-order finite-volume discretization schemes. Modified wall boundary conditions with respect to the conventional logarithmic-function approach were used, so as to make the local equilibrium hypothesis virtually ineffective.

Author(s):  
Mirko Baratta ◽  
Andrea E. Catania ◽  
Stefano d’Ambrosio

A general form of the stress-strain constitutive relation was introduced for the application of two nonlinear k-ε turbulence models, namely, the ARS model of Gatski and Speziale ([1]) and the Cubic model of Lien et al. ([2]), to the numerical analysis of flow fields in a test engine with flat-piston and bowl-in-piston arrangements, under swirling or no-swirling flow motored conditions. The model capabilities in capturing turbulent flow features were compared to those of the upgraded linear RNG k-ε model which was previously indicated as a good compromise between accuracy and computational cost ([3]). Evaluations were made on the basis of the predicted flow evolution throughout the whole engine cycle, as well as of the comparison between numerical and experimental results. Furthermore, the effect of the stress-strain relationship on the predicted averaged turbulence quantities and anisotropy invariant values were examined, in addition to the sensitivity of the nonlinear models to the mesh quality. Finally, prospects concerning possible improvements of turbulence Eddy Viscosity Models (EVM) were presented. The predictions were made by a newly developed CFD code embedding various accuracy-order finite-volume discretization schemes. Modified wall boundary conditions with respect to the conventional logarithmic-function approach were used, so as to give negligible importance to the local equilibrium hypothesis.


2019 ◽  
Vol 130 ◽  
pp. 01013
Author(s):  
Hariyo Priambudi Setyo Pratomo ◽  
Fandi Dwiputra Suprianto ◽  
Teng Sutrisno

Turbulence simulation remains one of the active research activities in computational engineering. Along with the increase in computing power and the prime motivation of improving the accuracy of statistical turbulence modeling approaches and reducing the expensive computational cost of both direct numerical and large turbulence scale- resolving simulations, various hybrid turbulence models being capable of capturing unsteadiness in the turbulence are now accessible. Nevertheless this introduces the daunting task to select an appropriate method for different cases as one can not know a priori the inherent nature of the turbulence. It is the aim of this paper to address recent progresses and further researches within a branch of the hybrid RANS-LES models examined by the first author as simple test cases but generating complex turbulent flows are available from experimentation. In particular, failure of a seamless hybrid formulation not explicitly dependent on the grid scale is discussed. From the literature, it is practical that at least one can go on with confidence when choosing a potential hybrid model by intuitively distinguishing between strongly and weakly unstable turbulent flows.


1995 ◽  
Vol 48 (4) ◽  
pp. 189-212 ◽  
Author(s):  
G. J. Brereton ◽  
R. R. Mankbadi

Turbulent flow which undergoes organized temporal unsteadiness is a subject of great importance to unsteady aerodynamic and thermodynamic devices. Of the many classes of unsteady flows, those bounded by rigid smooth walls are particularly amenable to fundamental studies of unsteady turbulence and its modeling. These flows are presently being given increased attention as interest grows in the prospect of predicting non-equilibrium turbulence and because of their relevance to turbulence–acoustics interactions, in addition to their importance as unsteady flows in their own right. It is therefore timely to present a review of recent advances in this area, with particular emphasis placed on physical understanding of the turbulent processes in these flows and the development of turbulence models to predict them. A number of earlier reviews have been published on unsteady turbulent flows, which have tended to focus on specific aspects of certain flows. This review is intended to draw together, from the diverse literature on the subject, information on fundamental aspects of these flows which are relevant to improved understanding and development of predictive models. Of particular relevance are issues of instability and transition to turbulence in reciprocating flows, the robustness of coherent structures in wall-bounded flows to forced perturbations (in contrast to the relative ease of manipulation in free shear flows), unsteady scalar transport, improved measurement technology, recent contributions to target data for model testing and the quasi-steady and non-steady rapid distortion approaches to turbulence modeling in these flows. The present article aims to summarize recent contributions to this research area, with a view to consolidating comprehension of the well-known basics of these flows, and drawing attention to critical gaps in information which restrict our understanding of unsteady turbulent flows.


2009 ◽  
Vol 131 (4) ◽  
Author(s):  
Pavel E. Smirnov ◽  
Florian R. Menter

A rotation-curvature correction suggested earlier by Spalart and Shur (1997, “On the Sensitization of Turbulence Models to Rotation and Curvature,” Aerosp. Sci. Technol., 1(5), pp. 297–302) for the one-equation Spalart–Allmaras turbulence model is adapted to the shear stress transport model. This new version of the model (SST-CC) has been extensively tested on a wide range of both wall-bounded and free shear turbulent flows with system rotation and/or streamline curvature. Predictions of the SST-CC model are compared with available experimental and direct numerical simulations (DNS) data, on the one hand, and with the corresponding results of the original SST model and advanced Reynolds stress transport model (RSM), on the other hand. It is found that in terms of accuracy the proposed model significantly improves the original SST model and is quite competitive with the RSM, whereas its computational cost is significantly less than that of the RSM.


Author(s):  
Dominik von Lavante ◽  
Eckart Laurien

With recent progress in high-temperature pebble-bed reactor programs research focus has started to include more ancillary engineering issues. One very important aspect for the realisability is the mixing of hot and colder helium in the reactor lower plenum. Under nominal operating conditions, depending on core design, the temperature of hot gas leaving the core can locally differ up to 210° C. Due to material limitations, these temperature differences have to be reduced to at least ±15° C. Several reduced-size air experiments have been performed on this problem, but their applicability to modern commercially sized reactors is not certain. With the rise in computing power CFD simulations can be performed in addition, but advanced turbulence modeling is necessary due to the highly swirling and turbulent nature of this flow. The presented work uses the geometry of the German HTR-Modul which consists of an annular mixing channel and radially arranged ribs. Using the commercial CFD code ANSYS CFX, we have made detailed analyses of the complex 3D vortical flow phenomena within this geometry. Several momentum transport turbulence models, e.g. the classical k-e model, advanced two-equation models and Reynolds-Stress Models were compared with respect to their accuracy for this particular flow. In addition, the full set of turbulent scalar flux transport equations was implemented for modeling the three components of turbulent transport of enthalpy seperately and were compared with the standard turbulent Prandtl number approach. As expected from previous work in related fields of turbulence modeling, the differences in predicting the mixing performance between models were significant. Only the full Reynolds-Stress model coupled with the scalar flux equations was able to reproduce the experimentally observed reduction of mixing efficiency with increasing Reynolds number. The correct scaling of mixing efficiencies demonstrates that the utilized turbulence models are able to reproduce the physics of the underlying flow. Hence they could be employed for the scaling and optimization of the lower plenum geometry. The results also showed that the original geometry used for the HTR-Modul is insufficient to provide adequate mixing, and that hence a not sufficiently mixed coolant for future reactor designs might be an issue. Based on this work, an optimization for future lower plenum geometries has become feasible.


2008 ◽  
Vol 130 (9) ◽  
Author(s):  
G. Dufour ◽  
J.-B. Cazalbou ◽  
X. Carbonneau ◽  
P. Chassaing

Rotation and curvature (RC) effects on turbulence are expected to impact losses and flow structure in turbomachines. This paper examines two recent eddy-viscosity-model corrections devised to account for these effects: the Spalart and Shur (1997, “On the Sensitization of Turbulence Models to Rotation and Curvature,” Aerosp. Sci. Technol., 1(5), pp. 297–302) correction to the model of Spalart and Allmaras (1994, “A One-Equation Turbulence Model for Aerodynamic Flows,” Rech. Aerosp., 1, pp. 5–21) and the correction of Cazalbou et al. (2005, “Two-Equation Modeling of Turbulent Rotating Flows,” Phys. Fluids., 17, p. 055110) to the (k,ϵ) model. The method of verification and validation is applied to assess the impact of these corrections on the computation of a centrifugal-compressor test case. First, a review of RC effects on turbulence as they apply to centrifugal compressors is made. The two corrected models are then presented. Second, the Radiver open test case (Ziegler K. U., Gallus, H. E., and Niehuis R., 2003, “A Study on Impeller Diffuser Interaction Part 1: Influence on the Performance,” ASME J. Turbomach, 125, pp. 173–182) is used as a basis for the assessment of the two corrections. After a physical-consistency analysis, the Richardson extrapolation is applied to quantify the numerical errors involved in all the calculations. Finally, experimental data are used to perform validation for both global and local predictions. The consistency analysis shows that both corrections lead to significant changes in the turbulent field, in perfect agreement with the underlying theoretical considerations. The uncertainty analysis shows that the predictions of the global performances are more sensitive to grid refinement than they are to RC turbulence modeling. However, the opposite conclusion is drawn with regard to the prediction of some local flow properties: Improvements are obtained with the RC corrections, the best results being observed for the RC-corrected (k,ϵ) model.


Author(s):  
O. Z. Mehdizadeh ◽  
L. Temmerman ◽  
B. Tartinville ◽  
Ch. Hirsch

Turbulence modeling remains an active CFD development front for turbomachinery as well as for general industrial applications. While DNS and even LES still seem out of reach within the typical industrial design cycle due to their high computational cost, RANS-based models remain the workhorse of CFD. Currently, the most widely used models are Linear Eddy-Viscosity Models (LEVM), despite their known limitations for certain flow complexities. Therefore, extending the reliability of eddy-viscosity models to more complex flows without significantly increasing the computational cost can immediately contribute to more reliable CFD results for wider range of applications. This, in turn, can further reduce the need for costly tests and consequently can reduce the product development cost. A promising approach to achieve this goal is using Explicit Algebraic Reynolds Stress Models (EARSM), obtained through a simplification of the full Differential Reynolds Stress Models (DRSM), and can be perceived as an extension of LEVMs by including the non-linear relation between the turbulence stress tensor, the mean-flow gradient and the turbulence scales. These models are thus less demanding than DRSM, yet capable of capturing more complex turbulence features, compared to LEVM, such as anisotropy in the normal stresses. This may be particularly important in corner flows, for instance, in the hub-blade regions or in diffusers. This work explores the application of EARSM models to a double diffuser and a high-performance centrifugal compressor stage (HPCC). The results are compared to available experimental data [1,2] showing the importance of including the anisotropy of turbulence in the model, particularly in presence of turbulent corner flows in a diffuser. Furthermore, the EARSM results are also compared to results from the commonly used SST turbulence model. The CFD comparison includes details of the flow structure in the diffuser, where the most noticeable impact from the use of EARSM turbulence models is expected.


2005 ◽  
Vol 73 (3) ◽  
pp. 405-412 ◽  
Author(s):  
Hermann F. Fasel ◽  
Dominic A. von Terzi ◽  
Richard D. Sandberg

A flow simulation Methodology (FSM) is presented for computing the time-dependent behavior of complex compressible turbulent flows. The development of FSM was initiated in close collaboration with C. Speziale (then at Boston University). The objective of FSM is to provide the proper amount of turbulence modeling for the unresolved scales while directly computing the largest scales. The strategy is implemented by using state-of-the-art turbulence models (as developed for Reynolds averaged Navier-Stokes (RANS)) and scaling of the model terms with a “contribution function.” The contribution function is dependent on the local and instantaneous “physical” resolution in the computation. This physical resolution is determined during the actual simulation by comparing the size of the smallest relevant scales to the local grid size used in the computation. The contribution function is designed such that it provides no modeling if the computation is locally well resolved so that it approaches direct numerical simulations (DNS) in the fine-grid limit and such that it provides modeling of all scales in the coarse-grid limit and thus approaches a RANS calculation. In between these resolution limits, the contribution function adjusts the necessary modeling for the unresolved scales while the larger (resolved) scales are computed as in large eddy simulation (LES). However, FSM is distinctly different from LES in that it allows for a consistent transition between RANS, LES, and DNS within the same simulation depending on the local flow behavior and “physical” resolution. As a consequence, FSM should require considerably fewer grid points for a given calculation than would be necessary for a LES. This conjecture is substantiated by employing FSM to calculate the flow over a backward-facing step and a plane wake behind a bluff body, both at low Mach number, and supersonic axisymmetric wakes. These examples were chosen such that they expose, on the one hand, the inherent difficulties of simulating (physically) complex flows, and, on the other hand, demonstrate the potential of the FSM approach for simulations of turbulent compressible flows for complex geometries.


Author(s):  
Lei-Yong Jiang ◽  
Ian Campbell

The flow field of a propane-air diffusion flame combustor with interior and exterior conjugate heat transfers was numerically investigated. Solutions obtained from four turbulence models together with a laminar flamelet combustion model, discrete ordinates radiation model and enhanced wall treatment are presented and discussed. The numerical results are compared, in detail, with a comprehensive database obtained from a series of experimental measurements. It is found that the Reynolds stress model (RSM), a second moment closure, illustrates superior performance over three popular two-equation eddy-viscosity models. Although the main flow features are captured by all four turbulence models, only the RSM is able to successfully predict the lengths of both recirculation zones and the turbulence kinetic energy distribution in the combustor chamber. In addition, it provides fairly good predictions for all Reynolds stress components, except for the circumferential normal stress at downstream sections. However, the superiority of the RSM is not so obvious for the temperature and species predictions in comparison with eddy-viscosity models, except for the standard k-ε model. This suggests that coupling between the RSM and combustion models needs to be further improved in order to enhance its applications in practical combustion systems.


Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 148 ◽  
Author(s):  
Chunhui Zhang ◽  
Charles Patrick Bounds ◽  
Lee Foster ◽  
Mesbah Uddin

In today’s road vehicle design processes, Computational Fluid Dynamics (CFD) has emerged as one of the major investigative tools for aerodynamics analyses. The age-old CFD methodology based on the Reynolds Averaged Navier–Stokes (RANS) approach is still considered as the most popular turbulence modeling approach in automotive industries due to its acceptable accuracy and affordable computational cost for predicting flows involving complex geometries. This popular use of RANS still persists in spite of the well-known fact that, for automotive flows, RANS turbulence models often fail to characterize the associated flow-field properly. It is even true that more often, the RANS approach fails to predict correct integral aerodynamic quantities like lift, drag, or moment coefficients, and as such, they are used to assess the relative magnitude and direction of a trend. Moreover, even for such purposes, notable disagreements generally exist between results predicted by different RANS models. Thanks to fast advances in computer technology, increasing popularity has been seen in the use of the hybrid Detached Eddy Simulation (DES), which blends the RANS approach with Large Eddy Simulation (LES). The DES methodology demonstrated a high potential of being more accurate and informative than the RANS approaches. Whilst evaluations of RANS and DES models on various applications are abundant in the literature, such evaluations on full-car models are relatively fewer. In this study, four RANS models that are widely used in engineering applications, i.e., the realizable k - ε two-layer, Abe–Kondoh–Nagano (AKN) k - ε low-Reynolds, SST k - ω , and V2F are evaluated on a full-scale passenger vehicle with two different front-end configurations. In addition, both cases are run with two DES models to assess the differences between the flow predictions obtained using RANS and DES.


Sign in / Sign up

Export Citation Format

Share Document