Turbulence Theory Research Articles

MHD turbulence: a biased review

This review of scaling theories of magnetohydrodynamic (MHD) turbulence aims to put the developments of the last few years in the context of the canonical time line (from Kolmogorov to Iroshnikov–Kraichnan to Goldreich–Sridhar to Boldyrev). It is argued that Beresnyak's (valid) objection that Boldyrev's alignment theory, at least in its original form, violates the Reduced-MHD rescaling symmetry can be reconciled with alignment if the latter is understood as an intermittency effect. Boldyrev's scalings, a version of which is recovered in this interpretation, and the concept of dynamic alignment (equivalently, local 3D anisotropy) are thus an example of a physical theory of intermittency in a turbulent system. The emergence of aligned structures naturally brings into play reconnection physics and thus the theory of MHD turbulence becomes intertwined with the physics of tearing, current-sheet disruption and plasmoid formation. Recent work on these subjects by Loureiro, Malletet al.is reviewed and it is argued that we may, as a result, finally have a reasonably complete picture of the MHD turbulent cascade (forced, balanced, and in the presence of a strong mean field) all the way to the dissipation scale. This picture appears to reconcile Beresnyak's advocacy of the Kolmogorov scaling of the dissipation cutoff (as$\mathrm {Re}^{3/4}$) with Boldyrev's aligned cascade. It turns out also that these ideas open the door to some progress in understanding MHD turbulence without a mean field – MHD dynamo – whose saturated state is argued to be controlled by reconnection and to contain, at small scales, a tearing-mediated cascade similar to its strong-mean-field counterpart (this is a new result). On the margins of this core narrative, standard weak-MHD-turbulence theory is argued to require some adjustment – and a new scheme for such an adjustment is proposed – to take account of the determining part that a spontaneously emergent 2D condensate plays in mediating the Alfvén-wave cascade from a weakly interacting state to a strongly turbulent (critically balanced) one. This completes the picture of the MHD cascade at large scales. A number of outstanding issues are surveyed: imbalanced turbulence (for which a new, tentative theory is proposed), residual energy, MHD turbulence at subviscous scales, and decaying MHD turbulence (where there has been dramatic progress recently, and reconnection again turned out to feature prominently). Finally, it is argued that the natural direction of research is now away from the fluid MHD theory and into kinetic territory – and then, possibly, back again. The review lays no claim to objectivity or completeness, focusing on topics and views that the author finds most appealing at the present moment.

Energy dynamics in the helical wakes of flapping systems

Hovering insects are adapted for high metabolism. Their physiology, morphology, and flight have similar scaling relations. However, there are differences in efficiency. For instance, bumblebees have superior efficiency over orchid bees. Here, we propose that unsteady aerodynamics can explain these differences. Inspired by helical wakes in the bumblebee flight, we study an unsteady viscous model to show that there is intermittently abundant energy in the wake which is available for extraction by flapping wings. We then investigate the transport phenomena in the model and identify the role of kinematic viscosity, nonlinearity, flow topology, and fluid dynamic current to justify the available energy. Wake energetics is applicable to other species such as hummingbirds whose responses rely on environmental conditions. The present work has implications in minimizing power consumption in bio-mimetic locomotion both in air and under water. Interdisciplinary applications include dynamic modeling in superfluid turbulence and magneto-fluid dynamic dynamo theory.

Wave turbulence in inertial electron magnetohydrodynamics

A wave turbulence theory is developed for inertial electron magnetohydrodynamics (IEMHD) in the presence of a relatively strong and uniform external magnetic field $\boldsymbol {B_0} = B_0 \hat {\boldsymbol {e}}_\|$ . This regime is relevant for scales smaller than the electron inertial length $d_e$ . We derive the kinetic equations that describe the three-wave interactions between inertial whistler or kinetic Alfvén waves. We show that for both invariants, energy and momentum, the transfer is anisotropic (axisymmetric) with a direct cascade mainly in the direction perpendicular ( $\perp$ ) to $\boldsymbol {B_0}$ . The exact stationary solutions (Kolmogorov–Zakharov spectra) are obtained for which we prove the locality. We also found the Kolmogorov constant $C_K \simeq 8.474$ . In the simplest case, the study reveals an energy spectrum in $k_\perp ^{-5/2} k_\|^{-1/2}$ (with k the wavenumber) and a momentum spectrum enslaved to the energy dynamics in $k_\perp ^{-3/2} k_\|^{-1/2}$ . These solutions correspond to a magnetic energy spectrum ${\sim }k_\perp ^{-9/2}$ , which is steeper than the EMHD prediction made for scales larger than $d_e$ . We conclude with a discussion on the application of the theory to space plasmas.

