Flow Driven by a Finite, Shrouded Spinning Disk With Suction

1985 ◽  
Vol 52 (4) ◽  
pp. 766-770 ◽  
Author(s):  
J. M. Hyun

Numerical solutions are presented for the flow driven by a spinning disk which forms an endwall of a finite, closed cylinder. The effects of imposing a uniform suction (or blowing) through the spinning disk in finite configuration are investigated. The Reynolds number is large and the cylinder aspect ratio is 0(1). Finite-difference techniques are employed to integrate the time-dependent Navier-Stokes equations. The initial state is taken to be a uniform axial motion. Integration is performed until an approximate steady state is attained. When there is no suction, the infinite disk model is shown to provide a qualitatively representative approximation to the flow in the central core region. As a suction (blowing) is imposed, the core rotation rate in the case of finite configuration becomes smaller (larger) than that for the case of no suction, which is in disagreement with the predictions of the infinite disk model. These significant discrepancies point to a fundamental difficulty of the infinite disk model to adequately describe the real flow infinite geometry when there is a mass flux across the system boundary. Plots showing the meridional stream function at various times are constructed. Details of the flow structure in the approximate steady state are analyzed. When there is a suction, a strong Ekman layer is present on the spinning disk but the Ekman layer on the stationary disk fades. When there is a blowing, a strong Ekman layer forms on the stationary disk. It is shown that the dynamic effects influencing the character of the flow are confined to these Ekman layers.

1992 ◽  
Vol 237 ◽  
pp. 413-434 ◽  
Author(s):  
Jae Min Hyun ◽  
Jun Sang Park

Spin-up flows of a compressible gas in a finite, closed cylinder from an initial state of rest are studied, The flow is characterized by small reference Ekman numbers, and the peripheral Mach number is O(1). Comprehensive numerical solutions have been obtained for the full, time-dependent compressible Navier-Stokes equations. The details of the flow, temperature, and density evolution are described. In the early phase of spin-up, owing to the thermoacoustic disturbances caused by the compressible Rayleigh effect, the flows are oscillatory, and this oscillatory behaviour is pronounced at higher Mach numbers. The principal dynamical role of the Ekman layer is dominant over moderate times of orders of the homogeneous spin-up timescales. Owing to the density stratification in the radial direction, the Ekman layer is thicker in the central region of the interior. The interior azimuthal flows are mainly uniform in the axial direction. As the Mach number increases, the rate of spin-up in the interior becomes slower, and the propagating shear front is more diffusive. Explicit comparisons with the results for an infinite cylinder are made to ascertain the contributions of the endwall disks. In contrast to the usual incompressible spin-up from rest, the viscous effects are relatively more important for the case of a compressible fluid.


1982 ◽  
Vol 104 (2) ◽  
pp. 191-197 ◽  
Author(s):  
T.-W. Kuo ◽  
F. V. Bracco

A scaling law for transient, turbulent, incompressible, round jets is reported. Numerical solutions of the Navier-Stokes equations were obtained using a k-ε model for turbulence. The constants of the k-ε model were optimized by comparing computed centerline velocity, mean radial velocity distribution, longitudinal kinetic energy distributions with those measured by other authors in steady round jets. The resulting constants are those also used in computations of steady planar jets except for the one that multiplies the source term in the ε-equation. After optimization, the agreement is satisfactory for all mean quantities but is still rather poor for the kinetic energy distribution. Parameteric studies of the transient were performed for 9•103 ≤ ReD ≤ 105. Then the definition was adopted that a jet reaches steady state between the nozzle and an axial location when, at that location, the centerline velocity achieves 70 percent of its steady state value, and characteristic steadying length and time scales (D•ReD0.053 and D•ReD0.053/u cL,0 respectively) were determined as well as a unique function that relates dimensionless steadying time to dimensionless steadying length. This function changes in a predictable way if a percent other than 70 is selected but the characteristic length and time scales do not. It is found that the 70 percent threshold is reached within the head vortex of the transient jet. Thus a transient jet, practically, is a steady jet except within its head vortex. This, in part, justifies our use of steady state k-ε constants in our transient computations. The computed jet tip arrival times are shown to compare favorably with measured ones.


