Flow and Heat Transfer Mechanisms in a Rotating Compressor Cavity Under Centrifugal Buoyancy-Driven Convection

This paper presents a systematic study of flow and heat transfer mechanisms in a compressor disc cavity with an axial throughflow under centrifugal buoyancy-driven convection, comparing with previously published experimental data. Wall-modelled large-eddy simulations are conducted for six operating conditions, covering a range of rotational Reynolds number (3.2x10^5 - 2.2x10^6), buoyancy parameter (0.11 - 0.26) and Rossby number (0.4 - 0.8). Numerical accuracy and computational efficiency of the simulations are considered. Wall heat transfer predictions are compared with measured data with a good level of agreement. A constant rothalpy core occurs at high Eckert number, appearing to reduce the driving buoyancy force. The flow in the cavity is turbulent with unsteady laminar Ekman layers observed on both discs except in the bore flow affected region on the downstream disc cob. The shroud heat transfer Nusselt number-Rayleigh number scaling agrees with that of natural convection under gravity for high Rayleigh numbers. Disc heat transfer is dominated by conduction across unsteady Ekman layers, except on the downstream disc cob. The disc bore heat transfer is close to a pipe flow forced convection correlation. The unsteady flow structure is investigated showing strong unsteadiness in the cavity that extends into the axial throughflow.

Flow and Heat Transfer Mechanisms in a Rotating Compressor Cavity Under Centrifugal Buoyancy-Driven Convection

This paper presents a systematic study of flow and heat transfer mechanisms in a compressor disc cavity with an axial throughflow under centrifugal buoyancy-driven convection, comparing with previously published experimental data. Wall-modelled large-eddy simulations are conducted for six operating conditions, covering a range of rotational Reynolds number (3.2 × 105 – 2.2 × 106), buoyancy parameter (0.11 – 0.26) and Rossby number (0.4 – 0.8). Numerical accuracy and computational efficiency of the simulations are considered. Wall heat transfer predictions are compared with measured data with a good level of agreement. A constant rothalpy core occurs at high Eckert number, appearing to reduce the driving buoyancy force. The flow in the cavity is turbulent with unsteady laminar Ekman layers observed on both discs except in the bore flow affected region on the downstream disc cob. The shroud heat transfer Nusselt number-Rayleigh number scaling agrees with that of natural convection under gravity for high Rayleigh numbers. Disc heat transfer is dominated by conduction across unsteady Ekman layers, except on the downstream disc cob. The disc bore heat transfer is close to a pipe flow forced convection correlation. The unsteady flow structure is investigated showing strong unsteadiness in the cavity that extends into the axial throughflow.

Influence of coastal upwelling on the surface current in the Bay of Bengal using HF radar and satellite observations

<p>Coastal Upwelling is a phenomenon in which cold and nutrient-enriched water from the Ekman layers reaches the surface enhancing the biological productivity of the upwelling region. In this work, an attempt is made to understand the influence of coastal upwelling on surface current variations during May 2018 to August 2018, when HF radar current observation (source: NIOT, India) is available. The wind-based Upwelling Index(UI<sub>wind</sub>) showed coastal upwelling throughout the study period. But the SST based upwelling index (UI<sub>sst</sub>) showed upwelling occurred only from May to the first week of June. Cross-shore components of HF radar-derived ocean surface current (CSSC)  showed strong similarity with UI<sub>sst</sub>. The first phase of upwelling from UI<sub>sst</sub> is observed to start on 5<sup>th</sup> May and lasts till 14<sup>th</sup> May with a maximum peak on around 10<sup>th</sup> May and having a horizontal extension of ~40 km. Then, there is a break period for about three days and after that, the second phase of upwelling starts on 17<sup>th</sup> May and lasts till 25<sup>th</sup> May with a maximum peak on around 20<sup>th</sup> May, but this time the horizontal extension is ~100 km which is much larger than during the first phase. A strong positive (from coast to offshore) CSSC is observed to start on around 5<sup>th</sup> May and lasts till 13<sup>th</sup> May with a maximum peak on around 10<sup>th</sup> May and having a horizontal extension of ~40 km, as observed from UIsst. A reversal of CSSC (towards coast) is noted on 14<sup>th</sup> May when the break of coastal upwelling is evident from UI<sub>sst</sub>. The CSSC then again started intensifying 15<sup>th</sup> May onwards and continued for ten days till 25<sup>th</sup> May, similar to UI<sub>sst</sub>.  The horizontal extension of the upwelling signature in the second phase of upwelling is ~70 km. Therefore, a 7-10 days of the coastal upwelling and its horizontal extension are identified in this study. This study suggests the use of high resolution (~6 km) HF radar current observation on the monitoring of coastal upwelling processes.</p>

Some Observations on the Computational Sensitivity of Rotating Cavity Flows

