scholarly journals Dynamical Simulation and Statistical Analysis of Velocity Fluctuations of a Turbulent Flow behind a Cube

2007 ◽  
Vol 2007 ◽  
pp. 1-28 ◽  
Author(s):  
T. F. Oliveira ◽  
R. B. Miserda ◽  
F. R. Cunha

A statistical approach for the treatment of turbulence data generated by computer simulations is presented. A model for compressible flows at large Reynolds numbers and low Mach numbers is used for simulating a backward-facing step airflow. A scaling analysis has justified the commonly used assumption that the internal energy transport due to turbulent velocity fluctuations and the work done by the pressure field are the only relevant mechanisms needed to model subgrid-scale flows. From the numerical simulations, the temporal series of velocities are collected for ten different positions in the flow domain, and are statistically treated. The statistical approach is based on probability averages of the flow quantities evaluated over several realizations of the simulated flow. We look at how long of a time average is necessary to obtain well-converged statistical results. For this end, we evaluate the mean-square difference between the time average and an ensemble average as the measure of convergence. This is an interesting question since the validity of the ergodic hypothesis is implicitly assumed in every turbulent flow simulation and its analysis. The ergodicity deviations from the numerical simulations are compared with theoretical predictions given by scaling arguments. A very good agreement is observed. Results for velocity fluctuations, normalized autocorrelation functions, power spectra, probability density distributions, as well as skewness and flatness coefficients are also presented.

Author(s):  
Xu Zhang ◽  
Dan Stanescu ◽  
Jonathan W. Naughton

This paper describes a turbulent flow simulation method, which is based on combination of spectral element and large eddy simulation (LES) technique. The robust, high-order discontinuous Galerkin (DG) spectral element method for large-eddy simulation of compressible flows allows for arbitrary order of accuracy and has excellent stability properties. A local spectral discretization in terms of Legendre polynomials is used on each element of the (possibly unstructured) mesh, which allows for high-accurate simulations of turbulent flows. Discontinuities across the interfaces of the elements are resolved using a Riemann solver. An isoparametric representation of the geometry is implemented, with boundaries of the domain discretized to the same order of accuracy as the solution, and explicit low-storage Runge-Kutta methods are used for time integration. Large eddy simulation has proven to be a valuable technique for the calculation of turbulent flows. An element based filtering technique is used in conjunction with the standard Smagorinsky eddy viscosity model to estimate the effect of sub-grid scales stresses in this paper. The recently developed nonlinear model [1] will also be added in the future. The final aim of this project is to use the LES methodology in swirling jet flow simulation. As a first step towards these simulations, simulations of compressible turbulent mixing layer and back-facing step are also performed to evaluate the robust method. Initial results based on both DNS and large eddy simulations are presented in this paper. Future work will be to validate the code.


1996 ◽  
Vol 07 (04) ◽  
pp. 543-561 ◽  
Author(s):  
WOLFGANG KALTHOFF ◽  
STEFAN SCHWARZER ◽  
GERALD RISTOW ◽  
HANS J. HERRMANN

We present a numerical method to deal efficiently with large numbers of particles in incompressible fluids. The interactions between particles and fluid are taken into account by a physically motivated ansatz based on locally defined drag forces. We demonstrate the validity of our approach by performing numerical simulations of sedimenting non-Brownian spheres in two spatial dimensions and compare our results with experiments. Our method reproduces qualitatively important aspects of the experimental findings, in particular the strong anisotropy of the hydrodynamic bulk self-diffusivities.


Author(s):  
J. Szydlowski ◽  
M. Costes

This paper presents numerical simulations of the flow around a NACA 0015 airfoil at static and dynamic stall. The treatment of these configurations is a very challenging task for CFD applications. The turbulent flow around the static and in pitch oscillation airfoil is computed using different approaches: 2D RANS, 3D RANS and DES methodologies and with finer and finer meshes in order to try to reach a space converged solution. The main conclusion of the paper is that the prediction of static and all the more dynamic stall is not mature with present modeling capabilities.


2021 ◽  
Vol 3 (8) ◽  
Author(s):  
Erik Flídr ◽  
Petr Straka ◽  
Milan Kladrubský ◽  
Tomáš Jelínek

AbstractThis contribution describes experimental and numerical research of an unsteady behaviour of a flow in an end-wall region of a linear nozzle cascade. Effects of compressibility ($$M_\mathrm {2,is}$$ M 2 , is ) and inlet flow angle ($$\alpha _1$$ α 1 ) were investigated. Reynolds number ($$Re_\mathrm {2,is}$$ R e 2 , is $$=8.5\times 10^5$$ = 8.5 × 10 5 ) was held constant for all tested cases. Unsteady pressure measurement was performed at the blade mid-span in the identical position $${\mathfrak {s}}$$ s to obtain reference data. Surface flow visualizations were performed as well as the steady pressure measurement to support conclusions obtained from the unsteady measurements. Comparison of the surface Mach number distributions obtained from the experiments and from the numerical simulations are presented. Flow visualizations are then compared with calculated limiting streamlines on the blade suction surface. It was shown, that the flow structures in the end-wall region were not affected by the primary flow at the blade mid-span, even when the shock wave formed. This conclusion was made from the experimental, numerical, steady as well as unsteady points of view. Three significant frequencies in the power spectra suggested that there was a periodical interaction between the vortex structures in the end-wall region. Based on the data analyses, anisotropic turbulence was observed in the cascade.


2021 ◽  
Vol 925 ◽  
Author(s):  
Kaushlendra Dubey ◽  
Sanjeev Sanghi ◽  
Amit Gupta ◽  
Supreet Singh Bahga

