Impact of heterogeneity on groundwater salinization due to coastal pumping

NASA Astrophysics Data System (ADS)

Yu, X.; Michael, H. A.

2017-12-01

Groundwater abstraction causes and accelerates seawater intrusion in many coastal areas. In heterogeneous aquifers, preferential flow paths can lead to fast intrusion, while low permeability layers can serve as barriers. The extent to which different types of heterogeneous aquifers are vulnerable to pumping-induced seawater intrusion has not been well studied. Here we show that the connectedness of pumping location and local boundary condition drive salinization patterns. Salinization patterns in homogeneous aquifers were relatively simple and only related to the hydraulic properties and pumping rate. The salinization rates and patterns in heterogeneous aquifers were much more complicated and related to pumping location, rate and depth, preferential flow path locations, and local boundary conditions. An intrusion classification approach was developed with three types in homogeneous aquifers and four types in heterogeneous aquifers. After classification the main factors of salinized areas, intrusion rates and salinization time could be identified. The ranges of these salinization assessment criteria suggested different aspect of groundwater vulnerability in each class. We anticipate the classification approach to be a starting point for more comprehensive groundwater abstraction vulnerability assessment (including consideration of pumping rates, locations and depths, connectivity, preferential flow paths, etc.), which is critical for coastal water resources management.

Flume experimental evaluation of the effect of rill flow path tortuosity on rill roughness based on the Manning–Strickler equation

USDA-ARS?s Scientific Manuscript database

Numerous soil erosion models compute concentrated flow hydraulics based on the Manning–Strickler equation (v = kSt R2/3 I1/2) even though the range of the application on rill flow is unclear. Unconfined rill morphologies generate local friction effects and consequently spatially variable rill roughn...

Path selection in the growth of rivers

DOE PAGES

Cohen, Yossi; Devauchelle, Olivier; Seybold, Hansjörg F.; ...

2015-11-02

River networks exhibit a complex ramified structure that has inspired decades of studies. But, an understanding of the propagation of a single stream remains elusive. In this paper, we invoke a criterion for path selection from fracture mechanics and apply it to the growth of streams in a diffusion field. We show that, as it cuts through the landscape, a stream maintains a symmetric groundwater flow around its tip. The local flow conditions therefore determine the growth of the drainage network. We use this principle to reconstruct the history of a network and to find a growth law associated withmore » it. Finally, our results show that the deterministic growth of a single channel based on its local environment can be used to characterize the structure of river networks.« less

Quasi-2D Unsteady Flow Procedure for Real Fluids (PREPRINT)

DTIC Science & Technology

2006-05-17

water /steam/ oil piping networks, refinery systems, gas-turbine secondary flow -path and cooling networks...friction factor, f, which is a function of the local Reynolds number and the wall surface roughness . For the viscous flow examples presented below, the...3.5 4 4.5 Time ( s ) V el oc ity (m / s ) Line 2 Inlet 25% 50% 75% Exit Velocity Figure 4. Water transient viscous pipe flow using

Generating Performance Models for Irregular Applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Friese, Ryan D.; Tallent, Nathan R.; Vishnu, Abhinav

2017-05-30

Many applications have irregular behavior --- non-uniform input data, input-dependent solvers, irregular memory accesses, unbiased branches --- that cannot be captured using today's automated performance modeling techniques. We describe new hierarchical critical path analyses for the \\Palm model generation tool. To create a model's structure, we capture tasks along representative MPI critical paths. We create a histogram of critical tasks with parameterized task arguments and instance counts. To model each task, we identify hot instruction-level sub-paths and model each sub-path based on data flow, instruction scheduling, and data locality. We describe application models that generate accurate predictions for strong scalingmore » when varying CPU speed, cache speed, memory speed, and architecture. We present results for the Sweep3D neutron transport benchmark; Page Rank on multiple graphs; Support Vector Machine with pruning; and PFLOTRAN's reactive flow/transport solver with domain-induced load imbalance.« less

Detection of small degree of nonuniformity in dialysate flow in hollow-fiber dialyzer using proton magnetic resonance imaging.

PubMed

Osuga, T; Obata, T; Ikehira, H

2004-04-01

A small degree of nonuniformity in dialysate flow in a hollow-fiber dialyzer was detected using proton magnetic resonance imaging (MRI). Since paramagnetic ions reduce the spin-lattice relaxation time of protons around them, MRI can detect Gd in water. An aqueous solution of a chelate compound of Gd was impulsively injected into the dialysate flow path at a flow rate of 500 cm(3) /m, which is that utilized in actual dialysis. Despite the apparent elimination of Gd from the dialysate flow path by the newly injected dialysate fluid after the injection of Gd was terminated, MRI revealed that Gd remained in the interior of the hollow fiber. The observed structure pattern of the Gd concentration profile revealed that the dialysate flow had a small degree of nonuniformity despite the currently established design to restrict channeling in dialysate flow. Local nonuniformity of the hollow-fiber density and vortex generation in the dialysate flow were considered to cause the nonuniformity in the dialysate flow.

Simulation of ground-water flow in the Cedar River alluvium, northwest Black Hawk County and southwest Bremer County, Iowa

USGS Publications Warehouse

Schaap, Bryan D.; Savoca, Mark E.; Turco, Michael J.

2003-01-01

In general, once high ground-water levels occur, either because of high Cedar River water Abstract levels or above normal local precipitation or both, ground-water in the central part of the study area along Highway 218 flows toward the south rather than following shorter flow paths to the Cedar River. Intermittent streams in the study area discharge substantial amounts of water from the ground-water flow system.

Coupling 3D groundwater modeling with CFC-based age dating to classify local groundwater circulation in an unconfined crystalline aquifer

NASA Astrophysics Data System (ADS)

Kolbe, Tamara; Marçais, Jean; Thomas, Zahra; Abbott, Benjamin W.; de Dreuzy, Jean-Raynald; Rousseau-Gueutin, Pauline; Aquilina, Luc; Labasque, Thierry; Pinay, Gilles

2016-12-01

Nitrogen pollution of freshwater and estuarine environments is one of the most urgent environmental crises. Shallow aquifers with predominantly local flow circulation are particularly vulnerable to agricultural contaminants. Water transit time and flow path are key controls on catchment nitrogen retention and removal capacity, but the relative importance of hydrogeological and topographical factors in determining these parameters is still uncertain. We used groundwater dating and numerical modeling techniques to assess transit time and flow path in an unconfined aquifer in Brittany, France. The 35.5 km2 study catchment has a crystalline basement underneath a ∼60 m thick weathered and fractured layer, and is separated into a distinct upland and lowland area by an 80 m-high butte. We used groundwater discharge and groundwater ages derived from chlorofluorocarbon (CFC) concentration to calibrate a free-surface flow model simulating groundwater flow circulation. We found that groundwater flow was highly local (mean travel distance = 350 m), substantially smaller than the typical distance between neighboring streams (∼1 km), while CFC-based ages were quite old (mean = 40 years). Sensitivity analysis revealed that groundwater travel distances were not sensitive to geological parameters (i.e. arrangement of geological layers and permeability profile) within the constraints of the CFC age data. However, circulation was sensitive to topography in the lowland area where the water table was near the land surface, and to recharge rate in the upland area where water input modulated the free surface of the aquifer. We quantified these differences with a local groundwater ratio (rGW-LOCAL), defined as the mean groundwater travel distance divided by the mean of the reference surface distances (the distance water would have to travel across the surface of the digital elevation model). Lowland, rGW-LOCAL was near 1, indicating primarily topographical controls. Upland, rGW-LOCAL was 1.6, meaning the groundwater recharge area is almost twice as large as the topographically-defined catchment for any given point. The ratio rGW-LOCAL is sensitive to recharge conditions as well as topography and it could be used to compare controls on groundwater circulation within or between catchments.

