scholarly journals On precessing flow in an oblate spheroid of arbitrary eccentricity

2014 ◽  
Vol 743 ◽  
pp. 358-384 ◽  
Author(s):  
Keke Zhang ◽  
Kit H. Chan ◽  
Xinhao Liao

AbstractWe consider a homogeneous fluid of viscosity $\nu $ confined within an oblate spheroidal cavity of arbitrary eccentricity $\mathcal{E}$ marked by the equatorial radius $d$ and the polar radius $d \sqrt{1-\mathcal{E}^2}$ with $0<\mathcal{E}<1$. The spheroidal container rotates rapidly with an angular velocity ${\boldsymbol{\Omega}}_0 $ about its symmetry axis and precesses slowly with an angular velocity ${\boldsymbol{\Omega}}_p$ about an axis that is fixed in space. It is through both topographical and viscous effects that the spheroidal container and the viscous fluid are coupled together, driving precessing flow against viscous dissipation. The precessionally driven flow is characterized by three dimensionless parameters: the shape parameter $\mathcal{E}$, the Ekman number ${\mathit{Ek}}=\nu /(d^2 \delimiter "026A30C {\boldsymbol{\Omega}}_0\delimiter "026A30C )$ and the Poincaré number ${\mathit{Po}}=\pm \delimiter "026A30C {\boldsymbol{\Omega}}_p\delimiter "026A30C / \delimiter "026A30C \boldsymbol{\Omega}_0\delimiter "026A30C $. We derive a time-dependent asymptotic solution for the weakly precessing flow in the mantle frame of reference satisfying the no-slip boundary condition and valid for a spheroidal cavity of arbitrary eccentricity at ${\mathit{Ek}}\ll 1$. No prior assumptions about the spatial–temporal structure of the precessing flow are made in the asymptotic analysis. We also carry out direct numerical simulation for both the weakly and the strongly precessing flow in the same frame of reference using a finite-element method that is particularly suitable for non-spherical geometry. A satisfactory agreement between the asymptotic solution and direct numerical simulation is achieved for sufficiently small Ekman and Poincaré numbers. When the nonlinear effect is weak with $\delimiter "026A30C {\mathit{Po}}\delimiter "026A30C \ll 1$, the precessing flow in an oblate spheroid is characterized by an azimuthally travelling wave without having a mean azimuthal flow. Stronger nonlinear effects with increasing $\delimiter "026A30C {\mathit{Po}}\delimiter "026A30C $ produce a large-amplitude, time-independent mean azimuthal flow that is always westward in the mantle frame of reference. Implications of the precessionally driven flow for the westward motion observed in the Earth’s fluid core are also discussed.

2012 ◽  
Vol 692 ◽  
pp. 420-445 ◽  
Author(s):  
Keke Zhang ◽  
Kit H. Chan ◽  
Xinhao Liao

