The Effect of Non-Circular Bearing Shapes in Hydrodynamic Journal Bearings on the Vibration Behavior of Turbocharger Structures

Increasing quality demands of combustion engines require, amongst others, improvements of the engine’s acoustics and all (sub)components mounted to the latter. A significant impact to the audible tonal noise spectrum results from the vibratory motions of fast-rotating turbocharger rotor systems in multiple hydrodynamic bearings such as floating bearing rings. Particularly, the study of self-excited non-linear vibrations of the rotor-bearing systems is crucial for the understanding, prevention or reduction of the noise and, consequently, for a sustainable engine acoustics development. This work presents an efficient modeling approach for the investigation, optimization, and design improvement of complex turbocharger rotors in hydrodynamic journal bearings, including floating bearing rings with circular and non-circular bearing geometries. The capability of tonal non-synchronous vibration prevention using non-circular bearing shapes is demonstrated with dynamic run-up simulations of the presented model. These findings and the performance of our model are compared and validated with results of a classical Laval/Jeffcott rotor-bearing model and a specific turbocharger model found in the literature. It is shown that the presented simulation method yields fast and accurate results and furthermore, that non-circular bearing shapes are an effective measure to reduce or even prevent self-excited tonal noise.


Introduction
The recent development of combustion engines shows not only a clear trend towards "downsizing", i.e., the reduction of quantity and/or volumetric size of the cylinders but also an increased demand of higher specific power of such new engines together with lower carbon dioxide emissions.These requirements may only be achieved by using highly efficient turbocharger systems.
In addition, the representative customer of the premium segment vehicle demands significantly higher expectations regarding vibration and acoustic comfort.Engine noises have been continuously reduced, such that the noise emitted by turbocharging systems has become a focus of the engine sound design.A special case of disturbing turbocharger noises are tonal noises which result from imbalance-excited or self-excited resonances of the rotor-bearing system.In [1] it was demonstrated how the subjective perception of tonal turbocharger noises can be assessed objectively.
Turbochargers are usually supported by hydrodynamic floating bearing rings, a bearing concept consisting of a freely rotating journal bearing bushing with separate oil films between the bushing and both the shaft and the center housing.Depending on the size and the load of the application other bearing concepts exist using a non-rotating bearing ring.These structures allow very high rotational speeds, but come with the challenge of acoustically perceivable self-excited shaft vibrations.Literature often refers to these vibrations as "oil whirl", "constant tone" or "subsynchronous whining" [2][3][4][5][6][7].
A good understanding of the cause and impact of self-excited vibratory motions in exhaust gas turbocharger systems is crucial to enhance the comfort of premium segment vehicles.To understand the dynamics of complex turbocharger rotors, one possible approach is the study of reduced models such as the Laval/Jeffcott rotor.In an earlier publication [7], the self-excited shaft whirls of customized Laval/Jeffcott rotor simulation models were identified and a sensitivity study of rotor dimensions, bearing geometry and tribological properties was performed.Furthermore, it was shown that based on the Laval/Jeffcott rotor model, a system optimization can be conducted regarding reductions of self-excited shaft motion durations as well as displacement amplitudes.In the literature measures of how to reduce displacement amplitudes of plain hydrodynamic bearing self-excited vibrations can be found, however, the complete removal of the addressed vibrations of single and multiple bearing assemblies, including their modeling, is rarely covered.The stabilizing effect of asymmetrical or even non-circular bearing shapes is known in the literature, e.g., lemon bore bearings or tilting pad bearings [8,9], however, extended research and especially a validated and efficient simulation method of non-circular hydrodynamic journal bearings is hardly reported to the best of the authors' knowledge.Using a complete turbocharger model, the present work shows the possibility of completely preventing self-excited motions using non-circular bearing designs.For comparison and model validation, the non-circular shape was introduced into the abovementioned customized Laval rotor and into a turbocharger model from the literature, [10].Analyzed orbits and spectrograms of these models show detailed insights about the sensitivities of onset rotational speeds, durations, and frequency characteristics of self-excited motions.As the outlook, a detailed comparison of simulation results and advanced measurements of a turbocharger during run-ups on a component hot gas test rig are subject of ongoing and future work.

Rotordynamic Model for Efficient Run-Up Simulations in HOTINT
In the present work we utilize the simulation tool HOTINT (High Order Time Integration, [11]) which provides a versatile framework for the efficient numerical modeling and simulation of general flexible multibody systems, and furthermore, includes functionality for simulator coupling, visualization, post-processing, and optimization (for further information, see, for instance, [12]).HOTINT is available online (www.hotint.org)both as maintained freeware release and open source version.
Used as simulation tool for a broad range of applications in mechanics and mechatronics, as well as coupled multi-physics problems, HOTINT features a modular structure which allows for an efficient implementation of extensions and additional functionality.In the context of the present work, the software was extended by hydrodynamic bearing elements based on analytical models (short and long bearing theory) and a numerical solution of the Reynolds equation.
The rotordynamic models used in the present work were defined in terms of the HOTINT script language which allows for a fully parametrized and generic model setup.Further details on the general structure of multibody simulation models in HOTINT, the employed model components, and the hydrodynamic bearing modeling are presented in the following sections.

