ALBERT

Description: We reconstructed the late Holocene relative sea-level (RSL) evolution of the ancient harbour of Naples, one of the largest coastal conurbations in the Mediterranean. We carried out multiproxy investigations, coupling archaeological evidence with biological indicators. Our data robustly constrain 2000 yr of non-monotonic changes in sea level, chiefly controlled by the complex volcano-tectonic processes that characterize the area. Between ∼200 BC and AD ∼0, a subsidence rate of more than ∼1.5 mm/yr enhanced the postglacial RSL rise, while negligible or moderate land uplift 〈 ∼0.5 mm/yr triggered a RSL stabilization during the Roman period (first five centuries AD). This stabilization was followed by a post-Roman enhancement of the sea-level rise when ground motion was negative, attested by a subsidence rate of ∼0.5 to ∼1 mm/yr. Our analysis seems to indicate very minor impacts of this nonmonotonic RSL evolution on the activities of the ancient harbour of Naples, which peaked from the third century BC to the second century AD. After this period, the progressive silting of the harbour basin made it impossible to safely navigate within the basin, leading to the progressive decline of the harbour.

Description: Published

Description: 284-298

Description: 1V. Storia eruttiva

Description: JCR Journal

Description: Late Quaternary landscapes of unglaciated Beringia were largely shaped by ice-wedge polygon tundra. Ice Complex (IC) strata preserve such ancient polygon formations. Here we report on the Yukagir IC from Bol’shoy Lyakhovsky Island in northeastern Siberia and suggest that new radioisotope disequilibria (230Th/U) dates of the Yukagir IC peat confirm its formation during the MIS 7a-c Interglacial. The preservation of the ice-rich Yukagir IC proves its resilience to Last Interglacial and lateglacial-Holocene warming. This study compares the Yukagir IC to IC strata of MIS 5, MIS 3 and MIS 2 ages exposed on Bol’shoy Lyakhovsky Island. Besides high intrasedimental ice content and syngenetic ice wedges intersecting silts, sandy silts, the Yukagir IC is characterized by high organic matter (OM) accumulation and low OM decomposition of a distinctive Drepanocladus moss-peat. The Yukagir IC pollen data reveals grassshrub-moss tundra indicating rather wet summer conditions similar to modern ones. The stable isotope composition of Yukagir IC wedge ice is similar to those of the MIS 5 and MIS 3 ICs pointing to similar atmospheric moisture generation and transport patterns in winter. Ice Complex data from glacial and interglacial periods provide insights into permafrost and climate dynamics since about 200 ka.

Description: The points where the horizontal component of the geomagnetic field vanishes are located in polar areas, far away from the geomagnetic (analytic) poles and the poles of rotation of the Earth and, differently from the geomagnetic poles, can be found experimentally with a magnetic survey to determine where the field is vertical. The experimental determination of the area where the total field is perfectly vertical, commonly known as dip pole, is not simple, due to the remoteness and harsh climatic conditions; another difficulty is related to the short term geomagnetic field variations, due to the interaction with the external solar wind, which causes the magnetospheric dynamics, particularly evident at high latitude, and as a consequence a displacement of the dip pole. Actually, the study of the dip pole displacements over short time scales can be an important tool for monitoring the magnetospheric dynamics at high latitude. In this study we present the updated location of the the dip poles, using data from the Swarm ESA’s constellation of satellites along their almost polar orbits. We also analyse the spatial shift of these areas during different seasons and interplanetary magnetic field orientations.

Description: Published

Description: 135-138

Description: 1A. Geomagnetismo e Paleomagnetismo

Description: N/A or not JCR

Description: We use low frequency geomagnetic field measurements at two Antarctic stations to statistically investigate the longitudinal location of the polar cusp. The two stations are both located in the polar cap at a geomagnetic latitude close to the cusp latitude; they are separated by one hour in magnetic local time. At each station the Pc5 power maximizes when the station approaches the cusp, i.e. around magnetic local noon. The comparison between the Pc5 power at the two stations allows to determine the longitudinal location of the cusp. Our analysis is conducted considering separately different orientation of the interplanetary magnetic field. The results, which indicate longitudinal shifts of the polar cusp depending on the selected conditions, are discussed in relation to previous studies of the polar cusp location based on polar magnetospheric satellite data.

Description: Published

Description: 139-141

Description: 1A. Geomagnetismo e Paleomagnetismo

Description: N/A or not JCR

