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: ABSTRACT A reactive transport modeling framework is presented that allows simultaneous assessment of groundwater flow, water quality evolution including δ 13 C, and 14 C activity or “age”. Through application of this framework, simulated 14 C activities can be directly compared with measured 14 C activities. This bypasses the need for interpretation of a 14 C age prior to flow simulation through factoring out processes other than radioactive decay, which typically involves simplifying assumptions regarding spatial and temporal variability in reactions, flow, and mixing. The utility of the approach is demonstrated for an aquifer system with spatially variable carbonate mineral distribution, multiple organic carbon sources, and transient boundary conditions for 14 C activity in the recharge water. In this case the simulated 14 C age was shown to be relatively insensitive to isotopic fractionation during DOC oxidation and variations in assumed DOC degradation behaviour. We demonstrate that the model allows quantitative testing of hypotheses regarding controls on groundwater age and water quality evolution for all three carbon isotopes. The approach also facilitates incorporation of multiple environmental tracers and combination with parameter optimization techniques. This article is protected by copyright. All rights reserved.

Description: Water scarcity is likely to increase in the coming years, making improvements in irrigation efficiency increasingly important. An emerging technology that promises to increase irrigation efficiency substantially is a wireless irrigation sensor network that uploads sensor data into irrigation management software, creating an integrated system that allows real-time monitoring and control of moisture status that has been shown in experimental settings to reduce irrigation costs, lower plant loss rates, shorten production times, decrease pesticide application, and increase yield, quality, and profit. We use an original survey to investigate likely initial acceptance, ceiling adoption rates, and profitability of this new sensor network technology in the nursery and greenhouse industry. We find that adoption rates for a base system and demand for expansion components are decreasing in price, as expected. The price elasticity of the probability of adoption suggests that sensor networks are likely to diffuse at a rate somewhat greater than that of drip irrigation. Adoption rates for a base system and demand for expansion components are increasing in specialization in ornamental production: Growers earning greater shares of revenue from greenhouse and nursery operations are willing to pay more for a base system and are willing to purchase larger numbers of expansion components at any given price. We estimate that growers who are willing to purchase a sensor network expect investment in this technology to generate significant profit, consistent with findings from experimental studies. 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: 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: This technical note presents a useful methodology for studying how the variance of hydraulic and/or reactive attributes of an aquifer are linked to the multi-scaled and hierarchical sedimentary architecture of the aquifer. A new recursive equation is derived which quantitatively describes how the variance is related to sedimentary facies defined at all scales across an entire stratal hierarchy. As compared to prior published equations that emphasize differences in means among facies populations within a hierarchical level, it emphasizes differences across levels. Because of the hierarchical relationships among the terms of the equation, we find it to be useful for conducting a holistic analysis of the relative contributions to the variance arising from all facies types defined across all scales. The methodology is demonstrated using appropriate field data, and is shown to be useful in defining parsimonious classification systems.

Description: ABSTRACT A primary concern for geologic carbon storage is the potential for leakage of stored carbon dioxide (CO 2 ) into the shallow subsurface where it could degrade the quality of groundwater and surface water. In order to predict and mitigate the potentially negative impacts of CO 2 leakage, it is important to understand the physical processes that CO 2 will undergo as it moves through naturally heterogeneous porous media formations. Previous studies have shown that heterogeneity can enhance the evolution of gas phase CO 2 in some cases, but the conditions under which this occurs have not yet been quantitatively defined, nor tested through laboratory experiments. This study quantitatively investigates the effects of geologic heterogeneity on the process of gas phase CO 2 evolution in shallow aquifers through an extensive set of experiments conducted in a column that was packed with layers of various test sands. Soil moisture sensors were utilized to observe the formation of gas phase near the porous media interfaces. Results indicate that the conditions under which heterogeneity controls gas phase evolution can be successfully predicted through analysis of simple parameters, including the dissolved CO 2 concentration in the flowing water, the distance between the heterogeneity and the leakage location, and some fundamental properties of the porous media. Results also show that interfaces where a less permeable material overlies a more permeable material affect gas phase evolution more significantly than interfaces with the opposite layering.

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: Fluvial sediment loads are frequently calculated with rating curves fit to measured sediment transport rates. Rating curves are often treated as statistical representations in which the fitted parameters have little or no physical meaning. Such models, however, may produce large errors when extrapolation is needed, and they provide no insight into the sediment transport process. It is shown that log-linear least squares, the usual method for fitting rating curves, does not generally produce physically meaningful parameter values. In addition, it cannot accommodate data that include zero-transport samples. Alternative fitting methods based non-linear least squares and on maximum likelihood parameter estimation are described and evaluated. The maximum likelihood approach is shown to fit synthetic data better than linear or non-linear least squares, and to perform well with data that include zero-transport samples. In contrast, non-linear least squares methods produce large errors in the parameter estimates when zero-transport samples are present or when the variance structure of the data is incorrectly specified. Analyses with fractional bedload data from a mountain stream suggest that bedload transport rates are gamma distributed, that the arrivals of bedload particles in a sampler conform to a Poisson distribution, and that the variance of non-zero samples can be expressed as a power function of the mean. Preliminary physical interpretations of variations in the rating curve parameters fit to fractional bedload data with the maximum likelihood method are proposed, and their relation to some previous interpretations of rating curve parameters are briefly discussed. This article is protected by copyright. All rights reserved.

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: 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: 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: This paper presents the results of a comprehensive model-based analysis of a uranyl [U(VI)] tracer test conducted at the U.S. DOE Hanford 300 Area (300A) IFRC. Despite the highly complex field conditions the numerical three-dimensional multi-component reactive transport model was able to capture most of the spatiotemporal variations of the observed U(VI) concentrations. A multi-model analysis was performed to interrogate the relative importance of various processes and factors for controlling field-scale reactive transport during the uranyl tracer test. The results indicate that multi-rate sorption/desorption, surface complexation reactions, and initial concentrations were the most important processes and factors controlling U(VI) migration. On the other hand, cation exchange reactions, the choice of the surface complexation model, and dual-domain mass transfer processes played less important roles under the prevailing field-test condition. Further analysis of the modeling results demonstrates that these findings are conditioned to the relatively stable groundwater chemistry and the selected length of the field experimental duration (16 days). The model analysis also revealed the crucial role of the intraborehole flow that occurred within the long-screened monitoring wells and thus affected both field measurements and simulated U(VI) concentrations as a combined effect of aquifer heterogeneity and dynamic flow conditions. This study provides the first highly data-constrained uranium transport simulations under highly dynamic flow conditions. It illustrates the value of reactive transport modeling for elucidating the relative importance of individual processes in controlling uranium transport under specific field-scale conditions.

Description: The dynamics of drying processes from porous media are critically influenced by the intensity of an adjacent free flow and by processes at the interface between free flow and the porous medium. In this paper, the influence of hydraulic properties of a porous medium and of the interaction between fluids and porous medium on the drying dynamics during the capillary-flow dominated stage-1 and transition to the diffusion-dominated stage-2 are studied using a coupled free-flow - porous-medium flow model on the REV scale. We present a detailed model concept that considers mass balance equations, an energy balance equation and the coupling to the adjacent free flow. Key microscale processes are identified and incorporated in the macroscale description of the evaporation process. Own experimental results are used to illustrate main features of the modeling framework. We demonstrate that the use of a homogeneous distribution of soil parameters without consideration of pore-scale induced nonlinearities in the numerical simulations results in a rather constant drying rate in stage-1, which was not observed for the high evaporative demand in the experiments. To account for the dependency of the drying rate on the surface moisture content, special conditions based on the work of Haghighi et al. [2013] and Schlünder [1988] are analyzed for their applicability on the REV scale. Typical features of a drying process, such as different stages of the drying rate, could be reproduced.

