ALBERT

Description: This investigation analyzes the effect of vortex wakes on the Lagrangian displacement of particles induced by the passage of an obstacle in a two-dimensional incompressible and inviscid fluid. In addition to the trajectories of individual particles, we also study their drift and the corresponding total drift areas in the Föppl and Kirchhoff potential flow models. Our findings, which are obtained numerically and in some regimes are also supported by asymptotic analysis, are compared to the wakeless potential flow which serves as a reference. We show that in the presence of the Föppl vortex wake, some of the particles follow more complicated trajectories featuring a second loop. The appearance of an additional stagnation point in the Föppl flow is identified as a source of this effect. It is also demonstrated that, while the total drift area increases with the size of the wake for large vortex strengths, it is actually decreased for small circulation values. On the other hand, the Kirchhoff flow model is shown to have an unbounded total drift area. By providing a systematic account of the wake effects on the drift, the results of this study will allow for more accurate modeling of hydrodynamic stirring.

Description: The slow motion of a circular cylinder in a plane Poiseuille flow in a microchannel is analyzed for a wide range of cylinder radii and positions across the channel. The cylinder translates parallel to the channel walls and rotates about its axis. The Stokes approximation is used and the problem is solved analytically using the Papkovich-Fadle eigenfunction expansion and the least-squares method. The stream function and the pressure distribution of the flow field are obtained as results. The force and moment exerted on the cylinder, and the pressure change far from the cylinder, are calculated and shown as functions of the size and location of the cylinder. The results confirm some reciprocal relations exactly. In particular, the translational and rotational velocities of the drifting cylinder in the existing Poiseuille flow are determined. The induced pressure change, when the cylinder drifts in the Poiseuille flow, is also calculated. Some typical streamline patterns, depending on the size and location of the cylinder, are shown and discussed. When the cylinder translates and/or rotates in the channel blocked at infinity, a series of Moffatt eddies appears far from the cylinder in the channel, as expected.

Description: Interactions between capillary and elastic effects are relevant to a variety of applications from micro- and nano-scale manufacturing to biological systems. In this work, we investigate capillary flows in flexible, millimeter-scale cylindrical elastic tubes. We demonstrate that surface tension can cause sufficiently flexible tubes to collapse and coalesce spontaneously through non-axisymmetric buckling, and develop criteria for the initial deformation and complete collapse of a circular tube. Experimental results for capillary rise and evaporation of a liquid in a flexible tube are presented, and several regimes are seen for the equilibrium state of a flexible tube deforming under capillary pressure. Deformations of the tube walls are measured in different regimes and compared with a shell theory model. Analysis and experimental results show that despite the complex and non-axisymmetric deformed shapes of cylindrical structures, the elastocapillary length used in previous literature for flat plates and sheets can be used to predict the behavior of flexible tubes.

Description: We present laboratory experimental results demonstrating that librational forcing of an ellipsoidal container of water can produce intense motions through the mechanism of a libration driven elliptical instability (LDEI). These libration studies are conducted using an ellipsoidal acrylic container filled with water. A particle image velocimetry method is used to measure the 2D velocity field in the equatorial plane over hundreds libration cycles for a fixed Ekman number, E = 2 × 10 −5 . In doing so, we recover the libration induced base flow and a time averaged zonal flow. Further, we show that LDEI in non-axisymmetric container geometries is capable of driving both intermittent and saturated turbulent motions in the bulk fluid. Additionally, we measure the growth rate and amplitude of the LDEI induced excited flow in a fully ellipsoidal container at more extreme parameters than previously studied [Noir et al. , “Experimental study of libration-driven flows in nonaxisymmetric containers,” Phys. Earth Planet. Inter. 204-205 , 1 (2012); Cébron et al. , Phys. Fluids 24 , 061703, “Libration driven elliptical instability,” (2012)]. Excitation of bulk filling turbulence by librational forcing provides a mechanism for transferring rotational energy into turbulent fluid motion and thus can play an important role in the thermal evolution, interior dynamics, and magneto-hydrodynamics of librating bodies, as appear to be common in solar system settings [e.g., Comstock and Bills, “A solar system survey of forced librations in longitude,” J. Geophys. Res. Planets 108 , 1 (2003)].

Description: We numerically study the displacement flow of two iso-viscous Newtonian fluids in an inclined two-dimensional channel, formed by two parallel plates. The results are complementary to our previous studies on displacement flows in pipes and channels. The heavier displacing fluid moves the lighter displaced fluid in the downward direction. Three dimensionless groups largely describe these flows: the densimetric Froude number ( Fr ), the Reynolds number ( Re ), and the duct inclination (β). As a first order approximation, we are able to classify different flow regimes phenomenologically in a two-dimensional ( Fr ; Re cosβ/ Fr )-plane and provide leading order expressions for the transitions between different regimes. The stabilizing and/or de-stabilizing effects of the imposed mean flow on buoyant exchange flows (zero imposed velocity) are described for a broad range of dimensionless parameters.

Description: Compressible granular materials are involved in many applications, some of them being related to energetic porous media. Gas permeation effects are important during their compaction stage, as well as their eventual chemical decomposition. Also, many situations involve porous media separated from pure fluids through two-phase interfaces. It is thus important to develop theoretical and numerical formulations to deal with granular materials in the presence of both two-phase interfaces and gas permeation effects. Similar topic was addressed for fluid mixtures and interfaces with the Discrete Equations Method (DEM) [R. Abgrall and R. Saurel, “Discrete equations for physical and numerical compressible multiphase mixtures,” J. Comput. Phys. 186 (2), 361-396 (2003)] but it seemed impossible to extend this approach to granular media as intergranular stress [K. K. Kuo, V. Yang, and B. B. Moore, “Intragranular stress, particle-wall friction and speed of sound in granular propellant beds,” J. Ballist. 4 (1), 697-730 (1980)] and associated configuration energy [J. B. Bdzil, R. Menikoff, S. F. Son, A. K. Kapila, and D. S. Stewart, “Two-phase modeling of deflagration-to-detonation transition in granular materials: A critical examination of modeling issues,” Phys. Fluids 11 , 378 (1999)] were present with significant effects. An approach to deal with fluid-porous media interfaces was derived in Saurel et al. [“Modelling dynamic and irreversible powder compaction,” J. Fluid Mech. 664 , 348-396 (2010)] but its validity was restricted to weak velocity disequilibrium only. Thanks to a deeper analysis, the DEM is successfully extended to granular media modelling in the present paper. It results in an enhanced version of the Baer and Nunziato [“A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials,” Int. J. Multiphase Flow 12 (6), 861-889 (1986)] model as symmetry of the formulation is now preserved. Several computational examples are shown to validate and illustrate method’s capabilities.

Description: We perform a theoretical and numerical study of the Coulomb-driven electroconvection flow of a dielectric liquid between two coaxial cylinders. The specific case, where the inner to outer diameter ratio is 0.5, is analyzed. A strong unipolar injection of ions either from the inner or outer cylinder is considered to introduce free charge carriers into the system. A finite volume method is used to solve all governing equations including Navier-Stokes equations and a simplified set of Maxwell’s equations. The flow is characterized by a subcritical bifurcation in the finite amplitude regime. A linear stability criterion and a nonlinear one that correspond to the onset and stop of the flow motion, respectively, are linked with a hysteresis loop. In addition, we also explore the behavior of the system for higher values of the stability parameter. For inner injection, we observe a transition between the patterns made of 7 and 8 cells, before an oscillatory regime is attained. Such a transition leads to a second finite amplitude stability criterion. A simple modal analysis reveals that the competition of different modes is at the origin of this behavior. The charge density, as well as velocity field distributions is provided to help understand the bifurcation behavior.

Description: Spur and groove (SAG) formations are found on the forereefs of many coral reefs worldwide. Modeling results have shown that SAG formations together with shoaling waves induce a nearshore Lagrangian circulation pattern of counter-rotating circulation cells, but these have never been observed in the field. We present results from two separate field studies of SAG formations on Palmyra Atoll which show their effect on waves to be small, but reveal a persistent order 1 cm/s depth-averaged Lagrangian offshore flow over the spur and onshore flow over the grooves. This circulation was stronger for larger, directly-incident waves and low alongshore flow conditions, consistent with predictions from modeling. Favorable forcing conditions must be maintained on the order of one hour to accelerate and develop the SAG circulation cells. The primary cross- and alongshore depth-averaged momentum balances were between the pressure gradient, radiation stress gradient and nonlinear convective terms, and the bottom drag was similar to values found on other reefs. The vertical structure of these circulation cells was previously unknown and the results show a complex horizontal offshore Lagrangian flow over the spurs near the surface driven by alongshore variability in radiation stress gradients. Vertical flow was downward over the spur and upward over the groove, likely driven by alongshore differences in bottom stress and not by vortex forcing. This article is protected by copyright. All rights reserved.

Description: Shipboard current measurements in the equatorial Indian Ocean in October and November of 2011 revealed oscillations in the meridional velocity with amplitude ~ 0.10m/s. These were clearest in a layer extending from ~300 to 600 m depth and had periods near 3 weeks. Phase propagation was upward. Measurements from two sequential time series at the equator, four meridional transects and one zonal transect are used to identify the oscillation as a Yanai wave packet and to establish its dominant frequency and vertical wavelength. The Doppler shift is accounted for, so that measured wave properties are translated into the reference frame of the mean zonal flow. We take advantage of the fact that, in the depth range where the wave signal was clearest, the time-averaged current and buoyancy frequency were nearly uniform with depth, allowing application of the classical theoretical representation of vertically propagating plane waves. Using the theory, we estimate wave properties that are not directly measured, such as the group velocity and the zonal wavelength and phase speed. The theory predicts a vertical energy flux that is comparable to that carried by midlatitude near-inertial waves. We also quantify the wave-driven meridional heat flux and the Stokes drift. This article is protected by copyright. All rights reserved.

Description: Understanding variability in the chlorophyll-specific absorption of marine phytoplankton, a ph * Chl (λ) , is essential for primary production modelling, calculation of underwater light field characteristics, and development of algorithms for remote sensing of chlorophyll concentrations. Previous field and laboratory studies have demonstrated significant apparent variability in a ph * Chl (λ) for natural samples and algal cultures. However, the potential impact of measurement uncertainties on derived values of a ph * Chl (λ) has received insufficient study. This study presents an analysis of measurement uncertainties for a data set collected in the Ligurian Sea in Spring and assesses the impact on estimates of a ph * Chl (λ) . It is found that a large proportion of apparent variability in this set of a ph * Chl (λ) can be attributed to measurement errors. Application of the same analysis to the global NOMAD data set suggests that a significant fraction of variability in a ph * Chl (λ) may also be due to measurement errors. This article is protected by copyright. All rights reserved.

Description: The South Atlantic is an important pathway for the inter-basin exchanges of heat and freshwater with strong influence on the global meridional overturning stability and variability. Along the 34°S parallel, a quarterly, high resolution XBT transect (AX18) samples the temperature structure in the upper ocean. The AX18 transect has been shown to be a useful component of a meridional overturning monitoring system of the region. However, a feasible, cost-effective design for an XBT-based system has not yet been developed. Here we use a high-resolution ocean assimilation product to simulate an XBT-based observational system across the South Atlantic. The sensitivity of the meridional heat transport, meridional overturning circulation, and geostrophic velocities to key observational and methodological assumptions is studied. Key assumptions taken into account are horizontal and temporal sampling of the transect, salinity and deep temperature inference, as well as the level of reference for geostrophic velocities. With the current sampling strategy, the largest errors in the meridional overturning and heat transport estimations are the reference (barotropic) velocity and the western boundary resolution. We show how altimetry can be used along with hydrography to resolve the barotropic component of the flow. We use the results obtained by the state estimation under observational assumptions to make recommendations for potential improvements in the AX18 transect implementation. This article is protected by copyright. All rights reserved.

Description: Phytoplankton patchiness, namely the heterogeneous distribution of microalgae over multiple spatial scales, dramatically impacts marine ecology. A spectacular example of such heterogeneity occurs in thin phytoplankton layers (TPLs), where large numbers of photosynthetic microorganisms are found within a small depth interval. Some species of motile phytoplankton can form TPLs by gyrotactic trapping due to the interplay of their particular swimming style (directed motion biased against gravity) and the transport by a flow with shear along the direction of gravity. Here we consider gyrotactic swimmers in numerical simulations of the Kolmogorov shear flow, both in laminar and turbulent regimes. In the laminar case, we show that the swimmer motion is integrable and the formation of TPLs can be fully characterized by means of dynamical systems tools. We then study the effects of rotational Brownian motion or turbulent fluctuations (appearing when the Reynolds number is large enough) on TPLs. In both cases, we show that TPLs become transient, and we characterize their persistence.

Description: Observations of the nearbed velocity field over a rippled sediment bed under asymmetric wave forcing conditions were collected using a submersible particle image velocimetry (PIV) system. To examine the role of bedform-induced dynamics in the total momentum transfer, a double-averaging technique was implemented on the two-dimensional time-dependent velocity field by means of the full momentum equation. This approach allows for direct determination of the bedform-induced stresses, ie. stresses that arise due to the presence of bedforms, which are zero in flat bed conditions. This analysis suggests that bedform-induced stresses are closely related to the presence of coherent motions and may be partitioned from the turbulent stresses. Inferences of stress provided by a bedload transport model suggest that total momentum transfer obtained from the double-averaging technique is capable of reproducing bedform mobilization. Comparisons between the total momentum transfer and stress estimates obtained from local velocity profiles show significant variability across the ripple, and suggest that an array of sensors is necessary to reproduce bedform evolution. The imbalance of momentum obtained by resolving the different terms constituting the near-bed momentum balance (i.e. acceleration deficit, stress gradient, and bedform-induced skin friction) provides an estimate of the bedform-induced pressure that is consistent with flow separation. This analysis reveals three regions in the flow: the free-stream, where all terms are relatively balanced; the near-bed, where momentum imbalance is significant during flow weakening; and below ripple crests, where bedform-induced pressure is the leading order mechanism.

