ALBERT

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: 〈b〉Ensemble forecasts of air quality in eastern China – Part 1: Model description and implementation of the MarcoPolo–Panda prediction system, version 1〈/b〉〈br〉 Guy P. Brasseur, Ying Xie, Anna Katinka Petersen, Idir Bouarar, Johannes Flemming, Michael Gauss, Fei Jiang, Rostislav Kouznetsov, Richard Kranenburg, Bas Mijling, Vincent-Henri Peuch, Matthieu Pommier, Arjo Segers, Mikhail Sofiev, Renske Timmermans, Ronald van der A, Stacy Walters, Jianming Xu, and Guangqiang Zhou〈br〉 Geosci. Model Dev., 12, 33-67, https://doi.org/10.5194/gmd-12-33-2019, 2019〈br〉 An operational multi-model forecasting system for air quality provides daily forecasts of ozone, nitrogen oxides, and particulate matter for 37 urban areas in China. The paper presents an intercomparison of the different forecasts performed during a specific period of time and highlights recurrent differences between the model output. Pathways to improve the forecasts by the multi-model system are suggested.

Description: 〈b〉RandomFront 2.3: a physical parameterisation of fire spotting for operational fire spread models – implementation in WRF-SFIRE and response analysis with LSFire+〈/b〉〈br〉 Andrea Trucchia, Vera Egorova, Anton Butenko, Inderpreet Kaur, and Gianni Pagnini〈br〉 Geosci. Model Dev., 12, 69-87, https://doi.org/10.5194/gmd-12-69-2019, 2019〈br〉 Wildfires are a concrete problem and impact on human life, property and the environment. An extremely dangerous phenomenon is so-called 〈q〉fire spotting〈/q〉, i.e., the generation of secondary ignitions responsible for dangerous flare-ups during wildfires. The aim of this research was to improve the tools used for risk management through the inclusion of fire spotting in operational wildfire simulators used by forest service agencies.

Description: 〈b〉Vertical characterization of aerosol optical properties and brown carbon in winter in urban Beijing, China〈/b〉〈br〉 Conghui Xie, Weiqi Xu, Junfeng Wang, Qingqing Wang, Dantong Liu, Guiqian Tang, Ping Chen, Wei Du, Jian Zhao, Yingjie Zhang, Wei Zhou, Tingting Han, Qingyun Bian, Jie Li, Pingqing Fu, Zifa Wang, Xinlei Ge, James Allan, Hugh Coe, and Yele Sun〈br〉 Atmos. Chem. Phys., 19, 165-179, https://doi.org/10.5194/acp-19-165-2019, 2019〈br〉 We present the first simultaneous real-time online measurements of aerosol optical properties at ground level and at 260 m on a meteorological tower in urban Beijing in winter. The vertical similarities and differences in scattering and absorption coefficients were characterized. The increases in MAC of BC were mainly associated with the coating materials on rBC. Coal combustion was the dominant source contribution of brown carbon followed by biomass burning and SOA in winter in Beijing.

Description: 〈b〉Dynamic changes in optical and chemical properties of tar ball aerosols by atmospheric photochemical aging〈/b〉〈br〉 Chunlin Li, Quanfu He, Julian Schade, Johannes Passig, Ralf Zimmermann, Daphne Meidan, Alexander Laskin, and Yinon Rudich〈br〉 Atmos. Chem. Phys., 19, 139-163, https://doi.org/10.5194/acp-19-139-2019, 2019〈br〉 〈p〉Following wood pyrolysis, tar ball aerosols were laboratory generated from wood tar separated into polar and nonpolar phases. Chemical information of fresh tar balls was obtained from a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) and single-particle laser desorption/resonance enhanced multiphoton ionization mass spectrometry (SP-LD-REMPI-MS). Their continuous refractive index (RI) between 365 and 425 nm was retrieved using a broadband cavity enhanced spectroscopy (BBCES). Dynamic changes in the optical and chemical properties for the nonpolar tar ball aerosols in 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical process were investigated in an oxidation flow reactor (OFR). Distinct differences in the chemical composition of the fresh polar and nonpolar tar aerosols were identified. Nonpolar tar aerosols contain predominantly high-molecular weight unsubstituted and alkyl-substituted polycylic aromatic hydrocarbons (PAHs), while polar tar aerosols consist of a high number of oxidized aromatic substances (e.g., methoxy-phenols, benzenediol) with higher O : C ratios and carbon oxidation states. Fresh tar balls have light absorption characteristics similar to atmospheric brown carbon (BrC) aerosol with higher absorption efficiency towards the UV wavelengths. The average retrieved RI is 〈span〉1.661+0.020〈i〉i〈/i〉〈/span〉 and 〈span〉1.635+0.003〈i〉i〈/i〉〈/span〉 for the nonpolar and polar tar aerosols, respectively, with an absorption Ångström exponent (AAE) between 5.7 and 7.8 in the detected wavelength range. The RI fits a volume mixing rule for internally mixed nonpolar/polar tar balls. The RI of the tar ball aerosols decreased with increasing wavelength under photochemical oxidation. Photolysis by UV light (254 nm), without strong oxidants in the system, slightly decreased the RI and increased the oxidation state of the tar balls. Oxidation under varying OH exposure levels and in the absence of 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 diminished the absorption (bleaching) and increased the O : C ratio of the tar balls. The photobleaching via OH radical initiated oxidation is mainly attributed to decomposition of chromophoric aromatics, nitrogen-containing organics, and high-molecular weight components in the aged particles. Photolysis of nitrous oxide (〈span〉N〈sub〉2〈/sub〉O〈/span〉) was used to simulate 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical aging of tar balls in the OFR. Under high-〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 conditions with similar OH exposure, photochemical aging led to the formation of organic nitrates, and increased both oxidation degree and light absorption for the aged tar ball aerosols. These observations suggest that secondary organic nitrate formation counteracts the bleaching by OH radical photooxidation to eventually regain some absorption of the aged tar ball aerosols. The atmospheric implication and climate effects from tar balls upon various oxidation processes are briefly discussed.〈/p〉

Description: 〈b〉Atmospheric new particle formation in China〈/b〉〈br〉 Biwu Chu, Veli-Matti Kerminen, Federico Bianchi, Chao Yan, Tuukka Petäjä, and Markku Kulmala〈br〉 Atmos. Chem. Phys., 19, 115-138, https://doi.org/10.5194/acp-19-115-2019, 2019〈br〉 The characteristics of new particle formation (NPF) in China, including frequency, formation rate, and particle growth rate, were summarized comprehensively and were compared among observations in different environments. The interactions between air pollution and NPF are discussed, as well as the possible reasons for more frequent NPF under heavy pollution conditions than in our current understanding. Significant and future research directions for NPF studies in China are also summarized.

