Modified Reynolds Equation for Squeeze-Film Air Damping of Slotted Plates in MEMS Devices

A differential equation for calculating squeeze-film air damping in slotted plates is developed by modifying the Reynolds equation. A term is added to account for the effect of airflow through the slots on the air damping of the plate. The end effect of the airflow in the slots is also treated by substituting an effective channel length for the geometric channel length （i. e. the thickness of the plate）..The damping pressure distribution, damping force, and damping force coefficient of the slotted plates can be found by solving the equation under appropriate boundary conditions. With restrictions on the thickness and the lateral dimensions of the slotted plate removed,the equation provides a useful tool for analysing the squeeze-film air damping effect of slotted plates with finite thickness and finite lateral dimensions. For a typical slotted plate structure, the damping force coefficient obtained by this equation agrees well with that generated by ANSYS.

MODELING THE INTERACTION OF SOLITARY WAVES AND SEMI-CIRCULAR BREAKWATERS BY USING UNSTEADY REYNOLDS EQUATIONS

A vertical 2-D numerical wave model was developed based on unsteady Reynolds equations. In this model, the k-epsilon models were used to close the Reynolds equations, and volume of fluid(VOF) method was used to reconstruct the free surface. The model was verified by experimental data. Then the model was used to simulate solitary wave interaction with submerged, alternative submerged and emerged semi-circular breakwaters. The process of velocity field, pressure field and the wave surface near the breakwaters was obtained. It is found that when the semi-circular breakwater is submerged, a large vortex will be generated at the bottom of the lee side wall of the breakwater; when the still water depth is equal to the radius of the semi-circular breakwater, a pair of large vortices will be generated near the shoreward wall of the semi-circular breakwater due to wave impacting, but the velocity near the bottom of the lee side wall of the breakwater is always relatively small. When the semi-circular breakwater is emerged, and solitary wave cannot overtop it, the solitary wave surface will run up and down secondarily during reflecting from the breakwater. It can be further used to estate the diffusing and transportation of the contamination and transportation of suspended sediment.

A GENERALIZED REYNOLDS EQUATION BASED ON NON-NEWTONIAN FLOW IN LUBRICATION MECHANICS

A generalized Reynolds equation based on non-Newtonian flow is derived in this paper.This equation is suitable for a number of non-Newtonian flow models and can be solved numerically to obtain pressure fields in thermalhydrodynamically or elastohydrodynamically lubricated fluid films.A mathematical ap- proach is given for solving simultaneously the shearing stress,shearing rate,velocity and equivalent viscosity.To show the application of this equation,two rheological models which have been widely used in lubrication mechnaics are incorporated into this equation to obtain numerical solutions to the line contact thermal elastohydrodynamic lubrication problem.

A new method to solve the Reynolds equation including mass-conserving cavitation by physics informed neural networks(PINNs)with both soft and hard constraints

In this work,a new method to solve the Reynolds equation including mass-conserving cavitation by using the physics informed neural networks(PINNs)is proposed.The complementarity relationship between the pressure and the void fraction is used.There are several difficulties in problem solving,and the solutions are provided.Firstly,the difficulty for considering the pressure inequality constraint by PINNs is solved by transferring it into one equality constraint without introducing error.While the void fraction inequality constraint is considered by using the hard constraint with the max-min function.Secondly,to avoid the fluctuation of the boundary value problems,the hard constraint method is also utilized to apply the boundary pressure values and the corresponding functions are provided.Lastly,for avoiding the trivial solution the limitation for the mean value of the void fraction is applied.The results are validated against existing data,and both the incompressible and compressible lubricant are considered.Good agreement can be found for both the domain and domain boundaries.

Homogenization of Reynolds Equation by Two-Scale Convergence

To increase the hydrodynamic performance in different machine elements, as e.g. journal bearings and thrust bearings, during lubrication it is important to understand the influence of surface roughness. In this connection one encounters homogenization of the incompressible Reynolds equation, where the roughness of the lubricated surface is assumed to be periodic. This problem has recently been studied in more engineering- oriented papers by using the formal method of multiple scale expansion. In this paper, we rigorously prove both homogenization and corrector results by using two-scale convergence, which may be regarded as a justification of the formal multiple scale expansion method described above. Moreover, some numerical illustrations and results are presented.

Water Profile over an Ogee-Shaped Crest and the Reynolds Fractional Equations