Navigating the Covid 19 Turbulence in Higher Education: Evidence from Turkish Faculty Members

Covid19 was the first pandemic of the modern era to strike with such virulence. We sought to understand this recent phenomenon and contribute to the empirical findings on the expectations from HEI leadership and management in Turkey. Drawing on the Turbulence Theory, we explored how the academic staff experienced the initial phase of the pandemic in Turkey and how they perceived the HE leaders’ navigation of the crisis at the selected universities. Within qualitative phenomonology, data from semi-structured interviews with a convenient sample of 10 academic staff in five public and five private universities in Turkey, was analysed through content analysis. Findings highlighted the opportunities and challenges of the pandemic for the faculty at personal and organizational level in an intersectional pattern. Moreover, the ways HEI leaders navigated the crisis created binaries in the form of experience vs. inexperience and trust vs. distrust. The challenges derived from the rapid but ineffective decision-making processes and the heightened surveillance mechanisms over the academic staff; which in some cases resulted in lack of trust. Hence, the turbulence level was shaped by how the universities and their leaders addressed it. In such cases, practices of building trustworthy connections, more distributive forms of leadership and robust communication; which would help the leaders to navigate the turbulence at times of crises are significant. Further recommendations are provided for research, policy and practice.

The role played by relational turbulence in managing agency problems among value chain partners in the sharing economy: A review of the antecedents, benefits, risks, and boundary conditions

The purpose of this paper is to present a detailed theoretical review based on the “relational turbulence model,” revealing current research gaps and demonstrating its applicability to the management of agency problems among value chain partners by using antecedents, benefits, risks, and boundary conditions. “Relational turbulence theory” proposes a communication perspective on interpersonal relationships that explains how relationship transitions polarize emotions and cognitive assessments, disrupting partner communication in the sharing economy. Relational turbulence theory has recently gained traction in value chain literature in relation to the management of agency problems between partners within the sharing economy context; however, there is a surprising lack of agreement on the antecedents, benefits, risks, and boundary conditions of such relationships. To address this gap, we reviewed 52 peer-reviewed publications published between 2011 and 2021. Furthermore, this paper reveals the current research gaps and future research topics by recommending the application of relational turbulence theory with respect to antecedents, benefits, risks, and boundary conditions in the management of agency problems among value chain partners.

A Study of Large-eddy Simulation using Statistical and Machine Learning Techniques

The numerical solution of Navier-Stokes (N-S) equations has been found useful in various disciplines during its development, especially in recent years. However, a large-eddy simulation method has been developed to model the subgrid-scale dissipation rate by closing the Navier-Stokes equations. Because the instantaneous and time-averaged statistic characteristics of the subgrid-scale turbulent kinetic energy and dissipation have been studied by large eddy simulation. The purpose of this study is to check the statistical and machine learning of the subgrid-scale energy dissipation. As we know that the current turbulence theory states that the vortex stretching mechanism transports energy from large to small scales and leads to a high energy dissipation rate in a turbulent flow. Hence, a vortex-stretching-based subgrid-scale model is considered regarding the square of the velocity gradient to detect the playing role of the vortex stretching mechanism. The study in this article has shown a two-step process. Considering a posteriori statistic of the velocity gradient is analyzed through higher-order statistics and joint probability density function. Secondly, a machine learning approach is studied on the same data. The results of the vortex-stretching-based subgrid-scale model are then compared with the other two dynamic subgrid models, such as the localized dynamic kinetic energy equation model and the TKE-based Deardorff model. The results suggest that the vortex-stretching-based model can detect the significant subgrid-scale dissipation of small-scale motions and predict satisfactory turbulence statistics of the velocity gradient tensor.

Mode vortex and turbulence in ventilated cavitation over hydrofoils

