Computing Blunt Body Flows on Coarse Grids Using Vorticity Confinement

2002 ◽  
Vol 124 (4) ◽  
pp. 876-885 ◽  
Author(s):  
M. Fan ◽  
Y. Wenren ◽  
W. Dietz ◽  
M. Xiao ◽  
J. Steinhoff

Over the last few years, a new flow computational methodology, vorticity confinement, has been shown to be very effective in treating concentrated vortical regions. These include thin vortex filaments which can be numerically convected over arbitrary distances on coarse Eulerian grids, while requiring only ∼2 grid cells across their cross section. They also include boundary layers on surfaces “immersed” in nonconforming uniform Cartesian grids, with no requirement for grid refinement or complex logic near the surface. In this paper we use vorticity confinement to treat flow over blunt bodies, including attached and separating boundary layers, and resulting turbulent wakes. In the wake it serves as a new, simple effective large-eddy simulation (LES). The same basic idea is applied to all of these features: At the smallest scales (∼2 cells) the vortical structures are captured and treated, effectively, as solitary waves that are solutions of nonlinear discrete equations on the grid. The method does not attempt to accurately discretize the Euler/Navier-Stokes partial differential equations (pde’s) for these small scales, but, rather, serves as an implicit, nonlinear model of the structures, directly on the grid. The method also allows the boundary layer to be effectively “captured.” In the turbulent wake, where there are many scales, small structures represent an effective small scale energy sink. However, they do not have the unphysical spreading due to numerical diffusion at these scales, which is present in conventional computational methods. The basic modeling idea is similar to that used in shock capturing, where intrinsically discrete equations are satisfied in thin, modeled regions. It is argued that, for realistic high Reynolds number flows, this direct, grid-based modeling approach is much more effective than first formulating model pde’s for the small scale, turbulent vortical regions and then discretizing them. Results are presented for three-dimensional flows over round and square cylinders and a realistic helicopter landing ship. Comparisons with experimental data are given. Finally, a new simpler formulation of vorticity confinement is given together with a related formulation for confinement of passive scalar fields.

1990 ◽  
Vol 216 ◽  
pp. 1-34 ◽  
Author(s):  
Rahul R. Prasad ◽  
K. R. Sreenivasan

The three-dimensional turbulent field of a passive scalar has been mapped quantitatively by obtaining, effectively instantaneously, several closely spaced parallel two-dimensional images; the two-dimensional images themselves have been obtained by laser-induced fluorescence. Turbulent jets and wakes at moderate Reynolds numbers are used as examples. The working fluid is water. The spatial resolution of the measurements is about four Kolmogorov scales. The first contribution of this work concerns the three-dimensional nature of the boundary of the scalar-marked regions (the ‘scalar interface’). It is concluded that interface regions detached from the main body are exceptional occurrences (if at all), and that in spite of the large structure, the randomness associated with small-scale convolutions of the interface are strong enough that any two intersections of it by parallel planes are essentially uncorrelated even if the separation distances are no more than a few Kolmogorov scales. The fractal dimension of the interface is determined directly by box-counting in three dimensions, and the value of 2.35 ± 0.04 is shown to be in good agreement with that previously inferred from two-dimensional sections. This justifies the use of the method of intersections. The second contribution involves the joint statistics of the scalar field and the quantity χ* (or its components), χ* being the appropriate approximation to the scalar ‘dissipation’ field in the inertial–convective range of scales. The third aspect relates to the multifractal scaling properties of the spatial intermittency of χ*; since all three components of χ* have been obtained effectively simultaneously, inferences concerning the scaling properties of the individual components and their sum have been possible. The usefulness of the multifractal approach for describing highly intermittent distributions of χ* and its components is explored by measuring the so-called singularity spectrum (or the f(α)-curve) which quantifies the spatial distribution of various strengths of χ*. Also obtained is a time sequence of two-dimensional images with the temporal resolution on the order of a few Batchelor timescales; this enables us to infer features of temporal intermittency in turbulent flows, and qualitatively the propagation speeds of the scalar interface. Finally, a few issues relating to the resolution effects have been addressed briefly by making point measurements with the spatial and temporal resolutions comparable with the Batchelor lengthscale and the corresponding timescale.


2013 ◽  
Vol 723 ◽  
pp. 404-428 ◽  
Author(s):  
D. A. Philips ◽  
R. Rossi ◽  
G. Iaccarino

