A Methodology for Simulating Compressible Turbulent Flows

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):  
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 modelling 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 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 modelling if the computation is locally well resolved so that it approaches a DNS in the fine-grid limit and such that it provides modelling of all scales in the coarsegrid limit and thus approaches an unsteady RANS calculation. In between these resolution limits, the contribution function adjusts the necessary modelling for the unresolved scales while the larger (resolved) scales are computed as in traditional LES. However, FSM is distinctly different from LES in that it allows for a consistent transition between (unsteady) 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 traditional LES. This conjecture is substantiated by employing FSM to calculate the flow over a backward-facing step at low Mach number and a supersonic, axisymmetric baseflow. 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 a wide range of compressible flows.


2002 ◽  
Vol 124 (4) ◽  
pp. 933-942 ◽  
Author(s):  
H. F. Fasel ◽  
J. Seidel ◽  
S. Wernz

A new flow simulation methodology (FSM) for computing turbulent shear flows is presented. The development of FSM was initiated in close collaboration with C. Speziale (then at Boston University). The centerpiece of FSM is a strategy to provide the proper amount of modeling of the subgrid scales. The strategy is implemented by use of a “contribution function” which is dependent on the local and instantaneous “physical” resolution in the computation. This physical resolution is obtained 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 the computation approaches a direct numerical simulation in the fine grid limit, or provides modeling of all scales in the coarse grid limit and thus approaches an unsteady 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 traditional large-eddy simulations (LES). However, a LES that is based on the present strategy is distinctly different from traditional LES in that the required amount of modeling is determined by physical considerations, and that state-of-the-art turbulence models (as developed for Reynolds-averaged Navier-Stokes) can be employed for modeling of the unresolved scales. Thus, in contrast to traditional LES based on the Smagorinsky model, with FSM a consistent approach (in the local sense) to the coarse grid and fine grid limits is possible. As a consequence of this, FSM should require much fewer grid points for a given calculation than traditional LES or, for a given grid size, should allow computations for larger Reynolds numbers. In the present paper, the fundamental aspects of FSM are presented and discussed. Several examples are provided. The examples were chosen such that they expose, on the one hand, the inherent difficulties of simulating complex wall bounded flows, and on the other hand demonstrate the potential of the FSM approach.


Author(s):  
Juan B. V. Wanderley ◽  
Gisele H. B. Souza ◽  
Carlos Levi

Author’s previous work Wanderley [1] presented an efficient numerical method to investigate VIV phenomenon on circular cylinders. The numerical model solves the unsteady Reynolds Average Navier–Stokes equations for slightly compressible flows using the Beam–Warming implicit factored scheme. In the present work, the effect of the turbulence model on the results is evaluated for both Baldwin Lomax and k-ε models. To demonstrate the quality of the numerical method, results for the transversal oscillation of a cylinder laterally supported by spring and damper are compared with experimental data. The application of the turbulence models showed the much better agreement of the k-ε model with the experimental results.


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.


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.


Author(s):  
M. S. Anand ◽  
R. Eggels ◽  
M. Staufer ◽  
M. Zedda ◽  
J. Zhu

A general-purpose combustion Computational Fluid Dynamics (CFD) design analysis tool has been developed. The method is pressure-based and applicable to both incompressible and compressible flows. The unstructured finite-volume approach used can take arbitrary shapes of mesh cells to resolve complicated geometries. Turbulence is simulated either by Reynolds-Averaged Navier-Stokes (RANS) or by Large Eddy Simulation (LES) approaches. Combustion is modeled by various combinations of combustion chemistry and combustion-turbulence models including transport probability density function (PDF) model. A Lagrangian approach is used to simulate fuel spray droplet. The resulting tool has being used in routine combustor simulations for a variety of commercial and military combustor development programs. Application examples presented include simulations of several combustors and comparisons with available rig data.


1990 ◽  
Vol 43 (5S) ◽  
pp. S240-S245 ◽  
Author(s):  
N. Aubry

The proper orthogonal decomposition (POD), also called Karhunen-Loe`ve expansion, which extracts ‘coherent structures’ from experimental data, is a very efficient tool for analyzing and modeling turbulent flows. It has been shown that it converges faster than any other expansion in terms of kinetic energy (Lumley 1970). First, the POD is applied to the chaotic solution of the Lorenz equations. The dynamics of the Lorenz attractor is reconstructed by only the first three POD modes. In the second part of this paper, we show how the POD can be used in turbulence modeling. The particular case studied is the wall region of a turbulent boundary layer. In this flow, the velocity field is expanded into POD modes in the normal direction and Fourier modes in the streamwise and spanwise directions. Dynamical systems are obtained by Galerkin projections of the Navier Stokes equations onto the different modes. Aubry et al. (1988) applied the technique to derive and study a ten dimensional representation which reproduced qualitatively the bursting event experimentally observed. It is shown that streamwise modes, absent in Aubry et al.’s model, participate to the bursting events. This agrees remarkably well with experimental observations. In both examples, the dynamics of the original system is very well recovered from the contribution of only a few modes.


Author(s):  
Omid Abouali ◽  
Mohammad M. Alishahi ◽  
Homayoon Emdad ◽  
Goodarz Ahmadi

A 3-D Thin Layer Navier-Stokes (TLNS) code for solving viscous supersonic flows is developed. The new code uses several numerical algorithms for space and time discretization together with appropriate turbulence modeling. Roe’s method is used for discretizing the convective terms and the central differencing scheme is employed for the viscous terms. An explicit time marching technique and a finite volume space discretization are used. The developed computational model can handle both laminar and turbulent flows. The Baldwin-Lomax model and Degani-Schiff modifications are used for turbulence modeling. The computational model is applied to a hypersonic laminar flow at Mach 7.95 around a cone at different incidence angles. The circumferential pressure distribution is compared with the experimental data. The cross-sectional Mach number contours are also presented. It is shown that in addition to the outer shock, a cross-flow shock wave is also present in the flow field. The cases of supersonic turbulent flows with Mach number 3 around a tangent-ogive with incidence angles of 6° and a secant-ogive with incidence angles of 10° are also studied. The circumferential pressure distributions are compared with the experimental data and the Euler code results and good agreement is obtained. The cross-sectional Mach number contours are also presented. It is shown that in this case also in addition to the outer shock, a cross-flow shock wave is also present at the incidence angle of 10°.


2005 ◽  
Author(s):  
Sowjanya Vijiapurapu ◽  
Jie Cui

The Reynolds averaged Navier-Stokes (RANS) equations were solved along with three turbulence models, namely κ-ε, κ-ω, and Reynolds stress models (RSM), to study the fully developed turbulent flows in circular pipes roughened by repeated square ribs. The spacing between the ribs was varied to form three representative types of surface roughness; d–type, intermediate, and k–type. Solutions of these flows at two Reynolds numbers were obtained using the commercial computational fluid dynamics (CFD) software Fluent. The numerical results were validated against experimental measurements and other numerical data published in literature. Extensive investigation of effects of rib spacing and Reynolds number on the pressure and friction resistance, flow and turbulence distribution was presented. The performance of three turbulence models was also compared and discussed.


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.


Sign in / Sign up

Export Citation Format

Share Document