In this paper, the complex interaction of gas-vortex-turbulence in ventilated cavitation over a NACA0015 hydrofoil is investigated by using large eddy simulation. The proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are used for analyzing the three-dimensional fluctuation modes of gas-liquid interfaces induced by turbulent coherent structures in ventilated cavitation. The results show that among the first four dominant POD modes possessing 99.5% modal kinetic energy in the ventilated cavitating flow, and the turbulent POD mode 2 contains the multi-scale coherent structures and acts as a bridge for the energy scale transition from mode 1 to mode 3 through the gas-liquid interface oscillations. The DMD modal structures are reconstructed from the roller-shaped energetic structure captured by POD mode 2, which further reveal the process of breaking down near the hydrofoil leading edge, changing into the hairpin shape and finally stretching into the horseshoe shape above the hydrofoil trailing edge for the turbulent coherent structures on the cavity interfaces. Such low-frequency wake structures are energetic but sustained by means of the high-frequency turbulent fluctuations, with the interfacial oscillation amplitude changed through the formation and transition of the multi-scale coherent structures. Further, the modal turbulent kinetic energy transport theory is introduced to describe the ventilated coherent structures. It is worth noting that the large-scale energetic structures in POD mode 2 are mainly drawing energy from mean flows through the Reynolds stress, and the modal turbulent energy stored by shedding vortex and air cavities in the ventilated wake flow is redistributed and balanced due to the turbulent and pressure diffusion effects. Besides, the uneven spatial distribution of turbulent energy caused by the vortex structures of different modal scales is proved to promote the spatial oscillation level of the ventilated cavity interfaces. This study is helpful to understand the interaction mechanism between turbulence and gas-liquid interfacial fluctuations for the ventilated wake flow.

How to Extract a Spectrum from Hydrodynamic Equations

The practical results gained from statistical theories of turbulence usually appear in the form of an inertial range energy spectrum $${\mathcal {E}}(k)\sim k^{-q}$$ and a cutoff wavenumber $$k_{c}$$ . For example, the values $$q=5/3$$ and $$\ell k_{c}\sim \mathrm {Re}^{3/4}$$ are intimately associated with Kolmogorov’s 1941 theory. To extract such spectral information from the Navier–Stokes equations, Doering and Gibbon (Phys. D 165, 163–175, 2020) introduced the idea of forming a set of dynamic wavenumbers $$\kappa _n(t)$$ from ratios of norms of solutions. The time averages of the $$\kappa _n(t)$$ can be interpreted as the 2nth moments of the energy spectrum. They found that $$1 < q \leqslant 8/3$$ , thereby confirming the earlier work of Sulem and Frisch (J. Fluid Mech. 72, 417–423, 1975) who showed that when spatial intermittency is included, no inertial range can exist in the limit of vanishing viscosity unless $$q \leqslant 8/3$$ . Since the $$\kappa _n(t)$$ are based on Navier–Stokes weak solutions, this approach connects empirical predictions of the energy spectrum with the mathematical analysis of the Navier–Stokes equations. This method is developed to show how it can be applied to many hydrodynamic models such as the two dimensional Navier–Stokes equations (in both the direct- and inverse-cascade regimes), the forced Burgers equation and shell models.

Physical regimes of electrostatic wave-wave nonlinear interactions generated by an electron beam propagating in a background plasma.

Electron-beam plasma interaction has long been a topic of great interest. Despite the success of the quasilinear and weak turbulence theories, their validities are limited by the requirements of a sufficiently dense mode spectrum and a small wave amplitude. In this paper, we extensively study the collective processes of a mono-energetic electron beam emitted from a thermionic cathode propagating through a cold plasma by performing high-resolution two-dimensional particle-in-cell simulations and using analytical theories. We confirm that, during the initial stage of two-stream instability between the beam and background cold electrons, it is saturated due to the well-known wave-trapping mechanism. Further evolution occurs due to strong wave-wave nonlinear processes. We show that the beam-plasma interaction can be classified into four different physical regimes in the parameter space for the plasma and beam parameters. The differences between the regimes are analyzed in detail. We identify a new regime in the strong Langmuir turbulence featured by what we call electron modulational instability (EMI) that could create a local Langmuir wave packet growing faster than the ion plasma frequency. Ions do not have time to respond to EMI in the initial growing stage. On a longer timescale, the action of the ponderomotive force produces very strong ion density perturbations, and eventually, the beam-plasma wave interaction stops being resonant due to the strong ion density perturbations. Consequently, in this EMI regime, electron beam-plasma interaction occurs in a repetitive (intermittent) process. The beam electrons are strongly scattered by waves, and the Langmuir wave spectrum is significantly broadened, which in turn gives rise to strong heating of bulk electrons. Associated energy transfer from the beam to the background plasma electrons has been studied. A resulting kappa (κ) distribution and a wave-energy spectrum E^{2}(k)∼k^{-5} are observed in the strong turbulent regime.

