An Analysis of Journal Orbits for Nonlinear Dynamic Bearing Systems

An analysis of journal centre orbits is presented in this paper based on a non-isothermal non-Newtonian fluid model for dynamically loaded bearing systems. A spectral element approach is used to solve a full set of coupled equations (kinematics and constitutive) governing the flow of the lubricant, and an operator-splitting spectral element technique is used to evaluate the dynamic energy equation. The motion of the journal is calculated on the basis of Newtonian mechanics incorporated with a simple cavitation model. The stability of the journal orbits is investigated under a wide range of the rotation speeds of journal. The unstable orbits arise as a sub-harmonic motion when the journal rotation speed is increased beyond a critical value. The influences of the oscillation speeds of the applied loads on the journal orbits are examined. The numerical simulations demonstrate that both the rotation speed of the journal and the oscillation speed of the applied load play an important role in determining the pattern of the journal orbits. The effects of square-wave and rotating applied loads on the journal orbits are also investigated. Received 22 April 1998 and accepted 26 May 1999     

基于ALE方法的3D充填流动模拟

基于任意拉格朗日-欧拉方法发展了三维充填流动的数值模拟方案.该方案采用ALE方法准确地追踪移动自由面的位置并避免了网格扭曲;基于移动最小二乘曲面拟合方法提出了移动自由面上网格节点重定位方法,将充填流动的网格更新过程简化为自由面附近的局部网格重划分过程,并通过分级多面体三角剖分实现,减小了网格划分的计算量,实现了实时网格生成.给出的数值算例结果表明了该数值模型对三维充填流动模拟的有效性.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

A numerical study of the propulsive efficiency of a flapping hydrofoil

A computational fluid dynamics study of the swimming efficiency of a two‐dimensional flapping hydrofoil at a Reynolds number of 1100 is presented. The model accounts fully for viscous effects that are particularly important when flow separation occurs. The model uses an arbitrary Lagrangian–Eulerian (ALE) method to track the moving boundaries of oscillatory and flapping bodies. A parametric analysis is presented of the variables that affect the motion of the hydrofoil as it moves through the flow along with flow visualizations in an attempt to quantify and qualify the effect that these variables have on the performance of the hydrofoil. Copyright © 2003 John Wiley & Sons, Ltd.     

A fully discrete ALE method over untwisted time–space control volumes

In this paper, a fully discrete high‐resolution arbitrary Lagrangian–Eulerian (ALE) method is developed over untwisted time–space control volumes. In the framework of the finite volume method, 2D Euler equations are discretized over untwisted moving control volumes, and the resulting numerical flux is computed using the generalized Riemann problem solver. Then, the fluid flows between meshes at two successive time steps can be updated without a remapping process in the classic ALE method. This remapping‐free ALE method directly couples the mesh motion into a physical variable update to reflect the temporal evolution in the whole process. An untwisted moving mesh is generated in terms of the vorticity‐free part of the fluid velocity according to the Helmholtz theorem. Some typical numerical tests show the competitive performance of the current method. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of liquid sloshing in curved‐wall containers with arbitrary Lagrangian–Eulerian method

There are many challenges in the numerical simulation of liquid sloshing in horizontal cylinders and spherical containers using the finite element method of arbitrary Lagrangian–Eulerian (ALE) formulation: tracking the motion of the free surface with the contact points, defining the mesh velocity on the curved wall boundary and updating the computational mesh. In order to keep the contact points slipping along the curved side wall, the shape vector in each time advancement is defined to modify the kinematical boundary conditions on the free surface. A special function is introduced to automatically smooth the nodal velocities on the curved wall boundary based on the liquid nodal velocities. The elliptic partial differential equation with Dirichlet boundary conditions can directly rezone the inner nodal velocities in more than a single freedom. The incremental fractional step method is introduced to solve the finite element liquid equations. The numerical results that stemmed from the algorithm show good agreement with experimental phenomena, which demonstrates that the ALE method provides an efficient computing scheme in moving curved wall boundaries. This method can be extended to 3D cases by improving the technique to compute the shape vector. Copyright © 2007 John Wiley & Sons, Ltd.     

High‐order time integrators for front‐tracking finite‐element analysis of viscous free‐surface flows