Constraining Path-Dependent Processes During Basalt-CO2 Interactions with Observations From Flow-Through and Batch Experiments

NASA Astrophysics Data System (ADS)

Thomas, D.; Garing, C.; Zahasky, C.; Harrison, A. L.; Bird, D. K.; Benson, S. M.; Oelkers, E. H.; Maher, K.

2017-12-01

Predicting the timing and magnitude of CO2 storage in basaltic rocks relies partly on quantifying the dependence of reactivity on flow path and mineral distribution. Flow-through experiments that use intact cores are advantageous because the spatial heterogeneity of pore space and reactive phases is preserved. Combining aqueous geochemical analyses and petrologic characterization with non-destructive imaging techniques (e.g. micro-computed tomography) constrains the relationship between irreversible reactions, pore connectivity and accessible surface area. Our work enhances these capabilities by dynamically imaging flow through vesicular basalts with Positron Emission Tomography (PET) scanning. PET highlights the path a fluid takes by detecting photons produced during radioactive decay of an injected radiotracer (FDG). We have performed single-phase, CO2-saturated flow-through experiments with basaltic core from Iceland at CO2 sequestration conditions (50 °C; 76-90 bar Ptot). Constant flow rate and continuous pressure measurements at the inlet and outlet of the core constrain permeability. We monitor geochemical evolution through cation and anion analysis of outlet fluid sampled periodically. Before and after reaction, we perform PET scans and characterize the core using micro-CT. The PET scans indicate a discrete, localized flow path that appears to be a micro-crack connecting vesicles, suggesting that vesicle-lining minerals are immediately accessible and important reactants. Rapid increases in aqueous cation concentration, pH and HCO3- indicate that the rock reacts nearly immediately after CO2 injection. After 24 hours the solute release decreases, which may reflect a transition to reaction with phases with slower kinetic dissolution rates (e.g. zeolites and glasses to feldspar), a decrease in available reactive surface area or precipitation. We have performed batch experiments using crushed material of the same rock to elucidate the effect of flow path geometry and mineral accessibility on geochemical evolution. Interestingly, surface area-normalized dissolution rates as evinced by SiO2 release in all experiments approach similar values ( 10-15 mol/cm2/s). Our experiments show how imaging techniques are helpful in interpreting path-dependent processes in open systems.

Geochemical modeling of iron, sulfur, oxygen and carbon in a coastal plain aquifer

USGS Publications Warehouse

Brown, C.J.; Schoonen, M.A.A.; Candela, J.L.

2000-01-01

Fe(III) reduction in the Magothy aquifer of Long Island, NY, results in high dissolved-iron concentrations that degrade water quality. Geochemical modeling was used to constrain iron-related geochemical processes and redox zonation along a flow path. The observed increase in dissolved inorganic carbon is consistent with the oxidation of sedimentary organic matter coupled to the reduction of O2 and SO4/2- in the aerobic zone, and to the reduction of SO4/2- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55 x 10-4 to 48.6 x 10-4 mmol 1-1 yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO4/2- (1.31 x 10-4 to 15 x 10-4 mmol 1-1 yr-1). The overall increase in SO4/2- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO4/2- loss through microbial reduction is exceeded by SO4/2- gain through diffusion from sediments and through the oxidation of FeS2. Geochemichal and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO4/2- -rich zones have localized SO4/2- -reducing zones in which the formation of iron disulfides been depleted by microbial reduction and resulted in decreases dissolved iron concentrations. These localized zones of SO4/2- reduction, which are important for assessing zones of low dissolved iron for water-supply development, could be overlooked by aquifer studies that rely only on groundwater data from well-water samples for geochemical modeling. (C) 2000 Elsevier Science B.V.Fe(III) reduction in the Magothy aquifer of Long Island, NY, results in high dissolved-iron concentrations that degrade water quality. Geochemical modeling was used to constrain iron-related geochemical processes and redox zonation along a flow path. The observed increase in dissolved inorganic carbon is consistent with the oxidation of sedimentary organic matter coupled to the reduction of O2 and SO42- in the aerobic zone, and to the reduction of SO42- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55??10-4 to 48.6??10-4mmol l-1yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO42- (1.31??10-4 to 15??10-4mmol l-1yr-1). The overall increase in SO42- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO42- loss through microbial reduction is exceeded by SO42- gain through diffusion from sediments and through the oxidation of FeS2. Geochemical and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO42--rich zones have been depleted by microbial reduction and resulted in localized SO42--reducing zones in which the formation of iron disulfides decreases dissolved iron concentrations. These localized zones of SO42- reduction, which are important for assessing zones of low dissolved iron for water-supply development, could be overlooked by aquifer studies that rely only on groundwater data from well-water samples for geochemical modeling.

Acoustic concentration of particles in fluid flow

DOEpatents

Ward, Michael D.; Kaduchak, Gregory

2010-11-23

An apparatus for acoustic concentration of particles in a fluid flow includes a substantially acoustically transparent membrane and a vibration generator that define a fluid flow path therebetween. The fluid flow path is in fluid communication with a fluid source and a fluid outlet and the vibration generator is disposed adjacent the fluid flow path and is capable of producing an acoustic field in the fluid flow path. The acoustic field produces at least one pressure minima in the fluid flow path at a predetermined location within the fluid flow path and forces predetermined particles in the fluid flow path to the at least one pressure minima.

A unified gas-kinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Chang, E-mail: cliuaa@ust.hk; Xu, Kun, E-mail: makxu@ust.hk; Sun, Quanhua, E-mail: qsun@imech.ac.cn

Fluid dynamic equations are valid in their respective modeling scales, such as the particle mean free path scale of the Boltzmann equation and the hydrodynamic scale of the Navier–Stokes (NS) equations. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. Even though the Boltzmann equation is claimed to be valid in all scales, many Boltzmann solvers, including direct simulation Monte Carlo method, require the cell resolution to the order of particle mean free path scale. Therefore, they are still single scale methods. In order to study multiscale flow evolution efficiently, themore » dynamics in the computational fluid has to be changed with the scales. A direct modeling of flow physics with a changeable scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS) is a direct modeling method in the mesh size scale, and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path. The cell size of UGKS is not limited by the particle mean free path. With the variation of the ratio between the numerical cell size and local particle mean free path, the UGKS recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale. The previous UGKS is mostly constructed from the evolution solution of kinetic model equations. Even though the UGKS is very accurate and effective in the low transition and continuum flow regimes with the time step being much larger than the particle mean free time, it still has space to develop more accurate flow solver in the region, where the time step is comparable with the local particle mean free time. In such a scale, there is dynamic difference from the full Boltzmann collision term and the model equations. This work is about the further development of the UGKS with the implementation of the full Boltzmann collision term in the region where it is needed. The central ingredient of the UGKS is the coupled treatment of particle transport and collision in the flux evaluation across a cell interface, where a continuous flow dynamics from kinetic to hydrodynamic scales is modeled. The newly developed UGKS has the asymptotic preserving (AP) property of recovering the NS solutions in the continuum flow regime, and the full Boltzmann solution in the rarefied regime. In the mostly unexplored transition regime, the UGKS itself provides a valuable tool for the non-equilibrium flow study. The mathematical properties of the scheme, such as stability, accuracy, and the asymptotic preserving, will be analyzed in this paper as well.« less

Ray-based acoustic localization of cavitation in a highly reverberant environment.

PubMed

Chang, Natasha A; Dowling, David R

2009-05-01