Description: Making useful predictions in ungauged basins is an incredibly difficult task given the limitations of hydrologic models to represent physical processes appropriately across the heterogeneity within and among different catchments. Here, we introduce a new method for this challenge, Bayes empirical Bayes, that allows for the statistical pooling of information from multiple donor catchments and provides the ability to transfer parametric distributions rather than single parameter sets to the ungauged catchment. Further, the methodology provides an efficient framework with which to formally assess predictive uncertainty at the ungauged catchment. We investigated the utility of the methodology under both synthetic and real data conditions, and with respect to its sensitivity to the number and quality of the donor catchments used. This study highlighted the ability of the hierarchical Bayes empirical Bayes approach to produce expected outcomes in both the synthetic and real data applications. The method was found to be sensitive to the quality (hydrologic similarity) of the donor catchments used. Results were less sensitive to the number of donor catchments, but indicated that predictive uncertainty was best constrained with larger numbers of donor catchments (but still adequate with fewer donors)

Description: Autumn is a season of dynamic change in forest streams of the northeastern USA due to effects of leaf fall on both hydrology and biogeochemistry. Few studies have explored how interactions of biogeochemical transformations, various nitrogen sources, and catchment flowpaths affect stream nitrogen variation during autumn. To provide more information on this critical period, we studied 1) the timing, duration, and magnitude of changes to stream nitrate, dissolved organic nitrogen (DON), and ammonium concentrations; 2) changes in nitrate sources and cycling; and 3) source areas of the landscape that most influence stream nitrogen. We collected samples at higher temporal resolution for a longer duration than typical studies of stream nitrogen during autumn. This sampling scheme encompassed the patterns and extremes that occurred during baseflow and stormflow events of autumn. Baseflow nitrate concentrations decreased by an order of magnitude from 5.4 to 0.7 μmol L -1 during the week when most leaves fell from deciduous trees. Changes to rates of biogeochemical transformations during autumn baseflow explained the low nitrate concentrations; in-stream transformations retained up to 72% of the nitrate that entered a stream reach. A decrease of in-stream nitrification coupled with assimilatory nitrate uptake was a primary factor in the seasonal nitrate decline. The period of low nitrate concentrations ended with a storm event in which stream nitrate concentrations increased by 25 fold. In the ensuing weeks, stormflow nitrate concentrations progressively decreased over closely-spaced, yet similarly sized events. Most stormflow nitrate originated from nitrification in near-stream areas with occasional, large inputs of unprocessed atmospheric nitrate, which has rarely been reported for non-snowmelt events. A maximum input of 33% unprocessed atmospheric nitrate to the stream occurred during one event. The large inputs of unprocessed atmospheric nitrate show direct and rapid effects on forest streams that may be widespread, although undocumented, throughout nitrogen-polluted temperate forests. In contrast to a week-long nitrate decline during peak autumn litterfall, baseflow DON concentrations increased after leaf fall and remained high for two months. Dissolved organic nitrogen was hydrologically flushed to the stream from riparian soils during stormflow. In contrast to distinct seasonal changes in baseflow nitrate and DON concentrations, ammonium concentrations were typically at or below detection limit, similar to the rest of the year. Our findings reveal couplings among catchment flow paths, nutrient sources and transformations that control seasonal extremes of stream nitrogen in forested landscapes.

Description: During the recent years there has been an increasing interest in multivariate frequency analysis of hydrological variables, e.g. those describing extreme events like rainfall, floods or droughts. The multivariate analysis provides a better understanding of the phenomena under investigation and an additional insight about the interrelationships between the different variables (e.g. peak, volume and duration of the flood), exploiting the complete structure of the problem and making a full use of the available data. However, while the developments on multivariate analysis of hydrological data has produced a large body of literature, a clear assessment of the use of these methods in the design and risk assessment of hydraulic structures is still a matter of debate. In the present work we illustrate a general, structure-based framework for the design and/or risk assessment of hydraulic structures in a bivariate environment; we also compare it to recently proposed methods which are based on the assumption of hydrological design events (as is customary in the univariate context). For illustration purposes, both the structure-based and the design event-based approaches are applied to the design of an idealized structure, thus exploring the differences among the methods as function of the parameters involved. Our work highlights that the return period of structure failure in a multivariate environment strictly depends on the particular structure under design, and in most cases the design of an hydraulic structure cannot be based on a single, hydrological multivariate design event. This acts as a warning for practitioners against the use of design methods based on single hydrological events, as usually done in the context of univariate hydrology, thus neglecting the interplay between the structure and the hydrological loads acting on it.

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: Proliferation of evapotranspiration (ET) products warrants comparison of these products. The study objective was to assess uncertainty in ET output from four land surface models (LSMs), Noah, Mosaic, VIC, and SAC in NLDAS-2, two remote sensing-based products, MODIS and AVHRR, and GRACE-inferred ET from a water budget with precipitation from PRISM, monitored runoff, and total water storage change (TWSC) from GRACE satellites. The three cornered hat method, which does not require a priori knowledge of the true ET value, was used to estimate ET uncertainties. In addition, TWSC or total water storage anomaly (TWSA) from GRACE was compared with water budget estimates of TWSC from a flux-based approach or TWSA from a storage-based approach. The analyses were conducted using data from three regions (humid – arid) in the South Central US as case studies. Uncertainties in ET are lowest in LSM ET (~5 mm/month), moderate in MODIS- or AVHRR-based ET (10 – 15 mm/month), and highest in GRACE-inferred ET (20 – 30 mm/month). There is a tradeoff between spatial resolution and uncertainty, with lower uncertainty in the coarser-resolution LSM ET (~14 km) relative to higher uncertainty in the finer-resolution (~ 1 ‒ 8 km) RS ET. Root-mean-square (RMS) of uncertainties in water budget estimates of TWSC is about half of RMS of uncertainties in GRACE-derived TWSC for each of the regions. Future ET estimation should consider a hybrid approach that integrates strengths of LSMs and satellite-based products to constrain uncertainties.

Description: It is well documented that deforestation results in an increase in landslide frequency due to the control that forest roots have on slope stability. The loss of forest vegetation leads to a reduction in soil cohesion and a decrease in the shear strength of the soil profile. As a result, the slope becomes more susceptible to landsliding and the return time of landslides decreases. When a landslide removes the soil profile, there may not be adequate time for seedlings to grow and enhance soil stability. In this study, we investigate whether bistable dynamics emerge from the interaction of forest vegetation with the formation and accumulation of colluvial deposits in soil-mantled landscapes. To that end, we develop deterministic and stochastic models of landslide occurrence with a dynamic vegetation component. Results show that bistability exists for the deterministic case for both steep and shallow hollows under event and supply limited conditions. However, for the stochastic case, the randomness of landslide occurrence largely changed the states of the system such that the system only exhibited one stable state, which was the fully vegetated condition. Examining different management practices under stochastic conditions showed that the system eventually recovered; however, management practices influenced the recovery time of the forest. Thus, different management practices could render the land in a state of low vegetation over economically significant time periods.