Path Large Deviations for the Kinetic Theory of Weak Turbulence

We consider a generic Hamiltonian system of nonlinear interacting waves with 3-wave interactions. In the kinetic regime of wave turbulence, which assumes weak nonlinearity and large system size, the relevant observable associated with the wave amplitude is the empirical spectral density that appears as the natural precursor of the spectral density, or spectrum, for finite system size. Following classical derivations of the Peierls equation for the moment generating function of the wave amplitudes in the kinetic regime, we propose a large deviation estimate for the dynamics of the empirical spectral density, where the number of admissible wavenumbers, which is proportional to the volume of the system, appears as the natural large deviation parameter. The large deviation stochastic Hamiltonian that quantifies the minus of the log probability of a trajectory is computed within the kinetic regime which assumes the Random Phase approximation for weak nonlinearity. We compare this Hamiltonian with the one for a system of modes interacting in a mean-field way with the empirical spectrum. Its relationship with the Random Phase and Amplitude approximation is discussed. Moreover, for the specific case when no forces and dissipation are present, a few fundamental properties of the large deviation dynamics are investigated. We show that the latter conserves total energy and momentum, as expected for a 3-wave interacting systems. In addition, we compute the equilibrium quasipotential and check that global detailed balance is satisfied at the large deviation level. Finally, we discuss briefly some physical applications of the theory.

Two-Dimensional Wavenumber Spectra on the Horizontal Submesoscale and Vertical Finescale

Abstract Horizontal and vertical wavenumbers (kx, kz) immediately below the Ozmidov wavenumber (N3/ε)1/2 are spectrally distinct from both isotropic turbulence (kx, kz &gt; 1 cpm) and internal waves as described by the Garrett–Munk (GM) model spectrum (kz &lt; 0.1 cpm). A towed CTD chain, augmented with concurrent Electromagnetic Autonomous Profiling Explorer (EM-APEX) profiling float microstructure measurements and shipboard ADCP surveys, are used to characterize 2D wavenumber (kx, kz) spectra of isopycnal slope, vertical strain, and isopycnal salinity gradient on horizontal wavelengths from 50 m to 250 km and vertical wavelengths of 2–48 m. For kz &lt; 0.1 cpm, 2D spectra of isopycnal slope and vertical strain resemble GM. Integrated over the other wavenumber, the isopycnal slope 1D kx spectrum exhibits a roughly +1/3 slope for kx &gt; 3 × 10−3 cpm, and the vertical strain 1D kz spectrum a −1 slope for kz &gt; 0.1 cpm, consistent with previous 1D measurements, numerical simulations, and anisotropic stratified turbulence theory. Isopycnal salinity gradient 1D kx spectra have a +1 slope for kx &gt; 2 × 10−3 cpm, consistent with nonlocal stirring. Turbulent diapycnal diffusivities inferred in the (i) internal wave subrange using a vertical strain-based finescale parameterization are consistent with those inferred from finescale horizonal wavenumber spectra of (ii) isopycnal slope and (iii) isopycnal salinity gradients using Batchelor model spectra. This suggests that horizontal submesoscale and vertical finescale subranges participate in bridging the forward cascade between weakly nonlinear internal waves and isotropic turbulence, as hypothesized by anisotropic turbulence theory.

Seasonality and Spatial Dependence of Mesoscale and Submesoscale Ocean Currents from Along-Track Satellite Altimetry

