ALBERT

Description: Amplifications of flow past a backward-facing step with respect to optimal inflow and initial perturbations are investigated at Reynolds number 500. Two mechanisms of receptivity to inflow noise are identified: the bubble-induced inflectional point instability and the misalignment effect downstream of the secondary bubble. Further development of the misalignment results in decay of perturbations from x=28 onwards (the step is located at x=0), as has been observed in previous non-normality studies (Blackburn et al., J. Fluid Mech., vol. 603, 2008, pp. 271-304), and eventually limits the receptivity. The receptivity is found to be maximized at an inflow perturbation frequency of ω=0.50 and a spanwise wavenumber of β=0, where the inflow noise takes full advantage of both mechanisms and is amplified over two orders of magnitude in terms of the velocity magnitude. In direct numerical simulations (DNS) of the flow perturbed by optimal or random inflow noise, vortex shedding, flapping of bubbles, three-dimensionality and turbulence are observed in succession as the magnitude of the inflow noise increases. Similar features of linear and nonlinear receptivity are observed at higher Reynolds numbers. The Strouhal number of the bubble flapping is 0.08, at which the receptivity to inflow noise reaches a maximum. This Strouhal number is close to reported values extracted from DNS or large eddy simulations (LES) at larger Reynolds numbers (Le et al., J. Fluid Mech., vol. 330, 1997, pp. 349-374; Kaiktsis et al., J. Fluid Mech., vol. 321, 1996, pp. 157-187; Métais, New Trends in Turbulence, 2001, Springer; Wee et al., Phys. Fluids, vol. 16, 2004, pp. 3361-3373). Methods to further clarify the mechanisms of receptivity and to suppress the noise amplifications by modifying the base flow using a linearly optimal body force are proposed. It is observed that the mechanisms of optimal noise amplification are fully revealed by the distribution of the base flow modification, which weakens the bubble instabilities and misalignment effects and subsequently reduces the receptivity significantly. Comparing the base flow modifications with respect to amplifications of inflow and initial perturbations, it is found that the maximum receptivity to initial perturbations is highly correlated with the receptivity to inflow noise at the optimal frequency ω=0.50, and the correlation reduces as the inflow frequency deviates from this optimal value. © 2015 Cambridge University Press.

Description: Linear and nonlinear transient growths of perturbations on a vortex ring up to Reynolds number ( circulation/viscosity) are studied. For short time intervals, perturbations around the ring axis undergo the strongest linear transient growth and lead to secondary structures in the form of ringlets, owing to the Orr mechanism and an inviscid vorticity-amplification mechanism: in contrast to the well-reported instabilities and lobe structures along the vortex ring core. These secondary ringlet structures induce a tertiary group of ringlets through similar transient perturbation growth. This cascade of ringlets lead to the breakup of the main ring even before activation of the vortex-core instabilities. Such a cascade scenario is also observed in the development of a vortex ring perturbed by random disturbance in the axis region. These new modes and mechanisms for the generation and breakup of vortex ring structures bring insights into the dynamics and control of vortex ring flows. © 2017 Cambridge University Press.

Description: Flow past a NACA 65 blade at chord-based Reynolds number 138 500 is studied using stability analysis, generalized (spatially weighted) transient growth analysis and direct numerical simulations (DNS). The mechanisms of transition on various sections of the blade observed in previous work by Zaki et al. (J. Fluid Mech., vol. 665, 2010, pp. 57-98) are examined, with a focus on the pressure side around the leading edge. In this region, the linearly most energetic perturbation has spanwise wavenumber 40π (five boundary-layer thicknesses) and is tilted against the mean shear to take advantage of the Orr mechanism. In a DNS, the nonlinear development of this optimal perturbation induces Λ structures, which are further stretched to hairpin vortices before breaking down to turbulence. At higher spanwise wavenumber, e.g. 120π, a free-stream optimal perturbation is obtained upstream of the leading edge, in the form of streamwise vortices. During its nonlinear evolution, this optimal perturbation tilts the mean shear and generates spanwise periodic high- and low-speed streaks. Then through a nonlinear lift-up mechanism, the low-speed streaks are lifted above the high-speed ones. This layout of streaks generates a mean shear with two inflectional points and activates secondary instabilities, namely inner and outer instabilities previously reported in the literature. © 2017 Cambridge University Press.

