USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Fluids Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Investigating the mixing efficiencies of liquid-to-liquid chemical injection manifolds for aquatic invasive species management

Aquatic invasive species (AIS) have spread throughout the United States via major rivers and tributaries. Locks and dams positioned along affected waterways, specifically lock chambers, are being evaluated as potential management sites to prevent further expansion into new areas. Recent research has shown that infusion of chemicals (e.g., carbon dioxide) into water can block or kill several invasive organisms and could be a viable option at navigational structures such as lock chambers because chemical infusion would not interfere with vessel passage or lock operation. Chemical treatments near lock structures will require large-scale fluid-mechanic systems and significant energy. Mixing must extend to all stagnation regions within a lock structure to prevent the passage of an invasive fish. This work describes the performance of both wall- and floor-based CO 2 -infused-water to water injection manifolds targeted for lock structures in terms of mixing time, mixing homogeneity, injection efficiency, and operational power requirements. Both systems have strengths and weaknesses so selection recommendations are given for applications such as open systems and closed systems.

Journal of Fluids Engineering

Groundwater

Introduction Groundwater represents the terrestrial subsurface component of the hydrologic cycle. As such, groundwater is generally in motion, moving from elevated areas of recharge to lower areas of discharge. Groundwater usually moves in accordance with Darcy’s law (Dalmont, Paris: Les Fontaines Publiques de la Ville de Dijon, 1856). Groundwater residence times can be under a day in small upland catchments to over a million years in subcontinental-sized desert basins. The broadest definition of groundwater includes water in the unsaturated zone, considered briefly here. Water chemically bound to minerals, as in gypsum (CaSO 4 • 2H 2 O) or hydrated clays, cannot flow in response to gradients in total hydraulic head (pressure head plus elevation head); such water is thus usually excluded from consideration as groundwater. In 1940, M. King Hubbert showed Darcy’s law to be a special case of thermodynamically based potential field equations governing fluid motion, thereby establishing groundwater hydraulics as a rigorous engineering science ( Journal of Geology 48, pp. 785–944). The development of computer-enabled numerical methods for solving the field equations with real-world approximating geometries and boundary conditions in the mid-1960s ushered in the era of digital groundwater modeling. An estimated 30 percent of global fresh water is groundwater, compared to 0.3 percent that is surface water, 0.04 percent atmospheric water, and 70 percent that exists as ice, including permafrost ( Shiklomanov and Rodda 2004 , cited under Groundwater Occurrence ). Groundwater thus constitutes the vast majority—over 98 percent—of the unfrozen fresh-water resources of the planet, excluding surface-water reservoirs. Environmental dimensions of groundwater are equally large, receiving attention on multiple disciplinary fronts. Riparian, streambed, and spring-pool habitats can be sensitively dependent on the amount and quality of groundwater inputs that modulate temperature and solutes, including nutrients and dissolved oxygen. Groundwater withdrawals can negatively impact riparian habitats by depriving ecosystems of adequate fresh water and fragmenting communities when streams go dry. Biochemical reactions in shallow groundwater can remove anthropogenically elevated nitrogen compounds and reduce—but only to a point—the greening of waterways and shorelines with periphyton and harmful algal blooms. Groundwater extraction for beneficial use is increasingly limited by water-quality constraints imposed by naturally occurring and introduced substances. Overdrafting can cause land-surface subsidence, damaging buildings and roads and disrupting canals, sewers, and other gravity-flow conveyances. Increases in groundwater levels can cause soil salinization in dry regions and erosive sapping and flooding in wet regions. Coastal saltwater intrusion, groundwater flooding, salinization associated with groundwater-irrigated agriculture, induced seismicity from injected wastes, and the detrimental impacts of groundwater depletion are among the major environmental challenges of our time.

Book chapter

1D numerical model of muddy subaqueous and subaerial debris flows