Description: Seepage flux from ephemeral streams can be an important component of the water balance in arid and semi-arid regions. An emerging technique for quantifying this flux involves the measurement and simulation of a flood wave as it moves along an initially dry channel. This study investigates the usefulness of including surface water and groundwater data to improve model calibration when using this technique. We trialed this approach using a controlled flow event along a 1387 m reach of artificial stream channel. Observations were then simulated using a numerical model that combines the diffusion wave approximation of the Saint-Vénant equations for streamflow routing, with Philips’ infiltration equation and the groundwater flow equation. Model estimates of seepage flux for the upstream segments of the study reach, where streambed hydraulic conductivities were approximately 10 1 m d -1 , were on the order of 10 -4 m 3 d -1 m -2 . In the downstream segments, streambed hydraulic conductivities were generally much lower but highly variable (~10 -3 – 10 -7 m d -1 ). The Latin Hypercube Monte Carlo sensitivity analysis showed that the flood front timing, surface water stage, groundwater heads and the predicted stream flow seepage were most influenced by specific yield. Furthermore, inclusion of groundwater data resulted in a higher estimate of total seepage estimates than if the flood front timing were used alone.

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: We present a probabilistic sediment cascade model to simulate sediment transfer in a mountain basin (Illgraben, Switzerland) where sediment is produced by hillslope landslides and rockfalls and exported out of the basin by debris flows and floods. The model conceptualizes the fluvial system as a spatially lumped cascade of connected reservoirs representing hillslope and channel storages where sediment goes through cycles of storage and remobilization by surface runoff. The model includes all relevant hydrological processes that lead to runoff formation in an Alpine basin, such as precipitation, snow accumulation, snow melt, evapotranspiration, and soil water storage. Although the processes of sediment transfer and debris flow generation are described in a simplified manner, the model produces complex sediment discharge behavior which is driven by the availability of sediment and antecedent wetness conditions (system memory) as well as the triggering potential (climatic forcing). The observed probability distribution of debris flow volumes and their seasonality in 2000-2009 are reproduced. The stochasticity of hillslope sediment input is important for reproducing realistic sediment storage variability, although many details of the hillslope landslide triggering procedures are filtered out by the sediment transfer system. The model allows us to explicitly quantify the division into transport and supply-limited sediment discharge events. We show that debris flows may be generated for a wide range of rainfall intensities because of variable antecedent basin wetness and snowmelt contribution to runoff, which helps to understand the limitations of methods based on a single rainfall threshold for debris flow initiation in Alpine basins.

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: Characterizing the complex geometries and the heterogeneity of the deposits in meandering river systems is a long-standing issue for the 3D modeling of alluvial formations. Such deposits are important sources of accessible groundwater in alluvial aquifers throughout the world and also play a major role as hydrocarbons reservoirs. In this paper we present a method to generate meandering river centerlines that are stochastic, geologically realistic, connected and conditioned to local observations or global geomorphological characteristics. The method is based on fast 1D multiple-point statistics in a transformed curvilinear domain: the succession in directions observed in a real world meandering river (the analog) is considered as statistical model for multiple-point statistics simulation. The integration of local data is accomplished by an inverse procedure ensuring that the channels pass through a given set of locations while conserving the high-order spatial characteristics of an analog. The methodology is applied on seven real world case studies. This work demonstrates the flexibility and the applicability of multiple-point statistics outside the standard paradigm that considers the simulation of a 2D or 3D variable with spatial coordinates.

Description: Large asymmetric bedforms commonly develop in rivers. The turbulence associated with flow separation that develops over their steep lee side is responsible for the form shear stress which can represent a substantial part of total shear stress in rivers. This paper uses the Delft3D modeling system to investigate the effects of bedform geometry and forcing conditions on flow separation length and associated turbulence, and bedform shear stress over angle-of-repose (30° lee side angle) bedforms. The model was validated with laboratory measurements that showed sufficient agreement to be used for a systematic analysis. The influence of flow velocity, bed roughness, relative height (bedform height / water depth) and aspect ratio (bedform height / length) on the variations of the normalized length of the flow separation zone, the extent of the wake region (where the turbulent kinetic energy (TKE) was more than 70% of the maximum TKE), the average TKE within the wake region and the form shear stress were investigated. Form shear stress was found not to scale with the size of the flow separation zone but to be related to the product of the normalized extent of the wake region (extent of the wake region / extent of water body above the bedform) and the average TKE within the wake region. The results add to understanding of the hydrodynamics of bedforms and may be used for the development of better parameterizations of small-scale processes for application in large-scale studies.

Description: We explore the bankfull width ( W bf ) vs. drainage area ( A da ) relationship across a range of climatic and geologic environments, and ask (1) is the relationship between ln ( W bf ) and ln ( A da ) best described by a linear function and (2) can a reliable relationship be developed for predicting W bf with A da as the only independent variable. The principal dataset for this study was compiled from regional curve studies and other reports that represent 1,018 sites (1 m ≤ W bf ≤ 110 m and 0.50 km 2 ≤ A da ≤ 22,000 km 2 ) in the continental U.S. Two additional datasets were used for validation. After dividing the data into small-, medium-, and large-size basins which, respectfully, correspond to A da 〈 4.95 km 2 , 4.95 km 2 ≤ A da 〈 337 km 2 , and A da ≥ 337 km 2 , regression lines from each dataset were compared using one-way analysis of covariance (ANCOVA). A second ANCOVA was performed to determine if mean annual precipitation ( P ) is an extraneous factor in the W bf vs. A da relationship. The ANCOVA results reveal that using A da alone does not yield a reliable W bf vs. A da relationship that is applicable across a wide range of environments and that P is a significant extraneous factor in the relationship. Considering data for very small basins ( A da ≤ 0.49 km 2 ) and very large basins ( A da ≥ 1.0×10 5 km 2 ) we conclude that a two-segment linear model is the most probable form of the ln ( W bf ) vs. ln ( A da ) relationship. This study provides useful information for building complex multivariate models for predicting W bf .

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: 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: Recent studies show that multimodel combinations improve hydroclimatic predictions by reducing model uncertainty. Given that climate forecasts are available from multiple climate models, which could be ingested with multiple watershed models, what is the best strategy to reduce the uncertainty in streamflow forecasts? To address this question, we consider three possible strategies: (1) reduce the input uncertainty first by combining climate models and then use the multimodel climate forecasts with multiple watershed models (MM-P) (2) ingest the individual climate forecasts (without multimodel combination) with various watershed models and then combine the streamflow predictions that arise from all possible combinations of climate and watershed models (MM-Q)(3) combine the streamflow forecasts obtained from multiple watershed models based on strategy (1) to develop a single streamflow prediction that reduces uncertainty in both climate forecasts and watershed models (MM-PQ) . For this purpose, we consider synthetic schemes that generate streamflow and climate forecasts, for comparing the performance of three strategies with the true streamflow generated by a given hydrologic model. Results from the synthetic study show that reducing input uncertainty first ( MM-P ) by combining climate forecasts results in reduced error in predicting the true streamflow compared to the error of multimodel streamflow forecasts obtained by combining streamflow forecasts from all-possible combination of individual climate model with various hydrologic models ( MM-Q ). Since the true hydrologic model structure is unknown, it is desirable to consider MM-PQ as an alternate choice that reduces both input uncertainty and hydrologic model uncertainty. Application on two watersheds in NC also indicates that reducing the input uncertainty first is critical before reducing the hydrologic model uncertainty.

