Multigrid Computation of Incompressible Flows Using Two-Equation Turbulence Models: Part II—Applications

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):  
Dimitrios A. Inglezakis ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

CFD (Computational Fluid Dynamics) solvers have become nowadays an integral part of the aerospace manufacturing process and product design, as their implementation allows for the prediction of the aerodynamic behavior of an aircraft in a relatively short period of time. Such an in-house academic solver, named Galatea, is used in this study for the prediction of the flow over the ARA (Aircraft Research Association) M151/1 aircraft model. The proposed node-centered finite-volume solver employs the RANS (Reynolds-Averaged Navier-Stokes) equations, combined with appropriate turbulence models, to account for the simulation of compressible turbulent flows on three-dimensional hybrid unstructured grids, composed of tetrahedral, prisms, and pyramids. A brief description of Galatea’s methodology is included, while attention is mainly directed toward the accurate prediction of pressure distribution on the wings’ surfaces of the aforementioned airplane, an uncommon combat aircraft research model with forward swept wings and canards. In particular, two different configurations of M151/1 were examined, namely, with parallel and expanding fuselage, while the obtained results were compared with those extracted with the commercial CFD software ANSYS CFX. A very good agreement is reported, demonstrating the proposed solver’s potential to predict accurately such demanding flows over complex geometries.


2003 ◽  
Vol 478 ◽  
pp. 227-235 ◽  
Author(s):  
J. D. GIBBON ◽  
Charles R. DOERING

Dissipation-range intermittency was first observed by Batchelor & Townsend (1949) in high Reynolds number turbulent flows. It typically manifests itself in spatio-temporal binary behaviour which is characterized by long, quiescent periods in the signal which are interrupted by short, active ‘events’ during which there are large excursions away from the average. It is shown that Leray's weak solutions of the three-dimensional incompressible Navier–Stokes equations can have this binary character in time. An estimate is given for the widths of the short, active time intervals, which decreases with the Reynolds number. In these ‘bad’ intervals singularities are still possible. However, the average width of a ‘good’ interval, where no singularities are possible, increases with the Reynolds number relative to the average width of a bad interval.


2009 ◽  
Vol 62 (4) ◽  
Author(s):  
Giancarlo Alfonsi

The approach of Reynolds-averaged Navier–Stokes equations (RANS) for the modeling of turbulent flows is reviewed. The subject is mainly considered in the limit of incompressible flows with constant properties. After the introduction of the concept of Reynolds decomposition and averaging, different classes of RANS turbulence models are presented, and, in particular, zero-equation models, one-equation models (besides a half-equation model), two-equation models (with reference to the tensor representation used for a model, both linear and nonlinear models are considered), stress-equation models (with reference to the pressure-strain correlation, both linear and nonlinear models are considered) and algebraic-stress models. For each of the abovementioned class of models, the most widely-used modeling techniques and closures are reported. The unsteady RANS approach is also discussed and a section is devoted to hybrid RANS/large methods.


1987 ◽  
Vol 109 (4) ◽  
pp. 345-352 ◽  
Author(s):  
M. Reggio ◽  
R. Camarero

A numerical procedure to solve three-dimensional incompressible flows in arbitrary shapes is presented. The conservative form of the primitive-variable formulation of the time-dependent Navier-Stokes equations written for a general curvilinear coordiante system is adopted. The numerical scheme is based on an overlapping grid combined with opposed differencing for mass and pressure gradients. The pressure and the velocity components are stored at the same location: the center of the computational cell which is used for both mass and the momentum balance. The resulting scheme is stable and no oscillations in the velocity or pressure fields are detected. The method is applied to test cases of ducting and the results are compared with experimental and numerical data.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


Author(s):  
N. Parameswara Rao ◽  
K. Arul Prakash

Numerical simulation of complex three-dimensional flow through the spiral casing has been studied using a finite element method. An explicit Eulerian velocity correction scheme has been employed to solve the Reynolds averaged Navier-stokes equations. The simulation has been performed to describe the flow in high Reynolds number (106) regime and two k-ε turbulence models (standard k-ε and RNG k-ε) have been used for computing the turbulent flow. A streamline upwind Petrov Galerkin technique has been used for spatial discretisation. The velocity field and the pressure distribution inside the spiral casing has been studied. It has been observed that very strong secondary flow is evolved on the cross-stream planes.


2013 ◽  
Vol 721 ◽  
pp. 58-85 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall ◽  
Andrew Walton