Description: An analysis is presented for the spatial and intensity distributions of North Atlantic extreme atmospheric events crossing the buoyant Amazon-Orinoco freshwater plume. The sea surface cooling amplitude in the wake of an ensemble of storm tracks travelling in that region is estimated from satellite products for the period 1998-2012. For the most intense storms, cooling is systematically reduced by ~50% over the plume area compared to surroundings open ocean waters. Historical salinity and temperature observations from in situ profiles indicate that salt-driven vertical stratification, enhanced oceanic heat content and barrier-layer presence within the plume waters are likely key oceanic factors to explain these results. Satellite SMOS surface salinity data combined with in situ observations are further used to detail the oceanic response to Category 4 hurricane Igor in 2010. Argo and satellite measurements confirm the haline stratification impact on the cooling inhibition as the hurricane crossed the river plume. Over this region, the SSS mapping capability is further tested and demonstrated to monitor the horizontal distribution of the vertical stratification parameter. SMOS SSS data can thus be used to consistently anticipate the cooling inhibition in the wake of TCs travelling over the Amazon-Orinoco plume region.

Description: The mixed layer salinity (MLS) budget of the tropical Pacific is investigated using results from a model of the Consortium for Estimating the Circulation and Climate of the Ocean (ECCO). The results focusing on the western Pacific freshwater pool indicate that the long-term averaged surface freshwater flux is well balanced by ocean dynamics, in which the subsurface processes account for the major part. The MLS budget shows significant seasonal and interannual variability, as a consequence of interplay among surface freshwater flux, advection, mixing, and vertical entrainment. On seasonal time scale, both the MLS and mixed layer depth are largely controlled by surface freshwater flux. The opposite phase between the subsurface processes and the barrier layer thickness confirms the important influence of the barrier layer on vertical mixing and entrainment from below. On interannual time scale, all the MLS budget terms show significant ENSO signal, which in turn is highly correlated with the salinity front and barrier layer thickness in the equatorial Pacific.

Description: To study the relationship of solar heat input into the Arctic open water and the variations of sea ice extent, improved satellite based estimates of shortwave radiative (SWR) fluxes and most recent observations of ice extent are used. The SWR flux estimates are based on observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and from the Advanced Very High Resolution Radiometer (AVHRR) for the period of 1984-2009. Ice extent information at 25-km resolution comes from Nimbus-7 SMMR and DMSP SSM/I Passive Microwave Data as generated with the NASA Team algorithm developed by the Oceans and Ice Branch, Laboratory for Hydrospheric Processes, NASA Goddard Space Flight Center. The trends of the solar heat input into the ocean and the open water fraction for 1984-2009 are found to be positive: 0.3%/year and 0.8% / year respectively at a 99% confidence level. There is an obvious transition region separating the 26 years into two periods: one with moderate change: 1984-2002; one with an abrupt growth in both solar heat input and open water fraction: 2003-2009. The impact of the observed changes on the reduction of winter ice growth in 2007 is estimated to be about 44 cm, and a delay in fall freeze-up as about 10~36 days.

Description: Numerous studies have shown the primary importance of wind stress curl in coastal upwelling dynamics. The main goal of this new analysis is to describe the QuikSCAT surface wind stress curl at various scales in the Benguela and Canary upwelling systems. The dominant spatial pattern is characterized by cyclonic curl near continental boundaries and anticyclonic curl offshore, in association with equatorward alongshore (upwelling favorable) wind stress. At a smaller scale, we demonstrate the sensitivity of the QuikSCAT wind stress curl to coastal processes related to sea surface temperature (SST) mesoscale fluctuations by presenting a linear relationship between the curl and crosswind SST gradients. Despite the spatial and temporal sensitivity of the underlying thermal coupling coefficient, a local analysis of the fraction of the curl ascribed to SST variability shows that SST is a main driver of the wind stress curl variability and magnitude over the upwelling extension zone (~100 to 300 km from the coast) in both the Canary and Benguela systems. Closer to the shore, the curl patterns derived from QuikSCAT observations are only loosely related to SST-wind interactions. As a working hypothesis, they can also be associated with the coastline geometry and orographic effects that are likely to play an important role in local cooling processes.

Description: To provide a probable explanation on the field observed rapid sedimentation process near river mouths, we investigate the convective sedimentation in stably stratified saltwater using 3D numerical simulations. Guided by the linear stability analysis (Yu et al. 2013, J. Geophys. Res. Oceans, 118, 256-272), this study focuses on the nonlinear interactions of several mechanisms, which lead to various sediment finger patterns, and the effective settling velocity for sediment ranging from clay (single-particle settling velocity V 0 = 0.0036 and 0.0144 mm/s, or particle diameter d = 2 and 4 μ m) to silt ( V 0 =0.36 mm/s, or d =20 μ m). For very fine sediment with V 0 =0.0036 mm/s, the convective instability is dominated by double diffusion, characterized by millimeter-scale fingers. Gravitational settling slightly increases the growth rate; however, it has notable effect on the downward development of vertical mixing shortly after the sediment interface migrates below the salt interface. For sediment with V 0 = 0.0144 mm/s, Rayleigh-Taylor instabilities become dominant before double-diffusive modes grow sufficiently large. Centimeter-scale and highly asymmetric sediment fingers are obtained due to nonlinear interactions between different modes. For sediment with V 0 = 0.36 mm/s, Rayleigh-Taylor mechanism dominates and the resulting centimeter-scale sediment fingers show a plume-like structure. The flow pattern is similar to that without ambient salt stratification. Rapid sedimentation with effective settling velocity on the order of 1 cm/s is likely driven by convective sedimentation for sediment with V 0 greater than 0.1 mm/s at concentration greater than 10 to 20 g/L.

Description: The surface layer of the southeast Pacific Ocean (SEP) requires an input of cold, fresh water to balance heat gain and evaporation from air-sea fluxes. Models typically fail to reproduce the cool sea surface temperatures (SST) of the SEP, limiting our ability to understand the variability of this climatically important region. We estimate the annual heat budget of the SEP for the period 2004 - 2009, using data from the upper 250 m of the Stratus mooring, located at 85°W 20°S, and from Argo floats. The surface buoy measures meteorological conditions and air-sea fluxes; the mooring line is heavily instrumented, measuring temperature, salinity, and velocity at more than 15 depth levels. We use a new method for estimating the advective component of the heat budget that combines Argo profiles and mooring velocity data, allowing us to calculate monthly profiles of heat advection. Averaged over the 6 year study period, we estimate a cooling advective heat flux of -41 ± 29 W m-2, accomplished by a combination of the mean gyre circulation, Ekman transport, and eddies. This compensates for warming fluxes of 32 ± 4 W m-2 due to air-sea fluxes and 7 ± 9 W m-2 due to vertical mixing and Ekman pumping. A salinity budget exhibits a similar balance, with advection of freshwater (-60 psu m) replenishing the freshwater lost through evaporation (47 psu m) and Ekman pumping (14 psu m).

Description: The model of gas bubble growth in high-viscous gas-saturated magmatic melt, subjected to rapid decompression, is presented in the current study. It is shown that consideration of unsteady character of the process is extremely important in a wide range of supersaturation. The analytical solution is found for the profile of dissolved gas concentration and the rate of bubble growth. The model of kinetics of overall degassing is developed. This model is based on distinguishing the so-called “forbidden” zone in the melt volume with suppressed formation of the new nucleation sites. The simple analytical dependences of the number of nucleating bubbles and typical nucleation time on the value of initial decompression were derived together with time dependence of volumetric concentration of the gas phase. Our results match the available experimental data.

Description: We investigate the effect of viscosity contrast on the stability of gravitationally unstable, diffusive layers in porous media. Our analysis helps evaluate experimental observations of various diffusive (boundary) layer models that are commonly used to study the sequestration of CO 2 in brine aquifers. We evaluate the effect of viscosity contrast for two basic models that are characterized with respect to whether or not the interface between CO 2 and brine is allowed to move. We find that diffusive layers are in general more unstable when viscosity decreases with depth within the layer compared to when viscosity increases with depth. This behavior is in contrast to the one associated with the classical displacement problem of gravitationally unstable diffusive layers that are subject to mean flow. For the classical problem, a greater instability is associated with the displacement of a more viscous, lighter fluid along the direction of gravity by a less viscous, heavier fluid. We show that the contrasting behavior highlighted in this study is a special case of the classical displacement problem that depends on the relative strength of the displacement and buoyancy velocities. We demonstrate the existence of a critical viscosity ratio that determines whether the flow is buoyancy dominated or displacement dominated. We explain the new behaviors in terms of the interaction of vorticity components related to gravitational and viscous effects.

Description: Acoustic-gravity waves are compression-type waves propagating with amplitudes governed by the restoring force of gravity. They are generated, among others, by wind-wave interactions, surface waves interactions, and submarine earthquakes. We show that acoustic-gravity waves contribute to deep ocean currents and circulation; they cause chaotic flow trajectories of individual water parcels, which can be transported by a few centimetres per second.

Description: The development of a round liquid jet under the influence of a confined coaxial flow of an immiscible liquid of comparable density (central to annular flow density ratio of 8:10) was investigated in the vicinity of the nozzle exit. Two flow regimes were considered; one where the annular flow is faster than the central jet, so the central liquid jet is accelerated and one where the annular flow is slower, so the central liquid jet is decelerated. The central jet was visualised by high speed photography. Three modes of jet development were identified and classified in terms of the Reynolds number, Re, of the central jet which was in the range of 525 〈 Re 〈 2725, a modified definition of the Weber number, We, which allows the distinction between accelerating and deceleration flows and was in the range of −22 〈 We 〈 67 and the annular to central Momentum Ratio, MR, of the two streams which was in the range of 3.6 〈 MR 〈 91. By processing the time resolved jet images using Proper Orthogonal Decomposition (POD), it was possible to reduce the description of jet morphology to a small number of spatial modes, which isolated the most significant morphologies of the jet development. In this way, the temporal and spatial characteristics of the instabilities on the interface were clearly identified which highlights the advantages of POD over direct observation of the images. Relationships between the flow parameters and the interfacial waves were established. The wavelength of the interfacial instability was found to depend on the velocity of the fastest moving stream, which is contrary to findings for fluids with large density differences.

Description: Interaction of a vortex ring impinging on multiple permeable screens orthogonal to the ring axis was studied to experimentally investigate the persistence and decay of vortical structures inside the screen array using digital particle image velocimetry in a refractive index matched environment. The permeable screens had porosities (open area ratios) of 83.8%, 69.0%, and 55.7% and were held by a transparent frame that allowed the screen spacing to be changed. Vortex rings were generated using a piston-cylinder mechanism at nominal jet Reynolds numbers of 1000, 2000, and 3000 with piston stroke length-to-diameter ratios of 2 and 3. The interaction of vortex rings with the porous medium showed a strong dependence of the overall flow evolution on the screen porosity, with a central flow being preserved and vortex ring-like structures (with smaller diameter than the primary vortex ring) being generated near the centerline. Due to the large rod size used in the screens, immediate reformation of the transmitted vortex ring with size comparable to the primary ring (as has been observed with thin screens) was not observed in most cases. Since the screens have lower complexity and high open area ratios, centerline vortex ring-like flow structures formed with comparable size to the screen pore size and penetrated through the screens. In the case of low porosity screens (55.7%) with large screen spacing, re-emergence of large scale (large separation), weak vortical structures/pairs (analogous to a transmitted vortex ring) was observed downstream of the first screen. Additional smaller scale vortical structures were generated by the interaction of the vortex ring with subsequent screens. The size distribution of the generated vortical structures were shown to be strongly affected by porosity, with smaller vortical structures playing a stronger role as porosity decreased. Finally, porosity significantly affected the decay of total energy, but the effect of screen spacing decreased as porosity decreased.

Description: Vortex cavitation forming in the leading-edge vortices of a delta wing was examined to determine how the individual cavitation bubbles incepted, grew, interacted with the underlying vortical flow and produced acoustic tones. The non-cavitating vortical flow over the delta wing was chosen to be similar to those previously reported in the literature. It was found that vortex breakdown was unaffected by the presence of incipient and developed vortex cavitation bubbles in the vortex core. While some cavitation bubbles incepted, grew, and collapsed relatively quickly, others reached an equilibrium position wherein the bubble tip was stationary in the laboratory frame at a particular location along the vortex axis. For a given attack angle, the equilibrium location moved upstream with a reduction in free stream cavitation number. It is shown that the existence of these stationary vortex bubbles is possible when there is a balance between the axial growth of the bubble along the vortex axis and the opposite motion of the axial jetting flow in the vortex core, and only a single equilibrium position is possible along the axially evolving vortex for a given free stream cavitation number. These transient and stationary vortex bubbles emit significant cavitation noise upon inception, growth, and collapse. The spectral content of the noise produced was expected to be related to the interaction of the bubble with the surrounding vortical flow in a manner similar to that reported in previous studies, where sustained tones were similar to the underlying vortex frequency. However, in the present study, the dominant frequency and higher harmonics of the tones occur at a higher frequency than that of the underlying vortex. Hence, it is likely that the highly elongated stationary bubbles have higher-order volume oscillations compared to the two-dimensional radial mode of the vortex cores of vortex cavitation bubbles with much smaller diameter-to-length ratios.

Description: Laminar flow over a periodic array of cylindrical surface roughness elements is simulated with an immersed boundary spectral method both to validate the method for subsequent studies and to examine how persistent streamwise vortices are introduced by a low Reynolds number roughness element. Direct comparisons are made with prior studies at a roughness-based Reynolds number Re k (= U ( k ) k / ν ) of 205 and a diameter to spanwise spacing ratio d / λ of 1/3. Downstream velocity contours match present and past experiments very well. The shear layer developed over the top of the roughness element produces the downstream velocity deficit. Upstream of the roughness element, the vortex topology is found to be consistent with juncture flow experiments, creating three cores along the recirculation line. Streamtraces stemming from these upstream cores, however, have unexpectedly little effect on the downstream flowfield as lateral divergence of the boundary layer quickly dissipates their vorticity. Long physical relaxation time of the recirculating wake behind the roughness remains a prominent issue for simulating this type of flowfield.