Abstract Along-track wavenumber spectral densities of sea surface height (SSH) are estimated from Jason-2 altimetry data as a function of spatial location and calendar month to understand the seasonality of meso- and submesoscale balanced dynamics across the global ocean. Regions with significant mode-1 and mode-2 baroclinic tides are rejected, restricting the analysis to the extratropics. Where balanced motion dominates, the SSH spectral density is averaged over all pass segments in a region for each calendar month and is fit to a four-parameter model consisting of a flat plateau at low wavenumbers, a transition at wavenumber k0 to a red power law spectrum k−s, and a white spectrum at high wavenumbers that models the altimeter noise. The monthly time series of the model parameters are compared to the evolution of the mixed layer. The annual mode of the spectral slope s reaches a minimum after the mixed layer deepens, and the annual mode of the bandpassed kinetic energy in the ranges [2k0, 4k0] and [k0, 2k0] peak ∼2 and ∼4 months, respectively, after the maximum of the annual mode of the mixed layer depth. This analysis is consistent with an energization of the submesoscale by a winter mixed layer instability followed by an inverse cascade of kinetic energy to the mesoscale, in agreement with prior modeling studies and in situ measurements. These results are compared to prior modeling, in situ, and satellite investigations of specific regions and are broadly consistent with them within measurement uncertainties. Significance Statement This paper uses satellite observations to understand the source of ocean dynamics at the 1–100-km scales at which vertical motion becomes important and which are thus relevant for biology and for the exchange of heat and carbon with the atmosphere. The observations are consistent with a seasonal variation of dynamics at these scales, predicted by a specific theory of upper-ocean turbulence and confirmed by modeling studies and regional observations. We update prior satellite-based studies by excluding regions with competing effects, by our treatment of the noise, and by our characterization of the seasonality. This work provides a template for analyzing data from the upcoming Surface Water and Ocean Topography (SWOT) satellite.

Spontaneously Stochastic Arnold’s Cat

We propose a simple model for the phenomenon of Eulerian spontaneous stochasticity in turbulence. This model is solved rigorously, proving that infinitesimal small-scale noise in otherwise a deterministic multi-scale system yields a large-scale stochastic process with Markovian properties. Our model shares intriguing properties with open problems of modern mathematical theory of turbulence, such as non-uniqueness of the inviscid limit, existence of wild weak solutions and explosive effect of random perturbations. Thereby, it proposes rigorous, often counterintuitive answers to these questions. Besides its theoretical value, our model opens new ways for the experimental verification of spontaneous stochasticity, and suggests new applications beyond fluid dynamics.

Study on the distribution law of airflow velocity in rectangular and semicircular roadway sections

AbstractReliable ventilation is the cornerstone of safe production in coal mines, and accurate monitoring of ventilation parameters is the fundamental guarantee of ventilation technology decision‐making. To improve the accuracy of airflow velocity monitoring in coal mine roadways, theoretical analysis, numerical simulation, and field tests were utilized to study the distribution law of airflow velocity in typical roadway sections. First, the calculation model of position of mean airflow velocity line in rectangular and semicircular arch roadway is established based on Boussinesq theory and Pelant turbulence theory. Then, to verify the correctness of the theoretical model, 25 groups of numerical simulation tests were conducted by using COMSOL‐Multiphysics 3.5a software. The errors between theoretical analysis and numerical simulation are all less than 4%. In addition, the numerical results also show that the contour line of airflow velocity in roadway section is consistent with the shape of roadway section, and the isoline of airflow velocity is basically parallel to the roadway wall. In addition, the closer to the roadway wall, the denser the airflow velocity isoline, indicating that the airflow velocity gradient near the wall is larger. And the thickness of the boundary layer decreases with the increase of airflow inlet velocity. Finally, field tests have been conducted in Chongqing Research Institute and Sima Coal Mine to further verify the correctness of the calculation model and numerical simulation results. The measured distribution law of airflow velocity is consistent with the numerical simulation. And the errors between theoretical analysis and field tests are all less than 4%.

Theoretical calculation of the laminar–turbulent transition in the round tube on the basis of stochastic theory of turbulence and equivalence of measures

Theoretical calculation of the laminar–turbulent transition in the round tube on the basis of stochastic theory of turbulence and equivalence of measures

Intermittency of gravity wave turbulence on the surface of an infinitely deep fluid: Numerical experiment

