A Computational Fluid Dynamics Simulation of Oil–Air Flow Between the Cage and Inner Race of an Aero-engine Bearing

In aeroengines the shafts are supported on bearings that carry the radial and axial loads. A ball bearing is made up of an inner-race, an outer-race and a cage which contains the balls, these together comprise the bearing elements. The bearings require oil for lubrication and cooling. The design of the bearing studied in this work is such that the oil is fed to the bearing through holes/slots in the inner race. At each axial feed location the oil is fed through a number of equispaced feedholes/slots but there is a different number of holes at each location. Once the oil has passed through the bearing it sheds outwards from both sides into compartments known as the bearing chambers. A number of studies have been carried out on the dynamics of bearings. Most of the analyses consider the contributions of fluid forces as small relative to the interaction of the bearing elements. One of the most sophisticated models for a cage-raceway analysis is based on the work of Ashmore et al. [1], where the cage-raceway is considered to be a short journal bearing divided into sectors by the oil feeds. It is further assumed that the oil exits from the holes and forms a continuous block of oil that exits outwards on both sides of the cage-raceway. In the model, the Reynolds equation is used to estimate the oil dynamics. Of interest in this current work is the behaviour of the oil and air within the space bounded by the cage and inner race. The aim is to determine whether oil feed to the bearing can be modelled as coming from a continuous slot or if the discrete entry points must be modelled. A Volume of Fluid Computational Fluid Dynamics approach is applied. A sector of a ball bearing is modelled with a fine mesh and the detailed simulations show the flow behaviour for different oil splits to the three feed locations of the bearing thus providing information useful to understanding oil shedding into the bearing chambers. The work shows that different flow behaviour is predicted by models where the oil inlets through a continuous slot compared to discrete entry holes. The form and speed of oil shedding from the bearing is found to depend strongly on shaft speed with the shedding speed being slightly higher than the cage linear speed. The break-up pattern of oil on the cage inner surface suggests smaller droplets will be shed at higher shaft speed

Comparison of Computational Results with a Low-g, Nitrogen Slosh and Boiling Experiment

This paper compares a fluid/thermal simulation, in Fluent, with a low-g, nitrogen slosh and boiling experiment. In 2010, the French Space Agency, CNES, performed cryogenic nitrogen experiments in a low-g aircraft campaign. From one parabolic flight, a low-g interval was simulated that focuses on low-g motion of nitrogen liquid and vapor with significant condensation, evaporation, and boiling. The computational results are compared with high-speed video, pressure data, heat transfer, and temperature data from sensors on the axis of the cylindrically shaped tank. These experimental and computational results compare favorably. The initial temperature stratification is in good agreement, and the two-phase fluid motion is qualitatively captured. Temperature data is matched except that the temperature sensors are unable to capture fast temperature transients when the sensors move from wet to dry (liquid to vapor) operation. Pressure evolution is approximately captured, but condensation and evaporation rate modeling and prediction need further theoretical analysis

Jet pumps for thermoacoustic applications: design guidelines based on a numerical parameter study

The oscillatory flow through tapered cylindrical tube sections (jet pumps) is characterized by a numerical parameter study. The shape of a jet pump results in asymmetric hydrodynamic end effects which cause a time-averaged pressure drop to occur under oscillatory flow conditions. Hence, jet pumps are used as streaming suppressors in closed-loop thermoacoustic devices. A two-dimensional axisymmetric computational fluid dynamics model is used to calculate the performance of a large number of conical jet pump geometries in terms of time-averaged pressure drop and acoustic power dissipation. The investigated geometrical parameters include the jet pump length, taper angle, waist diameter and waist curvature. In correspondence with previous work, four flow regimes are observed which characterize the jet pump performance and dimensionless parameters are introduced to scale the performance of the various jet pump geometries. The simulation results are compared to an existing quasi-steady theory and it is shown that this theory is only applicable in a small operation region. Based on the scaling parameters, an optimum operation region is defined and design guidelines are proposed which can be directly used for future jet pump design.Comment: The following article has been accepted by the Journal of the Acoustical Society of America. After it is published, it will be found at http://scitation.aip.org/JAS