Element and Model Description
A turbocharger rotor system consists of several rotating and non-rotating parts.For the simulation of the vibratory motion of the turbocharger, however, the rotating shaft with turbine and compressor wheel as well as the radial bearing elements (here, hydrodynamic bearings) are the most important components.Other elements such as the thrust bearing, spacer elements, or turbocharger housing, are not modeled in the present investigations or accounted for in a reduced way.For instance, the thrust bearing and spacer are represented by effective shaft stiffnesses, but not modeled as separate parts.Axial thrust is not considered in the model.Regarding the bearings, since we are primarily interested in the radial movement of the shafts and the bearing rings (in case of floating bearing rings), only the radial hydrodynamic forces are taken into account.The axial motion is not considered (to be exact, kinematically constrained) in the present investigations, since the axial hydrodynamic bearing forces are negligible compared to the radial forces and the coupling between axial and the radial motion is assumed to be small due to the small tilt angles of the shaft.Note that, furthermore, no tilting motion of the floating rings is taken into account.
Figure 1 shows a sketch of the basic object structure in HOTINT.In general, a model is composed of a set of rigid or flexible bodies which may be independent of each other or linked together by means of some kind of connectors (kinematic constraints, contact, etc.).The bodies may be subject to external loads, and furthermore, so-called IO elements (Input-Output elements) may be used to implement control or dynamically modify model parameters during a simulation.The run-up profile, for instance, is implemented by a set of IO elements which compute the angular frequency according to a specified mathematical expression in dependency of time and prescribe the result by modifying the offset of a kinematic velocity constraint with respect to the corresponding degree of freedom of a beam element node.Finally, sensors are used to measure physical quantities such as (generalized) displacements, velocities, and forces which in turn may be used either for post-processing purposes or as input for IO elements.instance, the thrust bearing and spacer are represented by effective shaft stiffnesses, but not modeled as separate parts.Axial thrust is not considered in the model.Regarding the bearings, since we are primarily interested in the radial movement of the shafts and the bearing rings (in case of floating bearing rings), only the radial hydrodynamic forces are taken into account.The axial motion is not considered (to be exact, kinematically constrained) in the present investigations, since the axial hydrodynamic bearing forces are negligible compared to the radial forces and the coupling between axial and the radial motion is assumed to be small due to the small tilt angles of the shaft.Note that, furthermore, no tilting motion of the floating rings is taken into account.
Figure 1 shows a sketch of the basic object structure in HOTINT.In general, a model is composed of a set of rigid or flexible bodies which may be independent of each other or linked together by means of some kind of connectors (kinematic constraints, contact, etc.).The bodies may be subject to external loads, and furthermore, so-called IO elements (Input-Output elements) may be used to implement control or dynamically modify model parameters during a simulation.The run-up profile, for instance, is implemented by a set of IO elements which compute the angular frequency according to a specified mathematical expression in dependency of time and prescribe the result by modifying the offset of a kinematic velocity constraint with respect to the corresponding degree of freedom of a beam element node.Finally, sensors are used to measure physical quantities such as (generalized) displacements, velocities, and forces which in turn may be used either for post-processing purposes or as input for IO elements.The rotordynamic model employed in the present work relies on rotor beam elements based on a linear beam formulation, including the Rayleigh model to represent material damping [13], which are used to discretize the turbocharger's shaft along its rotor axis combined with rigid disk elements which correctly represent gyroscopic effects (e.g., of the turbine or compressor wheel).The hydrodynamic bearings are implemented as special constraints which act on the respective nodes of the rotor beam elements.Further details on the implemented bearing formulation are given in Section 2.4.Imbalances can be modeled either by point masses connected to a local position on beam or disk elements, or by centrifugal forces  P,  P, (2πΩ S ) 2 , where  P, denotes the k-th imbalance mass,  P, the respective radial distance, and Ω  the rotor's angular frequency.The rotordynamic model employed in the present work relies on rotor beam elements based on a linear beam formulation, including the Rayleigh model to represent material damping [13], which are used to discretize the turbocharger's shaft along its rotor axis combined with rigid disk elements which correctly represent gyroscopic effects (e.g., of the turbine or compressor wheel).The hydrodynamic bearings are implemented as special constraints which act on the respective nodes of the rotor beam elements.Further details on the implemented bearing formulation are given in Section 2.4.Imbalances can be modeled either by point masses connected to a local position on beam or disk elements, or by centrifugal forces m P,k r P,k (2πΩ S ) 2 , where m P,k denotes the k-th imbalance mass, r P,k the respective radial distance, and Ω S the rotor's angular frequency.

Equations of Motion and Time Integration
The equations of motion for the considered rotordynamic system within the framework of multibody system dynamics may be written as: where q is the vector of generalized coordinates, M denotes the mass matrix, D is a matrix which includes both damping as well as the gyroscopic terms, K is the stiffness matrix, and F is a vector of (generalized) external forces which may depend on the generalized coordinates and velocities.C(q) denotes a set of algebraic equations which represent kinematic constraints, and λ is the corresponding vector of Lagrange multipliers.
The global system of Equations ( 1) is assembled from element-wise contributions.Each body contributes a mass matrix, flexible bodies additionally contribute an element-wise stiffness matrix, and damping and gyroscopic terms are accounted for in D. Connector elements can either be described by algebraic equations and associated Lagrange multipliers, or, in terms of a Penalty formulation, included in the force vector F q, .q .The latter is the case for the employed hydrodynamic bearing elements.The integration of the above governing equations with respect to time is performed numerically in HOTINT by means of stable implicit integrators based on Lobatto-IIIA or Radau-IIA formulas [12,14].The algebraic equations are typically reduced to index two (velocity-level) to enhance stability.Note that the efficient discretization of the rotor by means of beam elements and a fast implementation of the bearing elements within the discussed multibody dynamics framework (cf.also the following section) allow for efficient dynamic analyses of the non-linear rotor-bearing system.A typical computational effort for the simulation of a rotor run-up within 10 s of physical time, including the numerical solution of the Reynolds Equation (3) for each bearing, amounts to a few hours CPU time, which is still sufficient to reasonably conduct parameter studies.
To perform eigenvalue analyses, all algebraic constraints are replaced by stiff penalty constraints, such that no algebraic equations C(q) remain, and the system is linearized about a configuration q 0 , .q 0 which may either be its initial state or the current solution resulting from a dynamic simulation.In any case, the generalized eigenvalue problem M .. q + D + ∂F ∂ .q q 0 , .q 0 .q + K + ∂F ∂q q 0 , .q 0 q = 0 with q(t) = q k e λ k t and is solved, where q k denotes the k-th eigenvector corresponding to the complex eigenvalue Here, the imaginary part λ Im k represents the eigenmode's angular frequency, whereas a negative real part λ Re k corresponds to its damping.