Description: 〈b〉Dynamic changes in optical and chemical properties of tar ball aerosols by atmospheric photochemical aging〈/b〉〈br〉 Chunlin Li, Quanfu He, Julian Schade, Johannes Passig, Ralf Zimmermann, Daphne Meidan, Alexander Laskin, and Yinon Rudich〈br〉 Atmos. Chem. Phys., 19, 139-163, https://doi.org/10.5194/acp-19-139-2019, 2019〈br〉 〈p〉Following wood pyrolysis, tar ball aerosols were laboratory generated from wood tar separated into polar and nonpolar phases. Chemical information of fresh tar balls was obtained from a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) and single-particle laser desorption/resonance enhanced multiphoton ionization mass spectrometry (SP-LD-REMPI-MS). Their continuous refractive index (RI) between 365 and 425 nm was retrieved using a broadband cavity enhanced spectroscopy (BBCES). Dynamic changes in the optical and chemical properties for the nonpolar tar ball aerosols in 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical process were investigated in an oxidation flow reactor (OFR). Distinct differences in the chemical composition of the fresh polar and nonpolar tar aerosols were identified. Nonpolar tar aerosols contain predominantly high-molecular weight unsubstituted and alkyl-substituted polycylic aromatic hydrocarbons (PAHs), while polar tar aerosols consist of a high number of oxidized aromatic substances (e.g., methoxy-phenols, benzenediol) with higher O : C ratios and carbon oxidation states. Fresh tar balls have light absorption characteristics similar to atmospheric brown carbon (BrC) aerosol with higher absorption efficiency towards the UV wavelengths. The average retrieved RI is 〈span〉1.661+0.020〈i〉i〈/i〉〈/span〉 and 〈span〉1.635+0.003〈i〉i〈/i〉〈/span〉 for the nonpolar and polar tar aerosols, respectively, with an absorption Ångström exponent (AAE) between 5.7 and 7.8 in the detected wavelength range. The RI fits a volume mixing rule for internally mixed nonpolar/polar tar balls. The RI of the tar ball aerosols decreased with increasing wavelength under photochemical oxidation. Photolysis by UV light (254 nm), without strong oxidants in the system, slightly decreased the RI and increased the oxidation state of the tar balls. Oxidation under varying OH exposure levels and in the absence of 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 diminished the absorption (bleaching) and increased the O : C ratio of the tar balls. The photobleaching via OH radical initiated oxidation is mainly attributed to decomposition of chromophoric aromatics, nitrogen-containing organics, and high-molecular weight components in the aged particles. Photolysis of nitrous oxide (〈span〉N〈sub〉2〈/sub〉O〈/span〉) was used to simulate 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical aging of tar balls in the OFR. Under high-〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 conditions with similar OH exposure, photochemical aging led to the formation of organic nitrates, and increased both oxidation degree and light absorption for the aged tar ball aerosols. These observations suggest that secondary organic nitrate formation counteracts the bleaching by OH radical photooxidation to eventually regain some absorption of the aged tar ball aerosols. The atmospheric implication and climate effects from tar balls upon various oxidation processes are briefly discussed.〈/p〉

Description: 〈b〉Watershed classification for the Canadian prairie〈/b〉〈br〉 Jared D. Wolfe, Kevin R. Shook, Chris Spence, and Colin J. Whitfield〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-625,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 Watershed classification can identify regions expected to respond similarly to disturbance. Methods should extend beyond hydrology to include other environmental questions, such as ecology and water quality. We developed a classification for the Canadian Prairie, and identified seven classes defined by watershed characteristics, including elevation, climate, wetland density, and surficial geology. Results provide a basis for evaluating watershed response to land management and climate condition.

Description: 〈b〉Implications of water management representations for watershed hydrologic modeling in the Yakima River basin〈/b〉〈br〉 Jiali Qiu, Qichun Yang, Xuesong Zhang, Maoyi Huang, Jennifer C. Adam, and Keyvan Malek〈br〉 Hydrol. Earth Syst. Sci., 23, 35-49, https://doi.org/10.5194/hess-23-35-2019, 2019〈br〉 Complex water management activities challenge hydrologic modeling. We evaluated how different representations of reservoir operation and agricultural irrigation affect streamflow simulations in the Yakima River basin. Results highlight the importance of the inclusion of reliable reservoir and irrigation information in watershed models for improving watershed hydrology modeling. Models used here are public and hold the promise to benefit water assessment and management in other basins.

Description: 〈b〉Aircraft-based measurements of High Arctic springtime aerosol show evidence for vertically varying sources, transport and composition〈/b〉〈br〉 Megan D. Willis, Heiko Bozem, Daniel Kunkel, Alex K. Y. Lee, Hannes Schulz, Julia Burkart, Amir A. Aliabadi, Andreas B. Herber, W. Richard Leaitch, and Jonathan P. D. Abbatt〈br〉 Atmos. Chem. Phys., 19, 57-76, https://doi.org/10.5194/acp-19-57-2019, 2019〈br〉 The vertical distribution of Arctic aerosol is an important driver of its climate impacts. We present vertically resolved measurements of aerosol composition and properties made in the High Arctic during spring on an aircraft platform. We explore how aerosol properties are related to transport history and show evidence of vertical trends in aerosol sources, transport mechanisms and composition. These results will help us to better understand aerosol–climate interactions in the Arctic.

Description: 〈b〉Estimating the soil N〈sub〉2〈/sub〉O emission intensity of croplands in northwest Europe〈/b〉〈br〉 Vasileios Myrgiotis, Mathew Williams, Robert M. Rees, and Cairistiona F. E. Topp〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-490,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 This study focuses on a northwest European cropland region and shows that the type of crop growing on a soil has notable effects on the emission of nitrous oxide (N〈sub〉2〈/sub〉O – a greenhouse gas) from that soil. It was found that N〈sub〉2〈/sub〉O emissions from soils under oilseed cultivation are significantly higher than soils under cereal cultivation. This variation is mostly explained by the fact that oilseeds require more nitrogen (fertiliser) than cereals especially at early crop growth stages.

Description: 〈b〉A Double ITCZ Phenomenology of Wind Errors in the Equatorial Atlantic in Seasonal Forecasts with ECMWF Models〈/b〉〈br〉 Jonathan K. P. Shonk, Teferi D. Demissie, and Thomas Toniazzo〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1316,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 Modern climate models are affected by systematic biases that harm their ability to produce reliable seasonal forecasts and climate projections. In this study, we investigate causes of biases in wind patterns over the tropical Atlantic during northern spring in three related models. We find that the wind biases are associated with an increase in excess rainfall and convergence in the tropical western Atlantic at the start of April, leading to the redirection of trade winds away from the equator.

Description: 〈b〉Spiky Fluctuations and Scaling in High-Resolution EPICA Ice Core Dust Fluxes〈/b〉〈br〉 Shaun Lovejoy and Fabrice Lambert〈br〉 Clim. Past Discuss., https//doi.org/10.5194/cp-2018-171,2019〈br〉 〈b〉Manuscript under review for CP〈/b〉 (discussion: open, 0 comments)〈br〉 We analyze the statistical properties of the eight past glacial-interglacial cycles as well as subsections of a generic glacial cycle using the high-resolution dust flux dataset from the Antarctic EPICA Dome C ice core. We show that the high southern latitude climate during glacial maxima, interglacial, and glacial inception is generally more stable but more drought-prone than during mid-glacial conditions.

Description: 〈b〉Nutrient distribution and nitrogen and oxygen isotopic composition of nitrate in water masses of the subtropical South Indian Ocean〈/b〉〈br〉 Natalie C. Harms, Niko Lahajnar, Birgit Gaye, Tim Rixen, Kirstin Dähnke, Markus Ankele, Ulrich Schwarz-Schampera, and Kay-Christian Emeis〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-511,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 The Indian Ocean subtropical gyre is a large oligotrophic area that is likely to adjust to continued warming by increasing stratification, reduced nutrient supply, and decreasing biological production. In this study, we investigated concentrations of nutrients and stable isotopes of nitrate. We determine the lateral influence of water masses entering the gyre from the northern Indian Ocean and from the Southern Ocean and quantify the input of nitrogen by N〈sub〉2〈/sub〉-fixation into the surface layer.

Description: 〈b〉New particle formation in the marine atmosphere during seven cruise campaigns〈/b〉〈br〉 Yujiao Zhu, Kai Li, Yanjie Shen, Yang Gao, Xiaohuan Liu, Yang Yu, Huiwang Gao, and Xiaohong Yao〈br〉 Atmos. Chem. Phys., 19, 89-113, https://doi.org/10.5194/acp-19-89-2019, 2019〈br〉 In this paper, we investigate new particle formation (NPF) events during seven cruises. NPF events were observed on 25 days and were most likely associated with the long-range transport of anthropogenic air pollutants. The relationship between the net generated amount of new particles and their apparent formation rate is established and explained in terms of the roles of different vapor precursors. The survival probability of new particles to CCN size is also discussed.