Description: Transient energy growth of disturbances to co-rotating pairs of vortices with axial core flows is investigated in an analysis where vortex core expansion and vortex merging are included by adopting a time-evolving base flow. The dynamics of pairs are compared with those of individual vortices in order to highlight the effect of vortex interaction. Three typical vortex pair cases are studied, with the pairs comprised respectively of individually inviscidly unstable vortices at the streamwise wavenumber that maximizes the individual instabilities, viscously unstable vortices also at the streamwise wavenumber maximizing the individual instabilities and asymptotically stable vortices at streamwise wavenumber zero. For the inviscidly unstable case, the optimal perturbation takes the form of a superposition of two individual helical unstable modes and the optimal energy growth is similar to that predicted for an individual inviscid unstable vortex, while where the individual vortices are viscously unstable, the optimal disturbances within each core have similar spatial distributions to the individually stable case. For both of these cases, time horizons considered are much lower than those required for the merger of the undisturbed vortices. However, for the asymptotically stable case, large linear transient energy growth of optimal perturbations occurs for time horizons corresponding to vortex merging. Linear transient disturbance energy growth exhibited by pairs in this stable case is two to three orders of magnitude larger than that for a corresponding individual vortex. The superposition of the perturbation and the base flow shows that the perturbation has a displacement effect on the vortices in the base flow. Direct numerical simulations of stable pairs seeded by optimal initial perturbations have been carried out and acceleration/delay of vortex merging associated with a dual vortex meandering and vortex breakup related to axially periodic acceleration and delay of vortex merging are observed. For axially invariant cases, the sign of perturbation has an effect, as well as magnitude; the sign dependence relates to whether or not the perturbation adds to or subtracts from the swirl of the base flow. For a two-dimensional perturbation that adds to the swirl of the base flow, seeding with the linear optimal disturbance at a relative energy level 1 × 10 -4 induces the pair to move towards each other and approximately halves the time required for merger. Direct numerical simulation shows that the optimal three-dimensional perturbation can induce the vortex system to break up before merging occurs, since the two-dimensional nature of vortex merging is broken by the development of axially periodic perturbations. © 2012 Cambridge University Press.

Description: We determine optimal inflow boundary perturbations to steady flow through a straight inflexible tube with a smooth axisymmetric stenosis at a bulk-flow Reynolds number Re= 400, for which the flow is asymptotically stable. The perturbations computed produce an optimal gain, i.e. kinetic energy in the domain at a given time horizon normalized by a measure of time-integrated energy on the inflow boundary segment. We demonstrate that similarly to the optimal initial condition problem, the gain can be interpreted as the leading singular value of the forward linearized operator that evolves the boundary conditions to the final state at a fixed time. In this investigation we restrict our attention to problems where the temporal profile of the perturbations examined is a product of a Gaussian bell and a sinusoid, whose frequency is selected to excite axial wavelengths similar to those of the optimal initial perturbations in the same geometry. Comparison of the final state induced by the optimal boundary perturbation with that induced by the optimal initial condition demonstrates a close agreement for the selected problem. Previous works dealing with optimal boundary perturbation considered a prescribed spatial structure and computed an optimal temporal variation of a wall-normal velocity component, whereas in this paper we consider the problem of a prescribed temporal structure and compute the optimal spatial variation of velocity boundary conditions over a one-dimensional inflow boundary segment. The methodology is capable of optimizing boundary perturbations in general non-parallel two-and three-dimensional flows. © 2011 Cambridge University Press.

Description: The spectrum of the Batchelor vortex can be broadly split into a discrete spectrum, a potential spectrum and a free-stream spectrum where, since the last two spectra are both continuous, they can also be considered as one continuous spectrum. The discrete spectrum has been extensively studied but the continuous spectrum has received limited attention in the context of vortex flow. A local transient growth study is conducted and the contribution of the discrete spectrum and the continuous spectrum to the transient growth is separated by constructing optimal perturbations on the discrete or continuous sub-eigenspaces separately. It is found that the significant transient growth is mainly due to the non-normality of the continuous eigenmodes/spectrum whilst the discrete eigenmodes/spectrum have little contribution to the transient energy growth. A matrix-free method, which reduces to the local analysis when appropriate periodic boundary conditions are imposed, is also applied to investigate the transient growth in both a plane of constant azimuthal angle and a plane constant axial location. Previously studies by other authors have demonstrated that at zero azimuthal wavenumber the transient growth reaches infinitely large values over infinite time intervals while the optimal perturbations are located far from the vortex core. Therefore we limited our scope to small values of the time horizon so as to obtain reasonably strong transient effects stemming from physically relevant optimal perturbations. Two mechanisms of transient growth are observed: namely a redistribution of the azimuthal velocity to the azimuthal vorticity and interaction between out-of-vortex-core structures with those within the vortex core. A direct numerical simulation (DNS) of the vortex perturbed by optimal perturbations is conducted to investigate the nonlinear development of the optimal perturbations. In the azimuthally constant decomposed case, it is found that the optimal perturbation induces a string of bubble structures to be generated as a consequence of the non-orthogonality of continuous eigenmodes and the breakdown bubble is induced by viscous diffusion, while in the axially constant decomposition transient growth analysis, it is observed that the optimal perturbations associated with the continuous eigenmodes drive the vortex to vibrate around the initial vortex centre before eventually returning to its original position at larger times. This transient effect provides a mechanism for the 'vortex meandering' observed in previous experimental and numerical studies. These optimal perturbations associated with the continuous spectrum with out-of-vortex-core structures are observed to be activated by anisotropic inflow perturbations in the potential region. © 2012 Cambridge University Press.