Description: The understanding of reasons leading to non-uniqueness of soil erosion susceptibility is still inadequate, yet indispensable for establishing general relations between runoff volume and sediment yield. To obtain relevant insights, we performed a series of numerical simulations with a detailed hydrodynamic model using synthetic storms of varying intensity, duration, and lag time between events as representations of different hydrologic response conditions in a zero-order catchment. The design targeted to generate surface flow and ‘perturb’ soil substrate by a first rainfall event, creating a set of initial conditions in terms of flow and deposited sediment prior to the onset of a subsequent rainfall event. Due to the differential effect of (re)detachment and (re)entrainment processes on soil particles of varying sizes, the deposited sediment mass formed shielding layer. One of the essential results is that unless the initial condition of flow and sediment is identical, the same volume of runoff can generate different total sediment yields and their variation can reach up to ~200%. The effect is attributed to two major conflicting effects exerted by the deposited ‘initialization’ (soil antecedent condition) sediment mass: erosion enhancement, because of supply of highly erodible sediment, and erosion impediment, because of constrain on the availability of lighter particles by heavier sediment. Consistently with this inference, long-term simulations with continuous rainfall show that a peculiar feature of sediment yield series is the existence of maximum before the steady-state is reached. The two characteristic time scales, the time to peak and the time to steady-state, separate three characteristic periods that correspond to flow-limited, source-limited, and steady-state regimes. These time scales are log-linearly and negatively related to the spatially averaged Shields parameter: the smaller the rainfall input and the heavier a given particle is, the larger the two scales are. The results provide insights on how the existence of shield operates on erosion processes, possibly implying that accurate short-term predictions of geomorphic events from headwater areas may never become a tractable problem: the latter would require a detailed spatial characterization of particle size distribution prior to precipitation events.

Description: Aquifer hydraulic properties such as hydraulic conductivity ( K ) are ubiquitously heterogeneous and typically only a statistical characterization can be sought. Additionally statistical anisotropy at typical characterization scales is the rule. Thus, regardless of the processes governing solute transport at the local (pore) scale, transport becomes non-Fickian. Mass-transfer models provide an efficient tool that reproduces observed anomalous transport; in some cases though, these models lack predictability as model parameters cannot readily be connected to the physical properties of aquifers. In this study we focus on a multi-rate mass-transfer model (MRMT), and in particular the apparent capacity coefficient (β), which is a strong indicator of the potential of immobile zones to capture moving solute. We aim to find if the choice of an apparent β can be phenomenologically related to measures of statistical anisotropy. We analyzed an ensemble of random simulations of three-dimensional log-transformed multi-Gaussian permeability fields with stationary anisotropic correlation under convergent flow conditions. It was found that apparent β also displays an anisotropic behavior, physically controlled by the aquifer directional connectivity, which in turn is controlled by the anisotropic correlation model. A high hydraulic connectivity results in large β values. These results provide new insights into the practical use of mass-transfer models for predictive purposes.

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: In this paper, an operational algorithm is proposed for the mapping of surface moisture over the northern and central parts of Tunisia, in North Africa. A change detection approach is applied, using 160 multi-incidence Envisat ASAR Wide Swath images acquired in the horizontal polarization over a 7-year period. Parameterization of this algorithm is considered for three classes of vegetation cover density (NDVI〈0.25, 0.25〈NDVI〈0.5 and NDVI〉0.5), retrieved from SPOT-VGT decadal images. A relative soil moisture index, ranging between 0 (for the driest surfaces) and 1 (for saturated soils), is proposed for each date, with a resolution of 1 km. The retrieved soil moistures are validated by means of ground measurements based on continuous thetaprobe measurements, as well as low resolution (25 km) ERS and ASCAT soil moisture products from the Vienna University of Technology (TU Wien). A qualitative relationship between spatio-temporal variations of moisture and precipitation is also discussed.

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: 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: This paper presents the results of a data based comparative study of several hundred catchments across continental United States belonging to the MOPEX dataset to systematically explore the connection between the flood frequency curve and mean annual water balance. Mean annual water balance is expressed in terms of two similarity measures: (i) the climatic aridity index, AI , which is a measure of the competition between energy and water availability; and (ii) the baseflow index, BFI , which is a measure of total runoff partitioning into surface and subsurface components at the annual time scale. The data analyses showed that the aridity index, AI , has a first order control on the shape of the flood frequency curve (also known as the growth curve), as expressed in terms of both the mean and coefficient of variation ( C v ) of the annual maximum floods, once normalized by catchment size (i.e., specific flood discharge) While the mean annual (specific) flood discharge decreases with increasing aridity, C v increases with increasing aridity. On the other hand, the BFI was found to be a second order control on the flood frequency curve. Higher BFI , meaning higher contributions of subsurface flow to total streamflow, leads to a decrease of the mean annual (specific) flood discharge, and vice versa. The statistically significant relationship between AI and the flood frequency curve and the consistent shift of the growth curves with AI support the use of AI as a similarity measure for regionalization of flood frequency.

Description: The performance of glacio-hydrological models which simulate catchment response to climate variability depends to a large degree on the data used to force the models. The forcing data become increasingly important in high elevation, glacierised catchments where the interplay between extreme topography, climate and the cryosphere is complex. It is challenging to generate a reliable forcing dataset that captures this spatial heterogeneity. In this paper, we analyze the results of a one year field campaign focusing on air temperature and precipitation observations in the Langtang Valley in the Nepalese Himalayas. We use the observed time series to characterize both temperature lapse rates (LRs) and precipitation gradients (PGs). We study their spatial and temporal variability, and we attempt to identify possible controlling factors. We show that very clear LRs exist in the valley and that there are strong seasonal differences related to the water vapor content in the atmosphere. Results also show that the LRs are generally shallower than the commonly used environmental lapse rates. The analysis of the precipitation observations reveal that there is great variability in precipitation over short horizontal distances. A uniform valley wide PG cannot be established and several scale-dependent mechanisms may explain our observations. We complete our analysis by showing the impact of the observed LRs and PGs on the outputs of the TOPKAPI-ETH glacio-hydrological model. We conclude that LRs and PGs have a very large impact on the water balance composition and that short-term monitoring campaigns have the potential to improve model quality considerably.

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: Many physically-based hydrological/hydrogeological models used for predicting groundwater seepage areas, including topography-based index models such as TOPMODEL, rely on the Dupuit assumption. To ensure the sound use of these simplified models, knowledge of the conditions under which they provide a reasonable approximation is critical. In this study, a Dupuit solution for the seepage length in hillslope cross-sections is tested against a full-depth solution of saturated groundwater flow. In homogeneous hillslopes with horizontal impervious base and constant-slope topography, the comparison reveals that the validity of the Dupuit solution depends not only on the ratio of depth to hillslope length d/L (as might be expected), but also on the ratio of hydraulic conductivity to recharge K/R and on the topographic slope s . The validity of the Dupuit solution is shown to be in fact a unique function of another ratio, the ratio of depth to seepage length d/L S . For d/L S 〈 0.2, the relative difference between the two solutions is quite small (〈 14% for the wide range of parameter values tested), whereas for d/L S 〉 0.2, it increases dramatically. In practice, this criterion can be used to test the validity of Dupuit solutions. When d/L S increases beyond that cut-off, the ratio of seepage length to hillslope length L S /L given by the full-depth solution tends towards a non-zero asymptotic value. This asymptotic value is shown to be controlled by (and in many cases equal to) the parameter R/sK . Generalization of the findings to cases featuring heterogeneity, non-horizontal impervious base and variable-slope topography is discussed.