Description: 〈b〉Aircraft-based measurements of High Arctic springtime aerosol show evidence for vertically varying sources, transport and composition〈/b〉〈br〉 Megan D. Willis, Heiko Bozem, Daniel Kunkel, Alex K. Y. Lee, Hannes Schulz, Julia Burkart, Amir A. Aliabadi, Andreas B. Herber, W. Richard Leaitch, and Jonathan P. D. Abbatt〈br〉 Atmos. Chem. Phys., 19, 57-76, https://doi.org/10.5194/acp-19-57-2019, 2019〈br〉 The vertical distribution of Arctic aerosol is an important driver of its climate impacts. We present vertically resolved measurements of aerosol composition and properties made in the High Arctic during spring on an aircraft platform. We explore how aerosol properties are related to transport history and show evidence of vertical trends in aerosol sources, transport mechanisms and composition. These results will help us to better understand aerosol–climate interactions in the Arctic.

Description: 〈b〉Simulated single-layer forest canopies delay Northern Hemisphere snowmelt〈/b〉〈br〉 Markus Todt, Nick Rutter, Christopher G. Fletcher, and Leanne M. Wake〈br〉 The Cryosphere Discuss., https//doi.org/10.5194/tc-2018-270,2019〈br〉 〈b〉Manuscript under review for TC〈/b〉 (discussion: open, 0 comments)〈br〉 Vegetation is often represented by a single layer in global land models. Studies have found deficient simulation of thermal radiation beneath forest canopies when represented by single-layer vegetation. This study corrects thermal radiation in forests for a global land model using single-layer vegetation in order to assess the effect of deficient thermal radiation on snow cover and snowmelt. Results indicate that single-layer vegetation causes snow in forests being too cold and melting too late.

Description: 〈b〉Brief Communication: Early season snowpack loss and implications for oversnow vehicle recreation travel planning〈/b〉〈br〉 Benjamin J. Hatchett and Hilary G. Eisen〈br〉 The Cryosphere, 13, 21-28, https://doi.org/10.5194/tc-13-21-2019, 2019〈br〉 We examine the timing of early season snowpack relevant to oversnow vehicle (OSV) recreation over the past 3 decades in the Lake Tahoe region (USA). Data from two independent data sources suggest that the timing of achieving sufficient snowpack has shifted later by 2 weeks. Increasing rainfall and more dry days play a role in the later onset. Adaptation strategies are provided for winter travel management planning to address negative impacts of loss of early season snowpack for OSV usage.

Description: 〈b〉Aerosol radiative effects with MACv2〈/b〉〈br〉 Stefan Kinne〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-949,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 〈p〉onthly global maps for aerosol properties of the MACv2 climatology are applied in an off-line radiative transfer model to determine aerosol radiative effects. For details beyond global averages in most cases global maps are presented to visualize regional and seasonal details. Aside from the direct radiative (aerosol presence) effect, including those for aerosol components as extracted from MACv2 aerosol optics, also the major aerosol indirect radiative effect is covered. Hereby, the impact of smaller drops in water clouds due to added anthropogenic aerosol was simulated by applying a satellite retrieval based fit from locally associations between aerosol and drop concentrations over oceans. Present-day anthropogenic aerosols of MACv2 – on a global average basis – reduce the radiative net-fluxes at the top of the atmosphere (TOA) by −1.0 W/m〈sup〉2〈/sup〉 and at the surface by −2.1 W/m〈sup〉2〈/sup〉. Direct cooling contributions are only about half of indirect contributions (−.35 vs −.65) at TOA, but about twice at the surface (−1.45 vs −.65), as solar absorption of the direct effect warms the atmosphere by +1.1 W/m〈sup〉2〈/sup〉. Natural aerosols are on average less absorbing (for a relatively larger solar TOA cooling) and larger in size (now contributing with IR greenhouse warming). Thus, average TOA direct forcing efficiencies for total and anthropogenic aerosol happen to be similar: −11 W/m〈sup〉2〈/sup〉/AOD at all-sky and −24 W/m〈sup〉2〈/sup〉/AOD at clear-sky conditions. The present-day direct impact by all soot (BC) is globally averaged +0.55W/m〈sup〉2〈/sup〉 and at least half of it should be attributed to anthropogenic sources. Hereby any accuracy of anthropogenic impacts, not just for soot, suffers from the limited access to a pre-industrial reference. Anthropogenic uncertainty has a particular strong impact on aerosol indirect effects, which dominate the (TOA) forcing. Accounting for uncertainties in the anthropogenic definition, present-day aerosol forcing is estimated to stay within the −0.7 to −1.6 W/m〈sup〉2〈/sup〉 range, with a best estimate at −1 W/m〈sup〉2〈/sup〉. Calculations with model predicted temporal changes to anthropogenic AOD indicate that qualitatively the anthropogenic aerosol forcing has not changed much over the last decades and is not likely to increase over the next decades, despite strong regional shifts. These regional shifts explain most solar insolation (brightening or dimming) trends that have been observed by ground-based radiation data.〈/p〉

Description: 〈b〉Thermal structure of the mesopause region during the WADIS-2 rocket campaign〈/b〉〈br〉 Raimund Wörl, Boris Strelnikov, Timo P. Viehl, Josef Höffner, Pierre-Dominique Pautet, Michael J. Taylor, Yucheng Zhao, and Franz-Josef Lübken〈br〉 Atmos. Chem. Phys., 19, 77-88, https://doi.org/10.5194/acp-19-77-2019, 2019〈br〉 Simultaneous temperature measurements during the WADIS-2 rocket campaign are used to investigate the thermal structure of the mesopause region. Vertically and horizontally resolved in situ and remote measurements are in good agreement and show dominating long-term and large-scale waves with periods of 24 h and higher tidal harmonics. Only a few gravity waves with periods shorter than 6 h and small amplitudes are there.

Description: 〈b〉Verification of anthropogenic VOC emission inventory through ambient measurements and satellite retrievals〈/b〉〈br〉 Jing Li, Yufang Hao, Maimaiti Simayi, Yuqi Shi, Ziyan Xi, and Shaodong Xie〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1133,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 We established an emission inventory of anthropogenic VOCs in the Beijing–Tianjin–Hebei region of China. The developed emission inventory was evaluated through ambient measurements and satellite retrieval results. To obtain a more accurate emission inventory, we propose the investigation of the household coal consumption, the adjustment of EFs based on the latest pollution control policies, and the verification of the source profiles of OVOCs and halocarbons.

Description: 〈b〉Classification of aerosol population type and cloud condensation nuclei properties in a coastal California littoral environment using an unsupervised cluster model〈/b〉〈br〉 Samuel A. Atwood, Sonia M. Kreidenweis, Paul J. DeMott, Markus D. Petters, Gavin C. Cornwell, Andrew C. Martin, and Kathryn A. Moore〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1297,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 This paper presents measurements of aerosol particles at a coastal location. The particles were classified into distinct aerosol types using both microphysical measurements and meteorological information, allowing rapid changes between the aerosol types to be reliably identified. These particles can alter cloud and precipitation processes, and inclusion of the differences between types can improve atmospheric models and remote sensing retrievals in littoral zones.

Description: 〈b〉Development of a balloon-borne instrument for CO〈sub〉2〈/sub〉 vertical profile observations in the troposphere〈/b〉〈br〉 Mai Ouchi, Yutaka Matsumi, Tomoki Nakayama, Kensaku Shimizu, Takehiko Sawada, Toshinobu Machida, Hidekazu Matsueda, Yousuke Sawa, Isamu Morino, Osamu Uchino, Tomoaki Tanaka, and Ryoichi Imasu〈br〉 Atmos. Meas. Tech. Discuss., https//doi.org/10.5194/amt-2018-376,2019〈br〉 〈b〉Manuscript under review for AMT〈/b〉 (discussion: open, 0 comments)〈br〉 A novel, practical observation system for measuring tropospheric carbon dioxide (CO〈sub〉2〈/sub〉) concentrations carried by a small helium-filled balloon (CO〈sub〉2〈/sub〉 sonde), has been developed for the first time. The low-cost CO〈sub〉2〈/sub〉 sondes can potentially be used for frequently measurements of vertical profiles of CO〈sub〉2〈/sub〉 in any parts of the world providing as useful information to understand the global and regional carbon budgets by replenishing the present sparse observation coverage.

