Modelling and simulation of a sphere rolling on the inside of a rough vertical cylinder

A classical problem of nonholonomic system dynamics—the motion of a sphere on the inside of a rough vertical cylinder—is extended to rolling friction. The case study is modelled in independent coordinates. Due to the nonholonomic constraints imposed on the sphere, the governing equations arise as a set of differential-algebraic equations. The results of numerical simulations show the transition of the sphere from a sinusoid path on the vertical cylinder surface to a fall with slip. The physics of the ‘paradoxical’ motion is explained in detail.     

The boundary layer on a rotating sphere

Zusammenfassung Zur Berechnung der laminaren Grenzschicht auf einer rotierenden Kugel ist ein Differenzenverfahren entwickelt worden. Man kann sich dieses Verfahrens bedienen, um die Lösung schrittweise vom Pol zum Äquator zu erhalten. Die hier durchgeführten Berechnungen werden mit anderen Resultaten verglichen, welche man durch Reihenentwicklungen erhält.     

Note on the boundary layer on a rotating sphere

Zusammenfassung Die gleichmässige Strömung in der Grenzschicht bei einer gleichmässig rotierenden Kugel wird behandelt. Die Grenzschicht entsteht an den Polen und entwickelt sich in beiden Hemisphären nach dem Äquator zu. An den Polen verhält sich die Kugel wie eine rotierende Scheibe, und Flüssigkeit strömt hier ein. Am Äquator verhält sich die Kugel wie ein rotierender Zylinder, und die Flüssigkeit strömt in dessen Bereich aus.     

A note on a boundary-layer collision on a rotating sphere

A recent study of the unsteady boundary layer near the equator of a spinning sphere provided strong numerical evidence that the solution develops a singularity a finite time after the motion is initiated. In this paper we reformulate and complete the asymptotic structure proposed to describe this singularity. It is then self-consistent and provides convincing evidence of the accuracy of the computations and of the existence of the singularity.

---

Zusammenfassung Eine neuere Untersuchung über die nicht-stationäre Grenzschicht nahe am Äquator einer rotierenden Kugel ergab starke numerische Evidenz für das Auftreten einer Singularität in der Lösung nach einer endlichen Zeit, vom Beginn der Bewegung gerechnet. In dieser Arbeit wird die asymptotische Struktur zur Beschreibung der Singularität neu formuliert und ergänzt. Sie wird selbstkonsistent und gibt ein wichtiges Zeugnis für die Genauigkeit der Berechnungen und für die Existenz der Singularität.

     

Creeping flow of a sphere nearby a cylinder

We present a three-dimensional solution of a sphere nearby an infinite cylinder at low Reynolds number. We utilize the Lamb’s general solution based on spherical harmonics and develop a framework based on cylindrical harmonics to solve the flow field around the sphere and outside the cylinder, respectively. The solution is solved semi-analytically by considering geometrical parameters, including sphere radius, sphere velocity, separation distance and cylinder radius. The drag force coefficients of the sphere which are dependent on the distance between the cylinder surface and the sphere, as well as the velocity contours in the vicinity of the sphere, are analyzed. We also provide an analytical formula to calculate the drag force. The analytical formula has good quantitative agreement with the semi-analytical solution when the radius of the cylinder is smaller than the sphere. Such analysis can give insights into the details of the complex interaction between the sphere and cylinder.     

The thermal laminar boundary layer on a rotating sphere

Résumé Dans cet article, on étudie par un procédé numérique la couche limite laminaire thermique sur une sphère en rotation. En supposant que la sphère est maintenue à une température constante, on obtient la distribution des températures, ainsi que le transport de chaleur. On y discute aussi le problème de la sphère en rotation dans un cas plus général où la répartition de la température sur la surface de la sphère n'est plus uniforme.     

Latitudinally deforming rotating sphere

In this study, we consider a sphere with a surface that is fully covered by a stretchable elastic material. The radius of the sphere is fixed and it is also rotating about its radial axis. We investigate how the axisymmetric motion of a triggered fluid flow around the sphere is affected by the presence of both sphere rotation and latitudinal stretching. Considering that the deformation over the sphere commences at the pole, the problem is formulated such that the fluid flow near the pole is similar to the induced flow due to a linearly stretchable rotating disk, which has been described well in previous studies. When the rotation is omitted, the flow develops two-dimensionally under the action of pure stretching; otherwise, a three-dimensional axisymmetric fluid flow occurs, which is computed at each latitudinal angle both numerically and using a perturbation approach. The solution with wall deformation is different from the traditional character of the solution due to a solely rotating sphere. This solution is then used to compute the surface shears due to the physical drag and torque acting over the sphere. The contribution of wall stretching reduces the drag, whereas high rotation suppresses the effects of stretching to enhance the drag. More torque is required to rotate the sphere when both stretching and rotation mechanisms are in action.     

A note on Laplace’s equation inside a cylinder

Two difficulties connected with the solution of Laplace’s equation around an object inside an infinite circular cylinder are resolved. One difficulty is the non-convergence of Fourier transforms used, in earlier publications, to obtain the general solution, and the second difficulty concerns the existence of apparently different expressions for the solution. By using a Green’s function problem as an easily analyzed model problem, we show that, in general, Fourier transforms along the cylinder axis exist only in the sense of generalized functions, but when interpreted as such, they lead to correct solutions. We demonstrate the equivalence of the corrected solution to a different general solution, also previously published, but we point out that the two solutions have different numerical properties.     

