On Gravity Currents over a Porous Floor in Power-Law Cross-Section Channels

2020 ◽  
Author(s):  
T. Zemach
Keyword(s):  
Cross Section ◽  
Power Law ◽  
Gravity Currents

Critical regime of gravity currents flowing in non-rectangular channels with density stratification

2018 ◽  
Vol 840 ◽  
pp. 579-612
Author(s):  
L. Chiapponi ◽  
M. Ungarish ◽  
S. Longo ◽  
V. Di Federico ◽  
F. Addona
Keyword(s):  
Cross Section ◽  
Power Law ◽  
Rectangular Cross Section ◽  
Gravity Current ◽  
Gravity Currents ◽  
Critical Conditions ◽  
Uniform Density ◽  
Ambient Fluid ◽  
Real Numbers ◽  
Positive Real

We present theoretical and experimental analyses of the critical condition where the inertial–buoyancy or viscous–buoyancy regime is preserved in a uniform-density gravity current (which propagates over a horizontal plane) of time-variable volume ${\mathcal{V}}=qt^{\unicode[STIX]{x1D6FF}}$ in a power-law cross-section (with width described by $f(z)=bz^{\unicode[STIX]{x1D6FC}}$, where $z$ is the vertical coordinate, $b$ and $q$ are positive real numbers, and $\unicode[STIX]{x1D6FC}$ and $\unicode[STIX]{x1D6FF}$ are non-negative real numbers) occupied by homogeneous or linearly stratified ambient fluid. The magnitude of the ambient stratification is represented by the parameter $S$, with $S=0$ and $S=1$ describing the homogeneous and maximum stratification cases respectively. Earlier theoretical and experimental results valid for a rectangular cross-section ($\unicode[STIX]{x1D6FC}=0$) and uniform ambient fluid are generalized here to a power-law cross-section and stratified ambient. Novel time scalings, obtained for inertial and viscous regimes, allow a derivation of the critical flow parameter $\unicode[STIX]{x1D6FF}_{c}$ and the corresponding propagation rate as $Kt^{\unicode[STIX]{x1D6FD}_{c}}$ as a function of the problem parameters. Estimates of the transition length between the inertial and viscous regimes are also derived. A series of experiments conducted in a semicircular cross-section ($\unicode[STIX]{x1D6FC}=1/2$) validate the critical values $\unicode[STIX]{x1D6FF}_{c}=2$ and $\unicode[STIX]{x1D6FF}_{c}=9/4$ for the two cases $S=0$ and $1$. The ratio between the inertial and viscous forces is determined by an effective Reynolds number proportional to $q$ at some power. The threshold value of this number, which enables a determination of the regime of the current (inertial–buoyancy or viscous–buoyancy) in critical conditions, is determined experimentally for both $S=0$ and $S=1$. We conclude that a very significant generalization of the insights and results from two-dimensional (rectangular cross-section channel) gravity currents to power-law cross-sections is available.

The propagation of gravity currents in a V-shaped triangular cross-section channel: experiments and theory

2014 ◽  
Vol 754 ◽  
pp. 232-249 ◽  
Author(s):  
Marius Ungarish ◽  
Catherine A. Mériaux ◽  
Cathy B. Kurz-Besson
Keyword(s):  
Reynolds Number ◽  
Cross Section ◽  
Viscosity Coefficient ◽  
Gravity Currents ◽  
Quantitative Agreement ◽  
Flat Bottom ◽  
Special Configuration ◽  
Theoretical Predictions ◽  
Speed Of Propagation ◽  
High Reynolds