Open Physics ◽  
2011 ◽  
Vol 9 (1) ◽  
Author(s):  
Erik Sweet ◽  
Kuppalapalle Vajravelu ◽  
Robert Gorder

AbstractIn this paper we investigate the three-dimensional magnetohydrodynamic (MHD) rotating flow of a viscous fluid over a rotating sphere near the equator. The Navier-Stokes equations in spherical polar coordinates are reduced to a coupled system of nonlinear partial differential equations. Self-similar solutions are obtained for the steady state system, resulting from a coupled system of nonlinear ordinary differential equations. Analytical solutions are obtained and are used to study the effects of the magnetic field and the suction/injection parameter on the flow characteristics. The analytical solutions agree well with the numerical solutions of Chamkha et al. [31]. Moreover, the obtained analytical solutions for the steady state are used to obtain the unsteady state results. Furthermore, for various values of the temporal variable, we obtain analytical solutions for the flow field and present through figures.


1985 ◽  
Vol 107 (1) ◽  
pp. 92-96 ◽  
Author(s):  
Jae Min Hyun

Spin-up from rest in a cylinder with top and bottom endwall disks rotating in opposite directions (ΩT and ΩB are the respective rotation rate, but S[≡ ΩT/ΩB] < 0) is investigated. The sidewall is fixed to the faster-rotating disk. A finite-difference numerical model is adopted to integrate the unsteady Navier-Stokes equations. We consider a cylinder of aspect ratio 0(1) and minute Ekman numbers. Numerical solutions are presented to show the transient azimuthal flow structures, axial vorticity profiles, and meridional flow patterns. An azimuthal velocity front, which separates the rotating from the nonrotating fluid, propagates radially inward from the sidewall. The appearance of the front is similar to the front propagation in spin-up in a rigid cylinder. As S decreases from zero, the direction of rotation in the bulk of the interior fluid becomes the same as that of the faster-rotating disk. The azimuthal velocities are still vertically uniform in the bulk of the interior. The scaled time to reach the steady state decreases. The angular velocities of the interior fluid near the central axis become very small. Under counter-rotation, the meridional circulation forms a two-cell structure. A stagnation point appears on the slower-rotating disk. During spin-up, the stagnation point moves from the sidewall to its steady-state position. As counter-rotation increases, the radial distance traveled by the stagnation point decreases.


1973 ◽  
Vol 59 (2) ◽  
pp. 337-368 ◽  
Author(s):  
J. S. Allen

Numerical solutions of finite-difference approximations to the Navier–Stokes equations have been obtained for the axisymmetric motion of a Boussinesq liquid in a rigidly bounded rotating annulus. For most of the cases studied, a temperature difference is maintained between the top and bottom surfaces such that essentially a basic stable density stratification is imposed on the fluid. The side walls are thermally insulated and the motion is driven by a differential rotation of the top surface. Approximate steady-state solutions are obtained for various values of the Rossby number ε and the stratification parameter S = N2/Ω2, where N is the Brunt–Väisälä frequency and Ω is the rotational frequency. The changes in the flow field with the variation of these parameters is studied. Particular attention is given to an investigation of the meridional, or up welling, circulation and its dependence on the stratification parameter. The effects on the flow of different boundary conditions, such as an applied stress driving, specified temperature at the side walls and an applied heat flux at the top and bottom surfaces, are also investigated.


2002 ◽  
Vol 454 ◽  
pp. 345-364 ◽  
Author(s):  
E. MINKOV ◽  
M. UNGARISH ◽  
M. ISRAELI