Description: The Berau Continental Shelf is located close to the Equator in the Indonesian Archipelago, hosting a complex of coral reefs along its oceanic edge. The Berau coral reefs have a very high biodiversity, but the area is under serious risk due to river-derived nutrients and sediments. The region is characterized by weak winds, moderate tides and almost absent Coriolis forcing. Existing knowledge about river plume behaviour in tropical environments is limited. The aim of this paper is to investigate the influence of the subtle physical forcing on the dynamics of the Berau river plume. A three-dimensional model (ECOMSED) was calibrated with observational data. The model was forced by freshwater input from the Berau river distributaries, tides at the open boundaries and measured hourly wind. The model reproduces the freshwater dynamics on the shelf adequately and highlights that the river plume spreads symmetrically for river forcing only. Tides cause vertical mixing and suppress the cross-shelf spreading of the river plume. However, the spreading of the river plume over the shelf is mainly controlled by the weak monsoonal winds, resulting in a seasonal development. During the Southeast Monsoon, the southerly winds push the plume northeastward and cause a stratified water column in the northern part of the continental shelf. Northerly winds during the Northwest Monsoon disperse the plume to the south, promoting a vertically well-mixed water column. The results can be used to predict the possible impact of land-use changes in the steadily developing Berau region on coral reef health. This article is protected by copyright. All rights reserved.

Description: This study investigates the effect of adaptive (or targeted) observation on improving the mid-range (30 days) forecast skill of ocean state of the South China Sea (SCS). A region associated with the South China Sea Western Boundary Current (SCSWBC) is chosen as the “target” of the adaptive observation. The Conditional Nonlinear Optimal Perturbation (CNOP) approach is applied to a 3-dimensional ocean model and its adjoint model for determining the sensitive region. Results show that the initial errors in the sensitive region determined by the CNOP approach have significant impacts on the forecast of ocean state in the target region; thus, reducing these initial errors through adaptive observation can lead to a better 30-day prediction of ocean state in the target region. Our results suggest that implementing adaptive observation is an effective and cost-saving way to improve an ocean model's forecast skill over the SCS. This article is protected by copyright. All rights reserved.

Description: We present a year-round assessment of the hydrographic variability within the East Greenland Coastal Current on the Greenland shelf from five synoptic crossings and four years of moored hydrographic data. From the five synoptic sections the current is observed as a robust, surface intensified flow with a total volume transport of 0.66 ± 0.18 Sv and a freshwater transport of 42 ± 12 mSv. The moorings showed heretofore unobserved variability in the abundance of Polar and Atlantic water masses in the current on synoptic scales. This is exhibited as large vertical displacement of isotherms (often greater than 100 m). Seasonally, the current is hemmed into the coast during the fall by a full depth Atlantic Water layer that has penetrated onto the inner shelf. The Polar Water layer in the current then thickens through the winter and spring seasons increasing the freshwater content in the current; the timing implies that this is probably driven by the seasonally varying export of freshwater from the Arctic and not the local runoff from Greenland. The measured synoptic variability is enhanced during the winter and spring period due to a lower halocline and a concurrent enhancement in the along-coast wind speed. The local winds force much of the high-frequency variability in a manner consistent with downwelling, but variability distinct from downwelling is also visible. This article is protected by copyright. All rights reserved.

Description: The performance of several numerical ocean models is assessed with respect to their simulation of sea surface height (SSH) in the Arctic Ocean, and the main patterns of SSH variability and their causes over the past 40 years (1970-2009) are analyzed. In comparison to observations, all tested models broadly reproduce the mean SSH in the Arctic and reveal a good correlation with both tide gauge data and SSH anomalies derived from satellite observations. Although the models do not represent the positive Arctic SSH trend observed over the last two decades, their interannual-to-decadal SSH variability is in reasonable agreement with available measurements. Focusing on results from one of the models for a detailed analysis it is shown that the decadal-scale SSH variability over shelf areas and deep parts of the Arctic Ocean have pronounced differences that are determined mostly by salinity variations. A further analysis of the three time periods 1987-1992, 1993-2002 and 2003-2009, corresponding to the transition times between cyclonic and anticyclonic regimes of the atmospheric circulation over the Arctic, revealed an unusual increase of SSH in the Amerasian basin during 2003-2009. Results from this model support the recent finding that the increase is caused mainly by changes in freshwater content brought about by the freshwater export through the Canadian Arctic Archiplago and increased Ekman pumping in the Amerasian basin and partly by lateral freshwater transport changes, leading to a re-distribution of low-salinity shelf water. Overall we show that present day models can be used for investigating the reasons for low-frequency SSH variability in the region. This article is protected by copyright. All rights reserved.

Description: In this study, a neural network-based four-band model (NNFM) for the global oceanic and coastal waters has been developed in order to retrieve the total absorption coefficients a ( λ ). The applicability of the quasi-analytical algorithm (QAA) and NNFM models is evaluated by five independent datasets. Based on the comparison of a ( λ ) predicted by these two models with the field measurements taken from the global oceanic and coastal waters, it was found that both the QAA and NNFM models had good performances in deriving a ( λ ), but that the NNFM model works better than the QAA model. The results of the QAA model-derived a ( λ ), especially in highly turbid waters with strong backscattering properties of optical activity, was found to be lower than the field measurements. The QAA and NNFM models-derived a ( λ ) could be obtained from the MODIS data after atmospheric corrections. When compared with the field measurements, the NNFM model decreased by a 0.86 to 24.15% uncertainty (root mean square relative error) of the estimation from the QAA model in deriving a ( λ ) from the Bohai, Yellow, and East China seas. Finally, the NNFM model was applied to map the global climatological seasonal mean a (443) for the time range of July, 2002 to May, 2014. As expected, the a (443) value around the coastal regions was always larger than the open ocean around the equator. Viewed on a global scale, the oceans at a high latitude exhibited higher a (443) values than those at a low latitude. This article is protected by copyright. All rights reserved.

Description: In this study, we investigate the interannual variability of the sea surface temperature (SST) in the South China Sea (SCS) associated with two types of El Niño, namely, the eastern Pacific (EP) El Niño and the central Pacific (CP) El Niño. First, double warm peaks can occur during both types of El Niño events in the SCS. However, the strong warm basin mode can only develop in the EP El Niño, while the warm semi-basin mode exists during the CP El Niño. Associated with an anomalous positive (negative) net surface heat flux in the EP (CP) El Niño, along with a shallower thermocline with weaker (stronger) northeasterly wind anomalies, the SST anomalies become warmer (cooler) in the developing autumn. Over the background of cooling SST in autumn of CP El Niño, therefore, only a weak warming can occur in the subsequent years, which is limited in the western boundary area under the forcing of warm ocean advection. Second, the SST oscillation periods are different in these two types of El Niño. The SST evolution in the EP El Niño is negative-positive with a quasi-biennial oscillation, but that in the CP El Niño is positive-negative-positive-negative with an annual oscillation. It seems that the double cooling in the CP El Niño is phase-locked to the late autumn season. This article is protected by copyright. All rights reserved.

Description: Following the idea that analysis of in-situ information in the salt budget could be used as a surrogate for global “ocean rain gauge”, the annual mean oceanic net freshwater flux (E-P) was estimated from the Argo profiles and the wind stress data on a global scale. The comparison between the independent E-P estimation from Argo and the E-P product sets, including the combination of precipitation from TRMM, GPCP, CMAP and evaporation from OAFlux, GSSTF3 and IFREMER and E-P set from NEWS formed from satellite, generally show similar spatial patterns, particularly on the larg scale. However, there are differences among the different satellite-based E-P estimates and between satellite estimates and independent in-situ estimates. Based on the pattern correlation and the RMSD, the evaporation and precipitation from OAFlux and TRMM agrees best with the E-P estimated from the independent Argo-based estimates.

Description: The rapidly changing East Siberian Arctic Shelf (ESAS) receives large amounts of terrestrial organic carbon (OC) from coastal erosion and Russian-Arctic rivers. Climate warming increases thawing of coastal Ice Complex Deposits (ICD) and can change both the amount of released OC, as well as its propensity to be converted to greenhouse gases (fueling further global warming) or to be buried in coastal sediments. This study aimed to unravel the susceptibility to degradation, and transport and dispersal patterns of OC delivered to the ESAS. Bulk and molecular radiocarbon analyses on surface particulate matter (PM), sinking PM and underlying surface sediments illustrate the active release of old OC from coastal permafrost. Molecular tracers for recalcitrant soil OC showed ages of 3.4-13 14 C-ky in surface PM and 5.5-18 14 C-ky in surface sediments. The age difference of these markers between surface PM and surface sediments is larger (i) in regions with low OC accumulation rates, suggesting a weaker exchange between water column and sediments, and (ii) with increasing distance from the Lena River, suggesting preferential settling of fluvially-derived old OC nearshore. A dual-carbon end-member mixing model showed that (i) contemporary terrestrial OC is dispersed mainly by horizontal transport while being subject to active degradation, (ii) marine OC is most affected by vertical transport and also actively degraded in the water column, and (iii) OC from ICD settles rapidly and dominates surface sediments. Preferential burial of ICD-OC released into ESAS coastal waters might therefore lower the suggested carbon cycle climate feedback from thawing ICD permafrost. This article is protected by copyright. All rights reserved.

Description: The Arctic sea ice cover undergoes large changes over an annual cycle. In winter and spring the ice cover consists of large, snow-covered plate-like ice floes, with very little open water. By the end of summer the snow cover is gone and the large floes have broken into a complex mosaic of smaller, rounded floes surrounded by a lace of open water. This evolution strongly affects the distribution and fate of the solar radiation deposited in the ice-ocean system and consequently the heat budget of the ice cover. In particular, increased floe perimeter can result in enhanced lateral melting. We attempt to quantify the floe evolution process through the concept of a floe size distribution that is modified by lateral melting and floe breaking. A time series of aerial photographic surveys made during the SHEBA field experiment is analyzed to determine evolution of the floe size distribution from spring through summer. Based on earlier studies, we assume the floe size cumulative distribution could be represented by a power law D -α where D is the floe diameter. The exponent α as well as the number density of floes N tot are estimated from measurements of total ice area and perimeter. As summer progressed, there was an increase in α as the size distribution shifted towards smaller floes and the number of floes increased. Lateral melting causes the distribution to deviate from a power law for small floe sizes. This article is protected by copyright. All rights reserved.

Description: The responses of the ocean planktonic ecosystem to finite-amplitude perturbations are investigated using an ocean planktonic ecosystem model. Through changing the higher predation rate on zooplankton, multiple equilibria of the model, namely “high-nutrient” and “low-nutrient” states, are obtained under certain parameter values. Based on these states, the perturbations with maximum nonlinear growth are determined using the conditional nonlinear optimal perturbation (CNOP) method. The linear and nonlinear evolutions of the CNOP perturbation are compared. The results show that the nonlinear evolution of the perturbation leads to predator–prey oscillations with larger amplitude than the linear evolution. Besides, after the perturbation amplitude exceeds a critical value, the nonlinear evolution of the perturbation will induce the linearly stable ecosystem state to lose the stability and become nonlinearly unstable. This implies that nonlinear processes have important impacts on the stability of the ecosystem. Specifically, we identify the nonlinear processes related to zooplankton grazing to impact the stability most for the high-nutrient state, while for the low-nutrient state the main nonlinear process affecting the stability is the uptake process. These results help to improve our understanding of the sensitivity of the oceanic ecosystem model to finite-amplitude perturbations and the underlying nonlinear stability properties. This article is protected by copyright. All rights reserved.

Description: We evaluated the impacts of summertime mesoscale convective systems (MCS) on the heat balance and diel surface mixed layer (SML) dynamics of the Brazilian Amazon's Tucuruí Hydroelectric Reservoir (THR). We used a synergistic approach that combines in situ data, remote sensing data and three-dimensional (3D) modeling to investigate the typical behavior of the components of the heat balance and the SML dynamics. During the study period (the austral summer of 2012-2013), 22 days with MCS activity were identified. These events occurred approximately every 4 days, and they were most frequent during January (50% of the observations). An analysis of local meteorological data showed that when MCS occur, the environmental conditions at THR change significantly ( p-value 〈 0.01). The net longwave flux, which was the heat balance component most strongly impacted by MCS, increased more than 32% on days with MCS activity. The daily integrated heat balance became negative (-54 W m -2 ) on MCS days, while the balance was positive (19 W m -2 ) on non MCS days. In response to the changes in the heat balance, the SML dynamics changed when a MCS was over the THR. The SML depth was typically 28% higher on the days with MCS (~1.6 m) compared with the days without MCS (~1.3 m). The results indicate that MCS are one of the main meteorological disturbances driving the heat balance and the mixing dynamics of Amazonian hydroelectric reservoirs during the summer. These events may have implications for the water quality and greenhouse gas emissions of Amazonian reservoirs. This article is protected by copyright. All rights reserved.

Description: Solutal Marangoni instability (SMI) is investigated both in 2D and 3D using a combined Cahn-Hilliard and Navier-Stokes model in a finite system. Fe-Sn is chosen as a representative alloy system since the phase diagram reveals a region with a miscibility gap, where two liquid phases, namely, the Fe-rich phase L 1 and the Sn-rich phase L 2 , are in chemical equilibrium. In 3D, considering a perturbed liquid cylinder ( L 2 phase) with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H × H (length × width × height) surrounded by the phase L 1 , we find that the perturbation induced Marangoni flow is either clockwise or anti-clockwise depending on the mean curvature difference between the convex and concave regions which is affected by the ratio of λ/ R 0 . The critical ratio of λ/ R 0 for SMI is shown to be invariant for different Marangoni numbers as well as independent of the geometrical properties of the L 1 phase. In 2D, a perturbed liquid pipe with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H (length × height) is taken into account. Due to different curvature constitution, the critical ratio of λ/ R 0 for SMI depends on the height of the L 1 phase.

