Relaxation processes in hypersonic rarefied binary gas mixture flow around a cylinder

Using the direct simulation Monte Carlo method, the hypersonic flow of a binary gas mixture around a cylinder is investigated over a wide rarefaction range: from an almost continuum regime (at the Knudsen number Kn = 0.01) to free-molecular flow. The effect of a small admixture of heavy diatomic particles in a light gas flow on the relaxation processes near the cylinder and the heat flux is studied.     

Investigation of slow monatomic gas flows past a circular cylinder

Slow low-Knudsen-number monatomic-gas flow past a circular cylinder is numerically investigated on the basis of a model kinetic equation. The gas flow is described by a new kinetic equation, from which the continuum equations for slow nonisothermal gas flows containing temperature stresses follow rigorously. It is shown that a closed convective-flow region arises near a nonuniformly heated cylinder in a slow gas flow if the flow impinges on the hot side of its surface. Using a new model of the Boltzmann equation makes it possible to study gas flows both in continuum and rarefied flow regimes.     

Oblique shock/bow shock interference in rarefied flow past a cylinder

Direct statistical simulation is employed to study the flow of a rarefied diatomic gas past a cylinder in the presence of an incident oblique shock. The distinctive features of the formation of a high-pressure compressed-gas jet in the case of interference between the oblique shock and the bow shock are studied for different Reynolds numbers. The variation of the pressure and the heat transfer to the surface with the shock position relative to the center of the cylinder, the Reynolds number, and the surface temperature is analyzed. The results obtained are compared with the experimental data and the results of the numerical solutions of the Euler and boundary layer equations. Free-molecular-to-continuum flow transition is demonstrated with reference to the example of interference-free flow past a cylinder.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 171–180. Original Russian Text Copyright © 2004 by Gusev and Erofeev.     

Viscous oscillatory flow around a circular cylinder at low Keulegan–Carpenter numbers and frequency parameters

The results of a numerical study of the viscous oscillating flow around a circular cylinder at low Keulegan–Carpenter numbers (KC) and frequency parameters (β) are presented in this paper. The finite element method was used for the solution of the Navier–Stokes equations in the formulation where the streamfunction and vorticity are the field variables. The computation was conducted at Keulegan–Carpenter numbers extending up to KC=15 and frequency parameters ranging between β=6 and 100. At low values of the Keulegan–Carpenter number the flow remains symmetrical. As the Keulegan–Carpenter number is increased over a certain value which depends also on the frequency parameter, asymmetries appear in the flow which are eventually amplified and lead finally to complex vortex-shedding patterns, some of which are markedly different from those observed at higher frequency parameters. The solution revealed that although for certain values of KC and β the shedding of vortices is periodic, there also exists a complicated flow regime in which the flow is not periodic but switches between different modes in consecutive cycles of flow oscillation. For the various flow cases examined, the traces of the hydrodynamic forces are presented and the hydrodynamic coefficients and RMS values of the in-line force are compared with experimental evidence. © 1998 John Wiley & Sons, Ltd.     

Supersonic transverse rarefied gas flow past a strip

The supersonic flow of a monatomic gas consisting of hard spherical particles past a flat strip normal to the flow is investigated using the direct simulation Monte-Carlo (DSMC) method. The calculations are performed over the Knudsen and Mach number ranges 0.015–5 and 1.8–15, respectively. The structure of the compressed layer and the aerodynamic characteristics are systematically studied for the Mach number 5 and various Knudsen numbers. The dependences of the compressed-layer thickness in molecular free paths are found. The nonequilibrium processes in the neighborhood of the strip are described on the basis of the data on the temperature anisotropy with respect to three coordinates.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 159–167. Original Russian Text Copyright © 2005 by Maltsev and Rebrov.     

Unsteady characteristics of the cavity flow around a rough circular cylinder

Effects of surface roughness on the unsteady cavitating flow around a two-dimensional circular cylinder were experimentally investigated at Reynolds numbers from 1.36 × 10⁵ to 1.78 × 10⁵. Two patterns of surface roughness were investigated, a double-cut pattern and a single-cut pattern. The cavity elongates with an increase of the surface roughness, especially in supercavitating flow. However, for some roughness parameters tested, the cavity length exhibits an extreme decrease. In a particular case of the double-cut pattern, there exists the minimum cavity behind the cylinder.     

