Treatment of Natural Geometry in Finite Volume River Flow Computations

A method is proposed for the treatment of irregular bathymetry in one-dimensional finite volume computations of open-channel flow. The strategy adopted is based on a reformulation of the Saint-Venant equations. In contrast with the usual treatment of topography effects as source terms, the method accounts for slope and nonprismaticity by modifying the momentum flux. This makes it possible to precisely balance the hydrostatic pressure contributions associated with variations in valley geometry. The characteristic method is applied to the revised equations, yielding topographic corrections to the numerical fluxes of an upwind scheme. Further adaptations endow the scheme with an ability to capture transcritical sections and wetting fronts in channels of abrupt topography. To test the approach, the scheme is first applied to idealized benchmark problems. The method is then used to route a severe flood through a complex river system: the Tanshui in Northern Taiwan. Computational results compare favorably with gauge records. Discrepancies in water stage represent no more than a fraction of the magnitude of typical bathymetry variations.     

Total Load Transport Formula for Flow in Alluvial Channels

A user-friendly total bed-material load transport formula for flow in alluvial channels under equilibrium transport conditions has been developed based on dimensional analysis. The main advantages of this formula are its ease of computation, accuracy in prediction, and the wide range of application. The total sediment discharge gt is computed directly and is linearly related to the new total load transport parameter, TT. The latter involves variables that can be easily measured in field conditions, i.e., flow depth, mean flow velocity, energy slope, median sediment size and density, and water temperature. The factor of proportionality k in the formula has been checked for a wide range of hydraulic conditions and it remains a constant equal to 12.5. Comparisons between the computed and measured total sediment discharge indicate that the predictions are good.     

Discontinuous Galerkin Finite-Element Method for Transcritical Two-Dimensional Shallow Water Flows

A total variation diminishing Runge Kutta discontinuous Galerkin finite-element method for two-dimensional depth-averaged shallow water equations has been developed. The scheme is well suited to handle complicated geometries and requires a simple treatment of boundary conditions and source terms to obtain high-order accuracy. The explicit time integration, together with the use of orthogonal shape functions, makes the method for the investigated flows computationally as efficient as comparable finite-volume schemes. For smooth parts of the solution, the scheme is second order for linear elements and third order for quadratic shape functions both in time and space. Shocks are usually captured within only two elements. Several steady transcritical and transient flows are investigated to confirm the accuracy and convergence of the scheme. The results show excellent agreement with analytical solutions. For investigating a flume experiment of supercritical open-channel flow, the method allows very good decoupling of the numerical and mathematical model, resulting in a nearly grid-independent solution. The simulation of an actual dam break shows the applicability of the scheme to nontrivial bathymetry and wave propagation on a dry bed.     

Effect of Stilling Basin Geometry on Clear Water Scour Morphology Downstream of a Block Ramp

The scour mechanism downstream of a block ramp in clear water conditions is quite a complex phenomenon that depends on several parameters. Majors role are played by ramp configuration, hydraulic conditions (downstream tailwater level and discharge), material granulometry, and stilling basin geometry. Previous studies analyzed both the scour phenomenon and the effects of all the parameters involved except the last one (i.e., the case in which the ramp has the same width as the downstream stilling basin) and the case of symmetrically expanding stilling basins. This last basin configuration represents an optimal equilibrium between the necessity to dissipate energy and to create a natural pool for the biological species, thus it has a prominent ecological value. The present paper aims to assess the effect of both the width and length of the downstream stilling basin on scour features and flow pattern in clear water conditions. Different scour morphology types are distinguished and classified according to hydraulic and geometric conditions. Simple novel relationships are proposed to evaluate the scour depth and length and the maximum water level in the stilling basin. The results are valid for unsubmerged block ramps.     

Preserving Steady-State in One-Dimensional Finite-Volume Computations of River Flow

When using finite-volume methods and the conservative form of the Saint Venant equations in one-dimensional flow computations, it is important to establish the correct balance between the discretized flux vector and the geometric source terms. Over the last few years various improvements to numerical schemes have been presented to achieve this correct balance, focusing on the capability to simulate water at rest on irregular geometries (C-property). In this paper it is shown that common schemes can lead to energy-violating solutions in the case of steady flow. We present developments based on the Roe TVD finite-volume scheme for one-dimensional Saint Venant equations, which results in a method that not only satisfies the C-property, but also preserves the correct steady flow when stationary boundary conditions are used. We also present a totally irregular channel test case for the verification of the method.     

