Direct numerical simulations based on the incompressible nonlinear Navier–Stokes equations of the flow over the surface of a rotating disk have been conducted. An impulsive disturbance was introduced and its development as it travelled radially outwards and ultimately transitioned to turbulence has been analysed. Of particular interest was whether the nonlinear stability is related to the linear stability properties. Specifically three disk-edge conditions were considered; (i) a sponge region forcing the flow back to laminar flow, (ii) a disk edge, where the disk was assumed to be infinitely thin and (iii) a physically realistic disk edge of finite thickness. This work expands on the linear simulations presented by Appelquist et al. (J. Fluid. Mech., vol. 765, 2015, pp. 612–631), where, for case (i), this configuration was shown to be globally linearly unstable when the sponge region effectively models the influence of the turbulence on the flow field. In contrast, case (ii) was mentioned there to be linearly globally stable, and here, where nonlinearity is included, it is shown that both cases (ii) and (iii) are nonlinearly globally unstable. The simulations show that the flow can be globally linearly stable if the linear wavepacket has a positive front velocity. However, in the same flow field, a nonlinear global instability can emerge, which is shown to depend on the outer turbulent region generating a linear inward-travelling mode that sustains a transition front within the domain. The results show that the front position does not approach the critical Reynolds number for the local absolute instability, $R=507$. Instead, the front approaches $R=583$ and both the temporal frequency and spatial growth rate correspond to a global mode originating at this position.
Numerical simulations of the flow developing on the surface of a rotating disk are presented based on the linearized incompressible Navier-Stokes equations. The boundary-layer flow is perturbed by an impulsive disturbance within a linear global framework, and the effect of downstream turbulence is modelled by a damping region further downstream. In addition to the outward-travelling modes, inward-travelling disturbances excited at the radial end of the simulated linear region, r end , by the modelled turbulence are included within the simulations, potentially allowing absolute instability to develop. During early times the flow shows traditional convective behaviour, with the total energy slowly decaying in time. However, after the disturbances have reached r end , the energy evolution reaches a turning point and, if the location of r end is at a Reynolds number larger than approximately R = 594 (radius non-dimensionalized by √ ν/Ω * , where ν is the kinematic viscosity and Ω * is the rotation rate of the disk), there will be global temporal growth. The global frequency and mode shape are clearly imposed by the conditions at r end . Our results suggest that the linearized Ginzburg-Landau model by Healey (J. Fluid Mech., vol. 663, 2010, pp. 148-159) captures the (linear) physics of the developing rotating-disk flow, showing that there is linear global instability provided the Reynolds number of r end is sufficiently larger than the critical Reynolds number for the onset of absolute instability.
This paper proposes a resolution to the conundrum of the roles of convective and absolute instability in transition of the rotating-disk boundary layer. It also draws some comparison with swept-wing flows. Direct numerical simulations based on the incompressible Navier–Stokes equations of the flow over the surface of a rotating disk with modelled roughness elements are presented. The rotating-disk flow has been of particular interest for stability and transition research since the work by Lingwood (J. Fluid Mech., vol. 299, 1995, pp. 17–33) where an absolute instability was found. Here stationary disturbances develop from roughness elements on the disk and are followed from the linear stage, growing to saturation and finally transitioning to turbulence. Several simulations are presented with varying disturbance amplitudes. The lowest amplitude corresponds approximately to the experiment by Imayama et al. (J. Fluid Mech., vol. 745, 2014a, pp. 132–163). For all cases, the primary instability was found to be convectively unstable, and secondary modes were found to be triggered spontaneously while the flow was developing. The secondary modes further stayed within the domain, and an explanation for this is a proposed globally unstable secondary instability. For the low-amplitude roughness cases, the disturbances propagate beyond the threshold for secondary global instability before becoming turbulent, and for the high-amplitude roughness cases the transition scenario gives a turbulent flow directly at the critical Reynolds number for the secondary global instability. These results correspond to the theory of Pier (J. Engng Maths, vol. 57, 2007, pp. 237–251) predicting a secondary absolute instability. In our simulations, high temporal frequencies were found to grow with a large amplification rate where the secondary global instability occurred. For smaller radial positions, low-frequency secondary instabilities were observed, tripped by the global instability.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.