This paper proposes implicit Runge–Kutta (IRK) time integrators to improve the accuracy of a front‐tracking finite‐element method for viscous free‐surface flow predictions. In the front‐tracking approach, the modeling equations must be solved on a moving domain, which is usually performed using an arbitrary Lagrangian–Eulerian (ALE) frame of reference. One of the main difficulties associated with the ALE formulation is related to the accuracy of the time integration procedure. Indeed, most formulations reported in the literature are limited to second‐order accurate time integrators at best. In this paper, we present a finite‐element ALE formulation in which a consistent evaluation of the mesh velocity and its divergence guarantees satisfaction of the discrete geometrical conservation law. More importantly, it also ensures that the high‐order fixed mesh temporal accuracy of time integrators is preserved on deforming grids. It is combined with the use of a family of L‐stable IRK time integrators for the incompressible Navier–Stokes equations to yield high‐order time‐accurate free‐surface simulations. This is demonstrated in the paper using the method of manufactured solution in space and time as recommended in Verification and Validation. In particular, we report up to fifth‐order accuracy in time. The proposed free‐surface front‐tracking approach is then validated against cases of practical interest such as sloshing in a tank, solitary waves propagation, and coupled interaction between a wave and a submerged cylinder. Copyright © 2015 John Wiley & Sons, Ltd.     

Provably stable and time‐accurate extensions of Runge–Kutta schemes for CFD computations on moving grids

Extending fixed‐grid time integration schemes for unsteady CFD applications to moving grids, while formally preserving their numerical stability and time accuracy properties, is a nontrivial task. A general computational framework for constructing stability‐preserving ALE extensions of Eulerian multistep time integration schemes can be found in the literature. A complementary framework for designing accuracy‐preserving ALE extensions of such schemes is also available. However, the application of neither of these two computational frameworks to a multistage method such as a Runge–Kutta (RK) scheme is straightforward. Yet, the RK methods are an important family of explicit and implicit schemes for the approximation of solutions of ordinary differential equations in general and a popular one in CFD applications. This paper presents a methodology for filling this gap. It also applies it to the design of ALE extensions of fixed‐grid explicit and implicit second‐order time‐accurate RK (RK2) methods. To this end, it presents the discrete geometric conservation law associated with ALE RK2 schemes and a method for enforcing it. It also proves, in the context of the nonlinear scalar conservation law, that satisfying this discrete geometric conservation law is a necessary and sufficient condition for a proposed ALE extension of an RK2 scheme to preserve on moving grids the nonlinear stability properties of its fixed‐grid counterpart. All theoretical findings reported in this paper are illustrated with the ALE solution of inviscid and viscous unsteady, nonlinear flow problems associated with vibrations of the AGARD Wing 445.6. Copyright © 2011 John Wiley & Sons, Ltd.     

Calculation of capillary jet

The Arbitrary Lagrangian Eulerian (ALE) framework coupled with some boundary tracking techniques is proven to be an effective method for simulation of free‐surface flows. In this paper, a special ALE framework is derived with clarification of three velocities, the notion of mesh‐frozen and field‐frozen, and the notion of tentatively inertial coordinates. A weighted integral ALE governing equations are formulated on generic coordinates and discretized with a finite element method and linear implicit time scheme. The system is solved with a discrete operator splitting technique and superposition‐based logistic parallelization. The formulation and implementation are verified through several fixed‐geometry problems and a reasonably good parallel performance is observed. Capillary jet flow is the main problem of the paper and the numerical techniques for boundary tracking are elaborated, which include an indirect boundary tracking of flux method and an iterative direct boundary tracking method. Also, a high‐order compact scheme for dynamic boundary condition and a squeeze technique for kinematic boundary condition are adopted. The axisymmetric jet breakup is studied in detail and numerical results match with the published data very well. Numerical accuracy and sensitivity are studied, including effects of element type, time scheme, compact scheme, and boundary tracking techniques. Copyright © 2010 John Wiley & Sons, Ltd.     

Topology optimization of a flexible multibody system with variable-length bodies described by ALE–ANCF

Recent years have witnessed the application of topology optimization to flexible multibody systems (FMBS) so as to enhance their dynamic performances. In this study, an explicit topology optimization approach is proposed for an FMBS with variable-length bodies via the moving morphable components (MMC). Using the arbitrary Lagrangian–Eulerian (ALE) formulation, the thin plate elements of the absolute nodal coordinate formulation (ANCF) are used to describe the platelike bodies with variable length. For the thin plate element of ALE–ANCF, the elastic force and additional inertial force, as well as their Jacobians, are analytically deduced. In order to account for the variable design domain, the sets of equivalent static loads are reanalyzed by introducing the actual and virtual design domains so as to transform the topology optimization problem of dynamic response into a static one. Finally, the novel MMC-based topology optimization method is employed to solve the corresponding static topology optimization problem by explicitly evolving the shapes and orientations of a set of structural components. The effect of the minimum feature size on the optimization of an FMBS is studied. Three numerical examples are presented to validate the accuracy of the thin plate element of ALE–ANCF and the efficiency of the proposed topology optimization approach, respectively.     