A note on surface waves created by a rotating circular cylinder

A circular cylinder, placed vertically in a liquid with a free surface, makes oscillatory rotations about its axis. We find that small amplitude, radially progressive waves are created on the surface of the liquid that move a finite distance away from the cylinder.

---

Résumé Nous considérons les oscillations axiales d'un cylindre placé vérticalement dans une liquide à surface libre, et démontrons l'existence des ondes radiales à faible amplitude sur la surface s'eloignant à une distance finie du cylindre.

     

Optimal control of a cylinder rotating in a viscous fluid

Under consideration is the axisymmetric problem of optimal boundary control of a mechanical system consisting of two coaxial cylinders and an incompressible viscous fluid filling the region between them. The control parameter is the angular velocity of the outer cylinder. The goal is to stop the interior cylinder at a prescribed time with minimal energy expense. We prove that the problem is uniquely solvable and obtain the optimality system.     

The flow past a circular cylinder in a rotating frame

Summary In the present paper the flow past a circular cylinder in a uniformly rotating frame is investigated. A linear solution is given for the vertical shear layers surrounding the cylinder. The non-linear modifications of theE ^1/4 layer are considered and a criterion for separation obtained. The theoretical result is compared with experiment.

---

Résumé Cet étude s'occupe de l'écoulement d'un fluide au-dessus d'un cylindre circulaire dans une système qui tourne avec une vélocité uniforme. Une solution linéaire est presentée pour les couches verticales qui entourent le cylindre. Les modifications non-linéaires de la coucheE ^1/4 sont examinés et un critérium pour séparation est obtenu. Le résultat théorique est comparé avec l'experimentation.

     

The calculation of the thermal laminar boundary layer on a rotating sphere

Sommaire La communication expose une méthode permettant de déterminer la couche limite thermique d'une sphère en rotation dans un fluide incompressible à écoulement uniforme. La méthode se base sur l'idée deFrössling, selon laquelle la répartition des températures s'exprime en tant que développement en série de puissances, et utilise les résultats déjà obtenus parHoskin pour la répartition des vitesses dans la couche limite d'une sphère en rotation. L'auteur donne des résultats numériques obtenus sur une machine à calculer IBM 704.

---

Nomenclature x distance along surface of body in a meridian plane from forward stagnation point - y distance measured normal to surface - U velocity of main flow - U velocity of undisturbed stream - u component of velocity inx-direction - v component of velocity iny-direction - w transverse component of velocity due to spin - p pressure - density - T temperature - R distance of surface from axis of symmetry measured normal to axis of symmetry - coefficient of viscosity - kinematic viscosity - a radius of sphere - k thermal conductivity - angular velocity about axis of symmetry - Re Reynolds number - Pr Prandtl number - Nu Nusselt number - a constant - S area - Q quantity of heat transferred in unit time across areaS - d a characteristic length     

Approach to equilibrium of a rotating sphere in a Stokes flow

We study the unsteady rotary motion of a sphere immersed in a Stokes fluid. The equation of motion for the sphere leads to an integro-differential equation, and we are interested in the asymptotic behavior in time of the solution. Preparing initially the system (sphere + fluid) as a stationary state, we prove that the angular velocity of the sphere slows down with a law t ^−3/2 if no other forces than the one exerted by the fluid act on the sphere, while if the sphere is subject also to an elastic torque the asymptotic behavior of the angular position of the sphere is t ^−γ , with γ = 5/2 if the initial angular velocity is zero, γ = 3/2 otherwise. This behavior is due to the memory effect of the surrounding fluid. We discuss briefly other initial preparations of the system.     

Numerical studies of slow viscous rotating fluid past a sphere

Navier-Stokes equations for steady, viscous rotating fluid, rotating about the zaxis with angular velocity ω are linearized using Stokes approximation. The linearized Navier-Stokes equations governing the axisymmetric flow can be written as three coupled partial differential equations for the stream function, vorticity and rotational velocity component. Only one parameterR _eω =2ωa ²/v enters the resulting equations. Even the linearized equations are difficult to solve analytically and the method of matched asymptotic expansions is to be applied. Central differences are applied to the two-dimensional partial differential equations and are solved by the Peaceman-Rachford ADI method. The resulting algebraic equations are solved by successive over relaxation method. Streamlines are plotted for Ψ=0·01, 0·05, and 0·25 andR _eω =0·1, 0·3, 0·5.     

High Reynolds number turbulent flow past a rotating cylinder

High Reynolds-number flow over a rotating cylinder is investigated by two-dimensional numerical computations. The Reynolds-Averaged Navier-Stokes (RANS) equations are solved via the finite-volume method and they are closed by a modified k-ε turbulence model. The spin ratio a is defined as the ratio between the cylinder’s circumferential speed velocity to the free-stream varies from 2 to 8. The flow is examined at Reynolds numbers from 5 × 10⁵ to 5 × 10⁶, which is considered to be an interesting range for industrial flows. Available experimental and numerical data were used to verify the validity of the implemented procedure.The results revealed stabilization of the acting forces at high spin rates, thus indicating a flowfield with suppressed vortex-shedding activity, as it is expected, in accordance with theoretical considerations in previous studies. Load coefficients were found to be inversely proportional to the Reynolds number for most of the examined rotational rates.