Hydrodynamic Bearing Elements with Finite Length and Non-Circular Geometry
The starting point for the analysis and modeling of hydrodynamic bearings with circular geometry without a tilt between shaft and shell is the Reynold's equation for the hydrodynamic pressure field p(φ, x), where φ and x denote the cirumferential and axial coordinate of the oil film, h the oil film thickness, R i and R a the inner and outer bearing radius (i.e., the radius of shaft and shell), ω i and ω a the respective angular velocities, and η the dynamic viscosity of the lubricant [15].Here, the local oil film thickness h only depends on φ and the relative position of the centers of shaft and shell.For the most general case, i.e., non-circular bearing shapes including axial and radial grooves, as well as a tilt of the shaft, Now, the oil film thickness h is a function of both coordinates φ and x, as well as of the relative position and orientation between shaft and shell.Note that R i and R a here are functions of φ as well as the respective angular positions of shaft and shell and represent any general non-circular bearing geometry.Whereas Equation (3) can be simplified for the limit cases of infinitely short and long bearings and solved analytically (based on the short and long bearing theory, cf.[15] or [16]), the solution of (3) for finite-length bearings and of the generalized form (4) can only be obtained numerically.
In the present work, a finite differences scheme is used to discretize the governing equations.Depending on the dominance of the convective parts of Equation (4), i.e., the terms proportional to ∂h ∂φ ∂p ∂φ and ∂h ∂x ∂p ∂x , the scheme is switched from second-order central differences to a first-order upwind scheme (forward or backward differences) to preserve stability.Figure 2 shows a sketch of the computational domain with the regular finite difference grid of the [ϕ 1 , ϕ end ]×[x 1 , x end ], with x = 2x/L and the bearing length L. As to the boundary conditions of the differential equation, first of all, the periodic boundary conditions require p(ϕ + 2π, x) = p(ϕ, x), and specifically, p(ϕ 1 ) = p(ϕ end ), for the pressure field.Furthermore, the inlet pressure needs to be specified at the bearing's ends, i.e., p(x 1 ) = p in,1 and p(x end ) = p in,end , respectively, and can be additionally set on an oil inlet zone [ϕ sup,1 , ϕ sup,end × x sup,1 , x sup,end (cf. the orange area in the sketch in Figure 2 as well as Figure 3 for illustration).The additional inlet area does not need to conform to the underlying finite difference grid and is accounted for in a consistent way by adapted finite difference formulas in adjacent grid cells.For illustration of this approach, the second-order difference formulas used to approximate first and second derivatives with respect to x in a grid point i adjacent to a cell partly lying inside the inlet area (cut cell) are given by respectively.Here, ∆x denotes the regular (equidistant) resolution of the finite difference grid (i.e., the cell size), and α i is the fraction of a cut cell (in the direction of x) that lies outside the inlet area.Similar expressions can be derived for first-order forward or backward finite difference formulas for the upwinding scheme.Finally, since no cavitation model is used in the present investigations, the pressure is cut off at p = 0 in order to prevent the unphysical occurrence of negative pressure values.
The transition to the zero-pressure cut-off is made smooth, i.e., no discontinuity in the pressure field nor its gradient occurs, which is consistent with the so-called Reynolds boundary condition [15].
The system of linear equations resulting from the finite difference discretization is finally solved by SOR (Successive Over Relaxation) Gauss-Seidel iterations for both a preconditioning phase as well as the actual solution phase [17], which yields the discrete pressure field at all points of the finite difference grid.See Figure 2 for the schematic of the computational domain and Figure 3 for the illustration of a computed pressure field with an oil inlet area.Based on the pressure field , the velocity field , and subsequently, the hydrodynamic stress tensor consisting of the pressure contribution   = − and the viscous contribution   = η(∇ + ∇ T ) can be computed (for further details on the relation between pressure and velocity field cf.again [15]).Consequently, the bearing forces  and moment  about the rotor axis acting on the shaft and the corresponding counterforces and moments on the shell are given by integration over the bearing surface  as Figure 2. Schematic of the computational domain of the pressure field p(ϕ, x, t).A rectangular zone for a pressure supply area is marked in orange.
Lubricants 2017, 5, 6 6 of 19 The system of linear equations resulting from the finite difference discretization is finally solved by SOR (Successive Over Relaxation) Gauss-Seidel iterations for both a preconditioning phase as well as the actual solution phase [17], which yields the discrete pressure field at all points of the finite difference grid.See Figure 2 for the schematic of the computational domain and Figure 3 for the illustration of a computed pressure field with an oil inlet area.Based on the pressure field , the velocity field , and subsequently, the hydrodynamic stress tensor consisting of the pressure contribution   = − and the viscous contribution   = η(∇ + ∇ T ) can be computed (for further details on the relation between pressure and velocity field cf.again [15]).Consequently, the bearing forces  and moment  about the rotor axis acting on the shaft and the corresponding counterforces and moments on the shell are given by integration over the bearing surface  as The system of linear equations resulting from the finite difference discretization is finally solved by SOR (Successive Over Relaxation) Gauss-Seidel iterations for both a preconditioning phase as well as the actual solution phase [17], which yields the discrete pressure field at all points of the finite difference grid.See Figure 2 for the schematic of the computational domain and Figure 3 for the illustration of a computed pressure field with an oil inlet area.
Based on the pressure field p, the velocity field v, and subsequently, the hydrodynamic stress tensor consisting of the pressure contribution σ P = −pI and the viscous contribution σ visc = η ∇v + ∇v T can be computed (for further details on the relation between pressure and velocity field cf.again [15]).
Consequently, the bearing forces F and moment M about the rotor axis acting on the shaft and the corresponding counterforces and moments on the shell are given by integration over the bearing surface S as respectively, where n = (0, cos φ, sin φ) denotes the surface normal, and the rotor axis is aligned with the x-direction of the coordinate system.The dominating contribution to the bearing force (8) originates from the pressure terms.Note that S is always assumed as cylindrical bearing surface with a constant radius of R and its axis aligned with the x-axis of the coordinate system.For circular geometries, R is set to R a , in case of general bearing geometries, a reasonable value, e.g., the circumferential average of R a is chosen.The errors due to integration over that aligned cylindrical bearing surface S instead of the actual surface, possibly with additional tilt between shaft and shell, are on the order O(c) 1, with the bearing clearance c R i , R a , and O ∂h ∂φ 1, and thus, usually negligible in agreement with the underlying assumptions for the Reynolds Equations ( 3) and (4).Only in the vicinity of grooves, the term ∂h ∂φ may be large, which violates the above assumptions, yet introduces only local errors over a small range of φ.
The integrals (8) and ( 9) are finally evaluated numerically according to the discrete finite difference solution.The pressure contribution of the bearing force in (8), for instance, is approximated by with ∆ϕ = 2π/n ϕ and ∆x = L/n x .The resulting forces and moments are applied to the corresponding bodies in the multibody system and accounted for in the non-linear force vector F q, .q of the equation of motion (1).In case of full-floating bearing rings (without rotation constraints), where the floating ring is modeled as rigid body, the sum of the moments (9) from both the inner and outer oil film determines the ring speed about the rotor axis in the dynamic simulations.

Post Processing Methods
Post processing was done in Matlab 2015 [18].For the computation of the spectrograms based on data of 10 s run-up simulations, a fast Fourier transform (FFT) time window with a length of ∆t FFT = 0.11 s was used, which was shifted from 0 to 10 s consecutively by a time span of ∆t FFT,Shift = 0.6 ∆t FFT .The FFT frequency resolution is determined by the FFT time window length as ∆ f FFT = ∆t −1 FFT .The maximum analyzable frequency is given by the Nyquist criterion as the inverse of the time resolution (time step size) of the numerical solution data.

Results
This section shows results of simulated run-ups of three different rotor-bearing assembly models.Detailed descriptions of the model parameters are given in Appendix A.

Customized Laval/Jeffcott-Rotor in Plain Hydrodynamic Bearings
In [5], Schweizer presented a customized, extended version of the well-known Laval/Jeffcott-rotor [19], where the usually linear elastic and damped bearings were replaced by plain hydrodynamic bearings with circular geometry which introduce non-linear bearing forces.In this section such a rotor model is investigated using the simulation tool HOTINT [11].In a first step, for model validation a comparison between the results of Schweizer and the HOTINT implementation is drawn for circular-shaped bearings.Then, the model is extended by non-circular journal bearings in order to study the potential of preventing self-excited vibrations.Figure 4 shows a sketch of the investigated rotor system.