Description: 〈b〉Comparison of Different Sequential Assimilation Algorithms for Satellite-derived Leaf Area Index Using the Data Assimilation Research Testbed (lanai)〈/b〉〈br〉 Xiao-Lu Ling, Cong-Bin Fu, Zong-Liang Yang, and Wei-Dong Guo〈br〉 Geosci. Model Dev. Discuss., https//doi.org/10.5194/gmd-2018-232,2019〈br〉 〈b〉Manuscript under review for GMD〈/b〉 (discussion: open, 0 comments)〈br〉 Both observation and simulation can provide the temporal and spatial variation of vegetation characteristic, while they are not satisfactory for understanding the mechanism of the exchange between ecosystems and atmosphere. Data assimilation (DA) can combine the observation and models via mathematical statistical analysis.The results show that the Ensemble Adjust Kalman Filter (EAKF) is the optical algorithm. In addition, models perform better when the DA accept more proportion of observation.

Description: 〈b〉The importance of crystalline phases in ice nucleation by volcanic ash〈/b〉〈br〉 Elena C. Maters, Donald B. Dingwell, Corrado Cimarelli, Dirk Müller, Thomas F. Whale, and Benjamin J. Murray〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1326,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 This experimental study investigates the influence of volcanic ash chemical composition, crystallinity, and mineralogy on its ability to promote freezing of supercooled water. The results indicate that crystals in ash play a key role in this process, and suggest that alkali and plagioclase feldspars in ash may be highly ice-active. These findings contribute to improving understanding of the potential of ash emissions from different explosive eruptions to impact ice formation in the atmosphere.

Description: 〈b〉Large increases in N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉 together with a nucleation-modeparticle pool over the northwestern Pacific Ocean in the spring of 2014〈/b〉〈br〉 Juntao Wang, Yanjie Shen, Kai Li, Yang Gao, Huiwang Gao, and Xiaohong Yao〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1089,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 In this paper, we studied the spatiotemporal variability of N〈sub〉cn〈/sub〉 and particle number size distributions, as well as N〈sub〉ccn〈/sub〉 and CCN activities over the NWPO in the spring of 2014. Several new findings have been revealed and discussed, e.g., large increases of N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉, against historical data and small contributions of dust and BB aerosols to N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉 on the monthly time scale, etc. This work may help understand the influence of the current outflow from the Asian continent on the climate.

Description: 〈b〉Analyses of temperature and precipitation in the Indian Jammu and Kashmir region for the 1980–2016 period: implications for remote influence and extreme events〈/b〉〈br〉 Sumira Nazir Zaz, Shakil Ahmad Romshoo, Ramkumar Thokuluwa Krishnamoorthy, and Yesubabu Viswanadhapalli〈br〉 Atmos. Chem. Phys., 19, 15-37, https://doi.org/10.5194/acp-19-15-2019, 2019〈br〉 This paper is of first of its kind for the Jammu and Kashmir (western Himalayas) region, India. It shows the clear relation between the upper tropospheric Rossby wave activity (potential vorticity at the 350 K potential temperature and 200 mb level surface pressure) and the surface weather parameters (e.g., precipitation) over a period of 3 decades during 1980–2016. This indicates that the climatic weather pattern over the Kashmir region is influenced mostly by global climate change processes.

Description: 〈b〉Secondary organic aerosol formation from smoldering and flaming combustion of biomass: a box model parametrization based on volatility basis set〈/b〉〈br〉 Giulia Stefenelli, Jianhui Jiang, Amelie Bertrand, Emily A. Bruns, Simone M. Pieber, Urs Baltensperger, Nicolas Marchand, Sebnem Aksoyoglu, André S. H. Prévôt, Jay G. Slowik, and Imad El Haddad〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1308,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 Box model simulations, based on the volatility basis set (VBS) approach, of smog chamber wood combustion experiments conducted at different temperatures (−10 °C, 2 °C, 15 °C), emission loads, combustion conditions (flaming and smoldering) and residential stoves fabricated in the last two decades. Novel parameterization methods based on a genetic algorithm approach allowed estimation of precursor class contributions to SOA and evaluation of the effect of emission variability on SOA yield predictions.

Description: 〈b〉Impacts of Horizontal Resolution and Air–Sea Flux Parameterization on the Intensity and Structure of simulated Typhoon Haiyan (2013)〈/b〉〈br〉 Mien-Tze Kueh, Wen-Mei Chen, Yang-Fan Sheng, Simon C. Lin, Tso-Ren Wu, Eric Yen, Yu-Lin Tsai, and Chuan-Yao Lin〈br〉 Nat. Hazards Earth Syst. Sci. Discuss., https//doi.org/10.5194/nhess-2018-333,2019〈br〉 〈b〉Manuscript under review for NHESS〈/b〉 (discussion: open, 0 comments)〈br〉 In this study, we show that both the model horizontal resolution and air-sea flux parameterization can exert large influence on tropical cyclone intensity simulation but with different impacts on wind structures. We highlight the intensification and contraction of TC eyewall in response to the reduction of grid spacing. We also suggest that a well-developed eyewall is more conducive to the positive effect of flux formulas on TC development.

Description: 〈b〉Primary and secondary sources of ambient formaldehyde in the Yangtze River Delta based on OMPS observation〈/b〉〈br〉 Wenjing Su, Cheng Liu, Qihou Hu, Shaohua Zhao, Youwen Sun, Wei Wang, Yizhi Zhu, Jianguo Liu, and Jhoon Kim〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1192,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 For a better understanding of the pollution and atmospheric chemistry of HCHO, we evaluated its primary and secondary sources in the Yangtze River Delta based on HCHO column density from Ozone Mapping and Profiler Suite (OMPS) and combining with surface in-situ measurements. We found that secondary formation contributed most to ambient HCHO for a long time, but primary emission could be dominant in the winter. The usability of total HCHO as the proxy of VOCs reactivity depends on time scale.

Description: 〈b〉Local and regional contributions to fine particulate matter in the 18 cities of Sichuan Basin, southwestern China〈/b〉〈br〉 Xue Qiao, Hao Guo, Ya Tang, Pengfei Wang, Wenye Deng, Xing Zhao, Jianlin Hu, Qi Ying, and Hongliang Zhang〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1214,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 A source-oriented version of the CMAQ model was used to quantify contributions from nine regions to PM〈sub〉2.5〈/sub〉 and its components in the 18 cities within the Sichuan Basin. Non-local emissions contribute to 39~66 % and 25~52 % to the citywide average PM〈sub〉2.5〈/sub〉 concentrations of 45~126 and 14~31 μg m〈sup〉−3〈/sup〉 in the winter and summer, respectively. This study demonstrates the importance of joint emission control efforts among cities within the SCB and neighbor regions to the east.

Description: 〈b〉MAX-DOAS measurements of tropospheric NO〈sub〉2〈/sub〉 and HCHO in Nanjing and the comparison to OMI observations〈/b〉〈br〉 Ka Lok Chan, Zhuoru Wang, Aijun Ding, Klaus-Peter Heue, Yicheng Shen, Jing Wang, Feng Zhang, Nan Hao, and Mark Wenig〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1266,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 The manuscript presents long term observations of atmospheric nitrogen dioxide (NO〈sub〉2〈/sub〉) and formaldehyde (HCHO) in Nanjing using a MAX-DOAS instrument. The measurements were performed from April 2013 to February 2017. The MAX-DOAS measurements of NO〈sub〉2〈/sub〉 and HCHO are used to validate OMI satellite observations and to investigate the influences of region transport of air pollutants on the air quality in Nanjing.