AbstractWe investigate the motion of high-Reynolds-number gravity currents (GCs) in a horizontal channel of V-shaped cross-section combining lock-exchange experiments and a theoretical model. While all previously published experiments in V-shaped channels were performed with the special configuration of the full-depth lock, we present the first part-depth experiment results. A fixed volume of saline, that was initially of length $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}x_0$ and height $h_0$ in a lock and embedded in water of height $H_0$ in a long tank, was released from rest and the propagation was recorded over a distance of typically $ 30 x_0$. In all of the tested cases the current displays a slumping stage of constant speed $u_N$ over a significant distance $x_S$, followed by a self-similar stage up to the distance $x_V$, where transition to the viscous regime occurs. The new data and insights of this study elucidate the influence of the height ratio $H = H_0/h_0$ and of the initial Reynolds number ${\mathit{Re}}_0 = (g^{\prime }h_0)^{{{1/2}}} h_0/ \nu $, on the motion of the triangular GC; $g^{\prime }$ and $\nu $ are the reduced gravity and kinematic viscosity coefficient, respectively. We demonstrate that the speed of propagation $u_N$ scaled with $(g^{\prime } h_0)^{{{1/2}}}$ increases with $H$, while $x_S$ decreases with $H$, and $x_V \sim [{\mathit{Re}}_0(h_0/x_0)]^{{4/9}}$. The initial propagation in the triangle is 50 % more rapid than in a standard flat-bottom channel under similar conditions. Comparisons with theoretical predictions show good qualitative agreements and fair quantitative agreement; the major discrepancy is an overpredicted $u_N$, similar to that observed in the standard flat bottom case.

Non-Boussinesq gravity currents and surface waves generated by lock release in a circular-section channel: theoretical and experimental investigation

2019 ◽  
Vol 869 ◽  
pp. 610-633 ◽  
Author(s):  
L. Chiapponi ◽  
M. Ungarish ◽  
D. Petrolo ◽  
V. Di Federico ◽  
S. Longo
Keyword(s):  
Free Surface ◽  
Surface Waves ◽  
Cross Section ◽  
Cross Sections ◽  
Density Ratio ◽  
Circular Cross Section ◽  
Gravity Currents ◽  
Numerical Code ◽  
Long Distance ◽  
Front Speed

We present a combined theoretical and experimental study of lock-release inertial gravity currents (GCs) propagating in a horizontal channel of circular cross-section with open-top surface in the non-Boussinesq regime. A two-layer shallow-water (SW) model is developed for a generic shape of the cross-section with open top, and then implemented in a finite difference numerical code for the solution in a circular-cross-section channel of the type used in the experiments. The model predicts propagation with (almost) constant speed for a fairly long distance, accompanied by a depression of the ambient free open-top surface behind the front of the current. Sixteen experiments were conducted with a density ratio $r=0.587{-}0.939$ in full-depth and part-depth release conditions, measuring the front speed and the free-surface time series at four cross-sections. The channel was a circular tube 409 cm long, with a radius of 9.5 cm; the lengths of the locks were 52 and 103.5 cm. Density contrast was obtained by adding sodium chloride and dipotassium phosphate to fresh water. The theoretical values of the front speed and of the depression overestimate the experimental values, but they predict correctly their trend for varying parameters and provide reliable insights into the underlying mechanisms. In particular, we demonstrate that the circular cross-section increases the speed of propagation as compared to the standard rectangular cross-section case (for the same initial height and density ratio). The discrepancies between the SW predictions and the present experiments are of the same order of magnitude as those of previously published results for simpler systems (Boussinesq, rectangular). In addition to the depression, which is a wave bound to, and following the front of, the GC, the system also displays two kinds of free-surface waves, namely the initial bump (its amplitude is of the same order as the depression) and some short-length and low-amplitude waves in the tail of the bump. These free waves propagate with a celerity well predicted by the ‘fast’ eigenvalues of the mathematical model. Comparison is provided with the celerity of a solitary wave. It is expected that discrepancies between theory and experiments can be partly attributed to the presence of these waves. The reported insights and SW prediction method can be applied to a variety of cross-sections of practical interest (triangles, trapezoids, etc.).

A Numerical Study of Dean Instability in Non-Newtonian Fluids

2005 ◽  
Vol 128 (1) ◽  
pp. 34-41 ◽  
Author(s):  
H. Fellouah ◽  
C. Castelain ◽  
A. Ould El Moctar ◽  
H. Peerhossaini
Keyword(s):  
Aspect Ratio ◽  
Cross Section ◽  
Power Law ◽  
Numerical Study ◽  
Rectangular Cross Section ◽  
Newtonian Fluids ◽  
Dean Number ◽  
Bingham Number ◽  
Curved Channels ◽  
Power Law Index