Description: For several years, a promising Plasma Synthetic Jet actuator for high-speed flow control has been under development at ONERA. So far, its confined geometry and small space-time scales at play have prevented its full experimental characterization. Complementary accurate numerical simulations are then considered in this study in order to provide a complete aerothermodynamic description of the actuator. Two major obstacles have to be overcome with this approach: the modeling of the energy deposited by the electric arc and the accurate computation of the transient response of the cavity generating the pulsed jet. To solve the first problem, an Euler solver coupled with an electric circuit model was used to evaluate the energy deposition in the cavity. Such a coupling is performed by considering the electric field between the two electrodes. The second issue was then addressed by injecting these source terms in large Eddy simulations of the entire actuator. Aerodynamic results were finally compared with Schlieren visualizations. Using the proposed methodology, the temporal evolution of the jet front is remarkably well predicted.

Description: Snapshot and classical proper orthogonal decomposition (POD) are used to examine the large-scale, energetic motions in fully developed turbulent pipe flow at Re D = 47,000 and 93,000. The snapshot POD modes come in pairs, representing the same azimuthal mode number but with a simple phase shift. The first 10 snapshot POD modes, associated with the very large scale motions (VLSMs), contribute 43% of the average Reynolds shear stress, and for first 80 modes u ′ and v ′ are anti-correlated so that they all contribute to positive shear stress events. The attached motions are contained in the lower order modes, and detached motions do not appear until snapshot POD mode numbers ≥15. We find that snapshot POD can introduce mode mixing, which is avoided in classical POD. Classical POD also gives frequency information, confirming that the low order modes capture well the behavior of the very large scale motions.

Description: Deep water circulation in the Luzon Strait, which connects the Pacific deep circulation with the South China Sea throughflow, is investigated using a set of oceanographic observations combined with results from three numerical experiments. Both the in situ observations and the model show a deep water overflow in the Luzon Strait. Their results suggest that the deep Pacific water first flows into the Luzon Strait through the Bashi Channel (1.2 Sv, 1 Sv = 1×10 6 m 3 s -1 ) and the Taltung Canyon (0.4 Sv), then turns southward along the Luzon Trough, and finally enters the South China Sea primarily through two gaps in the Heng-Chun Ridge. Overall, the mean transport of the Luzon Strait overflow is about 1.5 Sv. Results from numerical experiments suggest that strong diapycnal mixing in the South China Sea and Luzon Strait, which sustains the baroclinic pressure gradient across the Luzon Strait, is the primary driving mechanism of the deep circulation in the Luzon Strait.

Description: Earlier studies indicate that during El Niño events the iceberg concentration increases in the east of the Pacific sector and in the west of the Atlantic sector of Southern Ocean, but decreases in the center of the Pacific sector. During La Niña the pattern of the iceberg concentration anomalies in these regions reverses. This iceberg redistribution is explained by anomalous winds and currents around an extensive positive atmospheric pressure anomaly that typically develops in the South-East Pacific during the warm El Niño-Southern Oscillation (ENSO) phase. In this study the results of iceberg observations during two cruises of the r/v “Akademik Fedorov” in Antarctica in January-February of 2008 (La Niña) and 2010 (El Niño) have been used to examine the consistency of changes in the iceberg distribution in the Southern Ocean related to El Niño events. The analysis of these observations has shown that in the Pacific Sector of Antarctica changes in the iceberg distribution between 2008 and 2010 followed the scenario outlined above and thus could be associated with the ENSO phase change. Contrary to earlier observations, the iceberg concentration in the Atlantic sector of Antarctica did not increase during 2010 El Niño. The latter is explained by a non-canonical type of 2010 El Niño, El Niño Modoki and associated atmospheric circulation pattern different from the canonical El Niño. Further analysis has shown that a more frequent occurrence of El Niño Modoki in recent years have resulted in weaker links between El Niño events and the Antarctic iceberg distribution.

Description: We explore the instabilities developed in a fluid in which viscosity depends on temperature. In particular, we consider a dependency that models a very viscous (and thus rather rigid) lithosphere over a convecting mantle. To this end, we study a 2D convection problem in which viscosity depends on temperature by abruptly changing its value by a factor of 400 within a narrow temperature gap. We conduct a study which combines bifurcation analysis and time-dependent simulations. Solutions such as limit cycles are found that are fundamentally related to the presence of symmetry. Spontaneous plate-like behaviors that rapidly evolve towards a stagnant lid regime emerge sporadically through abrupt bursts during these cycles. The plate-like evolution alternates motions towards either the right or the left, thereby introducing temporary asymmetries on the convecting styles. Further time-dependent regimes with stagnant and plate-like lids are found and described.

Description: We study the shock wave structure in a rarefied polyatomic gas based on a simplified model of extended thermodynamics in which the dissipation is due only to the dynamic pressure. In this case the differential system is very simple because it is a variant of Euler system with a new scalar equation for the dynamic pressure [T. Arima, S. Taniguchi, T. Ruggeri, and M. Sugiyama, Phys. Lett. A376, 2799–2803 (2012)]. It is shown that this theory is able to describe the three types of the shock wave structure observed in experiments: the nearly symmetric shock wave structure (Type A, small Mach number), the asymmetric structure (Type B, moderate Mach number), and the structure composed of thin and thick layers (Type C, large Mach number).

Description: In this paper, the scaling property of the inverse energy cascade and forward enstrophy cascade of the vorticity filed ω( x , y ) in two-dimensional (2D) turbulence is analyzed. This is accomplished by applying a Hilbert-based technique, namely Hilbert-Huang transform, to a vorticity field obtained from a 8192 2 grid-points direct numerical simulation of the 2D turbulence with a forcing scale k f = 100 and an Ekman friction. The measured joint probability density function p ( C , k ) of mode C i ( x ) of the vorticity ω and instantaneous wavenumber k ( x ) is separated by the forcing scale k f into two parts, which correspond to the inverse energy cascade and the forward enstrophy cascade. It is found that all conditional probability density function p ( C | k ) at given wavenumber k has an exponential tail. In the inverse energy cascade, the shape of p ( C | k ) does collapse with each other, indicating a nonintermittent cascade. The measured scaling exponent ζ ω I ( q ) is linear with the statistical order q , i.e., ζ ω I ( q ) = − q / 3 , confirming the nonintermittent cascade process. In the forward enstrophy cascade, the core part of p ( C | k ) is changing with wavenumber k , indicating an intermittent forward cascade. The measured scaling exponent ζ ω F ( q ) is nonlinear with q and can be described very well by a log-Poisson fitting: ζ ω F ( q ) = 1 3 q + 0.45 1 − 0 . 43 q . However, the extracted vorticity scaling exponents ζ ω ( q ) for both inverse energy cascade and forward enstrophy cascade are not consistent with Kraichnan's theory prediction. New theory for the vorticity field in 2D turbulence is required to interpret the observed scaling behavior.

Description: We report experimental observations of the controlled deformation of a dielectric liquid jet subjected to a local high-voltage electrostatic field in the direction normal to the jet. The jet deforms to the shape of an elliptic cylinder upon application of a normal electrostatic field. As the applied electric field strength is increased, the elliptic cylindrical jet deforms permanently into a flat sheet, and eventually breaks-up into droplets. We interpret this observation—the stretch of the jet is in the normal direction to the applied electric field—qualitatively using the Taylor-Melcher leaky dielectric theory, and develop a simple scaling model that predicts the critical electric field strength for the jet-to-sheet transition. Our model shows a good agreement with experimental results, and has a form that is consistent with the classical drop deformation criterion in the Taylor-Melcher theory. Finally, we statistically analyze the resultant droplets from sheet breakup, and find that increasing the applied electric field strength improves droplet uniformity and reduces droplet size.

Description: Normal impingement of a single droplet on a thin liquid film is investigated numerically solving the axisymmetric Navier-Stokes equations. Gravity and viscosity are taken into account whereas compressibility effects are neglected. Two phases are tracked by means of volume of fluid method and adaptive mesh refinement is used to increase accuracy of the interface. Numerical results are validated both qualitatively and quantitatively using experimental measurements. Effects of gas density, gas viscosity, and film thickness on the crown behavior are studied. Influence of droplet deviation from spherical shape on the crown behavior is investigated. It is shown that increasing the gas density leads to reduction of crown radius evolution rate, while gas viscosity does not affect the rate of crown radius evolution. Development rate of crown height decreases by increasing the gas density. Reynolds number and splashing regime can change the effect of gas viscosity on the crown height evolution. Deviation of droplet from sphere can change behavior of crown completely as result of change in droplet mass center position. Difference between numerical results and experimental ones is justified using different droplet shapes.

Description: We present an analytical study of peak mode isotachophoresis (ITP), and provide closed form solutions for sample distribution and electric field, as well as for leading-, trailing-, and counter-ion concentration profiles. Importantly, the solution we present is valid not only for the case of fully ionized species, but also for systems of weak electrolytes which better represent real buffer systems and for multivalent analytes such as proteins and DNA. The model reveals two major scales which govern the electric field and buffer distributions, and an additional length scale governing analyte distribution. Using well-controlled experiments, and numerical simulations, we verify and validate the model and highlight its key merits as well as its limitations. We demonstrate the use of the model for determining the peak concentration of focused sample based on known buffer and analyte properties, and show it differs significantly from commonly used approximations based on the interface width alone. We further apply our model for studying reactions between multiple species having different effective mobilities yet co-focused at a single ITP interface. We find a closed form expression for an effective-on rate which depends on reactants distributions, and derive the conditions for optimizing such reactions. Interestingly, the model reveals that maximum reaction rate is not necessarily obtained when the concentration profiles of the reacting species perfectly overlap. In addition to the exact solutions, we derive throughout several closed form engineering approximations which are based on elementary functions and are simple to implement, yet maintain the interplay between the important scales. Both the exact and approximate solutions provide insight into sample focusing and can be used to design and optimize ITP-based assays.

Description: A novel waveform modified from the standard-sinusoidal function is adopted to enhance the virtual aeroshaping effect of the synthetic jets positioned at the front stagnation point of a circular cylinder. The waveform is characterized by a control parameter, namely, the suction duty cycle factor k , which is the ratio of the time duration of the suction cycle to that of the blowing cycle. The strength of the synthetic jet vortex pair is enhanced by increasing the suction duty cycle factor. The periodic closed envelope forms upstream of the circular cylinder for k ≤ 1.00, while the quasi-steady open envelope forms for k ≥ 2.00, acting the virtual aeroshaping effect. As a result, both the statistical characteristics and the vortex dynamics of the near-wake flow field change with the suction duty cycle factor. The recirculation region downstream of the circular cylinder becomes smaller or even disappears, and thus, the drag coefficient over the circular cylinder is reduced by increasing the suction duty cycle factor to k ≥ 1.00. The statistical mean and fluctuating velocities show corresponding changes in the near wake with the different wake patterns. For k ≤ 0.50, the wake vortex shows the antisymmetric shedding mode which is similar with the natural case. For 1.00 ≤ k ≤ 2.00, the wake vortex shows the bistable state mode, where vortex sheds with symmetric or antisymmetric mode; the antisymmetric shedding mode dominates the global flow field for k = 1.00, while it is the symmetric shedding mode that dominates the flow field for k = 2.00. For k = 4.00, it shows the antisymmetric shedding mode with a shorter vortex formation length than the natural case. The above findings indicate that the virtual aeroshaping effect of the synthetic jets can be enhanced by increasing the suction duty cycle factor so as to increase the momentum coefficient while keeping other control parameters unchanged, providing us another way for effective flow control.

Description: The dynamics of vapor-liquid interface are important because interfacial instability determines bubble growth, detachment frequency, waiting time, shape of bubbles, and the interrelationship between bubble formation sites. In this study, a detailed numerical simulation has been performed to understand the transition in bubble release pattern and multimode bubble formation in saturated pool boiling. The interfaces drop down alternatively at the nodes and antinodes of the wavelengths dictated by Rayleigh-Taylor instability and Taylor-Helmholtz instability. Due to higher degrees of superheat, vapor jets emanate from nodes and antinodes. An attempt has been made to predict the maximum and minimum heat fluxes during saturated pool boiling.

Description: We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech.662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.

Description: Janus droplets are compound droplets that consist of two adhering drops of different fluids that are suspended in a third fluid. We use the Shan-Chen lattice Boltzmann method for multicomponent mixtures to simulate Janus droplets at rest and in shear. In this simulation model, interfacial tensions are not known a priori from the model parameters and must be determined using numerical experiments. We show that interfacial tensions obtained with the Young-Laplace law are consistent with those measured from the equilibrium geometry. The regimes of adhering, separated, and engulfing droplets were explored. Two different adhesion geometries were considered for two-dimensional simulations of Janus droplets in shear. The first geometry resembles two adhering circles with small overlap. In the second geometry, the two halves are semicircular. For both geometries, the rotation rate of the droplet depends on its orientation. The width of the periodic simulation domain also affects the rotation rate of both droplet types up to an aspect ratio of 6:1 (width:height). While the droplets with the first geometry oscillated about the middle of the domain, the droplets of the second geometry did not translate while rotating. A four-pole vortex structure inside droplets of the second geometry was found. These simulations of single Janus droplets reveal complex behaviour that implies a rich range of possibilities for the rheology of Janus emulsions.