Description: Isothermal compositional flow models require coupling transient compressible flows and advective transport systems of various chemical species in subsurface porous media. Building such numerical models is quite challenging and may be subject to many sources of uncertainties because of possible incomplete representation of some geological parameters that characterize the system's processes. Advanced data assimilation methods, such as the ensemble Kalman filter (EnKF), can be used to calibrate these models by incorporating available data. In this work, we consider the problem of estimating reservoir permeability using information about phase pressure as well as the chemical properties of fluid components. We carry out state-parameter estimation experiments using joint and dual updating schemes in the context of the EnKF with a two-dimensional single-phase compositional flow model (CFM). Quantitative and statistical analyses are performed to evaluate and compare the performance of the assimilation schemes. Our results indicate that including chemical composition data significantly enhances the accuracy of the permeability estimates. In addition, composition data provide more information to estimate system states and parameters than do standard pressure data. The dual state-parameter estimation scheme provides about 10% more accurate permeability estimates on average than the joint scheme when implemented with the same ensemble members, at the cost of twice more forward model integrations. At similar computational cost, the dual approach becomes only beneficial after using large enough ensembles.

Description: Groundwater-fed irrigation has been shown to deplete groundwater storage, decrease surface water runoff and increase evapotranspiration. Here we simulate soil moisture dependent groundwater-fed irrigation with an integrated hydrologic model. This allows for direct consideration of feedbacks between irrigation demand and groundwater depth. Special attention is paid to system dynamics in order to characterized spatial variability in irrigation demand and response to increased irrigation stress. A total of 80 years of simulation are completed for the Little Washita Basin in Southwestern Oklahoma, USA spanning a range of agricultural development scenarios and management practices. Results show regionally aggregated irrigation impacts consistent with other studies. However, here a spectral analysis reveals that groundwater-fed irrigation is also shown to amplify the annual streamflow cycle while dampening longer-term cyclical behavior with increased irrigation during climatological dry periods. Feedbacks between the managed and natural system are clearly observed with respect to both irrigation demand and utilization when water table depths are within a critical range. Although the model domain is heterogeneous with respect to both surface and subsurface parameters, relationships between irrigation demand, water table depth and irrigation utilization are consistent across space and between scenarios. Still, significant local heterogeneities are observed both with respect to transient behavior and response to stress. Spatial analysis of transient behavior shows that farms with groundwater depths within a critical depth range are most sensitive to management changes. Differences in behavior highlight the importance of groundwater's role in system dynamics in addition to water availability.

Description: Including positive feedbacks in hydrological models has recently been shown to result in complex behavior with multiple steady states. When a large disturbance, say a major drought, is simulated within such models the hydrology changes. Once the disturbance ends the hydrology does not return to that prior to the disturbance, but rather, persists within an alternate state. These multiple steady states (henceforth attractors ) exist for a single model parameterization and cause the system to have a finite resilience to such transient disturbances. A limitation of past hydrological resilience studies is that multiple attractors have been identified using mean annual or mean monthly forcing. Considering that most hydrological fluxes are subject to significant forcing stochasticity and do not operate at such large time scales, it remains an open question whether multiple hydrological attractors can exist when a catchment is subject to stochastic daily forcing. This question is the focus of this paper and it needs to be addressed prior to searching for multiple hydrological attractors in the field. To investigate this, a previously developed semi-distributed hill-slope ecohydrological model was adopted which exhibited multiple steady states under average monthly climate forcing. In this paper, the ecohydrological model was used to explore if feedbacks between the vegetation and a saline water table result in two attractors existing under daily stochastic forcing. The attractors and the threshold between them (henceforth repellor ) were quantified using a new limit cycle continuation technique that up-scaled climate forcing from daily to monthly (model and limit cycle code is freely available). The method was used to determine the values of saturated lateral hydraulic conductivity at which multiple attractors exist. These estimates were then assessed against time-integration estimates, which they agreed with. Overall, multiple attractors where found to exist under stochastic daily forcing. However, changing the climate forcing from monthly to daily did significantly reduce the parameter range over which two attractors existed. This suggests fewer catchments may have multiple attractors than previously considered.

Description: The companion paper showed that multiple steady state groundwater levels can exist within a hill-slope Boussinesq-vegetation model under daily stochastic forcing. Using a numerical limit-cycle continuation algorithm, the steady states (henceforth attractors ) and the threshold between them (henceforth repellor ) were quantified at a range of saturated lateral conductivity values, . This paper investigates if stochastic daily forcing can switch the catchment between both of the attractors. That is, an attractor may exist under average forcing conditions but can stochastic forcing switch the catchment into and out of each of the attractor basins?; i.e. making the attractor emerge . This was undertaken using the model of the companion paper and by completing daily time-integration simulations at six values of the saturated lateral hydraulic conductivity, ; three having two attractors and three having only a deep water table attractor. By graphically analyzing the simulations, and comparing against simulations from a model modified to have only one attractor, multiple attractors were found to emerge under stochastic daily forcing. However, the emergence of attractors was significantly more subtle and complex than that suggested by the companion paper. That is, an attractor may exist but never emerge; both attractors may exist and both may emerge but identifying the switching between attractors was often ambiguous; and only one attractor may exist and but a second temporary attractor may exist and emerge during periods of high precipitation. This subtle and complex emergence of attractors was explained using continuation analysis of the climate forcing rate, and not a model parameter such as . It showed that the temporary attractor existed over a large range of values and this suggests that more catchments may have multiple attractors than suggested by the companion paper. By combining this continuation analysis with the time-integration simulations, hydrological signatures indicative of a switch of multiple attractors were proposed. These signatures may provide a means for identifying actual catchments that have switched between multiple attractors.