A 1D numerical model of the downslope flow and deposition of muddy subaerial and subaqueous debris flows is presented. The model incorporates the Herschel-Bulkley and bilinear rheologies of viscoplastic fluid. The more familiar Bingham model is integrated into the Herschel-Bulkley rheological model. The conservation equations of mass and momentum of single-phase laminar debris flow are layer-integrated using the slender flow approximation. They are then expressed in a Lagrangian framework and solved numerically using an explicit finite difference scheme. Starting from a given initial shape, a debris flow is allowed to collapse and propagate over a specified topography. Comparison between the model predictions and laboratory experiments shows reasonable agreement. The model is used to study the effect of the ambient fluid density, initial shape of the failed mass, and rheological model on the simulated propagation of the front and runout characteristics of muddy debris flows. It is found that initial failure shape influence the front velocity but has little bearing on the final deposit shape. In the Bingham model, the excess of shear stress above the yield strength is proportional to the strain rate to the first power. This exponent is free to vary in the Herschel-Bulkley model. When it is set at a value lower than unity, the resulting final deposits are thicker and shorter than in the case of the Bingham rheology. The final deposit resulting from the bilinear model is longer and thinner than that from the Bingham model due to the fact that the debris flow is allowed to act as a Newtonian fluid at low shear rate in the bilinear model.

Journal of Hydraulic Engineering

Verification of vertically rotating flume using non-newtonian fluids

Three tests on non-Newtonian fluids were used to verify the use of a vertically rotating flume (VRF) for the study of the rheological properties of debris flow. The VRF is described and a procedure for the analysis of results of tests made with the VRF is presented. The major advantages of the VRF are a flow field consistent with that found in nature, a large particle-diameter threshold, inexpensive operation, and verification using several different materials; the major limitations are a lack of temperature control and a certain error incurred from the use of the Bingham plastic model to describe a more complex phenomenon. Because the VRF has been verified with non-Newtonian fluids as well as Newtonian fluids, it can be used to measure the rheological properties of coarse-grained debris-flow materials.

Journal of Hydraulic Engineering

Management of fluid mud in estuaries, bays, and lakes. II: Measurement, modeling, and management

Techniques for measurement, modeling, and management of fluid mud are available, but research is needed to improve them. Fluid mud can be difficult to detect, measure, or sample, which has led to new instruments and new ways of using existing instruments. Multifrequency acoustic fathometers sense neither density nor viscosity and are, therefore, unreliable in measuring fluid mud. Nuclear density probes, towed sleds, seismic, and drop probes equipped with density meters offer the potential for accurate measurements. Numerical modeling of fluid mud requires solving governing equations for flow velocity, density, pressure, salinity, water surface, plus sediment submodels. A number of such models exist in one-, two-, and three-dimensional form, but they rely on empirical relationships that require substantial site-specific validation to observations. Management of fluid mud techniques can be classified as those that accomplish: Source control, formation control, and removal. Nautical depth, a fourth category, defines the channel bottom as a specific fluid mud density or alternative parameter as safe for navigation. Source control includes watershed management measures to keep fine sediment out of waterways and in-water measures such as structures and traps. Formation control methods include streamlined channels and structures plus other measures to reduce flocculation and structures that train currents. Removal methods include the traditional dredging and transport of dredged material plus agitation that contributes to formation control and/or nautical depth. Conditioning of fluid mud by dredging and aerating offers the possibility of improved navigability. Two examples—the Atchafalaya Bar Channel and Savannah Harbor—illustrate the use of measurements and management of fluid mud.

Journal of Hydraulic Engineering

Evaluation of gas production potential from gas hydrate deposits in National Petroleum Reserve Alaska using numerical simulations

An evaluation of the gas production potential of Sunlight Peak gas hydrate accumulation in the eastern portion of the National Petroleum Reserve Alaska (NPRA) of Alaska North Slope (ANS) is conducted using numerical simulations, as part of the U.S. Geological Survey (USGS) gas hydrate Life Cycle Assessment program. A field scale reservoir model for Sunlight Peak is developed using Advanced Processes & Thermal Reservoir Simulator (STARS) that approximates the production design and response of this gas hydrate field. The reservoir characterization is based on available structural maps and the seismic-derived hydrate saturation map of the study region. A 3D reservoir model, with heterogeneous distribution of the reservoir properties (such as porosity, permeability and vertical hydrate saturation), is developed by correlating the data from the Mount Elbert well logs. Production simulations showed that the Sunlight Peak prospect has the potential of producing 1.53 × 10 9 ST m 3 of gas in 30 years by depressurization with a peak production rate of around 19.4 × 10 4 ST m 3 /day through a single horizontal well. To determine the effect of uncertainty in reservoir properties on the gas production, an uncertainty analysis is carried out. It is observed that for the range of data considered, the overall cumulative production from the Sunlight Peak will always be within the range of ±4.6% error from the overall mean value of 1.43 × 10 9 ST m 3 . A sensitivity analysis study showed that the proximity of the reservoir from the base of permafrost and the base of hydrate stability zone (BHSZ) has significant effect on gas production rates. The gas production rates decrease with the increase in the depth of the permafrost and the depth of BHSZ. From the overall analysis of the results it is concluded that Sunlight Peak gas hydrate accumulation behaves differently than other Class III reservoirs (Class III reservoirs are composed of a single layer of hydrate with no underlying zone of mobile fluids) due to its smaller thickness and high angle of dip.

