Coupling of the Reynolds Fluid-Film Equation with the 2D Navier-Stokes Equations

The pressure field in thin fluid films can quite precisely be calculated by Reynolds fluid-film equation. In some problems, it may be useful to couple thin fluid-films with general 2D or 3D fluid flows. In the current work, we analyze the fluid flow, pressure and temperature field in a hydrodynamic journal bearing with a rectangular oil groove. Pressure and temperature in the fluid gap are calculated by means of the Reynolds equation and the 2D energy equation. Cavitation effects are taken into account by incorporating a 2-phase cavitation approach. In order to calculate the velocity and pressure field in the oil groove, the 2D Navier-Stokes equations are used; the temperature distribution in the oil groove is computed by means of the 2D energy equation. Appropriate coupling conditions for velocity, pressure and temperature are formulated in order to couple the flow in the fluid gap with the flow in the oil groove. Thermal expansion of journal shaft and bearing housing are also taken into account, since the bearing clearance changes with increasing temperature. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-isothermal lubrication problem with Tresca fluid–solid interface law. Part II Asymptotic behavior of weak solutions

We study in this paper the asymptotic analysis of an incompressible Newtonian and non-isothermal problem, when one dimension of the fluid domain tends to zero. We prove the strong convergence of the unknowns which are the temperature, the velocity and the pressure of the fluid, we obtain the limit problem with the specific Reynolds equation, and we also prove the uniqueness of the limit temperature velocity and pressure distributions.     

Fully Coupled Thermo-Hydrodynamic Simulation Model for Journal Bearings

The isothermal form of Reynolds fluid film equation is used to predict the pressure generation in hydrodynamic journal bearings if temperature effects are neglected. Often, however, temperature effects may be important and cannot be neglected, because oil viscosity significantly varies with temperature. Also, thermal expansion of journal shaft and bearing housing must be taken into account since the bearing clearance changes with increasing temperature. Hence, the Reynolds pressure field equation, the energy equation for the fluid film and the heat transfer equations for journal and bearing housing have to be solved simultaneously. The coupled thermo-hydrodynamic fluid flow problem is mathematically defined by a system of nonlinear integro-differential equations. The governing equations are discretized and solved by a finite element approach. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large Time Behavior for a Simplified 1D Model of Fluid–Solid Interaction†

Abstract

In this article we consider a simple model in one space dimension for the interaction between a fluid and a solid represented by a point mass. The fluid is governed by the viscous Burgers equation and the solid mass, which shares the velocity of the fluid, is accelerated by the difference of pressure at both sides of it. We describe the asymptotic behavior of solutions for integrable data using energy estimates and scaling techniques. We prove that the asymptotic profile of the fluid is a self-similar solution of the Burgers equation with an appropriate total mass, and we describe the parabolic trajectory of the point mass. We also prove that, asymptotically, the difference of pressure to both sides of the point mass vanishes.     

A model for fluid flow in non-fully filled tribological interfaces – Part 1: Basics and Numerics

In literature, most contributions on starved lubrication focus on the occurring pressures in macroscopic devices. Hereby, usually the Reynolds equation is modified in different ways. In contrast to this proceeding, this paper's intention is the general investigation of this tribological regime to get a fundamental comprehension on the transition from boundary lubrication to mixed lubrication. The respective model describes the flow of the fluid through two rough surfaces moving relative to each other. The lack of fluid is regarded by the fact that elements may not be fully filled with the fluid. Only elements where the fluid fully fills the gap, generate a pressure. This effect is considered by a type of unilateral constraint in combination with a penalty function. The fluid flow is computed according to the Navier-Stokes equation. In combination with the continuity equation, a set of implicit nonlinear equations has to be solved. Its potential and basic application fields are finally discussed. A further paper will show applications of the algorithm towards different scenarios. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

How the Reynolds equation is related to the Stokes equations

Applied Mathematics & Optimization - The Reynolds equation, which gives the pressure distribution in a thin layer of lubricating fluid, is usually motivated in a rather heuristic way. We make...     

Reduction Approaches for Thermogasdynamic Lubrication Problems

Gas thrust bearings are often used in low-load applications, e.g. in air cycle machines, in micro gas turbines or in rotor systems for fuel cell applications, to support a shaft in axial direction. The pressure and temperature distribution in a gas thrust bearing pad are described by the generalized Reynolds equation according to Dowson and the 3D energy equation. In this paper, two different approaches are presented in order to reduce the dimension of the governing nonlinear integro-differential equation system and in order to stabilize the solution process. In the first reduction approach, the temperature in the fluid is averaged across the fluid film according to Lee and Kim. In the second approach, Legendre polynomials are used to approximate temperature, density and fluidity across the fluid film according to Elrod, Brewe and Moraru. The reduction techniques are compared with respect to numerical efficiency, accuracy and convergence behavior. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling and simulation of injection molding process of polymer melt by a robust SPH method