Description: ABSTRACT A large-time analytical solution is proposed for the spatial variance and coefficient of variation of the depth-averaged concentration due to instantaneous, cross-sectionally uniform solute sources in pseudo-rectangular open channel flows. The mathematical approach is based on the use of the Green functions and on the Fourier decomposition of the depth-averaged velocities, coupled with the method of the images. The variance spatial trend is characterized by a minimum at the center of the mass and two mobile, decaying symmetrical peaks which, at very large times, are located at the inflexion points of the average Gaussian distribution. The coefficient of variation, which provides an estimate of the expected percentage deviation of the depth-averaged point concentrations about the section-average, exhibits a minimum at the center which decays like t -1 and only depends on the river diffusive time-scale. The defect of cross-sectional mixing quickly increases with the distance from the center, and almost linearly at large times. Accurate numerical Lagrangian simulations were performed to validate the analytical results in pre-asymptotic and asymptotic conditions, referring to a particularly representative sample case for which cross-sectional depth and velocity measurements were known from a field survey. In addition, in order to discuss the practical usefulness of computing large-time concentration spatial moments in river flows, and resorting to directly measured input data, the order of magnitude of section-averaged concentrations and corresponding coefficients of variation was estimated in field conditions and for hypothetical contamination scenarios, considering a unit normalized mass impulsively injected across the transverse section of 81 U.S. rivers.

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: Two major causes of salt marsh loss are vertical drowning, when sediment accumulation on the platform cannot keep vertical pace with sea level rise, and horizontal retreat, associated with wave-induced marsh boundary erosion. Despite these processes having been extensively documented and modeled, is unclear which loss modality dominates given a set of environmental parameters. A three-point dynamic model was developed to predict marsh loss as a function of sea level rise, allochthonous sediment supply, wind regime, tidal range, and marsh bank and mudflat erodability. Marsh horizontal and vertical evolution was found to respond in opposing ways to wave induced erosion processes. Marsh horizontal retreat was triggered by large mudflats, strong winds, high erodability of marsh bank and mudflat, whereas the opposite conditions acted to reduce the sediment supply to the marsh platform, promoting marsh loss to drowning. With low and moderate rates of sea level rise (~ 5 mm/yr), retreat was found to be a more likely marsh loss modality than drowning. However, conditions associated with marsh retreat also increase the system resilience by transferring sediment on the marsh platform and preventing drowning. Our results suggest the use of a modular strategy for short-term marsh management: selectively protect extensive salt marsh regions by maintaining healthy vegetation on the platform, while allowing other areas to retreat, leveraging the natural resilience embedded in the lateral loss of marsh extent.

Description: Observational data and the Princeton Urban Canopy Model, with its detailed representation of urban heterogeneity and hydrological processes, are combined to study evaporation and turbulent water vapor transport over urban areas. The analyses focus on periods before and after precipitation events, at two sites in the Northeastern United States. Our results indicate that while evaporation from concrete pavements, building rooftops and asphalt surfaces is discontinuous and intermittent, overall these surfaces accounted for nearly 18% of total latent heat fluxes (LE) during a relatively wet 10-day period. More importantly, these evaporative fluxes have a significant impact on the urban surface energy balance, particularly during the 48 hours following a rain event when impervious evaporation is the highest. Thus, their accurate representation in urban models is critical. Impervious evaporation after rainfall is also shown to correlate the sources of heat and water at the earth surface, resulting in a conditional scalar transport similarity over urban terrain following rain events.

Description: The spatial and temporal dynamics of seasonal snow covers play a critical role for many hydrological, ecological, and climatic processes. This paper presents a new, innovative approach to continuously monitor these dynamics using numerous low-cost, standalone snow monitoring stations (SnoMoS). These stations provide snow and related meteorological data with a high temporal and spatial resolution. Data collected by SnoMoS include: snow depth, surface temperature, air temperature and humidity, total precipitation, global radiation, wind speed, and barometric pressure. A total of 99 sensors were placed over the winters 2010/11 and 2011/12 at multiple locations within three 40 - 180 km² basins in the Black Forest region of Southern Germany. The locations were chosen to cover a wide range of slopes, elevations, and expositions in a stratified sampling design. Furthermore, “paired stations” located in close proximity to each other, one in the open and one underneath various forest canopies, were set up to investigate the influence of vegetation on snow dynamics. The results showed that considerable differences in snow depth and, therefore, snow water equivalent (SWE) are present within the study area despite its moderate temperatures and medium elevation range (400 - 1500 m). The relative impact of topographical factors like elevation, aspect, and of different types of forest vegetation were quantified continuously and were found to change considerably over the winter period. The recorded differences in SWE and snow cover duration were large enough that they should be considered in hydrologic and climate models.

Description: Integrated land surface-groundwater models are valuable tools in simulating the terrestrial hydrologic cycle as a continuous system and exploring the extent of land surface–subsurface interactions from catchment to regional scales. However, the fidelity of model simulations is impacted not only by the vegetation and subsurface parameterizations, but also by the antecedent condition of model state variables, such as the initial soil moisture, depth to groundwater and ground temperature. In land surface modeling, a given model is often run repeatedly over a single year of forcing data until it reaches an equilibrium state: the point at which there is minimal artificial drift in the model state or prognostic variables (most often the soil moisture). For more complex coupled and integrated systems, where there is an increased computational cost of simulation and the number of variables sensitive to initialization is greater than in traditional uncoupled land surface modeling schemes, the challenge is to minimize the impact of initialization while using the smallest spin-up time possible. In this study, multi-criteria analysis was performed to assess the spin-up behavior of the ParFlow.CLM integrated groundwater-surface water-land surface model over a 208 km 2 sub-catchment of the Ringkobing Fjord catchment in Denmark. Various measures of spin-up performance were computed for model state variables such as the soil moisture and groundwater storage, as well as for diagnostic variables such as the latent and sensible heat fluxes. The impacts of initial conditions on surface water–groundwater interactions were then explored. Our analysis illustrates that the determination of an equilibrium state depends strongly on the variable and performance measure used. Choosing an improper initialization of the model can generate simulations that lead to a misinterpretation of land surface-subsurface feedback processes and result in large biases in simulated discharge. Estimated spin-up time from a series of spin-up functions revealed that 20 (or 21) years of simulation were sufficient for the catchment to equilibrate according to at least one criterion at the 0.1% (0.01%) threshold level. Amongst a range of convergence metrics examined, percentage changes in monthly values of groundwater and unsaturated zone storages produced a slow system convergence to equilibrium, whereas criteria based on ground temperature allowed a more rapid spin-up. Slow convergence of unsaturated and saturated zone storages is a result of the dynamic adjustment of the water table in response to a physically arbitrary or inconsistent initialization of a spatially uniform water table. Achieving equilibrium in subsurface storage ensured equilibrium across a spectrum of other variables, hence providing a good measure of system-wide equilibrium. Overall, results highlight the importance of correctly identifying the key variable affecting model equilibrium and also the need to use a multi-criteria approach to achieve a rapid and stable model spin-up.

Description: In this paper, we perform a detailed regional analysis of the link between meteorological drought indices and streamflow for a comprehensive Austrian dataset of 47 small to medium-size catchments in humid-temperate climate. Four drought indices considering different components of the catchment water balance are tested. We assess the quality of the link using rank correlation analysis, and the probability of detecting low flow events by hit-scores. Overall, correlations range between 0.4 and 0.8 and differ significantly between regions. A stratified analysis shows that the link is much stronger (i) for summer low flows and droughts than for anomalies in general, and (ii) for more humid than more arid conditions. Under more humid conditions streamflow droughts of small to medium-size catchments are to a large extent generated by climate forcing and therefore well represented by a simple meteorological index. Under increasingly dry conditions, the climate signal gets less predictive and storage properties of the underground become more important. A simple soil moisture accounting scheme (such as those of the Palmer index) can considerably improve the correlations. Overall, we conclude there is a significant link between meteorological drought and streamflow drought, except for catchments where groundwater storage and snow processes are important. The results are encouraging and provide a wealth of information which can profitably be used to set up statistical prediction models to estimate low flows from meteorological time series.