Weir crest must have the correct shape in the concave side of an ogee-shaped crest to diminish erosion. This shape can be obtained using an approximation of the fractional Reynolds equations when the water interacts with the surface. A model is introduced for the Reynolds stresses complemented with a closure relation of fractional origin. A power type solution is obtained for the main velocity and stress. Velocity profile is found based on the assumption of a steady flow and the energy conservation equation. A Froude number and the cubic equation of the weir are built. The dimensionless upstream velocity head is also determined which allow graphically showing the exponent and coefficient of the water-profile over an ogee-shaped crest. It is possible to estimate the occupied-space index trough an exponents＇ ratio of profile over the velocity head.

Derivation of a second-order model for Reynolds stress using renormalization group analysis and the two-scale expansion technique

With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples：turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder （GL） model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.

Improved finite difference method for pressure distribution of aerostatic bearing

An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses （less than 2 Win）are considered, the improved calculation method still yields a result and converges fast.

Numerical Analysis of a Spiral-groove Dry-gas Seal Considering Micro-scale Effects

A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication processes,while forming a high pressure gas film between two sealing faces due to the deceleration of the gas pumped in or out.There is little research into the effects and the influence on seal performance,if the grooves and the gas film are at the micro-scale.This paper investigates the micro-scale effects on spiral-groove dry-gas seal performance in a numerical solution of a corrected Reynolds equation.The Reynolds equation is discretized by means of the finite difference method with the second order scheme and solved by the successive-over-relaxation（SOR） iterative method.The Knudsen number of the flow in the sealing gas film is changed from 0.005 to 0.120 with a variation of film depth and sealing pressure.The numerical results show that the average pressure in the gas film and the sealed gas leakage increase due to micro-scale effects.The open force is enlarged,while the gas film stiffness is significantly decreased due to micro-scale effects.The friction torque and power consumption remain constant,even in low sealing pressure and spin speed conditions.In this paper,the seal performance at different rotor face spin speeds is also described.The proposed research clarifies the micro-scale effects in a spiral-groove dry-gas seal and their influence on seal performance,which is expected to be useful for the improvement of the design of dry-gas seal systems operating in the slip flow regime.

THERMAL ANALYSIS OF FRICTIONAL DISK IN SPEEDING WET CLUTCH

Multi-frictional disks are employed to transmit the torque in speeding wetclutch, and the oil thickness within frictional disks could be adjusted for practical outputspeeding. As oil combined with alpha-hydrocarbon or polyester is getting widely used as lubricantand the speeding wet clutch works within hydrodynamic lubrication, mixture lubrication, boundarylubrication and contact situation, established the thermal analysis model for investigating thebehavior of frictional disks in speeding wet clutch, which covers the power-law fluid model,Patir-Cheng average flow model, GT asperity contact model, oil film inertia and heat effects, meanenergy equation, and heat conduction equation. The formulas in the model are deduced. The numeralcalculation and analysis for hydrodynamic lubrication and mixture lubrication are executed.Relationships and variations among transmitting torque, mean push pressure, output speed, and loadare presented. Thermal effect should be considered during the hydrodynamic lubrication and largertransmitting torque makes earlier entrance to mixture lubrication.

Analysis of Thermoelastohydrodynamic Performance of Journal Misaligned Engine Main Bearings

To understand the engine main bearings＇ working condition is important in order to improve the performance of engine. However, thermal effects and thermal effect deformations of engine main bearings are rarely considered simultaneously in most studies. A typical finite element model is selected and the effect of thermoelastohydrodynamic（TEHD） reaction on engine main bearings is investigated. The calculated method of main bearing＇s thermal hydrodynamic reaction and journal misalignment effect is finite difference method, and its deformation reaction is calculated by using finite element method. The oil film pressure is solved numerically with Reynolds boundary conditions when various bearing characteristics are calculated. The whole model considers a temperature-pressure-viscosity relationship for the lubricant, surface roughness effect, and also an angular misalignment between the journal and the bearing. Numerical simulations of operation of a typical I6 diesel engine main bearing is conducted and importance of several contributing factors in mixed lubrication is discussed. The performance characteristics of journal misaligned main bearings under elastohydrodynamic（EHD） and TEHD loads of an I6 diesel engine are received, and then the journal center orbit movement, minimum oil film thickness and maximum oil film pressure of main bearings are estimated over a wide range of engine operation. The model is verified through the comparison with other present models. The TEHD performance of engine main bearings with various effects under the influences of journal misalignment is revealed, this is helpful to understand EHD and TEHD effect of misaligned engine main bearings.

