scholarly journals Effect of Local Turbulent Structures on the Hot Jet in Transitional Flow

2021 ◽  
Vol 2119 (1) ◽  
pp. 012009
Author(s):  
A Sakhnov ◽  
V V Lukashov

Abstract Turbulent parts localized in flow direction may arise in a pipe with transitional regime of the stable laminar Poiseuille flow. A key condition for occurrence of such structures is a pipe with rather long length relative to its diameter. Our paper presents numerical modelling of the hot air jet flowing from the long pipe into the cold open volume at Re=2426. The modelling was performed in OpenFOAM software on the basis of the large eddy simulation (LES) method. The WALE (Wall-adapting local eddy-viscosity) model was used for closure of Navier-Stokes equations on subgrid scales. We demonstrated that local turbulent structures have a weak effect on the hot jet at flowing into the cold open volume.

Author(s):  
Leopold Sˇkerget ◽  
Jure Ravnik

Detached turbulent flows are difficult to predict numerically and often serve as benchmark cases for developing new numerical schemes and new turbulent models. Turbulent flow over periodic hills is one such examples, since the flow exhibits separation and reattachment on a smoothly and/or sharp curved geometry, strong pressure gradients and fluctuation of the separation point in time. These cases have been chosen by many authors for testing different turbulence simulation approaches. When the bottom wall is heated, the complexity of the problem increased, since convective heat transfer is defined by small scale turbulent structures close to the wall. We developed a Reynolds-Averaged Navier-Stokes and Large Eddy Simulation solver based on the velocity-vorticity formulation of Navier Stokes equations. RANS equations are coupled by a low-Reynolds number turbulent model, while Smagorinsky subgrid model is used for LES. The governing equations are solved with a numerical solution algorithm, which is based on the boundary element method. The pressure field is computed in a post processing step by solving a Poisson equation. The single domain as well as domain decomposition approaches are applied. The developed method was validated using flow over periodic hills test case.


2011 ◽  
Vol 134 (2) ◽  
Author(s):  
Paul Tucker ◽  
Simon Eastwood ◽  
Christian Klostermeier ◽  
Richard Jefferson-Loveday ◽  
James Tyacke ◽  
...  

Unlike Reynolds-averaged Navier–Stokes (RANS) models that need calibration for different flow classes, LES (where larger turbulent structures are resolved by the grid and smaller modeled in a fashion reminiscent of RANS) offers the opportunity to resolve geometry dependent turbulence as found in complex internal flows—albeit at substantially higher computational cost. Based on the results for a broad range of studies involving different numerical schemes, large eddy simulation (LES) models and grid topologies, an LES hierarchy and hybrid LES related approach is proposed. With the latter, away from walls, no LES model is used, giving what can be termed numerical LES (NLES). This is relatively computationally efficient and makes use of the dissipation present in practical industrial computational fluid dynamics (CFD) programs. Near walls, RANS modeling is used to cover over numerous small structures, the LES resolution of which is generally intractable with current computational power. The linking of the RANS and NLES zones through a Hamilton–Jacobi equation is advocated. The RANS-NLES hybridization makes further sense for compressible flow solvers, where, as the Mach number tends to zero at walls, excessive dissipation can occur. The hybrid strategy is used to predict flow over a rib roughened surface and a jet impinging on a convex surface. These cases are important for blade cooling and show encouraging results. Further results are presented in a companion paper.


2011 ◽  
Vol 133 (5) ◽  
Author(s):  
F. P. P. Tan ◽  
N. B. Wood ◽  
G. Tabor ◽  
X. Y. Xu

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y+ less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y+ less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.


Author(s):  
Hongsik Im ◽  
Xiangying Chen ◽  
Gecheng Zha

Detached eddy simulation of an aeroelastic self-excited instability, flutter in NASA Rotor 67 is conducted using a fully coupled fluid/structre interaction. Time accurate compressible 3D Navier-Stokes equations are solved with a system of 5 decoupled modal equations in a fully coupled manner. The 5th order WENO scheme for the inviscid flux and the 4th order central differencing for the viscous flux are used to accurately capture interactions between the flow and vibrating blades with the DES (detached eddy simulation) of turbulence. A moving mesh concept that can improve mesh quality over the rotor tip clearance was implemented. Flutter simulations were first conducted from choke to stall using 4 blade passages. Stall flutter initiated at rotating stall onset, grows dramatically with resonance. The frequency analysis shows that resonance occurs at the first mode of the rotor blade. Before stall, the predicted responses of rotor blades decayed with time, resulting in no flutter. Full annulus simulation at peak point verifies that one can use the multi-passage approach with periodic boundary for the flutter prediction.


2013 ◽  
Vol 135 (7) ◽  
Author(s):  
A. Ghidoni ◽  
A. Colombo ◽  
S. Rebay ◽  
F. Bassi

In the last decade, discontinuous Galerkin (DG) methods have been the subject of extensive research efforts because of their excellent performance in the high-order accurate discretization of advection-diffusion problems on general unstructured grids, and are nowadays finding use in several different applications. In this paper, the potential offered by a high-order accurate DG space discretization method with implicit time integration for the solution of the Reynolds-averaged Navier–Stokes equations coupled with the k-ω turbulence model is investigated in the numerical simulation of the turbulent flow through the well-known T106A turbine cascade. The numerical results demonstrate that, by exploiting high order accurate DG schemes, it is possible to compute accurate simulations of this flow on very coarse grids, with both the high-Reynolds and low-Reynolds number versions of the k-ω turbulence model.


