共查询到20条相似文献，搜索用时 31 毫秒
1.
Ashes Banerjee Srinivas Pasupuleti Mritunjay Kumar Singh G.N. Pradeep Kumar 《Acta Geophysica》2018,66(1):81-91
Complexity of the pore geometry and the random nature of flow velocity make it difficult to predict and represent post laminar flow through porous media. Present study experimentally investigates the applicability of Forchheimer and Wilkins equations for post laminar flow where Darcy’s law is invalid due to predominant inertial effect. It is observed that both porosity and media size have significant influence over the coefficients of the Forchheimer coefficients. To incorporate the effect of porosity and media size, behaviour of Forchheimer coefficients are investigated with hydraulic radius as characteristic length. An inversely proportional variation trend is found for all the present and earlier reported data. A new empirical relation between Forchheimer coefficients and hydraulic radius is obtained which can be universally applicable for all media size and porosity. Coefficients of the Wilkins equation are found to be non-deviating for different hydraulic radius in the present study and in the reported literature validating its applicability in predicting the non laminar flow through porous media. Further the Wilkins equation is modified after incorporating the correction factors for better applicability on the field. 相似文献
2.
A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW‐2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells. 相似文献
3.
The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters. 相似文献
4.
《Advances in water resources》2005,28(10):1048-1056
For steady two-dimensional free surface flow over a horizontal impervious base, the Dupuit–Forchheimer theory assumes that the vertical component of velocity is zero, even for non-zero accretion rate at the free surface. This is improved by assuming that the vertical velocity component is zero at the base, and is proportional to height above the base. This requires the piezometric head to depend linearly on the square of the height, and the two parameters in this relation can be fitted to the two boundary conditions at the free surface, to give an expression for the free surface slope in terms of accretion, free surface height, and the pressure integral. For problems in which the pressure integral is known explicitly, this first order of ordinary differential equation for the free surface height can be solved numerically. The solutions are more accurate than the Dupuit–Forchheimer expressions for the free surface, and much easier to calculate than numerical solutions to the full two-dimensional problem. Four examples are given, leading to some simple analytical approximations for quantities of interest. 相似文献
5.
In this work we propose upscaling method for nonlinear Forchheimer flow in heterogeneous porous media. The generalized Forchheimer law is considered for incompressible and slightly-compressible single-phase flows. We use recently developed analytical results (Aulisa et al., 2009) [1] and formulate the resulting system in terms of a degenerate nonlinear flow equation for the pressure with the nonlinearity depending on the pressure gradient. The coarse scale parameters for the steady state problem are determined so that the volumetric average of velocity of the flow in the domain on fine scale and on coarse scale are close. A flow-based coarsening approach is used, where the equivalent permeability tensor is first evaluated following streamline methods for linear cases, and modified in order to take into account the nonlinear effects. Compared to previous works (Garibotti and Peszynska, 2009) [2], (Durlofsky and Karimi-Fard) [3], this approach can be combined with rigorous mathematical upscaling theory for monotone operators, (Efendiev et al., 2004) [4], using our recent theoretical results (Aulisa et al., 2009) [1]. The developed upscaling algorithm for nonlinear steady state problems is effectively used for variety of heterogeneities in the domain of computation. Direct numerical computations for average velocity and productivity index justify the usage of the coarse scale parameters obtained for the special steady state case in the fully transient problem. For nonlinear case analytical upscaling formulas in stratified domain are obtained. Numerical results were compared to these analytical formulas and proved to be highly accurate. 相似文献
6.
7.
Optimization of groundwater and other subsurface resources requires analysis of multiple‐well systems. The usual modeling approach is to apply a linear flow equation (e.g., Darcy's law in confined aquifers). In such conditions, the composite response of a system of wells can be determined by summating responses of the individual wells (the principle of superposition). However, if the flow velocity increases, the nonlinear losses become important in the near‐well region and the principle of superposition is no longer valid. This article presents an alternative method for applying analytical solutions of non‐Darcy flow for a single‐ to multiple‐well systems. The method focuses on the response of the central injection well located in an array of equally spaced wells, as it is the well that exhibits the highest pressure change within the system. This critical well can be represented as a single well situated in the center of a closed square domain, the width of which is equal to the well spacing. It is hypothesized that a single well situated in a circular region of the equivalent plan area adequately represents such a system. A test case is presented and compared with a finite‐difference solution for the original problem, assuming that the flow is governed by the nonlinear Forchheimer equation. 相似文献
8.
Melina G. Sidiropoulou Konstadinos N. Moutsopoulos Vassilios A. Tsihrintzis 《水文研究》2007,21(4):534-554
This study focuses on the determination of the Forchheimer equation coefficients a and b for non‐Darcian flow in porous media. Original theoretical equations are evaluated and empirical relations are proposed based on an investigation of available data in the literature. The validity of these equations is checked using existing experimental data, and their accuracy versus existing approaches is studied. On the basis of this analysis, some insight into the physical background of the phenomenon is also provided. The dependence of the coefficients a and b on the Reynolds number is also detected, and potential future research areas, e.g. investigation of inertial effects for consolidated porous media, are pointed out. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
《Continental Shelf Research》2006,26(17-18):2029-2049
A field experiment was carried out in Massachusetts Bay in August 1998 to assess the role of large-amplitude internal waves (LIWs) in resuspending bottom sediments. The field experiment consisted of a four-element moored array extending from just west of Stellwagen Bank (90-m water depth) across Stellwagen Basin (85- and 50-m water depth) to the coast (24-m water depth). The LIWs were observed in packets of 5–10 waves, had periods of 5–10 min and wavelengths of 200–400 m, and caused downward excursions of the thermocline of as much as 30 m. At the 85-m site, the current measured 1 m above bottom (mab) typically increased from near 0 to 0.2 m/s offshore in a few minutes upon arrival of the LIWs. At the 50-m site, the near-bottom offshore flow measured 6 mab increased from about 0.1 to 0.4–0.6 m/s upon arrival of the LIWs and remained offshore in the bottom layer for 1–2 h. The near-bottom currents associated with the LIWs, in concert with the tidal currents, were directed offshore and sufficient to resuspend the bottom sediments at both the 50- and 85-m sites. When LIWs are present, they may resuspend sediments for as long as 5 hours each tidal cycle as they travel westward across Stellwagen Basin. At 85-m water depth, resuspension associated with LIWs is estimated to occur for about 0.4 days each summer, about the same amount of time as caused by surface waves. 相似文献
10.
Observations from two SOund Detection And Ranging (SODAR) units, a 10 m micrometeorological tower and five Automated Surface Observing Stations (ASOS) were examined during several synoptic scale flow regimes over New York City after the World Trade Center disaster on September 11, 2001. An ARPS model numerical simulation was conducted to explore the complex mesoscale boundary layer structure over New York City. The numerical investigation examined the urban heat island, urban roughness effect and sea breeze structure over the New York City region. Estimated roughness lengths varied from 0.7 m with flow from the water to 4 m with flow through Manhattan. A nighttime mixed layer was observed over lower Manhattan, indicating the existence of an urban heat island. The ARPS model simulated a sea-breeze front moving through lower Manhattan during the study period consistent with the observations from the SODARs and the 10-m tower observations. Wind simulations showed a slowing and cyclonic turning of the 10-m air flow as the air moved over New York City from the ocean. Vertical profiles of simulated TKE and wind speeds showed a maximum in TKE over lower Manhattan during nighttime conditions. It appears that this TKE maximum is directly related to the influences of the urban heat island. 相似文献
11.
Vic Kelson 《Ground water》2012,50(6):918-926
Groundwater flow models are commonly used to design new wells and wellfields. As the spatial scale of the problem is large and much local‐scale detail is not needed, modelers often utilize two‐dimensional (2D) or quasi three‐dimensional models based on the Dupuit‐Forchheimer assumption. Dupuit models offer a robust set of tools for simulating regional groundwater flow including interactions with surface waters, the potential for well interference, and varying aquifer properties and recharge rates. However, given an assumed operating water level or drawdown at a well screen, Dupuit models systematically overpredict well yields. For design purposes, this discrepancy is unacceptable, and a method for predicting accurate well yields is needed. While published methods exist for vertical wells, little guidance is available for predicting yields in horizontal screens or collector wells. In plan view, a horizontal screen has a linear geometry, and will likely extend over several neighboring cells that may not align with rows or columns in a numerical model. Furthermore, the model must account for the effects of converging three‐dimensional (3D) flow to the well screens and hydraulic interference among the well screens; these all depend on the design of a specific well. This paper presents a new method for simulating the yield of angled or horizontal well screens in numerical groundwater flow models, specifically using the USGS code MODFLOW. The new method is compared to a detailed, 3D analytic element model of a collector well in a field of uniform flow. 相似文献
12.
L.K. Coachman 《Continental Shelf Research》1982,1(1):1-14
The broad (～500 km) southeastern Bering Sea continental shelf contains three fronts; outer (shelf break, ～170-m depth), middle (～100-m depth), and inner (～50-m depth). The shelf break and inner fronts appear to be analogous to similar fronts reported from other mid-latitude continental shelves; extensively studied examples are from the mid-Atlantic bight, off Nova Scotia, and around the British Isles. The middle front may have counterparts on the broad North Sea and East China Sea shelves.One-month current and temperature records from either side of the middle front, ～150 km landward from the shelf break, showed convergence in the layers deeper than 30 m in both the cross-shelf flow field and heat flux. The convergence was ～3 cm s?1, so an average upwelling at ～1 × 10?3cm s?1 and divergence in the surface layer were required to maintain continuity. Variations in the degree of convergence arose primarily from 1 to 5-day fluctuations in sub-tidal flow across the outer shelf domain seaward of the front.Diffusive landward heat flux was dominated by tidal scales. Horizontal eddy conductivities describing the flux were ～1 ? 106 on the landward side and ～5 × 106cm2 s?1 on the seaward side, and were less in the layers above the bottom layer. Advective flux by the mean flow was the same order as diffusive flux, but landward in the bottom layer and seaward in the mid-water column layers, in agreement with deductions from water mass analyses. Frontal effects reduced the net cross-shelf heat flux beneath ～30 m by about 50%. The observation of a flow convergence in the middle of a broad, flat continental shelf poses an important question of dynamics. 相似文献
13.
The nonhorizontal‐model‐layer (NHML) grid system is more accurate than the horizontal‐model‐layer grid system to describe groundwater flow in an unconfined sloping aquifer on the basis of MODFLOW‐2000. However, the finite‐difference scheme of NHML was based on the Dupuit‐Forchheimer assumption that the streamlines were horizontal, which was acceptable for slope less than 0.10. In this study, we presented a new finite‐difference scheme of NHML based on the Boussinesq assumption and developed a new package SLOPE which was incorporated into MODFLOW‐2000 to become the MODFLOW‐SP model. The accuracy of MODFLOW‐SP was tested against solution of Mac Cormack (1969). The differences between the solutions of MODFLOW‐2000 and MODFLOW‐SP were nearly negligible when the slope was less than 0.27, and they were noticeable during the transient flow stage and vanished in steady state when the slope increased above 0.27. We established a model considering the vertical flow using COMSOL Multiphysics to test the robustness of constrains used in MODFLOW‐SP. The results showed that streamlines quickly became parallel with the aquifer base except in the narrow regions near the boundaries when the initial flow was not parallel to the aquifer base. MODFLOW‐SP can be used to predict the hydraulic head of an unconfined aquifer along the profile perpendicular to the aquifer base when the slope was smaller than 0.50. The errors associated with constrains used in MODFLOW‐SP were small but noticeable when the slope increased to 0.75, and became significant for the slope of 1.0. 相似文献
14.
Frederick A. Cook Edward R. Decker Scott B. Smithson 《Earth and Planetary Science Letters》1978,40(3):316-326
A positive gravity anomaly, anomalous seismic structure, and abnormally high reduced heat flow occur in a possible extension of the Rio Grande rift in southern New Mexico. The reduced heat flux is explained by the transient effects of basaltic intrusions with dimensions based on interpretations from gravity and seismic data. Ages of volcanic activity in the rift area are used to determine times of thermal events. The effect of latent heat and sharp boundaries of the reduced heat flow anomaly also restrict modeling. The thermal model consists of three intrusions at the base of the crust at 35 m.y., 10 m.y. and 5 m.y. The initial temperatures of these intrusions would differ from those in the surrounding crust and upper mantle by about +800°C, +800°C, and +600°C, respectively. The 35-m.y. intrusion is consistent with Oligocene volcanism and tectonic activity in southwestern New Mexico. This mass also contributes to the Bouguer gravity high in the rift area, as do the intrusions at 10 and 5 m.y. The intrusions at 10 and 5 m.y. explain the excess reduced flux, with the 5-m.y. pulse accounting for most of the curvature in the sharp boundaries of the anomaly. The idealized model is consistent with late Tertiary andesites, the periods of little extrusive activity, and Pliocene-Quaternary basalts in the rift area. 相似文献
15.
SHEN Zhi-Gang 《国际泥沙研究》1999,(2)
IINTRODUCTIONRiversareanimportantcomponentofournatUralworld.ForthePurposeoffloodprotection,inlandshippingandlandutilization,manyrivershavebeenregulatedornarrowedinaschematic,monotonousmanner.Theriverbedandbanksaresmoothedinordertoreducetheflowresistanceandtoincreasetheflowcapacity.Thesetraditionalrivertrainingmeasures,whichcontradictnlderalevolutionlawofrivers,exertharmfuleffectsontheriverecosystemand"thealenvironmentalongtherivers.Theprimaryequilibriumofthebedloadtransportofriversaredam… 相似文献
16.
Flume experiments were conducted in a 6-m flume to determine the role of turbulence in the scour of pools. Paired results from constricted-flow experiments with and without a wake zone formed behind obstructions to flow show that pools are deeper and shorter when vertical free-shear layers are present. Although non-streamlined obstructions initially present more resistance to flow, channel-bed scour develops a pool morphology that lowers mechanical energy losses to levels below those in pools with streamlined obstructions. Scour primarily enlarges the cross-sectional area in the constricted section. Feedback between pool geometry and localized turbulence production play a secondary role in total scour but still exert a major influence on final energy slopes. The experiment shows that pools with wake zones evolve to reduce longitudinal energy expenditure with an associated reduction in total turbulence production associated with the obstruction. 相似文献
17.
The primary goals for most ground water capture systems (i.e., pump-and-treat systems) are that (1) all contaminants within zones of interest will eventually be captured and (2) the extraction and reinjection wells are best located and operated at optimal flow rates, creating hydraulically efficient flow systems. A new tool, MODular ALLocation (MODALL), is presented to aid in the design and assessment of capture systems. MODALL uses the MODFLOW-calculated cell-by-cell flow terms to evaluate internodal flow balances to determine the percentage of flow in each cell which has either originated from a given source(s) or flows to a specified sink(s). Output from MODALL can be easily displayed in isopleths of "capture fraction" (CF) to indicate the certainty or strength of capture in various areas. MODALL results are compared to the results from an analytical solution, a pathline analysis using MODPATH, and solute transport simulation with MT3DMS. A brief case study is also presented where MODALL is used to optimize an existing pump-and-treat system to more effectively and more efficiently contain a 5000-m long plume. 相似文献
18.
A model for the air–sea interface, based on the coupled pair of similarity relations for “aerodynamically” rough flow in both
fluids, is presented, which is applied to fetch-limited and high wind speed conditions which occur, for example, in hurricanes.
It is shown that the specification of the maximum 10-m drag coefficient and the 10-m wind speed and the peak wave speed at
which it occurs are sufficient to uniquely determine the drag law, which asymptotes at low wind speeds to a Charnock constant
similar to that for the fully developed wind wave sea and is almost independent of the peak wave speed at the maximum in drag
coefficient. A feature of the drag law is that it is of Charnock form, almost independent of the wave age, consistent with
the transfer of momentum to the wave spectrum being due to the smaller rather than the dominant wavelengths. The analysis
is also applied to a variable sea state in which either the surface wind or the surface Stokes drift vary, but the peak wave
speed is kept constant. The corresponding variability in the Charnock constant is in general accord with observations. 相似文献
19.
Hydraulic conductivity (K) and specific storage (S(s)) are required parameters when designing transient groundwater flow models. The purpose of this study was to evaluate the ability of commonly used hydrogeologic characterization approaches to accurately delineate the distribution of hydraulic properties in a highly heterogeneous glaciofluvial deposit. The metric used to compare the various approaches was the prediction of drawdown responses from three separate pumping tests. The study was conducted at a field site, where a 15 m × 15 m area was instrumented with four 18-m deep Continuous Multichannel Tubing (CMT) wells. Each CMT well contained seven 17 cm × 1.9 cm monitoring ports equally spaced every 2 m down each CMT system. An 18-m deep pumping well with eight separate 1-m long screens spaced every 2 m was also placed in the center of the square pattern. In each of these boreholes, cores were collected and characterized using the Unified Soil Classification System, grain size analysis, and permeameter tests. To date, 471 K estimates have been obtained through permeameter analyses and 270 K estimates from empirical relationships. Geostatistical analysis of the small-scale K data yielded strongly heterogeneous K fields in three-dimensions. Additional K estimates were obtained through slug tests in 28 ports of the four CMT wells. Several pumping tests were conducted using the multiscreen and CMT wells to obtain larger scale estimates of both K and S(s). The various K and S(s) estimates were then quantitatively evaluated by simulating transient drawdown data from three pumping tests using a 3D forward numerical model constructed using HydroGeoSphere (Therrien et al. 2005). Results showed that, while drawdown predictions generally improved as more complexity was introduced into the model, the ability to make accurate drawdown predictions at all CMT ports was inconsistent. 相似文献
20.
Paul Morgan 《Earth and Planetary Science Letters》1975,26(2):253-262
A standard core analysis technique has been modified to estimate porosities from measurements on rock fragments. For the range of rocks tested, chip-determined fractional porosities were within ±0.025 of the values measured on solid-core samples. This has enabled thermal conductivity measurements on rock fragments to be corrected for the effect of porosity, yielding agreement with conductivity determinations on solid core generally to better than ± 10%. The application of this is illustrated by the determination of heat flow in a 300-m borehole in western Cyprus (latitude 34°54′N, longitude 32°34′E, elevation 82 m). A decrease in temperature gradient with depth is almost completely compensated for by increasing thermal conductivity, and the best value for heat flow at this site is 23 ± 4 mW m?2. 相似文献