An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Hydrology and geomorphology of the Upper White Nile lakes and their relevance for water resources management in the Nile basin

Remote sensing data and digital elevation models were utilized to extract the catchment hydrological parameters and to delineate storage areas for the Ugandan Equatorial Lakes region. Available rainfall/discharge data are integrated with these morphometric data to construct a hydrological model that simulates the water balance of the different interconnected basins and enables the impact of potential management options to be examined. The total annual discharges of the basins are generally very low (less than 7% of the total annual rainfall). The basin of the shallow (5 m deep) Lake Kioga makes only a minor hydrological contribution compared with other Equatorial Lakes, because most of the overflow from Lake Victoria basin into Lake Kioga is lost by evaporation and evapotranspiration. The discharge from Lake Kioga could be significantly increased by draining the swamps through dredging and deepening certain channel reaches. Development of hydropower dams on the Equatorial Lakes will have an adverse impact on the annual water discharge downstream, including the occasional reduction of flow required for filling up to designed storage capacities and permanently increasing the surface areas of water that is exposed to evaporation. On the basis of modelling studies, alternative sites are proposed for hydropower development and water storage schemes

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

Two new species of Loricariidae (Teleostei: Siluriformes) from main channels of the upper and middle Amazon Basin, with discussion of deep water specialization in loricariids

