Global linear stability of the boundary-layer flow over a rotating sphere

We consider the linear global stability of the boundary-layer flow over a rotating sphere. Our results suggest that a self-excited linear global mode can exist when the sphere rotates sufficiently fast, with properties fixed by the flow at latitudes between approximately 55°-65° from the pole (depending on the rotation rate). A neutral curve for global linear instabilities is presented with critical Reynolds number consistent with existing experimentally measured values for the appearance of turbulence. The existence of an unstable linear global mode is in contrast to the literature on the rotating disk, where it is expected that nonlinearity is required to prompt the transition to turbulence. Despite both being susceptible to local absolute instabilities, we conclude that the transition mechanism for the rotating-sphere flow may be different to that for the rotating disk. © 2014 Elsevier Masson SAS. All rights reserved.

On the laminar-turbulent transition of the rotating-disk flow: The role of absolute instability

© 2014 Cambridge University Press. This paper describes a detailed experimental study using hot-wire anemometry of the laminar-turbulent transition region of a rotating-disk boundary-layer flow without any imposed excitation of the boundary layer. The measured data are separated into stationary and unsteady disturbance fields in order to elaborate on the roles that the stationary and the travelling modes have in the transition process. We show the onset of nonlinearity consistently at Reynolds numbers, R, of ∼ 510, i.e. at the onset of Lingwood's (J. Fluid Mech., vol. 299, 1995, pp. 17-33) local absolute instability, and the growth of stationary vortices saturates at a Reynolds number of ∼ 550. The nonlinear saturation and subsequent turbulent breakdown of individual stationary vortices independently of their amplitudes, which vary azimuthally, seem to be determined by well-defined Reynolds numbers. We identify unstable travelling disturbances in our power spectra, which continue to grow, saturating at around R=585, whereupon turbulent breakdown of the boundary layer ensues. The nonlinear saturation amplitude of the total disturbance field is approximately constant for all considered cases, i.e. different rotation rates and edge Reynolds numbers. We also identify a travelling secondary instability. Our results suggest that it is the travelling disturbances that are fundamentally important to the transition to turbulence for a clean disk, rather than the stationary vortices. Here, the results appear to show a primary nonlinear steep-fronted (travelling) global mode at the boundary between the local convectively and absolutely unstable regions, which develops nonlinearly interacting with the stationary vortices and which saturates and is unstable to a secondary instability. This leads to a rapid transition to turbulence outward of the primary front from approximately R=565 to 590 and to a fully turbulent boundary layer above 650.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.

Initiation of immersed granular avalanches.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.