Well-Balanced Scheme between Flux and Source Terms for Computation of Shallow-Water Equations over Irregular Bathymetry

Although many numerical techniques such as approximate Riemann solvers can be used to analyze subcritical and supercritical flows modeled by hyperbolic-type shallow-water equations, there are some difficulties in practical applications due to the numerical unbalance between source and flux terms. In this study, a revised surface gradient method is proposed that balances source and flux terms. The new numerical model employs the MUSCL–Hancock scheme and the HLLC approximate Riemann solver. Several verifications are conducted, including analyses of transcritical steady-state flows, unsteady dam break flows on a wet and dry bed, and flows over an irregular bathymetry. The model consistently returns accurate and reasonable results comparable to those obtained through analytical methods and laboratory experiments. The revised surface gradient method may be a simple but robust numerical scheme appropriate for solving hyperbolic-type shallow-water equations over an irregular bathymetry.     

Flow Computation and Mass Balance in Galerkin Finite-Element Groundwater Models

In most groundwater modeling studies, quantification of the flow rates at domain and subdomain boundaries is as important as the computation of the groundwater heads. The computation of these flow rates is not a trivial task when a finite-element method is chosen to solve the groundwater equation. Generally, it is believed that finite-element methods do not conserve mass locally. In this paper, a postprocessing technique is developed to compute mass-conserving flow rates at element faces. It postprocesses the groundwater head field obtained by the Galerkin finite-element method, and the calculated flow rates conserve mass locally and globally. The only requirement for the postprocessor to be applicable is the irrotationality of the flow field, i.e., the curl of the Darcy flux should be zero. The accuracy and the mass conservation properties of the new postprocessor are demonstrated using several test problems that include one-, two-, and three-dimensional flow systems in both homogeneous and heterogeneous aquifer conditions.     

Obsessive-compulsive disorder among schizophrenic patients: an exploratory study using functional magnetic resonance imaging data

The case of a 70-year-old woman with cerebral amyloid angiopathy (CAA) is presented. MRI of the head showed widespread miliary foci of haemorrhage within the cerebrum and cerebellum, with some additional linear lesions within the cerebral cortex and patchy lesions in the white matter. This is in contrast to the more usual pattern of intracranial haemorrhage in CAA, i.e., a lobar haematoma.     

Simple Optimal Downstream Feedback Canal Controllers: ASCE Test Case Results

In a companion paper, a class of downstream-water-level feedback canal controllers was described. Within this class, a particular controller is chosen by selecting which controller coefficients to optimize (tune), the remaining coefficients being set to zero. These controllers range from a series of simple proportional-integral (PI) controllers to a single centralized controller that considers lag times. In this paper, several controllers within this class were tuned with the same quadratic performance criteria (i.e., identical penalty functions for optimization). The resulting controllers were then tested through unsteady-flow simulation with the ASCE canal automation test cases for canal 1. Differences between canal and gate properties, as simulated and as assumed for tuning, reduced controller performance in terms of both water-level errors and gate movements. The test case restrictions placed on minimum gate movement caused water levels to oscillate around their set points. This resulted in steady-state errors and much more gate movement (hunting). More centralized controllers handle unscheduled flow changes better than a series of local PI controllers. Controllers that explicitly account for pool wave travel times did not improve control as much as expected. Sending control actions within a given pool to upstream pools improved performance, but caused oscillations in some cases, unless control signals were also sent downstream. A good compromise between controller performance and complexity is provided by controllers that pass feedback from a given water level to the check structure at the upstream end of its pool (i.e., that used for downstream control of an individual pool) and to all upstream and one downstream check structures.     

In disperse solution, "osmotic stress" is a restricted case of preferential interactions

In the practice of "osmotic stress," the effect of excluded cosolvents on a biochemical equilibrium is interpreted as the number of water molecules participating in the reaction. This action is attributed to lowering of solvent water activity by the cosolvent. This concept of osmotic stress in disperse solution is erroneous: (i) A cosolvent cannot be both excluded and inert, i.e., noninteracting, because exclusion requires a positive free energy change; (ii) a decrease in water activity alone by addition of solute cannot affect an equilibrium when the reacting surface is in contact with the solvent; and (iii) osmotic stress in disperse solution is a restricted case of preferential interactions; the reaction is driven by the free energy of cosolvent exclusion, and the derived number of water molecules is solely a measure of the mutual perturbations of the chemical potentials of the cosolvent and the protein.     