Effect of the accommodation coefficient on transfer processes in a supersonic rarefied gas flow around a transversely aligned cylinder

A supersonic flow around a cylinder is studied by the direct statistical Monte Carlo method in a wide range of rarefaction: from regimes close to continuum to free-molecular flow. The effect of the accommodation coefficient on the flow near the cylinder and on heat transfer between the gas and the cylinder is examined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 64–72, May–June, 2008.     

Investigation of scale effect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

Circulating rarefied gas flow in a rotating cylinder with a stationary upper end face

In order to analyze the intensity of the circulating flow in the rarefaction zone near the axis of a rapidly rotating cylinder, the problem of gas flow in a cylinder with a stationary upper end face is solved for intermediate Knudsen numbers by the direct statistical modeling method. The effect of the rarefaction of the gas on the intensity of the secondary flow is investigated. The contribution of the self-excited thermal circulation of the gas to the total circulating flow is analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 186–188, May–June, 1990.     

Numerical investigation of steady suction control of flow around a circular cylinder

This paper numerically investigates the effectiveness of the control of steady suction on a stationary circular cylinder with several isolated suction holes on the surface at a subcritical Reynolds number. The control effectiveness as a function of the azimuthal position, spanwise spacing and suction flow rate of the suction holes on the control of the aerodynamic forces on the cylinder and the suppression of alternate vortex shedding are taken into account. The study of the azimuthal location of the suction holes indicates that azimuthal angles of θ=90° and 270°, which are close to the separation point, provide the most substantial decreases in the aerodynamic forces. When restricted to the most effective azimuthal angle, a remarkable control effectiveness can be achieved when the axial spacing between two neighboring suction holes is less than a minimum value even under a small suction momentum coefficient. However, if the axial spacing exceeds the minimum spacing, the control effectiveness will not be saturated even under a very large suction momentum coefficient. Thus, the cause of the effective aerodynamic force control is suggested to be a result of obvious three-dimensional phenomenon in the near wake, which is characterized by the generation of a convergent flow between two neighboring suction hole sections and a stronger, larger three-dimensional vortex pair adjacent to the convergent flow. It has been suggested that this strongly three-dimensional flow pattern is induced by the strong interaction between two neighboring but counter-rotating three-dimensional vortices separately produced by two neighboring suction holes. Moreover, the effects of such three-dimensional flow patterns are investigated in detail based on variations in the flow field and sectional aerodynamic forces in different cross sections. Finally, the upper limit of the axial spacing between two neighboring suction holes to form such a three-dimensional flow pattern is suggested to be between 0.75 D and 1.5 D when the suction flow rate exceeds a certain value.     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

Buffeting lift forces and local air–water flow aspects around a rigid cylinder

A new experimental programme is conducted in order to relate the characteristics of two-phase flow around a rigid cylinder with the resulting lift forces. The local characteristics of air–water flow measured in the vicinity of the cylinder provide a useful source of information about the effects of flow on the excitation mechanisms. In particular, a selection of relevant parameters has been identified which, with the help of a standard dimensional analysis, may explain the energetic contents of buffeting forces. Among the parameters effective in reducing the data are the flow regime, bubble frequency and gravity forces. In addition, in the range of bubbly regimes, the magnitude of the random forces is found to be related to the local fluctuations of void fraction. Finally, a new formulation is proposed to collapse the dimensionless spectra of the buffeting lift forces in a single characteristic curve. This analysis shows a marked improvement over the collapse of data in comparison with previous normalized models.     

Very slow flow of Bingham viscoplastic fluid around a circular cylinder

Numerical simulations have been used to study the flow of a Bingham viscoplastic fluid around a circular cylinder in an infinite medium with negligible inertia effects. Papanastasiou's regularisation technique has been adopted to approximate the model. The case corresponding to preponderant plasticity effects has been particularly studied and convergence of the solutions examined in detail. The flow kinematics and stresses have been determined. The rigid zones have been identified and characterised. At large Oldroyd numbers, when plasticity effects become preponderant, a viscoplastic boundary layer appears around the cylinder. The characteristics of this viscoplastic boundary layer are quantified. The results are compared with existing theoretical results, concerning particularly the predictions of the viscoplastic boundary layer theory and the plasticity theory.     

