Turbulence in Hydraulic Jumps: Experimental Measurements

A hydraulic jump is the transition from a supercritical open channel flow to a subcritical regime. It is characterised by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air-water flow measurements were performed in hydraulic jump flows for a range of inflow Froude numbers. The experiments were conducted in a large-size facility using two types of phase-detection intrusive probes: i.e., single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The present study was focused on the turbulence characteristics of hydraulic jumps with partially-developed inflow conditions. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The present results highlighted some influence of the inflow Froude number onto the air bubble entrainment process. At the largest Froude numbers, the advected air bubbles were more thoroughly dispersed vertically, and larger amount of air bubbles were detected in the turbulent shear layer. In the air-water mixing layer, the maximum void fraction and bubble count rate data showed some longitudinal decay function in the flow direction. Such trends were previously reported in the literature. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The present data suggested some longitudinal decay of the turbulence intensity. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of CHANSON (2006,2007). The results provided the integral turbulent time and length scales of the eddy structures advecting the air bubbles in the developing shear layer. The experimental data showed that the auto-correlation time scale Txx was larger than the transverse cross-correlation time scale Txz. The integral turbulence length scale Lxz was a function of the inflow conditions, of the streamwise position (x-x1)/d1 and vertical elevation y/d1. Herein the dimensionless integral turbulent length scale Lxz/d1 was closely related to the inflow depth: i.e., Lxz/d1 = 0.2 to 0.8, with Lxz increasing towards the free-surface. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. A linear relationship was found between the normalized maximum free-surface fluctuation and the inflow Froude number.

Dynamic Similarity and Scale Effects Affecting Air Bubble Entrainment in Hydraulic Jumps

In an open channel, the transition from super- to sub-critical flow is a flow singularity (the hydraulic jump) characterised by a sharp rise in free-surface elevation, strong turbulence and air entrainment in the roller. A key feature of the hydraulic jump flow is the strong free-surface aeration and air-water flow turbulence. In the present study, similar experiments were conducted with identical inflow Froude numbers Fr1 using a geometric scaling ratio of 2:1. The results of the Froude-similar experiments showed some drastic scale effects in the smaller hydraulic jumps in terms of void fraction, bubble count rate and bubble chord time distributions. Void fraction distributions implied comparatively greater detrainment at low Reynolds numbers yielding some lesser aeration of the jump roller. The dimensionless bubble count rates were significantly lower in the smaller channel, especially in the mixing layer. The bubble chord time distributions were quantitatively close in both channels, and they were not scaled according to a Froude similitude. Simply the hydraulic jump remains a fascinating two-phase flow motion that is still poorly understood.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Implications of specimen preparation and of surface contamination for the measurement of the grain boundary carbon concentration of steels using x-ray microanalysis in an UHVFESTEM

The purpose of the present investigation was to gain an understanding of the nature of the carbon contamination on the surface of standard steel transmission electron spectroscopy (TEM) specimens, the effect of exposure of a clean specimen to normal laboratory air, and the efficacy of plasma-cleaning treatments. This knowledge is a necessary prerequisite to the development of appropriate specimen preparation and/or specimen cleaning methods. X-ray photoelectron spectroscopy in combination with argon ion beam profiling was used to characterize the specimen surfaces of X65 steel and 316 stainless steel. The only clean carbon-free surface obtained was that during argon etching of the sample in the surface analysis chamber. Any exposure of a previously cleaned sample to laboratory air resulted in a rapid carbon (hydrocarbon) contamination of the sample surface and the development of surface oxidation, Plasma cleaning with subsequent exposure of the specimen to the laboratory air also resulted in a carbon-contaminated surface. This suggests that procedures of preparation of TEM specimens of steels outside an ultrahigh vacuum chamber are unlikely to result in the lowering of contamination rates on specimens to levels where measurements for carbon in the grain boundaries are possible. What is needed is a cleaning system as an integral part of the specimen insertion system into the field-emission scanning transmission electron microscope. This cleaning could be carried out by argon ion etching. Copyright (C) 2000 John Wiley & Sons, Ltd.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Influences of surface and solvent on retention of HEMA/mixture components after evaporation