A numerical investigation on the vortex formation and flow separation of the oscillatory flow in jet pumps

A two-dimensional computational fluid dynamics model is used to predict the oscillatory flow through a tapered cylindrical tube section (jet pump) placed in a larger outer tube. Due to the shape of the jet pump, there will exist an asymmetry in the hydrodynamic end effects which will cause a time-averaged pressure drop to occur that can be used to cancel Gedeon streaming in a closed-loop thermoacoustic device. The performance of two jet pump geometries with different taper angles is investigated. A specific time-domain impedance boundary condition is implemented in order to simulate traveling acoustic wave conditions. It is shown that by scaling the acoustic displacement amplitude to the jet pump dimensions, similar minor losses are observed independent of the jet pump geometry. Four different flow regimes are distinguished and the observed flow phenomena are related to the jet pump performance. The simulated jet pump performance is compared to an existing quasi-steady approximation which is shown to only be valid for small displacement amplitudes compared to the jet pump length.Comment: The following article has been accepted by the Journal of the Acoustical Society of America. After it is published, it will be found at: http://scitation.aip.org/JAS

A Reduced Nonlinear Model for the Simulation of Two Phase Flow in a Horizontal Pipe.

In the last 10 years many 3D numerical schemes have been developed for the study the flow of a mixture of liquid and gas in a pipeline (Frank, Numerical simulation of slug flow regime for an air-water two-phase flow in horizontal pipes. In: The 11th international topical meeting on nuclear reactor thermal-hydraulics (NURETH-11), Avignon, 2005; Vallée et al., Nucl Eng Des 238(3):637–646, 2008; Höhne, Experiments and numerical simulations of horizontal two-phase flow regimes. In: Proceeding of the seventh international conference on CFD in the minerals and process industries, Melbourne, 2009; Bartosiewicz et al., Nucl Eng Des 240(9):2375–2381, 2010) but although they offer a very good accuracy, they are rarely fit for modelling a long pipe, due to the high computational costs. Then one is usually led to consider 1D models, see e.g. the works of Issa and his group (Issa and Kempf, Int J Multiphase Flow 29(1):69–95, 2003). Such models offer much faster simulations than 3D schemes, on the other hand they almost completely miss the dynamics in the transversal direction. Here we present a model able of representing the full 3D dynamics, but with the computational cost typical of 1D simulation. The main feature of our model consists in describing the dynamical variables in the direction transversal to the pipe by means of a family of functions depending on a set of parameters. The model is then solved by a standard finite volume scheme

Numerical modeling of oscillating Taylor bubbles

In this study, computational fluid dynamics (CFD) modeling is used to simulate Taylor bubbles rising in vertical pipes. Experiments indicate that in large diameter (0.29 m) pipes for an air–water system, the bubbles can rise in a oscillatory manner, depending on the method of air injection. The CFD models are able to capture this oscillatory behavior because the air phase is modeled as a compressible ideal gas. Insights into the flow field ahead and behind the bubble during contraction and expansion are shown. For a bubble with an initial pressure equal to the hydrostatic pressure at its nose, no oscillations are seen in the bubble as it rises. If the initial pressure in the bubble is set less than or greater than the hydrostatic pressure then the length of the bubble oscillates with an amplitude that depends on the magnitude of the initial bubble pressure relative to the hydrostatic pressure. The frequency of the oscillations is inversely proportional to the square root of the head of water above the bubble and so the frequency increases as the bubble approaches the water surface. The predicted frequency also depends inversely on the square root of the average bubble length, in agreement with experimental observations and an analytical model that is also presented. In this model, a viscous damping term due to the presence of a Stokes boundary layer for the oscillating cases is introduced for the first time and used to assess the effect on the oscillations of increasing the liquid viscosity by several orders of magnitude