Description: A numerical model of the northern California Current System along the coasts of Washington and British Columbia is used to quantify the impact of submarine canyons on upwelling from the continental slope onto the shelf. Comparisons with an extensive set of observations show that the model adequately represents the seasonal development of near-bottom density, as well as along-shelf currents that are critical in governing shelf-slope exchange. Additional model runs with simplified coastlines and bathymetry are used to isolate the effects of submarine canyons. Near submarine canyons, equatorward flow over the outer shelf is correlated with dense water at canyon heads and subsequent formation of closed cyclonic eddies, which are both associated with cross-shelf ageostrophic forces. Lagrangian particles tracked from the slope to mid shelf show that canyons are associated with upwelling from depths of ~140-260 m. Source depths for upwelling are shallower than 150 m at locations away from canyons and in a model run with bathymetry that is uniform in the along-shelf direction. Water upwelled through canyons is more likely to be found near the bottom over the shelf. Onshore fluxes of relatively saline water through submarine canyons are large enough to increase volume-averaged salinity over the shelf by 0.1--0.2 psu during the early part of the upwelling season. The nitrate input from the slope to the Washington shelf associated with canyons is estimated to be 30–60% of that upwelled to the euphotic zone by local wind-driven upwelling.

Description: The flow north of warm subtropical water though the northeastern Atlantic is known to have many pathways that vary over time. Here we use a combination of upper ocean current measurements between Greenland and Scotland near 60°N and satellite altimetry to examine the space-time variability of poleward transport. The high-resolution scans of currents in the top 400 m show that the Reykjanes Ridge serves as a very effective separator of flow towards the Nordic and Labrador Seas, respectively. Whereas the Labrador Sea branch exhibits two mean flows to the north on the western slope of the Reykjanes Ridge, the eastern branch flows north in roughly equal amounts over the deep Maury channel and east of Hatton Bank including the Slope Current. There is also a well-defined southward flow along the eastern slope of the Reykjanes Ridge. The satellite altimetric sea surface height (SSH) data show good overall agreement with geostrophically determined sea level difference from the repeat ADCP sections (1999-2002), but are unable to resolve the fine structure of the topographically defined mean circulation. The altimetric data show that variations in poleward flow west and east of the Reykjanes Ridge are strongly anticorrelated. They further reveal that the two eastern sub-branches also exhibit anticorrelated variability, but offset in time with respect to the Labrador Sea branch. Remarkably, all these variations cancel out for the entire Greenland-Scotland section leaving a gradual decrease in sea level difference of about 0.06 m over the 1993 to the end of 2010 observation period.

Description: Understanding the physics of water evaporation from saline porous media is important in many natural and engineering applications such as durability of building materials and preservation of monuments, water quality, and mineral-fluid interactions. We applied synchrotron x-ray micro-tomography to investigate the pore-scale dynamics of dissolved salt distribution in a three dimensional drying saline porous media using a cylindrical plastic column (15 mm in height and 8 mm in diameter) packed with sand particles saturated with CaI 2 solution (5% concentration by mass) with a spatial and temporal resolution of 12 μ m and 30 min, respectively. Every time the drying sand column was set to be imaged, two different images were recorded using distinct synchrotron x-rays energies immediately above and below the K-edge value of Iodine. Taking the difference between pixel gray values enabled us to delineate the spatial and temporal distribution of CaI 2 concentration at pore scale. Results indicate that during early stages of evaporation, air preferentially invades large pores at the surface while finer pores remain saturated and connected to the wet zone at bottom via capillary-induced liquid flow acting as evaporating spots. Consequently, the salt concentration increases preferentially in finer pores where evaporation occurs. Higher salt concentration was observed close to the evaporating surface indicating a convection-driven process. The obtained salt profiles were used to evaluate the numerical solution of the convection-diffusion equation (CDE). Results show that the macro-scale CDE could capture the overall trend of the measured salt profiles but fail to produce the exact slope of the profiles. Our results shed new insight on the physics of salt transport and its complex dynamics in drying porous media and establish synchrotron x-ray tomography as an effective tool to investigate the dynamics of salt transport in porous media at high spatial and temporal resolution.

Description: The behaviour of low Reynolds number, non-Boussinesq fountains from four different nozzle geometries (one circular and three rectangular nozzles) are studied. High speed laser schlieren imaging is used to study the fountain behaviour (frequency and penetration height). Bi-orthogonal decomposition and dynamic mode decomposition (DMD) are used to understand the unsteady characteristics of fountains. The flow regimes of fountains are classified as steady, flapping, and flapping-bobbing type. The DMD technique successfully separates the bobbing oscillation from the combined flapping-bobbing oscillation of the fountain. The frequency of the flapping oscillation, and the frequency of the bobbing oscillation in the flapping-bobbing regime scales as St h Fr h = C 1 and S t h F r h 2 = C 2 , respectively, where the characteristic length scale is the smallest dimension ( h ) of the nozzle. The mean steady state penetration heights ( Z s / h ) of “forced” low Reynolds number non-Boussinesq fountains are independent of nozzle shape (circular and rectangular), and scales linearly with the Froude number.

Description: The large-scale properties of self-similar unstably stratified homogeneous (USH) turbulence are investigated using an eddy-damped quasi-normal markovianized approximation of the nonlinear term. This analysis shows that a special role is played by the wave vectors contained in the equatorial plane, i.e., the plane perpendicular to gravity. It is indeed in this plane that turbulent spectra reach their maxima and evolve linearly from their initial condition when their initial infrared exponent is smaller than 4. At other angles, this property is not satisfied and turbulent spectra eventually undergo an evolution dominated by nonlinear backscattering processes. The self-similar evolution of USH turbulence is also shown to be related to the properties of large scales. In particular, the asymptotic growth rate of the mixing length depends on the initial infrared exponent in the equatorial plane. Besides, the self-similar asymptotic values of the concentration and velocity correlations also depend on the properties of large scales. This allows to derive relations between the correlations and the growth rate parameter.

Description: Secondary flow cells are commonly observed in straight laboratory channels, where they are often associated with duct corners. Here, we present velocity measurements acquired with an acoustic Doppler current profiler in a straight reach of the Seine river (France). We show that a remarkably regular series of stationary flow cells spans across the entire channel. They are arranged in pairs of counter-rotating vortices aligned with the primary flow. Their existence away from the river banks contradicts the usual interpretation of these secondary flow structures, which invokes the influence of boundaries. Based on these measurements, we use a depth-averaged model to evaluate the momentum transfer by these structures, and find that it is comparable with the classical turbulent transfer.

Description: Bubble size distributions ranging from 0.5 to 125 μm radius were measured optically during high winds of 13 m s -1 and large-scale wave breaking as part of the Southern Ocean Gas Exchange Experiment. Very small bubbles with radii less than 60 µm were measured at 6-9 m depth using optical measurements of the near forward volume scattering function and critical scattering angle for bubbles (˜80°). The bubble size distributions generally followed a power law distribution with mean slope values ranging from 3.6 to 4.6. The steeper slopes measured here were consistent with what would be expected near the base of the bubble plume. Bubbles, likely stabilized with organic coatings, were present for time periods on the order of 10-100 s at depths of 6-9 m. Here, relatively young seas, with an inverse wave age of approximately 0.88 and shorter characteristic wave scales, produced lower bubble concentrations, shallower bubble penetration depths, and steep bubble size distribution slopes. Conversely, older seas, with an inverse wave age of 0.70 and longer characteristic wave scales, produced relatively higher bubble concentrations penetrating to 15 m depth, larger bubble sizes and shallower bubble size distribution slopes. When extrapolated to 4 m depth using a previously published bubble size distribution, our estimates suggest that the deeply-penetrating small bubbles measured in the Southern Ocean supplied ˜36% of the total void fraction and likely contributed to the transfer and supersaturation of low solubility gases.

Description: Analysis of fluxes across the turbulent/non-turbulent interface (TNTI) of turbulent boundary layers is performed using data from two-dimensional particle image velocimetry (PIV) obtained at high Reynolds numbers. The interface is identified with an iso-surface of kinetic energy, and the rate of change of total kinetic energy ( K ) inside a control volume with the TNTI as a bounding surface is investigated. Features of the growth of the turbulent region into the non-turbulent region by molecular diffusion of K , viscous nibbling, are examined in detail, focussing on correlations between interface orientation, viscous stress tensor elements, and local fluid velocity. At the level of the ensemble (Reynolds) averaged Navier-Stokes equations (RANS), the total kinetic energy K is shown to evolve predominantly due to the turbulent advective fluxes occurring through an average surface which differs considerably from the local, corrugated, sharp interface. The analysis is generalized to a hierarchy of length-scales by spatial filtering of the data as used commonly in Large-Eddy-Simulation (LES) analysis. For the same overall entrainment rate of total kinetic energy, the theoretical analysis shows that the sum of resolved viscous and subgrid-scale advective flux must be independent of scale. Within the experimental limitations of the PIV data, the results agree with these trends, namely that as the filter scale increases, the viscous resolved fluxes decrease while the subgrid-scale advective fluxes increase and tend towards the RANS values at large filter sizes. However, a definitive conclusion can only be made with fully resolved three-dimensional data, over and beyond the large dynamic spatial range presented here. The qualitative trends from the measurement results provide evidence that large-scale transport due to the energy-containing eddies determines the overall rate of entrainment, while viscous effects at the smallest scales provide the physical mechanism ultimately responsible for entrainment. Data spanning over a decade in Reynolds number suggest that the fluxes (or the entrainment velocity) scale with the friction velocity (or equivalently the local turbulent fluctuating velocity), whereas Taylor microscale and boundary-layer thickness are the appropriate length scales at small and large filter sizes, respectively.

Description: In a recent direct numerical simulation (DNS) study [P. K. Yeung and K. R. Sreenivasan, “ Spectrum of passive scalars of high molecular diffusivity in turbulent mixing,” J. Fluid Mech.716, R14 (2013)] with Schmidt number as low as 1/2048, we verified the essential physical content of the theory of Batchelor, Howells, and Townsend [“Small-scale variation of convected quantities like temperature in turbulent fluid. 2. The case of large conductivity,” J. Fluid Mech.5, 134 (1959)] for turbulent passive scalar fields with very strong diffusivity, decaying in the absence of any production mechanism. In particular, we confirmed the existence of the −17/3 power of the scalar spectral density in the so-called inertial-diffusive range. In the present paper, we consider the DNS of the same problem, but in the presence of a uniform mean gradient, which leads to the production of scalar fluctuations at (primarily) the large scales. For the parameters of the simulations, the presence of the mean gradient alters the physics of mixing fundamentally at low Peclet numbers. While the spectrum still follows a −17/3 power law in the inertial-diffusive range, the pre-factor is non-universal and depends on the magnitude of the mean scalar gradient. Spectral transfer is greatly reduced in comparison with those for moderately and weakly diffusive scalars, leading to several distinctive features such as the absence of dissipative anomaly and a new balance of terms in the spectral transfer equation for the scalar variance, differing from the case of zero gradient. We use the DNS results to present an alternative explanation for the observed scaling behavior, and discuss a few spectral characteristics in detail.

Description: Wave interferometry is a remote sensing technique, which is increasingly employed in helioseismology, seismology, and acoustics to retrieve parameters of the propagation medium from two-point cross-correlation functions of random wave fields. Here, we apply interferometry to yearlong records of seafloor pressure at 28 locations off New Zealand's South Island to investigate propagation and directivity properties of infragravity waves away from shore. A compressed cross-correlation function technique is proposed to make the interferometry of dispersive waves more robust, decrease the necessary noise averaging time, and simplify retrieval of quantitative information from noise cross-correlations. The emergence of deterministic wave arrivals from cross-correlations of random wave fields is observed up to the maximum range of 692 km between the pressure sensors in the array. Free, linear waves with a strongly anisotropic distribution of power flux density are found to be dominant in the infragravity wave field. Lowest-frequency components of the infragravity wave field are largely isotropic. The anisotropy has its maximum in the middle of the spectral band and decreases at the high-frequency end of the spectrum. Highest anisotropy peaks correspond to waves coming from portions of the New Zealand's shoreline. Significant contributions are also observed from waves propagating along the coastline and probably coming from powerful sources in the northeast Pacific. Infragravity wave directivity is markedly different to the east and to the west of the South Island. The northwest coast of the South Island is found to be a net source of the infragravity wave energy.

Description: The Indonesian Throughflow (ITF) is the only open pathway for inter-ocean exchange between the Pacific and Indian Ocean basins at tropical latitudes. A proxy time series of ITF transport variability is developed using remotely-sensed altimeter data. The focus is on the three outflow passages of Lombok, Ombai and Timor that collectively transport the entire ITF into the Indian Ocean, and where direct velocity measurements are available to help ground-truth the transport algorithm. The resulting 18-year proxy time series shows strong interannual ITF variability. Significant trends of increased transport are found in the upper layer of Lombok Strait, and over the full depth in Timor Passage that are likely related to enhanced Pacific trade winds since the early 1990s. The partitioning of the total ITF transport through each of the major outflow passage varies according to the phase of the Indian Ocean Dipole (IOD) or El Niño-Southern Oscillation (ENSO). In general, Pacific ENSO variability is strongest in Timor Passage, most likely through the influence of planetary waves transmitted from the Pacific along the Northwest Australian shelf pathway. Somewhat surprisingly, concurrent El Niño and positive IOD episodes consistently show contradictory results from those composites constructed for purely El Niño episodes. This is particularly evident in Lombok and Ombai Straits, but also at depth in Timor Passage. This suggests that Indian Ocean dynamics likely win out over Pacific Ocean dynamics in gating the transport through the outflow passages during concurrent ENSO and IOD events.

Description: Time series of ice draft from 2003-2012 from moored sonar data are used to investigate variability and describe the reduction of the perennial sea ice cover in the Beaufort Gyre (BG), culminating in the extreme minimum in 2012. Negative trends in median ice drafts and most ice fractions are observed, while open water and thinnest ice fractions (〈0.3 m) have increased, attesting to the ablation or removal of the older sea ice from the BG over the nine year period. Monthly anomalies indicate a shift occurred toward thinner ice after 2007, in which the thicker ice evident at the northern stations was reduced. Differences in the ice characteristics between all of the stations also diminished, so that the ice cover throughout the region became statistically homogenous. The moored data are used in a relationship with satellite radiometer data to estimate ice volume changes throughout the BG. Summer solid fresh water content decreased drastically in consecutive years from 730 km 3 in 2006 to 570 km 3 in 2007, and to 240 km 3 in 2008. After a short rebound, solid fresh water fell below 220 km 3 in 2012. Meanwhile hydrographic data indicate that liquid fresh water in the BG in summer increased 5410 km 3 from 2003 to 2010 and decreased at least 210 km 3 by 2012. The reduction of both solid and liquid fresh water components indicates a net export of approximately 320 km 3 of fresh water from the region occurred between 2010 and 2012, suggesting that the anticyclonic atmosphere-ocean circulation has weakened.