Description: 〈b〉Detection of Outflow of Formaldehyde and Glyoxal from the African continent to the Atlantic Ocean with a MAX-DOAS Instrument〈/b〉〈br〉 Lisa K. Behrens, Andreas Hilboll, Andreas Richter, Enno Peters, Leonardo M. A. Alvarado, Anna B. Kalisz Hedegaard, Folkard Wittrock, John P. Burrows, and Mihalis Vrekoussis〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1286,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 MAX-DOAS measurements were conducted on the research vessel Maria S. Merian during a cruise from the Azores to South Africa in October 2016. The measurements indicate enhanced levels of HCHO and CHOCHO over the remote Atlantic Ocean, which is unexpected due to their short lifetime. Precursors of these gases or gas-aerosol combinations might be transported. Model simulations indicate potential source regions over the African continent, probably related to biomass burning or biogenic emissions.

Description: 〈b〉Develop a coupled agent-based modeling approach for uncertain water management decisions〈/b〉〈br〉 Jin-Young Hyun, Shih-Yu Huang, Y. C. Ethan Yang, Vincent Tidwell, and Jordan Macknick〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-555,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 This study applies a 〈q〉two-way〈/q〉 coupled agent-based model (ABM) with a River-Reservoir management model (RiverWare) to address the interaction between human and natural systems using Bayesian Inference (BI) mapping joined with Cost-Loss (CL). The calibration results show that this methodology can capture the historical pattern of both human activities and natural dynamics and outperforms those without using BI and CL.

Description: 〈b〉Potentials and pitfalls of permafrost active layer monitoring using the HVSR method: a case study in Svalbard〈/b〉〈br〉 Andreas Köhler and Christian Weidle〈br〉 Earth Surf. Dynam., 7, 1-16, https://doi.org/10.5194/esurf-7-1-2019, 2019〈br〉 The uppermost part of permanently frozen ground can thaw during summer and refreeze during winter. We use a method based on naturally generated seismic waves to continuously monitor these changes close to the research settlement of Ny-Ålesund in Svalbard between April and August 2016. Our results reveal some potential pitfalls when interpreting temporal variations in the data. However, we show that a careful data analysis makes this method a very useful tool for long-term permafrost monitoring.

Description: 〈b〉Long-profile evolution of transport-limited gravel-bed rivers〈/b〉〈br〉 Andrew D. Wickert and Taylor F. Schildgen〈br〉 Earth Surf. Dynam., 7, 17-43, https://doi.org/10.5194/esurf-7-17-2019, 2019〈br〉 Rivers can raise or lower their beds by depositing or eroding sediments. We combine equations for flow, channel/valley geometry, and gravel transport to learn how climate and tectonics shape down-valley profiles of river-bed elevation. Rivers steepen when they receive more sediment (relative to water) and become straighter with tectonic uplift. Weathering and breakdown of gravel is needed to produce gradually widening river channels with concave-up profiles that are often observed in the field.

Description: 〈b〉Definition differences and internal variability affect the simulated Arctic sea ice melt season〈/b〉〈br〉 Abigail Smith and Alexandra Jahn〈br〉 The Cryosphere, 13, 1-20, https://doi.org/10.5194/tc-13-1-2019, 2019〈br〉 Here we assessed how natural climate variations and different definitions impact the diagnosed and projected Arctic sea ice melt season length using model simulations. Irrespective of the definition or natural variability, the sea ice melt season is projected to lengthen, potentially by as much as 4–5 months by 2100 under the business as usual scenario. We also find that different definitions have a bigger impact on melt onset, while natural variations have a bigger impact on freeze onset.

Description: 〈b〉Three decades of simulated global terrestrial carbon fluxes from a data assimilation system confronted to different periods of observations〈/b〉〈br〉 Karel Castro-Morales, Gregor Schürmann, Christoph Köstler, Christian Rödenbeck, Martin Heimann, and Sönke Zaehle〈br〉 Biogeosciences Discuss., doi:10.5194/bg-2018-517,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 To obtain nearly thirty years of global terrestrial carbon fluxes, we simultaneously incorporated in a land surface model three different time periods of two observational data sets: absorbed photosynthetic active radiation and atmospheric CO〈sub〉2〈/sub〉 concentrations. One decade of data is enough to improve the modeled long-term trends and seasonal amplitudes of the assimilated variables, particularly in boreal regions. This model has the potential to provide short-term predictions of land carbon fluxes.

Description: 〈b〉Ostracods as ecological and isotopic indicators of lake water salinity changes: The Lake Van example〈/b〉〈br〉 Jeremy McCormack, Finn Viehberg, Derya Akdemir, Adrian Immenhauser, and Ola Kwiecien〈br〉 Biogeosciences Discuss., doi:10.5194/bg-2018-476,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 We juxtapose changes in ostracod taxonomy, morphology (noding) and oxygen (δ〈sup〉18〈/sup〉O) and carbon (δ〈sup〉13〈/sup〉C) isotopic composition for the last 150 kyr with independent low-resolution salinity proxies. We demonstrate that for Lake Van, salinity is the most important factor influencing the composition of the ostracod assemblage and the formation of nodes on the valves of limnocytherinae species. Ostracod δ〈sup〉18〈/sup〉O shows a higher sensibility towards climatic and hydrological variations than the bulk isotopy.

Description: 〈b〉Organic signatures in Pleistocene cherts from Lake Magadi (Kenya), analogs for early Earth hydrothermal deposits〈/b〉〈br〉 Manuel Reinhardt, Walter Goetz, Jan-Peter Duda, Christine Heim, Joachim Reitner, and Volker Thiel〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-513,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 Organic matter in Archean hydrothermal cherts may contain molecular traces of early life. Alteration processes during and after deposition, however, may have obliterated potential biosignatures. Our results from modern analog samples (Pleistocene cherts from Lake Magadi, Kenya) show that biomolecules can survive early hydrothermal destruction in the macromolecular fraction of the organic matter. A conservation of molecular biosignatures in Archean hydrothermal cherts therefore seems possible.

Description: 〈b〉Detection of Outflow of Formaldehyde and Glyoxal from the African continent to the Atlantic Ocean with a MAX-DOAS Instrument〈/b〉〈br〉 Lisa K. Behrens, Andreas Hilboll, Andreas Richter, Enno Peters, Leonardo M. A. Alvarado, Anna B. Kalisz Hedegaard, Folkard Wittrock, John P. Burrows, and Mihalis Vrekoussis〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1286,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 MAX-DOAS measurements were conducted on the research vessel Maria S. Merian during a cruise from the Azores to South Africa in October 2016. The measurements indicate enhanced levels of HCHO and CHOCHO over the remote Atlantic Ocean, which is unexpected due to their short lifetime. Precursors of these gases or gas-aerosol combinations might be transported. Model simulations indicate potential source regions over the African continent, probably related to biomass burning or biogenic emissions.

Description: 〈b〉Primary and secondary sources of ambient formaldehyde in the Yangtze River Delta based on OMPS observation〈/b〉〈br〉 Wenjing Su, Cheng Liu, Qihou Hu, Shaohua Zhao, Youwen Sun, Wei Wang, Yizhi Zhu, Jianguo Liu, and Jhoon Kim〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1192,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 For a better understanding of the pollution and atmospheric chemistry of HCHO, we evaluated its primary and secondary sources in the Yangtze River Delta based on HCHO column density from Ozone Mapping and Profiler Suite (OMPS) and combining with surface in-situ measurements. We found that secondary formation contributed most to ambient HCHO for a long time, but primary emission could be dominant in the winter. The usability of total HCHO as the proxy of VOCs reactivity depends on time scale.