ANALYSIS ON INFILTRATIVE CHARACTERISTICS OF DEFORMABLE POROUS MATAL RUBBER MATERIAL AND PARAMATERS IDENTIFICATION

The effect of deformation of porous material on infiltrative performance is investigated. Based on Darcy theory and Boit principle, the Reynolds equation and mathematical expression of deformable metal rubber （MR） material under laminar flow are obtained according to the change of porosity of metal rubber. It is shown that the throttle of MR material is dependent on its porosity and diameter of metal wires. It will be of great value for the application of MR in throttle field.

INVESTIGATION OF RATTLING THRESHOLD IN GEAR SYSTEM

A lumped mass gear-rattling model with backlash is established by considering the time varying mesh stiffness, composite transmission error and the torque fluctuation. Based on the principle that no separation will occur if the response amplitude is not larger than the static response, the threshold is obtained by theoretical derivation. The validity of the theoretical derivation is verified by the numerical method. The results show that the time-varying mesh stiffness has little effect on the threshold. When the exciting frequency is less than the system natural frequency, the theoretical agrees well with the numerical one. The higher the exciting frequency or the smaller the average input torque is, the easier the separation will occur. When the exciting frequency is larger than the natural frequency by a certain value, the theoretical threshold does not agree with the numerical results. The numerical results show that the motion state will change from the meshing state to the separating state directly. The phenomenon that the gear pair sometimes separates and sometimes meshes occurs only when the exciting frequency is smaller.

THE INTERACTIONS BETWEEN WAVE-CURRENTS AND OFFSHORE STRUCTURES WITH CONSIDERATION OF FLUID VISCOSITY

Study of the how held around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the how held around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averaged k-epsilon turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds number Re = 10(5), a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies.

ASYMPTOTIC CONVERGENCE OF A GENERALIZED NON-NEWTONIAN FLUID WITH TRESCA BOUNDARY CONDITIONS

The goal of this article is to study the asymptotic analysis of an incompressible Herschel-Bulkley fluid in a thin domain with Tresca boundary conditions.The yield stress and the constant viscosity are assumed to vary with respect to the thin layer parameterε.Firstly,the problem statement and variational formulation are formulated.We then obtained the existence and the uniqueness result of a weak solution and the estimates for the velocity field and the pressure independently of the parameterε.Finally,we give a specific Reynolds equation associated with variational inequalities and prove the uniqueness.

Influence of Oil Film Inertia on Piston Skirt Lubrication in High Speed Engine

By an iteration method, the influence of oil film inertia on piston skirt lubrication in a high speed engine is investigated. By alternately solving the Navier-Stocks equations and the Reynolds equation, the new iteration method can trace the variations of velocity field and pressure field with time. Based on this, the mixed-lubrication model suitable for the piston skirt of high engines is proposed. By introducing the inertia coefficient, the new lubrication model includes the inertia term in oil film. The model can be also used to solve for the lubrication performances of a piston skirt in low or medium speed engines and for lubrication problem in general excluding the inertia term of oil film , when the inertia coefficient is put equal to zero. The calculation results show that the influence of oil film inertia on the friction force increases with the ratio of the piston skirt's length to its diameter, the inertia coefficient and the eccentricities of the lower and the upper piston skirt, with other conditions kept constant, while the influence of it on the load capacity of oil film is small.

A solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation

This letter presents an analytical solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation.The solution is found by solving the Reynolds equation with Reynolds boundary condition for cavitation.The cavitation boundary is shown to be straight lines directed 108.4°against the sliding direction.The result is experimentally verified in the limit of large values of viscosity,sliding velocity and radius of a spherical ball.The solution raises questions about the coupling between cavitation and film rupture and can be used as an independent check on the validity of numerical solutions.

Flow Dynamics of a Spiral-groove Dry-gas Seal

The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.

Numerical Simulation on Hydraulic Performances of Quarter Circular Breakwater

Quarter circular breakwater （QCB） is a new-type breakwater developed from senti-circular breakwater （SCB）. The superstructure of QCB is composed of a quarter circular front wall, a horizontal base slab and a vertical rear wall. The width of QCB＇ s base slab is about half that of SCB, which makes QCB suitable to be used on relatively finn soil foundation. The numerical wave flume based on the Reynolds averaged Navier-Stokes equations for impressible viscosity fluid is adopted in this paper to simulate the hydraulic performances of QCB. Since the geometry of both breakwaters is similar and SCB has been studied in depth, the hydraulic performances of QCB are given in comparison with those of SCB.