Flow Between Eccentric Rotating Cylinders

1972 ◽  
Vol 94 (3) ◽  
pp. 266-274 ◽  
Author(s):  
R. C. DiPrima ◽  
J. T. Stuart

The flow between eccentric rotating cylinders is studied using a modified bipolar coordinate system. The Sommerfeld pressure distribution and associated flow in a journal bearing are derived from the Navier-Stokes equations by a straightforward and systematic expansion in the clearance ratio (assumed small). First order corrections for curvature and inertial effects are given. There is no restriction on the eccentricity.

2016 ◽  
pp. 90-92
Author(s):  
A. G. Obukhov ◽  
R. E. Volkov

It is proved that complex flows of the viscous compressible heat-conducting gas, arising during heating the vertical field, have a pronounced axial symmetry. Therefore, for the numerical solution of the full Navier-Stokes equations for description of such gas flows it are advisable to use a cylindrical coordinate system. This paper describes the transformation of the first projection of the equation of motion of the full Navier-Stokes equations system. The result of the transformation is a record of the first projection of the equation of a continuous medium motion in the cylindrical coordinate system.


2006 ◽  
Vol 128 (6) ◽  
pp. 1394-1399 ◽  
Author(s):  
Donghyun You ◽  
Meng Wang ◽  
Rajat Mittal ◽  
Parviz Moin

A novel structured grid approach which provides an efficient way of treating a class of complex geometries is proposed. The incompressible Navier-Stokes equations are formulated in a two-dimensional, generalized curvilinear coordinate system complemented by a third quasi-curvilinear coordinate. By keeping all two-dimensional planes defined by constant third coordinate values parallel to one another, the proposed approach significantly reduces the memory requirement in fully three-dimensional geometries, and makes the computation more cost effective. The formulation can be easily adapted to an existing flow solver based on a two-dimensional generalized coordinate system coupled with a Cartesian third direction, with only a small increase in computational cost. The feasibility and efficiency of the present method have been assessed in a simulation of flow over a tapered cylinder.


2016 ◽  
Vol 798 ◽  
pp. 643-681 ◽  
Author(s):  
E. D. Fernández-Nieto ◽  
J. Garres-Díaz ◽  
A. Mangeney ◽  
G. Narbona-Reina

In this work we present a multilayer shallow model to approximate the Navier–Stokes equations with the ${\it\mu}(I)$-rheology through an asymptotic analysis. The main advantages of this approximation are (i) the low cost associated with the numerical treatment of the free surface of the modelled flows, (ii) the exact conservation of mass and (iii) the ability to compute two-dimensional profiles of the velocities in the directions along and normal to the slope. The derivation of the model follows Fernández-Nieto et al. (J. Comput. Phys., vol. 60, 2014, pp. 408–437) and introduces a dimensional analysis based on the shallow flow hypothesis. The proposed first-order multilayer model fully satisfies a dissipative energy equation. A comparison with steady uniform Bagnold flow – with and without the sidewall friction effect – and laboratory experiments with a non-constant normal profile of the downslope velocity demonstrates the accuracy of the numerical model. Finally, by comparing the numerical results with experimental data on granular collapses, we show that the proposed multilayer model with the ${\it\mu}(I)$-rheology qualitatively reproduces the effect of the erodible bed on granular flow dynamics and deposits, such as the increase of runout distance with increasing thickness of the erodible bed. We show that the use of a constant friction coefficient in the multilayer model leads to the opposite behaviour. This multilayer model captures the strong change in shape of the velocity profile (from S-shaped to Bagnold-like) observed during the different phases of the highly transient flow, including the presence of static and flowing zones within the granular column.


Author(s):  
Quentin Desbonnets ◽  
Daniel Broc