Alaska

Fines classification based on sensitivity to pore-fluid chemistry

The 75-μm particle size is used to discriminate between fine and coarse grains. Further analysis of fine grains is typically based on the plasticity chart. Whereas pore-fluid-chemistry-dependent soil response is a salient and distinguishing characteristic of fine grains, pore-fluid chemistry is not addressed in current classification systems. Liquid limits obtained with electrically contrasting pore fluids (deionized water, 2-M NaCl brine, and kerosene) are combined to define the soil “electrical sensitivity.” Liquid limit and electrical sensitivity can be effectively used to classify fine grains according to their fluid-soil response into no-, low-, intermediate-, or high-plasticity fine grains of low, intermediate, or high electrical sensitivity. The proposed methodology benefits from the accumulated experience with liquid limit in the field and addresses the needs of a broader range of geotechnical engineering problems.

Journal of Geotechnical and Geoenvironmental Engin

Solubility of NaCl in aqueous electrolyte solutions from 10 to 100°C

The solubilities of NaCl in aqueous KCl, MgCl2, CaCl2, and mixed CaCl2-KCl solutions have been determined from 10 to 100??C. The data were fit to an equation, and the equation was used to calculate values of the change in solubility of NaCl, ???[NaCl]/???T. These values are required for calculations of the rate of migration of fluids in a thermal gradient in rock salt. The data obtained here indicate that the values of ???[NaCl]/???T are 36-73% greater for solutions containing divalent ions than for the NaCl-H2O system.

Journal of Chemical and Engineering Data

Well predictive performance of play-wide and Subarea Random Forest models for Bakken productivity

In recent years, geologists and petroleum engineers have struggled to clearly identify the mechanisms that drive productivity in horizontal, hydraulically-fractured oil wells producing from the middle member of the Bakken formation. This paper fills a gap in the literature by showing how this play’s heterogeneity affects factors that drive well productivity. It is important because understanding the relative strength of productivity drivers and how predictors vary spatially facilitates best-practices for well site selection and well completion design. The paper describes an application of the Random Forest (RF) machine learning technique to identify these mechanisms and to evaluate their importance across 9 subareas of the North Dakota portion of the Bakken play. The study examined productivity of 7311 wells initiating production from 2010 through 2017. Well productivity varied considerably across the 9 subareas within the play, so it was not surprising that the dominant predictors, the initial 180-day water cut and the 30-day initial gas production, vary spatially to mirror local conditions that strongly affect well productivity. The relative importance of well completion predictor variables, that is, the numbers of fractures stages per well, volume of injected proppant per stage, volume of injected fluids per stage, and lateral length, varied considerably across the subareas. Statistical permutation tests are presented that generally confirm the importance rankings. Subarea Random Forest models explained from 50 percent to 82 percent of the variation in productivity test samples while the play-wide model explained 73 percent of the test sample well productivity. Weakness in the predictive ability of the Random Forest models are traced to the limited variability in the training data. Implications of the empirical findings regarding the Bakken play for operators and for research and government institutions are discussed in the concluding section.

Montana, North Dakota, South Dakota

Granular-flow rheology: Role of shear-rate number in transition regime