Asymptotic solutions of the two-dimensional equations for a thin viscous shock layer in a rarefied gas

Two-dimensional hypersonic rarefied gas flow around blunt bodies is investigated for the continuum to free-molecular transition regime. In [1], as a result of an asymptotic analysis, three rarefied gas flow regimes, depending on the relationship between the problem parameters, were detected and one of these regimes was investigated. In the present study, asymptotic solutions of the thin viscous shock layer equations at small Reynolds numbers are obtained for the other two flow regimes. Analytical expressions for the heat transfer, friction and pressure coefficients are obtained as functions of the incident flow parameters and the body geometry and temperature. As the Reynolds number tends to zero, the values of these coefficients approach their values in free-molecular flow. The scaling parameters of hypersonic rarefied gas flow around bodies are determined for different regimes. The asymptotic solutions are compared with the results of direct Monte Carlo simulation.     

Numerical solutions for impulsively started and decelerated viscous flow past a circular cylinder

A finite difference study of the unsteady two-dimensional flow past a circular cylinder has been conducted using vorticity and streamfunction as the dependent variables. The two cases considered were impulsively started and decelerated flows. The impulsively started problem was considered to validate the method and has yielded results which agree quite closely with existing results from both calculations and experiments. The decelerated flow analysis produced results which can be explained in terms of induced velocity effects from existing wake vortices for both suddenly stopped and uniformly decelerated flows.     

Simulation of the rarefied gas flow around a perpendicular disk

The present work contributes to define the domain of validity of two continuum approaches, based on Navier–Stokes (NS) and quasi gas-dynamic (QGD) equations, respectively. Results obtained using each method are compared with those obtained using a direct simulation Monte Carlo (DSMC) method considered as a reference. QGD equations differ from NS ones by the presence of additional dissipative terms. The present paper includes a brief presentation of QGD equations and DSMC procedures used here. The rarefied flow around a perpendicular disk has been considered for a freestream Mach number varying from 2 to 20, a Knudsen number equal to 0.1 and two levels of wall temperature.     

Numerical simulation of a flow around an impulsively started radially deforming circular cylinder

The development of a two‐dimensional viscous incompressible flow generated by a deformable circular cylinder impulsively started into rectilinear motion is studied numerically for the Reynolds numbers equal to 550 and 3000. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the streamfunction is solved by a Fourier–Galerkin method in the direction of the flow that is assumed to remain symmetrical and a second‐order finite difference for the radial direction. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The computed results are compared qualitatively with experimental and numerical results done before in the particular non‐deformable case. The comparison is found to be satisfactory. The influence of the deformation of the cylinder on the flow structure and the drag coefficient is then analyzed. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical investigation from rarefied flow to continuum by solving the Boltzmann model equation

Based on the Bhatnagar–Gross–Krook (BGK) Boltzmann model equation, the unified simplified velocity distribution function equation adapted to various flow regimes can be presented. The reduced velocity distribution functions and the discrete velocity ordinate method are developed and applied to remove the velocity space dependency of the distribution function, and then the distribution function equations will be cast into hyperbolic conservation laws form with non‐linear source terms. Based on the unsteady time‐splitting technique and the non‐oscillatory, containing no free parameters, and dissipative (NND) finite‐difference method, the gas kinetic finite‐difference second‐order scheme is constructed for the computation of the discrete velocity distribution functions. The discrete velocity numerical quadrature methods are developed to evaluate the macroscopic flow parameters at each point in the physical space. As a result, a unified simplified gas kinetic algorithm for the gas dynamical problems from various flow regimes is developed. To test the reliability of the present numerical method, the one‐dimensional shock‐tube problems and the flows past two‐dimensional circular cylinder with various Knudsen numbers are simulated. The computations of the related flows indicate that both high resolution of the flow fields and good qualitative agreement with the theoretical, DSMC and experimental results can be obtained. Copyright © 2003 John Wiley & Sons, Ltd.