Numerical finite-difference results from the full axisymmetric incompressible Navier–Stokes equations are presented for the problem of the slow axial motion of a disk particle in an incompressible, rotating fluid in a long cylindrical container. The governing parameters are the Ekman number, E = ν*/(Ω*a*2), Rossby number, Ro = W*/(Ω*a*), and the dimensionless height of the container, 2H (the scaling length is the radius of the particle, a*; Ω* is the container angular velocity, W* is the particle axial velocity and ν* the kinematic viscosity). The study concerns the flow field for small values of E and Ro while HE is of order unity, and hence the appearance of a free Taylor column (slug) of fluid ‘trapped’ at the particle is expected. The numerical results are compared with predictions of previous analytical approximate studies. First, developed (quasi-steady-state) cases are considered. Excellent agreement with the exact linear (Ro = 0) solution of Ungarish & Vedensky (1995) is obtained when the computational Ro = 10−4. Next, the time-development for both an impulsive start and a start under a constant axial force is considered. A novel unexpected behaviour has been detected: the flow field first attains and maintains for a while the steady-state values of the unbounded configuration, and only afterwards adjusts to the bounded container steady state. Finally, the effects of the nonlinear momentum advection terms are investigated. It is shown that when Ro increases then the dimensionless drag (scaled by μ*a*W*) decreases, and the Taylor column becomes shorter, this effect being more pronounced in the rear region (μ* is the dynamic viscosity). The present results strengthen and extend the validity of the classical drag force predictions and therefore the issue of the large discrepancy between theory and experiments (Maxworthy 1970) concerning this force becomes more acute.


1985 ◽  
Vol 150 ◽  
pp. 121-138 ◽  
Author(s):  
Harold R. Vaughn ◽  
William L. Oberkampf ◽  
Walter P. Wolfe

The incompressible three-dimensional Navier–Stokes equations are solved numerically for a fluid-filled cylindrical cannister that is spinning and nutating. The motion of the cannister is characteristic of that experienced by spin-stabilized artillery projectiles. Equations for the internal fluid motion are derived in a non-inertial aeroballistic coordinate system. Steady-state numerical solutions are obtained by an iterative finite-difference procedure. Flow fields and liquid induced moments have been calculated for viscosities in the range of 0.9 × 104−1 × 109 cSt. The nature of the three-dimensional fluid motion inside the cylinder is discussed, and the moments generated by the fluid are explained. The calculated moments generally agree with experimental measurements.


2004 ◽  
Vol 20 (1) ◽  
pp. 13-25 ◽  
Author(s):  
Uzu- Kuei Hsu ◽  
Chang- Hsien Tai ◽  
Chien- Hsiung Tsai

ABSTRACTThe improved numerical approach is implemented with preconditioned Navier-Stokes solver on arbitrary three-dimensional (3-D) structured multi-block complex flowfield. With the successful application of time-derivative preconditioning, present hybrid finite volume solver is performed to obtain the steady state solutions in compressible and incompressible flows. This solver which combined the adjective upwind splitting method (AUSM) family of low-diffusion flux-splitting scheme with an optimally smoothing multistage scheme and the time-derivative preconditioning is used to solve both the compressible and incompressible Euler and Navier-Stokes equations. In addition, a smoothing procedure is used to provide a mechanism for controlling the numerical implementation to avoid the instability at stagnation and sonic region. The effects of preconditioning on accuracy and convergence to the steady state of the numerical solutions are presented. There are two validation cases and three complex cases simulated as shown in this study. The numerical results obtained for inviscid and viscous two-dimensional flows over a NACA0012 airfoil at free stream Mach number ranging from 0.1 to 1.0E-7 indicates that efficient computations of flows with very low Mach numbers are now possible, without losing accuracy. And it is effectively to simulate 3-D complex flow phenomenon from compressible flow to incompressible by using the advanced numerical methods.


2018 ◽  
Vol 39 (4) ◽  
pp. 045709 ◽  
Author(s):  
Mário S M N F Gomes ◽  
Pablo Martín-Ramos ◽  
Pedro S Pereira da Silva ◽  
Manuela Ramos Silva

Sign in / Sign up

Export Citation Format

Share Document