DSMC simulation of transition and turbulent flow in a Lid-driven cavity at high Mach number

Pradhan, S.K. and Kumaran, V. (2013) DSMC simulation of transition and turbulent flow in a Lid-driven cavity at high Mach number. In: ETC 2013 - 14th European Turbulence Conference, 1 September 2013 through 4 September 2013, France.

PDF ETC 2013.pdf - Published Version Restricted to Registered users only Download (10MB) | Request a copy

Abstract

The flow in a 2D and 3D lid-driven cavity, with two opposite walls moving in opposite directions with equal velocity, simulated using the direct simulation Monte Carlo (DSMC) method, has been used as a test bed for examining different aspects of transition and turbulence 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 at high Mach M = (Uw /sqrt(kB Tw /m)) and Reynolds numbers Re = (�av Uw Ly/μw). Here, Ly is the smallest box dimension, Uw and Tw are the wall velocity and temperature, �av is the volume-averaged gas density, and μw is the gas viscosity at the temperature corresponding to the wall temperature. The transition is found to be highly sub-critical in 2D even at high Mach number, with well-separated lower and upper critical Reynolds numbers for the transition from turbulent-laminar and laminar-turbulent transitions. The transition Reynolds number increases faster than linearly with Mach number. This implies that the Knudsen number at transition (ratio of mean-free-path and system size, also proportional to the ratio of Mach and Reynolds numbers) passes through a maximum as the Mach number is increased. This maximum value is small, less than 0.025, indicating that transition is a continuum phenomenon even at high Mach numbers. The transition to turbulence is sub-critical in 3D as well. The transition Reynolds number does increase with Mach number, and the Knudsen number increases monotonically with Mach number over the parameter range studied here. In the present analysis of high Mach turbulent flows using DSMC method, wall slip in the temperature and the velocities are found to be significant. Slip occurs because the temperature/velocity of the molecules incident on the wall could be very different from that of the the wall, even though the temperature/velocity of the reflected molecules is equal to that of the wall (Pradhan & Kumaran, J. Fluid Mech-2011). There is slip even in the mean velocity as well as the intensity of the turbulent velocity fluctuations tangential to the wall. Due to the non-zero variation in the tangential velocities, it is found that the wall-normal fluctuating velocity increases proportional to the distance from the wall, in contrast to an incompressible flow with no-slip boundary conditions where the velocity increases as the square of the distance from the wall. We find that the amplitudes of the tangential fluctuating velocities increase close to the one-third power of the distance from the wall. In a compressible turbulent flow, we examine the result that the ratio of the mean free path (λ) and Kolmogorov scale (η) increases proportional as (M/Re1/4), and it increases asymptotically with Mach number in the high Mach number limit. The simulation show that the ratio does decrease as Re�1/4, but it does not increase linearly with Mach number. The resolution suggested by our simulation is that even though the Mach number based on the wall velocity and temperature is large, the local Mach number based on the local dissipation velocity in regions of high shear decreases due to an increase in temperature. Due to this, the ratio of the mean free path and Kolmogorov scale appears to taper off in the high Mach number limit. An important finding is that the ratio of the mean free path and Kolmogorov scale shows very little variation across the domain, even though the mean free path and Kolmogorov scale individually show larger variations. The ratio of the mean free path and Kolmogorov scale was also shown to be equal to the local Mach number based on the local dissipation velocity, and can also be interpreted as the square root of the ratio of the strain rate and collision frequency. All of these quantities are remarkably invariant across the domain, indicating a coupling between the local temperature and the dissipation rate in a high Mach number turbulent flow. Copyright © ETC 2013 - 14th European Turbulence Conference.All rights reserved.

Item Type:	Conference Paper
Publication:	ETC 2013 - 14th European Turbulence Conference
Publisher:	Zakon Group LLC
Additional Information:	cited By 0; Conference of 14th European Turbulence Conference, ETC 2013 ; Conference Date: 1 September 2013 Through 4 September 2013; Conference Code:157594
Keywords:	Aerodynamics; Compressible flow; Computational complexity; Density of gases; Incompressible flow; Molecules; Monte Carlo methods; Pipe flow; Prandtl number; Reynolds number; Shear flow; Strain rate; Turbulence; Turbulent flow; Velocity, Direct simulation Monte Carlo method; DSMC simulation; Laminar turbulent transitions; Lid-driven cavities; No-slip boundary conditions; Transition Reynolds number; Transition to turbulence; Turbulent velocity fluctuations, Mach number
Department/Centre:	Division of Mechanical Sciences > Chemical Engineering
Date Deposited:	02 Sep 2020 07:24
Last Modified:	02 Sep 2020 07:24
URI:	http://eprints.iisc.ac.in/id/eprint/64887

Actions (login required)

View Item