AbstractResults from large-eddy simulations of short-range dispersion of a passive scalar from a point source release in an urban-like canopy are presented. The computational domain is that of a variable height array of buildings immersed in a pressure-driven, turbulent flow with a roughness Reynolds number ${\mathit{Re}}_{\tau } = 433$. A comparative study of several cases shows the changes in plume behaviour for different mean flow directions and source locations. The analysis of the results focuses on utilizing the high-fidelity datasets to examine the three-dimensional flow field and scalar plume structure. The detailed solution of the flow and scalar fields within the canopy allows for a direct assessment of the impact of local features of the building array geometry. The staggered, skewed and aligned arrangements of the buildings with respect to the oncoming flow were shown to affect plume development. Additional post-processing quantified this development through parameters fundamental to reduced-order Gaussian dispersion models. The parameters include measures of concentration decay with distance from the source as well as plume trajectory and spread. The horizontal plume trajectory and width were found to be more sensitive to source location variations, and hence local geometric features, than vertical plume parameters.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


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.


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%.


2017 ◽  
Vol 14 (03) ◽  
pp. 1750021 ◽  
Author(s):  
A. Niktash ◽  
B. P. Huynh

A windcatcher is a structure for providing natural ventilation using wind power; it is usually fitted on the roof of a building to exhaust the inside stale air to the outside and supplies the outside fresh air into the building interior space working by pressure difference between outside and inside of the building. In this paper, the behavior of free wind flow through a three-dimensional room fitted with a centered position two-canal bottom shape windcatcher model is investigated numerically, using a commercial computational fluid dynamics (CFD) software package and LES (Large Eddy Simulation) CFD method. The results have been compared with the obtained results for the same model but using RANS (Reynolds Averaged Navier–Stokes) CFD method. The model with its surrounded space has been considered in both method. It is found that the achieved results for the model from LES method are in good agreement with RANS method’s results for the same model.


Author(s):  
Hayder Schneider ◽  
Dominic von Terzi ◽  
Hans-Jo¨rg Bauer ◽  
Wolfgang Rodi

Reynolds-Averaged Navier-Stokes (RANS) calculations and Large-Eddy Simulations (LES) of the flow in two asymmetric three-dimensional diffusers were performed. The numerical setup was chosen to be in compliance with previous experiments. The aim of the present study is to find the least expensive method to compute reliably and accurately the impact of geometric sensitivity on the flow. RANS calculations fail to predict both the extent and location of the three-dimensional separation bubble. In contrast, LES is able to determine the amount of reverse flow and the pressure coefficient within the accuracy of experimental data.


Processes ◽  
2020 ◽  
Vol 8 (12) ◽  
pp. 1687
Author(s):  
Chao Yu ◽  
Xiangyao Xue ◽  
Kui Shi ◽  
Mingzhen Shao ◽  
Yang Liu

This paper compares the performances of three Computational Fluid Dynamics (CFD) turbulence models, Reynolds Average Navier-Stokes (RANS), Detached Eddy Simulation (DES), and Large Eddy Simulation (LES), for simulating the flow field of a wheel loader engine compartment. The distributions of pressure fields, velocity fields, and vortex structures in a hybrid-grided engine compartment model are analyzed. The result reveals that the LES and DES can capture the detachment and breakage of the trailing edge more abundantly and meticulously than RANS. Additionally, by comparing the relevant calculation time, the feasibility of the DES model is proved to simulate the three-dimensional unsteady flow of engine compartment efficiently and accurately. This paper aims to provide a guiding idea for simulating the transient flow field in the engine compartment, which could serve as a theoretical basis for optimizing and improving the layout of the components of the engine compartment.


1987 ◽  
Vol 184 ◽  
pp. 207-243 ◽  
Author(s):  
Ralph W. Metcalfe ◽  
Steven A. Orszag ◽  
Marc E. Brachet ◽  
Suresh Menon ◽  
James J. Riley

The three-dimensional stability of two-dimensional vortical states of planar mixing layers is studied by direct numerical integration of the Navier-Stokes equations. Small-scale instabilities are shown to exist for spanwise scales at which classical linear modes are stable. These modes grow on convective timescales, extract their energy from the mean flow and exist at moderately low Reynolds numbers. Their growth rates are comparable with the most rapidly growing inviscid instability and with the growth rates of two-dimensional subharmonic (pairing) modes. At high amplitudes, they can evolve into pairs of counter-rotating, streamwise vortices, connecting the primary spanwise vortices, which are very similar to the structures observed in laboratory experiments. The three-dimensional modes do not appear to saturate in quasi-steady states as do the purely two-dimensional fundamental and subharmonic modes in the absence of pairing. The subsequent evolution of the flow depends on the relative amplitudes of the pairing modes. Persistent pairings can inhibit three-dimensional instability and, hence, keep the flow predominantly two-dimensional. Conversely, suppression of the pairing process can drive the three-dimensional modes to more chaotic, turbulent-like states. An analysis of high-resolution simulations of fully turbulent mixing layers confirms the existence of rib-like structures and that their coherence depends strongly on the presence of the two-dimensional pairing modes.


Sign in / Sign up

Export Citation Format

Share Document