AbstractWe consider a homogeneous fluid of viscosity $\nu $ confined within an oblate spheroidal cavity, ${x}^{2} / {a}^{2} + {y}^{2} / {a}^{2} + {z}^{2} / ({a}^{2} (1\ensuremath{-} {\mathscr{E}}^{2} ))= 1$, with eccentricity $0\lt \mathscr{E}\lt 1$. The spheroidal container rotates rapidly with an angular velocity ${\mbit{\Omega} }_{0} $, which is fixed in an inertial frame and defines a small Ekman number $E= \nu / ({a}^{2} \vert {\mbit{\Omega} }_{0} \vert )$, and undergoes weak latitudinal libration with frequency $\hat {\omega } \vert {\mbit{\Omega} }_{0} \vert $ and amplitude $\mathit{Po}\vert {\mbit{\Omega} }_{0} \vert $, where $\mathit{Po}$ is the Poincaré number quantifying the strength of Poincaré force resulting from latitudinal libration. We investigate, via both asymptotic and numerical analysis, fluid motion in the spheroidal cavity driven by latitudinal libration. When $\vert \hat {\omega } \ensuremath{-} 2/ (2\ensuremath{-} {\mathscr{E}}^{2} )\vert \gg O({E}^{1/ 2} )$, an asymptotic solution for $E\ll 1$ and $\mathit{Po}\ll 1$ in oblate spheroidal coordinates satisfying the no-slip boundary condition is derived for a spheroidal cavity of arbitrary eccentricity without making any prior assumptions about the spatial–temporal structure of the librating flow. In this case, the librationally driven flow is non-axisymmetric with amplitude $O(\mathit{Po})$, and the role of the viscous boundary layer is primarily passive such that the flow satisfies the no-slip boundary condition. When $\vert \hat {\omega } \ensuremath{-} 2/ (2\ensuremath{-} {\mathscr{E}}^{2} )\vert \ll O({E}^{1/ 2} )$, the librationally driven flow is also non-axisymmetric but latitudinal libration resonates with a spheroidal inertial mode that is in the form of an azimuthally travelling wave in the retrograde direction. The amplitude of the flow becomes $O(\mathit{Po}/ {E}^{1/ 2} )$ at $E\ll 1$ and the role of the viscous boundary layer becomes active in determining the key property of the flow. An asymptotic solution for $E\ll 1$ describing the librationally resonant flow is also derived for an oblate spheroidal cavity of arbitrary eccentricity. Three-dimensional direct numerical simulation in an oblate spheroidal cavity is performed to demonstrate that, in both the non-resonant and resonant cases, a satisfactory agreement is achieved between the asymptotic solution and numerical simulation at $E\ll 1$.


2018 ◽  
Author(s):  
Liqiang Dong ◽  
Wenxian Lin ◽  
Steven W. Armfield ◽  
Michael P. Kirkpatrick ◽  
Nicholas Williamson ◽  
...  

2001 ◽  
Vol 438 ◽  
pp. 67-100 ◽  
Author(s):  
N. A. PATANKAR ◽  
P. Y. HUANG ◽  
T. KO ◽  
D. D. JOSEPH

In this paper we study the lift-off to equilibrium of a single circular particle in Newtonian and viscoelastic fluids by direct numerical simulation. A particle heavier than the fluid is driven forward on the bottom of a channel by a plane Poiseuille flow. After a certain critical Reynolds number, the particle rises from the wall to an equilibrium height at which the buoyant weight just balances the upward thrust from the hydrodynamic force. The aim of the calculation is the determination of the critical lift-off condition and the evolution of the height, velocity and angular velocity of the particle as a function of the pressure gradient and material and geometric parameters. The critical Reynolds number for lift-off is found to be larger for a heavier particle whereas it is lower for a particle in a viscoelastic fluid. A correlation for the critical shear Reynolds number for lift-off is obtained. The equilibrium height increases with the Reynolds number, the fluid elasticity and the slip angular velocity of the particle. Simulations of single particle lift-off at higher Reynolds numbers in a Newtonian fluid by Choi & Joseph (2001) but reported here show multiple steady states and hysteresis loops. This is shown here to be due to the presence of two turning points of the equilibrium solution.


2012 ◽  
Vol 709 ◽  
pp. 610-621 ◽  
Author(s):  
Xinhao Liao ◽  
Keke Zhang

AbstractWe investigate, through both asymptotic and numerical analysis, precessionally driven flow of a homogeneous fluid confined in a fluid-filled circular cylinder that rotates rapidly about its symmetry axis and precesses slowly about a different axis that is fixed in space. After demonstrating that the inviscid approximation is always divergent even far away from resonance, we derive a general asymptotic solution for an asymptotically small Ekman number in the rotating frame of reference describing the weakly precessing flow that satisfies the no-slip boundary condition and that is valid at or near or away from resonance. Numerical analysis of the same problem using the Galerkin method in terms of a Chebyshev polynomial expansion is also carried out, showing satisfactory agreement between the general asymptotic solution and the corresponding numerical solution at or near or away from resonance.


Sign in / Sign up

Export Citation Format

Share Document