Objectives: This study examined the retention of solvents within experimental HEMA/solvent primers after two conditions for solvent evaporation: from a free surface or from dentine surface. Methods: Experimental primers were prepared by mixing 35% HEMA with 65% water, methanol, ethanol or acetone (v/v). Aliquots of each primer (50 mu l) were placed on glass wells or they were applied to the surface of acid-etched dentine cubes (2 mm x 2 mm x 2 mm) (n = 5). For both conditions (i.e. from free surface or dentine cubes), change in primers mass due to solvent evaporation was gravimetrically measured for 10 min at 51% RH and 21 degrees C. The rate of solvent evaporation was calculated as a function of loss of primers mass (%) over time. Data were analysed by two-way ANOVA and Student-Newman-Keuls (p < 0.05). Results: There were significant differences between solvent retention(%) and evaporation rate (%/min) depending on the solvent present in the primer and the condition for evaporation (from free surface or dentine cubes) (p < 0.05). For both conditions, the greatest amount of retained solvent was observed for HEMA/water primer. The rate of solvent evaporation for HEMA/acetone primer was almost 2- to 10-times higher than for HEMA/water primer depending whether evaporation occurred, respectively, from a free surface or dentine cubes. The rate of solvent evaporation varied with time, being in general highest at the earliest periods. Conclusions: The rate of solvent evaporation and its retention into HEMA/solvent primers was influenced by the type of the solvent and condition allowed for their evaporation. (C) 2009 Elsevier Ltd. All rights reserved.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.

Caractéristiques diphasiques des écoulements sur les coursiers en marches d'escalier

Stepped cascades and spillways have been used for more than 3,500 years. The recent regain of interest for the stepped chute design is associated with the introduction of new construction techniques, the development of new design techniques and newer applications. Stepped chute flows are characterised by significant free-surface aeration that cannot be neglected. Two-phase flow measurements were conducted in a large-size model (h =0.10 m, α = 22o). Experimental observations demonstrate the existence of a transition flow regime for a relatively wide range of flow rates. Detailed air-water flow measurements were conducted for both skimming flows and transition flows. Skimming flows exhibit gradual variations of the air-water flow properties, whereas transition flows are characterised by rapid flow redistributions between adjacent steps.

Estudo numérico do escoamento e transporte de sedimentos em torno de pilares de pontes

Dissertação de mestrado integrado em Engenharia Civil

Bound states of 3He at the edge of a 4He drop on a cesium surface

We show that small amounts of 3He atoms, added to a 4He drop deposited on a flat cesium surface at zero temperature, populate bound states localized at the contact line. These edge states show up for drops large enough to develop well defined surface and bulk regions together with a contact line, and they are structurally different from the well-known Andreev states that appear at the free surface and at the liquid-solid interface of films. We illustrate the one-body density of 3He in a drop with 1000 4He atoms, and show that for a sufficiently large number of impurities the density profiles spread beyond the edge, coating both the curved drop surface and its flat base and eventually isolating it from the substrate.

Film-free laser forward printing of transparent and weakly absorbing liquids

A laser-based technique for printing transparent and weakly absorbing liquids is developed. Its principle of operation relies in the tight focusing of short laser pulses inside the liquid and close to its free surface, in such a way that the laser radiation is absorbed in a tiny volume around the beam waist, with practically no absorption in any other location along the beam path. If the absorbed energy overcomes the optical breakdown threshold, a cavitation bubble is generated, and its expansion results in the propulsion of a small fraction of liquid which can be collected on a substrate, leading to the printing of a microdroplet for each laser pulse. The technique does not require the preparation of the liquid in thin film form, and its forward mode of operation imposes no restriction concerning the optical properties of the substrate. These characteristics make it well suited for printing a wide variety of materials of interest in diverse applications. We demonstrate that the film-free laser forward printing technique is capable of printing microdroplets with good resolution, reproducibility and control, and analyze the influence of the main process parameter, laser pulse energy. The mechanisms of liquid printing are also investigated: time-resolved imaging provides a clear picture of the dynamics of liquid transfer which allows understanding the main features observed in the printed droplets.

Formation and distribution of point defects on a disclination line near a free nematic interface