Acoustic detection and localization of cavitation have inherent advantages over optical techniques because cavitation bubbles are natural sound sources, and acoustic transduction of cavitation sounds does not require optical access to the region of cavitating flow. In particular, near cavitation inception, cavitation bubbles may be visually small and occur infrequently, but may still emit audible sound pulses. In this investigation, direct-path acoustic recordings of cavitation events are made with 16 hydrophones mounted on the periphery of a water tunnel test section containing a low-cavitation-event-rate vortical flow. These recordings are used to localize the events in three dimensions via cross correlations to obtain arrival time differences. Here, bubble localization is hindered by reverberation, background noise, and the fact that both the pulse emission time and waveform are unknown. These hindrances are partially mitigated by a signal-processing scheme that incorporates straight-ray acoustic propagation and Monte-Carlo techniques for compensating ray-path, sound-speed, and hydrophone-location uncertainties. The acoustic localization results are compared to simultaneous optical localization results from dual-camera high-speed digital-video recordings. For 53 bubbles and a peak-signal to noise ratio frequency of 6.7 kHz, the root-mean-square spatial difference between optical and acoustic bubble location results was 1.94 cm. Parametric dependences in acoustic localization performance are also presented.

Acoustic concentration of particles in fluid flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ward, Michael W.; Kaduchak, Gregory

Disclosed herein is a acoustic concentration of particles in a fluid flow that includes a substantially acoustically transparent membrane and a vibration generator that define a fluid flow path therebetween. The fluid flow path is in fluid communication with a fluid source and a fluid outlet and the vibration generator is disposed adjacent the fluid flow path and is capable of producing an acoustic field in the fluid flow path. The acoustic field produces at least one pressure minima in the fluid flow path at a predetermined location within the fluid flow path and forces predetermined particles in the fluidmore » flow path to the at least one pressure minima.« less

Supersonic compressor

DOEpatents