It is well known that a fluid may strongly influence the dynamic behaviour of a structure. Many different physical phenomena may take place, depending on the conditions: fluid flow, fluid at rest, little or high displacements of the structure. Inertial effects can take place, with lower vibration frequencies, dissipative effects also, with damping, instabilities due to the fluid flow (Fluid Induced Vibration). In this last case the structure is excited by the fluid. Tube bundles structures are very common in the nuclear industry. The reactor cores and the steam generators are both structures immersed in a fluid which may be submitted to a seismic excitation or an impact. In this case the structure moves under an external excitation, and the movement is influence by the fluid. The main point in such system is that the geometry is complex, and could lead to very huge sizes for a numerical analysis. Homogenization models have been developed based on the Euler equations for the fluid. Only inertial effects are taken into account. A next step in the modelling is to build models based on the homogenization of the Navier-Stokes equations. The papers presents results on an important step in the development of such model: the analysis of the fluid flow in a oscillating tube bundle. The analysis are made from the results of simulations based on the Navier-Stokes equations for the fluid. Comparisons are made with the case of the oscillations of a single tube, for which a lot of results are available in the literature. Different fluid flow pattern may be found, depending in the Reynolds number (related to the velocity of the bundle) and the Keulegan-Carpenter number (related to the displacement of the bundle). A special attention is paid to the quantification of the inertial and dissipative effects, and to the forces exchanges between the bundle and the fluid. The results of such analysis will be used in the building of models based on the homogenization of the Navier-Stokes equations for the fluid.


2002 ◽  
Vol 465 ◽  
pp. 213-235 ◽  
Author(s):  
D. R. GRAHAM ◽  
J. J. L. HIGDON

Oscillatory forcing of a porous medium may have a dramatic effect on the mean flow rate produced by a steady applied pressure gradient. The oscillatory forcing may excite nonlinear inertial effects leading to either enhancement or retardation of the mean flow. Here, in Part 1, we consider the effects of non-zero inertial forces on steady flows in porous media, and investigate the changes in the flow character arising from changes in both the strength of the inertial terms and the geometry of the medium. The steady-state Navier–Stokes equations are solved via a Galerkin finite element method to determine the velocity fields for simple two-dimensional models of porous media. Two geometric models are considered based on constricted channels and periodic arrays of circular cylinders. For both geometries, we observe solution multiplicity yielding both symmetric and asymmetric flow patterns. For the cylinder arrays, we demonstrate that inertial effects lead to anisotropy in the effective permeability, with the direction of minimum resistance dependent on the solid volume fraction. We identify nonlinear flow phenomena which might be exploited by oscillatory forcing to yield a net increase in the mean flow rate. In Part 2, we take up the subject of unsteady flows governed by the full time-dependent Navier–Stokes equations.


Author(s):  
Changhu Xing ◽  
Minel J. Braun

Dynamic coefficients are very important for the stability of a hydrodynamic journal bearing and therefore for its design. In order to determine the stiffness, damping and added mass coefficients of the hydrodynamic bearing, the finite perturbation method around its stabilization position was employed. Based on the Reynolds equation with Gumbel cavitation algorithm, the maximum magnitude of the perturbation was judged by comparing results from finite perturbation (numerical way) to those from infinitesimal perturbation (additional analytical equations need to be derived based on order analysis), as well as theoretical analysis. Using the determined perturbation amplitude, the full three-dimensional Navier-Stokes equations in CFD-ACE+ were used to evaluate coefficients from an actual lubricant and compare to those obtained with Reynolds equation. Finally, a homogeneous gaseous cavitation algorithm is coupled with the Navier-Stokes equation to establish the pressure distribution in the bearing. When gas concentration was varied, the pressure distribution as well as the dynamic coefficients changed significantly.


2006 ◽  
Vol 129 (2) ◽  
pp. 221-230 ◽  
Author(s):  
Alex de Kraker ◽  
Ron A. J. van Ostayen ◽  
A. van Beek ◽  
Daniel J. Rixen