AbstractThe recently understood relationship between high-Reynolds-number vortex–wave interaction theory and computationally generated self-sustaining processes provides a possible route to an understanding of some of the underlying structures of fully turbulent flows. Here vortex–wave interaction (VWI) theory is used in the long streamwise wavelength limit to continue the development found at order-one wavelengths by Hall & Sherwin (J. Fluid Mech., vol. 661, 2010, pp. 178–205). The asymptotic description given reduces the Navier–Stokes equations to the so-called boundary-region equations, for which we find equilibrium states describing the change in the VWI as the wavelength of the wave increases from $O(h)$ to $O(Rh)$, where $R$ is the Reynolds number and $2h$ is the depth of the channel. The reduced equations do not include the streamwise pressure gradient of the perturbation or the effect of streamwise diffusion of the wave–vortex states. The solutions we calculate have an asymptotic error proportional to ${R}^{- 2} $ when compared to the full Navier–Stokes equations. The results found correspond to the minimum drag configuration for VWI states and might therefore be of relevance to the control of turbulent flows. The key feature of the new states discussed here is the thickening of the critical layer structure associated with the wave part of the flow to completely fill the channel, so that the roll part of the flow is driven throughout the flow rather than as in Hall & Sherwin as a stress discontinuity across the critical layer. We identify a critical streamwise wavenumber scaling, which, when approached, causes the flow to localize and take on similarities with computationally generated or experimentally observed turbulent spots. In effect, the identification of this critical wavenumber for a given value of the assumed high Reynolds number fixes a minimum box length necessary for the emergence of localized structures. Whereas nonlinear equilibrium states of the Navier–Stokes equations are thought to form a backbone on which turbulent flows hang, our results suggest that the localized states found here might play a related role for turbulent spots.


Author(s):  
Y Wang ◽  
S Komori

A pressure-based finite volume procedure developed previously for incompressible flows is extended to predict the three-dimensional compressible flow within a centrifugal impeller. In this procedure, the general curvilinear coordinate system is used and the collocated grid arrangement is adopted. Mass-averaging is used to close the instantaneous Navier-Stokes equations. The covariant velocity components are used as the main variables for the momentum equations, making the pressure-velocity coupling easier. The procedure is successfully applied to predict various compressible flows from subsonic to supersonic. With the aid of the k-ɛ turbulence model, the flow details within a centrifugal impeller are obtained using the present procedure. Predicted distributions of the meridional velocity and the static pressure are reasonable. Calculated radial velocities and flow angles are favourably compared with the measurements at the exit of the impeller.


2000 ◽  
Vol 123 (3) ◽  
pp. 680-685 ◽  
Author(s):  
L. He ◽  
K. Sato

A three-dimensional incompressible viscous flow solver of the thin-layer Navier-Stokes equations was developed for the unsteady turbomachinery flow computations. The solution algorithm for the unsteady flows combines the dual time stepping technique with the artificial compressibility approach for solving the incompressible unsteady flow governing equations. For time accurate calculations, subiterations are introduced by marching the equations in the pseudo-time to fully recover the incompressible continuity equation at each real time step, accelerated with a multi-grid technique. Computations of test cases show satisfactory agreements with corresponding theoretical and experimental results, demonstrating the validity and applicability of the present method to unsteady incompressible turbomachinery flows.


Author(s):  
Yannis Kallinderis ◽  
Hyung Taek Ahn

Numerical prediction of vortex-induced vibrations requires employment of the unsteady Navier-Stokes equations. Current Navier-Stokes solvers are quite expensive for three-dimensional flow-structure applications. Acceptance of Computational Fluid Dynamics as a design tool for the offshore industry requires improvements to current CFD methods in order to address the following important issues: (i) stability and computation cost of the numerical simulation process, (ii) restriction on the size of the allowable time-step due to the coupling of the flow and structure solution processes, (iii) excessive number of computational elements for 3-D applications, and (iv) accuracy and computational cost of turbulence models used for high Reynolds number flow. The above four problems are addressed via a new numerical method which employs strong coupling between the flow and the structure solutions. Special coupling is also employed between the Reynolds-averaged Navier-Stokes equations and the Spalart-Allmaras turbulence model. An element-type independent spatial discretization scheme is also presented which can handle general hybrid meshes consisting of hexahedra, prisms, pyramids, and tetrahedral.


Sign in / Sign up

Export Citation Format

Share Document