An effective mesh strategy for CFD modelling of polymer electrolyte membrane fuel cells

Computational fluid dynamics (CFD) is a major tool in PEM fuel cell research. Typical three-dimensional PEM fuel cell models involve more than 106 mesh elements. This makes the computation very intense and necessitates a methodology to mesh the computational domain that can keep the number of elements to a minimum while maintaining good accuracy. In this study, the effect of computational mesh in each direction on the accuracy of the solution is investigated in a systematic way. It is found that the mesh in different directions has a different degree of influence on the solution suggesting that the mesh in one direction can be coarser than the other. The proposed mesh strategy is capable of greatly reducing the number of mesh elements, hence computation time, while preserving the characteristics of important flow-field variables. Moreover, it is applicable to a wide range of cell sizes and flow-field configurations and should be used as a guideline for mesh generation

An analytical and experimental assessment of flexible road ironwork support structures

This paper describes work undertaken to investigate the mechanical performance of road ironwork installations in highways, concentrating on the chamber construction. The principal aim was to provide the background research which would allow improved designs to be developed to reduce the incidence of failures through improvements to the structural continuity between the installation and the surrounding pavement. In doing this, recycled polymeric construction materials (Jig Brix) were studied with a view to including them in future designs and specifications. This paper concentrates on the Finite Element (FE) analysis of traditional (masonry) and flexible road ironwork structures incorporating Jig Brix. The global and local buckling capacity of the Jig Brix elements was investigated and results compared well with laboratory measurements. FE models have also been developed for full-scale traditional (masonry) and flexible installations in a surrounding flexible (asphalt) pavement structure. Predictions of response to wheel loading were compared with full-scale laboratory measurements. Good agreement was achieved with the traditional (masonry) construction but poorer agreement for the flexible construction. Predictions from the FE model indicated that the use of flexible elements significantly reduces the tensile horizontal strain on the surface of the surrounding asphaltic material which is likely to reduce the incidence of surface cracking

Prediction of cavitation and induced erosion inside a high-pressure fuel pump

The operation of a high-pressure, piston-plunger fuel pump, oriented for use in the common rail circuit of modern Diesel engines for providing fuel to the injectors is investigated in the present study from a numerical perspective. Both the suction and pressurization phases of the pump stroke were simulated with the overall flow time be-ing in the order of 12•10-3 s. The topology of the cavitating flow within the pump con-figuration was captured through the use of an Equation of State (EoS) implemented in the framework of a barotropic, homogeneous equilibrium model. Cavitation was found to set in within the pressure chamber as early as 0.2•10-3 s in the operating cycle, while the minimum liquid volume fraction detected was in the order of 60% during the sec-ond period of the valve opening. Increase of the in-cylinder pressure during the final stages of the pumping stroke lead to the collapse of the previously arisen cavitation structures and three layout locations, namely the piston edge, the valve/valve-seat re-gion and the outlet orifice, were identified as vulnerable to cavitation-induced erosion through the use of cavitation-aggressiveness indicators

A hybrid actuator disc - full rotor CFD methodology for modelling the effects of wind turbine wake interactions on performance

The performance of individual wind turbines is crucial for maximum energy yield, however, their performance is often reduced when turbines are placed together in an array. The wake produced by the rotors interacts with downstream turbines, resulting in a reduction in power output. In this paper, we demonstrate a new and faster modelling technique which combines actuator disc theory, modelled using wind tunnel validated Computational Fluid Dynamics (CFD), and integrated into full rotor CFD simulations. This novel hybrid of techniques results in the ability to analyse performance when simulating various array layouts more rapidly and accurately than using either method on its own. It is shown that there is a significant power reduction from a downstream turbine that is subjected to the wake of an upstream turbine, and that this is due to both a reduction in power in the wind and also due to changes in the aerodynamics. Analysis of static pressure along the blade showed that as a result of wake interactions, a large reduction in the suction peak along the leading edge reduced the lift generated by the rotor and so reduced the torque production and the ability for the blade to extract energy from the wind