In the present study, the injection molding process of polymer melt based on the generalized Newtonian fluid model is investigated by a robust smoothed particle hydrodynamics (SPH) method. The numerical method is proposed by introducing a Rusanov flux into the continuity equation to improve the prediction of the pressure distribution and employing a corrected kernel gradient to improve the computational accuracy. In addition, a robust treatment of solid boundary is presented and verified by the spin-down problem. The merits of the robust SPH method are firstly illustrated by 2D dam breaking flow. Then the numerical method is extended to deal with the flow phenomena related to injection molding process of polymer melt. A number of numerical examples including 2D injection moldings of a thin plate mold, a circular disc with core, a ring-shaped channel, and a S-shaped cavity, and 3D injection moldings of a Z-shaped cavity and a four-legged fork are conducted. The numerical results are in agreement with the experiments, which demonstrate that the SPH method proposed here is capable of handling with injection molding process of polymer melt in a robust manner. Moreover, the robust SPH method allows to recover the fluctuations-free pressure and velocity fields which in most cases cannot be easily obtained by the traditional SPH method.     

Multigrid preconditioned Krylov subspace methods for three-dimensional numerical solutions of the incompressible Navier–Stokes equations

We consider numerical methods for the incompressible Reynolds averaged Navier–Stokes equations discretized by finite difference techniques on non-staggered grids in body-fitted coordinates. A segregated approach is used to solve the pressure–velocity coupling problem. Several iterative pressure linear solvers including Krylov subspace and multigrid methods and their combination have been developed to compare the efficiency of each method and to design a robust solver. Three-dimensional numerical experiments carried out on scalar and vector machines and performed on different fluid flow problems show that a combination of multigrid and Krylov subspace methods is a robust and efficient pressure solver. This revised version was published online in June 2006 with corrections to the Cover Date.     

Expansions at small Reynolds numbers for the flow past a porous circular cylinder

The purpose of this article is to use the method of matched asymptotic expansions (MMAE) in order to study the two-dimensional steady low Reynolds number flow of a viscous incompressible fluid past a porous circular cylinder. We assume that the flow inside the porous body is described by the continuity and Brinkman equations, and the velocity and boundary traction fields are continuous across the interface between the fluid and porous media. Formal expansions for the corresponding stream functions are used. We show that the force exerted by the exterior flow on the porous cylinder admits an asymptotic expansion with respect to low Reynolds numbers, whose terms depend on the characteristics of the porous cylinder. In addition, by considering Darcy's law for the flow inside the porous circular cylinder, an asymptotic formula for the force on the cylinder is obtained. Also, a porous circular cylinder with a rigid core inside is considered with Brinkman equation inside the porous region. Stress jump condition is used at the porous–liquid interface together with the continuity of velocity components and continuity of normal stress. Some particular cases, which refer to the low Reynolds number flow past a solid circular cylinder, have also been investigated.     

微极流体润滑在径向轴承中应用

本文从微极流体场方程出发,在润滑层的通常假设下,把它化简为两个独立的常微分方程组,并求得速度、微转动角速度的解析表达式.推导了微极流体润滑的雷诺方程,把它应用于有限长径向轴承的求解.通过数值计算得到了微极效应对各种动力参数、几何参数下轴承的压力分布、承载力、流量系数和摩擦系数的影响,并析了它的实际意义,使微极流体理论应用到工程问题又接近了一步.     

Boundary layer due to the axisymmetric spreading of material on the surface of a fluid

The effect of the axisymmetric spreading of a layer of material (oil or solid particles) on the surface of a viscous fluid is studied. Assuming high Reynolds numbers, the boundary layer equation is derived and solved for general power law surface velocities. The composite streamlines show sharp turns near the surface.     

Numerical simulation of flow past a sphere in vertical motion within a stratified fluid

In the present study we consider a viscous fluid, stratified by a diffusive saline agent and compute numerically the flow produced by a solid sphere moving vertically and uniformly. The governing equations describing this situation are solved on a variational grid. The results show the dependence of the boundary-layer separation point and the vanishing of vortices behind the sphere as the stratification increases at moderate Reynolds number flows. Details of the flow, density and pressure fields near the sphere are also shown. Important quantities for engineering use (drags, pressure and skin coefficients) are also computed and displayed in the Richardson vs. Reynolds number space. Comparison with experimental evidence shows and excellent agreement.     

Forces acting on particles in separated wall‐bounded shear flow

Trajectories of solid particles in laminar and turbulent flow over a backward‐facing step (BFS) were numerically computed by integrating the equation of motion for particles. The various forces acting on the particles [5],[6] were calculated for a variety of flow Reynolds numbers and for different particle characteristics such as the Stokes number and the particle‐to‐fluid density ratio. The investigation was conducted for the distinct flow regimes of the BFS flow separately. Generally, the drag and gravitation were found to be the most significant forces. The lift and history force were the next most important, mostly two orders of magnitude smaller, but in some cases closing up to the other two in importance. The pressure and virtual mass effects were very small for the majority of cases. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

三维PTT黏弹性液滴撞击固壁面问题的改进SPH模拟