Description: Water resources in the western United States are contingent on interannual variations in snowpack. Interannual snowpack variability has been attributed to large-scale climate patterns including the El Niño-Southern Oscillation (ENSO), however the contribution of snowfall frequency and extreme snowfall events to this variability are less well quantified. Long term records from Snowpack Telemetry and Cooperative Observer Program stations in the eleven western states were used to investigate these relationships by considering the number of snowfall days and snowfall water equivalent (SFE) of extreme snowfall events. The top decile of snowfall events contributed 20-38% of annual SFE, depending on the region. An average of 65% and 69% of the interannual variability in annual SFE was explained by snowfall days and SFE of top decile snowfall events, respectively, with extreme events being a more significant predictor at most stations. The latitudinal dipole in SFE during ENSO phases results from changes in snowfall frequency and extreme events. In the Pacific Northwest, above normal SFE during La Niña winters was a product of both larger contributions from extremes and more snowfall days, while below normal SFE during El Niño winters was primarily associated with a substantial reduction in extremes. Conversely, annual SFE during ENSO phases in the mountains of Arizona was more closely linked to fluctuations in snowfall days than extremes. Results indicate the importance of extreme snowfall events in shaping interannual variability in water resources and suggest that improved predictive ability may inform better water resource management now and in the coming decades.

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: Cell-centered finite volume methods are prevailing in numerical simulation of flow in porous media. However, due to the lack of cell-centered finite volume methods for mechanics, coupled flow and deformation is usually treated either by coupled finite-volume-finite element discretizations, or within a finite element setting. The former approach is unfavorable as it introduces two separate grid structures, while the latter approach loses the advantages of finite volume methods for the flow equation. Recently we proposed a cell-centered finite volume method for elasticity. Herein we explore the applicability of this novel method to provide a compatible finite volume discretization for coupled hydro-mechanic flows in porous media. We detail in particular the issue of coupling terms, and show how this is naturally handled. Furthermore, we observe how the cell-centered finite volume framework naturally allows for modeling fractured and fracturing porous media through internal boundary conditions. We support the discussion with a set of numerical examples: The convergence properties of the coupled scheme are first investigated; Secondly, we illustrate the practical applicability of the method both for fractured and heterogeneous media.

Description: The spatiotemporal distribution of Land Surface Temperature (LST) is linked to the partitioning of the coupled surface water and energy budgets. In watersheds with a strong seasonality in precipitation and vegetation cover, the temporal evolution of LST patterns are a signature of the interactions between the land surface and atmosphere. Nevertheless, few studies have sought to understand the topographical and ecohydrological controls on LST in regions of complex terrain. Numerical watershed models, tested against spatially-distributed field and remote sensing data, can aid in linking the seasonal evolution of LST to meteorology, terrain, soil and vegetation. In this study, we use a distributed hydrologic model to explore LST patterns in a semiarid mountain basin during the transition from a dry spring to the wetter North American monsoon (NAM). By accounting for vegetation greening through remotely-sensed parameters, the model reproduces LST and surface soil moisture observations derived from ground, aircraft and satellite platforms with good accuracy at individual sites and as spatial basin patterns. Distributed simulations reveal how LST varies with elevation, slope and aspect and the role played by the seasonal vegetation canopy in cooling the land surface and increasing the spatial variability in LST. As a result, LST is shown to track well with ecosystem-specific changes in vegetation cover, evapotranspiration and soil moisture during the NAM. Furthermore, vegetation greening is shown to modulate the spatial heterogeneity of LST during the NAM that should be considered in subsequent atmospheric studies in regions of complex terrain.

Description: Models of landscape evolution or hydrological processes typically depend on the accurate determination of upslope drainage area from digital elevation data, but such calculations can be very computationally demanding when applied to high-resolution topographic data. To overcome this limitation, we propose calculating drainage area in an implicit, iterative manner using linear solvers. The basis of this method is a recasting of the flow routing problem as a sparse system of linear equations, which can be solved using established computational techniques. This approach is highly parallelizable, enabling data to be spread over multiple computer processors. Good scalability is exhibited, rendering it suitable for contemporary high-performance computing architectures with many processors, such as graphics processing units (GPUs). In addition, the iterative nature of the computational algorithms we use to solve the linear system creates the possibility of accelerating the solution by providing an initial guess, making the method well suited to iterative calculations such as numerical landscape evolution models. We compare this method with a previously proposed parallel drainage area algorithm and present several examples illustrating its advantages, including a continent-scale flow routing calculation at 3 arcsecond resolution, improvements to models of fluvial sediment yield, and acceleration of drainage area calculations in a landscape evolution model. We additionally describe a modification that allows the method to be used for parallel basin delineation.

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: In this article, we discuss flows in shallow, stratified horizontal layers of two immiscible fluids. The top layer is an electrolyte which is electromagnetically driven and the bottom layer is a dielectric fluid. Using a quasi-two-dimensional approximation, which assumes a horizontal flow whose direction is independent of the vertical coordinate, we derive a generalized two-dimensional vorticity equation describing the evolution of the horizontal flow. Also, we derive an expression for the vertical profile of the horizontal velocity field. Measuring the horizontal velocity fields at the electrolyte-air and electrolyte-dielectric interfaces using particle image velocimetry, we validate the theoretical predictions of the horizontal velocity and its vertical profile for steady as well as for freely decaying Kolmogorov-like flows. Our analysis shows that by increasing the viscosity of the electrolyte relative to that of the dielectric, one may significantly improve the uniformity of the flow in the electrolyte, yielding excellent agreement between the analytical predictions and the experimental measurements.

Description: ABSTRACT Rising development costs and growing concerns over environmental impacts have led many communities to explore more diversified water management strategies. These 'portfolio'-style approaches integrate existing supply infrastructure with other options such as conservation measures or water transfers. Diversified water supply portfolios have been shown to reduce the capacity and costs required to meet demand, while also providing greater adaptability to changing hydrologic conditions. However, this additional flexibility can also cause unexpected reductions in revenue (from conservation) or increased costs (from transfers). The resulting financial instability can act as a substantial disincentive to utilities seeking to implement more innovative water management techniques. This study seeks to design portfolios that employ financial tools (e.g. contingency funds, index insurance) to reduce fluctuations in revenues and costs, allowing these strategies to achieve improved performance without sacrificing financial stability. This analysis is applied to the development of coordinated regional supply portfolios in the 'Research Triangle' region of North Carolina, an area comprising four rapidly growing municipalities. The actions of each independent utility become interconnected when shared infrastructure is utilized to enable inter-utility transfers, requiring the evaluation of regional tradeoffs in up to five performance and financial objectives. Diversified strategies introduce significant tradeoffs between achieving reliability goals and introducing burdensome variability in annual revenues and/or costs. Financial mitigation tools can mitigate the impacts of this variability, allowing for an alternative suite of improved solutions. This analysis provides a general template for utilities seeking to navigate the tradeoffs associated with more flexible, portfolio-style management approaches.