Description: 〈div data-abstract-type="normal"〉〈p〉The effectiveness of streak modes in controlling the oblique-type breakdown in a supersonic boundary-layer at Mach 2.0 is investigated using direct numerical simulations. Investigations in the literature have shown the effectiveness of streak modes in delaying the onset of transition dominated by two-dimensional waves, but in oblique breakdown, three-dimensional waves and a strong streak mode dominate the transition process. Paredes 〈span〉et al.〈/span〉 (〈span〉J. Fluid Mech.〈/span〉, vol. 831, 2017, pp. 524–553) discussed the possible stabilization of supersonic boundary layers by optimally growing streaks using parabolized stability equations. However, no study has as yet been reported regarding direct nonlinear control of oblique breakdown. This study deals with the effects of large-amplitude decaying streak modes generated by a blowing–suction strip at the wall to control full breakdown in a reference case. Modes with four to five times the fundamental wavenumber are found to be beneficial for controlling the transition. In the first region after the control-mode forcing, the beneficial mean-flow distortion (MFD), generated by inducing the control mode, is solely responsible for hampering the growth of the fundamental-mode. On the whole, the MFD and the three-dimensional part of the control contribute equally towards controlling the oblique breakdown. The results show significant suppression of transition, and substantial improvements have been observed in the levels of the skin-friction coefficient and wall-temperature in comparison to the uncontrolled case. Moreover, refreshing the control using an additional downstream control strip increases the gain. However, the forcing amplitude must be carefully chosen in order not to introduce a generalized inflection point in the spanwise averaged mean flow invoking enhanced disturbance growth.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The problem of interaction between disturbances and shock waves was solved by a theoretical approach called linear interaction analysis in the mid-twentieth century. More recently, great progress has been made in analysing shock–turbulence interactions by direct numerical simulation. However, an unsolved theoretical problem remains: What happens when no acoustic waves are stimulated behind the shock wave? The concept of a damped wave is introduced, which is a type of excited plane wave. Based on this, the dispersion and amplitude relationships between any incident plane wave and resulting stimulated waves are constructed analytically, systematically and comprehensively. The physical essence of damped waves and the existence of critical angles are clarified. It is demonstrated that a damped wave is a complex number space solution to the acoustic dispersion relationship under certain conditions. It acts as a bridge connecting fast and slow acoustic waves at the position where the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628130807083-0762:S0022112019004385:S0022112019004385_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 component of the group velocity is zero. There are two critical angles that can excite fast and slow acoustic waves, which determine the conditions that stimulate a damped wave. Our results show good agreement with theoretical and simulation results. The contribution of each excited wave to the transmission coefficient is evaluated, the distribution of the transmission coefficient is analysed and application to an engineering wedge model is performed.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The onset of thermal convection in a rapidly rotating spherical shell is studied by linear stability analysis based on the fully compressible Navier–Stokes equations. Compressibility is quantified by the number of density scale heights 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, which measures the intensity of density stratification of the motionless, polytropic base state. The nearly adiabatic flow with polytropic index 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is considered, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the adiabatic polytropic index. By investigating the stability of the base state with respect to the disturbance of specified wavenumber, the instability process is found to be sensitive to the Prandtl number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. For large 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the quasi-geostrophic columnar mode loses stability first; while for relatively small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 a new quasi-geostrophic compressible mode is identified, which becomes unstable first under strong density stratification. The inertial mode can also occur first for relatively small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and a certain intensity of density stratification in the parameter range considered. Although the Rayleigh numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the onsets of the quasi-geostrophic compressible mode and columnar mode are different by several orders of magnitude, we find that they follow very similar scaling laws with the Taylor number. The critical 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for convection onset is found to be always positive, in contrast with previous results based on the widely used anelastic model that convection can occur at negative 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. By evaluating the relative magnitude of the time derivative of density perturbation in the continuity equation, we show that the anelastic approximation in the present system cannot be applied in the small-〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and large-〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628092857401-0357:S0022112019004361:S0022112019004361_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 regime.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We study the breakup of confined fluid threads at low flow rates to understand instability mechanisms. To determine the critical conditions between the earlier quasi-stable necking stage and the later unstable collapse stage, simulations and experiments are designed to operate at an extremely low flow rate. The critical mean radii at the neck centres are identified by the stop-flow method for elementary microfluidic configurations. Two distinct origins of capillary instabilities are revealed for different confinement situations. One is the gradient of capillary pressure induced by the confinements of geometry and external flow, whereas the other is the competition between the capillary pressure and internal pressure determined by the confinements.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We report on a combined experimental and numerical study of convective heat transfer along ratchet surfaces in vertical natural convection (VC). Due to the asymmetry of the convection system caused by the asymmetric ratchet-like wall roughness, two distinct states exist, with markedly different orientations of the large-scale circulation roll (LSCR) and different heat transport efficiencies. Statistical analysis shows that the heat transport efficiency depends on the strength of the LSCR. When a large-scale wind flows along the ratchets in the direction of their smaller slopes, the convection roll is stronger and the heat transport is larger than the case in which the large-scale wind is directed towards the steeper slope side of the ratchets. Further analysis of the time-averaged temperature profiles indicates that the stronger LSCR in the former case triggers the formation of a secondary vortex inside the roughness cavity, which promotes fluid mixing and results in a higher heat transport efficiency. Remarkably, this result differs from classical Rayleigh–Bénard convection (RBC) with asymmetric ratchets (Jiang 〈span〉et al.〈/span〉, 〈span〉Phys. Rev. Lett.〈/span〉, vol. 120, 2018, 044501), wherein the heat transfer is stronger when the large-scale wind faces the steeper side of the ratchets. We reveal that the reason for the reversed trend for VC as compared to RBC is that the flow is less turbulent in VC at the same 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628072522266-0966:S0022112019004464:S0022112019004464_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Thus, in VC the heat transport is driven primarily by the coherent LSCR, while in RBC the ejected thermal plumes aided by gravity are the essential carrier of heat. The present work provides opportunities for control of heat transport in engineering and geophysical flows.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We explore the dynamics of inclined temporal gravity currents using direct numerical simulation, and find that the current creates an environment in which the flux Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, gradient Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and turbulent flux coefficient 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are constant across a large portion of the depth. Changing the slope angle 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 modifies these mixing parameters, and the flow approaches a maximum Richardson number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 at which the entrainment coefficient 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The turbulent Prandtl number remains 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for all slope angles, demonstrating that 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is not caused by a switch-off of the turbulent buoyancy flux as conjectured by Ellison (〈span〉J. Fluid Mech.〈/span〉, vol. 2, 1957, pp. 456–466). Instead, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 occurs as the result of the turbulence intensity going to zero as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, due to the flow requiring larger and larger shear to maintain the same level of turbulence. We develop an approximate model valid for small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 which is able to predict accurately 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline16.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and their maximum attainable values. The model predicts an entrainment law of the form 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190627132031350-0413:S0022112019004300:S0022112019004300_inline17.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, which is in good agreement with the simulation data. The simulations and model presented here contribute to a growing body of evidence that an approach to a marginally or critically stable, relatively weakly stratified equilibrium for stratified shear flows may well be a generic property of turbulent stratified flows.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Ice scallops are a small-scale (5–20 cm) quasi-periodic ripple pattern that occurs at the ice–water interface. Previous work has suggested that scallops form due to a self-reinforcing interaction between an evolving ice-surface geometry, an adjacent turbulent flow field and the resulting differential melt rates that occur along the interface. In this study, we perform a series of laboratory experiments in a refrigerated flume to quantitatively investigate the mechanisms of scallop formation and evolution in high resolution. Using particle image velocimetry, we probe an evolving ice–water boundary layer at sub-millimetre scales and 15 Hz frequency. Our data reveal three distinct regimes of ice–water interface evolution: a transition from flat to scalloped ice; an equilibrium scallop geometry; and an adjusting scallop interface. We find that scalloped-ice geometry produces a clear modification to the ice–water boundary layer, characterized by a time-mean recirculating eddy feature that forms in the scallop trough. Our primary finding is that scallops form due to a self-reinforcing feedback between the ice-interface geometry and shear production of turbulent kinetic energy in the flow interior. The length of this shear production zone is therefore hypothesized to set the scallop wavelength.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Direct numerical simulation (DNS) is performed for two wall-bounded flow configurations: laminar Couette flow at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and turbulent channel flow at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the shear stress at the wall. The top wall is smooth and the bottom wall is a realistically rough superhydrophobic surface (SHS), generated from a three-dimensional surface profile measurement. The air–water interface, which is assumed to be flat, is simulated using the volume-of-fluid (VOF) approach. The two flow cases are studied with varying interface heights 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to understand its effect on slip and drag reduction (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). For the laminar Couette flow case, the presence of the surface roughness is felt up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 of the channel height in the wall-normal direction. Nonlinear dependence of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 on 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is observed with three distinct regions. A nonlinear curve fit is obtained for gas fraction 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 determines the amount of slip area exposed to the flow. A power law fit is obtained from the data for the effective slip length as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and is compared to those derived for structured geometry. For the turbulent channel flow, statistics of the flow field are compared to that of a smooth wall to understand the effects of roughness and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Four cases are simulated ranging from fully wetted to fully covered and two intermediate regions in between. Scaling laws for slip length, slip velocity, roughness function and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are obtained for different penetration depths and are compared to past work for structured geometry. 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190628070320625-0031:S0022112019004191:S0022112019004191_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is shown to depend on a competing effect between slip velocity and turbulent losses due to the Reynolds shear stress contribution. Presence of trapped air in the cavities significantly alters near-wall flow physics where we examine near-wall structures and propose a physical mechanism for their behaviour. The fully wetted roughness increases the peak value of turbulent intensities, whereas the presence of the interface suppresses them. The pressure fluctuations have competing contributions between turbulent pressure fluctuations and stagnation due to asperities, the near-wall structure is altered and breaks down with increasing slip. Overall, there exists a competing effect between the interface and the asperities, the interface suppresses turbulence whereas the asperities enhance them. The present work demonstrates DNS over a realistic multiphase SHS for the first time, to the best of our knowledge.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We developed a numerical method for the set of equations governing fully compressible convection in the limit of infinite Prandtl numbers. Reduced models have also been analysed, such as the anelastic approximation and the anelastic liquid approximation. The tests of our numerical schemes against self-consistent criteria have shown that our numerical simulations are consistent from the point of view of energy dissipation, heat transfer and entropy budget. The equation of state of an ideal gas has been considered in this work. Specific effects arising because of the compressibility of the fluid are studied, like the scaling of viscous dissipation and the scaling of the heat flux contribution due to the mechanical power exerted by viscous forces. We analysed the solutions obtained with each model (fully compressible model, anelastic and anelastic liquid approximations) in a wide range of dimensionless parameters and determined the errors induced by each approximation with respect to the fully compressible solutions. Based on a rationale on the development of the thermal boundary layers, we can explain reasonably well the differences between the fully compressible and anelastic models, in terms of both the heat transfer and viscous dissipation dependence on compressibility. This could be mostly an effect of density variations on thermal diffusivity. Based on the different forms of entropy balance between exact and anelastic models, we find that a necessary condition for convergence of the anelastic results to the exact solutions is that the product 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 must be small compared to unity, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the ratio of the superadiabatic temperature difference to the adiabatic difference, and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the ratio of the superadiabatic heat flux to the heat flux conducted along the adiabat. The same condition seems also to be associated with a convergence of the computed heat fluxes. Concerning the anelastic liquid approximation, we confirm previous estimates by Anufriev 〈span〉et al.〈/span〉 (〈span〉Phys. Earth Planet. Inter.〈/span〉, vol. 152, 2005, pp. 163–190) and find that its results become generally close to those of the fully compressible model when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is small compared to unity, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the isobaric thermal expansion coefficient, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the temperature (here 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for an ideal gas) and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190625082702251-0920:S0022112019004208:S0022112019004208_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the dissipation number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor–Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, three-dimensional (3-D) roughness has not been previously considered in depth. Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (〈span〉Phys. Fluids〈/span〉, vol. 18, 031701, 2006) has been modified to simulate a random rough surface of monodisperse sand grains. Taylor numbers range from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉(corresponding to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). We focus on the influence of the roughness height 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. We analyse the global response of the system, expressed both by the dimensionless angular velocity transport 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and by the friction factor 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. An increase in friction with increasing roughness height is accompanied with enhanced plume ejection from the inner cylinder. Subsequently, we investigate the local response of the fluid flow over the rough surface. The equivalent sand grain roughness 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is calculated to be 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the size of the sand grains. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behaviour to that of the Nikuradse (〈span〉VDI-Forsch.〈/span〉, vol. 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer. Furthermore, we find that the dynamical effects of the sand grains are contained to the roughness sublayer 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 with 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621104521411-0582:S0022112019003768:S0022112019003768_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Geometric, kinematic and dynamic properties of focusing deep-water surface gravity wave packets are examined in a simplified model with the intent of deriving a wave breaking threshold parameter. The model is based on the spatial modified nonlinear Schrödinger equation of Dysthe (〈span〉Proc. R. Soc. Lond.〈/span〉 A, vol. 369 (1736), 1979, pp. 105–114). The evolution of initially narrow-banded and weakly nonlinear chirped Gaussian wave packets are examined, by means of a trial function and a variational procedure, yielding analytic solutions describing the approximate evolution of the packet width, amplitude, asymmetry and phase during focusing. A model for the maximum free surface gradient, as a function of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 the linear prediction of the maximum slope at focusing and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 the non-dimensional packet bandwidth, is proposed and numerically examined, indicating a quasi-self-similarity of these focusing events. The equations of motion for the fully nonlinear potential flow equations are then integrated to further investigate these predictions. It is found that a model of this form can characterize the bulk partitioning of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621105221417-0357:S0022112019004282:S0022112019004282_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 phase space, between non-breaking and breaking waves, serving as a breaking criterion. Application of this result to better understanding air–sea interaction processes is discussed.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We conduct direct numerical simulations (DNS) of the Cahn–Hilliard–Navier–Stokes (CHNS) equations to investigate the statistical properties of a turbulent phase-separating symmetric binary-fluid mixture. Turbulence causes an arrest of the phase separation which leads to the formation of a statistically steady emulsion. We characterise turbulent velocity fluctuations in an emulsion for different values of the Reynolds number and the Weber number. Our scale-by-scale kinetic energy budget analysis shows that the interfacial terms in the CHNS equations provide an alternative route for the kinetic energy transfer. By studying the probability distribution function (p.d.f.) of the energy dissipation rate, the vorticity magnitude and the joint-p.d.f. of the velocity-gradient invariants we show that the statistics of the turbulent fluctuations do not change with the Weber number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉With the aim to characterize the near-wall flow structures and their interaction with large-scale motions in the log-law region, time-resolved planar and volumetric flow field measurements were performed in the near-wall and log-law region of an adverse pressure gradient turbulent boundary layer following a zero pressure gradient turbulent boundary layer at a friction Reynolds number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Due to the high spatial and temporal resolution of the measurements, it was possible to resolve and identify uniform-momentum zones in the region 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 or 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621113954055-0112:S0022112019004087:S0022112019004087_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and to relate them with well known coherent flow motions near the wall. The space–time results confirm that the turbulent superstructures have a strong impact even on the very near-wall flow motion and also their alternating appearance in time and intensity could be quantified over long time sequences. Using the time record of the velocity field, rare localized separation events appearing in the viscous sublayer were also analysed. By means of volumetric particle tracking velocimetry their three-dimensional topology and dynamics could be resolved. Based on the results, a conceptual model was deduced that explains their rare occurrence, topology and dynamics by means of a complex interaction process between low-momentum turbulent superstructures, near-wall low-speed streaks and tilted longitudinal and spanwise vortices located in the near-wall region.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We experimentally investigate the extensional flow of a sheet – or curtain – of viscoelastic liquid falling freely from a slot at constant flow rate under gravity. Extruded liquids are aqueous solutions of flexible polyethylene oxide (PEO) and of semi-rigid partially hydrolysed polyacrylamide (HPAM) with low shear viscosities. Velocimetry measurements reveal that the mean velocity field 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the distance from the slot exit) does not reduce to a free fall. More precisely, we show that the liquid falls initially with sub-gravitational accelerations up to a distance from the slot which scales as 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is gravity and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the extensional relaxation time of the liquid) due to the stretching of polymer molecules. Beyond this elastic length, inertia dominates and the local acceleration reaches the asymptotic free-fall value 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The length of the sub-gravitational part of the curtain is shown to be much larger than the equivalent viscous length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for Newtonian liquids of density 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and dynamic viscosity 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 which is usually small compared to the curtain length. By analogy with Newtonian curtains, we show that the velocity field 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621102720940-0917:S0022112019003896:S0022112019003896_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 rescales on a master curve. Besides, the flow is shown to be only weakly affected by the history of polymer deformations in the die upstream of the curtain. Furthermore, investigations on the curtain stability reveal that polymer addition reduces the minimum flow rate required to maintain a continuous sheet of liquid.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Previous experiments have revealed that shock waves driven through dissipative media may become unstable, for example, in granular gases, and in molecular gases undergoing strong relaxation effects. The current paper addresses this problem of shock stability at the Euler and Navier–Stokes continuum levels in a system of disks (two-dimensional) undergoing activated inelastic collisions. The dynamics of shock formation and stability is found to be in very good agreement with earlier molecular dynamic simulations (Sirmas & Radulescu, 〈span〉Phys. Rev.〈/span〉 E, vol. 91, 2015, 023003). It was found that the modelling of shock instability requires the introduction of molecular noise for its development and sustenance. This is confirmed in two stability problems. In the first, the evolution of shock formation dynamics is monitored without noise, with only initial noise and with continuous molecular noise. Only the latter reproduces the results of shock instability of molecular dynamics simulations. In the second problem, the steady travelling wave solution is obtained for the shock structure in the inviscid and viscous limits and its nonlinear stability is studied with and without molecular fluctuations, again showing that instability can be sustained only in the presence of fluctuations. The continuum results show that instability takes the form of a rippled front of a wavelength comparable with the relaxation thickness of the steady shock wave, at scales at which molecular fluctuations become important, in excellent agreement with the molecular dynamic simulations.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We performed two-way coupled direct numerical simulations of turbulent channel flow with Lagrangian tracking of small, heavy spheres at a dimensionless gravitational acceleration of 0.077 in wall units, which is based on the flow condition in the experiment by Gerashchenko 〈span〉et al.〈/span〉 (〈span〉J. Fluid Mech.〈/span〉, vol. 617, 2008, pp. 255–281). We removed deposited particles after several collisions with the lower wall and then released new particles near the upper wall to observe direct interactions between particles and coherent structures of near-wall turbulence during gravitational settling through the mean shear. The results indicate that when the Stokes number is approximately 1 on the basis of the Kolmogorov time scale of the flow (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190624134045701-0463:S0022112019004002:S0022112019004002_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉), the so-called preferential sweeping occurs in association with coherent streamwise vortices, while the effect of crossing trajectories becomes significant for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190624134045701-0463:S0022112019004002:S0022112019004002_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Consequently, in either case, the settling particles deposit on the wall without strong accumulation in low-speed streaks in the viscous sublayer. When particles settle through near-wall turbulence from the upper wall, more small-scale vortical structures are generated in the outer layer as low-speed fluid is pulled farther in the direction of gravity, while the opposite is true near the lower wall.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Streaks have been found to be an important part of wall-turbulence dynamics. In this paper, we extend the analysis for unbounded shear flows, in particular a Mach 0.4 round jet, using measurements taken using dual-plane, time-resolved, stereoscopic particle image velocimetry (PIV) taken at pairs of jet cross-sections, allowing the evaluation of the cross-spectral density of streamwise velocity fluctuations resolved into azimuthal Fourier modes. From the streamwise velocity results, two analyses are performed: the evaluation of wavenumber spectra (assuming Taylor’s hypothesis for the streamwise coordinate) and a spectral proper orthogonal decomposition (SPOD) of the velocity field using PIV planes in several axial stations. The methods complement each other, leading to the conclusion that large-scale streaky structures are also present in turbulent jets where they experience large growth in the streamwise direction, energetic structures extending up to eight diameters from the nozzle exit. Leading SPOD modes highlight the large-scale, streaky shape of the structures, whose aspect ratio (streamwise over azimuthal length) is approximately 15. The data were further analysed using SPOD, resolvent and transient growth analyses, good agreement being observed between the models and the leading SPOD mode for the wavenumbers considered. The models also indicate that the lift-up mechanism is active in turbulent jets, with streamwise vortices leading to streaks. The results show that large-scale streaks are a relevant part of the jet dynamics.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We study single-phase and particle-laden turbulent channel flows bounded by two incompressible hyper-elastic walls with different deformability at bulk Reynolds number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The solid volume fraction of finite-size neutrally buoyant rigid spherical particles considered is 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The elastic walls are assumed to be of a neo-Hookean material. A fully Eulerian formulation is employed to model the elastic walls together with a direct-forcing immersed boundary method for the coupling between the fluid and the particles. The data show a significant drag increase and the enhancement of the turbulence activity with growing wall elasticity for both the single-phase and particle-laden flows when compared with the single-phase flow over rigid walls. Drag reduction and turbulence attenuation is obtained, on the other hand, with highly elastic walls when comparing the particle-laden flow with the single-phase flow for the same wall properties; the opposite effect, drag increase, is observed upon adding particles to the flow over less elastic walls. This is explained by investigating the near-wall turbulence, where the strong asymmetry in the magnitude of the wall-normal velocity fluctuations (favouring positive 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉), is found to push the particles towards the channel centre. The particle layer close to the wall contributes to turbulence production by increasing the wall-normal velocity fluctuations, so that in the absence of this layer, smaller wall deformations and in turn turbulence attenuation is observed. For a moderate wall elasticity, we increase the particle volume fraction up to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621111752900-0025:S0022112019004130:S0022112019004130_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and find that particle migration away from the wall is the cause of turbulence attenuation with respect to the flow over rigid walls. However, for this higher volume fractions, the particle induced stress compensates for the decreasing Reynolds shear stress, resulting in a higher overall drag for the case with elastic walls. The effect of the wall elasticity on the overall drag reduces significantly with increasing particle volume fraction.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present large-eddy simulation (LES) of flow past different airfoils with 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, based on the free-stream velocity and airfoil chord length, ranging from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. To avoid the challenging resolution requirements of the near-wall region, we develop a virtual wall model in generalized curvilinear coordinates and incorporate the non-equilibrium effects via proper treatment of the momentum equations. It is demonstrated that the wall model dynamically captures the instantaneous skin-friction vector field on arbitrary curved surfaces at the resolved scale. By combining the present wall model with the stretched-vortex subgrid-scale model, we apply the wall-modelled LES approach to three different airfoil cases, spanning different geometrical parameters, different attack angles and low to high 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The numerical results are verified with direct numerical simulation (DNS) at low 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, and validated with experiment data at higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, including typical aerodynamic properties such as pressure coefficient distributions, velocity components and also more challenging measurements such as skin-friction coefficient and Reynolds stresses. All comparisons show reasonable agreement, providing a measure of validity that enables us to further probe simulation results into aspects of flow physics that are not available from experiments. Two techniques to quantify hitherto unexplored physics of flows past airfoils are employed: one is the construction of the anisotropy invariant map, and the second is skin-friction portraits with emphasis on flow transition and unsteady separation along the airfoil surface. The anisotropy maps for all three 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 cases, show clearly that a portion of the flow field is aligned along the axisymmetric expansion line, corresponding to the turbulent boundary layer log-law behaviour and the appearance of turbulent transition. The instantaneous skin-friction portraits reveal a monotonic shrinking of the near wall structure scale. At 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the interaction between the primary separation bubble and the secondary separation bubble contributes to turbulent transition, similar to the case of flow past a cylinder. At higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the primary separation breaks into several small separation bubbles. At even higher 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190621103821391-0158:S0022112019003604:S0022112019003604_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, near the turbulent separation, the skin-friction lines show small-scale reversal flows that are similar to those observed in DNS of the flat plate turbulent separation. A notable feature of turbulent separation in flow past an airfoil is the appearance of turbulence structures and small-scale reversal flows in the spanwise direction due to the vortex shedding behaviour.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In the framework of the complete formulation of the conjugate problem, the liquid–gas flow structure arising upon local heating using thermal sources is investigated numerically. The two-layer system is confined by solid impermeable walls. The Navier–Stokes equations in the Boussinesq approximation in the ‘streamfunction–vorticity’ variables are used to describe the media motion. The dynamic conditions at the interface are formulated in terms of the tangential and normal velocities, while the temperature conditions at the external boundaries of the system take into account the presence of local heaters. The influence of the number of heaters and heating modes on the dynamics and character of the appearing convective regimes is analysed. The steady and commutated heating modes for one and two heaters arranged at the lower boundary are investigated. The heating initiates convective and thermocapillary mechanisms causing the fluid motion. Transient regimes with the successive formation of two-vortex, quadruple-vortex and two-vortex flows are observed before the stabilization of the system in the uniform heating mode. A stable thermocapillary deflection appears at the interface above the heater. The commutated mode of heating entails oscillations of the interface with a change in the deflection form and the formation of travelling vortices in the fluids. The impact of particular mechanisms on the flow patterns is analysed. The paper presents typical distributions of the velocity and temperature fields in the system and the position of the interface for the considered cases.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In the context of dynamic wetting, wall slip is often treated as a microscopic effect for removing viscous stress singularity at a moving contact line. In most drop spreading experiments, however, a considerable amount of slip may occur due to the use of polymer liquids such as silicone oils, which may cause significant deviations from the classical Tanner–de Gennes theory. Here we show that many classical results for complete wetting fluids may no longer hold due to wall slip, depending crucially on the extent of de Gennes’s slipping ‘foot’ to the relevant length scales at both the macroscopic and microscopic levels. At the macroscopic level, we find that for given liquid height 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and slip length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the apparent dynamic contact angle 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 can change from Tanner’s law 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to the strong-slip law 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the capillary number and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the macroscopic length scale. Such a no-slip-to-slip transition occurs at the critical capillary number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, accompanied by the switch of the ‘foot’ of size 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 from the inner scale to the outer scale with respect to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. A more generalized dynamic contact angle relationship is also derived, capable of unifying Tanner’s law and the strong-slip law under 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. We not only confirm the two distinct wetting laws using many-body dissipative particle dynamics simulations, but also provide a rational account for anomalous departures from Tanner’s law seen in experiments (Chen, 〈span〉J. Colloid Interface Sci〈/span〉., vol. 122, 1988, pp. 60–72; Albrecht 〈span〉et al.〈/span〉, 〈span〉Phys. Rev. Lett.〈/span〉, vol. 68, 1992, pp. 3192–3195). We also show that even for a common spreading drop with small macroscopic slip, slip effects can still be microscopically strong enough to change the microstructure of the contact line. The structure is identified to consist of a strongly slipping precursor film of length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 followed by a mesoscopic ‘foot’ of width 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 ahead of the macroscopic wedge, where 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline16.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is the molecular length. It thus turns out that it is the ‘foot’, rather than the film, contributing to the microscopic length in Tanner’s law, in accordance with the experimental data reported by Kavehpour 〈span〉et al.〈/span〉 (〈span〉Phys. Rev. Lett.〈/span〉, vol. 91, 2003, 196104) and Ueno 〈span〉et al.〈/span〉 (〈span〉Trans. ASME J. Heat Transfer〈/span〉, vol. 134, 2012, 051008). The advancement of the microscopic contact line is still led by the film whose length can grow as the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline17.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 power of time due to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619140754267-0248:S0022112019003525:S0022112019003525_inline18.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, as supported by the experiments of Ueno 〈span〉et al.〈/span〉 and Mate (〈span〉Langmuir〈/span〉, vol. 28, 2012, pp. 16821–16827). The present work demonstrates that the behaviour of a moving contact line can be strongly influenced by wall slip. Such slip-mediated dynamic wetting might also provide an alternative means for probing slippery surfaces.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We numerically investigated the unsteady dynamics of a two-dimensional airfoil undergoing a continuous, prescribed pitch-up motion and freely translating as a response to aerodynamic forces and the gravity field. The pitch-up motion was applied about an axis located 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190619142124772-0727:S002211201900421X:S002211201900421X_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 chord away from the leading edge and was parameterized using the shape change number, with a Reynolds number set to 2000. It was shown that the minimum kinetic energy reached by the airfoil depends stochastically and asymptotically on shape change numbers for values below and above 1, respectively. Very low kinetic energy levels (close to zero) can be reached in both stochastic and asymptotic regions but high shape change numbers are accompanied by significant gain in altitude which may be undesirable from a practical perspective. Rather, shape change numbers in the range [0.1–0.3] allow us to reach relatively low levels of kinetic energy for close perching locations. We showed that highly nonlinear fluid–structure interactions induced by massive flow separations and strong vortices are conducive to low kinetic energy, but responsible for the stochastic dependence of kinetic energy to shape change number, which can make perching manoeuvres hardly controllable for flying vehicles.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present a dynamic decomposition analysis of the wake flow in fluid–structure interaction (FSI) systems under both laminar and turbulent flow conditions. Of particular interest is to provide the significance of low-dimensional wake flow features and their interaction dynamics to sustain the free vibration of a square cylinder at a relatively low mass ratio. To obtain the high-dimensional data, we employ a body-conforming variational FSI solver based on the recently developed partitioned iterative scheme and the dynamic subgrid-scale turbulence model for a moderate Reynolds number (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉). The snapshot data from high-dimensional FSI simulations are projected to a low-dimensional subspace using the proper orthogonal decomposition (POD). We utilize each corresponding POD mode to detect features of the organized motions, namely, the vortex street, the shear layer and the near-wake bubble. We find that the vortex shedding modes contribute solely to the lift force, while the near-wake and shear layer modes play a dominant role in the drag force. We further examine the fundamental mechanism of this dynamical behaviour and propose a force decomposition technique via low-dimensional approximation. To elucidate the frequency lock-in, we systematically analyse the decomposed modes and their dynamical contributions to the force fluctuations for a range of reduced velocity at low Reynolds number laminar flow. These quantitative mode energy contributions demonstrate that the shear layer feeds the vorticity flux to the wake vortices and the near-wake bubble during the wake–body synchronization. Based on the decomposition of wake dynamics, we suggest an interaction cycle for the frequency lock-in during the wake–body interaction, which provides the interrelationship between the high-amplitude motion and the dominating wake features. Through our investigation of wake–body synchronization below critical 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 range, we discover that the bluff body can undergo a synchronized high-amplitude vibration due to flexibility-induced unsteadiness. Owing to the wake turbulence at a moderate Reynolds number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327055129472-0988:S002211201900140X:S002211201900140X_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, a distorted set of POD modes and the broadband energy distribution are observed, while the interaction cycle for the wake synchronization is found to be valid for the turbulent wake flow.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We use resolvent analysis to design active control techniques for separated flows over a NACA 0012 airfoil. Spanwise-periodic flows over the airfoil at a chord-based Reynolds number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and a free-stream Mach number of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are considered at two post-stall angles of attack of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325174522297-0677:S0022112019001630:S0022112019001630_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Near the leading edge, localized unsteady thermal actuation is introduced in an open-loop manner with two tunable parameters of actuation frequency and spanwise wavelength. To provide physics-based guidance for the effective choice of these control input parameters, we conduct global resolvent analysis on the baseline turbulent mean flows to identify the actuation frequency and wavenumber that provide large perturbation energy amplification. The present analysis also considers the use of a temporal filter to limit the time horizon for assessing the energy amplification to extend resolvent analysis to unstable base flows. We incorporate the amplification and response mode from resolvent analysis to provide a metric that quantifies momentum mixing associated with the modal structure. This metric is compared to the results from a large number of three-dimensional large-eddy simulations of open-loop controlled flows. With the agreement between the resolvent-based metric and the enhancement of aerodynamic performance found through large-eddy simulations, we demonstrate that resolvent analysis can predict the effective range of actuation frequency as well as the global response to the actuation input. We believe that the present resolvent-based approach provides a promising path towards mean flow modification by capitalizing on the dominant modal mixing.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In this work we present and demonstrate the reliability of a theoretical framework for the study of thermally driven turbulence. It consists of scale-by-scale budget equations for the second-order velocity and temperature structure functions and their limiting cases, represented by the turbulent kinetic energy and temperature variance budgets. This framework represents an extension of the classical Kolmogorov and Yaglom equations to inhomogeneous and anisotropic flows, and allows for a novel assessment of the turbulent processes occurring at different scales and locations in the fluid domain. Two relevant characteristic scales, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the velocity field and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 for the temperature field, are identified. These variables separate the space of scales into a quasi-homogeneous range, characterized by turbulent kinetic energy and temperature variance cascades towards dissipation, and an inhomogeneity-dominated range, where the production and the transport in physical space are important. This theoretical framework is then extended to the context of large-eddy simulation to quantify the effect of a low-pass filtering operation on both resolved and subgrid dynamics of turbulent Rayleigh–Bénard convection. It consists of single-point and scale-by-scale budget equations for the filtered velocity and temperature fields. To evaluate the effect of the filter length 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 on the resolved and subgrid dynamics, the velocity and temperature fields obtained from a direct numerical simulation are split into filtered and residual components using a spectral cutoff filter. It is found that when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is smaller than the minimum values of the cross-over scales given by 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the resolved processes correspond to the exact ones, except for a depletion of viscous and thermal dissipations, and the only role of the subgrid scales is to drain turbulent kinetic energy and temperature variance to dissipate them. On the other hand, the resolved dynamics is much poorer in the near-wall region and the effects of the subgrid scales are more complex for filter lengths of the order of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327124217270-0217:S0022112019001198:S0022112019001198_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 or larger. This study suggests that classic eddy-viscosity/diffusivity models employed in large-eddy simulation may suffer from some limitations for large filter lengths, and that alternative closures should be considered to account for the inhomogeneous processes at subgrid level. Moreover, the theoretical framework based on the filtered Kolmogorov and Yaglom equations may represent a valuable tool for future assessments of the subgrid-scale models.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Motivated by recent studies that have revealed the existence of trapped acoustic waves in subsonic jets (Towne 〈span〉et al.〈/span〉, 〈span〉J. Fluid Mech.〈/span〉, vol. 825, 2017, pp. 1113–1152), we undertake a more general exploration of the physics associated with acoustic modes in jets and wakes, using a double vortex-sheet model. These acoustic modes are associated with eigenvalues of the vortex-sheet dispersion relation; they are discrete modes, guided by the vortex sheet; they may be either propagative or evanescent; and under certain conditions they behave in the manner of acoustic-duct modes. By analysing these modes we show how jets and wakes may both behave as waveguides under certain conditions, emulating ducts with soft or hard walls, with the vortex-sheet impedance providing effective ‘wall’ conditions. We consider, in particular, the role that upstream-travelling acoustic modes play in the dispersion-relation saddle points that underpin the onset of absolute instability. The analysis illustrates how departure from duct-like behaviour is a necessary condition for absolute instability, and this provides a new perspective on the stabilising and destabilising effects of reverse flow, temperature ratio and compressibility; it also clarifies the differing symmetries of jet (symmetric) and wake (antisymmetric) instabilities. An energy balance, based on the vortex-sheet impedance, is used to determine stability conditions for the acoustic modes: these may become unstable in supersonic flow due to an energy influx through the shear layers. Finally, we construct the impulse response of flows with zero and finite shear-layer thickness. This allows us to show how the long-time wavepacket behaviour is indeed determined by interaction between Kelvin–Helmholtz and acoustic modes.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Following the recent observation that turbulent pipe flow can be relaminarised by a relatively simple modification of the mean velocity profile, we here carry out a quantitative experimental investigation of this phenomenon. Our study confirms that a flat velocity profile leads to a collapse of turbulence and in order to achieve the blunted profile shape, we employ a moving pipe segment that is briefly and rapidly shifted in the streamwise direction. The relaminarisation threshold and the minimum shift length and speeds are determined as a function of Reynolds number. Although turbulence is still active after the acceleration phase, the modulated profile possesses a severely decreased lift-up potential as measured by transient growth. As shown, this results in an exponential decay of fluctuations and the flow relaminarises. While this method can be easily applied at low to moderate flow speeds, the minimum streamwise length over which the acceleration needs to act increases linearly with the Reynolds number.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We investigate the gravity-driven flow of a thin film of liquid metal on a conducting conical substrate in the presence of a strong toroidal magnetic field (transverse to the flow and parallel to the substrate). We solve the leading-order governing equations in a physically relevant asymptotic limit to find the free-surface profile. We find that the leading-order fluid flow rate is a non-monotonic bounded function of the film height, and this can lead to singularities in the free-surface profile. We perform a detailed stability analysis and identify values of the relevant geometric, hydrodynamic and magnetic parameters such that the flow is stable.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We present a local stability analysis to investigate the effects of differential diffusion between momentum and density (quantified by the Schmidt number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) on the three-dimensional, short-wavelength instabilities in planar vortices with a uniform stable stratification along the vorticity axis. Assuming small diffusion in both momentum and density, but arbitrary values for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, we present a detailed analytical/numerical analysis for three different classes of base flows: (i) an axisymmetric vortex, (ii) an elliptical vortex and (iii) the flow in the neighbourhood of a hyperbolic stagnation point. While a centrifugally stable axisymmetric vortex remains stable for any 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, it is shown that 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 can have significant effects in a centrifugally unstable axisymmetric vortex: the range of unstable perturbations increases when 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is taken away from unity, with the extent of increase being larger for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 than for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. Additionally, for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, we report the possibility of oscillatory instability. In an elliptical vortex with a stable stratification, 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is shown to non-trivially influence the three different inviscid instabilities (subharmonic, fundamental and superharmonic) that have been previously reported, and also introduce a new branch of oscillatory instability that is not present at 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The unstable parameter space for the subharmonic (instability IA) and fundamental (instability IB) inviscid instabilities are shown to be significantly increased for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline12.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, respectively. Importantly, for sufficiently small and large 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline13.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, respectively, the maximum growth rate for instabilities IA and IB occurs away from the inviscid limit. The new oscillatory instability (instability III) is shown to occur only for sufficiently small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline14.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the signature of which is nevertheless present with zero growth rate in the inviscid limit. The Schmidt number is then shown to play no role in the evolution of transverse perturbations on the flow around a hyperbolic stagnation point with a stable stratification. We conclude by discussing the physical length scales associated with the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190327085207695-0702:S0022112019001472:S0022112019001472_inline15.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 instabilities, and their potential relevance in various realistic settings.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Instabilities and flow characteristics in the far wake of a circular cylinder are examined through direct numerical simulations. The transitions to the two-layered and secondary vortex streets are quantified by a new method based on the time-averaged transverse velocity field. Two processes for the transition to the secondary vortex street are observed: (i) the merging of two same-sign vortices over a range of low Reynolds numbers (〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) between 200 and 300, and (ii) the pairing of two opposite-sign vortices, followed by the merging of the paired vortices into subsequent vortices, over a range of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 between 400 and 1000. Single vortices may be generated between the merging cycles due to mismatch of the vortices. The irregular merging process results in flow irregularity and an additional frequency signal 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 (in addition to the primary vortex shedding frequency 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) in the two-layered and secondary vortex streets. In particular, a gradual energy transfer from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 with distance downstream is observed in the two-layered vortex street prior to the merging. The frequency spectra of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 are broad-band for 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉–300 but become increasingly sharp-peaked with increasing 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 because the vortex merging process becomes increasingly regular. The ratio of the sharp-peaked frequencies 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190326132244592-0183:S0022112019001678:S0022112019001678_inline11.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 is equal to the ratio of the numbers of vortices observed after and before the merging. The general conclusions drawn from a circular cylinder are expected to be applicable to other bluff bodies.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Results are presented from a laboratory study on the free-surface signal generated over an array of submerged circular cylinders, representative of submerged aquatic vegetation. We aim to understand whether aquatic ecosystems generate a surface signature that is indicative of both what is beneath the water surface as well as how it is altering the flow. A shear layer forms over the canopy, generating coherent vortex structures which eventually manifest in the free-surface slope field. We connect the vortex properties measured at the surface with measurements of the bulk flow, and show that correlations between these quantities are adequate to create a parameterized model in which the interior velocity profile can be predicted solely from measurements taken at the free surface. Experimental surface observations yield a Strouhal number that is twice the most amplified mode predicted by linear stability theory, suggesting that vortices may evolve between generation at the canopy height and their manifestation at the water surface. Additionally, the surface signal continues evolving with distance downstream, with vortices becoming spread farther apart and the passage frequency gradually decreasing. By the trailing edge of the canopy, surface-impacting boils emerge for canopies with higher submergence ratios, suggesting a transition from coherent two-dimensional rollers to transversely varying structures.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Results of experimental and numerical investigations of a supersonic flow around a cylinder with a frontal gas-permeable insert made of a high-porosity cellular material are presented. The measurements are performed in a T-327 supersonic blowdown wind tunnel at the free-stream Mach numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325181423329-0720:S0022112019001654:S0022112019001654_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, 7 and 21 in the range of the unit Reynolds numbers 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/resource/id/urn:cambridge.org:id:binary:20190325181423329-0720:S0022112019001654:S0022112019001654_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉. The drag coefficients for a cylinder with an aerospike and a cylinder with a frontal gas-permeable porous insert are obtained. For the cylinder with the frontal gas-permeable porous insert, variations of the insert length, cylinder diameter and pore size are considered, and the mechanism of drag reduction is found, which includes two supplementary processes: attenuation of the bow shock wave in a system of weaker shock waves, and formation of an effective pointed body. The experiments are accompanied by numerical simulations of the flow around the cylinder with the frontal high-porosity insert: the fields of parameters of the external flow and the flow inside the porous insert are obtained, the drag coefficients are calculated, and the shape of the effective body for the examined model is found. The structure of the high-porosity material is modelled by a system of staggered rings of different diameters aligned in the radial and longitudinal directions (skeleton model of a porous medium). Numerical simulations of the problem are performed by means of solving two-dimensional Reynolds-averaged Navier–Stokes equations written in an axisymmetric form. The experimental and numerical data reveal significant drag reduction in a wide range of supersonic flow conditions. The results obtained on the drag coefficient for the cylinder are generalized with the use of a parameter which includes the ratio of the cylinder diameter to the pore diameter in the insert and the Mach number. This parameter is proposed as a similarity criterion for the problem of a supersonic flow around a cylinder with a frontal high-porosity insert.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Decay of honeycomb-generated turbulence in a duct with a static transverse magnetic field is studied via direct numerical simulations. The simulations follow the revealing experimental study of Sukoriansky 〈span〉et al.〈/span〉 (〈span〉Exp. Fluids〈/span〉, vol. 4 (1), 1986, pp. 11–16), in particular the paradoxical observation of high-amplitude velocity fluctuations, which exist in the downstream portion of the flow when the strong transverse magnetic field is imposed in the entire duct including the honeycomb exit, but not in other configurations. It is shown that the fluctuations are caused by the large-scale quasi-two-dimensional structures forming in the flow at the initial stages of the decay and surviving the magnetic suppression. Statistical turbulence properties, such as the energy decay curves, two-point correlations and typical length scales are computed. The study demonstrates that turbulence decay in the presence of a magnetic field is a complex phenomenon critically depending on the state of the flow at the moment the field is introduced.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Faults in the Roer Valley Rift System (RVRS) act as barriers to horizontal groundwater flow. This causes steep cross-fault groundwater level steps (hydraulic head differences). An overview of the size and distribution of fault-related groundwater level steps and associated fault zone permeabilities is thus far lacking. Such an overview would provide useful insights for nature restoration projects in areas with shallow groundwater levels (〈span〉wijstgronden〈/span〉) on the foot wall of fault zones. In this review study, data on fault zone permeabilities and cross-fault hydraulic head differences were compiled from 39 sources of information, consisting of literature (starting from 1948), internal reports (e.g. from research institutes and drinking water companies), groundwater models, a geological database and new fieldwork. The data are unevenly distributed across the RVRS. Three-quarters of the data sources are related to the Peel Boundary Fault zone (PBFZ). This bias is probably caused by the visibility of fault scarps and fault-adjacent wet areas for the PBFZ, with the characteristic iron-rich groundwater seepage. Most data demonstrate a cross-fault phreatic groundwater level step of 1.0 to 2.5 m. Data for the Feldbiss Fault zone (FFZ) show phreatic cross-fault hydraulic head differences of 1.0 to 2.0 m. 〈span〉In situ〈/span〉 measured hydraulic conductivity data (K) are scarce. Values vary over three orders of magnitude, from 0.013 to 22.1 m d〈span〉−1〈/span〉, are non-directional and do not take into account heterogeneity caused by fault zones. The hydraulic conductivity (and hydraulic resistance) values used in three different groundwater models are obtained by calibration using field measurements. They also cover a large range, from 0.001 to 32 m d〈span〉−1〈/span〉 and from 5 to 100,000 days. Heterogeneity is also not taken into account in these models. The overview only revealed locations with a clear cross-fault groundwater level step, and at many locations the faults are visible on aerial photographs as cropmarks or as soil moisture contrasts at the surface. Therefore, it seems likely that all faults have a reduced permeability, which determines the size of the groundwater level steps. In addition, our results show that cross-fault hydraulic head gradients also correlate with topographic, fault-induced offsets, for both the Peel Boundary and the Feldbiss fault zone.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈img orientation="portrait" mimesubtype="gif" mimetype="image" position="float" type="simple" href="S0022112019008000_figAb" src="http://static.cambridge.org/content/id/urn%3Acambridge.org%3Aid%3Aarticle%3AS0022112019008000/resource/name/S0022112019008000_figAb.gif?pub-status=live"〉〈/p〉〈/div〉 〈div data-abstract-type="normal"〉〈p〉We consider a Stokeslet applied to a viscous fluid next to an infinite, flat wall, or in between two parallel walls. We calculate the forces exerted by the resulting flow on the confining boundaries, and use the results obtained to estimate the hydrodynamic contribution to the pressure exerted on boundaries by force-free self-propelled particles.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈img orientation="portrait" mimesubtype="gif" mimetype="image" position="float" type="simple" href="S0022112019007973_figAb" src="http://static.cambridge.org/content/id/urn%3Acambridge.org%3Aid%3Aarticle%3AS0022112019007973/resource/name/S0022112019007973_figAb.gif?pub-status=live"〉〈/p〉〈/div〉 〈div data-abstract-type="normal"〉〈p〉By means of three-dimensional direct numerical simulations, we investigate the influence of the regular roughness of heated and cooled plates on the mean heat transport in a cylindrical Rayleigh–Bénard convection cell of aspect ratio one. The roughness is introduced by a set of isothermal obstacles, which are attached to the plates and have a form of concentric rings of the same width. The considered Prandtl number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline1.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 equals 1, the Rayleigh number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline2.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 varies from 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline3.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 to 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline4.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, the number of rings on each plate is 1, 2, 4, 8 or 10, the height of the rings is varied from 1.5 % to 49 % of the cylinder height and the gap between the rings is varied from 1.5 % to 18.8 % of the cell diameter. Totally, 135 different cases are analysed. Direct numerical simulations show that with small 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline5.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 and wide roughness rings, a small reduction of the mean heat transport (the Nusselt number 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline6.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉) is possible, but, in most cases, the presence of the heated and cooled obstacles generally leads to an increase of 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline7.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉, compared to the case of classical Rayleigh–Bénard convection with smooth plates. When the rings are very tall and the gaps between them are sufficiently wide, the effective mean heat flux can be several times larger than in the smooth case. For a fixed geometry of the obstacles, the scaling exponent in the 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline8.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 versus 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline9.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 scaling first increases with growing 〈span〉〈span〉〈img data-mimesubtype="gif" data-type="simple" src="http://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20191105072504660-0909:S0022112019007973:S0022112019007973_inline10.gif"〉 〈span data-mathjax-type="texmath"〉 〈/span〉 〈/span〉〈/span〉 up to approximately 0.5, but then smoothly decreases back towards the exponent in the no-obstacle case.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉A stacked aeolian sequence with intercalated soils is presented from the southern Netherlands, which fully covers the Late Weichselian and Holocene periods. An integrated sedimentological (sedimentary structures, grain size), palynological (pollen) and dating approach (radiocarbon, optically stimulated luminescence (OSL)) was applied to unravel climatic and human forcing factors. The dating results of soils and sediments are compatible, and no large hiatuses between the radiocarbon-dated top of the soils and OSL-dated overlying sands were observed. It is argued that the peaty top of wet-type podzols can be used for reliable radiocarbon dating. This study reveals more phases than previously known of landscape stability (Usselo Soil and two podzol soils) and instability (Younger Coversand I and II, two drift-sand units) that are related to Late Weichselian climate change and Holocene human occupation. Regional aeolian deposition in source-bordering (river) dunes (Younger Coversand II) took place in the second part of the Younger Dryas, after 12.3 ka cal. BP, implying a delayed response to Younger Dryas cooling, vegetation cover decline and river pattern change of the Scheldt. The onset of podzolisation and development of ericaceous vegetation occurred prior to the introduction of Neolithic farming, which is earlier than previously assumed. Early podzolisation was followed by a short phase of local drift-sand deposition, at 〈span〉c〈/span〉.5500 cal. BP, that possibly relates to agriculture. Strong human impact on the landscape by deforestation and agriculture resulted in a second phase of widespread drift-sand deposition covering the younger podzol soil after AD 1000.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Early-maturing provitamin A (PVA) quality protein maize (QPM) hybrids with combined drought and low soil nitrogen (low-N) tolerance are needed to address malnutrition and food security problems in sub-Saharan Africa (SSA). The current study's objectives were to (i) examine combining ability of selected early maturing PVA-QPM inbreds for grain yield and other agronomic traits under drought, low-N, optimal environments and across environments, (ii) determine gene action conditioning PVA accumulation under optimal environments, (iii) classify inbreds into heterotic groups and identify testers and (iv) assess yield and stability of hybrids across environments. Ninety-six hybrids generated from 24 inbred lines using the North Carolina Design II together with four commercial hybrid controls were evaluated under drought, low-N and optimal environments in Nigeria in 2016 and 2017. Fifty-four selected hybrids were assayed for PVA carotenoid and tryptophan content. Additive genetic effects were greater than non-additive effects for grain yield and most agronomic traits under each and across environments. The gene action conditioning accumulation of PVA carotenoids under optimal growing conditions followed a pattern similar to that of grain yield and other yield-related traits. The inbred lines were categorized into four heterotic groups consistent with the pedigree records and with TZEIORQ 29 identified as the best male and female tester for heterotic group IV. No tester was found for the other groups. Hybrid TZEIORQ 24 × TZEIORQ 41 was the highest yielding and most stable across environments and should be further tested for consistent performance for commercialization in SSA.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉〈span〉Hordeum marinum〈/span〉 commonly known as sea barley is a salinity-tolerant species of grass. In the current study, 150 lines from ten populations of 〈span〉H. marinum〈/span〉 ssp. 〈span〉marinum〈/span〉 collected from five Tunisian bioclimatic sites were screened for polymorphism with 13 selected random amplified polymorphic DNA primers. Results exhibited a high level of polymorphism (160 polymorphic bands with an average of 12.46 per primer) and a high level of genetic diversity in all the studied populations (on average 〈span〉UHe〈/span〉 = 0.247 and 〈span〉I〈/span〉 = 0.358). High discrimination capacity was found for the 13 primers and a combination of three allowed assignation of a unique profile for each of the 150 lines. The partition of genetic diversity with Analysis of Molecular Variance suggested that the majority of genetic variation (67%) was within populations. The components between-populations within ecoregions and between-ecoregions explained 21 and 12%, respectively, of the total genetic variance. There was no significant association of population differentiation (Ф〈span〉PT〈/span〉) with geographical distance or altitudinal difference. Results also showed that the 150 lines grouped into three clusters with no respect to geographic origin. A sub-set of 13 lines was identified, which captured the maximum genetic diversity of the entire collection. The genetic variation found in this collection of 〈span〉H. marinum〈/span〉 is deemed to be useful in formulating conservation strategies for this species.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Potential planting area for tuber mustard was simulated using the Maxent model under current and future conditions based on 591 coordinates and 22 environmental layers. Model accuracy was excellent, with area under the receiving operator curve values of 0.967 and 0.958 for model training and testing, respectively. Dominant factors were mean diurnal range, mean temperature of the coldest quarter, annual mean temperature and minimum temperature of the coldest month, with thresholds of 6.5–7.5, 5.5–9, 16–19 and 2.0–6.5 °C, respectively. Under current conditions, suitable habitat areas (2.16% of total land in China) were concentrated mainly in Central, Southwest and East China, which can be defined as three occurrence and diffusion centres. In the 2050s and 2070s, suitable habitat areas are predicted to change to 3.72 and 3.92%, and 3.60 and 3.73% under scenarios RCP4.5 and RCP6.0, respectively, indicating that suitable habitat areas will increase slightly. However, future distribution of tuber mustard was predicted to differ among provinces or cities, i.e. predicted suitable habitat areas in Sichuan Province increased up to the 2050s but remained relatively unchanged between the 2050s and 2070s; in Chongqing city they first increased and then decreased; in Hunan, Anhui, Jiangsu, Zhejiang and Fujian Provinces they increased continuously; and in Guizhou, Hubei, Jiangxi Provinces and Shanghai city they first decreased, and then increased. The results from the current study provide useful information for management decisions of tuber mustard.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉A completely randomized experiment was designed to evaluate the effects of 〈span〉α〈/span〉-amylase (AMY) and glucoamylase (GLU) on total losses, fermentative profile, chemical composition and amylolytic activity of rehydrated maize. Eighty-four experimental silos of rehydrated maize [0.33 litres/kg ground maize, 4-mm theoretical particle size, and 625 g/kg dry matter (DM)] were assigned to the following treatments: (1) control (CON), no enzyme addition; (2) GLU added at 300 µl/kg of ground maize (as-fed); and (3) AMY added at 300 µl/kg of ground maize. Seven silos from each treatment were opened after 7, 14, 21 and 28 days. Differences among treatments were evaluated through orthogonal contrasts (CON 〈span〉v.〈/span〉 enzymes, and AMY 〈span〉v.〈/span〉 GLU). Time effects were decomposed using polynomial regression. Glucoamylase silage exhibited greater total losses than AMY. Enzymes increased acetate and lactic acid concentrations and decreased ethanol concentration. Regardless of treatment, gas, effluent and total fermentative losses linearly increased, whereas DM recovery linearly decreased with higher storage length. Glucoamylase silage had lower ammonia nitrogen and higher lactic acid concentrations than AMY. Enzyme treatments decreased silage neutral detergent fibre content and increased 〈span〉in vitro〈/span〉 DM degradation. Glucoamylase silage exhibited a more moderate starch content and greater 〈span〉in vitro〈/span〉 DM degradation than AMY. Storage time linearly decreased DM, starch and fibre content of rehydrated maize. 〈span〉In vitro〈/span〉 degradation of DM linearly increased as the storage length increased. This study showed evidence that enzymes with amylolytic activity, particularly GLU, improve the fermentative profile and DM degradation of rehydrated maize silage.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We argue in this article that the social and economic conditions in the Angkorian society of the tenth century or earlier contributed to the decline in status of some middle-level officials, as is evident from the mid-eleventh century. Many Angkorian inscriptions written between the late ninth and late twelfth centuries record purchases and donations of lands acquired for religious foundations. The texts often contain details of transactions and disputes seeking to validate title to these holdings. The buyers include middle-ranking 〈span〉loñ〈/span〉 and 〈span〉vāp〈/span〉, and increasingly, higher-ranking officials. An analysis of the roles and activities of the officials reveals something of their relative status and helps explain the disappearance of 〈span〉vāp〈/span〉 from the inscriptions in the eleventh century, and the relegation of 〈span〉loñ〈/span〉 to temple roles by the twelfth century. The transfer of communal lands and lands owned by these officials to elites is attributed to hierarchical restrictions on land purchases, a reduction in fiscal immunities, and the need for taxes to be paid to the centre with high-value goods in Angkor's moneyless economy.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The importance of spring grass in the diet of early lactation dairy cows has been widely acknowledged. Numerous studies completed on a plot/paddock basis have identified methods of increasing herbage availability in spring, but little focus was placed on how this impacts animal production. The aim of the current study was to investigate the impact of opening farm cover (OFC; grass availability on farm at turnout in spring; 1036 (high), 748 (medium) and 544 (low) kg DM/ha) and spring rotation length (fast – 56 days and normal – 63 days) on animal performance in early lactation and herbage production and quality. Spring rotation length had little effect on animal performance, herbage variables or sward composition. High OFC increased pre-grazing herbage mass, allowing for increased daily herbage allowance (DHA) compared to medium and low OFC. There was a reduced proportion of leaf in the sward of the high OFC compared to the low OFC, resulting in lower organic matter digestibility. Despite the reduction in sward quality observed as a consequence of achieving high OFC in spring, the greater DHA available to animals increased milk production (+1.4 kg milk/cow/day). Additionally, animals grazing a medium or low OFC had a greater requirement for silage supplementation in spring (+1.3 kg DM/cow/day). The benefits of the higher DHA highlighted in the current study suggest that autumn grazing management must be adapted to increase herbage availability in spring. However, the benefits observed in milk production did not persist beyond the first grazing rotation.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Lord Reid played a vital role as chairman in a Commonwealth commission in framing the Malayan Independence constitution between 1956 and 1957. The Scottish Lord of Appeal sought to ensure the commission's impartiality and to achieve a fair balance between the demands of the various interest groups. The Federation of Malaya was a complex emerging nation-state with a diverse population and the framers had to manage competing interests and demands. This article, through a close examination of the primary constitutional documents, considers Reid's influence on the framing of the Malayan (and hence, Malaysian) federal constitution. The article begins with a brief discussion of Lord Reid's appointment to head the commission and then considers in some depth areas where his influence on the framing of the draft constitution is evident. The article argues that Reid was the main playmaker and moderator during the constitution-framing process and played a critical role in ensuring a balance was achieved between the competing demands of the federal government and the states, safeguarding the fundamental rights of the citizens against the state, and in moderating the various communal demands.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Successive wars and the establishment of a border between the kingdom of Burma and British India in the nineteenth century challenged Burmese conceptions of sovereignty and political space. This essay investigates how European, and more specifically Anglo-American, notions of race, nation, and consular protection to nationals, progressively informed the Burmese concepts of ‘categories of people’ 〈span〉(lumyo)〈/span〉 and ‘subject’ 〈span〉(kyun)〈/span〉. First, I present the semantic evolution of these concepts in the 1820s–1830s, following the annexation of the western Burmese province of Arakan by British India in 1824. Then, I argue that the Burmese concept of 〈span〉lumyo〈/span〉 was progressively associated with the European concept of ‘nations’ in the 1850s–1860s, following the annexation of Lower Burma in 1852. Finally, I uncover developments in the 1870s, when British consular protection extended to several freshly categorised ‘nations’, such as Shan, Karenni, and Kachin.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉This volume contains an introduction and eight papers presented at an international symposium ‘Let's Talk about Trees’, which was organised by Ritsuko Kikusawa and hosted by the National Museum of Ethnology of Osaka, Japan, in February 2013. The stated purpose of the meeting was to evaluate the pros and cons of the classic tree model of historical linguistics in describing the order of splits within a language family. Because the problem of modelling relationships of descent is common to other disciplines, contributors were invited from a range of academic disciplines, including not only linguistics, but also what is described on page one as ‘cladistics’, ‘biology’ and ‘genetics’, although cladistics is clearly a part of biological taxonomy, and not an independent discipline.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The dynamics of a marine ice sheet's grounding lines determine the rate of ice discharge from the grounded part of ice sheet into surrounding oceans. In many locations in West Antarctica ice flows into ice shelves through ice streams experiencing low driving stress. However, existing simple theories of marine ice sheets are developed under the assumption of high basal and driving stress. Here we analyze the grounding line behavior of marine ice streams experiencing low basal shear and driving stress. We find that in this regime, the ice flux at the grounding line is a complex function of the geometry of the ice-stream bed, net accumulation rate and gradient of the net accumulation rate. Our analysis shows that the stability of distinct steady states is determined by the same parameters, suggesting a more complex (in)stability criterion than what is commonly referred to within the context of the ‘marine ice-sheet instability hypothesis’. We also determine characteristic timescales (〈span〉e〈/span〉-folding time) of ice-sheet configurations perturbed from their steady states. These timescales can be used to determine whether particular configurations can be considered in isolation from other components of the climate system or whether their effects and feedbacks between the ice sheet and the rest of the climate system have to be taken into account.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Although archaeological studies focusing on 19th-century sealing have been performed over the past 30 years, its history and sites have traditionally had low visibility in Antarctic narratives and the Antarctic Treaty System policymaking on heritage. Researchers face the challenge of increasing the visibility of sealers’ history and public awareness of the importance of conserving the oldest sites of Antarctica. In this paper, we propose that identifying patterns of tourism activity in the South Shetland Islands, specifically in their temporal and spatial dimensions, could help protect these sites and engage visitors with the early history of Antarctica. Data collected by the International Association of Antarctica Tour Operators were used to calculate landing point usage trends over time and the frequency of passenger landings from 2003–2004 to 2015–2016. We defined six different visitation patterns with temporal tendencies of passenger landings that varied from increasing, constant, or decreasing trends over time, differing in the magnitude and intensity of visitation. This information was used to assess the situation of particular sites located in the vicinity of tourism landing points. We set priorities for their conservation and management decisions and highlighted their relative potential to engage visitors with the stories of 19th-century sealing in Antarctica.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Glaciological ablation is computed from point-scale data at a few ablation stakes that are usually regressed as a function of elevation and averaged over the area-elevation distribution of a glacier. This method is contingent on a tight control of elevation on local ablation. However, in debris-covered glaciers, systematic and random spatial variations of debris thickness modify the ablation rates. We propose and test a method to compute sub-debris ablation where stake data are interpolated as a function of debris-thickness alone and averaged over the debris-thickness distribution at different parts of the glacier. We apply this method on Satopanth Glacier located in Central Himalaya utilising ~1000 ablation measurements obtained from a network of up to 56 stakes during 2015–2017. The estimated mean sub-debris ablation ranges between 1.5±0.2 to 1.7±0.3 cm d〈span〉−1〈/span〉. We show that the debris-thickness-dependent regression describes the spatial variability of the sub-debris ablation better than the elevation dependent regression. The uncertainties in ablation estimates due to the corresponding uncertainties in the measurement of ablation and debris-thickness distribution, and those due to interpolation procedures are estimated using Monte Carlo methods. Possible biases due to a finite number of stakes used are also investigated.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We explore the difficulty of achieving equity for women in two forest and livelihood restoration (FLR) pilot projects, one each in Papua New Guinea (PNG) and the Philippines. We use institutional bricolage as a framework to explain the context and background of stakeholders’ decision-making and the consequent impact on equity and benefit distribution. In the Philippines, material and institutional support was initially successful in assisting participants to establish small-scale tree plantations. A structured approach to institutional development has successfully evolved to meet the needs of women, even though corruption has re-emerged as a destabilizing influence. In PNG, despite success in establishing trees and crops, the participation of women was subjugated to traditional customs and norms that precluded them from engaging in land management decisions. The capacity-building and gender-equity principles of FLR consequently became compromised. We conclude that in some patriarchal societies achieving equity for women will be difficult and progress will be contingent on a detailed understanding of the effects of traditional customs and norms on participation and decision-making.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The terminal Ediacaran Dengying Formation (〈span〉c〈/span〉. 551.1–538.8 Ma) in South China is one of two successions where Ediacara-type macrofossils are preserved in carbonate facies along with skeletal fossils and bilaterian animal traces. Given the remarkable thickness of carbonate-bearing strata deposited in less than 12.3 million years, the Dengying Formation holds the potential for construction of a relatively continuous chemostratigraphic profile for the terminal Ediacaran Period. In this study, a detailed sedimentological and chemostratigraphic (δ〈span〉13〈/span〉C〈span〉carb〈/span〉, δ〈span〉18〈/span〉O〈span〉carb〈/span〉, δ〈span〉13〈/span〉C〈span〉org〈/span〉, δ〈span〉34〈/span〉S〈span〉pyrite〈/span〉, and 〈span〉87〈/span〉Sr/〈span〉86〈/span〉Sr) investigation was conducted on the Dengying Formation at the Gaojiashan section, Ningqiang County of southern Shaanxi Province, South China. Sedimentological results reveal an overall shallow-marine depositional environment. Carbonate breccia, void-filling botryoidal precipitates and aragonite crystal fans are common in the Algal Dolomite Member of the Dengying Formation, suggesting that peritidal facies were repeatedly karstified. The timing of karstification was likely early, probably soon after the deposition of the dolomite sediments. The presence of authigenic aragonite cements suggests high alkalinity in the terminal Ediacaran ocean. Geochemical analysis of micro-drilled samples shows that distinct compositions are registered in different carbonate phases, which should be considered when constructing chemostratigraphic profiles representative of true temporal variations in seawater chemistry. Integrated chemostratigraphic data suggest enhanced burial of organic carbon and pyrite, and the occurrence of extensive marine anoxia (at least in the Gaojiashan Member). Rapid basinal subsidence and carbonate accumulation during a time of elevated seawater alkalinity and increased rates of pyrite burial may have facilitated the evolutionary innovation of early biomineralizing metazoans.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Understanding drivers of deforestation is essential for developing any successful intervention to reduce forest degradation or loss, yet there remains relatively little consensus or clarity on how drivers should be identified and classified. To capture the full range of values and mediating factors that may contribute to land-use behaviours, an approach derived from a shared values perspective that includes a range of values associated with whole landscapes and ecosystems is required. We developed a model that combines behavioural theory with the Capability Approach as a conceptual framework through which to investigate the value–action gap. We used exploratory factor analysis (EFA) of Likert-scale responses to belief statements in order to identify land users’ shared values in the Sarstún Motagua region of Guatemala. We then qualify and quantify the role of capabilities in mediating between the shared values of different cultural groups of land users (Q’eqchi Maya and Ladinos) by comparing their factor scores with their self-reported forest cover change behaviours. Our results indicate that Maya and Ladinos share a set of values, but hold different value orientations that predict their behavioural intentions. We find that their different value orientations reflect behavioural intentions, but an understanding of the capabilities available to different groups is also necessary to fill the value–action gap. These findings have implications for behavioural theory, providing empirical links between shared values, capabilities and behaviour and identification of the role of value orientations, as well as demonstrating a useful approach for decision-makers seeking to understand drivers of change at landscape and whole-ecosystem levels.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Species distribution models (SDMs) are statistical tools used to develop continuous predictions of species occurrence. ‘Integrated SDMs’ (ISDMs) are an elaboration of this approach with potential advantages that allow for the dual use of opportunistically collected presence-only data and site-occupancy data from planned surveys. These models also account for survey bias and imperfect detection through the use of a hierarchical modelling framework that separately estimates the species–environment response and detection process. This is particularly helpful for conservation applications and predictions for rare species, where data are often limited and prediction errors may have significant management consequences. Despite this potential importance, ISDMs remain largely untested under a variety of scenarios. We performed an exploration of key modelling decisions and assumptions on an ISDM using the endangered Baird’s tapir (〈span〉Tapirus bairdii〈/span〉) as a test species. We found that site area had the strongest effect on the magnitude of population estimates and underlying intensity surface and was driven by estimates of model intercepts. Selecting a site area that accounted for the individual movements of the species within an average home range led to population estimates that coincided with expert estimates. ISDMs that do not account for the individual movements of species will likely lead to less accurate estimates of species intensity (number of individuals per unit area) and thus overall population estimates. This bias could be severe and highly detrimental to conservation actions if uninformed ISDMs are used to estimate global populations of threatened and data-deficient species, particularly those that lack natural history and movement information. However, the ISDM was consistently the most accurate model compared to other approaches, which demonstrates the importance of this new modelling framework and the ability to combine opportunistic data with systematic survey data. Thus, we recommend researchers use ISDMs with conservative movement information when estimating population sizes of rare and data-deficient species. ISDMs could be improved by using a similar parameterization to spatial capture–recapture models that explicitly incorporate animal movement as a model parameter, which would further remove the need for spatial subsampling prior to implementation.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We analyze the impact of earthquakes on nighttime lights at a sub-national level, i.e., on grids of different size. We argue that existing studies on the impact of natural disasters on economic development have several important limitations, both at the level of the outcome variable as well as at the level of the independent variable, e.g., the timing of an event and the measuring of its intensity. We aim to overcome these limitations by using geophysical event data on earthquakes together with satellite nighttime lights. Using panel fixed effects regressions covering the entire world for the period 1992–2013, we find that earthquakes reduce both light growth rates and light levels significantly. The effects persist for approximately 5 years, but we find no long-run effects. Effects are stronger the smaller the area of a unit of observation. National institutions and economic conditions are relevant moderating factors.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We recorded the ice motion and icequakes on the floating part of Langhovde Glacier in East Antarctica to better understand the dynamic behavior of ice shelves and floating tongues. Diurnal and semi-diurnal variations in ice motion and seismicity were simultaneously observed at all four global navigation satellite system and three seismic stations over 2 weeks. The short-term along-flow ice motion is explained by the elastic response of the glacier to ocean tide-induced hydrostatic stress variations, which decayed at a rate of 0.8 km〈span〉−1〈/span〉 toward the grounding line. We observed a large number of icequakes during mid-rising and high tides that covered a broad frequency range and formed two major groups of events centered at 10 and 120 Hz, respectively. The hourly occurrence rates were ~500 events h〈span〉−1〈/span〉, with the observed seismicity consistent with fracture due to floating tongue bending. We also observed minor secondary peaks at high ice speeds, which could reflect surface cracking due to stretching or basal friction. Our observation demonstrates that tidal-modulation was the main factor to fracture the floating tongue of Langhovde Glacier.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In Mali's current context, where the crops sector is particularly exposed and vulnerable to agricultural drought, this study assesses the economy-wide impacts of such events and the potential effectiveness of some adaptation strategies. Using a dynamic computable general equilibrium model, we conduct counterfactual simulations of various scenarios accounting for different levels of intensity and frequency of droughts over a 15-year period. We first show how mild, moderate and intense droughts currently experienced by the country affect its economic performance and considerably degrade the welfare of its households. We also show how these negative impacts could be aggravated in the future by the likely increased number of intense droughts threatened by global climate change. However, we finally show that there appears to be some room for Mali to maneuver in terms of drought-risk management policies, such as fostering the use of drought-tolerant crop varieties, improving drought early warning systems or extending irrigation capacities.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Changing rates of water input can affect both the flow of glaciers and ice sheets and their propensity to crevasse. Here we examine geodetic and seismic observations during two substantial (10–18-times background velocity) rain-induced glacier accelerations at Haupapa/Tasman Glacier, New Zealand. Changes in rain rate result in glacier acceleration and associated uplift, which propagate down-glacier. This pattern of acceleration results in a change to the strain rate field, which correlates with an order of magnitude increase in the apparent seismicity rate and an overall down-glacier migration in located seismicity. After each acceleration event the apparent seismicity rate decreases to below the pre-acceleration rate for 3 days. This suggests that seismic events associated with surface crevasse growth occur early during phases of glacier acceleration due to elevated extensional stresses, and then do not occur again until stresses recover.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We compared elastic moduli in polar firn derived from diving wave refraction seismic velocity analysis, firn-core density measurements and microstructure modelling based on firn-core data. The seismic data were obtained with a small electrodynamic vibrator source near Kohnen Station, East Antarctica. The analysis of diving waves resulted in velocity–depth profiles for different wave types (P-, SH- and SV-waves). Dynamic elastic moduli of firn were derived by combining P- and S-wave velocities and densities obtained from firn-core measurements. The structural finite-element method (FEM) was used to calculate the components of the elastic tensor from firn microstructure derived from X-ray tomography of firn-core samples at depths of 10, 42, 71 and 99 m, providing static elastic moduli. Shear and bulk moduli range from 0.39 to 2.42 GPa and 0.68 to 2.42 GPa, respectively. The elastic moduli from seismic observations and the structural FEM agree within 8.5% for the deepest achieved values at a depth of 71 m, and are within the uncertainty range. Our observations demonstrate that the elastic moduli of the firn can be consistently obtained from two independent methods which are based on dynamic (seismic) and static (tomography and FEM) observations, respectively, for deeper layers in the firn below ~10 m depth.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The Yuanjiacun banded iron formation (BIF) is hosted in lower Proterozoic metamorphic strata, and its structures are dominated by bands or streaks. Based on their differences in mineral compositions, the iron ores can be subdivided into haematite quartzite, magnetite quartzite, stilpnomelane magnetite quartzite and stilpnomelane haematite quartzite. The geochemical characteristics of the surrounding rocks show that the protoliths consisted of argillaceous and arenaceous sedimentary rocks. The predominant provenance was a high-maturity felsic sedimentary terrane. The absence of syn-depositional igneous rocks and the tectonic setting discrimination diagrams indicate that the Yuanjiacun BIF formed in a passive continental margin setting. Negligible terrigenous materials were involved in the precipitation of the Yuanjiacun BIF. The precipitation of the Yuanjiacun BIF was predominantly controlled by the mixing of seawater and hydrothermal fluids. Its metallogenic material originated from the leaching of mafic oceanic crust by hydrothermal fluids. The observed Ce anomaly deficiency and heavy Fe isotope enrichment indicate that the Yuanjiacun BIF formed in an anoxic marine environment. In a redox-stratified palaeo-ocean, the Yuanjiacun BIF formed in reducing seawater below the oxidation–reduction transition zone. The Si and O isotope compositions of quartz suggest that the formation of the Yuanjiacun BIF was closely related to submarine hydrothermal activity. The Si and Fe erupted from the seafloor and precipitated by supersaturation and biological oxidation under anoxic conditions, respectively.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉This paper examines and deconstructs the Rio Scale, focusing primarily on the recently published Rio Scale 2.0 concept, from the perspective of a social scientist. I argue that although there is value in developing tools to help astronomers and other scientists communicate their perceptions about the significance of a contact event to the media and the general public, the Rio Scale 2.0 remains problematic conceptually and, thus, does not represent a robust method for assessing or communicating the import of a valid contact. Therefore, it should not be used as a method for informing the media or the general public about scenarios that involve the detection of valid signals suggesting the existence of extraterrestrial intelligence.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉About 400 subglacial lakes are known from Antarctica. The question of whether life unique of subglacial lakes exists has been paramount since their discovery. Despite frequent evidence of microbial life mostly from accretion ice, subglacial lakes are characterized by physiologically hostile conditions to metazoan life, as we know it. Pure water (salinity ≤0.4–1.2%), extreme cold (−3°C), high hydrostatic pressure, areas of limited or no oxygen availability and permanent darkness altogether require physiological adaptations to these harsh conditions. The record of gene sequences including some associated with hydrothermal vents does foster the idea of metazoan life in Lake Vostok. Here, we synthesize the physico-chemical environment surrounding sub-glacial lakes and potential sites of hydrothermal activity and advocate that the physico-chemical stability found at these sites may be the most likely sites for metazoan life to exist. The unique conditions presented by Lake Vostok may also offer an outlook on life to be expected in extra-terrestrial subglacial environments, such as on Jupiter's moon Europa or Saturn's moon Enceladus.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The maximum height trees can grow on Earth is around 122–130 m. The height is constrained by two factors: the availability of water, and where water is not limiting, the pressure available to drive the column of water along the xylem vessels against the pull of gravity (cohesion tension). In turn the height of trees impacts the biodiversity of the environment in a number of ways. On Earth the largest trees are found in maritime temperate environments along the Pacific Northwest coasts of northern California and southern Oregon. These forests provide a large number of secondary habitats for species and serve as moisture pumps that return significant volumes of water to the lower atmosphere. In this work, we apply simple mathematical rules to illustrate how super-terran planets will have significantly smaller trees, with concomitant effects on the habitability of the planet. We also consider the impact of varying tree height on climate models.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Is our current concept of sustainable development too focused on our planet? Do we need a new conception of sustainability? The recent call for an 'ethics of planetary sustainability' by A. Losch may be understood as promoting an affirmative answer to this question. This essay analyses and assesses the introduction of a new concept of sustainability from the perspective of conceptual engineering. The central question is whether this new concept, which I call ‘transplanetary sustainability’, may improve our thinking, and, indirectly, our practices. I argue that a new notion of transplanetary sustainability advantageously points to considerations that matter from a moral point of view. It may also help us to be more exact and consistent in our moral thinking. At the same time, there are serious doubts as to how fundamental the concept of sustainability as such is from a theoretical perspective since it does not figure in prominent moral theories. Furthermore, in view of possible extraterrestrial beings that deserve to be taken into account from a moral point of view, the proposed revision of the concept may reach less than required. Nevertheless, since sustainability has had an impressive career in international politics, it is practically speaking important that sustainability be conceived such that outer space is taken into account.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The frost susceptibility of Australian commercial cereal crops, in particular wheat and barley, has become an economically devastating issue for growers. The relative risk to frost damage of the currently available varieties is obtained through testing varieties in a series of field experiments at locations susceptible to frost events (FEs). The experimental design, measurement protocols and resultant data from these frost expression experiments (FEEs) are complex due to the unpredictability of the timing and severity of FEs, and the maturity of the plants at the time of the events. Design and protocol complexities include the use of multiple sowing dates and the recording of plant maturity. Data difficulties include a high degree of unbalance, and in the instance of multiple frosts in a FEE, there is a longitudinal aspect. A linear mixed model analysis was adopted to accommodate these characteristics of individual FEEs and the multi-environment trial analysis of 17 FEEs. Finally, an approach is demonstrated for dissemination of results that are of use to both growers and breeders.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉This article provides an outline of the historical construction of Timorese (East Timorese and Indonesian West Timorese) geo-bodies and communal identities from the mid-nineteenth century to the present time, thereby reconstructing the origins of many national imaginings amongst the Timorese people. Since the controversial annexation of Portuguese Timor by Indonesia in 1976, (East) Timor has been constructed as a place of two territorial identities: Timor as a part of Indonesia and East Timor as a homogeneous nation distinct from Indonesia. However, representations of Timor had been much more fluid and inconsistent in preceding ages. This article studies various communities’ representations of Timor to reveal dialectic relations between diverse colonial and post-colonial representations of the Timorese spaces and their senses of belonging. Thereby, it problematises the political role of global and regional place-making in a contested Southeast Asian locale.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉This article examines 〈span〉gotong royong〈/span〉 as social citizenship in the second half of the twentieth century in Indonesia. 〈span〉Gotong royong〈/span〉 was traditionally understood as a collective spirit among neighbours to strengthen economic and social resilience. However, the institutionalisation of the system through massive development programmes since the 1940s has influenced the common people's perception and practice of it. This article argues that the co-option of 〈span〉gotong royong〈/span〉 as part of the discourse of nation-state building transformed the popular view of labour and capital mobilisation from an openly critical stance to apparent submission. 〈span〉Gotong royong〈/span〉 became a means of social engineering and an ingenious linguistic strategy by which state elites orchestrated control over the making of citizens. Still, Indonesians have remained vigilant of their participation in 〈span〉gotong royong〈/span〉, and about the potential for its misuse by corrupt officials.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Professional scientist-geographer Erich von Drygalski led the first German expedition to Antarctica in 1901–1903. The expedition saw itself as purely scientific, which turned out to be at odds with the expectations of Imperial Germany at the time. It was one of the first to use photography extensively and effectively to document and record scientific activities and to shape the public’s image of the work that was being done in this remote and unknown part of the world. Ice was the leitmotif of Drygalski’s life. He had prior experience in the Arctic, and the year spent in Antarctica confirmed his nuanced way of viewing the ice: on the one hand, and foremost, scholarly and objective, while still appreciating its aesthetic qualities; on the other, infused with feelings of human vulnerability. Using discourse analysis, this article examines Drygalski’s published work and photographs he chose to illustrate it, in order to investigate what the ice meant to him. In his writings, it was the scholarly, objective attitude which predominated and this may have contributed to the generally lacklustre reception of his Antarctic achievements. The photographs he chose to illustrate his published work, however, were many and varied, often capturing the awe-inspiring beauty of the ice and contributing to good sales of his narrative of the South Polar Expedition.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Global warming and the associated glacier retreat recently revealed the entrance to an ice–rock tunnel, at an altitude of ~3600 m a.s.l., in the uppermost portion of the Forni Glacier in the Central Italian Alps. The tunnel served as an entrance to an Austro-Hungarian cableway station excavated in the rocks during the Great War just behind the frontline. A comprehensive geophysical survey, based on seismic and ground-penetrating radar profiling, was then undertaken to map other possible World War I (WWI) remains still embedded in the ice. The ice–rock interface was reconstructed over the entire saddle and in the uppermost portion of the glacier. A prominent linear reflector was surprisingly similar to the common response of buried pipes. The reflector orientation, almost longitudinal to the slope, does not seem to be compatible with a glacial conduit or with other natural features. Numerical simulations of a series of possible targets constrained interpretation to a partly water-filled rounded shape cavity. The presence of a preserved WWI tunnel connecting Mount Vioz and Punta Linke could be considered a realistic hypothesis. The Forni glacier could be still considered polythermal and comprised of cold ice without basal sliding in its top portion.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Analysis of a recent surge of Morsnevbreen, Svalbard, is used to test predictions of the enthalpy balance theory of surging. High-resolution time series of velocities, ice thickness and crevasse distribution allow key elements of the enthalpy (internal energy) budget to be quantified for different stages of the surge cycle. During quiescence (1936–1990), velocities were very low, and geothermal heat slowly built-up enthalpy at the bed. Measurable mass transfer and frictional heating began in 1990–2010, then positive frictional heating-velocity feedbacks caused gradual acceleration from 2010 to 2015. Rapid acceleration occurred in summer 2016, when extensive crevassing and positive air temperatures allowed significant surface to bed drainage. The surge front reached the terminus in October 2016, coincident with a drop in velocities. Ice plumes in the fjord are interpreted as discharge of large volumes of supercooled water from the bed. Surge termination was prolonged, however, indicating persistence of an inefficient drainage system. The observations closely match predictions of the theory, particularly build-up of enthalpy from geothermal and frictional heat, and surface meltwater, and the concomitant changes in ice-surface elevation and velocity. Additional characteristics of the surge reflect spatial processes not represented in the model, but can be explained with respect to enthalpy gradients.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We determined ice velocities for the Academy of Sciences Ice Cap, Severnaya Zemlya, Russian Arctic, during November 2016–November 2017, by feature-tracking 54 pairs of Sentinel-1 synthetic-aperture radar images. Seasonal velocity variations with amplitudes up to 10% of the yearly-averaged velocity were observed. Shorter-term (−1. The difference from an estimate of ~ 1.4 Gt a〈span〉−1〈/span〉 for 2003–2009 was attributed to the initiation of ice stream flow in Basin BC. The total geodetic mass balance over 2012–2016 was − 1.72 ± 0.67 Gt a〈span〉−1〈/span〉 (− 0.31 ± 0.12 m w.e. a〈span〉−1〈/span〉). The climatic mass balance was not significantly different from zero, at 0.21 ± 0.68 Gt a〈span〉−1〈/span〉 (0.04 ± 0.12 m w.e. a〈span〉−1〈/span〉), and has remained near zero at decadal-scale for the last four decades. Therefore, the total mass balance has been controlled largely by variations in ice discharge, whose long-term changes do not appear to have responded to environmental changes but to the intrinsic characteristics of the ice cap governing tidewater glacier dynamics.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We compared systematic and random survey techniques to estimate breeding population sizes of burrow-nesting petrel species on Marion Island. White-chinned (〈span〉Procellaria aequinoctialis〈/span〉) and blue (〈span〉Halobaena caerulea〈/span〉) petrel population sizes were estimated in systematic surveys (which attempt to count every colony) in 2009 and 2012, respectively. In 2015, we counted burrows of white-chinned, blue and great-winged (〈span〉Pterodroma macroptera〈/span〉) petrels within 52 randomized strip transects (25 m wide, total 144 km). Burrow densities were extrapolated by Geographic Information System-derived habitat attributes (geology, vegetation, slope, elevation, aspect) to generate island-wide burrow estimates. Great-winged petrel burrows were found singly or in small groups at low densities (2 burrows ha〈span〉−1〈/span〉); white-chinned petrel burrows were in loose clusters at moderate densities (3 burrows ha〈span〉−1〈/span〉); and blue petrel burrows were in tight clusters at high densities (13 burrows ha〈span〉−1〈/span〉). The random survey estimated 58% more white-chinned petrels but 42% fewer blue petrels than the systematic surveys. The results suggest that random transects are best suited for species that are widely distributed at low densities, but become increasingly poor for estimating population sizes of species with clustered distributions. Repeated fixed transects provide a robust way to monitor changes in colony density and area, but might fail to detect the formation/disappearance of new colonies.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉The Southern Patagonia Icefield (SPI) withdrawal in recent decades shows contrasting behaviours between adjacent basins. One of the basins with highest volumetric losses is located at northernmost SPI. We refer to Jorge Montt tidewater glacier (48〈span〉°〈/span〉 30′S/73〈span〉°〈/span〉 30′W, 445 km〈span〉2〈/span〉 in 2018), which retreated 2.7 km between 2011 and 2018 and thinned at rates of up to 21 m a〈span〉−1〈/span〉 over this period. Based on the retreat record, remote-sensing imagery, field data, a mass-balance model and a calving parameterisation, we attempted to differentiate climatic-induced changes (i.e. surface mass balance) and dynamic responses (i.e. calving fluxes). The surface mass balance reached −4.15 km〈span〉3〈/span〉 w.e. a〈span〉−1〈/span〉 between 2012 and 2017. When frontal ablation is included, the net mass balance is −17.79 km〈span〉3〈/span〉 w.e. a〈span〉−1〈/span〉. This represents a change of trend compared with modelling estimations of positive surface mass balance prior to 2010. This shift is attributed to higher ablation rates given that accumulation is known to have increased between 1980 and 2015. The available evidence, therefore, indicates that frontal ablation is the main factor, supported by observed rates at Jorge Montt as high as 3.81 km〈span〉3〈/span〉 w.e. a〈span〉−1〈/span〉 in 2015, with ice velocities peaking at 11 km a〈span〉−1〈/span〉.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉In an effort to improve the reliability of Arctic sea-ice predictions, an ensemble-based Arctic Ice Ocean Prediction System (ArcIOPS) has been developed to meet operational demands. The system is based on a regional Arctic configuration of the Massachusetts Institute of Technology general circulation model. A localized error subspace transform ensemble Kalman filter is used to assimilate the weekly merged CryoSat-2 and Soil Moisture and Ocean Salinity sea-ice thickness data together with the daily Advanced Microwave Scanning Radiometer 2 (AMSR2) sea-ice concentration data. The weather forecasts from the Global Forecast System of the National Centers for Environmental Prediction drive the sea ice–ocean coupled model. The ensemble mean sea-ice forecasts were used to facilitate the Chinese National Arctic Research Expedition in summer 2017. The forecasted sea-ice concentration is evaluated against AMSR2 and Special Sensor Microwave Imager/Sounder sea-ice concentration data. The forecasted sea-ice thickness is compared to the in-situ observations and the Pan-Arctic Ice-Ocean Modeling and Assimilation System. These comparisons show the promising potential of ArcIOPS for operational Arctic sea-ice forecasts. Nevertheless, the forecast bias in the Beaufort Sea calls for a delicate parameter calibration and a better design of the assimilation system.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉We used fossil Chironomidae assemblages and the transfer function approach to reconstruct summer air temperatures over the past 300 years from a High Arctic lake in Hornsund, Svalbard. Our aims were to compare reconstructed summer temperatures with observed (last 100 years) seasonal temperatures, to determine a potential climate warming break point in the temperature series and to assess the significance and rate of the climate warming trend at the study site. The reconstructed temperatures were consistent with a previous proxy record from Svalbard and showed good correlation with the meteorological observations from Bjørnøya and Longyearbyen. From the current palaeoclimate record, we found a significant climate warming threshold in the 1930s, after which the temperatures rapidly increased. We also found that the climate warming trend was strong and statistically significant. Compared with the reconstructed Little Ice Age temperatures in late eighteenth century cooling culmination, the present day summer temperatures are 〉4°C higher and the temperature increase since the 1930s has been 0.5°C per decade. These results highlight the exceptionally rapid recent warming of southern Svalbard and add invaluable information on the seasonality of High Arctic climate change and Arctic amplification.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉Here we describe new microfossil assemblages for the Miocene Hobbs Glacier Formation and the first possibly indigenous assemblages for the Plio-Pleistocene Weddell Sea Formation on Seymour Island, West Antarctica. The assemblages are composed mainly of foraminifers, but radiolarians, calcitarchs and poriferan sclerites are also present. For the Hobbs Glacier Formation, we report the foraminifers 〈span〉Bolivina〈/span〉 sp., 〈span〉Oolina globosa〈/span〉 and 〈span〉Rosalina〈/span〉 cf. 〈span〉globularis〈/span〉; and for the Weddell Sea Formation, we report 〈span〉Favulina hexagona〈/span〉, 〈span〉Globigerinita uvula〈/span〉, 〈span〉Globocassidulina〈/span〉 cf. 〈span〉subglobosa〈/span〉 and 〈span〉Psammosphaera fusca〈/span〉. The low abundance and diversity of microfossils, allied with the complex taphonomical processes that prevailed in Antarctic glacial–marine palaeoenvironments, make it impossible to define whether the assemblages are composed of a mixture of indigenous and re-elaborated specimens or exclusively of re-elaborated remains. Nevertheless, the indigenous nature of some specimens is suggested by their inherent fragility, excellent preservation and/or taxonomic association with indigenous assemblages from correlated strata. The taxonomic compositions are not directly comparable with other Antarctic assemblages, although most of the species were previously reported from pre-Quaternary or modern deposits of both West and East Antarctica. This lack of correspondence is probably due to preservation biases, but any further significance is hidden by the complex taphonomy of the deposits.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉This study evaluates the effects of vegetative soil conservation practices (afforestation and/or bamboo planting) on farm profit and its components, revenue and variable cost. Since farmers self-select themselves as adopters of conservation measures, there could be a problem of selection bias in evaluating their soil conservation practices. We address the selection bias by using propensity score matching. We also check if there exists spatial spillover in adoption of vegetative conservation measures and how it affects matching. We use primary survey data from the Darjeeling district of the Eastern Himalayan region for the year 2013. Our results suggest strong spatial correlation. We find that the propensity score estimated from the spatial model provides better matches than the non-spatial model. While the results show that vegetative soil conservation can lead to significant gains in revenue, it also increases costs so that no significant gains in profit accrue to farmers.〈/p〉〈/div〉

Description: 〈div data-abstract-type="normal"〉〈p〉A survey of Antarctic toothfish (〈span〉Dissostichus mawsoni〈/span〉) was conducted in the northern Ross Sea region during the winter of 2016 to document the timing and location of spawning activity, to collect biological information about reproductive status during the spawning season and to look for temporal signals in biological data from 〈span〉D. mawsoni〈/span〉 that may indicate a spawning migration of mature toothfish from the continental slope region to the northern Ross Sea region. The 58 day survey showed that spawning of 〈span〉D. mawsoni〈/span〉 began on some seamounts by early July. No changes were detected between winter and summer in length, age, sex ratio or condition factor distributions for 〈span〉D. mawsoni〈/span〉 in the northern Ross Sea as hypothesized following a spawning migration from the slope to the northern Ross Sea region. These results suggest that the distribution of 〈span〉D. mawsoni〈/span〉 in the Ross Sea is mainly accomplished through ontogenetic migration and not annual return spawning migrations.〈/p〉〈/div〉