This paper examines the rationale behind the semiempirical formulation of a generalized viscoplastic fluid (GVF) model in the light of the Reiner-Rivlin constitutive theory and the viscoplastic theory, thereby identifying the parameters that control the rheology of granular flow. The shear-rate number ( N ) proves to be among the most significant parameters identified from the GVF model. As N → 0 and N →∞, the GVF model can reduce asymptotically to the theoretical stress versus shear-rate relations in the macroviscous and grain-inertia regimes, respectively, where the grain concentration ( C ) also plays a major role in the rheology of granular flow. Using available data obtained from the rotating-cylinder experiments of neutrally buoyant solid spheres dispersing in an interstitial fluid, the shear stress for granular flow in transition between the two regimes proves dependent on N and C in addition to some material constants, such as the coefficient of restitution. The insufficiency of data on rotating-cylinder experiments cannot presently allow the GVF model to predict how a granular flow may behave in the entire range of N ; however, the analyzed data provide an insight on the interrelation among the relevant dimensionless parameters.

Journal of Engineering Mechanics

Impact of pore fluid chemistry on fine-grained sediment fabric and compressibility

Fines, defined here as grains or particles, less than 75 μm in diameter, exist nearly ubiquitously in natural sediment, even those classified as coarse. Macroscopic sediment properties, such as compressibility, which relates applied effective stress to the resulting sediment deformation, depend on the fabric of fines. Unlike coarse grains, fines have sizes and masses small enough to be more strongly influenced by electrical interparticle forces than by gravity. These electrical forces acting through pore fluids are influenced by pore fluid chemistry changes. Macroscopic property dependence on pore fluid chemistry must be accounted for in sediment studies involving subsurface flow and sediment stability analyses, as well as in engineered flow situations such as groundwater pollutant remediation, hydrocarbon migration, or other energy resource extraction applications. This study demonstrates how the liquid limit‐based electrical sensitivity index can be used to predict sediment compressibility changes due to pore fluid chemistry changes. Laboratory tests of electrical sensitivity, sedimentation, and compressibility illustrate mechanisms linking microscale and macroscale processes for selected pure, end‐member fines. A specific application considered here is methane extraction via depressurization of gas hydrate‐bearing sediment, which causes a dramatic pore water salinity drop concurrent with sediment being compressed by the imposed effective stress increase.

Journal of Geophysical Research

Generalized viscoplastic modeling of debris flow

Various concepts have been proposed or used in the development of Theological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold's concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non‐Newtonian fluid mechanics. Although the generality of Bagnold's model is still at issue, debris‐flow modelers in Japan have generally accepted Takahashi's formulas derived from Bagnold's model. Some efforts have recently been made by theoreticians in non‐Newtonian fluid mechanics to modify or improve Bagnold's concept or model. A viable rheological model should consist both of a rate‐independent part and a ratedependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris‐flow modeling. In fact, Bagnold's model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold's simplified assumption of constant grain concentration.

Journal of Hydraulic Engineering

Application of a coupled vegetation competition and groundwater simulation model to study effects of sea level rise and storm surges on coastal vegetation

Global climate change poses challenges to areas such as low-lying coastal zones, where sea level rise (SLR) and storm-surge overwash events can have long-term effects on vegetation and on soil and groundwater salinities, posing risks of habitat loss critical to native species. An early warning system is urgently needed to predict and prepare for the consequences of these climate-related impacts on both the short-term dynamics of salinity in the soil and groundwater and the long-term effects on vegetation. For this purpose, the U.S. Geological Survey’s spatially explicit model of vegetation community dynamics along coastal salinity gradients (MANHAM) is integrated into the USGS groundwater model (SUTRA) to create a coupled hydrology–salinity–vegetation model, MANTRA. In MANTRA, the uptake of water by plants is modeled as a fluid mass sink term. Groundwater salinity, water saturation and vegetation biomass determine the water available for plant transpiration. Formulations and assumptions used in the coupled model are presented. MANTRA is calibrated with salinity data and vegetation pattern for a coastal area of Florida Everglades vulnerable to storm surges. A possible regime shift at that site is investigated by simulating the vegetation responses to climate variability and disturbances, including SLR and storm surges based on empirical information.

Journal of Marine Science and Engineering

General solutions For viscoplastic debris flow