We present an experimental and numerical investigation of electrokinetic instability (EKI) in microchannel flow with streamwise conductivity gradients, such as those observed during sample stacking in capillary electrophoresis. A plug of a low-conductivity electrolyte solution is initially sandwiched between two high-conductivity zones in a microchannel. This spatial conductivity gradient is subjected to an external electric field applied along the microchannel axis, and for sufficiently strong electric fields an instability sets in. We have explored the physics of this EKI through experiments and numerical simulations, and supplemented the results using scaling analysis. We performed EKI experiments at different electric field values and visualised the flow using a passive fluorescent tracer. The experimental data were analysed using the proper orthogonal decomposition technique to obtain a quantitative measure of the threshold electric field for the onset of instability, along with the corresponding coherent structures. To elucidate the physical mechanism underlying the instability, we performed high-resolution numerical simulations of ion transport coupled with fluid flow driven by the electric body force. Simulations reveal that the non-uniform electroosmotic flow due to axially varying conductivity field causes a recirculating flow within the low-conductivity region, and creates a new configuration wherein the local conductivity gradients are orthogonal to the applied electric field. This configuration leads to EKI above a threshold electric field. The spatial features of the instability predicted by the simulations and the threshold electric field are in good agreement with the experimental observations and provide useful insight into the underlying mechanism of instability.


Author(s):  
Veeraraghava R Hasti ◽  
Prithwish Kundu ◽  
Sibendu Som ◽  
Jay P Gore

The turbulent flow field in a practical gas turbine combustor is very complex because of the interactions between various flows resulting from components like multiple types of swirlers, dilution holes, and liner effusion cooling holes. Numerical simulations of flows in such complex combustor configurations are challenging. The challenges result from (a) the complexities of the interfaces between multiple three-dimensional shear layers, (b) the need for proper treatment of a large number of tiny effusion holes with multiple angles, and (c) the requirements for fast turnaround times in support of engineering design optimization. Both the Reynolds averaged Navier–Stokes simulation (RANS) and the large eddy simulation (LES) for the practical combustor geometry are considered. An autonomous meshing using the cut-cell Cartesian method and adaptive mesh refinement (AMR) is demonstrated for the first time to simulate the flow in a practical combustor geometry. The numerical studies include a set of computations of flows under a prescribed pressure drop across the passage of interest and another set of computations with all passages open with a specified total flow rate at the plenum inlet and the pressure at the exit. For both sets, the results of the RANS and the LES flow computations agree with each other and with the corresponding measurements. The results from the high-resolution LES simulations are utilized to gain fundamental insights into the complex turbulent flow field by examining the profiles of the velocity, the vorticity, and the turbulent kinetic energy. The dynamics of the turbulent structures are well captured in the results of the LES simulations.


2021 ◽  
Vol 931 ◽  
Author(s):  
D. Li ◽  
J. Komperda ◽  
A. Peyvan ◽  
Z. Ghiasi ◽  
F. Mashayek

The present paper uses the detailed flow data produced by direct numerical simulation (DNS) of a three-dimensional, spatially developing plane free shear layer to assess several commonly used turbulence models in compressible flows. The free shear layer is generated by two parallel streams separated by a splitter plate, with a naturally developing inflow condition. The DNS is conducted using a high-order discontinuous spectral element method (DSEM) for various convective Mach numbers. The DNS results are employed to provide insights into turbulence modelling. The analyses show that with the knowledge of the Reynolds velocity fluctuations and averages, the considered strong Reynolds analogy models can accurately predict temperature fluctuations and Favre velocity averages, while the extended strong Reynolds analogy models can correctly estimate the Favre velocity fluctuations and the Favre shear stress. The pressure–dilatation correlation and dilatational dissipation models overestimate the corresponding DNS results, especially with high compressibility. The pressure–strain correlation models perform excellently for most pressure–strain correlation components, while the compressibility modification model gives poor predictions. The results of an a priori test for subgrid-scale (SGS) models are also reported. The scale similarity and gradient models, which are non-eddy viscosity models, can accurately reproduce SGS stresses in terms of structure and magnitude. The dynamic Smagorinsky model, an eddy viscosity model but based on the scale similarity concept, shows acceptable correlation coefficients between the DNS and modelled SGS stresses. Finally, the Smagorinsky model, a purely dissipative model, yields low correlation coefficients and unacceptable accumulated errors.


Author(s):  
Peng Zhang ◽  
Yu Rao ◽  
Yanlin Li

This paper presents a numerical study on turbulent flow and heat transfer in the channels with a novel hybrid cooling structure with miniature V-shaped ribs and dimples on one wall. The heat transfer characteristics, pressure loss and turbulent flow structures in the channels with the rib-dimples with three different rib heights of 0.6 mm, 1.0 mm and 1.5 mm are obtained for the Reynolds numbers ranging from 18,700 to 60,000 by numerical simulations, which are also compared with counterpart of a pure dimpled and pure V ribbed channel. The results show that the overall Nusselt numbers of the V rib-dimple channel with the rib height of 1.5 mm is up to 70% higher than that of the channels with pure dimples. The numerical simulations show that the arrangement of the miniature V rib upstream each dimple induces complex secondary flow near the wall and generates downwashing vortices, which intensifies the flow mixing and turbulent kinetic energy in the dimple, resulting in significant improvement in heat transfer enhancement and uniformness.


2009 ◽  
Vol 630 ◽  
pp. 1-4 ◽  
Author(s):  
IVAN MARUSIC

Turbulent flows near walls have been the focus of intense study since their first description by Ludwig Prandtl over 100 years ago. They are critical in determining the drag and lift of an aircraft wing for example. Key challenges are to understand the physical mechanisms causing the transition from smooth, laminar flow to turbulent flow and how the turbulence is then maintained. Recent direct numerical simulations have contributed significantly towards this understanding.


Sign in / Sign up

Export Citation Format

Share Document