In this paper a multiscale method is presented that includes surface texture in a mixed lubrication journal bearing model. Recent publications have shown that the pressure generating effect of surface texture in bearings that operate in full film conditions may be the result of micro-cavitation and/or convective inertia. To include inertia effects, the Navier–Stokes equations have to be used instead of the Reynolds equation. It has been shown in earlier work (de Kraker et al., 2006, Tribol. Trans., in press) that the coupled two-dimensional (2D) Reynolds and 3D structure deformation problem with partial contact resulting from the soft EHL journal bearing model is not easy to solve due to the strong nonlinear coupling, especially for soft surfaces. Therefore, replacing the 2D Reynolds equation by the 3D Navier–Stokes equations in this coupled problem will need an enormous amount of computing power that is not readily available nowadays. In this paper, the development of a micro–macro multiscale method is described. The local (micro) flow effects for a single surface pocket are analyzed using the Navier–Stokes equations and compared to the Reynolds solution for a similar smooth piece of surface. It is shown how flow factors can be derived and added to the macroscopic smooth flow problem, that is modeled by the 2D Reynolds equation. The flow factors are a function of the operating conditions such as the ratio between the film height and the pocket dimensions, the surface velocity, and the pressure gradient over a surface texture unit cell. To account for an additional pressure buildup in the texture cell due to inertia effects, a pressure gain is introduced at macroscopic level. The method also allows for microcavitation. Microcavitation occurs when the pressure variation due to surface texture is larger than the average pressure level at that particular bearing location. In contrast with the work of Patir and Cheng (1978, J. Lubrication Technol., 78, pp. 1–10), where the microlevel is solved by the Reynolds equation, and the Navier–Stokes equations are used at the microlevel. Depending on the texture geometry and film height, the Reynolds equation may become invalid. A second pocket effect occurs when the pocket is located in the moving surface. In mixed lubrication, fluid can become trapped inside a pocket and squeezed out when the pocket is running into an area with higher contact load. To include this effect, an additional source term that represents the average fluid inflow due to the deformation of the surface around the pocket is added to the Reynolds equation at macrolevel. The additional inflow is computed at microlevel by numerical solution of the surface deformation for a single pocket that is subject to a contact load. The pocket volume is a function of the contact pressure. It must be emphasized that before ready-to-use results can be presented, a large number of simulations to determine the flow factors and pressure gain as a function of the texture parameters and operating conditions have yet to be done. Before conclusions can be drawn, regarding the dominanant mechanism(s), the flow factors and pressure gain have to be added to the macrobearing model. In this paper, only a limited number of preliminary illustrative simulation results, calculating the flow factors for a single 2D texture geometry, are shown to give insight into the method.


1989 ◽  
Vol 111 (3) ◽  
pp. 545-552 ◽  
Author(s):  
R. Nordmann ◽  
F. J. Dietzen ◽  
H. P. Weiser

The compressible flow in a seal can be described by the Navier-Stokes equations in connection with a turbulence model (k–ε model) and an energy equation. By introducing a perturbation analysis in these differential equations we obtain zeroth order equations for the centered position and first order equations for small motions of the shaft about the centered position. These equations are solved by a finite difference technique. The zeroth order equations describe the leakage flow. Integrating the pressure solution of the first order equations yields the fluid forces and the rotordynamic coefficients, respectively.


2017 ◽  
Vol 818 ◽  
pp. 344-365 ◽  
Author(s):  
Dominik Dierkes ◽  
Martin Oberlack

The present contribution is a significant extension of the work by Kelbin et al. (J. Fluid Mech., vol. 721, 2013, pp. 340–366) as a new time-dependent helical coordinate system has been introduced. For this, Lie symmetry methods have been employed such that the spatial dependence of the originally three independent variables is reduced by one and the remaining variables are: the cylindrical radius $r$ and the time-dependent helical variable $\unicode[STIX]{x1D709}=(z/\unicode[STIX]{x1D6FC}(t))+b\unicode[STIX]{x1D711}$, $b=\text{const.}$ and time $t$. The variables $z$ and $\unicode[STIX]{x1D711}$ are the usual cylindrical coordinates and $\unicode[STIX]{x1D6FC}(t)$ is an arbitrary function of time $t$. Assuming $\unicode[STIX]{x1D6FC}=\text{const.}$, we retain the classical helically symmetric case. Using this, and imposing helical invariance onto the equation of motion, leads to a helically symmetric system of Euler and Navier–Stokes equations with a time-dependent pitch $\unicode[STIX]{x1D6FC}(t)$, which may be varied arbitrarily and which is explicitly contained in all of the latter equations. This has been conducted both for primitive variables as well as for the vorticity formulation. Hence a significantly extended set of helically invariant flows may be considered, which may be altered by an external time-dependent strain along the axis of the helix. Finally, we sought new conservation laws which can be found from the helically invariant Euler and Navier–Stokes equations derived herein. Most of these new conservation laws are considerable extensions of existing conservation laws for helical flows at a constant pitch. Interestingly enough, certain classical conservation laws do not admit extensions in the new time-dependent coordinate system.


Sign in / Sign up

Export Citation Format

Share Document