A Methodology for Simulations of Complex Turbulent 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.

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.


2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.


2006 ◽  
Author(s):  
Francisco Elizalde-Blancas ◽  
Ismail Celik ◽  
Suryanarayana Pakalapati

In this study numerical solutions are presented for a steady state, incompressible, 2-D turbulent flow near a wall. For this specific problem a manufactured (exact) solution was provided by the organizers of the 2006 Lisbon Workshop [6]. With the help of manufactured solution, assessment of the true error and other relevant uncertainty measures are possible. The calculations were performed using the commercial flow solver FLUENT along with some user defined functions to define source terms and velocity profiles at boundaries. Though the flow regime is turbulent; the numerical solution is carried out for pseudo-laminar flow. This was done in order to avoid the errors implicit in turbulence models. The transformation from turbulent to laminar flow was done by defining a momentum source term which precludes the pressure gradient term. A detailed grid convergence analysis was performed. Using three-grid triplets the limiting values of the variables solved as the grid size tends to zero were calculated using different extrapolations. The L2 norms of the true error obtained from various extrapolations are assessed. These results exhibit solution convergence as the grid size decreases. It was also shown that cubic spline extrapolation perform the best among the methods considered.


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.


1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


Author(s):  
Lokesh Kalyan Gutti ◽  
◽  
Bhupendra Singh Chauhan ◽  
Hee-Chang Lim ◽  
◽  
...  

For incompressible flow simulation, it is commonly accepted to use uniform meshes to solve the governing equation of turbulent boundary layer. It follows the laws of conservation stabilizing the flow field in the domain and preventing odd-even decoupling in the pressure field. In this study, Large Eddy Simulation (LES) has been conducted in a long channel. In order to calculate the turbulent boundary layer in the channel, the unsteady Navier-Stokes equations has been adopted at a Reynolds number =180, which is based on mean centerline velocity and the half-width of the channel. The mesh used in this study was based on both stretch and uniform mesh having grid points, which is corresponding to . Turbulence statistics were also calculated to compare to the existing results. In the results, the turbu lent boundary layer was fully developed at around . In addition, fully developed channel flow was achieved at the non-dimensional time of .


2017 ◽  
Vol 139 (4) ◽  
Author(s):  
Renfang Huang ◽  
Xianwu Luo ◽  
Bin Ji ◽  
Qingfeng Ji

A modified partially averaged Navier–Stokes model (MPANS) is proposed by treating the standard k–ε model as the parent model and formulating the unresolved-to-total kinetic energy ratio fk as a function of the local grid size and turbulence length scale. Flows over a backward facing step are used to evaluate the performance of MPANS mode. Computations of the standard k–ε model, the constant fk partially averaged Navier–Stokes (PANS) models (fk = 0.6, 0.7), and the two-stage PANS model are carried out for comparisons. Based on the detailed analyses of calculated results and experimental data, the MPANS model performs better to predict the reattachment length together with the corner vortex and provides overall improved statistics of skin frictions, pressures, velocity profiles, and Reynolds stresses, demonstrating its promising applications in industrial turbomachines that often encounter with flow separations.


1998 ◽  
Author(s):  
Thomas Hildebrandt ◽  
Leonhard Fottner

A thorough numerical study was conducted to simulate the flow field inside a highly loaded linear turbine cascade. The numerical investigation was focused on the secondary flow field as well as on the prediction of the overall design goals within reasonable accuracy limits. The influence of grid resolution was investigated in order to obtain detailed information about the requirements of a grid independent solution. Three different two-equation turbulence models were applied to two numerical grids of different resolution. Emphasis was laid on separating the influences of grid resolution and turbulence models. The Mach- and Reynolds numbers as well as the level of freestream turbulence were set to values typical of turbomachinery conditions. The computational study was carried out using a 3D state-of-the-art blockstructured Navier-Stokes solver. The comparison of the numerical results with experiments clearly revealed the different degree of agreement between simulation and measurement. This paper describes the application of a modern flow solver to a testcase which is relevant for practical turbomachinery design purposes. The agreement between the experiments and the results of the numerical study is good and in most cases well within the accuracy limits proposed by Strazisar and Denton (1995). It was found out that the main effect on the quality of the computations is the resolution of the numerical grid. The finest grid used reached over one million points halfspan, showing clearly superior results compared with a coarser, though still fine grid. The influence of different turbulence models on the numerically obtained flow field was relatively small in comparison with the the grid influence.


Author(s):  
Armin Zirkel ◽  
Guido Doebbener ◽  
Eckart Laurien

This paper presents the current state of an ongoing analysis and validation of turbulence models for three dimensional numerical simulations (CFD-simulations) within containments of nuclear reactors. A severe accident flow inside a containment could be caused by a leak in the primary circuit of the reactor. It is characterized by different velocities, mass transport, the anisotropy and time dependency of the turbulence field as well as the transition between laminar and turbulent flows. Another issue is the complex geometry of a containment with different rooms and obstacles. The analysis is using the experimental data of the THAI model containment, using the TH18 experiment as the reference case. In order to investigate the numerical error, different advection schemes were used and a grid-dependency study was carried out within a half model of the geometry. On a sufficiently fine grid transient simulations were performed by using the Shear Stress Transport and the Reynolds Stress turbulence models. The results of the simulations are showing different deviations from the experiment. Along with the results, a guide for future work is discussed.


Sign in / Sign up

Export Citation Format

Share Document