2002 ◽  
Vol 124 (2) ◽  
pp. 413-423 ◽  
Author(s):  
L. S. Hedges ◽  
A. K. Travin ◽  
P. R. Spalart

The flow around a generic airliner landing-gear truck is calculated using the methods of Detached-Eddy Simulation, and of Unsteady Reynolds-Averaged Navier-Stokes Equations, with the Spalart-Allmaras one-equation model. The two simulations have identical numerics, using a multi-block structured grid with about 2.5 million points. The Reynolds number is 6×105. Comparison to the experiment of Lazos shows that the simulations predict the pressure on the wheels accurately for such a massively separated flow with strong interference. DES performs somewhat better than URANS. Drag and lift are not predicted as well. The time-averaged and instantaneous flow fields are studied, particularly to determine their suitability for the physics-based prediction of noise. The two time-averaged flow fields are similar, though the DES shows more turbulence intensity overall. The instantaneous flow fields are very dissimilar. DES develops a much wider range of unsteady scales of motion and appears promising for noise prediction, up to some frequency limit.


Author(s):  
Jagannath Mahato ◽  
Dhananjay Kumar Srivastava ◽  
Dinesh Kumar Chandraker ◽  
Rajaram Lakkaraju

Abstract Investigations on flow dynamics of a compound droplet have been carried out in a two-dimensional fully-developed Poiseuille flow by solving the Navier-Stokes equations with the evolution of the droplet using the volume of fluid method with interface compression. The outer droplet undergoes elongation similar to a simple droplet of same size placed under similar ambient condition in the flow direction, but, the inner droplet evolves in compressed form. The compound droplet is varied starting from the centerline towards the walls of the channel. The simulations showed that on applying an offset, asymmetric slipper-like shapes are observed as opposed to symmetric bullet-like shapes through the centerline. Temporal dynamics, deformation patterns, and droplet shell pinch-off mode vary with the offset, with induction of lateral migration. Also, investigations are done on the effect of various parameters like droplet size, Capillary number, and viscosity ratio on the deformation magnitude and lateral migration.


2018 ◽  
Author(s):  
Jiajun Chen ◽  
Yue Sun ◽  
Hang Zhang ◽  
Dakui Feng ◽  
Zhiguo Zhang

Mixing in pipe junctions can play an important role in exciting force and distribution of flow in pipe network. This paper investigated the cross pipe junction and proposed an improved plan, Y-shaped pipe junction. The numerical study of a three-dimensional pipe junction was performed for calculation and improved understanding of flow feature in pipe. The filtered Navier–Stokes equations were used to perform the large-eddy simulation of the unsteady incompressible flow in pipe. From the analysis of these results, it clearly appears that the vortex strength and velocity non-uniformity of centerline, can be reduced by Y-shaped junction. The Y-shaped junction not only has better flow characteristic, but also reduces head loss and exciting force. The results of the three-dimensional improvement analysis of junction can be used in the design of pipe network for industry.


2014 ◽  
Author(s):  
P. Bigay ◽  
A. Bardin ◽  
G. Oger ◽  
D. Le Touzé

In order to efficiently address complex problems in hydrodynamics, the advances in the development of a new method are presented here. This method aims at finding a good compromise between computational efficiency, accuracy, and easy handling of complex geometries. The chosen method is an Explicit Cartesian Finite Volume method for Hydrodynamics (ECFVH) based on a compressible (hyperbolic) solver, with a ghost-cell method for geometry handling and a Level-set method for the treatment of biphase-flows. The explicit nature of the solver is obtained through a weakly-compressible approach chosen to simulate nearly-incompressible flows. The explicit cell-centered resolution allows for an efficient solving of very large simulations together with a straightforward handling of multi-physics. A characteristic flux method for solving the hyperbolic part of the Navier-Stokes equations is used. The treatment of arbitrary geometries is addressed in the hyperbolic and viscous framework. Viscous effects are computed via a finite difference computation of viscous fluxes and turbulent effects are addressed via a Large-Eddy Simulation method (LES). The Level-Set solver used to handle biphase flows is also presented. The solver is validated on 2-D test cases (flow past a cylinder, 2-D dam break) and future improvements are discussed.


Inventions ◽  
2020 ◽  
Vol 5 (3) ◽  
pp. 33
Author(s):  
Horng-Wen Wu ◽  
Tang-Hong Chen ◽  
Nugroho-Putra Kelana ◽  
De-An Huang

This study analyzes transient turbulent modeling of three-dimensional multiple dimpled fin array using large eddy simulation (LES). The Navier–Stokes equations as well as the energy equation were constructed by the finite volume method and then discretized to form algebraic equations, which were solved by semi-implicit method for pressure-linked equation (SIMPLE). The solutions of temperature and velocity were obtained by iterating computation until it converged within each step. This simulation places nine fins on the bottom surface of a channel and changes the height of the dimple (0.4, 0.8, and 1.2 mm) with three different levels of Reynolds number (Re) (3500, 5000, and 6500) to investigate the temperature and flow field without gravity in forced convection. The results indicate that the dimpled fin array can generate vortices between the convex/concave dimples and the fin base and increase the influences of the height of the dimple on the flow field around the fin array. The averaged time-mean of the Nusselt number (Nu) for the dimple height of 0.8 mm is higher than that of the no-dimple case up to 14.4%, while the averaged time-mean Nu for the dimple height of 1.2 mm is lower than that of the no-dimple case up to 11.6%.


Sign in / Sign up

Export Citation Format

Share Document