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.

Analytic Benchmark Solutions for Open-Channel Flows

An important test of the quality of a computational model is its ability to reproduce standard test cases or benchmarks. For steady open–channel flow based on the Saint Venant equations some benchmarks exist for simple geometries from the work of Bresse, Bakhmeteff and Chow but these are tabulated in the form of standard integrals. This paper provides benchmark solutions for a wider range of cases, which may have a nonprismatic cross section, nonuniform bed slope, and transitions between subcritical and supercritical flow. This makes it possible to assess the underlying quality of computational algorithms in more difficult cases, including those with hydraulic jumps. Several new test cases are given in detail and the performance of a commercial steady flow package is evaluated against two of them. The test cases may also be used as benchmarks for both steady flow models and unsteady flow models in the steady limit.

Upwind solution of singular differential equations arising from steady channel flows

We study ordinary nonlinear singular differential equations which arise from steady conservation laws with source terms. An example of steady conservation laws which leads to those scalar equations is the Saint–Venant equations. The numerical solution of these scalar equations is sought by using the ideas of upwinding and discretisation of source terms. Both the Engquist–Osher scheme and the Roe scheme are used with different strategies for discretising the source terms.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

Extraction of tidal channel networks from airborne scanning laser altimetry and aerial photography

The study of the morphodynamics of tidal channel networks is important because of their role in tidal propagation and the evolution of salt-marshes and tidal flats. Channel dimensions range from tens of metres wide and metres deep near the low water mark to only 20-30cm wide and 20cm deep for the smallest channels on the marshes. The conventional method of measuring the networks is cumbersome, involving manual digitising of aerial photographs. This paper describes a semi-automatic knowledge-based network extraction method that is being implemented to work using airborne scanning laser altimetry (and later aerial photography). The channels exhibit a width variation of several orders of magnitude, making an approach based on multi-scale line detection difficult. The processing therefore uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels using a distance-with-destination transform. Breaks in the networks are repaired by extending channel ends in the direction of their ends to join with nearby channels, using domain knowledge that flow paths should proceed downhill and that any network fragment should be joined to a nearby fragment so as to connect eventually to the open sea.

Inhibition of synaptic transmission and G protein modulation by synthetic CaV2.2 Ca2+ channel peptides

Abstract: Modulation of presynaptic voltage-dependent Ca+ channels is a major means of controlling neurotransmitter release. The CaV 2.2 Ca2+ channel subunit contains several inhibitory interaction sites for Gβγ subunits, including the amino terminal (NT) and I–II loop. The NT and I–II loop have also been proposed to undergo a G protein-gated inhibitory interaction, whilst the NT itself has also been proposed to suppress CaV 2 channel activity. Here, we investigate the effects of an amino terminal (CaV 2.2[45–55]) ‘NT peptide’ and a I–II loop alpha interaction domain (CaV 2.2[377–393]) ‘AID peptide’ on synaptic transmission, Ca2+ channel activity and G protein modulation in superior cervical ganglion neurones (SCGNs). Presynaptic injection of NT or AID peptide into SCGN synapses inhibited synaptic transmission and also attenuated noradrenaline-induced G protein modulation. In isolated SCGNs, NT and AID peptides reduced whole-cell Ca2+ current amplitude, modified voltage dependence of Ca2+ channel activation and attenuated noradrenaline-induced G protein modulation. Co-application of NT and AID peptide negated inhibitory actions. Together, these data favour direct peptide interaction with presynaptic Ca2+ channels, with effects on current amplitude and gating representing likely mechanisms responsible for inhibition of synaptic transmission. Mutations to residues reported as determinants of Ca2+ channel function within the NT peptide negated inhibitory effects on synaptic transmission, Ca2+ current amplitude and gating and G protein modulation. A mutation within the proposed QXXER motif for G protein modulation did not abolish inhibitory effects of the AID peptide. This study suggests that the CaV 2.2 amino terminal and I–II loop contribute molecular determinants for Ca2+ channel function; the data favour a direct interaction of peptides with Ca2+ channels to inhibit synaptic transmission and attenuate G protein modulation. Non-technical summary: Nerve cells (neurones) in the body communicate with each other by releasing chemicals (neurotransmitters) which act on proteins called receptors. An important group of receptors (called G protein coupled receptors, GPCRs) regulate the release of neurotransmitters by an action on the ion channels that let calcium into the cell. Here, we show for the first time that small peptides based on specific regions of calcium ion channels involved in GPCR signalling can themselves inhibit nerve cell communication. We show that these peptides act directly on calcium channels to make them more difficult to open and thus reduce calcium influx into native neurones. These peptides also reduce GPCR-mediated signalling. This work is important in increasing our knowledge about modulation of the calcium ion channel protein; such knowledge may help in the development of drugs to prevent signalling in pathways such as those involved in pain perception.

Natural ventilation assessment in typical open an semi-open urban environments under various wind directions

Semi-open street roofs protect pedestrians from intense sunshine and rains. Their effects on natural ventilation of urban canopy layers (UCL) are less understood. This paper investigates two idealized urban models consisting of 4(2×2) or 16(4×4) buildings under a neutral atmospheric condition with parallel (0°) or non-parallel (15°,30°,45°) approaching wind. The aspect ratio (building height (H) / street width (W)) is 1 and building width is B=3H. Computational fluid dynamic (CFD) simulations were first validated by experimental data, confirming that standard k-ε model predicted airflow velocity better than RNG k-ε model, realizable k–ε model and Reynolds stress model. Three ventilation indices were numerically analyzed for ventilation assessment, including flow rates across street roofs and openings to show the mechanisms of air exchange, age of air to display how long external air reaches a place after entering UCL, and purging flow rate to quantify the net UCL ventilation capacity induced by mean flows and turbulence. Five semi-open roof types are studied: Walls being hung above street roofs (coverage ratio λa=100%) at z=1.5H, 1.2H, 1.1H ('Hung1.5H', 'Hung1.2H', 'Hung1.1H' types); Walls partly covering street roofs (λa=80%) at z=H ('Partly-covered' type); Walls fully covering street roofs (λa=100%) at z=H ('Fully-covered' type).They basically obtain worse UCL ventilation than open street roof type due to the decreased roof ventilation. 'Hung1.1H', 'Hung1.2H', 'Hung1.5H' types are better designs than 'Fully-covered' and 'Partly-covered' types. Greater urban size contains larger UCL volume and requires longer time to ventilate. The methodologies and ventilation indices are confirmed effective to quantify UCL ventilation.

Identifying the Open-Closed Field Line Boundary

Of all the various definitions of the polar cap boundary that have been used in the past, the most physically meaningful and significant is the boundary between open and closed field lines. Locating this boundary is very important as it defines which regions and phenomena are on open field lines and which are on closed. This usually has fundamental implications for the mechanisms invoked. Unfortunately, the open-closed boundary is usually very difficult to identify, particularly where it maps to an active reconnection site. This paper looks at the topological reconnection classes that can take place, both at the magnetopause and in the cross-tail current sheet and discusses the implications for identifying the open-closed boundary when reconnection is giving velocity filter dispersion of signatures. On the dayside, it is shown that the dayside boundary plasma sheet and low-latitude boundary layer precipitations are well explained as being on open field lines, energetic ions being present because of reflection of central plasma sheet ions off the two Alfvén waves launched by the reconnection site (the outer one of which is the magnetopause). This also explains otherwise anomalous features of the dayside convection pattern in the cusp region. On the nightside, similar considerations place the open-closed boundary somewhat poleward of the velocity-dispersed ion structures which are a signature of the plasma sheet boundary layer ion flows in the tail.