The intermittency of gravity wave turbulence on the surface of an infinitely deep fluid has been studied numerically. We estimated the scaling exponent of surface elevation structure-functions directly from the equations of motion, e.g., by numerical solution of the underlying equations of motion, without any additional closure assumptions (e.g., phase randomisation) and reduction to the kinetic equations as in the conventional framework of wave turbulence theory. With our high-resolution numerical methods for accurately modelling the nonlinear water surface, we were able to evaluate the scaling exponents of the structure-function up to relatively high orders (p=15), which has not been previously reported. We investigated the effect of wave steepness and forcing on intermittency and compared our results with other analytical and numerical studies, including results on the intermittency of wave turbulence of a different nature (turbulence of bending waves on a plate and magnetohydrodynamic turbulence). Our results support the conjecture of universality for intermittency phenomena in wave turbulence.

Swift generator for three-dimensional magnetohydrodynamic turbulence.

Magnetohydrodynamic turbulence is central to laboratory and astrophysical plasmas, and is invoked for interpreting many observed scalings. Verifying predicted scaling law behavior requires extreme-resolution direct numerical simulations (DNS), with needed computing resources excluding systematic parameter surveys. We here present an analytic generator of realistically looking turbulent magnetic fields, that computes three-dimensional (3D) O(1000^{3}) solenoidal vector fields in minutes to hours on desktop computers. Our model is inspired by recent developments in 3D incompressible fluid turbulence theory, where a Gaussian white noise vector subjected to a nonlinear transformation results in an intermittent, multifractal random field. Our B×C model has only few parameters that have clear geometric interpretations. We directly compare a (costly) DNS with a swiftly B×C-generated realization, in terms of its (1) characteristic sheetlike structures of current density, (2) volume-filling aspects across current intensity, (3) power-spectral behavior, (4) probability distribution functions of increments for magnetic field and current density, structure functions, and spectra of exponents, and (5) partial variance of increments. The model even allows to mimic time-evolving magnetic and current density distributions and can be used for synthetic observations on 3D turbulent data cubes.

Interplay of Thermalization and Strong Disorder: Wave Turbulence Theory, Numerical Simulations, and Experiments in Multimode Optical Fibers.

We address the problem of thermalization in the presence of a time-dependent disorder in the framework of the nonlinear Schrödinger (or Gross-Pitaevskii) equation with a random potential. The thermalization to the Rayleigh-Jeans distribution is driven by the nonlinearity. On the other hand, the structural disorder is responsible for a relaxation toward the homogeneous equilibrium distribution (particle equipartition), which thus inhibits thermalization (energy equipartition). On the basis of the wave turbulence theory, we derive a kinetic equation that accounts for the presence of strong disorder. The theory unveils the interplay of disorder and nonlinearity. It unexpectedly reveals that a nonequilibrium process of condensation and thermalization can take place in the regime where disorder effects dominate over nonlinear effects. We validate the theory by numerical simulations of the nonlinear Schrödinger equation and the derived kinetic equation, which are found in quantitative agreement without using any adjustable parameter. Experiments realized in multimode optical fibers with an applied external stress evidence the process of thermalization in the presence of strong disorder.

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

Performance of rough surface hydrodynamic circular and multi-lobe journal bearings in turbulent regimes

This article deals with the steady and dynamic performance analysis of rough surface circular and multi-lobe journal bearings operating with non-Newtonian lubricant in turbulent regimes. The numerical solution of Reynold's equation was obtained using the finite element method. To account for surface roughness, the Patir and Cheng model was used, and fluid turbulent behavior was modeled using the turbulent theory as proposed by Ng and Pan. The non-Newtonian behavior of the lubricant is presented by the Rabinowitsch fluid model, and JFO boundary conditions are used to solve Reynold's equation, considering gaseous cavitation. The bearing with micro-roughness has been considered in numerical analysis and transverse roughness has been reported to improve bearing performance indices. Among the different bearing design configurations, the two-lobe bearing configuration provides significantly improves the minimum film thickness, hydrodynamic pressure, and dynamic parameters. The presence of turbulent flow and transverse surface roughness produces a synergistic effect that enhances the minimum film thickness, stiffness, and damping parameters of journal bearings. This study provides a detailed comparison between various bearing designs and recommends the use of two-lobe bearings for high-speed applications.