We present a numerical study of Dean instability for non-Newtonian fluids in a laminar 180deg curved-channel flow of rectangular cross section. A methodology based on the Papanastasiou model (Papanastasiou, T. C., 1987, J. Rheol., 31(5), pp. 385–404) was developed to take into account the Bingham-type rheological behavior. After validation of the numerical methodology, simulations were carried out (using FLUENT CFD code) for Newtonian and non-Newtonian fluids in curved channels of square or rectangular cross section and for a large aspect and curvature ratios. A criterion based on the axial velocity gradient was defined to detect the instability threshold. This criterion was used to optimize the grid geometry. The effects of curvature and aspect ratio on the Dean instability are studied for all fluids, Newtonian and non-Newtonian. In particular, we show that the critical value of the Dean number decreases with increasing curvature ratio. The variation of the critical Dean number with aspect ratio is less regular. The results are compared to those for Newtonian fluids to emphasize the effect of the power-law index and the Bingham number. The onset of Dean instability is delayed with increasing power-law index. The same delay is observed in Bingham fluids when the Bingham number is increased.

Converging gravity currents over a permeable substrate

2015 ◽  
Vol 778 ◽  
pp. 669-690 ◽  
Author(s):  
Zhong Zheng ◽  
Sangwoo Shin ◽  
Howard A. Stone
Keyword(s):  
Power Law ◽  
Gravity Current ◽  
Gravity Currents ◽  
Numerical Predictions ◽  
Front Position ◽  
Dependent Position ◽  
Self Similar ◽  
Vertical Leakage ◽  
Similar Solutions ◽  
Lock Gate

We study the propagation of viscous gravity currents along a thin permeable substrate where slow vertical drainage is allowed from the boundary. In particular, we report the effect of this vertical fluid drainage on the second-kind self-similar solutions for the shape of the fluid–fluid interface in three contexts: (i) viscous axisymmetric gravity currents converging towards the centre of a cylindrical container; (ii) viscous gravity currents moving towards the origin in a horizontal Hele-Shaw channel with a power-law varying gap thickness in the horizontal direction; and (iii) viscous gravity currents propagating towards the origin of a porous medium with horizontal permeability and porosity gradients in power-law forms. For each of these cases with vertical leakage, we identify a regime diagram that characterizes whether the front reaches the origin or not; in particular, when the front does not reach the origin, we calculate the final location of the front. We have also conducted laboratory experiments with a cylindrical lock gate to generate a converging viscous gravity current where vertical fluid drainage is allowed from various perforated horizontal substrates. The time-dependent position of the propagating front is captured from the experiments, and the front position is found to agree well with the theoretical and numerical predictions when surface tension effects can be neglected.

Numerical Simulation of Laminar Flow of Yield-Power-Law Fluids in Conduits of Arbitrary Cross-Section

1993 ◽  
Vol 115 (4) ◽  
pp. 710-716 ◽  
Author(s):  
Idir Azouz ◽  
Siamack A. Shirazi ◽  
Ali Pilehvari ◽  
J. J. Azar
Keyword(s):  
Laminar Flow ◽  
Cross Section ◽  
Power Law ◽  
Annular Flow ◽  
Flow Behavior ◽  
Arbitrary Cross Section ◽  
Curvilinear Coordinates ◽  
Eccentric Annulus ◽  
Power Law Fluids ◽  
Eccentric Annuli

A numerical model has been developed to simulate laminar flow of Power-law and Yield-Power law fluids in conduits of arbitrary cross-section. The model is based on general, nonorthogonal, boundary-fitted, curvilinear coordinates, and represents a new approach to the solution of annular flow problems. The use of an effective viscosity in the governing equation of the flow allows the study of the flow behavior of any fluid for which the shear stress is a function of shear rate only. The model has been developed primarily to simulate annular flow of fluids used in drilling and completion operations of oil or gas wells. Predicted flow rates versus pressure gradient for laminar flow of Newtonian fluids in concentric and eccentric annuli, and Power-law fluids in concentric annuli compare very well with results derived from analytical expressions. Moreover, the predictions for laminar flow of Power-law and Yield-Power-law fluids in eccentric annuli are in excellent agreement with numerical and experimental data published in the literature. The model was also successfully applied to the case of laminar flow of Power-law fluids in an eccentric annulus containing a stationary bed of drilled cuttings and the results are presented herein.

Sign in / Sign up

Export Citation Format

Share Document