Theoretical velocity profile and theoretical pressure and concentration distributions for (steady) uniform debris flow in wide channels are derived from a generalized viscoplastic fluid (GVF) model without imposing R. A. Bagnold's assumption of constant grain concentration. Good agreement between the theoretical velocity profile and the experimental data of Japanese scientists strongly supports the validity of both the GVF model and the proposed method of solution from the model. It is shown that both E. C. Bingham and Bagnold versions (or submodels) of the GVF model can be used to simulate debris flow at the dynamic state. Although Bagnold's dilatant submodel appears to fit the Japanese data better than the Bingham submodel for flow of noncohesive grains, the choice between them is by no means clear-cut.

Journal of Hydraulic Engineering

Methane gas hydrate effect on sediment acoustic and strength properties

To improve our understanding of the interaction of methane gas hydrate with host sediment, we studied: (1) the effects of gas hydrate and ice on acoustic velocity in different sediment types, (2) effect of different hydrate formation mechanisms on measured acoustic properties (3) dependence of shear strength on pore space contents, and (4) pore pressure effects during undrained shear. A wide range in acoustic p-wave velocities (V p ) were measured in coarse-grained sediment for different pore space occupants. V p ranged from less than 1 km/s for gas-charged sediment to 1.77–1.94 km/s for water-saturated sediment, 2.91–4.00 km/s for sediment with varying degrees of hydrate saturation, and 3.88–4.33 km/s for frozen sediment. V p measured in fine-grained sediment containing gas hydrate was substantially lower (1.97 km/s). Acoustic models based on measured V p indicate that hydrate which formed in high gas flux environments can cement coarse-grained sediment, whereas hydrate formed from methane dissolved in the pore fluid may not. The presence of gas hydrate and other solid pore-filling material, such as ice, increased the sediment shear strength. The magnitude of that increase is related to the amount of hydrate in the pore space and cementation characteristics between the hydrate and sediment grains. We have found, that for consolidation stresses associated with the upper several hundred meters of sub-bottom depth, pore pressures decreased during shear in coarse-grained sediment containing gas hydrate, whereas pore pressure in fine-grained sediment typically increased during shear. The presence of free gas in pore spaces damped pore pressure response during shear and reduced the strengthening effect of gas hydrate in sands.

Journal of Petroleum Science and Engineering

Criteria for incipient motion of spherical sediment particles

Initiation of bed-load transport of uniform spherical sediment particles on a horizontal bed in an open-channel flow is studied. On the basis of micromechanical and fluid dynamical considerations, two separate criteria for the initiation of motion are derived: one for rolling and one for lifting. Fluid forces such as drag, shear lift, Magnus lift, and lift due to centrifugal force are included in the derivation. The formulation of the lift force is theoretical. No empirical coefficient is used other than the drag coefficient C d , which is well established in the literature. In the low particle Reynolds number regime, where u * d/ν ≤ 1, the dimensionless stress required to initiate lifting is found to be much higher than that of rolling, the former substantially above the Shields' curve while the latter is substantially below it. The theory is compared with Vanoni's 1964 data, with meaningful results. For higher particle Reynolds number, where u * d /ν> 10, the dimensionless stress needed to initiate lifting is closer to that of rolling. The Shields' curve for the most part lies between the two theoretical thresholds.

Journal of Hydraulic Engineering

Distributed shear of subglacial till due to Coulomb slip

In most models of the flow of glaciers on till beds, it has been assumed that till behaves as a viscoplastic fluid, despite contradictory evidence from laboratory studies. In accord with this assumption, displacement profiles measured in subglacial till have been fitted with viscoplastic models by estimating the stress distribution. Here we present a model that illustrates how observed displacement profiles can result from till deformation resisted solely by Coulomb friction. Motion in the till bed is assumed to be driven by brief departures from static equilibrium caused by fluctuations in effective normal stress. These fluctuations result from chains of particles that support intergranular forces that are higher than average and that form and fail at various depths in the bed during shearing. Newton’s second law is used to calculate displacements along slip planes and the depth to which deformation extends in the bed. Consequent displacement profiles are convex upward, similar to those measured by Boulton and colleagues at Breidamerkurjökull, Iceland. The model results, when considered together with the long-term and widespread empirical support for Coulomb models in soils engineering, indicate that efforts to fit viscoplastic flow models to till displacement profiles may be misguided.

Journal of Glaciology