Intrusion within a Simulated Water Distribution System due to Hydraulic Transients. I: Description of Test Rig and Chemical Tracer Method

A pilot-scale test rig was designed and constructed to simulate intrusion behavior associated with hydraulic transients initiated by sudden events such as rapid valve closure or uncontrolled change in on/off pump status (“pump trip”) in a water distribution system. After establishing steady state flow conditions, intrusion volumes were determined for 3.2 (1/8-in.) and 6.4 mm (1/4-in.) diam orifices overlaid with a column of water that provided 91 mm (3 ft) of external head. Average intrusion volumes associated with hydraulic transient events were determined by mass balance calculations using cesium as the tracer chemical. Average intrusion volumes were 11.4 and 71.2 mL for the two diameters, respectively. Given similar conditions in a water distribution system (i.e., an available pathway and favorable pressures), pathogens in the soil and water surrounding a water main potentially can intrude into the pipe during short-term pressure transient events.     

Effect of ecological viewing conditions on the Ames' distorted room illusion.

Ecological theory asserts that the Ames' distorted room illusion (DRI) occurs as a result of the artificial restriction of information pickup. According to J. J. Gibson (1966, 1979), the illusion is eliminated when binocular vision and/or head movement are allowed. Exp I, with 144 undergraduates, used size-matching technique employing discs placed within an Ames' distorted room to measure the DRI. Ss viewed the distorted room or a control apparatus under 4 viewing conditions (i.e., restricted or unrestricted head movement), using monocular and binocular vision. In Exp II, 20 Ss viewed binocularly and were instructed to move freely while making judgments. Findings show that the DRI decreased with increases in viewing access, although it persisted under all viewing conditions. The persistence of the illusion is seen as contradicting Gibson's position. (12 ref) (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

Experimental Evidence of Hysteresis in the Head-Discharge Relationship for a Leak in a Polyethylene Pipe

The relationship between leak outflow from a damaged pipe and flow condition inside the pipe plays a crucial role in pressurized pipe systems management. As an example, this relationship is used in leakage reduction techniques based on pressure control and in leak detection techniques based on inverse analysis. To explore the relationship between total head inside the pipe and leak outflow for a single leak in a polyethylene pipe, tests were carried out at the Water Engineering Laboratory of the University of Perugia. These tests point out that the viscoelastic nature of the pipe material gives rise to a hysteretic behavior of the investigated relationship, i.e., the outflow depends not only on the synchronous total head but also on the total head time history and variation rate.     

Blocking of subsequent and antecedent events.

Stimulus competition (e.g., blocking) has been observed between antecedent events (i.e., conditioned stimuli or potential causes), but recent evidence within the human causal learning literature suggests that it could also be obtained between subsequent events (i.e., unconditioned stimuli or potential effects). The present research tested this hypothesis with rat subjects. to avoid confounding the antecedent versus subsequent variable with the affective value of the events involved (i.e., unconditioned stimuli are ordinarily of greater affective value than conditioned stimuli), a preparation was used in which antecedent and subsequent events all lacked affective value during the blocking phases of the study. This was achieved through the use of sensory preconditioning. Blocking of subsequent events as well as antecedent events were observed. The challenge to most associative theories that is provided by blocking of subsequent events is discussed. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

Upwind Conservative Scheme for the Saint Venant Equations

An upwind conservative scheme with a weighted average water-surface-gradient approach is proposed to compute one-dimensional open channel flows. The numerical scheme is based on the control volume method. The intercell flux is computed by the one-sided upwind method. The water surface gradient is evaluated by the weighted average of both upwind and downwind gradients. The scheme is tested with various examples, including dam-break problems in channels with rectangular and triangular cross-sections, hydraulic jump, partial dam-break problem, overtopping flow, a steady flow over bump with hydraulic jump, and a dam-break flood case in a natural river valley. Comparisons between numerical and exact solutions or experimental data demonstrated that the proposed scheme is capable of accurately reproducing various open channel flows, including subcritical, supercritical, and transcritical flows. The scheme is inherently robust, stable, and monotone. The scheme does not require any special treatment, such as artificial viscosity or front tracking technique, to capture steep gradients or discontinuities in the solution.     

Cartesian Cut Cell Two-Fluid Solver for Hydraulic Flow Problems

