Modeling Flow into Horizontal Wells in a Dupuit-Forchheimer Model

Horizontal wells or radial collector wells are used in shallow aquifers to enhance water withdrawal rates. Groundwater flow patterns near these wells are three-dimensional (3D), but difficult to represent in a 3D numerical model because of the high degree of grid refinement needed. However, for the purpose of designing water withdrawal systems, it is sufficient to obtain the correct production rate of these wells for a given drawdown. We developed a Cauchy boundary condition along a horizontal well in a Dupuit-Forchheimer model. Such a steady-state 2D model is not only useful for predicting groundwater withdrawal rates but also for capture zone delineation in the context of source water protection. A comparison of our Dupuit-Forchheimer model for a radial collector well with a 3D model yields a nearly exact production rate. Particular attention is given to horizontal wells that extend underneath a river. A comparison of our approach with a 3D solution for this case yields satisfactory results, at least for moderate-to-large river bottom resistances.     

Velocity Distribution and the Moments of Turbulent Flow over a Sand-Gravel Mixture Bed

Water Resources - The current study performed the laboratory flume to examine the velocity distribution and the velocity moments of turbulent flow over a sand-gravel mixture bed. The 3D...     

Velocity Distribution and the Moments of Turbulent Flow over a Sand-Gravel Mixture Bed

Water Resources - The current study performed the laboratory flume to examine the velocity distribution and the velocity moments of turbulent flow over a sand-gravel mixture bed. The 3D...     

Experimental and Numerical Flow Simulation over Weirs

Water Resources - This paper reports a series of experimental studies done to simulate the flow behavior over crump and ogee type of weirs. The transition of subcritical to supercritical flow, as...     

Experimental and Numerical Investigation of the Flow Hydraulic in Gradual Transition Open Channels

Channel expansions are common in both natural and artificial open channels. With increasing cross-sectional dimensions in an expansion, the flow decelerates. Due to separation of flow and subsequent eddy formation, a significant head loss is occurred along the transition. This study presents the results of experimental investigations on subcritical flow along the expansive transition of rectangular to trapezoidal channels. Also, a numerical simulation was developed using the finite volume method with Reynolds Stress turbulent model. Water surface profiles and velocity distributions of flow through the transition were measured experimentally and compared with the numerical results. Also, hydraulic efficiency of the transition and coefficient of energy head loss were calculated. The results show that with increasing the upstream Froude number, hydraulic efficiency of the transition and coefficient of energy head loss are decreased and increased, respectively. The results also showed the ability of numerical simulation for simulating the flow separation zones and bed shear stress along the transition for different inlet discharges and inflow Froude numbers.     

Flow structure over a wavy bed with vegetation cover

On the basis of experiments carried out in flume with a wavy bed with vegetation cover, flow velocity, turbulence intensities and Reynolds stress distributions are investigated. The wavy bed was similar to dune in this study. The fixed artificial dunes were constructed over the bed and artificial vegetation put over them in a laboratory flume. An Acoustic Doppler Velocimeter and spatially-averaged method were applied to determine turbulent flow components and shear velocity. Results were compared with a gravel bedform. It was observed that vegetation cover influences considerably the flow structure and displays clearly the flow separation and reattachment point. The law of the wall was not valid within the vegetation cover, but it was fitted well to the zone above the vegetation cover within the inner layer. For a wavy bed having the same dimensions, shear velocity and friction factor over vegetation cover are 1.7 and 2.6 times of those for the gravel bedform, respectively. The results of laboratory study were compared with those of river study.     

Numerical Generation of Flow Nets — The FLOWNS Model

Abstract. A simple-to-use computer program (FLOWNS) has been developed for generating flow nets for any saturated rectangular domain with any combination of constant head or constant flux (including zero) boundary conditions. The program approximates with discrete values the continuous distributions of potential and stream function using finite-difference approximations of Laplace's equation. The hydraulic conductivity distribution may be anisotropic and/or heterogeneous. A contouring program is required to generate the final stream and equipotential lines.     

Flow dynamics over a foredune scarp

The nature of wind flow over a small, 0.6 m high foredune scarp is investigated on the Sir Richard Peninsula, South Australia during a variety of incident wind directions and speeds. The study provides additional supporting evidence that the presence of the scarp and the dune exerts a strong influence on a landwards trending reduction in wind velocity and an increase in turbulence, with the greatest area of turbulence occurring near and at the foredune scarp base. For an incident oblique wind, an alongshore helicoidal flow is formed within a separation region along the scarp basal region. In this region, the coefficient of variation (CV) of wind speed is high and displays significant fluctuations. The flow at the scarp crest is compressed, streamlined and accelerated, turbulence is suppressed, and local jets may occur depending on the incident wind approach angle. Jets are more likely where the incident flow is perpendicular or nearly so. A flow separation region does not develop downwind of the scarp crest where the morphology of the foredune stoss slope downwind of the scarp is more convex (as in this case) rather than relatively flat, and possibly due to the presence of vegetation at the scarp crest. A tentative model of the flow regions developed across a backshore–scarp–foredune region during oblique incident flow is provided. © 2018 John Wiley & Sons, Ltd.     

Groundwater Table Formation in the Coastal Zone over a Bed with Arbitrary Shape

The impact of tidal wave on a permeable shore consisting of porous material results in groundwater table lying above the average sea level. This effect, stemming from the nonlinearity of the process of seawater flow in the soil, is referred to as pumping effect.     

Open Channel Waves that Form after the Removal of a Screen before Uneven Bed

The results of an experimental study of the gravity waves that form in an open channel at the decomposition of an initial discontinuity in the free surface level above an uneven bed in the form of a shelf are discussed. Self-similar solutions obtained in the context of the first approximation of the shallow-water theory are tested. It is shown that more complex mathematical models must be used to adequately describe the waves of level drop in the upper pool and undular hydraulic jumps.     

Prediction of Flow around a Vertical Circular Pier in a Discordant Bed Channel Confluence

Water Resources - Understanding the flow dynamics around the hydraulic structure caused by complex junction flow is still a major challenge. In order to gain a better insight of impact of junction...     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Numerical simulation of a mesolow over a Gulf Stream filament

A three-dimensional mesoscale numerical model is used to investigate mesoscale circulation over a Gulf Stream filament. Two numerical experiments are performed with different initial uniform ambient wind speeds (U=0.1 m s^–1, 3.5 m s^–1 and 7 m s^–1) for a typical winter day. It is found that for both low and moderate winds, a closed mesoscale circulation forms over the Gulf Stream filament. When the Gulf Stream filament was removed, the model did not predict a mesoscale circulation. The modeled circulation over the filament is in agreement with the observations, suggesting that the atmospheric circulations over the filaments may be an important mechanism in the U.S. East Coast cyclogenesis.     

Numerical modeling of EM fields over local anomalies with a vertical axis of symmetry

A method for numerical modeling of quasi-stationary electromagnetic fields in axially symmetric media is proposed. It is based on the direct finite-element method and the use of special basis functions. Assuming cylindrical coordinates r, φ, z, the three-dimensional solution is presented as a superposition of fields with the azimuthal dependence exp(inφ). For each case this results in a system of two equations of elliptical type in two scalar functions in the (r, z) plane.

The discretization leads to the conservative nine-points difference scheme. The system of linear equations is solved by means of the LU-decomposition technique, the band structure of the matrix being taken into account.

The program is tested using analytical results (DC asymptote) for a near-surface inhomogeneity. Comparison also is made with 2D results (H-polarization) for the model of a local well conducting inclusion in a three-layered Earth.