基于光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法,对三维Phan-Thien Tanner(PTT)黏弹性液滴撞击固壁面问题进行了数值模拟.为了有效地防止粒子穿透固壁,且缩减三维数值模拟所消耗的计算时间,提出了一种适合三维数值模拟的改进固壁边界处理方法.为了消除张力不稳定性问题,采用一种简化的人工应力技术.应用改进SPH方法对三维PTT黏弹性液滴撞击固壁面问题进行了数值模拟,精细地捕捉了液滴在不同时刻的自由面,讨论了PTT黏弹性液滴不同于Newton(牛顿)液滴的流动特征,分析了PTT拉伸参数对液滴宽度、高度和弹性收缩比等的影响.模拟结果表明,改进SPH方法能够有效而准确地描述三维PTT黏弹性液滴撞击固壁面问题的复杂流变特性和自由面变化特征.     

On a non-isothermal,non-Newtonian lubrication problem with Tresca law: Existence and the behavior of weak solutions

We consider a problem describing the motion of an incompressible, non-isothermal, and non-Newtonian fluid in a three-dimensional thin domain. We first establish an existence result for weak solutions of this problem. Then we study the asymptotic analysis when one dimension of the fluid domain tends to zero. A specific weak Reynolds equation, the limit of Tresca fluid–solid boundary conditions, and the limit boundary conditions for the temperature are obtained. The uniqueness result for the limit problem is also proved.     

Modulation equations and Reynolds averaging for finite-amplitude non-linear waves in an incompressible fluid

A formal perturbation scheme is developed to determine originalmodulation equations for laminar finite-amplitude non-linearwaves in an incompressible fluid. Three idealized problems areanalysed. The modulation equations comprise conservation ofwaves, averaged conditions for conservation of mass, momentum,kinetic energy and angular momentum and the averaged projectionof the Navier–Stokes equations onto the vorticity vector.The last of these modulation equations, which is related tovortex stretching, only appears in 3D problems. The techniqueof Reynolds averaging is also employed to obtain equations forthe mean velocities and pressure. The Reynolds-averaged Navier–Stokesequations correspond to the modulation equations for conservationof mass and momentum. However, the Reynolds stress transportequations are shown to be inconsistent with the other necessarymodulation equations. In two further idealized problems, exactsolutions of the Navier–Stokes equations are obtainedby employing the modulation equations.     

A new more consistent Reynolds model for piezoviscous hydrodynamic lubrication problems in line contact devices

Hydrodynamic lubrication problems in piezoviscous regime are usually modeled by the classical Reynolds equation combined with a suitable law for the pressure dependence of viscosity. For the case of pressure–viscosity dependence in the Stokes equation, a new Reynolds equation in the thin film limit has been proposed by Rajagopal and Szeri. However, these authors consider some additional simplifications. In the present work, avoiding these simplifications and starting from a Stokes equation with pressure dependence of viscosity through Barus law, a new Reynolds model for line contact lubrication problems is deduced, in which the cavitation phenomenon is also taken into account. Thus, the new complete model consists of a nonlinear free boundary problem associated to the proposed new Reynolds equation.Moreover, the classical model, the one proposed by Rajagopal and Szeri and the here proposed one are simulated through the development of some numerical algorithms involving finite elements method, projected relaxation techniques, duality type numerical strategies and fixed point iteration techniques. Finally, several numerical tests are performed to carry out a comparative analysis among the different models.     

Numerical and analytical simulation of peristaltic flow of a Jeffrey-six constant fluid

This paper describes the Peristaltic flow of a Jeffrey-six constant fluid in an endoscope. The two-dimensional equation of Jeffrey-six constant fluid is simplified by making the assumptions of long wave length and low Reynolds number. The reduced momentum equations are solved with three methods, namely (i) Perturbation method, (ii) Homotopy analysis method, and (iii) shooting method. The comparison of the three solutions shows a very good agreement between the three results. The expressions for pressure rise and frictional forces per wave length have been also computed numerically. Finally, the pressure rise, frictional forces are plotted for different parameters of interest.     

Numerical simulation of three-dimensional incompressible fluid in a box flow passage considering fluid–structure interaction by differential quadrature method

A numerical simulation scheme of 3D incompressible viscous fluid in a box flow passage is developed to solve Navier–Stokes (N–S) equations by firstly taking fluid–structure interaction (FSI) into account. This numerical scheme with FSI is based on the polynomial differential quadrature (PDQ) approximation technique, in which motions of both the fluid and the solid boundary structures are well described. The flow passage investigated consists of four rectangular plates, of which two are rigid, while another two are elastic. In the simulation the elastic plates are allowed to vibrate subjected to excitation of the time-dependent dynamical pressure induced by the unsteady flow in the passage. Meanwhile, the vibrating plates change the flow pattern by producing many transient sources and sinks on the plates. The effects of FSI on the flow are evaluated by running numerical examples with the incoming flow’s Reynolds numbers of 3000, 7000 and 10,000, respectively. Numerical computations show that FSI has significant influence on both the velocity and pressure fields, and the DQ method developed here is effective for modelling 3D incompressible viscous fluid with FSI.