A two-fluid solver which can be applied to a variety of hydraulic flow problems has been developed. The scheme is based on the solution of the incompressible Euler equations for a variable density fluid system using the artificial compressibility method. The computational domain encompasses both water and air regions and the interface between the two fluids is treated as a contact discontinuity in the density field which is captured automatically as part of the solution using a high resolution Godunov-type scheme. A time-accurate solution has been achieved by using an implicit dual-time iteration technique. The complex geometry of the solid boundary arising in the real flow problems is represented using a novel Cartesian cut cell technique, which provides a boundary fitted mesh without the need for traditional mesh generation techniques. A number of test cases including the classical low amplitude sloshing tank and dam-break problems, as well as a collapsing water column hitting a downstream obstacle have been calculated using the present approach and the results compare very well with other theoretical and experimental results. Finally, a test case involving regular waves interacting with a sloping beach is also calculated to demonstrate the applicability of the method to real hydraulic problems.     

Flow of Steel in Mold Region During Continuous Casting

 The particle image velocimetry (PIV) technique was used to study the fluid flow phenomena that occurred during continuous casting, using a water model with dimensions of 1 840 mm×280 mm. Two types of solidified shells, ie, the smooth type and the coarse type, were used to characterize the dendrite in order to simulate different liquid solid interfacial conditions. The influence of the nozzle angle and the immersion depth of nozzle, as well as the casting speed on the flow behavior was investigated quantitatively. The results were as follows: (1) There are two large recirculations above and below the fluid jet in the mold, respectively, under the smooth interface condition. However, in the case of the dendrite solidified shell, it was found that the flow velocity of the fluid decreased and more smaller vortices appeared in the upper region of the mold. (2) The angle and the immersion depth of nozzle are two important factors affecting the flow pattern, and they are also capable of bringing about the change in the flow direction. (3) The higher the casting speed, the higher are the jet stream and the impacting point on the narrow face. However, the high casting speed causes serious fluctuation of the meniscus, and correspondingly leads to various defects.     

In vitro training model for endovascular embolization of cerebral aneurysms

We present a new in vitro model for the training of endovascular embolizations of cerebral aneurysms. The technique and laboratory set-up described allow the embolization of thin-walled, transparent silicone aneurysm models under "semi-physiologic" conditions, i.e., in a closed circulation with pulsatile flow and pressure, simulating the hemodynamic conditions encountered in biologic systems.     

Scouring at Bed Sills as a Response to Flash Floods

The temporal development of clear-water local scour depth at bed sills in uniform gravel beds is considered. Experiments are presented on the development of scour holes under unsteady hydraulic conditions, with the triangular-shaped hydrographs tested being of different durations and different rates of flow variation. Based on the experimental results and a theoretical framework, a method is given for the definition and prediction of the scouring process under unsteady flows in terms of a dimensionless temporal parameter. A “flash flood” is here defined as an event for which the scour doesn’t attain its potential magnitude, i.e., the equilibrium value for the peak hydrograph flow rate. This flood nature is dependent on both the characteristics of the flood event itself and the characteristics of the stream. A quantitative measure of what constitutes a flash flood is given in terms of the identified temporal parameter. Results show that the ratio between the final scour depth and the potential scour depth at a bed sill for a given hydrograph can be estimated as a function of the identified temporal parameter.     

Application of Three-Dimensional Hydrodynamic Model for Lake Okeechobee

A 3D hydrodynamic and heat transport model was developed for Lake Okeechobee. Continuity, momentum, and temperature transport equations were solved. Dynamically coupled transport equations for turbulent kinetic energy and turbulent scale also were solved. The numerical scheme used spatial finite differencing and a three-time-level, external-internal mode splitting procedure. A 28-day calibration was conducted, using measured bathymetry, rainfall, relative humidity, total solar radiation, wind velocity, inflow, and outflow data. During the calibration period, little rainfall occurred, and lake water levels receded. Water surface elevation, horizontal velocities, and temperature were computed. Agreement between observed and simulated values was based on graphical comparisons, minimizing mean absolute and root-mean-square errors, and spectral analysis. Comparisons showed that the model reproduced general observed trends and short-term fluctuations. The model's heat transport and turbulence closure schemes behaved as expected with regard to water column stratification and mixing. Simulation accuracy may potentially be improved by adding wind-wave and vegetation resistance algorithms to the model.