Description: The effects of water level variations on breaking wave setup over fringing reefs are assessed using field measurements obtained at three study sites in the Republic of the Marshall Islands and the Mariana Islands in the western tropical Pacific Ocean. At each site, reef flat setup varies over the tidal range with weaker setup at high tide and stronger setup at low tide for a given incident wave height. The observed water level dependence is interpreted in the context of radiation stress gradients specified by an idealized point break model generalized for non-normally incident waves. The tidally varying setup is due in part to depth-limited wave heights on the reef flat, as anticipated from previous reef studies, but also to tidally dependent breaking on the reef face. The tidal dependence of the breaking is interpreted in the context of the point break model in terms of a tidally varying wave height to water depth ratio at breaking. Implications for predictions of wave-driven setup at reef-fringed island shorelines are discussed.

Description: The distribution and thickness of sea ice in the Arctic is changing rapidly, resulting in changes to Arctic marine ecosystems. Seabirds are widely regarded as indicators of marine environmental change, and understanding their distribution patterns can serve as a tool to monitor and elucidate biological changes in the Arctic seas. We examined the at-sea distribution of seabirds in the North American Arctic in July and August, 2007-2012, and marine areas of high density were identified based on bird densities for four foraging guilds. Short-tailed shearwaters ( Puffinus tenuirostris ) were the most abundant species observed. Northern fulmars ( Fulmarus glacialis ), thick-billed murres ( Uria lomvia ) and dovekies ( Alle alle ) were also sighted in large numbers. Few birds were sighted between Dolphin and Union Strait and King William Island. Areas of high density over multiple years were found throughout the entire western portion of the study area (Bering Sea, Bering Strait, Chukchi Sea), Lancaster Sound, Baffin Bay, Davis Strait and the low-arctic waters off Newfoundland. These waters are characterized by high primary productivity. This study is the first to document the marine distribution of seabirds across the entire North American Arctic within the same time period, providing a critical baseline for monitoring the distribution and abundance of Arctic seabirds in a changing Arctic seascape.

Description: A non-equilibrium wall-model based on unsteady 3D Reynolds-averaged Navier-Stokes (RANS) equations has been implemented in an unstructured mesh environment. The method is similar to that of the wall-model for structured mesh described by Wang and Moin [Phys. Fluids14, 2043–2051 (2002)], but is supplemented by a new dynamic eddy viscosity/conductivity model that corrects the effect of the resolved Reynolds stress (resolved turbulent heat flux) on the skin friction (wall heat flux). This correction is crucial in predicting the correct level of the skin friction. Unlike earlier models, this eddy viscosity/conductivity model does not have a stress-matching procedure or a tunable free parameter, and it shows consistent performance over a wide range of Reynolds numbers. The wall-model is validated against canonical (attached) transitional and fully turbulent flows at moderate to very high Reynolds numbers: a turbulent channel flow at Re τ = 2000, an H-type transitional boundary layer up to Re θ = 3300, and a high Reynolds number boundary layer at Re θ = 31 000. Application to a separated flow over a NACA4412 airfoil operating close to maximum lift is also considered to test the performance of the wall-model in complex non-equilibrium flows.

Description: Conventional shallow water theory successfully reproduces many key features of the Jovian atmosphere: a mixture of coherent vortices and stable, large-scale, zonal jets whose amplitude decreases with distance from the equator. However, both freely decaying and forced-dissipative simulations of the shallow water equations in Jovian parameter regimes invariably yield retrograde equatorial jets, while Jupiter itself has a strong prograde equatorial jet. Simulations by Scott and Polvani [“Equatorial superrotation in shallow atmospheres,” Geophys. Res. Lett.35, L24202 (2008)] have produced prograde equatorial jets through the addition of a model for radiative relaxation in the shallow water height equation. However, their model does not conserve mass or momentum in the active layer, and produces mid-latitude jets much weaker than the equatorial jet. We present the thermal shallow water equations as an alternative model for Jovian atmospheres. These equations permit horizontal variations in the thermodynamic properties of the fluid within the active layer. We incorporate a radiative relaxation term in the separate temperature equation, leaving the mass and momentum conservation equations untouched. Simulations of this model in the Jovian regime yield a strong prograde equatorial jet, and larger amplitude mid-latitude jets than the Scott and Polvani model. For both models, the slope of the non-zonal energy spectra is consistent with the classic Kolmogorov scaling, and the slope of the zonal energy spectra is consistent with the much steeper spectrum observed for Jupiter. We also perform simulations of the thermal shallow water equations for Neptunian parameter values, with a radiative relaxation time scale calculated for the same 25 mbar pressure level we used for Jupiter. These Neptunian simulations reproduce the broad, retrograde equatorial jet and prograde mid-latitude jets seen in observations. The much longer radiative time scale for the colder planet Neptune explains the transition from a prograde to a retrograde equatorial jet, while the broader jets are due to the deformation radius being a larger fraction of the planetary radius.

Description: The hydrodynamic interaction of two deformable vesicles in shear flow induces a net displacement, in most cases an increase of their distance in the transverse direction. The statistical average of these interactions leads to shear-induced diffusion in the suspension, both at the level of individual particles which experience a random walk made of successive interactions, and at the level of suspension where a nonlinear down-gradient diffusion takes place, an important ingredient in the structuring of suspension flows. We make an experimental and computational study of the interaction of a pair of lipid vesicles in shear flow by varying physical parameters, and investigate the decay of the net lateral displacement with the distance between the streamlines on which the vesicles are initially located. This decay and its dependency upon vesicle properties can be accounted for by a simple model based on the well established law for the lateral drift of a vesicle in the vicinity of a wall. In the semi-dilute regime, a determination of self-diffusion coefficients is presented.

Description: Radiative transfer in sea ice is subject to anisotropic, multiple scattering. The impact of anisotropy on the light field under sea ice was found to be substantial and has been quantified. In this study, a large dataset of irradiance and radiance measurements under sea ice has been acquired with a Remotely Operated Vehicle (ROV) in the central Arctic. Measurements are interpreted in the context of numerical radiative transfer calculations, laboratory experiments, and microstructure analysis. The ratio of synchronous measurements of transmitted irradiance to radiance shows a clear deviation from an isotropic under-ice light field. We find that the angular radiance distribution under sea-ice is more downward directed than expected for an isotropic light field. This effect can be attributed to the anisotropic scattering coefficient within sea ice. Assuming an isotropic radiance distribution under sea ice leads to significant errors in light-field modeling and the interpretation of radiation measurements. Quantification of the light field geometry is crucial for correct conversion of radiance data acquired by Autonomous Underwater Vehicles (AUVs) and ROVs.

Description: In comparison to the deep ocean, the upper mixed layer is a region typically characterized by substantial vertical gradients in water properties. Within the Tropics, the rich variability in the vertical shapes and forms that these structures can assume through variation in the atmospheric forcing results in a differential effect in terms of the temperature and salinity stratification. Rather than focusing on the strong halocline above the thermocline, commonly referred to as the salinity barrier layer, the present study takes into account the respective thermal and saline dependencies in the Brunt-Väisälä frequency (N 2 ) in order to isolate the specific role of the salinity stratification in the layers above the main pycnocline. We examine daily vertical profiles of temperature and salinity from an ocean reanalysis over the period 2001-2007. We find significant seasonal variations in the Brunt-Väisälä frequency profiles are limited to the upper 300 m depth. Based on this we determine the ocean salinity stratification (OSS) to be defined as the stabilizing effect (positive values) due to the haline part of N 2 averaged over the upper 300m. In many regions of the tropics the OSS contributes 40 to 50% to N 2 as compared to the thermal stratification and, in some specific regions, exceeds it for a few months of the seasonal cycle. Away from the tropics, for example near the centers of action of the subtropical gyres, there are regions characterized by the permanent absence of OSS. In other regions previously characterized with salinity barrier layers the OSS obviously shares some common variations: however, we show that where temperature and salinity are mixed over the same depth, the salinity stratification can be significant. In addition, relationships between the OSS and the sea surface salinity are shown to be well defined and quasi-linear in the tropics, providing some indication that in the future, analyses that consider both satellite surface salinity measurements at the surface and vertical profiles at depth will result in a better determination of the role of the salinity stratification in climate prediction systems.

Description: Impacts of extreme Barents Sea air-sea exchanges are examined using the HadCM3 coupled ocean-atmosphere model. Variability in the Barents Sea winter air-sea density flux is found to be a potentially significant factor in determining changes in the southward transport of dense water through Fram Strait. The density flux variability is primarily driven by the thermal term, F T , due to heat loss to the atmosphere. The other two terms (haline flux and ice formation) play a relatively minor role. The difference in ocean circulation between winters with extreme strong and weak Barents Sea surface density flux anomalies is analysed. This reveals an increase in strong winters of both the north-westwards intermediate depth flow out of the basin and the east-west deep flows north of Spitsbergen and south through the Fram Strait. A linear fit yields a Fram Strait southward transport increase of 1.22 Sv for an increase in F T of 1x10 -6 kg m -2 s -1 . For the ten strongest Barents Sea surface density flux winters, the Fram Strait winter southward transport increases by 2.4 Sv. This compares with a reduction of 1.0 Sv for the corresponding weakest winters. Furthermore, the properties of the southwards flowing water are modified in strong density flux winters. In such winters, the Fram Strait water below 250 m is colder by up to 0.5 ° C and fresher by 0.05 than the climatological winter mean.

Description: In viscous-plastic (VP) sea-ice models, small deformations are approximated by irreversible viscous deformations, introducing a non-physical energy sink. As the spatial resolution and the degree of numerical convergence of the models increase, linear kinematic features (LKFs) are better resolved and more states of stress lie in the viscous regime. Energy dissipation in this non-physical viscous regime therefore increases. We derive a complete kinetic energy (KE) balance for sea ice, including plastic and viscous energy sinks to study energy dissipation. The main KE balance is between the energy input by the wind and the dissipation by the water drag and the internal stresses (dissipating respectively 87% and 13% of the energy input on an annual average). The internal stress term is mostly important in winter when ice-ice interactions are dominant. The energy input that is not dissipated locally is redistributed laterally by the internal stresses in regions of dissipation by small scale deformations (LKFs). Of the 13% dissipated annually by the internal stress term, 93% is dissipated in friction along LKFs (14% in ridging, 79% in shearing) and 7% is stored as potential energy in ridges. For all time and spatial scales tested, the frictional viscous dissipation is negligible in the KE balance. This conclusion remains valid when the spatial resolution and the numerical convergence of the simulations are increased. Overall, the results confirm the validity, from an energetical point of view, of the VP approximation.

Description: The magnetohydrodynamic Richtmyer-Meshkov instability is investigated for the case where the initial magnetic field is unperturbed and aligned with the mean interface location. For this initial condition, the magnetic field lines penetrate the perturbed density interface, forbidding a tangential velocity jump and therefore the presence of a vortex sheet. Through simulation, we find that the vorticity distribution present on the interface immediately after the shock acceleration breaks up into waves traveling parallel and anti-parallel to the magnetic field, which transport the vorticity. The interference of these waves as they propagate causes the perturbation amplitude of the interface to oscillate in time. This interface behavior is accurately predicted over a broad range of parameters by an incompressible linearized model derived presently by solving the corresponding impulse driven, linearized initial value problem. Our use of an equilibrium initial condition results in interface motion produced solely by the impulsive acceleration. Nonlinear compressible simulations are used to investigate the behavior of the transverse field magnetohydrodynamic Richtmyer-Meshkov instability, and the performance of the incompressible model, over a range of shock strengths, magnetic field strengths, perturbation amplitudes and Atwood numbers.

Description: We investigate experimentally the statistical properties of bedload transport induced by a steady, uniform, and laminar flow. We focus chiefly on lateral transport. The analysis is restricted to experiments where the flow-induced shear stress is just above the threshold for sediment transport. We find that, in this regime, the concentration of moving particles is low enough to neglect interactions between themselves. We can therefore represent bedload as a thin layer of independent walkers travelling over the bed surface. In addition to their downstream motion, the particles show significant fluctuations of their cross-stream velocity, likely due to the roughness of the underlying sediment bed. This causes particles to disperse laterally. Based on thousands of individual trajectories, we show that this lateral spreading is the manifestation of a random walk. The experiments are entirely consistent with Fickian diffusion.

Description: We study the behavior of heavy inertial particles in the flow field of two like-signed vortices. In a frame co-rotating with the two vortices, we find that stable fixed points exist for these heavy inertial particles; these stable frame-fixed points exist only for particle Stokes number St 〈 St cr . We estimate St cr and compare this with direct numerical simulations, and find that the addition of viscosity increases the St cr slightly. We find that the rate at which particles fall into the fixed points increases until the fixed points disappear at St = St cr . These frame-fixed points are between fixed points and limit cycles in character.

Description: The paper reports a new phenomenon—vortex flows in isothermal magnetic fluids in the vicinity of the localized source of magnetic field (magnetized iron sphere) induced by the drift of drop-like aggregates. Although the observed magnetic precipitation of drop-like aggregates resembles an ordinary rainfall in the Earth atmosphere, its origin and nature are quite different. In magnetic fluids this “rain” is induced by the non-uniform magnetic field and occurs at the scale of 1 mm, not at the scale of several kilometers as in the Earth atmosphere. The reason of this phenomenon is that the applied magnetic field initiates phase transition of “gas-liquid” type which is accompanied by formation of condensed phase represented by drop-like aggregates with the characteristic dimension of about tens of micrometers elongated along the field lines. Inhomogeneous spatial distribution of drop-like aggregates leads to deviation of the ponderomotive force, which is responsible for the formation of vortex flows in the fluid. The “rain” is the primary reason for the vortex flows and it lasts until all magnetic particles capable of condensing into drop-like aggregates precipitate at the surface of the condensation core (iron sphere). Thus, vortex flows induced by drop-like aggregate magnetophoresis represent one variant of “gas-liquid” phase transition. Hydrodynamic flows intensify mass transfer in vicinity of magnetic condensation core and considerably speed it up.

