Drag Coefficients of Viscous Spheres at Intermediate and High Reynolds Numbers

2001 ◽  
Vol 123 (4) ◽  
pp. 841-849 ◽  
Author(s):  
Zhi-Gang Feng ◽  
Efstathios E. Michaelides

A finite-difference scheme is used to solve the Navier-Stokes equations for the steady flow inside and outside viscous spheres in a fluid of different properties. Hence, the hydrodynamic force and the steady-state drag coefficient of the spheres are obtained. The Reynolds numbers of the computations range between 0.5 and 1000 and the viscosity ratio ranges between 0 (inviscid bubble) and infinity (solid particle). Unlike the numerical schemes previously implemented in similar studies (uniform grid in a stretched coordinate system) the present method introduces a two-layer concept for the computational domain outside the sphere. The first layer is a very thin one [ORe−1/2] and is positioned at the interface of the sphere. The second layer is based on an exponential function and covers the rest of the domain. The need for such a double-layered domain arises from the observation that at intermediate and large Reynolds numbers a very thin boundary layer appears at the fluid-fluid interface. The computations yield the friction and the form drag of the sphere. It is found that with the present scheme, one is able to obtain results for the drag coefficient up to 1000 with relatively low computational power. It is also observed that both the Reynolds number and the viscosity ratio play a major role on the value of the hydrodynamic force and the drag coefficient. The results show that, if all other conditions are the same, there is a negligible effect of the density ratio on the drag coefficient of viscous spheres.

2000 ◽  
Author(s):  
Zhigang Feng ◽  
Efstathios E. Michaelides

Abstract We have used a finite-difference scheme to solve the Navier-Stokes equations for the steady flow inside and outside viscous spheres in a fluid of different properties. Hence, we obtained the hydrodynamic forces and the steady-state drag coefficient of the spheres. The computational technique we have used enables us to extend the results to Reynolds numbers between 0 and 1,000. The viscosity ratio of the computations ranges between 0 (inviscid bubble) and infinity (solid particle). The method presented here makes use of a two-layer concept for the computational domain outside the sphere. The first layer is a very thin one [O(Re−1/2)] and is positioned at the interface of the sphere. The second layer is based on an exponential function and covers the rest of the domain. The computations yield the friction and the form drag of the sphere. It is observed that both the Reynolds number and the viscosity ratio play a major role on the value of the hydrodynamic force and the drag coefficient. If all other conditions are the same, there is a negligible effect of the density ratio on the drag coefficient of viscous spheres.


2021 ◽  
Vol 61 (SI) ◽  
pp. 155-162
Author(s):  
Petr Sváček

This paper is interested in the mathematical modelling of the voice production process. The main attention is on the possible closure of the glottis, which is included in the model with the concept of a fictitious porous media and using the Hertz impact force The time dependent computational domain is treated with the aid of the Arbitrary Lagrangian-Eulerian method and the fluid motion is described by the incompressible Navier-Stokes equations coupled to structural dynamics. In order to overcome the instability caused by the dominating convection due to high Reynolds numbers, stabilization procedures are applied and numerically analyzed for a simplified problem. The possible distortion of the computational mesh is considered. Numerical results are shown.


Author(s):  
Carl E. Rathmann

For well over 150 years now, theoreticians and practitioners have been developing and teaching students easily visualized models of fluid behavior that distinguish between the laminar and turbulent fluid regimes. Because of an emphasis on applications, perhaps insufficient attention has been paid to actually understanding the mechanisms by which fluids transition between these regimes. Summarized in this paper is the product of four decades of research into the sources of these mechanisms, at least one of which is a direct consequence of the non-linear terms of the Navier-Stokes equation. A scheme utilizing chaotic dynamic effects that become dominant only for sufficiently high Reynolds numbers is explored. This paper is designed to be of interest to faculty in the engineering, chemistry, physics, biology and mathematics disciplines as well as to practitioners in these and related applications.


2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.


1969 ◽  
Vol 11 (5) ◽  
pp. 445-453 ◽  
Author(s):  
A. K. Runchal ◽  
M. Wolfshtein

A procedure is proposed for the integration of the full Navier-Stokes equations for constant-property two-dimensional flows. In contrast with earlier procedures, the present one is capable of dealing with cases of very high Reynolds number. The power of the new procedure is demonstrated in two cases: (1) the square recirculating eddy, for which no solution was previously available for Reynolds numbers larger than about 400, and (2) an impinging jet, for which no solution was available previously. The procedure has also been applied to the square cavity at Reynolds numbers below 400; it gives results of an accuracy comparable with that of previous solutions, but with a smaller computing time.


Sign in / Sign up

Export Citation Format

Share Document