Customized Laval/Jeffcott-Rotor in Plain Hydrodynamic Bearings
In [5], Schweizer presented a customized, extended version of the well-known Laval/Jeffcott-rotor [19], where the usually linear elastic and damped bearings were replaced by plain hydrodynamic bearings with circular geometry which introduce non-linear bearing forces.In this section such a rotor model is investigated using the simulation tool HOTINT [11].In a first step, for model validation a comparison between the results of Schweizer and the HOTINT implementation is drawn for circular-shaped bearings.Then, the model is extended by non-circular journal bearings in order to study the potential of preventing self-excited vibrations.Figure 4 shows a sketch of the investigated rotor system.The simulated structure consists of an elastic Laval/Jeffcott rotor with a rigid disk mass centered on a massless shaft.Identical single hydrodynamic journal bearings on the left and right end support the assembly.A point mass  p is fixed to the disk at a given radial distance  p from the center of the shaft to account for the centrifugal imbalance force (Ω S ) =  P  P (2Ω S ()) 2 , where Ω S is the shaft's rotational speed and  is the time variable.To simulate a smooth run-up Ω S () = const.for  < 2 s and Figure 5 shows a comparison of results for the Laval/Jeffcott rotor model using HOTINT and the solution presented by Schweizer [5].The same parameter setup as in [5] was used, for details see the Appendix A. A spectrogram of the central node's  displacement (node 2), i.e., the disk displacement, is shown in Figure 5a.The abscissa corresponds to the vibration frequency, the ordinate depicts the time of the simulated run-up and the color displays the displacement value within one FFT-window.Clearly two kinds of displacements can be observed: (1) results from imbalance and (2) represents a self-excited, also called subsynchronous vibration, starting with Ω SW,Onset ≈ 400 s -1 , with an initial frequency of  SW,Onset ≈ 155 Hz ( SW,Onset = 4.5 s).The frequency of this phenomenon varies only slightly during the run-up, which is why the acoustical effect of this motion is often called "constant tone" in technical literature [7,[20][21][22][23][24][25][26][27][28].However, the term might be misleading as the frequency is not exactly constant over time.Since the addressed vibration is a consequence of a resonant shaft displacement, in the present work we refer to it as shaft whirl instead.The onset of the self-excited motion occurs slightly earlier in the authors' results than in the article by Schweizer.The difference can be accounted to of a cavitation model [29] in [5], which is not used in the HOTINT model.The difference in onset rotational speeds is also small.The simulated structure consists of an elastic Laval/Jeffcott rotor with a rigid disk mass centered on a massless shaft.Identical single hydrodynamic journal bearings on the left and right end support the assembly.A point mass m p is fixed to the disk at a given radial distance r p from the center of the shaft to account for the centrifugal imbalance force U(Ω S ) = m P r P (2πΩ S (t)) 2 , where Ω S is the shaft's rotational speed and t is the time variable.To simulate a smooth run-up Ω S (t) = const.for t < 2 s and Ω S (t) = const.+ A(1 − cos( t−2 s 10 s−2 s π)) for t ≥ 2 s is used.Figure 5 shows a comparison of results for the Laval/Jeffcott rotor model using HOTINT and the solution presented by Schweizer [5].The same parameter setup as in [5] was used, for details see the Appendix A. A spectrogram of the central node's y displacement (node 2), i.e., the disk displacement, is shown in Figure 5a.The abscissa corresponds to the vibration frequency, the ordinate depicts the time of the simulated run-up and the color displays the displacement value within one FFT-window.Clearly two kinds of displacements can be observed: (1) results from imbalance and (2) represents a self-excited, also called subsynchronous vibration, starting with Ω SW,Onset ≈ 400 s −1 , with an initial frequency of f SW,Onset ≈ 155 Hz (t SW,Onset = 4.5 s).The frequency of this phenomenon varies only slightly during the run-up, which is why the acoustical effect of this motion is often called "constant tone" in technical literature [7,[20][21][22][23][24][25][26][27][28].However, the term might be misleading as the frequency is not exactly constant over time.Since the addressed vibration is a consequence of a resonant shaft displacement, in the present work we refer to it as shaft whirl instead.The onset of the self-excited motion occurs slightly earlier in the authors' results than in the article by Schweizer.The difference can be accounted to of a cavitation model [29] in [5], which is not used in the HOTINT model.The difference in onset rotational speeds is also small.Figure 6 shows results for shaft displacements at the central node (Figure 6a) and the bearing eccentricities (Figure 6b).Again, the agreement is satisfactory.The good agreement between the obtained results and reference solution from literature with respect to frequency characteristics, displacement amplitudes, and bearing eccentricities confirms the validity of the implementation in HOTINT.
In view of practical applications, in some cases, self-excited motions of turbochargers can be audible and unwanted.Reduction methods mostly employ geometrical changes towards a certain asymmetry in the bearing design, which usually only results in a reduction of the amplitude or change of the onset speed of subsynchronous effects.In the following, we show that non-circular bearing shapes may allow a complete prevention of self-excited vibrations.To this end, a three-lobe sinusoidal clearance shape, as depicted in the cross-section in Figure 4, with a variation of 10% in the radial clearance with respect to the nominal clearance, ∆ = 0.1  NOM , is used instead of the circular bearing geometry.This modification is sufficient to completely suppress the occurrence of selfexcited motions as shown by the comparison of the results for circular vs. non-circular bearing geometries in Figure 7a,b.Further insights in the systems' characteristics may be obtained by a linearized eigenvalue analysis.Figure 7c,d show the frequencies   Im ,  ∈ {1,2,3,4}, of the lowest eigenmodes of a linearization of the system around a transient steady-state solution at a given constant rotational speed (see Equation ( 2)) for further details) after 0.5 s simulated time.Hence, we do not compute the linearization around a specifically chosen point.However, as can be seen in Figures 6 and 8, before the onset of the self-excited vibration, the orbit radius is small, and beyond the onset, the orbits in the bearings are almost circular such that the self-excited forces dominate the Figure 6 shows results for shaft displacements at the central node (Figure 6a) and the bearing eccentricities (Figure 6b).Again, the agreement is satisfactory.Figure 6 shows results for shaft displacements at the central node (Figure 6a) and the bearing eccentricities (Figure 6b).Again, the agreement is satisfactory.The good agreement between the obtained results and the reference solution from literature with respect to frequency characteristics, displacement amplitudes, and bearing eccentricities confirms the validity of the implementation in HOTINT.
In view of practical applications, in some cases, self-excited motions of turbochargers can be audible and unwanted.Reduction methods mostly employ geometrical changes towards a certain asymmetry in the bearing design, which usually only results in a reduction of the amplitude or change of the onset speed of subsynchronous effects.In the following, we show that non-circular bearing shapes may allow a complete prevention of self-excited vibrations.To this end, a three-lobe sinusoidal clearance shape, as depicted in the cross-section in Figure 4, with a variation of 10% in the radial clearance with respect to the nominal clearance, ∆ = 0.1  NOM , is used instead of the circular bearing geometry.This modification is sufficient to completely suppress the occurrence of selfexcited motions as shown by the comparison of the results for circular vs. non-circular bearing geometries in Figure 7a,b.Further insights in the systems' characteristics may be obtained by a linearized eigenvalue analysis.Figure 7c,d show the frequencies   Im ,  ∈ {1,2,3,4}, of the lowest eigenmodes of a linearization of the system around a transient steady-state solution at a given constant rotational speed (see Equation ( 2)) for further details) after 0.5 s simulated time.Hence, we do not compute the linearization around a specifically chosen point.However, as can be seen in Figures 6 and 8, before the onset of the self-excited vibration, the orbit radius is small, and beyond the onset, the orbits in the bearings are almost circular such that the self-excited forces dominate the The good agreement between the obtained results and the reference solution from literature with respect to frequency characteristics, displacement amplitudes, and bearing eccentricities confirms the validity of the implementation in HOTINT.
In view of practical applications, in some cases, self-excited motions of turbochargers can be audible and unwanted.Reduction methods mostly employ geometrical changes towards a certain asymmetry in the bearing design, which usually only results in a reduction of the amplitude or change of the onset speed of subsynchronous effects.In the following, we show that non-circular bearing shapes may allow a complete prevention of self-excited vibrations.To this end, a three-lobe sinusoidal clearance shape, as depicted in the cross-section in Figure 4, with a variation of 10% in the radial clearance with respect to the nominal clearance, ∆c = 0.1 c NOM , is used instead of the circular bearing geometry.This modification is sufficient to completely suppress the occurrence of self-excited motions as shown by the comparison of the results for circular vs. non-circular bearing geometries in Figure 7a,b.Further insights in the systems' characteristics may be obtained by a linearized eigenvalue analysis.Figure 7c,d show the frequencies λ Im k , k ∈ {1, 2, 3, 4}, of the lowest eigenmodes of a linearization of the system around a transient steady-state solution at a given constant rotational speed (see Equation ( 2)) for further details) after 0.5 s simulated time.Hence, we do not compute the linearization around a specifically chosen point.However, as can be seen in Figures 6 and 8, before the onset of the self-excited vibration, the orbit radius is small, and beyond the onset, the orbits in the bearings are almost circular such that the self-excited forces dominate the dynamic behavior.Therefore, the exact choice of the point of linearization is not crucial.The corresponding modes are a shaft bending mode at an almost constant frequency of around 240 Hz, a cylindrical rigid-body whirl mode at half the current rotational frequency of the shaft (i.e., where the effective hydrodynamic velocity amounts to zero, cf., for instance, [5,15,16]), and two coupled modes.The lowest coupled mode, which starts at zero at a shaft frequency of 0 Hz and locks in the bending mode at higher shaft speeds, corresponds to the self-excited subsynchronous motion as shown in Figure 7a.In the linearized analysis, the other modes show stronger damping (negative real parts λ Re k of the respective eigenvalues, cf.again Equation ( 2)) and thus are not excited during the dynamic run-up procedure.The lowest mode, in contrast, exhibits significantly smaller damping values and even instable behavior indicated by a positive real part of the complex eigenvalue in case of the circular bearing geometry, cf. Figure 7e.Note that the rotor speed at which this instability occurs at the first time corresponds well with the onset of the self-excited motion in the of the run-up simulation in Figure 7a.From this perspective, the main effect of the non-circular bearing geometry can be described as a shift of the damping of the instable mode to the stable region, which is clearly visible in Figure 7e.The eigenfrequencies, in contrary, do not change significantly, cf. Figure 7c,d.Physically, the additional damping, and thus, increased stability of the bearings with non-circular geometry is related to an additional pressure build-up in the oil film.It should be noted that the linearized eigenvalue analysis may only give indications on the actual run-up dynamics, however, it does not fully resemble the non-linearity of the system.Particularly, it does not include information about the energy transfer to (or from) self-excited vibrations which are characterized by the interplay between the dissipation and excitation.Stable limit cycle oscillations are defined by the equilibrium between transferred and dissipated energy (see, e.g., [30]).Having a closer look at the comparison between Figure 7a,e again, as discussed above, the onset of the subsynchronous motion corresponds well with the point where the corresponding eigenmode of the linearized systems becomes unstable (i.e., exhibits a positive real part of the eigenvalue).With increasing rotor speed, the eigenmode becomes stable again (negative real parts of the eigenvalue), which corresponds to a stable limit cycle oscillation where the energy transfer to the mode is compensated by its damping.Starting with a rotor speed of around 1500 s −1 , however, the amplitude of the self-excited oscillation starts to decrease, despite a decrease, in some points even negative damping of the respective eigenmode.These results lead to the conclusion that, in that range of rotor speeds, the energy transfer to the self-excited vibration gradually vanishes.
In literature, the abrupt change of displacement or eccentricities at t SW,Onset , c.f. Figure 6, is often described as a jump-like phenomenon of the shaft orbits to so-called limit cycles, see, e.g., [6,[31][32][33].Before the occurrence of the self-excited vibration, the shaft rotational axis performs small motions around a quasi-equilibrium state of the rotor at the particular rotational speed.At Ω SW,Onset the disk and bearing node orbits jump to the their maximum displacements (limit cycle) as indicated in the spectrograms and move around the bearing with an orbital frequency of f SW,Onset .The orbit sizes decrease for increasing shaft rotational speeds until the self-excited motion completely fades out.The shaft orbit within a circular-shaped bearing is depicted in Figure 8a demonstrating the sharp increase of the shaft displacements in the limit cycle beyond the onset speed.One orbit recirculation resulting from the self-excited vibration takes a time-period of 1/ f SW |Ω S =1250 s −1 and is modulated by the significantly smaller imbalance-excited motion with a period of 1/1250 s −1 .
the eigenmode becomes stable again (negative real parts of the eigenvalue), which corresponds to a stable limit cycle oscillation where the energy transfer to the mode is compensated by its damping.Starting with a rotor speed of around 1500 s −1 , however, the amplitude of the self-excited oscillation starts to decrease, despite a decrease, in some points even negative damping of the respective eigenmode.These results lead to the conclusion that, in that range of rotor speeds, the energy transfer to the self-excited vibration gradually vanishes.Figure 8b shows the bearing node shaft orbits in the non-circular bearings described above for ten shaft revolutions at speeds of 250 and 1250 s −1 in comparison with the maximum clearance in the bearings.The orbit size increases with rotational speed but is much smaller compared to a Laval/Jeffcott rotor with circular bearing shapes which illustrates the suppression of self-excited resonance phenomena.
To conclude, the validity of the employed simulation approach was shown by comparison with reference results from the literature.Furthermore, the potential of non-circularly shaped bearings for the prevention of self-excited shaft vibrations was investigated and demonstrated.In literature, the abrupt change of displacement or eccentricities at  SW,Onset , c.f. Figure 6, is often described as a jump-like phenomenon of the shaft orbits to so-called limit cycles, see, e.g., [6,[31][32][33].Before the occurrence of the self-excited vibration, the shaft rotational axis performs small motions around a quasi-equilibrium state of the rotor at the particular rotational speed.At Ω SW,Onset the disk and bearing node orbits jump to the their maximum displacements (limit cycle) as indicated in the spectrograms and move around the bearing with an orbital frequency of  SW,Onset .The orbit sizes decrease for increasing shaft rotational speeds until the self-excited motion completely fades out.The shaft orbit within a circular-shaped bearing is depicted in Figure 8a demonstrating the sharp increase of the shaft displacements in the limit cycle beyond the onset speed.One orbit recirculation resulting from the self-excited vibration takes a time-period of 1/ SW |Ω  =1250 s −1 and is modulated by the significantly smaller imbalance-excited motion with a period of 1/1250 s −1 .

Reduced Turbocharger Model with Floating-Bearing Rings
This section shows results of numerical rotor run-ups of a reduced turbocharger model presented in [10].There, simulated results for rotors with circular and lobe-shaped full-floating bearing rings are compared accordingly with experimental results from run-ups on a turbocharger balancing rig.Using the described validated simulation technique the reduced turbocharger model is reproduced and the results are compared with [10].
The reduced turbocharger model consists of two rigid disks representing the turbine and compressor wheel mounted to a rotating, elastic shaft modeled by twelve beam elements.Independently rotating floating rings support the structure at nodes 5 and 9, see also the sketch of the cross-section in Figure 9.The floating ring's inner and outer contours consist of six circumferentially in-phase sinusoidal lobes with an amplitude of 25% of the nominal radial clearance.
Imbalances with opposite orientation are added by mass points fixed to both disks.The oil inlets are represented as rectangular areas with constant oil pressure.For further details on the implementation and exact model data the reader is referred to Section 2 and Appendix A. Figure 8b shows the bearing node shaft orbits in the non-circular bearings described above for ten shaft revolutions at speeds of 250 and 1250 s −1 in comparison with the maximum clearance in the bearings.The orbit size increases with rotational speed but is much smaller compared to a Laval/Jeffcott rotor with circular bearing shapes which illustrates the suppression of self-excited resonance phenomena.
To conclude, the validity of the employed simulation approach was shown by comparison with reference results from the literature.Furthermore, the potential of non-circularly shaped bearings for the prevention of self-excited shaft vibrations was investigated and demonstrated.

Reduced Turbocharger Model with Floating-Bearing Rings
This section shows results of numerical rotor run-ups of a reduced turbocharger model presented in [10].There, simulated results for rotors with circular and lobe-shaped full-floating bearing rings are compared accordingly with experimental results from run-ups on a turbocharger balancing rig.Using the described validated simulation technique the reduced turbocharger model is reproduced and the results are compared with [10].
The reduced turbocharger model consists of two rigid disks representing the turbine and compressor wheel mounted to a rotating, elastic shaft modeled by twelve beam elements.Independently rotating floating rings support the structure at nodes 5 and 9, see also the sketch of the cross-section in Figure 9.The floating ring's inner and outer contours consist of six circumferentially in-phase sinusoidal lobes with an amplitude of 25% of the nominal radial clearance.
Imbalances with opposite orientation are added by mass points fixed to both disks.The oil inlets are represented as rectangular areas with constant oil pressure.For further details on the implementation and exact model data the reader is referred to Section 2 and Appendix A. Unlike as in [10] the hydrodynamic pressure communication in the oil connecting channels of the floating rings is not considered in the HOTINT model.However, as it will be shown, the effects on the bearing ring speeds are small.In both studies some simplifications, such as neglected eddy flow within the oil connecting channels as well as the assumption of constant oil viscosities and bearing clearances during the run-up, are introduced.Unlike as in [10] the hydrodynamic pressure communication in the oil connecting channels of the floating rings is not considered in the HOTINT model.However, as it will be shown, the effects on the bearing ring speeds are small.In both studies some simplifications, such as neglected eddy flow within the oil connecting as well as the assumption of constant oil viscosities and bearing clearances during the run-up, are introduced.
Calculated y displacements of the shaft in the turbine-sided bearing (node 5) are depicted for the floating ring in Figure 10a and the relative displacement between shaft and ring in Figure 10b, respectively.Compared to the analysis of the shaft displacement y S only, which is commonly done in the literature, the investigation of the absolute difference y FBR − y S between the y-displacements of the floating ring y FBR and the shaft makes the difference between instabilities in the inner and outer oil film, here denoted as ring whirl and shaft whirl, more obvious.Shaft whirls can be observed in both spectrograms of y FBR and y FBR − y S , whereas ring whirls only in the y FBR spectrogram.Note that in a y S spectrogram the ring whirl would be dominant, as discussed in the following.In general, corresponding displacements in Figure 10a The shaft whirl and the ring whirl become even more apparent in Figure 11 in the orbits of  FBR −  S and of  FBR .From Ω S,Start = 333 s −1 to Ω S,Bifurcation = 1470 s −1 (Figure 11a) both the bearing orbit and the relative shaft orbit are characterized by the shaft whirl.After the onset of the ring whirl the orbits are shaped by the resonant bearing displacements (Figure 11b).Both the bearing and the shaft displacements are now dominated by the ring whirl (blue solid line) such that the vector difference between the two is small (black solid line).The shaft starts to perform a self-excited motion beneath the minimum rotor speed.After a bifurcation at a shaft rotational speed of Ω S,Bifurcation = 1470 s −1 an additional self-excited motion comes into play.The sharp increase of the compressor ring speed indicates that the bifurcation can be assigned to a shaft whirl at the compressor side.Conversely, a higher shaft eccentricity leads to an increased acceleration of the floating ring.At a shaft speed of Ω S = 1685 s −1 the motion of the floating ring becomes resonant and a ring whirl occurs.Note that, since the floating rings carry the shaft, the ring whirl dominates the overall shaft motion (i.e., the displacement amplitudes).
The shaft whirl and the ring whirl become even more apparent in Figure 11 in the orbits of r FBR − r S and of r FBR .From Ω S,Start = 333 s −1 to Ω S,Bifurcation = 1470 s −1 (Figure 11a) both the bearing orbit and the relative shaft orbit are characterized by the shaft whirl.After the onset of the ring whirl the orbits are shaped by the resonant bearing displacements (Figure 11b).Both the bearing and the shaft displacements are now dominated by the ring whirl (blue solid line) such that the vector difference between the two is small (black solid line).
An introduction of a non-spherical bearing shape as depicted in the sketch of the cross-section of Figure 9 leads to a strikingly different characteristic of the shaft's and bearing rings' orbital motion.Figure 12 compares spectrograms of the vectorial sum of the bearing forces in the y direction for the turbocharger model with circular vs. non-circular bearing shapes.Furthermore, a comparison with spectrograms of accelerations obtained from measurements of the real turbocharger mounted on a balancing machine [10] is given.However, since the transfer function between bearing forces and accelerations is unknown, only the frequencies and relative values of the traces have been extracted.Additionally, [10] neither gives details about the exact rotor or bearing geometry nor about the mounting on the balancing machine.The agreement between our simulation results and the measurement data can therefore be seen as satisfactory, particularly the bifurcation of the shaft whirl corresponds well.It should be noted that the comparison between the simulation results and the measurement data needs to be taken care, since the model still contains simplifications (e.g., oil connecting channels in floating rings not modeled, constant temperature and viscosity assumed during run-up, housing and axial thrust not considered).Thus, the focus here is mainly on the correct representation of the qualitative dynamic behavior.The shaft whirl and the ring whirl become even more apparent in Figure 11 in the orbits of  FBR −  S and of  FBR .From Ω S,Start = 333 s −1 to Ω S,Bifurcation = 1470 s −1 (Figure 11a) both the bearing orbit and the relative shaft orbit are characterized by the shaft whirl.After the onset of the ring whirl the orbits are shaped by the resonant bearing displacements (Figure 11b).Both the bearing and the shaft displacements are now dominated by the ring whirl (blue solid line) such that the vector difference between the two is small (black solid line).An introduction of a non-spherical bearing shape as depicted in the sketch of the cross-section of Figure 9 leads to a strikingly different characteristic of the shaft's and bearing rings' orbital motion.Figure 12 compares spectrograms of the vectorial sum of the bearing forces in the  direction for the turbocharger model with circular vs. non-circular bearing shapes.Furthermore, a comparison with spectrograms of accelerations obtained from measurements of the real turbocharger mounted on a balancing machine [10] is given.However, since the transfer function between bearing forces and accelerations is unknown, only the frequencies and relative values of the traces have been extracted.Additionally, [10] neither gives details about the exact rotor or bearing geometry nor about the mounting on the balancing machine.The agreement between our simulation results and the measurement data can therefore be seen as satisfactory, particularly the bifurcation of the shaft whirl corresponds well.It should be noted that the comparison between the simulation results and the measurement data needs to be taken with care, since the model still contains simplifications (e.g., oil connecting channels in floating rings not modeled, constant temperature and viscosity assumed during run-up, housing and axial thrust not considered).Thus, the focus here is mainly on the correct representation of the qualitative dynamic behavior.The simulation of the turbocharger with non-circular bearings indicates that shaft whirls are no longer existent.The dominant ring whirl-with respect to the displacements-can be recognized by the deceleration of the turbine ring speed.It still occurs in case of non-circular bearings, however, the bifurcation and its onset shaft rotational speed is higher than in the turbocharger with circular shapes and its contribution to the bearing forces is small.A similar conclusion can be drawn from the experimental results of the balancing machine.Self-excited vibrations are still existent, but their The simulation of the turbocharger with non-circular bearings indicates that shaft whirls are no longer existent.The dominant ring whirl-with respect to the displacements-can be recognized by the deceleration of the turbine ring speed.It still occurs in case of non-circular bearings, however, the bifurcation and its onset shaft rotational speed is higher than in the turbocharger with circular shapes and its contribution to the bearing forces is small.A similar conclusion can be drawn from experimental results of the balancing machine.Self-excited vibrations are still existent, but their acceleration amplitudes are minor and perhaps result from production imperfection.
Comparing the simulated displacement spectrograms, Figure 10a,b, with the acceleration spectrograms, Figure 12b,d, it should be noted that the highest displacement value, the ring whirl, is responsible for the large orbit motion, the largest force and acceleration amplitudes, however, are found in the shaft whirls.For this reason, the ring whirl is not visible in the measured sensor acceleration.
In conclusion, in the simulation models as well as in agreement with measurements, non-circular bearing shapes have again proven to be an effective measure to (fully) suppress self-excited vibrations.

Conclusions and Outlook
An earlier published efficient hydrodynamic and rotordynamic modeling technique was used for the investigation of self-excited vibratory shaft and bearing ring motions of Laval/Jeffcott and turbocharger rotors.A detailed study of academic and real rotordynamic structures, validated by comparison to results from literature, was performed with special focus on the impact of non-spherical bearing shapes onto the shaft motion and the resulting bearing forces during rotor run-ups.Non-spherical bearing shapes such as lobes around the bearing's circumferential direction are shown to be promising as a means against self-excited shaft motions.In order to further qualify the simulation models as a tool for the early development stage of both Otto and Diesel combustion engines in parallel with their turbochargers, orbital measurements of shaft and bearing rings are important.A comparison of modeled and advanced measurement results is subject to future work-again with a focus on the effect and design of circular and non-circular bearing shapes (for a sketch see Figure 13, first simulation results are shown in Figure 14) with the perspective on the optimization of engine acoustics.
Lubricants 2017, 5, 6 15 of 19 found in the shaft whirls.For this reason, the ring whirl is not visible in the measured sensor acceleration.
In conclusion, in the simulation models as well as in agreement with measurements, non-circular bearing shapes have again proven to be an effective measure to (fully) suppress self-excited vibrations.

Conclusions and Outlook
An earlier published efficient hydrodynamic and rotordynamic modeling technique was used for the investigation of self-excited vibratory shaft and bearing ring motions of Laval/Jeffcott and turbocharger rotors.A detailed study of academic and real rotordynamic structures, validated by comparison to results from literature, was performed with special focus on the impact of non-spherical bearing shapes onto the shaft motion and the resulting bearing forces during rotor runups.Non-spherical bearing shapes such as lobes around the bearing's circumferential direction are shown to be promising as a means against self-excited shaft motions.In order to further qualify the simulation models as a tool for the early development stage of both Otto and Diesel combustion engines in parallel with their turbochargers, orbital measurements of shaft and bearing rings are important.A comparison of modeled and advanced measurement results is subject to future workagain with a focus on the effect and design of circular and non-circular bearing shapes (for a sketch see Figure 13, first simulation results are shown in Figure 14) with the perspective on the optimization of engine acoustics.Accelerations and kinematical data (e.g., floating ring motion and shaft orbits) will be obtained from experimental run-up studies with an advanced measurement setup on a turbocharger hot gas test rig.References [34,35] show that the influence of viscosity models can be large.In the future study, the effects of changing viscosity and clearance during the run-up will also be discussed.The combination of validated simulation models and state-of-the-art measurement techniques will lead to detailed insights into the rotordynamic characteristics of the investigated turbocharger model specifically, and the effect of non-circular bearing designs in general.Moreover, well-known results from the literature, such as distinct mode shapes of ring and shaft whirl (often referred to as conical vs. cylindrical motions, see, e.g., [3,20,23,[36][37][38][39]) will be studied based on measurement data.For the consideration of a non-spherical bearing geometry, the inner ring radius follows r B,in = r B,in,0 + 0.5[c FBR,in (1 − sin(3ϕ))].
The pressures at the axial boundaries of the hydrodynamic bearings is atmospheric pressure, p Left = p Right = 1 bar.
The oil inlet is a rectangular domain on top of the outside hydrodynamic FDM grid with a constant (absolute) pressure, p Inlet = 4.0 bar.It extends from ϕ sup,1 = −1.047rad to ϕ sup,end = 1.047 rad and from x sup,1 = −0.5 × 10 −3 m to x sup,end = 0.5 × 10 −3 m (cf.again the sketch of the domains in Figure 4) symmetrically around the center of the outer bearings.
A small numerical damping of d Turbine = d Compressor = 5 Ns•m −1 is added to both disks.The run-up is done numerically using a penalty constraint for the shaft's angular velocity with a stiffness of k Penalty = 1 × 10 4 Nms•rad −1 in order to prescribe the rotational speed between Ω S (t = 0) = 20 × 10 3 RPM to Ω S (t = 10s) = 130 × 10 3 RPM according to the formula The shaft's and the ring's axial motion is fixed at the first bearing node for numerical reasons.The tilting motion of the rings is also constrained.
Sensors measure the xand y-displacement of the turbine node, the compressor node, the bearing nodes, and accordingly the displacements of both floating rings.For the floating rings an inner and an outer eccentricity is evaluated.The static and dynamic eccentricity gives rise to inner and outer bearing forces and inner and outer bearing torques, which are observed as well.The rotational speed is stored for the shaft and for each ring individually.

Figure 2 .
Figure 2. Schematic of the computational domain of the pressure field (, , ).A rectangular zone for a pressure supply area is marked in orange.

Figure 3 .
Figure 3. Computed pressure field of a generalized shaft model supported in hydrodynamic radial bearings including a pressure supply area.Screenshots at different instants of time during a run-up simulation.

Figure 2 .
Figure 2. Schematic of the computational domain of the pressure field (, , ).A rectangular zone for a pressure supply area is marked in orange.

Figure 3 .
Figure 3. Computed pressure field of a generalized shaft model supported in hydrodynamic radial bearings including a pressure supply area.Screenshots at different instants of time during a run-up simulation.

Figure 3 .
Figure 3. Computed pressure field of a generalized shaft model supported in hydrodynamic radial bearings including a pressure supply area.Screenshots at different instants of time during a run-up simulation.

Figure 4 .
Figure 4. Model of a customized Laval/Jeffcott rotor using single hydrodynamic bearings.The solid and dash-dotted contours in the cross-sectional view on the right hand side correspond to a non-circular and circular shaped bearing, respectively.Finite element nodes are numbered 1-3.

Figure 4 .
Figure 4. Model of a customized Laval/Jeffcott rotor using single hydrodynamic bearings.The solid and dash-dotted contours in the cross-sectional view on the right hand side correspond to a non-circular and circular shaped bearing, respectively.Finite element nodes are numbered 1-3.

Figure 5 .
Figure 5.Comparison of results obtained with HOTINT and from [5].(a) Displacement spectrogram of the central node 2 in  direction; (b) Vibration frequencies from the spectrograms in [5] (cf.again Figure 4 for node numbers).Vibrations ( 1 ) and (2) result from imbalance-and self-excitation, respectively.

Figure 6 .
Figure 6.Comparison of results obtained with HOTINT (brown) against reference data from [5] (black dashed lines show envelopes).(a) Displacements of the central shaft node 2 in  direction; (b) Bearing eccentricities in the nodes 1 and 3 (cf.again Figure 4 for node numbers).

Figure 5 .
Figure 5.Comparison of results obtained with HOTINT and from [5].(a) Displacement spectrogram of the central node 2 in y direction; (b) Vibration frequencies from the spectrograms in [5] (cf.again Figure 4 for node numbers).Vibrations (1) and (2) result from imbalance-and self-excitation, respectively.

Figure 5 .
Figure 5.Comparison of results obtained with HOTINT and from [5].(a) Displacement spectrogram of the central node 2 in  direction; (b) Vibration frequencies from the spectrograms in [5] (cf.again Figure 4 for node numbers).Vibrations ( 1 ) and (2) result from imbalance-and self-excitation, respectively.

Figure 6 .
Figure 6.Comparison of results obtained with HOTINT (brown) against reference data from [5] (black dashed lines show envelopes).(a) Displacements of the central shaft node 2 in  direction; (b) Bearing eccentricities in the nodes 1 and 3 (cf.again Figure 4 for node numbers).

Figure 6 .
Figure 6.Comparison of results obtained with HOTINT (brown) against reference data from [5] (black dashed lines show envelopes).(a) Displacements of the central shaft node 2 in y direction; (b) Bearing eccentricities in the nodes 1 and 3 (cf.again Figure 4 for node numbers).

Figure 7 .
Figure 7. Spectrograms of the bearing nodes' (1 or 3 in Figure 4) -displacement of a run-up of a customized Laval/Jeffcott rotor in two identical single hydrodynamic bearings using (a) circular shaped and (b) non-circular shaped bearings; (1) and (2) indicate the imbalance-excited and the self-excited motions, respectively; (c,d) show the lowest four eigenfrequencies (first: black, second: blue, third: red, fourth: green) of the system linearized about a transient steady-state solution with a given constant rotational speed (see Equation (2)) for circular and non-circular bearing shapes; (e,f) show the damping parameter (again Equation (2)) of the first eigenmode (cf.black eigenfrequency in c,d)).

Figure 7 .
Figure 7. Spectrograms of the bearing nodes' (1 or 3 in Figure 4) y-displacement of a run-up of a customized Laval/Jeffcott rotor in two identical single hydrodynamic bearings using (a) circular shaped and (b) non-circular shaped bearings; (1) and (2) indicate the imbalance-excited and the self-excited motions, respectively; (c,d) show the lowest four eigenfrequencies (first: black, second: blue, third: red, fourth: green) of the system linearized about a transient steady-state solution with a given constant rotational speed (see Equation (2)) for circular and non-circular bearing shapes; (e,f) show the damping parameter (again Equation (2)) of the first eigenmode (cf.black eigenfrequency in c,d)).

Figure 8 .
Figure 8. Shaft motion at the central node (black, 2) and at the bearing nodes (blue, 1 or 3) during a run-up.The depicted orbits results from ten shaft revolutions at rotational speeds of 250 and 1250 s −1 .(a) Circular bearing shape; (b) non-circular bearing shape.

Figure 8 .
Figure 8. Shaft motion at the central node (black, 2) and at the bearing nodes (blue, 1 or 3) during a run-up.The depicted orbits results from ten shaft revolutions at rotational speeds of 250 and 1250 s −1 .(a) Circular bearing shape; (b) non-circular bearing shape.

Figure 9 .
Figure 9. Model of a reduced turbocharger model using floating bearing rings according to [10].Finite element nodes are numbered 1-13.

Figure 9 .
Figure 9. Model of a reduced turbocharger model using floating bearing rings according to [10].Finite element nodes are numbered 1-13 .

Figure 10 .
Figure 10.Results from simulated run-ups using the reduced turbocharger model in full-floating bearing rings.(a) Spectrogram of the turbine side floating ring displacement in  direction,  S ; (b) spectrogram of the shaft displacement in  direction at the turbine side bearing node relative to the bearing ring displacement,  FBR −  S ; The calculated turbine and compressor ring speeds are laid on top.

Figure 11 .
Figure 11.Reduced turbocharger model with circular bearing shape.Maximum inner and outer clearance and turbine side bearing ring orbit and shaft orbits at the turbine bearing node at shaft rotational speeds of (a) Ω S = 1250 s −1 ; (b) Ω S = 2000 s −1 .

Figure 10 .
Figure 10.Results from simulated run-ups using the reduced turbocharger model in full-floating bearing rings.(a) Spectrogram of the turbine side floating ring displacement in y direction, y S ; (b) spectrogram of the shaft displacement in y direction at the turbine side bearing node relative to the bearing ring displacement, y FBR − y S ; The calculated turbine and compressor ring speeds are laid on top.

Figure 10 .
Figure 10.Results from simulated run-ups using the reduced turbocharger model in full-floating bearing rings.(a) Spectrogram of the turbine side floating ring displacement in  direction,  S ; (b) spectrogram of the shaft displacement in  direction at the turbine side bearing node relative to the bearing ring displacement,  FBR −  S ; The calculated turbine and compressor ring speeds are laid on top.

Figure 11 .
Figure 11.Reduced turbocharger model with circular bearing shape.Maximum inner and outer clearance and turbine side bearing ring orbit and shaft orbits at the turbine bearing node at shaft rotational speeds of (a) Ω S = 1250 s −1 ; (b) Ω S = 2000 s −1 .

Figure 11 .
Figure 11.Reduced turbocharger model with circular bearing shape.Maximum inner and outer clearance and turbine side bearing ring orbit and shaft orbits at the turbine bearing node at shaft rotational speeds of (a) Ω S = 1250 s −1 ; (b) Ω S = 2000 s −1 .

Figure 12 .
Figure 12.Comparison of a reduced turbocharger model using floating bearing rings with (a) circular bearing shape from the present work; (b) circular bearing shape from data by Eling et al. [10]; (c) noncircular bearing shape from the present work; (d) non-circular bearing shape from data by Eling et al. [10].

Figure 12 .
Figure 12.Comparison of a reduced turbocharger model using floating bearing rings with (a) circular bearing shape from the present work; (b) circular bearing shape from data by Eling et al. [10]; (c) non-circular bearing shape from the present work; (d) non-circular bearing shape from data by Eling et al. [10].

Figure 13 .
Figure 13.Model of a turbocharger with floating-bearing rings.The sketch on the right-hand side shows a cross-section through a bearing which features an oil-inlet groove in the outer oil film and a non-circular geometry on the inside of the floating ring.

Figure 13 .
Figure 13.Model of a turbocharger with floating-bearing rings.The sketch on the right-hand side shows a cross-section through a bearing which features an oil-inlet groove in the outer oil film and a non-circular geometry on the inside of the floating ring.

Figure 13 .Figure 14 .
Figure 13.Model of a turbocharger with floating-bearing rings.The sketch on the right-hand side shows a cross-section through a bearing which features an oil-inlet groove in the outer oil film and a non-circular geometry on the inside of the floating ring.

Figure 14 .
Figure 14.Comparison of simulation results of a reduced turbocharger model using floating bearing rings with (a) circular inner bearing shape; (b) non-circular inner bearing shape.