Description: A novel algorithm is developed to calculate the nonlinear optimal boundary perturbations in three-dimensional incompressible flow. An optimal step length in the optimization loop is calculated without any additional calls to the Navier-Stokes equations. The algorithm is applied to compute the optimal inflow eddies for the flow around a wind turbine to clarify the mechanisms behind wake meandering, a phenomenon usually observed in wind farms. The turbine is modelled as an actuator disc using an immersed boundary method with the loading prescribed as a body force. At Reynolds number (based on free-stream velocity and turbine radius), the most energetic inflow perturbation has a frequency-2, and is in the form of an azimuthal wave with wavenumber and the same radius as the actuator disc. The inflow perturbation is amplified by the strong shear downstream of the edge of the disc and then tilts the rolling-up vortex rings to induce wake meandering. This mechanism is verified by studying randomly perturbed flow at. At five turbine diameters downstream of the disc, the axial velocity oscillates at a magnitude of more than 60Â % of the free-stream velocity when the magnitude of the inflow perturbation is 6Â % of the free-stream wind speed. The dominant Strouhal number of the wake oscillation is 0.16 at and keeps approximately constant at higher. This Strouhal number agrees well with previous experimental findings. Overall the observations indicate that the well-observed stochastic wake meandering phenomenon appearing far downstream of wind turbines is induced by large-scale (the same order as the turbine rotor) and low-frequency free-stream eddies. © 2018 Cambridge University Press.

Description: This study is focused on two- and three-dimensional incompressible flow past a circular cylinder for Reynolds number Re ≤ 1000. To gain insight into the mechanisms underlying the suppression of unsteadiness for this flow we determine the nonlinear optimal open-loop control driven by surface-normal wall transpiration. The spanwise-constant wall transpiration is allowed to oscillate in time, although steady forcing is determined to be most effective. At low levels of control cost, defined as the square integration of the control, the sensitivity of unsteadiness with respect to wall transpiration is a good approximation of the optimal control. The distribution of this sensitivity suggests that the optimal control at small magnitude is achieved by applying suction upstream of the upper and lower separation points and blowing at the trailing edge. At high levels of wall transpiration, the assumptions underlying the linearized sensitivity calculation become invalid since the base flow is eventually altered by the size of the control forcing. The large-magnitude optimal control is observed to spread downstream of the separation point and draw the shear layer separation towards the rear of the cylinder through suction, while blowing along the centreline eliminates the recirculation bubble in the wake. We further demonstrate that it is possible to completely suppress vortex shedding in two- and three-dimensional flow past a circular cylinder up to Re 1000, accompanied by 70% drag reduction when a nonlinear optimal control of moderate magnitude (with root-mean-square value 8 % of the free-stream velocity) is applied. This is confirmed through linearized stability analysis about the steady-state solution when the nonlinear optimal wall transpiration is applied. While continuously distributed wall transpiration is not physically realizable, the study highlights localized regions where discrete control strategies could be further developed. It also highlights the appropriate range of application of linear and nonlinear optimal control to this type of flow problem. © 2015 Cambridge University Press.

Description: SUMMARYThe production of winter wheat (Triticum aestivum L.) is affected by crop population structures and field microclimates. This 3-year study assessed the effect of different precision planting patterns and irrigation conditions on relative humidity (RH), air and soil temperature within the canopy, intercepted photosynthetically active radiation (iPAR), evapotranspiration (ET), water productivity (WP) and grain yields. Field experiments were conducted from 2011 to 2014 on a two-factor split-plot design with three replicates. The experiments involved three precision planting patterns (single row, alternating single and twin rows [hereafter ‘single–twin’] and twin row) and three irrigation treatments (0 mm (I0), 90 mm (I90) and 180 mm (I180)). Planting patterns and irrigation treatments exerted a significant effect on RH, air and soil temperature, iPAR, ET, WP and grain yield. The lowest RH and iPAR levels were detected in the single row pattern. When the irrigation treatment was identical, the highest soil and air temperatures were detected in the single row pattern, followed by the single–twin row and twin row patterns. Compared with the single row, the single–twin and twin row patterns increased ET by 0·3 and 1·4, WP by 4·7 and 5·7% and yields by 6·0 and 7·9%, respectively. Compared with I0, the I90 and I180 irrigation treatments increased ET by 0·3 and 1·4%, and WP by 4·7 and 5·7%, respectively. The grain yields of the twin row pattern were 5·8 and 1·7% higher than those of the single row and single–twin row patterns, respectively. Compared with I0, I90 increased yield by 19·3%. The twin row pattern improved crop structure and farmland microclimate by increasing RH and iPAR, and reducing soil and air temperatures, thus increasing grain yield. These results indicated that a twin row pattern effectively improved grain yield at I0. On the basis of iPAR, WP and grain yield, it was concluded that a twin row pattern combined with an I90 irrigation treatment provided optimal cropping conditions for the North China plain.