Description: 〈b〉Greenland Ice Sheet late-season melt: Investigating multi-scale drivers of K-transect events〈/b〉〈br〉 Thomas J. Ballinger, Thomas L. Mote, Kyle Mattingly, Angela C. Bliss, Edward Hanna, Dirk van As, Melissa Prieto, Saeideh Gharehchahi, Xavier Fettweis, Brice Noël, Paul C. J. P. Smeets, Mads H. Ribergaard, and John Cappelen〈br〉 The Cryosphere Discuss., https//doi.org/10.5194/tc-2018-285,2019〈br〉 〈b〉Manuscript under review for TC〈/b〉 (discussion: open, 0 comments)〈br〉 Studies have questioned links between Arctic marginal sea open water duration and Greenland Ice Sheet (GrIS) surface balance changes, namely melt events. Through analyses involving observations and climate models, we show that late summer through autumn “unseasonal” melt events are primarily driven by the northward movement of warm, moist air masses across the western ice sheet edge, while near-surface, off-ice winds block heat transfer off nearby Baffin Bay.

Description: 〈b〉Snow cover variability across glaciers in Nordenskiöldland (Svalbard) from point measurements in 2014–2016〈/b〉〈br〉 Marco Möller and Rebecca Möller〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-158,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 0 comments)〈br〉 Snow water equivalents (SWE) are useful for calibration and validation of different snow cover-related modeling efforts (e.g. glacier mass balance or snow drift). They form a measure of snow accumulation better than snow depth as they are independent of density differences and thus comparable. A dataset is presented which provides point data of snow depth and density measurements as well as of calculated SWE from 109 locations on glaciers in the central region of the Arctic archipelago Svalbard.

Description: 〈b〉Fifty years of recorded hillslope runoff on seasonally-frozen ground: The Swift Current, Saskatchewan, Canada dataset〈/b〉〈br〉 Anna E. Coles, Jeffrey J. McDonnell, and Brian G. McConkey〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-126,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 1 comment)〈br〉 〈p〉Long records of hillslope runoff and nutrients are rare – on seasonally frozen ground they are almost non-existent. The Swift Current hillslopes at the Swift Current Research and Development Centre on the Canadian Prairies provide such a long-term hydrological record. Runoff, runoff nutrient concentration, snowpack, soil moisture, and soil nutrient concentration were monitored on the three 5 ha hillslopes over a 50-year period (1962–2011). Runoff from the hillslopes was generated episodically during snowmelt and occasional rainfall events. Edge-of-field runoff was measured with a 0.61 m H-flume. Daily runoff nutrient concentration data are available for nitrate (March 1971–April 2011), ammonium (February 1996–April 2011), and orthophosphate (March–April 1971; June 1991–April 2011). Snowpack data (snowpack depth, density and water equivalent) were determined via manual snow surveys carried out several times each winter, between January and March, between 1965 and 2011 Gravimetric soil moisture content was measured in October and April each year between 1971 and 2011 at five depth intervals (0–15, 15–30, 30–60, 60–90, and 90–120 cm) at nine points on each hillslope. We summarize these hillslope data in two publically-available repositories: 1) 1962–2011 data on runoff, runoff nutrients, snowpack, soil moisture, soil nutrients, and crop and tillage practices at 〈a href="https://doi.org/10.23684/hhn5-rz52" target="_blank"〉https://doi.org/10.23684/hhn5-rz52〈/a〉; and 2) digital elevation data at 〈a href="https://doi.org/10.20383/101.011" target="_blank"〉https://doi.org/10.20383/101.011〈/a〉. Complete climate data recorded at a Environment and Climate Change Canada meteorological station located 390 m from the three hillslopes are publically-available at 〈a href="http://climate.weather.gc.ca/" target="_blank"〉http://climate.weather.gc.ca/〈/a〉.〈/p〉

Description: 〈b〉GemPy 1.0: open-source stochastic geological modeling and inversion〈/b〉〈br〉 Miguel de la Varga, Alexander Schaaf, and Florian Wellmann〈br〉 Geosci. Model Dev., 12, 1-32, https://doi.org/10.5194/gmd-12-1-2019, 2019〈br〉 〈i〉GemPy〈/i〉 is an open-source Python-based 3-D structural geological modeling software, which allows the implicit (i.e. automatic) creation of complex geological models from interface and orientation data. 〈i〉GemPy〈/i〉 is implemented in the programming language Python, making use of a highly efficient underlying library, Theano, for efficient code generation that performs automatic differentiation. This enables the link to probabilistic machine-learning and Bayesian inference frameworks.

Description: 〈b〉Comparison of the GRUAN data products for Meisei RS-11G and Vaisala RS92-SGP radiosondes at Tateno (36.06° N, 140.13° E), Japan〈/b〉〈br〉 Eriko Kobayashi, Shunsuke Hoshino, Masami Iwabuchi, Takuji Sugidachi, Kensaku Shimizu, and Masatomo Fujiwara〈br〉 Atmos. Meas. Tech. Discuss., https//doi.org/10.5194/amt-2018-416,2019〈br〉 〈b〉Manuscript under review for AMT〈/b〉 (discussion: open, 0 comments)〈br〉 The authors carried out dual flights of RS-11G and RS92-SGP radiosondes and investigated the differences in the performance of the radiosondes to help characterize GRUAN data products. A novel aspect of GRUAN data products is that vertically resolved uncertainty estimates and metadata are provided for each sounding and comparison of GRUAN data products is important in securing the temporal homogeneity of climate data records.

Description: 〈b〉WHU-SGCC: A novel approach for blending daily satellite (CHIRP) and precipitation observations over Jinsha River Basin〈/b〉〈br〉 Gaoyun Shen, Nengcheng Chen, Wei Wang, and Zeqiang Chen〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-150,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 0 comments)〈br〉 The development of effective methods for high-accuracy precipitation estimates over complex terrain and on a daily scale is important for mountainous hydrological applications. This study offers a novel approach called WHU-SGCC by blending rain gauged and satellite data, to estimate daily precipitation at 0.05° resolution over Jinsha River Basin, the complicated mountainous terrain with sparse rain gauge data, considering the spatial correlation and the historical precipitation characteristics.

Description: 〈b〉GemPy 1.0: open-source stochastic geological modeling and inversion〈/b〉〈br〉 Miguel de la Varga, Alexander Schaaf, and Florian Wellmann〈br〉 Geosci. Model Dev., 12, 1-32, https://doi.org/10.5194/gmd-12-1-2019, 2019〈br〉 〈i〉GemPy〈/i〉 is an open-source Python-based 3-D structural geological modeling software, which allows the implicit (i.e. automatic) creation of complex geological models from interface and orientation data. 〈i〉GemPy〈/i〉 is implemented in the programming language Python, making use of a highly efficient underlying library, Theano, for efficient code generation that performs automatic differentiation. This enables the link to probabilistic machine-learning and Bayesian inference frameworks.

Description: 〈b〉Merits of novel high-resolution estimates and existing long-term estimates of humidity and incident radiation in a complex domain〈/b〉〈br〉 Helene Birkelund Erlandsen, Lena Merete Tallaksen, and Jørn Kristiansen〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-121,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 0 comments)〈br〉 Robust estimates of runoff, snow, and evaporation rely on high quality estimates of incoming solar and thermal radiation at the surface, and near surface humidity. Taking advantage of the physical soundness of a numerical weather reanalyisis and the preciseness and spatial resolution of a national gridded temperature data set, new estimates of these variables are presented for Norway. Further, existing data sets and observations are compared, emphasising daily correlation, trends, and gradients.