Lawlor, Shawn P [Bellevue, WA; Novaresi, Mark A [San Diego, CA; Cornelius, Charles C [Kirkland, WA

2008-02-26

A gas compressor based on the use of a driven rotor having an axially oriented compression ramp traveling at a local supersonic inlet velocity (based on the combination of inlet gas velocity and tangential speed of the ramp) which forms a supersonic shockwave axially, between adjacent strakes. In using this method to compress inlet gas, the supersonic compressor efficiently achieves high compression ratios while utilizing a compact, stabilized gasdynamic flow path. Operated at supersonic speeds, the inlet stabilizes an oblique/normal shock system in the gasdyanamic flow path formed between the gas compression ramp on a strake, the shock capture lip on the adjacent strake, and captures the resultant pressure within the stationary external housing while providing a diffuser downstream of the compression ramp.

Modeling of ground-water flow in subsurface Austin Chalk and Taylor marl in Ellis County, Texas, near the superconducting super collider site

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mace, R.E.

1993-02-01

Numerical models are useful tools for developing an understanding of ground-water flow in sparsely characterized low-permeability aquifers. Finite-difference, cross-sectional models of Cretaceous chalk and marl formations near the Superconducting Super Collider (SSC) were constructed using MODFLOW to evaluate ground-water circulation paths and travel times. Weathered and fractured zones with enhanced permeability were included to assess the effect these features had on flow paths and times. Pump tests, slug tests, packer tests, core tests, and estimates were used to define hydraulic properties for model input. The model was calibrated with water-level data from monitor wells and from wire-line piezometers near amore » test shaft excavated by the SSC project. A ratio of vertical-to-horizontal permeability of 0.0085 was estimated through model calibration. A chalk-to-marl permeability ratio of 18 was needed to reproduce artesian head in a well completed in chalk beneath marl. Hydraulic head distributions and ground-water flow paths reflected local, intermediate, and regional flow systems with recharge beneath upland surface-water divides and discharge in valleys. Most of the flow (99%) occurred in the weathered zone, with average residence times of 5 to 10 years. Residence time in unweathered chalk bedrock was substantially longer, at an average of 1.7 Ma. As expected, the model demonstrated that deep and rapid ground-water circulation might occur in fracture zones. Particle paths calculated using MODPATH showed that ground-water travel times from recharge areas to the SSC subsurface facilities might be 20 to 60 years where flow is through fracture zones.« less

Augmenting two-dimensional hydrodynamic simulations with measured velocity data to identify flow paths as a function of depth on Upper St. Clair River in the Great Lakes basin

USGS Publications Warehouse

Holtschlag, D.J.; Koschik, J.A.

2005-01-01

Upper St. Clair River, which receives outflow from Lake Huron, is characterized by flow velocities that exceed 7 feet per second and significant channel curvature that creates complex flow patterns downstream from the Blue Water Bridge in the Port Huron, Michigan, and Sarnia, Ontario, area. Discrepancies were detected between depth-averaged velocities previously simulated by a two-dimensional (2D) hydrodynamic model and surface velocities determined from drifting buoy deployments. A detailed ADCP (acoustic Doppler current profiler) survey was done on Upper St. Clair River during July 1–3, 2003, to help resolve these discrepancies. As part of this study, a refined finite-element mesh of the hydrodynamic model used to identify source areas to public water intakes was developed for Upper St. Clair River. In addition, a numerical procedure was used to account for radial accelerations, which cause secondary flow patterns near channel bends. The refined model was recalibrated to better reproduce local velocities measured in the ADCP survey. ADCP data also were used to help resolve the remaining discrepancies between simulated and measured velocities and to describe variations in velocity with depth. Velocity data from ADCP surveys have significant local variability, and statistical processing is needed to compute reliable point estimates. In this study, velocity innovations were computed for seven depth layers posited within the river as the differences between measured and simulated velocities. For each layer, the spatial correlation of velocity innovations was characterized by use of variogram analysis. Results were used with kriging to compute expected innovations within each layer at applicable model nodes. Expected innovations were added to simulated velocities to form integrated velocities, which were used with reverse particle tracking to identify the expected flow path near a sewage outfall as a function of flow depth. Expected particle paths generated by use of the integrated velocities showed that surface velocities in the upper layers tended to originate nearer the Canadian shoreline than velocities near the channel bottom in the lower layers. Therefore, flow paths to U.S. public water intakes located on the river bottom are more likely to be in the United States than withdrawals near the water surface. Integrated velocities in the upper layers are generally consistent with the surface velocities indicated by drifting-buoy deployments. Information in the 2D hydrodynamic model and the ADCP measurements was insufficient to describe the vertical flow component. This limitation resulted in the inability to account for vertical movements on expected flow paths through Upper St. Clair River. A three dimensional hydrodynamic model would be needed to account for these effects.

Thermal transient anemometer

DOEpatents

Bailey, James L.; Vresk, Josip

1989-01-01

A thermal transient anemometer having a thermocouple probe which is utilized to measure the change in temperature over a period of time to provide a measure of fluid flow velocity. The thermocouple probe is located in the fluid flow path and pulsed to heat or cool the probe. The cooling of the heated probe or the heating of the cooled probe from the fluid flow over a period of time is measured to determine the fluid flow velocity. The probe is desired to be locally heated near the tip to increase the efficiency of devices incorporating the probe.

Mitigation of hazards from future lahars from Mount Merapi in the Krasak River channel near Yogyakarta, central Java

USGS Publications Warehouse

Ege, John R.; ,

1983-01-01

Procedures for reducing hazards from future lahars and debris flows in the Krasak River channel near Yogyakarta, Central Java, Indonesia, include (1) determining the history of the location, size, and effects of previous lahars and debris flows, and (2) decreasing flow velocities. The first may be accomplished by geologic field mapping along with acquiring information by interviewing local residents, and the second by increasing the cross sectional area of the river channel and constructing barriers in the flow path.

An implicit scheme with memory reduction technique for steady state solutions of DVBE in all flow regimes

NASA Astrophysics Data System (ADS)

Yang, L. M.; Shu, C.; Yang, W. M.; Wu, J.

2018-04-01

High consumption of memory and computational effort is the major barrier to prevent the widespread use of the discrete velocity method (DVM) in the simulation of flows in all flow regimes. To overcome this drawback, an implicit DVM with a memory reduction technique for solving a steady discrete velocity Boltzmann equation (DVBE) is presented in this work. In the method, the distribution functions in the whole discrete velocity space do not need to be stored, and they are calculated from the macroscopic flow variables. As a result, its memory requirement is in the same order as the conventional Euler/Navier-Stokes solver. In the meantime, it is more efficient than the explicit DVM for the simulation of various flows. To make the method efficient for solving flow problems in all flow regimes, a prediction step is introduced to estimate the local equilibrium state of the DVBE. In the prediction step, the distribution function at the cell interface is calculated by the local solution of DVBE. For the flow simulation, when the cell size is less than the mean free path, the prediction step has almost no effect on the solution. However, when the cell size is much larger than the mean free path, the prediction step dominates the solution so as to provide reasonable results in such a flow regime. In addition, to further improve the computational efficiency of the developed scheme in the continuum flow regime, the implicit technique is also introduced into the prediction step. Numerical results showed that the proposed implicit scheme can provide reasonable results in all flow regimes and increase significantly the computational efficiency in the continuum flow regime as compared with the existing DVM solvers.

Topographic Effects on Geologic Mass Movements

NASA Technical Reports Server (NTRS)

Baloga, Stephen M.; Frey, Herbert (Technical Monitor)

2000-01-01

This report describes research directed toward understanding the response of volcanic lahars and lava flows to changes in the topography along the path of the flow. We have used a variety of steady-state and time-dependent models of lahars and lava flows to calculate the changes in flow dynamics due to variable topography. These models are based on first-order partial differential equations for the local conservation of volume. A global volume conservation requirement is also imposed to determine the extent of the flow as a function of time and the advance rate. Simulated DEMs have been used in this report.

Postphloem, Nonvascular Transfer in Citrus

PubMed Central

Koch, Karen E.; Avigne, Wayne T.

1990-01-01

Postphloem, nonvascular assimilate transport occurs over an unusually long area in citrus fruit and thus facilitates investigation of this process relative to sugar entry into many sink structures. Labeled photosynthates moving into juice tissues of grapefruit (Citrus paradisi Macf.) slowed dramatically after entering the postphloem transport path (parenchyma cells, narrow portions of segment epidermis, and hair-like, parenchymatous stalks of juice sacs). Kinetic, metabolic, and compositional data indicated that transfer through the nonvascular area was delayed many hours by temporary storage and/or equilibration with sugars in compartments along the postphloem path. Labeled assimilates were generally recovered as sucrose throughout the path, and extent of hexose formation enroute bore no apparent relationship to the assimilate transfer process. Even after 24 hours, radiolabel was restricted to discrete, highly localized areas directly between vascular bundles and juice sacs. Postphloem transfer occurred against an ascending sucrose concentration gradient in young fruit, whereas a descending gradient (favoring diffusion/cytoplasmic streaming) developed only later in maturation. Involvement of a postphloem bulk flow is complicated in the present instance by the extremely limited water loss from juice sacs either via transpiration or fluid backflow. Nonetheless, tissue expansion can account for a collective water inflow of at least 1.0 milliliter per day throughout the majority of juice sac development, thus providing a modest, but potentially important means of nonvascular solution flow. Overall, data indicate postphloem transfer (a) can follow highly localized paths through sizable nonvascular areas (up to 3.0 centimeters total), (b) appears to involve temporary storage and/or equilibration with compartmentalized sugars enroute, (c) can occur either against an overall up-hill sugar gradient (young tissues) or along a descending gradient (near full expansion), and (d) appears to involve at least some contribution by nonvascular mass flow accommodated by tissue expansion. Images Figure 1 Figure 4 PMID:16667632

Flow paths in the Edwards aquifer, northern Medina and northeastern Uvalde counties, Texas, based on hydrologic identification and geochemical characterization and simulation

USGS Publications Warehouse

Clark, Allan K.; Journey, Celeste A.

2006-01-01

The U.S. Geological Survey, in cooperation with the San Antonio Water System, conducted a 4-year study during 2001– 04 to identify major ground-water flow paths in the Edwards aquifer in northern Medina and northeastern Uvalde Counties, Texas. The study involved use of geologic structure, surfacewater and ground-water data, and geochemistry to identify ground-water flow paths. Relay ramps and associated faulting in northern Medina County appear to channel ground-water flow along four distinct flow paths that move water toward the southwest. The northwestern Medina flow path is bounded on the north by the Woodard Cave fault and on the south by the Parkers Creek fault. Water moves downdip toward the southwest until the flow encounters a cross fault along Seco Creek. This barrier to flow might force part or most of the flow to the south. Departure hydrographs for two wells and discharge departure for a streamflow-gaging station provide evidence for flow in the northwestern Medina flow path. The north-central Medina flow path (northern part) is bounded by the Parkers Creek fault on the north and the Medina Lake fault on the south. The adjacent north-central Medina flow path (southern part) is bounded on the north by the Medina Lake fault and on the south by the Diversion Lake fault. The north-central Medina flow path is separated into a northern and southern part because of water-level differences. Ground water in both parts of the northcentral Medina flow path moves downgradient (and down relay ramp) from eastern Medina County toward the southwest. The north-central Medina flow path is hypothesized to turn south in the vicinity of Seco Creek as it begins to be influenced by structural features. Departure hydrographs for four wells and Medina Lake and discharge departure for a streamflow-gaging station provide evidence for flow in the north-central Medina flow path. The south-central Medina flow path is bounded on the north by the Seco Creek and Diversion Lake faults and on the south by the Haby Crossing fault. Because of bounding faults oriented northeast-southwest and adjacent flow paths directed south by other geologic structures, the south-central Medina flow path follows the configuration of the adjacent flow paths—oriented initially southwest and then south. Immediately after turning south, the south-central Medina flow path turns sharply east. Departure hydrographs for four wells and discharge departure for a streamflow-gaging station provide evidence for flow in the south-central Medina flow path. Statistical correlations between water-level departures for 11 continuously monitored wells provide additional evidence for the hypothesized flow paths. Of the 55 combinations of departure dataset pairs, the stronger correlations (those greater than .6) are all among wells in the same flow path, with one exception. Simulations of compositional differences in water chemistry along a hypothesized flow path in the Edwards aquifer and between ground-water and surface-water systems near Medina Lake were developed using the geochemical model PHREEQC. Ground-water chemistry for samples from five wells in the Edwards aquifer in the northwestern Medina flow path were used to evaluate the evolution of ground-water chemistry in the northwestern Medina flow path. Seven simulations were done for samples from pairs of these wells collected during 2001–03; three of the seven yielded plausible models. Ground-water samples from 13 wells were used to evaluate the evolution of ground-water chemistry in the north-central Medina flow path (northern and southern parts). Five of the wells in the most upgradient part of the flow path were completed in the Trinity aquifer; the remaining eight were completed in the Edwards aquifer. Nineteen simulations were done for samples from well pairs collected during 1995–2003; eight of the 19 yielded plausible models. Ground-water samples from seven wells were used to evaluate the evolution of ground-water chemistry in the south-central Medina flow path. One well was the Trinity aquifer end-member well upgradient from all flow paths, and another was a Trinity aquifer well in the most upgradient part of the flow path; all other wells were completed in the Edwards aquifer. Nine simulations were done for samples from well pairs collected during 1996–2003; seven of the nine yielded plausible models. The plausible models demonstrate that the four hypothesized flow paths can be partially supported geochemically. 

Fast interactive exploration of 4D MRI flow data

NASA Astrophysics Data System (ADS)

Hennemuth, A.; Friman, O.; Schumann, C.; Bock, J.; Drexl, J.; Huellebrand, M.; Markl, M.; Peitgen, H.-O.

2011-03-01

1- or 2-directional MRI blood flow mapping sequences are an integral part of standard MR protocols for diagnosis and therapy control in heart diseases. Recent progress in rapid MRI has made it possible to acquire volumetric, 3-directional cine images in reasonable scan time. In addition to flow and velocity measurements relative to arbitrarily oriented image planes, the analysis of 3-dimensional trajectories enables the visualization of flow patterns, local features of flow trajectories or possible paths into specific regions. The anatomical and functional information allows for advanced hemodynamic analysis in different application areas like stroke risk assessment, congenital and acquired heart disease, aneurysms or abdominal collaterals and cranial blood flow. The complexity of the 4D MRI flow datasets and the flow related image analysis tasks makes the development of fast comprehensive data exploration software for advanced flow analysis a challenging task. Most existing tools address only individual aspects of the analysis pipeline such as pre-processing, quantification or visualization, or are difficult to use for clinicians. The goal of the presented work is to provide a software solution that supports the whole image analysis pipeline and enables data exploration with fast intuitive interaction and visualization methods. The implemented methods facilitate the segmentation and inspection of different vascular systems. Arbitrary 2- or 3-dimensional regions for quantitative analysis and particle tracing can be defined interactively. Synchronized views of animated 3D path lines, 2D velocity or flow overlays and flow curves offer a detailed insight into local hemodynamics. The application of the analysis pipeline is shown for 6 cases from clinical practice, illustrating the usefulness for different clinical questions. Initial user tests show that the software is intuitive to learn and even inexperienced users achieve good results within reasonable processing times.

Estimation of distribution algorithm with path relinking for the blocking flow-shop scheduling problem

NASA Astrophysics Data System (ADS)

Shao, Zhongshi; Pi, Dechang; Shao, Weishi

2018-05-01

This article presents an effective estimation of distribution algorithm, named P-EDA, to solve the blocking flow-shop scheduling problem (BFSP) with the makespan criterion. In the P-EDA, a Nawaz-Enscore-Ham (NEH)-based heuristic and the random method are combined to generate the initial population. Based on several superior individuals provided by a modified linear rank selection, a probabilistic model is constructed to describe the probabilistic distribution of the promising solution space. The path relinking technique is incorporated into EDA to avoid blindness of the search and improve the convergence property. A modified referenced local search is designed to enhance the local exploitation. Moreover, a diversity-maintaining scheme is introduced into EDA to avoid deterioration of the population. Finally, the parameters of the proposed P-EDA are calibrated using a design of experiments approach. Simulation results and comparisons with some well-performing algorithms demonstrate the effectiveness of the P-EDA for solving BFSP.

Simulation and validation of concentrated subsurface lateral flow paths in an agricultural landscape

NASA Astrophysics Data System (ADS)

Zhu, Q.; Lin, H. S.

2009-08-01

The importance of soil water flow paths to the transport of nutrients and contaminants has long been recognized. However, effective means of detecting concentrated subsurface flow paths in a large landscape are still lacking. The flow direction and accumulation algorithm based on single-direction flow algorithm (D8) in GIS hydrologic modeling is a cost-effective way to simulate potential concentrated flow paths over a large area once relevant data are collected. This study tested the D8 algorithm for simulating concentrated lateral flow paths at three interfaces in soil profiles in a 19.5-ha agricultural landscape in central Pennsylvania, USA. These interfaces were (1) the interface between surface plowed layers of Ap1 and Ap2 horizons, (2) the interface with subsoil water-restricting clay layer where clay content increased to over 40%, and (3) the soil-bedrock interface. The simulated flow paths were validated through soil hydrologic monitoring, geophysical surveys, and observable soil morphological features. The results confirmed that concentrated subsurface lateral flow occurred at the interfaces with the clay layer and the underlying bedrock. At these two interfaces, the soils on the simulated flow paths were closer to saturation and showed more temporally unstable moisture dynamics than those off the simulated flow paths. Apparent electrical conductivity in the soil on the simulated flow paths was elevated and temporally unstable as compared to those outside the simulated paths. The soil cores collected from the simulated flow paths showed significantly higher Mn content at these interfaces than those away from the simulated paths. These results suggest that (1) the D8 algorithm is useful in simulating possible concentrated subsurface lateral flow paths if used with appropriate threshold value of contributing area and sufficiently detailed digital elevation model (DEM); (2) repeated electromagnetic surveys can reflect the temporal change of soil water storage and thus is a useful indicator of possible subsurface flow path over a large area; and (3) observable Mn distribution in soil profiles can be used as a simple indicator of water flow paths in soils and over the landscape; however, it does require sufficient soil sampling (by excavation or augering) to possibly infer landscape-scale subsurface flow paths. In areas where subsurface interface topography varies similarly with surface topography, surface DEM can be used to simulate potential subsurface lateral flow path reasonably so the cost associated with obtaining depth to subsurface water-restricting layer can be minimized.

A characterization of the coupled evolution of grain fabric and pore space using complex networks: Pore connectivity and optimized flows in the presence of shear bands

NASA Astrophysics Data System (ADS)

Russell, Scott; Walker, David M.; Tordesillas, Antoinette

2016-03-01

A framework for the multiscale characterization of the coupled evolution of the solid grain fabric and its associated pore space in dense granular media is developed. In this framework, a pseudo-dual graph transformation of the grain contact network produces a graph of pores which can be readily interpreted as a pore space network. Survivability, a new metric succinctly summarizing the connectivity of the solid grain and pore space networks, measures material robustness. The size distribution and the connectivity of pores can be characterized quantitatively through various network properties. Assortativity characterizes the pore space with respect to the parity of the number of particles enclosing the pore. Multiscale clusters of odd parity versus even parity contact cycles alternate spatially along the shear band: these represent, respectively, local jamming and unjamming regions that continually switch positions in time throughout the failure regime. Optimal paths, established using network shortest paths in favor of large pores, provide clues on preferential paths for interstitial matter transport. In systems with higher rolling resistance at contacts, less tortuous shortest paths thread through larger pores in shear bands. Notably the structural patterns uncovered in the pore space suggest that more robust models of interstitial pore flow through deforming granular systems require a proper consideration of the evolution of in situ shear band and fracture patterns - not just globally, but also inside these localized failure zones.

Characterizing Reactive Flow Paths in Fractured Cement

NASA Astrophysics Data System (ADS)

Wenning, Q. C.; Huerta, N. J.; Hesse, M. A.; Bryant, S. L.

2011-12-01

Geologic carbon sequestration can be a viable method for reducing anthropogenic CO2 flux into the atmosphere. However, the technology must be economically feasible and pose acceptable risk to stakeholders. One key risk is CO2 leakage out of the storage reservoir. Potential driving forces for leakage are the overpressure due to CO2 injection and the buoyancy of free phase CO2. Potential hazards of leakage are contamination of Underground Sources of Drinking Water or the atmosphere and would be deemed an unacceptable risk. Wells potentially provide a fast path for leakage from the reservoir. While the well's cement casing is reactive with CO2 and CO2-saturated brine, the low cement matrix permeability and slow diffusion rate make it unlikely that CO2 will escape through a properly constructed wellbore. However, highly permeable fractures with micrometer scale apertures can occur in cement casings. Reactions that occur in the flow in these fractures can either be self-limiting or self-enhancing. Therefore, understanding the reactive flow is critical to understanding of leakage evolution through these fractures. The goal of our work is to characterize the modification of the flow paths in the fracture due to reaction with acidic brine. With this aim we have characterized both the initial flow path of un-reactive flow and the final flow path after introduction of low-pH acid along the same fracture. Class H cement cores 3-6 cm in length and 2.5 cm diameter are created and a single natural and unique fracture is produced in each core using the Brazilian method. Our experimental fluid is injected at a constant rate into the cement core housed in a Hassler Cell under confining pressure. A solution of red dye and deionized water is pumped through the fracture to stain the un-reactive flow paths. Deionized water is then pumped through the core to limit diffusion of the dye into non-flowing portions of the fracture. After staining the initial flow path, low pH water due to hydrochloric acid (HCL), is pumped through the core at the same rate as the dye. The low pH water is used as a proxy for acidic CO2-saturated brine. Both staining from the un-reactive dye and acid produce visible permanent color alterations on the cement fracture plane. Results show that nearly the entire fracture width is stained by the red dye, with only a few asperities un-dyed. However the low pH HCl forms restricted reacted channels that are a subset of the area open to un-reactive flow, occupying only 10-50% of the entire fracture width. Low pH HCl is believed to be the driving force for the reaction that causes channeling. As acid flows through the fracture, calcium is stripped from the low pH high velocity flow front and precipitates along of the edges of the channel where pH is higher due to the lower flow velocities outside the channel. It is hypothesized that this mineral precipitation restricts the flow into localized channels within the plane of fractures having apertures of tens of micrometers. Reactions restrict the flow path to a smaller fraction of the surface, which may be an indication of self-limiting behavior.

Push-pull tracer tests: Their information content and use for characterizing non-Fickian, mobile-immobile behavior: INFORMATION CONTENT OF PUSH-PULL TESTS

DOE PAGES

Hansen, Scott K.; Berkowitz, Brian; Vesselinov, Velimir V.; ...

2016-12-01

Path reversibility and radial symmetry are often assumed in push-pull tracer test analysis. In reality, heterogeneous flow fields mean that both assumptions are idealizations. In this paper, to understand their impact, we perform a parametric study which quantifies the scattering effects of ambient flow, local-scale dispersion, and velocity field heterogeneity on push-pull breakthrough curves and compares them to the effects of mobile-immobile mass transfer (MIMT) processes including sorption and diffusion into secondary porosity. We identify specific circumstances in which MIMT overwhelmingly determines the breakthrough curve, which may then be considered uninformative about drift and local-scale dispersion. Assuming path reversibility, wemore » develop a continuous-time-random-walk-based interpretation framework which is flow-field-agnostic and well suited to quantifying MIMT. Adopting this perspective, we show that the radial flow assumption is often harmless: to the extent that solute paths are reversible, the breakthrough curve is uninformative about velocity field heterogeneity. Our interpretation method determines a mapping function (i.e., subordinator) from travel time in the absence of MIMT to travel time in its presence. A mathematical theory allowing this function to be directly “plugged into” an existing Laplace-domain transport model to incorporate MIMT is presented and demonstrated. Algorithms implementing the calibration are presented and applied to interpretation of data from a push-pull test performed in a heterogeneous environment. A successful four-parameter fit is obtained, of comparable fidelity to one obtained using a million-node 3-D numerical model. In conclusion, we demonstrate analytically and numerically how push-pull tests quantifying MIMT are sensitive to remobilization, but not immobilization, kinetics.« less

Push-pull tracer tests: Their information content and use for characterizing non-Fickian, mobile-immobile behavior: INFORMATION CONTENT OF PUSH-PULL TESTS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hansen, Scott K.; Berkowitz, Brian; Vesselinov, Velimir V.

Path reversibility and radial symmetry are often assumed in push-pull tracer test analysis. In reality, heterogeneous flow fields mean that both assumptions are idealizations. In this paper, to understand their impact, we perform a parametric study which quantifies the scattering effects of ambient flow, local-scale dispersion, and velocity field heterogeneity on push-pull breakthrough curves and compares them to the effects of mobile-immobile mass transfer (MIMT) processes including sorption and diffusion into secondary porosity. We identify specific circumstances in which MIMT overwhelmingly determines the breakthrough curve, which may then be considered uninformative about drift and local-scale dispersion. Assuming path reversibility, wemore » develop a continuous-time-random-walk-based interpretation framework which is flow-field-agnostic and well suited to quantifying MIMT. Adopting this perspective, we show that the radial flow assumption is often harmless: to the extent that solute paths are reversible, the breakthrough curve is uninformative about velocity field heterogeneity. Our interpretation method determines a mapping function (i.e., subordinator) from travel time in the absence of MIMT to travel time in its presence. A mathematical theory allowing this function to be directly “plugged into” an existing Laplace-domain transport model to incorporate MIMT is presented and demonstrated. Algorithms implementing the calibration are presented and applied to interpretation of data from a push-pull test performed in a heterogeneous environment. A successful four-parameter fit is obtained, of comparable fidelity to one obtained using a million-node 3-D numerical model. In conclusion, we demonstrate analytically and numerically how push-pull tests quantifying MIMT are sensitive to remobilization, but not immobilization, kinetics.« less

Thermal transient anemometer

DOEpatents

Bailey, J.L.; Vresk, J.

1989-07-18

A thermal transient anemometer is disclosed having a thermocouple probe which is utilized to measure the change in temperature over a period of time to provide a measure of fluid flow velocity. The thermocouple probe is located in the fluid flow path and pulsed to heat or cool the probe. The cooling of the heated probe or the heating of the cooled probe from the fluid flow over a period of time is measured to determine the fluid flow velocity. The probe is desired to be locally heated near the tip to increase the efficiency of devices incorporating the probe. 12 figs.

Local Magnetic Measurements of Trapped Flux Through a Permanent Current Path in Graphite

NASA Astrophysics Data System (ADS)

Stiller, Markus; Esquinazi, Pablo D.; Quiquia, José Barzola; Precker, Christian E.

2018-04-01

Temperature- and field-dependent measurements of the electrical resistance of different natural graphite samples suggest the existence of superconductivity at room temperature in some regions of the samples. To verify whether dissipationless electrical currents are responsible for the trapped magnetic flux inferred from electrical resistance measurements, we localized them using magnetic force microscopy on a natural graphite sample in remanent state after applying a magnetic field. The obtained evidence indicates that at room temperature a permanent current flows at the border of the trapped flux region. The current path vanishes at the same transition temperature T_c≈ 370 K as the one obtained from electrical resistance measurements on the same sample. This sudden decrease in the phase is different from what is expected for a ferromagnetic material. Time-dependent measurements of the signal show the typical behavior of flux creep of a permanent current flowing in a superconductor. The overall results support the existence of room-temperature superconductivity at certain regions in the graphite structure and indicate that magnetic force microscopy is suitable to localize them. Magnetic coupling is excluded as origin of the observed phase signal.

Geologic, hydrologic, and geochemical identification of flow paths in the Edwards Aquifer, northeastern Bexar and southern Comal Counties, Texas

USGS Publications Warehouse

Otero, Cassi L.

2007-01-01

The U.S. Geological Survey, in cooperation with the San Antonio Water System, conducted a 4-year study during 2002?06 to identify major flow paths in the Edwards aquifer in northeastern Bexar and southern Comal Counties (study area). In the study area, faulting directs ground water into three hypothesized flow paths that move water, generally, from the southwest to the northeast. These flow paths are identified as the southern Comal flow path, the central Comal flow path, and the northern Comal flow path. Statistical correlations between water levels for six observation wells and between the water levels and discharges from Comal Springs and Hueco Springs yielded evidence for the hypothesized flow paths. Strong linear correlations were evident between the datasets from wells and springs within the same flow path and the datasets from wells in areas where flow between flow paths was suspected. Geochemical data (major ions, stable isotopes, sulfur hexafluoride, and tritium and helium) were used in graphical analyses to obtain evidence of the flow path from which wells or springs derive water. Major-ion geochemistry in samples from selected wells and springs showed relatively little variation. Samples from the southern Comal flow path were characterized by relatively high sulfate and chloride concentrations, possibly indicating that the water in the flow path was mixing with small amounts of saline water from the freshwater/saline-water transition zone. Samples from the central Comal flow path yielded the most varied major-ion geochemistry of the three hypothesized flow paths. Central Comal flow path samples were characterized, in general, by high calcium concentrations and low magnesium concentrations. Samples from the northern Comal flow path were characterized by relatively low sulfate and chloride concentrations and high magnesium concentrations. The high magnesium concentrations characteristic of northern Comal flow path samples from the recharge zone in Comal County might indicate that water from the Trinity aquifer is entering the Edwards aquifer in the subsurface. A graph of the relation between the stable isotopes deuterium and delta-18 oxygen showed that, except for samples collected following an unusually intense rain storm, there was not much variation in stable isotope values among the flow paths. In the study area deuterium ranged from -36.00 to -20.89 per mil and delta-18 oxygen ranged from -6.03 to -3.70 per mil. Excluding samples collected following the intense rain storm, the deuterium range in the study area was -33.00 to -20.89 per mil and the delta-18 oxygen range was -4.60 to -3.70 per mil. Two ground-water age-dating techniques, sulfur hexafluoride concentrations and tritium/helium-3 isotope ratios, were used to compute apparent ages (time since recharge occurred) of water samples collected in the study area. In general, the apparent ages computed by the two methods do not seem to indicate direction of flow. Apparent ages computed for water samples in northeastern Bexar and southern Comal Counties do not vary greatly except for some very young water in the recharge zone in central Comal County.

Joint Use of ERT, Tracer, and Numerical Techniques to Image Preferential Flow Paths in a Fractured Granite Aquifer

NASA Astrophysics Data System (ADS)

Sanaga, S.; Vijay, S.; Kbvn, P.; Peddinti, S. R.; P S L, S.

2017-12-01

Fractured geologic media poses formidable challenges to hydrogeologists due of the strenuous mapping of fracture-matrix system and quantification of flow and transport processes. In this research, we demonstrated the efficacy of tracer-ERT studies coupled with numerical simulations to delineate preferential flow paths in a fractured granite aquifer of Deccan traps in India. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well located inside the IIT Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. Dynamic changes in sub-surface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements (R2=0.74). Fracture geometry and hydraulic properties derived from ERT and pumping tests were then used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that a dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by decrease in RMSE from 199 mg/l to 65 mg/l). A sensitivity analysis of the model parameters reveals that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations. Keywords: saline tracer, ERT, fractured granite, groundwater, preferential flow, numerical simulation

Evaluation of Flow Paths and Confluences for Saltwater Intrusion and Its Influence on Fish Species Diversity in a Deltaic River Network

NASA Astrophysics Data System (ADS)

Shao, X.; Cui, B.; Zhang, Z.; Fang, Y.; Jawitz, J. W.

2016-12-01

Freshwater in a delta is often at risk of saltwater intrusion, which has been a serious issue in estuarine deltas all over the world. Salinity gradients and hydrologic connectivity in the deltas can be disturbed by saltwater intrusion, which can fluctuate frequently and locally in time and space to affect biotic processes and then to affect the distribution patterns of the riverine fishes throughout the river network. Therefore, identifying the major flow paths or locations at risk of saltwater intrusion in estuarine ecosystems is necessary for saltwater intrusion mitigation and fish species diversity conservation. In this study, we use the betweenness centrality (BC) as the weighted attribute of the river network to identify the critical confluences and detect the preferential flow paths for saltwater intrusion through the least-cost-path algorithm from graph theory approach. Moreover, we analyse the responses of the salinity and fish species diversity to the BC values of confluences calculated in the river network. Our results show that the most likely location of saltwater intrusion is not a simple gradient change from sea to land, but closely dependent on the river segments' characteristics. In addition, a significant positive correlation between the salinity and the BC values of confluences is determined in the Pearl River Delta. Changes in the BC values of confluences can produce significant variation in the fish species diversity. Therefore, the dynamics of saltwater intrusion are a growing consideration for understanding the patterns and subsequent processes driving fish community structure. Freshwater can be diverted into these major flow paths and critical confluences to improve river network management and conservation of fish species diversity under saltwater intrusion.

Experimental investigation of turbine disk cavity aerodynamics and heat transfer

NASA Technical Reports Server (NTRS)

Daniels, W. A.; Johnson, B. V.

1993-01-01

An experimental investigation of turbine disk cavity aerodynamics and heat transfer was conducted to provide an experimental data base that can guide the aerodynamic and thermal design of turbine disks and blade attachments for flow conditions and geometries simulating those of the space shuttle main engine (SSME) turbopump drive turbines. Experiments were conducted to define the nature of the aerodynamics and heat transfer of the flow within the disk cavities and blade attachments of a large scale model simulating the SSME turbopump drive turbines. These experiments include flow between the main gas path and the disk cavities, flow within the disk cavities, and leakage flows through the blade attachments and labyrinth seals. Air was used to simulate the combustion products in the gas path. Air and carbon dioxide were used to simulate the coolants injected at three locations in the disk cavities. Trace amounts of carbon dioxide were used to determine the source of the gas at selected locations on the rotors, the cavity walls, and the interstage seal. The measurements on the rotor and stationary walls in the forward and aft cavities showed that the coolant effectiveness was 90 percent or greater when the coolant flow rate was greater than the local free disk entrainment flow rate and when room temperature air was used as both coolant and gas path fluid. When a coolant-to-gas-path density ratio of 1.51 was used in the aft cavity, the coolant effectiveness on the rotor was also 90 percent or greater at the aforementioned condition. However, the coolant concentration on the stationary wall was 60 to 80 percent at the aforementioned condition indicating a more rapid mixing of the coolant and flow through the rotor shank passages. This increased mixing rate was attributed to the destabilizing effects of the adverse density gradients.

Analysis of Lagrangian stretching in turbulent channel flow using a database task-parallel particle tracking approach

NASA Astrophysics Data System (ADS)

Meneveau, Charles; Johnson, Perry; Hamilton, Stephen; Burns, Randal

2016-11-01

An intrinsic property of turbulent flows is the exponential deformation of fluid elements along Lagrangian paths. The production of enstrophy by vorticity stretching follows from a similar mechanism in the Lagrangian view, though the alignment statistics differ and viscosity prevents unbounded growth. In this paper, the stretching properties of fluid elements and vorticity along Lagrangian paths are studied in a channel flow at Reτ = 1000 and compared with prior, known results from isotropic turbulence. To track Lagrangian paths in a public database containing Direct Numerical Simulation (DNS) results, the task-parallel approach previously employed in the isotropic database is extended to the case of flow in a bounded domain. It is shown that above 100 viscous units from the wall, stretching statistics are equal to their isotropic values, in support of the local isotropy hypothesis. Normalized by dissipation rate, the stretching in the buffer layer and below is less efficient due to less favorable alignment statistics. The Cramér function characterizing cumulative Lagrangian stretching statistics shows that overall the channel flow has about half of the stretching per unit dissipation compared with isotropic turbulence. Supported by a National Science Foundation Graduate Research Fellowship Program under Grant No. DGE-1232825, and by National Science Foundation Grants CBET-1507469, ACI-1261715, OCI-1244820 and by JHU IDIES.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Khinkis, Mark J.; Kozlov, Aleksandr P.

A radiant, non-catalytic recuperative reformer has a flue gas flow path for conducting hot exhaust gas from a thermal process and a reforming mixture flow path for conducting a reforming mixture. At least a portion of the reforming mixture flow path is positioned adjacent to the flue gas flow path to permit heat transfer from the hot exhaust gas to the reforming mixture. The reforming mixture flow path contains substantially no material commonly used as a catalyst for reforming hydrocarbon fuel (e.g., nickel oxide, platinum group elements or rhenium), but instead the reforming mixture is reformed into a higher calorificmore » fuel via reactions due to the heat transfer and residence time. In a preferred embodiment, a portion of the reforming mixture flow path is positioned outside of flue gas flow path for a relatively large residence time.« less

Moisture Risk in Unvented Attics Due to Air Leakage Paths

DOE Office of Scientific and Technical Information (OSTI.GOV)

Prahl, D.; Shaffer, M.

2014-11-01

IBACOS completed an initial analysis of moisture damage potential in an unvented attic insulated with closed-cell spray polyurethane foam. To complete this analysis, the research team collected field data, used computational fluid dynamics to quantify the airflow rates through individual airflow (crack) paths, simulated hourly flow rates through the leakage paths with CONTAM software, correlated the CONTAM flow rates with indoor humidity ratios from Building Energy Optimization software, and used Wärme und Feuchte instationär Pro two-dimensional modeling to determine the moisture content of the building materials surrounding the cracks. Given the number of simplifying assumptions and numerical models associated withmore » this analysis, the results indicate that localized damage due to high moisture content of the roof sheathing is possible under very low airflow rates. Reducing the number of assumptions and approximations through field studies and laboratory experiments would be valuable to understand the real-world moisture damage potential in unvented attics.« less

Moisture Risk in Unvented Attics Due to Air Leakage Paths

DOE Office of Scientific and Technical Information (OSTI.GOV)

Prahl, D.; Shaffer, M.

2014-11-01

IBACOS completed an initial analysis of moisture damage potential in an unvented attic insulated with closed-cell spray polyurethane foam. To complete this analysis, the research team collected field data, used computational fluid dynamics to quantify the airflow rates through individual airflow (crack) paths, simulated hourly flow rates through the leakage paths with CONTAM software, correlated the CONTAM flow rates with indoor humidity ratios from Building Energy Optimization software, and used Warme und Feuchte instationar Pro two-dimensional modeling to determine the moisture content of the building materials surrounding the cracks. Given the number of simplifying assumptions and numerical models associated withmore » this analysis, the results indicate that localized damage due to high moisture content of the roof sheathing is possible under very low airflow rates. Reducing the number of assumptions and approximations through field studies and laboratory experiments would be valuable to understand the real-world moisture damage potential in unvented attics.« less

Effects of Stochastic Traffic Flow Model on Expected System Performance

DTIC Science & Technology

2012-12-01

NSWC-PCD has made considerable improvements to their pedestrian flow modeling . In addition to the linear paths, the 2011 version now includes...using stochastic paths. 2.2 Linear Paths vs. Stochastic Paths 2.2.1 Linear Paths and Direct Maximum Pd Calculation Modeling pedestrian traffic flow...as a stochastic process begins with the linear path model . Let the detec- tion area be R x C voxels. This creates C 2 total linear paths, path(Cs

Sound source localization method in an environment with flow based on Amiet-IMACS

NASA Astrophysics Data System (ADS)

Wei, Long; Li, Min; Qin, Sheng; Fu, Qiang; Yang, Debin

2017-05-01

A sound source localization method is proposed to localize and analyze the sound source in an environment with airflow. It combines the improved mapping of acoustic correlated sources (IMACS) method and Amiet's method, and is called Amiet-IMACS. It can localize uncorrelated and correlated sound sources with airflow. To implement this approach, Amiet's method is used to correct the sound propagation path in 3D, which improves the accuracy of the array manifold matrix and decreases the position error of the localized source. Then, the mapping of acoustic correlated sources (MACS) method, which is as a high-resolution sound source localization algorithm, is improved by self-adjusting the constraint parameter at each irritation process to increase convergence speed. A sound source localization experiment using a pair of loud speakers in an anechoic wind tunnel under different flow speeds is conducted. The experiment exhibits the advantage of Amiet-IMACS in localizing a more accurate sound source position compared with implementing IMACS alone in an environment with flow. Moreover, the aerodynamic noise produced by a NASA EPPLER 862 STRUT airfoil model in airflow with a velocity of 80 m/s is localized using the proposed method, which further proves its effectiveness in a flow environment. Finally, the relationship between the source position of this airfoil model and its frequency, along with its generation mechanism, is determined and interpreted.

Quantification and Control of Wall Effects in Porous Media Experiments

NASA Astrophysics Data System (ADS)

Roth, E. J.; Mays, D. C.; Neupauer, R.; Crimaldi, J. P.

2017-12-01

Fluid flow dynamics in porous media are dominated by media heterogeneity. This heterogeneity can create preferential pathways in which local seepage velocities dwarf system seepage velocities, further complicating an already incomplete understanding of dispersive processes. In physical models of porous media flows, apparatus walls introduce preferential flow paths (i.e., wall effects) that may overwhelm other naturally occurring preferential pathways within the apparatus, leading to deceptive results. We used planar laser-induced fluorescence (PLIF) in conjunction with refractive index matched (RIM) porous media and pore fluid to observe fluid dynamics in the porous media, with particular attention to the region near the apparatus walls in a 17 cm x 8 cm x 7 cm uniform flow cell. Hexagonal close packed spheres were used to create an isotropic, homogenous porous media field in the interior of the apparatus. Visualization of the movement of a fluorescent dye revealed the influence of the wall in creating higher permeability preferential flow paths in an otherwise homogenous media packing. These preferential flow paths extended approximately one half of one sphere diameter from the wall for homogenously packed regions, with a quickly diminishing effect on flow dynamics for homogenous media adjacent to the preferential pathway, but with major influence on flow dynamics for adjoining heterogeneous regions. Multiple approaches to mitigate wall effects were investigated, and a modified wall was created such that the fluid dynamics near the wall mimics the fluid dynamics within the homogenous porous media. This research supports the design of a two-dimensional experimental apparatus that will simulate engineered pumping schemes for use in contaminant remediation. However, this research could benefit the design of fixed bed reactors or other engineering challenges in which vessel walls contribute to unwanted preferential flow.

Extensive electron transport and energization via multiple, localized dipolarizing flux bundles

NASA Astrophysics Data System (ADS)

Gabrielse, Christine; Angelopoulos, Vassilis; Harris, Camilla; Artemyev, Anton; Kepko, Larry; Runov, Andrei

2017-05-01

Using an analytical model of multiple dipolarizing flux bundles (DFBs) embedded in earthward traveling bursty bulk flows, we demonstrate how equatorially mirroring electrons can travel long distances and gain hundreds of keV from betatron acceleration. The model parameters are constrained by four Time History of Events and Macroscale Interactions during Substorms satellite observations, putting limits on the DFBs' speed, location, and magnetic and electric field magnitudes. We find that the sharp, localized peaks in magnetic field have such strong spatial gradients that energetic electrons ∇B drift in closed paths around the peaks as those peaks travel earthward. This is understood in terms of the third adiabatic invariant, which remains constant when the field changes on timescales longer than the electron's drift timescale: An energetic electron encircles a sharp peak in magnetic field in a closed path subtending an area of approximately constant flux. As the flux bundle magnetic field increases the electron's drift path area shrinks and the electron is prevented from escaping to the ambient plasma sheet, while it continues to gain energy via betatron acceleration. When the flux bundles arrive at and merge with the inner magnetosphere, where the background field is strong, the electrons suddenly gain access to previously closed drift paths around the Earth. DFBs are therefore instrumental in transporting and energizing energetic electrons over long distances along the magnetotail, bringing them to the inner magnetosphere and energizing them by hundreds of keV.