Description: Based on a long term simulation of an ocean-biogeochemical coupled model, we investigate the biogeochemical response to the two types of El Niño events, a Cold Tongue (CT)-El Niño and a Warm Pool (WP)-El Niño, in which a local maximum of anomalous sea surface temperature (SST) is located in the eastern and central tropical Pacific. Our model is able to reasonably simulate the characteristics of the biological variables in a way comparable to the observations. During the developing period, anomalous low chlorophyll appears in the eastern Pacific, while it appears in the central Pacific in the WP-El Niño. The difference in the spatial-temporal response of chlorophyll for the two types of El Niño events is mainly due to the eastward zonal advection of upper ocean currents, which plays a role in bringing nutrient-poor water from the western Pacific. During the decaying period of the WP-El Niño, anomalous high chlorophyll appears concurrently with anomalous low SST in the eastern Pacific. Conversely, anomalous high chlorophyll appears in the central Pacific prior to the decaying period of the CT-El Niño. In particular, the anomalous low sea level from the northwestern Pacific shifts to the southern equatorial region during the decaying period of the CT-El Niño. This drives anti-cyclonic boundary currents which enhance the Equatorial Undercurrent, playing a role in the supply of nutrients to the central equatorial Pacific, resulting in an increase in chlorophyll concentration in the same region.

Description: The barrier layer, the layer between the bottom of the density-defined mixed layer and the isothermal layer in the upper ocean, may play a role in air-sea dynamics. In the present study, data from Argo profiling floats in the tropical Indian Ocean and a mooring at 90°E, 0°N are used to examine subseasonal variations in upper ocean salinity and barrier-layer thickness (BLT) during boreal winter. In the eastern equatorial Indian Ocean, subseasonal variations in BLT are energetic. However, composites used to isolate the Madden-Julian Oscillation (MJO) component of the subseasonal signal reveal that, on average, the MJO anomaly in BLT is negligible despite large swings in both the mixed-layer depth and the isothermal-layer depth. This discrepancy is likely due to (a) noise from other subseasonal processes; and (b) the diversity of individual MJO events: the thickness of the mixed layer and the isothermal layer are sensitive to wind and rain forcing, so even subtle differences in the phasing and strength of MJO-related atmospheric anomalies can produce a very different effect on upper ocean stratification and hence on the thickness of the barrier layer. The effect of the barrier layer on the upper ocean response to MJO forcing is also evaluated. When the barrier layer is thick, entrainment cooling during the MJO is reduced, so the MJO drives a weaker sea surface temperature anomaly. This suggests that modulation of BLT can have significant consequences for the response of the upper ocean to the MJO, and hence, potentially, for feedbacks of the ocean onto the atmosphere on MJO timescales.

Description: The variability of the Antarctic Circumpolar Current (ACC) system is largely linked to atmospheric forcing. The objective of this work is to assess the link between local wind forcing mechanisms and the variability of the upper-ocean temperature and the dynamics of the different fronts in the ACC region south of South Africa. To accomplish this, in situ and satellite-derived observations are used between 1993-2010. The main finding of this work is that meridional changes in the westerlies linked with the Southern Annular Mode (SAM) drive temperature anomalies in the Ekman layer and changes in the Subantarctic Front (SAF) and Antarctic Polar Front (APF) transports through Ekman dynamics. The development of easterly anomalies between 35 o S-45 o S during positive SAM is linked to reduced (increased) SAF (APF) transports and a warmer mixed layer in the ACC. The link between the changes in the wind stress and the SAF and APF transport variations occurs through the development of Ekman pumping anomalies near the frontal boundaries, driving an opposite response on the SAF and APF transports. The observed wind-driven changes in the frontal transports suggest small changes to the net ACC transport. In addition, observations indicate that the SAF and APF locations in this region are not linked to the local wind forcing, emphasizing the importance of other factors (e.g. baroclinic instabilities generated by bottom topography) to changes in the frontal location. Results obtained here highlight the importance of repeat XBT temperature sections and their combined analysis with other in situ and remote sensing observations.

Description: Compressive sampling is well-known to be a useful tool used to resolve the energetic content of signals that admit a sparse representation. The broadband temporal spectrum acquired from point measurements in wall-bounded turbulence has precluded the prior use of compressive sampling in this kind of flow, however it is shown here that the frequency content of flow fields that have been Fourier transformed in the homogeneous spatial (wall-parallel) directions is approximately sparse, giving rise to a compact representation of the velocity field. As such, compressive sampling is an ideal tool for reducing the amount of information required to approximate the velocity field. Further, success of the compressive sampling approach provides strong evidence that this representation is both physically meaningful and indicative of special properties of wall turbulence. Another advantage of compressive sampling over periodic sampling becomes evident at high Reynolds numbers, since the number of samples required to resolve a given bandwidth with compressive sampling scales as the logarithm of the dynamically significant bandwidth instead of linearly for periodic sampling. The combination of the Fourier decomposition in the wall-parallel directions, the approximate sparsity in frequency, and empirical bounds on the convection velocity leads to a compact representation of an otherwise broadband distribution of energy in the space defined by streamwise and spanwise wavenumber, frequency, and wall-normal location. The data storage requirements for reconstruction of the full field using compressive sampling are shown to be significantly less than for periodic sampling, in which the Nyquist criterion limits the maximum frequency that can be resolved. Conversely, compressive sampling maximizes the frequency range that can be recovered if the number of samples is limited, resolving frequencies up to several times higher than the mean sampling rate. It is proposed that the approximate sparsity in frequency and the corresponding structure in the spatial domain can be exploited to design simulation schemes for canonical wall turbulence with significantly reduced computational expense compared with current techniques.

Description: We investigate the role of deep-ocean topography in scattering energy from the large spatial scales of the low mode internal tide to the smaller spatial scales of higher modes. The complete Green function method, which is not subject to the restrictions of the WKB approximation, is used for the first time to study the two-dimensional scattering of a mode-1 internal tide incident on subcritical and supercritical topography of any form in arbitrary stratifications. For an isolated Gaussian ridge in a uniform stratification, large amplitude critical topography is the most efficient at mode-1 scattering and small amplitude topography scatters with an efficiency on the order of 5-10%. In a nonuniform stratification with a pycnocline, the results are qualitatively the same as for a constant stratification, albeit with the key features shifted to larger height ratios. Having validated these results by direct comparison with the results of nonlinear numerical simulations, and in the process demonstrated that WKB results are not appropriate for reasonable ocean predictions, we proceed to use the Green function approach to quantify the role of topographic scattering for the region of the Pacific ocean surrounding the Hawaiian Islands chain. To the south, the Line Islands ridge is found to scatter ~ 40% of a mode-1 internal tide coming from the Hawaiian Ridge. To the north, realistic, small-amplitude, rough topography scatters ~ 5-10% of the energy out of mode-1 for transects of length 1000-3000km. A significant finding is that compared to large extents of small-amplitude, rough topography a single large topographic feature along the path of a mode-1 internal tide plays the dominant role in scattering the internal tide.

Description: The downslope descent of dense shelf water in the northwestern part of the Sea of Japan is investigated from a dynamical point of view, paying attention to the formation of bottom water in the Winter of 2001. It is supposed that before 1980, the shelf water in Peter the Great Bay sometimes descended far down the continental slope, at least partly reaching depths in excess of 3000 m (the foot of the continental slope). After 1980, however, the shelf water did not descend as far; it either descended only moderately or not at all. In Winter 2001, however, the dense shelf water again descended to depths greater than 3000 m, resulting in the formation of bottom water. Descents of more than 3000 m are due to low temperatures coupled with high salinities, whereas the moderate descents of the late 20th century were purely related to the low temperature of the shelf water. It is estimated that over the continental slope the shelf water becomes mixed with the ambient water in a ratio of about 1:9 for deep descents, whereas the ratio is about 1:5 for moderate descents. The formation of bottom water is greatly influenced by interannual atmospheric variability; thus in Winter 2001 a combination of the strengthened Siberian High (especially in its northern part) and the Aleutian Low advected very cold air into northeast Asia, producing dense shelf water and resulting in the formation of bottom water.

Description: The transports of water, heat and salt between the northwestern shelf and deep interior of the Black Sea are investigated using a high-resolution three-dimensional primitive equation model. From April to August, 2005, both onshore and offshore cross-shelf break transports in the top 20 m were 0.24 Sv on average, which is equivalent to the replacement of 60% of the volume of surface shelf waters (0 – 20 m) per month. Two main exchange mechanisms are studied: Ekman transport, and transport by mesoscale eddies and associated meanders of the Rim Current. The Ekman drift causes nearly uniform onshore or offshore flow over a large section of the shelf break, but it is confined to the upper layers. In contrast, eddies and meanders penetrate deep down to the bottom, but they are restricted laterally. During the strong wind events of April 15 – 22 and July 1 – 4, some 0.66×10 12 and 0.44×10 12 m 3 of water were removed from the northwestern shelf respectively. In comparison, the single long-lived Sevastopol Eddy generated a much larger offshore transfer of 2.84×10 12 m 3 over the period April 23 to June 30, which is equivalent to 102% of the volume of northwestern shelf waters. Over the study period, salt exchanges increased the average density of the shelf waters by 0.67 kg m -3 and reduced the density contrast between the shelf and deep sea, while lateral heat exchanges reduced the density of the shelf waters by 0.16 kg m -3 and sharpened the shelf break front.

Description: This paper presents extensive acoustic measurements on jets impinging on surfaces of various surface roughness values. Besides surface roughness, the effects of nozzle-to-plate spacing distance and nozzle pressure ratio are also investigated. Turbulent mixing noise and tonal noise are explained using far-field wall-jet velocity and impingement region temperature fields. The results demonstrate that roughness of the impingement plate widens the staging region of impingement noise. In general, high speed jet impinging on a rough plate generates less noise compared to a smooth plate. When tones are removed from the spectra, it is found that acoustic power monotonically decreases with increasing surface roughness. Thermal imaging in the stagnation region indicates that whenever tones are present, the temperature at the stagnation region is high. Further, sound pressure directivity pattern of impingement noise is constructed by superimposing a wall-jet and a free jet in the appropriate orientations.

Description: We investigate the elasticity of an isolated, threefold junction of soap films (Plateau border), which displays static undulations when liquid rapidly flows into it. By analyzing the shape of the Plateau border (thickness R and transverse displacement) as a function of the liquid flow rate Q , we show experimentally and theoretically that the elasticity of the Plateau border is dominated by the bending of the soap films pulling on the Plateau border. In this asymptotic regime, the undulation wavelength obeys the scaling law ∼ Q 2   R −2 and the decay length ∼ Q 2   R −4 .

Description: The case of a turbulent round jet impinging perpendicularly onto a rotating, heated disc is investigated, in order to understand the mechanisms at the origin of the influence of rotation on the radial wall jet and associated heat transfer. The present study is based on the complementary use of an analysis of the orders of magnitude of the terms of the mean momentum and Reynolds stress transport equations, available experiments, and dedicated Reynolds-averaged Navier–Stokes computations with refined turbulence models. The Reynolds number Re j = 14 500, the orifice-to-plate distance H = 5 D , where D is the jet-orifice diameter, and the four rotation rates were chosen to match the experiments of Minagawa and Obi [“Development of turbulent impinging jet on a rotating disk,” Int. J. Heat Fluid Flow25, 759–766 (2004)] and comparisons are made with the Nusselt number distribution measured by Popiel and Boguslawski [“Local heat transfer from a rotating disk in an impinging round jet,” J. Heat Transfer108, 357–364 (1986)], at a higher Reynolds number. The overestimation of turbulent mixing in the free-jet before the impact on the disk is detrimental to the prediction of the impingement region, in particular of the Nusselt number close to the symmetry axis, but the self-similar wall jet developing along the disk is correctly reproduced by the models. The analysis, experiments, and computations show that the rotational effect do not directly affect the outer layer, but only the inner layer of the wall jet. A noteworthy consequence is that entrainment at the outer edge of the wall jet is insensitive to rotation, which explains the dependence of the wall-jet thickness on the inverse of the non-dimensional rotation rate, observed in the experiments and the Reynolds stress model computations, but not reproduced by the eddy-viscosity models, due to the algebraic dependence to the mean flow. The analysis makes moreover possible the identification of a scenario for the appearance of rotational effects when the rotation rate is gradually increased. For weak rotation rates, the rotation-induced boundary layer appears but does not break the self-similar solution observed for the case without rotation. For intermediate rotation rates, the production of the azimuthal Reynolds stress becomes much stronger than other components, leading to a complete modification of the turbulence anisotropy which is reproduced only by Reynolds stress models. For strong rotation rates, centrifugal effects dominate, leading to an acceleration and thinning of the jet, and consequently an increase of turbulent production and heat transfer, reproduced by all the turbulence models.

Description: Microfluidic channels are powerful means of control of minute volumes such as droplets. These droplets are usually conveyed at will in an externally imposed flow which follows the geometry of the micro-channel. It has recently been pointed out by Dangla et al. [“Trapping microfluidic drops in wells of surface energy,” Phys. Rev. Lett.107(12), 124501 (2011)] that the motion of transported droplets may also be stopped in the flow, when they are anchored to grooves which are etched in the channels top wall. This feature of the channel geometry explores a direction that is usually uniform in microfluidics. Herein, this anchoring effect exploiting the three spatial directions is studied combining a depth averaged fluid description and a geometrical model that accounts for the shape of the droplet in the anchor. First, the presented method is shown to enable the capture and release droplets in numerical simulations. Second, this tool is used in a numerical investigation of the physical mechanisms at play in the capture of the droplet: a localized reduced Laplace pressure jump is found on its interface when the droplet penetrates the groove. This modified boundary condition helps the droplet cope with the linear pressure drop in the surrounding fluid. Held on the anchor the droplet deforms and stretches in the flow. The combination of these ingredients leads to recover the scaling law for the critical capillary number at which the droplets exit the anchors C a ★ ∝ h 2 / R 2 where h is the channel height and R the droplet undeformed radius.