Dynamic analysis of baffled fuel‐storage tanks using the ALE finite element method

This paper is concerned with the parametric investigation on the structural dynamic response of moving fuel‐storage tanks with baffles. Since the structural dynamic behaviour is strongly coupled with interior liquid motion, the design of a fuel‐storage tank securing the structural stability becomes the appropriate suppression of the flow motion, which is in turn related to the baffle design. In order to numerically investigate the parametric dynamic characteristics of moving tanks, we employ the arbitrary Lagrangian–Eulerian (ALE) finite element method that is widely being used to deal with the problems with free surface, moving boundary, large deformation and interface contact. Following the theoretical and numerical formulations of fluid‐structure interaction problems, we present parametric numerical results of a cylindrical fuel‐storage tank moving with uniform vertical acceleration, with respect to the baffle number and location, and the inner‐hole diameter. Copyright © 2003 John Wiley & Sons, Ltd.     

An anisothermal,compressible, piezoviscous model for journal‐bearing lubrication

Compressibility plays a significant role in the load‐bearing capacity of a journal bearing. This paper offers more realistic modelling of the lubricant than presented in an earlier paper (Int. J. Numer. Meth. Fluids 2007; 55 (11):1091–1120) by including variable sound speed, piezoviscosity and both temperature and shear thinning. The load‐bearing capacity of the journal bearing is sensitive to all of these attributes of the model, but piezoviscosity is found to be the most significant. The equations of motion are adapted to a moving frame to explore the stability of the journal in a more dynamic setting and it is found that a free journal using this model will spiral outward exhibiting half speed whirl. The model is discretized in 2D, semi‐implicitly in time and using the spectral element method in space. Numerical results are presented that highlight the contributions of the different elements in the model to journal stability. Copyright © 2008 John Wiley & Sons, Ltd.     

有限长轴承非稳态油膜力自由-移动边界问题的有限条解法

讨论了不可压缩流体润滑的动载径向滑动轴承油膜压力分布的自由移动边界问题的有限条解法．将自由边界问题转化为全域(矩形域)的具有不等式约束的微分方程的边值问题，进一步化为具有不等式约束的泛函极值问题。借助有限条法在矩形域上离散这个泛函，得到了一个特殊的二次泛函的规化问题。通过变量平移变换，使其化为标准的二次规划问题。然后借助于牛顿非光滑算法，迭代求解非线性的互补方程。给出了有限长轴承真实的油膜应力分布。对于所求解方程的系数矩阵的高度稀疏性。给出了紧缩存储算法。节省了存储空间和减少了计算量。算例表明该方法是有效的。     

An arbitrary Lagrangian Eulerian (ALE) formulation for free surface flows using the characteristic‐based split (CBS) scheme

An arbitrary Lagrangian Eulerian (ALE) method for non‐breaking free surface flow problems is presented. The characteristic‐based split (CBS) scheme has been employed to solve the ALE equations. A simple mesh smoothing procedure based on coordinate averaging (Laplacian smoothing) is employed in the calculations. The mesh velocity is calculated at each time step and incorporated as part of the scheme. Results presented show an excellent agreement with the available experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

An accurate finite element scheme with moving meshes for computing 3D‐axisymmetric interface flows

An accurate finite element scheme for computing 3D‐axisymmetric incompressible free surface and interface flows is proposed. It is based on the arbitrary Lagrangian Eulerian (ALE) approach using free surface/interface‐resolved moving meshes. Key features like the surface force, consisting of surface tension and the local curvature, and jumps in the density and viscosity over different fluid phases are precisely incorporated in the finite element formulation. The local curvature is approximated by using the Laplace–Beltrami operator technique combined with a boundary approximation by isoparametric finite elements. A new approach is used to derive the 3D‐axisymmetric form from the variational form in 3D‐Cartesian coordinates. Several test examples show the high accuracy and the robustness of the proposed scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

成型充填过程的ALE有限元模拟

在ALE框架中提出了一个用于成型充填过程有限元数值模拟的模型。应用ALE参考构形及ALE参考粒子速度描写充填过程中的熔体质量运动。摒弃了Hele-Shaw近似假定,因而所提出的模型能用于非薄壁型腔中高分子材料充填过程的数值模拟。应用基于时域分步算法的Taylor-Galerkin方法,对控制成型充填过程的守恒方程建立了弱形式。对移动自由面附近的充填材料区构造了网格生成算法与网格重划分方案。给出了在几种不同形状的典型腔体中充填过程的数值模拟结果,表明了所提出的ALE有限元模型模拟充填过程的有效性。