Point defects of opposite signs can alternately nucleate on the -1/2 disclination line that forms near the free surface of a confined nematic liquid crystal. We show the existence of metastable configurations consisting of periodic repetitions of such defects. These configurations are characterized by a minimal interdefect spacing that is seen to depend on sample thickness and on an applied electric field. The time evolution of the defect distribution suggests that the defects attract at small distances and repel at large distances.

Formation of disclinations line near a free nematic interface

We have studied the nucleation and the physical properties of a -1/2 wedge disclination line near the free surface of a confined nematic liquid crystal. The position of the disclination line has been related to the material parameters (elastic constants, anchoring energy, and favored anchoring angle of the molecules at the free surface). The use of a planar model for the structure of the director field (whose predictions have been contrasted to those of a fully three-dimensional model) has allowed us to relate the experimentally observed position of the disclination line to the relevant properties of the liquid crystals. In particular, we have been able to observe the collapse of the disclination line due to a temperature-induced anchoring-angle transition, which has allowed us to rule out the presence of a real disclination line near the nematic/isotropic front in directional growth experiments.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

A High Order Finite Volume Scheme for the 2D Shallow Water Equations Including Topography

Inhalt dieser Arbeit ist ein Verfahren zur numerischen Lösung der zweidimensionalen Flachwassergleichung, welche das Fließverhalten von Gewässern, deren Oberflächenausdehnung wesentlich größer als deren Tiefe ist, modelliert. Diese Gleichung beschreibt die gravitationsbedingte zeitliche Änderung eines gegebenen Anfangszustandes bei Gewässern mit freier Oberfläche. Diese Klasse beinhaltet Probleme wie das Verhalten von Wellen an flachen Stränden oder die Bewegung einer Flutwelle in einem Fluss. Diese Beispiele zeigen deutlich die Notwendigkeit, den Einfluss von Topographie sowie die Behandlung von Nass/Trockenübergängen im Verfahren zu berücksichtigen. In der vorliegenden Dissertation wird ein, in Gebieten mit hinreichender Wasserhöhe, hochgenaues Finite-Volumen-Verfahren zur numerischen Bestimmung des zeitlichen Verlaufs der Lösung der zweidimensionalen Flachwassergleichung aus gegebenen Anfangs- und Randbedingungen auf einem unstrukturierten Gitter vorgestellt, welches in der Lage ist, den Einfluss topographischer Quellterme auf die Strömung zu berücksichtigen, sowie in sogenannten \glqq lake at rest\grqq-stationären Zuständen diesen Einfluss mit den numerischen Flüssen exakt auszubalancieren. Basis des Verfahrens ist ein Finite-Volumen-Ansatz erster Ordnung, welcher durch eine WENO Rekonstruktion unter Verwendung der Methode der kleinsten Quadrate und eine sogenannte Space Time Expansion erweitert wird mit dem Ziel, ein Verfahren beliebig hoher Ordnung zu erhalten. Die im Verfahren auftretenden Riemannprobleme werden mit dem Riemannlöser von Chinnayya, LeRoux und Seguin von 1999 gelöst, welcher die Einflüsse der Topographie auf den Strömungsverlauf mit berücksichtigt. Es wird in der Arbeit bewiesen, dass die Koeffizienten der durch das WENO-Verfahren berechneten Rekonstruktionspolynome die räumlichen Ableitungen der zu rekonstruierenden Funktion mit einem zur Verfahrensordnung passenden Genauigkeitsgrad approximieren. Ebenso wird bewiesen, dass die Koeffizienten des aus der Space Time Expansion resultierenden Polynoms die räumlichen und zeitlichen Ableitungen der Lösung des Anfangswertproblems approximieren. Darüber hinaus wird die wohlbalanciertheit des Verfahrens für beliebig hohe numerische Ordnung bewiesen. Für die Behandlung von Nass/Trockenübergangen wird eine Methode zur Ordnungsreduktion abhängig von Wasserhöhe und Zellgröße vorgeschlagen. Dies ist notwendig, um in der Rechnung negative Werte für die Wasserhöhe, welche als Folge von Oszillationen des Raum-Zeit-Polynoms auftreten können, zu vermeiden. Numerische Ergebnisse die die theoretische Verfahrensordnung bestätigen werden ebenso präsentiert wie Beispiele, welche die hervorragenden Eigenschaften des Gesamtverfahrens in der Berechnung herausfordernder Probleme demonstrieren.