Description: 〈b〉Fifty years of recorded hillslope runoff on seasonally-frozen ground: The Swift Current, Saskatchewan, Canada dataset〈/b〉〈br〉 Anna E. Coles, Jeffrey J. McDonnell, and Brian G. McConkey〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-126,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 1 comment)〈br〉 〈p〉Long records of hillslope runoff and nutrients are rare – on seasonally frozen ground they are almost non-existent. The Swift Current hillslopes at the Swift Current Research and Development Centre on the Canadian Prairies provide such a long-term hydrological record. Runoff, runoff nutrient concentration, snowpack, soil moisture, and soil nutrient concentration were monitored on the three 5 ha hillslopes over a 50-year period (1962–2011). Runoff from the hillslopes was generated episodically during snowmelt and occasional rainfall events. Edge-of-field runoff was measured with a 0.61 m H-flume. Daily runoff nutrient concentration data are available for nitrate (March 1971–April 2011), ammonium (February 1996–April 2011), and orthophosphate (March–April 1971; June 1991–April 2011). Snowpack data (snowpack depth, density and water equivalent) were determined via manual snow surveys carried out several times each winter, between January and March, between 1965 and 2011 Gravimetric soil moisture content was measured in October and April each year between 1971 and 2011 at five depth intervals (0–15, 15–30, 30–60, 60–90, and 90–120 cm) at nine points on each hillslope. We summarize these hillslope data in two publically-available repositories: 1) 1962–2011 data on runoff, runoff nutrients, snowpack, soil moisture, soil nutrients, and crop and tillage practices at 〈a href="https://doi.org/10.23684/hhn5-rz52" target="_blank"〉https://doi.org/10.23684/hhn5-rz52〈/a〉; and 2) digital elevation data at 〈a href="https://doi.org/10.20383/101.011" target="_blank"〉https://doi.org/10.20383/101.011〈/a〉. Complete climate data recorded at a Environment and Climate Change Canada meteorological station located 390 m from the three hillslopes are publically-available at 〈a href="http://climate.weather.gc.ca/" target="_blank"〉http://climate.weather.gc.ca/〈/a〉.〈/p〉

Description: 〈b〉Calibration and analysis of the uncertainty in downscaling global land use and land cover projections from GCAM〈/b〉〈br〉 Min Chen, Chris R. Vernon, Maoyi Huang, Katherine V. Calvin, and Ian P. Kraucunas〈br〉 Geosci. Model Dev. Discuss., https//doi.org/10.5194/gmd-2018-248,2019〈br〉 〈b〉Manuscript under review for GMD〈/b〉 (discussion: open, 0 comments)〈br〉 Demeter is a community spatial downscaling model that disaggregates land use and land cover changes projected by integrated human-Earth system models. However, Demeter has not been intensively calibrated, and we still lack a good knowledge about its sensitivity to key parameters and the parameter uncertainties. This paper aims to solve this problem.

Description: 〈b〉WHU-SGCC: A novel approach for blending daily satellite (CHIRP) and precipitation observations over Jinsha River Basin〈/b〉〈br〉 Gaoyun Shen, Nengcheng Chen, Wei Wang, and Zeqiang Chen〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-150,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 0 comments)〈br〉 The development of effective methods for high-accuracy precipitation estimates over complex terrain and on a daily scale is important for mountainous hydrological applications. This study offers a novel approach called WHU-SGCC by blending rain gauged and satellite data, to estimate daily precipitation at 0.05° resolution over Jinsha River Basin, the complicated mountainous terrain with sparse rain gauge data, considering the spatial correlation and the historical precipitation characteristics.

Description: 〈b〉On the fine vertical structure of the low troposphere over the coastal margins of East Antarctica〈/b〉〈br〉 Étienne Vignon, Olivier Traullé, and Alexis Berne〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1197,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 The sea-level rise associated to climate change will depend on how much the Antarctic ice sheet gain – viz precipitation – or loose mass. The precipitation simulation by climate models used for projections depends on the representation of the atmospheric circulation over and around Antarctica. Using daily measurements from balloon soundings at nine stations, this study characterizes the structure of the atmosphere over the Antarctic coast and assess its representation in a climate model.

Description: 〈b〉Insights into the O : C dependent mechanisms controlling the evaporation of α-pinene secondary organic aerosol particles〈/b〉〈br〉 Angela Buchholz, Andrew T. Lambe, Arttu Ylisirniö, Zijun Li, Olli-Pekka Tikkanen, Celia Faiola, Eetu Kari, Liqing Hao, Olli Luoma, Wei Huang, Claudia Mohr, Douglas R. Worsnop, Sergey A. Nizkorodov, Taina Yli-Juuti, Siegfried Schobesberger, and Annele Virtanen〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1305,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 We studied the evaporation of α-pinene secondary organic aerosol particles in clean air to derive their volatility from the observed size changes. We found that the particles became more resilient to evaporation with increased oxidative age, possibly increasing their lifetime in the atmosphere. Also, increased relative humidity increased the particle evaporation. Mass spectrometry measurements of the particles at different stages of evaporation revealed some water-induced composition changes.

Description: 〈b〉Insights into the O : C dependent mechanisms controlling the evaporation of α-pinene secondary organic aerosol particles〈/b〉〈br〉 Angela Buchholz, Andrew T. Lambe, Arttu Ylisirniö, Zijun Li, Olli-Pekka Tikkanen, Celia Faiola, Eetu Kari, Liqing Hao, Olli Luoma, Wei Huang, Claudia Mohr, Douglas R. Worsnop, Sergey A. Nizkorodov, Taina Yli-Juuti, Siegfried Schobesberger, and Annele Virtanen〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1305,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 We studied the evaporation of α-pinene secondary organic aerosol particles in clean air to derive their volatility from the observed size changes. We found that the particles became more resilient to evaporation with increased oxidative age, possibly increasing their lifetime in the atmosphere. Also, increased relative humidity increased the particle evaporation. Mass spectrometry measurements of the particles at different stages of evaporation revealed some water-induced composition changes.

Description: 〈b〉Quantitative precipitation estimation with weather radar using a data- and information-based approach〈/b〉〈br〉 Malte Neuper and Uwe Ehret〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-606,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 In our study we apply a data-driven approach to estimate precipitation quantitatively with the use of weather radar data. The method is based on Information theory concepts. It uses predictive relations expressed by empirical discrete probability distributions directly derived from data rather than the standard deterministic functions.

Description: 〈b〉Develop a coupled agent-based modeling approach for uncertain water management decisions〈/b〉〈br〉 Jin-Young Hyun, Shih-Yu Huang, Y. C. Ethan Yang, Vincent Tidwell, and Jordan Macknick〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-555,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 This study applies a 〈q〉two-way〈/q〉 coupled agent-based model (ABM) with a River-Reservoir management model (RiverWare) to address the interaction between human and natural systems using Bayesian Inference (BI) mapping joined with Cost-Loss (CL). The calibration results show that this methodology can capture the historical pattern of both human activities and natural dynamics and outperforms those without using BI and CL.

Description: 〈b〉Stochastic modeling of flow and conservative transport in three-dimensional discrete fracture networks〈/b〉〈br〉 I-Hsien Lee, Chuen-Fa Ni, Fang-Pang Lin, Chi-Ping Lin, and Chien-Chung Ke〈br〉 Hydrol. Earth Syst. Sci., 23, 19-34, https://doi.org/10.5194/hess-23-19-2019, 2019〈br〉 Few studies focused on the direct solution of the ADE for 3-D DFNs. The study is the first to solve the ADE and focus on assessing the velocity uncertainty in 3-D DFNs. The velocity uncertainty shows a limited range of influence close to the mean diameter of a fracture. The information is useful for engineering designs at sites with fractured rocks. We quantified that the tracer test in wells might lead to the overestimation of mean concentration and induce high uncertainty in fractured media.

Description: 〈b〉Time-variability of the fraction of young water in a small headwater catchment〈/b〉〈br〉 Michael P. Stockinger, Heye R. Bogena, Andreas Lücke, Christine Stumpp, and Harry Vereecken〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-604,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 Precipitation moves through the soil to become streamwater. The fraction of precipitation that became streamwater after 3 months (Fyw) can be calculated with the stable isotopes of water. Previously this was done for all the isotope data available, e.g., for several years. We used one year of data to calculate Fyw and moved this calculation time window over the time series. Results highlight that Fyw varies in time. Comparison studies of different regions should take this into account.