Description: We present a numerical study of nanosecond pulsed dielectric barrier discharge (DBD) actuator operating in quiescent air at atmospheric condition. Our study concentrates on plasma discharge induced fluid dynamics and on exploration of parametric space of interest for voltage pulse in an attempt to shed some light into elucidation of the mechanisms whereby the generated shock wave propagates through and affects the external flow. Specifically, a one-dimensional, self-similar, local ionization kinetic model recently developed to predict key parameters of nanosecond pulsed plasma discharge is coupled with the compressible Navier-Stokes equations possibly for the first time. Within the considered range of parameters of the plasma model which is justified for the modeling of surface nanosecond pulsed discharge at atmospheric pressure, our coupled method is able to provide satisfactory prediction of the shock structure generated by the actuator for comparison with experiment, not only in the qualitative shock wave shape but also in quantitative shock front displacement. We provide a comprehensive analysis of the gas heating, shock wave initiation and evolution processes. For example, the characteristic time of the rapid localized heating responsible for shock wave generation, which is yet to be quantified experimentally, is found to be ∼350 ns. We conduct a parametric investigation by varying the peak voltage from 10 kV to 50 kV and rise time from 5 ns to 150 ns. The pressure wave whose behavior is found to be dominated by input voltage amplitude, introduces highly transient, localized disturbance to the quiescent air. In addition, the vortex induced by the shock passage is relatively weak. The interplay of the induced flows by a few successive plasma discharges operating at continuous mode does not appear to be significant, especially at low voltage amplitude.

Description: In this paper, we investigate the decay of incompressible homogeneous isotropic turbulence in a variable viscosity fluid. The viscosity coefficient is assumed to depend linearly on a scalar, representing either a temperature or a concentration, and obeying a simple advection-diffusion equation. At high Reynolds numbers, Direct Numerical Simulations (DNS) allow us to confirm the validity of Taylor's postulate that the dissipation is independent from the viscosity and its fluctuations. At low Reynolds numbers, we report the presence of extra energy at small scales due to these variable viscosity effects. This implies that the turbulent kinetic energy decreases less rapidly as a function of time in variable viscosity fluids. In order to explain this phenomenon and quantify its importance on the turbulent flow, we propose a statistical approach based on an eddy-damped quasi-normal Markovian (EDQNM) spectral closure which takes into account the nonlinearity introduced by variable viscosity. It is shown that this latter additional term is of constant sign in the energy spectrum equation and reduces the dissipation of the flow as observed. Also, by assuming the dominance of distant interactions between wave numbers, we can propose a simple formula expressing that variable viscosity effects lead to an effective reduction of the mean viscosity proportional to the variance of viscosity fluctuations.

Description: Liquid droplets flowing through a rectangular microfluidic channel develop a vortical flow field due to the presence of shear forces from the surrounding fluid. In this paper, we present an experimental and computational study of droplet velocities and internal flow patterns in a rectangular pressure-driven flow for droplet diameters ranging from 0.1 to 2 times the channel height. Our study shows excellent agreement with asymptotic predictions of droplet and interfacial velocities for infinitesimally small droplets. As the droplet diameter nears the size of the channel height, the droplet velocity slows significantly, and the changing external flow field causes a qualitative change in the location of internal vortices. This behavior is relevant for future studies of mass transfer in microfluidic devices.

Description: Laboratory experiments are performed to examine the formation of a crater in sediment by an impinging vertical turbulent jet. Light attenuation and a “depositometer,” which records conductivity through the bed from an array of electrodes, are used to measure the crater depth as a function of space and time. The onset of crater formation and deepening is best characterized in terms of the Rouse number, Rs (proportional to the particle settling speed divided by the centerline jet speed), rather than Shields number, Sh (proportional to the stress divided by the particle weight per unit area). The critical Rouse number, Rs c , is found to increase with the particle Reynolds number, Re p , as a power law with exponent 0.45 ± 0.03 for Re p ranging between 0.6 and 160. For smaller Rs, the crater is observed to deepen at a near-constant speed, while the crater radius remains constant. Bedload transport, measured in terms of the crater deepening speed, is determined to increase as Re p times the difference between Rs c and Rs.

Description: We report the experimental studies of the statistical and scaling properties of the fully developed turbulent regime in von Karman swirling flow between counter-rotating disks with and without blades using the only global measurements of the spatially averaged torque Γ and pressure p fluctuations in water and water-sugar solutions of different viscosities in the same cell geometry. We show that for all fluids under investigation probability distribution functions (PDFs) of the torque fluctuations δΓ/Γ rms are Gaussian in both the laminar and turbulent regimes and for the both types of the stirrers. On the contrary, PDFs of the pressure fluctuations change from Gaussian in the laminar regime into the skewed shape with the exponential tails toward low-pressure events for both the entrainment methods. Both the friction coefficient C f and normalized rms of the pressure fluctuations C p are independent of Re in the fully developed turbulent regime for all fluids under study and found in a good quantitative agreement with the previous results. We also observe that the internal flow variables such as the normalized torque Γ ¯ / V p r m s versus the “internal” Reynolds number Re rms = ( p rms /ρ) 1/2 R ρ/η instead of the global variables C f , C p versus Re show sharp transition into the well developed turbulent regime. We find that the scaling exponents of the fundamental characteristics based only on Γ and p measurements in the range of fully developed turbulent flow, namely, the integral, Taylor, and Kolmogorov dissipation lengths, as well as the Taylor-based Reynolds number R λ , are in rather fair agreement with the predictions. We would like to emphasize that scaling of the main turbulent parameters R λ , λ, η d obtained via the global variables is a very non-trivial result. It is not obvious that measurements based on the global quantities will provide the predicted scaling relations. The result on such scaling obtained previously strongly disagrees with the scaling predictions. Indeed, both Γ ¯ and p rms are averaged over the cell volume as well as all spatial scales, whereas the swirling flow is neither isotropic nor homogeneous. So the global variables being averaged over all spatial scales get contributions from the scales larger and smaller than those from the inertial range of scales. And finally, the normalized characteristic frequencies f p / f rot found in both the torque and pressure frequency power spectra in the fully developed turbulent regime have close values, are independent of Re , and associated with either the rotation or oscillation frequency of the main vortex.

Description: Numerical computations are presented to study the effect of soluble surfactant on the deformation and breakup of an axisymmetric drop or bubble stretched by an imposed linear strain flow in a viscous fluid. At the high values of bulk Peclet number Pe in typical fluid-surfactant systems, there is a thin transition layer near the interface in which the surfactant concentration varies rapidly. The large surfactant gradients are resolved using a fast and accurate “hybrid” numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the free boundary problem. The method is used to investigate the dependence of drop deformation on parameters that characterize surfactant solubility. We also compute resolved examples of tipstreaming, and investigate its dependence on parameters such as flow rate and bulk surfactant concentration.

Description: Multi-scale analysis is widely adopted in turbulence research for studying flow structures corresponding to specific length scales in the Kolmogorov spectrum. In the present work, a new methodology based on novel optimization techniques for scale decomposition is introduced, which leads to a bandpass filter with prescribed properties. With this filter, we can efficiently perform scale decomposition using Fourier transform directly while adequately suppressing Gibbs ringing artifacts. Both 2D and 3D scale decomposition results are presented, together with qualitative and quantitative analysis. The comparison with existing multi-scale analysis technique is conducted to verify the effectiveness of our method. Validation of this decomposition technique is demonstrated both qualitatively and quantitatively. The advantage of the proposed methodology enables a precise specification of continuous length scales while preserving the original structures. These unique features of the proposed methodology may provide future insights into the evolution of turbulent flow structures.

Description: Direct numerical simulations (DNS) are conducted for a Mach 2.75 turbulent boundary layer interacting with an impinging shock at three different shock incidence angles. The accuracies of DNS calculations are established by checking the convergence of flow statistics for various grids, by comparing the generated results with those in the literature and also by the balance of contributing terms in the turbulent kinetic energy equation. Instantaneous flow visualizations show the significant effect of shock on turbulence structure in the shock-boundary layer interaction zone and also in the flow downstream of the interaction region. The separation bubbles exhibit highly unsteady and three-dimensional behavior and are larger for stronger shocks but the maximum probability of flow separation is found to be independent of the shock strength. The differences between Reynolds- and Favre-averaged quantities are also observed to be small and largely independent of the shock intensity. The turbulent kinetic energy is amplified across the shock, mainly by the production term in the turbulent kinetic energy equation. The amplification of enstrophy across the shock zone is found to be due to the vortex stretching term in the enstrophy transport equation. A detailed examination of the terms in the turbulent kinetic equation shows a strong coupling between the mean and turbulent fields in the interaction region with energy being continuously exchanged from one field to another. However, the compressibility-related terms in the transport equations for turbulent kinetic energy and enstrophy are found to be small for the simulated flows.

Description: Lotus-type porous metals are a promising alternative for compact heat transfer applications. In lotus-type porous fins, jet impingement and transverse mixing play important roles for heat transfer: jets emerging from the pores impinge on the following fin and enhance heat transfer performance, while the transverse fluid motion advects heat away from the fin surface. By means of magnetic resonance imaging we have performed mean flow and scalar transport measurements through scaled-up replicas of two kinds of lotus-type porous fins: one with a deterministic hole pattern and staggered alignment, and one with a random hole pattern, but the same porosity and mean pore diameter. The choice of geometric parameters (fin spacing, thickness, porosity, and hole diameter) is based on previous thermal studies. The Reynolds number based on the mean pore diameter and inner velocity ranges from 80 to 3800. The measurements show that in the random hole pattern the jet characteristic length scale is substantially larger with respect to the staggered hole pattern. The random geometry also produces long coherent vortices aligned with the streamwise direction, which improves the transverse mixing. The random hole distribution causes the time mean streamlines to meander in a random-walk manner, and the diffusivity coefficient associated to the mechanical dispersion (which is nominally zero in the staggered hole configuration) is several times larger than the fluid molecular diffusivity at the higher Reynolds numbers. From the trends in maximum streamwise velocity, streamwise vorticity, and mechanical diffusivity, it is inferred that the flow undergoes a transition to an unsteady/turbulent regime around Reynolds number 300. This is supported by the measurements of concentration of an isokinetic non-buoyant plume of scalar injected upstream of the stack of fins. The total scalar diffusivity for the fully turbulent regime is found to be 22 times larger than the molecular diffusivity, but only 6 times higher than the mechanical diffusivity, indicating that the latter plays a significant role for heat transfer and mixing.

Description: Identifying sources and sinks of N 2 O can illuminate N cycling processes in marine systems, particularly where changes in dissolved O 2 can lead to changes in N cycling pathways (i.e., nitrification vs. denitrification). We measured N 2 O and NO 3 - concentration and their stable isotope ratios (δ 15 N and δ 18 O) in the water column and sediments of the oxygen minimum zone in the near-shore eastern subtropical North Pacific (23º to 34º N). Atmospheric efflux of N 2 O ranged from 2.2 to 17.9 μmol m -2 d -1 , or about 2 to 20 times higher than in oxygenated regions of the north Pacific. Surface waters were a source of 15 N-depleted and 18 O-enriched N 2 O to the atmosphere, indicating a bacterial, not archaeal, nitrification N 2 O source. Stable isotopes indicated that nitrification in both surface and intermediate waters (~ 0 to 200 m) was the major source of N 2 O in this study area, with denitrification acting as a small N 2 O sink in strongly O 2 -depleted waters. Denitrification had a larger impact on observed patterns of N 2 O and NO 3 - concentrations and isotope ratios in the southern oxygen minimum zone. Sediments were generally neutral or a weak sink for N 2 O, with only one site (Soledad basin) showing a positive efflux of +3.5 ± 1.0 μmol N 2 O-N m -2 d -1 . Sediment fluxes of N 2 O at all sites were several orders of magnitude smaller than fluxes of dinitrogen, nitrate, and ammonium measured in previous studies and did not appear to impact water column N 2 O concentrations. N 2 O was less than 0.1% of the N 2 efflux from sedimentary denitrification.

Description: Typhoon Megi passed between two subsurface moorings in the northern South China Sea in October 2010 and the upper ocean thermal and dynamical response with strong internal tides present was examined in detail. The entire observed water column (60-360 m) was cooled due to strong Ekman-pumped upwelling (up to 50 m in the thermocline) by Megi, with maximum cooling of 4.2 °C occurring in thermocline. A relatively weak (maximum amplitude of 0.4 m s -1 ) and quickly damped (e-folding timescale of 2 inertial periods) near-inertial oscillation (NIO) was observed in the mixed layer. Power spectrum and wavelet analyses both indicated an energy peak appearing at exactly the sum frequency fD1 (with maximum amplitude up to 0.2 m s -1 ) of NIO ( f ) and diurnal tide ( D1 ), indicating enhanced nonlinear wave-wave interaction between f and D1 during and after typhoon. Numerical experiments suggested that energy transfer from NIO to fD1 via nonlinear interaction between f and D1 may have limited the growth and accelerated the damping of mixed layer NIO generated by Megi. The occurrence of fD1 had a high correlation with NIO; the vertical nonlinear momentum term, associated with the vertical shear of NIO and vertical velocity of D1 or vertical shear of D1 and vertical velocity of NIO, was more than 10 times larger than the horizontal terms and was responsible for forcing fD1 . After Megi, surface-layer diurnal energy was enhanced by up to 100%, attributed to the combined effect of the increased surface-layer stratification and additional Megi-forced diurnal current.