Heat transfer inside rotating cavities plays an important role in gas turbine engineering. Flows in both compressors and turbine internal flow cavities exhibit self-generated large-scale inertial wave structures, and buoyancy effects are often important. Across the open literature on the topic, there seems to be no clear consensus on what the most suitable modelling fidelity is — although it is a widely held opinion that URANS approaches are less suitable than LES, many authors have succeeded in getting reasonable heat transfer results with URANS. There is also little knowledge of the validity of hybrid URANS/LES type approaches (such as DES) when it comes to predicting the heat transfer in these flows, and furthermore, on the sensitivity of the flow model validity to local driving aerothermal mechanisms in different parts of the cavity. This paper presents the results of a systematic investigation of a rotating cavity with axial throughflow at a Grashof number of 3 × 109. It is found that, for the case investigated, the disk Ekman layers remain laminar. This causes the disk heat transfer to be relatively insensitive to the modelling fidelity used with URANS, DES, and LES giving similar results. The effect of the disk thermal boundary condition is also investigated — it is found to have a significant effect on the direction of the near-wall flow at high radii, despite the large-scale flow structure within the cavity remaining essentially unchanged. This feedback of the disk heat transfer to the near-disk aerodynamics implies that conjugate heat transfer computations of rotating cavities may be worth investigating. On the shroud, URANS fails to resolve the heat transfer enhancement from small-scale buoyancy driven streaks, whilst these are captured by LES. DES also captures these streaks, as the URANS layer within which they are located returns a very small eddy viscosity, and behaves in a similar manner to LES.

Observations of Diurnal Coastal-Trapped Waves with a Thermocline-Intensified Velocity Field

Using 18 days of field observations, we investigate the diurnal (D1) frequency wave dynamics on the Tasmanian eastern continental shelf. At this latitude, the D1 frequency is subinertial and separable from the highly energetic near-inertial motion. We use a linear coastal-trapped wave (CTW) solution with the observed background current, stratification, and shelf bathymetry to determine the modal structure of the first three resonant CTWs. We associate the observed D1 velocity with a superimposed mode-zero and mode-one CTW, with mode one dominating mode zero. Both the observed and mode-one D1 velocity was intensified near the thermocline, with stronger velocities occurring when the thermocline stratification was stronger and/or the thermocline was deeper (up to the shelfbreak depth). The CTW modal structure and amplitude varied with the background stratification and alongshore current, with no spring–neap relationship evident for the observed 18 days. Within the surface and bottom Ekman layers on the shelf, the observed velocity phase changed in the cross-shelf and/or vertical directions, inconsistent with an alongshore propagating CTW. In the near-surface and near-bottom regions, the linear CTW solution also did not match the observed velocity, particularly within the bottom Ekman layer. Boundary layer processes were likely causing this observed inconsistency with linear CTW theory. As linear CTW solutions have an idealized representation of boundary dynamics, they should be cautiously applied on the shelf.

Stratified Ekman layers evolving under a finite-time stabilizing buoyancy flux

Stratified flow in nocturnal boundary layers is studied using direct numerical simulation (DNS) of the Ekman layer, a model problem that is useful to understand atmospheric boundary-layer (ABL) turbulence. A stabilizing buoyancy flux is applied for a finite time to a neutral Ekman layer. Based on previous studies and the simulations conducted here, the choice of $L_{\mathit{cri}}^{+}=Lu_{\ast }/\unicode[STIX]{x1D708}\approx 700$ ($L$ is the Obukhov length scale and $u_{\ast }$ is the friction velocity) provides a cooling flux that is sufficiently strong to cause the initial collapse of turbulence. The turbulent kinetic energy decays over a time scale of $4.06L/u_{\ast }$ during the collapse. The simulations suggest that imposing $L_{\mathit{cri}}^{+}\approx 700$ on the neutral Ekman layer results in turbulence collapse during the initial transient, independent of Reynolds number, $Re_{\ast }$. However, the long-time state of the flow, i.e. turbulent with spatial intermittency or non-turbulent, is found to depend on the initial value of $Re_{\ast }$ since the cooling flux and resultant stratification increase with $Re_{\ast }$ for a given $L^{+}$. The lower-$Re_{\ast }$ cases have sustained turbulence with shear and stratification profiles that evolve in a manner such that the gradient Richardson number, $Ri_{g}$, in the near-surface layer, including the low-level jet, remains subcritical. The highest $Re_{\ast }$ case has supercritical $Ri_{g}$ in the low-level jet and turbulence does not recover. A theoretical discussion is performed to infer that the bulk Richardson number, $Ri_{b}$, is more suitable than $L^{+}$ to determine the fate of stratified Ekman layers at late time. DNS results support the implications of $Ri_{b}$ for the effect of initial $Re_{\ast }$ and $L^{+}$ on the flow.

Energetics of Bottom Ekman Layers during Buoyancy Arrest

Turbulent bottom Ekman layers are among the most important energy conversion sites in the ocean. Their energetics are notoriously complex, in particular near sloping topography, where the feedback between cross-slope Ekman transports, buoyancy forcing, and mixing affects the energy budget in ways that are not well understood. Here, the authors attempt to clarify the energy pathways and different routes to mixing, using a combined theoretical and modeling approach. The analysis is based on a newly developed energy flux diagram for turbulent Ekman layers near sloping topography that allows for an exact definition of the different energy reservoirs and energy pathways. Using a second-moment turbulence model, it is shown that mixing efficiencies increase for increasing slope angle and interior stratification, but do not exceed the threshold of 5% except for very steep slopes, where the canonical value of 20% may be reached. Available potential energy generated by cross-slope advection may equal up to 70% of the energy lost to dissipation for upwelling-favorable flow, and up to 40% for downwelling-favorable flow.