Hydraulic fracture model comparison study: Complete results
Warpinski, N.R.; Abou-Sayed, I.S.; Moschovidis, Z.; Parker, C.
1993-02-01
Large quantities of natural gas exist in low permeability reservoirs throughout the US. Characteristics of these reservoirs, however, make production difficult and often economic and stimulation is required. Because of the diversity of application, hydraulic fracture design models must be able to account for widely varying rock properties, reservoir properties, in situ stresses, fracturing fluids, and proppant loads. As a result, fracture simulation has emerged as a highly complex endeavor that must be able to describe many different physical processes. The objective of this study was to develop a comparative study of hydraulic-fracture simulators in order to provide stimulation engineers with the necessary information to make rational decisions on the type of models most suited for their needs. This report compares the fracture modeling results of twelve different simulators, some of them run in different modes for eight separate design cases. Comparisons of length, width, height, net pressure, maximum width at the wellbore, average width at the wellbore, and average width in the fracture have been made, both for the final geometry and as a function of time. For the models in this study, differences in fracture length, height and width are often greater than a factor of two. In addition, several comparisons of the same model with different options show a large variability in model output depending upon the options chosen. Two comparisons were made of the same model run by different companies; in both cases the agreement was good. 41 refs., 54 figs., 83 tabs.
Air-water analogy and the study of hydraulic models
NASA Technical Reports Server (NTRS)
Supino, Giulio
1953-01-01
The author first sets forth some observations about the theory of models. Then he established certain general criteria for the construction of dynamically similar models in water and in air, through reference to the perfect fluid equations and to the ones pertaining to viscous flow. It is, in addition, pointed out that there are more cases in which the analogy is possible than is commonly supposed.
NASA Astrophysics Data System (ADS)
Kidson, R. L.; Richards, K. S.; Carling, P. A.
2006-02-01
Palaeoflood reconstructions based on stage evidence are typically conducted in data-poor field settings. Few opportunities exist to calibrate the hydraulic models used to estimate discharge from this evidence. Consequently, an important hydraulic model parameter, the roughness coefficient (e.g. Manning's n), is typically estimated by a range of approximate techniques, such as visual estimation and semi-empirical equations. These techniques contribute uncertainty to resulting discharge estimates, especially where the study reach exhibits sensitivity in the discharge-Manning's n relation. We study this uncertainty within a hydraulic model for a large flood of known discharge on the Mae Chaem River, northern Thailand. Comparison of the calibrated Manning's n with that obtained from semi-empirical equations indicates that these underestimate roughness. Substantial roughness elements in the extra-channel zone, inundated during large events, contribute significant additional sources of flow resistance that are captured neither by the semi-empirical equations, nor by existing models predicting stage-roughness variations. This bedrock channel exhibits a complex discharge-Manning's n relation, and reliable estimates of the former are dependent upon realistic assignment of the latter. Our study demonstrates that a large recent flood can provide a valuable opportunity to constrain this parameter, and this is illustrated when we model a palaeoflood event in the same reach, and subsequently examine the magnitude-return period consequences of discharge uncertainty within a flood frequency analysis, which contributes its own source of uncertainty.
Vermeyen, T.
1995-07-01
Bureau of Reclamation conducted this hydraulic model study to provide Pacific Gas and Electric Company with an evaluation of several selective withdrawal structures that are being considered to reduce intake flow temperatures through the Prattville Intake at Lake Almanor, California. Release temperature control using selective withdrawal structures is being considered in an effort to improve the cold-water fishery in the North Fork of the Feather River.
Hydraulic modelling for flood mapping and prevention: the case study of Cerfone River
NASA Astrophysics Data System (ADS)
Di Francesco, Silvia; Venturi, Sara; Manciola, Piergiorgio
2016-04-01
The research focuses on the hydraulic risk evaluation and danger estimation for different extreme flood events, in order to correctly implement mitigation measures in an anthropized basin. The Cerfone River (Tuscany, Italy), due to the several floods that have affected the neighbouring villages in recent years, is selected as case of study. A finite volume numerical model that solves the shallow water equations all over the computational domain, was used to simulate the unsteady evolution of the maximum extent of flooded areas for different scenarios. The one - dimensional approach (still widespread in engineering projects) can be inaccurate in complex flows, which are often two or three dimensional and sometimes does not manage to capture the flood spatial extents in terms of flow depth and velocity. The use of a two-dimensional numerical model seems to be the suitable instrument in terms of computational efficiency and adequacy of results. In fact it overcomes the limits of a one-dimensional modeling in terms of prediction of hydraulic variables with a less computational effort respect to a full 3d model. An accurate modeling of the river basin leads to the evaluation of the present hydraulic risk. Structural and non- structural measures are then studied, simulated and compared in order to define the optimal risk reduction plan for the area of study. At this aim, different flooding scenarios were simulated through the 2D mathematical model: i) existing state of the river and floodplain areas; ii) design of a levee to protect the most vulnerable populated areas against the flooding risk; iii) use of off - stream detention basins that strongly amplify the lamination capacity of floodplains. All these scenarios were simulated for different return periods: 50, 100, 200 and 500 years. The inputs of the hydraulic models are obtained in accordance with the legislative requirement of Tuscany Region; in particular discharge hydrographs are evaluate through the ALTo
NASA Astrophysics Data System (ADS)
Rajaram, H.; Birdsell, D.; Lackey, G.; Karra, S.; Viswanathan, H. S.; Dempsey, D.
2015-12-01
The dramatic increase in the extraction of unconventional oil and gas resources using horizontal wells and hydraulic fracturing (fracking) technologies has raised concerns about potential environmental impacts. Large volumes of hydraulic fracturing fluids are injected during fracking. Incidents of stray gas occurrence in shallow aquifers overlying shale gas reservoirs have been reported; whether these are in any way related to fracking continues to be debated. Computational models serve as useful tools for evaluating potential environmental impacts. We present modeling studies of hydraulic fracturing fluid and gas migration during the various stages of well operation, production, and subsequent plugging. The fluid migration models account for overpressure in the gas reservoir, density contrast between injected fluids and brine, imbibition into partially saturated shale, and well operations. Our results highlight the importance of representing the different stages of well operation consistently. Most importantly, well suction and imbibition both play a significant role in limiting upward migration of injected fluids, even in the presence of permeable connecting pathways. In an overall assessment, our fluid migration simulations suggest very low risk to groundwater aquifers when the vertical separation from a shale gas reservoir is of the order of 1000' or more. Multi-phase models of gas migration were developed to couple flow and transport in compromised wellbores and subsurface formations. These models are useful for evaluating both short-term and long-term scenarios of stray methane release. We present simulation results to evaluate mechanisms controlling stray gas migration, and explore relationships between bradenhead pressures and the likelihood of methane release and transport.
NASA Astrophysics Data System (ADS)
Gerke, Kirill; Khirevich, Siarhei; Sizonenko, Timofey; Karsanina, Marina; Umarova, Aminat; Korost, Dmitry; Matthai, Stephan; Mallants, Dirk
2016-04-01
To verify pore-scale modelling approach for determination of soil saturated hydraulic conductivity properties we scanned three cylindrical soil samples taken from A, Ah and B horizons using X-ray microtomography method. Resulting 3D soil images with resolutions of 15.25-20.96 μm were segmented into pores and solids and their maximum inscribed cube subvolumes were used as input data for three major pore-scale modelling methods to simulate saturated flow - lattice-Boltzmann method, finite-difference solution of the Stokes problem, and pore-network model. Provided that imaging resolution is high enough to capture the backbone of effective porosity and the main conducting pores all three methods resulted in simulated soil permeabilities close to experimental values for Ah and B samples. The resolution of A sample was not enough for an accurate modelling and we concluded that this soil requires multi-scale imaging to cover all relevant heterogeneities. We demonstrate that popular SWV method to choose segmentation threshold resulted in oversegmentation and order of magnitude higher permeability values. Careful manual thresholding combined with local segmentation algorithm provided much more accurate results. Detailed analysis of water retention curves showed that air-filled porosity at relevant pressure stages cannot be used for verification of the segmentation results. Representativity analysis by simulating flow in increasing soil volume up to 2.8 cm3 revealed no representative elementary volume (REV) within Ah sample and non-uniqueness of REV for B sample. The latter was explained by soil structure non-stationarity. We further speculate that structures soil horizons can exhibit no REV at all. We discuss numerous advantages of coupled imaging and pore-scale modelling approach and show how it can become a successor of the conventional soil coring method to parametrize large scale continuum models.
Hydraulic Redistribution: A Modeling Perspective
NASA Astrophysics Data System (ADS)
Daly, E.; Verma, P.; Loheide, S. P., III
2014-12-01
Roots play a key role in the soil water balance. They extract and transport water for transpiration, which usually represents the most important soil water loss in vegetated areas, and can redistribute soil water, thereby increasing transpiration rates and enhancing root nutrient uptake. We present here a two-dimensional model capable of describing two key aspects of root water uptake: root water compensation and hydraulic redistribution. Root water compensation is the ability of root systems to respond to the reduction of water uptake from areas of the soil with low soil water potential by increasing the water uptake from the roots in soil parts with higher water potential. Hydraulic redistribution is a passive transfer of water through the root system from areas of the soil with greater water potential to areas with lower water potential. Both mechanisms are driven by gradients of water potential in the soil and the roots. The inclusion of root water compensation and hydraulic redistribution in models can be achieved by describing root water uptake as a function of the difference in water potential between soil and root xylem. We use a model comprising the Richards equation for the water flow in variably saturated soils and the Darcy's equation for the water flow in the xylem. The two equations are coupled via a sink term, which is assumed to be proportional to the difference between soil and xylem water potentials. The model is applied in two case studies to describe vertical and horizontal hydraulic redistribution and the interaction between vegetation with different root depths. In the case of horizontal redistribution, the model is used to reproduce the fluxes of water across the root system of a tree subjected to uneven irrigation. This example can be extended to situations when only part of the root system has access to water, such as vegetation near creeks, trees at the edge of forests, and street trees in urban areas. The second case is inspired by recent
Cleary, M.P.
1994-02-01
This paper provides comments to a companion journal paper on predictive modeling of hydraulic fracturing patterns (N.R. Warpinski et. al., 1994). The former paper was designed to compare various modeling methods to demonstrate the most accurate methods under various geologic constraints. The comments of this paper are centered around potential deficiencies in the former authors paper which include: limited actual comparisons offered between models, the issues of matching predictive data with that from related field operations was lacking or undocumented, and the relevance/impact of accurate modeling on the overall hydraulic fracturing cost and production.
Anh Bui; Nam Dinh; Brian Williams
2013-09-01
In addition to validation data plan, development of advanced techniques for calibration and validation of complex multiscale, multiphysics nuclear reactor simulation codes are a main objective of the CASL VUQ plan. Advanced modeling of LWR systems normally involves a range of physico-chemical models describing multiple interacting phenomena, such as thermal hydraulics, reactor physics, coolant chemistry, etc., which occur over a wide range of spatial and temporal scales. To a large extent, the accuracy of (and uncertainty in) overall model predictions is determined by the correctness of various sub-models, which are not conservation-laws based, but empirically derived from measurement data. Such sub-models normally require extensive calibration before the models can be applied to analysis of real reactor problems. This work demonstrates a case study of calibration of a common model of subcooled flow boiling, which is an important multiscale, multiphysics phenomenon in LWR thermal hydraulics. The calibration process is based on a new strategy of model-data integration, in which, all sub-models are simultaneously analyzed and calibrated using multiple sets of data of different types. Specifically, both data on large-scale distributions of void fraction and fluid temperature and data on small-scale physics of wall evaporation were simultaneously used in this work’s calibration. In a departure from traditional (or common-sense) practice of tuning/calibrating complex models, a modern calibration technique based on statistical modeling and Bayesian inference was employed, which allowed simultaneous calibration of multiple sub-models (and related parameters) using different datasets. Quality of data (relevancy, scalability, and uncertainty) could be taken into consideration in the calibration process. This work presents a step forward in the development and realization of the “CIPS Validation Data Plan” at the Consortium for Advanced Simulation of LWRs to enable
NASA Astrophysics Data System (ADS)
Zehe, E.; Becker, R.; Schädel, W.
A dynamic system left without external disturbances, will always tend to a stable equilibrium state that is consistent with the internal physics. For natural soils such an equilibrium state is reached when the gradients of the total hydraulic potential tend to zero. This statement is still valid for heterogeneous soils, because the hydraulic po- tential is an intensive state variable and therefore continuous at discontinuities of the pore space. In contrary the soil water content is as an extensive property discontinu- ous at discontinuities of the pore space. Hence, a small scale soil moisture pattern that persists if the soil state tends to hydraulic equilibrium, reflects the lateral small scale variability of the pore space. The objectives of our study are to show a) whether and how we could use TDR observations to identify the small scale variability of the pore space. For that purpose we analyse artificial TDR measurements, taken from physi- cally based simulations of soil water dynamics in heterogeneous media. b) We want to introduce a new TDR technology which we call "Spatial TDR", that is suitable for that purposes. To produce the artificial TDR-datasets we generate random fields of soil porosity and saturated hydraulic conductivity with different statistical properties based on field data in a Luvisol and simulate artificial water dynamics in this model soil based on Richards-equation. Within this model soil we define several hypothetical "Spatial TDR" clusters, that differ in the lateral spacing and the number of the probes, in the temporal resolution of the hypothetical measurements and in the assumed mea- surement accuracy. If the model soil approaches hydraulic equilibrium, the remaining soil moisture pattern will be dominated by the statistical properties of the porosity. In contrary the variability of the hydraulic conductivity will dominate the soil moisture patterns during infiltration events. The hypothetical Spatial TDR measurements within the
Dynamics of Model Hydraulic Fracturing Liquid Studied by Two-Dimensional Infrared Spectroscopy
NASA Astrophysics Data System (ADS)
Daley, Kim; Kubarych, Kevin J.
2014-06-01
The technique of two-dimensional infrared (2DIR) spectroscopy is used to expose the chemical dynamics of various concentrations of polymers and their monomers in heterogeneous mixtures. An environmentally relevant heterogeneous mixture, which inspires this study, is hydraulic fracturing liquid (HFL). Hydraulic fracking is a technique used to extract natural gas from shale deposits. HFL consists of mostly water, proppant (sand), an emulsifier (guar), and other chemicals specific to the drilling site. Utilizing a metal carbonyl as a probe, we observe the spectral dynamics of the polymer, guar, and its monomer, mannose, and compare the results to see how hydration dynamics change with varying concentration. Another polymer, Ficoll, and its monomer, sucrose, are also compared to see how polymer size affects hydration dynamics. The two results are as follows: (1) Guar experiences collective hydration at high concentrations, where as mannose experiences independent hydration; (2) no collective hydration is observed for Ficoll in the same concentration range as guar, possibly due to polymer shape and size. HFL experiences extremely high pressure during natural gas removal, so future studies will focus on how increased pressure affects the hydration dynamics of polymers and monomers.
Helical coil thermal hydraulic model
NASA Astrophysics Data System (ADS)
Caramello, M.; Bertani, C.; De Salve, M.; Panella, B.
2014-11-01
A model has been developed in Matlab environment for the thermal hydraulic analysis of helical coil and shell steam generators. The model considers the internal flow inside one helix and its associated control volume of water on the external side, both characterized by their inlet thermodynamic conditions and the characteristic geometry data. The model evaluates the behaviour of the thermal-hydraulic parameters of the two fluids, such as temperature, pressure, heat transfer coefficients, flow quality, void fraction and heat flux. The evaluation of the heat transfer coefficients as well as the pressure drops has been performed by means of the most validated literature correlations. The model has been applied to one of the steam generators of the IRIS modular reactor and a comparison has been performed with the RELAP5/Mod.3.3 code applied to an inclined straight pipe that has the same length and the same elevation change between inlet and outlet of the real helix. The predictions of the developed model and RELAP5/Mod.3.3 code are in fairly good agreement before the dryout region, while the dryout front inside the helical pipes is predicted at a lower distance from inlet by the model.
NASA Astrophysics Data System (ADS)
McGuire, Daniel
A numerical tool for the simulation of the thermal dynamics of pipe networks with heat transfer has been developed with the novel capability of modeling supercritical fluids. The tool was developed to support the design and deployment of two thermal-hydraulic loops at Carleton University for the purpose of heat transfer studies in supercritical and near-critical fluids. First, the system was characterized based on its defining features; the characteristic length of the flow path is orders of magnitude larger than the other characteristic lengths that define the system's geometry; the behaviour of the working fluid in the supercritical thermodynamic state. An analysis of the transient thermal behaviour of the model's domains is then performed to determine the accuracy and range of validity of the modeling approach for simulating the transient thermal behaviour of a thermal-hydraulic loop. Preliminary designs of three test section geometries, for the purpose of heat transfer studies, are presented in support of the overall design of the Carleton supercritical thermal-hydraulic loops. A 7-rod-bundle, annular and tubular geometries are developed with support from the new numerical tool. Materials capable of meeting the experimental requirements while operating in supercritical water are determined. The necessary geometries to satisfy the experimental goals are then developed based on the material characteristics and predicted heat transfer behaviour from previous simulation results. An initial safety analysis is performed on the test section designs, where they are evaluated against the ASME Boiler, Pressure Vessel, and Pressure Piping Code standard, required for safe operation and certification.
An analytical model for hydraulic fracturing in shallow bedrock formations.
dos Santos, José Sérgio; Ballestero, Thomas Paul; Pitombeira, Ernesto da Silva
2011-01-01
A theoretical method is proposed to estimate post-fracturing fracture size and transmissivity, and as a test of the methodology, data collected from two wells were used for verification. This method can be employed before hydrofracturing in order to obtain estimates of the potential hydraulic benefits of hydraulic fracturing. Five different pumping test analysis methods were used to evaluate the well hydraulic data. The most effective methods were the Papadopulos-Cooper model (1967), which includes wellbore storage effects, and the Gringarten-Ramey model (1974), known as the single horizontal fracture model. The hydraulic parameters resulting from fitting these models to the field data revealed that as a result of hydraulic fracturing, the transmissivity increased more than 46 times in one well and increased 285 times in the other well. The model developed by dos Santos (2008), which considers horizontal radial fracture propagation from the hydraulically fractured well, was used to estimate potential fracture geometry after hydrofracturing. For the two studied wells, their fractures could have propagated to distances of almost 175 m or more and developed maximum apertures of about 2.20 mm and hydraulic apertures close to 0.30 mm. Fracturing at this site appears to have expanded and propagated existing fractures and not created new fractures. Hydraulic apertures calculated from pumping test analyses closely matched the results obtained from the hydraulic fracturing model. As a result of this model, post-fracturing geometry and resulting post-fracturing well yield can be estimated before the actual hydrofracturing. PMID:20572875
NASA Astrophysics Data System (ADS)
Abdelbaki, Chérifa; Benchaib, Mohamed Mouâd; Benziada, Salim; Mahmoudi, Hacène; Goosen, Mattheus
2016-04-01
For more effective management of water distribution network in an arid region, Mapinfo GIS (8.0) software was coupled with a hydraulic model (EPANET 2.0) and applied to a case study region, Chetouane, situated in the north-west of Algeria. The area is characterized not only by water scarcity but also by poor water management practices. The results showed that a combination of GIS and modeling permits network operators to better analyze malfunctions with a resulting more rapid response as well as facilitating in an improved understanding of the work performed on the network. The grouping of GIS and modeling as an operating tool allows managers to diagnosis a network, to study solutions of problems and to predict future situations. The later can assist them in making informed decisions to ensure an acceptable performance level for optimal network operation.
Studying Transonic Gases With a Hydraulic Analog
NASA Technical Reports Server (NTRS)
Wagner, W.; Lepore, F.
1986-01-01
Water table for hydraulic-flow research yields valuable information about gas flow at transonic speeds. Used to study fuel and oxidizer flow in high-pressure rocket engines. Method applied to gas flows in such equipment as furnaces, nozzles, and chemical lasers. Especially suitable when wall contours nonuniform, discontinuous, or unusually shaped. Wall shapes changed quickly for study and evaluated on spot. Method used instead of computer simulation when computer models unavailable, inaccurate, or costly to run.
ENHANCING HSPF MODEL CHANNEL HYDRAULIC REPRESENTATION
The Hydrological Simulation Program - FORTRAN (HSPF) is a comprehensive watershed model, which employs depth-area-volume-flow relationships known as hydraulic function table (FTABLE) to represent stream channel cross-sections and reservoirs. An accurate FTABLE determination for a...
Hydraulic modeling of heat dispersion in large lakes
NASA Technical Reports Server (NTRS)
Silberman, E.; Stefan, H.
1972-01-01
Case studies of hydraulic models are described for four major generating plants, including information and maps of thermal distribution. Information is of interest to agencies involved in thermal pollution control.
Hydraulic model of the Chesapeake Bay
NASA Technical Reports Server (NTRS)
Robinson, A. E., Jr.
1978-01-01
Preliminary planning for the formulation of the first year of hydraulic studies on the Chesapeake Bay model was recently completed. The primary purpose of this initial effort was to develop a study program that is both responsive to problems of immediate importance and at the same time ensure that from the very beginning of operation maximum economical use is made of the model. The formulation of this preliminary study plan involved an extensive analysis of the environmental, economic, and social aspects of a series of current problems in order to establish a priority listing of their importance. The study program that evolved is oriented towards the analysis of the effects of some of the works of man on the Chesapeake Bay estuarine environment.
Analytical Modeling of Shale Hydraulic Fracturing and Gas Production
NASA Astrophysics Data System (ADS)
Xu, W.
2012-12-01
Shale gas is abundant all over the world. Due to its extremely low permeability, extensive stimulation of a shale reservoir is always required for its economic production. Hydraulic fracturing has been the primary method of shale reservoir stimulation. Consequently the design and optimization of a hydraulic fracturing treatment plays a vital role insuring job success and economic production. Due to the many variables involved and the lack of a simple yet robust tool based on fundamental physics, horizontal well placement and fracturing job designs have to certain degree been a guessing game built on previous trial and error experience. This paper presents a method for hydraulic fracturing design and optimization in these environments. The growth of a complex hydraulic fracture network (HFN) during a fracturing job is equivalently represented by a wiremesh fracturing model (WFM) constructed on the basis of fracture mechanics and mass balance. The model also simulates proppant transport and placement during HFN growth. Results of WFM simulations can then be used as the input into a wiremesh production model (WPM) constructed based on WFM. WPM represents gas flow through the wiremesh HFN by an elliptic flow and the flow of gas in shale matrix by a novel analytical solution accounting for contributions from both free and adsorbed gases stored in the pore space. WPM simulation is validated by testing against numerical simulations using a commercially available reservoir production simulator. Due to the analytical nature of WFM and WPM, both hydraulic fracturing and gas production simulations run very fast on a regular personal computer and are suitable for hydraulic fracturing job design and optimization. A case study is presented to demonstrate how a non-optimized hydraulic fracturing job might have been optimized using WFM and WPM simulations.Fig. 1. Ellipsoidal representation of (a) stimulated reservoir and (b) hydraulic fracture network created by hydraulic
ERIC Educational Resources Information Center
Decker, Robert L.; Kirby, Klane
This curriculum guide contains a course in hydraulics to train entry-level workers for automotive mechanics and other fields that utilize hydraulics. The module contains 14 instructional units that cover the following topics: (1) introduction to hydraulics; (2) fundamentals of hydraulics; (3) reservoirs; (4) lines, fittings, and couplers; (5)…
Recent advances in modeling of well hydraulics
NASA Astrophysics Data System (ADS)
Yeh, Hund-Der; Chang, Ya-Chi
2013-01-01
Well hydraulics is a discipline to understand the process of flow to the well in an aquifer which is regarded as a source of groundwater. A variety of analytical and numerical models have been developed over the last few decades to provide a framework for understanding and quantifying the flow behavior in aquifer systems. In this review, we first briefly introduce the background of the theory of well hydraulics and the concepts, methodologies, and applications of analytical, semi-analytical, numerical and approximate methods in solving the well-hydraulic problems. We then address the subjects of current interests such as the incorporation of effects of finite well radius, wellbore storage, well partial penetration, and the presence of skin into various practical problems of groundwater flow. Furthermore, we also summarize recent developments of flow modeling such as the flow in aquifers with horizontal wells or collector wells, the capture zone delineation, and the non-Darcian flow in porous media and fractured formations. Finally, we present a comprehensive review on the numerical calculations for five well functions frequently appearing in well-hydraulic literature and suggest some topics in groundwater flow for future research.
Hydraulic modelling of sewage exfiltration.
Karpf, Christian; Traenckner, Jens; Krebs, Peter
2009-01-01
Exfiltration of waste water in sewer networks represents a potential danger for the soil and the aquifer. Various modelling approaches have been proposed to quantify sewerage exfiltration and its spatial and temporal variation. Common models are based on the law of Darcy, extended by a more or less detailed consideration of the expansion of leaks, the characteristics of the soil and the colmation layer. In the paper investigations are introduced, which are focused on the actual water content of the soil and its influence on exfiltration rates. Modelling results with HYDRUS 1D show, that under unsaturated conditions initial exfiltration rates are increased compared to saturated conditions. In experiments it was found, that the matrix potential increases the tightness of the colmation layer. Further a colmation model was deduced, which allows the calculation of the thickness and conductivity of the colmation layer. PMID:19403969
ERIC Educational Resources Information Center
Engelbrecht, Nancy; And Others
These instructional materials provide an orientation to hydraulics for use at the postsecondary level. The first of 12 sections presents an introduction to hydraulics, including discussion of principles of liquids, definitions, liquid flow, the two types of hydraulic fluids, pressure gauges, and strainers and filters. The second section identifies…
ENHANCING HYDROLOGICAL SIMULATION PROGRAM - FORTRAN MODEL CHANNEL HYDRAULIC REPRESENTATION
The Hydrological Simulation Program– FORTRAN (HSPF) is a comprehensive watershed model that employs depth-area - volume - flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross-sections and reservoirs. ...
Subgrid spatial variability of soil hydraulic functions for hydrological modelling
NASA Astrophysics Data System (ADS)
Kreye, Phillip; Meon, Günter
2016-07-01
State-of-the-art hydrological applications require a process-based, spatially distributed hydrological model. Runoff characteristics are demanded to be well reproduced by the model. Despite that, the model should be able to describe the processes at a subcatchment scale in a physically credible way. The objective of this study is to present a robust procedure to generate various sets of parameterisations of soil hydraulic functions for the description of soil heterogeneity on a subgrid scale. Relations between Rosetta-generated values of saturated hydraulic conductivity (Ks) and van Genuchten's parameters of soil hydraulic functions were statistically analysed. An universal function that is valid for the complete bandwidth of Ks values could not be found. After concentrating on natural texture classes, strong correlations were identified for all parameters. The obtained regression results were used to parameterise sets of hydraulic functions for each soil class. The methodology presented in this study is applicable on a wide range of spatial scales and does not need input data from field studies. The developments were implemented into a hydrological modelling system.
Development of a hydraulic model of the human systemic circulation
NASA Technical Reports Server (NTRS)
Sharp, M. K.; Dharmalingham, R. K.
1999-01-01
Physical and numeric models of the human circulation are constructed for a number of objectives, including studies and training in physiologic control, interpretation of clinical observations, and testing of prosthetic cardiovascular devices. For many of these purposes it is important to quantitatively validate the dynamic response of the models in terms of the input impedance (Z = oscillatory pressure/oscillatory flow). To address this need, the authors developed an improved physical model. Using a computer study, the authors first identified the configuration of lumped parameter elements in a model of the systemic circulation; the result was a good match with human aortic input impedance with a minimum number of elements. Design, construction, and testing of a hydraulic model analogous to the computer model followed. Numeric results showed that a three element model with two resistors and one compliance produced reasonable matching without undue complication. The subsequent analogous hydraulic model included adjustable resistors incorporating a sliding plate to vary the flow area through a porous material and an adjustable compliance consisting of a variable-volume air chamber. The response of the hydraulic model compared favorably with other circulation models.
ERIC Educational Resources Information Center
Decker, Robert L.
Designed for use in courses where students are expected to become proficient in the area of hydraulics, including diesel engine mechanic programs, this curriculum guide is comprised of fourteen units of instruction. Unit titles include (1) Introduction, (2) Fundamentals of Hydraulics, (3) Reservoirs, (4) Lines, Fittings, and Couplers, (5) Seals,…
NASA Astrophysics Data System (ADS)
Wollschläger, U.; Gerhards, H.; Schneider, S.; Roth, K.
2007-12-01
Estimating field scale hydraulic properties is still a challenge in hydrology. Most classical methods require undisturbed soil samples that have to be excavated during time consuming and labour intensive field work which is often followed by tedious measurements of hydraulic properties in the laboratory. Since these methods can only be applied with a limited number of samples, often only a few point measurements need to be used to characterize field scale hydraulic properties while layer geometry has to be derived from interpolation of these values and additional drilling. The combination of geophysical measurement techniques and hydraulic modeling offers an attractive alternative to bridge the gap between i) few accurate point measurements that are used to infer local hydraulic properties and ii) spatial mapping of the respective layers over large scales. We use a time series of water contents measured in a soil profile with time domain reflectometry to estimate hydraulic properties of the different soil layers with a 1D hydraulic inverse model. Here, hydraulic properties are estimated from \\it in situ \
NASA Astrophysics Data System (ADS)
Fuchsluger, Martin; Götzl, Gregor
2014-05-01
In general most aquifers have a much larger lateral extent than vertical. This fact leads to the application of the Dupuit-Forchheimer assumptions to many groundwater problems, whereas a two dimensional simulation is considered sufficient. By coupling transient fluid flow modeling with heat transport the 2D aquifer approximation is in many cases insufficient as it does not consider effects of the subjacent and overlying aquitards on heat propagation as well as the impact of surface climatic effects on shallow aquifers. A shallow Holocene aquifer in Vienna served as a case study to compare different modeling approaches in two and three dimensions in order to predict the performance and impact of a thermal aquifer utilization for heating (1.3 GWh) and cooling (1.4 GWh) of a communal building. With the assumption of a 6 doublets well field, the comparison was realized in three steps: At first a two dimensional model for unconfined flow was set up, assuming a varying hydraulic conductivity as well as a varying top and bottom elevation of the aquifer (gross - thickness). The model area was chosen along constant hydraulic head at steady state conditions. A second model was made by mapping solely the aquifer in three dimensions using the same subdomain and boundary conditions as defined in step one. The third model consists of a complete three dimensional geological build-up including the aquifer as well as the overlying and subjacent layers and additionally an annually variable climatic boundary condition at the surface. The latter was calibrated with measured water temperature at a nearby water gauge. For all three models the same annual operating mode of the 6 hydraulic doublets was assumed. Furthermore a limited maximal groundwater temperature at a range between 8 and 18 °C as well as a constrained well flow rate has been given. Finally a descriptive comparison of the three models concerning the extracted thermal power, drawdown, temperature distribution and Darcy
NASA Astrophysics Data System (ADS)
Kobor, J. S.; O'Connor, M. D.; Sherwood, M. N.
2013-12-01
Effective floodplain management and restoration requires a detailed understanding of floodplain processes not readily achieved using standard one-dimensional hydraulic modeling approaches. The application of more advanced numerical models is, however, often limited by the relatively high costs of acquiring the high-resolution topographic data needed for model development using traditional surveying methods. The increasing availability of LiDAR data has the potential to significantly reduce these costs and thus facilitate application of multi-dimensional hydraulic models where budget constraints would have otherwise prohibited their use. The accuracy and suitability of LiDAR data for supporting model development can vary widely depending on the resolution of channel and floodplain features, the data collection density, and the degree of vegetation canopy interference among other factors. More work is needed to develop guidelines for evaluating LiDAR accuracy and determining when and how best the data can be used to support numerical modeling activities. Here we present two recent case studies where LiDAR datasets were used to support floodplain and sediment transport modeling efforts. One LiDAR dataset was collected with a relatively low point density and used to study a small stream channel in coastal Marin County and a second dataset was collected with a higher point density and applied to a larger stream channel in western Sonoma County. Traditional topographic surveying was performed at both sites which provided a quantitative means of evaluating the LiDAR accuracy. We found that with the lower point density dataset, the accuracy of the LiDAR varied significantly between the active stream channel and floodplain whereas the accuracy across the channel/floodplain interface was more uniform with the higher density dataset. Accuracy also varied widely as a function of the density of the riparian vegetation canopy. We found that coupled 1- and 2-dimensional hydraulic
23 CFR 650.111 - Location hydraulic studies.
Code of Federal Regulations, 2010 CFR
2010-04-01
... § 650.111 (c) and (d) shall be summarized in environmental review documents prepared pursuant to 23 CFR... 23 Highways 1 2010-04-01 2010-04-01 false Location hydraulic studies. 650.111 Section 650.111... BRIDGES, STRUCTURES, AND HYDRAULICS Location and Hydraulic Design of Encroachments on Flood Plains §...
NASA Astrophysics Data System (ADS)
Kim, J.; Schumann, G.; Neal, J. C.; Lin, S.
2013-12-01
Earth is the only planet possessing an active hydrological system based on H2O circulation. However, after Mariner 9 discovered fluvial channels on Mars with similar features to Earth, it became clear that some solid planets and satellites once had water flows or pseudo hydrological systems of other liquids. After liquid water was identified as the agent of ancient martian fluvial activities, the valley and channels on the martian surface were investigated by a number of remote sensing and in-suit measurements. Among all available data sets, the stereo DTM and ortho from various successful orbital sensor, such as High Resolution Stereo Camera (HRSC), Context Camera (CTX), and High Resolution Imaging Science Experiment (HiRISE), are being most widely used to trace the origin and consequences of martian hydrological channels. However, geomorphological analysis, with stereo DTM and ortho images over fluvial areas, has some limitations, and so a quantitative modeling method utilizing various spatial resolution DTMs is required. Thus in this study we tested the application of hydraulics analysis with multi-resolution martian DTMs, constructed in line with Kim and Muller's (2009) approach. An advanced LISFLOOD-FP model (Bates et al., 2010), which simulates in-channel dynamic wave behavior by solving 2D shallow water equations without advection, was introduced to conduct a high accuracy simulation together with 150-1.2m DTMs over test sites including Athabasca and Bahram valles. For application to a martian surface, technically the acceleration of gravity in LISFLOOD-FP was reduced to the martian value of 3.71 m s-2 and the Manning's n value (friction), the only free parameter in the model, was adjusted for martian gravity by scaling it. The approach employing multi-resolution stereo DTMs and LISFLOOD-FP was superior compared with the other research cases using a single DTM source for hydraulics analysis. HRSC DTMs, covering 50-150m resolutions was used to trace rough
Effect of turbulence models on the submerged hydraulic jump simulation
NASA Astrophysics Data System (ADS)
Shekari, Y.; Javan, M.; Eghbalzadeh, A.
2015-05-01
This study presents a numerical investigation and prediction of the flow field in threedimensional submerged hydraulic jumps. The volume of fluid (VOF) method is used to simulate the free surface. The turbulent structure is simulated by using different turbulence models, such as the standard k-ɛ model, RNG k-ɛ model, realizable k-ɛ model, and Reynolds-stress model (RSM) closure schemes. The capabilities of the turbulence models are investigated with the standard wall functions and enhanced wall treatment methods. A comparison between the numerical and experimental results shows that the numerical model is adequate for predicting the flow pattern and free surface of submerged hydraulic jumps. The RNG k-ɛ turbulence model with the enhanced wall treatment method ensures the highest accuracy in the water surface simulation. Near the channel bed of a fully developed region, the RSM model with the enhanced wall treatment method shows better agreement with the experimental longitudinal velocity than the other turbulence models. The standard k-ɛ model predicts the longitudinal velocity more accurately than the RNG and realizable k-ɛ models.
Representing plant hydraulics in a global Earth system model.
NASA Astrophysics Data System (ADS)
Kennedy, D.; Gentine, P.
2015-12-01
Earth system models need improvement to reproduce observed seasonal and diurnal cycles of photosynthesis and respiration. Model water stress parameterizations lag behind the plant physiology literature. A plant hydraulics model is developed and deployed in a global Earth system model (NCAR CESM 1.2.2 with CLM 4.5). Assimilation and transpiration are attenuated according to literature cavitation curves. Water stress is evaluated based on plant functional type hydraulic parameters forced by soil moisture and atmospheric conditions. Resolving the plant water status allows for modelling divergent strategies for water stress. The case of isohydric versus anisohydric species is presented, showing that including plant hydraulic traits alter modelled photosynthesis and transpiration.
Design, test and model of a hybrid magnetostrictive hydraulic actuator
NASA Astrophysics Data System (ADS)
Chaudhuri, Anirban; Yoo, Jin-Hyeong; Wereley, Norman M.
2009-08-01
The basic operation of hybrid hydraulic actuators involves high frequency bi-directional operation of an active material that is converted to uni-directional motion of hydraulic fluid using valves. A hybrid actuator was developed using magnetostrictive material Terfenol-D as the driving element and hydraulic oil as the working fluid. Two different lengths of Terfenol-D rod, 51 and 102 mm, with the same diameter, 12.7 mm, were used. Tests with no load and with load were carried out to measure the performance for uni-directional motion of the output piston at different pumping frequencies. The maximum no-load flow rates were 24.8 cm3 s-1 and 22.7 cm3 s-1 with the 51 mm and 102 mm long rods respectively, and the peaks were noted around 325 Hz pumping frequency. The blocked force of the actuator was close to 89 N in both cases. A key observation was that, at these high pumping frequencies, the inertial effects of the fluid mass dominate over the viscous effects and the problem becomes unsteady in nature. In this study, we also develop a mathematical model of the hydraulic hybrid actuator in the time domain to show the basic operational principle under varying conditions and to capture phenomena affecting system performance. Governing equations for the pumping piston and output shaft were obtained from force equilibrium considerations, while compressibility of the working fluid was taken into account by incorporating the bulk modulus. Fluid inertia was represented by a lumped parameter approach to the transmission line model, giving rise to strongly coupled ordinary differential equations. The model was then used to calculate the no-load velocities of the actuator at different pumping frequencies and simulation results were compared with experimental data for model validation.
Modeling soil detachment capacity by rill flow using hydraulic parameters
NASA Astrophysics Data System (ADS)
Wang, Dongdong; Wang, Zhanli; Shen, Nan; Chen, Hao
2016-04-01
The relationship between soil detachment capacity (Dc) by rill flow and hydraulic parameters (e.g., flow velocity, shear stress, unit stream power, stream power, and unit energy) at low flow rates is investigated to establish an accurate experimental model. Experiments are conducted using a 4 × 0.1 m rill hydraulic flume with a constant artificial roughness on the flume bed. The flow rates range from 0.22 × 10-3 m2 s-1 to 0.67 × 10-3 m2 s-1, and the slope gradients vary from 15.8% to 38.4%. Regression analysis indicates that the Dc by rill flow can be predicted using the linear equations of flow velocity, stream power, unit stream power, and unit energy. Dc by rill flow that is fitted to shear stress can be predicted with a power function equation. Predictions based on flow velocity, unit energy, and stream power are powerful, but those based on shear stress, especially on unit stream power, are relatively poor. The prediction based on flow velocity provides the best estimates of Dc by rill flow because of the simplicity and availability of its measurements. Owing to error in measuring flow velocity at low flow rates, the predictive abilities of Dc by rill flow using all hydraulic parameters are relatively lower in this study compared with the results of previous research. The measuring accuracy of experiments for flow velocity should be improved in future research.
Modeling the effects of hydraulic stimulation on geothermal reservoirs
NASA Astrophysics Data System (ADS)
De Simone, Silvia; Vilarrasa, Victor; Carrera, Jesús; Alcolea, Andrés; Meier, Peter
2013-04-01
Geothermal energy represents a huge power source that can provide clean energy in potentially unlimited supply. When designing geothermal energy production from deep hot rocks, permeability is considered to control the economic efficiency of the heat extraction operations. In fact, a high permeability heat exchanger is required to achieve a cost-competitive power generation. The typical procedure entails intercepting naturally fractured rocks and enhancing their permeability by means of stimulation. Hydraulic stimulation is the most widely used method. It involves the massive injection of a large volume of water at high flow rates to increase the downhole pore pressure. This overpressure reduces the effective stresses, which tends to induce shearing along the fracture planes. In this way permeability is enhanced due to dilatancy, especially in the direction perpendicular to shear. These processes usually trigger microseismic events, which are sometimes of sufficient magnitude to be felt by the local population. This causes a negative impact on the local population and may compromise the continuation of the project. Hence, understanding the mechanisms triggering these induced micro-earthquakes is important to properly design and manage geothermal stimulation and operations so as to prevent them. We analyzed the thermo-hydro-mechanical response of a fractured deep rock mass subjected to hydraulic stimulation. Considering that seismicity is triggered when failure condition are reached, we studied the variation of the stress regime due to the hydraulic and thermal perturbations during fluid injection. Starting with a simplified model with constant permeability fault zones, more sophisticated schemes are considered to simulate the behavior of the discontinuity zones, including permeability variation associated to temperature, pressure and stress regime changes. Numerical simulations are performed using the finite element numerical code CODE_BRIGHT, which allows to solve
Mathematical modeling of bent-axis hydraulic piston motors
NASA Technical Reports Server (NTRS)
Bartos, R. D.
1992-01-01
Each of the DSN 70-m antennas uses 16 bent-axis hydraulic piston motors as part of the antenna drive system. On each of the two antenna axes, four motors are used to drive the antenna and four motors provide counter torque to remove the backlash in the antenna drive train. This article presents a mathematical model for bent-axis hydraulic piston motors. The model was developed to understand the influence of the hydraulic motors on the performance of the DSN 70-m antennas' servo control system.
We have conducted numerical simulation studies to assess the potential for injection-induced fault reactivation and notable seismic events associated with shale-gas hydraulic fracturing operations. The modeling is generally tuned toward conditions usually encountered in the Marce...
Evaluating models for predicting hydraulic characteristics of layered soils
NASA Astrophysics Data System (ADS)
Mavimbela, S. S. W.; van Rensburg, L. D.
2012-01-01
Soil water characteristic curve (SWCC) and unsaturated hydraulic conductivity (K-coefficient) are critical hydraulic properties governing soil water activity on layered soils. Sustainable soil water conservation would not be possible without accurate knowledge of these hydraulic properties. Infield rainwater harvesting (IRWH) is one conservation technique adopted to improve the soil water regime of a number of clay soils found in the semi arid areas of Free State province of South Africa. Given that SWCC is much easier to measure, most soil water studies rely on SWCC information to predict in-situ K-coefficients. This work validated this practice on the Tukulu, Sepane and Swartland layered soil profiles. The measured SWCC was first described using Brooks and Corey (1964), van Genuchten (1980) and Kasugi (1996) parametric models. The conductivity functions of these models were then required to fit in-situ based K-coefficients derived from instantaneous profile method (IPM). The same K-coefficient was also fitted by HYDRUS 1-D using optimised SWCC parameters. Although all parametric models fitted the measured SWCC fairly well their corresponding conductivity functions could not do the same when fitting the in-situ based K-coefficients. Overestimates of more than 2 orders of magnitude especially at low soil water content (SWC) were observed. This phenomenon was pronounced among the upper horizons that overlaid a clayey horizon. However, optimized α and n parameters using HYDRUS 1-D showed remarkable agreement between fitted and in-situ K-coefficient with root sum of squares error (RMSE) recording values not exceeding unity. During this exercise the Brooks and Corey was replaced by modified van Genuchten model (Vogel and Cislerova, 1988) since it failed to produce unique inverse solutions. The models performance appeared to be soil specific with van Genuchten-Mualem (1980) performing fairly well on the Orthic and neucutanic horizons while its modified form fitted very
Jones, Joseph L.; Haluska, Tana L.; Kresch, David L.
2001-01-01
A method of updating flood inundation maps at a fraction of the expense of using traditional methods was piloted in Washington State as part of the U.S. Geological Survey Urban Geologic and Hydrologic Hazards Initiative. Large savings in expense may be achieved by building upon previous Flood Insurance Studies and automating the process of flood delineation with a Geographic Information System (GIS); increases in accuracy and detail result from the use of very-high-accuracy elevation data and automated delineation; and the resulting digital data sets contain valuable ancillary information such as flood depth, as well as greatly facilitating map storage and utility. The method consists of creating stage-discharge relations from the archived output of the existing hydraulic model, using these relations to create updated flood stages for recalculated flood discharges, and using a GIS to automate the map generation process. Many of the effective flood maps were created in the late 1970?s and early 1980?s, and suffer from a number of well recognized deficiencies such as out-of-date or inaccurate estimates of discharges for selected recurrence intervals, changes in basin characteristics, and relatively low quality elevation data used for flood delineation. FEMA estimates that 45 percent of effective maps are over 10 years old (FEMA, 1997). Consequently, Congress has mandated the updating and periodic review of existing maps, which have cost the Nation almost 3 billion (1997) dollars. The need to update maps and the cost of doing so were the primary motivations for piloting a more cost-effective and efficient updating method. New technologies such as Geographic Information Systems and LIDAR (Light Detection and Ranging) elevation mapping are key to improving the efficiency of flood map updating, but they also improve the accuracy, detail, and usefulness of the resulting digital flood maps. GISs produce digital maps without manual estimation of inundated areas between
NASA Astrophysics Data System (ADS)
WöHling, Thomas; Vrugt, Jasper A.
2008-12-01
Most studies in vadose zone hydrology use a single conceptual model for predictive inference and analysis. Focusing on the outcome of a single model is prone to statistical bias and underestimation of uncertainty. In this study, we combine multiobjective optimization and Bayesian model averaging (BMA) to generate forecast ensembles of soil hydraulic models. To illustrate our method, we use observed tensiometric pressure head data at three different depths in a layered vadose zone of volcanic origin in New Zealand. A set of seven different soil hydraulic models is calibrated using a multiobjective formulation with three different objective functions that each measure the mismatch between observed and predicted soil water pressure head at one specific depth. The Pareto solution space corresponding to these three objectives is estimated with AMALGAM and used to generate four different model ensembles. These ensembles are postprocessed with BMA and used for predictive analysis and uncertainty estimation. Our most important conclusions for the vadose zone under consideration are (1) the mean BMA forecast exhibits similar predictive capabilities as the best individual performing soil hydraulic model, (2) the size of the BMA uncertainty ranges increase with increasing depth and dryness in the soil profile, (3) the best performing ensemble corresponds to the compromise (or balanced) solution of the three-objective Pareto surface, and (4) the combined multiobjective optimization and BMA framework proposed in this paper is very useful to generate forecast ensembles of soil hydraulic models.
Implementation of diverse tree hydraulics in a land surface model
NASA Astrophysics Data System (ADS)
Wolf, A.; Shevliakova, E.; Malyshev, S.; Weng, E.; Pacala, S. W.
2013-12-01
Increasing attention has been devoted to the occurence of drought kill in forests worldwide. These mortality events are significant disruptions to the terrestrial carbon cycle, but the mechanisms required to represent drought kill are not represented in terrestrial carbon cycle models. In part, this is due to the challenge of representing the diversity of hydraulic strategies, which include stomatal sensitivity to water deficit and woody tissue vulnerability to cavitation at low water potential. In part, this is due to the challenge of representing this boundary value problem numerically, because the hydraulic components determine water potential at the leaf, but the stomatal conductance on the leaf also determines the hydraulic gradients within the plant. This poster will describe the development of a land surface model parameterization of diverse tree hydraulic strategies.
Bayesian Model Averaging of Artificial Intelligence Models for Hydraulic Conductivity Estimation
NASA Astrophysics Data System (ADS)
Nadiri, A.; Chitsazan, N.; Tsai, F. T.; Asghari Moghaddam, A.
2012-12-01
This research presents a Bayesian artificial intelligence model averaging (BAIMA) method that incorporates multiple artificial intelligence (AI) models to estimate hydraulic conductivity and evaluate estimation uncertainties. Uncertainty in the AI model outputs stems from error in model input as well as non-uniqueness in selecting different AI methods. Using one single AI model tends to bias the estimation and underestimate uncertainty. BAIMA employs Bayesian model averaging (BMA) technique to address the issue of using one single AI model for estimation. BAIMA estimates hydraulic conductivity by averaging the outputs of AI models according to their model weights. In this study, the model weights were determined using the Bayesian information criterion (BIC) that follows the parsimony principle. BAIMA calculates the within-model variances to account for uncertainty propagation from input data to AI model output. Between-model variances are evaluated to account for uncertainty due to model non-uniqueness. We employed Takagi-Sugeno fuzzy logic (TS-FL), artificial neural network (ANN) and neurofuzzy (NF) to estimate hydraulic conductivity for the Tasuj plain aquifer, Iran. BAIMA combined three AI models and produced better fitting than individual models. While NF was expected to be the best AI model owing to its utilization of both TS-FL and ANN models, the NF model is nearly discarded by the parsimony principle. The TS-FL model and the ANN model showed equal importance although their hydraulic conductivity estimates were quite different. This resulted in significant between-model variances that are normally ignored by using one AI model.
Fluid-solid interaction model for hydraulic reciprocating O-ring seals
NASA Astrophysics Data System (ADS)
Liao, Chuanjun; Huang, Weifeng; Wang, Yuming; Suo, Shuangfu; Liu, Ying
2013-01-01
Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. Only elastic deformations of hydraulic reciprocating seals were discussed, and hydrodynamic effects were neglected in many studies. The physical process of the fluid-solid interaction effect did not be clearly presented in the existing fluid-solid interaction models for hydraulic reciprocating O-ring seals, and few of these models had been simultaneously validated through experiments. By exploring the physical process of the fluid-solid interaction effect of the hydraulic reciprocating O-ring seal, a numerical fluid-solid interaction model consisting of fluid lubrication, contact mechanics, asperity contact and elastic deformation analyses is constructed with an iterative procedure. With the SRV friction and wear tester, the experiments are performed to investigate the elastohydrodynamic lubrication characteristics of the O-ring seal. The regularity of the friction coefficient varying with the speed of reciprocating motion is obtained in the mixed lubrication condition. The experimental result is used to validate the fluid-solid interaction model. Based on the model, The elastohydrodynamic lubrication characteristics of the hydraulic reciprocating O-ring seal are presented respectively in the dry friction, mixed lubrication and full film lubrication conditions, including of the contact pressure, film thickness, friction coefficient, liquid film pressure and viscous shear stress in the sealing zone. The proposed numerical fluid-solid interaction model can be effectively used to analyze the operation characteristics of the hydraulic reciprocating O-ring seal, and can also be widely used to study other hydraulic reciprocating seals.
NASA Astrophysics Data System (ADS)
Tsai, Frank T.-C.; Li, Xiaobao
2008-09-01
This study proposes a Bayesian model averaging (BMA) method to address parameter estimation uncertainty arising from nonuniqueness in parameterization methods. BMA is able to incorporate multiple parameterization methods for prediction through the law of total probability and to obtain an ensemble average of hydraulic conductivity estimates. Two major issues in applying BMA to hydraulic conductivity estimation are discussed. The first problem is using Occam's window in usual BMA applications to measure approximated posterior model probabilities. Occam's window only accepts models in a very narrow range, tending to single out the best method and discard other good methods. We propose a variance window to replace Occam's window to cope with this problem. The second problem is the Kashyap information criterion (KIC) in the approximated posterior model probabilities, which tends to prefer highly uncertain parameterization methods by considering the Fisher information matrix. With sufficient amounts of observation data, the Bayesian information criterion (BIC) is a good approximation and is able to avoid controversial results from using KIC. This study adopts multiple generalized parameterization (GP) methods such as the BMA models to estimate spatially correlated hydraulic conductivity. Numerical examples illustrate the issues of using KIC and Occam's window and show the advantages of using BIC and the variance window in BMA application. Finally, we apply BMA to the hydraulic conductivity estimation of the "1500-foot" sand in East Baton Rouge Parish, Louisiana.
Refilling of a Hydraulically Isolated Embolized Xylem Vessel: Model Calculations
VESALA, TIMO; HÖLTTÄ, TEEMU; PERÄMÄKI, MARTTI; NIKINMAA, EERO
2003-01-01
When they are hydraulically isolated, embolized xylem vessels can be refilled, while adjacent vessels remain under tension. This implies that the pressure of water in the refilling vessel must be equal to the bubble gas pressure, which sets physical constraints for recovery. A model of water exudation into the cylindrical vessel and of bubble dissolution based on the assumption of hydraulic isolation is developed. Refilling is made possible by the turgor of the living cells adjacent to the refilling vessel, and by a reflection coefficient below 1 for the exchange of solutes across the interface between the vessel and the adjacent cells. No active transport of solutes is assumed. Living cells are also capable of importing water from the water‐conducting vessels. The most limiting factors were found to be the osmotic potential of living cells and the ratio of the volume of the adjacent living cells to that of the embolized vessel. With values for these of 1·5 MPa and 1, respectively, refilling times were in the order of hours for a broad range of possible values of water conductivity coefficients and effective diffusion distances for dissolved air, when the xylem water tension was below 0·6 MPa and constant. Inclusion of the daily pattern for xylem tension improved the simulations. The simulated gas pressure within the refilling vessel was in accordance with recent experimental results. The study shows that the refilling process is physically possible under hydraulic isolation, while water in surrounding vessels is under negative pressure. However, the osmotic potentials in the refilling vessel tend to be large (in the order of 1 MPa). Only if the xylem water tension is, at most, twice atmospheric pressure, the reflection coefficient remains close to 1 (0·95) and the ratio of the volume of the adjacent living cells to that of the embolized vessel is about 2, does the osmotic potential stay below 0·4 MPa. PMID:12588721
Coupled Hydrological and Hydraulic Modeling for Flood Mapping
NASA Astrophysics Data System (ADS)
Drobot, Radu; Draghia, Aurelian
2014-05-01
The delineation of the flooded areas involves both hydrological and hydraulic modeling. Usually, the hydrological and hydraulic processes are separately treated. In the proposed methodology, the coupled modeling of the hydrological and hydraulic processes is used. The calibration and validation of the hydrological parameters is undertaken based on historical floods using the corresponding precipitations for the same period. The calibration process was more complicated in the presence of reservoirs, when not only the discharges downstream but also the water level in the reservoirs had to be accurately reproduced. The time step for precipitation is 1 hour, corresponding to the concentration time of the smallest catchments. The maximum annual precipitation for different time steps (1; 3; 6; 24 hours) were statistically processed and based on these results the cumulative rainfall curves and the synthetic hyetographs were derived. The rainfall duration is depending on the concentration time. Mike 11 with UHM module based on SCS model was used for coupled hydrological and hydraulic modeling. The coupled hydrological and hydraulic simulation for the scaled precipitation leads both at the computation of the components which contribute to the generation of the P% flood at the Hydrometric stations as well as to the determination of the discharge hydrograph along the main river. Based on these results the flood hazard maps were obtained using a DTM based on Lidar data. The methodology was applied for a river basin in Romania of 12500 km2.
Hydraulic fracturing: insights from field, lab, and numerical studies
NASA Astrophysics Data System (ADS)
Walsh, S. D.; Johnson, S.; Fu, P.; Settgast, R. R.
2011-12-01
Hydraulic fracturing has become an increasingly important technique in stimulating reservoirs for gas, oil, and geothermal energy production. In use commercially since the 1950's, the technique has been widely lauded, when combined with other techniques, for enabling the development of shale gas resources in the United States, providing a valuable and extensive source of domestic energy. However, the technique has also drawn a degree of notoriety from high-profile incidents involving contamination of drinking water associated with gas extraction operations in the Marcellus shale region. This work highlights some of the insights on the behavior of subsurface hydraulic fracturing operations that have been derived from field and laboratory observations as well as from numerical simulations. The sensitivity of fracture extent and orientation to parameters such as matrix material heterogeneity, presence and distribution of discontinuities, and stress orientation is of particular interest, and we discuss this in the context of knowledge derived from both observation and simulation. The limitations of these studies will also be addressed in terms of resolution, uncertainty, and assumptions as well as the balance of fidelity to cost, both in computation time (for numerical studies) and equipment / operation cost (for observational studies). We also identify a number of current knowledge gaps and propose alternatives for addressing those gaps. We especially focus on the role of numerical studies for elucidating key concepts and system sensitivities. The problem is inherently multi-scale in both space and time as well as highly coupled hydromechanically, and, in several applications, thermally as well. We will summarize the developments to date in analyzing these systems and present an approach for advancing the capabilities of our models in the short- to long-term and how these advances can help provide solutions to reduce risk and improve efficiency of hydraulic fracturing
Effects of model layer simplification using composite hydraulic properties
Kuniansky, Eve L.; Sepulveda, Nicasio
2011-01-01
Groundwater provides much of the fresh drinking water to more than 1.5 billion people in the world (Clarke et al., 1996) and in the United States more that 50 percent of citizens rely on groundwater for drinking water (Solley et al., 1998). As aquifer systems are developed for water supply, the hydrologic system is changed. Water pumped from the aquifer system initially can come from some combination of inducing more recharge, water permanently removed from storage, and decreased groundwater discharge. Once a new equilibrium is achieved, all of the pumpage must come from induced recharge and decreased discharge (Alley et al., 1999). Further development of groundwater resources may result in reductions of surface water runoff and base flows. Competing demands for groundwater resources require good management. Adequate data to characterize the aquifers and confining units of the system, like hydrologic boundaries, groundwater levels, streamflow, and groundwater pumping and climatic data for recharge estimation are to be collected in order to quantify the effects of groundwater withdrawals on wetlands, streams, and lakes. Once collected, three-dimensional (3D) groundwater flow models can be developed and calibrated and used as a tool for groundwater management. The main hydraulic parameters that comprise a regional or subregional model of an aquifer system are the hydraulic conductivity and storage properties of the aquifers and confining units (hydrogeologic units) that confine the system. Many 3D groundwater flow models used to help assess groundwater/surface-water interactions require calculating ?effective? or composite hydraulic properties of multilayered lithologic units within a hydrogeologic unit. The calculation of composite hydraulic properties stems from the need to characterize groundwater flow using coarse model layering in order to reduce simulation times while still representing the flow through the system accurately. The accuracy of flow models with
EPA Study of Hydraulic Fracturing and Drinking Water Resources
In its FY2010 Appropriations Committee Conference Report, Congress directed EPA to study the relationship between hydraulic fracturing and drinking water, using: • Best available science • Independent sources of information • Transparent, peer-reviewed process • Consultatio...
Empirical flow parameters : a tool for hydraulic model validity
Asquith, William H.; Burley, Thomas E.; Cleveland, Theodore G.
2013-01-01
The objectives of this project were (1) To determine and present from existing data in Texas, relations between observed stream flow, topographic slope, mean section velocity, and other hydraulic factors, to produce charts such as Figure 1 and to produce empirical distributions of the various flow parameters to provide a methodology to "check if model results are way off!"; (2) To produce a statistical regional tool to estimate mean velocity or other selected parameters for storm flows or other conditional discharges at ungauged locations (most bridge crossings) in Texas to provide a secondary way to compare such values to a conventional hydraulic modeling approach. (3.) To present ancillary values such as Froude number, stream power, Rosgen channel classification, sinuosity, and other selected characteristics (readily determinable from existing data) to provide additional information to engineers concerned with the hydraulic-soil-foundation component of transportation infrastructure.
Effects of model layer simplification using composite hydraulic properties
Sepulveda, Nicasio; Kuniansky, Eve L.
2010-01-01
The effects of simplifying hydraulic property layering within an unconfined aquifer and the underlying confining unit were assessed. The hydraulic properties of lithologic units within the unconfined aquifer and confining unit were computed by analyzing the aquifer-test data using radial, axisymmetric two-dimensional (2D) flow. Time-varying recharge to the unconfined aquifer and pumping from the confined Upper Floridan aquifer (USA) were simulated using 3D flow. Conceptual flow models were developed by gradually reducing the number of lithologic units in the unconfined aquifer and confining unit by calculating composite hydraulic properties for the simplified lithologic units. Composite hydraulic properties were calculated using either thickness-weighted averages or inverse modeling using regression-based parameter estimation. No significant residuals were simulated when all lithologic units comprising the unconfined aquifer were simulated as one layer. The largest residuals occurred when the unconfined aquifer and confining unit were aggregated into a single layer (quasi-3D), with residuals over 100% for the leakage rates to the confined aquifer and the heads in the confining unit. Residuals increased with contrasts in vertical hydraulic conductivity between the unconfined aquifer and confining unit. Residuals increased when the constant-head boundary at the bottom of the Upper Floridan aquifer was replaced with a no-flow boundary.
Influence of Rock Fabric on Hydraulic Fracture Propagation: Laboratory Study
NASA Astrophysics Data System (ADS)
Stanchits, S. A.; Desroches, J.; Burghardt, J.; Surdi, A.; Whitney, N.
2014-12-01
Massive hydraulic fracturing is required for commercial gas production from unconventional reservoirs. These reservoirs are often highly fractured and heterogeneous, which may cause significant fracture complexity and also arrest propagation of hydraulic fractures, leading to production decrease. One of the goals of our study was to investigate the influence of rock fabric features on near-wellbore fracture geometry and complexity. We performed a series of laboratory tests on Niobrara outcrop shale blocks with dimensions of 30 x 30 x 36 inches in a true-triaxial loading frame. Acoustic Emission (AE) technique was applied to monitor hydraulic fracture initiation and dynamics of fracture propagation. After the tests, the shape of the created hydraulic fracture was mapped by goniometry technique. To estimate fracture aperture, particles of different sizes were injected with fracturing fluid. In all tests, AE analysis indicated hydraulic fracture initiation prior to breakdown or the maximum of wellbore pressure. In most tests, AE analysis revealed asymmetrical hydraulic fracture shapes. Post-test analysis demonstrated good correspondence of AE results with the actual 3D shape of the fracture surface map. AE analysis confirmed that in some of these tests, the hydraulic fracture approached one face of the block before the maximum wellbore pressure had been reached. We have found that in such cases the propagation of hydraulic fracture in the opposite direction was arrested by the presence of mineralized interfaces. Mapping the distribution of injected particles confirmed the creation of a narrow-width aperture in the vicinity of pre-existing interfaces, restricting fracture conductivity. Based on the results of our study, we concluded that the presence of planes of weakness, such as mineralized natural fractures, can result in the arrest of hydraulic fracture propagation, or in poor fracture geometries with limited aperture, that in turn could lead to high net pressure
EPA releases progress report on hydraulic fracturing study
NASA Astrophysics Data System (ADS)
Showstack, Randy
2013-01-01
The U.S. Environmental Protection Agency (EPA) provided a 21 December progress report on its ongoing national study about the potential impacts of hydraulic fracturing on drinking water resources. The agency said that a draft of the congressionally requested study will be released in 2014 for public and peer review and that its progress report does not draw conclusions about the potential impacts of hydraulic fracturing, often referred to as fracking.
Mechanical stability of propped hydraulic fractures: A numerical study
Asgian, M.I.; Cundall, P.A.; Brady, B.H.
1995-03-01
Proppant is sometimes produced along with hydrocarbons in hydraulically fractured petroleum wells. Sometimes 10% to 20% of the proppant is backproduced, which can lead to damaged equipment and downtime. Furthermore, proppant flowback can lead to a substantial loss of fracture conductivity. A numerical study was conducted to help understand what conditions are likely to lead to proppant flowback. In the simulations, the mechanical interaction of a larger number (several thousand) individual proppant grains was modeled with a distinct-element-type code. The numerical simulations show that hydraulic fractures propped with cohesionless, unbonded proppant fail under closure stress at a critical ratio of mean grain diameter to fracture width. This is consistent with published laboratory studies. The simulations identify the mechanism (arch failure) that triggers the mechanical instability and also show that the primary way that drawdowns (less than {approx} 75 psi/ft) affect proppant flowback is to transport loose proppant grains in front of the stable arch to the wellbore. Drawdowns > 75 psi/ft are sufficient to destabilize the arch and to cause progressive failure of the propped fractures.
NASA Astrophysics Data System (ADS)
Hsiao, Kai-Wen; Hsu, Yu-Chao; Jan, Chyan-Deng; Su, Yu-Wen
2016-04-01
The inclined rectangular chute construction is a common structure used in hydraulic engineering for typical reasons such as the increase of bottom slope, the transition from side channel intakes to tunnel spillways, the drainage construction, and the reduction of chute width due to bridges, flood diversion structures or irrigation systems. The converging vertical sidewalls of a chute contraction deflect the supercritical flow to form hydraulic shock waves. Hydraulic shock waves have narrow and locally extreme wavy surfaces, which commonly results in the requirement of higher height of sidewalls. Therefore, predicting the possible height and position of maximum hydraulic shock wave are necessary to design the required height of sidewalls to prevent flow overtopping. In this study, we used a three-dimensional computation fluid dynamics model (i.e., FLOW-3D) to simulate the characteristics of hydraulic shock waves in an inclined chute contraction. For this purpose, the parameters of simulated hydraulic shock wave, such as the shock angle, maximum shock wave height and maximum shock wave position in various conditions are compared with those calculated by the empirical relations obtained from literatures. We showed that the simulated results are extremely close to the experimental results. The numerical results validated the applicability of these empirical relations and extend their applicability to higher approach Froude numbers from 3.51 to 7.27. Furthermore, we also applied the Yuan-Shan-Tsu flood diversion channel under 200-year peak flow condition to FLOW-3D model to simulate the hydraulic shock waves and validate the effect of the installation of a diversion pier in the channel on promoting the stability of flow fluid. The results revealed that a diversion pier installed in the Yuan-Shan-Tsu flood diversion channel is helpful for improving the stability of flow field. In summary, this study demonstrates that FLOW-3D model can be used to simulate the
Coupled thermal-hydraulic-chemical modelling of enhanced geothermal systems
NASA Astrophysics Data System (ADS)
Bächler, D.; Kohl, T.
2005-05-01
The study investigates thermal-, hydraulic- and chemically coupled processes of enhanced geothermal systems (EGS). On the basis of the two existing numerical codes, the finite element program FRACTURE and the geochemical module of CHEMTOUGH, FRACHEM was developed, to simulate coupled thermal-hydraulic-chemical (THC) processes, accounting for the Soultz specific conditions such as the high salinity of the reservoir fluid and the high temperatures. The finite element part calculates the thermal and hydraulic field and the geochemical module the chemical processes. According to the characteristics of the Soultz EGS reservoir, the geochemical module was modified. (i) The Debye-Huckel approach was replaced by the Pitzer formalism. (ii) New kinetic laws for calcite, dolomite, quartz and pyrite were implemented. (iii) The porosity-permeability relation was replaced by a new relation for fractured rock. (iv) The possibility of re-injecting the produced fluid was implemented. The sequential non-iterative approach (SNIA) was used to couple transport and reactions. Sensitivity analyses proved the proper functionality of FRACHEM, but highlighted the sensitivity of the SNIA approach to time steps. To quantify the FRACHEM results, a comparative simulation with the code SHEMAT was conducted, which validated FRACHEM. Coupled THC processes in a fractured zone in the Soultz reservoir at 3500 m (T0= 165 °C), which occur as a result of the injection of fluid (Tinj= 65 °C) at one end of the zone and the production at the other end, were modelled for 2 yr. Calcite is the most reactive mineral and therefore the porosity and permeability evolution results from the calcite reactions: near the injection point, porosity and permeability increase and near the production well they decrease. After 2 yr, the system seems to be very close to steady-state. Therefore, mineral dissolution and precipitation during the circulation of the fluid in the reservoir do not represent a limiting factor on
Simplified hydraulic model of French vertical-flow constructed wetlands.
Arias, Luis; Bertrand-Krajewski, Jean-Luc; Molle, Pascal
2014-01-01
Designing vertical-flow constructed wetlands (VFCWs) to treat both rain events and dry weather flow is a complex task due to the stochastic nature of rain events. Dynamic models can help to improve design, but they usually prove difficult to handle for designers. This study focuses on the development of a simplified hydraulic model of French VFCWs using an empirical infiltration coefficient--infiltration capacity parameter (ICP). The model was fitted using 60-second-step data collected on two experimental French VFCW systems and compared with Hydrus 1D software. The model revealed a season-by-season evolution of the ICP that could be explained by the mechanical role of reeds. This simplified model makes it possible to define time-course shifts in ponding time and outlet flows. As ponding time hinders oxygen renewal, thus impacting nitrification and organic matter degradation, ponding time limits can be used to fix a reliable design when treating both dry and rain events. PMID:25225940
LVP modeling and dynamic characteristics prediction of a hydraulic power unit in deep-sea
NASA Astrophysics Data System (ADS)
Cao, Xue-peng; Ye, Min; Deng, Bin; Zhang, Cui-hong; Yu, Zu-ying
2013-03-01
A hydraulic power unit (HPU) is the driving "heart" of deep-sea working equipment. It is critical to predict its dynamic performances in deep-water before being immerged in the seawater, while the experimental tests by simulating deep-sea environment have many disadvantages, such as expensive cost, long test cycles, and difficult to achieve low-temperature simulation, which is only used as a supplementary means for confirmatory experiment. This paper proposes a novel theoretical approach based on the linear varying parameters (LVP) modeling to foresee the dynamic performances of the driving unit. Firstly, based on the varying environment features, dynamic expressions of the compressibility and viscosity of hydraulic oil are derived to reveal the fluid performances changing. Secondly, models of hydraulic system and electrical system are accomplished respectively through studying the control process and energy transfer, and then LVP models of the pressure and flow rate control is obtained through the electro-hydraulic models integration. Thirdly, dynamic characteristics of HPU are obtained by the model simulating within bounded closed sets of varying parameters. Finally, the developed HPU is tested in a deep-sea imitating hull, and the experimental results are well consistent with the theoretical analysis outcomes, which clearly declare that the LVP modeling is a rational way to foresee dynamic performances of HPU. The research approach and model analysis results can be applied to the predictions of working properties and product designs for other deep-sea hydraulic pump.
European hydraulic geometries for continental SCALE environmental modelling
NASA Astrophysics Data System (ADS)
Pistocchi, A.; Pennington, D.
2006-10-01
SummaryThis paper presents a geo-exploratory analysis carried out using publicly accessible data to estimate river channel width and mean lake depth. These two parameters, combined with other more readily available information, are key in many environmental assessment applications such as sedimentological and ecological studies, and the prediction of surface water residence times. For river width, runoff estimates have been used to produce a European map of mean annual discharge based on the flow accumulation computed for the GTOPO30 digital elevation model (DEM). A regression equation is fitted to predict river width as a function of river discharge. For lake depth, a landscape roughness index has been calculated on the GTOPO30 DEM. Based on this index, a regression equation predicts lake depth, allowing to populate the existing geographic databases with information on estimated depth. Using these models for river width and lake depth it is possible to compute "representative" values of velocity and depth in rivers, as well as residence time in lakes at a pan-European level. These parameters define the hydraulic geometry of European inland waters for typical conditions with errors on the predictions generally within a factor of 2. Errors are generally unbiased, and geographic patterns (such as macroscopic differentiation of catchments and physiographic environments) are reproduced in a qualitatively correct way, so that overall representativeness of the values for continental scale modelling is further corroborated.
NASA Astrophysics Data System (ADS)
Williams, M. D.; Thorne, P. D.; Bergeron, M. P.; Vermeul, V. R.; Ward, D. L.
2006-12-01
A three dimensional groundwater flow and transport model was calibrated against observations of both hydraulic head and tritium plume concentrations measured in wells. Hydraulic parameters were estimated with a transient inverse process using UCODE, a universal inverse modeling code developed jointly by the U.S. Geological Survey and the International Groundwater Modeling Center at the Colorado School of Mines. Previous groundwater models at the site had been calibrated using hydraulic head data in the transient inverse calibration process. The resulting models were good at fitting the hydraulic head data, but did not perform well in replicating the movement of contaminant plumes over this period. A separate transient inverse calibration effort used only tritium measurements collected from wells at the site over the operational period, along with estimates of the water volume and tritium mass discharged to the aquifer, to estimate the hydraulic properties. The resulting model did a better job of replicating the overall shape and development of the tritium plume, but did not do as well in matching the hydraulic heads. Both the hydraulic head and tritium concentration data sets were used jointly in the transient inverse process for this study. These data included 47,739 measurements of hydraulic head from 543 wells and 37,802 measurements of tritium concentrations from 1,201 wells. The transient inverse process estimated hydraulic conductivity for 18 facies-based zones in the main sand and gravel units in the unconfined aquifer. A simplified weighting scheme for the hydraulic head and tritium data was developed so that the overall sum-of-squared residuals for the inverse runs were roughly equally weighted for the two data sets. Preliminary simulation results from this combined calibration dataset show a good fit for both the evolving tritium plume and hydraulic head measurements over the operational period.
Thermal hydraulics modeling of the US Geological Survey TRIGA reactor
NASA Astrophysics Data System (ADS)
Alkaabi, Ahmed K.
The Geological Survey TRIGA reactor (GSTR) is a 1 MW Mark I TRIGA reactor located in Lakewood, Colorado. Single channel GSTR thermal hydraulics models built using RELAP5/MOD3.3, RELAP5-3D, TRACE, and COMSOL Multiphysics predict the fuel, outer clad, and coolant temperatures as a function of position in the core. The results from the RELAP5/MOD3.3, RELAP5-3D, and COMSOL models are similar. The TRACE model predicts significantly higher temperatures, potentially resulting from inappropriate convection correlations. To more accurately study the complex fluid flow patterns within the core, this research develops detailed RELAP5/MOD3.3 and COMSOL multichannel models of the GSTR core. The multichannel models predict lower fuel, outer clad, and coolant temperatures compared to the single channel models by up to 16.7°C, 4.8°C, and 9.6°C, respectively, as a result of the higher mass flow rates predicted by these models. The single channel models and the RELAP5/MOD3.3 multichannel model predict that the coolant temperatures in all fuel rings rise axially with core height, as the coolant in these models flows predominantly in the axial direction. The coolant temperatures predicted by the COMSOL multichannel model rise with core height in the B-, C-, and D-rings and peak and then decrease in the E-, F-, and G-rings, as the coolant tends to flow from the bottom sides of the core to the center of the core in this model. Experiments at the GSTR measured coolant temperatures in the GSTR core to validate the developed models. The axial temperature profiles measured in the GSTR show that the flow patterns predicted by the COMSOL multichannel model are consistent with the actual conditions in the core. Adjusting the RELAP5/MOD3.3 single and multichannel models by modifying the axial and cross-flow areas allow them to better predict the GSTR coolant temperatures; however, the adjusted models still fail to predict accurate axial temperature profiles in the E-, F-, and G-rings.
Studies investigate effects of hydraulic fracturing
NASA Astrophysics Data System (ADS)
Balcerak, Ernie
2012-11-01
The use of hydraulic fracturing, also known as fracking, to enhance the retrieval of natural gas from shale has been increasing dramatically—the number of natural gas wells rose about 50% since 2000. Shale gas has been hailed as a relatively low-cost, abundant energy source that is cleaner than coal. However, fracking involves injecting large volumes of water, sand, and chemicals into deep shale gas reservoirs under high pressure to open fractures through which the gas can travel, and the process has generated much controversy. The popular press, advocacy organizations, and the documentary film Gasland by Josh Fox have helped bring this issue to a broad audience. Many have suggested that fracking has resulted in contaminated drinking water supplies, enhanced seismic activity, demands for large quantities of water that compete with other uses, and challenges in managing large volumes of resulting wastewater. As demand for expanded domestic energy production intensifies, there is potential for substantially increased use of fracking together with other recovery techniques for "unconventional gas resources," like extended horizontal drilling.
Application study of magnetic fluid seal in hydraulic turbine
NASA Astrophysics Data System (ADS)
Yu, Z. Y.; Zhang, W.
2012-11-01
The waterpower resources of our country are abundant, and the hydroelectric power is developed, but at present the main shaft sealing device of hydraulic turbine is easy to wear and tear and the leakage is great. The magnetic fluid seal has the advantages of no contact, no wear, self-healing, long life and so on. In this paper, the magnetic fluid seal would be used in the main shaft of hydraulic turbine, the sealing structure was built the model, meshed the geometry, applied loads and solved by using MULTIPHYSICS in ANSYS software, the influence of the various sealing structural parameters such as tooth width, height, slot width, sealing gap on the sealing property were analyzed, the magnetic fluid sealing device suitable for large-diameter shaft and sealing water was designed, the sealing problem of the hydraulic turbine main shaft was solved effectively which will bring huge economic benefits.
Nonlinear mathematical modeling and sensitivity analysis of hydraulic drive unit
NASA Astrophysics Data System (ADS)
Kong, Xiangdong; Yu, Bin; Quan, Lingxiao; Ba, Kaixian; Wu, Liujie
2015-09-01
The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displacement changes of the piston are ignored, even experiment verification is not conducted. Therefore, in view of deficiencies above, a nonlinear mathematical model is established in this paper, including dynamic characteristics of servo valve, nonlinear characteristics of pressure-flow, initial displacement of servo cylinder piston and friction nonlinearity. The transfer function block diagram is built for the hydraulic drive unit closed loop position control, as well as the state equations. Through deriving the time-varying coefficient items matrix and time-varying free items matrix of sensitivity equations respectively, the expression of sensitivity equations based on the nonlinear mathematical model are obtained. According to structure parameters of hydraulic drive unit, working parameters, fluid transmission characteristics and measured friction-velocity curves, the simulation analysis of hydraulic drive unit is completed on the MATLAB/Simulink simulation platform with the displacement step 2 mm, 5 mm and 10 mm, respectively. The simulation results indicate that the developed nonlinear mathematical model is sufficient by comparing the characteristic curves of experimental step response and simulation step response under different constant load. Then, the sensitivity function time-history curves of seventeen parameters are obtained, basing on each state vector time-history curve of step response characteristic. The maximum value of displacement variation percentage and the sum of displacement variation absolute values in the sampling time are both taken as sensitivity indexes. The sensitivity indexes values above are calculated and shown visually in histograms under different working conditions, and change rules are analyzed. Then the sensitivity
Modeling of Propagation of Interacting Cracks Under Hydraulic Pressure Gradient
Huang, Hai; Mattson, Earl Douglas; Podgorney, Robert Karl
2015-04-01
A robust and reliable numerical model for fracture initiation and propagation, which includes the interactions among propagating fractures and the coupling between deformation, fracturing and fluid flow in fracture apertures and in the permeable rock matrix, would be an important tool for developing a better understanding of fracturing behaviors of crystalline brittle rocks driven by thermal and (or) hydraulic pressure gradients. In this paper, we present a physics-based hydraulic fracturing simulator based on coupling a quasi-static discrete element model (DEM) for deformation and fracturing with conjugate lattice network flow model for fluid flow in both fractures and porous matrix. Fracturing is represented explicitly by removing broken bonds from the network to represent microcracks. Initiation of new microfractures and growth and coalescence of the microcracks leads to the formation of macroscopic fractures when external and/or internal loads are applied. The coupled DEM-network flow model reproduces realistic growth pattern of hydraulic fractures. In particular, simulation results of perforated horizontal wellbore clearly demonstrate that elastic interactions among multiple propagating fractures, fluid viscosity, strong coupling between fluid pressure fluctuations within fractures and fracturing, and lower length scale heterogeneities, collectively lead to complicated fracturing patterns.
The EPA's Study on the Potential Impacts of Hydraulic Fracturing on Drinking Water Resources
NASA Astrophysics Data System (ADS)
Burden, Susan
2013-03-01
Natural gas plays a key role in our nation's clean energy future. The United States has vast reserves of natural gas that are commercially viable as a result of advances in horizontal drilling and hydraulic fracturing technologies, which enable greater access to gas in rock formations deep underground. These advances have spurred a significant increase in the production of both natural gas and oil across the country. However, as the use of hydraulic fracturing has increased, so have concerns about its potential human health and environmental impacts, especially for drinking water. In response to public concern, the US Congress requested that the US Environmental Protection Agency (EPA) conduct scientific research to examine the relationship between hydraulic fracturing and drinking water resources. In 2011, the EPA began research to assess the potential impacts of hydraulic fracturing on drinking water resources, if any, and to identify the driving factors that may affect the severity and frequency of such impacts. The study is organized around the five stages of the hydraulic fracturing water cycle, from water acquisition through the mixing of chemicals and the injection of fracturing fluid to post-fracturing treatment and/or disposal of wastewater. EPA scientists are using a transdisciplinary research approach involving laboratory studies, computer modeling, toxicity assessments, and case studies to answer research questions associated with each stage of the water cycle. This talk will provide an overview of the EPA's study, including a description of the hydraulic fracturing water cycle and a summary of the ongoing research projects.
Technical Review of the UNET2D Hydraulic Model
Perkins, William A.; Richmond, Marshall C.
2009-05-18
The Kansas City District of the US Army Corps of Engineers is engaged in a broad range of river management projects that require knowledge of spatially-varied hydraulic conditions such as velocities and water surface elevations. This information is needed to design new structures, improve existing operations, and assess aquatic habitat. Two-dimensional (2D) depth-averaged numerical hydraulic models are a common tool that can be used to provide velocity and depth information. Kansas City District is currently using a specific 2D model, UNET2D, that has been developed to meet the needs of their river engineering applications. This report documents a tech- nical review of UNET2D.
Sensitivity study on hydraulic well testing inversion using simulated annealing
Nakao, Shinsuke; Najita, J.; Karasaki, Kenzi
1997-11-01
For environmental remediation, management of nuclear waste disposal, or geothermal reservoir engineering, it is very important to evaluate the permeabilities, spacing, and sizes of the subsurface fractures which control ground water flow. Cluster variable aperture (CVA) simulated annealing has been used as an inversion technique to construct fluid flow models of fractured formations based on transient pressure data from hydraulic tests. A two-dimensional fracture network system is represented as a filled regular lattice of fracture elements. The algorithm iteratively changes an aperture of cluster of fracture elements, which are chosen randomly from a list of discrete apertures, to improve the match to observed pressure transients. The size of the clusters is held constant throughout the iterations. Sensitivity studies using simple fracture models with eight wells show that, in general, it is necessary to conduct interference tests using at least three different wells as pumping well in order to reconstruct the fracture network with a transmissivity contrast of one order of magnitude, particularly when the cluster size is not known a priori. Because hydraulic inversion is inherently non-unique, it is important to utilize additional information. The authors investigated the relationship between the scale of heterogeneity and the optimum cluster size (and its shape) to enhance the reliability and convergence of the inversion. It appears that the cluster size corresponding to about 20--40 % of the practical range of the spatial correlation is optimal. Inversion results of the Raymond test site data are also presented and the practical range of spatial correlation is evaluated to be about 5--10 m from the optimal cluster size in the inversion.
Modeling Reactor Coolant Systems Thermal-Hydraulic Transients
Energy Science and Technology Software Center (ESTSC)
1999-10-05
RELAP5/MOD3.2* is used to model reactor coolant systems during postulated accidents. The code models the coupled behavior of the reactor coolant system and the core for loss-of-coolant accidents and operational transients such as anticipated transients without scram, loss of offsite power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits simulating a variety of thermal-hydraulic systems. Control system and secondary system components are included to allow modeling of themore » plant controls, turbines, condensers, and secondary feedwater systems.« less
Advanced geothermal hydraulics model -- Phase 1 final report, Part 2
W. Zheng; J. Fu; W. C. Maurer
1999-07-01
An advanced geothermal well hydraulics model (GEODRIL) is being developed to accurately calculate bottom-hole conditions in these hot wells. In Phase 1, real-time monitoring and other improvements were added to GEODRIL. In Phase 2, GEODRIL will be integrated into Marconi's Intelligent Drilling Monitor (IDM) that will use artificial intelligence to detect lost circulation, fluid influxes and other circulation problems in geothermal wells. This software platform has potential for significantly reducing geothermal drilling costs.
NASA Astrophysics Data System (ADS)
Ghysels, Gert; Huysmans, Marijke
2015-04-01
Characterization of groundwater-surface water exchange fluxes is important for assessing riparian ecology, determining quantity and quality of pumped groundwater close to rivers, modeling groundwater flow, predicting flood peaks and low flows, and river water quality. The exchange fluxes between river and aquifer are strongly influenced by the hydraulic conductivity of the riverbed which can vary several orders of magnitude and shows a strong spatial variation. Direct measurement of riverbed hydraulic conductivity is cumbersome and therefore often indirect data such as temperature data or calibration of groundwater models are used to constrain riverbed hydraulic conductivity. In these approaches, the riverbed is usually represented as a homogeneous geological structure and the spatial variation of riverbed hydraulic conductivity is thus neglected. However, neglecting this spatial variation can lead to systematic underestimation of net river-aquifer exchange fluxes and may have important implications for the estimation of peak mass flows, for the hydrochemistry of streambed sediments, nutrient cycling and biogeochemical gradients. The MODFLOW software is the most wide-spread package used for groundwater modelling. In MODFLOW rivers are usually modelled using the River-package. However, in this package no distinction can be made between horizontal and vertical riverbed hydraulic conductivity and the riverbed cannot be subdivided into layers with different hydraulic characteristics. Riverbed sediments are strongly layered and thus another approach is advised. Different ways of introducing heterogeneous riverbeds in MODFLOW groundwater flow models are explored and compared. The influence of heterogeneous riverbeds on groundwater-surface water exchange fluxes is analyzed for two case studies: the Aa River in the Nete catchment and a stretch of the Dijle River near the nature reserve 'de Doode Bemde' (Belgium). For both cases fine-scale distributed local groundwater flow
Examining Alpine Meadow Restoration Techniques through Hydraulic Modeling
NASA Astrophysics Data System (ADS)
Moore, C. E.; Deems, J. S.; Loheide, S. P.; Lundquist, J. D.
2009-12-01
Stream stage controls groundwater fluctuations and, hence, vegetation communities more dramatically than stream discharge. A fluctuating climate has impacted stream discharge in the snow dominated alpine regions of the West Coast. Just how much these fluctuations in discharge affect stream stage depends on specific channel characteristics, such as slope, roughness, and morphology. Here we use a one dimensional hydraulic model, HEC-RAS (Hydraulic Engineering Center - River Analysis System) to model stream stage along the Tuolumne River, given a time series of stream discharge. Extensive hydroclimatic monitoring since 2001, and groundwater monitoring since 2006, make Tuolumne Meadows, in Yosemite National Park, California a unique location to test new model approaches, techniques, and linkages. Tuolumne, like many meadows in the Sierra, has seen degrading human impacts. In order to determine the most plausible, efficient and effective strategy of restoring impacted meadows, different management scenarios are modeled. Model results from different restoration scenarios, created by changing channel characteristics (ex: pristine conditions to vegetation enhancement or engineered structures, such as weirs or log jams), will demonstrate possible changes to river stage. HEC-RAS modeling then provides critical boundary conditions for groundwater table modeling. This paired technique of using surface and groundwater modeling to determine vegetation communities under future climates is easily transferable to other snow dominated basins and allows non-invasive investigation of multiple restoration options. Presented are the initial results of HEC-RAS modeling at Tuolumne Meadows as a collaborative project with surface, groundwater, and vegetation modeling in a changing climate.
Hydrological and hydraulic models for determination of flood-prone and flood inundation areas
NASA Astrophysics Data System (ADS)
Aksoy, Hafzullah; Sadan Ozgur Kirca, Veysel; Burgan, Halil Ibrahim; Kellecioglu, Dorukhan
2016-05-01
Geographic Information Systems (GIS) are widely used in most studies on water resources. Especially, when the topography and geomorphology of study area are considered, GIS can ease the work load. Detailed data should be used in this kind of studies. Because of, either the complication of the models or the requirement of highly detailed data, model outputs can be obtained fast only with a good optimization. The aim in this study, firstly, is to determine flood-prone areas in a watershed by using a hydrological model considering two wetness indexes; the topographical wetness index, and the SAGA (System for Automated Geoscientific Analyses) wetness index. The wetness indexes were obtained in the Quantum GIS (QGIS) software by using the Digital Elevation Model of the study area. Flood-prone areas are determined by considering the wetness index maps of the watershed. As the second stage of this study, a hydraulic model, HEC-RAS, was executed to determine flood inundation areas under different return period-flood events. River network cross-sections required for this study were derived from highly detailed digital elevation models by QGIS. Also river hydraulic parameters were used in the hydraulic model. Modelling technology used in this study is made of freely available open source softwares. Based on case studies performed on watersheds in Turkey, it is concluded that results of such studies can be used for taking precaution measures against life and monetary losses due to floods in urban areas particularly.
Floodplain Mapping Using Hydraulic Simulation Model in GIS
NASA Astrophysics Data System (ADS)
Salimi, Shokoufeh; Ghanbarpour, M. Reza; Solaimani, Karim; Ahmadi, Mirkhalegh Z.
In this research, a methodology was applied to integrate hydraulic simulation model, HEC-RAS and GIS analysis for delineation of flood extents and depths within a selected reach of Zaremroud River in Iran. Floodplain modeling is a recently new and applied method in river engineering discipline and is essential for prediction of flood hazards. It is necessary to simulate complicated hydraulic behavior of the river in a more simple way, for the purpose of managing and performing all river training practices. In this research, steady flow was simulated along 3 km end of Zaremroud River, upstream of the Tajan River in North of Iran. Floodplain zonation maps were derived using integrating of HEC-RAS and GIS analysis. Delineation of flood extents and depths within the floodplain were conducted in different return periods. Critical flooding area along the river was distinguished based on the grid layer of flood depths. The results indicated that hydraulic simulation by integrating with GIS analysis could be effective for various kinds of floodplain management and different scenarios for river training practices and flood mitigation planning.
Views on the future of thermal hydraulic modeling
Ishii, M.
1997-07-01
It is essential for the U.S. NRC to sustain the highest level of the thermal-hydraulics and reactor safety research expertise and continuously improve their accident analysis capability. Such expertise should span over four different areas which are strongly related to each other. These are: (1) Reactor Safety Code Development, (2) Two-phase Flow Modeling, (3) Instrumentation and Fundamental Experimental Research, and (4) Separate Effect and Integral Test. The NRC is already considering a new effort in the area of advanced thermal-hydraulics effort. Its success largely depends on the availability of a significantly improved two-phase flow formulation and constitutive relations supported by detailed experimental data. Therefore, it is recommended that the NRC start significant research efforts in the areas of two-phase flow modeling, instrumentation, basic and separate effect experiments which should be pursued systematically and with clearly defined objectives. It is desirable that some international program is developed in this area. This paper is concentrated on those items in the thermal-hydraulic area which eventually determine the quality of future accident analysis codes.
Statistical-physical model of the hydraulic conductivity
NASA Astrophysics Data System (ADS)
Usowicz, B.; Marczewski, W.; Usowicz, J. B.; Lukowski, M. I.
2012-04-01
The water content in unsaturated subsurface soil layer is determined by processes of exchanging mass and energy between media of soil and atmosphere, and particular members of layered media. Generally they are non-homogeneous on different scales, considering soil porosity, soil texture including presence of vegetation elements in the root zone, and canopy above the surface, and varying biomass density of plants above the surface in clusters. That heterogeneity determines statistically effective values of particular physical properties. This work considers mainly those properties which determine the hydraulic conductivity of soil. This property is necessary for characterizing physically water transfer in the root zone and access of nutrient matter for plants, but it also the water capacity on the field scale. The temporal variability of forcing conditions and evolutionarily changing vegetation causes substantial effects of impact on the water capacity in large scales, bringing the evolution of water conditions in the entire area, spanning a possible temporal state in the range between floods and droughts. The dynamic of this evolution of water conditions is highly determined by vegetation but is hardly predictable in evaluations. Hydrological models require feeding with input data determining hydraulic properties of the porous soil which are proposed in this paper by means of the statistical-physical model of the water hydraulic conductivity. The statistical-physical model was determined for soils being typical in Euroregion Bug, Eastern Poland. The model is calibrated on the base of direct measurements in the field scales, and enables determining typical characteristics of water retention by the retention curves bounding the hydraulic conductivity to the state of water saturation of the soil. The values of the hydraulic conductivity in two reference states are used for calibrating the model. One is close to full saturation, and another is for low water content far
Pen Branch Delta and Savannah River Swamp Hydraulic Model
Chen, K.F.
1999-05-13
The proposed Savannah River Site (SRS) Wetlands Restoration Project area is located in Barnwell County, South Carolina on the southwestern boundary of the SRS Reservation. The swamp covers about 40.5 km2 and is bounded to the west and south by the Savannah River and to the north and east by low bluffs at the edge of the Savannah River floodplain. Water levels within the swamp are determined by stage along the Savannah River, local drainage, groundwater seepage, and inflows from four tributaries, Beaver Dam Creek, Fourmile Branch, Pen Branch, and Steel Creek. Historic discharges of heated process water into these tributaries scoured the streambed, created deltas in the adjacent wetland, and killed native vegetation in the vicinity of the delta deposits. Future releases from these tributaries will be substantially smaller and closer to ambient temperatures. One component of the proposed restoration project will be to reestablish indigenous wetland vegetation on the Pen Branch delta that covers about 1.0 km2. Long-term predictions of water levels within the swamp are required to determine the characteristics of suitable plants. The objective of the study was to predict water levels at various locations within the proposed SRS Wetlands Restoration Project area for a range of Savannah River flows and regulated releases from Pen Branch. TABS-MD, a United States Army Corps of Engineer developed two-dimensional finite element open channel hydraulic computer code, was used to model the SRS swamp area for various flow conditions.
Technology Transfer Automated Retrieval System (TEKTRAN)
Concentrated flow is often the dominant source of water erosion following disturbance on rangeland. Because of the lack of studies that explain the hydraulics of concentrated flow on rangelands, cropland-based equations have typically been used for rangeland hydrology and erosion modeling, leading t...
NASA Astrophysics Data System (ADS)
Rosenzweig, Ravid; Furman, Alex; Dosoretz, Carlos; Shavit, Uri
2014-07-01
Biofilm effects on water flow in unsaturated environments have largely been ignored in the past. However, intensive engineered systems that involve elevated organic loads such as wastewater irrigation, effluent recharge, and bioremediation processes make understanding how biofilms affect flow highly important. In the current work, we present a channel-network model that incorporates water flow, substrate transport, and biofilm dynamics to simulate the alteration of soil hydraulic properties, namely water retention and conductivity. The change in hydraulic properties due to biofilm growth is not trivial and depends highly on the spatial distribution of the biofilm development. Our results indicate that the substrate mass transfer coefficient across the water-biofilm interface dominates the spatiotemporal distribution of biofilm. High mass transfer coefficients lead to uncontrolled biofilm growth close to the substrate source, resulting in preferential clogging of the soil. Low mass transfer coefficients, on the other hand, lead to a more uniform biofilm distribution. The first scenario leads to a dramatic reduction of the hydraulic conductivity with almost no change in water retention, whereas the second scenario has a smaller effect on conductivity but a larger influence on retention. The current modeling approach identifies key factors that still need to be studied and opens the way for simulation and optimization of processes involving significant biological activity in unsaturated soils.
Development of two-phase pipeline hydraulic analysis model based on Beggs-Brill correlation
NASA Astrophysics Data System (ADS)
Waluyo, Joko; Hermawan, Achilleus; Indarto
2016-06-01
The hydraulic analysis is an important stage in a reliable pipeline design. In the implementation, fluid distribution from a source to the sinks often occurs on parallel pipeline networks. The solution to the problem is complicated because of the iterative technique requirement. Regarding its solution effectiveness, there is a need for analysis related to the model and the solution method. This study aims to investigate pipeline hydraulic analysis on distributing of two-phase fluids flow. The model uses Beggs-Brill correlation to converse mass flow rates into pressure drops. In the solution technique, the Newton-Raphson iterative method is utilized. The iterative technique is solved using a computer program. The study is carried out using a certain pipeline network. The model is validated by comparing between Beggs-Brill towards Mukherjee-Brill correlation. The result reveals that the computer program enables solving of iterative calculation on the parallel pipeline hydraulic analysis. Convergence iteration is achieved by 50 iterations. The main results of the model are mass flow rate and pressure drop. The mass flow rate is obtained in the deviation up to 2.06%, between Beggs-Brill and Mukherjee-Brill correlation. On the other hand, the pressure gradient deviation is achieved on a higher deviation due to the different approach of the two correlations. The model can be further developed in the hydraulic pipeline analysis for two-phase flow.
Using Hydraulic Network Models to Teach Electric Circuit Principles
NASA Astrophysics Data System (ADS)
Jones, Irvin; EERC (Engineering Education Research Center) Collaboration
2013-11-01
Unlike other engineering disciplines, teaching electric circuit principles is difficult for some students because there isn't a visual context to rely on. So concepts such as electric potential, current, resistance, capacitance, and inductance have little meaning outside of their definition and the derived mathematical relationships. As a work in progress, we are developing a tool to support teaching, learning, and research of electric circuits. The tool will allow the user to design, build, and operate electric circuits in the form of hydraulic networks. We believe that this system will promote greater learning of electric circuit principles by visually realizing the conceptual and abstract concepts of electric circuits. Furthermore, as a teaching and learning tool, the hydraulic network system can be used to teach and improve comprehension of electrical principles in K through 12 classrooms and in cross-disciplinary environments such as Bioengineering, Mechanical Engineering, Industrial Engineering, and Aeronautical Engineering. As a research tool, the hydraulic network can model and simulate micro/nano bio-electro-chemical systems. Organization within the Swanson School of Engineering at the University of Pittsburgh.
Modelling of hydraulic fracture propagation in inhomogeneous poroelastic medium
NASA Astrophysics Data System (ADS)
Baykin, A. N.; Golovin, S. V.
2016-06-01
In the paper a model for description of a hydraulic fracture propagation in inhomogeneous poroelastic medium is proposed. Among advantages of the presented numerical algorithm, there are incorporation of the near-tip analysis into the general computational scheme, account for the rock failure criterion on the base of the cohesive zone model, possibility for analysis of fracture propagation in inhomogeneous reservoirs. The numerical convergence of the algorithm is verified and the agreement of our numerical results with known solutions is established. The influence of the inhomogeneity of the reservoir permeability to the fracture time evolution is also demonstrated.
Comparison of Soil Hydraulic Parameterizations for Mesoscale Meteorological Models.
NASA Astrophysics Data System (ADS)
Braun, Frank J.; Schädler, Gerd
2005-07-01
Soil water contents, calculated with seven soil hydraulic parameterizations, that is, soil hydraulic functions together with the corresponding parameter sets, are compared with observational data. The parameterizations include the Campbell/Clapp-Hornberger parameterization that is often used by meteorologists and the van Genuchten/Rawls-Brakensiek parameterization that is widespread among hydrologists. The observations include soil water contents at several soil depths and atmospheric surface data; they were obtained within the Regio Klima Projekt (REKLIP) at three sites in the Rhine Valley in southern Germany and cover up to 3 yr with 10-min temporal resolution. Simulations of 48-h episodes, as well as series of daily simulations initialized anew every 24 h and covering several years, were performed with the “VEG3D” soil-vegetation model in stand-alone mode; furthermore, 48-h episodes were simulated with the model coupled to a one-dimensional atmospheric model. For the cases and soil types considered in this paper, the van Genuchten/Rawls-Brakensiek model gives the best agreement between observed and simulated soil water contents on average. Especially during episodes with medium and high soil water content, the van Genuchten/Rawls-Brakensiek model performs better than the Campbell/Clapp-Hornberger model.
A Lagrangian Approach to Modelling Proppant Transport with Tip Screen-Out in KGD Hydraulic Fractures
NASA Astrophysics Data System (ADS)
Dontsov, E. V.; Peirce, A. P.
2015-11-01
This study introduces a continuum approach to model proppant transport in hydraulic fractures in a Lagrangian frame of reference. The model for the proppant transport is based on the recently obtained slurry flow solution inside a channel, where the latter utilizes a phenomenological constitutive relationship for a slurry. This approach allows us to describe the transition from Poiseuille flow with an effective viscosity to Darcy flow as the particle concentration increases towards the maximum value. The algorithm is presented for the one-dimensional case, for which propagation of a symmetric Kristinovich-Zheltov-Geertsma-De Klerk fracture is considered. To examine the effectiveness of the Lagrangian approach for proppant transport modelling, a set of parameters, for which proppant particles reach the fracture tip and cause the development of a proppant plug is selected. In this situation, the coupling between the hydraulic fracture propagation and proppant transport is the most significant. To estimate the accuracy of the Lagrangian proppant transport model, the results are compared to the predictions of an Eulerian proppant transport model, which utilizes the same algorithm for hydraulic fracture propagation. It is shown that, although both approaches have the same convergence rate, the error of the Lagrangian approach is three to five times smaller, which depends on the number of proppant elements used in the Lagrangian approach. This permits us to use a coarser mesh for hydraulic fracture propagation, and to obtain results with similar accuracy up to a hundred times faster.
McGraw, D.; Oberlander, P.
2007-12-18
The purpose of this study is to report on the results of a preliminary modeling framework to investigate the causes of the large hydraulic gradient north of Yucca Mountain. This study builds on the Saturated Zone Site-Scale Flow and Transport Model (referenced herein as the Site-scale model (Zyvoloski, 2004a), which is a three-dimensional saturated zone model of the Yucca Mountain area. Groundwater flow was simulated under natural conditions. The model framework and grid design describe the geologic layering and the calibration parameters describe the hydrogeology. The Site-scale model is calibrated to hydraulic heads, fluid temperature, and groundwater flowpaths. One area of interest in the Site-scale model represents the large hydraulic gradient north of Yucca Mountain. Nearby water levels suggest over 200 meters of hydraulic head difference in less than 1,000 meters horizontal distance. Given the geologic conceptual models defined by various hydrogeologic reports (Faunt, 2000, 2001; Zyvoloski, 2004b), no definitive explanation has been found for the cause of the large hydraulic gradient. Luckey et al. (1996) presents several possible explanations for the large hydraulic gradient as provided below: The gradient is simply the result of flow through the upper volcanic confining unit, which is nearly 300 meters thick near the large gradient. The gradient represents a semi-perched system in which flow in the upper and lower aquifers is predominantly horizontal, whereas flow in the upper confining unit would be predominantly vertical. The gradient represents a drain down a buried fault from the volcanic aquifers to the lower Carbonate Aquifer. The gradient represents a spillway in which a fault marks the effective northern limit of the lower volcanic aquifer. The large gradient results from the presence at depth of the Eleana Formation, a part of the Paleozoic upper confining unit, which overlies the lower Carbonate Aquifer in much of the Death Valley region. The
Thermal hydraulic modeling of integrated cooling water systems
Niyogi, K.K.; Rathi, J.S.; Phan, T.Q.; Chaudhary, A.
1994-12-31
Thermal hydraulic modeling of cooling water systems has been extended to multiple system configurations with heat exchangers as interface components between systems. The computer program PC-TRAX has been used as the basic tool for the system simulation. Additional heat exchanger modules have been incorporated to accurately predict the thermal performance of systems for the design as well as off-design conditions. The modeling accommodates time-dependent changes in conditions, temperature and pressure controllers, and detailed physical parameters of the heat exchangers. The modeling has been illustrated with examples from actual plant systems. An integrated system consisting of Spent Fuel Pool, Primary Component Cooling Water, and Service Water System has been successfully modeled to predict their performance under normal operations and emergency conditions. System configurations are changed from the base model by using a command module.
Discrete modeling of hydraulic fracturing processes in a complex pre-existing fracture network
NASA Astrophysics Data System (ADS)
Kim, K.; Rutqvist, J.; Nakagawa, S.; Houseworth, J. E.; Birkholzer, J. T.
2015-12-01
Hydraulic fracturing and stimulation of fracture networks are widely used by the energy industry (e.g., shale gas extraction, enhanced geothermal systems) to increase permeability of geological formations. Numerous analytical and numerical models have been developed to help understand and predict the behavior of hydraulically induced fractures. However, many existing models assume simple fracturing scenarios with highly idealized fracture geometries (e.g., propagation of a single fracture with assumed shapes in a homogeneous medium). Modeling hydraulic fracture propagation in the presence of natural fractures and homogeneities can be very challenging because of the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and pre-existing natural fractures. In this study, the TOUGH-RBSN code for coupled hydro-mechanical modeling is utilized to simulate hydraulic fracture propagation and its interaction with pre-existing fracture networks. The simulation tool combines TOUGH2, a simulator of subsurface multiphase flow and mass transport based on the finite volume approach, with the implementation of a lattice modeling approach for geomechanical and fracture-damage behavior, named Rigid-Body-Spring Network (RBSN). The discrete fracture network (DFN) approach is facilitated in the Voronoi discretization via a fully automated modeling procedure. The numerical program is verified through a simple simulation for single fracture propagation, in which the resulting fracture geometry is compared to an analytical solution for given fracture length and aperture. Subsequently, predictive simulations are conducted for planned laboratory experiments using rock-analogue (soda-lime glass) samples containing a designed, pre-existing fracture network. The results of a preliminary simulation demonstrate selective fracturing and fluid infiltration along the pre-existing fractures, with additional fracturing in part
Use of hydraulic models to identify and resolve design isssues in FGD systems
Strock, T.W.; Gohara, W.F.
1995-06-01
The hydraulics within a wet flue gas desulfurization (FGD) scrubber involve several complex two-phase gas/liquid interactions that directly affect the scrubber pressure drop, mist elimination efficiency, and the mass transfer process of SO{sub 2} removal. Current industrial efforts to develop cost effective, high-efficiency wet FGD scrubbers are focusing, in part, on the hydraulics. The development of an experimental approach and test facility for understanding and optimizing wet scrubber flow characteristics has been completed. Hydraulic models simulate full-scale units and allow the designer to view the gas/liquid flow interactions. Modeling procedures for downsizing the wet scrubber for the laboratory have been developed and validated with field data comparisons. A one-eighth scale hydraulic model has been used to study several FGD scrubber design issues. Design changes to reduce capital and operating cost have been developed and tested. Recently, the model was used to design a commercial, uniform flow, high gas velocity absorber for the next generation of FGD systems.
Vrugt, Jasper A; Wohling, Thomas
2008-01-01
Most studies in vadose zone hydrology use a single conceptual model for predictive inference and analysis. Focusing on the outcome of a single model is prone to statistical bias and underestimation of uncertainty. In this study, we combine multi-objective optimization and Bayesian Model Averaging (BMA) to generate forecast ensembles of soil hydraulic models. To illustrate our method, we use observed tensiometric pressure head data at three different depths in a layered vadose zone of volcanic origin in New Zealand. A set of seven different soil hydraulic models is calibrated using a multi-objective formulation with three different objective functions that each measure the mismatch between observed and predicted soil water pressure head at one specific depth. The Pareto solution space corresponding to these three objectives is estimated with AMALGAM, and used to generate four different model ensembles. These ensembles are post-processed with BMA and used for predictive analysis and uncertainty estimation. Our most important conclusions for the vadose zone under consideration are: (1) the mean BMA forecast exhibits similar predictive capabilities as the best individual performing soil hydraulic model, (2) the size of the BMA uncertainty ranges increase with increasing depth and dryness in the soil profile, (3) the best performing ensemble corresponds to the compromise (or balanced) solution of the three-objective Pareto surface, and (4) the combined multi-objective optimization and BMA framework proposed in this paper is very useful to generate forecast ensembles of soil hydraulic models.
NASA Astrophysics Data System (ADS)
Stucchi Boschi, Raquel; Qin, Mingming; Gimenez, Daniel; Cooper, Miguel
2016-04-01
Modeling is an important tool for better understanding and assessing land use impacts on landscape processes. A key point for environmental modeling is the knowledge of soil hydraulic properties. However, direct determination of soil hydraulic properties is difficult and costly, particularly in vast and remote regions such as one constituting the Amazon Biome. One way to overcome this problem is to extrapolate accurately estimated data to pedologically similar sites. The van Genuchten (VG) parametric equation is the most commonly used for modeling SWRC. The use of a Bayesian approach in combination with the Markov chain Monte Carlo to estimate the VG parameters has several advantages compared to the widely used global optimization techniques. The Bayesian approach provides posterior distributions of parameters that are independent from the initial values and allow for uncertainty analyses. The main objectives of this study were: i) to estimate hydraulic parameters from data of pasture and forest sites by the Bayesian inverse modeling approach; and ii) to investigate the extrapolation of the estimated VG parameters to a nearby toposequence with pedologically similar soils to those used for its estimate. The parameters were estimated from volumetric water content and tension observations obtained after rainfall events during a 207-day period from pasture and forest sites located in the southeastern Amazon region. These data were used to run HYDRUS-1D under a Differential Evolution Adaptive Metropolis (DREAM) scheme 10,000 times, and only the last 2,500 times were used to calculate the posterior distributions of each hydraulic parameter along with 95% confidence intervals (CI) of volumetric water content and tension time series. Then, the posterior distributions were used to generate hydraulic parameters for two nearby toposequences composed by six soil profiles, three are under forest and three are under pasture. The parameters of the nearby site were accepted when
Modeling of Interaction of Hydraulic Fractures in Complex Fracture Networks
NASA Astrophysics Data System (ADS)
Kresse, O. 2; Wu, R.; Weng, X.; Gu, H.; Cohen, C.
2011-12-01
A recently developed unconventional fracture model (UFM) is able to simulate complex fracture network propagation in a formation with pre-existing natural fractures. Multiple fracture branches can propagate at the same time and intersect/cross each other. Each open fracture exerts additional stresses on the surrounding rock and adjacent fractures, which is often referred to as "stress shadow" effect. The stress shadow can cause significant restriction of fracture width, leading to greater risk of proppant screenout. It can also alter the fracture propagation path and drastically affect fracture network patterns. It is hence critical to properly model the fracture interaction in a complex fracture model. A method for computing the stress shadow in a complex hydraulic fracture network is presented. The method is based on an enhanced 2D Displacement Discontinuity Method (DDM) with correction for finite fracture height. The computed stress field is compared to 3D numerical simulation in a few simple examples and shows the method provides a good approximation for the 3D fracture problem. This stress shadow calculation is incorporated in the UFM. The results for simple cases of two fractures are presented that show the fractures can either attract or expel each other depending on their initial relative positions, and compares favorably with an independent 2D non-planar hydraulic fracture model. Additional examples of both planar and complex fractures propagating from multiple perforation clusters are presented, showing that fracture interaction controls the fracture dimension and propagation pattern. In a formation with no or small stress anisotropy, fracture interaction can lead to dramatic divergence of the fractures as they tend to repel each other. However, when stress anisotropy is large, the fracture propagation direction is dominated by the stress field and fracture turning due to fracture interaction is limited. However, stress shadowing still has a strong effect
A Simple Hydraulic Analog Model of Oxidative Phosphorylation.
Willis, Wayne T; Jackman, Matthew R; Messer, Jeffrey I; Kuzmiak-Glancy, Sarah; Glancy, Brian
2016-06-01
Mitochondrial oxidative phosphorylation is the primary source of cellular energy transduction in mammals. This energy conversion involves dozens of enzymatic reactions, energetic intermediates, and the dynamic interactions among them. With the goal of providing greater insight into the complex thermodynamics and kinetics ("thermokinetics") of mitochondrial energy transduction, a simple hydraulic analog model of oxidative phosphorylation is presented. In the hydraulic model, water tanks represent the forward and back "pressures" exerted by thermodynamic driving forces: the matrix redox potential (ΔGredox), the electrochemical potential for protons across the mitochondrial inner membrane (ΔGH), and the free energy of adenosine 5'-triphosphate (ATP) (ΔGATP). Net water flow proceeds from tanks with higher water pressure to tanks with lower pressure through "enzyme pipes" whose diameters represent the conductances (effective activities) of the proteins that catalyze the energy transfer. These enzyme pipes include the reactions of dehydrogenase enzymes, the electron transport chain (ETC), and the combined action of ATP synthase plus the ATP-adenosine 5'-diphosphate exchanger that spans the inner membrane. In addition, reactive oxygen species production is included in the model as a leak that is driven out of the ETC pipe by high pressure (high ΔGredox) and a proton leak dependent on the ΔGH for both its driving force and the conductance of the leak pathway. Model water pressures and flows are shown to simulate thermodynamic forces and metabolic fluxes that have been experimentally observed in mammalian skeletal muscle in response to acute exercise, chronic endurance training, and reduced substrate availability, as well as account for the thermokinetic behavior of mitochondria from fast- and slow-twitch skeletal muscle and the metabolic capacitance of the creatine kinase reaction. PMID:26807634
75 FR 35023 - Informational Public Meetings for Hydraulic Fracturing Research Study
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-21
... AGENCY Informational Public Meetings for Hydraulic Fracturing Research Study AGENCY: Environmental... between hydraulic fracturing and drinking water. The meetings are open to all interested parties and will... Hydraulic Fracturing Study informational meetings are as follows: July 8, 2010, from 6 p.m. to 10 p.m.,...
NASA Astrophysics Data System (ADS)
Lu, C.; Zhang, Y.; Shu, L.; Chen, X.; Chen, S.; Li, S.; Wang, G.; Li, J.
2015-05-01
The paper aims to evaluate the impacts of the average hydraulic conductivity of the heterogeneous aquifer on the estimated hydraulic conductivity using the observations from pumping tests. The results of aquifer tests conducted at a karst aquifer are first introduced. A MODFLOW groundwater flow model was developed to perform numerical pumping tests, and the heterogeneous hydraulic conductivity (K) field was generated using the Monte Carlo method. The K was estimated by the Theis solution for an unconfined aquifer. The effective hydraulic conductivity (Ke) was calculated to represent the hydraulic conductivity of a heterogeneous aquifer. The results of numerical simulations demonstrate that Ke increase with the mean of hydraulic conductivity (EK), and decrease with the coefficient of variation of the hydraulic conductivity (Cv). The impact of spatial variability of K on the estimated Ke at two observation wells with smaller EK is less significant compared to the cases with larger EK.
Sánchez, F; Viedma, A; Kaiser, A S
2016-09-15
Fluid dynamic behaviour plays an important role in wastewater treatment. An efficient treatment requires the inexistence of certain hydraulic problems such as dead zones or short-circuiting flows. Residence time distribution (RTD) analysis is an excellent technique for detecting these inefficiencies. However, many wastewater treatment installations include water or sludge recycling systems, which prevent us from carrying out a conventional tracer pulse experiment to obtain the RTD curve of the installation. This paper develops an RTD analysis of an activated sludge reactor with recycling system. A tracer experiment in the reactor is carried out. Three analytical models, derived from the conventional pulse model, are proposed to obtain the RTD curve of the reactor. An analysis of the results is made, studying which model is the most suitable for each situation. This paper is useful to analyse the hydraulic efficiency of reactors with recycling systems. PMID:27288672
Modeling hydraulic regenerative hybrid vehicles using AMESim and Matlab/Simulink
NASA Astrophysics Data System (ADS)
Lynn, Alfred; Smid, Edzko; Eshraghi, Moji; Caldwell, Niall; Woody, Dan
2005-05-01
This paper presents the overview of the simulation modeling of a hydraulic system with regenerative braking used to improve vehicle emissions and fuel economy. Two simulation software packages were used together to enhance the simulation capability for fuel economy results and development of vehicle and hybrid control strategy. AMESim, a hydraulic simulation software package modeled the complex hydraulic circuit and component hardware and was interlinked with a Matlab/Simulink model of the vehicle, engine and the control strategy required to operate the vehicle and the hydraulic hybrid system through various North American and European drive cycles.
HYDRAULIC STUDIES AND CLEANING EVALUATIONS OF ULTRAVIOLET DISINFECTION UNITS
Various types of operating ultraviolet disinfection reactor designs were evaluated for hydraulic characteristics and cleaning requirements. The fluorocarbon polymer tube designs promote plug-flow behavior because of their relatively high length-to-diameter ratio. Hydraulic evalua...
Assessing the Transferability of Hydraulic Habitat Models for Atlantic Salmon Fry
NASA Astrophysics Data System (ADS)
Millidine, K. J.; Malcolm, I.; Fryer, R. J.
2015-12-01
Hydraulic habitat models, which are logistically and technically challenging and expensive to produce, are frequently transferred between rivers without validation. Although this is known to be associated with problems, few studies have assessed the potential consequences for model predictions. This study investigated the local (within sub-catchment) transfer of hydraulic habitat models developed for Atlantic salmon (Salmo salar) fry. Detailed 2D hydraulic models were developed for two adjacent reaches, each containing pool, riffle, glide and run habitats where salmon fry were stocked at uniform saturated densities. Substrate and cover were characterised using transects. Generalised Additive Models (GAM's) were fitted to seasonal fry abundance data, with Froude number, dominant substrate and cover included as predictor variables. Despite attempts to select reaches with similar characteristics, the spatial distribution of Froude, dominant substrate and cover differed, with substrate and cover exhibiting the greatest inter-reach differences. Froude was the most important individual predictor of fry abundance, with the highest densities observed at moderate Froude across all seasons. When transferred between reaches, models which contained multiple predictor variables and their interactions transferred less well than models containing Froude alone potentially reflecting inter-reach differences in the distribution of substrate and cover. This study suggests that (1) habitat models should be developed at sites offering maximum environmental complexity at a local level (2) scientists and managers should avoid transferring models between locations with different environmental characteristics, especially in the absence of model validation (3) complex models should be avoided (4) the transferability of Froude only models should be further investigated, if predictions of habitat quality are to be made at new sites.
NASA Astrophysics Data System (ADS)
Huisman, J. A.; Rings, J.; Vrugt, J. A.; Sorg, J.; Vereecken, H.
2010-01-01
SummaryCoupled hydrogeophysical inversion aims to improve the use of geophysical data for hydrological model parameterization. Several numerical studies have illustrated the feasibility and advantages of a coupled approach. However, there is still a lack of studies that apply the coupled inversion approach to actual field data. In this paper, we test the feasibility of coupled hydrogeophysical inversion for determining the hydraulic properties of a model dike using measurements of electrical resistance tomography (ERT). Our analysis uses a two-dimensional (2D) finite element hydrological model (HYDRUS-2D) coupled to a 2.5D finite element electrical resistivity code (CRMOD), and includes explicit recognition of parameter uncertainty by using a Bayesian and multiple criteria framework with the DREAM and AMALGAM population based search algorithms. To benchmark our inversion results, soil hydraulic properties determined from ERT data are compared with those separately obtained from detailed in situ soil water content measurements using Time Domain Reflectometry (TDR). Our most important results are as follows. (1) TDR and ERT data theoretically contain sufficient information to resolve most of the soil hydraulic properties, (2) the DREAM-derived posterior distributions of the hydraulic parameters are quite similar when estimated separately using TDR and ERT measurements for model calibration, (3) among all parameters, the saturated hydraulic conductivity of the dike material is best constrained, (4) the saturation exponent of the petrophysical model is well defined, and matches independently measured values, (5) measured ERT data sufficiently constrain model predictions of water table dynamics within the model dike. This finding demonstrates an innate ability of ERT data to provide accurate hydrogeophysical parameterizations for flooding events, which is of particular relevance to dike management, and (6) the AMALGAM-derived Pareto front demonstrates trade-off in the
Modelling the hydraulic and geochemical evolution of a hillslope transect
NASA Astrophysics Data System (ADS)
Maier, U.; Flegr, M.; Rügner, H.; Selle, B.; Grathwohl, P.
2012-04-01
Especially the long-term response of subsurface-surface water systems and their water quality on pressures such as diffuse pollution from agricultural activities or atmospheric deposition are not well understood at the catchment scale. The goal of this study is to quantify the factors that influence the geochemical interactions between soils and seepage water, rocks and groundwater as well the subsurface/surface water interaction in a comprehensive way. Only very recently numerical codes became available which not only couple flow and reactive transport but also unsaturated and saturated zone allowing it to follow water quality from the infiltration into the soil until discharge in a spring or river, of which the code MIN3P was used. In order to assess groundwater quality, the origin of the recharge water and its chemical evolution pathways have to be known at the catchment scale. This involves precipitation, evapotranspiration, seepage water infiltration, interflow and perched water, and finally the groundwater as well as the effluent to rivers or springs in a watershed. Water quality is affected by rainwater pH and dissolved solids, leaching of potential pollutants from top-soils, release of CO2 from organic matter oxidation / microbial respiration in the unsaturated zone and the water-rock interaction in the subsurface. A promising approach to identify the principal processes is the selection of vertical profiles along streamlines across the area of interest. That way, numerical simulations requiring only short computational time could be utilized to describe the water flow and solute transport from elevated parts of the catchment to the receiving stream. Eventually, this approach can be extended to capture a watershed in a three-dimensional model. A geologic model representative for a typical valley scenario in a triassic sequence landscape composed of of sandstones and marls was set up, consisting of a i) sand and gravel aquifer, ii) underlying and hill
Empirical correction of a hydraulic model of mutual displacement of solutions in a porous bed
Kuznetskii, R.S.
1986-09-10
This paper seeks to experimentally verify and correct a hydraulic model for the forced displacement of a homogeneous solution in a porous or granular bed or in a pipe by calibrating the model against empirical data for the hydraulic behavior of sulfonated coal, an ion exchange resin, calcium carbonate, and diatomite.
NASA Astrophysics Data System (ADS)
Kanno, C.; McLaughlin, M.; Blotevogel, J.; Benson, D. A.; Borch, T.; McCray, J. E.
2015-12-01
Hydraulic fracturing has revolutionized the U.S.'s energy portfolio by making shale reservoirs productive and commercially viable. However, the public is concerned that the chemical constituents in hydraulic fracturing fluid, produced water, or natural gas itself could potentially impact groundwater or adjacent streams. Here, we conduct fate and transport simulations of surface spills, the most likely contamination pathway to occur during oil and gas production operations, to evaluate whether or not these spills pose risks to groundwater quality. We focus on the South Platte Alluvial Aquifer, which is located in the greater Denver metro area and overlaps a zone of high-density oil and gas development. The purpose of this work is to assess the mobility and persistence of chemical contaminants (e.g. biocides, friction reducers, surfactants, hydrocarbons, etc.) —based on sorption to soil, degradation potential, co-contaminant interactions, and spill conditions—and to understand the site characteristics and hydrologic conditions that would make a particular location prone to groundwater quality degradation in the event of an accidental release. We propose a coupled analytical-numerical approach that could be duplicated by environmental consultants. Results suggest that risk of groundwater pollution, based on predicted concentration at the groundwater table, is low in most areas of the South Platte system for the contaminants investigated under common spill conditions. However, substantial risk may exist in certain areas where the groundwater table is shallow. In addition, transport of certain contaminants is influenced by interactions with other constituents in produced or stimulation fluids. By helping to identify locations in the Front Range of Colorado that are at low or high risk for groundwater contamination due to a surface spill, it is our hope that this work will aid in improving prevention, mitigation, and remediation practices so that decision-makers can
A Bayesian Chance-Constrained Method for Hydraulic Barrier Design Under Model Structure Uncertainty
NASA Astrophysics Data System (ADS)
Chitsazan, N.; Pham, H. V.; Tsai, F. T. C.
2014-12-01
The groundwater community has widely recognized the model structure uncertainty as the major source of model uncertainty in groundwater modeling. Previous studies in the aquifer remediation design, however, rarely discuss the impact of the model structure uncertainty. This study combines the chance-constrained (CC) programming with the Bayesian model averaging (BMA) as a BMA-CC framework to assess the effect of model structure uncertainty in the remediation design. To investigate the impact of the model structure uncertainty on the remediation design, we compare the BMA-CC method with the traditional CC programming that only considers the model parameter uncertainty. The BMA-CC method is employed to design a hydraulic barrier to protect public supply wells of the Government St. pump station from saltwater intrusion in the "1,500-foot" sand and the "1-700-foot" sand of the Baton Rouge area, southeastern Louisiana. To address the model structure uncertainty, we develop three conceptual groundwater models based on three different hydrostratigraphy structures. The results show that using the traditional CC programming overestimates design reliability. The results also show that at least five additional connector wells are needed to achieve more than 90% design reliability level. The total amount of injected water from connector wells is higher than the total pumpage of the protected public supply wells. While reducing injection rate can be achieved by reducing reliability level, the study finds that the hydraulic barrier design to protect the Government St. pump station is not economically attractive.
Effects of various parameters on hydraulic fracture geometry. [Theoretical and experimental studies
Hanson, M.E.; Anderson, G.D.; Shaffer, R.J.
1980-05-01
Theoretical models have been applied to analyze some aspects of the dynamics of fracturing near material interfaces. Results of these calculations indicate that variation of material properties across a well bonded interface can cause dynamic material response resulting from the fracturing which could enhance propagation across the inerface. Effects of friction have also been analyzed theoretically; however, in the frictional calculations the wave mechanics have been ignored. These calculations have shown that frictional slip along the interface tends to draw a pressurized fracture toward the interface; this motion tends to reduce the chances of penetrating the material across the frictional interface. Small scale laboratory experiments are performed to study the effects of frictional characteristics on hydraulic fracture growth across unbonded interfaces in rocks. Various lubricants and mechanical preparation of the interface surfaces are used to vary the coefficients of friction on the interface surfaces. It is found that the frictional shear stress that the interface surface can support determines whether a hydraulically driven crack will cross the interface. Pre-existing cracks impede the propagation of the hydraulic fracture across the interface. These experimental results on the effects of friction on the interface and the effects of pre-existing cracks on hydraulic fracture penetration of interfaces are consistent with the predictions of the numerical model calculations. 11 figures.
Modelling Subduction Zone Magmatism Due to Hydraulic Fracture
NASA Astrophysics Data System (ADS)
Lawton, R.; Davies, J. H.
2014-12-01
The aim of this project is to test the hypothesis that subduction zone magmatism involves hydraulic fractures propagating from the oceanic crust to the mantle wedge source region (Davies, 1999). We aim to test this hypothesis by developing a numerical model of the process, and then comparing model outputs with observations. The hypothesis proposes that the water interconnects in the slab following an earthquake. If sufficient pressure develops a hydrofracture occurs. The hydrofracture will expand in the direction of the least compressive stress and propagate in the direction of the most compressive stress, which is out into the wedge. Therefore we can calculate the hydrofracture path and end-point, given the start location on the slab and the propagation distance. We can therefore predict where water is added to the mantle wedge. To take this further we have developed a thermal model of a subduction zone. The model uses a finite difference, marker-in-cell method to solve the heat equation (Gerya, 2010). The velocity field was prescribed using the analytical expression of cornerflow (Batchelor, 1967). The markers contained within the fixed grid are used to track the different compositions and their properties. The subduction zone thermal model was benchmarked (Van Keken, 2008). We used the hydrous melting parameterization of Katz et.al., (2003) to calculate the degree of melting caused by the addition of water to the wedge. We investigate models where the hydrofractures, with properties constrained by estimated water fluxes, have random end points. The model predicts degree of melting, magma productivity, temperature of the melt and water content in the melt for different initial water fluxes. Future models will also include the buoyancy effect of the melt and residue. Batchelor, Cambridge UP, 1967. Davies, Nature, 398: 142-145, 1999. Gerya, Cambridge UP, 2010. Katz, Geochem. Geophys. Geosy, 4(9), 2003 Van Keken et.al. Phys. Earth. Planet. In., 171:187-197, 2008.
Modelling of of hydraulic fractures trajectories in inhomogeneous stress field
NASA Astrophysics Data System (ADS)
Andreev, A. A.; Galybin, A.
2013-05-01
The paper examines an actual problem of oil and gas production -- modelling of the hydro-fracture trajectories depending on ihomogeneous distributions of pore pressure. The results could serve for improvement of the design of hydraulic fracturing in the oil/gas fields. The methods of the plane elasticity theory and fracture mechanics are employed. It is assumed, that in addition to the homogeneous field of natural stress the reservoir is also subjected to additional stresses caused by technological reasons, which makes the total stress field to be inhomogeneous. Therefore, the objective is to model a curvilinear crack path in an elastic inhomogeneous-loaded plane depending on the different mechanical parameters that control the stress state of the reservoir. For the simulation of the trajectory of a crack the method of boundary integral equation is used. The algorithms of step-by-step determination of the crack's trajectory development using the criterion of maximum tensile stresses at the end of the cracks have been developed. For the numerical realization of the solution we used a special modification of the method of mechanical quadratures providing effective and fast solution of the corresponding system of singular integral equation. The solution for the hydro-fracture path have been simulated for the case of inhomogeneous stress field due to presence of injection well for several physical models.
Sensitivity study on hydraulic well testing inversion using simulated annealing
Nakao, Shinsuke; Najita, J.; Karasaki, Kenzi
1999-10-01
Cluster variable aperture (CVA) simulated annealing has been used as an inversion technique to construct fluid flow models of fractured formations based on transient pressure data from hydraulic tests. A two-dimensional fracture network system is represented as a filled regular lattice of fracture elements. The algorithm iteratively changes element apertures for a cluster of fracture elements in order to improve the match to observed pressure transients. Aperture size is chosen randomly from a list of discrete apertures. The cluster size is held constant through the iterations. Since hydraulic inversion is inherently nonunique, it is important to use additional information. The authors investigated the relationship between the scale of heterogeneity and the optimal cluster size and shape to enhance convergence of the inversion and improve the results. In a spatially correlated transmissivity field, a cluster size corresponding to about 20% to 40% of the practical range of the spatial correlation is optimal. Inversion results of the Raymond test site data are also presented and based on an optimal cluster size; the practical range of the spatial correlation is estimated to be 5 to 10 m.
Sensitivity study on hydraulic well testing inversion using simulated annealing.
Nakao, S; Najita, J; Karasaki, K
1999-01-01
Cluster variable aperture (CVA) simulated annealing has been used as an inversion technique to construct fluid flow models of fractured formations based on transient pressure data from hydraulic tests. A two-dimensional fracture network system is represented as a filled regular lattice of fracture elements. The algorithm iteratively changes element apertures for a cluster of fracture elements in order to improve the match to observed pressure transients. Aperture size is chosen randomly from a list of discrete apertures. The cluster size is held constant throughout the iterations. Since hydraulic inversion is inherently nonunique, it is important to use additional information. We investigated the relationship between the scale of heterogeneity and the optimal cluster size and shape to enhance convergence of the inversion and improve the results. In a spatially correlated transmissivity field, a cluster size corresponding to about 20 % to 40 % of the practical range of the spatial correlation is optimal. Inversion results of the Raymond test site data are also presented and based on an optimal cluster size; the practical range of the spatial correlation is estimated to be 5 to 10 m. PMID:19125927
NASA Astrophysics Data System (ADS)
Yan, Binyan; Dickinson, Robert E.
2014-11-01
Hydraulic redistribution is the process of soil water transport through the low-resistance pathway provided by plant roots. It has been observed in field studies and proposed to be one of the processes that enable the Amazon rainforest to resist periodical dry spells without experiencing water limitations. How and to what extent hydraulic redistribution may increase vegetation resistance to longer or more severe droughts than seasonal dryness have not been investigated yet, which is the focus of this study. The artificially prolonged drought produced by the rainfall exclusion experiment is used as an example of long drought, and the 2005 drought is used as a severe drought. The parameterization of hydraulic redistribution proposed by Ryel et al. (2002) was incorporated into the Community Land Model version 4 (CLM4). Three paired numerical experiments were conducted, one set using the default model (CTL) and the other using the model with considerations of hydraulic redistribution (HR). Results show that the vegetation response (including evapotranspiration, biomass, and leaf area index (LAI)) to dryness of all the three types is better captured with hydraulic redistribution incorporated. Plants are more resistant to dryness when hydraulic redistribution increases plant water availability and thus facilitates their growth. When a drought is long lasting, the vegetation response is delayed by hydraulic redistribution. Therefore, if a drought ends earlier than permanent damage is made, the magnitude of vegetation response will be lowered by this mechanism, i.e., the vegetation will be more resistant to dryness.
Hydraulic Fracturing and Drinking Water Resources: Update on EPA Hydraulic Fracturing Study
Natural gas plays a key role in our nation's energy future and the process known as hydraulic fracturing (HF) is one way of accessing that resource. Over the past few years, several key technical, economic, and energy developments have spurred increased use of HF for gas extracti...
NASA Astrophysics Data System (ADS)
He, C.; Jin, X.; Zhang, L.; Zhang, X.
2014-12-01
Heterogeneity of soil hydraulic properties directly affects variations of hydrological processes at corresponding scales. Understanding spatial variation of soil hydraulic properties such as soil moisture is therefore fundamental for modeling watershed ecohydrological processes. As part of the National Science Foundation of China (NSFC) funded ''Integrated Ecohydrological Research Plan of the Heihe River Watershed'', this study established an observation network that consists of sampling points, zones, and tributaries to analyze spatial variations of soil hydraulic properties in the Upper Reach of the Heihe River Watershed, a second largest inland river (terminal lake) with a drainage area of over 128,000 km2 in Northwest China. Spatial heterogeneity of soil properties was analyzed based on the large number of soil sampling and in situ observations. The spatial clustering method, Full-Order-CLK was employed to derive five soil heterogeneous zones (Configuration 97, 80, 65, 40, and 20). Subsequently, SWAT model was used to quantify the impact of the spatial heterogeneity of soil hydraulic properties on hydrologic process in the study watershed. Results show the simulations by the SWAT model with the spatially clustered soil hydraulic information from the field sampling data had much better representation of the soil heterogeneity and more accurate performance than the model using the average soil property values for each soil type derived from the coarse soil datasets (Gansu Soil Handbook at 1:1,000,000 scale). Thus, incorporating detailed field sampling soil heterogeneity data greatly improves performance in hydrologic modeling.
Jones, Joseph L.; Fulford, Janice M.; Voss, Frank D.
2002-01-01
A system of numerical hydraulic modeling, geographic information system processing, and Internet map serving, supported by new data sources and application automation, was developed that generates inundation maps for forecast floods in near real time and makes them available through the Internet. Forecasts for flooding are generated by the National Weather Service (NWS) River Forecast Center (RFC); these forecasts are retrieved automatically by the system and prepared for input to a hydraulic model. The model, TrimR2D, is a new, robust, two-dimensional model capable of simulating wide varieties of discharge hydrographs and relatively long stream reaches. TrimR2D was calibrated for a 28-kilometer reach of the Snoqualmie River in Washington State, and is used to estimate flood extent, depth, arrival time, and peak time for the RFC forecast. The results of the model are processed automatically by a Geographic Information System (GIS) into maps of flood extent, depth, and arrival and peak times. These maps subsequently are processed into formats acceptable by an Internet map server (IMS). The IMS application is a user-friendly interface to access the maps over the Internet; it allows users to select what information they wish to see presented and allows the authors to define scale-dependent availability of map layers and their symbology (appearance of map features). For example, the IMS presents a background of a digital USGS 1:100,000-scale quadrangle at smaller scales, and automatically switches to an ortho-rectified aerial photograph (a digital photograph that has camera angle and tilt distortions removed) at larger scales so viewers can see ground features that help them identify their area of interest more effectively. For the user, the option exists to select either background at any scale. Similar options are provided for both the map creator and the viewer for the various flood maps. This combination of a robust model, emerging IMS software, and application
Burton, Taylour G; Rifai, Hanadi S; Hildenbrand, Zacariah L; Carlton, Doug D; Fontenot, Brian E; Schug, Kevin A
2016-03-01
Hydraulic fracturing operations have been viewed as the cause of certain environmental issues including groundwater contamination. The potential for hydraulic fracturing to induce contaminant pathways in groundwater is not well understood since gas wells are completed while isolating the water table and the gas-bearing reservoirs lay thousands of feet below the water table. Recent studies have attributed ground water contamination to poor well construction and leaks in the wellbore annulus due to ruptured wellbore casings. In this paper, a geospatial model of the Barnett Shale region was created using ArcGIS. The model was used for spatial analysis of groundwater quality data in order to determine if regional variations in groundwater quality, as indicated by various groundwater constituent concentrations, may be associated with the presence of hydraulically fractured gas wells in the region. The Barnett Shale reservoir pressure, completions data, and fracture treatment data were evaluated as predictors of groundwater quality change. Results indicated that elevated concentrations of certain groundwater constituents are likely related to natural gas production in the study area and that beryllium, in this formation, could be used as an indicator variable for evaluating fracturing impacts on regional groundwater quality. Results also indicated that gas well density and formation pressures correlate to change in regional water quality whereas proximity to gas wells, by itself, does not. The results also provided indirect evidence supporting the possibility that micro annular fissures serve as a pathway transporting fluids and chemicals from the fractured wellbore to the overlying groundwater aquifers. PMID:26745299
Implicit level set algorithms for modelling hydraulic fracture propagation.
Peirce, A
2016-10-13
Hydraulic fractures are tensile cracks that propagate in pre-stressed solid media due to the injection of a viscous fluid. Developing numerical schemes to model the propagation of these fractures is particularly challenging due to the degenerate, hypersingular nature of the coupled integro-partial differential equations. These equations typically involve a singular free boundary whose velocity can only be determined by evaluating a distinguished limit. This review paper describes a class of numerical schemes that have been developed to use the multiscale asymptotic behaviour typically encountered near the fracture boundary as multiple physical processes compete to determine the evolution of the fracture. The fundamental concepts of locating the free boundary using the tip asymptotics and imposing the tip asymptotic behaviour in a weak form are illustrated in two quite different formulations of the governing equations. These formulations are the displacement discontinuity boundary integral method and the extended finite-element method. Practical issues are also discussed, including new models for proppant transport able to capture 'tip screen-out'; efficient numerical schemes to solve the coupled nonlinear equations; and fast methods to solve resulting linear systems. Numerical examples are provided to illustrate the performance of the numerical schemes. We conclude the paper with open questions for further research. This article is part of the themed issue 'Energy and the subsurface'. PMID:27597787
Engine Hydraulic Stability. [injector model for analyzing combustion instability
NASA Technical Reports Server (NTRS)
Kesselring, R. C.; Sprouse, K. M.
1977-01-01
An analytical injector model was developed specifically to analyze combustion instability coupling between the injector hydraulics and the combustion process. This digital computer dynamic injector model will, for any imposed chamber of inlet pressure profile with a frequency ranging from 100 to 3000 Hz (minimum) accurately predict/calculate the instantaneous injector flowrates. The injector system is described in terms of which flow segments enter and leave each pressure node. For each flow segment, a resistance, line lengths, and areas are required as inputs (the line lengths and areas are used in determining inertance). For each pressure node, volume and acoustic velocity are required as inputs (volume and acoustic velocity determine capacitance). The geometric criteria for determining inertances of flow segments and capacitance of pressure nodes was set. Also, a technique was developed for analytically determining time averaged steady-state pressure drops and flowrates for every flow segment in an injector when such data is not known. These pressure drops and flowrates are then used in determining the linearized flow resistance for each line segment of flow.
Assessing the impact of different sources of topographic data on 1-D hydraulic modelling of floods
NASA Astrophysics Data System (ADS)
Ali, A. Md; Solomatine, D. P.; Di Baldassarre, G.
2015-01-01
Topographic data, such as digital elevation models (DEMs), are essential input in flood inundation modelling. DEMs can be derived from several sources either through remote sensing techniques (spaceborne or airborne imagery) or from traditional methods (ground survey). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), the light detection and ranging (lidar), and topographic contour maps are some of the most commonly used sources of data for DEMs. These DEMs are characterized by different precision and accuracy. On the one hand, the spatial resolution of low-cost DEMs from satellite imagery, such as ASTER and SRTM, is rather coarse (around 30 to 90 m). On the other hand, the lidar technique is able to produce high-resolution DEMs (at around 1 m), but at a much higher cost. Lastly, contour mapping based on ground survey is time consuming, particularly for higher scales, and may not be possible for some remote areas. The use of these different sources of DEM obviously affects the results of flood inundation models. This paper shows and compares a number of 1-D hydraulic models developed using HEC-RAS as model code and the aforementioned sources of DEM as geometric input. To test model selection, the outcomes of the 1-D models were also compared, in terms of flood water levels, to the results of 2-D models (LISFLOOD-FP). The study was carried out on a reach of the Johor River, in Malaysia. The effect of the different sources of DEMs (and different resolutions) was investigated by considering the performance of the hydraulic models in simulating flood water levels as well as inundation maps. The outcomes of our study show that the use of different DEMs has serious implications to the results of hydraulic models. The outcomes also indicate that the loss of model accuracy due to re-sampling the highest resolution DEM (i.e. lidar 1 m) to lower resolution is much less than the loss of model accuracy due
Yoshimura, K.; Sakashita, S.; Ando, K.; Bruines, P.; Blechschmidt, I.; Kickmaier, W.; Onishi, Y.; Nishiyama, S.
2007-07-01
The objective of this study is to establish a technique to obtain hydraulic conductivity distribution in granite rock masses using seismic tomography. We apply the characteristic that elastic wave velocity disperses in fully saturated porous media on frequency and this velocity dispersion is governed by the hydraulic conductivity - this characteristic has been confirmed in laboratory experiments. The feasibility and design of the field experiment was demonstrated in a first step with numerical simulations. In a second step we applied the technique to the fractured granite at the Grimsel Test Site in Switzerland. The emphasis of the field campaign was on the evaluation of the range of applicability of this technique. The field campaign was structured in three steps, each one corresponding to a larger spatial scale. First, the seismic tomography was applied to a small area - the two boreholes were located at a distance of 1.5 m. In the following step, we selected a larger area, in which the distance of the boreholes amounts to 10 m and the field corresponds to a more complex geology. Finally we applied the testing to a field where the borehole distance was of the order of 75 m. We also drilled a borehole to confirm hydraulic characteristic and reviewed hydraulic model in the 1.5 m cross-hole location area. The results from the field campaign are presented and their application to the various fields are discussed and evaluated. (authors)
NASA Astrophysics Data System (ADS)
Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.
2015-02-01
The rapid progress of Light Detection And Ranging (LiDAR) technology has made acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially with regards to the study of floodplain flow modeling. High-resolution DEM data include many redundant interpolation points, needs a high amount of calculation, and does not match the size of computational mesh. These disadvantages are a common problem for floodplain flow modeling studies. Two-dimensional (2-D) hydraulic modeling, a popular method of analyzing floodplain flow, offers high precision of elevation parameterization for computational mesh while ignoring much micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, enabling the calculation of flood water levels in DEM grid cells through local inverse distance weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of the connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha polder, a flood storage area of Dongting Lake, using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than DEM only simulations.
NASA Astrophysics Data System (ADS)
Rockhold, M. L.; Yarwood, R. R.; Niemet, M. R.; Bottomley, P. J.; Selker, J. S.
Bacterial-induced changes in the hydraulic properties of porous media are important in a variety of disciplines. Most of the previous research on this topic has focused on liquid-saturated porous media systems that are representative of aquifer sediments. Unsaturated or variably saturated systems such as soils require additional considerations that have not been fully addressed in the literature. This paper reviews some of the earlier studies on bacterial-induced changes in the hydraulic properties of saturated porous media, and discusses characteristics of unsaturated or variably saturated porous media that may be important to consider when modeling such phenomena in these systems. New data are presented from experiments conducted in sand-packed columns with initially steady unsaturated flow conditions that show significant biomass-induced changes in pressure heads and water contents and permeability reduction during growth of a Pseudomonas fluorescens bacterium.
Rockhold, Mark L.; Yarwood, R. R.; Niemet, Michael R.; Bottomley, Peter J.; Selker, John S.
2002-07-26
Bacterial-induced changes in the hydraulic properties of porous media are important in a variety of disciplines. Most of the pervious research on this topic has focused on liquid-saturated porous media systems that are representative of aquifer sediments. Unsaturated or variably saturated systems such as soils require additional considerations that have not been fully addressed in the literature. This paper reviews some of the earlier studies on bacterial-induced changes in the hydraulic properties of saturated porous media, and discusses characteristics of unsaturated or variably saturated porous media that may be important to consider when modeling such phenomena in these systems. New data are presented from experiments conducted in sand-packed columns with initially steady unsaturated flow conditions that show significant biomass-induced changes in pressure heads and water contents and permeability reduction during growth of a Pseudomonas fluorescens bacterium.
Assessment of uncertainties of the models used in thermal-hydraulic computer codes
NASA Astrophysics Data System (ADS)
Gricay, A. S.; Migrov, Yu. A.
2015-09-01
The article deals with matters concerned with the problem of determining the statistical characteristics of variable parameters (the variation range and distribution law) in analyzing the uncertainty and sensitivity of calculation results to uncertainty in input data. A comparative analysis of modern approaches to uncertainty in input data is presented. The need to develop an alternative method for estimating the uncertainty of model parameters used in thermal-hydraulic computer codes, in particular, in the closing correlations of the loop thermal hydraulics block, is shown. Such a method shall feature the minimal degree of subjectivism and must be based on objective quantitative assessment criteria. The method includes three sequential stages: selecting experimental data satisfying the specified criteria, identifying the key closing correlation using a sensitivity analysis, and carrying out case calculations followed by statistical processing of the results. By using the method, one can estimate the uncertainty range of a variable parameter and establish its distribution law in the above-mentioned range provided that the experimental information is sufficiently representative. Practical application of the method is demonstrated taking as an example the problem of estimating the uncertainty of a parameter appearing in the model describing transition to post-burnout heat transfer that is used in the thermal-hydraulic computer code KORSAR. The performed study revealed the need to narrow the previously established uncertainty range of this parameter and to replace the uniform distribution law in the above-mentioned range by the Gaussian distribution law. The proposed method can be applied to different thermal-hydraulic computer codes. In some cases, application of the method can make it possible to achieve a smaller degree of conservatism in the expert estimates of uncertainties pertinent to the model parameters used in computer codes.
NASA Astrophysics Data System (ADS)
Savage, James; Pianosi, Francesca; Bates, Paul; Freer, Jim; Wagener, Thorsten
2015-04-01
Predicting flood inundation extents using hydraulic models is subject to a number of critical uncertainties. For a specific event, these uncertainties are known to have a large influence on model outputs and any subsequent analyses made by risk managers. Hydraulic modellers often approach such problems by applying uncertainty analysis techniques such as the Generalised Likelihood Uncertainty Estimation (GLUE) methodology. However, these methods do not allow one to attribute which source of uncertainty has the most influence on the various model outputs that inform flood risk decision making. Another issue facing modellers is the amount of computational resource that is available to spend on modelling flood inundations that are 'fit for purpose' to the modelling objectives. Therefore a balance needs to be struck between computation time, realism and spatial resolution, and effectively characterising the uncertainty spread of predictions (for example from boundary conditions and model parameterisations). However, it is not fully understood how much of an impact each factor has on model performance, for example how much influence changing the spatial resolution of a model has on inundation predictions in comparison to other uncertainties inherent in the modelling process. Furthermore, when resampling fine scale topographic data in the form of a Digital Elevation Model (DEM) to coarser resolutions, there are a number of possible coarser DEMs that can be produced. Deciding which DEM is then chosen to represent the surface elevations in the model could also influence model performance. In this study we model a flood event using the hydraulic model LISFLOOD-FP and apply Sobol' Sensitivity Analysis to estimate which input factor, among the uncertainty in model boundary conditions, uncertain model parameters, the spatial resolution of the DEM and the choice of resampled DEM, have the most influence on a range of model outputs. These outputs include whole domain maximum
A statistical model for seismic hazard assessment of hydraulic-fracturing-induced seismicity
NASA Astrophysics Data System (ADS)
Hajati, T.; Langenbruch, C.; Shapiro, S. A.
2015-12-01
We analyze the interevent time distribution of hydraulic-fracturing-induced seismicity collected during 18 stages at four different regions. We identify a universal statistical process describing the distribution of hydraulic-fracturing-induced events in time. The distribution of waiting times between subsequently occurring events is given by the exponential probability density function of the homogeneous Poisson process. Our findings suggest that hydraulic-fracturing-induced seismicity is directly triggered by the relaxation of stress and pore pressure perturbation initially created by the injection. Therefore, compared to this relaxation, the stress transfer caused by the occurrence of preceding seismic events is mainly insignificant for the seismogenesis of subsequently occurring events. We develop a statistical model to compute the occurrence probability of hydraulic-fracturing-induced seismicity. This model can be used to assess the seismic hazard associated with hydraulic fracturing operations. No aftershock triggering has to be included in the statistical model.
Influence of aquifer geometry on karst hydraulics using different distributive modeling approaches
NASA Astrophysics Data System (ADS)
Oehlmann, Sandra; Geyer, Tobias; Licha, Tobias; Birk, Steffen
2013-04-01
The simulation of flow and transport processes in karst systems is a challenge due to the unknown location of highly conductive karst conduit networks. In this work, the influence of aquifer geometry, particularly the geometry of highly conductive discrete elements, on three-dimensional groundwater flow in a large-scale aquifer system is examined. The area of investigation comprises several springs on the Western Swabian Alb / Germany and has an area of approximately 150 km2. The largest spring therein is the Gallusquelle with an annual average discharge of 0.5 m3/s. Long-term spring hydrographs and hydraulic head measurements, as well as several tracer tests, are available from previous work and are used for model calibration. Four distributive continuum and discrete flow models with different degrees of complexity were set-up employing the finite element simulation software Comsol Multiphysics®. Stationary groundwater flow equations were implemented for single continuum and hybrid modeling. The aquifer geometry was modeled previously with the software Geological Objects Computer Aided Design® (GoCAD®) and transferred to Comsol® software. Simulation results show that not only the location of karst conduits but also their geometry has significant impact on the simulated spring discharge and hydraulic head distribution. A constant conduit radius leads to distorted hydraulic head contour lines and a conduit restrained flow regime close to the spring, while a linearly increasing radius towards the spring leads to evenly distributed contour lines. Models with such an increase in conduit diameters allow the simulation of annual discharge for several springs. This result is in agreement with synthetic karst genesis models, which suggest an increase of conduit diameters towards karst springs because of a positive correlation between flow rates and carbonate solution. The software Comsol Multiphysics®, while rarely used for groundwater flow modeling, was found to meet
Technology Transfer Automated Retrieval System (TEKTRAN)
We examined the hydraulics of concentrated flow using unconfined field experimental data over diverse rangeland landscapes, and developed new empirical prediction models of different rangeland concentrated flow hydraulic parameters, which can be applicable across a wide span of rangeland sites, soil...
Two-Dimensional Coupled Distributed Hydrologic-Hydraulic Model Simulation on Watershed
NASA Astrophysics Data System (ADS)
Cea, Miguel; Rodriguez, Martin
2016-03-01
The objective of this work is to develop a coupled distributed model that enables to analyze water movement in watershed as well as analyze the rainfall-runoff. More specifically, it allows to estimate the various hydrologic water cycle variables at each point of the watershed. In this paper, we have carried out a coupled model of a distributed hydrological and two-dimensional hydraulic models. We have incorporated a hydrological rainfall-runoff model calculated by cell based on the Soil Conservation Service (SCS) method to the hydraulic model, leaving it for the hydraulic model (GUAD2D) to conduct the transmission to downstream cells. The goal of the work is demonstrate the improved predictive capability of the coupled Hydrological-Hydraulic models in a watershed.
Coupled Numerical Study of Turbidity Currents, Internal Hydraulic Jump and Morphological Signatures
NASA Astrophysics Data System (ADS)
Hu, P.; Cao, Z.; He, Z.; Gareth, P.
2013-12-01
Abstract: The last two decades have seen intensive experimental and numerical studies of the occurrence condition of internal hydraulic jump in turbidity currents and the induced morphological signatures (Garcia and Parker 1989; Kostic and Parker 2006). Yet there are two critical issues that remain insufficiently or inappropriately addressed. First, depositional turbidity currents are imposed on steep slopes in both flume experiments and numerical cases, exclusively based on a configuration consisting of an upstream sloping portion and a downstream horizontal portion linked by a slope break. This appears physically counterintuitive as steep slope should favour self-accelerating erosional turbidity currents (Parker et al. 1986). The second issue concerns the numerical studies. There exist significant interactions among the current, sediment transport and bed topography. Due to the slope break in bed, the current may experience an internal hydraulic jump, leaving morphological signatures on the bed, which in turn affects the current evolution. Nevertheless, simplified decoupled models are exclusively employed in previous numerical investigations, in which the interactions are either partly or completely ignored without sufficient justification. The present paper aims to address the above-mentioned two issues relevant to the occurrence condition of the internal hydraulic jump and the induced morphological signatures. A recently developed well-balanced coupled numerical model for turbidity currents (Hu et al. 2012) is applied. In contrast to previous studies, erosional turbidity currents will be imposed at the upstream boundary, which is much more typical of the field. The effects of sediment size, bed slope decrease, and upstream and downstream boundary conditions are revealed in detail. In addition, the evolution of turbidity currents over a bed characterized by gradual decrease in slope is also discussed. References Garcia, M. H., and Parker, G. (1989). Experiments
Assessing the impact of different sources of topographic data on 1-D hydraulic modelling of floods
NASA Astrophysics Data System (ADS)
Ali, A. Md; Solomatine, D. P.; Di Baldassarre, G.
2014-07-01
Topographic data, such as digital elevation models (DEMs), are essential input in flood inundation modelling. DEMs can be derived from several sources either through remote sensing techniques (space-borne or air-borne imagery) or from traditional methods (ground survey). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), the Light Detection and Ranging (LiDAR), and topographic contour maps are some of the most commonly used sources of data for DEMs. These DEMs are characterized by different precision and accuracy. On the one hand, the spatial resolution of low-cost DEMs from satellite imagery, such as ASTER and SRTM, is rather coarse (around 30-90 m). On the other hand, LiDAR technique is able to produce a high resolution DEMs (around 1m), but at a much higher cost. Lastly, contour mapping based on ground survey is time consuming, particularly for higher scales, and may not be possible for some remote areas. The use of these different sources of DEM obviously affects the results of flood inundation models. This paper shows and compares a number of hydraulic models developed using HEC-RAS as model code and the aforementioned sources of DEM as geometric input. The study was carried out on a reach of the Johor River, in Malaysia. The effect of the different sources of DEMs (and different resolutions) was investigated by considering the performance of the hydraulic models in simulating flood water levels as well as inundation maps. The outcomes of our study show that the use of different DEMs has serious implications to the results of hydraulic models. The outcomes also indicates the loss of model accuracy due to re-sampling the highest resolution DEM (i.e. LiDAR 1 m) to lower resolution are much less compared to the loss of model accuracy due to the use of low-cost DEM that have not only a lower resolution, but also a lower quality. Lastly, to better explore the sensitivity of the hydraulic models
Comparison between InfoWorks hydraulic results and a physical model of an urban drainage system.
Rubinato, Matteo; Shucksmith, James; Saul, Adrian J; Shepherd, Will
2013-01-01
Urban drainage systems are frequently analysed using hydraulic modelling software packages such as InfoWorks CS or MIKE-Urban. The use of such modelling tools allows the evaluation of sewer capacity and the likelihood and impact of pluvial flood events. Models can also be used to plan major investments such as increasing storage capacity or the implementation of sustainable urban drainage systems. In spite of their widespread use, when applied to flooding the results of hydraulic models are rarely compared with field or laboratory (i.e. physical modelling) data. This is largely due to the time and expense required to collect reliable empirical data sets. This paper describes a laboratory facility which will enable an urban flood model to be verified and generic approaches to be built. Results are presented from the first phase of testing, which compares the sub-surface hydraulic performance of a physical scale model of a sewer network in Yorkshire, UK, with downscaled results from a calibrated 1D InfoWorks hydraulic model of the site. A variety of real rainfall events measured in the catchment over a period of 15 months (April 2008-June 2009) have been both hydraulically modelled and reproduced in the physical model. In most cases a comparison of flow hydrographs generated in both hydraulic and physical models shows good agreement in terms of velocities which pass through the system. PMID:23863430
Triaxial coreflood study of the hydraulic fracturing of Utica Shale
NASA Astrophysics Data System (ADS)
Carey, J. W.; Frash, L.; Viswanathan, H. S.
2015-12-01
One of the central questions in unconventional oil and gas production research is the cause of limited recovery of hydrocarbon. There are many hypotheses including: 1) inadequate penetration of fractures within the stimulated volume; 2) limited proppant delivery; 3) multiphase flow phenomena that blocks hydrocarbon migration; etc. Underlying any solution to this problem must be an understanding of the hydrologic properties of hydraulically fractured shale. In this study, we conduct triaxial coreflood experiments using a gasket sealing mechanism to characterize hydraulic fracture development and permeability of Utica Shale samples. Our approach also includes fracture propagation with proppants. The triaxial coreflood experiments were conducted with an integrated x-ray tomography system that allows direct observation of fracture development using x-ray video radiography and x-ray computed tomography at elevated pressure. A semi-circular, fracture initiation notch was cut into an end-face of the cylindrical samples (1"-diameter with lengths from 0.375 to 1"). The notch was aligned parallel with the x-ray beam to allow video radiography of fracture growth as a function of injection pressure. The proppants included tungsten powder that provided good x-ray contrast for tracing proppant delivery and distribution within the fracture system. Fractures were propagated at injection pressures in excess of the confining pressure and permeability measurements were made in samples where the fractures propagated through the length of the sample, ideally without penetrating the sample sides. Following fracture development, permeability was characterized as a function of hydrostatic pressure and injection pressure. X-ray video radioadiography was used to study changes in fracture aperture in relation to permeability and proppant embedment. X-ray tomography was collected at steady-state conditions to fully characterize fracture geometry and proppant distribution.
NASA Astrophysics Data System (ADS)
Favrel, A.; Landry, C.; Müller, A.; Avellan, F.
2012-11-01
Resonance in hydraulic systems is characterized by pressure fluctuations of high amplitude which can lead to undesirable and dangerous effects, such as noise, vibration and structural failure. For a Francis turbine operating at partial load, the cavitating vortex rope developing at the outlet of the runner induces pressure fluctuations which can excite the hydraulic system resonance, leading to undesirable large torque and power fluctuations. At resonant operating points, the prediction of amplitude pressure fluctuations by hydro-acoustic models breaks down and gives unreliable results. A more detailed knowledge of the eigenmodes and a better understanding of phenomenon occurring at resonance could allow improving the hydro-acoustic models prediction.This paper presents an experimental identification of a resonance observed in a close-looped hydraulic system with a Francis turbine reduced scale model operating at partial load. The resonance is excited matching one of the test rig eigenfrequencies with the vortex rope precession frequency. At this point, the hydro-acoustic response of the test rig is studied more precisely and used finally to reproduce the shape of the excited eigenmode.
Modeling water-table fluctuations in a sloping aquifer with random hydraulic conductivity
NASA Astrophysics Data System (ADS)
Srivastava, Kirti; Rai, S.; Singh, R.
2002-01-01
To prevent environmental problems like water logging and increase in soil salinity which are responsible for the degradation of the top productive soils, an optimum ditch drainage design is required. For this purpose a knowledge of the spatio-temporal variation of the water table is essential. In this study the spatio-temporal variation of the water table in a sloping ditch drainage system has been modeled from a stochastic point of view, incorporating randomness in hydraulic conductivity to get the expression for the mean and the standard deviation of the water-table height. The hydraulic conductivity has been considered to be a realization of a log-normal distribution. Application of these expressions in the prediction of mean water-table variation with the associated error bounds has been demonstrated with the help of a ditch drainage problem of a sloping aquifer. The sensitivity analysis has also been carried out to see the effect of variability in the hydraulic conductivity on the water-table fluctuations. The error bounds quantified on the water-table height will thus help in the decision-making process for proper drainage design.
NASA Astrophysics Data System (ADS)
Neuville, Amélie; Toussaint, Renaud; Schmittbuhl, Jean
2011-09-01
Natural open joints in rocks commonly present multiscale self-affine apertures. This geometrical complexity affects fluid transport and heat exchange between the flowing fluid and the surrounding rock. In particular, long range correlations of self-affine apertures induce strong channelling of the flow which influences both mass and heat advection. A key question is to find a geometrical model of the complex aperture that describes at best the macroscopic properties (hydraulic conductivity, heat exchange) with the smallest number of parameters. Solving numerically the Stokes and heat equations with a lubrication approximation, we show that a low pass filtering of the aperture geometry provides efficient estimates of the effective hydraulic and thermal properties (apertures). A detailed study of the influence of the bandwidth of the lowpass filtering on these transport properties is also performed. For instance, keeping the information of amplitude only of the largest Fourier length scales allows us to reach already an accuracy of 9 per cent on the hydraulic and the thermal apertures.
NASA Astrophysics Data System (ADS)
Liu, Z.; Rajib, M. A.; Jafarzadegan, K.; Merwade, V.
2015-12-01
Application of land surface/hydrologic models within an operational flood forecasting system can provide probable time of occurrence and magnitude of streamflow at specific locations along a stream. Creating time-varying spatial extent of flood inundation and depth requires the use of a hydraulic or hydrodynamic model. Models differ in representing river geometry and surface roughness which can lead to different output depending on the particular model being used. The result from a single hydraulic model provides just one possible realization of the flood extent without capturing the uncertainty associated with the input or the model parameters. The objective of this study is to compare multiple hydraulic models toward generating ensemble flood inundation extents. Specifically, relative performances of four hydraulic models, including AutoRoute, HEC-RAS, HEC-RAS 2D, and LISFLOOD are evaluated under different geophysical conditions in several locations across the United States. By using streamflow output from the same hydrologic model (SWAT in this case), hydraulic simulations are conducted for three configurations: (i) hindcasting mode by using past observed weather data at daily time scale in which models are being calibrated against USGS streamflow observations, (ii) validation mode using near real-time weather data at sub-daily time scale, and (iii) design mode with extreme streamflow data having specific return periods. Model generated inundation maps for observed flood events both from hindcasting and validation modes are compared with remotely sensed images, whereas the design mode outcomes are compared with corresponding FEMA generated flood hazard maps. The comparisons presented here will give insights on probable model-specific nature of biases and their relative advantages/disadvantages as components of an operational flood forecasting system.
NASA Astrophysics Data System (ADS)
Laine-Kaulio, H.; Karvonen, T.; Koivusalo, H.; Lauren, A.; Saastamoinen, S.
2009-04-01
Shallow till layers typically overlay bedrock in forested areas in the boreal region. In forested tills, preferential flowpaths related to the soil structure have a decisive influence on hydrogeological properties such as the soil hydraulic conductivity. Hydraulic conductivity is also proven to depend on the observation scale. Traditional soil core samples cannot capture the impact of soil structure on hillslope scale conductivities. Measurements and observations made at different scales, combined with simulation models, are essential for investigating conductivity properties and flow and transport processes in forest soils. This study combined a set of soil analyses and field experiments with physics-based modelling to investigate the hydraulic properties of a forested till slope in Finland. The main objective was to i) determine the saturated hydraulic conductivity in the study slope with methods related to different scales, and to ii) study the utilisation of the conductivity results in modelling flow and solute transport in the slope. Soil sampling, dye, and ion tracer experiments were conducted in a forested hillslope in Eastern Finland. In the 20 m long study section of the slope the mean slope was about 15 %. The haplic podsol profile above bedrock had a thickness of 0.8 m and was formed of sandy till. The soil was very stony and heterogeneous in terms of granularity and pore size distribution. Granularity, porosity and proportion of macropores reduced clearly with depth. Dye tracer experiments revealed three types of preferential flow routes in the slope: i) stone surfaces, ii) areas of coarse-grained soil material, and iii) decayed root channels. Both living roots and preferential flowpaths reached the transitional zone of the podsol at about 0.5 m depth, but living roots were not found to function unequivocally as preferential flowpaths. The saturated hydraulic conductivity was determined using three methods: i) from soil core samples in laboratory, ii
Thermal hydraulic characteristics study of prototype NET and CEA cable-in-conduit conductors (CICCs)
Maekawa, Ryuji
1995-10-31
The thermal hydraulic characteristics of low temperature helium in a Cable-in-Conduit Conductor (CICC) significantly affects the overall design and performance of the associated large scale superconducting magnet system. It is essential to understand the transient and steady state behavior of the helium in the conductor. Throughout the development of CICCs, the reduction of flow impedance has been one of the key factors to improving the overall pressure drop. The newly developed CICC for the ITER project has a hybrid cooling scheme: a central channel that is surrounded by bundles, for which the thermal hydraulic characteristics are not well understood. This thesis describes an experimental and analytical investigation of thermal hydraulic characteristics of low temperature helium in conventional and hybrid CICCS. Pressure drop measurements for both NET and CEA conductors have been conducted, using low temperature helium and liquid nitrogen to obtain a range of Reynolds numbers. The results are correlated with classical friction factor and Reynolds number analysis. The flow impedance reduction of the CEA conductor is described by measures of a developed flow model. Thermally induced flow in the CEA conductor has been studied with an inductive heating method. The induced velocity in the central channel is measured by a Pitot tube with steady state Reynolds number up to {approximately}7000. The transient pressure wave propagation has been recorded with pressure transducers placed equally along the conductor. The supercritical helium temperature in the central channel has been measured with the thermometer probe. However, the reduction of the central channel area significantly affects the overall thermal hydraulic characteristics of the conductor. The results suggest the importance of the central channel. A transient heat transfer experiment studied the.transverse heat transfer mechanism in the CEA conductor. The temperatures in the central channel and bundle region
NASA Astrophysics Data System (ADS)
Driba, D. L.; De Lucia, M.; Peiffer, S.
2014-12-01
Fluid-rock interactions in geothermal reservoirs are driven by the state of disequilibrium that persists among solid and solutes due to changing temperature and pressure. During operation of enhanced geothermal systems, injection of cooled water back into the reservoir disturbs the initial thermodynamic equilibrium between the reservoir and its geothermal fluid, which may induce modifications in permeability through changes in porosity and pore space geometry, consequently bringing about several impairments to the overall system.Modeling of fluid-rock interactions induced by injection of cold brine into Groß Schönebeck geothermal reservoir system situated in the Rotliegend sandstone at 4200m depth have been done by coupling geochemical modeling Code Phreeqc with OpenGeoSys. Through batch modeling the re-evaluation of the measured hydrochemical composition of the brine has been done using Quintessa databases, the results from the calculation indicate that a mineral phases comprising of K-feldspar, hematite, Barite, Calcite and Dolomite was found to match the hypothesis of equilibrium with the formation fluid, Reducing conditions are presumed in the model (pe = -3.5) in order to match the amount of observed dissolved Fe and thus considered as initial state for the reactive transport modeling. based on a measured composition of formation fluids and the predominant mineralogical assemblage of the host rock, a preliminary 1D Reactive transport modeling (RTM) was run with total time set to 30 years; results obtained for the initial simulation revealed that during this period, no significant change is evident for K-feldspar. Furthermore, the precipitation of calcite along the flow path in the brine results in a drop of pH from 6.2 to a value of 5.2 noticed over the simulated period. The circulation of cooled fluid in the reservoir is predicted to affect the temperature of the reservoir within the first 100 -150m from the injection well. Examination of porosity change in
Rutqvist, Jonny; Rinaldi, Antonio P.; Cappa, Frédéric; Moridis, George J.
2013-07-01
We have conducted numerical simulation studies to assess the potential for injection-induced fault reactivation and notable seismic events associated with shale-gas hydraulic fracturing operations. The modeling is generally tuned towards conditions usually encountered in the Marcellus shale play in the Northeastern US at an approximate depth of 1500 m (~;;4,500 feet). Our modeling simulations indicate that when faults are present, micro-seismic events are possible, the magnitude of which is somewhat larger than the one associated with micro-seismic events originating from regular hydraulic fracturing because of the larger surface area that is available for rupture. The results of our simulations indicated fault rupture lengths of about 10 to 20 m, which, in rare cases can extend to over 100 m, depending on the fault permeability, the in situ stress field, and the fault strength properties. In addition to a single event rupture length of 10 to 20 m, repeated events and aseismic slip amounted to a total rupture length of 50 m, along with a shear offset displacement of less than 0.01 m. This indicates that the possibility of hydraulically induced fractures at great depth (thousands of meters) causing activation of faults and creation of a new flow path that can reach shallow groundwater resources (or even the surface) is remote. The expected low permeability of faults in producible shale is clearly a limiting factor for the possible rupture length and seismic magnitude. In fact, for a fault that is initially nearly-impermeable, the only possibility of larger fault slip event would be opening by hydraulic fracturing; this would allow pressure to penetrate the matrix along the fault and to reduce the frictional strength over a sufficiently large fault surface patch. However, our simulation results show that if the fault is initially impermeable, hydraulic fracturing along the fault results in numerous small micro-seismic events along with the propagation, effectively
Growth model for large branched three-dimensional hydraulic crack system in gas or oil shale.
Chau, Viet T; Bažant, Zdeněk P; Su, Yewang
2016-10-13
Recent analysis of gas outflow histories at wellheads shows that the hydraulic crack spacing must be of the order of 0.1 m (rather than 1 m or 10 m). Consequently, the existing models, limited to one or several cracks, are unrealistic. The reality is 10(5)-10(6) almost vertical hydraulic cracks per fracking stage. Here, we study the growth of two intersecting near-orthogonal systems of parallel hydraulic cracks spaced at 0.1 m, preferably following pre-existing rock joints. One key idea is that, to model lateral cracks branching from a primary crack wall, crack pressurization, by viscous Poiseuille-type flow, of compressible (proppant-laden) frac water must be complemented with the pressurization of a sufficient volume of micropores and microcracks by Darcy-type water diffusion into the shale, to generate tension along existing crack walls, overcoming the strength limit of the cohesive-crack or crack-band model. A second key idea is that enforcing the equilibrium of stresses in cracks, pores and water, with the generation of tension in the solid phase, requires a new three-phase medium concept, which is transitional between Biot's two-phase medium and Terzaghi's effective stress and introduces the loading of the solid by pressure gradients of diffusing pore water. A computer program, combining finite elements for deformation and fracture with volume elements for water flow, is developed to validate the new model.This article is part of the themed issue 'Energy and the subsurface'. PMID:27597791
75 FR 36387 - Informational Public Meetings for Hydraulic Fracturing Research Study; Correction
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-25
..., 2010, in FR doc. 2010-14897, on page 35023, in the third Column, correct the Web site addresses shown... AGENCY Informational Public Meetings for Hydraulic Fracturing Research Study; Correction AGENCY... Hydraulic Fracturing Research Study. The document contained an incorrect EPA Web site address in two...
Final Plan to Study the Potential Impacts of Hydraulic Fracturing on Drinking Water Resources
The overall purpose of this study is to elucidate the relationship, if any, between hydraulic fracturing and drinking water resources. More specifically, the study has been designed to assess the potential impacts of hydraulic fracturing on drinking water resources and to identif...
A new empirical model for estimating the hydraulic conductivity of low permeability media
NASA Astrophysics Data System (ADS)
Qi, S.; Wen, Z.; Lu, C.; Shu, L.; Shao, J.; Huang, Y.; Zhang, S.; Huang, Y.
2015-05-01
Hydraulic conductivity (K) is one of the significant soil characteristics in terms of flow movement and solute transport. It has been recognized that K is statistically related to the grain-size distribution. Numerous models have been developed to reveal the relationship between K and the grain-size distribution of soil, but most of these are inappropriate for fine-grained media. Therefore, a new empirical model for estimating K of low permeability media was proposed in this study. In total, the values of K of 30 soil samples collected in the Jiangning District of Nanjing were measured using the single-ring infiltrometer method. The new model was developed using the percentages of sand, silt and clay-sized particles, and the first and the second rank moment of the grain-size through the moment method as predictor variables. Multivariate nonlinear regression analysis yielded a coefficient of determination (R2) of 0.75, indicating that this empirical model seems to provide a new approach for the indirect determination of hydraulic conductivity of low permeability media.
Experimental Study on Abrasive Waterjet Polishing of Hydraulic Turbine Blades
NASA Astrophysics Data System (ADS)
Khakpour, H.; Birglenl, L.; Tahan, A.; Paquet, F.
2014-03-01
In this paper, an experimental investigation is implemented on the abrasive waterjet polishing technique to evaluate its capability in polishing of surfaces and edges of hydraulic turbine blades. For this, the properties of this method are studied and the main parameters affecting its performance are determined. Then, an experimental test-rig is designed, manufactured and tested to be used in this study. This test-rig can be used to polish linear and planar areas on the surface of the desired workpieces. Considering the number of parameters and their levels, the Taguchi method is used to design the preliminary experiments. All experiments are then implemented according to the Taguchi L18 orthogonal array. The signal-to-noise ratios obtained from the results of these experiments are used to determine the importance of the controlled polishing parameters on the final quality of the polished surface. The evaluations on these ratios reveal that the nozzle angle and the nozzle diameter have the most important impact on the results. The outcomes of these experiments can be used as a basis to design a more precise set of experiments in which the optimal values of each parameter can be estimated.
Numerical simulation of the two-phase flows in a hydraulic coupling by solving VOF model
NASA Astrophysics Data System (ADS)
Luo, Y.; Zuo, Z. G.; Liu, S. H.; Fan, H. G.; Zhuge, W. L.
2013-12-01
The flow in a partially filled hydraulic coupling is essentially a gas-liquid two-phase flow, in which the distribution of two phases has significant influence on its characteristics. The interfaces between the air and the liquid, and the circulating flows inside the hydraulic coupling can be simulated by solving the VOF two-phase model. In this paper, PISO algorithm and RNG k-ɛ turbulence model were employed to simulate the phase distribution and the flow field in a hydraulic coupling with 80% liquid fill. The results indicate that the flow forms a circulating movement on the torus section with decreasing speed ratio. In the pump impeller, the air phase mostly accumulates on the suction side of the blades, while liquid on the pressure side; in turbine runner, air locates in the middle of the flow passage. Flow separations appear near the blades and the enclosing boundaries of the hydraulic coupling.
Control method and system for hydraulic machines employing a dynamic joint motion model
Danko, George
2011-11-22
A control method and system for controlling a hydraulically actuated mechanical arm to perform a task, the mechanical arm optionally being a hydraulically actuated excavator arm. The method can include determining a dynamic model of the motion of the hydraulic arm for each hydraulic arm link by relating the input signal vector for each respective link to the output signal vector for the same link. Also the method can include determining an error signal for each link as the weighted sum of the differences between a measured position and a reference position and between the time derivatives of the measured position and the time derivatives of the reference position for each respective link. The weights used in the determination of the error signal can be determined from the constant coefficients of the dynamic model. The error signal can be applied in a closed negative feedback control loop to diminish or eliminate the error signal for each respective link.
Mathematical modeling of hydraulic fracturing in coal seams
Olovyanny, A.G.
2005-02-01
Hydraulic fracturing of coal seam is considered as a process of development of discontinuities in rock mass elements due to change in hydrogeomechanical situation on filtration of fluid under pressure. Failure is associated with excess of the effective stresses over the rock tension strength. The problem on filtration and failure of massif is solved by the finite-element method using the procedure of fictitious nodal forces.
A summary of EPA's research relating to potential impacts of hydraulic fracturing on drinking water resources will be presented. Background about the study plan development will be presented along with an analysis of the water cycle as it relates to hydraulic fracturing processe...
Modeling and simulation of hydraulic vibration system based on bond graph and Matlab/Simulink
NASA Astrophysics Data System (ADS)
Lian, Hongzhen; Kou, Ziming
2008-10-01
The hydraulic vibration system controlled by wave exciter is a mechanic-electric-fluid integration system, and it has high dynamic characteristics. Modeling and simulation for it has come to professional's attention in the field of hydraulic vibration industry, because it is nonlinear and complex. In this paper, a method has been proposed. By using power bond graph method, the bond graph model for it can be established, meanwhile, it is proposed that controlled parameters are considered to join the model, in order to control power flow alternated; and the mathematical model(state equations) of this system can be built according to bond graph theory and controlled relations; then simulation model can be built by using Matlab/Simulink software, the model can intuitively express system's power flow direction and controlled relations. To the question that stiff equation appears easily in model of hydraulic system, we can choose the adapting algorithm offered by Matlab software to obtain the more precise simulation results.
Experimental studies of rock fracture behavior related to hydraulic fracture
NASA Astrophysics Data System (ADS)
Ma, Zifeng
The objective of this experimental investigation stems from the uncontrollable of the hydraulic fracture shape in the oil and gas production field. A small-scale laboratory investigation of crack propagation in sandstone was first performed with the objective to simulate the field fracture growth. Test results showed that the fracture resistance increased with crack extension, assuming that there was an interaction between crack faces (bridging, interlocking, and friction). An acoustic emission test was conducted to examine the existence of the interaction by locating AE events and analyzing waveform. Furthermore, the effects of confining stress, loading rate, stress field, and strength heterogeneous on the tortuosity of the fracture surface were experimentally investigated in the study. Finally, a test was designed and conducted to investigate the crack propagation in a stratified media with permeability contrast. Crack was observed to arrested in an interface. The phenomenon of delamination along an interface between layers with permeability contrast was observed. The delamination was proposed to be the cause of crack arrest and crack jump in the saturated stratified materials under confinement test.
NASA Astrophysics Data System (ADS)
Mahmud, K.; Mariethoz, G.; Baker, A.; Sharma, A.
2015-01-01
Hydraulic conductivity is one of the most critical and at the same time one of the most uncertain parameters in many groundwater models. One problem commonly faced is that the data are usually not collected at the same scale as the discretized elements used in a numerical model. Moreover, it is common that different types of hydraulic conductivity measurements, corresponding to different spatial scales, coexist in a studied domain, which have to be integrated simultaneously. Here we address this issue in the context of Image Quilting, one of the recently developed multiple-point geostatistics methods. Based on a training image that represents fine-scale spatial variability, we use the simplified renormalization upscaling method to obtain a series of upscaled training images that correspond to the different scales at which measurements are available. We then apply Image Quilting with such a multiscale training image to be able to incorporate simultaneously conditioning data at several spatial scales of heterogeneity. The realizations obtained satisfy the conditioning data exactly across all scales, but it can come at the expense of a small approximation in the representation of the physical scale relationships. In order to mitigate this approximation, we iteratively apply a kriging-based correction to the finest scale that ensures local conditioning at the coarsest scales. The method is tested on a series of synthetic examples where it gives good results and shows potential for the integration of different measurement methods in real-case hydrogeological models.
Testing sensitivity of the LISFLOOD subgrid hydraulic model to SAR image derived information
NASA Astrophysics Data System (ADS)
Wood, Melissa; Bates, Paul; Neal, Jeff; Hostache, Renaud; Matgen, Patrick; Chini, Marco; Giustarini, Laura
2013-04-01
There has been much interest in the use of Synthetic Aperture Radar (SAR) images to indirectly estimate flood extent and flood elevation to aid the understanding of fluvial flood inundation processes. SAR remote sensing satellites are capable of all-weather day/night observations that can discriminate between land and smooth open water surfaces over large scales. By combining SAR derived information with 2D hydraulic models and terrain data, the mechanisms of flooding can be better simulated therefore enabling more accurate and reliable flood forecasting. The objective of this study is to test the sensitivity of a LISFLOOD subgrid 2D model to its main parameters (i.e. roughness coefficient, river bathymetry) using SAR derived flood extent maps. Because of SAR imaging techniques and processing steps used to derive the flood information, any SAR-derived flood extent image will contain inherent uncertainty. We therefore use the uncertainty of the SAR information to obtain a range of plausible parameters to test sensitivity of the hydraulic model. LISFLOOD is a distributed 2D model developed at the University of Bristol and designed for use with larger ungauged river catchments. The version used employs a subgrid procedure which allows any size of river channel below that of the grid resolution to be represented. This procedure has been shown to improve hydraulic connectivity within the modelled flooded area and thus improve flood prediction for data sparse areas. A hydrodynamic LISFLOOD subgrid model of the River Severn at Tewkesbury covering a domain area of 50x70km and including the confluence with a major tributary (the River Avon) will be utilised. A complete storm hydrograph will be used as inflow to the model to simulate the full flood event. Surveyed cross section and gauged daily flows are also available for the River Severn. Therefore, the model results using variable parameters can be compared against results obtained from ground observations to further
NASA Astrophysics Data System (ADS)
Perkone, E.; Delina, A.; Saks, T.; Raga, B.; Jātnieks, J.; Klints, I.; Popovs, K.; Babre, A.; Bikše, J.; Kalvāns, A.; Retike, I.; Ukass, J.
2012-04-01
very important to take into account the fact that groundwater flow in carbonate aquifers is often almost entirely dependent on jointing and concomitant joint enlargement by dissolution. In this study pumping test results provide a wide range of hydraulic conductivity values, for example in Pļaviņu aquifer hydraulic conductivity varies from 0,03 - 266 m/day but in Daugava aquifer values range from 0,06 - 735 m/day. Pumping test results is provided by Latvian Environment, Geology and Meteorology Centre. Studying average values of hydraulic conductivity there exists a correlation between K and aquifer flat depth - Daugava aquifer, which in geological structure, is located above the Pļaviņu aquifer has higher average K value - 32 m/day, in Pļaviņu aquifer - 27 m/day. Correlative study of the depth and hydraulic conductivity allowed to characterize the mean values as function of the aquifer depth for the regional groundwater flow modelling. This study is supported by the European Social Fund project No. 2009/0212/1DP/1.1.1.2.0/09/APIA/VIAA/060
Analytical solutions in a hydraulic model of seepage with sharp interfaces
NASA Astrophysics Data System (ADS)
Kacimov, A. R.
2002-02-01
Flows in horizontal homogeneous porous layers are studied in terms of a hydraulic model with an abrupt interface between two incompressible Darcian fluids of contrasting density driven by an imposed gradient along the layer. The flow of one fluid moving above a resting finger-type pool of another is studied. A straight interface between two moving fluids is shown to slump, rotate and propagate deeper under periodic drive conditions than in a constant-rate regime. Superpropagation of the interface is related to Philip's superelevation in tidal dynamics and acceleration of the front in vertical infiltration in terms of the Green-Ampt model with an oscillating ponding water level. All solutions studied are based on reduction of the governing PDE to nonlinear ODEs and further analytical and numerical integration by computer algebra routines.
NASA Astrophysics Data System (ADS)
Fabris, L.; Malcolm, I.; Millidine, K. J.; Buddendorf, B.; Tetzlaff, D.; Soulsby, C.
2015-12-01
Wild Atlantic salmon populations in Scottish rivers constitute an important economic and recreational resource, as well as being a key component of biodiversity. Salmon have very specific habitat requirements at different life stages and their distribution is therefore strongly influenced by a complex suite of biological and physical controls. Previous research has shown that stream hydrodynamics and channel morphology have a strong influence on the distribution and density of juvenile salmon. Here, we utilise a unique 20 year data set of spatially distributed juvenile salmon densities derived from annual electro-fishing surveys in an upland Scottish river. We examine to what extent the spatial and temporal variability of in-stream hydraulics regulates the spatial and temporal variability in the performance and density of juvenile salmon. A 2-D hydraulic model (River2D) is used to simulate water velocity and water depth under different flow conditions for seven different electro-fishing sites. The selected sites represent different hydromorphological environments including plane-bed, step-pool and pool riffle reaches. The bathymetry of each site was characterised using a total station providing an accurate DTM of the bed, and hydraulic simulations were driven by 20 year stream flow records. Habitat suitability curves, based on direct observations during electro-fishing surveys, were produced for a range of hydraulic indices for juvenile salmon. The hydraulic simulations showed marked spatial differences in juvenile habitat quality both within and between reaches. They also showed marked differences both within and between years. This is most evident in extreme years with wet summers when salmon feeding opportunities may be constrained. Integration of hydraulic habitat models, with fish preference curves and the long term hydrological data allows us to assess whether long-term changes in hydroclimate may be affecting juvenile salmonid populations in the study stream
NASA Astrophysics Data System (ADS)
Yao, Yao
2012-05-01
Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a Perkins-Kern-Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone
The Effect of Soil Hydraulic Properties vs. Soil Texture in Land Surface Models
NASA Technical Reports Server (NTRS)
Gutmann, E. D.; Small, E. E.
2005-01-01
This study focuses on the effect of Soil Hydraulic Property (SHP) selection on modeled surface fluxes following a rain storm in a semi-arid environment. SHPs are often defined based on a Soil Texture Class (STC). To examine the effectiveness of this approach, the Noah land surface model was run with each of 1306 soils in a large SHP database. Within most STCs, the outputs have a range of 350 W/m2 for latent and sensible heat fluxes, and 8K for surface temperature. The average difference between STC median values is only 100 W/m2 for latent and sensible heat. It is concluded that STC explains 5-15% of the variance in model outputs and should not be used to determine SHPs.
NASA Astrophysics Data System (ADS)
Amano, Ayako; Sakuma, Taisuke; Kazama, So
This study evaluated waterborne infectious diseases risk and incidence rate around Phonm Penh in Cambodia. We use the hydraulic flood simulation, coliform bacterium diffusion model, dose-response model and outpatient data for quantitative analysis. The results obtained are as follows; 1. The incidence (incidence rate) of diarrhea as water borne diseases risk is 0.14 million people (9%) in the inundation area. 2. The residents in the inundation area are exposed up to 4 times as high risk as daily mean calculated by the integrated model combined in the regional scale. 3.The infectious disease risk due to floods and inundation indicated is effective as an element to explain the risk. The scenario explains 34% number of patient estimated by the outpatient data.
Moon, H.Y. ); Advani, S.H.; Lee, T.S. )
1992-11-01
Hydraulic fracturing plays a pivotal role in the enhancement of oil and gas production recovery from low permeability reservoirs. The process of hydraulic fracturing entails the generation of a fracture by pumping fluids blended with special chemicals and proppants into the payzone at high injection rates and pressures to extend and wedge fractures. The mathematical modeling of hydraulically induced fractures generally incorporates coupling between the formation elasticity, fracture fluid flow, and fracture mechanics equations governing the formation structural responses, fluid pressure profile, and fracture growth. Two allied unsymmetric elliptic fracture models are developed for fracture configuration evolutions in three-layered rock formations. The first approach is based on a Lagrangian formulation incorporating pertinent energy components associated with the formation structural responses and fracture fluid flow. The second model is based on a generalized variational principle, introducing an energy rate related functional. These models initially simulate a penny-shaped fracture, which becomes elliptic if the crack tips encounters (upper and/or lower) barriers with differential reservoir properties (in situ stresses, 16 elastic moduli, and fracture toughness-contrasts and fluid leak-off characteristics). The energy rate component magnitudes are determined to interpret the governing hydraulic fracture mechanisms during fracture evolution. The variational principle is extended to study the phenomenon and consequences of fluid lag in fractures. Finally, parametric sensitivity and energy rate investigations to evaluate the roles of controllable hydraulic treatment variables and uncontrollable reservoir property characterization parameters are performed. The presented field applications demonstrate the overall capabilities of the developed models. These studies provide stimulation treatment guidelines for fracture configuration design, control, and optimization.
NASA Astrophysics Data System (ADS)
Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.
2015-08-01
The rapid progress of lidar technology has made the acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially in regards to the study of floodplain flow. However, high-resolution DEM data pose several disadvantages for floodplain modeling studies; e.g., the data sets contain many redundant interpolation points, large numbers of calculations are required to work with data, and the data do not match the size of the computational mesh. Two-dimensional (2-D) hydraulic modeling, which is a popular method for analyzing floodplain flow, offers highly precise elevation parameterization for computational mesh while ignoring much of the micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, thus enabling the calculation of flood water levels in DEM grid cells through local inverse distance-weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha detention basin, which is a flood storage area of Dongting Lake in China, by using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than simulations only with DEM.
NASA Astrophysics Data System (ADS)
Wang, W. L.; Yu, D. S.; Zhou, Z.
2015-10-01
Due to the high-speed operation of modern rail vehicles and severe in-service environment of their hydraulic dampers, it has become important to establish more practical and accurate damper models and apply those models in high-speed transit problem studies. An improved full parametric model with actual in-service parameters, such as variable viscous damping, comprehensive stiffness and small mounting clearance was established for a rail vehicle's axle-box hydraulic damper. A subtle variable oil property model was built and coupled to the modelling process, which included modelling of the dynamic flow losses and the relief-valve system dynamics. The experiments validated the accuracy and robustness of the established full in-service parametric model and simulation which captured the damping characteristics over an extremely wide range of excitation speeds. Further simulations were performed using the model to uncover the effects of key in-service parameter variations on the nominal damping characteristics of the damper. The obtained in-service parametric model coupled all of the main factors that had significant impacts on the damping characteristics, so that the model could be useful in more extensive parameter effects analysis, optimal specification and product design optimisation of hydraulic dampers for track-friendliness, ride comfort and other high-speed transit problems.
Numerical modeling of the near-field hydraulics of water wells.
Houben, Georg J; Hauschild, Sarah
2011-01-01
Numerical flow models can be a useful tool for dimensioning water wells and to investigate the hydraulics in their near-field. Fully laminar flow can be assumed for all models calculated up to the screen. Therefore models can be used to predict--at least qualitatively, neglecting turbulent losses inside the well--the spatial distribution of inflow into the well and the overall hydraulic performance of different combinations of aquifer parameters and technical installations. Models for both horizontal (plan view) and vertical flow (cross section) to wells were calculated for a variety of setups. For the latter, this included variations of hydraulic conductivity of the screen, pump position, and aquifer heterogeneity. Models of suction flow control devices showed that they indeed can homogenize inflow, albeit at the cost of elevated entrance losses. PMID:20860691
NASA Astrophysics Data System (ADS)
Owen, Gareth; Wilkinson, Mark; Nicholson, Alex; Quinn, Paul; O'Donnell, Greg
2015-04-01
river stem and principal tributaries, it is possible to understand in detail how floods develop and propagate, both temporally and spatially. Traditional rainfall-runoff modelling involves the calibration of model parameters to achieve a best fit against an observed flow series, typically at a single location. The modelling approach adopted here is novel in that it directly uses the nested observed information to disaggregate the outlet hydrograph in terms of the source locations. Using a combination of local evidence and expert opinion, the model can be used to assess the impacts of distributed land use management changes and NFM on floods. These studies demonstrate the power of networks of observational instrumentation for constraining hydraulic and hydrologic models for use in prediction.
Hirano, Masashi
1997-07-01
This paper describes the results of a scoping study on seismically induced resonance of nuclear-coupled thermal-hydraulic instability in BWRs, which was conducted by using TRAC-BF1 within a framework of a point kinetics model. As a result of the analysis, it is shown that a reactivity insertion could occur accompanied by in-surge of coolant into the core resulted from the excitation of the nuclear-coupled instability by the external acceleration. In order to analyze this phenomenon more in detail, it is necessary to couple a thermal-hydraulic code with a three-dimensional nuclear kinetics code.
Physical Hydraulic Model of Side-Channel Spillway of Lambuk DAM, Bali
NASA Astrophysics Data System (ADS)
Harifa, A. C.; Sholichin, M.; Othman, F. B.
2013-12-01
The spillway is among the most important structures of a dam project. A spillway is designed to prevent overtopping of a dam at a place that is not designed for overtopping. Side-channel spillways are commonly used to release water flow from a reservoir in places where the sides are steep and have a considerable height above the dam. Experimental results were collected with a hydraulic model of the side-channel spillway for releasing the peak overflow of Lambuk Dam. This dam is, located on the Lambuk River, which is a tributary of the Yeh Hoo River ~ 34.6 km north of Denpasar on the island of Bali. The bituminous geomembrane faced dam is 24 m in height, with a 35-m wide spillway. The length of the side channel is 35 m long, with 58 m of transition channel, 67.37 m of chuteway channel and 22.71 m of stilling basin. The capacity of the spillway is 231.91 m3/s and the outlet works capacity is 165.28 m3/s. The reservoir is designed for irrigation and water supply. The purpose of this study was to optimize the designed of the structure and to ensure its safe operation. In hydraulic model may help the decision-makers to visualize the flow field before selecting a ';suitable' design. The hydraulic model study was performed to ensure passage of the maximum discharge at maximum reservoir capacity; to study the spillway approach conditions, water surface profiles, and flow patterns in the chuteway; and to reveal potential demerits of the proposed hydraulic design of various structures and explore solutions. The model was constructed at 1 : 40 scale, Reservoir topography was modeled using concrete, the river bed using sand and some gravel, the river berm using concrete, and the spillway and channel using Plexiglas. Water was measured using Rectangular contracted weir. Design floods (with return period in year) were Q2 = 111.40 m3/s, Q5 = 136.84 m3/s, Q10 = 159.32 m3/s, Q25 = 174.61 m3/s, Q50 = 185.13 m3/s, Q100 = 198.08 m3/s, Q200 = 210.55 m3/s, Q1000 = 231.91 m3/s and the
Copper River hydraulic study at Million Dollar Bridge, Alaska
Jones, Stanley H.; Barber, William F.
1980-01-01
The Copper River hydraulic conditions in the vicinity of the Million-Dollar Bridge, Alaska, at the outlet to Miles Lake are described. The water discharge, lake and river bed profiles, bathymetry, velocity, and direction of flow are presented. (Kosco-USGS)
NASA Astrophysics Data System (ADS)
Vogel, Tomas; Votrubova, Jana; Dusek, Jaromir; Dohnal, Michal
2016-02-01
In the context of soil water flow modeling, root water uptake is often evaluated based on water potential difference between the soil and the plant (the water potential gradient approach). Root water uptake rate is modulated by hydraulic resistance of both the root itself, and the soil in the root vicinity. The soil hydraulic resistance is a function of actual soil water content and can be assessed assuming radial axisymmetric water flow toward a single root (at the mesoscopic scale). In the present study, three approximate solutions of mesoscopic root water uptake - finite difference approximation, steady-state solution, and steady-rate solution - are examined regarding their ability to capture the pressure head variations in the root vicinity. Insignificance of their differences when implemented in the macroscopic soil water flow model is demonstrated using the critical root water uptake concept. Subsequently, macroscopic simulations of coupled soil water flow and root water uptake are presented for a forest site under temperate humid climate. Predicted soil water pressure heads and actual transpiration rates are compared with observed data. Scenario simulations illustrate uncertainties associated with estimates of root geometrical and hydraulic properties. Regarding the actual transpiration prediction, the correct characterization of active root system geometry and hydraulic properties seems far more important than the choice of a particular mesoscopic model.
NASA Astrophysics Data System (ADS)
Illman, Walter A.; Berg, Steven J.; Zhao, Zhanfeng
2015-05-01
The robust performance of hydraulic tomography (HT) based on geostatistics has been demonstrated through numerous synthetic, laboratory, and field studies. While geostatistical inverse methods offer many advantages, one key disadvantage is its highly parameterized nature, which renders it computationally intensive for large-scale problems. Another issue is that geostatistics-based HT may produce overly smooth images of subsurface heterogeneity when there are few monitoring interval data. Therefore, some may question the utility of the geostatistical inversion approach in certain situations and seek alternative approaches. To investigate these issues, we simultaneously calibrated different groundwater models with varying subsurface conceptualizations and parameter resolutions using a laboratory sandbox aquifer. The compared models included: (1) isotropic and anisotropic effective parameter models; (2) a heterogeneous model that faithfully represents the geological features; and (3) a heterogeneous model based on geostatistical inverse modeling. The performance of these models was assessed by quantitatively examining the results from model calibration and validation. Calibration data consisted of steady state drawdown data from eight pumping tests and validation data consisted of data from 16 separate pumping tests not used in the calibration effort. Results revealed that the geostatistical inversion approach performed the best among the approaches compared, although the geological model that faithfully represented stratigraphy came a close second. In addition, when the number of pumping tests available for inverse modeling was small, the geological modeling approach yielded more robust validation results. This suggests that better knowledge of stratigraphy obtained via geophysics or other means may contribute to improved results for HT.
NASA Astrophysics Data System (ADS)
Zhang, Donghua; Madsen, Henrik; Ridler, Marc E.; Refsgaard, Jens C.; Jensen, Karsten H.
2015-12-01
The ensemble Kalman filter (EnKF) is a popular data assimilation (DA) technique that has been extensively used in environmental sciences for combining complementary information from model predictions and observations. One of the major challenges in EnKF applications is the description of model uncertainty. In most hydrological EnKF applications, an ad hoc model uncertainty is defined with the aim of avoiding a collapse of the filter. The present work provides a systematic assessment of model uncertainty in DA applications based on combinations of forcing, model parameters, and state uncertainties. This is tested in a case where groundwater hydraulic heads are assimilated into a distributed and integrated catchment-scale model of the Karup catchment in Denmark. A series of synthetic data assimilation experiments are carried out to analyse the impact of different model uncertainty assumptions on the feasibility and efficiency of the assimilation. The synthetic data used in the assimilation study makes it possible to diagnose model uncertainty assumptions statistically. Besides the model uncertainty, other factors such as observation error, observation locations, and ensemble size are also analysed with respect to performance and sensitivity. Results show that inappropriate definition of model uncertainty can greatly degrade the assimilation performance, and an appropriate combination of different model uncertainty sources is advised.
NASA Astrophysics Data System (ADS)
Simpson, Scott C.; Meixner, Thomas
2012-02-01
Flood events can induce temporal changes in streambed elevation and particle-size composition, which may influence the bed's hydraulic properties and stream-aquifer fluxes during and after an event. This study combines a set of previously developed modeling approaches to create a synthetic flood event during which bed sediment is entrained and deposited as a function of hydraulic conditions and particle size. One simulated river reach in a state of approximate dynamic equilibrium is chosen to investigate the impacts of size-selective sediment transport on stream-aquifer interaction. Along this reach, the preferential entrainment of fine sediment during the flood's rising limb leads to overall bed coarsening, and increases in vertical hydraulic conductivity (Kbv) and downward fluxes of floodwater into the streambed. Progressively finer sediment layers are deposited during the event's falling limb, causing the redevelopment of a colmation (clogging) layer on the bed surface and a decline in overall Kbv by the event's conclusion. This reduction in Kbv leads to prolonged retention of event water in the streambed (after the reach reverts from losing to gaining river conditions) when compared with what is expected if pre-event Kbv values are used to estimate river-aquifer exchanges. This process of sequential bed coarsening and fining during a flood event provides a mechanistic explanation for the event size-and-duration threshold, inferred in some systems, that must be exceeded for significant amounts of flood recharge to occur. The major consequences of these processes—enhanced infiltration and prolonged floodwater retention—have potentially major implications for groundwater-surface water interactions, water quality, contaminant transport, and riparian biogeochemistry.
Numerical Simulation Study on the Hydraulic Behavior in Closed Fractures
NASA Astrophysics Data System (ADS)
Zhang, C.; Shu, L.; Wen, Z.; Wu, G.; Zhang, X.; Hu, B. X.
2015-12-01
As the main repositories for groundwater in karst systems, fractures involve the movement and storage of groundwater. Fundamentally, Navier-Stokes (NS) Equations is used to describe flow in fractures. However, due to the complexity of Navier-Stokes Equations, it is rarely applied to solve fracture flow problems. Thus, based on some simplifications, Stokes equations, Reynolds equations and Cubic Law (CL) are derived to describe fracture flow. The validity of the three simplified equations were extensively studies. Among the three simplified equations of NS, CL is the simplest and used to describe flow in open, smooth and paralleled fractures. In the previous work, most researchers focused on the open fractures. But it's the closed fractures exist widely in the field not the open fractures. The objective of this paper aims to check the validity of CL in closed fractures with different apertures and widths of fracture. After comparing the experimental results and simulations results from the COMSOL Multiphysics (FEM), this software was applied to solve the 3D or 2D NS equations in the closed fractures. The results obtained from NS simulation results and calculation results from CL were compared to indicate the degree of the validity of CL in application. A critical velocity was proposed to illustrate the validity of CL in closed fractures. Furthermore, the impacts of aperture size, width of fracture size, and velocity magnitude on both the hydraulic conductivity and velocity profile were also analyzed. The results showed the CL was capable of describing flow in closed fractures when the velocity was less than the critical velocity varying from 0.02 to 30.08cm/s. The ratio of NS results and CL results was between 0.9 and 2, with velocity varying from 0 to 40cm/s. The discrepancy between NS equation and CL increased with Reynolds number, increased with aperture size and decreased with width of fracture.
NASA Astrophysics Data System (ADS)
Blessent, Daniela; Therrien, René; Lemieux, Jean-Michel
2011-12-01
This paper presents numerical simulations of a series of hydraulic interference tests conducted in crystalline bedrock at Olkiluoto (Finland), a potential site for the disposal of the Finnish high-level nuclear waste. The tests are in a block of crystalline bedrock of about 0.03 km3 that contains low-transmissivity fractures. Fracture density, orientation, and fracture transmissivity are estimated from Posiva Flow Log (PFL) measurements in boreholes drilled in the rock block. On the basis of those data, a geostatistical approach relying on a transitional probability and Markov chain models is used to define a conceptual model based on stochastic fractured rock facies. Four facies are defined, from sparsely fractured bedrock to highly fractured bedrock. Using this conceptual model, three-dimensional groundwater flow is then simulated to reproduce interference pumping tests in either open or packed-off boreholes. Hydraulic conductivities of the fracture facies are estimated through automatic calibration using either hydraulic heads or both hydraulic heads and PFL flow rates as targets for calibration. The latter option produces a narrower confidence interval for the calibrated hydraulic conductivities, therefore reducing the associated uncertainty and demonstrating the usefulness of the measured PFL flow rates. Furthermore, the stochastic facies conceptual model is a suitable alternative to discrete fracture network models to simulate fluid flow in fractured geological media.
Thermal-hydraulic modeling needs for passive reactors
Kelly, J.M.
1997-07-01
The U.S. Nuclear Regulatory Commission has received an application for design certification from the Westinghouse Electric Corporation for an Advanced Light Water Reactor design known as the AP600. As part of the design certification process, the USNRC uses its thermal-hydraulic system analysis codes to independently audit the vendor calculations. The focus of this effort has been the small break LOCA transients that rely upon the passive safety features of the design to depressurize the primary system sufficiently so that gravity driven injection can provide a stable source for long term cooling. Of course, large break LOCAs have also been considered, but as the involved phenomena do not appear to be appreciably different from those of current plants, they were not discussed in this paper. Although the SBLOCA scenario does not appear to threaten core coolability - indeed, heatup is not even expected to occur - there have been concerns as to the performance of the passive safety systems. For example, the passive systems drive flows with small heads, consequently requiring more precision in the analysis compared to active systems methods for passive plants as compared to current plants with active systems. For the analysis of SBLOCAs and operating transients, the USNRC uses the RELAP5 thermal-hydraulic system analysis code. To assure the applicability of RELAP5 to the analysis of these transients for the AP600 design, a four year long program of code development and assessment has been undertaken.
NASA Astrophysics Data System (ADS)
Iden, Sascha C.; Peters, Andre; Durner, Wolfgang
2015-11-01
The prediction of unsaturated hydraulic conductivity from the soil water retention curve by pore-bundle models is a cost-effective and widely applied technique. One problem for conductivity predictions from retention functions with continuous derivatives, i.e. continuous water capacity functions, is that the hydraulic conductivity curve exhibits a sharp drop close to water saturation if the pore-size distribution is wide. So far this artifact has been ignored or removed by introducing an explicit air-entry value into the capillary saturation function. However, this correction leads to a retention function which is not continuously differentiable. We present a new parameterization of the hydraulic properties which uses the original saturation function (e.g. of van Genuchten) and introduces a maximum pore radius only in the pore-bundle model. In contrast to models using an explicit air entry, the resulting conductivity function is smooth and increases monotonically close to saturation. The model concept can easily be applied to any combination of retention curve and pore-bundle model. We derive closed-form expressions for the unimodal and multimodal van Genuchten-Mualem models and apply the model concept to curve fitting and inverse modeling of a transient outflow experiment. Since the new model retains the smoothness and continuous differentiability of the retention model and eliminates the sharp drop in conductivity close to saturation, the resulting hydraulic functions are physically more reasonable and ideal for numerical simulations with the Richards equation or multiphase flow models.
Upper and Middle Tiete River Basin dam-hydraulic system, travel time and temperature modeling
NASA Astrophysics Data System (ADS)
Devkota, Bishnu; Imberger, Jörg
2012-12-01
SummaryTiete River System in the State of Sao Paolo, Brazil is characterized by complex hydraulics and operational problems due to series of dams and point and diffuse inflows along the river. A one dimension Lagrangian river model was developed and applied to the 313 km reach of the Upper and Middle Tiete River Basin from the Penha Dam to the head water of Bara Bonita Reservoir, a stretch of river that includes six small to medium size dams (3.4-22 m high) including the Pirapora Reservoir and 26 inflows into the river (11 tributaries, 9 diffuse source areas, and discharges of 4 cities stormwater and 2 wastewater treatment plants. The conservative tracer transport and temperature model that accounts for the short and long wave radiation and heat transfers at the free surface was included and solved using the Crank-Nicholson scheme. The time variable catchment input to the model was the simulated output of the external hydrological model called Runoff Load Model which results were provided by CETESB. The numerical treatment of series of dams and spillway (that included uncontrolled overflow spillway, gate-controlled ogee spillway; and underflow gates and tunnels) and parameterisation of hydraulic jumps are described. Special attention was focused on the high spatial and temporal variation of flows in Tiete River Basin, a result of the large variation in catchment inflows and channel geometry due to dams and reservoirs along the river. Predicted and measured spatial and seasonal variation of flow and temperature profiles along the river show good agreement. The simulated travel time of conservative tracer is compared against the CETESB's 1982 and 1984 field study data in a 254 km reach of the Middle Tiete River that again shows good agreement. Being Lagrangian in construction, this new model is computationally efficient making it an ideal tool for long term simulation for water resource planning, management and operation decision making in a large and complex river
Bumgarner, Johnathan R; McCray, John E
2007-06-01
During operation of an onsite wastewater treatment system, a low-permeability biozone develops at the infiltrative surface (IS) during application of wastewater to soil. Inverse numerical-model simulations were used to estimate the biozone saturated hydraulic conductivity (K(biozone)) under variably saturated conditions for 29 wastewater infiltration test cells installed in a sandy loam field soil. Test cells employed two loading rates (4 and 8cm/day) and 3 IS designs: open chamber, gravel, and synthetic bundles. The ratio of K(biozone) to the saturated hydraulic conductivity of the natural soil (K(s)) was used to quantify the reductions in the IS hydraulic conductivity. A smaller value of K(biozone)/K(s,) reflects a greater reduction in hydraulic conductivity. The IS hydraulic conductivity was reduced by 1-3 orders of magnitude. The reduction in IS hydraulic conductivity was primarily influenced by wastewater loading rate and IS type and not by the K(s) of the native soil. The higher loading rate yielded greater reductions in IS hydraulic conductivity than the lower loading rate for bundle and gravel cells, but the difference was not statistically significant for chamber cells. Bundle and gravel cells exhibited a greater reduction in IS hydraulic conductivity than chamber cells at the higher loading rates, while the difference between gravel and bundle systems was not statistically significant. At the lower rate, bundle cells exhibited generally lower K(biozone)/K(s) values, but not at a statistically significant level, while gravel and chamber cells were statistically similar. Gravel cells exhibited the greatest variability in measured values, which may complicate design efforts based on K(biozone) evaluations for these systems. These results suggest that chamber systems may provide for a more robust design, particularly for high or variable wastewater infiltration rates. PMID:17449084
Otten, Alexander; van Vuuren, Wieke; Stienen, Arno; van Asseldonk, Edwin; Schouten, Alfred; van der Kooij, Herman
2011-01-01
Robotics used for diagnostic measurements on, e.g. stroke survivors, require actuators that are both stiff and compliant. Stiffness is required for identification purposes, and compliance to compensate for the robots dynamics, so that the subject can move freely while using the robot. A hydraulic actuator can act as a position (stiff) or a torque (compliant) actuator. The drawback of a hydraulic actuator is that it behaves nonlinear. This article examines two methods for controlling a nonlinear hydraulic actuator. The first method that is often applied uses an elastic element (i.e. spring) connected in series with the hydraulic actuator so that the torque can be measured as the deflection of the spring. This torque measurement is used for proportional integral control. The second method of control uses the inverse of the model of the actuator as a linearizing controller. Both methods are compared using simulation results. The controller designed for the series elastic hydraulic actuator is faster to implement, but only shows good performance for the working range for which the controller is designed due to the systems nonlinear behavior. The elastic element is a limiting factor when designing a position controller due to its low torsional stiffness. The model-based controller linearizes the nonlinear system and shows good performance when used for torque and position control. Implementing the model-based controller does require building and validating of the detailed model. PMID:22275654
NASA Astrophysics Data System (ADS)
Chen, Yi-Feng; Hu, Shao-Hua; Hu, Ran; Zhou, Chuang-Bing
2015-04-01
High-pressure packer test (HPPT) is an enhanced constant head packer test for characterizing the permeability of fractured rocks under high-pressure groundwater flow conditions. The interpretation of the HPPT data, however, remains difficult due to the transition of flow conditions in the conducting structures and the hydraulic fracturing-induced permeability enhancement in the tested rocks. In this study, a number of HPPTs were performed in the sedimentary and intrusive rocks located at 450 m depth in central Hainan Island. The obtained Q-P curves were divided into a laminar flow phase (I), a non-Darcy flow phase (II), and a hydraulic fracturing phase (III). The critical Reynolds number for the deviation of flow from linearity into phase II was 25-66. The flow of phase III occurred in sparsely to moderately fractured rocks, and was absent at the test intervals of perfect or poor intactness. The threshold fluid pressure between phases II and III was correlated with RQD and the confining stress. An Izbash's law-based analytical model was employed to calculate the hydraulic conductivity of the tested rocks in different flow conditions. It was demonstrated that the estimated hydraulic conductivity values in phases I and II are basically the same, and are weakly dependent on the injection fluid pressure, but it becomes strongly pressure dependent as a result of hydraulic fracturing in phase III. The hydraulic conductivity at different test intervals of a borehole is remarkably enhanced at highly fractured zone or contact zone, but within a rock unit of weak heterogeneity, it decreases with the increase of depth.
NASA Astrophysics Data System (ADS)
Bermudez, Maria; Neal, Jeffrey C.; Bates, Paul D.; Coxon, Gemma; Freer, Jim E.; Cea, Luis; Puertas, Jeronimo
2015-04-01
Flood inundation models require appropriate boundary conditions to be specified at the limits of the domain, which commonly include upstream flow rate and water stage at the downstream boundary of the reach being modelled. These data are usually acquired from gauging stations of the river network, where measured water levels are converted to discharge by means of a rating curve. Although generally treated as deterministic relationships, rating curves are subject to uncertainties which can limit their accuracy and applicability to any particular event, particularly for low-frequency events, where the rating curve relationship is sparsely sampled and needs to be extrapolated. Consequently, derived streamflow records can be particularly uncertain and/or problematic for high-return period estimates. Moreover, the limited number of gauges in reach-scale studies usually requires flow to be routed from upstream gauge stations to the boundary of the model domain, which introduces additional uncertainty. This is more likely to be the case with complex modelling approaches, in which only limited areas can be simulated due to computational cost and model setup requirements. In this study, a method to incorporate the rating curve uncertainty into the predictions of a reach-scale flood inundation model is proposed. Firstly, the uncertainty in the rating curves is quantified using a non-parametric local weighted regression approach and uncertainty bounds for discharge and water elevations at gauging locations are generated. A regional simplified-physics hydraulic model is then applied to combine these uncertainties and generate an ensemble of discharge and water elevation time series at the boundaries of a local-scale high complexity hydraulic model. Finally, the impact of this uncertainty on the local model performance is evaluated using flood extent data and measured water levels within the local model area. The local-scale model is applied to 7 km of the river Severn passing
Daly, Keith R.; Mooney, Sacha J.; Bennett, Malcolm J.; Crout, Neil M. J.; Roose, Tiina; Tracy, Saoirse R.
2015-01-01
Understanding the dynamics of water distribution in soil is crucial for enhancing our knowledge of managing soil and water resources. The application of X-ray computed tomography (CT) to the plant and soil sciences is now well established. However, few studies have utilized the technique for visualizing water in soil pore spaces. Here this method is utilized to visualize the water in soil in situ and in three-dimensions at successive reductive matric potentials in bulk and rhizosphere soil. The measurements are combined with numerical modelling to determine the unsaturated hydraulic conductivity, providing a complete picture of the hydraulic properties of the soil. The technique was performed on soil cores that were sampled adjacent to established roots (rhizosphere soil) and from soil that had not been influenced by roots (bulk soil). A water release curve was obtained for the different soil types using measurements of their pore geometries derived from CT imaging and verified using conventional methods, such as pressure plates. The water, soil, and air phases from the images were segmented and quantified using image analysis. The water release characteristics obtained for the contrasting soils showed clear differences in hydraulic properties between rhizosphere and bulk soil, especially in clay soil. The data suggest that soils influenced by roots (rhizosphere soil) are less porous due to increased aggregation when compared with bulk soil. The information and insights obtained on the hydraulic properties of rhizosphere and bulk soil will enhance our understanding of rhizosphere biophysics and improve current water uptake models. PMID:25740922
Daly, Keith R; Mooney, Sacha J; Bennett, Malcolm J; Crout, Neil M J; Roose, Tiina; Tracy, Saoirse R
2015-04-01
Understanding the dynamics of water distribution in soil is crucial for enhancing our knowledge of managing soil and water resources. The application of X-ray computed tomography (CT) to the plant and soil sciences is now well established. However, few studies have utilized the technique for visualizing water in soil pore spaces. Here this method is utilized to visualize the water in soil in situ and in three-dimensions at successive reductive matric potentials in bulk and rhizosphere soil. The measurements are combined with numerical modelling to determine the unsaturated hydraulic conductivity, providing a complete picture of the hydraulic properties of the soil. The technique was performed on soil cores that were sampled adjacent to established roots (rhizosphere soil) and from soil that had not been influenced by roots (bulk soil). A water release curve was obtained for the different soil types using measurements of their pore geometries derived from CT imaging and verified using conventional methods, such as pressure plates. The water, soil, and air phases from the images were segmented and quantified using image analysis. The water release characteristics obtained for the contrasting soils showed clear differences in hydraulic properties between rhizosphere and bulk soil, especially in clay soil. The data suggest that soils influenced by roots (rhizosphere soil) are less porous due to increased aggregation when compared with bulk soil. The information and insights obtained on the hydraulic properties of rhizosphere and bulk soil will enhance our understanding of rhizosphere biophysics and improve current water uptake models. PMID:25740922
NASA Astrophysics Data System (ADS)
Baroni, G.; Facchi, A.; Gandolfi, C.; Ortuani, B.; Horeschi, D.; van Dam, J. C.
2009-06-01
Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. These data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different outputs (i.e. evapotranspiration, water content in the root zone, fluxes through the bottom boundary of the root zone) of two hydrological models with different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental field. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. Three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek; d) HYPRES; and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June-October 2006. Results showed a wide range of soil hydraulic parameter values evaluated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter α of the Van Genuchten curve. This is reflected in a variability of the modeling results which is, as expected, different for
NASA Astrophysics Data System (ADS)
Baroni, G.; Facchi, A.; Gandolfi, C.; Ortuani, B.; Horeschi, D.; van Dam, J. C.
2010-02-01
Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. The data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different variables (i.e. evapotranspiration, average water content in the root zone, flux at the bottom boundary of the root zone) simulated by two hydrological models of different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental profile. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. The remaining three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek, d) HYPRES, and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June - October 2006. Results showed a wide range of soil hydraulic parameter values generated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter α of the van Genuchten curve. This is reflected in a variability of the modeling results which is
NASA Astrophysics Data System (ADS)
Takeda, M.; Hiratsuka, T.; Ito, K.
2008-12-01
Precise characterization and modeling of groundwater flow systems are necessary for realistic performance assessments of radioactive waste disposal. In groundwater flow modeling, the gravity potential is commonly assumed as the dominant driving force of regional groundwater flow. However, the gravity potential flow model may have a limited ability to reconstruct the excess fluid pressure distributions occasionally observed in low-permeability formations. To improve groundwater flow models, geologic processes such as compaction disequilibrium, tectonic forces and diagenetic reactions have been invoked to reconstruct excess pressures. On the other hand, chemical osmosis has recently been considered as one of the driving forces of groundwater flow and a factor causing excess pressures in clay-rich formations with vertical salinity gradients. If a formation medium acts as a semi-permeable membrane, chemical osmosis induces a fluid movement in the direction of increasing salinity. Consequently, fluid pressure could increase where salinity is high and decrease where it is low. Thus, chemical osmosis could induce a fluid flow countering the pressure-driven flow in the formation. When osmotic- and pressure-driven flows equilibrate, the net flow ceases while the fluid pressures remain in disequilibrium. This means that the direction of groundwater flow might be misinterpreted without differentiating osmotically-induced pressure from those induced by other causes. However, the formation media are not perfect membranes, as they allow solute diffusion that accord to the salinity gradients. As a result, osmotic pressure would dissipate as the solutes diffuse from high to low concentrations. That means the time period during which the osmotic pressures are held in the formation depends on hydraulic and diffusive properties other than the chemico-osmotic property of the formation media. The osmotic pressures have indeed been observed in natural formations, and the chemico
NASA Astrophysics Data System (ADS)
Gayler, Sebastian; Salima-Sultana, Daisy; Selle, Benny; Ingwersen, Joachim; Wizemann, Hans-Dieter; Högy, Petra; Streck, Thilo
2016-04-01
Soil water extraction by roots affects the dynamics and distribution of soil moisture and controls transpiration, which influences soil-vegetation-atmosphere feedback processes. Consequently, root water uptake requires close attention when predicting water fluxes across the land surface, e.g., in agricultural crop models or in land surface schemes of weather and climate models. The key parameters for a successful simultaneous simulation of soil moisture dynamics and evapotranspiration in Richards equation-based models are the soil hydraulic parameters, which describe the shapes of the soil water retention curve and the soil hydraulic conductivity curve. As measurements of these parameters are expensive and their estimation from basic soil data via pedotransfer functions is rather inaccurate, the values of the soil hydraulic parameters are frequently inversely estimated by fitting the model to measured time series of soil water content and evapotranspiration. It is common to simulate root water uptake and transpiration by simple stress functions, which describe from which soil layer water is absorbed by roots and predict when total crop transpiration is decreased in case of soil water limitations. As for most of the biogeophysical processes simulated in crop and land surface models, there exist several alternative functional relationships for simulating root water uptake and there is no clear reason for preferring one process representation over another. The error associated with alternative representations of root water uptake, however, contributes to structural model uncertainty and the choice of the root water uptake model may have a significant impact on the values of the soil hydraulic parameters estimated inversely. In this study, we use the agroecosystem model system Expert-N to simulate soil moisture dynamics and evapotranspiration at three agricultural field sites located in two contrasting regions in Southwest Germany (Kraichgau, Swabian Alb). The Richards
NASA Astrophysics Data System (ADS)
Wissmeier, L.; Barry, D. A.
2009-08-01
The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards' equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC's comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil's hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.
Penny, G.S.; Conway, M.W.
1996-02-01
The primary objective of this project is to provide laboratory data that is pertinent to designing hydraulic fracturing treatments for coalbed methane. Coal fluid interactions studies, fracture conductivity, fluid leak-off through cleats, rheology, and proppant transport are designed to respresent Black Warrior and San Juan treatments. A second objective is to apply the information learned in laboratory testing to actual hydraulic fracturing treatments in order to improve results. A final objective is to review methods currently used to catalog well performance following hydraulic fracturing for the purpose of placing the data in a useable database that can be accessed by users to determine the success of various treatment scenarios.
Numerical Hydraulic Study on Seawater Cooling System of Combined Cycle Power Plant
NASA Astrophysics Data System (ADS)
Kim, J. Y.; Park, S. M.; Kim, J. H.; Kim, S. W.
2010-06-01
As the rated flow and pressure increase in pumping facilities, a proper design against surges and severe cavitations in the pipeline system is required. Pressure surge due to start-up, shut-down process and operation failure causes the water hammer in upstream of the closing valve and the cavitational hammer in downstream of the valve. Typical cause of water hammer is the urgent closure of valves by breakdown of power supply and unexpected failure of pumps. The abrupt changes in the flow rate of the liquid results in high pressure surges in upstream of the valves, thus kinetic energy is transformed into potential energy which leads to the sudden increase of the pressure that is called as water hammer. Also, by the inertia, the liquid continues to flow downstream of the valve with initial speed. Accordingly, the pressure decreases and an expanding vapor bubble known as column separation are formed near the valve. In this research, the hydraulic study on the closed cooling water heat exchanger line, which is the one part of the power plant, is introduced. The whole power plant consists of 1,200 MW combined power plant and 220,000 m3/day desalination facility. Cooling water for the plant is supplied by sea water circulating system with a capacity of 29 m3/s. The primary focus is to verify the steady state hydraulic capacity of the system. The secondary is to quantify transient issues and solutions in the system. The circuit was modeled using a commercial software. The stable piping network was designed through the hydraulic studies using the simulation for the various scenarios.
NASA Astrophysics Data System (ADS)
Bour, O.; Klepikova, M.; Le Borgne, T.; De Dreuzy, J.
2013-12-01
Inverse modeling of hydraulic and geometrical properties of fractured media is a very challenging objective due to the spatial heterogeneity of the medium and the scarcity of data. Here we present a flow tomography approach that permits to characterize the location, the connectivity and the hydraulic properties of main flow paths in fractured media. The accurate characterization of the location, hydraulic properties and connectivity of major fracture zones is essential to model flow and solute transport in fractured media. Cross-borehole flowmeter tests, which consist of measuring changes in vertical borehole flows when pumping a neighboring borehole, were shown to be an efficient technique to provide information on the properties of the flow zones that connect borehole pairs [Paillet, 1998; Le Borgne et al., 2006]. The interpretation of such experiments may however be quite uncertain when multiple connections exist. In this study, we explore the potential of flow tomography (i.e., sequential cross-borehole flowmeter tests) for characterizing aquifer heterogeneity. We first propose a framework for inverting flow and drawdown data to infer fracture connectivity and transmissivities. Here we use a simplified discrete fracture network approach that highlights main connectivity structures. This conceptual model attempts to reproduce fracture network connectivity without taking fracture geometry (length, orientation, dip) into account. We then explore the potential of the method for simplified synthetic fracture network models and quantify the sensitivity of drawdown and borehole flow velocities to the transmissivity of the connecting flowpaths. Flow tomography is expected to be most effective if cross-borehole pumping induces large changes in vertical borehole velocities. The uncertainty of the transmissivity estimates increases for small borehole flow velocities. The uncertainty about the transmissivity of fractures that connect the main flowpath but not the boreholes
A novel approach to model hydraulic and electrical conductivity in fractal porous media
NASA Astrophysics Data System (ADS)
Ghanbarian, B.; Daigle, H.; Sahimi, M.
2014-12-01
Accurate prediction of conductivity in partially-saturated porous media has broad applications in various phenomena in porous media, and has been studied intensively since the 1940s by petroleum, chemical and civil engineers, and hydrologists. Many of the models developed in the past are based on the bundle of capillary tubes. In addition, pore network models have also been developed for simulating multiphase fluid flow in porous media and computing the conductivity in unsaturated porous media. In this study, we propose a novel approach using concepts from the effective-medium approximation (EMA) and percolation theory to model hydraulic and electrical conductivity in fractal porous media whose pore-size distributions exhibit power-law scaling. In our approach, the EMA, originally developed for predicting electrical conductivity of composite materials, is used to predict the effective conductivity, from complete saturation to some intermediate water content that represents a crossover point. Below the crossover water content, but still above a critical saturation (percolation threshold), a universal scaling predicted by percolation theory, a power law that expresses the dependence of the conductivity on the water content (less a critical water saturation) with an exponent of 2, is invoked to describe the effective conductivity. In order to evaluate the accuracy of the approach, experimental data were used from the literature. The predicted hydraulic conductivities for most cases are in excellent agreement with the data. In a few cases the theory underestimates the hydraulic conductivities, which correspond to porous media with very broad pore-size distribution in which the largest pore radius is more than 7 orders of magnitude greater than the smallest one. The approach is also used to predict the saturation dependence of the electrical conductivity for experiments in which capillary pressure data are available. The results indicate that the universal scaling of
Specification of Flow Conditions in the Mathematical Model of Hydraulic Damper
NASA Astrophysics Data System (ADS)
Svoboda, Rudolf; Škliba, Jan; Matějec, Radek
Hydraulic dampers represent one of the basic instruments for absorption of vibration in dynamic systems. The damper is substituted either by its hydraulic velocity characteristic, or directly by a mathematical model as a dynamic subsystem. The standard damper model does not provide satisfactory results especially those concerning the strokes of all four damper valves. To improve these results it is necessary to simulate the flows through valves more precisely and, last but not least, to set adequately correct values to all essential parameters of both mechanical and hydraulic parts of the damper. In the paper is presented a new, corrected formula for discharge flow coefficients based on the measurements of flow characteristics of throttle elements of the damper with constant as well as variable slot width. The experimental equipment used for identification process is described as well.
NASA Astrophysics Data System (ADS)
Ding, Fei; Han, Xu; Luo, Zhen; Zhang, Nong
2012-12-01
In this paper, a new hydraulically interconnected suspension (HIS) system is proposed for the implementation of a resistance control for the pitch and bounce modes of tri-axle heavy trucks. A lumped-mass half-truck model is established using the free-body diagram method. The equations of motion of a mechanical and hydraulic coupled system are developed by incorporating the hydraulic strut forces into the mechanical subsystem as externally applied forces. The transfer matrix method (TMM) is used to evaluate the impedance matrix of the hydraulic subsystem consisting of models of fluid pipes, damper valves, accumulators, and three-way junctions. The TMM is further applied to find the quantitative relationships between the hydraulic strut forces and boundary flow of the mechanical-fluid interactive subsystem. The modal analysis method is employed to perform the vibration analysis between the trucks with the conventional suspension and the proposed HIS. Comparison analysis focuses on free vibration with identified eigenvalues and eigenvectors, isolation vibration capacity, and force vibration in terms of the power spectrum density responses. The obtained results show the effectiveness of the proposed HIS system in reducing the pitch motion of sprung mass and simultaneously maintaining the ride comfort. The pitch stiffness is increased while the bounce stiffness is slightly softened. The peak values of sprung mass and wheel hop motions are greatly reduced, and the vibration decay rate of sprung mass is also significantly increased.
Nonlinear dynamic modeling for smart material electro-hydraulic actuator development
NASA Astrophysics Data System (ADS)
Larson, John P.; Dapino, Marcelo J.
2013-03-01
Smart material electro-hydraulic actuators use hydraulic rectification by one-way check valves to amplify the motion of smart materials, such as magnetostrictives and piezoelectrics, in order to create compact, lightweight actuators. A piston pump driven by a smart material is combined with a hydraulic cylinder to form a self-contained, power-by-wire actuator that can be used in place of a conventional hydraulic system without the need for hydraulic lines and a centralized pump. The performance of an experimental actuator driven by a 12.7 mm diameter, 114 mm length Terfenol-D rod is evaluated over a range of applied input frequencies, loads, and currents. The peak performance achieved is 37 W, moving a 220 N load at a rate of 17 cm/s and producing a blocked pressure of 12.5 MPa. Additional tests are conducted to quantify the dynamic behavior of the one-way reed valves using a scanning laser vibrometer to identify the frequency response of the reeds and the effect of the valve seat and fluid mass loading. A lumped-parameter model is developed for the system that includes valve inertia and fluid response nonlinearities, and the model results are compared with the experimental data.
Centrifugal slurry pump wear and hydraulic studies. Phase II report. Experimental studies
Mistry, D.; Cooper, P.; Biswas, C.; Sloteman, D.; Onuschak, A.
1983-01-01
This report describes the work performed by Ingersoll-Rand Research, Inc., under Phase II, Experimental Studies for the contract entitled, Centrifugal Slurry Pump Wear and Hydraulic Studies. This work was carried out for the US Department of Energy under Contract No. DE-AC-82PC50035. The basic development approach pursued this phase is presented, followed by a discussion on wear relationships. The analysis, which resulted in the development of a mathematical wear model relating pump life to some of the key design and operating parameters, is presented. The results, observations, and conclusions of the experimental investigation on small scale pumps that led to the selected design features for the prototype pump are discussed. The material investigation was performed at IRRI, ORNL and Battelle. The rationale for selecting the materials for testing, the test methods and apparatus used, and the results obtained are presented followed by a discussion on materials for a prototype pump. In addition, the prototype pump test facility description, as well as the related design and equipment details, are presented. 20 references, 53 figures, 13 tables.
NASA Astrophysics Data System (ADS)
Freedman, V. L.; Bacon, D. H.; Saripalli, K. P.; Meyer, P. D.
2001-12-01
Precipitation and dissolution of minerals in the subsurface can cause a significant reduction in porosity and permeability by plugging pore throats in aquifer and reservoir media. Changes in these two basic properties of the medium also result in significant changes in the remaining 'derived properties' (i.e., relative permeability, fluid-fluid and fluid-solid interfacial areas, pore and particle size distributions) and the constitutive relationships among these properties. Very few published works on modeling the influence of chemical reactions and fluid flow on porosity and permeability account for the spatial and temporal changes in the hydrologic properties on flow and transport. This study reports on the development of a methodology for modeling changes in permeability of unsaturated sub-surface media due to glass and mineral precipitation and dissolution reactions using a film depositional modeling approach. The model is based on the assumption that the mineral precipitate is deposited on the pore walls as a continuous film, causing a reduction in permeability. In this study, the film depositional model is developed for a discrete pore-size distribution, which is determined using the unsaturated hydraulic properties of the porous medium. This facilitates the process of dynamically updating the unsaturated hydraulic parameters used to describe fluid flow through the media. The resulting algorithms are implemented in the multiphase, multicomponent reactive transport code STORM (Sub-surface Transport over Reactive Multiphases). The modeling approach is tested using the Hanford's Immobilized Low Activity Waste (ILAW) repository, where the low-level waste from fuel fabrication activities is being vitrified and emplaced in the sub-surface. Results from simulation of the simultaneous dissolution of low-level glassified waste and secondary mineral precipitation show that the film depositional model based on the Mualem approach reasonably predicts permeability changes
Controls on peat hydraulic conductivity: bridging the gap between models and measurements
NASA Astrophysics Data System (ADS)
Morris, Paul J.; Baird, Andy J.; Belyea, Lisa R.
2015-04-01
Peat hydraulic properties such as saturated hydraulic conductivity (Ksat) are variable in space and time. A broadly accepted but poorly understood relationship between peat decomposition and Ksat is highly influential on the behaviour of models of peatland development, but existing data do not allow this relationship to be parameterised satisfactorily. Previous empirical studies have typically used simple metrics of time-integrated decomposition such as fibre content, von Post score or light transmission (humification) to predict Ksat. By contrast, peatland development models represent the state of peat decomposition as a more abstract fraction of remaining mass. As such, a gap exists between the requirements of simulation models and the available empirical data - a gap that we seek to address. We collected fourteen shallow (~ 0.5 m long, 0.1 m diameter) peat cores from a Swedish raised bog. Our sampling locations comprised two vegetation microhabitats, hummocks and hollows (n = 7 for each), combined factorially with two locations - the flat, treeless, central bog plateau (n = 8) and the treed, sloping bog margin (n = 6). In the laboratory we split the cores into 0.06 m depth intervals and measured horizontal Ksat and dry bulk peat carbon to nitrogen concentration ratios (C:N). Following a published method we used these C:N values to approximate the state of peat decomposition (fractional mass remaining) for depth intervals deeper than ~ 0.1 m. We used a step-up procedure to fit a linear mixed-effects (LME) model to the data, so as to predict log10 (Ksat) from two continuous predictor variables, depth and our C:N-derived approximation of fractional mass remaining; and two categorical predictors, microhabitat (hummock, hollow) and location (central plateau, bog margin). The LME model predicts that log10 (Ksat) declines linearly with increasing depth below the surface (p < 0.001) and with increasing peat decomposition (p = 0.012); and that it is greater in hummock
NASA Astrophysics Data System (ADS)
Idzikowski, A.; Salamon, S.
2013-06-01
A general characteristics of a car hydraulic braking system (CHBS) is presented in this publication. A graphical model of properties-component objects is developed for the above-mentioned system. Moreover, four mathematical models in terms of logic, the set theory and the Boolean algebra of Boolean functions are developed. The examination is ended with a general model of the CHBS for n - Boolean variables and the construction and mathematical-technical interpretation of this model is presented.
A Model for Hydraulic Properties Based on Angular Pores with Lognormal Size Distribution
NASA Astrophysics Data System (ADS)
Durner, W.; Diamantopoulos, E.
2014-12-01
Soil water retention and unsaturated hydraulic conductivity curves are mandatory for modeling water flow in soils. It is a common approach to measure few points of the water retention curve and to calculate the hydraulic conductivity curve by assuming that the soil can be represented as a bundle of capillary tubes. Both curves are then used to predict water flow at larger spatial scales. However, the predictive power of these curves is often very limited. This can be very easily illustrated if we measure the soil hydraulic properties (SHPs) for a drainage experiment and then use these properties to predict the water flow in the case of imbibition. Further complications arise from the incomplete wetting of water at the solid matrix which results in finite values of the contact angles between the solid-water-air interfaces. To address these problems we present a physically-based model for hysteretic SHPs. This model is based on bundles of angular pores. Hysteresis for individual pores is caused by (i) different snap-off pressures during filling and emptying of single angular pores and (ii) by different advancing and receding contact angles for fluids that are not perfectly wettable. We derive a model of hydraulic conductivity as a function of contact angle by assuming flow perpendicular to pore cross sections and present closed-form expressions for both the sample scale water retention and hydraulic conductivity function by assuming a log-normal statistical distribution of pore size. We tested the new model against drainage and imbibition experiments for various sandy materials which were conducted with various liquids of differing wettability. The model described both imbibition and drainage experiments very well by assuming a unique pore size distribution of the sample and a zero contact angle for the perfectly wetting liquid. Eventually, we see the possibility to relate the particle size distribution with a model which describes the SHPs.
NASA Astrophysics Data System (ADS)
Wissmeier, L. C.; Barry, D. A.
2008-12-01
Precipitation/dissolution induces changes in the pore radii of water-filled pores, and, consequently, affects flow in porous media. The selective radius shift model was developed to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties of unsaturated soils. The model considers the dependency of the amount of mineral precipitation/dissolution within a pore on the local pore volume. Furthermore, it accounts for precipitation/dissolution taking place only in the water-filled part of the pore space. The pore bundle concept was used to relate the pore-scale process of dissolution/precipitation to changes in macroscopic soil hydraulic properties. In the numerical model, the finite change in mineral volume at a discrete time step leads to a discontinuous pore-size distribution, because only the water-filled pores are affected. This pore-size distribution is converted back to a discontinuous soil moisture characteristic to which, at every time step, a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with the selective radius shift model a head-based solution of Richards' equation for aqueous phase flow was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC's comprehensive geochemical reactions. The model was applied to kinetic halite dissolution and calcite precipitation as a consequence of cation exchange in a variety of unsaturated flow situations. The applications showed marked changes in the soil's hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on the water content. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.
NASA Astrophysics Data System (ADS)
Zemenkova, M. Yu; Shipovalov, A. N.; Zemenkov, Yu D.
2016-04-01
The main technological equipment of pipeline transport of hydrocarbons are hydraulic machines. During transportation of oil mainly used of centrifugal pumps, designed to work in the “pumping station-pipeline” system. Composition of a standard pumping station consists of several pumps, complex hydraulic piping. The authors have developed a set of models and algorithms for calculating system reliability of pumps. It is based on the theory of reliability. As an example, considered one of the estimation methods with the application of graph theory.
Modeling and control of a hydraulically actuated flexible-prismatic link robot
Love, L.; Kress, R.; Jansen, J.
1996-12-01
Most of the research related to flexible link manipulators to date has focused on single link, fixed length, single plane of vibration test beds. In addition, actuation has been predominantly based upon electromagnetic motors. Ironically, these elements are rarely found in the existing industrial long reach systems. This manuscript describes a new hydraulically actuated, long reach manipulator with a flexible prismatic link at Oak Ridge National Laboratory (ORNL). Focus is directed towards both modeling and control of hydraulic actuators as well as flexible links that have variable natural frequencies.
Seismic studies of a massive hydraulic fracturing experiment
House, L.; Keppler, H.; Kaieda, H.
1985-01-01
During a massive hydraulic fracturing experiment carried out at Fenton Hill, New Mexico, 850 microearthquakes, ranging in magnitudes from -3 to 0, were located reliably using arrival times recorded at a set of 5 downhole geophone stations. A subset of these events were located using an upgraded hodogram technique. The seismicity defines a tabular zone with horizontal extent of 900 m, vertical extent of 800 m, and thickness of 150 m. This zone strikes N340/sup 0/E, and dips 75/sup 0/ to the east; its position indicates that no hydraulic connection between the two predrilled wells could be achieved by the fracturing. The distribution of locations obtained from arrival times shows good agreement with those derived from hodograms. Well constrained fault plane solutions were determined for 26 of the larger microearthquakes observed at a surface seismic net. Most solutions display one nearly vertical nodal plane that strikes close to N - S, and a T axis that trends roughly E - W, in agreement with regional indicators of the least principal stress direction. 9 refs., 6 figs.
NASA Astrophysics Data System (ADS)
Montanari, M.; Hostache, R.; Matgen, P.; Schumann, G.; Pfister, L.; Hoffmann, L.
2008-11-01
Two of the most relevant components of any flood forecasting system, namely the rainfall-runoff and flood inundation models, increasingly benefit from the availability of spatially distributed Earth Observation data. With the advent of microwave remote sensing instruments and their all weather capabilities, new opportunities have emerged over the past decade for improved hydrologic and hydraulic model calibration and validation. However, the usefulness of remote sensing observations in coupled hydrologic and hydraulic models still requires further investigations. Radar remote sensing observations are readily available to provide information on flood extent. Moreover, the fusion of radar imagery and high precision digital elevation models allows estimating distributed water levels. With a view to further explore the potential offered by SAR images, this paper investigates the usefulness of remote sensing-derived water stages in a modelling sequence where the outputs of hydrologic models (rainfall-runoff models) serve as boundary condition of flood inundation models. The methodology consists in coupling a simplistic 3-parameter conceptual rainfall-runoff model with a 1-D flood inundation model. Remote sensing observations of flooded areas help to identify and subsequently correct apparent volume errors in the modelling chain. The updating of the soil moisture module of the hydrological model is based on the comparison of water levels computed by the coupled hydrologic-hydraulic model with those estimated using remotely sensed flood flood extent. The potential of the proposed methodology is illustrated with data collected during a storm event of the Alzette River (Grand-Duchy of Luxembourg). The study contributes to assessing the value of remote sensing data for evaluating the saturation status of a river basin.
NASA Astrophysics Data System (ADS)
Montanari, M.; Hostache, R.; Matgen, P.; Schumann, G.; Pfister, L.; Hoffmann, L.
2009-03-01
Two of the most relevant components of any flood forecasting system, namely the rainfall-runoff and flood inundation models, increasingly benefit from the availability of spatially distributed Earth Observation data. With the advent of microwave remote sensing instruments and their all weather capabilities, new opportunities have emerged over the past decade for improved hydrologic and hydraulic model calibration and validation. However, the usefulness of remote sensing observations in coupled hydrologic and hydraulic models still requires further investigations. Radar remote sensing observations are readily available to provide information on flood extent. Moreover, the fusion of radar imagery and high precision digital elevation models allows estimating distributed water levels. With a view to further explore the potential offered by SAR images, this paper investigates the usefulness of remote sensing-derived water stages in a modelling sequence where the outputs of hydrologic models (rainfall-runoff models) serve as boundary condition of flood inundation models. The methodology consists in coupling a simplistic 3-parameter conceptual rainfall-runoff model with a 1-D flood inundation model. Remote sensing observations of flooded areas help to identify and subsequently correct apparent volume errors in the modelling chain. The updating of the soil moisture module of the hydrologic model is based on the comparison of water levels computed by the coupled hydrologic-hydraulic model with those estimated using remotely sensed flood extent. The potential of the proposed methodology is illustrated with data collected during a storm event on the Alzette River (Grand-Duchy of Luxembourg). The study contributes to assess the value of remote sensing data for evaluating the saturation status of a river basin.
Multiscale modelling of hydraulic conductivity in vuggy porous media
Daly, K. R.; Roose, T.
2014-01-01
Flow in both saturated and non-saturated vuggy porous media, i.e. soil, is inherently multiscale. The complex microporous structure of the soil aggregates and the wider vugs provides a multitude of flow pathways and has received significant attention from the X-ray computed tomography (CT) community with a constant drive to image at higher resolution. Using multiscale homogenization, we derive averaged equations to study the effects of the microscale structure on the macroscopic flow. The averaged model captures the underlying geometry through a series of cell problems and is verified through direct comparison to numerical simulations of the full structure. These methods offer significant reductions in computation time and allow us to perform three-dimensional calculations with complex geometries on a desktop PC. The results show that the surface roughness of the aggregate has a significantly greater effect on the flow than the microstructure within the aggregate. Hence, this is the region in which the resolution of X-ray CT for image-based modelling has the greatest impact. PMID:24511248
NASA Astrophysics Data System (ADS)
White, J. R.; Wang, H.; Jawitz, J. W.; Sees, M. D.
2004-12-01
The Orlando Easterly Wetland (OEW), the largest municipal treatment wetland in Florida, began operation in 1987 mainly for reducing nutrient loads in tertiary treated domestic wastewater produced by the city of Orlando. After more than ten years of operation, a decrease in total P removal effectiveness has occurred since 1999, even though the effluent concentration of the wetland has remained below the permitted limit of 0.2 mg/L,. Hydraulic inefficiency in the wetland, especially in the front-end cells of the north flow train, was identified as a primary cause of the reduced treatment effectiveness. In order to improve the hydraulic performance of the OEW and maintain its efficient phosphorus treatment, a rejuvenation program (including muck removal followed by re-vegetation) was initiated on the front-end cells of the north flow train in 2002. The effectiveness of this activity for the improvement of hydraulic performance was evaluated with a tracer test and subsequent moment and model analyses for the tracer resident time distribution (RTDs). Results were compared to similar tracer tests conducted prior to rejuvenation activities. The models included one-path tank-in-series (TIS), two-path TIS, one-dimensional transport with inflow and storage (OTIS), plug flow with dispersion (PFD), and plug flow with fractional dispersion (PFFD). The hydraulic performance was characterized by both wetland hydraulic efficiency and the spreading of tracers. The results demonstrated that the rejuvenation considerably improved the hydraulic performance in the restored area. Also presented is a comparison of the wetland response between both bromide and lithium tracers, and the determination of the complete moments of residence time distributions (RTD) in cell-network wetlands.
Thermal-hydraulic modeling of the Pennsylvania State University Breazeale Nuclear Reactor (PSBR)
NASA Astrophysics Data System (ADS)
Chang, Jong E.
2005-11-01
Earlier experiments determined that the Pennsylvania State University Breazeale Nuclear Reactor (PSBR) core is cooled, not by an axial flow, but rather by a strong cross flow due to the thermal expansion of the coolant. To further complicate the flow field, a nitrogen-16 (N-16) pump was installed above the PSBR core to mix the exiting core buoyant thermal plume in order to delay the rapid release of radioactive N-16 to the PSBR pool surface. Thus, the interaction between the N-16 jet flow and the buoyancy driven flow complicates the analysis of the flow distribution in the PSBR pool. The main objectives of this study is to model the thermal-hydraulic behavior of the PSBR core and pool. During this study four major things were performed including the Computational Fluid Dynamics (CFD) model for the PSBR pool, the stand-alone fuel rod model for a PSBR fuel rod, the velocity measurements in and around the PSBR core, and the temperature measurements in the PSBR pool. Once the flow field was predicted by the CFD model, the measurement devices were manufactured and calibrated based on the CFD results. The major contribution of this study is to understand and to explain the flow behavior in the PSBR subchannels and pool using the FLOW3D model. The stand-alone dynamic fuel rod model was developed to determine the temperature distribution inside a PSBR fuel rod. The stand-alone fuel rod model was coupled to the FLOW3D model and used to predict the temperature behavior during steady-state and pulsing. The heat transfer models in the stand-alone fuel rod code are used in order to overcome the disadvantage of the CFD code, which does not calculate the mechanical stress, the gap conductance, and the two phase heat transfer. (Abstract shortened by UMI.)
Ecohydrological streamflow distributions and hydraulic food chain models
NASA Astrophysics Data System (ADS)
Ceola, S.; Botter, G.; Bertuzzo, E.; Mari, L.; Rodriguez-Iturbe, I.; Rinaldo, A.
2010-12-01
A comprehensive probabilistic characterization of streamflow variability in river basins has noteworthy scientific and social implications due to the relevant impacts on in-stream biogeochemical processes, human exploitations of stream water and ecological services of riparian and riverine environments. To this aim, a comparative study of ecohydrological models of streamflow probability distributions (pdfs), p(Q), against field data gathered in different hydrological contexts is described. Streamflows measured in several catchments across various climatic regions of North-Eastern Italy and the United States are employed. The relevance of the work stems from the implied analytical predictive ability of hydrologic variability, whose role on stream and riparian ecological processes and large-scale management schemes is fundamental. The tools employed are analytical models of p(Q) (and of the related flow duration curve, D(Q)) derived by coupling suitable storage-discharge relations with a stochastic description of streamflow production through soil moisture dynamics, and are expressed as a function of few macroscopic rainfall, soil, vegetation and geomorphological parameters. The performances of a version of the model (which includes the effects of nonlinear subsurface storage-discharge relations) are compared with those provided by its linear version through the analysis of 13 test catchments belonging to various climatic and geomorphic contexts. A general agreement between predicted and observed daily streamflows pdfs is shown, though differences emerge between the linear and the nonlinear approaches. In particular, by including the effects of a nonlinear storage-discharge relation the model accuracy is shown to increase with respect to the linear scheme in most examined cases. We show that this is not simply attributable to the added parameter but corresponds to a proper likelihood increase. Inferences on the nonlinear character of the relation between subsurface
NASA Astrophysics Data System (ADS)
Vitillo, F.; Vitale Di Maio, D.; Galati, C.; Caruso, G.
2015-11-01
A CFD analysis has been carried out to study the thermal-hydraulic behavior of liquid metal coolant in a fuel assembly of triangular lattice. In order to obtain fast and accurate results, the isotropic two-equation RANS approach is often used in nuclear engineering applications. A different approach is provided by Non-Linear Eddy Viscosity Models (NLEVM), which try to take into account anisotropic effects by a nonlinear formulation of the Reynolds stress tensor. This approach is very promising, as it results in a very good numerical behavior and in a potentially better fluid flow description than classical isotropic models. An Anisotropic Shear Stress Transport (ASST) model, implemented into a commercial software, has been applied in previous studies, showing very trustful results for a large variety of flows and applications. In the paper, the ASST model has been used to perform an analysis of the fluid flow inside the fuel assembly of the ALFRED lead cooled fast reactor. Then, a comparison between the results of wall-resolved conjugated heat transfer computations and the results of a decoupled analysis using a suitable thermal wall-function previously implemented into the solver has been performed and presented.
Interpretation of hydraulic model outputs in supporting ecologially-led river restoration.
NASA Astrophysics Data System (ADS)
Gillies, E.; Moir, H. J.
2014-12-01
In river systems, hydrodynamic forces are a major driver for geomorphic change, and a major influence on aquatic habitat. Representative modelling of channel hydraulics is therefore an invaluable tool in ecologically-led river restoration design, enabling a quantitative and objective assessment of complex processes that are essential to achieve biophysical objectives. Hydraulic modelling can form part of an iterative design process, utilised to indicate 'design performance' through a wide range of descriptors (e.g. physical heterogeneity, hydraulic micro-habitat) afforded by each stage of the design. However, it is important that the limitations of any computational fluid dynamic approach (e.g., 2D depth-averaged simulations) are well communicated to other specialists (e.g., geomorphologists, ecologists), managers, regulators and clients. One aspect of hydraulics where this is specifically important is in the use of vorticity and coherent vortex structures as indicators of habitat suitability and the probability of regions of scour or deposition. It has long been recognised that eddies structures provide important physical and ecological function in rivers but calculating vorticity changes using hydraulic models is relatively new in river hydraulics. The use of such calculations is growing in both academia and industry, especially as restoration approaches such as engineered log jams introduce significant vorticity into the flow. However, the introduction of vorticity needs to be modelled with very high resolutions near solid boundaries, and the convected vortex structures themselves are inherently three dimensional. Neither of these is routinely captured in river hydraulic models. Specifically, the interpretation of vorticity patterns or coherent vortex structures from 2D depth averaged river models must be treated with caution, or provided with further interpretation. We present work that demonstrates how coherent vortex structures and vorticity have been used to
NASA Astrophysics Data System (ADS)
Matgen, P.; Montanari, M.; Hostache, R.; Pfister, L.; Hoffmann, L.; Plaza, D.; Pauwels, V. R. N.; de Lannoy, G. J. M.; de Keyser, R.; Savenije, H. H. G.
2010-09-01
With the onset of new satellite radar constellations (e.g. Sentinel-1) and advances in computational science (e.g. grid computing) enabling the supply and processing of multi-mission satellite data at a temporal frequency that is compatible with real-time flood forecasting requirements, this study presents a new concept for the sequential assimilation of Synthetic Aperture Radar (SAR)-derived water stages into coupled hydrologic-hydraulic models. The proposed methodology consists of adjusting storages and fluxes simulated by a coupled hydrologic-hydraulic model using a Particle Filter-based data assimilation scheme. Synthetic observations of water levels, representing satellite measurements, are assimilated into the coupled model in order to investigate the performance of the proposed assimilation scheme as a function of both accuracy and frequency of water level observations. The use of the Particle Filter provides flexibility regarding the form of the probability densities of both model simulations and remote sensing observations. We illustrate the potential of the proposed methodology using a twin experiment over a widely studied river reach located in the Grand-Duchy of Luxembourg. The study demonstrates that the Particle Filter algorithm leads to significant uncertainty reduction of water level and discharge at the time step of assimilation. However, updating the storages of the model only improves the model forecast over a very short time horizon. A more effective way of updating thus consists in adjusting both states and inputs. The proposed methodology, which consists in updating the biased forcing of the hydraulic model using information on model errors that is inferred from satellite observations, enables persistent model improvement. The present schedule of satellite radar missions is such that it is likely that there will be continuity for SAR-based operational water management services. This research contributes to evolve reactive flood management into
Global river flood hazard maps: hydraulic modelling methods and appropriate uses
NASA Astrophysics Data System (ADS)
Townend, Samuel; Smith, Helen; Molloy, James
2014-05-01
Flood hazard is not well understood or documented in many parts of the world. Consequently, the (re-)insurance sector now needs to better understand where the potential for considerable river flooding aligns with significant exposure. For example, international manufacturing companies are often attracted to countries with emerging economies, meaning that events such as the 2011 Thailand floods have resulted in many multinational businesses with assets in these regions incurring large, unexpected losses. This contribution addresses and critically evaluates the hydraulic methods employed to develop a consistent global scale set of river flood hazard maps, used to fill the knowledge gap outlined above. The basis of the modelling approach is an innovative, bespoke 1D/2D hydraulic model (RFlow) which has been used to model a global river network of over 5.3 million kilometres. Estimated flood peaks at each of these model nodes are determined using an empirically based rainfall-runoff approach linking design rainfall to design river flood magnitudes. The hydraulic model is used to determine extents and depths of floodplain inundation following river bank overflow. From this, deterministic flood hazard maps are calculated for several design return periods between 20-years and 1,500-years. Firstly, we will discuss the rationale behind the appropriate hydraulic modelling methods and inputs chosen to produce a consistent global scaled river flood hazard map. This will highlight how a model designed to work with global datasets can be more favourable for hydraulic modelling at the global scale and why using innovative techniques customised for broad scale use are preferable to modifying existing hydraulic models. Similarly, the advantages and disadvantages of both 1D and 2D modelling will be explored and balanced against the time, computer and human resources available, particularly when using a Digital Surface Model at 30m resolution. Finally, we will suggest some
Analytical Study of Cavitation Surge in a Hydraulic System.
Kang, Donghyuk; Yokota, Kazuhiko
2014-10-01
In order to clarify effects of an accumulator, pipe lengths and gradients of pressure and suction performances on cavitation surge, one-dimensional stability analyses of cavitation surge were performed in hydraulic systems consisting of an upstream tank, an inlet pipe, a cavitating pump, a downstream pipe, and a downstream tank. An accumulator located upstream or downstream of the cavitating pump was included in the analysis. Increasing the distance between the upstream accumulator and the cavitating pump enlarged the stable region. On the other hand, decreasing the distance between the downstream accumulator and the cavitating pump enlarged the stable region. Furthermore, the negative gradient of a suction performance curve and the positive gradient of a pressure performance curve cause cavitation surge. PMID:25278638
NASA Astrophysics Data System (ADS)
Małoszewski, Piotr; Wachniew, Przemysław; Czupryński, Piotr
2006-12-01
SummaryCombined use of tracers and mathematical modelling for evaluation of hydraulic characteristics of constructed wetlands is presented for the subsurface-flow system with Phragmites australis in Nowa Słupia (Poland). Instantaneously injected bromide and tritium tracers were used to obtain residence time distributions of wastewaters in three parallel inhomogeneous gravel cells of the wetland. The multi flow dispersion model, which assumes the existence of several flow-paths with different hydraulic properties was developed using the respective parallel combination of analytical solutions of the one-dimensional advection-dispersion equation. The model was used successfully to fit the experimental tracer breakthrough curves. Different flow components were identified and wastewater volumes, water-saturated porosity, mean wastewater travel times, longitudinal dispersivities as well as hydraulic conductivity of wetland cells were derived from model parameters. The variation in flow components and apparent hydraulic characteristics among wetland cells relate to the improper design and maintenance of the wetland. The single fissure dispersion model, which assumes possible diffusion of tracers into the zones with stagnant water during convective-dispersive flow through the mobile zone is adopted to the research conditions and used to model the TBC-s for one cell. The results show that this model can be calibrated with the satisfactory accuracy in that cell but yields unacceptable values of some parameters.
The Effect of Loading Rate on Hydraulic Fracturing in Synthetic Granite - a Discrete Element Study
NASA Astrophysics Data System (ADS)
Tomac, I.; Gutierrez, M.
2015-12-01
Hydraulic fracture initiation and propagation from a borehole in hard synthetic rock is modeled using the two dimensional Discrete Element Method (DEM). DEM uses previously established procedure for modeling the strength and deformation parameters of quasi-brittle rocks with the Bonded Particle Model (Itasca, 2004). A series of simulations of laboratory tests on granite in DEM serve as a reference for synthetic rock behavior. Fracturing is enabled by breaking parallel bonds between DEM particles as a result of the local stress state. Subsequent bond breakage induces fracture propagation during a time-stepping procedure. Hydraulic fracturing occurs when pressurized fluid induces hoop stresses around the wellbore which cause rock fracturing and serves for geo-reservoir permeability enhancement in oil, gas and geothermal industries. In DEM, a network of fluid pipes and reservoirs is used for mathematical calculation of fluid flow through narrow channels between DEM particles, where the hydro-mechanical coupling is fully enabled. The fluid flow calculation is superimposed with DEM stress-strain calculation at each time step. As a result, the fluid pressures during borehole pressurization in hydraulic fracturing, as well as, during the fracture propagation from the borehole, can be simulated. The objective of this study is to investigate numerically a hypothesis that fluid pressurization rate, or the fluid flow rate, influences upon character, shape and velocity of fracture propagation in rock. The second objective is to better understand and define constraints which are important for successful fracture propagation in quasi-brittle rock from the perspective of flow rate, fluid density, viscosity and compressibility relative to the rock physical properties. Results from this study indicate that not only too high fluid flow rates cause fracture arrest and multiple fracture branching from the borehole, but also that the relative compressibility of fracturing fluid and
NASA Astrophysics Data System (ADS)
Hu, R.; Brauchler, R.; Herold, M.; Bayer, P.; Sauter, M.
2009-04-01
the hydraulic gradient does not. By this trick, transient data can be analyzed with the computational efficiency of a steady state model, which proceeds hundreds of times faster than transient models. Finally, a specific storage distribution can be calculated from the diffusivity and hydraulic conductivity reconstructions derived from travel time and steady shape inversion. The groundwork of this study is the aquifer-analogue study from BAYER (1999), in which six parallel profiles of a natural sedimentary body with a size of 16m x 10m x 7m were mapped in high resolution with respect to structural and hydraulic parameters. Based on these results and using geostatistical interpolation methods, MAJI (2005) designed a three dimensional hydraulic model with a resolution of 5cm x 5cm x 5cm. This hydraulic model was used to simulate a large number of short term pumping tests in a tomographical array. The high resolution parameter reconstructions gained from the inversion of simulated pumping test data demonstrate that the proposed inversion scheme allows reconstructing the individual architectural elements and their hydraulic properties with a higher resolution compared to conventional hydraulic and geological investigation methods. Bayer P (1999) Aquifer-Analog-Studium in grobklastischen braided river Ablagerungen: Sedimentäre/hydrogeologische Wandkartierung und Kalibrierung von Georadarmessungen, Diplomkartierung am Lehrstuhl für Angewandte Geologie, Universität Tübingen, 25 pp. Maji, R. (2005) Conditional Stochastic Modelling of DNAPL Migration and Dissolution in a High-resolution Aquifer Analog, Ph.D. thesis at the University of Waterloo, 187 pp.
Impact of representation of hydraulic structures in modelling a Severn barrage
NASA Astrophysics Data System (ADS)
Bray, Samuel; Ahmadian, Reza; Falconer, Roger A.
2016-04-01
In this study, enhancements to the numerical representation of sluice gates and turbines were made to the hydro-environmental model Environmental Fluid Dynamics Code (EFDC), and applied to the Severn Tidal Power Group Cardiff-Weston Barrage. The extended domain of the EFDC Continental Shelf Model (CSM) allows far-field hydrodynamic impact assessment of the Severn Barrage, pre- and post-enhancement, to demonstrate the importance of accurate hydraulic structure representation. The enhancements were found to significantly affect peak water levels in the Bristol Channel, reducing levels by nearly 1 m in some areas, and even affect predictions as far-field as the West Coast of Scotland, albeit to a far lesser extent. The model was tested for sensitivity to changes in the discharge coefficient, Cd, used in calculating discharge through sluice gates and turbines. It was found that the performance of the Severn Barrage is not sensitive to changes to the Cd value, and is mitigated through the continual, rather than instantaneous, discharge across the structure. The EFDC CSM can now be said to be more accurately predicting the impacts of tidal range proposals, and the investigation of sensitivity to Cd improves the confidence in the modelling results, despite the uncertainty in this coefficient.
Investigation of approximations in thermal-hydraulic modeling of core conversions
Garner, Patrick L.; Hanan, Nelson A.
2008-07-15
Neutronics analyses for core conversions are usually fairly detailed, for example representing all 4 flats and all 4 corners of all 6 tubes of all 20 IRT-3M or -4M fuel assemblies in the core of the VVR-SM reactor in Uzbekistan. The coupled neutronics and thermal-hydraulic analysis for safety analysis transients is usually less detailed, for example modeling only a hot and an average fuel plate and the associated coolant. Several of the approximations have been studied using the RELAP5 and PARET computer codes in order to provide assurance that the lack of full detail is not important to the safety analysis. Two specific cases studied are (1) representation of a core of same- type fuel assemblies by a hot and an average assembly each having multiple channels as well as by merely a hot and average channel and (2) modeling a core containing multiple fuel types as the sum of fractional core models for each fuel type. (author)
Rock deformation models and fluid leak-off in hydraulic fracturing
NASA Astrophysics Data System (ADS)
Yarushina, Viktoriya M.; Bercovici, David; Oristaglio, Michael L.
2013-09-01
Fluid loss into reservoir rocks during hydraulic fracturing is modelled via a poro-elastoplastic pressure diffusion equation in which the total compressibility is a sum of fluid, rock and pore space compressibilities. Inclusion of pore compressibility and porosity-dependent permeability in the model leads to a strong pressure dependence of leak-off (i.e. drainage rate). Dilation of the matrix due to fluid invasion causes higher rates of fluid leak-off. The present model is appropriate for naturally fractured and tight gas reservoirs as well as for soft and poorly consolidated formations whose mechanical behaviour departs from simple elastic laws. Enhancement of the leak-off coefficient by dilation, predicted by the new model, may help explain the low percentage recovery of fracturing fluid (usually between 5 and 50 per cent) in shale gas stimulation by hydraulic fracturing.
Thermal hydraulic modeling of the mock fuel facility
NASA Astrophysics Data System (ADS)
Gardner, Jacob
The major focus of this thesis was to make improved three dimensional models of the Mock Fuel Facility. Three distinct experiment types run with the Mock Fuel Facility (MFF) were the main focus of this thesis. Two of the experiments were modeled and an in-depth analysis of the model results was performed to gain a better understanding of the Mock Fuel Facility. For the third experiment the process of creating a model was begun. There were multiple purposes for the work completed in this thesis. The work was done partially to gain a greater understanding of the UMass Lowell Research Reactor (UMLRR). There is minimal instrumentation within the UMLRR to measure localized temperatures within the UMLRR. It is hoped that the work done in this thesis will provide a basis for future modeling work which will give insight into the temperature profiles within the UMLRR. This work is also being done to gain insight into the capabilities of the COMSOL multiphysics modelling software and evaluate its potential for future modelling work. Finally this work is also being done for its potential as an educational tool. The MFF and COMSOL have potential to be used for experimental lab work by students to learn about computer modeling and validation.
NASA Astrophysics Data System (ADS)
Battiston, Stéphanie; Allenbach, Bernard
2010-05-01
The exceptional characteristics of the December 2003 Rhône flood event (particularly high water flows, extent of the affected area, important damages especially in the region of Arles) make it be considered as a reference flood episode of this French river and a very well-known event. During the crisis, the International Charter "Space and Major Disasters" was triggered by the French Civil Protection for the rapid mapping of the flooding using Earth Observation imagery in order to facilitate crisis operations. As a result, more than 60 satellite images covering the flood were acquired over a 10 days period following the peak flow. Using the opportunity provided by this incomparable data coverage, the French Ministry of the Environment ordered a study on the evaluation of remote sensing's potential benefits for flood management. One of the questions asked by the risk managers was: what type of flood information can be provided by the different remote sensing platforms? Elements of response were delivered mainly in the form of a comprehensive compilation of maps and illustrations, displaying the main hydraulic elements (static ones as well as dynamic ones), initially listed and requested by hydrologists (more precisely, by a regional engineering society specialised in hydraulics and hydrology and in charge of a field campaign during the event), observed on different optical images of the flood event having affected the plain between Tarascon (upstream) and Arles (downstream). It is seen that a careful mapping of all flood traces visible on remote sensing event imagery - apparent water, moisture traces, breaches, overflows, stream directions, impermeable boundaries … - delivers a valuable vision of the flood's occurrence combining accuracy and comprehensiveness. In fact, optical imagery offers a detailed vision of the event : moisture traces complete flood traces extent; the observation of draw-off directions through waterproof barriers reveals hydraulic
NASA Astrophysics Data System (ADS)
von Boetticher, Albrecht; Turowski, Jens M.
2013-04-01
Sediment transport plays the key role in the forming fluvial systems, and it has become a major issue for the management and operation of hydropower plants and shipping channels as well as for the understanding of flooding hazards. The interactions between grain, hydraulics, and local topography are complicated, and it is difficult to quantify them in the laboratory, in the field, or in simple theoretical models based on two-dimensional force balances. In many theoretical models of incipient motion, grains are assumed to have simple geometrical shapes (usually spheres), despite the fact that in the real world grain shapes can vary from platy to elongated to symmetrical, and from angular to rounded. In addition, grain shape may influence protrusion and angle of repose for a grain sitting in the same pocket, and both local grain environment and grain shape in turn influence the local flow field around the grain, and thus the forces acting on it. In this project a numerical model will be developed that is capable of resolving the interaction between the grain and the flow. The model will be used to study the influence of the grain environment and grain shape on initiation of motion. The aims of the project are three-fold. First, we want to develop a numerical model that is adapted to and can efficiently deal with the specific problems occurring in the study of initiation of bedload sediment motion in streams. Second, we want to identify which of the parameters describing grain shape and local topography (e.g., protrusion, angle of repose, grain angularity) dominantly influence the initiation of motion of an individual grain. Third, we want to clarify when and where simple analytical approximations of the sort that have been used in many previous models are applicable. In three tasks we will develop the model, validate it against laboratory and field data, and perform a number of systematic numerical experiments to better understand incipient motion.
Prototype Data Models and Data Dictionaries for Hanford Sediment Physical and Hydraulic Properties
Rockhold, Mark L.; Last, George V.; Middleton, Lisa A.
2010-09-30
The Remediation Decision Support (RDS) project, managed by the Pacific Northwest National Laboratory (PNNL) for the U.S. Department of Energy (DOE) and the CH2M HILL Plateau Remediation Company (CHPRC), has been compiling physical and hydraulic property data and parameters to support risk analyses and waste management decisions at Hanford. In FY09 the RDS project developed a strategic plan for a physical and hydraulic property database. This report documents prototype data models and dictionaries for these properties and associated parameters. Physical properties and hydraulic parameters and their distributions are required for any type of quantitative assessment of risk and uncertainty associated with predictions of contaminant transport and fate in the subsurface. The central plateau of the Hanford Site in southeastern Washington State contains most of the contamination at the Site and has up to {approx}100 m of unsaturated and unconsolidated or semi-consolidated sediments overlying the unconfined aquifer. These sediments contain a wide variety of contaminants ranging from organic compounds, such as carbon tetrachloride, to numerous radionuclides including technetium, plutonium, and uranium. Knowledge of the physical and hydraulic properties of the sediments and their distributions is critical for quantitative assessment of the transport of these contaminants in the subsurface, for evaluation of long-term risks and uncertainty associated with model predictions of contaminant transport and fate, and for evaluating, designing, and operating remediation alternatives. One of the goals of PNNL's RDS project is to work with the Hanford Environmental Data Manager (currently with CHPRC) to develop a protocol and schedule for incorporation of physical property and hydraulic parameter datasets currently maintained by PNNL into HEIS. This requires that the data first be reviewed to ensure quality and consistency. New data models must then be developed for HEIS that are
USER'S MANUAL FOR EXPLORE-I: A RIVER BASIN WATER QUALITY MODEL (HYDRAULIC MODULE ONLY)
EXPLORE-1 is a computer program that simulates the dynamic hydraulic and water quality characteristics of a river basin. It can be used to study the effects of various flow conditions, waste discharges and treatment schemes on the water quality conditions of lakes, reservoirs, an...
NASA Astrophysics Data System (ADS)
Müller, Daniel; Regenspurg, Simona; Milsch, Harald; Blöcher, Guido; Kranz, Stefan; Saadat, Ali
2014-05-01
In aquifer thermal energy storage (ATES) systems, large amounts of energy can be stored by injecting hot water into deep or intermediate aquifers. In a seasonal production-injection cycle, water is circulated through a system comprising the porous aquifer, a production well, a heat exchanger and an injection well. This process involves large temperature and pressure differences, which shift chemical equilibria and introduce or amplify mechanical processes. Rock-fluid interaction such as dissolution and precipitation or migration and deposition of fine particles will affect the hydraulic properties of the porous medium and may lead to irreversible formation damage. In consequence, these processes determine the long-term performance of the ATES system and need to be predicted to ensure the reliability of the system. However, high temperature and pressure gradients and dynamic feedback cycles pose challenges on predicting the influence of the relevant processes. Within this study, a reservoir model comprising a coupled hydraulic-thermal-chemical simulation was developed based on an ATES demonstration project located in the city of Berlin, Germany. The structural model was created with Petrel, based on data available from seismic cross-sections and wellbores. The reservoir simulation was realized by combining the capabilities of multiple simulation tools. For the reactive transport model, COMSOL Multiphysics (hydraulic-thermal) and PHREEQC (chemical) were combined using the novel interface COMSOL_PHREEQC, developed by Wissmeier & Barry (2011). It provides a MATLAB-based coupling interface between both programs. Compared to using COMSOL's built-in reactive transport simulator, PHREEQC additionally calculates adsorption and reaction kinetics and allows the selection of different activity coefficient models in the database. The presented simulation tool will be able to predict the most important aspects of hydraulic, thermal and chemical transport processes relevant to
Antiquity versus modern times in hydraulics - a case study
NASA Astrophysics Data System (ADS)
Stroia, L.; Georgescu, S. C.; Georgescu, A. M.
2010-08-01
Water supply and water management in Antiquity represent more than Modern World can imagine about how people in that period used to think about, and exploit the resources they had, aiming at developing and improving their society and own lives. This paper points out examples of how they handled different situations, and how they managed to cope with the growing number of population in the urban areas, by adapting or by improving their water supply systems. The paper tries to emphasize the engineering contribution of Rome and the Roman Empire, mainly in the capital but also in the provinces, as for instance the today territory of France, by analysing some aqueducts from the point of view of modern Hydraulic Engineering. A third order polynomial regression is proposed to compute the water flow rate, based on the flow cross-sectional area measured in quinaria. This paper also emphasizes on contradictory things between what we thought we knew about Ancient Roman civilization, and what could really be proven, either by a modern engineering approach, a documentary approach, or by commonsense, where none of the above could be used. It is certain that the world we live in is the heritage of the Greco-Roman culture and therefore, we are due to acknowledge their contribution, especially taking into account the lack of knowledge of that time, and the poor resources they had.
NASA Astrophysics Data System (ADS)
Nguyen, P.; Sorooshian, S.; Hsu, K.; AghaKouchak, A.; Sanders, B. F.
2013-12-01
Flash floods are considered one of the most hazardous natural disasters, which kills thousands of people and causes billions of US dollar economic damages annually world-wide. Forecasting flash floods to provide accurate warnings in a timely manner is still challenging. At the Center for Hydrometeorology and Remote Sensing (CHRS) at the University of California, Irvine, we have been developing a coupled high resolution distributed hydrologic-hydraulic system for flash flood modeling which has been successfully tested for some selected areas in the U.S. and has potential to be implemented in global scale. The system employs the National Weather Service's distributed hydrologic model (HL-RDHM) as a rainfall-runoff generator, and a high-resolution hydraulic model (BreZo) for simulating the channel and flood-plain flows realistically. In this research, we evaluate the system for flash flood warning using multiple precipitation sources (gauge, radar and satellite and forecast). A flash flood event occurring on June 11, 2010 in the Upper Little Missouri River watershed in Arkansas is used as a case study. The catchment was delineated into 123 sub-catchments based on the 10m Digital Elevation Model (DEM) topography data from USGS. From HL-RDHM surface runoff, 123 hydrographs can be derived and connected as inputs to BreZo. The system was calibrated using NEXRAD Stage IV radar-based rainfall by tuning the roughness parameter in BreZo to best match the USGS discharge observation at the catchment outlet. The results show good agreement with the USGS gauge flow measurement (Nash-Sutcliffe coefficient = 0.91) when using Stage IV data. The system is under investigation with satellite-based precipitation data, rain gauge and Global Forecast System (GFS) data and will be reported in the presentation.
NASA Astrophysics Data System (ADS)
Couvreur, Valentin; Kandelous, Maziar; Mairesse, Harmony; Baram, Shahar; Moradi, Ahmad; Pope, Katrin; Hopmans, Jan
2015-04-01
Groundwater quality is specifically vulnerable in irrigated agricultural lands in California and many other (semi-)arid regions of the world. The routine application of nitrogen fertilizers with irrigation water in California is likely responsible for the high nitrate concentrations in groundwater, underlying much of its main agricultural areas. To optimize irrigation/fertigation practices, it is essential that irrigation and fertilizers are applied at the optimal concentration, place, and time to ensure maximum root uptake and minimize leaching losses to the groundwater. The applied irrigation water and dissolved fertilizer, root nitrate and water uptake interact with soil and root properties in a complex manner that cannot easily be resolved. It is therefore that coupled experimental-modelling studies are required to allow for unravelling of the relevant complexities that result from typical variations of crop properties, soil texture and layering across farmer-managed fields. A combined field monitoring and modelling approach was developed to quantify from simple measurements the leaching of water and nitrate below the root zone. The monitored state variables are soil water content within the root zone, soil matric potential below the root zone, and nitrate concentration in the soil solution. Plant and soil properties of incremented complexity are optimized with the software HYDRUS in an inverse modelling scheme, which allows estimating leaching under constraint of hydraulic principles. Questions of optimal irrigation and fertilization timing can then be addressed using predictive results and global optimization algorithms.
Combined seismic and hydraulic method of modeling flow in fractured low permeability rocks
Witherspoon, P.A.; Long, J.C.S.; Majer, E.L.; Myer, L.R.
1987-06-01
Modeling flow of ground water in hard rocks where a network of fractures provides the dominant flow paths is a major problem. This paper summarizes a program of investigations currently underway in this laboratory to characterize the geometry of fractured rocks and develop methods of handling flow in such systems. Numerical models have been developed to investigate flow behavior in two- and three-dimensional fracture networks. The results demonstrate the insights that can be gained from modeling studies of fractured rocks. A key problem is gathering the necessary data on fracture geometry. Investigations have been started to determine how vertical seismic profiling (VSP) might be improved and applied to this problem. A VSP experiment in The Geysers geothermal field in northern California, where fracture orientation is known, produced shear wave splitting and velocity anisotropy in agreement with theory. The results suggest the potential application of 3-component, multi-source VSP data in determining fracture orientation and average spacing. We believe a combination of seismic and hydraulic methods can greatly enhance an understanding of fluid flow and transport in low permeability rock systems where fractures provide the dominant paths. 40 refs, 16 figs., 4 tabs.
NASA Astrophysics Data System (ADS)
Lerat, Julien; Perrin, Charles; Andréassian, Vazken; Loumagne, Cécile; Ribstein, Pierre
2012-02-01
SummaryAccurate prediction of a flood inundation area constitutes an essential part of a flood forecasting system. When a river reach receives significant lateral inflows, flood inundation modelling requires the joint application of a hydrological model to calculate lateral inflows and a hydraulic model to calculate water levels along the river reach. In this study, we compared different strategies to couple the GR4J lumped rainfall-runoff model and the linearised diffusive wave propagation model. These strategies introduced variations in the nature of the connections between the two models using combinations of point and uniformly distributed lateral inflows. The coupled model was then applied to the Illinois River case study for which 10 years of hourly data were made available within the Distributed Model Intercomparison Project, now in phase 2. The simulations were assessed at the downstream end of the reach and at two interior points considered to be ungauged during the calibration process. The results first show that including uniformly distributed inflows made the coupled model more robust and stable compared to only using point flow input. A similar level of performance was reached with models using point inflows only, but at the cost of more uncertain parameters and less stable model performance when changing test periods. Second, identifying the optimal number of tributaries to be modelled individually by the hydrological model was easier when a combination of uniformly distributed and point inflows was used. In this case, model performance was less sensitive to the number of tributaries used and the inclusion of two or three tributaries appeared sufficient to obtain satisfactory performance for the simulations on the main channel. Last, the results on the main channel differed from those obtained on tributaries: overall performance was better on the main channel and required a lower degree of lateral inflow resolution, which suggests that upstream flow
NASA Astrophysics Data System (ADS)
Ni, C.-F.; Huang, Y.-J.; Dong, J.-J.; Yeh, T.-C. J.
2015-12-01
The transient hydraulic tomography survey (THTS) is a conceptually improved technique that efficiently estimates detailed variations in aquifer parameters. Based on the concept of the THTS, we developed a geostatistical inverse model to characterize saturated hydraulic conductivity (K) and the specific yield (Sy) in transient and unconfined aquifer systems. In this study, a synthetic example was first used to assess the accuracy of the developed inverse model. Multiple random K and Sy realizations with different variances of natural logarithm of K (lnK) were generated and systematically compared to evaluate the effects of joint inversion on K estimations. The model was implemented in field-scale, cross-hole injection tests in a shallow and highly permeable unconfined aquifer near the middle reaches of the Wu River in central Taiwan. To assess the effect of constant head boundary conditions on the estimation results, two additional modeling domains were evaluated on the basis of the same field data from the injection tests. The results of the synthetic example showed that the proposed inverse model can effectively reproduce the predefined K patterns and magnitudes. However, slightly less detail was obtained for the Sy field based on the sampling data from sequential transient hydraulic stresses. The joint inversion by using transient head observations could slightly decrease the accuracy of K estimations. The model implementation for field-scale injection tests showed that the model can estimate K and Sy fields with detailed spatial variations. Estimation results showed a relatively homogeneous aquifer for the tested well field. Results based on the three modeling domains showed similar patterns and magnitudes of K and Sy near the well locations. These results indicated that the THTS is relatively insensitive to artificially drawn boundary conditions even under transient conditions.
Hydraulic modeling for lahar hazards at cascades volcanoes
Costa, J.E.
1997-01-01
The National Weather Service flood routing model DAMBRK is able to closely replicate field-documented stages of historic and prehistoric lahars from Mt. Rainier, Washington, and Mt. Hood, Oregon. Modeled time-of-travel of flow waves are generally consistent with documented lahar travel-times from other volcanoes around the world. The model adequately replicates a range of lahars and debris flows, including the 230 million km3 Electron lahar from Mt. Rainier, as well as a 10 m3 debris flow generated in a large outdoor experimental flume. The model is used to simulate a hypothetical lahar with a volume of 50 million m3 down the East Fork Hood River from Mt. Hood, Oregon. Although a flow such as this is thought to be possible in the Hood River valley, no field evidence exists on which to base a hazards assessment. DAMBRK seems likely to be usable in many volcanic settings to estimate discharge, velocity, and inundation areas of lahars when input hydrographs and energy-loss coefficients can be reasonably estimated.
Model and analysis of a cylindrical in-line hydraulic suppressor with a solid compressible liner
NASA Astrophysics Data System (ADS)
Marek, K. A.; Earnhart, N. E.; Cunefare, K. A.
2014-12-01
An in-line hydraulic noise suppressor with a lossy, compressible liner made of foamed polyurethane liner is introduced which is intended to provide an alternative to current in-line silencing devices using compressed nitrogen gas volumes. The liner is engineered to be compressible at elevated pressures, such that it can provide effective noise abatement for practical hydraulic systems. In support of such work, a multimodal model is developed to characterize the device and the liner material. Because the hydraulic system is pressurized after insertion of the liner, the model must address liner compression and the corresponding small gaps introduced in the expansion volume; additionally, both compression and shear wave propagation must be considered in the liner. Several mode matching solutions are investigated, and a pseudoinverse mode matching method is found to provide good convergence characteristics. The multimodal model is validated against a finite element model, and also used in an optimization algorithm to estimate the material properties of a prototype liner using experimental transmission loss data. Experimental results show broadband transmission loss performance at 2.8 MPa system pressure; transmission loss decreases with increasing system pressure, and data at 4.1 MPa system pressure produces about 4 dB less transmission loss than a similarly sized commercial device. The multimodal model with estimated material properties at 2.8 MPa achieves a root mean squared error of 1.7 dB or less for two different length devices over a frequency range of 50-2000 Hz.
Inverse modeling of a multistep outflow experiment fordetermining hysteretic hydraulic properties
Faybishenko, B.; Finsterle, S.; Sonnenborg, T.O.
1998-05-01
A new, closed-form hysteretic model of the capillary pressure-saturation and relative permeability-saturation relationship has been implemented into ITOUGH2. The hysteretic capillary pressure function is based on the van Genuchten model, with a modified version of the dependent domain model of Mualem to describe the scanning curves. Hysteresis in the relative permeability relations is considered to be mainly a result of nonwetting fluid entrap- ment. The hysteresis model was used in combination with inverse modeling techniques to examine the potential of a simple drainage- imbibition experiment to determine hysteretic hydraulic properties.
The use of sediment deposition maps as auxiliary data for hydraulic model calibration
NASA Astrophysics Data System (ADS)
Mukolwe, Micah; Di Baldassarre, Giuliano; Solomatine, Dimitri
2013-04-01
One aspect of the French disaster mitigation setup is the statutory Risk Prevention Plans (PPR, Plans de Prévention des Risques); i.e. spatial identification of potential disasters and mitigation measures. The maps are categorised into three zones depicting increasing disaster severity and potential mitigation measures (RTM, 1999). Taking the example of the city of Barcelonnette, in South France (French Alps), floods have been the most frequent occurring natural hazard (Flageollet et al., 1996). Consequently, a case is put forward for the need for accurate flood extent delineation to support the decision making process. For this study, the Barcelonnette case study was considered, whereby the last devastating flooding was in June 1957 (Weber, 1994). Contrary to the recent advances in the proliferation of data to support flood inundation studies (Bates, 2012; Bates, 2004; Di Baldassarre and Uhlenbrook, 2012; Schumann et al., 2009), constraints are faced when analysing flood inundation events that occurred before the 1970's. In absence of frequent flooding, the analysis of historical flood extents may play an important role in shaping the awareness of local stakeholders and support land-use and urban planning. This study is part of a probabilistic flood mapping (e.g. Di Baldassarre et al., 2010, Horrit, 2006) of the valley carried out in a Monte-Carlo framework, while taking into account the peak flow and the parametric uncertainty. The simulations were carried out using the sub-grid channel model extension of the LISFLOOD-FP hydraulic model (Bates et al, 2010; Neal et al., 2012). Sediment deposition maps (Lecarpentier, 1963) were used to analyse the model performance, additionally the graduation of the sediment deposition sizes showed the flood propagation and was used to analyse the model runs. However, there still remains the challenge of quantifying the uncertainty in the sediment deposition map and the actual flood extent.
NASA Technical Reports Server (NTRS)
Gutmann, Ethan D.; Small, Eric E.
2007-01-01
Soil hydraulic properties (SHPs) regulate the movement of water in the soil. This in turn plays an important role in the water and energy cycles at the land surface. At present, SHPS are commonly defined by a simple pedotransfer function from soil texture class, but SHPs vary more within a texture class than between classes. To examine the impact of using soil texture class to predict SHPS, we run the Noah land surface model for a wide variety of measured SHPs. We find that across a range of vegetation cover (5 - 80% cover) and climates (250 - 900 mm mean annual precipitation), soil texture class only explains 5% of the variance expected from the real distribution of SHPs. We then show that modifying SHPs can drastically improve model performance. We compare two methods of estimating SHPs: (1) inverse method, and (2) soil texture class. Compared to texture class, inverse modeling reduces errors between measured and modeled latent heat flux from 88 to 28 w/m(exp 2). Additionally we find that with increasing vegetation cover the importance of SHPs decreases and that the van Genuchten m parameter becomes less important, while the saturated conductivity becomes more important.
Hydraulic conductivity study of compacted clay soils used as landfill liners for an acidic waste.
Hamdi, Noureddine; Srasra, Ezzeddine
2013-01-01
Three natural clayey soils from Tunisia were studied to assess their suitability for use as a liner for an acid waste disposal site. An investigation of the effect of the mineral composition and mechanical compaction on the hydraulic conductivity and fluoride and phosphate removal of three different soils is presented. The hydraulic conductivity of these three natural soils are 8.5 × 10(-10), 2.08 × 10(-9) and 6.8 × 10(-10)m/s for soil-1, soil-2 and soil-3, respectively. Soil specimens were compacted under various compaction strains in order to obtain three wet densities (1850, 1950 and 2050 kg/m(3)). In this condition, the hydraulic conductivity (k) was reduced with increasing density of sample for all soils. The test results of hydraulic conductivity at long-term (>200 days) using acidic waste solution (pH=2.7, charged with fluoride and phosphate ions) shows a decrease in k with time only for natural soil-1 and soil-2. However, the specimens of soil-2 compressed to the two highest densities (1950 and 2050 kg/m(3)) are cracked after 60 and 20 days, respectively, of hydraulic conductivity testing. This damage is the result of a continued increase in the internal stress due to the swelling and to the effect of aggressive wastewater. The analysis of anions shows that the retention of fluoride is higher compared to phosphate and soil-1 has the highest sorption capacity. PMID:22980909
A compartmental model to describe hydraulics in a full-scale waste stabilization pond.
Alvarado, Andres; Vedantam, Sreepriya; Goethals, Peter; Nopens, Ingmar
2012-02-01
The advancement of experimental and computational resources has facilitated the use of computational fluid dynamics (CFD) models as a predictive tool for mixing behaviour in full-scale waste stabilization pond systems. However, in view of combining hydraulic behaviour with a biokinetic process model, the computational load is still too high for practical use. This contribution presents a method that uses a validated CFD model with tracer experiments as a platform for the development of a simpler compartmental model (CM) to describe the hydraulics in a full-scale maturation pond (7 ha) of a waste stabilization ponds complex in Cuenca (Ecuador). 3D CFD models were validated with experimental data from pulse tracer experiments, showing a sufficient agreement. Based on the CFD model results, a number of compartments were selected considering the turbulence characteristics of the flow, the residence time distribution (RTD) curves and the dominant velocity component at different pond locations. The arrangement of compartments based on the introduction of recirculation flow rate between adjacent compartments, which in turn is dependent on the turbulence diffusion coefficient, is illustrated. Simulated RTD's from a systemic tanks-in-series (TIS) model and the developed CM were compared. The TIS was unable to capture the measured RTD, whereas the CM predicted convincingly the peaks and lags of the tracer experiment using only a minimal fraction of the computational demand of the CFD model. Finally, a biokinetic model was coupled to both approaches demonstrating the impact an insufficient hydraulic model can have on the outcome of a modelling exercise. TIS and CM showed drastic differences in the output loads implying that the CM approach is to be used when modelling the biological performance of the full-scale system. PMID:22137448
The effect of dynamic changes in soil bulk density on hydraulic properties: modeling approaches
NASA Astrophysics Data System (ADS)
Assouline, Shmuel
2014-05-01
Natural and artificial processes, like rainfall-induced soil surface sealing or mechanical compaction, disturb the soil structure and enhance dynamic changes of the related pore size distribution. These changes may influence many aspects of the soil-water-plant-atmosphere system. One of the easiest measurable variables is the soil bulk density. Approaches are suggested that could model the effect of the change in soil bulk density on soil permeability, water retention curve (WRC) and unsaturated hydraulic conductivity function (HCF). The resulting expressions were calibrated and validated against experimental data corresponding to different soil types at various levels of compaction, and enable a relatively good prediction of the effect of bulk density on the soil hydraulic properties. These models allow estimating the impact of such changes on flow processes and on transport properties of heterogeneous soil profiles.
Modeling the hydraulics of root growth in three dimensions with phloem water sources.
Wiegers, Brandy S; Cheer, Angela Y; Silk, Wendy K
2009-08-01
Primary growth is characterized by cell expansion facilitated by water uptake generating hydrostatic (turgor) pressure to inflate the cell, stretching the rigid cell walls. The multiple source theory of root growth hypothesizes that root growth involves transport of water both from the soil surrounding the growth zone and from the mature tissue higher in the root via phloem and protophloem. Here, protophloem water sources are used as boundary conditions in a classical, three-dimensional model of growth-sustaining water potentials in primary roots. The model predicts small radial gradients in water potential, with a significant longitudinal gradient. The results improve the agreement of theory with empirical studies for water potential in the primary growth zone of roots of maize (Zea mays). A sensitivity analysis quantifies the functional importance of apical phloem differentiation in permitting growth and reveals that the presence of phloem water sources makes the growth-sustaining water relations of the root relatively insensitive to changes in root radius and hydraulic conductivity. Adaptation to drought and other environmental stresses is predicted to involve more apical differentiation of phloem and/or higher phloem delivery rates to the growth zone. PMID:19542299
NASA Astrophysics Data System (ADS)
Liu, C. H.
2015-12-01
Atmospheric stability has substantial effects on the flows and heat/mass transport processes. While extensive studies have been conducted for neutral and unstable stabilities, rather limited studies have been devoted to stable stratification. Major technical reason is the demanding spatio-temporal resolution required to solve the small scales in stratified turbulent flows. Instead of continuous density variation, we use the single-layer hydraulics model (analogous to shallow water equations for global dynamics), to simulate the stratified flows and turbulence structure over hypothetical urban areas. An array of identical ribs in cross flows is used to model an idealized urban surface and the aerodynamic resistance is controlled by the separation among the ribs. Two immiscible fluids (water and air) with a large density difference (three order of magnitude) are used to simulate the stratification. The key assumption is that the density in the (lower) single layer is uniform. As a result, the stratification is measured by the Froude number Fr (= U/(gH)1/2; where U is the flow speed, g the gravitational acceleration and H the single-layer depth). One of the characteristics of single-layer hydraulics model is hydraulic jump which occurs when the flows are slowing down from Fr > 1 (high-speed flows over smoother surfaces) to Fr < 1 (lower-speed flows over rougher surfaces). It is noteworthy that kinetic energy does not conserve across hydraulic jump that, unavoidably, cascades to turbulent kinetic energy (TKE). We thus hypotheses that the elevated TKE could modify the street-level ventilation mechanism in the stratified flows across an abrupt change in surface roughness entering urban areas. Large-eddy simulation and laboratory-scale water channel experiments are sought to improve our understanding of the occurrence of hydraulic jump and the associated street-level ventilation mechanism in the stratified flows over urban areas. Preliminary results, by comparing the
TWIST: a transient two-dimensional intra-subassembly thermal hydraulics model for LMFBRs
Khatib-Rahbar, M.; Cazzoli, E.G.
1984-06-03
Mathematical models and numerical methods for a two-dimensional porous body simulation of steady state and transient thermal-hydraulics conditions in LMFBR subassemblies resulting in the TWIST computer code are presented. Comparison of calculated results to steady state and transient out-of-pile sodium experiments show good agreement for cross-assembly temperature distributions for a wide range of heat transfer and flow conditions.
NASA Astrophysics Data System (ADS)
Yoon, Jeoung Seok; Zang, Arno; Zimmermann, Günter; Stephansson, Ove; Min, Ki-Bok
2015-04-01
This paper presents discrete element based numerical model which is applied to simulation of multiple stage hydraulic fracturing in crystalline granitic geothermal reservoir. Target site modeled locates in south of state of Saxony Germany. Particle Flow Code 2D (Itasca) is used in which fluid flow algorithm and moment tensor based seismicity computation algorithm are implemented. Crystalline rock layer to be stimulated locates at 4-6 km depth with relative low density of pre-existing joints and faults. Hydraulic stimulation is modeled with five stages of fluid injection with distance of several hundreds of meters. Hydraulic fracturing is done on the stages from toe to heel direction along a series of sub-horizontally drilled wellbore with constant rate of fluid injection. Fracture propagation paths and induced seismic events are documented based on their time of occurrence and their magnitude. In addition to the evolution of the fracture propagation path and distribution of the induced events, migration of the injected fluid is investigated in space and time. This is to see how the results relate to the fluid migration front in low permeability crystalline reservoir subjected to multiple stage hydraulic fracturing. Moreover, this paper addresses advantages and disadvantages of the inclined drilling of the wellbore in low permeability reservoir and multi-stage fracturing setting. We try to seek an optimum inclination of the drilling in relation to the gradients and magnitudes of the in situ stresses, which are horizontal minimum and vertical stresses. Preliminary modeling results show that inclination angle of the drilling has a significant effect on lowering of the stress shadow effect and level of induced seismicity in terms of total number, magnitudes and the Gutenberg-Richter relation.
NASA Technical Reports Server (NTRS)
Bartos, R. D.
1992-01-01
As the pointing accuracy and service life requirements of the DSN 70 meter antenna increase, it is necessary to gain a more complete understanding of the servo hydraulic system in order to improve system designs to meet the new requirements. A mathematical model is developed for the servovalve incorporated into the hydraulic system of the 70 meter antenna and uses experimental data to verify the validity of the model and to identify the model parameters.
Hydraulic-fracture propagation in layered rock: experimental studies of fracture containment
Teufel, L. W.; Clark, J. A.
1981-01-01
Fracture geometry is an important concern in the design of a massive hydraulic fracture treatment for improved natural gas recovery from tight gas sands. Possible prediction of vertical fracture growth and containment in layered rock requires an improved understanding of the parameters which may control fracture growth across layer interfaces. We have conducted laboratory hydraulic fracture experiments and elastic finite element studies which show that at least two distinct geologic conditions may inhibit or contain the vertical growth of hydraulic fractures in layered rock; (1) a weak interfacial shear strength of the layers and (2) a compressional increase in the minimum horizontal stress in the bounding layer. The second condition is more important and more likely to occur at depth. Variations in the horizontal stress can result from differences in elastic properties of individual layers in a layered rock sequence. A compressional increase in the minimum horizontal stress can occur in going from high shear modulus into low shear modulus layers.
Modelling cave flow hydraulics in the Notranjski Kras, Slovenia
NASA Astrophysics Data System (ADS)
Kaufmann, Georg; Gabrovsek, Franci
2015-04-01
The Notranjski Kras region is a karst region in western Slovenia, developed in Cretaceous limestone. The region is characterised by hilly relief, with peaks reaching 1300 m elevation. Several well-developed cave systems drain the karst aquifer, providing preferential flow pathes along two sections: The Pivka River, which sinks into Postojnska Jama and reappears in Planinska Jama, and the Stržen and Cerkniščica rivers, which sink into Karlovica Jama, flow through Zelške Jama and Tkalca Jama and also reappear in Planinska Jama. Both sub-surface flow pathes merge in Planinska Jama, providing water for the Unica river. The Unica river leaves Planinska Jama via a large karst srping and passes through Planinsko Polje, disappearing again through two groups of ponors, finally emerging in the Ljubljanka Springs at around 300 m asl. The sub-surface flow path through the Postojnska Jama cave system has been monitored with 7 stations distributed along the flow path, monitoring stage and temperature. We have used the stage data to model flow through the cave system with the program package SWMM, simulating the active parts of Postojnska Jama with simplified geometry. From the comparison of stage observations and predictions, we identified key sections in the cave, which control the sub-surface flow, such as passage constrictions, sumps and by-passes. Using a formal inverse procedure, we determined the geometry of this key sections by fitting predicted to observed stages, and we achieved a very high degree of correlation.
Drift trajectories of a floating human body simulated in a hydraulic model of Puget Sound.
Ebbesmeyer, C C; Haglund, W D
1994-01-01
After a young man jumped off a 221-foot (67 meters) high bridge, the drift of the body that beached 20 miles (32 km) away at Alki Point in Seattle, Washington was simulated with a hydraulic model. Simulations for the appropriate time period were performed using a small floating bead to represent the body in the hydraulic model at the University of Washington. Bead movements were videotaped and transferred to Computer Aided Drafting (AutoCAD) charts on a personal computer. Because of strong tidal currents in the narrow passage under the bridge (The Narrows near Tacoma, WA), small changes in the time of the jump (+/- 30 minutes) made large differences in the distance the body traveled (30 miles; 48 km). Hydraulic and other types of oceanographic models may be located by contacting technical experts known as physical oceanographers at local universities, and can be utilized to demonstrate trajectories of floating objects and the time required to arrive at selected locations. Potential applications for forensic death investigators include: to be able to set geographic and time limits for searches; determine potential origin of remains found floating or beached; and confirm and correlate information regarding entry into the water and sightings of remains. PMID:8113703