Hemiancistrus pankimpuju, new species, and Panaque bathyphilus, new species, are described from the main channel of the upper (Maranon) and middle (Solimoes)Amazon River, respectively. Both species are diagnosed by having a nearly white body, long filamentous extensions of both simple caudal-fin rays, small eyes, absence of an iris operculum and unique combinations of morphometrics. The coloration and morphology of these species, unique within Loricariidae, are hypothesized to be apomorphies associated with life in the dark, turbid depths of the Amazon mainstem. Extreme elongation of the caudal filaments in these and a variety of other main channel catfishes is hypothesized to have a mechanosensory function associated with predator detection.

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide`s primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

Combined Vertical And Lateral Channel Evolution Numerical Modeling

Vertical stream bed erosion has been studied routinely and its modeling is getting widespread acceptance. The same cannot be said with lateral stream bank erosion since its measurement or numerical modeling is very challenging. Bank erosion, however, can be important to channel morphology. It may contribute significantly to the overall sediment budget of a stream, is a leading cause of channel migration, and is the cause of major channel maintenance. However, combined vertical and lateral channel evolution is seldom addressed. In this study, a new geofluival numerical model is developed to simulate combined vertical and lateral channel evolution. Vertical erosion is predicted with a 2D depth-averaged model SRH-2D, while lateral erosion is simulated with a linear retreat bank erosion model developed in this study. SRH-2D and the bank erosion model are coupled together both spatially and temporally through a common mesh and the same time advancement. The new geofluvial model is first tested and verified using laboratory meander channels; good agreement are obtained between predicted bank retreat and measured data. The model is then applied to a 16-kilometer reach of Chosui River, Taiwan. Vertical and lateral channel evolution during a three-year period (2004 to 2007) is simulated and results are compared with the field data. It is shown that the geofluvial model correctly captures all major erosion and deposition patterns. The new model is shown to be useful for identifying potential erosion sites and providing information for river maintenance planning.

Model Configuration And Data Management In The Short-Term Water Information Forecasting Tools

The Short-term Water Information and Forecasting Tools (SWIFT) is a suite of tools for flood and short-term streamflow forecasting, consisting of a collection of hydrologic model components and utilities. Catchments are modeled using conceptual subareas and a node-link structure for channel routing. The tools comprise modules for calibration, model state updating, output error correction, ensemble runs and data assimilation. Given the combinatorial nature of the modelling experiments and the sub-daily time steps typically used for simulations, the volume of model configurations and time series data is substantial and its management is not trivial. SWIFT is currently used mostly for research purposes but has also been used operationally, with intersecting but significantly different requirements. Early versions of SWIFT used mostly ad-hoc text files handled via Fortran code, with limited use of netCDF for time series data. The configuration and data handling modules have since been redesigned. The model configuration now follows a design where the data model is decoupled from the on-disk persistence mechanism. For research purposes the preferred on-disk format is JSON, to leverage numerous software libraries in a variety of languages, while retaining the legacy option of custom tab-separated text formats when it is a preferred access arrangement for the researcher. By decoupling data model and data persistence, it is much easier to interchangeably use for instance relational databases to provide stricter provenance and audit trail capabilities in an operational flood forecasting context. For the time series data, given the volume and required throughput, text based formats are usually inadequate. A schema derived from CF conventions has been designed to efficiently handle time series for SWIFT.

A Study On The Flow Characteristics Influenced By Hydraulic Structure At A Channel Junction

Natural riversare consisting of various networks as junction andstreams. And sediment and erosion are occurred by specific stream condition. When flood season,large discharge flew in the river and river bed changed by high flow velocity. Especially junction area’s flow characteristics are very complex. The purpose of this study is to analyze the flow characteristics in channel junction, which are most influenced by large discharge like flooding and input water from tributary. We investigate the flow characteristics by using hydrodynamics and transport module in MIKE 3 FM. MIKE 3 FM model was helpful tool to analysis 3D hydrodynamics, erosion and sediment effect from channel bed. We analyze flow characteristics at channel junction. Also we consider hydraulic structures like a bridge pier which is influencing flow characteristics like a flow velocity, water level, erosion and scour depth in channel bed. In the model, we controlled discharge condition according to Froude Number and reflect various grain diameter size and flow ratio change in main stream and tributary. In the result, flow velocity, water level, erosion and sediment depth are analyzed. Additionally, we suggest a these result relationship with equations. This study will help the understand flow characteristics and influence of hydraulic structure in channel junction. Acknowledgments This research was supported by a grant (12-TI-C01) from Advanced Water Management Research Program funded by Ministry of Land, Infrastructure and Transport of Korean government.

L-Type calcium channels mediate water intake induced by angiotensin injection into median preoptic nucleus

Calcium ions are widely accepted as critically important in responses of neurons to a stimulus. We have show previously the central involvement of angiotensin II (ANGII) in water intake. This study determined whether voltage-dependent calcium channels are involved in ANGII-induced behavioral drinking implicating nitrergic mechanism. The antidipsogenic actions of L-type calcium channel antagonists nifedipine, on ANGII-induced drinking behavior were studied when it is injected into the median preoptic nucleus (MnPO). The influence of nitric oxide (NO) on nifedipine antidipsogenic action was also studied by utilizing the N-W-nitro-L-arginine methyl ester (L-NAME) a constitutive nitric oxide synthase inhibitor constitutive (cNOSI) and 7-nitroindazol (7-NIT) a specific neuronal nitric oxide synthase inhibitor (nNOSI) and L-arginine a NO donor. Rats 200-250 g, with cannulae implanted into MnPO, pre-treated into MnPO with either nifedipine, followed by ANGII, drank significantly less water than controls during the first 15 min after injection. However, L-NAME potentiated the dipsogenic effect of ANGII that is blocked by prior injection of nifedipine and L-arginine. 7-NIT injected prior to ANGII into MnPO also potentiated the dipsogenic effect of ANGII but with a less intensity than L-NAME that it is also blocked by prior injection of nifedipine. The results described in this paper provide evidence that calcium channels play important roles in the ANGII-induced behavioral water intake. The structures containing NO in the brain such as MnPO include both endothelial cells and neurons might be responsible for the influence of nifedipine on dipsogenic effect of ANGII. These data shows the correlation between L-type calcium channel and a free radical gas NO produced endogenously from amino acids L-arginine by endothelial and neuronal NO synthase in the control of ANGII-dipsogenic effect. This suggests that an L-type calcium channel participates in both short- and longer-term neuronal actions of ANGII by nitrergic way. (c) 2006 Elsevier B.V. All rights reserved.

On a possible dual role for the lateral septal area 5-HT1A receptor system in the regulation of water intake and urinary excretion

The 5-hydroxytryptamine (5-HT)(1A) receptor system plays a prominent role in a variety of physiological functions and behavior and regulation of this responsiveness of the receptor system has been implicated in the central regulation of water intake and urinary excretion. The lateral septal area (LSA) exhibits a high density of 5-HT1A receptors, as well as a subpopulation of oxytocin (OT) receptors. Here we report the effects of pMPPF (a selective 5-HT1A antagonist), d(CH2)(5)[Tyr(Me)(2)Thr(4), Orn(5), Tyr(NH2)(9)]-vasotocin (an OT antagonist), and that 5-HT1A receptor system is regulated as a consequence of activation of the Na+ channel by veratridine. Cannulae were implanted into the LSA of rats to enable the introduction of the drugs. Injections of 8-OH-DPAT (a 5-HT1A agonist) blocked water intake and increased urinary excretion, while pMPPF or the OT antagonist injected bilaterally before 8-OH-DPAT blocked its inhibitory effect on water intake and its diuretic effect. In contrast, increases in extracellular sodium levels induced by the sodium channel modulator, veratridine, enhanced 5-HT1A responsiveness for water intake and reduced the diuretic effects induced by 8-OH-DPAT. These trials demonstrated that the responsiveness of the 5-HT1A receptor system in the LSA can be enhanced or depressed as a consequence of an induced rise in extracellular sodium. (C) 2010 Elsevier B.V. All rights reserved.

Changes in abiotic characteristics of water in the Paranapanema River and three lateral lagoons at mouth zone of the Jurumirim Reservoir during the flood period, São Paulo, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide's primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.