Description: 〈b〉The significant contribution of HONO to secondary pollutants during a severe winter pollution event in southern China〈/b〉〈br〉 Xiao Fu, Tao Wang, Li Zhang, Qinyi Li, Zhe Wang, Men Xia, Hui Yun, Weihao Wang, Chuan Yu, Dingli Yue, Yan Zhou, Junyun Zheng, and Rui Han〈br〉 Atmos. Chem. Phys., 19, 1-14, https://doi.org/10.5194/acp-19-1-2019, 2019〈br〉 This study has identified the major contributors to the observed high HONO levels during a severe winter pollution episode and highlighted the importance of HONO chemistry in the combined photochemical and haze pollution in a subtropical region. It also highlighted the critical need to include and update HONO sources in regional air quality models in order to predict ozone and other secondary pollutants better during heavy pollution events in southern China and similar regions.

Description: 〈b〉A reconstruction of warm water inflow to Upernavik Isstrøm since AD 1925 and its relation to glacier retreat〈/b〉〈br〉 Flor Vermassen, Nanna Andreasen, David J. Wangner, Nicolas Thibault, Marit-Solveig Seidenkrantz, Rebecca Jackson, Sabine M. Schmidt, Kurt H. Kjær, and Camilla S. Andresen〈br〉 Clim. Past Discuss., https//doi.org/10.5194/cp-2018-174,2019〈br〉 〈b〉Manuscript under review for CP〈/b〉 (discussion: open, 0 comments)〈br〉 This study shows that warming of ocean waters is related to the retreat of Upernavik Isstrøm, a glacier in Northwest Greenland. We show that in the 1930s and after 2000 the waters in the fjord warmed and the glacier retreated. We found this by investigating microfossils from sediments in Upernavik Fjord; different species occur in response to warmer waters.

Description: 〈b〉The CMEMS GlobColour 〈i〉Chlorophyll-a〈/i〉 Product Based on Satellite Observation〈/b〉〈br〉 Philippe Garnesson, Antoine Mangin, Odile Fanton d'Andon, Julien Demaria, and Marine Bretagnon〈br〉 Ocean Sci. Discuss., https//doi.org/10.5194/os-2018-155,2019〈br〉 〈b〉Manuscript under review for OS〈/b〉 (discussion: open, 0 comments)〈br〉 〈p〉This work concerns the chlorophyll products based on Satellite Observation and disseminated in the frame of the Copernicus Marine Environmental Monitoring Service (CMEMS).〈/p〉 〈p〉 This work highlights the main advantages provided by the Copernicus Globcolour processor which is used to serve the CMEMS with a long time series from 1997 to present with level 3 & 4 products at Global level (4 km of spatial resolution) and for the Atlantic level 4 product (1 km).〈/p〉 〈p〉 It discusses the different ways to merge data coming from different sensors and it is shown that the GlobColour processor approach provide a better flexibility. At present, it is the only one CMEMS processor able to ingest the OLCI-S3A in the merged product (OLCI-S3A data are ingested in the operational CMEMS products since the April 2018 release).〈/p〉 〈p〉 Behind the merging, the flagging strategy to go from level 2 provided by spatial agencies to the level 3 CMEMS products is also discussed. A better spatial coverage is demonstrated, including the coastal area which is of particular interest for many users involved in the EU Water Framework and Marine Strategy Framework Directive.〈/p〉

Description: 〈b〉Induced seismicity in geologic carbon storage〈/b〉〈br〉 Víctor Vilarrasa, Jesus Carrera, Sebastià Olivella, Jonny Rutqvist, and Lyesse Laloui〈br〉 Solid Earth Discuss., https//doi.org/10.5194/se-2018-129,2019〈br〉 〈b〉Manuscript under review for SE〈/b〉 (discussion: open, 0 comments)〈br〉 To meet the goal of the Paris Agreement to limit temperature increase below 2 ºC, Geologic Carbon Storage (GCS) will be necessary at the gigatone scale. But to successfully deploy GCS, seismicity induced by CO〈sub〉2〈/sub〉 injection should be controlled and maintained below a threshold that does not generate nuisances to the population. We conclude that felt induced seismicity can be minimized provided that a proper site characterization, monitoring and pressure management are performed.

Description: 〈b〉Improved aerosol correction for OMI tropospheric NO〈sub〉2〈/sub〉 retrieval over East Asia: constraint from CALIOP aerosol vertical profile〈/b〉〈br〉 Mengyao Liu, Jintai Lin, K. Folkert Boersma, Gaia Pinardi, Yang Wang, Julien Chimot, Thomas Wagner, Pinhua Xie, Henk Eskes, Michel Van Roozendael, François Hendrick, Pucai Wang, Ting Wang, Yingying Yan, Lulu Chen, and Ruijing Ni〈br〉 Atmos. Meas. Tech., 12, 1-21, https://doi.org/10.5194/amt-12-1-2019, 2019〈br〉 China has become the world’s largest emitter of NO〈sub〉〈i〉x〈/i〉〈/sub〉, which mainly comes from vehicle exhaust, power plants, etc. However, there are no official ground-based measurements before 2013, so satellites have been widely used to monitor and analyze NO〈sub〉〈i〉x〈/i〉〈/sub〉 pollution here. Aerosol is the key factor influencing the accuracy of the satellite NO〈sub〉〈i〉x〈/i〉〈/sub〉 product. Our study provides a more accurate way to account for aerosol's influence compared to current widely used products.

Description: 〈b〉Comparison of the GRUAN data products for Meisei RS-11G and Vaisala RS92-SGP radiosondes at Tateno (36.06° N, 140.13° E), Japan〈/b〉〈br〉 Eriko Kobayashi, Shunsuke Hoshino, Masami Iwabuchi, Takuji Sugidachi, Kensaku Shimizu, and Masatomo Fujiwara〈br〉 Atmos. Meas. Tech. Discuss., https//doi.org/10.5194/amt-2018-416,2019〈br〉 〈b〉Manuscript under review for AMT〈/b〉 (discussion: open, 0 comments)〈br〉 The authors carried out dual flights of RS-11G and RS92-SGP radiosondes and investigated the differences in the performance of the radiosondes to help characterize GRUAN data products. A novel aspect of GRUAN data products is that vertically resolved uncertainty estimates and metadata are provided for each sounding and comparison of GRUAN data products is important in securing the temporal homogeneity of climate data records.

Description: 〈b〉Snow cover variability across glaciers in Nordenskiöldland (Svalbard) from point measurements in 2014–2016〈/b〉〈br〉 Marco Möller and Rebecca Möller〈br〉 Earth Syst. Sci. Data Discuss., https//doi.org/10.5194/essd-2018-158,2019〈br〉 〈b〉Manuscript under review for ESSD〈/b〉 (discussion: open, 0 comments)〈br〉 Snow water equivalents (SWE) are useful for calibration and validation of different snow cover-related modeling efforts (e.g. glacier mass balance or snow drift). They form a measure of snow accumulation better than snow depth as they are independent of density differences and thus comparable. A dataset is presented which provides point data of snow depth and density measurements as well as of calculated SWE from 109 locations on glaciers in the central region of the Arctic archipelago Svalbard.

Description: 〈b〉Analyses of temperature and precipitation in the Indian Jammu and Kashmir region for the 1980–2016 period: implications for remote influence and extreme events〈/b〉〈br〉 Sumira Nazir Zaz, Shakil Ahmad Romshoo, Ramkumar Thokuluwa Krishnamoorthy, and Yesubabu Viswanadhapalli〈br〉 Atmos. Chem. Phys., 19, 15-37, https://doi.org/10.5194/acp-19-15-2019, 2019〈br〉 This paper is of first of its kind for the Jammu and Kashmir (western Himalayas) region, India. It shows the clear relation between the upper tropospheric Rossby wave activity (potential vorticity at the 350 K potential temperature and 200 mb level surface pressure) and the surface weather parameters (e.g., precipitation) over a period of 3 decades during 1980–2016. This indicates that the climatic weather pattern over the Kashmir region is influenced mostly by global climate change processes.