Description: Despite significant advances during the last decades, there are still many processes related to non-equilibrium flow and transport in macroporous soil that are far from completely understood. The use of X-ray for imaging time-lapse 3-D solute transport has a large potential to help advance the knowledge in this field. We visualized the transport of potassium iodide (20 g iodide l -1 H 2 O) front through a small undisturbed soil column (height 3.8 cm, diameter 6.8 cm) under steady-state hydraulic conditions using an industrial X-ray scanner. In addition, the electrical conductivity was measured in the effluent solution during the experiment. We attained a series of seventeen 3-D difference images which we related to iodide concentrations using a linear calibration relationship. The solute transport through the soil mainly took place in two cylindrical macropores, by-passing more than 90% of the bulk soil volume during the entire experiment. From these macropores the solute diffused into the surrounding soil matrix. We illustrated the properties of the investigated solute transport by comparing it to a 1-D convective-dispersive transport and by calculating the temporal evolution of the dilution index. We furthermore showed that the tracer diffusion from one of the macropores into the surrounding soil matrix could not be exactly fitted with the cylindrical diffusion equation. We believe that similar studies will help establish links between soil structure and solute transport processes and lead to improvements in models for solute transport through undisturbed soil.

Description: A semi-analytical grid-free series solution method is presented for modeling 3-D steady-state free boundary groundwater-surface water exchange in geometrically complex stratified aquifers. Continuous solutions for pressure in the subsurface are determined semi-analytically, as is the location of the water table surface. Mass balance is satisfied exactly over the entire domain except along boundaries and interfaces between layers, where errors are shown to be acceptable. The solutions are derived and demonstrated on a number of test cases and the errors are assessed and discussed. This accurate and grid-free scheme can also be a helpful tool for providing insight into lake-aquifer and stream-aquifer interactions. Here, it is used to assess the impact of lake sediment geometry and properties on lake-aquifer interactions. Various combinations of lake sediment are considered and the appropriateness of the Dupuit-Forchheimer approximation for simulating lake bottom flux distribution is investigated. In addition, the method is applied to a test problem of surface seepage flows from a complex topographic surface; this test case demonstrated the method's efficacy for simulating physically realistic domains.

Description: This article is aimed at quantifying and discussing the relative magnitude of key terms of the equation of conservation of turbulent kinetic energy (TKE) in the inter-stem space of a flow within arrays of vertical cylinders simulating plant stems of emergent and rigid vegetation. The spatial distribution of turbulent quantities and mean flow variables are influenced by two fundamental space scales, the diameter of the stems and the local stem areal number-density. Both may vary considerably since the areal distribution of plant stems in natural systems is generally not homogeneous; they are often arranged in alternating sparse and dense patches. The magnitude of the terms of the budget of TKE in the inter-stem space has seldom been quantified experimentally and is currently not well-known. This work addresses this research need. New databases, consisting of three-component LDA velocity series and two-component PIV velocity maps, obtained in carefully controlled laboratory conditions, were used to calculate the terms of the TKE budget. The physical system comprises random arrays of rigid and emergent cylinders with longitudinally varying areal number-density. It is verified that the main source of TKE is vortex shedding from individual cylinders. The rates of production and dissipation are not in equilibrium. Regions with negative production, a previously unreported feature, are identified. Turbulent transport is particularly important along the von Kármán vortex street. Convective rate of change of TKE and pressure diffusion are most relevant in the vicinity of the cylinders.

Description: A major difficulty in modeling multiphase flow in porous media is the emergence of trapped phases. Our experiments demonstrate that gas can be trapped in either single-pores, multi-pores, or in large connected networks. These large connected clusters can comprise up to 8 grain-volumes and can contain up to 50% of the whole trapped gas volume. About 85% of the gas volume is trapped by multi-pore gas clusters. This variety of possible trapped gas clusters of different shape and volume will lead to a better process understanding of bubble-mediated mass transfer. Since multi-pore gas bubbles are in contact with the solid surface through ultra-thin adsorbed water films the interfacial area between trapped gas clusters and intergranular capillary water is only about 80% of the total gas surface. We could derive a significant (R 2 = 0.98) linear relationship between the gas-water-interface and gas saturation. We found no systematic dependency of the front velocity of the invading water phase in the velocity range from 0.1 to 0.6 cm/min corresponding to capillary numbers from 2×10 -7 to 10 -6 . Our experimental results indicate that the capillary trapping mechanism is controlled by the local pore structure and local connectivity and not by thermodynamics, i.e. by the minimum of the Free Energy , at least in the considered velocity range. Consistent with this physical picture is our finding that the trapping frequency (= bubble-size distribution) reflects the pore-size distribution for the whole range of pore radii, i.e. the capillary trapping process is determined by statistics and not by thermodynamics.

Description: We investigate the direct enstrophy cascade of two-dimensional decaying turbulence in a flowing soap film channel. We use a coarse-graining approach that allows us to resolve the nonlinear dynamics and scale-coupling simultaneously in space and in scale. From our data, we verify an exact relation due to Eyink [“Local energy flux and the refined similarity hypothesis,” J. Stat. Phys.78, 335–351 (1995); Eyink “Exact results on scaling exponents in the 2D enstrophy cascade,” Phys. Rev. Lett.74, 3800–3803 (1995)] between traditional 3rd-order structure function and the enstrophy flux obtained by coarse-graining. We also present experimental evidence that enstrophy cascades to smaller (larger) scales with a 60% (40%) probability, in support of theoretical predictions by Merilees and Warn [“On energy and enstrophy exchanges in two-dimensional non-divergent flow ,” J. Fluid Mech.69, 625–630 (1975)] which appear to be valid in our flow owing to the ergodic nature of turbulence. We conjecture that their kinematic arguments break down in quasi-laminar 2D flows. We find some support for these ideas by using an Eulerian coherent structure identification technique, which allows us to determine the effect of flow topology on the enstrophy cascade. A key finding is that “centers” are inefficient at transferring enstrophy between scales, in contrast to “saddle” regions which transfer enstrophy to small scales with high efficiency.

Description: Hydrological and geochemical processes controlling the pore water chemistry in a permafrost wetland, with loam overlain by sphagnum peat, were investigated. The vertical distributions of dissolved Cl, and of pore water δ 18 O, appeared unrelated to ion freeze-out and isotope ice-water fractionation processes, respectively, dismissing solute freeze-out as a main control on the water chemistry. However, concentrations of major ions, others than Cl, generally increased with depth into the active layer. A conceptual model for water and solute movement in the active layer was derived. The model indicates upwards diffusive transport of elements, released in the loam layer by mineral weathering, to the peat layer, in which lateral advective transport dominates. Active layer pore water and water of melted core sections of permafrost were of Ca-Mg-HCO 3 type (1:1:4 stoichiometry) and were subsaturated for calcite and dolomite. The results are consistent with an annual cycling of inorganic carbon species, Ca and Mg, via cryogenic carbonate precipitation during fall freeze-up and their re-dissolution following spring thaw. Similarly, elevated Fe 2+ concentrations appear to be related to cryogenic siderite formation. Pore water in the active layer showed high partial pressures CO 2 , indicating the feasibility of bubble ebullition as a greenhouse gas emission pathway from permafrost wetlands. Elevated concentrations of geogenic trace elements (Ni, Al and As) were observed, and the controlling geochemical processes are discussed. The conceptual model for water and solute movement was applied to quantify the contribution of released trace elements to a downstream lake in the permafrost catchment.